Abstract
The effect of the presence of invasive macroalgae Caulerpa cylindracea Sonder in the seagrass meadow Cymodocea nodosa (Ucria) Ascherson was studied by comparing the dynamics of biological and physicochemical parameters in the tissues of these two macrophytes and in the sediments underlying an invaded meadow (mixed settlement) and a C. cylindracea monospecific settlement. The study was conducted for 15 months, encompassing two summers during 2017 and 2018 (July to September) when maximum C. nodosa growth occurs. During 2017 C. cylindracea biomasses in the mixed settlement (79.5 ± 28.2 to 264.6 ± 65.1 g m–2) were lower than in its monospecific stands (113.4 ± 48.0 to 399.3 ± 56.3 g m–2). In the same period, less reducing conditions in the sediment underlying the mixed settlement were indicated by deeper redox transition depths (RTD: 8–12 mm) than those observed below the C. cylindracea monospecific settlement (RTD: 4–7 mm). In June 2018, C. cylindracea proliferated in both settlements reaching very similar biomasses that were maintained until September 2018 (mixed: 131.5 ± 23.0 to 172.5 ± 16.3 g m–2; monospecific: 162.8 ± 32.5 to 178.8 ± 30.0 g m–2). In parallel, a considerable lowering of RTD (5–7 mm) under the mixed settlement indicated the progression of stronger reducing conditions similar to those observed under the monospecific settlement (RTD: 0–7 mm). This alteration was followed by a decrease in C. nodosa below-ground biomass (89.3 ± 16.0 to 140.3 ± 24.3 g m–2), that became considerably lower than in the same period of 2017 (242.3 ± 44.3 to 346.9 ± 32.1 g m–2). At the same time, the above-ground biomass of C. nodosa (72.3 ± 14.8 to 110.3 ± 13.4 g m–2) showed no difference to the summer of 2017 (69.0 ± 15.4 to 116.0 ± 37.4 g m–2). The resulting increase of the above- to below-ground biomass ratio indicated the disruption of the meadow stability. More intense spawning of C. cylindracea in the mixed settlement during the summer 2018 hindered its expected proliferation in October 2018, while the below-ground biomass of C. nodosa increased concomitantly with the deepening of the RTD suggesting a possible recovery of the meadow stability.
Effects of the presence of the invasive macroalgae C. cylindracea in the seagrass meadow evidenced by substantial loss in below-ground biomass of C. nodosa and lowering of the redox transition depth in the sediment underlying the mixed settlement.
Introduction
Seagrasses are one of the most important habitat-forming species in Mediterranean coastal areas and play a key role in the preservation of marine biodiversity and carbon sequestration (; ). They can be threatened by nutrient loading, warming, sediment runoff, physical disturbance, intense grazing, algal blooms, and invasion by non-indigenous species (; ; , ; ). Generally, the introduction of the tropical invasive macroalgae of the genus Caulerpa in the Mediterranean Sea was recognized as a significant environmental concern due to its modification of benthic communities and reduction of biodiversity on large spatial scales (; ; ; ; ).
Invasive Caulerpa species can strongly compete with native species in temperate areas, leading to their decline and constituting monospecific beds (). In seagrass habitats, this alga overgrows their rhizomes; interacting in both the below- and above-ground tissues, affecting nutrient acquisition and light availability (). In contrast, seagrass provides a sub-canopy microhabitat to Caulerpa sp. that differs from the non-vegetated sediment being characterized by canopy shading (), slowing down water motion () and releasing secondary metabolites that can interfere with understory species and grazers. In temperate areas, the extensive spreading of Caulerpa was ascribed to the successful adaptation to seasonal changes in environmental conditions (; ; ) in addition to both modes of propagation, by vegetative and sexual reproduction (). Numerous short-term studies relate the presence of Caulerpa to deterioration processes of seagrass meadows, by competitive displacement or by altering the sediment quality which becomes more suitable for the invasive macroalgae growth (; ; ).
The sediments underlying seagrasses are generally less hypoxic, present lower levels of nitrogen fixation and lower sulfide concentrations than those settled by Caulerpa sp. (; ; ; ; ; ). The presence of Caulerpa modifies benthic microbial community composition (), as well as patterns of bacterially mediated organic matter diagenesis by stimulating N2 fixation due to the release of secondary metabolites through plant rhizoids (), and modifying other microbial processes in the sediments. Sulfate-reducing bacteria elevate the sulfate reduction rates and H2S production in such sediments. Hydrogen sulfide is highly toxic to plants, and its possible intrusion into their tissues may reduce the growth and survival of the seagrass ().
Caulerpa cylindracea Sonder was first identified in the year 2000 in the central Adriatic () and successively spread north reaching the coldest area of its expansion in the Mediterranean (). This invasive macroalgae colonized all kinds of substrates, from rocky bottoms to seagrass meadows (; ). Given the wide distribution of C. cylindracea along the northern Adriatic coast and the need for data on the ecological status of seagrass meadows in this area, a study was conducted to investigate the growth cycle of C. cylindracea in the seagrass meadow Cymodocea nodosa (Ucria) Ascherson and the possible effect of their interactions on the sediment’s physico-chemical properties. We hypothesized that the presence of C. cylindracea can disturb the stability of the C. nodosa meadow by modifying the sediment conditions that concurrently favor the growth of C. cylindracea and limit it for C. nodosa. An effect of the presence of C. cylindracea in the C. nodosa meadow was evaluated by comparing the physiological and growth state of the macrophytes, environmental conditions in the water column and sediment between the summer of 2017 and the summer of 2018, the periods of C. nodosa maximum growth. This study aims to offer further insights to the complex interactions between invasive and native macrophytes and their relative importance in the maintenance of the seagrass habitat stability.
Materials and Methods
Study Site and Samplings
Funtana Bay is located along the western coast of Istria (Croatia) in the northern Adriatic, 30 km northwest of Rovinj (45°10′42″N; 13°35′51″E, Figure 1). The coast of Funtana Bay is rocky, moderately indented, and partially exposed to waves coming from the west and northwest. The sea floor is partly overgrown with C. cylindracea, which can be found alone or mixed with the seagrass C. nodosa. A monthly field survey was carried out from July 2017 to October 2018. Samplings were performed at three proximate sites; in the C. nodosa meadow invaded by the green macroalgae C. cylindracea (mixed settlement, sampling depth: 2–2.5 m), in the C. cylindracea settlement (monospecific settlement, sampling depth: 1–1.5 m) and in a non-vegetated area (sampling depth: 1–1.5 m).
FIGURE 1
Seawater for nutrients, chlorophyll a (Chl a), particulate matter (PM) concentration, and prokaryotic abundance (PA) analysis was collected with plastic bottles (10 L). C. nodosa (with rhizomes and roots), C. cylindracea thalli, and macroalgae were collected by snorkeling using three square frames (20 × 20 cm) scattered randomly in the areas of highest vegetation presence. Sediment samples were collected within the vegetated (mixed and monospecific settlement) and in the non-vegetated sediment by divers using plastic corers (length, 15 cm; area 15.9 cm2). The cores for granulometric composition, total lipids, and PA were cut into 1 cm sections and lyophilized with the exception of sections used for determining the abundance of prokaryotes which were weighed (2 g) and fixed with formaldehyde (final conc. 4% v/v).
Environmental Parameters
Temperature (T) was continuously measured at 30 min. intervals by HOBO pendant Temp/Light Data Loggers (Onset, United States) placed in the mixed and monospecific settlement. Salinity (S) was measured on each sampling by a pIONeer 65 probe (Radiometer analytical, Copenhagen). Nutrients (NO3, NO2, NH4, PO4, and SiO4) were analyzed according to . Chl a was determined using a fluorometer after seawater filtration through a Whatman GF/F filter (). Particulate matter (PM) was determined by weighing after filtering 5 L of seawater on pre-weighed, combusted Whatman GF/F filters, and dried at 60°C.
Prokaryotic Abundance Determination
Prokaryotic abundance in seawater was determined after fixing 2 mL of samples with formaldehyde (final conc. 4% v/v) and staining with 4,6-diamidino-2-phenylindol (DAPI, 1 μg mL–1 final conc.) for 10 min (). Prokaryotes in sediments were detached by adding Tween 80 (50 μL) and ultrasonicating for 15 min (). After the sonication, DAPI (final conc. 5 μg/mL) was added to 1 mL of the supernatant. Samples stained with DAPI were then filtered through black polycarbonate filters (Whatman, Nuclepore, 0.22 μm) and prokaryotes were counted under an epifluorescence microscope (Zeiss Axio Imager Z1).
Biometry of Macrophytes
The macrophytes from each square were thoroughly washed with seawater to eliminate the associated sediment, and separated into three groups of macrophytes: C. cylindracea thalli, C. nodosa (leaves and rhizomes with roots), and macroalgae thalli. For C. nodosa, shoots were counted and the length of the longest leaf of each shoot was measured. Macroalgae from each sample were placed on the surface that has an area equal to the sample in situ (20 × 20 cm) and identified to species level. The cover of each thallus, representing the surface covered in orthogonal projection, was determined according to . Species of macroalgae were affiliated to the following morpho-functional groups: corticated, articulated calcareous, filamentous, foliose, leathery, and crustose (). Individual samples were rinsed again, weighed, dried at 60°C for 48 h, re-weighed, and the dry mass expressed in g m–2.
Sediment Granulometric Composition
Each sample was wet sieved through standard ASTM sieves (4-, 2-, 1-, 0.5-, 0.25-, 0.125-, 0.063-mm mesh size). The fraction passing through the 0.063-mm sieve was analyzed according to standard sedigraph procedures (). The retained particles were dried and weighed. Merging the data acquired by both techniques a continuous grain size range was obtained, analyzed with the statistic package Gradistat v 6.0, and classified according to . Sediment permeability was calculated following the empirical relation by based on the median grain size (dg): k = Dar⋅735⋅106⋅dg2, where k is the permeability in m2 and Dar is the conversion factor for unit Darcy into m2 (=9.869.10–13).
Oxygen (O2), Hydrogen Sulfide (H2S), and Redox Potential (Eh)
The vertical distribution of O2, H2S, and Eh was measured in sediment cores using a motorized micromanipulator (MMS9083) equipped with OX-100 and H2S-200 microsensor probes, and a redox microelectrode RD-200 coupled with the reference electrode REF-RM (Unisense A/S, Denmark). All microsensor probes were calibrated according to the manufacturer’s instructions. The OX-100 probe was calibrated with a two-point oxic–anoxic calibration, H2S-200 with eight-point calibration in fresh Na2S solutions [1–300 μM in a de-oxygenated calibration buffer (NaAc/HAc, pH < 4)], RD-200 coupled with REF-RM by a two-point calibration using synchronized immersion of electrodes in quinhydrone redox solutions in pH 4 and pH 7 buffers. During measurements, in situ temperature was maintained by keeping the sediment cores in a glass container filled with seawater from the sampling location. From July to October 2017, H2S in pore waters centrifuged from 5 mm sediment cores section was measured according to .
Total Lipids, Fatty Acid Composition, and Sulfur
The detailed procedure for determination of total lipids (TL), fatty acid composition [fatty acid methyl esters (FAME)], and elemental sulfur (S0) was reported previously (). Briefly, lyophilized samples were weighed and consecutively ultrasonically extracted into a solvent mixture of dichloromethane/methanol (DCM: MeOH, 2:1) and separated into layers. Lipid containing layers were evaporated to dryness and weighed. The fatty acid composition was determined from lipid extracts, after saponification, acidification, methylation, and extraction in DCM.
Fatty acid methyl esters were analyzed by an Agilent gas-liquid chromatography (GLC) 6890 N GC System with a 5973 Network Mass Selective Detector, capillary column (30 m × 0.3 mm × 0.25 μm; cross-linked 5% phenylmethylsiloxane) and ultra-high purity helium as the carrier gas. The GLC settings were: programmed column temperature increase from 145 to 215°C by 4°C/min, then by 1°C/min to 225°C, and by 4°C/min to 270°C at constant column pressure of 2.17 kPa. The concentration of elemental sulfur (S0) in FAME chromatograms was estimated using the calibration curve determined under the same GLC settings as FAME.
Data Analyses
Analysis of variance (ANOVA) was performed by GMAV5 for Windows. The effects of factors settlements (mixed and monospecific), summer (2017 and 2018), and sampling months (July, August, and September) on the biomasses of C. cylindracea and macroalgae were tested by three-way ANOVA. The effects of factors summer (2017 and 2018) and sampling months (July, August, and September) on the above-ground biomass, below-ground biomass, shoot density, and the longest leaf length of C. nodosa were tested by two-way ANOVA. In addition, the effects of factors sampling site (mixed, monospecific, and non-vegetated) and summer (2017 and 2018) on total lipid and prokaryotes concentration in the sediment were tested by two-way ANOVA. Homogeneity of variances was tested by Cochran’s test. Significant ANOVA results (p < 0.05) were followed with post hoc Student–Newman–Keuls (SNK) test. Correlations were tested using Pearson’s correlation coefficient (r) with a level of statistical significance, p < 0.05. Only sediment samples originating from the upper part of cores (0–4 cm), where most rhizomes and roots reside, were used for analyses. A multivariate principal component analysis (PCA, Primer 6) was used to recognize the most important variables which explained the differences between vegetated sediments. Correlation matrices were based on H2S, S0, Eh, O2, TL, and PAs as variables, which were normalized due to their different scales. Only the principal components with eigenvalues > 1 were considered.
Results
Environmental Variables
During the summer of 2017, daily means of sea-bottom temperature ranged between 26 and 28°C. During autumn and winter, seawater temperatures decreased. The coldest period was recorded at the end of February 2018, with minima of 8.49 and 8.26°C in the mixed and monospecific settlement, respectively, but it lasted only a few days. After a steep increase from April to June 2018, in August 2018, the temperature reached with large fluctuations a maximum of 29.34°C in the mixed and 28.84°C in the monospecific settlement (Figure 2A). The fluctuations in daily temperatures relative to the monthly averages over the sampling period July/August 2017 were considerably lower (up to 1°C) and more uniform than those obtained for July/August 2018 (up to 3°C) (Figure 2B). The temperature conditions in the mixed settlement were more variable than in the monospecific settlement. Half-hourly values were up to 1.34°C higher during the day and up to 1.12°C lower during the night in the mixed than in the monospecific settlement (Figure 2C).
FIGURE 2
Concentrations of inorganic nutrients (DIN: 1.07–3.40 μM; PO4: 0–0.08 μM; SiO4: 1.24–3.21 μM; and Chl a (0.26–0.85 μg L–1) in Funtana bay were generally low. The abundance of prokaryotes in the ambient water (2.8–8.9 × 105 cell mL–1) correlated with seawater temperature (r = 0.618; p < 0.05) and PM (4.94–15.47 mg L–1) with Chl a (r = 0.49, p < 0.05). The variations in salinity (S: 35.5–38.5) were irregular, with the lowest values between August and October 2018 (Supplementary Table S1).
Biometry
Cymodocea nodosa (Figure 3A) showed a regular seasonal fluctuation in biomass, leaf length, and shoot density which peaked in summer 2017 and 2018 (up to 116.0 ± 27.5 g m–2, 130.6 ± 20.6 mm, and 2043 ± 100 shoots m–2, respectively) and decreased in winter/early spring (down to 17.5 ± 4.0 g m–2, 52.6 ± 23.8 mm, and 42 ± 62 shoots m–2, respectively). In the summer of 2017, the C. nodosa above-ground biomass, the leaves length, and shoot density were not significantly different from those measured during summer 2018 (Supplementary Tables S2, S3). The biomasses of rhizomes and roots reached their maximum during summer 2017 (up to 346.9 ± 32.1 g m–2). After the winter minimum (down to 76.7 ± 45.0 g m–2), the below-ground biomass increased until May 2018 (254.8 ± 33.7 g m–2), then decreased until September 2018 (89.3 ± 16.0 g m–2), and began to increase again in October 2018 (122.3 ± 24.12 g m–2). During summer 2018, the below-ground biomass was significantly lower than in the summer of 2017 (Supplementary Table S2).
FIGURE 3
In the mixed settlement (Figure 3B), C. cylindracea biomass peaked in October 2017 (264.7 ± 65.1 g m–2). After a significant decline in December 2017 (107.7 ± 19.5 g m–2), biomass reached a minimum in May 2018 (18.5 ± 4.5 g m–2). Following the increase in June 2018 (131.5 ± 23.5 g m–2), the biomass moderately increased until October 2018 (177.7 ± 13.2 g m–2). In the monospecific settlement, C. cylindracea biomass rapidly increased in October 2017 (346.3 ± 69.7 g m–2), gradually reaching its maximum in December 2017 (399.2 ± 56.3 g m–2). After the minimum in March 2018 (57.5 ± 32.3 g m–2), the biomass increased in June 2018 (162.7 ± 32.5 g m–2), and remained fairly stable until a considerable increase in October 2018 (301.7 ± 24.5 g m–2). During the summer of 2017, C. cylindracea biomass in the mixed settlement was significantly lower than in the monospecific settlement, while in 2018, no significant difference in biomass between the two settlements was observed (Supplementary Table S4).
In the mixed and monospecific settlements, similar macroalgae taxa were present throughout the study period. The corticated Rytiphlaea tinctoria (Rhodophyta), foliose Udotea petiolata (Chlorophyta), and articulated calcareous Halimeda tuna (Chlorophyta) algae were almost always present in the both settlements (Supplementary Table S5). Filamentous macroalgae (Cladophora sp., Polysiphonia sp.) generally contributed more to the macroalgae cover of the mixed settlement, while their presence with a higher coverage in both settlements was observed during the summer 2018 (Supplementary Figure S1). The macroalgae biomass in the monospecific settlement was significantly lower than in the mixed settlement whereas no differences between two summers were observed (Supplementary Table S4). In both settlements, the highest macroalgae biomass was recorded in October 2017 (mixed: 173.5 ± 14.2 g m–2 and monospecific: 130.7 ± 9.8 g m–2). In the mixed settlement, a macroalgae biomass increase (169.3 ± 15.0 g m–2) was also observed in April 2018 with a presence and higher coverage of filamentous macroalgae, particularly Cladophora prolifera, in contrast to the monospecific settlement (Figure 3C and Supplementary Figure S1).
Total Lipid (TL) Concentrations and Fatty Acid Composition
TL in the C. nodosa leaves (11.5 ± 2.6 to 29.6 ± 3.1 mg g–1 dry weight, DW) showed two maxima, in March and August 2018, while the minimum was measured in November 2017 (Figure 4A). The major fatty acid components were palmitic (C16:0) among the saturated (SAT), oleic (C18:1n-9) in monounsaturated fatty acids (MUFA) and α-linolenic (C18:3 n-3, ALA) in polyunsaturated fatty acids (PUFA). ALA generated a pattern of PUFA and the unsaturation degree (UND) dynamics. PUFA and UND were high (61.5–65.9% and 5.29–6.43, respectively) during the summer of 2017 and decreased toward November 2017 (55.6%, 4.21). From December 2017 onward, PUFA and UND increased and reached a maximum in March 2018 (67.9%, 6.59, respectively). After a decrease in April 2018 and a gradual increase until July 2018, PUFA and UND decreased in August 2018 (59.3%, 4.56, respectively), and then increasing until October 2018 (Figure 4B and Supplementary Table S6).
FIGURE 4
TL extracted from the C. cylindracea from the mixed and monospecific settlement co-varied and was considerably higher in the thalli of the mixed settlement, except in August 2018. Generally, maximum values were reached in February (mixed: 22.5 ± 1.0 mg g–1 DW, monospecific: 19.5 ± 1.1 mg g–1 DW) and minimum in August 2018 (mixed: 7.2 ± 1.5 mg g–1 DW, monospecific: 10.2 ± 1.3 mg g–1 DW, Figure 4A). The percent distributions of dominant components among the fatty acids were seasonally stable. Saturates represented the major fraction of total fatty acids in the mixed (30.1–93.7%) and monospecific (29.7–70.3%) settlements, except for February 2018, and were always dominated by C16:0 and to a lesser extent, by C14:0. MUFA was mostly represented by C16:1(n-7) and C18:1(n-9). PUFA varied inversely to SAT, reaching the maximum in February and minimum in August 2018 in mixed (3.0–55.5%) and monospecific (19.4–54.6%) settlements. PUFA were dominated by C18:3(n-3) with an important contribution of C18:2(n-6), C16:3(n-3), C20:5(n-3), and C20:4(n-6) fatty acids (Supplementary Table S7). UND increased until February 2018, when it started to decline until March to May 2018. In C. cylindracea of both sites, a considerably lower UND were observed in July and August 2018 than during the rest of the study period (Figure 4B).
Sediment
Granulometric Composition
According to the granulometric composition, the vegetated (underlying mixed and monospecific settlements) and non-vegetated sediments were classified as slightly gravelly muddy sand. Given the median grain sizes (dg) and permeability (k), the vegetated sediments were fine-grained (dg < 165 μm) with low permeability (k < 2.10–11 m2), while the non-vegetated sediment was classified as medium-grained (165 μm < dg < 379 μm) and permeable (2.10–11 m2 < k < 10–10 m2). Generally, the dg distribution showed a shift toward coarser particles from the sediment surface until 4 cm and from the sediments of the mixed to monospecific settlements and to the non-vegetated sediment (Figure 5).
FIGURE 5
Oxygen (O2), Redox Potential (Eh), Hydrogen Sulfide (H2S), and Elemental Sulfur (S0)
Oxygen concentrations (O2) in the bottom water of the monospecific settlement (0.1 ± 0.2 to 126.8 ± 30.0 μM) were lower than in the mixed settlement (38.9 ± 25.1 to 195.8 ± 5.9 μM) or above non-vegetated sediment (111.5 ± 4.8 to 198.8 ± 5.1 μM). At all three sites, the variations of oxygen content relative to saturation level (O2/O2 sat) throughout the sampling period were generally similar. In the monospecific settlement, the bottom waters were hypoxic in July and August 2017 and in April, May, and July 2018, while in the mixed settlement in September and October 2018 (Figure 6A). Similar monthly variations in O2 concentrations were observed in the sediment (0–4 cm) at all stations reaching their maximum O2 concentrations during March and April 2018. Lower values and narrower ranges of O2 were found in the sediment below the monospecific settlement (0 to 17.6 ± 24.7 μM) than below mixed settlements (0.5 ± 0.5 to 68.6 ± 97.1 μM) and non-vegetated (1.5 ± 1.3 to 49.9 ± 39.6 μM) sediments.
FIGURE 6
The redox transition depth (RTD, Eh ≤ 0) showed a similar trend in all sediments, being generally shallower below the C. cylindracea monospecific settlement (Figure 6B). In sediments below mixed (8–28 mm), monospecific settlements (4–19 mm), and non-vegetated sediment (14–30 mm), RTDs gradually deepened from July 2017 to March 2018. Afterward it progressed abruptly toward the sediments surface in April to May 2018 followed by moderate fluctuations until September 2018. During summer 2018, in the sediment below the mixed settlement, RTD (5–7 mm) was shallower than during the summer of 2017 (8–9 mm). In contrast, RTDs were similar in both seasons below the C. cylindracea monospecific settlement (2017: 4–7 mm; 2018: 0–7 mm) and in non-vegetated sediment (2017: 15–17 mm; 2018: 18–20 mm). Eh differed between stations, being more negative in the sediment below the monospecific settlement than in the other sediments.
Concentrations of free H2S in sediment pore water generally increased with depth. Below the mixed settlement, the H2S concentrations (up to 29.49 ± 10.18 μM) showed a more pronounced and frequent formation of accumulation zones, between the surface and 2 cm of sediment depth. Such appearances were recorded in December 2017, and in May 2018 with H2S rising to the surface of sediment and intruding into bottom waters in September 2018, when the highest concentrations were measured (up to 90.8 ± 30.5 μM). For H2S concentrations in the monospecific settlement (up to 36.8 ± 10.6 μM), the elevation toward the sediment surface was recorded in April 2018 and then in July 2018 when it appeared in bottom waters (Figure 7). In non-vegetated sediment, H2S concentrations (up to 7.4 ± 3.1 μM) were generally low (data not shown). During summer 2018, in the sediment below the mixed settlement H2S concentrations were higher and accumulation fronts shallower than during summer 2017.
FIGURE 7
The ranges of approximated S0 concentrations in the sediment below the monospecific settlement (0.01–0.23 mg⋅g–1 DW) were higher than below mixed settlement (0.01–0.09 mg⋅g–1 DW) and in the non-vegetated sediment (0–0.03 mg⋅g–1 DW, data not shown). In all sediments, the concentrations of S0 generally varied in opposite mode to H2S concentrations, being lower in H2S accumulation zones.
Prokaryotic Abundance
Prokaryotic abundance in the sediment layers (0–4 cm) underlying the mixed settlement (15.8 ± 3.3 to 33.4 ± 2.7.107 cells g–1 fresh weight, FW) was similar to the sediment below the monospecific settlement (16.8 ± 1.1 to 29.9 ± 4.5.107 cells g–1 FW). Both vegetated sediments displayed a significantly higher abundance than the non-vegetated (12.9 ± 1.6 to 18.8 ± 1.9 .107 cells g–1 FW) sediment (Supplementary Table S8). In all sediments, PA was lower during the winter and spring (February to May 2018). After an increase in June 2018, high PA was maintained throughout the summer of 2018 in all sediments, being significantly higher than during the summer of 2017 (Supplementary Table S8 and Supplementary Figure S2A).
Total Lipids (TL)
The concentrations of total lipids in the surface sediments below the mixed (1.15 ± 0.08 to 1.62 ± 0.10 mg g–1 DW), monospecific (1.00 ± 0.07 to 1.22 ± 0.09 mg g–1 DW) settlements and non-vegetated sediments (0.38 ± 0.04 to 0.55 ± 0.14 mg g–1 DW) differed significantly, being ranked as mixed > monospecific > non-vegetated. During the summer of 2018, TL in the sediments at all sites were significantly lower in comparison to the summer of 2017 (Supplementary Table S8 and Supplementary Figure S2B).
Relationship Between Environmental Parameters in Vegetated Sediments
The relationships between H2S, Eh, O2, S0, TL, and PA in sediments underlying mixed and monospecific settlements were analyzed using a PCA. PC1 explained 51.1% and PC2 23.1% of the variability (Figure 8). PC1 separated most of the sediments below the mixed settlement (except the sediments from June to September 2018) according to higher concentrations of O2, TL, and more positive Eh from most of the sediments below the monospecific settlement (except the sediments in February and March 2018). The monospecific settlement was characterized by higher PA and S0 concentrations. On PC2, more negative Eh, higher H2S, and lower S0 concentrations caused the separation of sediments below the mixed (May to September 2018) and monospecific settlement (April and July 2018) from all other samples (Figure 8).
FIGURE 8
Discussion
C. nodosa in Control of the Growth of C. cylindracea
The presence of C. cylindracea in the meadow did not affect the shoot density and the above-ground biomass of C. nodosa throughout the study period. Invaded meadow of C. nodosa showed consistent seasonal fluctuations reaching the highest above-ground biomass in the summers of 2017 and 2018, the lowest in winter, and exhibited exponential growth in spring as commonly observed for this species in monospecific meadows (; ; ). Likewise, the seasonal cycle of C. cylindracea in both settlements, with increasing biomass from summer to autumn, a drastic regression in winter and a renewed growth in spring occurred as reported for other temperate areas (; ; ; ). However, our results showed that during 2017 C. cylindracea biomass in the mixed settlement was lower than in its monospecific stands. This finding could be ascribed to the shading provided by C. nodosa leaves and to the reduced availability of open substrate for C. cylindracea rhizoid attachment in the C. nodosa meadow. The contribution of the macroalgae to the effects of shading and open substrate reduction on C. cylindracea growth was larger in the mixed settlement, most probably due to the significantly higher biomass in comparison to the monospecific settlement.
Shading by reduction of incident light attenuates or even prevents photosynthesis, retarding the growth and decreasing the biomass production of autotrophs (). As a physiological response to a shaded habitat, C. cylindracea in the mixed settlement generally accumulated more lipids and showed a higher unsaturation of the thalli in comparison to the monospecific one predominantly during the intensive meadow growth. By increased unsaturation, membrane fluidity is increased and the electron transport in the photosystem is facilitated to approach optimal photosynthetic activity under low light conditions (; ; ).
Caulerpa cylindracea, like other rhizophytic algae, tends to bury their rhizoids in the sediment for nutrient acquisition, thus gaining an advantage over non-rhizophytic macroalgae relying only on the nutrients from the water column. In contrast with other rooted primary producers, such as seagrasses, C. cylindracea competes for nutrients from the water and from the sediment. The ability of Caulerpa to take up and translocate nutrients from the sediment supports its biomass production () and spreading over the bare sediment () and assists the growth and development of the monospecific settlement. This nutrient uptake ability of Caulerpa is critical since low-nutrient conditions generally characterize the coastal waters of western Istria (). In contrast, C. cylindracea in the mixed settlement, in addition to competition for nutrients, faced a reduced availability of free sediment for rhizoid burial due to the overwhelming presence of C. nodosa below-ground tissue. The reduced open substrate and therefore insufficient support with nutrients from the sediment most likely caused a decrease in biomass after substantial proliferation in October 2017, in contrast to the C. cylindracea settlement where high biomass was maintained until December 2017. Previous research showed that in healthy growing seagrass meadows, invasive macroalgae, regardless of their biomass, have little effect on sediment biogeochemistry or habitat modification (). Thus, while the resistance to the invasion of intact seagrass beds, such as those of Posidonia oceanica, has primarily been attributed to light limitation (; ), the success of the invasion might also be controlled by the limited availability of sediment substrate. Therefore, the lack of open substrate was crucial factor in limiting C. cylindracea growth within the meadow during 2017, since the shading provided by C. nodosa leaves and macroalgae continually affected this understory species during the entire study period.
The sediments underlying the mixed settlements showed a higher organic enrichment than those below the monospecific stands, suggesting that the combined effect of both macrophytes increased the sediment organic pools and enhanced the trapping of organic matter (). Their different granulometric composition additionally emphasizes this divergence in organic matter enrichment between the two sediments. The absorptive capacity for organic molecules and the degree of organic matter preservation was higher in the mixed area due to the smaller grain size of the sediment. Generally higher H2S concentrations and shallower depths of the sulfide front were found below the invaded meadow in comparison to the sediment below C. cylindracea monospecific stands (Figure 7). However, in the monospecific settlement, lower oxygen concentrations, more negative Eh coupled with RTD closer to the sediment surface were observed. The oxygen depletion in C. cylindracea sediment might be explained by a considerably higher S0 concentration likely derived from more intense sulfide oxidation. This process could indicate a rapid and efficient removal of the produced sulfide by the chemical reaction with free oxygen or more likely mediated by sulfide-oxidizing bacteria (; ). Additionally, the presence of densely packed thalli covering the sediment surface limited the oxygen diffusion and penetration into the sediment. By the simultaneous consumption of H2S and O2 highly reducing conditions near the sediment surface are created.
In contrast, in the sediment under the mixed settlement, the deeper RTD might be related to a below-ground tissue due to the vertical coupling between its biomass and positive redox anomaly (). A more extensive sediment layer with oxic conditions (Eh > 0) accompanied with a more positive redox potential was probably caused by the better oxygenation of the sediment due to radial oxygen leakage from the C. nodosa below-ground tissue along with O2 diffusion from bottom waters in contact with the meadow (). The difference between the RTDs of the monospecific and mixed settlements was more pronounced until December 2017 (Figure 6) when both macrophytes were declining in biomass. During the winter/early spring, H2S production decreased in both settlements, likely due to the reduced activity of sulfate-reducing prokaryotes at lower temperatures, and the sediments gradually shifting toward a more oxidized state. H2S detected within the rooted area of the invaded meadow could be attributed to the presence of reducing micro-niches where anaerobic metabolism could occur irrespective of surrounding redox conditions (; ). By rapid desaturation of increasing lipids, both macrophytes increased their cell membranes fluidity to achieve chilling resistance at lower temperatures (). In both settlements, C. cylindracea survived winter temperatures in the form of small thallus fragments (rhizoids, stolons, and propagules). Similar and less drastic reductions in biomass were also reported in warmer regions of the Mediterranean Sea (; ). Thus, the mild winter with short-lasting low temperatures (below 9°C) likely allowed a faster recovery and growth of C. cylindracea during the spring ().
The bloom of macroalgae within the invaded C. nodosa meadow in April 2018 (Figure 3C) most likely contributed to prolong the regression of C. cylindracea until May 2018, in contrast to the monospecific settlement where its growing phase started earlier (Figure 3B). Filamentous macroalgae, which in April 2018 contributed noticeably to the meadow cover (Supplementary Figure S1), generally grow better than other morpho-functional groups due to their higher nutrient uptake rates (; ). These macroalgae likely exploited this advantage in the competition for nutrients with C. cylindracea, that was still in dormant phase when the metabolic needs are at a minimum, and monopolized the available substrate. This short-term burst in macroalgal biomass followed by a steep decrease in biomass (Figure 3C) suggests that at least a part of this material became available for degradation.
C. cylindracea Proliferation Led to the Alteration of Sediment Conditions and Disturbance of the Stability of the C. nodosa Meadow
In April and May 2018, C. nodosa entered its growth phase. The shoot density along with below-ground tissue considerably increased. The sediment nutrient uptake by small seagrass species, such as C. nodosa, is closely dependent on nutrients regenerated from recently deposited organic matter (). In May 2018, an increase in temperature and availability of the macroalgal remains probably stimulated the activity of sulfate-reducing prokaryotes and the release of nutrients along with H2S production leading to the development of more reducing conditions. From June onward, the conditions in the invaded meadow sediment shifted to a new regime, characterized by expanding anoxia closer to the sediment surface, that become similar to that under the C. cylindracea monospecific settlement (Figures 6, 8). Such highly reducing sediment conditions had been likely stimulated by a production of new and very active rhizoid clusters buried in the sediments. Already in June 2018, C. cylindracea proliferated in both settlements and their biomasses became very similar. C. nodosa responded to the newly established sediment conditions by decreasing the below-ground tissue biomass while shoot density and above-ground biomass continued to increase. It seems very likely that anoxia driven by C. cylindracea proliferation severely reduced or even prevented the translocation of photoassimilates and the transport of photosynthetically produced oxygen from leaves and shoots to the rhizomes and roots (; ). In accordance, the energy production in the roots was impaired (). Since C. nodosa cannot tolerate such long-lasting anoxia without reserves imported from the above-ground tissue, the reduction of below-ground biomass is likely unavoidable. Such unbalanced growth in which C. nodosa kept a more substantial fraction of CO2 fixation in the photosynthetic compartment and was not stored as reserves lasted until October 2018 (Figure 3A). This resulted in considerably higher above- to below-ground biomass ratios than observed during the summer of 2017 and reported for healthy growing and stable C. nodosa meadows (; ). Therefore, the control of C. cylindracea over C. nodosa was supported by its alteration of the below-ground processes, as Caulerpa species are regularly abundant in disturbed seagrass beds and their decaying patches (; ; ; ).
During July and August 2018 C. cylindracea likely underwent spawning. The release of gametes, which has a high energetic cost, resulted in degraded thalli condition manifested by discoloration. Such conditions were evidenced by a considerably decreased UND of thalli (Figure 4B), which derived from the sharp depletion of PUFAs, most likely due to their incorporation into gametes enabling their long-term viability (). The spawning was likely triggered by large fluctuations in daily temperatures that C. cylindracea faced in July and August 2018 unlike during the same period in 2017 (Figure 2B) when only vegetative propagation occurred. However, the greater decrease in physiological indices of C. cylindracea in the meadow, matching to its poorer condition, suggests that a comparatively larger part of the population spawned than in the monospecific settlement. This difference could be ascribed to a wider range of day–night temperature oscillations (Figure 2C) combined with a milder hydrodynamic regime due to the flow reduction of the waters by the seagrass leaves, which is more favorable for spawning algae (). During the spawning, no decrease in C. cylindracea biomass was observed in either settlement. However, it is possible that energy invested in this process, especially in the mixed settlement became unavailable for growth in early autumn in contrast to the previous year. Unexpectedly, macroalgae displayed a lack of biomass increase in both settlements. While the temperature and its large oscillations might have also influenced their growth cycle, the impaired physiological condition of the C. cylindracea might have contributed as well, but the exact mechanism remains unclear. However, the resulting poor competitive ability of C. cylindracea due to exhaustion by spawning was likely exploited by filamentous macroalgae as their contribution in the cover of both settlements was strikingly higher than in the previous summer (Supplementary Figure S1). The temperature fluctuations did not have a major effect on the conditions of C. nodosa above-ground tissues. Only a moderate decrease in the degree of unsaturation was observed (Figure 4B), presenting a typical metabolic adjustment to a sudden increase in temperature (up to 29.34°C) in August 2018. A similar response of seagrass was observed during an experiment where C. nodosa and P. oceanica were exposed to a heat wave reaching a plateau of 29°C () which is comparable to our field conditions. C. nodosa tolerates heat stress by sustaining enhanced photosynthesis, an ability allowing its growth and persistence in shallow environments (; ).
In September 2018, the minimized or even disrupted root respiration of the reduced below-ground tissue biomass likely enabled a H2S penetration into the plant, causing a strong shedding of the shoots. During October 2018, the conditions in the water column and sediment returned to a sustainable state favorable for the recovery of rhizomes. The slight increase in shoot density indicated that C. nodosa might have regained its fitness likely due to the increased below-ground biomass and the re-establishment of more balanced growth. This revitalization accompanied with the lack of invasive macroalgae growth might have increased the stability of the meadow, thereby re-exerting its control over C. cylindracea growth and propagation.
Ecological Significance
Our results showed that the changes in the sediment conditions, caused by invasive C. cylindracea proliferation, appear to be the main factor that affected the decrease in C. nodosa below-ground biomass production, thus disrupting the stability of the seagrass meadow. The presence of C. cylindracea did not affect the C. nodosa above-ground biomass and shoot density. The invaded seagrass meadow appears to keep or even enhance its trapping of particles and thereby, increasing water clarity and providing protection against coastal erosion. However, the presence of C. cylindracea might change the biodiversity and secondary production of foraging benthic species since the invaded meadow does not provide the same trophic and nursery support as the monospecific C. nodosa habitats (; ). Therefore, changes in the structural complexity of the coastal ecosystem where these two macrophytes coexist could be expected.
Even though C. nodosa is a protected species () and an important biological element of coastal water quality included in European directives and monitoring programs (), its invaded meadows should receive more attention. The prolonged destabilization of the invaded meadows could increase their susceptibility to coastal modification and water quality degradation which are considered the main causes of C. nodosa decline and extinction as reported in European coastal areas including the northern Adriatic (; ; ; ). So far, containment of C. cylindracea was not successful and further spread along the Adriatic coast can be expected. For these reasons, invaded meadows should be considered as more vulnerable habitats and given a higher monitoring priority. In particular, the results of our study present an early warning to our local/regional authorities to implement policies and strategies to prevent further coastal destruction (e.g., touristic facilities, mechanical damage by dredging and anchoring) and mitigate the impact of existing anthropogenic activities (e.g., untreated and uncontrolled release of wastewaters) negatively affecting the quality of coastal waters. Further studies of C. nodosa and C. cylindracea interactions at larger spatial and temporal scales are needed to predict and evaluate the stability, conservation, and sustainability of the functions and services provided by invaded C. nodosa meadows to the shallow coastal areas.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author/s.
Author contributions
MN, MK, and GH contributed to conceptualization. MK, PP, MM, II, LI, IF, and MN contributed to investigation. MN contributed to formal analysis and writing—original draft. MK, GH, PP, LI, II, IF, and MM contributed to writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
The financial support was provided by the Croatian Science Foundation to MN (project IP-2016-06-7118, MICRO-SEAGRASS). GH was supported by grants of the Austrian Science Fund (ARTEMIS; #P28781-B21 and Microbial Nitrogen Cycling; # W1257-B20).
Acknowledgments
We sincerely thank J. Jakovčević and M. Buterer for nutrient and chlorophyll a determination, and A. Budiša and I. Haberle for occasional help during separation and biometry of plant material.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.602055/full#supplementary-material
References
1
AgostiniS.PergentG.MarchandB. (2003). Growth and primary production of Cymodocea nodosa in a coastal lagoon.Aquat. Bot.76185–193. 10.1016/s0304-3770(03)00049-4
2
AlcoverroT.ZimmermanR. C.KohrsD. G.AlberteR. S. (1999). Resource allocation and sucrose mobilization in light-limited eelgrass Zostera marina.Mar. Ecol. Prog. Ser.187121–131. 10.3354/meps187121
3
AlexandreA.SantosR. (2020). High nitrogen and phosphorous acquisition by belowground parts of Caulerpa prolifera (Chlorophyta) contribute to the species’ rapid spread in Ria Formosa lagoon, Southern Portugal.J. Phycol.56608–617. 10.1111/jpy.12988
4
Beca-CarreteroP.GuihéneufF.Marín-GuiraoL.Bernardeau-EstellerJ.García-MuñozR.StengelD. B.et al (2018). Effects of an experimental heat wave on fatty acid composition in two Mediterranean seagrass species.Mar. Pollut. Bull.13427–37. 10.1016/j.marpolbul.2017.12.057
5
Beca-CarreteroP.GuihéneufF.WintersG.StengelD. B. (2019). Depth-induced adjustment of fatty acid and pigment composition suggests high biochemical plasticity in the tropical seagrass Halophila stipulacea.Mar. Ecol. Prog. Ser.608105–117. 10.3354/meps12816
6
BlažinaM.IvešaL. J.NajdekM. (2009). Caulerpa racemosa: adaptive varieties studied by fatty acid composition (Northern Adriatic Sea, Vrsar, Croatia).Eur. J. Phycol.44183–189. 10.1080/09670260802428250
7
BoudouresqueC. F.BernardG.PergentG.ShiliA.VerlaqueM. (2009). Regression of Mediterranean seagrasses caused by natural processes and anthropogenic disturbances and stress: a critical review.Bot. Mar.52395–418. 10.1515/bot.2009.057
8
BulleriF.BalataD.BertocciI.TamburelloL.Benedetti-CecchiL. (2010). The seaweed Caulerpa racemosa on Mediterranean rocky reefs: from passenger to driver of ecological change.Ecology912205–2212. 10.1890/09-1857.1
9
CebriánJ.PedersenM. F.KroegerK. D.ValielaI. (2000). Fate of production of the seagrass Cymodocea nodosa in different stages of meadow formation.Mar. Ecol. Prog. Ser.204119–130. 10.3354/meps204119
10
CeccherelliG.CinelliF. (1997). Short-term effects of nutrient enrichment of the sediment and interactions between the seagrass Cymodocea nodosa and the introduced green alga Caulerpa taxifolia in a Mediterranean bay.J. Exp. Mar. Biol. Ecol.217165–177. 10.1016/s0022-0981(97)00050-6
11
CeccherelliG.CinelliF. (1999). Effects of Posidonia oceanica canopy on Caulerpa taxifolia size in a north-western Mediterranean bay.J. Exp. Mar. Biol. Ecol.24019–36. 10.1016/s0022-0981(99)00044-1
12
CeccherelliG.PiazziL.CinelliF. (2000). Response of the non-indigenous Caulerpa racemosa (Forsskål) J. Agardh to the native seagrass Posidonia oceanica (L.) Delile: effect of density of shoots and orientation of edges of meadows.J. Exp. Mar. Biol. Ecol.243227–240. 10.1016/s0022-0981(99)00122-7
13
CeccherelliG.PinnaS.CussedduV.BulleriF. (2014). The role of disturbance in promoting the spread of the invasive seaweed Caulerpa racemosa in seagrass meadows.Biol. Inv.162737–2745. 10.1007/s10530-014-0700-7
14
ChisholmJ. R. M.MoulinP. (2003). Stimulation of nitrogen fixation in refractory organic sediments by Caulerpa taxifolia (Chlorophyta).Limnol. Oceanogr.48787–794. 10.4319/lo.2003.48.2.0787
15
CliftonK. E. (1997). Mass spawning by green algae on coral reefs.Science2751116–1118. 10.1126/science.275.5303.1116
16
ClineJ. D. (1969). Spectrophotometric determination of hydrogen sulfide in natural waters.Limnol. Oceanogr.14454–458. 10.4319/lo.1969.14.3.0454
17
CormaciM.FurnariG.GiacconeG. (2004). “Macrophytobenthos,” inMediterranean Marine Benthos: A Manual of Methods for Its Sampling and Study, edsGambiM. C.DappianoM. (Genova: Società Italiana di Biologia Marina), 217–246.
18
CúcioC.EngelenA. H.CostaR.MuyzerG. (2016). Rhizosphere microbiomes of European seagrasses are selected by the plant, but are not species specific.Front. Microbiol.7:440. 10.3389/fmicb.2016.00440
19
de los SantosC. B.Krause-JensenD.AlcoverroT.MarbàN.DuarteC. M.van KatwijkM. M.et al (2019). Recent trend reversal for declining European seagrass meadows.Nat. Commun.10:3356. 10.1038/s41167-019-11340-4
20
De VillèleX.VerlaqueM. (1995). Changes and degradation in Posidonia oceanica bed invaded by the introduced tropical alga Caulerpa taxifolia in the northwestern Mediterranean.Bot. Mar.3879–87.
21
DuarteC. M.KennedyH.MarbàN.GaciaE.FourqureanJ. W.BegginsJ.et al (2013). Seagrass community metabolism: assessing the capacity of seagrass meadows for carbon burial: current limitations and future strategies.Ocean Coast. Manag.8332–38. 10.1016/j.ocecoaman.2011.09.001
22
DuarteC. M.MerinoM.AgawinN. S. R.UriJ.FortesM. D.GallegosM. E.et al (1998). Root production and belowground seagrass biomass.Mar. Ecol. Prog. Ser.17197–108. 10.3354/meps171097
23
EnríquezS.MarbàN.DuarteC. M.van TussenbroekB. I.Reyes-ZavalaG. (2001). Effects of seagrass Thalassia testudinum on sediment redox.Mar. Ecol. Prog. Ser.219149–158. 10.3354/meps219149
24
EpsteinS. S.RosselJ. (1995). Enumeration of sandy sediment bacteria: search for optimal protocol.Mar. Ecol. Prog. Ser.117289–298. 10.3354/meps117289
25
FolkR. L. (1954). The distinction between grain size and mineral composition in sedimentary-rock nomenclature.J. Geol.62344–359. 10.1086/626171
26
FrederiksenM. S.GludR. N. (2006). Oxygen dynamics in the rhizosphere of Zostera marina: a two-dimensional planar optode study.Limnol. Oceanogr.511072–1083. 10.4319/lo.2006.51.2.1072
27
GallucciF.HutchingsP.GribbenP.FonsecaG. (2012). Habitat alteration and community-level effects of an invasive ecosystem engineer: a case study along the coast of NSW, Australia.Mar. Ecol. Prog. Ser.44995–108. 10.3354/meps09547
28
GambiM. C.NowellA. R. M.JumarsP. A. (1990). Flume observations on flow dynamics in Zostera marina (eelgrass) beds.Mar. Ecol. Prog. Ser.61159–169. 10.3354/meps061159
29
GangiA. F. (1985). Permeability of unconsolidated sands and porous rocks.J. Geophys. Res. Solid903099–3104. 10.1029/jb090ib04p03099
30
Garcias-BonetN.MarbàN.HolmerM.DuarteC. M. (2008). Effects of sediment sulfides on seagrass Posidonia oceanica meristematic activity.Mar. Ecol. Prog. Ser.3721–6. 10.3354/meps07714
31
GlasbyT. M. (2013). Caulerpa taxifolia in seagrass meadows: killer or opportunistic weed?Biol. Inv.151017–1035. 10.1007/s10530-012-0347-1
32
GribbenP. E.ThomasT.PuscedduA.BonechiL.BianchelliS.BuschiE.et al (2018). Below-ground processes control the success of an invasive seaweed.J. Ecol.1062082–2095. 10.1111/1365-2745.12966
33
GuschinaI. A.HarwoodJ. (2009). “Algal lipids and effect of environment on their biochemistry,” inLipids in Aquatic Ecosystems, edsArtsM.BrettM.KainzM. (New York, NY: Springer), 1–24. 10.1007/978-0-387-89366-2_1
34
HolmerM.MarbàN.LamoteM.DuarteC. M. (2009). Deterioration of sediment quality in seagrass meadows (Posidonia oceanica) Invaded by Macroalgae (Caulerpa sp.).Estuar. Coast.32456–466. 10.1007/s12237-009-9133-4
35
HolmerM.PedersenO.IkejimaK. (2006). Sulfur cycling and sulfide intrusion in mixed Southeast Asian tropical seagrass meadows.Bot. Mar.4991–102. 10.1515/bot.2006.013
36
Holm-HansenO.LorenzenC. J.HolmesR. W.StricklandJ. D. H. (1965). Fluorometric determination of chlorophyll.J. Conseil.3013–15. 10.1093/icesjms/30.1.3
37
IvančićI.PaliagaP.PfannkuchenM.ÐakovacT.NajdekM.SteinerP.et al (2018). Seasonal variations in extracellular enzymatic activity in marine snow-associated microbial communities and their impact on the surrounding water.FEMS Microbiol. Ecol.94:fyi198. 10.1093/femsec/fiy198
38
IvešaL.BlažinaM.NajdekM. (2004). Seasonal variation in fatty acid composition of Caulerpa taxifolia (M. Vahl.) C. Ag. in the northern Adriatic Sea (Malinska, Croatia).Bot. Mar.47209–214. 10.1515/BOT.2004.021
39
IvešaL.DevescoviM. (2006). Seasonal vegetation patterns of the introduced Caulerpa racemosa (Caulerpales, Chlorophyta) in the northern Adriatic Sea (Vrsar, Croatia).Period. Biol.108111–116. 10.1127/nova.hedwigia/50/1990/111
40
IvešaL.DjakovacT.DevescoviM. (2015). Spreading patterns of the invasive Caulerpa cylindracea Sonder along the west Istrian Coast (northern Adriatic Sea, Croatia).Mar. Environ. Res.1071–7. 10.1016/j.marenvres.2015.03.008
41
JørgensenB. B. (1977). The sulfur cycle of a coastal marine sediment (Limfjorden, Denmark).Limnol. Oceanogr.22814–832. 10.4319/lo.1977.22.5.0814
42
KhotimchenkoS. V. (1995). Fatty acid composition of green algae of the genus Caulerpa.Bot. Mar.38509–512. 10.1515/botm.1995.38.1-6.509
43
LittlerM. M.LittlerD. S. (1984). Relationships between macroalgal functional form groups and substrata stability in a subtropical rocky-intertidal system.J. Exp. Mar. Biol. Ecol.7413–34. 10.1016/0022-0981(84)90035-2
44
Marín-GuiraoL.Bernardeau-EstellerJ.RuizJ. M.Sandoval-GilJ. M. (2015). Resistance of Posidonia oceanica seagrass meadows to the spread of the introduced green alga Caulerpa cylindracea: assessment of the role of light.Biol. Inv.171989–2009. 10.1007/s10530-015-0852-0
45
Marín-GuiraoL.RuizJ. M.DattoloE.Garcia-MunozR.ProcacciniG. (2016). Physiological and molecular evidence of differential short-term heat tolerance in Mediterranean seagrasses.Sci. Rep.6:28615. 10.1038/srep28615
46
McKinnonJ. G.GribbenP. E.DavisA. R.JolleyD. F.WrightJ. T. (2009). Differences in soft-sediment macrobenthic assemblages invaded by Caulerpa taxifolia compared to uninvaded habitats.Mar. Ecol. Prog. Ser.38059–71. 10.3354/meps07926
47
MICROMERITICS (2002). SediGraph 5100 Particle Size Analysis System Operator’ Manual.Norcross: Micromeritics Instrument Corporation.
48
NajdekM.KorlevićM.PaliagaP.MarkovskiM.IvančićI.IvešaL.et al (2020). Dynamics of environmental conditions during a decline of a Cymodocea nodosa meadow.Biogeosciences173299–3315. 10.5194/bg-17-3299-2020
49
NguyenH. M.KimM.RalphP. J.Marín-GuiraoL.PerniceM.ProcacciniG. (2020). Stress memory in seagrasses: first insight into the effects of thermal priming and the role of epigenetic modifications.Front. Plant Sci.11:494. 10.3389/fpls.2020.00494
50
Orlando-BonacaM.FrancéJ.MavričB.GregoM.LipejL.Flander PutrleV.et al (2015). A new index (MediSkew) for the assessment of the Cymodocea nodosa (Ucria) Ascherson meadow’s status.Mar. Environ. Res.110132–141. 10.1016/j.marenvres.2015.08.009
51
Orlando-BonacaM.FrancéJ.MavričB.LipejL. (2019). Impact of the port of Koper on Cymodocea nodosa meadow.Annales29187–194.
52
OrthR. J.CarruthersT. J. B.DennisonW. C.DuarteC. M.FourqureanJ. W.HeckK. L.Jr.et al (2006). A global crisis for seagrass ecosystems.BioScience56987–793.
53
PanayotidisP.ŽuljevićA. (2001). Sexual reproduction of the invasive green alga Caulerpa racemosa var. occidentalis in the Mediterranean Sea.Oceanol. Acta24199–203. 10.1016/s0399-1784(01)01142-2
54
Pardal-SouzaA. L.DiasG. M.JenkinsS. R.CiottiÁM.ChristofolettiR. A. (2017). Shading impacts by coastal infrastructure on biological communities from subtropical rocky shores.J. Appl. Ecol.54826–835. 10.1111/1365-2664.12811
55
PearsonG. A.SerraoE. A. (2006). Revisiting synchronous gamete release by fucoid algae in the intertidal zone: fertilization success and beyond?Int. Comp. Biol.46587–597. 10.1093/icb/ic1030
56
PedersenM. F.BorumJ. (1997). Nutrient control of estuarine macroalgae: growth strategy and the balance between nitrogen requirements and uptake.Mar. Ecol. Prog. Ser.161155–163. 10.3354/meps161155
57
PiazziL.BalataD. (2009). Invasion of alien macroalgae in different Mediterranean habitats.Biol. Inv.11193–204. 10.1007/s10530-008-9224-3
58
PiazziL.BalataD.ForesiL.CristaudoC.CinelliF. (2007). Sediment as a constituent of Mediterranean benthic communities dominated by Caulerpa racemosa var. cylindracea.Sci. Mar.71129–135. 10.3989/scimar.2007.71n1129
59
PiazziL.CeccherelliG.CinelliF. (2001). Threat to macroalgal diversity: effects of the introduced green alga Caulerpa racemosa in the Mediterranean.Mar. Ecol. Prog. Ser.210149–159. 10.3354/meps210149
60
PiazziL.CinelliF. (1999). Development and seasonal dynamics of a population of the tropical alga Caulerpa racemosa (Forsskâl) J. Agardh in the Mediterranean.Cryptogamie Algol.20295–300. 10.1016/j.aquabot.2005.02.008
61
PorterK. G.FeigY. S. (1980). The use of DAPI for identification and counting aquatic microflora.Limnol. Oceanogr.25943–984. 10.4319/lo.1980.25.5.0943
62
RasheedM. A.UnsworthR. K. F. (2011). Long-term climate-associated dynamics of a tropical seagrass meadow: implications for the future.Mar. Ecol. Prog. Ser.42293–103. 10.3354/meps08925
63
RizzoL.PuscedduA.StabiliL.AlifanoP.FraschettiS. (2017). Potential effects of invasive seaweed (Caulerpa cylindracea, Sonder) on sedimentary organic matter and microbial metabolic activities.Sci. Rep.7:12113. 10.1038/s41598-017-12556-4
64
Samper-VillarrealJ.LovelockC. E.SaundersM. I.RoelfsemaC.MumbyP. J. (2016). Organic carbon in seagrass sediment is influenced by seagrass canopy complexity, turbidity, wave height, and water depth.Limnol. Oceanogr.61938–952. 10.1002/lno.10262
65
SmithR. D.PregnallA. M.AlberteR. S. (1988). Effects of anaerobiosis on root metabolism of Zostera marina (eelgrass): implications for survival in reducing sediments.Mar. Biol.98131–141. 10.1007/bf00392668
66
StricklandJ. D. H.ParsonsT. R. (1972). A practical handbook of seawater analysis.Bull. Fish. Res. Board. Can.1671–310.
67
TamburelloL.MaggiE.Benedetti-CecchiL.BellistriG.RattrayA. J.RavaglioliC.et al (2015). Variation in the impact of non-native seaweeds along gradients of habitat degradation: a meta-analysis and an experimental test.Oikos1241121–1131. 10.1111/oik.02197
68
TaylorR. B.PeekJ. T. A.ReesT. A. V. (1998). Scaling of ammonium uptake by seaweeds to surface area:volume ratio: geographical variation and the role of uptake by passive diffusion.Mar. Ecol. Prog. Ser.169143–148. 10.3354/meps169143
69
TerradosJ.RosJ. D. (1992). Growth and primary production of Cymodocea nodosa (Ucria) Ascherson in a mediterranean coastal lagoon: the Mar Menor (SE Spain).Aquat. Bot.4363–74. 10.1016/0304-3770(92)90014-A
70
TuyaF.Ribeiro-LeiteL.Arto-CuestaN.CocaJ.HarounR.EspinoF. (2014). Decadal changes in the structure of Cymodocea nodosa seagrass meadows: natural vs. human influences.Estuar. Coast. Shelf Sci.13741–49. 10.1016/j.ecss.2013.11.026
71
UNEP-PAM-RAC/SPA (2012). Protocol Concerning Specially Protected Areas and Biological Diversity in the Mediterranean. Annex II. List of Endangered or Threatened Species.Paris: UNEP.
72
Vaquer-SunyerR.DuarteC. M. (2008). Thresholds of hypoxia for marine biodiversity.PNAS10515452–15457. 10.1073/pnas.0803833105
73
WackerA.PiephoM.HarwoodJ. L.GuschinaI. A.ArtsM. T. (2016). Light-induced changes in fatty acid profiles of specific lipid classes in several freshwater phytoplankton species.Front. Plant Sci.7:264. 10.3389/fpls.2016.00264
74
WilliamsS. L. (1984). Uptake of sediment ammonium and translocation in a marine green macroalga Caulerpa cupressoides.Limnol. Oceanogr.29374–379. 10.4319/lo.1984.29.2.0374
75
WilliamsS. L. (2007). Introduced species in seagrass ecosystems: status and concerns.J. Exp. Mar. Biol. Ecol.35089–110. 10.1016/j.jembe.2007.05.032
76
ZavodnikN.TraviziA.De RosaS. (1998). Seasonal variations in the rate of photosynthetic activity and chemical composition of the seagrass Cymodocea nodosa (Ucr.) Asch.Sci. Mar.62301–309. 10.3989/scimar.1998.62n4301
77
ZimmermanR. C.AlberteR. S. (1996). Effect of light/dark transition on carbon translocation in eelgrass Zostera marina seedlings.Mar. Ecol. Prog. Ser.136305–309. 10.3354/meps136305
78
ŽuljevićA.AntolićB.OnofriV. (2003). First record of Caulerpa racemosa (Caulerpales: Chlorophyta) in the Adriatic Sea.J. Mar. Biol. Assoc. UK83711–712. 10.1017/s0025315403007689h
Summary
Keywords
Caulerpa cylindracea, Cymodocea nodosa, invaded seagrass meadow, sediment biogeochemistry, northern Adriatic Sea, seagrass meadow stability
Citation
Najdek M, Korlević M, Paliaga P, Markovski M, Ivančić I, Iveša L, Felja I and Herndl GJ (2020) Effects of the Invasion of Caulerpa cylindracea in a Cymodocea nodosa Meadow in the Northern Adriatic Sea. Front. Mar. Sci. 7:602055. doi: 10.3389/fmars.2020.602055
Received
02 September 2020
Accepted
26 November 2020
Published
23 December 2020
Volume
7 - 2020
Edited by
Johan Schijf, University of Maryland Center for Environmental Science (UMCES), United States
Reviewed by
Zhijian Jiang, Chinese Academy of Sciences, China; Gloria Peralta, University of Cádiz, Spain; Jaime Bernardeau, Instituto Español de Oceanográfia, Spain
Updates
Copyright
© 2020 Najdek, Korlević, Paliaga, Markovski, Ivančić, Iveša, Felja and Herndl.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mirjana Najdek, najdek@cim.irb.hr
This article was submitted to Marine Biogeochemistry, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.