Next Article in Journal
Poor Correlation between Diamondback Terrapin Population Estimates Using Two New Estimation Methods
Next Article in Special Issue
Perception of Changes in Marine Benthic Habitats: The Relevance of Taxonomic and Ecological Memory
Previous Article in Journal
Phylogenetics of Molecular Regulators Contributing to Plant Stress Tolerance
Previous Article in Special Issue
Filling the Gap and Improving Conservation: How IUCN Red Lists and Historical Scientific Data Can Shed More Light on Threatened Sharks in the Italian Seas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity and Distribution Patterns of Hard Bottom Polychaete Assemblages in the North Adriatic Sea (Mediterranean)

1
Interdepartmental Research Centre for Environmental Sciences, University of Bologna, Via Sant’Alberto 163, 48123 Ravenna, Italy
2
Department of Biological and Environmental Sciences and Technologies, University of Salento, Via Provinciale Lecce-Monteroni, 73100 Lecce, Italy
3
Center for Marine Research, Ruđer Bošković Institute, Giordano Paliaga 5, 52210 Rovinj, Croatia
4
Stazione Zoologica Anton Dohrn, Integrative Marine Ecology Department, Villa Comunale, 80121 Naples, Italy
*
Author to whom correspondence should be addressed.
Submission received: 3 September 2020 / Revised: 11 October 2020 / Accepted: 19 October 2020 / Published: 21 October 2020

Abstract

:
The knowledge on the hard bottom polychaete assemblages in the Northern Adriatic Sea, a Mediterranean region strongly affected by environmental pressures, is scarce and outdated. The objective of this paper was to update the information on polychaete diversity and depict their patterns of natural spatial variation, in relation to changes in algal coverage at increasing depth. Hard bottom benthos was quantitatively sampled by scraping off the substrate from three stations at Sveti Ivan Island (North Adriatic) at three depths (1.5 m, 5 m and 25 m). Polychaete fauna comprised 107 taxa (the majority of them identified at species level) belonging to 22 families, with the family Syllidae ranking first in terms of number of species, followed by Sabellidae, Nereididae, Eunicidae and Serpulidae. Considering the number of polychaete species and their identity, the present data differed considerably from previous studies carried out in the area. Two alien species, Lepidonotus tenuisetosus, which represented a new record for the Adriatic Sea, and Nereis persica, were recorded. The highest mean abundance, species diversity and internal structural similarity of polychaete assemblages were found at 5 m depth, characterised by complex and heterogeneous algal habitat. The DISTLM forward analysis revealed that the distribution of several algal taxa as well as some algal functional-morphological groups significantly explained the observed distribution patterns of abundance and diversity of polychaete assemblages. The diversity of the North Adriatic hard bottom polychaete fauna is largely underestimated and needs regular updating in order to detect and monitor changes of benthic communities in the area.

Graphical Abstract

1. Introduction

The North Adriatic Sea is the northernmost sector of the Mediterranean Sea, with peculiar geomorphological, hydrographical and biogeographical characteristics. With an average depth of 35 m, this semienclosed basin represents the most extensive region of shallow water in the Mediterranean [1,2], being one of the most productive areas too, due to high amount of nutrients loaded by the Po River [3]. It is also the coldest Mediterranean sector, together with the Gulf of Lion and the North Aegean Sea, thus inhabited by species of boreal affinity [4,5]. The Adriatic Sea and in particular its northern part, exhibits the highest species richness of invertebrates in the Mediterranean basin [6]. Being also densely populated and thus under high anthropogenic pressures, the North Adriatic is a sensitive area currently undergoing severe environmental changes (climate change, fishing impacts, destruction of habitats, pollution, introduction of non indigenous species) that affect the benthic communities [6,7,8]. Updating the knowledge about benthic diversity and understanding patterns of benthic assemblages’ vertical and horizontal spatial variation are benchmarks for detecting and monitoring environmental changes, also according to the European Union Marine Strategy Framework Directive [9].
Polychaetes are among the most abundant and species-rich marine benthic groups, showing a wide functional diversity and adaptation to different environmental conditions [10,11]. Thus, they are often used as surrogates to estimate the state and dynamics of benthic communities [12,13,14,15]. North Adriatic soft bottom polychaete fauna is well known, e.g., [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], while the knowledge on hard bottom polychaete assemblages from natural substrates is scarce, i.e., [33,34,35,36,37,38,39,40,41,42] and their diversity might be largely underestimated as indicated by recent studies dealing with Syllidae and Sabellidae polychaete families [43,44].
The structure of hard bottom benthic assemblages is characterised by having high small- and middle-scale spatial variability (i.e., patchiness), both alongshore and at different depths, which is caused by various interplaying biological (e.g., predation, competition, recruitment) and physical-chemical (e.g., light intensity, temperature, salinity, hydrodynamics, sedimentation, habitat complexity) factors [45,46,47,48,49,50]. In particular, the distribution of hard-bottom polychaetes is strongly dependent on the bathymetric variation in algal composition and the associated changes in algal forms [51,52,53]. However, studies aiming at understanding the role of the above-mentioned factors in structuring polychaete assemblages have not been done in the North Adriatic so far.
The aims of our study were: (1) to update the knowledge of the faunal composition of the North Adriatic hard bottom polychaete assemblages, and (2) to assess variation of their spatial distribution in relation to changes in algal assemblages along a bathymetric gradient and according to substrate orientation.

2. Materials and Methods

The study area was in the vicinity of the city of Rovinj (Croatia, North Adriatic Sea) at Sveti Ivan Island. Islands of the Rovinj archipelago and the coastal area of up to 500 m from the coastline, were proclaimed by Rovinj Municipality a natural landscape reserve. The area is characterized by calcareous rocky shelf extending from 0 to about 25 m depth, with a gentle-medium slope. The submarine slopes of the Sveti Ivan Island are representing typical infralittoral environments of the North Adriatic Shelf.
Collecting surveys were carried out in June 2007, taking into consideration that the maximum development of macroalgal assemblages in the Northern Adriatic Sea occurs in the spring–early summer period [54,55]. Benthos was sampled using scuba diving at stations A (N 45° 02.69′, E 13° 37.18′) and B (N 45° 02.7′, E 13° 37.48′) on the southern side and station C (N 45° 02.87′, E 13° 37.34′) on the northern side of the Island (Figure 1). At each station, three depths (1.5 m, 5 m and 25 m) were appointed along a vertical transect and at each depth three replicates of 10 × 10 cm surface quadrats covered with macroalgae were randomly chosen (27 samples in total). Samples were collected by scraping off the substrate including the whole algal coverage present within the 10 × 10 cm quadrats using hammer and chisel. The scraped material was collected within plastic bags. Although ordinarily a sampling area of 20 × 20 cm is suggested for benthic studies in the Mediterranean [56], the small sample size (10 × 10 cm) in this research was chosen in order to minimize sampling impacts in the natural reserve. Furthermore, 10 × 10 cm replicate areas were already used in other studies dealing with polychaetes in the Mediterranean Sea and revealed to allow acceptable representation of polychaete diversity and distribution patterns [11,57,58]. Each replicate unit (10 × 10 cm surface quadrats) was photographed underwater to facilitate the determination of the associated algal taxa and their percent coverage.
In the laboratory, samples were fixed in 8% formaldehyde seawater solution and subsequently rinsed with fresh water and sieved through 0.5mm mesh and preserved in 70% ethanol. Polychaetes and algae were sorted and determined to the lowest taxonomic level possible using stereo- and light microscope. In some cases, it was not possible to identify the organisms at species level but only at higher taxonomic categories (e.g., genus, family). For this reason, we used the term taxa instead of species throughout the manuscript to indicate the recorded taxonomic entities. When specimens belonging to the same genus were clearly different from each other, however, we considered them belonging to different undetermined species (e.g., Nereis sp. 1 and Nereis sp. 2). The coverage of each algal thallus, representing the surface covered in an orthogonal projection, was determined according to Boudouresque [56] and Cormaci et al. [59]. Algae from each sample were placed on the surface that has an area equal to that sampled in situ (10 × 10 cm). Consequently, for each algal species the percentage of total quadrat area (100 cm2) covered by the projection of all the thalli was estimated. Value for total algal coverage can reach more than 100% in the presence of multilayered assemblages (as canopy forming algae) or epibiosis [56]. Additionally, algae were grouped in the following functional-morphological groups: Articulated calcareous, Corticated, Encrusting, Filamentous, Foliose and Leathery [60,61], and coverage of each functional-morfological group per sample was calculated as well. The polychaete specimens were deposited at the Center for Marine Research of the Ruđer Bošković Institute in Rovinj (Croatia) (IRB-CIM) and in the collection of the Natural History Museum in Rijeka (Croatia) (PMR).
For each replicate sample, polychaete assemblages were characterised by their respective abundance (N), species richness (S), Hill’s species diversity index (N1) and Hill’s evenness index (N10) [62]. In order to graphically represent trends in the number of species found within the collected samples, species accumulation curve of observed species (Sobs) was created. Moreover, to estimate the number of species potentially present in the area, curves of estimated number of species were calculated using the Jacknife 1, Jacknife 2 and Bootstrap methods [63].
Nonparametric distance-based permutational analysis of variance by permutation of residuals under a reduced model (PERMANOVA) [64,65] was used to test for differences in univariate indices (based on Euclidean distances of untransformed data) and in multivariate structure (based on Bray–Curtis similarity of untransformed data) of the polychaete assemblages between stations and depths. PERMANOVA design included two crossed factors: station (3 levels, random) and depth (3 levels, fixed). Post-hoc pair-wise comparisons allowed detecting the source of significant variations. For significant terms, a permutational analysis of multivariate dispersions (PERMDISP) [66] was used to test the homogeneity of samples dispersion from their group centroids. When the number of permutations was low (less than 1000), Monte Carlo probability (P(MC)) was considered instead of permutational probability (P(perm)). To calculate p values for PERMANOVA and PERMDISP, 9999 permutations were used.
Multivariate patterns were visualised by nonmetric multidimensional scaling ordination (nMDS). The similarity percentage routine (SIMPER) [67] (70% cut off), was used to detect the taxa most responsible for within-depth similarity and between-depth dissimilarity and, at each depth, the between-station dissimilarity. When analysing dissimilarity between stations, taxa were considered important if they exceeded an arbitrarily chosen threshold of 4% of dissimilarity between stations at each depth.
Potential relationships among structuring algal taxa, algal functional-morphological groups, depth and orientation of the sampling station in respect to the island geography (south/north)) and the variation of polychaete assemblages were assessed by nonparametric multiple regression analyses based on Bray–Curtis dissimilarities [65] exerting the distance-based multivariate analysis for a linear model using forward selection procedure (DISTLM), with 9999 permutations. Resemblance matrix produced by DISTLM analyses informed on the correlation among all pairs of explanatory variables to check for multicolinearity [68]. The predictor variables included depth, orientation and percent coverage of each structuring algal species (cut-off 5% of the total cover) in the first DISTLM analysis and depth, orientation and percent coverage of algal functional-morphological groups in the second one. Results of the forward selection procedure with the sequential tests (i.e., fitting each variable one at a time, conditional on the variables that are already included in the model) were presented. All analyses were done using PRIMER v.6 [69], with the add-on PERMANOVA+ [66].

3. Results

3.1. Algal Assemblages

Forty-eight algal taxa were recorded in the research area (Table S1). Total mean algal coverage was higher at 1.5 m and 5 m depths if compared to 25 m depth in all three stations (Figure 2). At 1.5 m depth, algal assemblages were characterised by the high coverage of the articulated calcareous algae (mostly Corallina officinalis Linnaeus, Haliptilon sp. and Jania spp.), the corticated algae (mostly Alsidium sp. and Laurencia obtusa (Hudson) J.V.Lamouroux) and the foliose algae (mostly Padina pavonica (Linnaeus) Thivy, Dictyota dichotoma (Hudson) J.V.Lamouroux and Dictyota dichotoma var. intricata (C.Agardh) Greville) (Figure 3a). Filamentous algae (mostly Ectocarpus sp. and Carradoriella elongata (Hudson) A.M.Savoie & G.W.Saunders) and encrusting algae (mostly Valonia utricularis (Roth) C.Agardh and Peyssonnelia rubra (Greville) J.Agardh) were present with low percent coverage at 1.5 m depth. At 5 m depth, encrusting algae (mostly P. rubra and P. heteromorpha (Zanardini) Athanasiadis) dominated, followed by foliose algae (mostly P. pavonica, Flabellia petiolata (Turra) Nizamuddin, 1987 and D. dichotoma), while articulated calcareous algae (mostly C. officinalis, Haliption sp. and Jania spp) were present with a lower coverage (Figure 3b). Corticated algae (mostly Gelidium spinosum (S.G.Gmelin) P.C.Silva) and filamentous algae (mostly Polysiphonia sp. and Cladophora spp.) had very low percent coverage at 5 m depth. At 25 m depth, filamentous algae (mostly Cladophora spp., Sphacelaria plumula Zanardini, and Polysiphonia sp.) were prevalent and followed by encrusting algae (mostly P. rubra and an unidentified rose coloured encrusting algae) (Figure 3c). Foliose algae (mostly Rhodymenia sp., D. dichotoma and unidentified foliose algae) showed a very low coverage at 25 m depth. At this depth, corticated algae (mostly Rodriguezella sp.) were present only at stations B and C, while articulated calcareous algae (only Halimeda tuna (J.Ellis & Solander) J.V.Lamouroux) were found only at station C, all of them with very low coverage. Cystoseira compressa (Esper) Gerloff & Nizamuddin, the single species of the leathery functional group, appeared with a 2% coverage only in one sample from 1.5 m depth (Figure 2).

3.2. Composition and Diversity of Polychaete Assemblages

As a whole, 1993 polychaete specimens from 107 taxa (94 at species level) belonging to 22 families were found (Table S2). The richest families in terms of number of species were Syllidae (39 species), Sabellidae (13), Nereididae (12), Eunicidae (7) and Serpulidae (7), while eleven families were represented by only one species. Altogether 62 species from 16 families were found at 1.5 m depth, 73 species from 15 families at 5 m depth and 66 species from 15 families at 25 m depth.
The mean number of individuals and species, as well as the mean Hill’s species diversity index, were the highest at 5 m depth, while no particular pattern could be observed in the Hill’s evenness index (Figure 4a,b). Significant differences in species richness and Hill’s species diversity index were revealed both horizontally (between stations) and vertically (between depths), while differences in polychaete abundance and Hill’s evenness index were significant only between stations (Table 1). In particular, the assemblages at 5 m depth had significantly higher species richness and Hill’s species diversity index than those of 1.5 m and 25 m depths, at most of the stations (Table S3). Significant differences between stations in all univariate diversity descriptors were revealed mostly at 1.5m depth and to a lesser extent also at 25 m depth (Table S3).
Syllidae was the dominant family at all three depths, both in abundance and species richness, followed by Nereididae, Sabellidae and Eunicidae at 1.5 m depth and, at 5m depth, by Sabellidae, Nereididae and Eunicidae considering abundance, and Nereididae, Eunicidae and Sabellidae considering species richness (Figure 5). At 25 m depth, Syllidae were particularly dominant, and followed by Sabellidae, Nereididae, Serpulidae and Eunicidae, both considering abundance and number of species. Abundance and diversity of Serpulidae increased with depth. In fact, at 1.5 m only one specimen of Vermiliopsis infundibulum (Philippi, 1844) was found, while the Serpulidae were represented by 18 specimens belonging to 4 species at 5 m depth and 23 specimens belonging to 6 species at 25 m depth.
Most species were locally rare, with 36 species found only in one sample (Figure S1), while only Sphaerosyllis pirifera Claparède, 1868, Amphiglena mediterranea (Leydig, 1851) and Platynereis dumerilii (Audouin & Milne Edwards, 1833) were widely distributed in the whole area in most of the samples (25, 21 and 21 respectively). Neither the accumulation, nor the estimator curves (Jacknife1, Jacknife2, Bootstrap) reached the asymptote (Figure S2), suggesting a potential higher number of species ranging from 121 to 167 in the area.

3.3. Patterns of Variation of Polychaete Assemblages Structure

There were significant differences in the structure of polychaete assemblages both among stations and depths, but also for the interaction term station × depth (Table 1), with the significant alongshore variation occurring only at 1.5 m depth, between stations B and A and between stations B and C (Table S3). Moreover, significant differences in the structure of polychaete assemblages were revealed at all stations between 1.5 m and 25 m depth, and between 5 m and 25 m depth, while differences between 1.5 m and 5 m depth were significant only at station B. Results of PERMDISP analyses confirmed that these differences were not barely due to differences in the dispersion of the samples (Table S3).
At all stations, difference in polychaete assemblages between three depths was clearly evident, with assemblages from 1.5 m and 5 m depths being more similar among each other and different from those at 25 m depth (Figure 6). Station A showed the highest scattered distribution among replicates at 25m depth, and the most homogeneous assemblages were those at 5m depth. In fact, the average similarity in species composition and abundance between samples was the highest at 5 m depth (45.22%), intermediate at 25 m depth (39.51%) and the lowest at 1.5 m depth (34.39%) (Table 2A). The most abundant species were the sabellid Amphiglena mediterranea, the nereidids Platynereis dumerilii, Nereis usticensis Cantone, Catalano & Badalamenti, 2003 and Nereis pulsatoria (Savigny, 1822) and the syllids Syllis rosea (Langerhans, 1879), S. pirifera and Exogone dispar (Webster, 1879) at 1.5m depth; the sabellid A. mediterranea, the syllids S. pirifera, E. dispar, Syllis variegata Grube, 1860, Syllis prolifera Krohn, 1852 and Syllis corallicola Verrill, 1900, and nereidids P. dumerilii, Nereis sp. 1 and unidentified juvenile nereidids at 5m depth; and the syllids Syllis armillaris (O.F. Müller, 1776), Syllis gracilis Grube, 1840, Syllis gerlachi (Hartmann-Schröder, 1960) and S. pirifera at 25 m depth (Table 2A). The lowest average dissimilarity in species composition and abundance was between 1.5 m and 5 m depth (68.67%), while it was higher between 1.5 m and 25 m (88.42%) and between 5 m and 25 m depth (79.59%) (Table 2B). Differences in abundance of the most abundant species of Syllidae (i.e., S. rosea, S. prolifera, S. armillaris, S. gracilis, S. gerlachi, S. pirifera), Nereididae (i.e., N. usticensis, P. dumerilii, Nereis sp. 1, Nereididae juv. indet.) and Sabellidae (i.e., A. mediterranea) at different depths, were mainly responsible for these dissimilarities (Table 2B). Amphiglena mediterranea had the highest abundance at 5 m depth, was slightly less abundant at 1.5 m depth and poorly represented at 25 m depth. Syllis armillaris and S. gracilis were very abundant at 25 m depth and poorly represented at 1.5 and 5 m depth. Only low percentages of dissimilarities were due to the differences in taxonomic composition (i.e., presence/absence of species). Namely, the sabellid Amphicorina rovignensis Mikac, Giangrande & Licciano, 2013 characterised 1.5 m depth but was absent at 5 m depth, N. usticensis and S. rosea characterised 1.5 m depth but were absent at 25 m depth, the syllid Odontosyllis ctenostoma Claparède, 1868 characterised 5 m depth, but was absent at 25 m depth, and the sabellid Hypsicomus stichophthalmos (Grube, 1863) characterised 25 m depth but was absent at 5 m depth (Table 2B).
The species that contributed with 4% or more to the dissimilarity between stations at each depth, were mostly the same ones that characterised that depth (Table S4): N. usticensis, N. pulsatoria, S. rosea, S. prolifera, A. mediterranea, S. pirifera, P. dumerilii, E. dispar, Ceratonereis (Composetia) costae (Grube, 1840) and Nereididae juv. indet. at 1.5m depth; A. mediterranea, A. rovignensis, S. prolifera, S. gerlachi, Nereis sp. 1, S. pirifera and P. dumerilii at 5 m depth; and S. armillaris, S. pirifera, H. stichophthalmos, S. gracilis, S. variegata, S. gerlachi and Spirobranchus triqueter (Linnaeus, 1758) at 25 m depth.
When considering the percent coverage of each algal taxon separately, 6 variables (over a total of 36) significantly explained 57.2% of the variation of the polychaete assemblages (Table 3A). Depth explained 28.4% of the variation, while the contribution of orientation was not significant. Overall, the percent coverage of algal taxa significantly explained 28.7% of the observed variation. Peyssonnelia rubra (7.8%), Haliption sp. (7.2%), one taxon of filamentous algae (5.74%), Padina pavonica (3.9%) and Codium efusum (4.1%) were significantly related to distribution of polychaete assemblages. When considering percent coverage of algal functional-morphological groups, 4 out of 9 variables significantly explained 46.7% of the variation of the polychaete assemblages (Table 3B). Among them, depth explained 28.4%, while algal functional-morphological groups all together explained 18.2% of the variation of polychaete assemblages. Orientation was again not significant. In particular, the encrusting (7.3%), foliose (6.7%) and articulated calcareous (4.2%) algae were significantly related to the distribution of polychaete assemblages.

4. Discussion

This is the first study of the hard bottom polychaete assemblages in the North Adriatic Sea over a wide bathymetric range also applying an experimental design that allows describing patterns of spatial distribution in relation to algal coverage. Previous hard bottom studies in the area considered mostly shallower benthic habitats, using qualitative sampling methods making results hardly comparable to ours. In fact, Zavodnik [40,41] reported 38 polychaete species from Cystoseira barbata (Stackhouse) C.Agardh assemblages sampled at 0.5–1 m depth and 69 species examining different brown algae from 0–3 m depths. Later, Amoureux and Katzmann [34] reported 52 species from rocky bottom algal assemblages at 6–9 m depth, while Katzmann [37] collected 93 species from Cystoseira assemblages sampled at different depths ranging from 0–2.5 m. Amoureux [33] found 55 polychaete species among Cystoseira and Peyssonelia algae sampled at 3 m depth. In general, we can notice higher species richness (107 taxa from 22 families) in our study, which can be explained especially by the wider depth range considered herein. Indeed, if we consider only shallower habitats, species richness found herein (62 species at 1.5 m depth and 73 species at 5 m depth) is more similar to what was previously reported for the area. Moreover, our sampling covers the most representative hard-bottom vegetated habitats in the area, supposedly allowing collecting a larger number of polychaete taxa. The macroalgal assemblages found at Sveti Ivan Island are characteristic for the western coast of the North Adriatic Sea, with photophilic algae typically dominating in shallow, and sciaphilic in deeper waters [55,70], and they can be considered as representative habitats to study the diversity of rocky bottom polychaetes in the region. The ecological status of Sv. Ivan island was previously determined using macroalgae grouped in two ecological status groups (ESG I and ESG II) according to Orfanidis et al. [71], and categorized as overall high, while moderate at 1 m depth and good at 3 m depth [54,55]. The algal assemblages herein observed corresponded to those found in the above-mentioned studies. The articulated calcareous, corticated and foliose functional-morphological groups were dominant at 1.5 m depth, the articulated calcareous, encrusting and foliose ones at 5 m depth and the encrusting and filamentous ones at 25 m depth. Due to the strong exposure, the leathery macroalgal group was poorly developed [55,72], and represented only by C. compressa, with very low abundance. From a faunal point of view the present data on polychaetes differ considerably from previous studies, with 23 species recently reported for the first time in the North Adriatic region, including four new records for the Adriatic Sea, as described elsewhere [43,44,73], while the sabellid Amphicorina rovignensis, collected at Sveti Ivan during our survey, was recently described as new to science [44]. Peculiarly, some of the newly reported species were among the most abundant ones in the examined habitats, in particular Nereis usticensis, Syllis corallicola, S. gerlachi, S. rosea and, the most frequent species in this research, S. pirifera. Results of faunal and biogeographic analyses of Syllidae from Sveti Ivan Island indicated possible “meridionalization” of the North Adriatic fauna, i.e., the process of establishment of native Mediterranean warm-water species typical from southern sectors in colder northern areas [43,74]. Two species recorded during our survey are considered alien in the Mediterranean being considered Lessepsian migrants, Nereis persica Fauvel, 1911 and Lepidonotus tenuisetosus (Gravier, 1902) [73,75] and are deposited in the collection of the Natural History Museum in Rijeka (Croatia). Nereis persica was previously reported in the Adriatic Sea only twice, in its northern part [76,77]; however, those records were considered questionable [75]. Elsewhere in the Mediterranean, it was reported along the coast of Israel and Turkey [78,79]. If previous records of this species in the Adriatic Sea are eventually erroneous [75], our finding could represent a north-western widening of the species areal in the Mediterranean. Lepidonotus tenuisetosus was until now reported only in the Eastern Mediterranean, along the coasts of Israel, Egypt, Turkey and Greece [75,79,80,81]. Our finding of L. tenuisetosus represents the first report for the Adriatic Sea and might indicate a north-western widening of the species’ distribution. However, present data are not fully comparable with previous studies from a qualitative point of view, because the knowledge on polychaetes taxonomy is continuously evolving and new species are recorded and newly described also in the Adriatic Sea [73]. New records of some species in our research (such as Nereis usticensis Cantone, Catalano & Badalamenti, 2001 and several syllid species) might arise from the fact that former studies were carried out before those species were scientifically described [43]. Moreover, polychaete systematics changed a lot meanwhile and some taxa are not considered valid at present, such as for example previously reported Pionosyllis serrata Southern, 1914, Syllis brevipennis (Grube, 1863) and Vermiliopsis richardi Fauvel, 1909 [33,37], herein listed as Nudisyllis pulligera (Krohn, 1852), Trypanosyllis (Trypanosyllis) coeliaca Claparède, 1868 and Vermiliopsis labiata (O. G. Costa, 1861) respectively. All this emphasizes the importance of research such as the present one, aiming at updating the knowledge of polychaetes biodiversity of an area.
Our data are more comparable to more recent studies conducted with similar experimental designs in the South Adriatic Sea by Giangrande et al. [82,83]. These authors found 152 polychaete taxa from 22 families in the first study [82] and 118 taxa from 18 families in the second one [83]. Species mostly characterising 5 m and 25 m depth in the south Adriatic [82] were mainly different from those characterising the same depths in our study, with only P. dumerilii being in common at 5 m and S. armillaris and S. gerlachi at 25 m depth. These differences could derive from local differences in the composition of algal assemblages, from different sampling periods (May and November in Giangrande et al. [82] and July in our research) and also from the biogeographic distribution of the Adriatic Sea polychaete fauna. Several species that characterised 5 m and 25 m depth in the south Adriatic were overall absent in our research, in particular Syllis pulvinata (Langerhans, 1881) that characterised 5 m and 25 m depth and S. golfonovensis (Hartmann-Schröder, 1962) and Kefersteinia cirrhata (Keferstein, 1862) that characterised 25 m depth. In fact, S. pulvinata and S. golfonovensis are species that were up to date found only in the southern part of the Adriatic, while it remains difficult to explain the absence of K. cirrata in our samples, since this species is reported as widely distributed in the whole Adriatic Sea [73]. This species is reported as Psamathe fusca Johnston, 1836 in the recent-most polychaete checklist of the Adriatic Sea, based on the synonymy proposed by Pleijel [84]. Although high, the number of taxa from Sveti Ivan Island is considerably lower than that reported in Giangrande et al. [82], possibly because of the smaller sampling surface and less extensive sampling period. Indeed, species area estimator curves suggested potentially higher species richness in our study area. Many species were rare and the analyses of distribution patterns indicating that additional sampling would probably yield more species, as well as it would presumably do sampling in different seasons, considering potential seasonal variability of hard bottom polychaete assemblages [51,58]. Thus, further studies encompassing different seasons and spatial scales should be done in order to upgrade the knowledge of diversity and spatial-temporal variation of hard bottom polychaete assemblages in the North Adriatic region.
The most abundant and species rich families found herein (particularly the Syllidae, but also the Sabellidae, Nereididae, Eunicidae and Serpulidae) are commonly reported as the most characteristic in the Mediterranean hard bottom polychaete assemblages [52,58,82,85,86].
The structure of the polychaete assemblages was highly variable alongshore, but only at the shallowest sites, while it clearly varied bathymetrically, which agrees with the most common trends previously identified, e.g., [11,82]. High variability in the shallowest habitats could be promoted by the high variability of environmental factors (temperature, salinity, hydrodynamics, light intensity etc.) [45]. However, the small sample size used herein might account for differences in presence/absence of several taxa that were represented by few individuals in the overall samples analysed and might have influenced the observed patterns of variation. The highest abundance and diversity and the highest similarity among polychaete assemblages was found at 5 m depth, as well as the highest number of species contributing to the similarity between samples, indicating that assemblages at this depth are the most diverse and structurally complex [87]. The high complexity of the algal coverage at intermediate depth, which can be considered an ecotone where photophilic and sciaphilic conditions coexist, together with a decrease in environmental variability, compared to the shallower habitat, possibly explain the observed increase in diversity [83,88]. In fact, at 5 m depth, the rich coverage of foliose algae (Padina pavonica, Flabelia petiolata and Dictyota dichotoma), structurally complex articulated calcareous algae (Corallina officinalis and Jania sp.) and encrusting algae (P. rubra, P. heteromorpha), with the last two forms being known to entrap considerable quantities of sediment, created altogether a complex and heterogeneous habitat suitable for diverse epifaunal and infaunal polychaete species. We expected the assemblages from 25 m depth to show the highest homogeneity as a consequence of the supposedly more stable environmental conditions. However, the within group similarity was lower at 25 m depth than at 5 m depth, likely because the sites at 25 m depth were situated at the end of the rocky slope, very close to the soft bottom, which may give rise to occasional sedimentation. Enhanced sedimentation combined with reduced light intensity, together with the simplification of the algal-habitat structure, could also contribute to explaining the lower abundance and species richness at 25m depth [48,89]. However, we cannot exclude an effect of the limited sampling effort in the observed patterns of polychaete assemblage variation.
There was a common trend along the whole research area, characterised by differences in the structure of polychaete assemblages between the shallower depths (1.5 m and 5 m) and 25 m depth as already reported elsewhere [82]. The assemblages from 1.5 m and 5 m depths showed significant differences in their structure only at one station. In fact, assemblages from these shallower habitats were commonly characterised by species typical of shallow photophilic environments reported within a variety of algal assemblages (such are Amphiglena mediterranea, Platynereis dumerilii, Syllis prolifera and Exogone dispar) [52,53,90,91], while those from 25 m depth were characterised by species usually also found in sciaphilic habitats (such are S. armillaris and S. gerlachi) [82,90]. The increase of species richness of Serpulidae with depth observed in our research was already reported elsewhere in the Mediterranean [57,92]. It is probably related to the combining effects of the increment of hard (both lithic and organogenic) substrata, low hydrodynamic energy and shadowing, which favour the development of underlying biogenic concretions hosting species with coralligenous affinity [92]. As a whole, the distribution pattern of polychaete assemblages in the examined area appeared related mostly to depth, which, per se, covaries with different environmental variables (temperature, light intensity, hydrodynamics, sedimentation, etc). It was also shown to be related to the algal distribution, being correlation higher considering single algal species than algal functional-morphological groups. The encrusting calcareous alga P. rubra and the articulated calcareous alga Haliption sp. ranked first among algal predictor variables of polychaete distribution, followed by one filamentous and two foliose species. The high percent coverage of articulated calcareous and encrusting algae certainly contributed to explain the high abundance and diversity of polychaetes, due to their high structural complexity that provided wide panoply of suitable microhabitats [88,93]. However, taxa with less complex morphology, i.e., the unidentified filamentous species, also contributed to significantly explain the variability of the polychaete assemblages. In particular, filamentous algal species characterised the habitat at 25 m depth. The spatial complexity of algal thalli surely represents an important factor influencing polychaete distribution, but other non-three-dimensional algal features (e.g., production of antagonistic metabolites, epiphyte colonization, palatability, capability to entrap sediment, etc.) could also be important and deserve further investigations [51,94].
In recent years, the North Adriatic macroalgal assemblages have suffered severe changes (such as reduction in algal coverage, variation in richness and species composition and simplification of the community structure), due to different natural and human driven impacts [95,96], which could likely provoke alterations of the resident polychaete assemblages. However, our algal-based predictor variables explained only part of the observed variability of the polychaete assemblages. We assume that other environmental variables, as well as biotic interactions among polychaetes and between them and other benthic invertebrates (e.g., competition for food or space, recruitment, predation, etc.) may contribute to explain part of the unexplained variability and, thus, should be considered in future studies [11,51,82].
Our results suggest that the diversity of the North Adriatic hard bottom polychaete fauna may be largely underestimated. Further faunal and ecological studies over larger spatial and temporal scale are thus welcome in order to implement our knowledge on diversity and distribution of polychaete assemblages in the area, which will serve as a necessary base to detect changes and predict consequences of natural and anthropogenic disturbances on benthic communities in this important and sensitive Mediterranean sector.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1424-2818/12/10/408/s1, Table S1: Mean coverage (MC) (%) (± standard error (±SE), in italics) of algal species and functional-morphological groups per station and depth. Most characteristic species are marked with asterisk (*); Table S2: Mean abundance (MA) (± standard error (±SE), in italics) of Polychaete taxa per station and depth. In bold values of most abundant taxa on each station and depth; Table S3. Results of PERMANOVA pair-wise and PERMDISP analyses testing for differences in abundance (N), species richness (S), Hill’s species diversity index (N1), Hill’s evenness index (N10) and structure (Stru) of polychaete assemblages: between stations, separately for each depth and between depths, separately for each station. Up, unique perms; t, t-test; P (perm), probability; P (MC), Monte Carlo probability. Significant P-values (p < 0.05) are given in bold; Table S4. Results of SIMPER analyses (cut-off 70%) used to identify taxa that mostly contribute to faunal dissimilarity between stations at each depth. Species contributing to dissimilarity with more than 4% are marked with asterisk. Av.Ab = mean abundance, Diss = mean dissimilarity, Diss/SD = dissimilarity/standard deviation, Contrib% = contribution relative to single taxon, Cum% = cumulative contribution; Figure S1. Distribution of species according to their frequency in the studied samples; Figure S2. Species area accumulation curve (Sobs, Species observed) and estimator curves (Jacknife1, Jacknife2 and Bootstrap).

Author Contributions

The study was conceived and designed by B.M., A.G. and L.M., B.M. and A.J. conducted sampling campaign. Analyses of benthic assemblages and polychaetes were done by B.M., A.G., M.L. and L.M., while Lj.I. analysed algae. Data analysis was performed by B.M. and L.M. The first draft of the manuscript was written by B.M. and all authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by THE GOVERNMENT OF THE REPUBLIC OF ITALY through a research grant for Croatian citizens for the academic year 2007/2008, awarded to the author BM for the execution of the project “Faunal and biogeographic analyses of the Adriatic Sea polychaetes: taxonomic comparison of the polychaetes of Istrian and Salento coasts”.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Buljan, M.; Zore-Armanda, M. Oceanographical properties of the Adriatic Sea. Oceanogr. Mar. Biol. Annu. Rev. 1976, 14, 11–98. [Google Scholar]
  2. McKinney, F.K. The Northern Adriatic Ecosystem: Deep Time in a Shallow Sea; Columbia University Press: New York, NY, USA, 2007; 328p. [Google Scholar]
  3. Degobbis, D.; Gilmartin, M. Nitrogen, phosphorus and biogenic silicon budgets for the northern Adriatic Sea. Oceanol. Acta 1990, 13, 31–45. [Google Scholar]
  4. Boero, F.; Bonsdorff, E. A conceptual framework for marine biodiversity and ecosystem functioning. Mar. Ecol Evol Perspect 2007, 28 (Suppl. 1), 134–145. [Google Scholar] [CrossRef]
  5. Mikac, B.; Semprucci, F.; Guidi, L.; Ponti, M.; Abbiati, M.; Balsamo, M.; Dovgal, I. Newly discovered associations between peritrich ciliates (Ciliophora: Peritrichia) and scale polychaetes (Annelida: Polynoidae and Sigalionidae) with a review of polychaete–peritrich epibiosis. Zool. J. Linn. Soc 2019, 20, 1–15. [Google Scholar] [CrossRef]
  6. Coll, M.; Piroddi, C.; Albouy, C.; Lasram, F.B.R.; Cheung, W.W.L.; Christensen, V.; Karpouzi, V.S.; Guilhaumon, F.; Mouillot, D.; Paleczny, M.; et al. The Mediterranean Sea under siege: Spatial overlap between marine biodiversity, cumulative threats and marine reserves. Glob. Ecol. Biogeogr. 2012, 21, 465–480. [Google Scholar] [CrossRef]
  7. Coll, M.; Piroddi, C.; Steenbeek, J.; Kaschner, K.; Ben Rais Lasram, F.; Aguzzi, J.; Ballesteros, E.; Bianchi, C.N.; Corbera, J.; Dailianis, T.; et al. The Biodiversity of the Mediterranean Sea: Estimates, patterns, and threats. PLoS ONE 2012, 5. [Google Scholar] [CrossRef] [Green Version]
  8. Micheli, F.; Halpern, B.S.; Walbridge, S.; Ciriaco, S.; Ferretti, F.; Fraschetti, S.; Lewison, R.; Nykjaer, L.; Rosenberg, A.A. Cumulative Human Impacts on Mediterranean and Black Sea Marine Ecosystems: Assessing Current Pressures and Opportunities. PLoS ONE 2013, 8, e79889. [Google Scholar] [CrossRef] [Green Version]
  9. European Commission. Directive 2008/56/EC of the European Parliament and of the Council of 17 June 2008 establishing a framework for Community action in the field of marine environmental policy (Marine Strategy Framework Directive). Off. J. Eur. Union 2008, L 164, 25.6.2008, 19–40.
  10. Knox, G.A. The role of Polychaetes in benthic soft-bottom communities. In Essays on Polychaetous Annelids in Memory of Olga Hartmann; Reish, D., Fauchald, K., Eds.; Allan Hancock Foundation: Los Angeles, CA, USA, 1977; pp. 547–604. [Google Scholar]
  11. Musco, L. Ecology and diversity of Mediterranean hard-bottom Syllidae (Annelida): A community-level approach. Mar. Ecol Prog. Ser. 2012, 461, 107–119. [Google Scholar] [CrossRef] [Green Version]
  12. Giangrande, A.; Licciano, M.; Musco, L. Polychaetes as environmental indicators revisited. Mar. Pollut Bull. 2005, 50, 1153–1162. [Google Scholar] [CrossRef]
  13. Aguado-Giménez, F.; Gairin, J.I.; Martinez-Garcia, E.; Fernandez-Gonzalez, V.; Ballester Molto, M.; Cerezo-Valverde, J.; Sanchez-Jerez, P. Application of taxocene surrogation and taxonomic sufficiency concepts to fish farming environmental monitoring. Comparison of BOPA index versus polychaete assemblage structure. Mar. Environ. Res. 2015, 103, 27–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Musco, L.; Terlizzi, A.; Licciano, M.; Giangrande, A. Taxonomic structure and the effectiveness of surrogates in environmental monitoring: A lesson from polychaetes. Mar. Ecol Prog. Ser. 2009, 383, 199–210. [Google Scholar] [CrossRef] [Green Version]
  15. Soares-Gomes, A.; Mendes, C.L.T.; Tavares, M.; Santi, L. Taxonomic sufficiency of polychaete taxocenes for estuary monitoring. Ecol Indic. 2012, 15, 149–156. [Google Scholar] [CrossRef]
  16. Aleffi, F.; Bettoso, N.; Solis-Weiss, V. Spatial distribution of soft–bottom polychaetes along western coast of the northern Adriatic Sea (Italy). Ann. Ser. Hist Nat. 2003, 13, 211–222. [Google Scholar]
  17. Amoureux, L. Inventaire d’une petit collection d’Annélides Polychètes des parages sud de Rovinj (Haute–Adriatique). Thalass. Jugosl. 1976, 12, 381–390. [Google Scholar]
  18. Amoureux, L. Annélides Polychètes recueillies par D. Zavodnik. Thalass. Jugosl. 1983, 19, 1–6. [Google Scholar]
  19. Castelli, A.; Prevedelli, D. Effetto del fenomeno delle mucillaggini dell’estate 1989 sul popolamento a policheti di un microhabitat salmastro preso Punta Marina (Ravenna). Biol. Mar. Meditsuppl. Al Not. S.I.B.M. 1993, 1, 35–38. [Google Scholar]
  20. Crema, R.; Prevedelli, D.; Valentini, A.; Castelli, A. Recovery of the macrozoobenthic community of the Comacchio lagoon system (Northern Adriatic Sea). Ophelia 2000, 52, 143–152. [Google Scholar] [CrossRef]
  21. Fauvel, P. Annélides Polychètes de Rovigno d’Istria. Thalassia 1934, 1, 1–78. [Google Scholar]
  22. Fauvel, P. Annelida Polychaeta della Laguna di Venezia. Bollettino. R. Com. Talassogr. Ital. 1938, 246, 1–27. [Google Scholar]
  23. Fauvel, P. Annélides Polychètes de la Haute Adriatique. Thalassia 1940, 4, 1–24. [Google Scholar]
  24. Gamulin-Brida, H.; Požar, A.; Zavodnik, D. Contributions aux recherches sur la bionomie benthique des fonds meubles de l’Adriatique du nord (II). Biol. Glas 1968, 21, 157–201. [Google Scholar]
  25. Gillet, P. Annélides Polychètes des fonds meubles du Canal de Lim prés de Rovinj (Yugoslavie). Thalassia Jugosl. 1986, 22/21, 127–138. [Google Scholar]
  26. Maggiore, F.; Keppel, E. Biodiversity and distribution of polychaetes and molluscs along the Dese estuary (Lagoon of Venice, Italy). Hydrobiologia 2007, 588, 189–203. [Google Scholar] [CrossRef]
  27. Mikac, B.; Musco, L.; Đakovac, T.; Giangrande, A.; Terlizzi, T. Long–term changes in North Adriatic soft–bottom polychaete assemblages following a dystrophic crisis. Ital. J. Zool. 2011, 78, 304–316. [Google Scholar] [CrossRef]
  28. Mistri, M.; Fano, A.E.; Ghion, F.; Rossi, R. Disturbance and community pattern of polychaetes inhabiting Valle Magnavacca (Valli di Comacchio, Northern Adriatic Sea, Italy). Mar. Ecol. Pszni 2002, 23, 31–49. [Google Scholar] [CrossRef]
  29. Occhipinti-Ambrogi, A.; Savini, D.; Forni, G. Macrobenthos community structural changes off Cesenatico coast (Emilia Romagna, Northern Adriatic), a six–year monitoring programme. Sci. Tot Environ. 2005, 353, 317–328. [Google Scholar] [CrossRef]
  30. Prevedelli, D.; Bellucci, L.G.; Simononi, R.; Ansaloni, I.; Frignani, M.; Ravaioli, M.; Castelli, A. Macrobenthos and environmental characteristics of the Venice lagoon. Atti. Soc. Nat. Mat. Modena 2007, 138, 151–161. [Google Scholar]
  31. Zahtila, E. Offshore polychaete fauna in the northern Adriatic with trophic characteristic. Period. Biol. 1997, 99, 213–217. [Google Scholar]
  32. Zavodnik, D.; Vidaković, J.; Amoureux, L. Contribution to sediment macrofauna in the area of Rovinj (North Adriatic Sea). Cah Biol. Mar. 1985, 26, 431–444. [Google Scholar]
  33. Amoureux, L. Annélides Polychètes de l’îlot Banjole (prés de Rovinj, haute–Adriatique). Cah Biol. Mar. 1975, 16, 231–244. [Google Scholar]
  34. Amoureux, L.; Katzmann, W. Note faunistique et écologique sur une collection d’Annélides Polychètes de substrats rocheux circalittoraux de la région de Rovinj (Yougoslavie). Zool Anz 1971, 186, 114–122. [Google Scholar]
  35. Casellato, S.; Masiero, L.; Sichirollo, E.; Soresi, S. Hidden secrets of the Northern Adriatic: ‘‘Tegnue’’, peculiar reefs. Centr. Eur. J. Biol. 2007, 2, 122–136. [Google Scholar] [CrossRef]
  36. Casellato, S.; Stefanon, A. Coralligenous habitat in the northern Adriatic Sea: An overview. Mar. Ecol. 2008, 29, 321–341. [Google Scholar] [CrossRef]
  37. Katzmann, W. Polychaeten (Errantier, Sedentarier) aus nord—Adriatischen Cystoseira—Beständen und deren Epiphyten. Oecologia 1971, 8, 31–51. [Google Scholar] [CrossRef] [PubMed]
  38. Pitacco, V.; Mavrič, B.; Orlando-Bonaca, M.; Lipej, L. Rocky macrozoobenthos mediolittoral community in the Gulf of Trieste (North Adriatic) along a gradient of hydromorphological modifications. Acta Adriat. 2013, 54, 67–86. [Google Scholar]
  39. Pitacco, V.; Orlando-Bonaca, M.; Mavrič, B.; Lipej, L. Macrofauna associated with a bank of Cladocora caespitosa (Anthozoa, scleractinia) in the Gulf of Trieste (northern Adriatic). Ann. Ser. Hist Nat. 2014, 24, 1–14. [Google Scholar]
  40. Zavodnik, D. Prispevk k poznavanju naselja Cystoseira barbata (Good. & Wood.) C. Ag. v severnem Jadranu. Biol. Vest 1965, 13, 87–101. [Google Scholar]
  41. Zavodnik, D. Dinamika litoralnega fitala na zahodnoistrski obali. Razpr. Slov. Akad. Znan. Umet. 1967, 10, 5–71. [Google Scholar]
  42. Zavodnik, D.; Legac, M.; Gluhak, T. An account of the marine fauna of Pag Island (Adriatic Sea, Croatia). Nat. Croat. 2006, 15, 65–107. [Google Scholar]
  43. Mikac, B.; Musco, L. Faunal and biogeographic analysis of Syllidae (Polychaeta) from Rovinj (Croatia, northern Adriatic Sea). Sci. Mar. 2010, 74, 353–370. [Google Scholar] [CrossRef] [Green Version]
  44. Mikac, B.; Liciano, M.; Giangrande, A. Sabellidae and Fabriciidae (Polychaeta) of the Adriatic Sea with particular retrospect to the Northern Adriatic and the description of two new species. J. Mar. Biol. Assoc. Uk 2013, 93, 1511–1524. [Google Scholar] [CrossRef]
  45. Benedetti-Cecchi, L.; Bulleri, F.; Cinelli, F. The interplay of physical and biological factors in maintaining midshore and low-shore assemblages on rocky coasts in the north-west Mediterranean. Oecologia 2000, 123, 406–417. [Google Scholar] [CrossRef] [PubMed]
  46. Çinar, M.E. Ecological features of Syllidae (Polychaeta) from shallow-water benthic environments of the Aegean Sea, eastern Mediterranean. J. Mar. Biol. Ass Uk 2003, 83, 737–745. [Google Scholar] [CrossRef]
  47. Fraschetti, S.; Terlizzi, A.; Benedetti-Cecchi, L. Patterns of distribution of rocky marine assemblages: Evidence of relevant scales of variation. Mar. Ecol. Prog. Ser. 2005, 296, 13–29. [Google Scholar] [CrossRef] [Green Version]
  48. Irving, A.D.; Connell, S.D. Sedimentation and light penetration interact to maintain heterogeneity of subtidal habitats: Algal versus invertebrate dominated assemblages. Mar. Ecol. Prog. Ser. 2002, 245, 83–91. [Google Scholar] [CrossRef] [Green Version]
  49. Terlizzi, A.; Anderson, M.J.; Fraschetti, S.; Benedetti-Cecchi, L. Scales of spatial variation in Mediterranean subtidal sessile assemblages at different depths. Mar. Ecol. Prog. Ser. 2007, 332, 25–39. [Google Scholar] [CrossRef] [Green Version]
  50. Underwood, A.J.; Chapman, M.G. Scales of spatial patterns of distribution of intertidal invertebrates. Oecologia 1996, 107, 212–224. [Google Scholar] [CrossRef]
  51. Dorgham, M.M.; Hamdy, R.; El-Rashidy, H.H.; Atta, M.M.; Musco, L. Distribution patterns of shallow water polychaetes (Annelida) along the coast of Alexandria, Egypt (eastern Mediterranean). Mediterr. Mar. Sci. 2014, 15, 635–649. [Google Scholar] [CrossRef] [Green Version]
  52. Giangrande, A. Polychaete zonation and its relation to algal distribution down a vertical cliff in the western Mediterranean (Italy): A structural analysis. J. Exp. Mar. Biol. Ecol. 1988, 120, 263–276. [Google Scholar] [CrossRef]
  53. Somaschini, A. Policheti della biocenosi ad alghe fotofile (Facies a Corallina elongata) nel Lazio settentrionale. Atti Sac. Toscana Sci. Nat. Mem Ser. B 1988, 95, 83–94. [Google Scholar]
  54. Iveša, L. Dinamika Popluacija Makrofitobentosa na Hridinastim Dnima uz Zapadnu Obalu Istre. Ph.D. Thesis, Faculty of Science, University of Zagreb, Zagreb, Croatia, 2005. [Google Scholar]
  55. Iveša, Lj.; Lyons, D.M.; Devescovi, M. Assessment of the ecological status of north-eastern Adriatic coastal waters (Istria, Croatia) using macroalgal assemblages for the European Union Water Framework Directive. Aquat. Conserv. Mar. Freshw. Ecosyst. 2009, 19, 14–23. [Google Scholar] [CrossRef]
  56. Boudouresque, C. Méthodes d’étude qualitative et quantitative du benthos (en particulier du phytobenthos). Tethys 1971, 3, 79–104. [Google Scholar]
  57. Casoli, E.; Bonifazi, A.; Ardizzone, G.; Gravina, M.F. How algae influence sessile marine organisms: The tubeworms case of study. Estuar. Coast. Shelf Sci. 2016, 178, 12–20. [Google Scholar] [CrossRef]
  58. Fraschetti, S.; Giangrande, A.; Terlizzi, A.; Miglietta, M.P.; Della Tommasa, L.; Boero, F. Spatio-temporal variation of hydroids and polychaetes associated with Cystoseira amentacea (Fucales: Phaeophyceae). Mar. Biol. 2002, 140, 949–957. [Google Scholar] [CrossRef]
  59. Cormaci, M.; Furnari, G.; Giaccone, G. Macrophytobenthos. In Mediterranean Marine Benthos: A Manual of Methods for Its Sampling and Study; Gambi, M.C., Dappiano, M., Eds.; Biol. Mar. Medit. 11 (1), Chapter 7; Società Italiana di Biologia Marina: Genova, Italy, 2004; pp. 217–246. [Google Scholar]
  60. Littler, M.M.; Littler, D.S. The evolution of thallus form and survival strateges in benthic marine macroalgae: Field and laboratory tests of a functional form model. Am. Nat. 1980, 116, 5–44. [Google Scholar] [CrossRef] [Green Version]
  61. Littler, M.M.; Littler, D.S. Relationships between macroalgal functional form groups and substrata stability in a subtropical rocky-intertidal system. J. Exp. Mar. Biol. Ecol. 1984, 74, 13–34. [Google Scholar] [CrossRef]
  62. Magurran, A.E. Measuring Biological Diversity; Blackwell Publishing: Oxford, UK, 2004; 256p. [Google Scholar]
  63. Colwell, R.K.; Coddington, J.A. Estimating terrestrial biodiversity through extrapolation. In Biodiversity: Measurement and Estimation; Hawksworth, D.L., Ed.; The Royal Society: London, UK, 1994; pp. 101–118. [Google Scholar]
  64. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Aust. J. Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  65. McArdle, B.H.; Anderson, M.J. Fitting multivariate models to community data: A comment on distance-based redundancy analysis. Ecology 2001, 82, 290–297. [Google Scholar] [CrossRef]
  66. Anderson, M.J.; Gorley, R.N.; Clarke, K.R. PERMANOVA+ for PRIMER: Guide to Software and Statistical Methods; PRIMER-E: Plymouth, UK, 2008. [Google Scholar]
  67. Clarke, K.R. Nonparametric multivariate analyses of changes in community structure. Aust. J. Ecol. 1993, 18, 17–143. [Google Scholar] [CrossRef]
  68. Anderson, M.J. DISTLM Forward: A FORTRAN Computer Program to Calculate a Distance-Based Multivariate Analysis for a Linear Model Using Forward Selection; Department of Statistics, University of Auckland: Auckland, New Zealand, 2003. [Google Scholar]
  69. Clarke, K.R.; Gorley, R.N. PRIMER v6: User Manual/Tutorial; PRIMER-E: Plymouth, UK, 2006; 192p. [Google Scholar]
  70. Pérès, J.M.; Gamulin-Brida, H. Biološka Oceanografija. Bentos. Bentoska Bionomija Jadranskog Mora; Školska knjiga: Zagreb, Croatia, 1973; 493p. [Google Scholar]
  71. Orfanidis, S.; Panayotidis, P.; Stamatis, N.V. An insight to ecological evaluation index (EEI). Ecol. Indic. 2003, 3, 27–33. [Google Scholar] [CrossRef]
  72. Munda, I.M. Changes in the benthic algal associations of the vicinity of Rovinj (Istrian coast, North Adriatic) caused by organic wastes. Acta Adriat 1980, 21, 299–332. [Google Scholar]
  73. Mikac, B. A sea of worms: Polychaete checklist of the Adriatic Sea. Zootaxa 2015, 3943, 1–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bianchi, C.N. Biodiversity issues for the forthcoming tropical Mediterranean Sea. Hydrobiologia 2007, 580, 7–21. [Google Scholar] [CrossRef]
  75. Zenetos, A.; Gofas, S.; Verlaque, M.; Çinar, M.E.; García Raso, J.E.; Bianchi, C.N.; Morri, C.; Azzurro, E.; Bilecenoglu, M.; Froglia, C.; et al. Alien species in the Mediterranean Sea by 2010. A contribution to the application of European Union’s Marine Strategy Framework Directive (MSFD). Part I. Spatial distribution. Medit Mar. Sci. 2010, 11/2, 381–493. [Google Scholar] [CrossRef]
  76. Amoureux, L. Annélides Polychètes de l’Adriatique. Thalass. Jugosl. 1983, 19, 15–21. [Google Scholar]
  77. Zavodnik, D.; Kovačić, M. Index of marine fauna in Rijeka Bay (Adriatic Sea, Croatia). Nat. Croat. 2000, 9, 297–379. [Google Scholar]
  78. Ben-Eliahu, M.N. Polychaeta errantia of the Suez Canal. Isr. J. Zool. 1972, 21, 189–237. [Google Scholar]
  79. Çinar, M.E. Alien polychaete species (Annelida: Polychaeta) on the southern coast of Turkey (Levantine Sea, eastern Mediterranean), with 13 new records for the Mediterranean Sea. J. Nat. Hist 2009, 43, 2283–2328. [Google Scholar] [CrossRef]
  80. Barnich, R.; Fiege, D. The Aphroditoidea (Annelida: Polychaeta) of the Mediterranean Sea. Abh. Der. Senckenbergischen Nat. Ges. Frankfut Am. Main 2003, 559, 1–170. [Google Scholar]
  81. Chatzigeorgiou, G.; Faulwetter, S. Polychaetes from Two Subtidal Rocky Shores of the North Coast of Crete, Collected for the NaGISA Project 2007–2008. Version 1.8, 2020, Sampling Event Dataset. Available online: https://www.gbif.org/occurrence/1705514079 (accessed on 13 April 2020). [CrossRef]
  82. Giangrande, A.; Delos, A.L.; Fraschetti, S.; Musco, L.; Licciano, M.; Terlizzi, A. Polychaete assemblages along a rocky shore on the South Adriatic coast (Mediterranean Sea): Patterns of spatial distribution. Mar. Biol. 2003, 143, 1109–1116. [Google Scholar] [CrossRef]
  83. Giangrande, A.; Delos, A.L.; Musco, L.; Licciano, M.; Pierri, C. Polychaete assemblages of rocky shore along the South Adriatic coast (Mediterranean Sea). Cah Biol. Mar. 2004, 45, 85–95. [Google Scholar]
  84. Pleijel, F. Phylogeny and classification of Hesionidae (Polychaeta). Zool Scr. 1998, 27, 89–163. [Google Scholar] [CrossRef]
  85. Çinar, M.E.; Gönlügür-Demirci, G. Polychaete assemblages on shallow-water benthic habitats along the Sinop Peninsula (Black Sea, Turkey). Cah Biol. Mar. 2005, 46, 253–263. [Google Scholar]
  86. Gambi, M.C.; Musco, L.; Giangrande, A.; Badalamenti, F.; Micheli, F.; Kroeker, K.J. Distribution and functional traits of polychaetes in a CO2 vent system: Winners and losers among closely related species. Mar. Ecol. Prog. Ser. 2016, 550, 121–134. [Google Scholar] [CrossRef] [Green Version]
  87. Dahl, L.; Dahl, K. Temporal, spatial and substrate-dependent variations of Danish hard-bottom macrofauna. Helgol. Mar. Res. 2002, 56, 159–168. [Google Scholar] [CrossRef]
  88. Serrano, A.; Preciado, I. Environmental factors structuring polychaete communities in shallow rocky habitats: Role of physical stress versus habitat complexity. Helgol Mar. Res. 2007, 61, 17–29. [Google Scholar] [CrossRef] [Green Version]
  89. Airoldi, L. The effects of sedimentation on rocky coast assemblages. Oceanogr Mar. Biol. Ann. Rev. 2003, 41, 161–236. [Google Scholar]
  90. Sardá, R. Polychaete communities related to plant covering in the midlittoral and infralittoral zones of the Balearic Islands (Western Mediterranean). Pszni Mar. Ecol. 1991, 12, 341–360. [Google Scholar] [CrossRef]
  91. López, E.; Viéitez, J.M. Polychaete assemblages on non-encrusting infralittoral algae from the Chafarinas Islands (SW Mediterranean). Cah Biol. Mar. 1999, 40, 375–384. [Google Scholar]
  92. Sanfilippo, R.; Rosso, A.; Sciuto, F.; Serio, D.; Catra, M.; Alongi, G.; Negri, M.P.; Leonardi, R.; Viola, A. Serpulid polychaetes from Cystoseira communities in the Ionian Sea, Mediterranean. Vie Et Milieu 2017, 67, 217–226. [Google Scholar]
  93. Tena, J.; Capaccioni-Azzati, R.; Torres-Gravila, F.J.; Garcia-Carrascosa, A.M. Polychaetes associated with different facies of the photophilic algal community in the Chafarinas archipelago (SW Mediterranean). Bull. Mar. Sci. 2000, 67, 55–72. [Google Scholar]
  94. Cacabelos, E.; Olabarria, C.; Incera, M.; Troncoso, J.S. Effects of habitat structure and tidal height on epifaunal assemblages associated with macroalgae. Estuar Coast. Shelf Sci. 2010, 89, 43–52. [Google Scholar] [CrossRef]
  95. Falace, A.; Alongi, G.; Cormaci, M.; Furnari, G.; Curiel, D.; Cecere, E.; Petrocelli, A. Changes in the benthic algae along the Adriatic Sea in the last three decades. Chem. Ecol. 2010, 26, 77–90. [Google Scholar] [CrossRef]
  96. Marchini, A.; Ragazzola, F.; Vasapollo, C.; Castelli, A.; Cerrati, G.; Gazzola, F.; Jiang, C.; Langeneck, J.; Manauzzi, M.C.; Musco, L.; et al. Intertidal Mediterranean Coralline Algae Habitat Is Expecting a Shift Toward a Reduced Growth and a Simplified Associated Fauna Under Climate Change. Front. Mar. Sci. 2019, 6, 106. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the study area and sampling stations (AC). Arrows indicating city of Rovinj area and Sveti Ivan Island.
Figure 1. Location of the study area and sampling stations (AC). Arrows indicating city of Rovinj area and Sveti Ivan Island.
Diversity 12 00408 g001
Figure 2. Mean total percent coverage of algae and algal functional-morphological groups, per stations and depths.
Figure 2. Mean total percent coverage of algae and algal functional-morphological groups, per stations and depths.
Diversity 12 00408 g002
Figure 3. Representative algal assemblages at: (a) 1.5 m depth; (b) 5 m depth; (c) 25 m depth.
Figure 3. Representative algal assemblages at: (a) 1.5 m depth; (b) 5 m depth; (c) 25 m depth.
Diversity 12 00408 g003
Figure 4. (a) Mean (±SE) total abundance (N; lines) and species richness (S; bars); (b) Mean (±SE) Hill’s species diversity index (N1; bars) and Hill’s evenness index (N10; lines) at each studies station and depth.
Figure 4. (a) Mean (±SE) total abundance (N; lines) and species richness (S; bars); (b) Mean (±SE) Hill’s species diversity index (N1; bars) and Hill’s evenness index (N10; lines) at each studies station and depth.
Diversity 12 00408 g004
Figure 5. Cumulative abundance (N; on the left side) and species richness (S; on the right side) of polychaetes per family at the three studied depths. CI—Cirratulidae, EN—Eunicidae, GL—Glyceridae, NE—Nereididae, PH—Phyllodocidae, PO—Polynoidae, SA—Sabellidae, SE—Serpulidae, SY—Syllidae, TE—Terebellidae, OTH—other families.
Figure 5. Cumulative abundance (N; on the left side) and species richness (S; on the right side) of polychaetes per family at the three studied depths. CI—Cirratulidae, EN—Eunicidae, GL—Glyceridae, NE—Nereididae, PH—Phyllodocidae, PO—Polynoidae, SA—Sabellidae, SE—Serpulidae, SY—Syllidae, TE—Terebellidae, OTH—other families.
Diversity 12 00408 g005
Figure 6. Nonmetric multidimensional scaling (MDS) ordination plot based on Bray–Curtis similarity of untransformed data, comparing structure of polychaete assemblages between samples. AC = stations.
Figure 6. Nonmetric multidimensional scaling (MDS) ordination plot based on Bray–Curtis similarity of untransformed data, comparing structure of polychaete assemblages between samples. AC = stations.
Diversity 12 00408 g006
Table 1. Results of PERMANOVA and PERMDISP analyses (untransformed data) testing for differences in abundance (N), species richness (S), Hill’s species diversity index (N1), Hill’s evenness index (N10) and structure (Stru) of polychaete assemblages between stations (st) and depths (de). df, degrees of freedom; SS, sum of squares; MS, mean squares; Up, unique perms; F, F-ratio; P (perm), probability. Significant p-values (p  <  0.05) are given in bold.
Table 1. Results of PERMANOVA and PERMDISP analyses (untransformed data) testing for differences in abundance (N), species richness (S), Hill’s species diversity index (N1), Hill’s evenness index (N10) and structure (Stru) of polychaete assemblages between stations (st) and depths (de). df, degrees of freedom; SS, sum of squares; MS, mean squares; Up, unique perms; F, F-ratio; P (perm), probability. Significant p-values (p  <  0.05) are given in bold.
PERMANOVAPERMDISP
SourcedfSSMSPseudo-FUpP (perm)df1df2FP (perm)
st25442.32721.14.753999530.01922242.43510.1038
Nde24607.22303.64.331660660.0995 --
stxde42127.3531.810.9290899560.4792 --
Res1810303572.41
st2178.7489.376.486699490.00762240.34130.8379
Sde2627.63313.8118.97638250.02542240.23660.7997
stxde466.14816.5371.200399500.3431 --
Res1824813.778
st2101.9550.9747.959499550.00282241.07830.4925
N1de2282.77141.389.794960910.04592243.26310.0885
stxde457.73714.4342.253899500.0995 --
Res18115.286.4043
st20.08090.04054.29399500.03042241.13210.4103
N10de20.00650.00320.2507560860.7786 --
stxde40.05170.01291.371899370.2792 --
Res180.16960.0094
st27282.23641.12.507799000.00012244.25130.0549
Strude226808134044.378461140.01762242.47350.1668
stxde4122463061.42.108498430.00018188.04710.0122
Res18261361452
Table 2. Results of SIMPER analyses (cut-off 70%) used to identify taxa that mostly contribute to (A) faunal similarity within depths, (B) faunal dissimilarity between depths. Abund = mean abundance, Sim% = mean similarity, Sim/SD = similarity/standard deviation, Contrib% = contribution relative to single taxon, Cum% = cumulative contribution, Av.Ab = mean abundance, Diss% = mean dissimilarity, Diss/SD = dissimilarity/standard deviation.
Table 2. Results of SIMPER analyses (cut-off 70%) used to identify taxa that mostly contribute to (A) faunal similarity within depths, (B) faunal dissimilarity between depths. Abund = mean abundance, Sim% = mean similarity, Sim/SD = similarity/standard deviation, Contrib% = contribution relative to single taxon, Cum% = cumulative contribution, Av.Ab = mean abundance, Diss% = mean dissimilarity, Diss/SD = dissimilarity/standard deviation.
GroupSpecies AbundSim%Sim/ SDContrib%Cum%
ASim%
1.5 mAmphiglena mediterranea 9.678.141.6323.6823.68
34.39% Platynereis dumerilii 6.115.411.4815.7239.40
Nereis usticensis 11.002.450.337.1346.54
Syllis rosea 7.002.451.117.1253.66
Sphaerosyllis pirifera 4.672.430.717.0660.72
Nereis pulsatoria 3.112.351.056.8567.56
Exogone dispar 2.891.620.654.7072.26
5 mAmphiglena mediterranea 17.5610.461.4323.1423.14
45.22% Sphaerosyllis pirifera 5.674.061.468.9732.11
Nereididae juv. indet. 4.003.042.116.7138.82
Nereis sp. 1 4.562.921.376.4645.28
Exogone dispar 3.672.792.136.1751.45
Syllis variegata 3.332.551.335.6357.08
Platynereis dumerilii 3.562.501.175.5262.60
Syllis prolifera 5.332.120.714.6967.29
Syllis corallicola 3.112.082.294.5971.88
25 mSyllis armillaris 13.2213.051.5533.0433.04
39.51% Syllis gracilis 6.678.381.8821.2154.25
Sphaerosyllis pirifera 8.785.930.9815.0169.27
Syllis gerlachi 4.113.651.449.2478.50
GroupsSpeciesAv. AbAv.AbDiss%Diss/SDContrib%Cum%
B Diss% Group 1.5 Group 5
1.5 and 5Amphiglena mediterranea9.6717.567.341.3510.6910.69
68.67%Nereis usticensis11.000.115.990.718.7219.41
Syllis rosea7.000.444.010.585.8325.24
Syllis prolifera5.445.333.891.165.6630.91
Sphaerosyllis pirifera4.675.673.381.134.9235.83
Platynereis dumerilii6.113.562.531.213.6839.51
Nereis sp. 11.114.562.391.223.4842.99
Nereididae juv. indet.2.004.002.161.473.1546.14
Syllis variegata0.333.332.011.442.9249.07
Exogone dispar2.893.671.951.292.8451.90
Syllis corallicola1.003.111.891.692.7554.65
Syllis gerlachi0.223.001.880.762.7457.39
Nereis pulsatoria3.111.001.621.182.3559.74
Amphicorina rovignensis0.003.331.460.412.1361.87
Nereis rava0.782.671.461.302.1364.00
Ceratonereis (Composetia) costae1.112.001.420.922.0666.07
Odontosyllis ctenostoma0.332.441.360.861.9868.05
Lysidice unicornis2.330.441.141.111.6669.71
Syllis armillaris1.001.891.110.981.6271.33
Group 1.5Group 25
1.5 and 25Syllis armillaris1.0013.229.061.7510.2410.24
88.42%Amphiglena mediterranea9.670.337.291.538.2418.48
Nereis usticensis11.000.007.200.708.1426.63
Sphaerosyllis pirifera4.678.786.061.136.8533.48
Syllis rosea7.000.005.260.615.9539.43
Syllis gracilis0.226.675.091.735.7545.18
Platynereis dumerilii6.110.444.621.265.2350.41
Syllis prolifera5.441.003.820.864.3254.73
Syllis gerlachi0.224.112.871.313.2457.97
Hypsicomus stichophthalmos0.113.112.470.432.7960.76
Nereis pulsatoria3.110.112.451.262.7763.54
Exogone dispar2.890.112.340.992.6566.18
Lysidice unicornis2.330.111.621.251.8368.01
Nereididae juv. indet.2.000.111.430.741.6269.63
Dodecaceria concharum1.670.221.300.981.4771.10
Group 5Group 25
5 and 25Amphiglena mediterranea17.560.3311.401.5214.3314.33
79.59%Syllis armillaris1.8913.227.941.629.9724.30
Sphaerosyllis pirifera5.678.785.071.236.3730.66
Syllis gracilis0.676.674.291.655.4036.06
Syllis prolifera5.331.003.241.134.0840.14
Nereis sp. 14.560.562.891.263.6343.77
Nereididae juv. indet.4.000.112.731.573.4347.20
Syllis gerlachi3.004.112.641.113.3250.52
Exogone dispar3.670.112.591.383.2553.77
Platynereis dumerilii3.560.442.461.123.0956.86
Hypsicomus stichophthalmos0.003.112.170.412.7259.59
Syllis variegata3.331.561.891.172.3861.96
Syllis corallicola3.110.441.751.432.2064.17
Amphicorina rovignensis3.330.221.670.442.1066.26
Odontosyllis ctenostoma2.440.001.560.851.9668.22
Nereis rava2.671.111.551.161.9570.17
Table 3. Results of DISTLM-forward analysis. (A) Variables: percent coverage of each algal taxa (5% cut-off), depth and orientation), (B) Variables: percent coverage of algal functional-morphological groups, depth and orientation). Only variables significantly contributing to explain variation of polychaete assemblages (p < 0.05) are presented. Prop.: proportion of explained variation; Cumul.: cumulative proportion of explained variation.
Table 3. Results of DISTLM-forward analysis. (A) Variables: percent coverage of each algal taxa (5% cut-off), depth and orientation), (B) Variables: percent coverage of algal functional-morphological groups, depth and orientation). Only variables significantly contributing to explain variation of polychaete assemblages (p < 0.05) are presented. Prop.: proportion of explained variation; Cumul.: cumulative proportion of explained variation.
VariableSS(trace)Pseudo-FPProp.Cumul.
Depth206179.93940.00010.284480.28448
APeyssonnelia rubra5681.72.95320.00050.0783980.36287
Haliptilon sp.5214.12.92790.00080.0719470.43482
Filamentous sp. 14161.52.4880.0010.0574230.49224
Padina pavonica2822.71.74470.02630.0389490.53119
Codium efusum2954.41.90480.02080.0407660.57196
Depth206179.93940.00010.284480.28448
BEncrusting5292.22.72770.00070.0730240.3575
Foliose4853.62.67640.00110.0669710.42447
Articulated calcareous3064.31.74440.03060.0422820.46675
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mikac, B.; Licciano, M.; Jaklin, A.; Iveša, L.; Giangrande, A.; Musco, L. Diversity and Distribution Patterns of Hard Bottom Polychaete Assemblages in the North Adriatic Sea (Mediterranean). Diversity 2020, 12, 408. https://0-doi-org.brum.beds.ac.uk/10.3390/d12100408

AMA Style

Mikac B, Licciano M, Jaklin A, Iveša L, Giangrande A, Musco L. Diversity and Distribution Patterns of Hard Bottom Polychaete Assemblages in the North Adriatic Sea (Mediterranean). Diversity. 2020; 12(10):408. https://0-doi-org.brum.beds.ac.uk/10.3390/d12100408

Chicago/Turabian Style

Mikac, Barbara, Margherita Licciano, Andrej Jaklin, Ljiljana Iveša, Adriana Giangrande, and Luigi Musco. 2020. "Diversity and Distribution Patterns of Hard Bottom Polychaete Assemblages in the North Adriatic Sea (Mediterranean)" Diversity 12, no. 10: 408. https://0-doi-org.brum.beds.ac.uk/10.3390/d12100408

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop