Wenchao Deng1,2
Shunyang Chen1,2,3Shiquan Chen4
Bingpeng Xing1,2
Zhuhua Chan1
Yao Zhang5*Bin Chen1,2,3
Guangcheng Chen1,2,3*- 1Third Institute of Oceanography, Ministry of Natural Resources, Xiamen, China
- 2Observation and Research Station of Coastal Wetland Ecosystem in Beibu Gulf, Ministry of Natural Resources, Beihai, China
- 3Key Laboratory of Marine Ecological Conservation and Restoration, Ministry of Natural Resources, Xiamen, China
- 4Hainan Academy of Ocean and Fisheries Sciences, Haikou, China
- 5State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
Introduction: Seagrass-associated microbial communities play a crucial role in the growth and health of seagrasses. However, like seagrass meadows, seagrass-associated microbial communities are often affected by eutrophication. It remains unclear how eutrophication influences the composition and function of microbial communities associated with different parts of seagrass.
Methods: We employed prokaryotic 16S rRNA gene high-throughput sequencing combining microbial community structure analysis and co-occurrence network analysis to investigate variances in microbial community compositions, potential functions and complexities across sediment, seagrass leaves, and seawater within different eutrophic areas of two adjacent seagrass meadows on Hainan Island, China.
Results: Our results indicated that microbial diversity on seagrass leaves was significantly lower than in sediment but significantly higher than in seawater. Both sediment and phyllosphere microbial diversity showed no significant difference between the highly eutrophic and less eutrophic sites in each lagoon. However, sediment microbial diversity was higher in the more eutrophic lagoon, while phyllosphere microbial diversity was higher in the less eutrophic lagoon. Heavy eutrophication increased the relative abundance of phyllosphere microorganisms potentially involved in anaerobic metabolic processes, while reducing those responsible for beneficial functions like denitrification. The main factor affecting microbial diversity was organic carbon in seawater and sediment, with high organic carbon levels leading to decreased microbial diversity. The co-occurrence network analysis revealed that heavy eutrophication notably reduced the complexity and internal connections of the phyllosphere microbial community in comparison to the sediment and seawater microbial communities. Furthermore, ternary analysis demonstrated that heavy eutrophication diminished the external connections of the phyllosphere microbial community with the sediment and seawater microbial communities.
Conclusion: The pronounced decrease in biodiversity and complexity of the phyllosphere microbial community under eutrophic conditions can lead to greater microbial functional loss, exacerbating seagrass decline. This study emphasizes the significance of phyllosphere microbial communities compared to sediment microbial communities in the conservation and restoration of seagrass meadows under eutrophic conditions.
1 Introduction
Seagrasses are considered one of the most important coastal habitats due to their support for a wide range of keystone marine species from different trophic levels. However, they are facing decline as a result of natural causes and human-induced activities (Bujang et al., 2018; Ramesh et al., 2019). Human activities, such as untreated sewage discharge and fertilizer runoff, which are rich in nutrients, contribute to the loss of seagrass (Bujang et al., 2006; Ugarelli et al., 2017), emphasizing the need to understand mechanisms affecting seagrass health under human influence to prevent further decline. Microbial communities are essential for maintaining a healthy ecosystem by playing important roles in several ecological processes, such as the oxidation of organic carbon, denitrification, nitrification, and sulfur oxidation and reduction (Nealson, 1997; Cardinale et al., 2012). Seagrass-associated microbial communities, found around the rhizosphere and phyllosphere, interact closely with seagrass, contributing to its growth (Ugarelli et al., 2017). For instance, rhizosphere microorganisms can provide inorganic nutrients to seagrass by degrading seagrass-derived organic matters (Evrard et al., 2005) and fixing nitrogen gas, especially in oligotrophic seagrass beds where nitrogen is limited (Patriquin and Knowles, 1972).
The phyllosphere refers to the above-ground surface of a plant leaves, serving as habitat for microorganisms (Vogel et al., 2021). Recent studies have revealed that the phyllosphere of seagrasses hosts rich and diverse microbial communities (Ugarelli et al., 2018; Vogel et al., 2021). These communities are often dominated by various groups capable of diverse metabolic functions, ranging from aerobic respiration to nitrogen fixation and fermentation (Agawin et al., 2016; Crump et al., 2018; Chen et al., 2022), thus supporting or inhibiting the health and growth of their hosts (Ugarelli et al., 2017; Crump et al., 2018). Though both sediment/rhizosphere and phylloshpere microbial communities are essential for the health and growth of seagrass, sediment/rhizosphere microbial communities have received enormous and widespread attention. Martin et al. (2020) have suggested that rhizosphere microorganisms could be treated as indicators of seagrass health. However, the microbial communities on seagrass leaves have been much less studied.
The physicochemical environments within a seagrass meadow ecosystem vary significantly among sediment, seagrass leaves, and seawater. Consequently, environmental changes are likely to have distinct effects on microbial communities in these different habitats. Previous investigations of seagrass-associated microbial communities have solely focused on the communities themselves (Liu et al., 2018; Zhang et al., 2020; Vogel et al., 2021), while minimal information is available about the differences in the impact of environmental changes on microbial communities across different habitats within seagrass ecosystems. Recent studies that simultaneously examined microbial communities associated with different habitats within a seagrass ecosystem found distinct differences among them due to environmental variations (Ugarelli et al., 2018; Rabbani et al., 2021; Banister et al., 2022). For instance, phyllosphere microbial communities were influenced by environmental factors such as temperature, water depth, salinity, and light (Vogel et al., 2020, 2021). In comparison, root and rhizosphere microbial communities in the seagrass meadow near a volcanic CO2 vent were more significantly affected by pH (Banister et al., 2022).
Human activities inducing eutrophication can affect the health of seagrass through various mechanisms such as light reduction, ammonium toxicity, water-column nitrate inhibition, and algal blooms (Burkholder et al., 2007; Dolbeth et al., 2011). Consequently, the health condition of seagrasses affects the community structure and function of microorganisms within seagrass meadows (Liu et al., 2018; Martin et al., 2020; Wang et al., 2020). For instance, the relative abundances of potential pathogen groups in seagrass meadow sediment near nutrient sources were significantly higher than in sediment further from these sources (Liu et al., 2018). However, the effects of eutrophication on the seagrass phyllosphere microbial community are largely unknown, leading to a knowledge gap about the distinct effects of eutrophication on sediment, phylloshpere, and seawater microbial community structures.
The Xincun (XC) and Lian (LA) bays are adjacent lagoons located in the southeastern part of Hainan Island, China (18°23′ to 18°26′, 109°58′ to 110°03′) (Figure 1). Both are natural, nearly closed lagoons with extensive, contiguous seagrass meadows distributed inside, serving as important fishing ports and mariculture areas in Hainan province. However, factors such as enclosed aquaculture, increased fishing activities and land-based wastewater pollution have degraded the ecological environment of both lagoons. Moreover, large amounts of bait are released from fish farms, which release organic matter and consume oxygen in the water, resulting in severe eutrophication of the water bodies. Additionally, weak water exchange between the lagoons exacerbates the eutrophication problem. This eutrophication poses a severe threat to the growth of seagrass populations in the lagoons (Li et al., 2022). To prevent further environmental degradation, fish farming has been banned in the LA lagoon, whereas it remains active in the XC lagoon. Extensive seagrass mortality has been observed along the southeast coast of the XC lagoon, with most of the seagrass leaves fallen off and in a state of decay.
Figure 1. Map of sampling sites in Xincun and Lian lagoons in Hainan Island, South China Sea.
In this study, we investigated the microbial community structure in sediment, seagrass leaves, and seawater habitats in the XC and LA lagoons, with the aim of exploring variations in the impact of eutrophication on phyllosphere microbial diversity, composition, and potential functions compared to sediment and seawater microbial communities within seagrass ecosystems. We also aimed to assess the complexity of the microbial community in each habitat through co-occurrence network analysis, and to examine the interactions of microbial communities from the three habitats under various eutrophication levels. We hypothesized that eutrophication may have greater effects on the diversity and complexity of phyllosphere microbial community than sediment microbial community, as phyllosphere microbial communities are more exposed to eutrophic substances in seawater compared to sediment microbial communities. Those, in turn, weaken the ability of the phyllosphere microbial community to help seagrass cope with environmental changes. This work will helps us understand the microbial mechanisms by which eutrophication affects seagrass health and enhances our recognition of the importance of phyllosphere microbial communities.
2 Materials and methods
2.1 Study area description and sampling
In the XC lagoon, seagrass is mainly distributed in the southern and eastern parts, forming continuous patches, while sporadic distribution is observed in other areas. In the LA lagoon, seagrass is primarily found in the western and southern regions, where it is extensively distributed. In the present study, two sampling sites in the lagoons of XC and LA were chosen, respectively (Figure 1). The sites XC1 and XC2 located in southern and northeastern areas of the XC lagoon, and LA1 and LA2 are in the northwestern area of the LA lagoon. The site2 (XC2 or LA2) was observed more polluted than their counterpart site1 (XC1 or LA1) in each lagoon, respectively. At each site, 100 g of sediment was collected from the surface of seagrass meadow in June 2022. Additionally, leaves of Enhalus acoroides, one of the dominant seagrass species in the two bays (Zhang et al., 2020), were collected near the sediment sampling location. Furthermore, 500 mL of seawater above 0.2 m of the seagrass leaves was collected. The microbes in seawater were filtered on a 0.22 μm-pore size polycarbonate filter (Millipore Corp., Merck KGaA, Germany). All samples were taken in triplicate. The sediment, seagrass leaves, and the polycarbonate filters were immediately flash-frozen in liquid nitrogen upon collection and subsequently stored at −80°C.
2.2 Physicochemical analysis
The seawater temperature, salinity, and pH were determined in situ using a thermometer, a portable refractometer, and a pH meter (WP-81, TPS, Banksia Scientific Co., Brisbane, Australia) respectively. Seawater NH4+, NO3−, and NO2− concentrations were measured according to methods described by Deng et al. (2021) with a four-channel continuous Flow Technicon AA3 Auto-Analyzer (Bran-Lube GmbH, Norderstedt, Germany); the concentration of labile phosphate (LP) was determined using the phosphomolybdenum blue spectrophotometric method with a UV–visible spectrophotometer Evolution™ 350 (Thermo Scientific, Waltham, United States) (He and Honeycutt, 2005); chemical oxygen demand (COD) and dissolved oxygen (DO) were measured using alkaline potassium permanganate method (Rho et al., 2018) and iodometric titration method, respectively, with a titrator; chlorophyll a (Chla) was determined using spectrophotometry method (Aminot and Rey, 2000). For sediment chemical parameters, organic carbon (OC) and total nitrogen (TN) contents were analyzed using a TOC analyser coupled with a nitrogen analyser (Vario TOC Cube, Elementar Analysensysteme, Langenselbold, Germany) after acidification to remove carbonates (Zheng et al., 2023); total phosphorus (TP) content was measured using spectrophotometer method with a UV–visible spectrophotometer Evolution™ 350 (Thermo Scientific, Waltham, United States) (Aspila et al., 1976); sulfide content was determined using iodometric titration method (Pawlak and Pawlak, 1999) with a titrator. All soil data were expressed based on their oven-dry weight (105°C).
2.3 DNA extraction and sequencing
DNA of microbes in sediment was extracted from about 0.5 g (wet weight) of sediments using Fast DNA® SPIN Kit for Soil Kits (MP Biomedicals, California, United States) according to the manufacturer’s specifications. Microbial DNA on polycarbonate filters was extracted using the phenol-chloroform-isoamyl alcohol method (Nercessian et al., 2005). The microbes attached to seagrass leaves were rinsed off with 1% sterile NaCl solution, and then the microbial DNA was extracted using Fast DNA® SPIN Kit for Soil Kits (MP Biomedicals, Irvine, United States). The purity and quality of the DNA was tested using Thermo Scientific NanoDrop 2000c (Thermo Scientific, Waltham, United States). The DNA was stored at −20°C for further analysis. Prokaryotic V3–V4 hypervariable regions in 16S rRNA genes were amplified using universal primers 515F (5′-GTGYCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACNVGGGTWTCTAAT-3′) (Pichler et al., 2018). The amplicon library was sequenced using the Illumina Nova 6000 platform to obtain paired-end sequences (Guangdong Magigene Biotechnology Co., Ltd. Guangzhou, China).
2.4 Sequences processing
The obtained raw paired-end sequences were analyzed using the QIIME version 1.9.1 pipeline (Caporaso et al., 2010). Briefly, the raw sequences were merged into full-length amplicon sequences using Usearch fastq_mergepairs (V10, http://www.drive5.com/usearch/) (Edgar and Flyvbjerg, 2015). The above-processed clean sequences were clustered into operational taxonomic units (OTUs) with a cut-off 0.03 using the UPARSE 11 pipeline (Edgar, 2013), and chimera checking was performed using the UCHIME algorithm (Edgar et al., 2011). The Ribosomal Database Project (RDP) classifier was used to assign OTUs to taxonomic groups based on the SILVA database (Quast et al., 2012) with a minimum confidence of 0.8. To obtain a final OTU matrix that contains the valid read counts of each OTU in each sample, the OTUs that were assigned to chloroplasts or mitochondria, as well as those were not assigned to any kingdom, were discarded.
2.5 Microbial community analysis
The number of prokaryotic sequences varied from 41,399 to 103,425 (average 76,668) per sample. To normalize the data, the sequences were rarefied and subsampled to the minimum number in all samples (Zhang et al., 2016). We were interested in assessing whether significant differences existed in α-diversity of the microbial communities between different habitats (sediment, phyllosphere, and seawater), lagoons (XC Bay vs. LA Bay) and levels of eutrophication (heavy pollution-site2 vs. less pollution-site1). For this purpose, we calculated the ACE, Chao1, Shannon2 (log base 2), and Simpson diversity indices based on the OTUs matrices. Two-way ANOVAs was performed to identify the main and interactive effects of the lagoons and sampling sites (eutrophication) on the α-diversity of the microbial communities, after which post hoc Tukey’s HSD test was applied to test the pairwise group differences. The normality of each index was determined using the Shapiro–Wilk test (Chao, p > 0.05; ACE, p > 0.05; Shannon, p > 0.05; Simpson, p > 0.05) (Royston, 1992).
To visualize the community dissimilarity between sites, principal coordinate analysis (PCoA) were applied based on Bray–Curtis dissimilarities of the OTU matrices using vegan package (Dixon, 2003) in R. A one-way permutational multivariate ANOVA (PERMANOVA) (Anderson, 2005) was applied to identify the differences of microbial community among the three habitats and among the four sampling sites of each habitat based on Bray–Curtis dissimilarities of the OTU matrices. The overall compositional differences between the lagoons and sites of each habitat were tested using two-way PERMANOVA. The permutation test of multivariate homogeneity of group dispersions was applied to evaluate the effect of group dispersion on the PERMANOVA results (Anderson, 2006). The linear discriminant analysis effect size (LEfSe), based on the relative abundance of OTUs at the genus level, was applied to identify the groups that display significant differences in relative abundance between XC and LA lagoons, and between site1 and site2 within each habitat (Segata et al., 2011). The bacterial functions were predicted by FAPROTAX based on the OTU matrices (Louca et al., 2016).
2.6 Correlation analysis
The correlations between microbial community compositions and physicochemical parameters were examined using Mantel test, with dissimilarity matrices of microbial communities based on Bray–Curtis distances between samples. Pearson tests were run to determine the correlations between the physicochemical parameters and α-diversities. The analyses were performed and the results were visualized by using the packages linkET (Huang, 2021) and dplyr (Wickham et al., 2021) in R. The Hmisc package in R was used to examine the Pearson correlations between OTUs and environmental variables (Harrell and Harrell, 2019). Only the top 100 OTUs in total relative abundance across all samples from each habitat are used for analysis. The interactive platform gephi was used to visualize networks (Bastian et al., 2009) with only correlations with absolute value of r > 0.6 and p < 0.05 were retained as valid co-occurrence events.
To access the complexity of the microbial community within the seagrass ecosystems, a co-occurrence network consisting of all samples at each site was constructed. To avoid potential spurious correlations, the top 500 dominant OTUs in each site were selected. The correlation among microbial OTUs was accessed based on Pearson’s correlation, and correlation coefficients >0.8 with a corresponding of p-value <0.01 were retained to generate the network. All networks were constructed using ggClusterNet package in R (Wen et al., 2022), and visualized in Fruchterman–Reingold layout using the interactive Gephi platform.1 To describe the topology and compare the complexity of each network, we calculated a set of metrics, including nodes, edges, average degree, network diameter, average path length, graph density, modularity, and average clustering coefficient, in Gephi. A higher number of nodes and edges, average degree, average clustering coefficient, and graph density, along with lower average path lengths, indicating greater potential complexity within networks (Wang et al., 2023).
3 Results
3.1 Physicochemical parameters
In seawater, the temperature (29.5–33.3°C), salinity (33.3–34.0), and pH (7.85–7.89) showed slight variations among the sites XC1, XC2, LA1, and LA2 (Table 1). The COD and the concentration of Chla and LP were several micrograms per liter at the four sites, with higher levels observed at site2 (XC2 or LA2) than at site1 (XC1 or LA1) in each lagoon. DO was lower at site2 compared to site1 in each lagoon. The TIN concentrations were high (7.87–12.18 μM) across all sites, with the concentrations of NH4+, NO3−, and NO2− were higher at XC1 than at XC2, but lower at LA1 than at LA2. The higher COD and LP content but lower inorganic nitrogen content at XC2 compared to XC1 might be because the higher concentrations of organic carbon and phosphorus nutrients stimulated the assimilation of inorganic nitrogen by heterotrophic microorganisms, thereby decreasing its concentration (Deng et al., 2021). In sediment, the content of all tested nutrients, including OC, sulfide, TP, and TN, especially the sulfide, were higher at site2 than at site1 in each lagoon (Table 1). These results indicate that the site2 (XC2 or LA2) was more eutrophic than site1 (XC1 or LA1) in each lagoon, and the XC lagoon was relative more eutrophic than the LA lagoon.
Table 1. Physicochemical parameters of the seawater and sediment at sampling sites.
3.2 α-diversity of microbial community
The α-diversity indices, including ACE, Chao1, and Shannon, showed a successive and significant decrease in microbial diversity from sediment to seagrass leaves to seawater (Supplementary Figure S1). The Simpson index also showed that microbial diversity in sediment and on seagrass leaves was significantly higher than in seawater (Supplementary Figure S1). For each habitat, the ACE (p < 0.001), Shannon (p < 0.005), and Simpson (p < 0.005) indices indicated that sediment microbial diversity in the XC lagoon was higher than in the LA lagoon. However, all α-diversity indexes (ACE, p < 0.01; Chao1, p < 0.01; Shannon, p < 0.05; Simpson, p < 0.05) suggested that the phyllosphere microbial diversity in the XC lagoon was lower than those in the LA lagoon (Figure 2 and Supplementary Table S1). Moreover, the α-diversity indices (Chao1, p < 0.01; Shannon, p < 0.05; Simpson, p < 0.05) mostly indicated that seawater microbial diversity was higher at the less eutrophic site compared to the highly eutrophic site in both lagoons (Figure 2 and Supplementary Table S1).
Figure 2. The values of α-diversity index ACE, Chao1, Shannon, and Simpson of sediment, phyllosphere, and seawater microbial communities are shown as bar plot. The bar color represents different sites.
3.3 Community composition of prokaryotic 16S rRNA gene among the habitats
PCoA showed a clear distinction of microbial community composition among the sediment, seagrass leaves, and seawater habitats (Figure 3A). The average microbial relative abundance at levels from phylum to family in the three habitats is shown in Figure 4. Alphaproteobacteria, Gammaproteobacteria, and Deltaproteobacteria were the three most dominant class in Proteobacteria, with Deltaproteobacteria being abundant only in sediment (9.3%), and Alphaproteobacteria having the highest relative abundance on seagrass leaves (57.0%) and in seawater (46.5%). Gammaproteobacteria was most abundant on seagrass leaves (14.7%) compared to sediment (8.8%) and seawater (4.6%). Vibrionales within Gammaproteobacteria was most abundant on seagrass leaves. The Deltaproteobacteria in sediment was mainly composed of Desulfobacterales. The relative abundance of the class Bacteroidia gradually decreased from sediment (15.9%) to seagrass leaves (12.7%) to seawater (5.4%). Within the phylum Actinobacteria, the class Actinobacteria (mainly Micrococcales) was detected only abundant in seawater (14.5%), while Acidimicrobiia was detected with the relative abundance ranging from 1 to 2.9% in all three habitats. Campylobacteria (phylum Epsilonbacteraeota) was the most abundant class in sediment (27.0%) and had a minor proportion on seagrass leaves (1.3%). The phyla Chloroflexi (5.0%) and Acidobacteria (2.3%) were exclusively abundant in sediment. Armatimonadetes was exclusively abundant in seawater (2.3%).
Figure 3. Principal coordinate analysis (PCoA) revealed the differences of the prokaryotic community composition among the habitats of sediment (red), seagrass leaves (green), and seawater (blue) (A), and among the four sampling sites in sediment (B), on seagrass leaves (C), and in seawater (D), respectively.
Figure 4. Taxonomic composition of microbial communities from the (A) sediment (Sed), (B) seagrass leaves (SL), and (C) seawater (SW). The pie area represents the average relative abundance of taxa in each habitat. The unlabeled area in family level are marked as the others with relative abundance <2%.
3.4 Comparison of prokaryotic community structure among sites in each habitat
PCoA analyses indicated that microbial community structure differed significantly among the four sampling sites (Figures 3B–D). Further two-way PERMANOVA analyses (Supplementary Table S2) showed that lagoons (XC vs. LA) played a major role in determining sediment microbial composition (F = 9.2, p = 0.001) compared to sampling sites (site1 vs. site2) (F = 3.8, p = 0.022; Figure 3B). In contrast, both lagoon (F = 13.9, p = 0.001) and sampling sites (F = 11.4, p = 0.001; Figure 3D) played similar roles in determining seawater microbial composition. The combined effect of lagoons and sampling sites was also statistically significant for both sediment (F = 4.3, p = 0.009) and phyllosphere (F = 5.9, p = 0.002) microbial community. For seagrass leaves, microbial composition was significantly affected by lagoons (F = 6.3, p = 0.001) but not by sampling sites (F = 1.3, p = 0.213; Figure 3C).
LEfSe results showed significant differences in the relative abundance of microbial groups (including class and order) between the XC and LA lagoons: 17 groups in sediment, 20 groups in the phyllosphere, and 10 groups in seawater (Figures 5A,C,E). However, there were only 4 groups in sediment, 0 groups in the phyllosphere, and 7 groups in seawater showing significant differences between highly and less eutrophic sites (Figures 5B,D,F). On seagrass leaves, Campylobacterales, Bacteroidales, Oceanospirillales, and Vibrionales were enriched in XC, whereas Microtrichales, Bacteroidia, and Clostridiales were enriched in LA (Figure 5C).
Figure 5. Relative abundance of sediment (A,B), phyllosphere (C,D), and seawater (E,F) prokaryotic taxa significantly differentiated between XC and LA (A,C,E) and between site1 and site2 (B,D,F), identified by linear discriminant analysis effect size (LEfSe). Only lineages with linear discriminate analysis (LDA) values >3.5 are displayed. The diameter of each small solid circle is proportional to abundance of the given taxon. The multiclass analysis strategy was all-against-all (more strict). Large circles indicate phylogenetic levels (from domain to genus) in reverse order. Classes and orders are labeled.
3.5 Correlation between microbial community and physicochemical parameters
Pearson correlation analysis showed that seawater and phyllosphere microbial α-diversity exhibited various correlations with seawater physicochemical parameters (Figures 6A,B). For the seawater microbial community, seawater salinity and DO were positively correlated with diversity, while COD, Chla and LP were negatively correlated with diversity, and nitrogen nutrients showed no obvious correlation with diversity (Figure 6A). For the phyllosphere microbial community, seawater salinity and DO were also strongly correlated with diversity, but the correlations of COD, Chla, and LP with diversity were weakened; however, nitrogen nutrients became have some weak correlations with diversity: NO2− was positively correlated, and NO3− and NH4+ were negatively correlated (Figure 6B). Regarding sediment physicochemical parameters, OC had a negative correlation, while TP and TN showed a positive correlation with sediment microbial diversity. These correlations were reversed for phyllosphere microbial diversity, where OC showed a positive correlation, and TP and TN showed negative correlations (Figures 6C,D). Sulfide showed no obvious correlation with either sediment or phyllosphere microbial diversity (Figures 6C,D).
Figure 6. Correlations of microbial community compositions and α-diversity (based on OTUs) with geochemical parameters of (A) seawater and (C) sediment, and the correlations of phyllosphere microbial compositions and α-diversity with geochemical parameters of (B) seawater and (D) sediment. The color of lines indicates the p-value of the Mantel test; the width of color lines indicates the r-value. Asterisks indicate Pearson’s p < 0.05 and the color bar is based on Pearson’s correlation coefficients. DO, dissolved oxygen; COD, chemical oxygen demand; TIN, total inorganic nitrogen; LP, labile phosphorus; OC, organic carbon; TP, total phosphorus; TN, total nitrogen.
Beside seagrass salinity and DO being significantly linked with phyllosphere microbial composition, the physicochemical parameters generally showed no significant correlation with microbial compositions (Figure 6). The co-occurrence network showed that microbial taxa (OTUs) were generally connected with only one or two physicochemical parameters, except that many OTUs of phyllosphere Rhodobacterales connected with two or more seawater physicochemical parameters (Supplementary Figure S2). Notably, among the seawater physicochemical parameters, NO2− had the highest number of links with seawater and phyllosphere microbial taxa. For sediment physicochemical parameters, OC had the highest number of links with phyllosphere microbial taxa.
3.6 Microbial community analysis by using co-occurrence networks and ternary plot
The topological features of the network varied considerably between the two lagoons and between the highly and less eutrophic sites (Figure 7A and Supplementary Table S3). Compared to the highly eutrophic site in each lagoon, microbial community at the less eutrophic site had a more complex network with higher edges, a higher average degree, density and average clustering coefficient, and a lower average path length and diameter. Five dominant ecological clusters were identified across the four sites (Figure 7A). Members in Cluster 1 and Cluster 2 accounted for the highest relative abundance in sediment (> 67%), Cluster 3 and Cluster 4 were dominated by taxa from seagrass leaves (>72%), and Cluster 5 was dominated by taxa from seawater (>73%) (Figure 7B). Among these clusters, Cluster 1 and Cluster 2 had the highest average degree, and Cluster 5 had the lowest average degree (Supplementary Table S4).
Figure 7. (A) Co-occurrence network of the sediment, phylloshpere, and seawater microbial community at each site. Only Pearson’s correlation coefficient (r > 0.8 or r < 0.8 significant at p < 0.01) is shown. The nodes and edges are colored according to cluster. Node size is proportional to the average degree of each OTU, and edge thickness is proportional to the weight of each correlation. (B) Bar plot shows the composition of the clusters in class level and their total relative abundance in each habitat at each site.
While the average degree of the sediment clusters (Clusters 1 and 2), seawater cluster (Cluster 5), and the phylloshpere Cluster 4 remained relatively stable across the four sites, the average degree of the phylloshpere Cluster 3 decreased greatly from lagoon LA to XC (Figure 7 and Supplementary Table S4). In addition, the percentages of the number of links (edges) between Cluster 1 and Cluster 2 relative to the total edges within these clusters were not higher in the XC lagoon (15.2%–17.2%) compared to the LA lagoon (10.7%–15.2%), whereas those values between Clusters 3 and 4 were much lower (0.3%–1.2%) in XC than in LA (3.2%–3.3%) (Figure 7 and Supplementary Table S5). Additionally, the links between Cluster 3 and Cluster 4 were much weaker at the highly eutrophic site (0.3%) compared to the less eutrophic site (1.2%) in XC, while those values were similar between the highly and less eutrophic sites in the LA (Supplementary Table S5).
Given the large variation in the relative abundance of each OTU among the three habitats and strict criteria (r > 0.8, p < 0.01) used for building co-occurrence networks, no connection was observed between clusters dominated by taxa from different habitats in the co-occurrence networks (Figure 7). Therefore, we also used a ternary plot to access the interconnections of microbial groups among the three habitats at each site (Figure 8). The results indicate that most genera were enriched in one of the three habitats (E_Sed, E_SG, and E_SW), with no or few genera sharing similar relative abundance across all three habitats (C_All) at each site. At site XC1, the majority of genera were located near or on the axis SG-Sed, with the fewest genera near or on the Sed-SW axis (Figure 8A). Conversely, at site XC2, most genera were located near or on the Sed-SW axis (Figure 8B). Similarly, the total relative abundance of the connecters between seagrass leaves and sediment (Sed-SG) and between seagrass leaves and seawater (SW-SG) was higher at XC1 than at XC2 (Supplementary Figure S3). However, those values did not differ obviously in the LA lagoon (Supplementary Figure S3).
Figure 8. Comparison of microbial composition from the habitats of sediment (Sed), seagrass leaves (SL), and seawater (SW) at (A) XC1, (B) XC2, (C) LA1, and (D) LA2 sites. Circle size represents the highest relative abundance of each genus among the three habitats at each site. If the difference in relative abundance of a genus was less than two-fold across the three habitats, it was noted as a core connector (C_All, gray circle). If the relative abundance of a genus in one habitat was at least two-fold higher than its abundance in the remaining two habitats, the genus was noted as an enriched (E.) genus (open circles) in that habitat. The remaining genera with higher relative abundance in two habitats compared to the other habitat were classified as “connector (C)” (closed color circles) between those two habitats. The relative abundance is the average relative abundance of the three parallels in each habitat of each site.
3.7 Functional prediction of the microbial community
Only 2,735 out of 16,212 OTUs (16.9%) could be assigned to at least one function by using FAPROTAX. Chemoheterotrophy/aerobic chemoheterotrophy were the most abundant functions across all three habitats (Figure 9). Sulfate respiration was the second most abundant function in sediment; other functions related to sulfur metabolism were also more abundant in sediment than on seagrass leaves and in seawater. Functions related to nitrogen metabolism (including nitrate respiration, denitrification, nitrogen fixation), as well as methanol oxidation, methylotrophy, chitinolysis, cellulolysis, and ureolysis, were commonly most abundant on LA seagrass leaves. In contrast, nitrate reduction and fermentation were commonly most abundant on XC seagrass leaves. Additionally, functions related to photoautotrophy had a higher relative abundance on LA seagrass leaves compared to XC seagrass leaves.
Figure 9. The functional composition of microbial community was predicted by using FAPROTAX. Heatmap represent the relative abundance of functional groups of prokaryotic communities.
4 Discussion
4.1 Differences of sediment, phyllosphere, and seawater microbial community structure within seagrass ecosystem
Previous investigations of microbial communities in coastal ecosystem have indicated that microbial diversity was much higher in sediment than in seawater, with distinct microbial composition observed between these two environments (Won et al., 2017; Cleary et al., 2018; Ul-Hasan et al., 2019). However, these studies did not specifically focus on the phyllosphere microbial community. Our results regarding sediment and seawater microbial diversity in line with these findings and further identify that phyllosphere microbial diversity is significantly higher than in seawater and significantly lower than in sediment (Supplementary Figure S1). The microbial community composition also varied obviously among the three habitats (Figures 3A, 4). For example, the Desulfobacterales, which comprises anaerobic sulfate-reducing microorganisms (Burow et al., 2014; Cúcio et al., 2016; Dyksma et al., 2018), and Sulfurovum, which primarily consists of sulfur-oxidizing chemoautotrophs (Inagaki et al., 2004; Giovannelli et al., 2016; Xie et al., 2021), were most abundant in sediment (Figure 9).
High abundance of both sulfur reducers and oxidizers have been reported in the sediment of various seagrass ecosystem (Zhang et al., 2020; Rabbani et al., 2021; Mohapatra et al., 2022). However, in non-seagrass coastal sediment, only sulfur reducers are typically found in high abundance (Cleary et al., 2018; Ul-Hasan et al., 2019; Chen et al., 2020). Previous studies have shown that high concentration of nutrients and organic matter stimulated the abundance of Campylobacterales in non-seagrass coastal sediment (Aires et al., 2019; Salem et al., 2019). Our co-occurrence network analysis also showed that the most abundant OUT of Campylobacterales positively connected with organic carbon content (OC) within sediment (Supplementary Figure S2C). These findings likely explain the high abundance of sulfur oxidizers in seagrass sediments: seagrass provides abundant organic matter to the sediments, and sulfur reducers can utilize organic matter as electron donors to produce reduced sulfide (Wu et al., 2021), which in turn stimulates the growth of sulfur oxidizers. Since reduced sulfur compounds, especially hydrogen sulfide, are known phytotoxins in eutrophic sediment (Goodman et al., 1995; Lamers et al., 2012), their oxidation by the abundant Campylobacterales generated a positive feedback loop of the organic matter provided by seagrass. This process helps mitigate the toxic effects of hydrogen sulfide and promotes seagrass adaptation in the two lagoons (Zhang et al., 2020).
Another function that notably differed among the three habitats was organic matter degradation. Taxa such as Clostridiales, Thalassobaculales, and Bacteroidales, which were predominantly abundant in sediment, play important roles in degrading recalcitrant organic matters, such as seagrass cell wall polysaccharides and humic-rich organic matter generated through the biochemical transformations of dead seagrass (Riviere et al., 2009; van der Lelie et al., 2012; Vanwonterghem et al., 2014; Aires et al., 2019; Grevesse et al., 2022). In contrast, the order Rhodobacterales, which was the most abundant group in seawater and on seagrass leaves, is well-known for its rapid response to low-molecular-weight organic matters (Taylor et al., 2014; Varela et al., 2020; Pontiller et al., 2022). Micrococcales, known for its importance in degrading refractory organic matters such as cellulose, lignin, and polyaromatic hydrocarbons (Liu et al., 2020), was the second most abundant taxon in seawater (Figure 4). Thus, Micrococcales in eutrophic seawater of seagrass ecosystem is likely an important degrader of organic matter derived from dead seagrass biomass.
As seagrass is a key species in seagrass ecosystem, it is crucial to highlight the functions that enriched in the phyllosphere microbial community. Several important functions, such as denitrification, nitrogen fixation, urelysis, methylotrophy, chitinolysis, and cellulolysis, were observed to be enriched on seagrass leaves (Figure 9). This is reasonable because cellulose, methanol, and chitin are exudates released by seagrass leaves (Venkatachalam et al., 2015; Sanders-Smith et al., 2020), which can obtain inorganic nitrogen with the assistance of microbial nitrogen fixation (Agawin et al., 2016) and the degradation of dissolved organic nitrogen (Tarquinio et al., 2018). While seagrass sediment is well-known as a hotspot for denitrification and nitrogen loss (Hemminga et al., 1991; Garcias-Bonet et al., 2018), phyllosphere denitrification has received comparatively less attention. Our findings suggest that denitrification on seagrass leaves may also serve as a potential pathway for nitrogen loss, thereby balancing coastal anthropogenic nitrogen inputs. Consequently, our results revealed that sediment, seagrass leaves, and seawater in seagrass ecosystem are occupied by distinct microbial groups dedicated to metabolizing various kinds of organic matters and nutrients. This diversity in microbial functions is essential for maintaining the health of seagrass ecosystems.
4.2 Differential impact of eutrophication on seagrass microbial communities compared to sediment and seawater microbial communities
High microbial diversity can enhance the stability and resilience of ecosystems, because microbial communities with high diversity are more likely to include species that can adapt to new environments when ecosystems are subjected to external disturbances (Girvan et al., 2004). In the present work, sediment microbial diversity was higher in the relative more eutrophic XC lagoon than in the LA lagoon, while phylloshpere microbial diversity was higher in LA than in XC (Figure 2). This result indicates that the diversity of phyllosphere microbial communities is more susceptible to reduction due to eutrophication compared to the diversity of sediment microbial communities, potentially leading to the decay of seagrass leaves under eutrophic conditions. This inference can be supported by the enrichment of representative Oceanospirillales and Vibrionales on seagrass leaves in the relative more eutrophic XC lagoon (Figures 4, 5C). Because, the Oceanospirillales and Vibrionales were also observed a significant increase in abundance on seaweed due to environmental degradation (Aires et al., 2018), and Vibrionales is a representative microorganism that commonly associated with stress and diseased seaweeds, coral, and animals (Meron et al., 2011; Vezzulli et al., 2016; Aires et al., 2018).
We also found that the relative harsher environment in the XC lagoon reduced the relative abundance of phyllosphere functional taxa related to important functions, such as denitrification, while increasing the relative abundance of anaerobic functional taxa involved in processes like fermentation and nitrate reduction (Figure 9). This shift may accelerate seagrass leaf decay and reduce denitrification, exacerbating eutrophication. Additionally, the phyllosphere microbial diversity and compositions between the highly and less eutrophic sites were more similar (Figure 3C) compared to those of seawater microbial compositions (Figures 3D, 5D). This indicates that seagrass leaf plays a pivotal role in shaping the microbial community and maintaining the stability of community structure by releasing exudates such as methanol, cellulose, and phenols (Papazian et al., 2019; Sanders-Smith et al., 2020) under different eutrophic conditions. Another possibility is that the phyllosphere microbial community at the less eutrophic site had been subjected to similar degree of impact by eutrophication as at the highly eutrophic site in each lagoon. The microbial-mediated plant communication theory is suitable to interpret this possibility: the secretion of signaling substances under certain conditions triggers the expression of plant microbiome-shape genes, thus establishing unique plant-associated microbial community; the influenced microbes might, in turn, regulate the plants themselves and other spatially proximate plants via released volatile organic compounds or quorum-sensing signaling molecules (Liu et al., 2023; Lv et al., 2024).
Correlation analysis indicated that COD and OC were the primary factors influencing the seawater and sediment microbial diversity, respectively, and high levels of COD and OC were associated with reduced diversity (Figures 6A,C). This is because a high content of organic carbon can lead to strong selection for specialist hydrocarbon-degrading microorganisms, such as the Rhodobacterales in this study, resulting in the disappearance of certain groups of microorganisms (Nogales et al., 2011). Unlike COD, which only influenced seawater microbial diversity, factors such as the salinity and DO, had similar influence on both seawater and phylloshpere microbial diversity (Figures 6A,B). However, this did not mean that organic carbon was not related to phylloshere microbial diversity. Because, the degradation of high levels of organic matter can reduce DO. The lack of an obvious correlation between seawater COD and phylloshere microbial diversity might be because the phylloshere microbial community can be more directly affected by seagrass-derived organic carbon (Ugarelli et al., 2017).
Microbial composition was not correlated with almost all physicochemical parameters (Figure 6). However, in some environments such as terrestrial grasslands and starved seawater macrocosm, microbial composition was showed strong correlation with environmental factors (Zhang et al., 2021; Cornell et al., 2023). Compared with those relative stable environments, the environmental factors in the XC and LA lagoons can be changed by various pathways, including human activities, seagrass, and tidal action, resulting in more complex and variable interactions between the microbial community and environmental factors. The fact that the majority of the top 100 abundant taxa have no correlation with any environmental factors (Supplementary Figure S2) also suggests that the factors influencing their relative abundance varied across different sites. Meanwhile, we still detected some taxa strongly correlated with environmental factors (Supplementary Figure S2), and some of those correlations are known reasonable. For example, Flavobacteriales, the well know degrader of high-molecular-weight organic matters (Bernardet, 2015), were positively corelated with COD in seawater (Supplementary Figure S2A), and Desulfobacterales, a sulfur reducer, were positively correlated with sulfide in sediment (Supplementary Figure S2C). The most unexpected result is that seawater NO2− had the most number of links with various seawater and phyllosphere taxa (Supplementary Figures S2A,B), despite its concentration being much lower than NO3− and NH4+ (Table 1). This work highlights the importance of paying attention to the influence of seawater NO2− on seawater and phyllosphere microbial composition.
4.3 Significant effects of eutrophication on the complexity of phyllosphere microbial community
The various effects of eutrophication on sediment, phyllosphere, and seawater microbial diversity and composition might influence the microorganisms-mediated flow of energy and elements within the seagrass ecosystem (Konopka, 2009; Ren et al., 2023), thereby impacting the health conditions of the seagrass ecosystem. In ecosystems, exchanges of energy and elements among different species can give rise to complex interactions (Montoya et al., 2006; Yuan et al., 2021), often represented as networks (Pržulj and Malod-Dognin, 2016). It is well-documented that there are links between network complexity and ecosystem functioning and stability (Landi et al., 2018; Yuan et al., 2021; Wang et al., 2023). The co-occurrence network of the integral microbial community from the three habitats at each sampling site revealed network complexity at the highly eutrophic site was lower than at the less eutrophic site in each lagoon and the difference in network properties between the highly and less eutrophic sites was greater in the more eutrophic XC lagoon than in the less eutrophic LA lagoon (Supplementary Table S3). Thus, although the microbial diversity and composition analysis, specifically in sediment and on seagrass leaves, did not find obvious differences between highly eutrophic and less eutrophic sites in each lagoon, the co-occurrence network analysis indicated that heavy eutrophication decreased the complexity of microbial communities and, consequently, the multi-functionality of seagrass leaves (Cheng et al., 2021; Wang et al., 2023).
Moreover, the greater decrease in the connections among phylloshpere microorganisms (Figure 7 and Supplementary Table S4) indicates that phylloshpere microbial complexity may be more susceptible to reduction than sediment and seawater microbial complexity under eutrophic conditions. In addition, there was a greater decrease in internal connections within the phyllosphere microbial community compared to the sediment microbial community under more eutrophic conditions (Figure 7 and Supplementary Table S5). The ternary plot also depicted a similar trend, showing that heavy eutrophication obviously decreased the interaction of phyllosphere microbial communities with sediment and seawater microbial communities, and this trend was more pronounced in the more eutrophic XC lagoon (Figure 8). As microorganisms can communicate with each other by releasing volatile organic compounds or quorum-sensing signaling molecules (Liu et al., 2023), the decrease in internal and external interactions of the phylloshpere microbial community might limit the response of seagrass leaves to environmental changes in sediment, seawater, and on seagrass leaves.
By examining the spatial network dynamics of sediment, phyllospere, and seawater microbial communities within seagrass ecosystem in response to eutrophication, this study provides important insights into the microbial mechanisms by which eutrophication threatens seagrass health. Although environmental changes are known to trigger complex interactive effects on seagrass-associated microbial community structure (Korlević et al., 2021; Vogel et al., 2021; Markovski et al., 2022), it is not clear whether eutrophication has different effects on the microbial community structures in the three habitats. This study documents that eutrophication significantly decreased the complexity and internal interactions of the phyllosphere microbial community more than those of the sediment microbial community, and also reduced the external interaction of the phyllosphere microbial community with sediment and seawater microbial communities. Microbial network complexity is closely associated with the functional structure of microbial communities and ecosystem functional processes (Yuan et al., 2021). Therefore, the phyllosphere microbial community is more susceptible to functional decline due to eutrophication, suggesting that seagrass decline is more likely to originate from the decline in the functions of the phyllosphere microbial communities rather than sediment microbial communities. Furthermore, since networked communities have strong linkages with ecosystem functioning, our work suggests that preserving seagrass phyllosphere microbial network structure should be prioritized over the sediment microbial network structure for future seagrass ecosystem conservation or restoration. In addition, similar to microbial biodiversity, which is dependent on both space and time, network features are also temporally dynamic (Yuan et al., 2021). Thus, future studies on microbial networks within seagrass ecosystem need to integrate the spatial dynamics along with their temporal dynamics for seagrass ecosystem conservation or restoration.
5 Conclusion
In summary, through the simultaneous investigation of sediment, phylloshpere, and seawater microbial communities within seagrass meadow ecosystems that suffered different degrees of eutrophication, we assessed the various effects of eutrophication on microbial communities from the three habitats. The results demonstrated that eutrophication had more pronounced effects on the decrease in biodiversity and complexity of the phyllosphere microbial community compared to those of sediment and seawater microbial communities. This can lead to greater functional loss in the phyllosphere microbial community, reducing their ability to effectively assist seagrass leaves in adapting to environmental changes (Cornell et al., 2023). This work contributes to the understanding of the microbial mechanisms underlying the effects of eutrophication on the health of seagrass ecosystems and highlights the importance of focusing on microbial communities on seagrass leaves during seagrass meadow conservation and restoration. High-resolution spatial and time-series monitoring of phyllosphere microbial communities may become a feasible approach to assist in seagrass conservation and restoration.
Data availability statement
Raw sequencing data of the prokaryotic 16S rRNA gene are available at National Center for Biotechnology Information (NCBI) Sequence Read Archive under BioProject accession number PRJNA1113396 with BioSample accession number SAMN41454142 – SAMN41454177.
Author contributions
WD: Conceptualization, Data curation, Formal analysis, Funding acquisition, Visualization, Writing – original draft. SuC: Investigation, Writing – original draft. SiC: Data curation, Investigation, Writing – original draft. BX: Investigation, Writing – original draft. ZC: Resources, Writing – original draft. YZ: Conceptualization, Writing – review & editing. BC: Writing – review & editing. GC: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by grants from the National Key Research and Development Program of China (2022YFF0802201), the National Natural Science Foundation of China (42188102) and Key Research and Development Program of Hainan Province, China (ZDYF2023SHFZ172). The Scientific Research Foundation of Third Institute of Oceanography, MNR (2023016), and the Xiamen Natural Science Foundation (3502Z202373080) also contributed.
Acknowledgments
The authors are grateful to Dr. Yong Yang, Dr. Danyun Ou, Ms. Limei Xie, Mr. Wenshuo An, and Hongyi Wu for their assistance with the field sampling and laboratory work.
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.
Publisher’s note
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1449545/full#supplementary-material
Footnotes
References
Agawin, N. S., Ferriol, P., Cryer, C., Alcon, E., Busquets, A., Sintes, E., et al. (2016). Significant nitrogen fixation activity associated with the phyllosphere of Mediterranean seagrass Posidonia oceanica: first report. Mar. Ecol. Prog. Ser. 551, 53–62. doi: 10.3354/meps11755
Aires, T., Muyzer, G., Serrão, E. A., and Engelen, A. H. (2019). Seaweed loads cause stronger bacterial community shifts in coastal lagoon sediments than nutrient loads. Front. Microbiol. 9:3283. doi: 10.3389/fmicb.2018.03283
Aires, T., Serebryakova, A., Viard, F., Serrão, E. A., and Engelen, A. H. (2018). Acidification increases abundances of Vibrionales and Planctomycetia associated to a seaweed-grazer system: potential consequences for disease and prey digestion efficiency. PeerJ 6:e4377. doi: 10.7717/peerj.4377
Aminot, A., and Rey, F. (2000). Standard procedure for the determination of chlorophyll a by spectroscopic methods. ICES 112, 1–2.
Anderson, M. J. (2005). Permutational multivariate analysis of variance. Auckland: Department of Statistics, University of Auckland, 32–46.
Anderson, M. J. (2006). Distance-based tests for homogeneity of multivariate dispersions. Biometrics 62, 245–253. doi: 10.1111/j.1541-0420.2005.00440.x
Aspila, K. I., Agemian, H., and Chau, A. (1976). A semi-automated method for the determination of inorganic, organic and total phosphate in sediments. Analyst 101, 187–197. doi: 10.1039/AN9760100187
Banister, R. B., Schwarz, M. T., Fine, M., Ritchie, K. B., and Muller, E. M. (2022). Instability and stasis among the microbiome of seagrass leaves, roots and rhizomes, and nearby sediments within a natural pH gradient. Microb. Ecol. 84, 703–716. doi: 10.1007/s00248-021-01867-9
Bastian, M., Heymann, S., and Jacomy, M. (2009). Gephi: an open source software for exploring and manipulating networks. ICWSM Conf. 3, 361–362. doi: 10.1609/icwsm.v3i1.13937
Bernardet, J.-F. (2015). Flavobacteriales ord. nov. Bergey’s manual of systematics of archaea and bacteria. Trujillo, M. E., Dedysh, S., DeVos, P., Hedlund, B., Kämpfer, P., Rainey, F. A., et al. (New York: Springer). 1–2
Bujang, J. S., Zakaria, M. H., and Arshad, A. (2006). Distribution and significance of seagrass ecosystems in Malaysia. Aquat. Ecosyst. Health Manag. 9, 203–214. doi: 10.1080/14634980600705576
Bujang, J. S., Zakaria, M. H., and Short, F. T. (2018). “Seagrass in Malaysia: Issues and challenges ahead”, in The wetland book II: Distribution, description, and conservation. Eds. C. M. Finlayson, G. R. Milton, R. C. Prentice, and N. C. Davidson. (Netherlands: Springer), 1875–1883.
Burkholder, J. M., Tomasko, D. A., and Touchette, B. W. (2007). Seagrasses and eutrophication. J. Exp. Mar. Biol. Ecol. 350, 46–72. doi: 10.1016/j.jembe.2007.06.024
Burow, L., Woebken, D., Marshall, I., Singer, S., Pett-Ridge, J., Prufert-Bebout, L., et al. (2014). Identification of Desulfobacterales as primary hydrogenotrophs in a complex microbial mat community. Geobiology 12, 221–230. doi: 10.1111/gbi.12080
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cardinale, B. J., Duffy, J. E., Gonzalez, A., Hooper, D. U., Perrings, C., Venail, P., et al. (2012). Biodiversity loss and its impact on humanity. Nature 486, 59–67. doi: 10.1038/nature11148
Chen, Q., Fan, J., Su, J., Ming, H., Sun, Z., Li, M., et al. (2020). Spatial distribution characteristics of bacterial community structure and gene abundance in sediments of the Bohai Sea. Acta Oceanol. Sin. 39, 69–78. doi: 10.1007/s13131-020-1554-8
Chen, J., Zang, Y., Yang, Z., Qu, T., Sun, T., Liang, S., et al. (2022). Composition and functional diversity of epiphytic bacterial and fungal communities on marine macrophytes in an intertidal zone. Front. Microbiol. 13:839465. doi: 10.3389/fmicb.2022.839465
Cheng, X., Yun, Y., Wang, H., Ma, L., Tian, W., Man, B., et al. (2021). Contrasting bacterial communities and their assembly processes in karst soils under different land use. Sci. Total Environ. 751:142263. doi: 10.1016/j.scitotenv.2020.142263
Cleary, D. F. R., Polónia, A. R. M., Becking, L. E., De Voogd, N. J., Purwanto,, Gomes, H., et al. (2018). Compositional analysis of bacterial communities in seawater, sediment, and sponges in the Misool coral reef system, Indonesia. Mar. Biodivers. 48, 1889–1901. doi: 10.1007/s12526-017-0697-0
Cornell, C. R., Zhang, Y., Ning, D., Xiao, N., Wagle, P., Xiao, X., et al. (2023). Land use conversion increases network complexity and stability of soil microbial communities in a temperate grassland. ISME J. 17, 2210–2220. doi: 10.1038/s41396-023-01521-x
Crump, B. C., Wojahn, J. M., Tomas, F., and Mueller, R. S. (2018). Metatranscriptomics and amplicon sequencing reveal mutualisms in seagrass microbiomes. Front. Microbiol. 9:388. doi: 10.3389/fmicb.2018.00388
Cúcio, C., Engelen, A. H., Costa, R., and Muyzer, G. (2016). Rhizosphere microbiomes of European seagrasses are selected by the plant, but are not species specific. Front. Microbiol. 7:440. doi: 10.3389/fmicb.2016.00440
Deng, W., Wang, S., Wan, X., Zheng, Z., Jiao, N., Kao, S. J., et al. (2021). Potential competition between marine heterotrophic prokaryotes and autotrophic picoplankton for nitrogen substrates. Limnol. Oceanogr. 66, 3338–3355. doi: 10.1002/lno.11883
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Dolbeth, M., Cardoso, P., and Pardal, M. Â. (2011). “Impact of eutrophication on the seagrass assemblages of the Mondego estuary (Portugal)” in Eutrophication: causes, consequences and control. eds. A. A. Ansari, S. Singh Gill, G. R. Lanza, and W. Rast (Netherlands: Springer), 225–246.
Dyksma, S., Pjevac, P., Ovanesov, K., and Mussmann, M. (2018). Evidence for H2 consumption by uncultured Desulfobacterales in coastal sediments. Environ. Microbiol. 20, 450–461. doi: 10.1111/1462-2920.13880
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. doi: 10.1038/nmeth.2604
Edgar, R. C., and Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 31, 3476–3482. doi: 10.1093/bioinformatics/btv401
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Evrard, V., Kiswara, W., Bouma, T. J., and Middelburg, J. J. (2005). Nutrient dynamics of seagrass ecosystems: 15N evidence for the importance of particulate organic matter and root systems. Mar. Ecol. Prog. Ser. 295, 49–55. doi: 10.3354/meps295049
Garcias-Bonet, N., Fusi, M., Ali, M., Shaw, D. R., Saikaly, P. E., Daffonchio, D., et al. (2018). High denitrification and anaerobic ammonium oxidation contributes to net nitrogen loss in a seagrass ecosystem in the central Red Sea. Biogeosciences 15, 7333–7346. doi: 10.5194/bg-15-7333-2018
Giovannelli, D., Chung, M., Staley, J., Starovoytov, V., Le Bris, N., and Vetriani, C. (2016). Sulfurovum riftiae sp. nov., a mesophilic, thiosulfate-oxidizing, nitrate-reducing chemolithoautotrophic epsilonproteobacterium isolated from the tube of the deep-sea hydrothermal vent polychaete Riftia pachyptila. Int. J. Syst. Evol. Microbiol. 66, 2697–2701. doi: 10.1099/ijsem.0.001106
Girvan, M. S., Bullimore, J., Ball, A. S., Pretty, J. N., and Osborn, A. M. (2004). Responses of active bacterial and fungal communities in soils under winter wheat to different fertilizer and pesticide regimens. Appl. Environ. Microbiol. 70, 2692–2701. doi: 10.1128/AEM.70.5.2692-2701.2004
Goodman, J. L., Moore, K. A., and Dennison, W. C. (1995). Photosynthetic responses of eelgrass (Zostera marina L.) to light and sediment sulfide in a shallow barrier island lagoon. Aquat. Bot. 50, 37–47. doi: 10.1016/0304-3770(94)00444-Q
Grevesse, T., Guéguen, C., Onana, V. E., and Walsh, D. A. (2022). Degradation pathways for organic matter of terrestrial origin are widespread and expressed in Arctic Ocean microbiomes. Microbiome 10:237. doi: 10.1186/s40168-022-01417-6
Harrell, F. E. Jr., and Harrell, M. F. E. Jr.. (2019). Package ‘hmisc’. CRAN2018. Available at: https://cran.r-project.org/web/packages/Hmisc/Hmisc.pdf. 235–236
He, Z., and Honeycutt, C. W. (2005). A modified molybdenum blue method for orthophosphate determination suitable for investigating enzymatic hydrolysis of organic phosphates. Commun. Soil Sci. Plant Anal. 36, 1373–1383. doi: 10.1081/CSS-200056954
Hemminga, M., Harrison, P., and van Lent, F. (1991). The balance of nutrient losses and gains in seagrass meadows. Mar. Ecol. Prog. Ser. 71, 85–96. doi: 10.3354/meps071085
Huang, H. (2021). LinkET: everything is linkable. R package version 0.0.3. https://github.com/Hy4m/linkET
Inagaki, F., Takai, K., Nealson, K. H., and Horikoshi, K. (2004). Sulfurovum lithotrophicum gen. nov., sp. nov., a novel sulfur-oxidizing chemolithoautotroph within the ε-Proteobacteria isolated from Okinawa Trough hydrothermal sediments. Int. J. Syst. Evol. Microbiol. 54, 1477–1482. doi: 10.1099/ijs.0.03042-0
Konopka, A. (2009). What is microbial community ecology? ISME J. 3, 1223–1230. doi: 10.1038/ismej.2009.88
Korlević, M., Markovski, M., Zhao, Z., Herndl, G. J., and Najdek, M. (2021). Seasonal dynamics of epiphytic microbial communities on marine macrophyte surfaces. Front. Microbiol. 12:671342. doi: 10.3389/fmicb.2021.671342
Lamers, L. P., van Diggelen, J. M., Op den Camp, H. J., Visser, E. J., Lucassen, E. C., Vile, M. A., et al. (2012). Microbial transformations of nitrogen, sulfur, and iron dictate vegetation composition in wetlands: a review. Front. Microbiol. 3:156. doi: 10.3389/fmicb.2012.00156
Landi, P., Minoarivelo, H. O., Brännström, Å., Hui, C., and Dieckmann, U. (2018). Complexity and stability of ecological networks: a review of the theory. Popul. Ecol. 60, 319–345. doi: 10.1007/s10144-018-0628-3
Li, Y., Bai, J., Zhang, L., and Yang, Z. (2022). Mapping and spatial variation of seagrasses in Xincun, Hainan Province, China, based on satellite images. Remote Sens. 14:2373. doi: 10.3390/rs14102373
Liu, Q., Cheng, L., Nian, H., Jin, J., and Lian, T. (2023). Linking plant functional genes to rhizosphere microbes: a review. Plant Biotechnol. J. 21, 902–917. doi: 10.1111/pbi.13950
Liu, S., Deng, Y., Jiang, Z., Wu, Y., Huang, X., and Macreadie, P. I. (2020). Nutrient loading diminishes the dissolved organic carbon drawdown capacity of seagrass ecosystems. Sci. Total Environ. 740:140185. doi: 10.1016/j.scitotenv.2020.140185
Liu, S., Jiang, Z., Deng, Y., Wu, Y., Zhang, J., Zhao, C., et al. (2018). Effects of nutrient loading on sediment bacterial and pathogen communities within seagrass meadows. Microbiologyopen 7:e00600. doi: 10.1002/mbo3.600
Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507
Lv, M. H., Shi, W. C., Li, M. C., Zhou, B., Liu, Y. X., and Gao, Z. (2024). Ms gene and Mr gene: microbial-mediated spatiotemporal communication between plants. iMeta 3:e210. doi: 10.1002/imt2.210
Markovski, M., Najdek, M., Herndl, G. J., and Korlević, M. (2022). Compositional stability of sediment microbial communities during a seagrass meadow decline. Front. Mar. Sci. 9:966070. doi: 10.3389/fmars.2022.966070
Martin, B. C., Alarcon, M. S., Gleeson, D., Middleton, J. A., Fraser, M. W., Ryan, M. H., et al. (2020). Root microbiomes as indicators of seagrass health. FEMS Microbiol. Ecol. 96:fiz201. doi: 10.1093/femsec/fiz201
Meron, D., Atias, E., Iasur, K. L., Elifantz, H., Minz, D., Fine, M., et al. (2011). The impact of reduced pH on the microbial community of the coral Acropora eurystoma. ISME J. 5, 51–60. doi: 10.1038/ismej.2010.102
Mohapatra, M., Manu, S., Dash, S. P., and Rastogi, G. (2022). Seagrasses and local environment control the bacterial community structure and carbon substrate utilization in brackish sediments. J. Environ. Manag. 314:115013. doi: 10.1016/j.jenvman.2022.115013
Montoya, J. M., Pimm, S. L., and Solé, R. V. (2006). Ecological networks and their fragility. Nature 442, 259–264. doi: 10.1038/nature04927
Nealson, K. H. (1997). Sediment bacteria: who’s there, what are they doing, and what’s new? Annu. Rev. Earth Planet. Sci. 25, 403–434. doi: 10.1146/annurev.earth.25.1.403
Nercessian, O., Noyes, E., Kalyuzhnaya, M. G., Lidstrom, M. E., and Chistoserdova, L. (2005). Bacterial populations active in metabolism of C1 compounds in the sediment of Lake Washington, a freshwater lake. Appl. Environ. Microbiol. 71, 6885–6899. doi: 10.1128/AEM.71.11.6885-6899.2005
Nogales, B., Lanfranconi, M. P., Piña-Villalonga, J. M., and Bosch, R. (2011). Anthropogenic perturbations in marine microbial communities. FEMS Microbiol. Rev. 35, 275–298. doi: 10.1111/j.1574-6976.2010.00248.x
Papazian, S., Parrot, D., Burýšková, B., Weinberger, F., and Tasdemir, D. (2019). Surface chemical defence of the eelgrass Zostera marina against microbial foulers. Sci. Rep. 9:3323. doi: 10.1038/s41598-019-39212-3
Patriquin, D., and Knowles, R. (1972). Nitrogen fixation in the rhizosphere of marine angiosperms. Mar. Biol. 16, 49–58. doi: 10.1007/BF00347847
Pawlak, Z., and Pawlak, A. S. (1999). Modification of iodometric determination of total and reactive sulfide in environmental samples. Talanta 48, 347–353. doi: 10.1016/S0039-9140(98)00253-7
Pichler, M., Coskun, Ö. K., Ortega-Arbulú, A. S., Conci, N., Wörheide, G., Vargas, S., et al. (2018). A 16S rRNA gene sequencing and analysis protocol for the Illumina MiniSeq platform. Microbiologyopen 7:e00611. doi: 10.1002/mbo3.611
Pontiller, B., Martínez-García, S., Joglar, V., Amnebrink, D., Pérez-Martínez, C., González, J. M., et al. (2022). Rapid bacterioplankton transcription cascades regulate organic matter utilization during phytoplankton bloom progression in a coastal upwelling system. ISME J. 16, 2360–2372. doi: 10.1038/s41396-022-01273-0
Pržulj, N., and Malod-Dognin, N. (2016). Network analytics in the age of big data. Science 353, 123–124. doi: 10.1126/science.aah3449
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Rabbani, G., Yan, B. C., Lee, N. L. Y., Ooi, J. L. S., Lee, J. N., Huang, D., et al. (2021). Spatial and structural factors shape seagrass-associated bacterial communities in Singapore and peninsular Malaysia. Front. Mar. Sci. 8:659180. doi: 10.3389/fmars.2021.659180
Ramesh, R., Banerjee, K., Paneerselvam, A., Raghuraman, R., Purvaja, R., and Lakshmi, A. (2019). Importance of seagrass management for effective mitigation of climate change. Coast. Manag. 14, 283–299. doi: 10.1016/B978-0-12-810473-6.00015-7
Ren, H., Wang, G., Ding, W., Li, H., Shen, X., Shen, D., et al. (2023). Response of dissolved organic matter (DOM) and microbial community to submerged macrophytes restoration in lakes: a review. Environ. Res. 231:116185. doi: 10.1016/j.envres.2023.116185
Rho, T., Choi, S.-H., Kim, E. S., Kang, N. Y., Cho, S. R., Khang, S. H., et al. (2018). Optimization of chemical oxygen demand determination in seawater samples using the alkaline potassium permanganate method. Ocean Sci. J. 53, 611–619. doi: 10.1007/s12601-018-0040-0
Riviere, D., Desvignes, V., Pelletier, E., Chaussonnerie, S., Guermazi, S., Weissenbach, J., et al. (2009). Towards the definition of a core of microorganisms involved in anaerobic digestion of sludge. ISME J. 3, 700–714. doi: 10.1038/ismej.2009.2
Royston, P. (1992). Approximating the Shapiro–Wilk W-test for non-normality. Stat. Comput. 2, 117–119. doi: 10.1007/BF01891203
Salem, F. B., Saïd, O. B., Cravo-Laureau, C., Mahmoudi, E., Bru, N., Monperrus, M., et al. (2019). Bacterial community assemblages in sediments under high anthropogenic pressure at Ichkeul Lake/Bizerte Lagoon hydrological system, Tunisia. Environ. Pollut. 252, 644–656. doi: 10.1016/j.envpol.2019.05.146
Sanders-Smith, R., Segovia, B. T., Forbes, C., Hessing-Lewis, M., Morien, E., Lemay, M. A., et al. (2020). Host-specificity and core taxa of seagrass leaf microbiome identified across tissue age and geographical regions. Front. Ecol. Evol. 8:605304. doi: 10.3389/fevo.2020.605304
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60–R18. doi: 10.1186/gb-2011-12-6-r60
Tarquinio, F., Bourgoure, J., Koenders, A., Laverock, B., Säwström, C., and Hyndes, G. A. (2018). Microorganisms facilitate uptake of dissolved organic nitrogen by seagrass leaves. ISME J. 12, 2796–2800. doi: 10.1038/s41396-018-0218-6
Taylor, J. D., Cottingham, S. D., Billinge, J., and Cunliffe, M. (2014). Seasonal microbial community dynamics correlate with phytoplankton-derived polysaccharides in surface coastal waters. ISME J. 8, 245–248. doi: 10.1038/ismej.2013.178
Ugarelli, K., Chakrabarti, S., Laas, P., and Stingl, U. (2017). The seagrass holobiont and its microbiome. Microorganisms 5:81. doi: 10.3390/microorganisms5040081
Ugarelli, K., Laas, P., and Stingl, U. (2018). The microbial communities of leaves and roots associated with turtle grass (Thalassia testudinum) and manatee grass (Syringodium filliforme) are distinct from seawater and sediment communities, but are similar between species and sampling sites. Microorganisms 7:4. doi: 10.3390/microorganisms7010004
Ul-Hasan, S., Bowers, R. M., Figueroa-Montiel, A., Licea-Navarro, A. F., Beman, J. M., Woyke, T., et al. (2019). Community ecology across bacteria, archaea and microbial eukaryotes in the sediment and seawater of coastal Puerto Nuevo, Baja California. PloS One 14:e0212355. doi: 10.1371/journal.pone.0212355
van der Lelie, D., Taghavi, S., McCorkle, S. M., Li, L.-L., Malfatti, S. A., Monteleone, D., et al. (2012). The metagenome of an anaerobic microbial community decomposing poplar wood chips. PLoS One 7:e36740. doi: 10.1371/journal.pone.0036740
Vanwonterghem, I., Jensen, P. D., Dennis, P. G., Hugenholtz, P., Rabaey, K., and Tyson, G. W. (2014). Deterministic processes guide long-term synchronised population dynamics in replicate anaerobic digesters. ISME J. 8, 2015–2028. doi: 10.1038/ismej.2014.50
Varela, M. M., Rodríguez-Ramos, T., Guerrero-Feijóo, E., and Nieto-Cid, M. (2020). Changes in activity and community composition shape bacterial responses to size-fractionated marine DOM. Front. Microbiol. 11:586148. doi: 10.3389/fmicb.2020.586148
Venkatachalam, A., Govinda, R., Thirunavukkarasu, N., and Suryanarayanan, T. (2015). Endophytic fungi of marine algae and seagrasses: a novel source of chitin modifying enzymes. Mycosphere 6, 345–355. doi: 10.5943/mycosphere/6/3/10
Vezzulli, L., Grande, C., Reid, P. C., Hélaouët, P., Edwards, M., Höfle, M. G., et al. (2016). Climate influence on Vibrio and associated human diseases during the past half-century in the coastal North Atlantic. Proc. Natl. Acad. Sci U.S.A. 113, E5062–E5071. doi: 10.1073/pnas.1609157113
Vogel, M. A., Mason, O. U., and Miller, T. E. (2020). Host and environmental determinants of microbial community structure in the marine phyllosphere. PLoS One 15:e0235441. doi: 10.1371/journal.pone.0235441
Vogel, M. A., Mason, O. U., and Miller, T. E. (2021). Composition of seagrass phyllosphere microbial communities suggests rapid environmental regulation of community structure. FEMS Microbiol. Ecol. 97:fiab013. doi: 10.1093/femsec/fiab013
Wang, L., Tomas, F., and Mueller, R. S. (2020). Nutrient enrichment increases size of Zostera marina shoots and enriches for sulfur and nitrogen cycling bacteria in root-associated microbiomes. FEMS Microbiol. Ecol. 96:fiaa129. doi: 10.1093/femsec/fiaa129
Wang, X., Zhang, Q., Zhang, Z., Li, W., Liu, W., Xiao, N., et al. (2023). Decreased soil multifunctionality is associated with altered microbial network properties under precipitation reduction in a semiarid grassland. iMeta 2:e106. doi: 10.1002/imt2.106
Wen, T., Xie, P., Yang, S., Niu, G., Liu, X., Ding, Z., et al. (2022). ggClusterNet: an R package for microbiome network analysis and modularity-based multiple network layouts. iMeta 1:e32. doi: 10.1002/imt2.32
Wickham, H., François, R., Henry, L., and Müller, K. (2021). dplyr: a grammar of data manipulation. R package version 1.0.7. Available at: https://CRAN.R-project.org/package=dplyr
Won, N.-I., Kim, K.-H., Kang, J. H., Park, S. R., and Lee, H. J. (2017). Exploring the impacts of anthropogenic disturbance on seawater and sediment microbial communities in Korean coastal waters using metagenomics analysis. Int. J. Environ. Res. Public Health 14:130. doi: 10.3390/ijerph14020130
Wu, B., Liu, F., Fang, W., Yang, T., Chen, G.-H., He, Z., et al. (2021). Microbial sulfur metabolism and environmental implications. Sci. Total Environ. 778:146085. doi: 10.1016/j.scitotenv.2021.146085
Xie, S., Wang, S., Li, D., Shao, Z., Lai, Q., Wang, Y., et al. (2021). Sulfurovum indicum sp. nov., a novel hydrogen-and sulfur-oxidizing chemolithoautotroph isolated from a deep-sea hydrothermal plume in the northwestern Indian Ocean. Int. J. Syst. Evol. Microbiol. 71:4748. doi: 10.1099/ijsem.0.004748
Yuan, M. M., Guo, X., Wu, L., Zhang, Y., Xiao, N., Ning, D., et al. (2021). Climate warming enhances microbial network complexity and stability. Nat. Clim. Change 11, 343–348. doi: 10.1038/s41558-021-00989-9
Zhang, L., Chen, M., Chen, X., Wang, J., Zhang, Y., Xiao, X., et al. (2021). Nitrifiers drive successions of particulate organic matter and microbial community composition in a starved macrocosm. Environ. Int. 157:106776. doi: 10.1016/j.envint.2021.106776
Zhang, Y., Deng, W., Xie, X., and Jiao, N. (2016). Differential incorporation of carbon substrates among microbial populations identified by field-based, DNA stable-isotope probing in South China Sea. PLoS One 11:e0157178. doi: 10.1371/journal.pone.0157178
Zhang, X., Zhao, C., Yu, S., Jiang, Z., Liu, S., Wu, Y., et al. (2020). Rhizosphere microbial community structure is selected by habitat but not plant species in two tropical seagrass beds. Front. Microbiol. 11:161. doi: 10.3389/fmicb.2020.00161
Keywords: seagrass ecosystem, microbial diversity, phyllosphere microbial community, co-occurrence network, eutrophication
Citation: Deng W, Chen S, Chen S, Xing B, Chan Z, Zhang Y, Chen B and Chen G (2024) Impacts of eutrophication on microbial community structure in sediment, seawater, and phyllosphere of seagrass ecosystems. Front. Microbiol. 15:1449545. doi: 10.3389/fmicb.2024.1449545
Edited by:
Yaodong He, Zhejiang Ocean University, ChinaReviewed by:
Mohan Bai, Chinese Academy of Agricultural Sciences, ChinaPeng Zhang, Guangdong Ocean University, China
Copyright © 2024 Deng, Chen, Chen, Xing, Chan, Zhang, Chen and Chen. 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: Guangcheng Chen, Z2MuY2hlbkB0aW8ub3JnLmNu; Yao Zhang, eWFvemhhbmdAeG11LmVkdS5jbg==