“…Here we considered all the exclusively meiofaunal phyla (Gnathostomulida, Kinorhyncha, Loricifera, Gastrotricha, and Tardigrada) and other metazoans that can be meiofaunal-sized during their life cycle (temporary meiofaunal taxa) ( Higgins & Thiel, 1988 ; Giere, 2009 ). After this filtering, the final dataset contained 10 phyla (Annelida, Cnidaria, Crustacea, Echinodermata, Gastrotricha, Mollusca, Nematoda, Nemertea, Platyhleminthes, and Rotifera) as previously implemented in other studies ( Brannock & Halanych, 2015 ; Bernardino et al, 2019 ; Fais et al, 2020 ; Bellisario et al, 2021 ; Castro et al, 2021 ; Coppo et al, 2023 ). Permutational analysis of variance (PERMANOVA; Anderson, Gorley & Clarke, 2008 ) was performed to compare environmental variables (rainfall, temperature, salinity, carbonate content, grain size, total organic matter, and biopolymeric composition), seascape coverage (the extent of an area that is encompassed within any of the MBON Marine Seascapes categories), and meiofaunal data (diversity metrics—Shannnon’s diversity index, phylogenetic diversity, and abundance of sequence reads) among sampled seasons (summer, autumn, winter, and spring) and sampled stations at Gramuté beach.…”