1. Introduction
Polyplacophora are mollusks typically found in midlitoral and shallow sublitoral boulder fields [
1]. They are considered grazers, possibly stabilizing the barrens created by sea urchins [
2]. Their bioerosive potential has also been established and is comparable to other important bioeroders [
3]. There is extensive literature on a variety of organisms inhabiting this zone in the eastern Mediterranean [
4,
5]. Data regarding the presence of chiton species in the eastern Mediterranean are rather scarce and no studies exist regarding the spatio-temporal variations of distribution and abundance [
6,
7].
The patterns of distribution have been shown to be substantially different in a spatial scale, ranging from some centimeters to hundreds of meters, regarding intertidal mollusks [
8]. Temporally, the dispersion pattern of these organisms can significantly change seasonally, but also over longer periods of time [
9,
10]. A number of studies have investigated the influence of a variety of abiotic and biotic factors, mainly along the Pacific and Atlantic coasts, where it was evident that more data would be essential for a better understanding of the ecological aspects of Polyplacophora [
11,
12].
Growth parameters estimates are ecologically important in that they provide insights regarding the trophic interactions of different trophic levels. Although most species of chitons are relatively long-lived and exhibit slow growth (e.g.,
Plaxiphora aurata in Argentina with a longevity of 6–7 years,
K = 0.359), these results may be biased towards studies of species with larger sizes [
13,
14,
15]. Growth, amongst other live history features, is strongly influenced by temperature [
16]. The plasticity of the growth parameters of
Chiton articulatus in response to global warming thermal events has already been observed [
17]. Other abiotic and biotic factors might also influence these parameters, such as depth, body condition and latitude, thus leading to differences between populations of the same species [
15,
18,
19,
20].
This is the first study to examine the spatio-temporal effects on dispersion and abundance and to estimate the growth parameters for three co-occurring intertidal and shallow subtidal Polyplacophora species of the Eastern Mediterranean Sea. The present study will serve as a steppingstone for further research, possibly on the density dependent effects of these organisms on midlittoral and shallow sublittoral ecosystems. This will be particularly important in the future, since these ecosystems are expected to be notably susceptible to ocean acidification.
4. Discussion
Previous studies in the eastern Mediterranean focused on the presence and species composition of Polyplacophora, rather than spatio-temporal changes in specific populations.
Rhyssoplax olivacea was characterized as the dominant chiton species in the eastern coasts of Turkey and the most common in Greece, as also shown in the present study [
32,
33]. This species was the only one to exhibit both spatial and temporal differences in population density and individual size. Region was the main factor affecting the densities and sizes for all species. The present study examined the temporal effect in two scales (season, month) and reported a different influence for each species. Monthly changes in weight may result from a series of reproductive processes taking place throughout the year, as shown by previous literature [
21,
34]. If this is the case for the species examined here, then the present results might infer some variation in reproductive strategies, although further research is necessary.
The relationship between the size of boulder and the number of individuals has been previously tested for species of the genus Ischnochiton in Australia, where no significance was found [
12]. In the present study,
Lepidopleurus cajetanus was the only species where a relationship between boulder surface and density was reported (
Table 4). However, the underlying causes as to why this species shows a preference for larger boulders is yet unclear. This species is reported to be mostly found under stones embedded in sand, the underside colors of these stones matching the yellow or whitish color of the animals [
32]. This has also been observed for a number of other species; for instance, the red colored chiton
Callochiton septemvalvis on red algae, the white
Stenosemus albus on white algae, while color morphs for
Lepidochitona cinerea and
Acanthochitona crinita were found in habitats with matching colors [
35,
36]. During field sampling, this tendency was also observed for
L. cajetanus in the present study (authors’ observations). The preference in boulder size combined with the hypothesized camouflage behavior might serve as a defensive mechanism against both predators and disturbance, since smaller boulders are more easily overturned through wave action.
Generally, chitons are characterized by a more-or-less oval body shape, which, however, can vary from broad oval to vermiform among the 1000 living species [
37]. Intra- and interspecific allometric shifts in body shape and size are sometimes related to a particular niche adaptation [
38]. Allometry in chitons is also closely related to life stage, as was shown in previous studies [
21]. Data on polyplacophoran allometry are mainly available for tropical environments, where both negative and positive allometries have been observed for different species (see [
18] for more details). For
Chiton articulatus, a latitudinal shift in weight–length has been observed, with increased b values reported with increased latitude, which might be explained by the ectothermic nature of this species [
19]. In the present study, several differences were observed regarding the allometric relationships among the examined morphometric variables for all species with the effect of region. Since the two sites are in close proximity with each other, intraspecific differences might be explained by other environmental or biotic factors, potentially resulting in some physiological variations. More importantly, the three species exhibited different allometric trends, with
Acanthochitona fascicularis showing isometry for the weight–length relationship (
Figure 6), while
L. cajetanus was the only species where positive allometry was reported (
Figure 8). Differences in allometry often result from environmental and habitat adaptations, while they might also infer changes in resource and energy allocation during different ontogenic or reproductive stages ([
18,
20,
21,
38,
39,
40]). This might indicate differences in terms of environmental influence and habitat adaptations, possibly leading to different ecological niche occupation among the three species. Future studies examining the shift in allometric trends for different life stages or for different reproductive periods will help in explaining the inter- and intraspecific differences in allometry observed here.
Mollusks on intertidal boulder fields often show aggregated patterns of dispersion [
40]. Many factors have been proposed to result in clustered patterns, some of those being stress as a result of desiccation, wave action or food availability [
41,
42,
43]. However, a study on the genus of Ischnochiton in southern Australia concluded that physical stress did not influence the dispersion of these species. Furthermore, the same study reported no clear relationships between the dispersion of chitons and other assemblages neither on nor under boulders [
12]. In the present study, the three examined species exhibited spatio-temporal aggregated patterns, except for
A. fascicularis in autumn and winter.
L. cajetanus exhibited the highest overdispersion with spatio-temporal fluctuations, whereas the overdispersion pattern of
R. olivacea seems to remain constant.
Published data regarding mostly large chitons from tropical waters hypothesize that these organisms are relatively long-lived, while exhibiting slow growth [
44,
45,
46]. For
Katharina tunicata, the growth parameters
L∞ and
K were estimated to be 114 mm and 0.217, with an observed maximum age 17 years. An examination of the growth of
Plaxiphora aurata in Argentina reported values of 53 mm and 0.359 for
L∞ and
K, respectively, while longevity was estimated at 6–7 years [
15]. For
Chiton articulatus in Mexico,
L∞ was estimated at 71.5 and 81 mm for cold and warm annual events, respectively, using Schnute models. The same study reported maximum ages of 2.6 and 1.3 years, although these results would correspond to 5.2 and 4.7 years if the von Bertalanffy equation was used [
17]. However, one study on
Acanthochitona rectrojecta, a small-sized chiton, indicated differences in the reproduction, growth and mortality of this species compared to other larger, co-occurring chitons [
9]. For
Acanthochitona crinite, lifespan was estimated at nearly 1 year [
47].
In the present study, the estimated values for
L∞ were substantially lower, while
K values appeared to be two to three times higher. Longevity ranged from 4.72 years for
R. olivacea to 3.4 years for
L. cajetanus and
A. fascicularis. These values for maximum age are similar with those reported for larger species in tropical environments. However, the present study might overestimate the longevity of the three species, due to limitations regarding the growth estimates induced by potential constant recruitment. A seasonal variation in growth rates is often reported for various species, predominately a decrease in winter and increase during spring and summer [
15,
48,
49]. Using the seasonal oscillation index and
ts, the present study reports a strong positive term in growth during spring for
L. cajetanus and during August for
R. olivacea and
A. fascicularis. Recently, it was found that
L. cajetanus might represent a complex of species occurring throughout the Mediterranean, rather than a single species [
50]. This might introduce additional variation regarding the growth parameters presented here.
This is the first study to report spatial and temporal variations in shape and size regarding three co-occurring shallow sublittoral Polyplacophora species from the Eastern Mediterranean Sea. Although overdispersion was observed for all species, the reasons leading to clustered aggregations for these organisms are yet largely unknown. The three species exhibited different allometric patterns, possibly due to differences in a combination of ecological adaptations (e.g., resistance to disturbance) and reproductive strategies (e.g., recruitment patterns), which might be related to differences in resource allocation. Finally, it is the first time that specific parameters for growth and age are reported for chitons native to the Mediterranean Sea. The ecological aspects of the Mediterranean Polyplacophora are largely understudied, leading to a gap in knowledge regarding the ecological zones, especially the shallow sublittoral, where they are most commonly found. Further research towards their interactions with other organisms, as well as their reproduction, needs to be conducted, especially regarding their ecological niches and potential competition for food or space.