Heterotrophic bacteria trigger transcriptome remodelling in the photosynthetic picoeukaryote Micromonas commoda.

Marine biogeochemical cycles are built on interactions between surface ocean microbes, particularly those connecting phytoplankton primary producers to heterotrophic bacteria. Details of these associations are not well understood, especially in the case of direct influences of bacteria on phytoplankton physiology. Here we catalogue how the presence of three marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14 and Polaribacter dokdonensis MED152) individually and uniquely impact gene expression of the picoeukaryotic alga Micromonas commoda RCC 299. We find a dramatic transcriptomic remodelling by M. commoda after 8 h in co-culture, followed by an increase in cell numbers by 56 h compared with the axenic cultures. Some aspects of the algal transcriptomic response are conserved across all three bacterial co-cultures, including an unexpected reduction in relative expression of photosynthesis and carbon fixation pathways. Expression differences restricted to a single bacterium are also observed, with the Flavobacteriia P. dokdonensis uniquely eliciting changes in relative expression of algal genes involved in biotin biosynthesis and the acquisition and assimilation of nitrogen. This study reveals that M. commoda has rapid and extensive responses to heterotrophic bacteria in ways that are generalizable, as well as in a taxon specific manner, with implications for the diversity of phytoplankton-bacteria interactions ongoing in the surface ocean.

Functional annotation and importance of marine bacterial transporters of plankton exometabolites.

Metabolite exchange within marine microbial communities transfers carbon and other major elements through global cycles and forms the basis of microbial interactions. Yet lack of gene annotations and concern about the quality of existing ones remain major impediments to revealing currencies of carbon flux. We employed an arrayed mutant library of the marine bacterium Ruegeria pomeroyi DSS-3 to experimentally annotate substrates of organic compound transporter systems, using mutant growth and compound drawdown analyses to link transporters to their cognate substrates. Mutant experiments verified substrates for thirteen R. pomeroyi transporters. Four were previously hypothesized based on gene expression data (taurine, glucose/xylose, isethionate, and cadaverine/putrescine/spermidine); five were previously hypothesized based on homology to experimentally annotated transporters in other bacteria (citrate, glycerol, N-acetylglucosamine, fumarate/malate/succinate, and dimethylsulfoniopropionate); and four had no previous annotations (thymidine, carnitine, cysteate, and 3-hydroxybutyrate). These bring the total number of experimentally-verified organic carbon influx transporters to 18 of 126 in the R. pomeroyi genome. In a longitudinal study of a coastal phytoplankton bloom, expression patterns of the experimentally annotated transporters linked them to different stages of the bloom, and also led to the hypothesis that citrate and 3-hydroxybutyrate were among the most highly available bacterial substrates. Improved functional annotation of the gatekeepers of organic carbon uptake is critical for deciphering carbon flux and fate in microbial ecosystems.

Ecological divergence of syntopic marine bacterial species is shaped by gene content and expression.

Identifying mechanisms by which bacterial species evolve and maintain genomic diversity is particularly challenging for the uncultured lineages that dominate the surface ocean. A longitudinal analysis of bacterial genes, genomes, and transcripts during a coastal phytoplankton bloom revealed two co-occurring, highly related Rhodobacteraceae species from the deeply branching and uncultured NAC11-7 lineage. These have identical 16S rRNA gene amplicon sequences, yet their genome contents assembled from metagenomes and single cells indicate species-level divergence. Moreover, shifts in relative dominance of the species during dynamic bloom conditions over 7 weeks confirmed the syntopic species' divergent responses to the same microenvironment at the same time. Genes unique to each species and genes shared but divergent in per-cell inventories of mRNAs accounted for 5% of the species' pangenome content. These analyses uncover physiological and ecological features that differentiate the species, including capacities for organic carbon utilization, attributes of the cell surface, metal requirements, and vitamin biosynthesis. Such insights into the coexistence of highly related and ecologically similar bacterial species in their shared natural habitat are rare.

A mutant fitness assay identifies bacterial interactions in a model ocean hot spot.

Bacteria that assemble in phycospheres surrounding living phytoplankton cells metabolize a substantial proportion of ocean primary productivity. Yet the type and extent of interactions occurring among species that colonize these micron-scale "hot spot" environments are challenging to study. We identified genes that mediate bacterial interactions in phycosphere communities by culturing a transposon mutant library of copiotrophic bacterium Ruegeria pomeroyi DSS-3 with the diatom Thalassiosira pseudonana CCMP1335 as the sole source of organic matter in the presence or absence of other heterotrophic bacterial species. The function of genes having significant effects on R. pomeroyi fitness indicated explicit cell-cell interactions initiated in the multibacterial phycospheres. We found that R. pomeroyi simultaneously competed for shared substrates while increasing reliance on substrates that did not support the other species' growth. Fitness outcomes also indicated that the bacterium competed for nitrogen in the forms of ammonium and amino acids; obtained purines, pyrimidines, and cofactors via crossfeeding; both initiated and defended antagonistic interactions; and sensed an environment with altered oxygen and superoxide levels. The large genomes characteristic of copiotrophic marine bacteria are hypothesized to enable responses to dynamic ecological challenges occurring at the scale of microns. Here, we discover >200 nonessential genes implicated in the management of fitness costs and benefits of membership in a globally significant bacterial community.

13C NMR metabolomics: J-resolved STOCSY meets INADEQUATE.

Robust annotation of metabolites is a challenging task in metabolomics. Among available applications, 13C NMR experiment INADEQUATE determines direct 13C-13C connectivity unambiguously, offering indispensable information on molecular structure. Despite its great utility, it is not always practical to collect INADEQUATE data on every sample in a large metabolomics study because of its relatively long experiment time. Here, we propose an alternative approach that maintains the quality of information but saves experiment time. In this approach, individual samples in a study are first screened by 13C homonuclear J-resolved experiment (JRES). Next, JRES data are processed by statistical total correlation spectroscopy (STOCSY) to extract peaks that behave similarly among samples. Finally, INADEQUATE is collected on one internal pooled sample to select STOCSY peaks that originate from the same compound. We tested this concept using the 13C-labeled endometabolome of a model marine diatom strain incubated under various settings, intending to cover a range of metabolites produced under different external conditions. This scheme was able to extract known diatom metabolites proline, 2,3-dihydroxypropane-1-sulfonate (DHPS), ß-1,3-glucan, choline, and glutamate. This pipeline also detected unknown compounds with structural information, which is valuable in metabolomics where a priori knowledge of metabolites is not always available. The ability of this scheme was seen even in sugar regions, which are usually challenging in 1H NMR due to severe peak overlap. JRES and INADEQUATE were highly complementary; INADEQUATE provided directly-bonded 13C networks, whereas JRES linked INADEQUATE networks within the same compound but broken by nitrogen or sulfur atoms, highlighting the advantage of this integrated approach.

Bacterial transcriptional response to labile exometabolites from photosynthetic picoeukaryote Micromonas commoda.

Dissolved primary production released into seawater by marine phytoplankton is a major source of carbon fueling heterotrophic bacterial production in the ocean. The composition of the organic compounds released by healthy phytoplankton is poorly known and difficult to assess with existing chemical methods. Here, expression of transporter and catabolic genes by three model marine bacteria (Ruegeria pomeroyi DSS-3, Stenotrophomonas sp. SKA14, and Polaribacter dokdonensis MED152) was used as a biological sensor of metabolites released from the picoeukaryote Micromonas commoda RCC299. Bacterial expression responses indicated that the three species together recognized 38 picoeukaryote metabolites. This was consistent with the Micromonas expression of genes for starch metabolism and synthesis of peptidoglycan-like intermediates. A comparison of the hypothesized Micromonas exometabolite pool with that of the diatom Thalassiosira pseudonana CCMP1335, analyzed previously with the same biological sensor method, indicated that both phytoplankton released organic acids, nucleosides, and amino acids, but differed in polysaccharide and organic nitrogen release. Future ocean conditions are expected to favor picoeukaryotic phytoplankton over larger-celled microphytoplankton. Results from this study suggest that such a shift could alter the substrate pool available to heterotrophic bacterioplankton.

Dissolved organic metabolite extraction from high-salt media.

We describe considerations and strategies for developing a nuclear magnetic resonance (NMR) sample preparation method to extract low molecular weight metabolites from high-salt spent media in a model coculture system of phytoplankton and marine bacteria. Phytoplankton perform half the carbon fixation and oxygen generation on Earth. A substantial fraction of fixed carbon becomes part of a metabolite pool of small molecules known as dissolved organic matter (DOM), which are taken up by marine bacteria proximate to phytoplankton. There is an urgent need to elucidate these metabolic exchanges due to widespread anthropogenic transformations on the chemical, phenotypic, and species composition of seawater. These changes are increasing water temperature and the amount of CO2 absorbed by the ocean at energetic costs to marine microorganisms. Little is known about the metabolite-mediated, structured interactions occurring between phytoplankton and associated marine bacteria, in part because of challenges in studying high-salt solutions on various analytical platforms. NMR analysis is problematic due to the high-salt content of both natural seawater and culture media for marine microbes. High-salt concentration degrades the performance of the radio frequency coil, reduces the efficiency of some pulse sequences, limits signal-to-noise, and prolongs experimental time. The method described herein can reproducibly extract low molecular weight DOM from small-volume, high-salt cultures. It is a promising tool for elucidating metabolic flux between marine microorganisms and facilitates genetic screens of mutant microorganisms.

Divergent gene expression responses in two Baltic Sea heterotrophic model bacteria to dinoflagellate dissolved organic matter.

Phytoplankton release massive amounts of dissolved organic matter (DOM) into the water column during recurring blooms in coastal waters and inland seas. The released DOM encompasses a complex mixture of both known and unknown compounds, and is a rich nutrient source for heterotrophic bacteria. The metabolic activity of bacteria during and after phytoplankton blooms can hence be expected to reflect the characteristics of the released DOM. We therefore investigated if bacterioplankton could be used as "living sensors" of phytoplankton DOM quantity and/or quality, by applying gene expression analyses to identify bacterial metabolisms induced by DOM. We used transcriptional analysis of two Baltic Sea bacterial isolates (Polaribacter sp. BAL334 [Flavobacteriia] and Brevundimonas sp. BAL450 [Alphaproteobacteria]) growing with DOM from axenic cultures of the dinoflagellate Prorocentrum minimum. We observed pronounced differences between the two bacteria both in growth and the expressed metabolic pathways in cultures exposed to dinoflagellate DOM compared with controls. Differences in metabolic responses between the two isolates were caused both by differences in gene repertoire between them (e.g. in the SEED categories for membrane transport, motility and photoheterotrophy) and the regulation of expression (e.g. fatty acid metabolism), emphasizing the importance of separating the responses of different taxa in analyses of community sequence data. Similarities between the bacteria included substantially increased expression of genes for Ton and Tol transport systems in both isolates, which are commonly associated with uptake of complex organic molecules. Polaribacter sp. BAL334 showed stronger metabolic responses to DOM harvested from exponential than stationary phase dinoflagellates (128 compared to 26 differentially expressed genes), whereas Brevundimonas sp. BAL450 responded more to the DOM from stationary than exponential phase dinoflagellates (33 compared to 6 differentially expressed genes). These findings suggest that shifts in bacterial metabolisms during different phases of phytoplankton blooms can be detected in individual bacterial species and can provide insights into their involvement in DOM transformations.

Oxidative Stress Regulates a Pivotal Metabolic Switch in Dimethylsulfoniopropionate Degradation by the Marine Bacterium Ruegeria pomeroyi.

Dimethylsulfoniopropionate (DMSP) is an abundant organic compound in marine surface water and source of dimethyl sulfide (DMS), the largest natural sulfur source to the upper atmosphere. Marine bacteria either mineralize DMSP through the demethylation pathway or transform it to DMS through the cleavage pathway. Factors that regulate which pathway is utilized are not fully understood. In chemostat experiments, the marine Roseobacter Ruegeria pomeroyi DSS-3 was exposed to oxidative stress either during growth with H2O2 or by mutation of the gene encoding catalase. Oxidative stress reduced expression of the genes in the demethylation pathway and increased expression of those encoding the cleavage pathway. These results are contrary to the sulfur demand hypothesis, which theorizes that DMSP metabolism is driven by sulfur requirements of bacterial cells. Instead, we find strong evidence consistent with oxidative stress control over the switch in DMSP metabolism from demethylation to DMS production in an ecologically relevant marine bacterium. IMPORTANCE Dimethylsulfoniopropionate (DMSP) is the most abundant low-molecular-weight organic compound in marine surface water and source of dimethyl sulfide (DMS), a climatically active gas that connects the marine and terrestrial sulfur cycles. Marine bacteria are the major DMSP consumers, either generating DMS or consuming DMSP as a source of reduced carbon and sulfur. However, the factors regulating the DMSP catabolism in bacteria are not well understood. Marine bacteria are also exposed to oxidative stress. RNA sequencing (RNA-seq) experiments showed that oxidative stress induced in the laboratory reduced expression of the genes encoding the consumption of DMSP via the demethylation pathway and increased the expression of genes encoding DMS production via the cleavage pathway in the marine bacterium Ruegeria pomeroyi. These results support a model where DMS production in the ocean is regulated in part by oxidative stress.

Microbial metabolites in the marine carbon cycle.

One-quarter of photosynthesis-derived carbon on Earth rapidly cycles through a set of short-lived seawater metabolites that are generated from the activities of marine phytoplankton, bacteria, grazers and viruses. Here we discuss the sources of microbial metabolites in the surface ocean, their roles in ecology and biogeochemistry, and approaches that can be used to analyse them from chemistry, biology, modelling and data science. Although microbial-derived metabolites account for only a minor fraction of the total reservoir of marine dissolved organic carbon, their flux and fate underpins the central role of the ocean in sustaining life on Earth.

Growth-stage-related shifts in diatom endometabolome composition set the stage for bacterial heterotrophy.

Phytoplankton-derived metabolites fuel a large fraction of heterotrophic bacterial production in the global ocean, yet methodological challenges have limited our understanding of the organic molecules transferred between these microbial groups. In an experimental bloom study consisting of three heterotrophic marine bacteria growing together with the diatom Thalassiosira pseudonana, we concurrently measured diatom endometabolites (i.e., potential exometabolite supply) by nuclear magnetic resonance (NMR) spectroscopy and bacterial gene expression (i.e., potential exometabolite uptake) by metatranscriptomic sequencing. Twenty-two diatom endometabolites were annotated, with nine increasing in internal concentration in the late stage of the bloom, eight decreasing, and five showing no variation through the bloom progression. Some metabolite changes could be linked to shifts in diatom gene expression, as well as to shifts in bacterial community composition and their expression of substrate uptake and catabolism genes. Yet an overall low match indicated that endometabolome concentration was not a good predictor of exometabolite availability, and that complex physiological and ecological interactions underlie metabolite exchange. Six diatom endometabolites accumulated to higher concentrations in the bacterial co-cultures compared to axenic cultures, suggesting a bacterial influence on rates of synthesis or release of glutamate, arginine, leucine, 2,3-dihydroxypropane-1-sulfonate, glucose, and glycerol-3-phosphate. Better understanding of phytoplankton metabolite production, release, and transfer to assembled bacterial communities is key to untangling this nearly invisible yet pivotal step in ocean carbon cycling.

Diel investments in metabolite production and consumption in a model microbial system.

Organic carbon transfer between surface ocean photosynthetic and heterotrophic microbes is a central but poorly understood process in the global carbon cycle. In a model community in which diatom extracellular release of organic molecules sustained growth of a co-cultured bacterium, we determined quantitative changes in the diatom endometabolome and the bacterial uptake transcriptome over two diel cycles. Of the nuclear magnetic resonance (NMR) peaks in the diatom endometabolites, 38% had diel patterns with noon or mid-afternoon maxima; the remaining either increased (36%) or decreased (26%) through time. Of the genes in the bacterial uptake transcriptome, 94% had a diel pattern with a noon maximum; the remaining decreased over time (6%). Eight diatom endometabolites identified with high confidence were matched to the bacterial genes mediating their utilization. Modeling of these coupled inventories with only diffusion-based phytoplankton extracellular release could not reproduce all the patterns. Addition of active release mechanisms for physiological balance and bacterial recognition significantly improved model performance. Estimates of phytoplankton extracellular release range from only a few percent to nearly half of annual net primary production. Improved understanding of the factors that influence metabolite release and consumption by surface ocean microbes will better constrain this globally significant carbon flux.

Bacterial responses to background organic pollutants in the northeast subarctic Pacific Ocean.

Thousands of man-made synthetic chemicals are released to oceans and compose the anthropogenic dissolved organic carbon (ADOC). Little is known about the effects of this chronic pollution on marine microbiome activities. In this study, we measured the pollution level at three sites in the Northeast Subarctic Pacific Ocean (NESAP) and investigated how mixtures of three model families of ADOC at different environmentally relevant concentrations affected naturally occurring marine bacterioplankton communities' structure and metabolic functioning. The offshore northernmost site (North) had the lowest concentrations of hydrocarbons, as well as organophosphate ester plasticizers, contrasting with the two other continental shelf sites, the southern coastal site (South) being the most contaminated. At North, ADOC stimulated bacterial growth and promoted an increase in the contribution of some Gammaproteobacteria groups (e.g. Alteromonadales) to the 16 rRNA pool. These groups are described as fast responders after oil spills. In contrast, minor changes in South microbiome activities were observed. Gene expression profiles at Central showed the coexistence of ADOC degradation and stress-response strategies to cope with ADOC toxicities. These results show that marine microbial communities at three distinct domains in NESAP are influenced by background concentrations of ADOC, expanding previous assessments for polar and temperate waters.

Quantification of Amine- and Alcohol-Containing Metabolites in Saline Samples Using Pre-extraction Benzoyl Chloride Derivatization and Ultrahigh Performance Liquid Chromatography Tandem Mass Spectrometry (UHPLC MS/MS).

Dissolved metabolites serve as nutrition, energy, and chemical signals for microbial systems. However, the full scope and magnitude of these processes in marine systems are unknown, largely due to insufficient methods, including poor extraction of small, polar compounds using common solid-phase extraction resins. Here, we utilized pre-extraction derivatization and ultrahigh performance liquid chromatography electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) to detect and quantify targeted dissolved metabolites in seawater and saline culture media. Metabolites were derivatized with benzoyl chloride by their primary and secondary amine and alcohol functionalities and quantified using stable isotope-labeled internal standards (SIL-ISs) produced from 13C6-labeled benzoyl chloride. We optimized derivatization, extraction, and sample preparation for field and culture samples and evaluated matrix-derived biases. We have optimized this quantitative method for 73 common metabolites, of which 50 cannot be quantified without derivatization due to low extraction efficiencies. Of the 73 metabolites, 66 were identified in either culture media or seawater and 45 of those were quantified. This derivatization method is sensitive (detection limits = pM to nM), rapid (∼5 min per sample), and high throughput.

Niche dimensions of a marine bacterium are identified using invasion studies in coastal seawater.

Niche theory is a foundational ecological concept that explains the distribution of species in natural environments. Identifying the dimensions of any organism's niche is challenging because numerous environmental factors can affect organism viability. We used serial invasion experiments to introduce Ruegeria pomeroyi DSS-3, a heterotrophic marine bacterium, into a coastal phytoplankton bloom on 14 dates. RNA-sequencing analysis of R. pomeroyi was conducted after 90 min to assess its niche dimensions in this dynamic ecosystem. We identified ~100 external conditions eliciting transcriptional responses, which included substrates, nutrients, metals and biotic interactions such as antagonism, resistance and cofactor synthesis. The peak bloom was characterized by favourable states for most of the substrate dimensions, but low inferred growth rates of R. pomeroyi at this stage indicated that its niche was narrowed by factors other than substrate availability, most probably negative biotic interactions with the bloom dinoflagellate. Our findings indicate chemical and biological features of the ocean environment that can constrain where heterotrophic bacteria survive.

Resource partitioning of phytoplankton metabolites that support bacterial heterotrophy.

The communities of bacteria that assemble around marine microphytoplankton are predictably dominated by Rhodobacterales, Flavobacteriales, and families within the Gammaproteobacteria. Yet whether this consistent ecological pattern reflects the result of resource-based niche partitioning or resource competition requires better knowledge of the metabolites linking microbial autotrophs and heterotrophs in the surface ocean. We characterized molecules targeted for uptake by three heterotrophic bacteria individually co-cultured with a marine diatom using two strategies that vetted the exometabolite pool for biological relevance by means of bacterial activity assays: expression of diagnostic genes and net drawdown of exometabolites, the latter detected with mass spectrometry and nuclear magnetic resonance using novel sample preparation approaches. Of the more than 36 organic molecules with evidence of bacterial uptake, 53% contained nitrogen (including nucleosides and amino acids), 11% were organic sulfur compounds (including dihydroxypropanesulfonate and dimethysulfoniopropionate), and 28% were components of polysaccharides (including chrysolaminarin, chitin, and alginate). Overlap in phytoplankton-derived metabolite use by bacteria in the absence of competition was low, and only guanosine, proline, and N-acetyl-D-glucosamine were predicted to be used by all three. Exometabolite uptake pattern points to a key role for ecological resource partitioning in the assembly marine bacterial communities transforming recent photosynthate.

Transcriptional activity differentiates families of Marine Group II Euryarchaeota in the coastal ocean.

Marine Group II Euryarchaeota (Candidatus Poseidoniales), abundant but yet-uncultivated members of marine microbial communities, are thought to be (photo)heterotrophs that metabolize dissolved organic matter (DOM), such as lipids and peptides. However, little is known about their transcriptional activity. We mapped reads from a metatranscriptomic time series collected at Sapelo Island (GA, USA) to metagenome-assembled genomes to determine the diversity of transcriptionally active Ca. Poseidoniales. Summer metatranscriptomes had the highest abundance of Ca. Poseidoniales transcripts, mostly from the O1 and O3 genera within Ca. Thalassarchaeaceae (MGIIb). In contrast, transcripts from fall and winter samples were predominantly from Ca. Poseidoniaceae (MGIIa). Genes encoding proteorhodopsin, membrane-bound pyrophosphatase, peptidase/proteases, and part of the ß-oxidation pathway were highly transcribed across abundant genera. Highly transcribed genes specific to Ca. Thalassarchaeaceae included xanthine/uracil permease and receptors for amino acid transporters. Enrichment of Ca. Thalassarchaeaceae transcript reads related to protein/peptide, nucleic acid, and amino acid transport and metabolism, as well as transcript depletion during dark incubations, provided further evidence of heterotrophic metabolism. Quantitative PCR analysis of South Atlantic Bight samples indicated consistently abundant Ca. Poseidoniales in nearshore and inshore waters. Together, our data suggest that Ca. Thalassarchaeaceae are important photoheterotrophs potentially linking DOM and nitrogen cycling in coastal waters.

Genome Sequences and Metagenome-Assembled Genome Sequences of Microbial Communities Enriched on Phytoplankton Exometabolites.

We report 11 bacterial draft genome sequences and 38 metagenome-assembled genomes (MAGs) from marine phytoplankton exometabolite enrichments. The genomes and MAGs represent marine bacteria adapted to the metabolite environment of phycospheres, organic matter-rich regions surrounding phytoplankton cells, and are useful for exploring functional and taxonomic attributes of phytoplankton-associated bacterial communities.

Ecological drivers of bacterial community assembly in synthetic phycospheres.

In the nutrient-rich region surrounding marine phytoplankton cells, heterotrophic bacterioplankton transform a major fraction of recently fixed carbon through the uptake and catabolism of phytoplankton metabolites. We sought to understand the rules by which marine bacterial communities assemble in these nutrient-enhanced phycospheres, specifically addressing the role of host resources in driving community coalescence. Synthetic systems with varying combinations of known exometabolites of marine phytoplankton were inoculated with seawater bacterial assemblages, and communities were transferred daily to mimic the average duration of natural phycospheres. We found that bacterial community assembly was predictable from linear combinations of the taxa maintained on each individual metabolite in the mixture, weighted for the growth each supported. Deviations from this simple additive resource model were observed but also attributed to resource-based factors via enhanced bacterial growth when host metabolites were available concurrently. The ability of photosynthetic hosts to shape bacterial associates through excreted metabolites represents a mechanism by which microbiomes with beneficial effects on host growth could be recruited. In the surface ocean, resource-based assembly of host-associated communities may underpin the evolution and maintenance of microbial interactions and determine the fate of a substantial portion of Earth's primary production.

Correction: Single-cell genomics shedding light on marine Thaumarchaeota diversification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.