Predicting responses of plankton to variations in essential nutrients is hampered by limited in situ measurements, a poor understanding of community composition, and the lack of reference gene catalogs for key taxa. Iron is a key driver of plankton dynamics and, therefore, of global biogeochemical cycles and climate. To assess the impact of iron availability on plankton communities, we explored the comprehensive bio‐oceanographic and bio‐omics data sets from Tara Oceans in the context of the iron products from two state‐of‐the‐art global scale biogeochemical models. We obtained novel information about adaptation and acclimation toward iron in a range of phytoplankton, including picocyanobacteria and diatoms, and identified whole subcommunities covarying with iron. Many of the observed global patterns were recapitulated in the Marquesas archipelago, where frequent plankton blooms are believed to be caused by natural iron fertilization, although they are not captured in large‐scale biogeochemical models. This work provides a proof of concept that integrative analyses, spanning from genes to ecosystems and viruses to zooplankton, can disentangle the complexity of plankton communities and can lead to more accurate formulations of resource bioavailability in biogeochemical models, thus improving our understanding of plankton resilience in a changing environment.

Nitrogen fixation has a critical role in marine primary production, yet our understanding of marine nitrogen-fixers (diazotrophs) is hindered by limited observations. Here, we report a quantitative image analysis pipeline combined with mapping of molecular markers for mining >2,000,000 images and >1300 metagenomes from surface, deep chlorophyll maximum and mesopelagic seawater samples across 6 size fractions (<0.2–2000 μm). We use this approach to characterise the diversity, abundance, biovolume and distribution of symbiotic, colony-forming and particle-associated diazotrophs at a global scale. We show that imaging and PCR-free molecular data are congruent. Sequence reads indicate diazotrophs are detected from the ultrasmall bacterioplankton (<0.2 μm) to mesoplankton (180–2000 μm) communities, while images predict numerous symbiotic and colony-forming diazotrophs (>20 µm). Using imaging and molecular data, we estimate that polyploidy can substantially affect gene abundances of symbiotic versus colony-forming diazotrophs. Our results support the canonical view that larger diazotrophs (>10 μm) dominate the tropical belts, while unicellular cyanobacterial and non-cyanobacterial diazotrophs are globally distributed in surface and mesopelagic layers. We describe co-occurring diazotrophic lineages of different lifestyles and identify high-density regions of diazotrophs in the global ocean. Overall, we provide an update of marine diazotroph biogeographical diversity and present a new bioimaging-bioinformatic workflow.