Cyanobacterial diazotrophs are considered to be the most important source of fixed N2 in the open ocean. Biological N2 fixation is catalyzed by the extremely O2-sensitive nitrogenase enzyme. In cyanobacteria without specialized N2-fixing cells (heterocysts), mechanisms such as decoupling photosynthesis from N2 fixation in space or time are involved in protecting nitrogenase from the intracellular O2 evolved by photosynthesis. However, it is not known how cyanobacterial cells limit O2 diffusion across their membranes to protect nitrogenase in ambient O2-saturated surface ocean waters. Here, we explored all known genomes of the major marine cyanobacterial lineages for the presence of hopanoid synthesis genes, since hopanoids are a class of lipids that might act as an O2 diffusion barrier. We found that, whereas all non-heterocyst-forming cyanobacterial diazotrophs had hopanoid synthesis genes, none of the marine Synechococcus, Prochlorococcus (non-N2-fixing), and marine heterocyst-forming (N2-fixing) cyanobacteria did. Finally, we conclude that hopanoid-enriched membranes are a conserved trait in non-heterocyst-forming cyanobacterial diazotrophs that might lower the permeability to extracellular O2 This membrane property coupled with high respiration rates to decrease intracellular O2 concentration may therefore explain how non-heterocyst-forming cyanobacterial diazotrophs can fix N2 in the fully oxic surface ocean.

Cyanobacterial diazotrophs are considered to be the most important source of fixed N₂ in the open ocean. Biological N₂ fixation is catalyzed by the extremely O₂-sensitive nitrogenase enzyme. In cyanobacteria without specialized N₂-fixing cells (heterocysts), mechanisms such as decoupling photosynthesis from N₂ fixation in space or time are involved in protecting nitrogenase from the intracellular O₂ evolved by photosynthesis. However, it is not known how cyanobacterial cells limit O₂ diffusion across their membranes to protect nitrogenase in ambient O₂-saturated surface ocean waters. Here, we explored all known genomes of the major marine cyanobacterial lineages for the presence of hopanoid synthesis genes, since hopanoids are a class of lipids that might act as an O₂ diffusion barrier. We found that, whereas all non−heterocyst-forming cyanobacterial diazotrophs had hopanoid synthesis genes, none of the marine Synechococcus, Prochlorococcus (non−N₂-fixing), and marine heterocyst-forming (N₂-fixing) cyanobacteria did. Finally, we conclude that hopanoid-enriched membranes are a conserved trait in non−heterocyst-forming cyanobacterial diazotrophs that might lower the permeability to extracellular O₂. This membrane property coupled with high respiration rates to decrease intracellular O₂ concentration may therefore explain how non−heterocyst-forming cyanobacterial diazotrophs can fix N₂ in the fully oxic surface ocean.

The multiple symbiotic partnerships between closely related species of the haptophyte algae Braarudosphaera bigelowii and the nitrogen-fixing cyanobacteria Candidatus Atelocyanobacterium thalassa (UCYN-A) contribute importantly to the nitrogen and carbon cycles in vast areas of the ocean. The diversity of the eukaryotic 18S rDNA phylogenetic gene marker has helped to identify some of these symbiotic haptophyte species, yet we still lack a genetic marker to assess its diversity at a finer scale. One of such genes is the ammonium transporter (amt) gene, which encodes the protein that might be involved in the uptake of ammonium from UCYN-A in these symbiotic haptophytes. Here, we designed three specific PCR primer sets targeting the amt gene of the haptophyte species (A1-Host) symbiotic with the open ocean UCYN-A1 sublineage, and tested them in samples collected from open ocean and near-shore environments. Regardless of the primer pair used at Station ALOHA, which is where UCYN-A1 is the pre-dominant UCYN-A sublineage, the most abundant amt amplicon sequence variant (ASV) was taxonomically classified as A1-Host. In addition, two out of the three PCR primer sets revealed the existence of closely-related divergent haptophyte amt ASVs (>95% nucleotide identity). These divergent amt ASVs had higher relative abundances than the haptophyte typically associated with UCYN-A1 in the Bering Sea, or co-occurred with the previously identified A1-Host in the Coral Sea, suggesting the presence of new diversity of closely-related A1-Hosts in polar and temperate waters. Therefore, our study reveals an overlooked diversity of haptophytes species with distinct biogeographic distributions partnering with UCYN-A, and provides new primers that will help to gain new knowledge of the UCYN-A/haptophyte symbiosis.

The unicellular N₂-fixing cyanobacteria UCYN-A live in symbiosis with haptophytes in the Braarudosphaera bigelowii lineage. Maintaining N₂-fixing symbioses between two unicellular partners requires tight coordination of multiple biological processes including cell growth and division and, in the case of the UCYN-A symbiosis, N₂ fixation of the symbiont and photosynthesis of the host. In this system, it is thought that the host photosynthesis supports the high energetic cost of N₂ fixation, and both processes occur during the light period. However, information on this coordination is very limited and difficult to obtain because the UCYN-A symbiosis has yet to be available in culture. Natural populations containing the UCYN-A2 symbiosis were manipulated to explore the effects of alterations of regular light and dark periods and inhibition of host photosynthesis on N₂ fixation (single cell N₂ fixation rates), nifH gene transcription, and UCYN-A2 cell division (fluorescent in situ hybridization and nifH gene abundances). The results showed that the light period is critical for maintenance of regular patterns of gene expression, N₂ fixation and symbiont replication and cell division. This study suggests a crucial role for the host as a producer of fixed carbon, rather than light itself, in the regulation and implementation of these cellular processes in UCYN-A.

The symbiotic unicellular cyanobacterium Candidatus Atelocyanobacterium thalassa (UCYN-A) is one of the most abundant and widespread nitrogen (N₂)-fixing cyanobacteria in the ocean. Although it remains uncultivated, multiple sublineages have been detected based on partial nitrogenase (nifH) gene sequences, including the four most commonly detected sublineages UCYN-A1, UCYN-A2, UCYN-A3 and UCYN-A4. However, very little is known about UCYN-A3 beyond the nifH sequences from nifH gene diversity surveys. In this study, single cell sorting, DNA sequencing, qPCR and CARD-FISH assays revealed discrepancies involving the identification of sublineages, which led to new information on the diversity of the UCYN-A symbiosis. 16S rRNA and nifH gene sequencing on single sorted cells allowed us to identify the 16S rRNA gene of the uncharacterized UCYN-A3 sublineage. We designed new CARD-FISH probes that allowed us to distinguish and observe UCYN-A2 in a coastal location (SIO Pier; San Diego) and UCYN-A3 in an open ocean location (Station ALOHA; Hawaii). Moreover, we reconstructed about 13% of the UCYN-A3 genome from Tara Oceans metagenomic data. Finally, our findings unveil the UCYN-A3 symbiosis in open ocean waters suggesting that the different UCYN-A sublineages are distributed along different size fractions of the plankton defined by the cell-size ranges of their prymnesiophyte hosts.

The unicellular cyanobacterium UCYN-A, one of the major contributors to nitrogen fixation in the open ocean, lives in symbiosis with single-celled phytoplankton. UCYN-A includes several closely related lineages whose partner fidelity, genome-wide expression and time of evolutionary divergence remain to be resolved. Here we detect and distinguish UCYN-A1 and UCYN-A2 lineages in symbiosis with two distinct prymnesiophyte partners in the South Atlantic Ocean. Both symbiotic systems are lineage specific and differ in the number of UCYN-A cells involved. Our analyses infer a streamlined genome expression towards nitrogen fixation in both UCYN-A lineages. Comparative genomics reveal a strong purifying selection in UCYN-A1 and UCYN-A2 with a diversification process ∼91 Myr ago, in the late Cretaceous, after the low-nutrient regime period occurred during the Jurassic. These findings suggest that UCYN-A diversified in a co-evolutionary process, wherein their prymnesiophyte partners acted as a barrier driving an allopatric speciation of extant UCYN-A lineages.

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H₂ oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Abstract. Marine diazotrophs convert dinitrogen (N2) gas into bioavailable nitrogen (N), supporting life in the global ocean. In 2012, the first version of the global oceanic diazotroph database (version 1) was published. Here, we present an updated version of the database (version 2), significantly increasing the number of in situ diazotrophic measurements from 13 565 to 55 286. Data points for N2 fixation rates, diazotrophic cell abundance, and nifH gene copy abundance have increased by 184 %, 86 %, and 809 %, respectively. Version 2 includes two new data sheets for the nifH gene copy abundance of non-cyanobacterial diazotrophs and cell-specific N2 fixation rates. The measurements of N2 fixation rates approximately follow a log-normal distribution in both version 1 and version 2. However, version 2 considerably extends both the left and right tails of the distribution. Consequently, when estimating global oceanic N2 fixation rates using the geometric means of different ocean basins, version 1 and version 2 yield similar rates (43–57 versus 45–63 Tg N yr−1; ranges based on one geometric standard error). In contrast, when using arithmetic means, version 2 suggests a significantly higher rate of 223±30 Tg N yr−1 (mean ± standard error; same hereafter) compared to version 1 (74±7 Tg N yr−1). Specifically, substantial rate increases are estimated for the South Pacific Ocean (88±23 versus 20±2 Tg N yr−1), primarily driven by measurements in the southwestern subtropics, and for the North Atlantic Ocean (40±9 versus 10±2 Tg N yr−1). Moreover, version 2 estimates the N2 fixation rate in the Indian Ocean to be 35±14 Tg N yr−1, which could not be estimated using version 1 due to limited data availability. Furthermore, a comparison of N2 fixation rates obtained through different measurement methods at the same months, locations, and depths reveals that the conventional 15N2 bubble method yields lower rates in 69 % cases compared to the new 15N2 dissolution method. This updated version of the database can facilitate future studies in marine ecology and biogeochemistry. The database is stored at the Figshare repository (https://doi.org/10.6084/m9.figshare.21677687; Shao et al., 2022).