Next Article in Journal
Phylogenomic, Biogeographic, and Evolutionary Research Trends in Arachnology
Next Article in Special Issue
Marine Protected Areas Management in the Mediterranean Sea—The Case of Croatia
Previous Article in Journal
Can Thailand Protect 30% of Its Land Area for Biodiversity, and Will This Be Enough?
Previous Article in Special Issue
Bycatch Dynamics from a Small-Scale Shrimp Trap Fishery in the Mediterranean Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insights into a Mediterranean Sea Benthic Habitat: High Diversity of Epiphytic Bryozoan Assemblages on Phyllophora crispa (Rhodophyta) Mats

1
Marine Ecology Department, Faculty of Biology and Chemistry, University of Bremen, 28359 Bremen, Germany
2
Department of Environmental Biology, Sapienza University of Rome, 00185 Rome, Italy
3
Animal Ecology, Department of Ecology and Genetics, Uppsala University, 75236 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(5), 346; https://doi.org/10.3390/d14050346
Submission received: 4 March 2022 / Revised: 12 April 2022 / Accepted: 28 April 2022 / Published: 28 April 2022
(This article belongs to the Special Issue Marine Biodiversity and Ecosystems Management)

Abstract

:
With its geographically isolated location and geological history, the Mediterranean Sea harbors well-known biodiversity hotspots, such as Posidonia oceanica seagrass meadows. Recently, long-living mats formed by the fleshy red alga Phyllophora crispa have been described to be associated with a high diversity of sessile invertebrates in the Tyrrhenian Sea. One of the key taxa among these sessile invertebrates are bryozoans: their abundance, diversity, and spatial distribution in P. crispa mats represent a gap in scientific knowledge. Thus, we conducted a pilot study on bryozoan assemblages associated with P. crispa mats around Giglio Island (Tuscan Archipelago, Italy) in 2018, followed by a comparative study on four sites distributed around the island in the subsequent year, 2019. We compared these findings to bryozoan abundance and diversity on P. oceanica shoots and leaves during the second expedition. The findings revealed more than 46 families, with a significantly higher number of taxa identified in P. crispa mats (33) than in P. oceanica meadows (29). The Shannon diversity index was similar between P. crispa and P. oceanica shoots, while Pielou’s evenness index was lower in P. crispa mats. The most abundant families reported across all habitats were Crisiidae, Aetidae, and Lichenoporidae; but the most abundant family on P. crispa was Chlidoniidae (Chlidonia pyriformis). The assemblages associated with P. crispa differed among sites, with higher abundances but lower diversity on the exposed southernmost site. The total bryozoan abundance was significantly higher on P. crispa (average 2.83 × 106 ± 1.99 × 106 colonies per m2 seafloor) compared to P. oceanica meadows (average 0.54 × 106 ± 0.34 × 106 colonies per m2 seafloor). Our results show a high diversity of bryozoans on P. crispa thalli compared to P. oceanica meadows, which was consistent throughout the study. These findings confirm the value of the red alga-generated habitat for associated bryozoans and may have implications for future biodiversity assessments and conservation measures.

1. Introduction

Understanding the patterns of variability of benthic assemblages represents one of the main goals for ecologists and has a pivotal role in managing and conserving marine habitats [1]. This information might help scientists predict or understand organisms’ responses to global environmental changes. Hotspots of biodiversity in the Mediterranean Sea are often created by key engineering species that provide structurally complex habitats for associated communities [2]. Some of the most-studied engineered habitats in the Mediterranean Sea are seagrass meadows of the endemic angiosperm Posidonia oceanica (L.) Delile, 1813 (Figure 1C), and coralligenous reefs, mainly formed by the accumulation of crustose Rhodophytes thalli [3,4,5,6]. The rooted, flowering plants of P. oceanica build up dense meadows consisting of two sub habitats: a canopy of up to 50 cm long leaves that grows from a dense network of stems and roots [7]. The crustose structures of coralligenous reefs provide a rigid substrate, characterized by systems of canals and crevices [5]. Both of these ecosystems promote high biodiversities of sessile invertebrates, as they provide different environmental gradients (e.g., light intensity, water movement, food availability), as well as shelter and space for larval settlement [8,9,10] (Figure 1D,E).
Bryozoans are filter feeders and form one of the most abundant and diverse groups of epiphytic invertebrates on host organisms, such as seagrass (i.e., P. oceanica) and macroalgae [10,11,12,13,14,15]. In particular, the richest bryozoan diversity in the Mediterranean Sea has been reported on coralligenous reefs and marine caves due to the availability of several microhabitats that enhance the presence of bryozoans characterized by different shapes and ecological traits [16,17,18]. Their typical colonial structures consist of often highly specialized zooids and may vary between thin crusts, erect and branched forms, or larger rigid structures [19,20]. Many species are considered bioindicators for environmental changes, as they often respond faster to environmental and human-mediated pressures [21,22,23]. Because of their calcium carbonate hulls, they are sensitive to ocean acidification [24,25,26] and hold an essential part in fossil records [27]. Some bryozoan species have been reported as habitat-forming organisms playing a pivotal role in promoting biodiversity [24], e.g., by overgrowing P. oceanica shoots and preventing the settlement of other species [28] (Figure 1D). Furthermore, bryozoans play an essential role as primary consumers by transferring particulate organic matter from the water column into the benthic community [29,30].
The red macroalgae Phyllophora crispa ((Hudson) P.S.Dixon, 1964) is known for forming dense mats and hosting a diverse community of epiphytic invertebrates in the Black Sea [31]. In the north-western Mediterranean Sea, P. crispa forms dense mats [32,33,34] (Figure 1A), which have recently been shown to host high diversities of invertebrate fauna. Especially epiphytic filter feeders (e.g., Bryozoa, Serpulidae), which benefit from the algal thalli as a substrate and accumulated food particles from the water column, have been identified to contribute to the associated biodiversity [34,35,36] (Figure 1B). However, little is known about the variations of bryozoan diversity inside P. crispa mats over space and time or their diversity compared to other Mediterranean habitats (e.g., P. oceanica meadows).
To address this knowledge gap, we carried out a comparative field study on the abundance and diversity of bryozoans at four locations of P. crispa mats around Giglio Island (Tuscan Archipelago, Italy, Tyrrhenian Sea) in two consecutive seasons (2018 and 2019) and on a P. oceanica meadow as a reference habitat (2019). We decided on P. oceanica as a reference habitat because its biodiversity has been well studied during the last decades, and the structure is more similar to the fleshy P. crispa thalli than, e.g., the calcareous substrate of coralligenous reefs. Additionally, extensive coralligenous reefs are generally found deeper, and other mat-forming macroalgae are not present at the same depth as P. crispa mats within the study area. With this work, we aim to answer the following research questions:
  • What are the abundances and diversity of bryozoans inside P. crispa mats compared to P. oceanica meadows?
  • Which are the most abundant families in the investigated habitats, and which families are unique to P. crispa mats?
  • What is the spatial variability of the bryozoan assemblages inside P. crispa mats?

2. Materials and Methods

2.1. Location and Sampling Procedure

The study area is located at the island of Giglio (42°21′19.4″ N 10°54′06.1″ E, Figure 2) and is characterized by steep granite slopes alternating with sandy bottoms. The infralittoral seabeds are colonized by Posidonia oceanica meadows, Phyllophora crispa mats, and coralligenous reefs [35]. Scientific SCUBA divers collected all samples at water depths of 30 ± 4 m at four sites around the island. Sampling took place between May and June 2018 (21 P. crispa samples on Site Mix) and between May and July 2019 (4 P. crispa samples per site, resulting in 16 total samples; 9 P. oceanica leaves, and 10 P. oceanica shoots on Site Mix). The sampling sites were chosen for their similar topography and occurrences of P. crispa mats of at least 90% coverage (by visual census; Appendix A Figure A1) at the target depth of 30 m. In 2018, we sampled P. crispa material from Site Mix. In 2019, all four sites (site PC1, PC2, PC3, and Site Mix) were sampled for P. crispa mats, and Site Mix was additionally sampled for P. oceanica material. Site Mix was the only site with a continuous P. oceanica meadow present at the target depth.
The P. crispa mats were sampled using a metal quadrat (size 30 × 30 cm) that was placed randomly inside a continuous mat of at least 5 cm thickness to define the exact sampling area. All material of these main samples was then carefully scraped off the rock surface with a spatula directly under the holdfast to avoid breaking or removing epiphytic organisms.
The P. oceanica leaves and shoots were cut with scissors, directly on the sheath or rhizome branching point, on meadows close to the sampling depth of P. crispa mats (30 ± 4 m). Additionally, we counted the number of P. oceanica shoots per m2 (n = 16 counts) and the number of leaves per shoot (n = 32 counts) for later extrapolating bryozoan colonies per m2 seafloor. These density measures were carried out by counting the number of shoots within a 40 × 40 cm frame. Shoots and leaves were treated as separate sub-habitats considering their different ecological traits, particularly regarding their longevity as a fundamental trait for larval settlement [37,38].
All main samples were carefully transferred into plastic jars immediately after sampling. Every jar contained approximately one-third of sampled material and two-thirds of seawater to avoid oxygen depletion during the transport to the holding facilities in the Institute for Marine Biology (IfMB, located in the near bay of Campese). The main samples were then kept in aerated seawater tanks at constant temperature (18 °C; equivalent to in situ temperature) before they were analyzed within three days after sampling.

2.2. Species Identification and Abundance Assessment

From the P. crispa main samples, subsamples between 20 and 100 g wet weight were extracted, roughly 20% of the respective main sample. The span of different wet weights resulted from the different amounts of main samples due to the randomly picked sampling area. The remaining material of the main samples was analyzed for other taxonomic groups that are not within the frame of this study. The shoots and leaves of P. oceanica were analyzed as a whole under stereo magnifiers (maximum 40× magnification). Bryozoan colonies were identified to the lowest possible taxonomic level using the relevant literature (Appendix A, Table A2). The abundance was assessed by counting the number of colonies. In the case of branching, stolonal taxa—e.g., Chlidonia pyriformis (Bertoloni, 1810)—all parts connected by stolons were considered one colony. We opted to work with the family level for further analysis to reduce observer bias and increase work efficiency. Previous studies have shown that this method may not result in a significant loss of information [39]. Furthermore, the taxonomic sufficiency hypothesis applied to Mediterranean peculiar habitats revealed that surrogate taxonomic levels higher than species could be used to highlight the diversity pattern of benthic assemblages [40].
The surface area of all subsamples was assessed as follows. For P. oceanica shoots, the length and diameter were measured, and surface area calculations were based on an assumed cylindrical shape. For P. oceanica leaves, the length and width were measured to calculate the rectangular surface (times two, to account for both sides of the leaf). The P. crispa subsamples were flattened with a glass pane on laminated graph paper before being photographed from above with a fixed tripod. The surface area was then determined using ImageJ (version 1.52o, https://imagej.nih.gov/ij/, accessed on 23 April 2019) and multiplied by two to account for both sides of the thalli. The wet weight of all P. crispa main samples and subsamples was assessed after shaking off excess water to extrapolate the abundances from the subsamples to the main sample and finally to the surface of the seafloor. Bryozoan abundance was then calculated as the number of colonies per m2 of seafloor ± standard deviation (SD) (Appendix A, Formulas (A1)–(A3)).

2.3. Diversity Descriptors and Statistical Analysis

Bryozoan diversity was assessed using four descriptors: total numbers of families per site and habitat, Shannon diversity index [41], and Pielou’s evenness index [42]. The descriptors were calculated as means per site and habitat and were reported with the respective SD.
Differences in the composition of bryozoan assemblages in P. crispa mats among sites were tested using multivariate permutational analysis of variance (PERMANOVA [43]). The source of significant results (p < 0.05) was tested using Tukey’s honestly significant difference (HSD) test. Based on the results of the multivariate analysis, the data of the northern sites were pooled for further comparison among habitats (Appendix A, Table A1).
Statistical differences in the diversity descriptors among sites and (sub-) habitats were assessed with pairwise Wilcoxon–Mann–Whitney tests. The comparison among habitats on the northern sites was visualized using non-metric multidimensional scaling (NMDS). In addition, the bryozoan assemblages were clustered with a Spearman ranked correlation (average linkage) on a family level to visualize the composition of the bryozoan assemblages using the software’ heatmapper’ [42]. Analyses and plots were made with R (version 4.0.5) [44].
To avoid potential variability between the observer and different sampling efforts in the two subsequent seasons, we decided against a direct statistical comparison on a temporal scale.

3. Results

3.1. Bryozoan Richness and Abundance

A total of 17,822 bryozoan colonies were found (2018 = 10,312; 2019 = 7510) and 46 taxa identified to at least the family level. The highest number of families was identified in Phyllophora crispa mats in 2019 (33), while Posidonia oceanica meadows held 29 families (Figure 3). A similar amount was found in P. crispa mats during the pilot study in 2018 (28). In the comparative study of 2019, the two habitats shared 21 families, while this number was higher for P. crispa and P. oceanica shoots (20) than for P. crispa and P. oceanica leaves (13). A significantly higher number of families (p < 0.05) was reported in P. crispa compared to P. oceanica sub-habitats (Figure 4F).
The most abundant families across all habitats were Crisiidae, Aetidae, and Lichenoporidae (Table 1). The highest density of bryozoan colonies was found in P. crispa mats in 2018 (average 2,827,762 ± 1,984,965 colonies per m2 seafloor). Among P. crispa sites, the density was highest on the southernmost site (site PC3) and lowest on the northernmost site (site PC2, Figure 4A). This observation was the only significant effect among sites in the diversity descriptors of P. crispa mats and confirmed the results of the multivariate analysis (Appendix A, Table A1). Based on these results, the data of the northern P. crispa sites were pooled for further comparison among (sub-) habitats. The lowest abundance was recorded on P. oceanica leaves (average 177,912 ± 104,999 colonies per m2 seafloor, Figure 4E). The most abundant family contributing to the measured densities was Chlidoniidae (Chlidonia pyriformis) in P. crispa mats and P. oceanica shoots, with a higher abundance on P. crispa. Candidae were abundant on P. crispa mats and P. oceanica shoots as well. On the leaves of P. oceanica, Haplopomidae and Tubuliporidae were most abundant (Table 1). It is also notable that P. crispa mats host additional growth types, such as petraliform (e.g., Beania hirtissima (Heller, 1867)), encrusting (e.g., Watersipora sp.), and creeping (e.g., Aeta sp.), while P. oceanica leaves mainly support encrusting forms (e.g., Haplopoma sp.).

3.2. Diversity Indices

The Shannon index was highest in P. crispa mats in 2019 at sites PC1 and PC2 ((mean 2.2 ± 0.1), Figure 4C), while in the direct comparison of the northern P. crispa sites to P. oceanica on Site Mix, the P. oceanica shoots showed slightly higher values (mean 2.3 ± 0.4) than P. crispa ((mean 2.1 ± 0.2) Figure 4C).
The evenness (Pielou’s index) was similar on P. oceanica leaves and shoots (mean 0.20 ± 0.01, and 0.18 ± 0.01, respectively), compared to P. crispa (mean 0.20 ± 0.01) on the northern sites. The northern sites showed slightly higher values, while the southern site was lower, without significant effects (Figure 4D).

3.3. Structure of Bryozoan Assemblages

The diversity inside the P. crispa mats consisted of many taxa unique to this habitat. Out of the 46 families identified during this study, 18 were exclusively found on P. crispa. In the P. oceanica samples, we found 3 families not present on P.crispa (Figure 5). Most of the taxa solely found on P. crispa belonged to Watersiporidae (212.608 ± 37.342 and 24.788 ± 5.333 colonies per m2 in 2018 and 2019, respectively).
The cluster analysis highlights the differences between P. oceanica sub-habitats and P. crispa and between sampling years (Figure 5). It also shows the dominance of C. pyriformis (the only species of Chlioniidae found) and Crisiidae (mainly Filicrisia geniculate and Crisia sp.) in P. crispa mats and on P. oceanica shoots. The family of Tubuliporidae was abundant across samples and habitats. Between the two years of P. crispa sampling, Watersiporidae and Candidae showed higher abundances in 2018, while Aetidae were more abundant in 2019 (Figure 5, Table 1).
The non-metric multidimensional scaling (NMDS) further highlights the differences in the composition of bryozoan families among the (sub-) habitats in 2019 (only northern sites; Figure 6). While all three (sub-) habitats form distinct clusters, the P. oceanica sub-habitats samples show a higher degree of scattering than the P. crispa samples.

4. Discussion

4.1. Differences in Bryozoan Abundances and Diversity between Phyllophora crispa Mats and Posidonia oceanica Sub-Habitats

Our results show differences among the three (sub-) habitats, with a clear trend of higher abundances and diversity of the bryozoan community inside the P. crispa mats (Figure 4E–H). The results have shown that the bryozoan assemblages in P. crispa mats and on P. oceanica shoots are similar concerning the chosen diversity descriptors (Figure 4G,H), whereas the distinct clusters in the NMDS analysis (Figure 6) highlight the differences according to the structure of the assemblages. While the quantitative diversity (abundance and number of taxa) was significantly higher in P. crispa mats, the diversity indices were similar to P. oceanica meadows. The relatively low indices for P. crispa compared to the indices on P. oceanica are likely due to the high abundance of one species (Chlidonia pyriformis, Chlioniidae), while the diversity on P. oceanica is more evenly distributed (Figure 4H; Figure 5). The erect colonies of C. pyriformis are typical for calm water conditions and are often associated with macro-algae [45]. High abundances of erect and branching colonies (e.g., C. pyriformis, Scrupocellaria sp., Crisia sp.) on P. crispa and P. oceanica shoots are likely related to water current gradients inside P. crispa mats [46], as demonstrated before for P. oceanica meadows [47,48,49]. This is further supported by an experimental study that has confirmed negative phototaxis for the larvae of some erect bryozoan species [50], since strong light gradients are present in both habitats [46]. These traits, accompanied by the higher longevity of P. crispa and P. oceanica shoots compared to the relatively short-lived P. oceanica leaves, could be the main reason for the higher similarity of the two habitats in bryozoan families’ composition (Table 1) and diversity (Figure 4G). A richer diversity on P. oceanica shoots compared to the leaves has also been reported before in the Mediterranean Sea [51]. While the leaves of P. oceanica mainly host encrusting forms, P. crispa hosted additional growth types (petraliform and creeping), further underlining the structural diversity of the bryozoan community. As demonstrated above for P. oceanica meadows, this structural diversity is likely linked to gradients in water currents and light intensity, which have recently been confirmed for P. crispa mats [46]. Further investigations are needed to identify the mechanisms of these gradients and their influence on the epiphytic community in this specific habitat.
The high amount of rare bryozoan taxa found exclusively on P. crispa during this study further emphasizes the relevance of P. crispa for the quantitative diversity and its role as an essential habitat besides classically known hotspots, such as P. oceanica meadows. Furthermore, a high diversity of filter feeders potentially supports a diverse and productive food web by transferring biomass from the water column into benthic communities [29]. Previous studies have shown that other sessile filter feeders are abundant inside P. crispa mats [34,35,36,39]. Bryozoans attract a wide range from invertebrates to fish and from incidental to specialized predators [52]. Some predators are highly specialized on bryozoans (e.g., nudibranchs [53,54] or pycnogonids [55]). We also thus expect to find high abundances and diversities of mobile benthic predators inside these mats in the future.
Some characteristic species for the P. oceanica leaf assemblages, such as the endemic Electra posidoniae (Electridae; Gautier, 1954), were not found. This absence can be explained by the depth limitations and seasonality of these species [56].

4.2. Spatio-Temporal Variability of the Bryozoan Community inside Phyllophora crispa Mats

The high diversity and abundance of bryozoans found in the pilot study in 2018 were confirmed in the comparative study in 2019 (Figure 5, Table 1). Because sampling took place within the same period (May–July) in both years, we can not evaluate seasonal changes during the year. The differences in the taxa composition (higher abundance of Watersiporidae and Candidae in 2018; higher abundance of Aetidae in 2019; Figure 5) are consistent with a previous study carried out on the eastern side of the island [34], where Watersipora sp. (Watersiporidae) was the most abundant species at depths between 25 and 35 m. Assemblages of cheilostome bryozoans (such as Candidae, Aeteidae, and Watersiporidae) have been shown to follow seasonal and depth-related variations [10]. Our study focused on the same depth range and season during both years to avoid impacts on the results. However, satellite data from both years show differences in the Mediterranean Sea surface temperature, with relatively high values for 2018 and lower values for the first half of 2019 [56]. This effect has potentially shifted the natural seasonal variations and could have caused the observed effect.
Among the sampling sites of P. crispa during the 2019 campaign, our analysis showed higher diversity but lower abundance values for the northern sites (PC1 and PC2, Figure 4A–D). The high abundance accompanied by low diversity indices at Site PC3 was mainly driven by an exceptionally high density of C. pyriformis. The larvae of cheilostomate bryozoans (such as C. pyriformis) are often selective for suitable settling grounds [57]. The observed differences at Site PC3 could result from alterations in the environmental conditions due to the more exposed location towards the prevalent southern currents in the area [58], and thus an enhanced larval supply of this generalistic species. It remains unclear to what extent these differences could be explained by the relatively sheltered western site (Site Mix) from prevalent southern currents [58], resulting in different hydrodynamic patterns inside the bay. In addition, Campese bay is known for extensive touristic usage during the summer months, which might result in changes in the water quality. Physical disturbances and changes in water quality have been reported to affect bryozoan abundance and diversity [22,59,60]. Recent studies highlight the effects of local temperature and salinity changes in intertidal areas [61,62]. In the context of reoccurring temperature anomalies [63,64,65] and impacts on the thermohaline circulation [63], these effects are potentially also becoming relevant for sublittoral coastal habitats in the Mediterranean Sea. Further research is needed to describe hydrodynamic patterns and potential disturbances on these algal mats inside the bay.

5. Conclusions

We conclude that P. crispa mats provide an essential habitat for bryozoan diversity, harboring a high density of bryozoan colonies and a large number of families that were not present in the P. oceanica reference habitat. The Mediterranean “bryodiversity” (bryozoan diversity) has been recently estimated to 556 species, representing 9.6% of global bryozoan diversity [17]. About 79% of the bryozoan species in the Mediterranean Sea colonize coralligenous and dark and semi-dark cave habitats (219 and 220 species, respectively). Nevertheless, as recently highlighted by [66] for the mesophotic reefs in the Adriatic Sea, the understanding of bryozoans’ diversity and ecological roles in Mediterranean habitats is still far from being thoroughly investigated. Our results further strengthen the significance of P. crispa mats as a habitat harboring an exceptional bryozoan diversity, along with previous studies on epiphytic epifauna [34,35,36,67]. Regional human impacts and climate change threaten biodiversity in the Mediterranean Sea [68,69,70]. Therefore, identifying, protecting, and enhancing highly diverse habitats have become central parts of conservation strategies [71]. We suggest further investigations on the distribution of P. crispa mats along the Mediterranean coastline to confirm previous results on a larger scale and gain knowledge on the distribution of significant P. crispa aggregations. Furthermore, it is essential to understand how P. crispa mats are threatened by local and regional environmental impacts.

Author Contributions

Conceptualization, F.I.R., J.P., N.S. and C.W.; methodology, F.I.R., J.P. and N.S.; validation, F.I.R., E.C. and C.W.; formal analysis, F.I.R.; investigation, F.I.R., J.P. and N.S.; resources, C.W.; data curation, F.I.R.; writing—original draft preparation, F.I.R.; writing—review and editing, F.I.R., E.C., J.P., N.S. and C.W.; visualization, F.I.R.; supervision, C.W.; project administration, C.W.; funding acquisition, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by baseline funding of the Marine Ecology Department, University of Bremen (Bremen, Germany), and the “Institut für Marine Biologie” (IfMB, Karlsruhe, Germany). J.P. and N.S. received funding via the ERASMUS+ program.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The dataset supporting the conclusions of this study is available in the “PANGAEA” online repository: doi.org/10.1594/PANGAEA.942472.

Acknowledgments

The authors would like to thank Jenny Tuček and Mischa Schwarzmeier (IfMB), as well as Reiner and Regina Krumbach (Campese Diving Center), for logistical support throughout our study. We are also thankful to Susann Roßbach for providing helpful feedback on the manuscript and Anette Reh for support in sampling activities.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Diver assessing Phyllophora crispa mat on site PC2. The picture shows the typical closed coverage of Phyllophora crispa across the rocky surface of all sampling sites (Picture: F. I. Rossbach).
Figure A1. Diver assessing Phyllophora crispa mat on site PC2. The picture shows the typical closed coverage of Phyllophora crispa across the rocky surface of all sampling sites (Picture: F. I. Rossbach).
Diversity 14 00346 g0a1
Table A1. Results of permutational multivariate analysis of variance (PERMANOVA) of bryozoan communities among P. crispa sites, and pairwise comparison (Tukey’s honestly significant difference (HSD) test) of sites and (sub-) habitats. Significant results (p < 0.05) are indicated in bold.
Table A1. Results of permutational multivariate analysis of variance (PERMANOVA) of bryozoan communities among P. crispa sites, and pairwise comparison (Tukey’s honestly significant difference (HSD) test) of sites and (sub-) habitats. Significant results (p < 0.05) are indicated in bold.
PERMANOVA broyzoan assemblages P. crispa all sites
SourceDfSSR2Fp
Site30.94410.400072.66740.018
Residual121.41580.59993
Total152.35991
PERMANOVA broyzoan assemblages P. crispa, northern sites
SourceDfSSR2Fp
Site20.353930.24971.49760.192
Residual91.063480.7503
Total111.417411
Pairwise comparison all sites
pairs DfSSFR2pp adj
SiteMixSitePC310.2588422.3915550.2849950.0970.582
SiteMixSitePC210.2954352.62520.3043640.1090.654
SiteMixSitePC110.1277571.118430.1571180.3591
SitePC3SitePC210.6856685.6326920.4842120.0290.174
SitePC3SitePC110.4128043.3446770.3579230.090.54
SitePC2SitePC110.1077010.8432120.1232190.4441
Pairwise comparison (sub-) habitats, northern sites
pairs DfSSFR2pp adj
P. oceanica shootP. oceanica leaf13.1058114.896580.2870440.0010.003
P. oceanica shootP. crispa mat11.6540079.5246270.240980.0010.003
P. oceanica leafP. crispa mat11.92349810.447260.2648410.0010.003
Table A2. Literature used for species identification.
Table A2. Literature used for species identification.
AuthorsYearTitle
Ryland, J. S. & Hayward, P. J.1977British Anascan Bryozoans
Hayward, P. J. & Ryland, J. S.1979British Ascophoran Bryozoans
Hayward, P. J. & Ryland, J. S.1985Cyclostome Bryozoans
Hayward, P. J. 1985Ctenostome Bryozoans
Zabala, M. & Maluquer, P.1988Treballs del museu de zoologia–illustrated keys for the classification of Mediterranean Bryozoa
Hayward, P. J. & Ryland, J. S.1995Handbook of the Marine Fauna of North-West Europe
Hayward, P. J. & Ryland, J. S.1998Cheilostomatous Bryozoa: Part 1 Aeteoidea-Cribrilinoidea
Hayward, P. J. & Ryland, J. S.1999Cheilostomatous Bryozoa: Part 2 Hippothooidae - Celleporoidae
Bedini, R.2003Gli animali delle praterie a Poseidonia oceanica: dai macroinvertebrati ai pesci
Formula (A1): Calculation of bryozoan colonies on P. crispa per m2 seafloor ( C o l S F ) from colonies per m2 substrate ( C o l S S ), using wet weights of the main sample ( W W M S ) and subsample ( W W S S ), and surface area of the subsample ( S A S S ) (0.09 m2 corresponds to the size of the sampling frame):
C o l S F = W W M S × S A S S W W S S S A S S × C o l S S × 1   m 2 0.09   m 2
Formula (A2): Calculation of bryozoan colonies on P. oceanica leaves ( C o l S F ) per m2 seafloor using the mean leaf surface area ( S A l e a f A V G ), the mean number of leaves per m2 (162), the surface area of investigated leaf sample ( S A l e a f S S ), and colony count per m2 substrate ( C o l S S ):
C o l S F = 162 × S A l e a f A V G S A l e a f S S × C o l S S
Formula (A3): Calculation of bryozoan colonies on P. oceanica shoots per m2 seafloor ( C o l S F ) using the shoot surface area ( S A s h o o t S S ), the mean number of shoots per m2 (40.5), and colony count per m2 substrate ( C o l S S ):
C o l S F = 40.5 × S A s h o o t S S S A s h o o t S S × C o l S S

References

  1. Gaston, K.J. Global patterns in biodiversity. Nature 2000, 405, 220–227. [Google Scholar] [CrossRef] [PubMed]
  2. Medail, F.; Quezel, P. Biodiversity Hotspots in the Mediterranean Basin: Setting Global Conservation Priorities. Conserv. Biol. 1999, 13, 1510–1513. [Google Scholar] [CrossRef]
  3. Boudouresque, C.F.; Bernard, G.; Bonhomme, P.; Charbonnel, E.; Diviacco, G.; Meinesz, A.; Pergent, G.; Pergent-Martini, C.; Ruitton, S.; Tunesi, L. Préservation et Conservation des Herbiers à Posidonia Oceanica; Ramoge: Marseille, France, 2006; ISBN 2905540303. [Google Scholar]
  4. Mazzella, L.; Buia, M.C.; Gambi, M.M.C.; Lorenti, M.; Russo, G.F.; Scipione, M.B.; Zupo, V. Plant-animal trophic relationships in the Posidonia oceanica ecosystem of the Mediterranean Sea: A review. Plant-Anim. Interact. Mar. Benthos 1992, 46, 165–187. [Google Scholar] [CrossRef]
  5. Ballesteros, E. Mediterranean coralligenous assemblages: A synthesis of present knowledge. In Oceanography and Marine Biology: An Annual Review; Gibson, R.N., Atkinson, R.J.A., Gordon, J.D.M., Eds.; Taylor & Francis: London, UK, 2006; pp. 123–195. [Google Scholar]
  6. Ingrosso, G.; Abbiati, M.; Badalamenti, F.; Bavestrello, G.; Belmonte, G.; Cannas, R.; Benedetti-Cecchi, L.; Bertolino, M.; Bevilacqua, S.; Bianchi, C.N.; et al. Mediterranean Bioconstructions Along the Italian Coast, 1st ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2018; Volume 79, ISBN 9780128151013. [Google Scholar]
  7. Borum, J.; Duarte, C.M.; Krause-Jensen, D.; Greve, T.M. (Eds.) European Seagrasses: An Introduction to Monitoring and Management; The M&MS Project: New York, NY, USA, 2004; ISBN 8789143213. [Google Scholar]
  8. Donnarumma, L.; Lombardi, C.; Cocito, S.; Gambi, M.C. Settlement pattern of Posidonia oceanica epibionts along a gradient of ocean acidification: An approach with mimics. Mediterr. Mar. Sci. 2014, 15, 498–509. [Google Scholar] [CrossRef] [Green Version]
  9. Cocito, S. Bioconstruction and biodiversity: Their mutual influence. Sci. Mar. 2004, 68, 137–144. [Google Scholar] [CrossRef] [Green Version]
  10. Lepoint, G.; Balancier, B.; Gobert, S. Seasonal and depth-related biodiversity of leaf epiphytic Cheilostome Bryozoa in a Mediterranean Posidonia oceanica meadow. Cah. Biol. Mar. 2014, 55, 57–67. [Google Scholar]
  11. Lepoint, G.; Gobert, S.; Bouquegneau, J.M.; Havelange, S.; Gobert, S.; Bouquegneau, J.M. Fauna vs. flora contribution to the leaf epiphytes biomass in a Posidonia oceanica seagrass bed (Revellata Bay, Corsica). Hydrobiologia 1999, 394, 63–67. [Google Scholar] [CrossRef]
  12. Pardi, G.; Piazzi, L.; Balata, D.; Papi, I.; Cinelli, F.; Benedetti-Cecchi, L. Spatial variability of Posidonia oceanica (L.) Delile epiphytes around the mainland and the islands of Sicily (Mediterranean Sea). Mar. Ecol. 2006, 27, 397–403. [Google Scholar] [CrossRef]
  13. Balata, D.; Nesti, U.; Piazzi, L. Patterns of spatial variability of seagrass epiphytes in the north-west Mediterranean Sea Patterns of spatial variability of seagrass epiphytes in the north-west Mediterranean Sea. Mar. Biol. 2007, 151, 2025–2035. [Google Scholar] [CrossRef]
  14. Nesti, U.; Piazzi, L.; Balata, D. Variability in the structure of epiphytic assemblages of the Mediterranean seagrass Posidonia oceanica in relation to depth. Mar. Ecol. 2009, 30, 276–287. [Google Scholar] [CrossRef]
  15. Gluhak, T.; Lewis, J.E.; Popijac, A. Bryozoan fauna of Green Island, Taiwan: First indications of biodiversity. Zool. Stud. 2007, 46, 397–426. [Google Scholar]
  16. Rosso, A.; Gerovasileiou, V.; Sanfilippo, R.; Guido, A. Bryozoan assemblages from two submarine caves in the Aegean Sea (Eastern Mediterranean). Mar. Biodivers. 2019, 49, 707–726. [Google Scholar] [CrossRef]
  17. Rosso, A.; Di Martino, E. Bryozoan diversity in the Mediterranean Sea: An update. Mediterr. Mar. Sci. 2016, 17, 567–607. [Google Scholar] [CrossRef]
  18. Harmelin, J.G. Bryozoan facies in the coralligenous community: Two assemblages with contrasting features at Port-Cros Archipelago (Port-Cros National Park, France, Mediterranean). Sci. Rep. Port-Cros. Natl. Park 2017, 31, 105–123. [Google Scholar]
  19. Buchsbaum, R.; Buchsbaum, M.; Pearse, M.; Pearse, V. Animals Without Backbones, 3rd ed.; University of Chicago Press: Chicago, IL, USA, 1987. [Google Scholar]
  20. McKinney, F.; Jackson, J. Bryozoan Evolution; University of Chicago Press: Chicago, IL, USA, 1989. [Google Scholar]
  21. Peterson, B.J.; Frankovich, T.A.; Zieman, J.C. Response of seagrass epiphyte loading to field manipulations of fertilization, gastropod grazing and leaf turnover rates. J. Exp. Mar. Biol. Ecol. 2007, 349, 61–72. [Google Scholar] [CrossRef]
  22. Casoli, E.; Nicoletti, L.; Mastrantonio, G.; Jona-Lasinio, G.; Belluscio, A.; Ardizzone, G.D. Scuba diving damage on coralligenous builders: Bryozoan species as an indicator of stress. Ecol. Indic. 2017, 74, 441–450. [Google Scholar] [CrossRef]
  23. Reverter-Gil, O.; Souto, J. Watersiporidae (Bryozoa) in Iberian waters: An update on alien and native species. Mar. Biodivers. 2019, 49, 2735–2752. [Google Scholar] [CrossRef] [Green Version]
  24. Lombardi, C.; Taylor, P.D.; Cocito, S. Bryozoan Constructions in a Changing Mediterranean Sea. In The Mediterranean Sea: Its History and Present Challenges; Goffredo, S., Dubinsky, Z., Eds.; Springer Science and Business Media, LLC: Dordrecht, The Netherlands, 2014; pp. 373–384. ISBN 9789400767041. [Google Scholar]
  25. Lombardi, C.; Gambi, M.C.; Vasapollo, C.; Taylor, P.; Cocito, S. Skeletal alterations and polymorphism in a Mediterranean bryozoan at natural CO2 vents. Zoomorphology 2011, 130, 135–145. [Google Scholar] [CrossRef]
  26. Lombardi, C.; Cocito, S.; Gambi, M.; Cisterna, B.; Flach, F.; Taylor, P.; Keltie, K.; Freer, A.; Cusack, M. Effects of ocean acidification on growth, organic tissue and protein profile of the Mediterranean bryozoan Myriapora truncata. Aquat. Biol. 2011, 13, 251–262. [Google Scholar] [CrossRef] [Green Version]
  27. Hageman, S.J.; Bone, Y.; Mcgowran, B.; James, N.P. Bryozoan colonial growth-forms as paleoenvironmental indicators: Evaluation of methodology. Palaios 1997, 12, 405–419. [Google Scholar] [CrossRef]
  28. Cigliano, M.; Cocito, S.; Gambi, M.C. Epibiosis of Calpensia nobilis (Esper) (Bryozoa: Cheilostomida) on Posidonia oceanica (L.) Delile rhizomes: Effects on borer colonization and morpho-chronological features of the plant. Aquat. Bot. 2007, 86, 30–36. [Google Scholar] [CrossRef]
  29. Gili, J.M.; Coma, R. Benthic suspension feeders: Their paramount role in littoral marine food webs. Trends Ecol. Evol. 1998, 13, 316–321. [Google Scholar] [CrossRef]
  30. Belloni, B.; Sartoretto, S.; Cresson, P.; Bouchoucha, M.; Guillou, G.; Lebreton, B.; Ruitton, S.; Harmelin-Vivien, M. Food Web Structure of a Mediterranean Coralligenous System. In Proceedings of the 3rd Mediterranean Symposium on the Conservation of Coralligenous & Other Calcareous Bio-Concretions, Antalya, Turkey, 15–16 January 2019; p. 30. [Google Scholar]
  31. Kostylev, E.F.; Tkachenko, F.P.; Tretiak, I.P. Establishment of “Zernov’s Phyllophora field” marine reserve: Protection and restoration of a unique ecosystem. Ocean Coast. Manag. 2010, 53, 203–208. [Google Scholar] [CrossRef]
  32. Navone, A.; Bianchi, C.N.; Orru, P.; Ulzega, A. Saggio di cartografia geomorfologica e bionomica nel parco marino di Tavolara-Capo Coda di Cavallo (Sardegna nord-orientale). Oebalia 1992, XVII, 469–478. [Google Scholar]
  33. Bianchi, C.N.; Morri, C.; Navone, A. I popolamenti delle scogliere rocciose sommerse dell’Area Marina Protetta di Tavolara Punta Coda Cavallo (Sardegna nord-orientale). Sci. Rep. Port-Cros Natl. Park 2010, 24, 39–85. [Google Scholar]
  34. Bonifazi, A.; Ventura, D.; Gravina, M.F.; Lasinio, G.J.; Belluscio, A.; Ardizzone, G.D. Unusual algal turfs associated with the rhodophyta Phyllophora crispa: Benthic assemblages along a depth gradient in the Central Mediterranean Sea. Estuar. Coast. Shelf Sci. 2017, 185, 77–93. [Google Scholar] [CrossRef]
  35. Casoli, E.; Bonifazi, A.; Ardizzone, G.; Gravina, M.F. How algae influence sessile marine organisms: The tube worms case of study. Estuar. Coast. Shelf Sci. 2016, 178, 12–20. [Google Scholar] [CrossRef]
  36. Rossbach, F.I.; Casoli, E.; Beck, M.; Wild, C. Mediterranean Red Macro Algae Mats as Habitat for High Abundances of Serpulid Polychaetes. Diversity 2021, 40, 265. [Google Scholar] [CrossRef]
  37. Kikuchi, T. Handbook of Seagrass Biology: An Ecosystem Perspective; Phillips, R.C., McRoy, C.P., Eds.; Garland STPM Press: New York, NY, USA, 1980. [Google Scholar]
  38. Kikuchi, T.; Pérès, J.M. Animal communities in seagrass beds: A review. In Seagrass Ecosystems: A Scientific Perspective; McRoy, C.P., Helfferich, C., Eds.; Marcel Dekker: New York, NY, USA, 1967; pp. 147–193. [Google Scholar]
  39. Rossbach, F.I.; Merk, B.; Wild, C. High Diversity and Abundance of Foraminifera Associated with Mediterranean Benthic Red Algae Mats. Diversity 2021, 14, 21. [Google Scholar] [CrossRef]
  40. Çinar, M.E.; Féral, J.P.; Arvanitidis, C.; David, R.; Taşkin, E.; Sini, M.; Dailianis, T.; Doğan, A.; Gerovasileiou, V.; Evcen, A.; et al. Coralligenous assemblages along their geographical distribution: Testing of concepts and implications for management. Aquat. Conserv. Mar. Freshw. Ecosyst. 2020, 30, 1578–1594. [Google Scholar] [CrossRef]
  41. Shannon, C.E. A mathematical theory of communication. Bell Syst. Technol. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  42. Pielou, E.C. The measurement of diversity in different types of biological collections. J. Theor. Biol. 1966, 13, 131–144. [Google Scholar] [CrossRef]
  43. 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]
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria; Available online: http://www.r-project.org (accessed on 24 June 2021).
  45. Pisano, E.; Boyer, M. Development pattern of an infralittoral bryozoan community in the western Mediterranean Sea. Mar. Ecol. Prog. Ser. 1985, 27, 195–202. [Google Scholar] [CrossRef]
  46. Schmidt, N.; El-khaled, Y.C.; Rossbach, F.I.; Wild, C. Fleshy red algae mats influence their environment in the Mediterranean Sea. Front. Mar. Sci. 2021, 8, 721626. [Google Scholar] [CrossRef]
  47. Terrados, J.; Duarte, C.M. Experimental evidence of reduced particle resuspension within a seagrass (Posidonia oceanica L.) meadow. J. Exp. Mar. Biol. Ecol. 2000, 243, 45–53. [Google Scholar] [CrossRef]
  48. Gacia, E.; Duarte, C.M. Sediment retention by a Mediterranean Posidonia oceanica meadow: The balance between deposition and resuspension. Estuar. Coast. Shelf Sci. 2001, 52, 505–514. [Google Scholar] [CrossRef]
  49. Gacia, E.; Granata, T.C.; Duarte, C.M. An approach to measurement of particle flux and sediment retention within seagrass (Posidonia oceanica) meadows. Aquat. Bot. 1999, 65, 255–268. [Google Scholar] [CrossRef]
  50. Wendt, D.E.; Woollacott, R.M. Ontogenies of Phototactic Behavior and Metamorphic Competence in Larvae of Three Species of Bugula (Bryozoa). Invertebr. Biol. 1999, 118, 75–84. [Google Scholar] [CrossRef] [Green Version]
  51. Kocak, F.; Balduzzi, A.; Benli, H.A. Epiphytic bryozoan community of Posidonia oceanica (L.) Delile meadow in the northern Cvorus (Eastern Mediterranean). Indian J. Mar. Sci. 2002, 31, 235–238. [Google Scholar]
  52. Mckinney, F.K.; Taylor, P.D.; Lidgard, S. Predation on Bryozoans and its Reflection in the Fossil Record. In Predator-Prey Interactions in the Fossil Record; Kelley, P., Kowalewski, M., Hansen, T.A., Eds.; Springer: New York, NY, USA, 2003; pp. 239–261. ISBN 978-1-4615-0161-9. [Google Scholar]
  53. Harvell, C.D. Why Nudibranchs are Partial Predators: Intracolonial Variation in Bryozoan Palatability. Ecology 1984, 65, 716–724. [Google Scholar] [CrossRef]
  54. Harvell, C.D. Predator-induced defense in a marine bryozoan. Science 1984, 224, 1357–1359. [Google Scholar] [CrossRef] [PubMed]
  55. Dietz, L.; Dömel, J.S.; Leese, F.; Lehmann, T.; Melzer, R.R. Feeding ecology in sea spiders (Arthropoda: Pycnogonida): What do we know? Front. Zool. 2018, 15, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Pisano, A.; Fanelli, C.; Nardelli, B.B. Mediterranean Sea Anomaly Time Series of Sea Surface Temperature. Available online: https://marine.copernicus.eu/de/node/6323 (accessed on 24 June 2021).
  57. Harmelin, J.G. Diversity of bryozoans in a Mediterranean sublittoral cave with bathyal-like conditions: Role of dispersal processes and local factors. Mar. Ecol. Prog. Ser. 1997, 153, 139–152. [Google Scholar] [CrossRef]
  58. Boero, F.; De Leo, F.; Fraschetti, S.; Ingrosso, G. The Cells of Ecosystem Functioning: Towards a Holistic Vision of Marine Space, 1st ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2019; Volume 82. [Google Scholar]
  59. Harmelin, J.G.; Capo, S. Effects of sewage on bryozoan diversity in Mediterranean rocky bottoms. In Proceedings of the Bryozoan Studies 2001: Proceedings of the 12th International Bryozoology Association Conference, Dublin, Ireland, 16–21 July 2002; Swets & Zeitlinger: Sassenheim, The Netherlands; pp. 151–158. [Google Scholar]
  60. Piazzi, L.; Bonaviri, C.; Castelli, A.; Ceccherelli, G.; Costa, G.; Curini-Galletti, M.; Langeneck, J.; Manconi, R.; Montefalcone, M.; Pipitone, C.; et al. Biodiversity in canopy-forming algae: Structure and spatial variability of the Mediterranean Cystoseira assemblages. Estuar. Coast. Shelf Sci. 2018, 207, 132–141. [Google Scholar] [CrossRef]
  61. Evseeva, O.Y.; Ishkulova, T.G.; Dvoretsky, A.G. Environmental Drivers of an Intertidal Bryozoan Community in the Barents Sea: A Case Study. Animals 2022, 12, 552. [Google Scholar] [CrossRef]
  62. Bračun, S.; Wagner, M.; Koblmüller, S. Spatio-temporal occurrence patterns of epibiota along the leaves of the seagrass Cymodocea nodosa in the Northern Adriatic Sea. Mar. Biol. Res. 2021, 17, 592–602. [Google Scholar] [CrossRef]
  63. Lejeusne, C.; Chevaldonné, P.; Pergent-Martini, C.; Boudouresque, C.F.; Pérez, T. Climate change effects on a miniature ocean: The highly diverse, highly impacted Mediterranean Sea. Trends Ecol. Evol. 2010, 25, 250–260. [Google Scholar] [CrossRef]
  64. Marbà, N.; Jordà, G.; Agustí, S.; Girard, C.; Duarte, C.M. Footprints of climate change on Mediterranean Sea biota. Front. Mar. Sci. 2015, 2, 56. [Google Scholar] [CrossRef]
  65. Rivetti, I.; Fraschetti, S.; Lionello, P.; Zambianchi, E.; Boero, F. Global warming and mass mortalities of benthic invertebrates in the Mediterranean Sea. PLoS ONE 2014, 9, e115655. [Google Scholar] [CrossRef] [Green Version]
  66. Giampaoletti, J.; Cardone, F.; Corriero, G.; Gravina, M.F.; Nicoletti, L. Sharing and Distinction in Biodiversity and Ecological Role of Bryozoans in Mediterranean Mesophotic Bioconstructions. Front. Mar. Sci. 2020, 7, 581292. [Google Scholar] [CrossRef]
  67. Casoli, E.; Bonifazi, A.; Giandomanico, A.; Gravina, M.F.; Russo, G.F.; Sandulli, R.; Donnarumma, L. Comparative Analysis of Mollusc Assemblages from Different Hard Bottom Habitats in the Central Tyrrhenian Sea. Diversity 2019, 11, 74. [Google Scholar] [CrossRef] [Green Version]
  68. Buonocore, E.; Donnarumma, L.; Appolloni, L.; Miccio, A.; Russo, G.F.; Franzese, P.P. Marine natural capital and ecosystem services: An environmental accounting model. Ecol. Modell. 2020, 424, 109029. [Google Scholar] [CrossRef]
  69. Verdura, J.; Linares, C.; Ballesteros, E.; Coma, R.; Uriz, M.J.; Bensoussan, N.; Cebrian, E. Biodiversity loss in a Mediterranean ecosystem due to an extreme warming event unveils the role of an engineering gorgonian species. Sci. Rep. 2019, 9, 5911. [Google Scholar] [CrossRef]
  70. Mazaris, A.D.; Kallimanis, A.; Gissi, E.; Pipitone, C.; Danovaro, R.; Claudet, J.; Rilov, G.; Badalamenti, F.; Stelzenmüller, V.; Thiault, L.; et al. Threats to marine biodiversity in European protected areas. Sci. Total Environ. 2019, 677, 418–426. [Google Scholar] [CrossRef]
  71. Miu, I.V.; Rozylowicz, L.; Popescu, V.D.; Anastasiu, P. Identification of areas of very high biodiversity value to achieve the EU biodiversity strategy for 2030 key commitments. PeerJ 2020, 8, e10067. [Google Scholar] [CrossRef]
Figure 1. Edge of a Phyllophora crispa mat (A) with details on epiphytic fauna (B), including bryozoans, serpulids, and foraminiferans. Overview of Posidonia oceanica meadow (C), with details on the leaves (D) and shoots (E). Pictures: E.C. (A); F.R. (B,C,E); N.S. (D)).
Figure 1. Edge of a Phyllophora crispa mat (A) with details on epiphytic fauna (B), including bryozoans, serpulids, and foraminiferans. Overview of Posidonia oceanica meadow (C), with details on the leaves (D) and shoots (E). Pictures: E.C. (A); F.R. (B,C,E); N.S. (D)).
Diversity 14 00346 g001
Figure 2. Locations of the sampling sites in the study area (Isola del Giglio, Italy). Red dots mark sampling sites of Phyllophora crispa in 2019, the yellow dot marks the site of comparative sampling in both years (2018 and 2019), and Posidonia oceanica reference habitat (2019).
Figure 2. Locations of the sampling sites in the study area (Isola del Giglio, Italy). Red dots mark sampling sites of Phyllophora crispa in 2019, the yellow dot marks the site of comparative sampling in both years (2018 and 2019), and Posidonia oceanica reference habitat (2019).
Diversity 14 00346 g002
Figure 3. Number of families found per habitat during the second sampling season (2019), showing shared families between the respective habitats.
Figure 3. Number of families found per habitat during the second sampling season (2019), showing shared families between the respective habitats.
Diversity 14 00346 g003
Figure 4. Comparison of abundance and diversity descriptors among sites and habitats. First row (AD): all four investigated Phyllophora crispa sites in 2019. Second row (EH): all three (sub-) habitats on the northern sites (Site Mix, site PC1, site PC2) in 2019. Columns show the density of colonies per m2 seafloor, number of identified taxa, Shannon index, and Pielou’s index. Statistical differences (p < 0.05; Wilcoxon–Mann–Whitney test) are indicated by small letters (a–c) where significant results occurred. Black dots resemble outliers.
Figure 4. Comparison of abundance and diversity descriptors among sites and habitats. First row (AD): all four investigated Phyllophora crispa sites in 2019. Second row (EH): all three (sub-) habitats on the northern sites (Site Mix, site PC1, site PC2) in 2019. Columns show the density of colonies per m2 seafloor, number of identified taxa, Shannon index, and Pielou’s index. Statistical differences (p < 0.05; Wilcoxon–Mann–Whitney test) are indicated by small letters (a–c) where significant results occurred. Black dots resemble outliers.
Diversity 14 00346 g004
Figure 5. Visualization of Spearman ranked correlation of the abundances (colonies per m2 seafloor) of bryozoan families. Samples of Posidonia oceanica leaves and shoots are indicated in green and grey, respectively. Phyllophora crispa samples are indicated in light red (2018) and dark red (2019). Bryozoa families highlighted in green were exclusive to Posidonia oceanica, highlighted in red were exclusive to Phyllophora crispa.
Figure 5. Visualization of Spearman ranked correlation of the abundances (colonies per m2 seafloor) of bryozoan families. Samples of Posidonia oceanica leaves and shoots are indicated in green and grey, respectively. Phyllophora crispa samples are indicated in light red (2018) and dark red (2019). Bryozoa families highlighted in green were exclusive to Posidonia oceanica, highlighted in red were exclusive to Phyllophora crispa.
Diversity 14 00346 g005
Figure 6. Non-metric multidimensional scaling (NMDS) plot comparing the bryozoan communities found on the three (sub-) habitats across the northern sites (Site Mix, site PC1, site PC2) during the second sampling season (2019).
Figure 6. Non-metric multidimensional scaling (NMDS) plot comparing the bryozoan communities found on the three (sub-) habitats across the northern sites (Site Mix, site PC1, site PC2) during the second sampling season (2019).
Diversity 14 00346 g006
Table 1. Most abundant families of every (sub-) habitat and both years of Phyllophora crispa sampling (top 5 indicated in bold for each habitat) and Posidonia oceanica sampling in 2019. Numbers show the mean number of colonies per m2 seafloor ± standard deviation (SD).
Table 1. Most abundant families of every (sub-) habitat and both years of Phyllophora crispa sampling (top 5 indicated in bold for each habitat) and Posidonia oceanica sampling in 2019. Numbers show the mean number of colonies per m2 seafloor ± standard deviation (SD).
P. crispa 2018P. crispa 2019P. oceanica Leaves 2019P. oceanica Shoots 2019
FamilyMeanSDMeanSDMeanSDMeanSD
Aetidae15,9765887112,14013,31252,53512,36045,0645863
Candidae268,36893,217518323430037,9165463
Chlidoniidae1,287,926276,7381,250,665445,41900131,95446,711
Crisiidae192,08922,353347,94854,64448,94514,28864,7439896
Haplopomidae4674165256,64020,818155,14936,10600
Tubuliporidae240,49533,541200,12731,55163,96619,21033,6567157
Watersiporidae212,60837,34224,78853330000
Unknown113,59350,098130,74222,592183,53636,54436281166
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rossbach, F.I.; Casoli, E.; Plewka, J.; Schmidt, N.; Wild, C. New Insights into a Mediterranean Sea Benthic Habitat: High Diversity of Epiphytic Bryozoan Assemblages on Phyllophora crispa (Rhodophyta) Mats. Diversity 2022, 14, 346. https://doi.org/10.3390/d14050346

AMA Style

Rossbach FI, Casoli E, Plewka J, Schmidt N, Wild C. New Insights into a Mediterranean Sea Benthic Habitat: High Diversity of Epiphytic Bryozoan Assemblages on Phyllophora crispa (Rhodophyta) Mats. Diversity. 2022; 14(5):346. https://doi.org/10.3390/d14050346

Chicago/Turabian Style

Rossbach, Felix Ivo, Edoardo Casoli, Julia Plewka, Neele Schmidt, and Christian Wild. 2022. "New Insights into a Mediterranean Sea Benthic Habitat: High Diversity of Epiphytic Bryozoan Assemblages on Phyllophora crispa (Rhodophyta) Mats" Diversity 14, no. 5: 346. https://doi.org/10.3390/d14050346

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