Environmental gradients and habitat specificity structure benthic microbial assemblages in a temperate seagrass ecosystem

Seagrass meadows are globally distributed coastal ecosystems that provide habitats for diverse species and facilitate ecosystem services, such as carbon storage and nutrient cycling. However, factors linked to human activities and climate change are causing seagrass decline worldwide, thus impacting benthic community structure and ecosystem functioning. Although macrofaunal communities in seagrass habitats are well studied, the structure and function of bacteria/archaea and microeukaryote assemblages remains underexplored. In this study, we used environmental DNA (eDNA) metabarcoding (16S and 18S rRNA genes) to analyze bacteria/archaea (16S), microeukaryotes (18S raw sediment), and meiofauna (18S Ludox) communities in seagrass and bare sediment habitats across three sites along an estuarine gradient in Bodega Harbor, California. We found that microbial and microeukaryote communities were primarily structured by site as a response to the environmental gradient (e.g., sediment properties), whereas meiofauna responded more strongly to the presence of seagrass. Beta diversity showed strong differentiation among estuarine sites, while alpha diversity varied within sites, with bare sediments often displaying higher diversity. Network analysis revealed significant co-occurrences between bacteria and nematodes, including chemosynthetic symbiont-bearing lineages in mudflats and stress-tolerant bioindicators. Our findings underscore the importance of environmental gradients and microhabitat features in shaping benthic biodiversity in temperate seagrass ecosystems. Furthermore, it highlights potential microbial-nematode associations, likely a response to environmental variation . By integrating multiple benthic components, our study advances understanding of biodiversity drivers in temperate estuarine seagrass habitats and provides new insights on microbial community structure in these ecosystems.