A distinct, high-affinity, alkaline phosphatase facilitates occupation of P-depleted environments by marine picocyanobacteria.

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus, the two most abundant phototrophs on Earth, thrive in oligotrophic oceanic regions. While it is well known that specific lineages are exquisitely adapted to prevailing in situ light and temperature regimes, much less is known of the molecular machinery required to facilitate occupancy of these low-nutrient environments. Here, we describe a hitherto unknown alkaline phosphatase, Psip1, that has a substantially higher affinity for phosphomonoesters than other well-known phosphatases like PhoA, PhoX, or PhoD and is restricted to clade III Synechococcus and a subset of high light I-adapted Prochlorococcus strains, suggesting niche specificity. We demonstrate that Psip1 has undergone convergent evolution with PhoX, requiring both iron and calcium for activity and likely possessing identical key residues around the active site, despite generally very low sequence homology. Interrogation of metagenomes and transcriptomes from TARA oceans and an Atlantic Meridional transect shows that psip1 is abundant and highly expressed in picocyanobacterial populations from the Mediterranean Sea and north Atlantic gyre, regions well recognized to be phosphorus (P)-deplete. Together, this identifies psip1 as an important oligotrophy-specific gene for P recycling in these organisms. Furthermore, psip1 is not restricted to picocyanobacteria and is abundant and highly transcribed in some α-proteobacteria and eukaryotic algae, suggesting that such a high-affinity phosphatase is important across the microbial taxonomic world to occupy low-P environments.

Long-read powered viral metagenomics in the oligotrophic Sargasso Sea.

Dominant microorganisms of the Sargasso Sea are key drivers of the global carbon cycle. However, associated viruses that shape microbial community structure and function are not well characterised. Here, we combined short and long read sequencing to survey Sargasso Sea phage communities in virus- and cellular fractions at viral maximum (80 m) and mesopelagic (200 m) depths. We identified 2,301 Sargasso Sea phage populations from 186 genera. Over half of the phage populations identified here lacked representation in global ocean viral metagenomes, whilst 177 of the 186 identified genera lacked representation in genomic databases of phage isolates. Viral fraction and cell-associated viral communities were decoupled, indicating viral turnover occurred across periods longer than the sampling period of three days. Inclusion of long-read data was critical for capturing the breadth of viral diversity. Phage isolates that infect the dominant bacterial taxa Prochlorococcus and Pelagibacter, usually regarded as cosmopolitan and abundant, were poorly represented.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size.

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10 k - 100 k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Increased cyanophage infection at the bottom of the euphotic zone, especially in the fall.

Picocyanobacteria contribute greatly to offshore primary production with cells extending through the deep euphotic zone. Literature indicates high viral infection of cyanobacteria in ocean transition zones. We postulate that the bottom of the euphotic zone is a transition zone, where communities transition from phototrophic to aphotic processes. We use single-copy core genes to examine cyanophage to cyanobacteria ratios in cellular metagenomes in the subtropical North Atlantic and Pacific. Cyanophage to cyanobacteria terL/rpoB ratios generally increase to >10 in the deep euphotic zone. As light levels decrease in the fall, Prochlorococcus in the deep euphotic zone experience reduced light levels. We find clear differences between spring (Geotraces GA02) and fall (GA03) in the North Atlantic, with terL/rpoB ratios increasing to >40 in the fall. When examining 23 months of the North Pacific Hawaii Ocean Timeseries, the depth of elevated cyanophage to cyanobacteria ratios in cellular metagenomes negatively correlated with surface photosynthetic radiation (PAR), particularly with the change in PAR, which reflected the season. In fall, all picocyanobacteria ecotypes were found at depths enriched with viruses, while in summer, only low light ecotypes were affected. Thus, we find high cyanophage infection both in the deep euphotic zone and during seasonal transitions.

Global scale phylogeography of functional traits and microdiversity in Prochlorococcus.

Prochlorococcus is the most numerically abundant photosynthetic organism in the surface ocean. The Prochlorococcus high-light and warm-water adapted ecotype (HLII) is comprised of extensive microdiversity, but specific functional differences between microdiverse sub-clades remain elusive. Here we characterized both functional and phylogenetic diversity within the HLII ecotype using Bio-GO-SHIP metagenomes. We found widespread variation in gene frequency connected to local environmental conditions. Metagenome-assembled marker genes and genomes revealed a globally distributed novel HLII haplotype defined by adaptation to chronically low P conditions (HLII-P). Environmental correlation analysis revealed different factors were driving gene abundances verses phylogenetic differences. An analysis of cultured HLII genomes and metagenome-assembled genomes revealed a subclade within HLII, which corresponded to the novel HLII-P haplotype. This work represents the first global assessment of the HLII ecotype's phylogeography and corresponding functional differences. These findings together expand our understanding of how microdiversity structures functional differences and reveals the importance of nutrients as drivers of microdiversity in Prochlorococcus.

Choanoflagellates alongside diverse uncultured predatory protists consume the abundant open-ocean cyanobacterium Prochlorococcus.

Prochlorococcus is a key member of open-ocean primary producer communities. Despite its importance, little is known about the predators that consume this cyanobacterium and make its biomass available to higher trophic levels. We identify potential predators along a gradient wherein Prochlorococcus abundance increased from near detection limits (coastal California) to >200,000 cells mL-1 (subtropical North Pacific Gyre). A replicated RNA-Stable Isotope Probing experiment involving the in situ community, and labeled Prochlorococcus as prey, revealed choanoflagellates as the most active predators of Prochlorococcus, alongside a radiolarian, chrysophytes, dictyochophytes, and specific MAST lineages. These predators were not appropriately highlighted in multiyear conventional 18S rRNA gene amplicon surveys where dinoflagellates and other taxa had highest relative amplicon abundances across the gradient. In identifying direct consumers of Prochlorococcus, we reveal food-web linkages of individual protistan taxa and resolve routes of carbon transfer from the base of marine food webs.

Chitin utilization by marine picocyanobacteria and the evolution of a planktonic lifestyle.

Marine picocyanobacteria Prochlorococcus and Synechococcus, the most abundant photosynthetic cells in the oceans, are generally thought to have a primarily single-celled and free-living lifestyle. However, while studying the ability of picocyanobacteria to supplement photosynthetic carbon fixation with the use of exogenous organic carbon, we found the widespread occurrence of genes for breaking down chitin, an abundant source of organic carbon that exists primarily as particles. We show that cells that encode a chitin degradation pathway display chitin degradation activity, attach to chitin particles, and show enhanced growth under low light conditions when exposed to chitosan, a partially deacetylated soluble form of chitin. Marine chitin is largely derived from arthropods, which underwent major diversifications 520 to 535 Mya, close to when marine picocyanobacteria are inferred to have appeared in the ocean. Phylogenetic analyses confirm that the chitin utilization trait was acquired at the root of marine picocyanobacteria. Together this leads us to postulate that attachment to chitin particles allowed benthic cyanobacteria to emulate their mat-based lifestyle in the water column, initiating their expansion into the open ocean, seeding the rise of modern marine ecosystems. Subsequently, transitioning to a constitutive planktonic life without chitin associations led to cellular and genomic streamlining along a major early branch within Prochlorococcus. Our work highlights how the emergence of associations between organisms from different trophic levels, and their coevolution, creates opportunities for colonizing new environments. In this view, the rise of ecological complexity and the expansion of the biosphere are deeply intertwined processes.

Understanding opposing predictions of Prochlorococcus in a changing climate.

Statistically derived species distribution models (SDMs) are increasingly used to predict ecological changes on a warming planet. For Prochlorococcus, the most abundant phytoplankton, an established statistical prediction conflicts with dynamical models as they predict large, opposite, changes in abundance. We probe the SDM at various spatial-temporal scales, showing that light and temperature fail to explain both temporal fluctuations and sharp spatial transitions. Strong correlations between changes in temperature and population emerge only at very large spatial scales, as transects pass through transitions between regions of high and low abundance. Furthermore, a two-state model based on a temperature threshold matches the original SDM in the surface ocean. We conclude that the original SDM has little power to predict changes when Prochlorococcus is already abundant, which resolves the conflict with dynamical models. Our conclusion suggests that SDMs should prove efficacy across multiple spatial-temporal scales before being trusted in a changing ocean.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies.

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Draft genomes of three closely related low light-adapted Prochlorococcus.

OBJECTIVES: The marine cyanobacterium Prochlorococcus is a critical part of warm ocean ecosystems and a model for studying microbial evolution and ecology. To expand the representation of this organism's vast wild diversity in sequence collections, we performed a set of isolation efforts targeting low light-adapted Prochlorococcus. Three genomes resulting from this larger body of work are described here. DATA DESCRIPTION: We present draft-quality Prochlorococcus genomes from enrichment cultures P1344, P1361, and P1363, sampled in the North Pacific. The genomes were built from Illumina paired reads assembled de novo. Supporting datasets of raw reads, assessments, and sequences from co-enriched heterotrophic marine bacteria are also provided. These three genomes represent members of the low light-adapted LLIV Prochlorococcus clade that are closely related, with 99.9% average nucleotide identity between pairs, yet vary in gene content. Expanding the powerful toolkit of Prochlorococcus genomes, these sequences provide an opportunity to study fine-scale variation and microevolutionary processes.

Growth dynamics and transcriptional responses of a Red Sea Prochlorococcus strain to varying temperatures.

Prochlorococcus play a crucial role in the ocean's biogeochemical cycling, but it remains controversial how they will respond to global warming. Here we assessed the response to temperature (22-30°C) of the growth dynamics and gene expression profiles of a Red Sea Prochlorococcus strain (RSP50) in a non-axenic culture. Both the specific growth rate (0.55-0.80 day-1 ) and cell size (0.04-0.07 µm3 ) of Prochlorococcus increased significantly with temperature. The primary production released extracellularly ranged from 20% to 34%, with humic-like fluorescent compounds increasing up to fivefold as Prochlorococcus reached its maximum abundance. At 30°C, genes involved in carbon fixation such as CsoS2 and CsoS3 and photosynthetic electron transport including PTOX were downregulated, suggesting a cellular homeostasis and energy saving mechanism response. In contrast, PTOX was found upregulated at 22°C and 24°C. Similar results were found for transaldolase, related to carbon metabolism, and citrate synthase, an important enzyme in the TCA cycle. Our data suggest that in spite of the currently warm temperatures of the Red Sea, Prochlorococcus can modulate its gene expression profiles to permit growth at temperatures lower than its optimum temperature (28°C) but is unable to cope with temperatures exceeding 30°C.

Production, homology modeling and mutagenesis studies on GlcH glucose transporter from Prochlorococcus sp. strain SS120.

The marine cyanobacterium Prochlorococcus is one of the main primary producers on Earth, which can take up glucose by using the high affinity, multiphasic transporter GlcH. We report here the overexpression of glcH from Prochlorococcus marinus strain SS120 in Escherichia coli. Modeling studies of GlcH using the homologous MelB melibiose transporter from Salmonella enterica serovar Typhimurium showed high conservation at the overall fold. We observed that an important structural interaction, mediated by a strong hydrogen bond between D8 and R141, is conserved in Prochlorococcus, although the corresponding amino acids in MelB from Salmonella are different. Biased docking studies suggested that when glucose reaches the pocket of the transporter and interacts with D8 and R141, the hydrogen bond network in which these residues are involved could be disrupted, favoring a conformational change with the subsequent translocation of the glucose molecule towards the cytoplasmic region of the pmGlcH structure. Based on these theoretical predictions and on the conservation of N117 and W348 in other MelB structures, D8, N117, R141 and W348 were mutated to glycine residues. Their key role in glucose transport was evaluated by glucose uptake assays. N117G and W348G mutations led to 17 % decrease in glucose uptake, while D8G and R141G decreased the glucose transport by 66 % and 92 % respectively. Overall, our studies provide insights into the Prochlorococcus 3D-structure of GlcH, paving the way for further analysis to understand the features which are involved in the high affinity and multiphasic kinetics of this transporter.

Single-cell measurements and modelling reveal substantial organic carbon acquisition by Prochlorococcus.

Marine phytoplankton are responsible for about half of the photosynthesis on Earth. Many are mixotrophs, combining photosynthesis with heterotrophic assimilation of organic carbon, but the relative contribution of these two lifestyles is unclear. Here single-cell measurements reveal that Prochlorococcus at the base of the photic zone in the Eastern Mediterranean Sea obtain only ~20% of carbon required for growth by photosynthesis. This is supported by laboratory-calibrated calculations based on photo-physiology parameters and compared with in situ growth rates. Agent-based simulations show that mixotrophic cells could grow tens of metres deeper than obligate photo-autotrophs, deepening the nutricline by ~20 m. Time series from the North Atlantic and North Pacific indicate that, during thermal stratification, on average 8-10% of the Prochlorococcus cells live without enough light to sustain obligate photo-autotrophic populations. Together, these results suggest that mixotrophy underpins the ecological success of a large fraction of the global Prochlorococcus population and its collective genetic diversity.

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.

Cyanophages from a less virulent clade dominate over their sister clade in global oceans.

Environmental virus communities are highly diverse. However, the infection physiology underlying the evolution of diverse phage lineages and their ecological consequences are largely unknown. T7-like cyanophages are abundant in nature and infect the marine unicellular cyanobacteria, Synechococcus and Prochlorococcus, important primary producers in the oceans. Viruses belonging to this genus are divided into two distinct phylogenetic clades: clade A and clade B. These viruses have narrow host-ranges with clade A phages primarily infecting Synechococcus genotypes, while clade B phages are more diverse and can infect either Synechococcus or Prochlorococcus genotypes. Here we investigated infection properties (life history traits) and environmental abundances of these two clades of T7-like cyanophages. We show that clade A cyanophages have more rapid infection dynamics, larger burst sizes and greater virulence than clade B cyanophages. However, clade B cyanophages were at least 10-fold more abundant in all seasons, and infected more cyanobacteria, than clade A cyanophages in the Red Sea. Models predicted that steady-state cyanophage abundances, infection frequency, and virus-induced mortality, peak at intermediate virulence values. Our findings indicate that differences in infection properties are reflected in virus phylogeny at the clade level. They further indicate that infection properties, together with differences in subclade diversity and host repertoire, have important ecological consequences with the less aggressive, more diverse virus clade having greater ecological impacts.

Microbial Sources of Exocellular DNA in the Ocean.

Exocellular DNA is operationally defined as the fraction of the total DNA pool that passes through a membrane filter (0.1 µm). It is composed of DNA-containing vesicles, viruses, and free DNA and is ubiquitous in all aquatic systems, although the sources, sinks, and ecological consequences are largely unknown. Using a method that provides separation of these three fractions, we compared open ocean depth profiles of DNA associated with each fraction. Pelagibacter-like DNA dominated the vesicle fractions for all samples examined over a depth range of 75 to 500 m. Viral DNA consisted predominantly of myovirus-like and podovirus-like DNA and contained the highest proportion of unannotated sequences. Euphotic zone free DNA (75 to 125 m) contained primarily bacterial and viral sequences, with bacteria dominating samples from the mesopelagic zone (500 to 1,000 m). A high proportion of mesopelagic zone free DNA sequences appeared to originate from surface waters, including a large amount of DNA contributed by high-light Prochlorococcus ecotypes. Throughout the water column, but especially in the mesopelagic zone, the composition of free DNA sequences was not always reflective of cooccurring microbial communities that inhabit the same sampling depth. These results reveal the composition of free DNA in different regions of the water column (euphotic and mesopelagic zones), with implications for dissolved organic matter cycling and export (by way of sinking particles and/or migratory zooplankton) as a delivery mechanism. IMPORTANCE With advances in metagenomic sequencing, the microbial composition of diverse environmental systems has been investigated, providing new perspectives on potential ecological dynamics and dimensions for experimental investigations. Here, we characterized exocellular free DNA via metagenomics, using a newly developed method that separates free DNA from cells, viruses, and vesicles, and facilitated the independent characterization of each fraction. The fate of this free DNA has both ecological consequences as a nutrient (N and P) source and potential evolutionary consequences as a source of genetic transformation. Here, we document different microbial sources of free DNA at the surface (0 to 200 m) versus depths of 250 to 1,000 m, suggesting that distinct free DNA production mechanisms may be present throughout the oligotrophic water column. Examining microbial processes through the lens of exocellular DNA provides insights into the production of labile dissolved organic matter (i.e., free DNA) at the surface (likely by viral lysis) and processes that influence the fate of sinking, surface-derived organic matter.

Diverse Subclade Differentiation Attributed to the Ubiquity of Prochlorococcus High-Light-Adapted Clade II.

Prochlorococcus is the key primary producer in marine ecosystems, and the high-light-adapted clade II (HLII) is the most abundant ecotype. However, the genomic and ecological basis of Prochlorococcus HLII in the marine environment has remained elusive. Here, we show that the ecologically coherent subclade differentiation of HLII corresponds to genomic and ecological characteristics on the basis of analyses of 31 different strains of HLII, including 12 novel isolates. Different subclades of HLII with different core and accessory genes were identified, and their distribution in the marine environment was explored using the TARA Oceans metagenome database. Three major subclade groups were identified, viz., the surface group (HLII-SG), the transition group (HLII-TG), and the deep group (HLII-DG). These subclade groups showed different temperature ranges and optima for distribution. In regression analyses, temperature and nutrient availability were identified as key factors affecting the distribution of HLII subclades. A 35% increase in the relative abundance of HLII-SG by the end of the 21st century was predicted under the Representative Concentration Pathway 8.5 scenario. Our results show that the ubiquity and distribution of Prochlorococcus HLII in the marine environment are associated with the differentiation of diverse subclades. These findings provide insights into the large-scale shifts in the Prochlorococcus community in response to future climate change. IMPORTANCEProchlorococcus is the most abundant oxygenic photosynthetic microorganism on Earth, and high-light-adapted clade II (HLII) is the dominant ecotype. However, the factors behind the dominance of HLII in the vast oligotrophic oceans are still unknown. Here, we identified three distinct groups of HLII subclades, viz., the surface group (HLII-SG), the transition group (HLII-TG), and the deep group (HLII-DG). We further demonstrated that the ecologically coherent subclade differentiation of HLII corresponds to genomic and ecological characteristics. Our study suggests that the differentiation of diverse subclades underlies the ubiquity and distribution of Prochlorococcus HLII in the marine environment and provides insights into the shifts in the Prochlorococcus community in response to future climate change.

Siderophores as an iron source for picocyanobacteria in deep chlorophyll maximum layers of the oligotrophic ocean.

Prochlorococcus and Synechococcus are the most abundant photosynthesizing organisms in the oceans. Gene content variation among picocyanobacterial populations in separate ocean basins often mirrors the selective pressures imposed by the region's distinct biogeochemistry. By pairing genomic datasets with trace metal concentrations from across the global ocean, we show that the genomic capacity for siderophore-mediated iron uptake is widespread in Synechococcus and low-light adapted Prochlorococcus populations from deep chlorophyll maximum layers of iron-depleted regions of the oligotrophic Pacific and S. Atlantic oceans: Prochlorococcus siderophore consumers were absent in the N. Atlantic ocean (higher new iron flux) but constituted up to half of all Prochlorococcus genomes from metagenomes in the N. Pacific (lower new iron flux). Picocyanobacterial siderophore consumers, like many other bacteria with this trait, also lack siderophore biosynthesis genes indicating that they scavenge exogenous siderophores from seawater. Statistical modeling suggests that the capacity for siderophore uptake is endemic to remote ocean regions where atmospheric iron fluxes are the smallest, especially at deep chlorophyll maximum and primary nitrite maximum layers. We argue that abundant siderophore consumers at these two common oceanographic features could be a symptom of wider community iron stress, consistent with prior hypotheses. Our results provide a clear example of iron as a selective force driving the evolution of marine picocyanobacteria.

Phosphonate production by marine microbes: Exploring new sources and potential function.

SignificancePhosphonates are a class of phosphorus metabolites characterized by a highly stable C-P bond. Phosphonates accumulate to high concentrations in seawater, fuel a large fraction of marine methane production, and serve as a source of phosphorus to microbes inhabiting nutrient-limited regions of the oligotrophic ocean. Here, we show that 15% of all bacterioplankton in the surface ocean have genes phosphonate synthesis and that most belong to the abundant groups Prochlorococcus and SAR11. Genomic and chemical evidence suggests that phosphonates are incorporated into cell-surface phosphonoglycoproteins that may act to mitigate cell mortality by grazing and viral lysis. These results underscore the large global biogeochemical impact of relatively rare but highly expressed traits in numerically abundant groups of marine bacteria.

Broad phylogenetic and functional diversity among mixotrophic consumers of Prochlorococcus.

Small eukaryotic phytoplankton are major contributors to global primary production and marine biogeochemical cycles. Many taxa are thought to be mixotrophic, but quantitative studies of phagotrophy exist for very few. In addition, little is known about consumers of Prochlorococcus, the abundant cyanobacterium at the base of oligotrophic ocean food webs. Here we describe thirty-nine new phytoplankton isolates from the North Pacific Subtropical Gyre (Station ALOHA), all flagellates ~2-5 µm diameter, and we quantify their ability to graze Prochlorococcus. The mixotrophs are from diverse classes (dictyochophytes, haptophytes, chrysophytes, bolidophytes, a dinoflagellate, and a chlorarachniophyte), many from previously uncultured clades. Grazing ability varied substantially, with specific clearance rate (volume cleared per body volume) varying over ten-fold across isolates and six-fold across genera. Slower grazers tended to create more biovolume per prey biovolume consumed. Using qPCR we found that the haptophyte Chrysochromulina was most abundant among the isolated mixotrophs at Station ALOHA, with 76-250 cells mL-1 across depths in the upper euphotic zone (5-100 m). Our results show that within a single ecosystem the phototrophs that ingest bacteria come from many branches of the eukaryotic tree, and are functionally diverse, indicating a broad range of strategies along the spectrum from phototrophy to phagotrophy.