Structural variability of protist assemblages in surface sediments across Italian Mediterranean marine subregions

Marine sediments host heterogeneous protist communities consisting of both living benthic microorganisms and planktonic resting stages. Despite their key functions in marine ecosystem processes and biogeochemical cycles, their structure and dynamics are largely unknown. In the present study, with a spatially intensive sampling design we investigated benthic protist diversity and function of surface sediment samples from three subregions of the Mediterranean Sea, through an environmental DNA metabarcoding approach targeting the 18S V4 region of rRNA gene. Protists were characterized at the taxonomic level and trophic function, both in terms of alpha diversity and community composition, testing for potential differences among marine subregions and bathymetric groups. Overall, Alveolata and Stramenopiles were the two divisions that dominated the communities. These dominant groups exhibited significant differences among the three Mediterranean subregions in the alpha diversity estimates based on the detected ASVs, for all computed indices (ASV richness, Shannon and Simpson indices). Protist communities were also found to be significantly different in terms of composition at the order rank in the three subregions p-value < 0.01). These differences were mainly driven by Anoecales, Peridiniales, Borokales, Paraliales and Gonyaulacales, which together contributed almost 80% of the average dissimilarity. Anoecales was the dominant order in the Ionian – Central Mediterranean and Adriatic Sea, but with considerably different relative abundances (52% and 36%, respectively), while Borokales was the dominant order in the Western Mediterranean Sea (33%). Similarly, significant differences among the three marine subregions were also highlighted when protist assemblages were examined in terms of trophic function, both in terms of alpha diversity (calculated on the ASVs for each trophic group) and community composition p-value < 0.01. In particular, the Adriatic Sea stood out for having the highest relative abundance of autotrophic/mixotrophic components in the surface sediments analyzed. Conversely, no significant differences in protist assemblages were found among depth groups. This study provided new insights into the taxonomic and trophic composition of benthic protist communities found in Mediterranean surface sediments, revealing geographical differences among regional seas. The results were discussed in relation to the Mediterranean environmental features that could generate the differences among benthic protist communities.