Global maps of the Ocean microbiome: Presence probability (0-1) of 53 Operational Taxonomic Units of the phylum Picozoa (pOTUs) in the epipelagic layer (0-200 m), projected under contemporary (2005-2017) environmental conditions, using an ensemble of habitat suitability models.

This collection of global maps provides global projections of habitat suitability index (HSI) estimated for each of 53 different picozoa operational taxonomic units (pOTUs), as described in Huber et al. 2024. For each pOTU-level map, mean, max, min, stdev, and coefficient of variation of HSI values computed across several algorithms are provided. We developed the single pOTU HSI maps from 2,394 eDNA samples relative to the pico-size fraction (i.e. 0.2 to 5 μm), which were retrieved from the EukBank database. Picozoa samples were classified based on their depth into epipelagic (0-200 m), mesopelagic (200-1000 m) and bathypelagic (>1000 m) samples. To estimate the HSI, we used species distribution models based on a presence/absence design; among all retrieved samples of pico-eukaryotes, georeferenced samples where a specific pOTU was detected were treated as presences, while sampling locations where that pOTU was not detected were treated as (pseudo-) absences. For each pOTU, we contrasted average multi-year depth-specific environmental conditions at presence vs absence points using an ensemble of modelling algorithms, including Generalized Linear Models (GLM), Artificial Neural Network (ANN), Boosted Regression Trees (BRT) and Random Forest (RF). Models were calibrated using multi-year average conditions across epipelagic, mesopelagic and bathypelagic layers. HSI maps were projected over epipelagic conditions only, due to lack of sufficient data coverage for deeper layers. Environmental predictors were retrieved from the World Ocean Atlas 2018. 
This collection, and other collections of Global maps of the Ocean microbiome are described in the AtlantECO-MAPS publication (https://doi.org/10.5281/zenodo.18696236).