Biological nitrogen fixation is a major factor contributing to microbial primary productivity in the open ocean. The current view depicts a few cyanobacterial diazotrophs as the most relevant marine nitrogen fixers, whereas heterotrophic diazotrophs are more diverse and considered to have lower impacts on the nitrogen balance. Here, we used 891 Tara Oceans metagenomes to create a manually curated, non-redundant genomic database corresponding to free-living, as well as filamentous, colony-forming, particle-attached and symbiotic bacterial and archaeal populations occurring in the surface of five oceans and two seas. Notably, the database provided the genomic content of eight cyanobacterial diazotrophs including Trichodesmium populations and a newly discovered population similar to Richelia, as well as 40 heterotrophic bacterial diazotrophs organized into three main functional groups that considerably expand the known diversity of abundant marine nitrogen fixers compared to previous genomic surveys. Critically, these 48 populations may account for more than 90% of cells containing known nifH genes and occurring in the sunlit ocean, suggesting that the genomic characterization of the most abundant marine diazotrophs may be nearing completion. The newly identified heterotrophic bacterial diazotrophs are widespread, express their nifH genes in situ, and co-occur under nitrate-depleted conditions in large size fractions where they might form aggregates providing the low-oxygen microenvironments required for nitrogen fixation. Most significantly, we found heterotrophic bacterial diazotrophs to be more abundant than cyanobacterial diazotrophs in most metagenomes from the open oceans and seas. This large-scale environmental genomic survey emphasizes the considerable potential of heterotrophs in the marine nitrogen balance.

Biological nitrogen fixation sustains ~50% of ocean primary production. However, our understanding of marine N2-fixers (diazotrophs) is hindered by limited observations. Here, we developed a quantitative image analysis pipeline in concert with mapping of molecular markers for mining >2,000,000 images and >1,300 metagenomes from Tara Oceans, covering surface, deep chlorophyll maximum and mesopelagic layers across 6 organismal size fractions (0-2000 μm). Imaging and molecular data were remarkably congruent. Diazotrophs were detected from ultrasmall bacterioplankton (<0.2 μm) to mesoplankton (180 to 2000 μm). We identified several new high density regions of diazotrophs. Distributional and abundance patterns support the previous canonical view that larger sized diazotrophs (>10 μm) dominate the tropical belts, while unicellular diazotrophs were found in surface and mesopelagic samples. Multiple co-occurring diazotrophic lineages were frequently encountered, suggesting that complex overlapping niches are common. Overall, this work provides an updated global snapshot of marine diazotroph biogeographical diversity and highlights new sources and sinks of diazotroph-fueled new production.

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison withTaraOceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.