Novel lineages of Prochlorococcus thrive within the oxygen minimum zone of the eastern tropical South Pacific

Novel isolates expand the physiological diversity of Prochlorococcus and illuminate its macroevolution

Prochlorococcus is a diverse picocyanobacterial genus and the most abundant phototroph on Earth. Its photosynthetic diversity divides it into high-light (HL)- or low-light (LL)-adapted groups representing broad phylogenetic grades-each composed of several monophyletic clades. Here, we physiologically characterize four new Prochlorococcus strains isolated from below the deep chlorophyll maximum in the North Pacific Ocean. We combine these physiological properties with genomic analyses to explore the evolution of photosynthetic antennae and discuss potential macroevolutionary implications. The isolates belong to deeply branching low-light-adapted clades that have no other cultivated representatives and display some unusual characteristics. For example, despite its otherwise low-light-adapted physiological characteristics, strain MIT1223 has low chl b2 content similar to high-light-adapted strains. Isolate genomes revealed that each strain contains a unique arsenal of pigment biosynthesis and binding alleles that have been horizontally acquired, contributing to the observed physiological diversity. Comparative genomic analysis of all picocyanobacteria reveals that Pcb, the major pigment carrying protein in Prochlorococcus, greatly increased in copy number and diversity per genome along a branch that coincides with the loss of facultative particle attachment. Collectively, these observations support a recently developed macroevolutionary model, in which niche-constructing radiations allowed ancestral lineages of picocyanobacteria to transition from a particle-attached to planktonic lifestyle and broadly colonize the euphotic zone.IMPORTANCEThe marine cyanobacterium, Prochlorococcus, is among the Earth's most abundant organisms, and much of its genetic and physiological diversity remains uncharacterized. Although field studies help reveal the scope of diversity, cultured isolates allow us to link genomic potential to physiological processes, illuminate eco-evolutionary feedbacks, and test theories arising from comparative genomics of wild cells. Here, we report the isolation and characterization of novel low-light (LL)-adapted Prochlorococcus strains that fill in multiple evolutionary gaps. These new strains are the first cultivated representatives of the LLVII and LLVIII paraphyletic grades of Prochlorococcus, which are broadly distributed in the lower regions of the ocean euphotic zone. Each of these grades is a unique, highly diverse section of the Prochlorococcus tree that separates distinct ecological groups: the LLVII grade branches between monophyletic clades that have facultatively particle-associated and constitutively planktonic lifestyles, whereas the LLVIII grade lies along the branch that leads to all high-light (HL)-adapted clades. Characterizing strains and genomes from these grades yields insights into the large-scale evolution of Prochlorococcus. The new LLVII and LLVIII strains are adapted to growth at very low irradiance levels and possess unique light-harvesting gene signatures and pigmentation. The LLVII strains represent the most basal Prochlorococcus group with a major expansion in photosynthetic antenna genes. Furthermore, a strain from the LLVIII grade challenges the paradigm that all LL-adapted Prochlorococcus exhibit high ratios of chl b:a2. These findings provide insights into the photophysiological evolution of Prochlorococcus and redefine what it means to be a low- vs high-light-adapted Prochlorococcus cell.

Prochlorococcus marinus responses to light and oxygen

Prochlorococcus marinus, the smallest picocyanobacterium, comprises multiple clades occupying distinct niches, currently across tropical and sub-tropical oligotrophic ocean regions, including Oxygen Minimum Zones. Ocean warming may open growth-permissive temperatures in new, poleward photic regimes, along with expanded Oxygen Minimum Zones. We used ocean metaproteomic data on current Prochlorococcus marinus niches, to guide testing of Prochlorococcus marinus growth across a matrix of peak irradiances, photoperiods, spectral bands and dissolved oxygen. MED4 from Clade HLI requires greater than 4 h photoperiod, grows at 25 μmol O2 L-1 and above, and exploits high cumulative diel photon doses. MED4, however, relies upon an alternative oxidase to balance electron transport, which may exclude it from growth under our lowest, 2.5 μmol O2 L-1, condition. SS120 from clade LLII/III is restricted to low light under full 250 μmol O2 L-1, shows expanded light exploitation under 25 μmol O2 L-1, but is excluded from growth under 2.5 μmol O2 L-1. Intermediate oxygen suppresses the cost of PSII photoinactivation, and possibly the enzymatic production of H2O2 in SS120, which has limitations on genomic capacity for PSII and DNA repair. MIT9313 from Clade LLIV is restricted to low blue irradiance under 250 μmol O2 L-1, but exploits much higher irradiance under red light, or under lower O2 concentrations, conditions which slow photoinactivation of PSII and production of reactive oxygen species. In warming oceans, range expansions and competition among clades will be governed not only by light levels. Short photoperiods governed by latitude, temperate winters, and depth attenuation of light, will exclude clade HLI (including MED4) from some habitats. In contrast, clade LLII/III (including SS120), and particularly clade LLIV (including MIT9313), may exploit higher light niches nearer the surface, under expanding OMZ conditions, where low O2 relieves the stresses of oxidation stress and PSII photoinhibition.

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.

Associations between picocyanobacterial ecotypes and cyanophage host genes across ocean basins and depth

Background

Cyanophages, viruses that infect cyanobacteria, are globally abundant in the ocean's euphotic zone and are a potentially important cause of mortality for marine picocyanobacteria. Viral host genes are thought to increase viral fitness by either increasing numbers of genes for synthesizing nucleotides for virus replication, or by mitigating direct stresses imposed by the environment. The encoding of host genes in viral genomes through horizontal gene transfer is a form of evolution that links viruses, hosts, and the environment. We previously examined depth profiles of the proportion of cyanophage containing various host genes in the Eastern Tropical North Pacific Oxygen Deficient Zone (ODZ) and at the subtropical North Atlantic (BATS). However, cyanophage host genes have not been previously examined in environmental depth profiles across the oceans.

Methodology

We examined geographical and depth distributions of picocyanobacterial ecotypes, cyanophage, and their viral-host genes across ocean basins including the North Atlantic, Mediterranean Sea, North Pacific, South Pacific, and Eastern Tropical North and South Pacific ODZs using phylogenetic metagenomic read placement. We determined the proportion of myo and podo-cyanophage containing a range of host genes by comparing to cyanophage single copy core gene terminase (terL). With this large dataset (22 stations), network analysis identified statistical links between 12 of the 14 cyanophage host genes examined here with their picocyanobacteria host ecotypes.

Results

Picyanobacterial ecotypes, and the composition and proportion of cyanophage host genes, shifted dramatically and predictably with depth. For most of the cyanophage host genes examined here, we found that the composition of host ecotypes predicted the proportion of viral host genes harbored by the cyanophage community. Terminase is too conserved to illuminate the myo-cyanophage community structure. Cyanophage cobS was present in almost all myo-cyanophage and did not vary in proportion with depth. We used the composition of cobS phylotypes to track changes in myo-cyanophage composition.

Conclusions

Picocyanobacteria ecotypes shift with changes in light, temperature, and oxygen and many common cyanophage host genes shift concomitantly. However, cyanophage phosphate transporter gene pstS appeared to instead vary with ocean basin and was most abundant in low phosphate regions. Abundances of cyanophage host genes related to nutrient acquisition may diverge from host ecotype constraints as the same host can live in varying nutrient concentrations. Myo-cyanophage community in the anoxic ODZ had reduced diversity. By comparison to the oxic ocean, we can see which cyanophage host genes are especially abundant (nirA, nirC, and purS) or not abundant (myo psbA) in ODZs, highlighting both the stability of conditions in the ODZ and the importance of nitrite as an N source to ODZ endemic LLV Prochlorococcus.

Quantification of Marine Picocyanobacteria on Water Column Particles and in Sediments Using Real-Time PCR Reveals Their Role in Carbon Export

Picocyanobacteria are the most abundant primary producers in the ocean and play a fundamental role in marine carbon cycling. Quantification of picocyanobacteria on sinking particles and in sediments is essential to understanding their contribution to the biological carbon pump. We designed a primer set targeting the 16S-23S rRNA internal transcribed spacer (ITS) sequence of cyanobacteria and established a quantitative PCR (qPCR) method for quantifying the ITS sequence abundance. High-throughput sequencing confirmed that this primer set can cover broad diversities of marine picocyanobacteria and avoid amplification of other marine cyanobacteria such as Trichodesmium and Crocosphaera. Amplification efficiencies were slightly different when seven marine Synechococcus and Prochlorococcus strains were assayed. The qPCR results were comparable with flow cytometry for water samples. Using this method, we found that, in the dark ocean, picocyanobacterial ITS sequence abundances were 10 to 100 copies/mL in the size fraction of 0.2 to 3 μm, which were 1 to 3 orders of magnitude more abundant than on the >3-μm particles. We also found that picocyanobacterial ITS abundance in sediment ranged from 10⁵ to 10⁷ copies/g along two nearshore-to-offshore transects in the northern South China Sea. These results further explain the important role of picocyanobacteria in carbon export. Collectively, we provide a qPCR method quantifying the total abundance of marine picocyanobacteria on water column particles and in sediments. Moreover, this newly designed primer set can be also applied to investigate the community of picocyanobacteria via high-throughput sequencing. IMPORTANCE Picocyanobacteria are the most abundant primary producers in the ocean. However, quantification of picocyanobacteria on the sinking particles and in sediments remains challenging using flow cytometry or epifluorescence microscopy. Here, we developed a real-time PCR method to quantify picocyanobacteria using a newly designed primer set specifically targeting the 16S-23S rRNA ITS sequence of cyanobacteria. We showed that in the dark ocean, picocyanobacteria are 1 to 3 orders of magnitude more abundant in small particles (0.2 to 3 μm) than in larger particles (>3 μm). This result supports the important role of direct sinking free-living picocyanobacteria cells in the carbon export to deep ocean. We also found that the picocyanobacterial ITS sequence abundance were 10⁵ to 10⁷ copies per gram in sediments, suggesting significant accumulation of sinking picocyanobacteria in the benthic ecosystem. This qPCR method can be used to quantify the contribution of picocyanobacteria to the biological carbon pump.

The 'oxygen' in oxygen minimum zones

Aerobic processes require oxygen, and anaerobic processes are typically hindered by it. In many places in the global ocean, oxygen is completely removed at mid-water depths forming anoxic oxygen minimum zones (A-OMZs). Within the oxygen gradients linking oxygenated waters with A-OMZs, there is a transition from aerobic to anaerobic microbial processes. This transition is not sharp and there is an overlap between processes using oxygen and those using other electron acceptors. This review will focus on the oxygen control of aerobic and anaerobic metabolisms and will explore how this overlap impacts both the carbon and nitrogen cycles in A-OMZ environments. We will discuss new findings on non-phototrophic microbial processes that produce oxygen, and we focus on how oxygen impacts the loss of fixed nitrogen (as N₂ ) from A-OMZ waters. There are both physiological and environmental controls on the activities of microbial processes responsible for N₂ loss, and the environmental controls are active at extremely low levels of oxygen. Understanding how these controls function will be critical to understanding and predicting how fixed-nitrogen loss in the oceans will respond to future global warming.

Physical-chemical factors influencing the vertical distribution of phototrophic pico-nanoplankton in the Oxygen Minimum Zone (OMZ) off Northern Chile: The relative influence of low pH/low O2 conditions

The vertical distribution of phytoplankton is of fundamental importance in the structure, dynamic, and biogeochemical pathways in marine ecosystems. Nevertheless, what are the main factors determining this distribution remains as an open question. Here, we evaluated the relative influence of environmental factors that might control the coexistence and vertical distribution of pico-nanoplankton associated with the OMZ off northern Chile. Our results showed that in the upper layer Synechococcus-like cells were numerically important at all sampling stations. Pico-nano eukaryotes and phototrophic nanoflagellates (PNF) also showed high abundances in the upper layer decreasing in abundance down to the upper oxycline, while only Prochlorococcus showed high abundances under oxycline and within the oxygen-depleted layer. Statistical analyses evidenced that temperature, oxygen, and carbonate chemistry parameters (pH and dissolved inorganic carbon, DIC) influenced significantly the vertical distribution of phototrophic pico-nanoplankton. Additionally, we experimentally-evaluated the combined effect of low pH/low O₂ conditions on a nanophytoplankton species, the haptophyte Imantonia sp. Under control conditions (pH = 8.1; O₂ = 287.5 μM, light = 169.6 μEm^-2s^-1), Imantonia sp. in vivo fluorescence increased over fifty times, inducing supersaturated O₂ conditions (900 μM) and an increasing pH (8.5), whereas upon an experimental treatment mimicking OMZ conditions (pH = 7.5; O₂ = 55.6 μM; light = 169.6 μEm^-2s^-1), in vivo fluorescence declined dramatically, suggesting that Imantonia sp. did not survive. Although preliminary, our study provides evidence about the role of low pH/low O₂ conditions on the vertical distribution of nanophytoplankton, which deserve future attention through both fieldwork and more extended experimental experiences.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

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.

An analysis of protists in Pacific oxygen deficient zones: implications for Prochlorococcus and N2 -producing bacteria

Ocean oxygen deficient zones (ODZs) host 30%-50% of marine N₂ production. Cyanobacteria photosynthesizing in the ODZ create a secondary chlorophyll maximum and provide organic matter to N₂ -producing bacteria. This chlorophyll maximum is thought to occur due to reduced grazing in anoxic waters. We first examine ODZ protists with long amplicon reads. We then use non-primer-based methods to examine the composition and relative abundance of protists in metagenomes from the Eastern Tropical North and South Pacific ODZs and compare these data to the oxic Hawaii Ocean Time-series (HOT) in the North Pacific. We identify and quantify protists in proportion to the total microbial community. From metagenomic data, we see a large drop in abundance of fungi and protists such as choanoflagellates, radiolarians, cercozoa and ciliates in the ODZs but not in the oxic mesopelagic at HOT. Diplonemid euglenozoa were the only protists that increased in the ODZ. Dinoflagellates and foraminifera reads were also present in the ODZ though less abundant compared to oxic waters. Denitrification has been found in foraminifera but not yet in dinoflagellates. DNA techniques cannot separate dinoflagellate cells and cysts. Metagenomic analysis found taxonomic groups missed by amplicon sequencing and identified trends in abundance.

The Enzymology of Ocean Global Change

A small subset of marine microbial enzymes and surface transporters have a disproportionately important influence on the cycling of carbon and nutrients in the global ocean. As a result, they largely determine marine biological productivity and have been the focus of considerable research attention from microbial oceanographers. Like all biological catalysts, the activity of these keystone biomolecules is subject to control by temperature and pH, leaving the crucial ecosystem functions they support potentially vulnerable to anthropogenic environmental change. We summarize and discuss both consensus and conflicting evidence on the effects of sea surface warming and ocean acidification for five of these critical enzymes [carbonic anhydrase, ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO), nitrogenase, nitrate reductase, and ammonia monooxygenase] and one important transporter (proteorhodopsin). Finally, we forecast how the responses of these few but essential biocatalysts to ongoing global change processes may ultimately help to shape the microbial communities and biogeochemical cycles of the future greenhouse ocean.

Microbial community structures and important taxa across oxygen gradients in the Andaman Sea and eastern Bay of Bengal epipelagic waters

In oceanic oxygen minimum zones (OMZs), the abundances of aerobic organisms significantly decrease and energy shifts from higher trophic levels to microorganisms, while the microbial communities become critical drivers of marine biogeochemical cycling activities. However, little is known of the microbial ecology of the Andaman Sea and eastern Bay of Bengal (BoB) OMZs. In the present study, a total of 131 samples which from the Andaman Sea and eastern BoB epipelagic waters were analyzed. The microbial community distribution patterns across oxygen gradients, including oxygenic zones (OZs, dissolved oxygen [DO] ≥ 2 mg/L), oxygen limited zones (OLZs, 0.7 mg/L < DO < 2 mg/L), and OMZs (DO ≤ 0.7 mg/L), were investigated. Mantel tests and Spearman's correlation analysis revealed that DO was the most important driver of microbial community structures among several environmental factors. Microbial diversity, richness, and evenness were highest in the OLZs and lowest in the OZs. The microbial community compositions of OZ and OMZ waters were significantly different. Random forest analysis revealed 24 bioindicator taxa that differentiated OZ, OLZ, and OMZ water communities. These bioindicator taxa included Burkholderiaceae, HOC36, SAR11 Clade IV, Thioglobaceae, Nitrospinaceae, SAR86, and UBA10353. Further, co-occurrence network analysis revealed that SAR202, AEGEAN-169, UBA10353, SAR406, and Rhodobacteraceae were keystone taxa among the entire interaction network of the microbial communities. Functional prediction further indicated that the relative abundances of microbial populations involved in nitrogen and sulfur cycling were higher in OMZs. Several microbial taxa, including the Thioglobaceae, Nitrospinaceae, SAR202, SAR406, WPS-2, UBA10353, and Woeseiaceae, may be involved in nitrogen and/or sulfur cycling, while also contributing to oxygen consumption in these waters. This study consequently provides new insights into the microbial community structures and potentially important taxa that contribute to oxygen consumption in the Andaman Sea and eastern BoB OMZ.

Substantial oxygen consumption by aerobic nitrite oxidation in oceanic oxygen minimum zones

Oceanic oxygen minimum zones (OMZs) are globally significant sites of biogeochemical cycling where microorganisms deplete dissolved oxygen (DO) to concentrations <20 µM. Amid intense competition for DO in these metabolically challenging environments, aerobic nitrite oxidation may consume significant amounts of DO and help maintain low DO concentrations, but this remains unquantified. Using parallel measurements of oxygen consumption rates and ¹⁵N-nitrite oxidation rates applied to both water column profiles and oxygen manipulation experiments, we show that the contribution of nitrite oxidation to overall DO consumption systematically increases as DO declines below 2 µM. Nitrite oxidation can account for all DO consumption only under DO concentrations <393 nM found in and below the secondary chlorophyll maximum. These patterns are consistent across sampling stations and experiments, reflecting coupling between nitrate reduction and nitrite-oxidizing Nitrospina with high oxygen affinity (based on isotopic and omic data). Collectively our results demonstrate that nitrite oxidation plays a pivotal role in the maintenance and biogeochemical dynamics of OMZs.

Dynamics and Distribution of Marine Synechococcus Abundance and Genotypes during Seasonal Hypoxia in a Coastal Marine Ranch

Marine Synechococcus are an ecologically important picocyanobacterial group widely distributed in various oceanic environments. Little is known about the dynamics and distribution of Synechococcus abundance and genotypes during seasonal hypoxia in coastal zones. In this study, an investigation was conducted in a coastal marine ranch along two transects in Muping, Yantai, where hypoxic events (defined here as the dissolved oxygen concentration <3 mg L−1) occurred in the summer of 2015. The hypoxia occurred in the bottom waters from late July and persisted until late August. It was confined at nearshore stations of the two transects, one running across a coastal ranch and the other one outside. During this survey, cell abundance of Synechococcus was determined with flow cytometry, showing great variations ranging from 1 × 104 to 3.0 × 105 cells mL−1, and a bloom of Synechococcus occurred when stratification disappeared and hypoxia faded out outside the ranch. Regression analysis indicated that dissolved oxygen, pH, and inorganic nutrients were the most important abiotic factors in explaining the variation in Synechococcus cell abundance. Diverse genotypes (mostly belonged to the sub-clusters 5.1 and 5.2) were detected using clone library sequencing and terminal restriction fragment length polymorphism analysis of the 16S–23S rRNA internal transcribed spacer region. The richness of genotypes was significantly related to salinity, temperature, silicate, and pH, but not dissolved oxygen. Two environmental factors, temperature and salinity, collectively explained 17% of the variation in Synechococcus genotype assemblage. With the changes in population composition in diverse genotypes, the Synechococcus assemblages survived in the coastal hypoxia event and thrived when hypoxia faded out.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

The cyanobacterium Prochlorococcus has divergent light-harvesting antennae and may have evolved in a low-oxygen ocean

The marine unicellular cyanobacterium Prochlorococcus is the most abundant photosynthetic organism on Earth. Members of this genus are classically thought to be adapted to high-oxygen and nutrient-poor ocean conditions, with a principle divergence between high-light and low-light ecotypes. We show that the most basal Prochlorococcus lineages are adapted to the low-oxygen, low-light, and high-nutrient conditions found in the dimly illuminated waters of anoxic marine zones. The most basal lineages have retained phycobilisomes as light-harvesting antennae—a characteristic of most other cyanobacteria—whose loss was thought to define all Prochlorococcus. As oxygenic photosynthesis drove ocean oxidation in the ancient Earth, oxygen appears to have played as much a role as light and nutrients in driving Prochlorococcus evolution.

Marine picocyanobacteria of the genus Prochlorococcus are the most abundant photosynthetic organisms in the modern ocean, where they exert a profound influence on elemental cycling and energy flow. The use of transmembrane chlorophyll complexes instead of phycobilisomes as light-harvesting antennae is considered a defining attribute of Prochlorococcus. Its ecology and evolution are understood in terms of light, temperature, and nutrients. Here, we report single-cell genomic information on previously uncharacterized phylogenetic lineages of this genus from nutrient-rich anoxic waters of the eastern tropical North and South Pacific Ocean. The most basal lineages exhibit optical and genotypic properties of phycobilisome-containing cyanobacteria, indicating that the characteristic light-harvesting antenna of the group is not an ancestral attribute. Additionally, we found that all the indigenous lineages analyzed encode genes for pigment biosynthesis under oxygen-limited conditions, a trait shared with other freshwater and coastal marine cyanobacteria. Our findings thus suggest that Prochlorococcus diverged from other cyanobacteria under low-oxygen conditions before transitioning from phycobilisomes to transmembrane chlorophyll complexes and may have contributed to the oxidation of the ancient ocean.

Genome-resolved viral ecology in a marine oxygen minimum zone

Oxygen minimum zones (OMZs) are critical to marine nitrogen cycling and global climate change. While OMZ microbial communities are relatively well-studied, little is known about their viruses. Here, we assess the viral community ecology of 22 deeply sequenced viral metagenomes along a gradient of oxygenated to anoxic waters (<0.02 μmol/l O₂ ) in the Eastern Tropical South Pacific (ETSP) OMZ. We identified 46 127 viral populations (≥5 kb), which augments the known viruses from ETSP by 10-fold. Viral communities clustered into six groups that correspond to oceanographic features. Oxygen concentration was the predominant environmental feature driving viral community structure. Alpha and beta diversity of viral communities in the anoxic zone were lower than in surface waters, which parallels the low microbial diversity seen in other studies. ETSP viruses were largely endemic, with the majority of shared viruses (87%) also present in other OMZ samples. We detected 543 putative viral-encoded auxiliary metabolic genes (AMGs), of which some have a distribution that reflects physico-chemical characteristics across depth. Together these findings provide an ecological baseline for viral community structure, drivers and population variability in OMZs that will help future studies assess the role of viruses in these climate-critical environments.

Spatial Distribution Patterns of Bacterioplankton in the Oxygen Minimum Zone of the Tropical Mexican Pacific

Microbial communities within oxygen minimum zones (OMZs) are crucial drivers of marine biogeochemical cycles; however, we still lack an understanding of how these communities are distributed across an OMZ. We explored vertical (from 5 to 500 m depth) and horizontal (coast to open ocean) distribution of bacterioplankton and its relationships with the main oceanographic conditions in three transects of the tropical Mexican Pacific OMZ. The distribution of the microbial diversity and the main clades changed along the transition from oxygen-rich surface water to the OMZ core, demonstrating the sensitivity of key bacterial groups to deoxygenation. The euphotic zone was dominated by Synechococcales, followed by Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales, whereas the OMZ core was dominated by SAR11, followed by SAR406, SAR324, SAR202, UBA10353 marine group, Thiomicrospirales and Nitrospinales. The marked environmental gradients along the water column also supported a high potential for niche partitioning among OMZ microorganisms. Additionally, in the OMZ core, bacterial assemblages from the same water mass were more similar to each other than those from another water mass. There were also important differences between coastal and open-ocean communities: Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales were more abundant in coastal areas, while Synechococcales, SAR406, SAR324, SAR202, UBA10353 marine group, and Thiomicrospirales were more abundant in the open ocean. Our results suggest a biogeographic structure of the bacterioplankton in this OMZ region, with limited community mixing across water masses, except in upwelling events, and little dispersion of the community by currents in the euphotic zone.

Cyanophage host-derived genes reflect contrasting selective pressures with depth in the oxic and anoxic water column of the Eastern Tropical North Pacific

Cyanophages encode host-derived genes that may increase their fitness. We examined the relative abundance of 18 host-derived cyanophages genes in metagenomes and viromes along depth profiles from the Eastern Tropical North Pacific Oxygen Deficient Zone (ETNP ODZ) where Prochlorococcus dominates a secondary chlorophyll maximum within the ODZ. Cyanophages at the oxic primary chlorophyll maximum encoded genes related to light and phosphate stress (psbA, psbD and pstS in T4-like and psbA in T7-like), but the proportion of cyanophage with these genes decreased with depth. The proportion of cyanophage with purine biosynthesis genes increased with depth in T4-like, but not T7-like cyanophages. No additional host-derived genes were found in deep T7-like cyanophages, suggesting that T4-like and T7-like cyanophages have different host-derived gene acquisition strategies, possibly linked to their different genome packaging mechanisms. In contrast to the ETNP, in the oxic North Atlantic T4-like cyanophages encoded psbA and pstS throughout the euphotic zone. Differences in pstS between the ETNP and the North Atlantic stations were consistent with differences in phosphate concentrations in those regimes. We suggest that the low proportion of cyanophage with psbA within the ODZ reflects the stably stratified low-light conditions occupied by their hosts, a Prochlorococcus ecotype endemic to ODZs.

Biogeochemistry and hydrography shape microbial community assembly and activity in the eastern tropical North Pacific Ocean oxygen minimum zone

Oceanic oxygen minimum zones (OMZs) play a pivotal role in biogeochemical cycles due to extensive microbial activity. How OMZ microbial communities assemble and respond to environmental variation is therefore essential to understanding OMZ functioning and ocean biogeochemistry. Sampling along depth profiles at five stations in the eastern tropical North Pacific Ocean (ETNP), we captured systematic variations in dissolved oxygen (DO) and associated variables (nitrite, chlorophyll, and ammonium) with depth and between stations. We quantitatively analysed relationships between oceanographic gradients and microbial community assembly and activity based on paired 16S rDNA and 16S rRNA sequencing. Overall microbial community composition and diversity were strongly related to regional variations in density, DO, and other variables (regression and redundancy analysis r²  = 0.68-0.82), displaying predictable patterns with depth and between stations. Although similar factors influenced the active community, diversity was substantially lower within the OMZ. We also identified multiple active microbiological networks that tracked specific gradients or features - particularly subsurface ammonium and nitrite maxima. Our findings indicate that overall microbial community assembly is consistently shaped by hydrography and biogeochemistry, while active segments of the community form discrete networks inhabiting distinct portions of the water column, and that both are tightly tuned to environmental conditions in the ETNP.

Microbial niches in marine oxygen minimum zones

In the ocean's major oxygen minimum zones (OMZs), oxygen is effectively absent from sea water and life is dominated by microorganisms that use chemicals other than oxygen for respiration. Recent studies that combine advanced genomic and chemical detection methods are delineating the different metabolic niches that microorganisms can occupy in OMZs. Understanding these niches, the microorganisms that inhabit them, and their influence on marine biogeochemical cycles is crucial as OMZs expand with increasing seawater temperatures.

Diversity and Spatial Distribution of Chromophytic Phytoplankton in the Bay of Bengal Revealed by RuBisCO Genes (rbcL)

Phytoplankton are the basis of primary production and play important roles in regulating energy export in marine ecosystems. Compared to other regions, chromophytic phytoplankton are considerably understudied in the Bay of Bengal (BOB). Here, we investigated community structure and spatial distribution of chromophytic phytoplankton in the BOB by using RuBisCO genes (Form ID rbcL). High throughput sequencing of rbcL genes revealed that diatoms, cyanobacteria (Cyanophyceae), Pelagophyceae, Haptophyceae, Chrysophyceae, Eustigamatophyceae, Xanthophyceae, Cryptophyceae, Dictyochophyceae, and Pinguiophyceae were the most abundant groups recovered in the BOB. Abundances and distribution of diatoms and Pelagophyceae were further verified using quantitative PCR analyses which showed the dominance of these groups near the Equator region (p < 0.01) where upwelling was likely the source of nutrients. Further, redundancy analysis demonstrated that temperature was an important environmental driver in structuring distributions of Cyanophyceae and dominant chromophytic phytoplankton. Morphological identification and quantification confirmed the dominance of diatoms, and also detected other cyanobacteria and dinoflagellates that were missing in our molecular characterizations. Pearson’s correlations of these morphologically identified phytoplankton with environmental gradients also indicated that nutrients and temperature were key variables shaping community structure. Combination of molecular characterization and morphological identification provided a comprehensive overview of chromophytic phytoplankton. This is the first molecular study of chromophytic phytoplankton accomplished in the BOB, and our results highlight a combination of molecular analysis targeting rbcL genes and microscopic detection in examining phytoplankton composition and diversity.

Cyanobacteria and cyanophage contributions to carbon and nitrogen cycling in an oligotrophic oxygen-deficient zone

Up to half of marine N losses occur in oxygen-deficient zones (ODZs). Organic matter flux from productive surface waters is considered a primary control on N₂ production. Here we investigate the offshore Eastern Tropical North Pacific (ETNP) where a secondary chlorophyll a maximum resides within the ODZ. Rates of primary production and carbon export from the mixed layer and productivity in the primary chlorophyll a maximum were consistent with oligotrophic waters. However, sediment trap carbon and nitrogen fluxes increased between 105 and 150 m, indicating organic matter production within the ODZ. Metagenomic and metaproteomic characterization indicated that the secondary chlorophyll a maximum was attributable to the cyanobacterium Prochlorococcus, and numerous photosynthesis and carbon fixation proteins were detected. The presence of chemoautotrophic ammonia-oxidizing archaea and the nitrite oxidizer Nitrospina and detection of nitrate oxidoreductase was consistent with cyanobacterial oxygen production within the ODZ. Cyanobacteria and cyanophage were also present on large (>30 μm) particles and in sediment trap material. Particle cyanophage-to-host ratio exceeded 50, suggesting that viruses help convert cyanobacteria into sinking organic matter. Nitrate reduction and anammox proteins were detected, congruent with previously reported N₂ production. We suggest that autochthonous organic matter production within the ODZ contributes to N₂ production in the offshore ETNP.

Strain specific differences in rates of Photosystem II repair in picocyanobacteria correlate to differences in FtsH protein levels and isoform expression patterns

Picocyanobacteria are the numerically dominant photoautotrophs of the oligotrophic regions of Earth’s oceans. These organisms are characterized by their small size and highly reduced genomes. Strains partition to different light intensity and nutrient level niches, with differing photosynthetic apparatus stoichiometry, light harvesting machinery and susceptibility to photoinactivation. In this study, we grew three strains of picocyanobacteria: the low light, high nutrient strain Prochlorococcus marinus MIT 9313; the high light, low nutrient Prochlorococcus marinus MED 4; and the high light, high nutrient marine Synechococcus strain WH 8102; under low and high growth light levels. We then performed matched photophysiology, protein and transcript analyses. The strains differ significantly in their rates of Photosystem II repair under high light and in their capacity to remove the PsbA protein as the first step in the Photosystem II repair process. Notably, all strains remove the PsbD subunit at the same rate that they remove PsbA. When grown under low light, MIT 9313 loses active Photosystem II quickly when shifted to high light, but has no measurable capacity to remove PsbA. MED 4 and WH 8102 show less rapid loss of Photosystem II and considerable capacity to remove PsbA. MIT 9313 has less of the FtsH protease thought to be responsible for the removal of PsbA in other cyanobacteria. Furthermore, by transcript analysis the predominant FtsH isoform expressed in MIT 9313 is homologous to the FtsH 4 isoform characterized in the model strain Synechocystis PCC 6803, rather than the FtsH 2 and 3 isoforms thought to be responsible for PsbA degradation. MED 4 on the other hand shows high light inducible expression of the isoforms homologous to FtsH 2 and 3, consistent with its faster rate of PsbA removal. MIT 9313 has adapted to its low light environment by diverting resources away from Photosystem II content and repair.

Comparison of photosynthetic performances of marine picocyanobacteria with different configurations of the oxygen-evolving complex

The extrinsic PsbU and PsbV proteins are known to play a critical role in stabilizing the Mn₄CaO₅ cluster of the PSII oxygen-evolving complex (OEC). However, most isolates of the marine cyanobacterium Prochlorococcus naturally miss these proteins, even though they have kept the main OEC protein, PsbO. A structural homology model of the PSII of such a natural deletion mutant strain (P. marinus MED4) did not reveal any obvious compensation mechanism for this lack. To assess the physiological consequences of this unusual OEC, we compared oxygen evolution between Prochlorococcus strains missing psbU and psbV (PCC 9511 and SS120) and two marine strains possessing these genes (Prochlorococcus sp. MIT9313 and Synechococcus sp. WH7803). While the low light-adapted strain SS120 exhibited the lowest maximal O₂ evolution rates (P_max per divinyl-chlorophyll a, per cell or per photosystem II) of all four strains, the high light-adapted strain PCC 9511 displayed even higher P^Chl_max and P^PSII_max at high irradiance than Synechococcus sp. WH7803. Furthermore, thermoluminescence glow curves did not show any alteration in the B-band shape or peak position that could be related to the lack of these extrinsic proteins. This suggests an efficient functional adaptation of the OEC in these natural deletion mutants, in which PsbO alone is seemingly sufficient to ensure proper oxygen evolution. Our study also showed that Prochlorococcus strains exhibit negative net O₂ evolution rates at the low irradiances encountered in minimum oxygen zones, possibly explaining the very low O₂ concentrations measured in these environments, where Prochlorococcus is the dominant oxyphototroph.

Genome Rearrangement Shapes Prochlorococcus Ecological Adaptation

Prochlorococcus, the most abundant and smallest known free-living photosynthetic microorganism, plays a key role in marine ecosystems and biogeochemical cycles. Prochlorococcus genome evolution is a fundamental issue related to how Prochlorococcus clades adapted to different ecological niches. Recent studies revealed that the gene gain and loss is crucial to the clade differentiation. The significance of our research is that we interpreted the Prochlorococcus genome evolution from the perspective of genome structure and associated the genome rearrangement with the Prochlorococcus clade differentiation and subsequent ecological adaptation.

Prochlorococcus is the most abundant and smallest known free-living photosynthetic microorganism and is a key player in marine ecosystems and biogeochemical cycles. Prochlorococcus can be broadly divided into high-light-adapted (HL) and low-light-adapted (LL) clades. In this study, we isolated two low-light-adapted clade I (LLI) strains from the western Pacific Ocean and obtained their genomic data. We reconstructed Prochlorococcus evolution based on genome rearrangement. Our results showed that genome rearrangement might have played an important role in Prochlorococcus evolution. We also found that the Prochlorococcus clades with streamlined genomes maintained relatively high synteny throughout most of their genomes, and several regions served as rearrangement hotspots. Backbone analysis showed that different clades shared a conserved backbone but also had clade-specific regions, and the genes in these regions were associated with ecological adaptations.

IMPORTANCEProchlorococcus, the most abundant and smallest known free-living photosynthetic microorganism, plays a key role in marine ecosystems and biogeochemical cycles. Prochlorococcus genome evolution is a fundamental issue related to how Prochlorococcus clades adapted to different ecological niches. Recent studies revealed that the gene gain and loss is crucial to the clade differentiation. The significance of our research is that we interpreted the Prochlorococcus genome evolution from the perspective of genome structure and associated the genome rearrangement with the Prochlorococcus clade differentiation and subsequent ecological adaptation.

Cryptic oxygen cycling in anoxic marine zones

Oxygen availability drives changes in microbial diversity and biogeochemical cycling between the aerobic surface layer and the anaerobic core in nitrite-rich anoxic marine zones (AMZs), which constitute huge oxygen-depleted regions in the tropical oceans. The current paradigm is that primary production and nitrification within the oxic surface layer fuel anaerobic processes in the anoxic core of AMZs, where 30-50% of global marine nitrogen loss takes place. Here we demonstrate that oxygenic photosynthesis in the secondary chlorophyll maximum (SCM) releases significant amounts of O2 to the otherwise anoxic environment. The SCM, commonly found within AMZs, was dominated by the picocyanobacteria Prochlorococcus spp. Free O2 levels in this layer were, however, undetectable by conventional techniques, reflecting a tight coupling between O2 production and consumption by aerobic processes under apparent anoxic conditions. Transcriptomic analysis of the microbial community in the seemingly anoxic SCM revealed the enhanced expression of genes for aerobic processes, such as nitrite oxidation. The rates of gross O2 production and carbon fixation in the SCM were found to be similar to those reported for nitrite oxidation, as well as for anaerobic dissimilatory nitrate reduction and sulfate reduction, suggesting a significant effect of local oxygenic photosynthesis on Pacific AMZ biogeochemical cycling.

Global genetic capacity for mixotrophy in marine picocyanobacteria

The assimilation of organic nutrients by autotrophs, a form of mixotrophy, has been demonstrated in the globally abundant marine picocyanobacterial genera Prochlorococcus and Synechococcus. However, the range of compounds used and the distribution of organic compound uptake genes within picocyanobacteria are unknown. Here we analyze genomic and metagenomic data from around the world to determine the extent and distribution of mixotrophy in these phototrophs. Analysis of 49 Prochlorococcus and 18 Synechococcus isolate genomes reveals that all have the transporters necessary to take up amino acids, peptides and sugars. However, the number and type of transporters and associated catabolic genes differ between different phylogenetic groups, with low-light IV Prochlorococcus, and 5.1B, 5.2 and 5.3 Synechococcus strains having the largest number. Metagenomic data from 68 stations from the Tara Oceans expedition indicate that the genetic potential for mixotrophy in picocyanobacteria is globally distributed and differs between clades. Phylogenetic analyses indicate gradual organic nutrient transporter gene loss from the low-light IV to the high-light II Prochlorococcus. The phylogenetic differences in genetic capacity for mixotrophy, combined with the ubiquity of picocyanobacterial organic compound uptake genes suggests that mixotrophy has a more central role in picocyanobacterial ecology than was previously thought.

Response of Prochlorococcus to varying CO2:O2 ratios

Carbon fixation has a central role in determining cellular redox poise, increasingly understood to be a key parameter in cyanobacterial physiology. In the cyanobacterium Prochlorococcus-the most abundant phototroph in the oligotrophic oceans-the carbon-concentrating mechanism is reduced to the bare essentials. Given the ability of Prochlorococcus populations to grow under a wide range of oxygen concentrations in the ocean, we wondered how carbon and oxygen physiology intersect in this minimal phototroph. Thus, we examined how CO2:O2 gas balance influenced growth and chlorophyll fluorescence in Prochlorococcus strain MED4. Under O2 limitation, per-cell chlorophyll fluorescence fell at all CO2 levels, but still permitted substantial growth at moderate and high CO2. Under CO2 limitation, we observed little growth at any O2 level, although per-cell chlorophyll fluorescence fell less sharply when O2 was available. We explored this pattern further by monitoring genome-wide transcription in cells shocked with acute limitation of CO2, O2 or both. O2 limitation produced much smaller transcriptional changes than the broad suppression seen under CO2 limitation and CO2/O2 co-limitation. Strikingly, both CO2 limitation conditions initially evoked a transcriptional response that resembled the pattern previously seen in high-light stress, but at later timepoints we observed O2-dependent recovery of photosynthesis-related transcripts. These results suggest that oxygen has a protective role in Prochlorococcus when carbon fixation is not a sufficient sink for light energy.

Aerobic Microbial Respiration In Oceanic Oxygen Minimum Zones

Oxygen minimum zones are major sites of fixed nitrogen loss in the ocean. Recent studies have highlighted the importance of anaerobic ammonium oxidation, anammox, in pelagic nitrogen removal. Sources of ammonium for the anammox reaction, however, remain controversial, as heterotrophic denitrification and alternative anaerobic pathways of organic matter remineralization cannot account for the ammonium requirements of reported anammox rates. Here, we explore the significance of microaerobic respiration as a source of ammonium during organic matter degradation in the oxygen-deficient waters off Namibia and Peru. Experiments with additions of double-labelled oxygen revealed high aerobic activity in the upper OMZs, likely controlled by surface organic matter export. Consistently observed oxygen consumption in samples retrieved throughout the lower OMZs hints at efficient exploitation of vertically and laterally advected, oxygenated waters in this zone by aerobic microorganisms. In accordance, metagenomic and metatranscriptomic analyses identified genes encoding for aerobic terminal oxidases and demonstrated their expression by diverse microbial communities, even in virtually anoxic waters. Our results suggest that microaerobic respiration is a major mode of organic matter remineralization and source of ammonium (~45-100%) in the upper oxygen minimum zones, and reconcile hitherto observed mismatches between ammonium producing and consuming processes therein.

Physiology and evolution of nitrate acquisition in Prochlorococcus

Prochlorococcus is the numerically dominant phototroph in the oligotrophic subtropical ocean and carries out a significant fraction of marine primary productivity. Although field studies have provided evidence for nitrate uptake by Prochlorococcus, little is known about this trait because axenic cultures capable of growth on nitrate have not been available. Additionally, all previously sequenced genomes lacked the genes necessary for nitrate assimilation. Here we introduce three Prochlorococcus strains capable of growth on nitrate and analyze their physiology and genome architecture. We show that the growth of high-light (HL) adapted strains on nitrate is ∼17% slower than their growth on ammonium. By analyzing 41 Prochlorococcus genomes, we find that genes for nitrate assimilation have been gained multiple times during the evolution of this group, and can be found in at least three lineages. In low-light adapted strains, nitrate assimilation genes are located in the same genomic context as in marine Synechococcus. These genes are located elsewhere in HL adapted strains and may often exist as a stable genetic acquisition as suggested by the striking degree of similarity in the order, phylogeny and location of these genes in one HL adapted strain and a consensus assembly of environmental Prochlorococcus metagenome sequences. In another HL adapted strain, nitrate utilization genes may have been independently acquired as indicated by adjacent phage mobility elements; these genes are also duplicated with each copy detected in separate genomic islands. These results provide direct evidence for nitrate utilization by Prochlorococcus and illuminate the complex evolutionary history of this trait.

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

The cyanobacteria Prochlorococcus and Synechococcus are important marine primary producers. We explored their distributions and covariance along a physico-chemical gradient from coastal to open ocean waters in the Northeastern Pacific Ocean. An inter-annual pattern was delineated in the dynamic transition zone where upwelled and eastern boundary current waters mix, and two new Synechococcus clades, Eastern Pacific Clade (EPC) 1 and EPC2, were identified. By applying state-of-the-art phylogenetic analysis tools to bar-coded 16S amplicon datasets, we observed higher abundance of Prochlorococcus high-light I (HLI) and low-light I (LLI) in years when more oligotrophic water intruded farther inshore, while under stronger upwelling Synechococcus I and IV dominated. However, contributions of some cyanobacterial clades were proportionally relatively constant, e.g. Synechococcus EPC2. In addition to supporting observations that Prochlorococcus LLI thrive at higher irradiances than other LL taxa, the results suggest LLI tolerate lower temperatures than previously reported. The phylogenetic precision of our 16S rRNA gene analytical approach and depth of bar-coded sequencing also facilitated detection of clades at low abundance in unexpected places. These include Prochlorococcus at the coast and Cyanobium-related sequences offshore, although it remains unclear whether these came from resident or potentially advected cells. Our study enhances understanding of cyanobacterial distributions in an ecologically important eastern boundary system.

Genomic potential for nitrogen assimilation in uncultivated members of Prochlorococcus from an anoxic marine zone

Cyanobacteria of the genus Prochlorococcus are the most abundant photosynthetic marine organisms and key factors in the global carbon cycle. The understanding of their distribution and ecological importance in oligotrophic tropical and subtropical waters, and their differentiation into distinct ecotypes, is based on genetic and physiological information from several isolates. Currently, all available Prochlorococcus genomes show their incapacity for nitrate utilization. However, environmental sequence data suggest that some uncultivated lineages may have acquired this capacity. Here we report that uncultivated low-light-adapted Prochlorococcus from the nutrient-rich, low-light, anoxic marine zone (AMZ) of the eastern tropical South Pacific have the genetic potential for nitrate uptake and assimilation. All genes involved in this trait were found syntenic with those present in marine Synechococcus. Genomic and phylogenetic analyses also suggest that these genes have not been aquired recently, but perhaps were retained from a common ancestor, highlighting the basal characteristics of the AMZ lineages within Prochlorococcus.

Deoxygenation alters bacterial diversity and community composition in the ocean's largest oxygen minimum zone

Oceanic oxygen minimum zones (OMZs) have a central role in biogeochemical cycles and are expanding as a consequence of climate change, yet how deoxygenation will affect the microbial communities that control these cycles is unclear. Here we sample across dissolved oxygen gradients in the oceans' largest OMZ and show that bacterial richness displays a unimodal pattern with decreasing dissolved oxygen, reaching maximum values on the edge of the OMZ and decreasing within it. Rare groups on the OMZ margin are abundant at lower dissolved oxygen concentrations, including sulphur-cycling Chromatiales, for which 16S rRNA was amplified from extracted RNA. Microbial species distribution models accurately replicate community patterns based on multivariate environmental data, demonstrate likely changes in distributions and diversity in the eastern tropical North Pacific Ocean, and highlight the sensitivity of key bacterial groups to deoxygenation. Through these mechanisms, OMZ expansion may alter microbial composition, competition, diversity and function, all of which have implications for biogeochemical cycling in OMZs.

Distribution and diversity of Prochlorococcus ecotypes in the Red Sea

Photosynthetic prokaryotes of the genus Prochlorococcus play a major role in global primary production in the world's oligotrophic oceans. A recent study on pelagic bacterioplankton communities in the northern and central Red Sea indicated that the predominant cyanobacterial 16S rRNA gene sequence types were from Prochlorococcus cells belonging to a high-light-adapted ecotype (HL II). In this study, we analyzed microdiversity of Prochlorococcus sp. at multiple depths within and below the euphotic zone in the northern, central, and southern regions of the Red Sea, as well as in surface waters in the same locations, but in a different season. Prochlorococcus dominated the communities in clone libraries of the amplified 16S-23S rRNA internal transcribed spacer (ITS) region. Almost no differences were found between samples from coastal or open-water sites, but a high diversity of Prochlorococcus ecotypes was detected at 100-meter depth in the water column. In addition, an unusual dominance of HL II-related sequences was observed in deeper waters. Our results indicate that the Red Sea harbors diverse Prochlorococcus lineages, but no novel ecotypes, despite its unusual physicochemical properties.

Application of pyrosequencing method for investigating the diversity of synechococcus subcluster 5.1 in open ocean

Synechococcus are distributed throughout the world's oceans and are composed of diverse genetic lineages. However, as they are much less abundant than Prochlorococcus in oligotrophic open oceans, their in-depth genetic diversity cannot be investigated using commonly used primers targeting both Prochlorococcus and Synechococcus. Thus, in this study, we designed a primer specific to the 16S-23S rRNA internal transcribed spacer (ITS) of the Synechococcus subcluster 5.1. Using the primer, we could selectively amplify Synechococcus sequences in oligotrophic seawater samples. Further, we showed that a barcoded amplicon pyrosequencing method could be applicable to investigate Synechococcus diversity using sequences retrieved in GenBank and obtained from environmental samples. Allowing sequence analyses of a large number of samples, this high-throughput method would be useful to study global biodiversity and biogeographic patterns of Synechococcus in marine environments.

Giant hydrogen sulfide plume in the oxygen minimum zone off Peru supports chemolithoautotrophy

In Eastern Boundary Upwelling Systems nutrient-rich waters are transported to the ocean surface, fuelling high photoautotrophic primary production. Subsequent heterotrophic decomposition of the produced biomass increases the oxygen-depletion at intermediate water depths, which can result in the formation of oxygen minimum zones (OMZ). OMZs can sporadically accumulate hydrogen sulfide (H₂S), which is toxic to most multicellular organisms and has been implicated in massive fish kills. During a cruise to the OMZ off Peru in January 2009 we found a sulfidic plume in continental shelf waters, covering an area >5500 km², which contained ∼2.2×10⁴ tons of H₂S. This was the first time that H₂S was measured in the Peruvian OMZ and with ∼440 km³ the largest plume ever reported for oceanic waters. We assessed the phylogenetic and functional diversity of the inhabiting microbial community by high-throughput sequencing of DNA and RNA, while its metabolic activity was determined with rate measurements of carbon fixation and nitrogen transformation processes. The waters were dominated by several distinct γ-, δ- and ε-proteobacterial taxa associated with either sulfur oxidation or sulfate reduction. Our results suggest that these chemolithoautotrophic bacteria utilized several oxidants (oxygen, nitrate, nitrite, nitric oxide and nitrous oxide) to detoxify the sulfidic waters well below the oxic surface. The chemolithoautotrophic activity at our sampling site led to high rates of dark carbon fixation. Assuming that these chemolithoautotrophic rates were maintained throughout the sulfidic waters, they could be representing as much as ∼30% of the photoautotrophic carbon fixation.

Postulated changes such as eutrophication and global warming, which lead to an expansion and intensification of OMZs, might also increase the frequency of sulfidic waters. We suggest that the chemolithoautotrophically fixed carbon may be involved in a negative feedback loop that could fuel further sulfate reduction and potentially stabilize the sulfidic OMZ waters.

Oceanographic and biological effects of shoaling of the oxygen minimum zone

Long-term declines in oxygen concentrations are evident throughout much of the ocean interior and are particularly acute in midwater oxygen minimum zones (OMZs). These regions are defined by extremely low oxygen concentrations (<20-45 μmol kg(-1)), cover wide expanses of the ocean, and are associated with productive oceanic and coastal regions. OMZs have expanded over the past 50 years, and this expansion is predicted to continue as the climate warms worldwide. Shoaling of the upper boundaries of the OMZs accompanies OMZ expansion, and decreased oxygen at shallower depths can affect all marine organisms through multiple direct and indirect mechanisms. Effects include altered microbial processes that produce and consume key nutrients and gases, changes in predator-prey dynamics, and shifts in the abundance and accessibility of commercially fished species. Although many species will be negatively affected by these effects, others may expand their range or exploit new niches. OMZ shoaling is thus likely to have major and far-reaching consequences.

Microbial oceanography of anoxic oxygen minimum zones

Vast expanses of oxygen-deficient and nitrite-rich water define the major oxygen minimum zones (OMZs) of the global ocean. They support diverse microbial communities that influence the nitrogen economy of the oceans, contributing to major losses of fixed nitrogen as dinitrogen (N(2)) and nitrous oxide (N(2)O) gases. Anaerobic microbial processes, including the two pathways of N(2) production, denitrification and anaerobic ammonium oxidation, are oxygen-sensitive, with some occurring only under strictly anoxic conditions. The detection limit of the usual method (Winkler titrations) for measuring dissolved oxygen in seawater, however, is much too high to distinguish low oxygen conditions from true anoxia. However, new analytical technologies are revealing vanishingly low oxygen concentrations in nitrite-rich OMZs, indicating that these OMZs are essentially anoxic marine zones (AMZs). Autonomous monitoring platforms also reveal previously unrecognized episodic intrusions of oxygen into the AMZ core, which could periodically support aerobic metabolisms in a typically anoxic environment. Although nitrogen cycling is considered to dominate the microbial ecology and biogeochemistry of AMZs, recent environmental genomics and geochemical studies show the presence of other relevant processes, particularly those associated with the sulfur and carbon cycles. AMZs correspond to an intermediate state between two "end points" represented by fully oxic systems and fully sulfidic systems. Modern and ancient AMZs and sulfidic basins are chemically and functionally related. Global change is affecting the magnitude of biogeochemical fluxes and ocean chemical inventories, leading to shifts in AMZ chemistry and biology that are likely to continue well into the future.

ProPortal: a resource for integrated systems biology of Prochlorococcus and its phage

ProPortal (http://proportal.mit.edu/) is a database containing genomic, metagenomic, transcriptomic and field data for the marine cyanobacterium Prochlorococcus. Our goal is to provide a source of cross-referenced data across multiple scales of biological organization—from the genome to the ecosystem—embracing the full diversity of ecotypic variation within this microbial taxon, its sister group, Synechococcus and phage that infect them. The site currently contains the genomes of 13 Prochlorococcus strains, 11 Synechococcus strains and 28 cyanophage strains that infect one or both groups. Cyanobacterial and cyanophage genes are clustered into orthologous groups that can be accessed by keyword search or through a genome browser. Users can also identify orthologous gene clusters shared by cyanobacterial and cyanophage genomes. Gene expression data for Prochlorococcus ecotypes MED4 and MIT9313 allow users to identify genes that are up or downregulated in response to environmental stressors. In addition, the transcriptome in synchronized cells grown on a 24-h light–dark cycle reveals the choreography of gene expression in cells in a ‘natural’ state. Metagenomic sequences from the Global Ocean Survey from Prochlorococcus, Synechococcus and phage genomes are archived so users can examine the differences between populations from diverse habitats. Finally, an example of cyanobacterial population data from the field is included.

Novel lineages of Prochlorococcus and Synechococcus in the global oceans

Picocyanobacteria represented by Prochlorococcus and Synechococcus have an important role in oceanic carbon fixation and nutrient cycling. In this study, we compared the community composition of picocyanobacteria from diverse marine ecosystems ranging from estuary to open oceans, tropical to polar oceans and surface to deep water, based on the sequences of 16S-23S rRNA internal transcribed spacer (ITS). A total of 1339 ITS sequences recovered from 20 samples unveiled diverse and several previously unknown clades of Prochlorococcus and Synechococcus. Six high-light (HL)-adapted Prochlorococcus clades were identified, among which clade HLVI had not been described previously. Prochlorococcus clades HLIII, HLIV and HLV, detected in the Equatorial Pacific samples, could be related to the HNLC clades recently found in the high-nutrient, low-chlorophyll (HNLC), iron-depleted tropical oceans. At least four novel Synechococcus clades (out of six clades in total) in subcluster 5.3 were found in subtropical open oceans and the South China Sea. A niche partitioning with depth was observed in the Synechococcus subcluster 5.3. Members of Synechococcus subcluster 5.2 were dominant in the high-latitude waters (northern Bering Sea and Chukchi Sea), suggesting a possible cold-adaptation of some marine Synechococcus in this subcluster. A distinct shift of the picocyanobacterial community was observed from the Bering Sea to the Chukchi Sea, which reflected the change of water temperature. Our study demonstrates that oceanic systems contain a large pool of diverse picocyanobacteria, and further suggest that new genotypes or ecotypes of picocyanobacteria will continue to emerge, as microbial consortia are explored with advanced sequencing technology.