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.

Diel variations in the estimated refractive index of bulk oceanic particles.

The index of refraction (n) of particles is an important parameter in optical models that aims to extract particle size and carbon concentrations from light scattering measurements. An inadequate choice of n can critically affect the characterization and interpretation of optically-derived parameters, including those from satellite-based models which provide the current view of how biogeochemical processes vary over the global ocean. Yet, little is known about how n varies over time and space to inform such models. Particularly, in situ estimates of n for bulk water samples and at diel-resolving time scales are rare. Here, we demonstrate a method to estimate n using simultaneously and independently collected particulate beam attenuation coefficients, particle size distribution data, and a Mie theory model. We apply this method to surface waters of the North Pacific Subtropical Gyre (NPSG) at hourly resolution. Clear diel cycles in n were observed, marked by minima around local sunrise and maxima around sunset, qualitatively consistent with several laboratory-based estimates of n for specific phytoplankton species. A sensitivity analysis showed that the daily oscillation in n amplitude was somewhat insensitive to broad variations in method assumptions, ranging from 11.3 ± 4.3% to 16.9 ± 2.9%. Such estimates are crucial for improvement of algorithms that extract the particle size and production from bulk optical measurements, and could potentially help establish a link between n variations and changes in cellular composition of in situ particles.

Empirical relationship between nifH gene abundance and diazotroph cell concentration in the North Pacific Subtropical Gyre.

Cyanobacterial N2 -fixing microorganisms (diazotrophs) play a critical role in nitrogen and carbon cycling in the oceans; hence, accurate measurements of diazotroph abundance are imperative for understanding ocean biogeochemistry. Marine diazotroph abundances are often assessed using qPCR of the nifH gene, a sensitive, taxa-specific, and time/cost-efficient method. However, the validity of nifH abundance as a proxy for cell concentration has recently been questioned. Here, we compare nifH gene abundances to cell counts for four diazotroph taxa (Trichodesmium, Crocosphaera, Richelia, and Calothrix) on two cruises to the North Pacific Subtropical Gyre, one of the largest habitats for marine diazotrophs. nifH:cell relationships were strong and significant for Crocosphaera, Richelia, and Calothrix (nifH:cell 1.51-2.58; R2  = 0.89-0.96) but were not significant for Trichodesmium, despite previous studies reporting significant nifH:cell relationships for this organism. Limited available data suggest that empirical nifH:cell can vary among studies but that relationships are usually significantly linear and >1:1. Our study indicates that nifH gene abundance, while not a direct measure of cells, is a useful quantitative proxy for diazotroph abundance.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean.

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

The dynamic trophic architecture of open-ocean protist communities revealed through machine-guided metatranscriptomics.

Intricate networks of single-celled eukaryotes (protists) dominate carbon flow in the ocean. Their growth, demise, and interactions with other microorganisms drive the fluxes of biogeochemical elements through marine ecosystems. Mixotrophic protists are capable of both photosynthesis and ingestion of prey and are dominant components of open-ocean planktonic communities. Yet the role of mixotrophs in elemental cycling is obscured by their capacity to act as primary producers or heterotrophic consumers depending on factors that remain largely uncharacterized. Here, we develop and apply a machine learning model that predicts the in situ trophic mode of aquatic protists based on their patterns of gene expression. This approach leverages a public collection of protist transcriptomes as a training set to identify a subset of gene families whose transcriptional profiles predict trophic mode. We applied our model to nearly 100 metatranscriptomes obtained during two oceanographic cruises in the North Pacific and found community-level and population-specific evidence that abundant open-ocean mixotrophic populations shift their predominant mode of nutrient and carbon acquisition in response to natural gradients in nutrient supply and sea surface temperature. Metatranscriptomic data from ship-board incubation experiments revealed that abundant mixotrophic prymnesiophytes from the oligotrophic North Pacific subtropical gyre rapidly remodeled their transcriptome to enhance photosynthesis when supplied with limiting nutrients. Coupling this approach with experiments designed to reveal the mechanisms driving mixotroph physiology provides an avenue toward understanding the ecology of mixotrophy in the natural environment.

Overlooked and widespread pennate diatom-diazotroph symbioses in the sea.

Persistent nitrogen depletion in sunlit open ocean waters provides a favorable ecological niche for nitrogen-fixing (diazotrophic) cyanobacteria, some of which associate symbiotically with eukaryotic algae. All known marine examples of these symbioses have involved either centric diatom or haptophyte hosts. We report here the discovery and characterization of two distinct marine pennate diatom-diazotroph symbioses, which until now had only been observed in freshwater environments. Rhopalodiaceae diatoms Epithemia pelagica sp. nov. and Epithemia catenata sp. nov. were isolated repeatedly from the subtropical North Pacific Ocean, and analysis of sequence libraries reveals a global distribution. These symbioses likely escaped attention because the endosymbionts lack fluorescent photopigments, have nifH gene sequences similar to those of free-living unicellular cyanobacteria, and are lost in nitrogen-replete medium. Marine Rhopalodiaceae-diazotroph symbioses are a previously overlooked but widespread source of bioavailable nitrogen in marine habitats and provide new, easily cultured model organisms for the study of organelle evolution.

Complex marine microbial communities partition metabolism of scarce resources over the diel cycle.

Complex assemblages of microbes in the surface ocean are responsible for approximately half of global carbon fixation. The persistence of high taxonomic diversity despite competition for a small suite of relatively homogeneously distributed nutrients, that is, 'the paradox of the plankton', represents a long-standing challenge for ecological theory. Here we find evidence consistent with temporal niche partitioning of nitrogen assimilation processes over a diel cycle in the North Pacific Subtropical Gyre. We jointly analysed transcript abundances, lipids and metabolites and discovered that a small number of diel archetypes can explain pervasive periodic dynamics. Metabolic pathway analysis of identified diel signals revealed asynchronous timing in the transcription of nitrogen uptake and assimilation genes among different microbial groups-cyanobacteria, heterotrophic bacteria and eukaryotes. This temporal niche partitioning of nitrogen uptake emerged despite synchronous transcription of photosynthesis and central carbon metabolism genes and associated macromolecular abundances. Temporal niche partitioning may be a mechanism by which microorganisms in the open ocean mitigate competition for scarce resources, supporting community coexistence.

The rate and fate of N2 and C fixation by marine diatom-diazotroph symbioses.

N2 fixation constitutes an important new nitrogen source in the open sea. One group of filamentous N2 fixing cyanobacteria (Richelia intracellularis, hereafter Richelia) form symbiosis with a few genera of diatoms. High rates of N2 fixation and carbon (C) fixation have been measured in the presence of diatom-Richelia symbioses. However, it is unknown how partners coordinate C fixation and how the symbiont sustains high rates of N2 fixation. Here, both the N2 and C fixation in wild diatom-Richelia populations are reported. Inhibitor experiments designed to inhibit host photosynthesis, resulted in lower estimated growth and depressed C and N2 fixation, suggesting that despite the symbionts ability to fix their own C, they must still rely on their respective hosts for C. Single cell analysis indicated that up to 22% of assimilated C in the symbiont is derived from the host, whereas 78-91% of the host N is supplied from their symbionts. A size-dependent relationship is identified where larger cells have higher N2 and C fixation, and only N2 fixation was light dependent. Using the single cell measures, the N-rich phycosphere surrounding these symbioses was estimated and contributes directly and rapidly to the surface ocean rather than the mesopelagic, even at high estimated sinking velocities (<10 m d-1). Several eco-physiological parameters necessary for incorporating symbiotic N2 fixing populations into larger basin scale biogeochemical models (i.e., N and C cycles) are provided.

Marine Community Metabolomes Carry Fingerprints of Phytoplankton Community Composition.

Phytoplankton transform inorganic carbon into thousands of biomolecules that represent an important pool of fixed carbon, nitrogen, and sulfur in the surface ocean. Metabolite production differs between phytoplankton, and the flux of these molecules through the microbial food web depends on compound-specific bioavailability to members of a wider microbial community. Yet relatively little is known about the diversity or concentration of metabolites within marine plankton. Here, we compare 313 polar metabolites in 21 cultured phytoplankton species and in natural planktonic communities across environmental gradients to show that bulk community metabolomes reflect the chemical composition of the phytoplankton community. We also show that groups of compounds have similar patterns across space and taxonomy, suggesting that the concentrations of these compounds in the environment are controlled by similar sources and sinks. We quantify several compounds in the surface ocean that represent substantial understudied pools of labile carbon. For example, the N-containing metabolite homarine was up to 3% of particulate carbon and is produced in high concentrations by cultured Synechococcus, and S-containing gonyol accumulated up to 2.5 nM in surface particles and likely originates from dinoflagellates or haptophytes. Our results show that phytoplankton composition directly shapes the carbon composition of the surface ocean. Our findings suggest that in order to access these pools of bioavailable carbon, the wider microbial community must be adapted to phytoplankton community composition.IMPORTANCE Microscopic phytoplankton transform 100 million tons of inorganic carbon into thousands of different organic compounds each day. The structure of each chemical is critical to its biological and ecosystem function, yet the diversity of biomolecules produced by marine microbial communities remained mainly unexplored, especially small polar molecules which are often considered the currency of the microbial loop. Here, we explore the abundance and diversity of small biomolecules in planktonic communities across ecological gradients in the North Pacific and within 21 cultured phytoplankton species. Our work demonstrates that phytoplankton diversity is an important determinant of the chemical composition of the highly bioavailable pool of organic carbon in the ocean, and we highlight understudied yet abundant compounds in both the environment and cultured organisms. These findings add to understanding of how the chemical makeup of phytoplankton shapes marine microbial communities where the ability to sense and use biomolecules depends on the chemical structure.

Particulate Metabolites and Transcripts Reflect Diel Oscillations of Microbial Activity in the Surface Ocean.

Light fuels photosynthesis and organic matter production by primary producers in the sunlit ocean. The quantity and quality of the organic matter produced influence community function, yet in situ measurements of metabolites, the products of cellular metabolism, over the diel cycle are lacking. We evaluated community-level biochemical consequences of oscillations of light in the North Pacific Subtropical Gyre by quantifying 79 metabolites in particulate organic matter from 15 m every 4 h over 8 days. Total particulate metabolite concentration peaked at dusk and represented up to 2% of total particulate organic carbon (POC). The concentrations of 55/79 (70%) individual metabolites exhibited significant 24-h periodicity, with daily fold changes from 1.6 to 12.8, often greater than those of POC and flow cytometry-resolvable biomass, which ranged from 1.2 to 2.8. Paired metatranscriptome analysis revealed the taxa involved in production and consumption of a subset of metabolites. Primary metabolites involved in anabolism and redox maintenance had significant 24-h periodicity and diverse organisms exhibited diel periodicity in transcript abundance associated with these metabolites. Compounds with osmotic properties displayed the largest oscillations in concentration, implying rapid turnover and supporting prior evidence of functions beyond cell turgor maintenance. The large daily oscillation of trehalose paired with metatranscriptome and culture data showed that trehalose is produced by the nitrogen-fixing cyanobacterium Crocosphaera, likely to store energy for nighttime metabolism. Together, paired measurements of particulate metabolites and transcripts resolve strategies that microbes use to manage daily energy and redox oscillations and highlight dynamic metabolites with cryptic roles in marine microbial ecosystems.IMPORTANCE Fueled by light, phytoplankton produce the organic matter that supports ocean ecosystems and carbon sequestration. Ocean change impacts microbial metabolism with repercussions for biogeochemical cycling. As the small molecule products of cellular metabolism, metabolites often change rapidly in response to environmental conditions and form the basis of energy and nutrient management and storage within cells. By pairing measurements of metabolites and gene expression in the stratified surface ocean, we reveal strategies of microbial energy management over the day-night cycle and hypothesize that oscillating metabolites are important substrates for dark respiration by phytoplankton. These high-resolution diel measurements of in situ metabolite concentrations form the basis for future work into the specific roles these compounds play in marine microbial communities.

A system of coordinated autonomous robots for Lagrangian studies of microbes in the oceanic deep chlorophyll maximum.

The deep chlorophyll maximum (DCM) layer is an ecologically important feature of the open ocean. The DCM cannot be observed using aerial or satellite remote sensing; thus, in situ observations are essential. Further, understanding the responses of microbes to the environmental processes driving their metabolism and interactions requires observing in a reference frame that moves with a plankton population drifting in ocean currents, i.e., Lagrangian. Here, we report the development and application of a system of coordinated robots for studying planktonic biological communities drifting within the ocean. The presented Lagrangian system uses three coordinated autonomous robotic platforms. The focal platform consists of an autonomous underwater vehicle (AUV) fitted with a robotic water sampler. This platform localizes and drifts within a DCM community, periodically acquiring samples while continuously monitoring the local environment. The second platform is an AUV equipped with environmental sensing and acoustic tracking capabilities. This platform characterizes environmental conditions by tracking the focal platform and vertically profiling in its vicinity. The third platform is an autonomous surface vehicle equipped with satellite communications and subsea acoustic tracking capabilities. While also acoustically tracking the focal platform, this vehicle serves as a communication relay that connects the subsea robot to human operators, thereby providing situational awareness and enabling intervention if needed. Deployed in the North Pacific Ocean within the core of a cyclonic eddy, this coordinated system autonomously captured fundamental characteristics of the in situ DCM microbial community in a manner not possible previously.

Evaluation of argon-induced hydrogen production as a method to measure nitrogen fixation by cyanobacteria.

The production of dihydrogen (H2 ) is an enigmatic yet obligate component of biological dinitrogen (N2 ) fixation. This study investigates the effect on H2 production by N2 fixing cyanobacteria when they are exposed to either air or a gas mixture consisting of argon, oxygen, and carbon dioxide (Ar:O2 :CO2 ). In the absence of N2 , nitrogenase diverts the flow of electrons to the production of H2 , which becomes a measure of Total Nitrogenase Activity (TNA). This method of argon-induced hydrogen production (AIHP) is much less commonly used to infer rates of N2 fixation than the acetylene reduction (AR) assay. We provide here a full evaluation of the AIHP method and demonstrate its ability to achieve high-resolution measurements of TNA in a gas exchange flow-through system. Complete diel profiles of H2 production were obtained for N2 fixing cyanobacteria despite the absence of N2 that broadly reproduced the temporal patterns observed by the AR assay. Comparison of H2 production under air versus Ar:O2 :CO2 revealed the efficiency of electron usage during N2 fixation and place these findings in the broader context of cell metabolism. Ultimately, AIHP is demonstrated to be a viable alternative to the AR assay with several additional merits that provide an insight into cell physiology and promise for successful field application.

Anthropogenic Asian aerosols provide Fe to the North Pacific Ocean.

Fossil-fuel emissions may impact phytoplankton primary productivity and carbon cycling by supplying bioavailable Fe to remote areas of the ocean via atmospheric aerosols. However, this pathway has not been confirmed by field observations of anthropogenic Fe in seawater. Here we present high-resolution trace-metal concentrations across the North Pacific Ocean (158°W from 25°to 42°N). A dissolved Fe maximum was observed around 35°N, coincident with high dissolved Pb and Pb isotope ratios matching Asian industrial sources and confirming recent aerosol deposition. Iron-stable isotopes reveal in situ evidence of anthropogenic Fe in seawater, with low δ56Fe (-0.23‰ > δ56Fe > -0.65‰) observed in the region that is most influenced by aerosol deposition. An isotope mass balance suggests that anthropogenic Fe contributes 21-59% of dissolved Fe measured between 35° and 40°N. Thus, anthropogenic aerosol Fe is likely to be an important Fe source to the North Pacific Ocean.

Diel variability of bulk optical properties associated with the growth and division of small phytoplankton in the North Pacific Subtropical Gyre.

Cross-platform observing systems are requisite to capturing the temporal and spatial dynamics of particles in the ocean. We present simultaneous observations of bulk optical properties, including the particulate beam attenuation (cp) and backscattering (bbp) coefficients, and particle size distributions collected in the North Pacific Subtropical Gyre. Clear and coherent diel cycles are observed in all bulk and size-fractionated optical proxies for particle biomass. We show evidence linking diurnal increases in cp and bbp to daytime particle growth and division of cells, with particles <7µm driving the daily cycle of particle production and loss within the mixed layer. Flow cytometry data reveal the nitrogen-fixing cyanobacterium Crocosphaera (∼4-7µm) to be an important driver of cp at the time of sampling, whereas Prochlorococcus dynamics (∼0.5µm) were essential to reproducing temporal variability in bbp. This study is a step towards improved characterization of the particle size range represented by in situ bulk optical properties and a better understanding of the mechanisms that drive variability in particle production in the oligotrophic open ocean.

Biological Control of Chromium Redox and Stable Isotope Composition in the Surface Ocean.

While chromium stable isotopes (δ53Cr) have received significant attention for their utility as a tracer of oxygen availability in the distant geological past, a mechanistic understanding of modern oceanic controls on Cr and δ53Cr is still lacking. Here we present total dissolved δ53Cr, concentrations of Cr (III) and total dissolved Cr, and net community productivity (NCP) from the North Pacific. Chromium concentrations show surface depletions in waters with elevated NCP, but not in lower productivity waters. Observed Cr deficits correspond well with calculated Cr export derived from NCP and Cr:C ratios of natural phytoplankton and marine particulates. Chromium (III) concentrations are stable over the diel cycle yet correlate with NCP, with maxima found in highly productive surface waters but not in lower productivity waters, indicating biological control on Cr (III). The relationship between Cr (III) and δ53Cr suggests that δ53Cr distributions may be controlled by the removal of isotopically light Cr (III) at an isotopic enrichment factor (∆53Cr) of -1.08‰ ± 0.25 relative to total dissolved δ53Cr, in agreement with the global δ53Cr-Cr fractionation factor (-0.82‰ ± 0.05). No perturbation to δ53Cr, Cr, or Cr (III) is observed in oxygen-depleted waters (~10 µmol/kg), suggesting no strong control by O2 availability, in agreement with other recent studies. Therefore, we propose that biological productivity is the primary control on Cr and δ53Cr in the modern ocean. Consequently, δ53Cr records in marine sediments may not faithfully record oxygen availability in the Late Quaternary. Instead, our data demonstrate that δ53Cr records may be a useful tracer for biological productivity.

Climate-driven oscillation of phosphorus and iron limitation in the North Pacific Subtropical Gyre.

The supply of nutrients is a fundamental regulator of ocean productivity and carbon sequestration. Nutrient sources, sinks, residence times, and elemental ratios vary over broad scales, including those resulting from climate-driven changes in upper water column stratification, advection, and the deposition of atmospheric dust. These changes can alter the proximate elemental control of ecosystem productivity with cascading ecological effects and impacts on carbon sequestration. Here, we report multidecadal observations revealing that the ecosystem in the eastern region of the North Pacific Subtropical Gyre (NPSG) oscillates on subdecadal scales between inorganic phosphorus (P i ) sufficiency and limitation, when P i concentration in surface waters decreases below 50-60 nmol⋅kg-1 In situ observations and model simulations suggest that sea-level pressure changes over the northwest Pacific may induce basin-scale variations in the atmospheric transport and deposition of Asian dust-associated iron (Fe), causing the eastern portion of the NPSG ecosystem to shift between states of Fe and P i limitation. Our results highlight the critical need to include both atmospheric and ocean circulation variability when modeling the response of open ocean pelagic ecosystems under future climate change scenarios.

Daily changes in phytoplankton lipidomes reveal mechanisms of energy storage in the open ocean.

Sunlight is the dominant control on phytoplankton biosynthetic activity, and darkness deprives them of their primary external energy source. Changes in the biochemical composition of phytoplankton communities over diel light cycles and attendant consequences for carbon and energy flux in environments remain poorly elucidated. Here we use lipidomic data from the North Pacific subtropical gyre to show that biosynthesis of energy-rich triacylglycerols (TAGs) by eukaryotic nanophytoplankton during the day and their subsequent consumption at night drives a large and previously uncharacterized daily carbon cycle. Diel oscillations in TAG concentration comprise 23 ± 11% of primary production by eukaryotic nanophytoplankton representing a global flux of about 2.4 Pg C yr-1. Metatranscriptomic analyses of genes required for TAG biosynthesis indicate that haptophytes and dinoflagellates are active members in TAG production. Estimates suggest that these organisms could contain as much as 40% more calories at sunset than at sunrise due to TAG production.