Cyanobacterial distributions along a physico-chemical gradient in the Northeastern Pacific Ocean

Recurring seasonality exposes dominant species and niche partitioning strategies of open ocean picoeukaryotic algae

Ocean spring phytoplankton blooms are dynamic periods important to global primary production. We document vertical patterns of a diverse suite of eukaryotic algae, the prasinophytes, in the North Atlantic Subtropical Gyre with monthly sampling over four years at the Bermuda Atlantic Time-series Study site. Water column structure was used to delineate seasonal stability periods more ecologically relevant than seasons defined by calendar dates. During winter mixing, tiny prasinophytes dominated by Class II comprise 46  ±  24% of eukaryotic algal (plastid-derived) 16S rRNA V1-V2 amplicons, specifically Ostreococcus Clade OII, Micromonas commoda, and Bathycoccus calidus. In contrast, Class VII are rare and Classes I and VI peak during warm stratified periods when surface eukaryotic phytoplankton abundances are low. Seasonality underpins a reservoir of genetic diversity from multiple prasinophyte classes during warm periods that harbor ephemeral taxa. Persistent Class II sub-species dominating the winter/spring bloom period retreat to the deep chlorophyll maximum in summer, poised to seed the mixed layer upon winter convection, exposing a mechanism for initiating high abundances at bloom onset. Comparisons to tropical oceans reveal broad distributions of the dominant sub-species herein. This unparalleled window into temporal and spatial niche partitioning of picoeukaryotic primary producers demonstrates how key prasinophytes prevail in warm oceans.

3D intrusions transport active surface microbial assemblages to the dark ocean

Subtropical oceans contribute significantly to global primary production, but the fate of the picophytoplankton that dominate in these low-nutrient regions is poorly understood. Working in the subtropical Mediterranean, we demonstrate that subduction of water at ocean fronts generates 3D intrusions with uncharacteristically high carbon, chlorophyll, and oxygen that extend below the sunlit photic zone into the dark ocean. These contain fresh picophytoplankton assemblages that resemble the photic-zone regions where the water originated. Intrusions propagate depth-dependent seasonal variations in microbial assemblages into the ocean interior. Strikingly, the intrusions included dominant biomass contributions from nonphotosynthetic bacteria and enrichment of enigmatic heterotrophic bacterial lineages. Thus, the intrusions not only deliver material that differs in composition and nutritional character from sinking detrital particles, but also drive shifts in bacterial community composition, organic matter processing, and interactions between surface and deep communities. Modeling efforts paired with global observations demonstrate that subduction can flux similar magnitudes of particulate organic carbon as sinking export, but is not accounted for in current export estimates and carbon cycle models. Intrusions formed by subduction are a particularly important mechanism for enhancing connectivity between surface and upper mesopelagic ecosystems in stratified subtropical ocean environments that are expanding due to the warming climate.

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

A single Prochlorococcus ecotype dominates the tropical Bay of Bengal with ultradian growth

The Bay of Bengal (BoB) spans >2.2 million km2 in the northeastern Indian Ocean and is bordered by dense populations that depend upon its resources. Over recent decades, a shift from larger phytoplankton to picoplankton has been reported, yet the abundance, activity, and composition of primary producer communities are not well-characterized. We analysed the BoB regions during the summer monsoon. Prochlorococcus ranged up to 3.14 × 105 cells mL-1 in the surface mixed layer, averaging 1.74 ± 0.46 × 105 in the upper 10 m and consistently higher than Synechococcus and eukaryotic phytoplankton. V1-V2 rRNA gene amplicon analyses showed the High Light II (HLII) ecotype formed 98 ± 1% of Prochlorococcus amplicons in surface waters, comprising six oligotypes, with the dominant oligotype accounting for 65 ± 4% of HLII. Diel sampling of a coherent water mass demonstrated evening onset of cell division and rapid Prochlorococcus growth between 1.5 and 3.1 div day-1, based on cell cycle analysis, as confirmed by abundance-based estimates of 2.1 div day-1. Accumulation of Prochlorococcus produced by ultradian growth was restricted by high loss rates. Alongside prior Arabian Sea and tropical Atlantic rates, our results indicate Prochlorococcus growth rates should be reevaluated with greater attention to latitudinal zones and influences on contributions to global primary production.

The Bay of Bengal exposes abundant photosynthetic picoplankton and newfound diversity along salinity-driven gradients

The Bay of Bengal (BoB) is a 2,600,000 km2 expanse in the Indian Ocean upon which many humans rely. However, the primary producers underpinning food chains here remain poorly characterized. We examined phytoplankton abundance and diversity along strong BoB latitudinal and vertical salinity gradients-which have low temperature variation (27-29°C) between the surface and subsurface chlorophyll maximum (SCM). In surface waters, Prochlorococcus averaged 11.7 ± 4.4 × 104  cells ml-1 , predominantly HLII, whereas LLII and 'rare' ecotypes, HLVI and LLVII, dominated in the SCM. Synechococcus averaged 8.4 ± 2.3 × 104  cells ml-1 in the surface, declined rapidly with depth, and population structure of dominant Clade II differed between surface and SCM; Clade X was notable at both depths. Across all sites, Ostreococcus Clade OII dominated SCM eukaryotes whereas communities differentiated strongly moving from Arabian Sea-influenced high salinity (southerly; prasinophytes) to freshwater-influenced low salinity (northerly; stramenopiles, specifically, diatoms, pelagophytes, and dictyochophytes, plus the prasinophyte Micromonas) surface waters. Eukaryotic phytoplankton peaked in the south (1.9 × 104  cells ml-1 , surface) where a novel Ostreococcus was revealed, named here Ostreococcus bengalensis. We expose dominance of a single picoeukaryote and hitherto 'rare' picocyanobacteria at depth in this complex ecosystem where studies suggest picoplankton are replacing larger phytoplankton due to climate change.

Differential Timing for Glucose Assimilation in Prochlorococcus and Coexistent Microbial Populations in the North Pacific Subtropical Gyre

The marine cyanobacterium Prochlorococcus can utilize glucose as a source of carbon. However, the relative importance of inorganic and organic carbon assimilation and the timing of glucose assimilation are still poorly understood in these numerically dominant cyanobacteria. Here, we investigated whole microbial community and group-specific primary production and glucose assimilation using incubations with radioisotopes combined with flow cytometry cell sorting. We also studied changes in the microbial community structure in response to glucose enrichments and analyzed the transcription of Prochlorocccus genes involved in carbon metabolism and photosynthesis. Our results showed a diel variation for glucose assimilation in Prochlorococcus, with maximum assimilation at midday and minimum at midnight (~2-fold change), which was different from that of the total microbial community. This suggests that the timing in glucose assimilation in Prochlorococcus is coupled to photosynthetic light reactions producing energy, it being more convenient for Prochlorococcus to show maximum glucose uptake precisely when the rest of microbial populations have their minimum glucose uptake. Many transcriptional responses to glucose enrichment occurred after 12- and 24-h periods, but community composition did not change. High-light Prochlorococcus strains were the most impacted by glucose addition, with transcript-level increases observed for genes in pathways for glucose metabolism, such as the pentose phosphate pathway, the Entner-Doudoroff pathway, glycolysis, respiration, and glucose transport. While Prochlorococcus C assimilation from glucose represented less than 0.1% of the bacterium's photosynthetic C fixation, increased assimilation during the day and glcH gene upregulation upon glucose enrichment indicate an important role of mixotrophic C assimilation by natural populations of Prochlorococcus. IMPORTANCE Several studies have demonstrated that Prochlorococcus, the most abundant photosynthetic organism on Earth, can assimilate organic molecules, such as amino acids, amino sugars, ATP, phosphonates, and dimethylsulfoniopropionate. This autotroph can also assimilate small amounts of glucose, supporting the hypothesis that Prochlorococcus is mixotrophic. Our results show, for the first time, a diel variability in glucose assimilation by natural populations of Prochlorococcus with maximum assimilation during midday. Based on our previous results, this indicates that Prochlorococcus could maximize glucose uptake by using ATP made during the light reactions of photosynthesis. Furthermore, Prochlorococcus showed a different timing of glucose assimilation from the total population, which may offer considerable fitness advantages over competitors "temporal niches." Finally, we observed transcriptional changes in some of the genes involved in carbon metabolism, suggesting that Prochlorococcus can use both pathways previously proposed in cyanobacteria to metabolize glucose.

Spatiotemporal Variations in Antarctic Protistan Communities Highlight Phytoplankton Diversity and Seasonal Dominance by a Novel Cryptophyte Lineage

The Andvord fjord in the West Antarctic Peninsula (WAP) is known for its productivity and abundant megafauna. Nevertheless, seasonal patterns of the molecular diversity and abundance of protistan community members underpinning WAP productivity remain poorly resolved. We performed spring and fall expeditions pursuing protistan diversity, abundance of photosynthetic taxa, and the connection to changing conditions. 18S rRNA amplicon sequence variant (ASV) profiles revealed diverse predatory protists spanning multiple eukaryotic supergroups, alongside enigmatic heterotrophs like the Picozoa. Among photosynthetic protists, cryptophyte contributions were notable. Analysis of plastid-derived 16S rRNA ASVs supported 18S ASV results, including a dichotomy between cryptophytes and diatom contributions previously reported in other Antarctic regions. We demonstrate that stramenopile and cryptophyte community structures have distinct attributes. Photosynthetic stramenopiles exhibit high diversity, with the polar diatom Fragilariopsis cylindrus, unidentified Chaetoceros species, and others being prominent. Conversely, ASV analyses followed by environmental full-length rRNA gene sequencing, electron microscopy, and flow cytometry revealed that a novel alga dominates the cryptophytes. Phylogenetic analyses established that TPG clade VII, as named here, is evolutionarily distinct from cultivated cryptophyte lineages. Additionally, cryptophyte cell abundance correlated with increased water temperature. Analyses of global data sets showed that clade VII dominates cryptophyte ASVs at Southern Ocean sites and appears to be endemic, whereas in the Arctic and elsewhere, Teleaulax amphioxeia and Plagioselmis prolonga dominate, although both were undetected in Antarctic waters. Collectively, our studies provide baseline data against which future change can be assessed, identify different diversification patterns between stramenopiles and cryptophytes, and highlight an evolutionarily distinct cryptophyte clade that thrives under conditions enhanced by warming.

Selective Uptake of Pelagic Microbial Community Members by Caribbean Reef Corals

We identify interactions between coral grazing behavior and the growth rates and cell abundances of pelagic microbial groups found surrounding a Caribbean reef. During incubation experiments with three reef corals, reductions in microbial cell abundance differed according to coral species and suggest specific coral or microbial mechanisms are at play.

Coral reefs are possible sinks for microbes; however, the removal mechanisms at play are not well understood. Here, we characterize pelagic microbial groups at the CARMABI reef (Curaçao) and examine microbial consumption by three coral species: Madracis mirabilis, Porites astreoides, and Stephanocoenia intersepta. Flow cytometry analyses of water samples collected from a depth of 10 m identified 6 microbial groups: Prochlorococcus, three groups of Synechococcus, photosynthetic eukaryotes, and heterotrophic bacteria. Minimum growth rates (μ) for Prochlorococcus, all Synechococcus groups, and photosynthetic eukaryotes were 0.55, 0.29, and 0.45 μ day⁻¹, respectively, and suggest relatively high rates of productivity despite low nutrient conditions on the reef. During a series of 5-h incubations with reef corals performed just after sunset or prior to sunrise, reductions in the abundance of photosynthetic picoeukaryotes, Prochlorococcus and Synechococcus cells, were observed. Of the three Synechococcus groups, one decreased significantly during incubations with each coral and the other two only with M. mirabilis. Removal of carbon from the water column is based on coral consumption rates of phytoplankton and averaged between 138 ng h⁻¹ and 387 ng h⁻¹, depending on the coral species. A lack of coral-dependent reduction in heterotrophic bacteria, differences in Synechococcus reductions, and diurnal variation in reductions of Synechococcus and Prochlorococcus, coinciding with peak cell division, point to selective feeding by corals. Our study indicates that bentho-pelagic coupling via selective grazing of microbial groups influences carbon flow and supports heterogeneity of microbial communities overlying coral reefs.

IMPORTANCE We identify interactions between coral grazing behavior and the growth rates and cell abundances of pelagic microbial groups found surrounding a Caribbean reef. During incubation experiments with three reef corals, reductions in microbial cell abundance differed according to coral species and suggest specific coral or microbial mechanisms are at play. Peaks in removal rates of Prochlorococcus and Synechococcus cyanobacteria appear highest during postsunset incubations and coincide with microbial cell division. Grazing rates and effort vary across coral species and picoplankton groups, possibly influencing overall microbial composition and abundance over coral reefs. For reef corals, use of such a numerically abundant source of nutrition may be advantageous, especially under environmentally stressful conditions when symbioses with dinoflagellate algae break down.

Seasonal and Geographical Transitions in Eukaryotic Phytoplankton Community Structure in the Atlantic and Pacific Oceans

Much is known about how broad eukaryotic phytoplankton groups vary according to nutrient availability in marine ecosystems. However, genus- and species-level dynamics are generally unknown, although important given that adaptation and acclimation processes differentiate at these levels. We examined phytoplankton communities across seasonal cycles in the North Atlantic (BATS) and under different trophic conditions in the eastern North Pacific (ENP), using phylogenetic classification of plastid-encoded 16S rRNA amplicon sequence variants (ASVs) and other methodologies, including flow cytometric cell sorting. Prasinophytes dominated eukaryotic phytoplankton amplicons during the nutrient-rich deep-mixing winter period at BATS. During stratification (‘summer’) uncultured dictyochophytes formed ∼35 ± 10% of all surface plastid amplicons and dominated those from stramenopile algae, whereas diatoms showed only minor, ephemeral contributions over the entire year. Uncultured dictyochophytes also comprised a major fraction of plastid amplicons in the oligotrophic ENP. Phylogenetic reconstructions of near-full length 16S rRNA sequences established 11 uncultured Dictyochophyte Environmental Clades (DEC). DEC-I and DEC-VI dominated surface dictyochophytes under stratification at BATS and in the ENP, and DEC-IV was also important in the latter. Additionally, although less common at BATS, Florenciella-related clades (FC) were prominent at depth in the ENP. In both ecosystems, pelagophytes contributed notably at depth, with PEC-VIII (Pelagophyte Environmental Clade) and (cultured) Pelagomonas calceolata being most important. Q-PCR confirmed the near absence of P. calceolata at the surface of the same oligotrophic sites where it reached ∼1,500 18S rRNA gene copies ml^–1 at the DCM. To further characterize phytoplankton present in our samples, we performed staining and at-sea single-cell sorting experiments. Sequencing results from these indicated several uncultured dictyochophyte clades are comprised of predatory mixotrophs. From an evolutionary perspective, these cells showed both conserved and unique features in the chloroplast genome. In ENP metatranscriptomes we observed high expression of multiple chloroplast genes as well as expression of a selfish element (group II intron) in the psaA gene. Comparative analyses across the Pacific and Atlantic sites support the conclusion that predatory dictyochophytes thrive under low nutrient conditions. The observations that several uncultured dictyochophyte lineages are seemingly capable of photosynthesis and predation, raises questions about potential shifts in phytoplankton trophic roles associated with seasonality and long-term ocean change.

Phytoplankton transcriptomic and physiological responses to fixed nitrogen in the California current system

Marine phytoplankton are responsible for approximately half of photosynthesis on Earth. However, their ability to drive ocean productivity depends on critical nutrients, especially bioavailable nitrogen (N) which is scarce over vast areas of the ocean. Phytoplankton differ in their preferences for N substrates as well as uptake efficiencies and minimal N requirements relative to other critical nutrients, including iron (Fe) and phosphorus. In this study, we used the MicroTOOLs high-resolution environmental microarray to examine transcriptomic responses of phytoplankton communities in the California Current System (CCS) transition zone to added urea, ammonium, nitrate, and also Fe in the late summer when N depletion is common. Transcript level changes of photosynthetic, carbon fixation, and nutrient stress genes indicated relief of N limitation in many strains of Prochlorococcus, Synechococcus, and eukaryotic phytoplankton. The transcriptomic responses helped explain shifts in physiological and growth responses observed later. All three phytoplankton groups had increased transcript levels of photosynthesis and/or carbon fixation genes in response to all N substrates. However, only Prochlorococcus had decreased transcript levels of N stress genes and grew substantially, specifically after urea and ammonium additions, suggesting that Prochlorococcus outcompeted other community members in these treatments. Diatom transcript levels of carbon fixation genes increased in response to Fe but not to Fe with N which might have favored phytoplankton that were co-limited by N and Fe. Moreover, transcription patterns of closely related strains indicated variability in N utilization, including nitrate utilization by some high-light adapted Prochlorococcus. Finally, up-regulation of urea transporter genes by both Prochlorococcus and Synechococcus in response to filtered deep water suggested a regulatory mechanism other than classic control via the global N regulator NtcA. This study indicated that co-existing phytoplankton strains experience distinct nutrient stresses in the transition zone of the CCS, an understudied region where oligotrophic and coastal communities naturally mix.

Small phytoplankton dominate western North Atlantic biomass

The North Atlantic phytoplankton spring bloom is the pinnacle in an annual cycle that is driven by physical, chemical, and biological seasonality. Despite its important contributions to the global carbon cycle, transitions in plankton community composition between the winter and spring have been scarcely examined in the North Atlantic. Phytoplankton composition in early winter was compared with latitudinal transects that captured the subsequent spring bloom climax. Amplicon sequence variants (ASVs), imaging flow cytometry, and flow-cytometry provided a synoptic view of phytoplankton diversity. Phytoplankton communities were not uniform across the sites studied, but rather mapped with apparent fidelity onto subpolar- and subtropical-influenced water masses of the North Atlantic. At most stations, cells < 20-µm diameter were the main contributors to phytoplankton biomass. Winter phytoplankton communities were dominated by cyanobacteria and pico-phytoeukaryotes. These transitioned to more diverse and dynamic spring communities in which pico- and nano-phytoeukaryotes, including many prasinophyte algae, dominated. Diatoms, which are often assumed to be the dominant phytoplankton in blooms, were contributors but not the major component of biomass. We show that diverse, small phytoplankton taxa are unexpectedly common in the western North Atlantic and that regional influences play a large role in modulating community transitions during the seasonal progression of blooms.

Diel transcriptional response of a California Current plankton microbiome to light, low iron, and enduring viral infection

Phytoplankton and associated microbial communities provide organic carbon to oceanic food webs and drive ecosystem dynamics. However, capturing those dynamics is challenging. Here, an in situ, semi-Lagrangian, robotic sampler profiled pelagic microbes at 4 h intervals over ~2.6 days in North Pacific high-nutrient, low-chlorophyll waters. We report on the community structure and transcriptional dynamics of microbes in an operationally large size class (>5 μm) predominantly populated by dinoflagellates, ciliates, haptophytes, pelagophytes, diatoms, cyanobacteria (chiefly Synechococcus), prasinophytes (chiefly Ostreococcus), fungi, archaea, and proteobacteria. Apart from fungi and archaea, all groups exhibited 24-h periodicity in some transcripts, but larger portions of the transcriptome oscillated in phototrophs. Periodic photosynthesis-related transcripts exhibited a temporal cascade across the morning hours, conserved across diverse phototrophic lineages. Pronounced silica:nitrate drawdown, a high flavodoxin to ferredoxin transcript ratio, and elevated expression of other Fe-stress markers indicated Fe-limitation. Fe-stress markers peaked during a photoperiodically adaptive time window that could modulate phytoplankton response to seasonal Fe-limitation. Remarkably, we observed viruses that infect the majority of abundant taxa, often with total transcriptional activity synchronized with putative hosts. Taken together, these data reveal a microbial plankton community that is shaped by recycled production and tightly controlled by Fe-limitation and viral activity.

Distribution of Abundant and Active Planktonic Ciliates in Coastal and Slope Waters Off New England

Despite their important role of linking microbial and classic marine food webs, data on biogeographical patterns of microbial eukaryotic grazers are limited, and even fewer studies have used molecular tools to assess active (i.e., those expressing genes) community members. Marine ciliate diversity is believed to be greatest at the chlorophyll maximum, where there is an abundance of autotrophic prey, and is often assumed to decline with depth. Here, we assess the abundant (DNA) and active (RNA) marine ciliate communities throughout the water column at two stations off the New England coast (Northwest Atlantic)—a coastal station 43 km from shore (40 m depth) and a slope station 135 km off shore (1,000 m). We analyze ciliate communities using a DNA fingerprinting technique, Denaturing Gradient Gel Electrophoresis (DGGE), which captures patterns of abundant community members. We compare estimates of ciliate communities from SSU-rDNA (abundant) and SSU-rRNA (active) and find complex patterns throughout the water column, including many active lineages below the photic zone. Our analyses reveal (1) a number of widely-distributed taxa that are both abundant and active; (2) considerable heterogeneity in patterns of presence/absence of taxa in offshore samples taken 50 m apart throughout the water column; and (3) three distinct ciliate assemblages based on position from shore and depth. Analysis of active (RNA) taxa uncovers biodiversity hidden to traditional DNA-based approaches (e.g., clone library, rDNA amplicon studies).

Synechococcus Assemblages across the Salinity Gradient in a Salt Wedge Estuary

Synechococcus are the most abundant and widely distributed picocyanobacteria in the ocean. The salt-wedge type of estuary possesses the complete horizontal and vertical gradient of salinity together with other physical and chemical parameters. In order to reveal whether such a complex environmental gradient harbors a high diversity of Synechococcus, we investigated the abundance, taxonomic composition and pigment genetic diversity of Synechococcus in surface and bottom waters across the salinity gradient in a salt-wedge estuary by flow cytometric analysis and pyrosequencing of the rpoC1 gene and cpcBA operon (encoding phycocyanin). Synechococcus were ubiquitously distributed in the studied region, with clear spatial variations both horizontally and vertically. The abundance and diversity of Synechococcus were low in the freshwater-dominated low salinity waters. By pyrosequencing of the rpoC1 gene, we have shown that with the increase of salinity, the dominant Synechococcus shifted from the freshwater Synechococcus to the combination of phylogenetic subcluster 5.2 and freshwater Synechococcus, and then the strictly marine subcluster 5.1 clade III. Besides, the composition of Synechococcus assemblage in the deep layer was markedly different from the surface in the stratified waters (dissimilarities: 40.32%-95.97%, SIMPER analysis). High abundance of clade III Synechococcus found in the brackish waters may revise our previous understanding that strains of this clade prefers oligotrophic environment. Our data also suggested that both the phylogenetic subcluster 5.3 Synechococcus, a lineage that was not well understood, and subcluster 5.1 clade I, a typical cold water lineage, were widely distributed in the bottom layer of the estuary. Clade I detected in the studied region was mainly contributed by subclade IG. Analysis of the cpcBA operon sequences revealed niche partitioning between type 1 and type 3 Synechococcus, with type 2 distributed broadly across the whole environmental gradients. Our results suggest that the salt wedge estuary provides various niches for different lineages of Synechococcus, making it an environment with high Synechococcus diversity compared with adjacent freshwater and shelf sea environments.

Genomic differentiation among wild cyanophages despite widespread horizontal gene transfer

Background

Genetic recombination is a driving force in genome evolution. Among viruses it has a dual role. For genomes with higher fitness, it maintains genome integrity in the face of high mutation rates. Conversely, for genomes with lower fitness, it provides immediate access to sequence space that cannot be reached by mutation alone. Understanding how recombination impacts the cohesion and dissolution of individual whole genomes within viral sequence space is poorly understood across double-stranded DNA bacteriophages (a.k.a phages) due to the challenges of obtaining appropriately scaled genomic datasets.

Results

Here we explore the role of recombination in both maintaining and differentiating whole genomes of 142 wild double-stranded DNA marine cyanophages. Phylogenomic analysis across the 51 core genes revealed ten lineages, six of which were well represented. These phylogenomic lineages represent discrete genotypic populations based on comparisons of intra- and inter- lineage shared gene content, genome-wide average nucleotide identity, as well as detected gaps in the distribution of pairwise differences between genomes. McDonald-Kreitman selection tests identified putative niche-differentiating genes under positive selection that differed across the six well-represented genotypic populations and that may have driven initial divergence. Concurrent with patterns of recombination of discrete populations, recombination analyses of both genic and intergenic regions largely revealed decreased genetic exchange across individual genomes between relative to within populations.

Conclusions

These findings suggest that discrete double-stranded DNA marine cyanophage populations occur in nature and are maintained by patterns of recombination akin to those observed in bacteria, archaea and in sexual eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3286-x) contains supplementary material, which is available to authorized users.

Variable but persistent coexistence of Prochlorococcus ecotypes along temperature gradients in the ocean's surface mixed layer

The vast majority of the phytoplankton communities in surface mixed layer of the oligotrophic ocean are numerically dominated by one of two ecotypes of Prochlorococcus, eMIT9312 or eMED4. In this study, we surveyed large latitudinal transects in the Atlantic and Pacific Ocean to determine if these ecotypes discretely partition the surface mixed layer niche, or if populations exist as a continuum along key environmental gradients, particularly temperature. Transitions of dominance occurred at approximately 19-21°C, with the eMED4 ecotype dominating the colder, and eMIT9312 ecotype dominating the warmer regions. Within these zones of regional dominance, however, the minority ecotype was not competed to extinction. Rather, a robust log-linear relationship between ecotype ratio and temperature characterized this stabilized coexistence: for every 2.5°C increase in temperature, the eMIT9312:eMED4 ratio increased by an order of magnitude. This relationship was observed in both quantitative polymerase chain reaction and in pyrosequencing assays. Water column stratification also contributed to the ecotype ratio along the basin-scale transects, but to a lesser extent. Finally, instances where the ratio of the eMED4 and eMIT9312 abundances did not correlate well with temperature were identified. Such occurrences are likely due to changes in water temperatures outpacing changes in community structure.

Abundance and Biogeography of Picoprasinophyte Ecotypes and Other Phytoplankton in the Eastern North Pacific Ocean

Eukaryotic algae within the picoplankton size class (≤2 μm in diameter) are important marine primary producers, but their spatial and ecological distributions are not well characterized. Here, we studied three picoeukaryotic prasinophyte genera and their cyanobacterial counterparts, Prochlorococcus and Synechococcus, during two cruises along a North Pacific transect characterized by different ecological regimes. Picoeukaryotes and Synechococcus reached maximum abundances of 1.44 × 10(5) and 3.37 × 10(5) cells · ml(-1), respectively, in mesotrophic waters, while Prochlorococcus reached 1.95 × 10(5) cells · ml(-1) in the oligotrophic ocean. Of the picoeukaryotes, Bathycoccus was present at all stations in both cruises, reaching 21,368 ± 327 18S rRNA gene copies · ml(-1). Micromonas and Ostreococcus clade OI were detected only in mesotrophic and coastal waters and Ostreococcus clade OII only in the oligotrophic ocean. To resolve proposed Bathycoccus ecotypes, we established genetic distances for 1,104 marker genes using targeted metagenomes and the Bathycoccus prasinos genome. The analysis was anchored in comparative genome analysis of three Ostreococcus species for which physiological and environmental data are available to facilitate data interpretation. We established that two Bathycoccus ecotypes exist, named here BI (represented by coastal isolate Bathycoccus prasinos) and BII. These share 82% ± 6% nucleotide identity across homologs, while the Ostreococcus spp. share 75% ± 8%. We developed and applied an analysis of ecomarkers to metatranscriptomes sequenced here and published -omics data from the same region. The results indicated that the Bathycoccus ecotypes cooccur more often than Ostreococcus clades OI and OII do. Exploratory analyses of relative transcript abundances suggest that Bathycoccus NRT2.1 and AMT2.2 are high-affinity NO3 (-) and low-affinity NH4 (+) transporters, respectively, with close homologs in multiple picoprasinophytes. Additionally, in the open ocean, where dissolved iron concentrations were low (0.08 nM), there appeared to be a shift to the use of nickel superoxide dismutases (SODs) from Mn/Fe/Cu SODs closer inshore. Our study documents the distribution of picophytoplankton along a North Pacific ecological gradient and offers new concepts and techniques for investigating their biogeography.

Ecophysiology of uncultivated marine euryarchaea is linked to particulate organic matter

Particles in aquatic environments host distinct communities of microbes, yet the evolution of particle-specialized taxa and the extent to which specialized microbial metabolism is associated with particles is largely unexplored. Here, we investigate the hypothesis that a widely distributed and uncultivated microbial group—the marine group II euryarchaea (MGII)—interacts with living and detrital particulate organic matter (POM) in the euphotic zone of the central California Current System. Using fluorescent in situ hybridization, we verified the association of euryarchaea with POM. We further quantified the abundance and distribution of MGII 16 S ribosomal RNA genes in size-fractionated seawater samples and compared MGII functional capacity in metagenomes from the same fractions. The abundance of MGII in free-living and >3 μm fractions decreased with increasing distance from the coast, whereas MGII abundance in the 0.8–3 μm fraction remained constant. At several offshore sites, MGII abundance was highest in particle fractions, indicating that particle-attached MGII can outnumber free-living MGII under oligotrophic conditions. Compared with free-living MGII, the genome content of MGII in particle-associated fractions exhibits an increased capacity for surface adhesion, transcriptional regulation and catabolism of high molecular weight substrates. Moreover, MGII populations in POM fractions are phylogenetically distinct from and more diverse than free-living MGII. Eukaryotic phytoplankton additions stimulated MGII growth in bottle incubations, providing the first MGII net growth rate measurements. These ranged from 0.47 to 0.54 d⁻¹. However, MGII were not recovered in whole-genome amplifications of flow-sorted picoeukaryotic phytoplankton and heterotrophic nanoflagellates, suggesting that MGII in particle fractions are not physically attached to living POM. Collectively, our results support a linkage between MGII ecophysiology and POM, implying that marine archaea have a role in elemental cycling through interactions with particles.