Next Article in Journal
A Review of the Impacts of Roads on Wildlife in Semi-Arid Regions
Next Article in Special Issue
Species-Specific Functional Morphology of Four US Atlantic Coast Dune Grasses: Biogeographic Implications for Dune Shape and Coastal Protection
Previous Article in Journal / Special Issue
The Influence of Canopy Cover on the Ecological Function of A Key Autogenic Ecosystem Engineer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea

1
P.P. Shirshov Institute of Oceanology of the Russian Academy of Sciences, Moscow 117997, Russia
2
Invertebrate Zoology Department, St.-Petersburg State University, Saint-Petersburg 199034, Russia
*
Author to whom correspondence should be addressed.
Diversity 2019, 11(5), 80; https://doi.org/10.3390/d11050080
Submission received: 14 April 2019 / Revised: 10 May 2019 / Accepted: 15 May 2019 / Published: 17 May 2019
(This article belongs to the Special Issue Diversity of Ecosystem Engineers in the World Coasts and Oceans)

Abstract

:
Facilitation by foundation species commonly structures terrestrial and marine communities. Intraspecific variation in individual properties of these strong facilitators can affect the whole suite of the dependent taxa. Marine macroalgae often act as ecosystem engineers, providing shelter and substrate for numerous associated organisms. Epibiosis of foliose red algae, however, remains underexplored, especially in the high latitudes. Here we studied sessile macrobenthic assemblages associated with a foliose red algae Phycodrys rubens in the White Sea (66° N) shallow subtidal, and the effect of individual plant properties on their structure. The blades of P. rubens develop annually, and it is possible to tell the young (usually larger) plant parts from the old ones. We hypothesized that epibenthic community structure depends on plant part age and size. We examined epibiosis on 110 plants at two sites, and the results generally supported our hypotheses. Old plant parts were several times smaller, and had higher total cover than young parts. Sponges strongly dominated the epibiosis on old parts, and young parts were dominated by polychaetes and bryozoans. Plant part surface area negatively correlated with total cover on young parts, while on old parts the relatioship was location-specific. On young parts the relative abundance of a polychaete Circeis armoricana increased with surface area, and the proportion of sponges decreased. The patterns indicate that epibenthic community structure is linked to the demography of an ecosystem engineer.

1. Introduction

Interspecific biotic interactions are powerful drivers of community structure. Though ecological research has long revolved around competition or predation, in recent decades the interest in positive interactions revealed the crucial role of facilitation [1,2]. Facilitation of multiple dependent species by a strong facilitator (‘foundation species’, [3]) is common in both terrestrial and aquatic ecosystems, with trees and corals being prominent examples.
Space is one of the most limiting resources in marine benthic communities [4]. Marine macrophytes act as secondary substrate and provide complex habitat architecture, functioning as ecosystem engineers [5]. By hosting numerous epibenthic taxa, they increase abundances and diversity in marine ecosystems (reviewed in reference [6]). Although species composition, spatial structure and functioning of the epibenthic assemblages associated with brown and green algae have been extensively explored (reviewed in reference [7]), macrobenthic epibiosis of red algae has been largely overlooked. The research has been focused on bacterial colonisation and algal anti-fouling activity (e.g., references [8,9,10]). Epibiotic assemblages has been described mostly for red coralline algae (e.g., [11,12]), whilst the epibiosis of foliose species is underexplored [13,14] and rarely mentioned in recent reviews (see reference [15]). At the same time rhodophytes are widespread and the most diverse group of multicellular algae [16] utilizing the largest depth range (e.g., Table 4 in reference [17]), and are commonly dominant or subdominant in lower subtidal areas in polar waters [18]. While the effect of biotic interactions on community structure is supposed to weaken with latitude, the evidence for this relationship is mostly based on comparisons of tropical and temperate ecosystems [19,20]. This knowledge gap makes the research on communities from polar regions structured by facilitation especially important.
Research on variation of the effects of foundation species on the associated organisms has mostly focused on interspecific differences between facilitators: e.g., different species of macroalgae [21,22], trees [23,24] or benthic invertebrates [25,26] hosting different dependent assemblages. However, intraspecific individual variation in foundation species has also been recently recognized as an important structuring agent. Individual properties of a facilitator which affect the dependent taxa include its genetic diversity [27], variation in behavioral patterns [28], morphological traits [29] and ontogenetic variation [30,31]. Also, an individual ecosystem engineer can provide a set of functionally different microhabitats, such as tree crowns and hollows or kelp blades and stipes, which develop different associated assemblages of dependent species [29,30,32].
Here we explored sessile macrobenthic assemblages associated with a foliose red algae Phycodrys rubens in the White Sea (66° N) shallow subtidal, and the effect of individual plant properties on their structure. P. rubens is a red foliose seaweed that is frequently found all around the North Atlantic [18], and is common in the White Sea subtidal. Though individual specimens can be found in the upper subtidal kelp forests, deeper P. rubens form a so-called «red algae belt» along with several other red algae. Life history, growth and reproduction of P. rubens were studied in various parts of its geographical range [33,34,35,36,37], including a single example of research conducted specifically in the White Sea [38]. Here, P. rubens lives for up to four years. The size of an individual plant is highly variable (length varies from 10 to 250 mm, weight—from 0.005 to 12 g [38]). New blades emerge mainly in late winter and early spring. P. rubens grows mostly in spring and autumn, and its blades are largest in November. Later under the ice cover (which lasts from November to early May [39]) most blades disintegrate, leaving only midribs which turn into side branches. The following spring these branches develop young blades that are easily distinguishable from older ones throughout the year [38]. Consequently, in summer and autumn an individual plant consists of relatively large, soft and intact young blades, and (for the plants which survived at least one winter) an old part formed by branches and rigid partly disintegrated blades from previous years. We hypothesized that due to the seasonal growth and disintegration patterns of P. rubens, epibiotic assemblages may differ on young and old parts of its thallus. We also hypothesized that plant size may affect epibiosis in terms of species composition and abundance, linking the community structure to the demography of an ecosystem engineer.

2. Materials and Methods

To determine patterns of abundance of sessile organisms associated with Phycodrys rubens, SCUBA divers sampled two subtidal sites with contrasting hydrological regimes in the Velikaya Salma Strait between the Karelian shore and Velikii Island in the western part of Kandalaksha Bay (the White Sea): ‘Site K’ (near Kamenukha island, 11 m deep, 66°33.028′ N 33°9.295′ E) and ‘Site V’ (near Velikii island, 12 m deep, 66°33.437′ N 33°6.877′ E). Strong currents in Velikaya Salma Strait are determined by semi-diurnal tides [39]. Site V was located in the narrowest part of the Strait (opposite the White Sea Biological Station of the Moscow State University). Site K was 1.7 km east from Site V in a wider part of the strait closer to the open sea. As a result, the surface current velocity during a flood tide was much higher at Site V, peaking at 1.33 m·s−1 compared to 0.75 m·s−1 at Site K (measured in July 2018 at 1 and 5 m depth with a Valeport 106 current meter). Both sites had a rocky bottom with cobble and gravel. Ice cover here typically lasts from December to May. Waters of Velikaya Salma show no thermal stratification because of tidal-mixing with water temperatures not exceeding 6–7 °C until July and raising up to 17–20 °C in August, while salinity stays closely around 25‰ throughout summer and fall [39,40].
We collected all red macroalgae from 0.25 m2 square frames randomly placed on the bottom 2–3 m from each other. Three frames per site were obtained on July 6–7, 2015, and two on September 24, 2015. The frames from the same site and month were sampled simultaneously. Two additional frames per site obtained on October 1, 2016 were only used to determine the proportion of Phycodrys rubens compared to other macroalgal species (see Figure S1 and Tables S1–S3). All plants collected were identified to the species level (and individually wet weighed in September 2015 and October 2016), but only the thalli of the dominant P. rubens from 2015 were used for further examination of epibenthic fauna. Individual P. rubens stipe with its blades and rhizoids was regarded as a separate sample. In July, when the covers of epibenthic organisms on P. rubens were exceedingly low, epibiosis was examined on 5 random individual P. rubens from each frame (30 plants in total). In September 2015 (when the covers were much higher) we examined epibiosis on 15 random individual P. rubens plus 5 largest of the remaining arrays from each frame (80 plants in total). These additional largest plants were added to compensate for their disproportionately high contribution in total substrate area and to explore the relationship between the structure of epibiosis and plant size. Per frame, the 5 largest plants together provided 569 ± 157 cm2 of the total surface area and 56 ± 11 cm2 of the area covered by epibionts compared to 245 ± 61 and 32 ± 10 cm2 (correspondingly) contributed by 15 random plants (sites pooled). All the plants were preserved in 10% buffered formalin on sea water prior to examination.
The blades of P. rubens develop annually so that the current year’s young blades, which constitute the major part of the thallus, are visually distinguishable from the older part remaining from the previous 1–3 years [38]. We examined the epibiosis separately on these two parts (hereafter ‘young’ and ‘old’, see Figure 1 for the scheme). Importantly, many plants examined had no old part, especially in September samples. For each young or old plant part we documented its surface area (accurate to 1 mm2), and identified to species level (except for sponges) and counted all sessile organisms attached. We determined the area of each sponge found (accurate to 0.1 mm2), calculated the number of units (zooids/polyps) in each hydrozoan or bryozoan colony, and individually measured the opercular diameter in serpulid polychaetes (accurate to 0.05 mm2) to estimate areas covered and percent covers according to the size-area allometric relationships established from subsamples (Supplementary File 1). Undetermined bryozoan ancestrulae with individual areas less than 0.15 mm2 were excluded from further analyses.
To assess the distribution of P. rubens by surface area and roughly estimate its contribution to substrate pool for epibenthic organisms, we approximated the surface area of the plants not preserved for examination from their wet weight (in grams) based on linear relationships established from measurements of the plants examined in September 2015. Sessile epibiota was not removed prior to weighing because of logistical constraints. A single outlier extremely heavily fouled by sponges (having, in contrast to all other plants, a sponge weighing more than a plant itself) was excluded. We used separate relationships for young (Syoung, cm2) and old (Sold, cm2) parts of the plant: Syoung = 27.189638 · Weight + 5.944277, R2 = 0.68, n = 79, and Sold = 6.942554 · Weight + 0.463375 R2 = 0.72, n = 79.
For the plants with both young and old parts, we analyzed the relationship between the total cover of epibenthic organisms and plant part age (young or old), sampling month (July or September), and site (K or V) with beta-regression (see below), using individual plant ID as a random blocking factor nested in site × month interaction. Total cover of epibiosis on the plants with only a young part was analyzed with a reduced model with only month and site as predictors. To assess the relationships between multivariate community structure and plant part age, sampling month and location, we applied permutational analyses of variance (PERMANOVA, see reference [41]) with the same factor sets. We used Bray-Curtis similarities calculated on fourth root transformed percent covers as a distance measure for PERMANOVAs, and visualized multivariate differences between the assemblages associated with plant parts of different ages in September (when both young and old parts were sufficiently covered) with principal coordinates ordination (PCO). We also compared average total cover and the covers of top abundant species between young and old parts of the plants with both young and old parts using the Wilcoxon matched pairs test, and between July and September—with the Mann-Whitney U-test (separately for young and old parts). The covers of top abundant species could not be analyzed in a similar fashion to total cover, since they had zero values, and the only available software implementation of zero-inflated beta-regression mixed models (i.e., supporting random effects) did not converge. Consequently, the covers of top abundant species were only compared between young and old plant parts by using pairwise tests as described above.
To assess the effects of location and individual size variation in the ecosystem engineer on the suite of dependent species, we examined the effects of site, P. rubens surface area and their interaction on total percent cover separately for young and old parts of P. rubens collected in September 2015 with beta-regression (see below). The effects of the same predictors on Shannon-Wiener species diversity and the number of species were analyzed with type III sum of squares ANCOVAs. We visualized multivariate differences between the assemblages associated with plants of different size and location with PCO based on the Bray-Curtis similarity matrix calculated on fourth root transformed percent covers. The effects of site, plant part surface area and their interaction on multivariate community structure were tested with PERMANOVA analyses [41] separately for young and old plant parts. The relationships between plant part surface area, sampling site and covers of 5 top abundant species were also analyzed with beta-regression. Since the two sites were intentionally selected to have contrasting hydrological regimes, the corresponding effect with two levels was treated as fixed in all the analyses. We did not combine assessing the effects of plant part surface area and age in a single model since many plants in September did not have an old part, and these were, on average, much smaller than the plants having both parts, which would lead to biased estimates.
For all the proportional data (total cover and relative covers of dominant species) where we fitted beta-regression models, we used zero-inflated distributions if zero observations occurred. We addressed heteroscedasticity by modeling the variances where this improved models’ quality based on the generalized Akaike Information Criterion (AIC). In one case (relative covers of a bryozoan Electra pilosa on old blades) a zero-inflated model could not converge, since all the zero observations were concentrated at one site. Separate models for the two sites were fitted in this case. We checked distributions of residuals using a detrended quantile plot (‘worm plot’, [42]). Analyses were conducted using the ‘gamlss’ package [43] for R [44]. To assess the effect of location and plant surface area on average individual size of epibenthic organisms we used type III sum of squares ANCOVAs with site, area and their interaction as predictors. Replicate 0.25 m2 frames were pooled by site in all the analyses since (i) their number (2–3 per sample) was too small to be used as a random blocking factor in the models (at least 5–6 levels required), and (ii) plant areas, the proportion of young/old plant parts’ surface areas, and total cover of epibenthic organisms showed no significant differences between the simultaneously sampled frames from the same site (Tables S4 and S5). In addition, we supplemented all the analyses supporting our main conclusions with alternative models where we included Frame as a fixed effect (see Tables S6–S10), which had no principal effect on the results. Means are reported ± standard error (SE). All the raw data needed to reproduce the analyses are available as supplementary data (Supplementary File 1).

3. Results

P. rubens was a dominant red algae in all the samples (Table S1). Smallest plants were most abundant, while largest ones contributed most substrate surface for epibionts (Table S2, Figure S1), in total approximately 0.6–1.7 m2 per 1 m2 of the bottom in September 2015 and 2.4–3.3 m2 in September 2016 (Table S3). In terms of surface area young parts of P. rubens in September were 2–3 times larger than old ones (Table S3). Also, in September 2015, 57% of individual plants had only a young part and no old part compared to 17% in July. In July 2015 we found, in total, 24 species of sessile epibionts on 30 P. rubens individuals examined. There were 11 bryozoans, 9 hydroids, 3 bivalves, 1 ascidian, 1 spirorbid polychaete and unidentified sponges. In September 2015 there were 41 species of sessile epibionts: 20 bryozoans, 11 hydroids, 3 bivalves, 3 ascidians, 2 epiphytic red algae, 1 barnacle, 1 spirorbid polychaete and unidentified sponges on 80 P. rubens individuals examined. Total cover (sites and frames pooled) was several times higher in September (10.7 ± 0.9%) vs. July (2.6 ± 0.4%).
Effects of plant part age and month on total cover of epibenthic organisms were both significant (Table 1, Table S6). Regardless of location and season total cover was much higher on old than young parts. The difference was significantly higher in July (15–20-fold) compared to September (3–5-fold) (Table 1 and Table 2). Total cover on young plant parts in September was dramatically higher than in July (Mann-Whitney U-test, Z = −7.374, p = 1.66·10−13, sites pooled). In July, in fact, young parts were almost clear of epibionts. The covers on old parts were also higher in September than in July (Mann-Whitney U-test, Z = −2.214, p = 0.0241, sites pooled). As a result, given the difference in size between young and old plant parts, in terms of substrate area covered in July, 25 ± 6% and 75 ± 6% of P. rubens epibionts occupied young and old parts, correspondingly (n = 6 frames, sites pooled). In September the proportion was 53 ± 9% and 47 ± 9% on young and old plant parts, correspondingly (n = 4 frames, sites pooled).
In July total and average numbers of species were lower on young plant parts than on old ones at both sites. In September, however, the number of species (being generally higher than in July) was lower on young parts only at Site K. While in July the Shannon-Wiener species diversity index was higher in the epibiosis on old parts compared to young ones, in September in contrast it became higher on young parts compared to old ones regardless of the location (Table 2).
The effects of plant part age, month and location on multivariate community structure were all significant (Table 3, Table S7). Young plant parts in July were nearly empty, in September being co-dominated by bryozoans Cribrilina annulata, Celleporella hyalina and Electra pilosa, a serpulid polychaete Circeis armoricana, and sponges. Old parts, in contrast, were strongly dominated by sponges (especially in September) with C. annulata, E. pilosa, C. armoricana and juvenile mytilid bivalves as subdominants (Table 4, Figure 1). In July percent covers of sponges and C. armoricana were significantly higher on old parts regardless of the location. C. annulata, C. hyalina and juvenile mytilids either had equally low covers on old and young plant parts, or were significantly more abundant on old ones, depending on the location. E. pilosa was nearly absent (Table 4). In September sponges and juvenile mytilids were significantly more abundant on old parts than on young ones. E. pilosa displayed a similar distribution, but only at Site V where it was relatively abundant. In contrast, C. armoricana and C. hyalina were either equally abundant on young and old parts or more abundant on young parts, depending on the location (Table 4). In September C. armoricana and juvenile mytilid bivalves were subdominant on old parts at both sites, while the subdominant bryozoans were C. annulata at Site K and E. pilosa at Site V. At the same time on young parts bryozoans comprised more than a half of total cover, with C. annulata dominating at Site K, and E. pilosa and C. hyalina at Site V. Sponges and C. armoricana were equally abundant here at both sites (Figure 1). Average size of an individual or colony was either larger on old parts or was not affected by substrate age depending on the species and location (Table 4). On a PCO plot plant parts sampled in September were grouped by age and site, while the old and young parts from same plants were generally shifted from each other along the vector of the sponges abundance (Figure 2).
Since in July we did not sample largest plants, which contribute the most substrate surface, the connection between the epibiosis structure and plant surface area was studied based on September samples. Total cover of epibionts on young and old plant parts had different relationships with their surface area (Table 5 and Table S8, Figure 3). The relationship was negative for the total cover on young plant parts at both sites. The smallest plants with the total area of the young part being less than 10 cm2 had young parts that were 3.1%–38.8% (9.7 ± 1.1% on average) covered, while the total cover of young parts in the largest plants (young part >50 cm2) never exceeded 8.6%. Total cover on young parts at Site K was generally lower than at Site V. In contrast, the effect of surface area on total cover of old plant parts was site-specific, strongly positive on Site V and almost absent at Site K (Table 5 and Table S8, Figure 3). The total cover on old plant parts at Site K was higher than at Site V.
On young plant parts, the Shannon-Wiener diversity index (H’, calculated from percent covers) and the number of species both significantly increased with surface area with a slope that was different between the sites (there was a significant Site × Area interaction). On old parts only the number of species was positively affected by the plant part surface area, while H’ was not affected neither by area nor by location (Table 6, Table S9). Plant surface area and location significantly affected the multivariate community structure both on old and young plant parts (Table 7, Table S10). On the PCO plots plant parts of both age grouped by sampling site (along the primary axis of variation) and along the surface area gradient (secondary axis of variation) (Figure 4, Figure 5).
Four of the five top abundant taxa (hereafter ‘dominants’) on old and young plant parts in September 2015 were common: Circeis armoricana, Cribrilina annulata, Electra pilosa and sponges. While Celleporella hyalina was the fifth one on young parts, old blades were co-dominated by juvenile mytilids. Together these 5 top abundant taxa contributed 79.9 ± 2.7% (n = 46) of total cover on old plant parts, and 79.6 ± 1.9% (n = 80) on young parts (sites and frames pooled). Composition of these dominants strongly changed with the surface area on young parts and much less on old ones. In terms of relative abundances, either plant part surface area or its interaction with Site affected C. armoricana, C. hyalina and sponges on young parts, and only C. annulata on old parts (Table 8, Figure 6). In addition, average colony size in C. hyalina on young parts was negatively correlated with area (Table S11). C. armoricana relative cover on young parts increased with area at both sites (but faster at Site K) (Table 8). Its individual size also slightly increased with area at Site V (Table S11). Sponges, in contrast, decreased their relative abundance (slightly slower at Site K) on larger young parts at both sites (Table 8). Overall, on young parts only C. armoricana of all the dominants showed similar constant increases of area covered with increasing substrate area at both sites (Figure 6).
On old plant parts, C. annulata was the only dominant species with relative abundance significantly affected by surface area of P. rubens. However, the effect was negative at Site V, where the covers of C. annulata were very low, and almost absent at Site K, where this bryozoan was abundant. In addition, while the non-zero values of C. annulata relative abundance decreased with increasing area, the frequency of its absence also decreased on larger plants (as indicated by a significant effect in Nu-model, see Table 8), further reducing the net effect.

4. Discussion

At both locations studied, Phycodrys rubens was the dominant seaweed and added substantial substrate space to a habitat. Consistent with our hypotheses plant part age and size both correlated with species composition and abundances. Old parts were dominated by sponges with total percent cover of epibionts differently affected by surface area depending on the location. Young parts were dominated by bryozoans and spirorbid polychaetes, and had lower total percent cover which was negatively affected by surface area regardless of the location.
Macroalgae as foundation species can host a rich epibiotic community, comprised mostly of epiphytic algae and sessile invertebrates [15]. The number of species and diversity of epibiosis associated with P. rubens were similar to other epibiotic assemblages in high latitudes, both on brown [30,45,46] and red algae [13,47]. Species composition was consistent with scarce data on red algae epibiosis in the White Sea, with bryozoans being the most abundant and species-rich group [13,48]. Particularly, bryozoans Electra pilosa, Celleporella hyalina, Cribrilina annulata, and serpulid sedentary polychaetes are common in boreal and arctic epiphytic assemblages associated with kelps and red algae [13,45,47]. However, according to our findings, the epibiosis on P. rubens lacked any abundant epiphytic algae, which tend to subdominate alongside bryozoans [45] or predominate [46] in comparable systems. The number of species was similar to the only estimate available for P. rubens epibiosis [13].
In marine benthic habitats substrate space is commonly in short supply [4,49]. It is either released from occupants by disturbance events, clearing ‘Type I patches’ sensu Sousa [50], or created anew by habitat engineers like seaweeds or bivalves, which serve as substrate for epibiosis (‘Type II patches’). P. rubens blades are typical Type II patches. Species composition and abundances of the residents of such a patch can generally result from the interplay of colonization history and further survival driven by interspecific interference and predation or grazing [51]. Colonization history, in turn, is most affected by propagule supply and larval selectivity, while interspecific interactions depend on space deficiency and relative competitive ability; substrate longevity should affect both.
Our results revealed the relationships between the plant part age, total cover and species composition of Phycodrys-associated epibiotic community. Substrate age can directly result in a higher number of species and total cover on older plant parts: these have been exposed to potential colonizers for longer (from 1 to 3 years) and could accumulate the results of recruitment, which often has strong interannual variation [52,53]. This, however, can only affect species with a long enough lifespan. The colony lifespan at least of the two dominant bryozoans, E. pilosa and C. annulata, is limited mostly by substrate lifetime and reaches up to 2–3 years. Both species definitely can make it through winter [54,55] as their zooids hibernate during the cold season and emerge again with rising water temperature. Though the data on longevity of sponges in our study region is scarce, the lifespan of common species, Myxilla incrustans and Iophon piceus, in the White Sea is reported to exceed 4 years [56]. Taking into account that sponges are generally longest-living animals [57], the species inhabiting Phycodrys (presumably Mycale sp. and Halichondria sitiens) are likely to last long and successfully overcome the cold season. Mytilid bivalves found on P. rubens (Musculus discors and Modiolus modiolus) are also perennial: M. discors has a lifespan of 2–5 years [58], and M. modiolus can last for 50 years or more [59]. However, these species do not necessarily stay always attached to a same substrate: like other mytilids they can use macroalgae as a transitional habitat for juveniles [60,61]. Another dominant bryozoan, C. hyalina, in contrast, has an annual lifespan [62]. Similarly, largest individuals of a serpulid polychaete C. armoricana are found on the thalli Laminaria that renew each year [63]. Interestingly, consistent with these species-specific differences in life histories, relatively short-living C. armoricana and C. hyalina in September on at least at one site had higher abundances on young plant parts compared to old ones. The dominants with longer lifespan were either more abundant on old parts or not affected by substrate age (Table 4).
Old parts of P.rubens are closer to the bottom than young ones, potentially causing the difference in flow velocity and illumination. However, possible consequences of this difference are complicated. Decreased flow beneath macroalgae resulting in increased sedimentation affects feeding performance and lowers growth rates in bryozoans, bivalves and cirripedes, at the same time being beneficial for serpulid polychaetes and sponges [64,65]. Flow intensity also affects larval settlement of sessile invertebrates (reviewed in reference [61]). While settling bryozoans Bugula turrita and Tubulipora sp. or a serpulid polychaete Pseudochitinopoma occidentalis prefer the zones of reduced flow velocity and turbulence, a hydroid Tubularia crocea is attracted by higher velocities [66,67]. The difference in lighting conditions resulting in microalgal turf development intensity may be beneficial for sponges or ascidians [64,68], but disadvantageous for serpulids [64]. On Fucus serratus several epiphytic sponges are more abundant on relatively rigid fibrous basal segments of Fucus serratus due to selective settling or differential early survival [7,69].
Heavier fouling itself can further facilitate some species to colonize old Phycodrys parts. Extensively covered by epibionts old parts become more rigid and rugged, which may favor other epibionts. Primary epibionts commonly provide more suitable substrate for the secondary ones than a basibiont, developing ‘fouling cascades’ [70]. Intensive fouling can reduce the amount of anti-fouling chemicals released by algae [14] which also can provide positive feedback for subsequent arrivals. On the other hand, settling larvae of many invertebrate species strictly prefer younger parts of marine macroalgae [71,72,73,74,75]. This mechanism allows the vulnerable postlarval stages to occupy the most long-living and unfouled part of the basiphyte. In particular, bryozoans Membranipora sp. and Scrupocellaria reptans prefer basal growing parts of kelp [71,73], while Alcyonidium hirsutum and Flustrellidra hispida settle mostly on the frond tips—the youngest parts of Fucus serratus [72]. A spirorbid polychaete Neodexiospira brasiliensis selects the youngest regions when settling on seagrasses species [75], while Spirorbis corallinae and S. pagenstecheri show the same pattern on Laminaria digitata [71].
With several times higher total percent cover (up to 36% in September on Site V, see Table 2) and much longer existence term, the epibiosis on old plant parts must experience the pressure of interspecific interference. While competition can hardly shape community structure on young parts having total cover within 10%, on old parts the importance of overgrowth interactions is quite expectable. In tropical epibioses species abundances and composition are commonly determined by the relative competitive strength of epibionts [76,77,78]. Almost nothing is known though of the competitive hierarchies among benthic sessile invertebrates in polar regions. Yet, the strong dominance of sponges, which is the most important feature of old Phycodrys blades compared to young ones (Figure 1, Figure 2) is seemingly caused by their highly competitive strength. Indeed, based on the studies in temperate and tropical waters higher taxonomic groups have the following competitive hierarchy: ascidians (mostly colonial) and sponges are superior in relation to bryozoans, and the latter are better competitors for space than cirripedians, polychaetes and hydroids [79].
In particular, on coral reefs at Rio Bueno, Jamaica, sponges are the strongest overgrowth competitors for space which results in their ultimate dominance on older regions of coral undersurfaces, with young coral edges left for inferior competitors including cheilostome bryozoans [80]. Alongside remarkable resistance to environmental stresses (e.g., references [81,82] sponges can exhibit strong allelochemical effects on their competitors [76,83,84].
Consistent with our observations on the Phycodrys-associated community, epibiosis on macroalgae commonly varies across a single thallus age-, structure- or function-dependently [30,74,85,86,87,88,89,90]. On Lobophora variegata (a brown alga) community structure strongly depends on the blade side. Upper- and undersides show different species composition, and undersides also harbor greater numbers of individuals. The supposed advantages of an underside include (i) providing potential shelter from herbivorous grazers, (ii) hosting several species of secondary basiphytes, (iii) containing lower amounts of defensive substances, which serve as anti-fouling and anti-grazing agents [85]. Epibiosis on a common brown kelps Saccharina latissima and Laminaria digitata show the decrease in species diversity on older parts of the lamina [30]. One of possible reasons for such a decline is that the distal (and oldest) part of the lamina is mostly exposed to physical stress and tissue decay. Algal epiphytes of a kelp Ecklonia radiata are, in contrast, most abundant on oldest parts of the thalli [86]. While the authors propose the length of time spent in water as a principal determinant of these age-related differences, the position in the water column and surface rugosity are also considered important.
In addition to plant parts’ age, epibiosis on P. rubens was also related to their surface area. At both sites, total percent cover on young plant parts dropped on larger plants, while the proportion of solitary organisms, number of species and diversity increased. On old plant parts (which were several times smaller than young parts) the effect on total cover was location-specific, species diversity was constant, while the number of species also increased with substrate surface area. Similarly, in the epibiosis of Fucus serratus, plant size affects abundances and within-plant distribution of sponges and tunicates [91]. In marine benthic communities there are also other Type II biogenic habitat patches where epibiosis is affected by substrate surface area. The relationship between patch size and community structure is best studied in the systems with strong competition for substrate space. There the larger the patch, the higher the probability that one would be colonized by poorly-recruiting strong competitors (tunicates or sponges), causing space monopolization [51]. Patch size thus shifts the principal structuring process from colonization to interspecific interference. On young P. rubens parts total percent cover is much lower, and the role of competition for space unlikely is important. Consistently, the relative abundances of potentially strong competitors (modular organisms, especially sponges) did not increase, but, on the contrary, decreased on larger plants (Table 8, Figure 6). Our results suggest that the effect of substrate size on the epibiosis of relatively large and lightly fouled young parts is driven by recruitment limitations of dominant species. It looks like only a polychaete Circeis armoricana of top abundant species had a strong enough propagule supply in 2015 to maintain nearly constant percent cover despite the substrate size and sampling location (Figure 6). Other dominants underperformed either at one of the sites (bryozoans) or on both (sponges). A bryozoan Celleporella hyalina had also smaller mean colony size on larger plants. This may either reflect that substrate size affects colony growth or results from seasonal variation in recruitment: largest plants could be colonized during the heaviest settlement peak, which might coincidentally be the latest leading to smaller (younger) colonies. Settlement and survival conditions on larger plants can be less advantageous compared to smaller ones: a major part of young blades in large plants is located further from the bottom than in small ones, and consequently is more exposed to the flow. Increased flow velocity may in turn negatively affect larval settlement, colony growth and feeding in many epibiotic species (reviewed in reference [7]).
The effects of substrate size observed on relatively small and close to the bottom (compared to young ones) old Phycodrys parts are less interpretable. Here, much higher total percent cover strongly increased with the substrate surface area at Site V and slightly decreased at Site K. Cribrilina annulata was the only dominant species with relative abundance affected, and the effect was also location-specific, but inverse. While sponges constituted more than a half of average total cover on old blades, their relative abundances exhibited no relationship with substrate size. Importantly, in contrast to young plant parts, which clearly have developed within the recent season, old parts comprise the blades of variable undetermined age, possibly different for the two sites. This, together with a plethora of other factors, including the location-specific difference in flow regime, could contribute to the observed inconsistency in plant size effect between the sampling locations. The relationship between the substrate area and the number of species was positive both on young and on old plant parts. This result is rather obvious. Indeed, larger plants have a higher probability of hosting rare species with limited propagule supply, and at the same time are more complex due to a higher number of blades, while patch size and habitat complexity are both known to increase the species number [92].
The proportion of young and old parts and the size of an individual plant reflects the ontogenetic stage of P. rubens and local growth conditions. In addition to the limited data on potential community-wide effects of ontogenetic variation in an ecosystem engineer for seaweeds [90], the examples from terrestrial habitats are available. For instance, species diversity and composition of lichens epiphytic on ash trees changes significantly with tree age [93]. In tropical rainforests species richness of ants and herbivore insects increases with host tree size, which also affects species composition [29]. Similarly, the observed effects of individual P. rubens properties on the associated epibiosis show how the demography of this ecosystem engineer can generate a spatial structure for the entire community.
The patterns documented here can result from multiple interacting processes. A survey of interannual variation in community structure is seemingly needed to interpret the observations on epibiosis covering old blades, where the traces of multiple previous years accumulate. While only manipulative experiments can reveal causal relationships behind the correlations we found, understudied systems like the one described first need observations to provide the context for further experiments [94]. Phycodrys rubens grows in easily accessible shallow subtidal, where the disturbance level is low. An individual plant commonly has a small rock or gravel as a substrate, thus offering a convenient unit for field and laboratory manipulations. Epibenthic community associated with P. rubens is thus potentially a perfect system for further experiments aimed towards isolating the contributions of colonization, interspecific interference and possible grazing to the patterns observed.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-2818/11/5/80/s1, Figure S1: Size structure of Phycodrys rubens and its proportion compared to other red macroalgal species in September 2015 and 2016 by site and frame, Table S1: Average number and wet weight of foliouse red algae per 1 m2 of the bottom in September-October by species and year, Table S2: Contribution of 20 largest Phycodrys rubens plants to approximated total area in a frame and total wet weight in September-October by year, site and frame, Table S3: Approximated surface area of Phycodrys rubens per 1 m2 in September-October by plant part age, site and year, frames pooled, Table S4: Mann-Whitney U-tests comparing total cover of the epibiosis on Phycodrys rubens in September 2015 between the two replicate frames by site and plant part age, Table S5: Kolmogorov–Smirnov and Mann-Whitney U-tests comparing total plant surface areas and the proportion of young blades’ surface of Phycodrys rubens between the pairs of replicate frames in September 2015 by sampling site, Table S6: Effects of plant part age, month, and location on total cover of Phycodrys rubens epibiosis (alternative analysis with Frame factor included), Table S7: Effects of plant part age, month, and location on multivariate community structure of epibiosis on Phycodrys rubens (alternative analysis with Frame factor included), Table S8: Effects of plant part size and location on total cover on young and old plant parts of Phycodrys rubens in September 2015 (alternative analysis with Frame factor included), Table S9: Effects of plant part size and location on the number of speicies and diversity on young and old plant parts of Phycodrys rubens in September 2015 (alternative analysis with Frame factor included), Table S10: Effects of plant part size and location on multivariate community structure of epibiosis on young and old plant parts of Phycodrys rubens in September 2015 (alternative analysis with Frame factor included), Table S11: Effects of and plant part size and location on mean individual size of 5 top abundant species on young and old plant parts of Phycodrys rubens in September 2015, Supplementary File 1: Raw data.

Author Contributions

Conceptualization, A.C., A.A, and E.Y.; Methodology, A.C., A.A, and E.Y.; Software, A.C. and E.Y.; Validation, E.Y. and A.A.; Formal analysis, A.C. and E.Y.; Investigation, A.C. and A.A.; Resources, A.C., A.A. and E.Y.; Data curation, E.Y.; Writing—original draft preparation, A.C.; Writing—review and editing, A.C. and E.Y.; Visualization, A.C.; supervision, E.Y.; Project administration, E.Y.; Funding acquisition, E.Y.

Funding

This research was funded by Russian Foundation for Basic Research (RFBR), grants 14-04-00972a and 17-04-00651a.

Acknowledgments

We thank our colleagues Marina Varfolomeeva and Vadim Mokievsky for fruitful discussion and comments on the earlier versions of the manuscript. We acknowledge the support of numerous field volunteers, especially Ksenia Shunkina, Vladimir Krapivin, Nikolay Neretin, Glafira Kolbasova, Anton Makarov and Vladimir Chava. We are extremely grateful to Alexander Tsetlin who gathered us at the White Sea Biological Station and thus conceived this research, while the staff of the Station provided accommodation and technical support, and the diving team backed our underwater affairs with assistance and equipment. The Valeport current meter was kindly offered by Natalia Kosevich. Our heartfelt thanks to Jorge Gutiérrez for invitation to participate in this volume and a chance to celebrate the New Year by meeting a deadline.

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. Stachowicz, J.J. Mutualism, facilitation, and the structure of ecological communities: Positive interactions play a critical, but underappreciated, role in ecological communities by reducing physical or biotic stresses in existing habitats and by creating new habitats on which many species depend. AIBS Bull. 2001, 51, 235–246. [Google Scholar]
  2. Bruno, J.F.; Stachowicz, J.J.; Bertness, M.D. Inclusion of facilitation into ecological theory. Trends Ecol. Evol. 2003, 18, 119–125. [Google Scholar] [CrossRef]
  3. Dayton, P.K. Toward an understanding of community resilience and the potential effects of enrichments to the benthos at McMurdo Sound, Antarctica. In Proceedings of the Colloquium on Conservation Problems in Antarctica; Dayton, P., Parker, B.C., Eds.; Allen Press: Lawrence, KS, USA, 1972; pp. 81–96. [Google Scholar]
  4. Dayton, P.K. Competition, disturbance, and community organization: The provision and subsequent utilization of space in a rocky intertidal community. Ecol. Monogr. 1971, 41, 351–389. [Google Scholar] [CrossRef]
  5. Jones, C.G.; Lawton, J.H.; Shachak, M. Organisms as ecosystem engineers. Oikos 1994, 69, 373–386. [Google Scholar] [CrossRef]
  6. Thornber, C.S.; Jones, E.; Thomsen, M.S. Epibiont-marine macrophyte assemblages. In Marine Macrophytes As Foundation Species; Olafsson, E., Ed.; CRC Press: Boca Raton, FL, USA, 2016; pp. 43–65. [Google Scholar]
  7. Seed, R.; O’Connor, R.J. Community organization in marine algal epifaunas. Annu. Rev. Ecol. Syst. 1981, 12, 49–74. [Google Scholar] [CrossRef]
  8. Maximilien, R.; de Nys, R.; Holmström, C.; Gram, L.; Givskov, M.; Crass, K.; Kjelleberg, S.; Steinberg, P.D. Chemical mediation of bacterial surface colonisation by secondary metabolites from the red alga Delisea pulchra. Aquat. Microb. Ecol. 1998, 15, 233–246. [Google Scholar] [CrossRef]
  9. Nylund, G.M.; Pavia, H. Chemical versus mechanical inhibition of fouling in the red alga Dilsea carnosa. Mar. Ecol. Prog. Ser. 2005, 299, 111–121. [Google Scholar] [CrossRef]
  10. Paul, N.A.; de Nys, R.; Steinberg, P.D. Chemical defence against bacteria in the red alga Asparagopsis armata: Linking structure with function. Mar. Ecol. Prog. Ser. 2006, 306, 87–101. [Google Scholar] [CrossRef]
  11. Amado-Filho, G.M.; Maneveldt, G.; Pereira-Filho, G.H.; Manso, R.C.C.; Bahia, R. Seaweed diversity associated with a Brazilian tropical rhodolith bed. Cienc. Mar. 2010, 36, 371–391. [Google Scholar] [CrossRef] [Green Version]
  12. Peña, V.; Bárbara, I. Maërl community in the north-western Iberian Peninsula: A review of floristic studies and long-term changes. Aquat. Conserv. Mar. Freshw. Ecosyst. 2008, 18, 339–366. [Google Scholar] [CrossRef]
  13. Grishankov, A.V. Preliminary observation on the structure of macroepifaunal communities associated with the surface of Phycodrys rubens (Rhodophyta) from the White Sea. Vestnik Leningradskogo Universiteta 2000, 2, 101–104. [Google Scholar]
  14. D’Antonio, C. Epiphytes on the rocky intertidal red alga Rhodomela larix (Turner) C. Agardh: Negative effects on the host and food for herbivores? J. Exp. Mar. Biol. Ecol. 1985, 86, 197–218. [Google Scholar] [CrossRef]
  15. Olafsson, E. Marine Macrophytes as Foundation Species; CRC Press: Boca Raton, FL, USA, 2016; pp. 43–65. [Google Scholar]
  16. Searles, R.B. The strategy of the red algal life history. Am. Nat. 1980, 115, 113–120. [Google Scholar] [CrossRef]
  17. Markager, S.; Sand-Jensen, K. Light requirements and depth zonation of marine macroalgae. Mar. Ecol.-Prog. Ser. 1992, 88, 83. [Google Scholar] [CrossRef]
  18. Wiencke, C.; Clayton, M.N.; Gómez, I.; Iken, K.; Lüder, U.H.; Amsler, C.D.; Karsten, U.; Hanelt, D.; Bischof, K.; Dunton, K. Life strategy, ecophysiology and ecology of seaweeds in polar waters. Rev. Environ. Sci. Biotechnol. 2006, 6, 95–126. [Google Scholar] [CrossRef] [Green Version]
  19. Schemske, D.W.; Mittelbach, G.G.; Cornell, H.V.; Sobel, J.M.; Roy, K. Is there a latitudinal gradient in the importance of biotic interactions? Annu. Rev. Ecol. Evol. Syst. 2009, 40, 245–269. [Google Scholar] [CrossRef]
  20. Freestone, A.L.; Osman, R.W.; Ruiz, G.M.; Torchin, M.E. Stronger predation in the tropics shapes species richness patterns in marine communities. Ecology 2011, 92, 983–993. [Google Scholar] [CrossRef]
  21. Parker, J.D.; Duffy, J.E.; Orth, R.J. Plant species diversity and composition: Experimental effects on marine epifaunal assemblages. Mar. Ecol. Prog. Ser. 2001, 224, 55–67. [Google Scholar] [CrossRef]
  22. Schneider, F.I.; Mann, K.H. Species specific relationships of invertebrates to vegetation in a seagrass bed. I. Correlational studies. J. Exp. Mar. Biol. Ecol. 1991, 145, 101–117. [Google Scholar] [CrossRef]
  23. Lee, P.Y.; Rotenberry, J.T. Relationships between bird species and tree species assemblages in forested habitats of eastern North America. J. Biogeogr. 2005, 32, 1139–1150. [Google Scholar] [CrossRef]
  24. Lambais, M.R.; Crowley, D.E. Bacterial diversity in tree canopies of the Atlantic forest. In Encyclopedia of Metagenomics; Nelson, K.E., Ed.; Springer: New York, NY, USA, 2015; pp. 49–54. [Google Scholar]
  25. Yakovis, E.L.; Artemieva, A.V.; Shunatova, N.N.; Varfolomeeva, M.A. Multiple foundation species shape benthic habitat islands. Oecologia 2008, 155, 785–795. [Google Scholar] [CrossRef]
  26. Morrow, K.M.; Moss, A.G.; Chadwick, N.E.; Liles, M.R. Bacterial associates of two Caribbean coral species reveal species-specific distribution and geographic variability. Appl. Environ. Microbiol. 2012, 78, 6438–6449. [Google Scholar] [CrossRef]
  27. Whitham, T.G.; DiFazio, S.P.; Schweitzer, J.A.; Shuster, S.M.; Allan, G.J.; Bailey, J.K.; Woolbright, S.A. Extending genomics to natural communities and ecosystems. Science 2008, 320, 492–495. [Google Scholar] [CrossRef]
  28. Pruitt, J.N.; Modlmeier, A.P. Animal personality in a foundation species drives community divergence and collapse in the wild. J. Anim. Ecol. 2015, 84, 1461–1468. [Google Scholar] [CrossRef]
  29. Campos, R.I.; Vasconcelos, H.L.; Ribeiro, S.P.; Neves, F.S.; Soares, J.P. Relationship between tree size and insect assemblages associated with Anadenanthera macrocarpa. Ecography 2006, 29, 442–450. [Google Scholar] [CrossRef]
  30. Carlsen, B.P.; Johnsen, G.; Berge, J.; Kuklinski, P. Biodiversity patterns of macro-epifauna on different lamina parts of Laminaria digitata and Saccharina latissima collected during spring and summer 2004 in Kongsfjorden, Svalbard. Polar Biol. 2007, 30, 939–943. [Google Scholar] [CrossRef]
  31. Taylor, A.; Burns, K. Epiphyte community development throughout tree ontogeny: An island ontogeny framework. J. Veg. Sci. 2015, 26, 902–910. [Google Scholar] [CrossRef]
  32. Christie, H.; Jørgensen, N.M.; Norderhaug, K.M.; Waage-Nielsen, E. Species distribution and habitat exploitation of fauna associated with kelp (Laminaria hyperborea) along the Norwegian coast. J. Mar. Biol. Assoc. UK 2003, 83, 687–699. [Google Scholar] [CrossRef]
  33. Lund, S. The marine algae of East Greenland I. Taxonomic part. Meddr. Gronland 1959, 156, 1–247. [Google Scholar]
  34. Kain, J.M. The reproductive phenology of nine species of Rhodophyta in the subtidal region of the Isle of Man. Br. Phycol. J. 1982, 17, 321–331. [Google Scholar] [CrossRef] [Green Version]
  35. Kain, J.M. Plant size and reproductive phenology of six species of Rhodophyta in subtidal Isle of Man. Br. Phycol. J. 1986, 21, 129–138. [Google Scholar] [CrossRef]
  36. Voskoboinikov, G.M.; Breeman, A.M.; Van den Hoek, C.; Makarov, V.N.; Schoschina, E.V. Influence of temperature and photoperiod on survival and growth of North East Atlantic isolates of Phycodrys rubens (Rhodophyta) from different latitudes. Bot. Mar. 1996, 39, 341–346. [Google Scholar] [CrossRef]
  37. Makarov, V.N.; Makarov, M.V.; Schoschina, E.V. Seasonal dynamics of growth in the Barents Sea seaweeds: Endogenous and exogenous regulation. Bot. Mar. 1999, 42, 43–49. [Google Scholar] [CrossRef]
  38. Schoschina, E.V. Seasonal and age dynamics of growth and reproduction of Phycodrys rubens (Rhodophyta) in the Barents and White Seas. Aquat. Bot. 1996, 55, 13–30. [Google Scholar] [CrossRef]
  39. Mileikovsky, S.A. Seasonal and daily dynamics in pelagic larvae of marine shelf bottom invertebrates in nearshore waters of Kandalaksha Bay (White Sea). Mar. Biol. 1970, 5, 180–194. [Google Scholar] [CrossRef]
  40. Tikhonenkov, D.V.; Mazei, Y.A.; Mylnikov, A.P. Species diversity of heterotrophic flagellates in White Sea littoral sites. Eur. J. Protistol. 2006, 42, 191–200. [Google Scholar] [CrossRef] [PubMed]
  41. Anderson, M.J. Permutation tests for univariate or multivariate analysis of variance and regression. Can. J. Fish. Aquat. Sci. 2001, 58, 626–639. [Google Scholar] [CrossRef]
  42. Buuren, S.V.; Fredriks, M. Worm plot: A simple diagnostic device for modeling growth reference curves. Stat. Med. 2001, 20, 1259–1277. [Google Scholar] [CrossRef] [PubMed]
  43. Stasinopoulos, D.M.; Rigby, R.A. Generalized additive models for location scale and shape (GAMLSS) in R. J. Stat. Softw. 2007, 23, 1–46. [Google Scholar] [CrossRef]
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 13 April 2019).
  45. Shunatova, N.; Nikishina, D.; Ivanov, M.; Berge, J.; Renaud, P.E.; Ivanova, T.; Granovitch, A. The longer the better: The effect of substrate on sessile biota in Arctic kelp forests. Polar Biol. 2018, 41, 993–1011. [Google Scholar] [CrossRef]
  46. Fredriksen, S.; Christie, H.; Andre Sæthre, B. Species richness in macroalgae and macrofauna assemblages on Fucus serratus L. (Phaeophyceae) and Zostera marina L.(Angiospermae) in Skagerrak, Norway. Mar. Biol. Res. 2005, 1, 2–19. [Google Scholar] [CrossRef]
  47. Lippert, H.; Iken, K.; Rachor, E.; Wiencke, C. Macrofauna associated with macroalgae in the Kongsfjord (Spitsbergen). Polar Biol. 2001, 24, 512–522. [Google Scholar]
  48. Mikhaylova, T.A.; Aristov, D.A.; Naumov, A.D.; Malavenda, S.S.; Savchenko, O.N.; Bijagov, K.L. Diversity and structure of epibenthic communities of the red algae zone at the White Sea. Polar Biol. 2019. [Google Scholar] [CrossRef]
  49. Wahl, M. Marine epibiosis. I. Fouling and antifouling: Some basic aspects. Mar. Ecol. Prog. Ser. 1989, 58, 175–189. [Google Scholar] [CrossRef]
  50. Sousa, W.P. Disturbance and patch dynamics on rocky intertidal shores. In The Ecology of Natural Disturbance and Patch Dynamics; White, P.S., Pickett, S.T.A., Eds.; Academic Press: New York, NY, USA, 1985; pp. 101–124. [Google Scholar]
  51. Keough, M.J. Effects of patch size on the abundance of sessile marine invertebrates. Ecology 1984, 65, 423–437. [Google Scholar] [CrossRef]
  52. Kendall, M.A.; Bowman, R.S.; Williamson, P.; Lewis, J.R. Annual variation in the recruitment of Semibalanus balanoides on the North Yorkshire coast 1969–1981. J. Mar. Biol. Assoc. UK 1985, 65, 1009–1030. [Google Scholar] [CrossRef]
  53. Yakovis, E.L.; Artemieva, A.V.; Fokin, M.V.; Varfolomeeva, M.A.; Shunatova, N.N. Synchronous annual recruitment variation in barnacles and ascidians in the White Sea shallow subtidal 1999–2010. Hydrobiologia 2013, 706, 69–79. [Google Scholar] [CrossRef]
  54. Shunkina, K.V. Reconstruction of the Life Cycle of Cheilostome Bryozoan Cribrilina Annulata (Gymnolaemata) Living on the Laminaria (The White Sea, Chupa Inlet). Master’s Thesis, St. Petersburg State University, Department of Invertebrate Zoology, St. Petersburg, Russia, 2010. [Google Scholar]
  55. Marcus, E. Beobachtungen und Versuche an lebeden Meeresbryozoen. Zool. Jahrb. Abt. Syst. Ökol. Geogr. Tiere. 1926, 52, 1–102. [Google Scholar]
  56. Ereskovsky, A.V. Reproduction cycles and strategies of the cold-water sponges Halisarca dujardini (Demospongiae, Halisarcida), Myxilla incrustans and Iophon piceus (Demospongiae, Poecilosclerida) from the White Sea. Biol. Bull. 2000, 198, 77–87. [Google Scholar] [CrossRef]
  57. Petralia, R.S.; Mattson, M.P.; Yao, P.J. Aging and longevity in the simplest animals and the quest for immortality. Ageing Res. Rev. 2014, 16, 66–82. [Google Scholar] [CrossRef] [Green Version]
  58. Thorson, G. Studies on the egg capsules and development of Arctic marine prosobranchs. Medd. Groenl. 1935, 100, 1–71. [Google Scholar]
  59. Anwar, N.A.; Richardson, C.A.; Seed, R. Age determination, growth rate and population structure of the horse mussel Modiolus modiolus. J. Mar. Biol. Assoc. UK 1990, 70, 441–457. [Google Scholar] [CrossRef]
  60. Bayne, B.L. Marine Mussels: Their Ecology and Physiology; Cambridge University Press: Cambridge, UK, 1976; 506p. [Google Scholar]
  61. Dobretsov, S.; Wahl, M. Larval recruitment of the blue mussel Mytilus edulis: The effect of flow and algae. J. Exp. Mar. Biol. Ecol. 2008, 355, 137–144. [Google Scholar] [CrossRef]
  62. Eggleston, D. Patterns of reproduction in the marine Ectoprocta of the Isle of Man. J. Nat. Hist. 1972, 6, 31–38. [Google Scholar] [CrossRef]
  63. Bergan, P. On the anatomy and reproduction biology in Spirorbis (Daudin). Nytt. Mag. Zool. 1953, 1, 1–26. [Google Scholar]
  64. Eckman, J.E.; Duggins, D.O. Life and death beneath macrophyte canopies: Effects of understory kelps on growth rates and survival of marine, benthic suspension feeders. Oecologia 1991, 87, 473–487. [Google Scholar] [CrossRef]
  65. Leichter, J.J.; Witman, J.D. Water flow over subtidal rock walls: Relation to distributions and growth rates of sessile suspension feeders in the Gulf of Maine Water flow and growth rates. J. Exp. Mar. Biol. Ecol. 1997, 209, 293–307. [Google Scholar] [CrossRef]
  66. Mullineaux, L.S.; Butman, C.A. Recruitment of encrusting benthic invertebrates in boundary-layer flows: A deep-water experiment on Cross Seamount. Limnol. Oceanogr. 1990, 35, 409–423. [Google Scholar] [CrossRef] [Green Version]
  67. Eckman, J.E.; Duggins, D.O. Larval settlement in turbulent pipe flows. J. Mar. Res. 1998, 56, 1285–1312. [Google Scholar] [CrossRef]
  68. Young, C.M.; Chia, F.S. Microhabitat-associated variability in survival and growth of subtidal solitary ascidians during the first 21 days after settlement. Mar. Biol. 1984, 81, 61–68. [Google Scholar] [CrossRef]
  69. Boaden, P.J.S.; O’Connor, R.J.; Seed, R. The composition and zonation of a Fucus serratus community in Strangford Lough, Co. Down. J. Exp. Mar. Biol. Ecol. 1975, 17, 111–136. [Google Scholar] [CrossRef]
  70. Gutiérrez, J.L.; Palomo, M.G. Increased algal fouling on mussels with barnacle epibionts: A fouling cascade. J. Sea Res. 2016, 112, 49–54. [Google Scholar] [CrossRef]
  71. Stebbing, A.R.D. Preferential settlement of a bryozoan and serpulid larvae on the younger parts of Laminaria fronds. J. Mar. Biol. Assoc. UK 1972, 52, 765–772. [Google Scholar] [CrossRef]
  72. Ryland, J.S. Physiology and ecology of marine bryozoans. Adv. Mar. Biol. 1977, 14, 285–443. [Google Scholar]
  73. Seed, R. Observations on the ecology of Membranipora (Bryozoa) and a major predator Doridella steinbergae (Nudibranchiata) along the fronds of Laminaria saccharina at Friday Harbor, Washington. J. Exp. Mar. Biol. Ecol. 1976, 24, 1–17. [Google Scholar] [CrossRef]
  74. Durante, K.M.; Chia, F.S. Epiphytism on Agarum fimbriatum: Can herbivore preferences explain distributions of epiphytic bryozoans? Mar. Ecol. Prog. Ser. Oldendorf 1991, 77, 279–287. [Google Scholar] [CrossRef]
  75. Hamamoto, K.; Mukai, H. Effects of larval settlement and post-settlement mortality on the distribution pattern and abundance of the spirorbid tube worm Neodexiospira brasiliensis (Grube)(Polychaeta) living on seagrass leaves. Mar. Ecol. 1999, 20, 251–272. [Google Scholar] [CrossRef]
  76. Jackson, J.B.C.; Buss, L.E.O. Alleopathy and spatial competition among coral reef invertebrates. Proc. Natl. Acad. Sci. USA 1975, 72, 5160–5163. [Google Scholar] [CrossRef] [Green Version]
  77. Buss, L.W.; Jackson, J.B.C. Competitive networks: Nontransitive competitive relationships in cryptic coral reef environments. Am. Nat. 1979, 113, 223–234. [Google Scholar] [CrossRef]
  78. McCook, L.; Jompa, J.; Diaz-Pulido, G. Competition between corals and algae on coral reefs: A review of evidence and mechanisms. Coral Reefs 2001, 19, 400–417. [Google Scholar] [CrossRef]
  79. Russ, G.R. Overgrowth in a marine epifaumal community: Competitive hierarchies and competitive networks. Oecologia 1982, 53, 12–19. [Google Scholar] [CrossRef]
  80. Jackson, J.B.C.; Winston, J.E. Ecology of cryptic coral reef communities. I. Distribution and abundance of major groups of encrusting organisms. J. Exp. Mar. Biol. Ecol. 1982, 57, 135–147. [Google Scholar] [CrossRef]
  81. Carballo, J.L. Effect of natural sedimentation on the structure of tropical rocky sponge assemblages. Ecoscience 2006, 13, 119–130. [Google Scholar] [CrossRef]
  82. Rützler, K. Impact of crustose clionid sponges on Caribbean reef corals. Acta Geol. Hisp. 2002, 37, 61–72. [Google Scholar]
  83. Porter, J.W.; Targett, N.M. Allelochemical interactions between sponges and corals. Biol. Bull. 1988, 175, 230–239. [Google Scholar] [CrossRef]
  84. De Voogd, N.J.; Becking, L.E.; Hoeksema, B.W.; Noor, A.; van Soest, R.W. Sponge interactions with spatial competitors in the Spermonde Archipelago. Bollettino Musei degli Istituti Biologici Universita Genova 2003, 68, 253–261. [Google Scholar]
  85. Fricke, A.; Titlyanova, T.V.; Nugues, M.M.; Bischof, K. Depth-related variation in epiphytic communities growing on the brown alga Lobophora variegata in a Caribbean coral reef. Coral Reefs 2011, 30, 967–973. [Google Scholar] [CrossRef]
  86. Jennings, J.G.; Steinberg, P.D. Phlorotannins versus other factors affecting epiphyte abundance on the kelp Ecklonia radiata. Oecologia 1997, 109, 461–473. [Google Scholar] [CrossRef]
  87. Arrontes, J. Composition, distribution on host, and seasonality of epiphytes on three intertidal algae. Bot. Mar. 1990, 33, 205–212. [Google Scholar] [CrossRef]
  88. Bernstein, B.B.; Jung, N. Selective pressures and coevolution in a kelp canopy community in southern California. Ecol. Monogr. 1979, 49, 335–355. [Google Scholar] [CrossRef]
  89. O’Connor, R.J.; Boaden, P.J.S.; Seed, R. Niche breadth in Bryozoa as a test of competition theory. Nature 1975, 256, 307. [Google Scholar] [CrossRef]
  90. Ryland, J.S. Observations on some epibionts of gulf-weed, Sargassum natans (L.) Meyen. J. Exp. Mar. Biol. Ecol. 1974, 14, 17–25. [Google Scholar] [CrossRef]
  91. Boaden, P.J.S.; O’Connor, R.J.; Seed, R. The fauna of a Fucus serratus L. community: Ecological isolation in sponges and tunicates. J. Exp. Mar. Biol. Ecol. 1976, 21, 249–267. [Google Scholar] [CrossRef]
  92. Matias, M.G.; Underwood, A.J.; Hochuli, D.F.; Coleman, R.A. Independent effects of patch size and structural complexity on diversity of benthic macroinvertebrates. Ecology 2010, 91, 1908–1915. [Google Scholar] [CrossRef] [Green Version]
  93. Johansson, P.; Rydin, H.; Thor, G. Tree age relationships with epiphytic lichen diversity and lichen life history traits on ash in southern Sweden. Ecoscience 2007, 14, 81–91. [Google Scholar] [CrossRef]
  94. Underwood, A.J.; Chapman, M.G.; Connell, S.D. Observations in ecology: You can’t make progress on processes without understanding the patterns. J. Exp. Mar. Biol. Ecol. 2000, 250, 97–115. [Google Scholar] [CrossRef]
Figure 1. Dominant species composition (on piecharts) and average total cover (in white boxes) of epibiosis on young and old plant parts of Phycodrys rubens in July and September 2015. Piechart diameter denotes log total percent covers (frames pooled). Po—Porifera, My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia). Ch—Celleporella hyalina (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa). Ca—Circeis armoricana (Polychaeta). Cs—Calycella syringa (Hydrozoa), Oi—Orthopyxis integra (Hydrozoa), O—all other species.
Figure 1. Dominant species composition (on piecharts) and average total cover (in white boxes) of epibiosis on young and old plant parts of Phycodrys rubens in July and September 2015. Piechart diameter denotes log total percent covers (frames pooled). Po—Porifera, My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia). Ch—Celleporella hyalina (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa). Ca—Circeis armoricana (Polychaeta). Cs—Calycella syringa (Hydrozoa), Oi—Orthopyxis integra (Hydrozoa), O—all other species.
Diversity 11 00080 g001
Figure 2. Principal coordinates ordination (PCO) of sessile species percent covers on young and old parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Dashed lines connect young and old parts of a single plant in those having both. Species with R2 (in brackets) > 0.2: Po—Porifera; My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae (Crisia eburnea and Crisiella producta) (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa); Cs—Calycella syringa (Hydrozoa), Lp—Lafoea pocillum (Hydrozoa); Ca—Circeis armoricana (Polychaeta); logS—log-transformed total surface area of the corresponding plant part (young or old).
Figure 2. Principal coordinates ordination (PCO) of sessile species percent covers on young and old parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Dashed lines connect young and old parts of a single plant in those having both. Species with R2 (in brackets) > 0.2: Po—Porifera; My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae (Crisia eburnea and Crisiella producta) (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa); Cs—Calycella syringa (Hydrozoa), Lp—Lafoea pocillum (Hydrozoa); Ca—Circeis armoricana (Polychaeta); logS—log-transformed total surface area of the corresponding plant part (young or old).
Diversity 11 00080 g002
Figure 3. Total percent cover of epibiosis on young and old parts of Phycodrys rubens in September 2015 by surface area of these parts. ‘L’ marks auxiliary non-randomly sampled 5 largest plants per frame. Fit lines plotted according to significant relationships (see Table 5 for details on beta-regression models).
Figure 3. Total percent cover of epibiosis on young and old parts of Phycodrys rubens in September 2015 by surface area of these parts. ‘L’ marks auxiliary non-randomly sampled 5 largest plants per frame. Fit lines plotted according to significant relationships (see Table 5 for details on beta-regression models).
Diversity 11 00080 g003
Figure 4. Principal coordinates ordination (PCO) of sessile species percent covers on old parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Labels (1 or 2) indicate the replicate frame number. Species with R2 (in brackets) > 0.2: Mo—Molgula sp. (Ascidiacea); My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae f. gen. spp. (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa), Ta—Tegella armifera (Bryozoa); Cs—Calycella syringa (Hydrozoa), Lp—Lafoea pocillum (Hydrozoa); Ca—Circeis armoricana; logS—log-transformed surface area of an old plant part.
Figure 4. Principal coordinates ordination (PCO) of sessile species percent covers on old parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Labels (1 or 2) indicate the replicate frame number. Species with R2 (in brackets) > 0.2: Mo—Molgula sp. (Ascidiacea); My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae f. gen. spp. (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa), Ta—Tegella armifera (Bryozoa); Cs—Calycella syringa (Hydrozoa), Lp—Lafoea pocillum (Hydrozoa); Ca—Circeis armoricana; logS—log-transformed surface area of an old plant part.
Diversity 11 00080 g004
Figure 5. Principal coordinates ordination (PCO) of sessile species percent covers on young parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Labels (1 or 2) indicate the replicate frame number. Species with R2 (in brackets) > 0.2: My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae (Crisia eburnea and Crisiella producta) (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa); Cs—Calycella syringa (Hydrozoa), Tm—Tiaropsis multicirrata (Hydrozoa); Ca—Circeis armoricana (Polychaeta); logS—log-transformed surface area of a young plant part.
Figure 5. Principal coordinates ordination (PCO) of sessile species percent covers on young parts of Phycodrys rubens in September 2015. Bray-Curtis similarity on fourth root transformed data. Labels (1 or 2) indicate the replicate frame number. Species with R2 (in brackets) > 0.2: My—juvenile Mytilidae (Musculus discors and Modiolus modiolus) (Bivalvia); Cc—Callopora craticula (Bryozoa), Ch—Celleporella hyalina (Bryozoa), Cr—Crisiidae (Crisia eburnea and Crisiella producta) (Bryozoa), Cra—Cribrilina annulata (Bryozoa), Ep—Electra pilosa (Bryozoa), Lv—Lichenopora verrucaria (Bryozoa); Cs—Calycella syringa (Hydrozoa), Tm—Tiaropsis multicirrata (Hydrozoa); Ca—Circeis armoricana (Polychaeta); logS—log-transformed surface area of a young plant part.
Diversity 11 00080 g005
Figure 6. Area covered and relative percent covers of dominant sessile species on young and old parts of Phycodrys rubens by plant part surface area in September 2015. Fit lines plotted where plant part area effect or its interaction with Site were significant (see Table 8 for details on models). Dots mark auxiliary non-randomly sampled 5 largest plants per frame (see Methods for details).
Figure 6. Area covered and relative percent covers of dominant sessile species on young and old parts of Phycodrys rubens by plant part surface area in September 2015. Fit lines plotted where plant part area effect or its interaction with Site were significant (see Table 8 for details on models). Dots mark auxiliary non-randomly sampled 5 largest plants per frame (see Methods for details).
Diversity 11 00080 g006
Table 1. Effects of plant part age (fixed, young or old), month (fixed, July or September), and location (fixed, sites K or V) on total cover of Phycodrys rubens epibiosis in 2015: beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. Auxiliary non-randomly sampled (in September only) five largest plants per frame not included (see Materials and Methods). Significant terms in a mean model highlighted in bold.
Table 1. Effects of plant part age (fixed, young or old), month (fixed, July or September), and location (fixed, sites K or V) on total cover of Phycodrys rubens epibiosis in 2015: beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. Auxiliary non-randomly sampled (in September only) five largest plants per frame not included (see Materials and Methods). Significant terms in a mean model highlighted in bold.
Plant AgeModelSource of VariationEstimateSEt-Valuep
Plants having both old and young parts (with Plant ID as a random blocking factor nested in Site × Month; n = 51 plants)Mean (Mu)Intercept−4.28180.2837−15.091<0.001
Site [level ‘K’]1.08150.51882.0850.040
Age [level ‘Old’]2.57560.36707.018<0.001
Month [level ‘September’])1.92480.30536.305<0.001
Site × Age1.14470.62801.8230.072
Site × Month0.91850.54361.6900.095
Age×Month1.33770.49002.7300.008
Site × Age × Month0.38320.79510.4820.631
Variance (Sigma)Intercept−2.00560.2469−8.121<0.001
Site [level ‘K’]−0.16870.3894−0.4330.666
Age [level ‘Old’]1.17520.35583.3030.001
Month [level ‘September’]−0.13650.3419-−0.3990.691
Site × Age0.34490.54090.6380.525
Site × Month0.30680.49990.6140.541
Age × Month0.73070.51781.4110.162
Site × Age × Month0.90980.73301.2410.218
Plants having only young parts (n = 39 plants)Mean (Mu)Intercept−4.47940.2277−19.674<0.001
Site [level ‘K’]0.56790.81580.6960.492
Month [level ‘September’]2.44660.28398.618<0.001
Site × Month−1.17460.8399−1.3980.172
Variance (Sigma)Intercept−3.34970.5212−6.427<0.001
Site [level ‘K’]1.87910.80052.3470.025
Month [level ‘September’]2.17930.55943.896<0.001
Site × Month−2.88020.8503−3.3870.002
Table 2. Average total cover, number of species and Shannon-Wiener diversity index (H’) of epibiosis on Phycodrys rubens by month, site and plant part age in 2015 (frames pooled).
Table 2. Average total cover, number of species and Shannon-Wiener diversity index (H’) of epibiosis on Phycodrys rubens by month, site and plant part age in 2015 (frames pooled).
ParameterMonthSite KSite V
YoungOldYoungOld
Average total cover, %July0.8 ± 0.416.2 ± 4.01.3 ± 0.415.4 ± 3.3
September7.1 ± 0.635.9 ± 8.710.4 ± 1.423.8 ± 6.8
Total number of species (average number of species per plant)July10
(3.9 ± 0.5)
15
(4.8 ± 0.6)
16
(6.7 ± 0.4)
23
(7.1 ± 0.8)
September26
(9.8 ± 0.5)
35
(11.3 ± 1.7)
27
(8.5 ± 0.5)
27
(8.8 ± 1.2)
Total H’ (average H’ per plant)July0.83
(0.67 ± 0.12)
1.20
(0.73 ± 0.14)
1.34
(0.86 ± 0.09)
1.67
(1.16 ± 0.12)
September1.63
(1.33 ± 0.04)
1.47
(1.08 ± 0.17)
2.09
(1.36 ± 0.07)
1.39
(1.29 ± 0.08)
Table 3. Effects of plant part age (young or old), month (July or September), and location (sites K or V) on multivariate community structure of epibiosis on Phycodrys rubens: type III sum of squares PERMANOVA for 9999 permutations on fourth root transformed covers of sessile epibionts, Bray-Curtis similarity, frames pooled. Significant terms highlighted in bold. [f]—fixed effect, [r]—random effect.
Table 3. Effects of plant part age (young or old), month (July or September), and location (sites K or V) on multivariate community structure of epibiosis on Phycodrys rubens: type III sum of squares PERMANOVA for 9999 permutations on fourth root transformed covers of sessile epibionts, Bray-Curtis similarity, frames pooled. Significant terms highlighted in bold. [f]—fixed effect, [r]—random effect.
Plant AgeSource of VariationdfSSMSPseudo-FpUnique Permutations
Plants having both old and young partsSite [f]127,41127,41118.560.00019941
Month [f]131,63531,63521.420.00019945
Age [f]124,17224,17225.840.00019930
S×M [f]1718371834.870.00019945
S×A [f]1292129213.120.00769940
M×A [f]116,01516,01517.120.00019955
S × M × A [f]1141614161.510.18509944
Plant ID (S×M) [r]4769,38914761.580.00029733
Error4743,960935
Plants having only young partsSite [f]1277927792.410.03619952
Month [f]111,99011,99010.400.00019949
S×M [f]1390639063.390.00559954
Error3540,3381153
Table 4. Mean percent covers and mean individual sizes (approximated area occupied) of top abundant species on young and old parts of the same Phycodrys rubens pairwise compared (Wilcoxon matched pairs test, P). K—Site K, V—Site V. n.v.—no variance or sample size too small. n.c.—no co-occurrence on old and young parts of a single plant. N—sample size. Significantly higher mean in a pair highlighted in bold.
Table 4. Mean percent covers and mean individual sizes (approximated area occupied) of top abundant species on young and old parts of the same Phycodrys rubens pairwise compared (Wilcoxon matched pairs test, P). K—Site K, V—Site V. n.v.—no variance or sample size too small. n.c.—no co-occurrence on old and young parts of a single plant. N—sample size. Significantly higher mean in a pair highlighted in bold.
Mean Percent Cover, %Individual Mean Area, mm2
SpeciesSiteYoung partsOld partsPNYoung partsOld partsPN
Circeis armoricanaK0.02±0.011.01±0.280.002120.044±0.0050.534±0.0420.00311
V0.03±0.011.06±0.280.006130.039±0.0010.512±0.0380.00510
Cribrilina annulataK0.00±0.005.33±2.170.004120.127±0.0002.793±0.4200.0187
V0.01±0.000.67±0.520.508130.140±0.0102.593±0.3130.0684
Electra pilosaK0.00±0.000.00±0.00n.v.12n.c.0
V0.00±0.000.17±0.130.68613n.c.0
Celleporella hyalinaK0.00±0.000.01±0.010.68612n.c.0
V0.39±0.071.83±0.580.013130.625±0.0552.133±0.2650.00411
PoriferaK0.37±0.358.72±4.370.0281211.369±8.06920.981±17.081n.v.2
V0.61±0.488.20±3.270.0171341.056±21.57324.138±7.3060.4654
Juvenile mytilidsK0.00±0.000.28±0.230.345121.73214.628n.v.1
V0.02±0.010.68±0.240.002130.629±0.1181.086±0.6550.4636
Circeis armoricanaK1.38±0.091.79±0.280.326250.250±0.0140.578±0.0200.00021
V1.31±0.140.55±0.160.003210.325±0.0110.451±0.0440.01016
Cribrilina annulataK2.78±0.333.38±0.620.300251.005±0.0691.510±0.1520.00222
V0.38±0.070.22±0.090.099211.229±0.2381.463±0.1960.2088
Electra pilosaK0.05±0.020.06±0.020.664250.639±0.1730.722±0.1470.50811
V1.50±0.333.61±0.910.016211.923±0.2872.755±0.6510.06321
Celleporella hyalinaK0.17±0.050.27±0.070.092250.697±0.0661.620±0.2890.01517
V2.92±0.400.82±0.140.000211.408±0.1521.599±0.1440.28717
PoriferaK1.34±0.4219.14±5.680.000258.216±1.30918.355±6.4900.14719
V0.72±0.3315.92±4.330.000218.683±1.84253.231±36.6580.27913
Juvenile mytilidsK0.19±0.112.23±0.690.000251.955±0.5487.417±1.4720.00216
V0.12±0.055.43±1.530.001212.170±0.8735.808±1.2770.06012
Table 5. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on total cover of epibiosis on young and old parts of Phycodrys rubens in September 2015: beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. Significant terms in a mean model highlighted in bold.
Table 5. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on total cover of epibiosis on young and old parts of Phycodrys rubens in September 2015: beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. Significant terms in a mean model highlighted in bold.
AgeModelSource of VariationEstimateSEt-Valuep
Young parts
(n = 80)
Mean (Mu)Intercept−2.18370.1084−20.15<0.001
Site [level ‘K’]−0.35420.1330−2.660.009
Area−0.00350.0014−2.590.012
Site × Area0.00050.00160.360.718
Variance (Sigma)Intercept−1.28910.1561−8.26<0.001
Site [level ‘K’]−0.78230.2157−3.630.001
Area−0.01820.0040−4.52<0.001
Site × Area0.01500.00453.400.001
Old parts
(n = 46)
Mean (Mu)Intercept−1.66140.3020−5.50<0.001
Site [level ‘K’]1.54200.42813.600.001
Area0.07540.02363.200.003
Site × Area−0.09080.0241−3.760.001
Variance (Sigma)Intercept−0.61860.3181−1.950.059
Site [level ‘K’]2.08410.40575.14<0.001
Area0.01180.02520.470.642
Site × Area−0.07100.0269−2.630.012
Table 6. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on the number of speicies and diversity (H’) on young and old plant parts of Phycodrys rubens in September 2015: type III sum of squares ANCOVAs, frames pooled. Arrows show the sign of the relationship where covariate effect or interaction is significant. Significant terms are highlighted in bold.
Table 6. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on the number of speicies and diversity (H’) on young and old plant parts of Phycodrys rubens in September 2015: type III sum of squares ANCOVAs, frames pooled. Arrows show the sign of the relationship where covariate effect or interaction is significant. Significant terms are highlighted in bold.
AgeParameterSource of VariationdfSum of SquaresF-ValuepSign
Young partsH’ (based on % cover)Site10.00220.0240.877
Area11.783819.375<0.001
Site×Area11.177912.7940.001K V
Error766.9969
Number of speciesSite119.22.96580.089
Area1671.4103.8066<0.001
Site × Area171.311.02850.001K V
Error76491.5
Old partsH’ (based on % cover)Site10.49891.94860.170
Area10.21830.85260.361
Site × Area10.61832.41520.128
Error4210.752
Number of speciesSite136.12.89090.096
Area1551.844.1654<0.001
Site × Area18.90.71510.403
Error42524.7
Table 7. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on the multivariate community structure of epibiosis on young and old plant parts of Phycodrys rubens in September 2015: type III sum of squares PERMANOVA for 9999 permutations on fourth root transformed covers of sessile epibionts, Bray-Curtis similarity, frames pooled. Significant terms are highlighted in bold.
Table 7. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on the multivariate community structure of epibiosis on young and old plant parts of Phycodrys rubens in September 2015: type III sum of squares PERMANOVA for 9999 permutations on fourth root transformed covers of sessile epibionts, Bray-Curtis similarity, frames pooled. Significant terms are highlighted in bold.
AgeSource of VariationdfSSMSPseudo-FpUnique Permutations
Young partsArea1549854986.620.00019952
Site113,84113,84116.690.00019936
Site × Area1159315931.920.06519937
Error7663,026829
Old partsArea1470947083.970.00019948
Site1978097808.250.00019936
Site × Area1204320431.720.07829946
Error4249,7821185
Table 8. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on relative abundances of 5 top abundant species on young and old plant parts of Phycodrys rubens in September 2015: zero-inflated (logit link) beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. I—intercept, Sources of variation: S—Site (K), A—plant part surface area. The model did not converge for E. pilosa on old blades since all the zero observations were at Site K, thus two separate models for the two sites were fitted. Significant terms in mean and zero models highlighted in bold.
Table 8. Effects of plant part size (surface Area of the corresponding part, cm2) and location (Sites K or V) on relative abundances of 5 top abundant species on young and old plant parts of Phycodrys rubens in September 2015: zero-inflated (logit link) beta-regression, mean (logit link) and variance (log link) modeled, frames pooled. I—intercept, Sources of variation: S—Site (K), A—plant part surface area. The model did not converge for E. pilosa on old blades since all the zero observations were at Site K, thus two separate models for the two sites were fitted. Significant terms in mean and zero models highlighted in bold.
SpeciesMean Model (Mu)Variance Model (Sigma)Zeros Model (Nu)
Estimate (SE)pEstimate (SE)pEstimate (SE)p
Young Parts (n = 80)
Cribrilina annulataI−2.4951
(0.1870)
<0.0012.0878
(0.2763)
<0.001−0.1898
(0.7094)
0.790
S2.1997
(0.2236)
<0.001−0.0523
(0.3811)
0.892−21.2419
(5354.2603)
0.997
A0.0007
(0.0028)
0.8060.0238
(0.0075)
0.002−0.1596
(0.1069)
0.140
S × A0.0006
(0.0030)
0.842−0.0180
(0.0083)
0.033−0.1596
(71.9445)
0.998
Electra pilosaI−1.7665
(0.1845)
<0.0011.6823
(0.2641)
<0.0012.0963
(2.0629)
0.313
S1.9029
(0.2571)
<0.0012.0667
(0.3928)
<0.001−2.3008
(2.1404)
0.286
A0.0030
(0.0035)
0.39580.0125
(0.0073)
0.089−1.1162
(0.7365)
0.134
S × A−0.0073
(0.0040)
0.0756−0.0075
(0.0081)
0.3591.0704
(0.7370)
0.151
Celleporella hyalinaI−0.8769
(0.1452)
<0.0011.5891
(0.2785)
<0.0011.9564
(2.0067)
0.333
S−2.3326
(0.2234)
<0.0011.7150
(0.3971)
<0.001−1.6762
(2.1876)
0.446
A0.0019
(0.0017)
0.2710.0334
(0.0076)
<0.001−1.0639
(0.7096)
0.138
S × A0.0059
(0.0025)
0.019−0.0264
(0.0084)
0.0020.9373
(0.7146)
0.194
Circeis armoricanaI−2.1489
(0.1509)
<0.001−0.8095
(0.1605)
<0.001no zeros
S0.8676
(0.1785)
<0.001−0.4126
(0.2279)
0.074
A0.0063
(0.0024)
0.012−0.0098
(0.0043)
0.025
S × A−0.0046
(0.0026)
0.0770.0065
(0.0048)
0.180
PoriferaI−0.8429
(0.2670)
0.0020.6823
(0.3210)
0.037−0.1903
(0.4041)
0.639
S0.0421
(0.3471)
0.9040.8275
(0.4501)
0.0700.4485
(0.6518)
0.494
A0.0221
(0.0049)
<0.0010.0341
(0.0089)
<0.001−0.0103
(0.0123)
0.405
S × A0.0116
(0.0056)
0.040−0.0303
(0.0095)
0.002−0.0272
(0.0232)
0.244
Old Parts (n = 46)
Cribrilina annulataI−1.5135
(0.7367)
0.0464.3858
(1.5457)
0.0076.0933
(2.6097)
0.024
S−0.0450
(0.7783)
0.954−3.0769
(1.5904)
0.0606.4730
(2.7646)
0.024
A0.1484
(0.0434)
0.0010.0246
(0.0779)
0.7540.4881
(0.1982)
0.018
S×A0.1505
(0.0439)
0.0010.0137
(0.0796)
0.8650.2536
(0.2567)
0.329
Electra pilosa
(separate models for the two sites)
Site V: I−1.1564 (0.4217)0.0140.1092 (0.4254)0.801no zeros
Site V: A−0.0509 (0.0290)0.117−0.0420 (0.0343)0.237
Site K: I−4.9332 (0.4195)<0.0014.9017 (0.6593)<0.0011.5588 (0.8440)0.080
Site K: A−0.0132 (0.0134)0.3390.0141 (0.0213)0.5140.1572 (0.074)0.048
Juvenile mytilidsI−1.4250
(0.3840)
<0.0010.5747
(0.5548)
0.3062.5499
(1.9499)
0.198
S−0.2370
(0.4135)
0.570−0.5368
(0.6891)
0.440−0.2620
(2.3573)
0.912
A−0.0089
(0.0195)
0.6500.1132
(0.0412)
0.009−0.8268
(0.5108)
0.113
S × A−0.0006
(0.0196)
0.9770.0372
(0.0447)
0.410−0.0138
(0.6892)
0.984
Circeis armoricanaI−3.2273
(0.3373)
<0.0013.7427
(0.5848)
<0.0010.7101
(1.0069)
0.485
S1.2267
(0.4094)
0.005−2.1756
(0.7120)
0.00451.2633
(299.8967)
0.865
A−0.0327
(0.0247)
0.1920.0015
(0.0392)
0.970−0.2601
(0.1598)
0.111
S × A0.0317
(0.0251)
0.2140.0577
(0.0428)
0.185−50.1948
(286.5729)
0.862
PoriferaI−0.7947
(0.3675)
0.0361.3691
(0.4779)
0.0071.0340
(2.8280)
0.716
S0.7172
(0.4500)
0.119−1.2206
(0.5771)
0.040−15.466
(225.1950)
0.946
A0.0426
(0.0294)
0.155−0.0216
(0.0355)
0.546−1.1330
(1.1680)
0.338
S × A−0.0469
(0.0298)
0.1230.0883
(0.0387)
0.0281.1330
(9.672)
0.907

Share and Cite

MDPI and ACS Style

Chava, A.; Artemieva, A.; Yakovis, E. Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea. Diversity 2019, 11, 80. https://doi.org/10.3390/d11050080

AMA Style

Chava A, Artemieva A, Yakovis E. Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea. Diversity. 2019; 11(5):80. https://doi.org/10.3390/d11050080

Chicago/Turabian Style

Chava, Alexandra, Anna Artemieva, and Eugeniy Yakovis. 2019. "Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea" Diversity 11, no. 5: 80. https://doi.org/10.3390/d11050080

APA Style

Chava, A., Artemieva, A., & Yakovis, E. (2019). Plant Part Age and Size Affect Sessile Macrobenthic Assemblages Associated with a Foliose Red Algae Phycodrys rubens in the White Sea. Diversity, 11(5), 80. https://doi.org/10.3390/d11050080

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