Unicellular eukaryotic phytoplankton, such as diatoms, rely on microbial communities for survival despite lacking specialized compartments to house microbiomes (e.g., animal gut). Microbial communities have been widely shown to benefit from diatom excretions that accumulate within the microenvironment surrounding phytoplankton cells, known as the phycosphere. However, mechanisms that enable diatoms and other unicellular eukaryotes to nurture specific microbiomes by fostering beneficial bacteria and repelling harmful ones are mostly unknown. We hypothesized that diatom exudates may tune microbial communities and employed an integrated multiomics approach using the ubiquitous diatom Asterionellopsis glacialis to reveal how it modulates its naturally associated bacteria. We show that A. glacialis reprograms its transcriptional and metabolic profiles in response to bacteria to secrete a suite of central metabolites and two unusual secondary metabolites, rosmarinic acid and azelaic acid. While central metabolites are utilized by potential bacterial symbionts and opportunists alike, rosmarinic acid promotes attachment of beneficial bacteria to the diatom and simultaneously suppresses the attachment of opportunists. Similarly, azelaic acid enhances growth of beneficial bacteria while simultaneously inhibiting growth of opportunistic ones. We further show that the bacterial response to azelaic acid is numerically rare but globally distributed in the world’s oceans and taxonomically restricted to a handful of bacterial genera. Our results demonstrate the innate ability of an important unicellular eukaryotic group to modulate select bacteria in their microbial consortia, similar to higher eukaryotes, using unique secondary metabolites that regulate bacterial growth and behavior inversely across different bacterial populations. Significance Phytoplankton are major primary producers in the marine environment that excrete a wide range of metabolites. These exudates support the growth of surrounding bacteria that in turn provide phytoplankton cells with resources and growth cofactors to proliferate. Unlike multicellular eukaryotes with dedicated structures housing microbiomes, mechanisms enabling unicellular eukaryotic phytoplankton hosts to modulate potential symbionts and opportunists within a natural microbial community are unknown. Exposure of a host phytoplankton cell to its natural microbial community triggers major transcriptional and metabolic reprogramming to release unique secondary metabolites that selectively enable growth and attachment of symbiotic taxa while simultaneously suppressing the colonization of nonsymbiont bacteria. These results suggest strong and highly selective microbiome-modulating strategies shared across the unicellular and multicellular eukaryotic lineages.

Coral reefs are threatened throughout the world. A major factor contributing to their decline is outbreaks and propagation of coral diseases. Due to the complexity of coral-associated microbe communities, little is understood in terms of disease agents, hosts and vectors. It is known that compromised health in corals is correlated with shifts in bacterial assemblages colonizing coral mucus and tissue. However, general disease patterns remain, to a large extent, ambiguous as comparative studies over species, regions, or diseases are scarce. Here, we compare bacterial assemblages of samples from healthy (HH) colonies and such displaying signs of White Plague Disease (WPD) of two different coral species (Pavona duerdeni and Porites lutea) from the same reef in Koh Tao, Thailand, using 16S rRNA gene microarrays. In line with other studies, we found an increase of bacterial diversity in diseased (DD) corals, and a higher abundance of taxa from the families that include known coral pathogens (Alteromonadaceae, Rhodobacteraceae, Vibrionaceae). In our comparative framework analysis, we found differences in microbial assemblages between coral species and coral health states. Notably, patterns of bacterial community structures from HH and DD corals were maintained over species boundaries. Moreover, microbes that differentiated the two coral species did not overlap with microbes that were indicative of HH and DD corals. This suggests that while corals harbor distinct species-specific microbial assemblages, disease-specific bacterial abundance patterns exist that are maintained over coral species boundaries.

Summary: Interactions between phytoplankton and bacteria play major roles in global biogeochemical cycles and oceanic nutrient fluxes. These interactions occur in the microenvironment surrounding phytoplankton cells, known as the phycosphere. Bacteria in the phycosphere use either chemotaxis or attachment to benefit from algal excretions. Both processes are regulated by quorum sensing (QS), a cell–cell signalling mechanism that uses small infochemicals to coordinate bacterial gene expression. However, the role of QS in regulating bacterial attachment in the phycosphere is not clear. Here, we isolated a Sulfitobacter pseudonitzschiae F5 and a Phaeobacter sp. F10 belonging to the marine Roseobacter group and an Alteromonas macleodii F12 belonging to Alteromonadaceae, from the microbial community of the ubiquitous diatom Asterionellopsis glacialis. We show that only the Roseobacter group isolates (diatom symbionts) can attach to diatom transparent exopolymeric particles. Despite all three bacteria possessing genes involved in motility, chemotaxis, and attachment, only S. pseudonitzschiae F5 and Phaeobacter sp. F10 possessed complete QS systems and could synthesize QS signals. Using UHPLC–MS/MS, we identified three QS molecules produced by both bacteria of which only 3‐oxo‐C _16:1 ‐HSL strongly inhibited bacterial motility and stimulated attachment in the phycosphere. These findings suggest that QS signals enable colonization of the phycosphere by algal symbionts.

Oceanic microbial diversity covaries with physicochemical parameters. Temperature, for example, explains approximately half of global variation in surface taxonomic abundance. It is unknown, however, whether covariation patterns hold over narrower parameter gradients and spatial scales, and extending to mesopelagic depths. We collected and sequenced 45 epipelagic and mesopelagic microbial metagenomes on a meridional transect through the eastern Red Sea. We asked which environmental parameters explain the most variation in relative abundances of taxonomic groups, gene ortholog groups, and pathways—at a spatial scale of <2000 km, along narrow but well-defined latitudinal and depth-dependent gradients. We also asked how microbes are adapted to gradients and extremes in irradiance, temperature, salinity, and nutrients, examining the responses of individual gene ortholog groups to these parameters. Functional and taxonomic metrics were equally well explained (75–79%) by environmental parameters. However, only functional and not taxonomic covariation patterns were conserved when comparing with an intruding water mass with different physicochemical properties. Temperature explained the most variation in each metric, followed by nitrate, chlorophyll, phosphate, and salinity. That nitrate explained more variation than phosphate suggested nitrogen limitation, consistent with low surface N:P ratios. Covariation of gene ortholog groups with environmental parameters revealed patterns of functional adaptation to the challenging Red Sea environment: high irradiance, temperature, salinity, and low nutrients. Nutrient-acquisition gene ortholog groups were anti-correlated with concentrations of their respective nutrient species, recapturing trends previously observed across much larger distances and environmental gradients. This dataset of metagenomic covariation along densely sampled environmental gradients includes online data exploration supplements, serving as a community resource for marine microbial ecology.

This study investigates the diversity of the most abundant primary producer in the under-explored Red Sea. While we found diverse Prochlorococcus lineages, no n

The marine picocyanobacteria Prochlorococcus represent a significant fraction of the global pelagic bacterioplankton community. Specifically, in the surface waters of the Red Sea, they account for around 91% of the phylum Cyanobacteria. Previous work suggested a widespread presence of high-light (HL)-adapted ecotypes in the Red Sea with the occurrence of low-light (LL)-adapted ecotypes at intermediate depths in the water column. To obtain a more comprehensive dataset over a wider biogeographical scope, we used a 454-pyrosequencing approach to analyze the diversity of the Prochlorococcus rpoC1 gene from a total of 113 samples at various depths (up to 500 m) from 45 stations spanning the Red Sea basin from north to south. In addition, we analyzed 45 metagenomes from eight stations using hidden Markov models based on a set of reference Prochlorococcus genomes to 1) estimate the relative abundance of Prochlorococcus based on 16S rRNA gene sequences, and 2) identify and classify rpoC1 sequences as an assessment of the community structure of Prochlorococcus in the northern, central and southern regions of the basin without amplification bias. Analyses of metagenomic data indicated that Prochlorococcus occurs at a relative abundance of around 9% in samples from surface waters (25, 50, 75 m), 3% in intermediate waters (100 m) and around 0.5% in deep-water samples (200 – 500 m). Results based on rpoC1 sequences using both methods showed that HL II cells dominate surface waters and were also present in deep-water samples. Prochlorococcus communities in intermediate waters (100 m) showed a higher diversity and co-occurrence of low-light and high-light ecotypes. Prochlorococcus communities at each depth range (surface, intermediate, deep sea) did not change significantly over the sampled transects spanning most of the Saudi waters in the Red Sea. Statistical analyses of rpoC1 sequences from metagenomes indicated that the vertical distribution of Prochlorococcus in the water column is correlated with physicochemical gradients: temperature and oxygen are positively correlated with HL II (R2= 0.71, p-value ≤ 0.05) while Chl a concentration, nutrient concentrations and salinity correlate with the prevalence of LL clades.