Skip to main content

Low impact of Zostera marina meadows on sediment and water microbiota under brackish conditions

Abstract

Background

Zostera marina is an important ecosystem engineer influencing shallow water environments and possibly shaping the microbiota in surrounding sediments and water. Z. marina is typically found in marine systems, but it can also proliferate under brackish conditions. Changes in salinity generally have a strong impact on the biota, especially at the salty divide between salinity 6 and 9. To better understand the impact of the salty divide on the interaction between Z. marina and the surrounding sediment and water microbiota, we investigated the effects of Z. marina meadows on the surrounding microbiota across a salinity range of 6–15 in the Baltic Sea during the summer using 16S and 18S rRNA gene amplicon sequencing.

Results

Salinity was the most important factor for structuring the microbiota within both water and sediment. The presence of Z. marina affected the composition of the bacterial and eukaryotic community and bacterial alpha diversity in the sediment. However, this effect was confined to alpha-mesohaline conditions (salinity 9–15). The impact of Z. marina below salinity 9 on water and sediment microbiota was insignificant.

Conclusions

Increasing salinity was associated with a longer leaf length of Z. marina, causing an increased canopy height, which affects the sediment microbiota through reduced water velocity. Hence, we propose that the canopy effect may be the major predictor explaining Z. marina’s interactions with the surrounding microbiota at salinity 9–15. These findings emphasize the importance of the physical effects of Z. marina meadow ecosystem services and have important implications for Z. marina management under brackish conditions in a changing climate.

Introduction

In shallow, sheltered coastal regions, Zostera marina (also called eelgrass) builds large meadows that provide ecosystem services, including sediment stabilization, food provision, shelter, and nutrient recycling. It maintains intimate ecological interactions with microbial consortia that live in association with plants and within the surrounding seawater and sediments [1, 2]. Despite their ecological importance, Z. marina meadows decline globally by 2–5% annually, primarily due to human pressures, including climate change [3, 4]. Z. marina is widespread along coastlines throughout the Northern Hemisphere and has been suggested as a model system for aquatic plants [5]. Because the leaves and roots are constantly submerged, microorganisms attached to Z. marina interact strongly with the surrounding water and sediment bacterial communities [2, 6]. For example, the release of organic molecules (e.g., amino acids and sugars) and gases (CO2, N2) by aquatic macrophytes is associated with the presence of specific sulfate-reducing and nitrogen-fixing bacterial communities [7,8,9,10]. However, interactions with the surrounding environment extend beyond directly attached bacteria and may also be recognized within the Z. marina meadow, e.g., by the bacterial community in the sediment surrounding the roots [11,12,13] as well as the edges of Z. marina meadow patches [5].

In contrast to the effects on sediment, the influence of Z. marina on bacteria in the water columns of its meadows is still debated and seems to depend on environmental conditions [14,15,16,17,18]. For instance, based on a reduction in colony-forming bacterial units in the intertidal regions of tropical islands, Lamb et al. [19] suggested that Z. marina acts as a natural water filtration system for allochthonous pathogenic bacteria. Similarly, other studies have found lower levels of potentially pathogenic bacteria in Z. marina meadows than in areas without Z. marina [20,21,22,23]. Tasdemir et al. [24] and Millan et al. [25] reported the presence of potentially antibiotic-producing bacteria in Z. marina leaves, suggesting an influence on the surrounding seawater. However, yet other studies did not provide conclusive evidence of pathogen reduction by Z. marina in coastal waters [14, 15, 26]. Hence, the interaction between Z. marina and the surrounding bacterial community is not fully understood [14,15,16, 18].

In this study, we characterized the water and sediment microbiota inside, on the edge of, and outside of Z. marina meadows at different salinity levels at Baltic Sea coastal sites. Z. marina meadows are present within a salinity range of 5–35. Within this salinity range, especially at the salty divide at salinities 6–9, the sediment and water microbiota change from marine-related to freshwater-related organisms [27,28,29,30,31]. The salinity-driven shift in the biota at the salty divide [29, 32, 33] is a fundamental concept in ecology that has effects on species richness and community composition in different biological groups [34]. Few studies have investigated the effects of the host and its surrounding microbiota among the salty divide [35, 36]. Previous studies investigating the impact of salinity on bacterial communities associated with macrophytes suggested that host and habitat were the most important factors in structuring their composition [37, 38], with the strongest changes occurring in the salty divide [39, 40]. We hypothesized that (1) the direct effect of Z. marina on the sediment and water microbiota differs depending on the salinity and that (2) the effect of Z. marina on the meadow water and sediment microbiota decreases with increasing distance from the Z. marina meadow. To address these hypotheses, we investigated Z. marina meadows in the Baltic Sea, where a stable salinity gradient with a characteristic sediment and water bacterial community exists [29, 30] and where Z. marina is highly abundant along the salinity gradient [41].

Materials and methods

Samples were collected as described by Riedinger et al. [26] in shallow coastal areas of the Baltic Sea (Fig. 1A-C). Each sampling station consisted of three subsamples taken within the Z. marina meadow (“inside”), approximately 15 m from the meadow (“edge”) and at least 100 m from the Z. marina meadow (“outside”) (Fig. 1D). Samples for DNA extraction from water and sediment were collected from each substation. For water, 100 mL syringes were used to sample ca. 5 cm from the Z. marina and 20 cm above the sediment surface. Sediment samples were collected using sterile plastic tubes by scraping the upper 1 cm at each station.

Fig. 1
figure 1

Sampling location (A) alpha-mesohaline stations colored by salinity; (B) overview of the Baltic Sea sites where the subgraphs are marked (C) beta-mesohaline stations colored by salinity. (D) Sampling scheme for each station. Samples were collected inside the Zostera marina field, on the edge (15 m), and outside (100 m) of the Z. marina field

The Z. marina densities were counted in triplicate 20 × 20 cm squares. Leaf length was measured for 30 plants in each meadow. All samples were immediately transferred to a 4 °C cooler and stored (maximally 8 h) until processing. Salinity, temperature, and depth were measured using a CTD48M (Sea & Sun Technology). The dissolved oxygen and pH were measured using a HQ40D portables 2-channel multimeter. Phosphate (PO43−), nitrate (NO3), nitrite (NO2), ammonium (NH4+), and silicate (SiO2) concentrations were measured using a Seal Analytical QuAAtro automated continuous flow analyzer (SEAL Analytical Ltd. Nordestedt. Germany). Chlorophyll-a (Chl-a) was measured fluorometrically using a 10- AU-005-CE fluorometer (Turner San Jose, USA).

To determine bacterial and eukaryotic community composition, water was filtered through 0.2 μm polyvinylidene fluoride membrane filters (Merck, Darmstadt, Germany) and shock frozen in liquid nitrogen. DNA was extracted after bead beating using a DNeasy PowerSoil Pro Kit Pro (Qiagen. Hilden. Germany) following the manufacturer’s protocol after grinding the filters. For the sediment samples, a subsample (300 mg) of the sediment was transferred into bead-beating tubes. The samples were placed on ice, sonicated twice for 7 min, and then subjected to bead-beating for 30 s at 4 m/s on a vortex adaptor. Subsequently, the manufacturer’s instructions were followed, and the DNA yield was quantified using the PicoGreen assay (Thermo Fisher, Waltham. USA).

The 16S rRNA genes were amplified using V3–V4 primers [29], and the 18S rRNA genes were amplified by targeting the eukaryotic V4 region [42]. The samples were sequenced on a MiSeq (Illumina Inc. San Diego. CA. US) for the 16S and 18S rRNA genes by SciLifeLab/NGI (Solna. Sweden). Amplicon sequencing data were processed using the workflow described by Riedinger et al. [26] (https://github.com/biodiversitydata-se/amplicon-multi-cutadapt). ASVs with similar taxonomic assignment were aggregated.

For the bacterial community, sequences assigned to Archaea, chloroplasts, and mitochondria were removed before analysis because the primers used covered only minor parts of these groups. After the removal of these groups, a total of 11,599,548 reads were available for the bacterial community (Appendix Fig. 1A). For the eukaryotic data, reads assigned to unclassified Opisthokonta, Fungi, Embryophyceae (including Zostera spp.), and Metazoa were excluded, resulting in a total of 20,702,520 reads (Appendix Fig. 1B).

Statistical analysis

The amplicon sequence variant (ASV) richness (SOBS) and Shannon index were estimated using Explicet [43], which performs rarefaction-based analysis via bootstrapping. For all stations, bootstrap resampling was conducted at the size of the smallest library at the rarefaction point (Bacteria: 10,257 excluding two samples; Eukarya 25,417 excluding two samples) to compare ASVs between libraries with equal sampling efforts. The Kruskal–Wallis test and a post hoc Tukey`s pairwise test was used to determine significant differences between the numbers of ASVs in the samples. A sequential Bonferroni correction was used for multiple comparisons. For the sample groups, p > 0.01 indicated that the variances were not significantly different from each other.

Due to the large spread between individual sampling sequencing depth (Bacteria: 8,720 − 138,054 reads, Eukarya: 9,699 − 213,379 reads), quality-trimmed sequencing reads were transformed using a centered log-ratio (CLR) transformation. This compositional data analysis approach [44] separates data variance according to differences rather than abundances. PCA on the pairwise Aitchison distances (Euclidean distance on CLR-transformed data [45]) was used as an exploratory analysis of the microbiota composition [46]. Then, the data were split into sub-datasets, and each micro-environment was tested separately to assess the impact of environmental factors. The effects of environmental variables on the community composition patterns of water and sediment were analyzed using PERMANOVA and two-way PERMANOVA based on Aitchison distances. The number of permutations was set to 9999, and p < 0.01 was considered significant.

For analyzing differences between microbiota at different water depths, the averages of Z. marina abundance and length at the different sampling stations were calculated, with values above the average categorized as high and below the average as low (Appendix Table 1).

Table 1 Physicochemical parameters measured at the stations (see Fig. 1A-C). The substations “inside, edge, and outside” (see Fig. 1D) were averaged because the differences were minimal. Zostera marina shoot length was determined from inside samples. (Chl a = Chlorophyll a, DOC = dissolved organic carbon, PON = particular organic matter, DN = dissolved nitrogen, NA = not available)

Results

Characterization of sampling sites

The samples were categorized following a modified Venice system [47], which takes the salty divide into account (Table 1): alpha-mesohaline conditions (salinity 9.3–15.5; stations BV5, BV6, BV7, BV8, BV10, BV22, BV24) and beta-mesohaline conditions (salinity 6.1–7.7; stations BV9, BV11, BV12, BV14, BV15, BV16, BV18, BV21). Both alpha- and beta-mesohaline sites were characterized by an average depth of ~ 2.6 m, low waves (~ 0.2 m), and wind speed (2.1 Bft). The average temperature and salinity were higher at alpha-mesohaline conditions compared to beta-mesohaline conditions. Chl-a, POC, PON, and SiO2 were, on average, more than twice as high under beta-mesohaline conditions than under alpha-mesohaline conditions. In particular, the beta-mesohaline station BV11 (Greifswalder Bodden) where higher compared to the average chemical measurements for Chl-a (19.7 mg m− 3), PON (34.9 µmol L− 1), POC (265.3 µg L− 1), DOC (614.8 µmol L− 1) and SiO2 (43.3 µmol L− 1) indicating the presence of a phytoplankton bloom.

Microbial patterns in sediment and water within and outside the Zostera marina meadows

PCA of the bacterial (Fig. 2A) and eukaryotic (Fig. 2B) community composition showed a clear separation between the water and sediment microbiota along the first component (two-way PERMANOVA: Bacteria F = 131, p < 0.01; Eukarya F = 87, p < 0.01), independent of all other factors measured. Along the second component, a separation between alpha- and beta-mesohaline samples was apparent (two-way PERMANOVA: Bacteria F = 19, p < 0.01; Eukarya F = 16, p < 0.01). The interaction between the sediment/water and salinity category was also significant (two-way PERMANOVA: Bacteria F = 15, p < 0.01; Eukarya F = 11, p < 0.01). Because the sediment and water microbiota were clearly separated, we analyzed them separately. Envfit analysis of the sediment and water microbiota showed that the community composition was strongest associated (highest r²) with salinity (Table 2).

Fig. 2
figure 2

Principal component analysis of the (A) Bacterial community and (B) Eukaryotic community based on Aitchison distances in sediment (brown) and water (blue). Dots indicate samples taken under alpha-mesohaline conditions (salinity 9.3–16), and circles samples taken under beta-mesohaline conditions (salinity 6-7.8)

Table 2 Results of the envfit analysis. PC1 and PC2 show the relationships of the environmental variables with the first and second PCA axes (based Aitchison distance). R² is the variation explained by the multiple regression model. Pr(> r) indicates the significance of the multiple regression. Significance codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05. (A) bacterial community water samples and (B) bacterial community sediment samples (C) eukaryotic community water samples and (D) eukaryotic community sediment samples

To examine the potential influence of the presence of Z. marina in water and sediment microbiota, samples were divided according to proximity to Z. marina meadows (inside, edge, and outside) and analyzed separately for alpha- and beta-mesohaline conditions and for water and sediment samples. PC analysis and PERMANOVA revealed no significant effect of proximity to Z. marina on the water bacterial and eukaryotic community composition under neither alpha (F = 0.9, p = 0.54) nor beta (F = 1.1, p = 0.23)-mesohaline conditions (Appendix Fig. 2C, D, G, H). The interaction effect between salinity and proximity to Z. marina was also nonsignificant (two-way PERMANOVA: F = 0.59, p = 0.72). Similarly, no significant differences in the water column alpha-diversity (neither SOBS nor Shannon index) in relation to Z. marina proximity were detected under either alpha-mesohaline or beta-mesohaline conditions (Kruskal-Wallis test: p > 0.01; Fig. 3, Appendix 3).

Fig. 3
figure 3

Alpha diversity (Sobs) of the (A) bacterial community and (B) eukaryotic community

In contrast, both the bacterial and eukaryotic community compositions of the sediment samples were significantly different inside the Z. marina meadow compared to on the edge or outside the meadows for the alpha-mesohaline samples (F = 2.7, p < 0.01; Fig. 4, Appendix 2 B, F). However, in the lower saline beta-mesohaline samples, the effect of Z. marina proximity was not observed on the microbiota (F = 1.0 p = 0.3; Appendix 2 A, E). The interaction effect between salinity and proximity to Z. marina was insignificant (two-way PERMANOVA: F = 0.84, p = 0.08). Similar to the effects on the community composition, the presence of Z. marina affected the alpha-diversity (SOBS and Shannon index) of the bacterial sediment community, only at alpha-mesohaline conditions, where a significantly higher SOBS was observed inside the Z. marina meadow than at the edge and outside (Kruskal-Wallis test: p < 0.01; Fig. 3). Interestingly, no significant differences were found in the SOBS and Shannon index between the inside-edge-outside samples of the eukaryotes for sediment or for the water samples (Kruskal-Wallis test: p > 0.01, Fig. 3, Appendix 3).

Fig. 4
figure 4

Effect of Zostera marina at alpha-mesohaline conditions on the sedimental (A) bacterial community and (B) sediment eukaryotic community

Z. marina leaf length differed significantly between alpha and beta-mesohaline conditions (p < 0.01, Appendix 4). Additionally, there was a significant difference in the bacterial and eukaryotic communities among the different Z. marina length categories (Table 3). In addition, Z. marina abundance had a significant influence (p < 0.01) on the microbiota in contrast to water depth.

Table 3 PERMANOVA of the bacterial (A) and eukaryotic (B) community compositions at water depth and of Z. marina abundance and leaf length categories. (A)

Differences in taxonomic composition

The most abundant (> 1%) phyla/classes of the bacterial communities in our study included Bacteroidetes (12–31%), Alphaproteobacteria (7–33%), and Gammaproteobacteria (10–22%; Fig. 5). The sediment samples contained higher relative abundances of Gammaproteobacteria (sediment 18% vs. 13% water), Planctomycetota (sediment 8% vs. 2% water), and Desulfobacterota (sediment 8% vs. 1% water). Alphaprotoeobacteria (sediment 7% vs. 27% water), Bacteroideta (sediment 17% vs. 24% water), Myxococcota (sediment 5% vs. > 1% water), Chloroflexi (sediment 4% vs. > 1% water), and Cyanobacteria (sediment 9% vs. 16% water) differed. Gammaproteobacteria (beta-mesohaline 15% vs. 21% alpha-mesohaline), Actinobacteria (beta-mesohaline 19% vs. 12% alpha-mesohaline) and Bacteroideta (beta-mesohaline 14% vs. 20% alpha-mesohaline) differed between the alpha and beta-mesohaline conditions. Within the Z. marina meadow under alpha-mesohaline conditions, a high abundance of Campylobacterota (4% inside vs. > 1% outside) was observed. The alpha-mesohaline and beta-mesohaline water bacterial community compositions also differed. Alphaproteobacteria (beta-mesohaline 21% vs. 32% alpha-mesohaline) and Bacteroidota (beta-mesohaline 18% vs. 30% alpha-mesohaline) were more abundant under alpha-mesohaline conditions, whereas Actinobacteria (beta-mesohaline 21% vs. 10% alpha-mesohaline) and Cyanobacteria (beta-mesohaline 22% vs. 9% alpha-mesohaline) were more abundant under beta-mesohaline conditions.

Fig. 5
figure 5

Bacterial phyla/classes (A) and eukaryotic phyla/classes (B) with abundances > 1% in the sediment and water at alpha-mesohaline and beta-mesohaline conditions

On a finer taxonomic level, Cyanobium was most abundant ASV in the water column (max = 21%) especially under beta-mesohaline conditions, but it was also abundant in all other habitats (Fig. 6A). The abundant bacterial groups typically found in the water column were SAR11 clades, NS3a, NS5, SAR86, and Planktomarina under alpha-mesohaline conditions, whereas hgcI (beta-mesohaline 9.3% vs. 2.1% alpha-mesohaline) and CL500-29 (beta-mesohaline 4.8% vs. 1.1% alpha-mesohaline) were more abundant under beta-mesohaline conditions. The bacterial communities in the sediments were dominated by Cyanobium, Illumatobacter, unclassified Saprospiraceae, and the desulfosarcinal Sva0081 group. The bacterial genus Campylobacterota with higher abundances inside the Z. marina meadow at alpha-mesohaline conditions was mostly assigned to Sulfurovorum. Despite having a low overall abundance (0.8%), the genus has an abundance of 4% within the meadow. Representatives of the genus Sulfurovorum are able to oxidize sulfur.

Fig. 6
figure 6

Heatmap of the relative abundances of (A) bacterial and (B) eukaryotic amplicon sequence variants (ASVs) > 1% in sediment and water at alpha-mesohaline and beta-mesohaline conditions. (unclas. = unclassified)

Similar to the bacterial community composition, the community of microbial eukaryotes differed among the sediment and water samples (Fig. 6B). In water, Picochlorum and Monodus were abundant genera at beta-mesohaline conditions. Rhizosolenia delicatula was abundant under alpha-mesohaline conditions in water and sediment (beta-mesohaline 0% vs. 4.7% and 13.2% alpha-mesohaline). Choricystis sp. were more abundant inside (inside 13.3% vs. 3.6% outside) the Z. marina meadow and Gymnodiniaceae, rather than outside the meadow (inside 0.3% vs. 12.0% outside).

Discussion

Z. marina meadows stabilize sediments, provide food and shelter, and recycle nutrients, providing critical aquatic ecosystem services. Z. marina has been suggested to impact its environment beyond the direct Z. marina leaves (phyllosphere) [14,15,16,17, 19] and roots (rhizosphere) [11,12,13]. Our results showed a limited effect of Z. marina on the surrounding water and sediment microbiota since only at alpha-mesohaline conditions in sediments a significant effect of Z. marina on the microbiota was observed. Therefore, the first hypothesis was partially supported. Previous studies also observed an influence of Z. marina under saltwater conditions [12,13,14, 16, 17, 48,49,50].

The missing effect of Z. marina meadows at lower salinity could be related to the symbiont (microbiota) or host (Z. marina). Among the mechanisms by which Z. marina influences the surrounding microbiome is the release of bioavailable carbon and the resulting accumulation of microbial grazers [51] and bacteria that excrete antimicrobial substances [52]. However, these effects appear to be local and may be detected only when the phyllosphere or the rhizosphere is directly sampled. The release of oxygen from the roots of Z. marina has also been shown to influence the surrounding microbiota [6]. However, since we sampled oxygenated top-layer sediments, the impact of released oxygen was likely minor. Another proposed mechanism is the absorption of pathogens by biofilms [53], which was not observed in in lower salinities [26].

Mechanisms of Z. marina influencing the surrounding sediment and water include a reduction of wave energy and increase sedimentation as a physical effect of Z. marina canopy structure (“canopy effect”) [54, 55]. Under beta-mesohaline conditions, Z. marina has a lower productivity that is connected with a significant reduced leaf length (Table 1, Appendix 4). Shorter leaves also reduce the canopy effect, which reduces the capacity of Z. marina to trap fine sediments, including particle-attached bacterial communities [56] from the water column. At higher salinities, longer leaves are more likely to trap particle-attached bacterial communities, resulting in a significant change in the meadow surface sediment bacterial community [18]. Accordingly, a correlation between the canopy structure of Z. marina and the microbiota was previously reported [5, 18, 55, 57]. However, other local factors, in addition to salinity, including sediment type, history, and wave exposure, can also influence leaf length [58, 59].

In contrast to previous studies, no specific bacterial community was found on the edge of a Z. marina meadow in the water and sediment [5, 18]. Therefore, the second hypothesis was rejected. This could be attributed to the fact that the edges of the Z. marina meadow in our study were predominantly patchy and sampled in a distance of 15 m from the main meadow. Samples closer to the meadow may show a different result. We also found higher sediment bacterial diversity within the Z. marina meadows than outside under alpha-mesohaline conditions. This increased diversity of microbiota in close proximity to plant roots due to exudates was previously described as the “rhizosphere effect” [60]. Other studies [5, 61, 62] in contrast suggested higher diversity in bulk sediment outside the Z. marina meadow because the rhizosphere selects microorganisms from the bare sediment and reduces overall microbial diversity. The discrepancy may be connected to the brackish conditions in our study compared with the saltwater conditions in other studies and the fact that we did not sample the rhizosphere, but rather surface sediment in the Z. marina meadow.

The missing “Z. marina meadow”-effect of the water bacterial and microeukaryotic communities is in line with studies on the San Diego coastal water, in which the location of the sampling site was identified as the major driver of the bacterial community composition rather than proximity to the Z. marina meadow [15]. Other studies have reported a reduction in colony-forming units in Z. marina beds compared with surrounding water [20,21,22,23]. However, these studies were conducted at higher salinity. The reduced canopy effect in our study may also have reduced the effects on the water bacterial community. In addition to the canopy effect, a potential mechanism by which Z. marina directly influences water microbial communities is the release of exudates as dissolved organic matter (DOM) [6, 63]. The high DOM concentration in the Baltic Sea [64], which is caused by terrestrial influences, especially in the littoral zone, may mask the effects of DOM release on aquatic microbial communities. Moreover, DOM can be confounded by salinity; therefore, more studies, including experiments under controlled conditions, are needed to support the observed trends in this study.

Other effects on the microbiota

The bacterial and eukaryotic communities in the sediment and water were significantly different (Fig. 2). This is in accordance with previous studies in which the effect of microenvironment type has been shown to have a stronger impact on determining the bacterial community composition than salinity [14, 17, 65,66,67]. The second most important factor influencing microbiota structure was salinity. Based on the salinity range covered (salinity 6–15), clustering of the samples according to alpha-mesohaline and beta-mesohaline conditions was expected [28,29,30, 34, 68]. The SOBS, which reflects microbial richness, was highest in the sediment (Fig. 3), which is consistent with previous studies [5, 29, 30]. In contrast to the species minimum at the salty divide (“horohalinicum”) observed for macroorganisms [34] and phytoplankton [68], the bacterial and microeukaryotic community species richness in the water was rather constant at different salinities (Fig. 3). No significant effects on bacterial and microeukaryotic SOBS have been reported in previous studies [29, 31, 69, 70]. In accordance with these studies, a shift was observed between the dominance of clades present in saltwater, including SAR11 clade I, NS3a, NS5, SAR86, and Planktomarina, and the dominance of clades typically found in freshwater, including HGC-I, SAR11-III, and CL500-29, at the salty divide. Most ASVs were assigned only to their taxonomic genus due to the short-read length of the method and the lack of species descriptions in current databases [71]. This resulted in a relatively low diversity of taxonomic genera and many undescribed taxonomic clades (Fig. 6).

For the water microbiota, a rather gradual change in community composition was observed with changing salinity, whereas the distinction between the alpha- and beta-mesohaline microbial communities in sediments was more pronounced (Fig. 2). This could be connected with a higher resistance of sediment bacterial communities to changes in salinity [66] and stronger mixing of water masses than in sediment, but also with differences in the sediment characteristics of the two salinity conditions.

Conclusion

A global effort is underway to restore coastal and estuarine ecosystems, and a major part of this effort is the reestablishment of Z. marina [72]. Our findings provide fundamental information about the impact of Z. marina on sediment and water microbiota under brackish conditions. We found that the influence of Z. marina on the microbiota was lower under brackish conditions, especially below a salinity of 9. The shift in salinity has a significant effect on both the microbiota and host appearance. Because leaf length can affect the impact on the sediment bacterial community in a Z. marina meadow, we suggest that the canopy effect is pivotal to the influence of Z. marina meadows on the surrounding environment. This finding supports the previous hypothesis that the effects of Z. marina are rather passive [18, 55]. The results should be considered when managing brackish Z. marina meadow ecosystem services in brackish environments, such as the Baltic Sea, Black Sea, and larger estuaries.

Data availability

The datasets supporting the conclusions of this article are available in the European Nucleotide Archive under accession number PRJEB68222, which is in compliance with the Minimum Information about any (X) Sequence (MIxS) standard through the brokerage service GFBio93. Environmental data are available from IOWMeta (doi.io-warnemuende.de/10.12754/data-2023-0010).

References

  1. Walters K. The effects of complex trophic interactions on a marine microbenthic community. Ecology. 1993;74(5):1475–89.

    Article  Google Scholar 

  2. Duarte CM, Holmer M, Marbà N. Plant–microbe interactions in seagrass meadows. Interact between macro-and Microorganisms Mar Sediments. 2005;60:31–60.

    Article  Google Scholar 

  3. Dunic JC, Brown CJ, Connolly RM, Turschwell MP, Côté IM. Long-term declines and recovery of meadow area across the world’s seagrass bioregions. Glob Change Biol. 2021;27(17):4096–109.

    Article  CAS  Google Scholar 

  4. Duarte CM. The future of seagrass meadows. Environ Conserv. 2002;29(2):192–206.

    Article  Google Scholar 

  5. Ettinger CL, Voerman SE, Lang JM, Stachowicz JJ, Eisen JA. Microbial communities in sediment from Zostera marina patches, but not the Z. marina leaf or root microbiomes, vary in relation to distance from patch edge. PeerJ. 2017;5:e3246.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Ugarelli K, Chakrabarti S, Laas P, Stingl U. The seagrass holobiont and its microbiome. Microorganisms. 2017;5(4):81.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Donnelly A, Herbert R. Bacterial interactions in the rhizosphere of seagrass communities in shallow coastal lagoons. J Appl Microbiol. 1998;85(S1):S151–60.

    Article  Google Scholar 

  8. Nielsen LB, Finster K, Welsh DT, Donelly A, Herbert RA, De Wit R, Lomstein BA. Sulphate reduction and nitrogen fixation rates associated with roots, rhizomes and sediments from Zostera noltii and Spartina maritima meadows. Environ Microbiol. 2001;3(1):63–71.

    Article  CAS  PubMed  Google Scholar 

  9. Pollard PC, Moriarty D. Organic carbon decomposition, primary and bacterial productivity, and sulphate reduction, in tropical seagrass beds of the Gulf of Carpentaria, Australia. Mar Ecol Progress Ser Oldendorf. 1991;69(1):149–59.

    Article  CAS  Google Scholar 

  10. Liu P, Zou S, Zhang H, Liu Q, Song Z, Huang Y, Hu X. Genome-resolved metagenomics provides insights into the microbialmediated sulfur and nitrogen cycling in temperate seagrass meadows. Front Mar Sci. 2023;10:1245288.

    Article  Google Scholar 

  11. Sun Y, Song Z, Zhang H, Liu P, Hu X. Seagrass vegetation affect the vertical organization of microbial communities in sediment. Mar Environ Res. 2020;162:105174.

    Article  CAS  PubMed  Google Scholar 

  12. Zheng P, Wang C, Zhang X, Gong J. Community structure and abundance of archaea in a Zostera marina meadow: a comparison between seagrass-colonized and bare sediment sites. Archaea. 2019;2019.

  13. Jensen SI, Kühl M, Prieme A. Different bacterial communities associated with the roots and bulk sediment of the seagrass Zostera marina. FEMS Microbiol Ecol. 2007;62(1):108–17.

    Article  CAS  PubMed  Google Scholar 

  14. Crump BC, Koch EW. Attached bacterial populations shared by four species of aquatic angiosperms. Appl Environ Microbiol. 2008;74(19):5948–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Webb SJ, Rabsatt T, Erazo N, Bowman JS. Impacts of Zostera eelgrasses on microbial community structure in San Diego coastal waters. Elementa: Science of the Anthropocene. 2019;7.

  16. Crump BC, Wojahn JM, Tomas F, Mueller RS. Metatranscriptomics and amplicon sequencing reveal mutualisms in seagrass microbiomes. Front Microbiol. 2018;9:388.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Fahimipour AK, Kardish MR, Lang JM, Green JL, Eisen JA, Stachowicz JJ. Global-scale structure of the eelgrass microbiome. Appl Environ Microbiol. 2017;83(12):e03391–16.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Iqbal MM, Nishimura M, Sano M, Yoshizawa S. Particle-attached microbes in Eelgrass Vegetation areas Differ in Community structure depending on the Distance from the Eelgrass Bed. Microbes Environ. 2023;38(3):ME23013.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Lamb JB, Van De Water JA, Bourne DG, Altier C, Hein MY, Fiorenza EA, et al. Seagrass ecosystems reduce exposure to bacterial pathogens of humans, fishes, and invertebrates. Science. 2017;355(6326):731–3.

    Article  CAS  PubMed  Google Scholar 

  20. Palazón A, Aragonés L, López I, López-Úbeda R, Saval Pérez JM. Determination of the most influential factors in the concentration of bacteria in coastal waters. 2018.

  21. Deng Y, Liu S, Feng J, Wu Y, Mao C. What drives putative bacterial pathogens removal within seagrass meadows? Mar Pollut Bull. 2021;166:112229.

    Article  CAS  PubMed  Google Scholar 

  22. Sun H, Wang T, Liu S, Tang X, Sun J, Liu X, et al. Novel insights into the rhizosphere and seawater microbiome of Zostera marina in diverse mariculture zones. Microbiome. 2024;12(1):27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Reusch TB, Schubert PR, Marten S-M, Gill D, Karez R, Busch K, Hentschel U. Lower Vibrio spp. abundances in Zostera marina leaf canopies suggest a novel ecosystem function for temperate seagrass beds. Mar Biol. 2021;168:1–6.

    Article  Google Scholar 

  24. Tasdemir D, Scarpato S, Utermann-Thüsing C, Jensen T, Blümel M, Wenzel-Storjohann A, et al. Epiphytic and endophytic microbiome of the seagrass Zostera marina: do they contribute to pathogen reduction in seawater? Sci Total Environ. 2024;908:168422.

    Article  CAS  PubMed  Google Scholar 

  25. Millan C. The effect of seagrass components on the associated bacterial diversity and antimicrobial activity in Zostera marina. 2022.

  26. Riedinger DJ, Fernández-Juárez V, Delgado LF, Sperlea T, Hassenrück C, Herlemann DP, et al. Control of Vibrio vulnificus proliferation in the Baltic Sea through eutrophication and algal bloom management. Commun Earth Environ. 2024;5(1):246.

    Article  Google Scholar 

  27. Crump BC, Armbrust EV, Baross JA. Phylogenetic analysis of particle-attached and free-living bacterial communities in the Columbia River, its estuary, and the adjacent coastal ocean. Appl Environ Microbiol. 1999;65(7):3192–204.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Lozupone CA, Knight R. Global patterns in bacterial diversity. Proceedings of the National Academy of Sciences. 2007;104(27):11436-40.

  29. Herlemann DP, Labrenz M, Jürgens K, Bertilsson S, Waniek JJ, Andersson AF. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011;5(10):1571–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Klier J, Dellwig O, Leipe T, Jürgens K, Herlemann DP. Benthic bacterial community composition in the oligohaline-marine transition of surface sediments in the Baltic Sea based on rRNA analysis. Front Microbiol. 2018;9:236.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Mazur-Marzec H, Andersson AF, Błaszczyk A, Dąbek P, Górecka E, Grabski M, et al. Biodiversity of microorganisms in the Baltic Sea: the power of novel methods in the identification of marine microbes. FEMS Microbiol Rev. 2024;48(5):fuae024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Paver SF, Muratore D, Newton RJ, Coleman ML. Reevaluating the salty divide: phylogenetic specificity of transitions between marine and freshwater systems. Msystems 2018;3(6). https://doiorg.publicaciones.saludcastillayleon.es/10.1128/msystems.00232-18.

  33. Kinne O. Salinity-invertebrates. Marine ecology. 1971.

  34. Remane A. Die Brackwasserfauna. Verhandlungen Der Deutschen Zoologischen Gesellschaft. 1934;36:34–74.

    Google Scholar 

  35. Kivistik C, Käiro K, Tammert H, Sokolova IM, Kisand V, Herlemann DP. Distinct stages of the intestinal bacterial community of Ampullaceana Balthica after salinization. Front Microbiol. 2022;13:767334.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Kivistik C, Knobloch J, Käiro K, Tammert H, Kisand V, Hildebrandt J-P, Herlemann DP. Impact of salinity on the gastrointestinal bacterial community of Theodoxus fluviatilis. Front Microbiol. 2020;11:683.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Hempel M, Blume M, Blindow I, Gross EM. Epiphytic bacterial community composition on two common submerged macrophytes in brackish water and freshwater. BMC Microbiol. 2008;8:1–10.

    Article  Google Scholar 

  38. Behera P, Mahapatra S, Mohapatra M, Kim JY, Adhya TK, Raina V, et al. Salinity and macrophyte drive the biogeography of the sedimentary bacterial communities in a brackish water tropical coastal lagoon. Sci Total Environ. 2017;595:472–85.

    Article  CAS  PubMed  Google Scholar 

  39. Van der Loos LM, D’hondt S, Engelen AH, Pavia H, Toth GB, Willems A, et al. Salinity and host drive Ulva-associated bacterial communities across the Atlantic–Baltic Sea gradient. Mol Ecol. 2023;32(23):6260–77.

    Article  PubMed  Google Scholar 

  40. Mohapatra M, Manu S, Dash SP, Rastogi G. Seagrasses and local environment control the bacterial community structure and carbon substrate utilization in brackish sediments. J Environ Manage. 2022;314:115013.

    Article  CAS  PubMed  Google Scholar 

  41. Asmala E, Gustafsson C, Krause-Jensen D, Norkko A, Reader H, Staehr PA, Carstensen J. Role of eelgrass in the coastal filter of contrasting Baltic Sea environments. Estuaries Coasts. 2019;42:1882–95.

    Article  CAS  Google Scholar 

  42. Hugerth LW, Muller EE, Hu YO, Lebrun LA, Roume H, Lundin D, et al. Systematic design of 18S rRNA gene primers for determining eukaryotic diversity in microbial consortia. PLoS ONE. 2014;9(4):e95567.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Robertson CE, Harris JK, Wagner BD, Granger D, Browne K, Tatem B, et al. Explicet: graphical user interface software for metadata-driven management, analysis and visualization of microbiome data. Bioinformatics. 2013;29(23):3100–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Gloor GB, Macklaim JM, Pawlowsky-Glahn V, Egozcue JJ. Microbiome datasets are compositional: and this is not optional. Front Microbiol. 2017;8:2224.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Aitchison J. Principal component analysis of compositional data. Biometrika. 1983;70(1):57–65.

    Article  Google Scholar 

  46. Ramette A. Multivariate analyses in microbial ecology. FEMS Microbiol Ecol. 2007;62(2):142–60. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/j.1574-6941.2007.00375.x

    Article  CAS  PubMed  Google Scholar 

  47. den Hartog C. Typologie Des Brackwassers. Helgoländer Wissenschaftliche Meeresuntersuchungen. 1964;10(1):377–90.

    Article  Google Scholar 

  48. Cúcio C, Engelen AH, Costa R, Muyzer G. Rhizosphere microbiomes of European seagrasses are selected by the plant, but are not species specific. Front Microbiol. 2016;7:440.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Sun F, Zhang X, Zhang Q, Liu F, Zhang J, Gong J. Seagrass (Zostera marina) colonization promotes the accumulation of diazotrophic bacteria and alters the relative abundances of specific bacterial lineages involved in benthic carbon and sulfur cycling. Appl Environ Microbiol. 2015;81(19):6901–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Markovski M, Najdek-Dragić M, Herndl GJ, Korlević M. Compositional stability of sediment microbial communities during a seagrass meadow decline. 2022.

  51. Vymazal J. Removal of enteric bacteria in constructed treatment wetlands with emergent macrophytes: a review. J Environ Sci Health. 2005;40(6–7):1355–67.

    Article  CAS  Google Scholar 

  52. Sundaravadivel M, Vigneswaran S. Constructed wetlands for wastewater treatment. Crit Rev Environ Sci Technol. 2001;31(4):351–409.

    Article  CAS  Google Scholar 

  53. Stott R, Tanner C. Influence of biofilm on removal of surrogate faecal microbes in a constructed wetland and maturation pond. Water Sci Technol. 2005;51(9):315–22.

    Article  CAS  PubMed  Google Scholar 

  54. Agawin NS, Duarte CM. Evidence of direct particle trapping by a tropical seagrass meadow. Estuaries. 2002;25:1205–9.

    Article  Google Scholar 

  55. Liu S, Wu Y, Luo H, Ren Y, Jiang Z, Zhang X, et al. Seagrass canopy structure mediates putative bacterial pathogen removal potential. Front Mar Sci. 2023;9:1076097.

    Article  Google Scholar 

  56. Rieck A, Herlemann DP, Jürgens K, Grossart H-P. Particle-associated differ from free-living bacteria in surface waters of the Baltic Sea. Front Microbiol. 2015;6:1297.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Iqbal MM, Nishimura M, Haider MN, Yoshizawa S. Microbial communities on eelgrass (Zostera marina) thriving in Tokyo Bay and the possible source of leaf-attached microbes. Front Microbiol. 2023;13:1102013.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Stevenson A, O Corcora TC, Hukriede W, Schubert PR, Reusch TB. Substantial seagrass blue carbon pools in the southwestern Baltic Sea include relics of terrestrial peatlands. Front Mar Sci. 2022;9:Art. Nr. 949101.

  59. Salo T, Pedersen MF, Boström C. Population specific salinity tolerance in eelgrass (Zostera marina). J Exp Mar Biol Ecol. 2014;461:425–9.

    Article  CAS  Google Scholar 

  60. Dotaniya M, Meena V. Rhizosphere effect on nutrient availability in soil and its uptake by plants: a review. Proceedings of the National Academy of Sciences, India Section B: Biological Sciences. 2015;85:1–12.

  61. Ling N, Wang T, Kuzyakov Y. Rhizosphere bacteriome structure and functions. Nat Commun. 2022;13(1):836.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Leach JE, Triplett LR, Argueso CT, Trivedi P. Communication in the phytobiome. Cell. 2017;169(4):587–96.

    Article  CAS  PubMed  Google Scholar 

  63. Chen PW-Y, Annabel CN, Olivia M, Chou W-C, Chen J-J, Mukhanov V et al. Investigation of the Growth and Mortality of Bacteria and Synechococcus spp. in Unvegetated and Seagrass Habitats. 2024.

  64. Seidel M, Manecki M, Herlemann DP, Deutsch B, Schulz-Bull D, Jürgens K, Dittmar T. Composition and transformation of dissolved organic matter in the Baltic Sea. Front Earth Sci. 2017;5:31.

    Article  Google Scholar 

  65. Han Z, Wang C, Lei B, Hui N, Yu Y, Shi Y, Zheng J. A limited overlap of interactions between the bacterial community of water and sediment in wetland ecosystem of the Yellow River floodplain. Front Microbiol. 2023;14:1193940.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Tammert H, Kivistik C, Kisand V, Käiro K, Herlemann DP. Resistance of freshwater sediment bacterial communities to salinity disturbance and implication for industrial salt discharge and climate change based salinization. Front Microbiomes. 2024;2:1232571.

  67. Gordon-Bradley N, Lymperopoulou DS, Williams HN. Differences in bacterial community structure on Hydrilla verticillata and Vallisneria americana in a freshwater spring. Microbes Environ. 2014;29(1):67–73.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Olli K, Ptacnik R, Klais R, Tamminen T. Phytoplankton species richness along coastal and estuarine salinity continua. Am Nat. 2019;194(2):E41–51.

    Article  PubMed  Google Scholar 

  69. Herlemann DP, Lundin D, Andersson AF, Labrenz M, Jürgens K. Phylogenetic signals of salinity and season in bacterial community composition across the salinity gradient of the Baltic Sea. Front Microbiol. 2016;7:1883.

    Article  PubMed  PubMed Central  Google Scholar 

  70. Telesh IV, Schubert H, Skarlato SO. Revisiting Remane’s concept: evidence for high plankton diversity and a protistan species maximum in the horohalinicum of the Baltic Sea. Mar Ecol Prog Ser. 2011;421:1–11.

    Article  Google Scholar 

  71. Glöckner FO, Yilmaz P, Quast C, Gerken J, Beccati A, Ciuprina A, et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J Biotechnol. 2017;261:169–76.

    Article  PubMed  Google Scholar 

  72. Unsworth RK, Jones BL, Coals L, Furness E, Inman I, Rees SC, Evans AJ. Overcoming ecological feedbacks in seagrass restoration. Restor Ecol. 2024;32(4):e14101.

    Article  Google Scholar 

Download references

Acknowledgements

The authors acknowledge support from the National Genomics Infrastructure in Genomics Application Stockholm funded by the Science for Life Laboratory, the Knut and Alice Wallenberg Foundation and the Swedish Research Council, and SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science for assistance with massively parallel sequencing and access to the UPPMAX computational infrastructure. We are grateful to the crew and captain of the R/V Elisabeth Mann Borgese (EMB283), the research divers of the Leibniz Institute for Baltic Sea Research, the University of Rostock, Klaipeda University, and the Estonian Marine Institute, to Jonas Nilsson for assisting with sampling in Sweden, and to Adam Woźniczka and Jarone Pinhassi for providing lab space. In addition, the Leibniz Institute for Baltic Sea Research Warnemünde is conducting research in the new research field “Shallow Water Processes and Transitions to the Baltic Scale.”

Funding

Open Access funding enabled and organized by Projekt DEAL. This work resulted from the BiodivERsA project “Pathogenic Vibrio bacteria in the current and future Baltic Sea waters: mitigating the problem” (BaltVib), funded by the European Union and the Estonian Research Council (grant T210076PKKH / P200028PKKH), the Federal Ministry of Education and Research, Germany (grant 16LC2022A), the Innovation Fund Denmark (grant 0156-00001B), the Research Council of Lithuania (grant S-BIODIVERSA-21-1), the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning FORMAS (grant 2020–02366), the Polish National Science Center (grant 2020/02/Y/NZ8/00009), and the Research Council of Finland (Academy of Finland; grant 344743). In addition, the Leibniz Institute for Baltic Sea Research Warnemünde is conducting research in the new research field “Shallow Water Processes and Transitions to the Baltic Scale.”

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization by DH, AA, CP, LR, MK, SK, MR, MB, ML; Data Curation by all authors; Formal Analysis by DH, LD, DR, VJ, KK, AA, CP, LR, MK, SK, ML; Funding acquisition by DH, AA, CP, LR, MK, TR, MR, ML. Project administration by SK, ML; Writing original draft by DH, DR, ML; Writing, review, and edition by all authors; Investigation by DH, DR, KK, AA, CP, LR, MB, GG, SK, GM, MR, ML; Visualization by DH; Software by AA, LD, VK, DH; Validation by SK.

Corresponding author

Correspondence to Daniel P. R. Herlemann.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have approved the final manuscript.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Herlemann, D.P.R., Delgado, L.F., Riedinger, D.J. et al. Low impact of Zostera marina meadows on sediment and water microbiota under brackish conditions. Environmental Microbiome 20, 2 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-024-00662-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40793-024-00662-6

Keywords