Fig. 3

PCO ordination plot of sampling stations based on the Bray–Curtis similarity calculated by taxa. Stations are labelled as “Season × Area” interaction levels (WiMar, Winter Margherita di Savoia; WiBar, Winter Barletta; WiLes, Winter Lesina; WiVar, Winter Varano; SuMar, Summer Margherita; SuBar, Summer Barletta; SuLes, Summer Lesina; SuVar, Summer Varano). Vectors indicate the correlation of FGs to the multidimensional space according to Pearson’s correlation coefficient (ρ > 0.4). Taxa are coded as: Acanthocardia tuberculata (Act), Acteon tornatilis (Acte), Astropecten spp. (Ass), Diogenes pugilator (Dip), Dosinia lupinus (Dol), Donax semistriatus (Dos), Echinocardium cordatum (Ecc), Glycimeris glycimeris (Glg), Glycera rouxii (Glr), Hyalinoecia tubicola (Hyt), Holothuria poli (Hop), Liocarcinus depurator (Lid), Maldanidae (Mald), Moerella donacina (Mod), Oestergrenia digitata (Oed), Owenia fusiformis (Owf), Peronidia albicans (Pea), Sipunculus nudus (Sin), Tritia mutabilis (Trm) Tritia reticulata (Trr)