Next Article in Journal
Metabarcoding Extends the Distribution of Porphyra corallicola (Bangiales) into the Arctic While Revealing Novel Species and Patterns for Conchocelis Stages in the Canadian Flora
Previous Article in Journal
What Does “ITS” Say about Hybridization in Lineages of Sarsia (Corynidae, Hydrozoa) from the White Sea?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plant Compartments Shape the Assembly and Network of Vallisneria natans-Associated Microorganisms

1
School of Resources and Environmental Engineering, Anhui University, Hefei 230601, China
2
School of Resources and Civil Engineering, Suzhou University, Suzhou 234000, China
*
Author to whom correspondence should be addressed.
Diversity 2023, 15(5), 676; https://doi.org/10.3390/d15050676
Submission received: 29 March 2023 / Revised: 28 April 2023 / Accepted: 11 May 2023 / Published: 18 May 2023

Abstract

:
The submerged plant Vallisneria natans can provide an attachment matrix and habitat for diverse microorganisms and plays an important role in maintaining the structure and function of the shallow lake ecosystem. However, little is known about how V. natans-related microorganism components, especially bacteria, adapt to specific plant compartments. In this study, we investigated the assembly and network of bacterial communities living in different plant compartments (sediment, rhizosphere, rhizoplane, root endosphere, and leaf endosphere) associated with V. natans by 16S rRNA gene sequencing. The results showed that the diversity and network complexity of the bacterial community in the sediment was significantly higher than that in other plant compartments. The bacterial community composition showed that the dominant phyla were Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Desulfobacterota, and Chloroflexi, among which Proteobacteria were extremely abundant in all samples, and there were notable differences in bacterial community composition related to plant compartments. Different networks based on sediment and plant compartments showed distinct co-occurrence patterns and exhibited distinct topological features. Additionally, functional predictions from FAPROTAX indicate that the predominant biogeochemical cycle function of the V. natans-related bacterial community is to participate in the carbon and nitrogen cycle. These results strongly suggested how the microbial community adapted to different plant compartments and provided theoretical and technical data for isolating beneficial bacteria from macrophytes in the future.

1. Introduction

Throughout the plant growth phase, the interrelationships between the plant and the microbial community hold significant importance for plants [1]. These microorganisms, including various mutualistic bacteria, colonize different plant compartments, such as the rhizosphere, rhizoplane, root endosphere, and leaf endosphere [2]. Plant-associated microorganisms play a crucial role in the biogeochemical processes of aquatic ecosystems, providing nutrient sources, enhancing disease resistance, and facilitating pollutant degradation [3,4]. Furthermore, plant-associated microorganisms are enriched by plant secretions and adapt to various ecological niches [5]. Some microorganisms are capable of recognizing signal molecules and settling in specific compartmentalized niches, while other microorganisms are filtered out [6]. This compartment enrichment in particular species is influenced by both biotic and abiotic factors, particularly in plant compartments more closely associated with microbial communities [7]. Previous studies have demonstrated that bacterial communities associated with macrophytes can promote plant growth, nutrient uptake, and tolerance to abiotic and biotic stresses [8]. Additionally, the complex co-associations between bacterial microbiota and macrophytes play a crucial role in improving water quality [9]. For example, Singh et al. used aquatic plants and root bacteria to remediate the chromium-containing wastewater, and after 72 h, each pollution parameter was significantly reduced [10]. In addition, macrophytes can affect sedimental biogeochemical cycling by increasing the activity of microbes such as methanogens [11].
Previous studies have demonstrated that different plant compartments are colonized by distinct microbial communities [12]. Soil habitats serve as an extremely rich microbial reservoir for host–microbial selection [13]. For the root zone, plant-associated microbial communities originate from the soil, where bacterial populations follow a hierarchical selection mechanism, initially colonizing the rhizosphere, adhering to the rhizoplane, and ultimately entering the root endosphere [14]. Root exudates (such as organic acids) can stimulate bacterial growth and extracellular enzyme activity, thus affecting the biogeochemical cycles of carbon, nitrogen, and other elements in the soil [15]. Therefore, it is important to understand the effect of plant roots on the composition of microbial communities in the rhizosphere and rhizoplane. After entering the plant interior, endophytic bacteria may spread throughout the plant and colonize the leaf tissue. However, due to the physiological requirements necessary to occupy these plant ecological niches, only a few bacteria can colonize the leaf tissue of their host plant [16]. The movement of bacteria inside plants is supported by bacterial flagellum and the plant transpiration stream. Migration along that intercellular space requires the secretion of cell-wall-degrading enzymes such as cellulases and pectinases [17]. Therefore, the bacteria that migrate to the leaf tissue are well adapted to this unique endophytic ecological niche [18]. Microbial endophytes in the leaf are also important to the host plant, promoting plant growth and controlling plant diseases [19]. A better understanding of endophytic bacteria may help elucidate their function and potential role in improving plant performance.
Vallisneria natans, a perennial submerged aquatic herbaceous plant, is one of the most important primary producers in shallow lakes and can provide attachment substrates and habitats for a variety of microorganisms [9]. Its biological characteristics help maintain the structure and function of shallow lake ecosystems [20]. The V. natans not only absorbs and accumulates nutrients, organic pollutants, and heavy metal ions in the water through its root system and leaf organs but also secretes oxygen into the sediment to improve the redox environment of the sediment. This species is often selected as an important aquatic plant species when constructing water purification systems and is widely used in the restoration of shallow lake ecosystems and the construction of benign ecosystems [21]. Its sensitivity to high concentrations of environmental pollutants, weak resistance to non-biological stress, and low biological utilization of ions and compounds all restrict its restoration efficiency [22].
In recent years, important progress has been made in the research on the combined remediation technology for submerged plants and microorganisms. However, the application effect of functional bacteria is unstable, and further exploration of functional microbial germplasm resources is needed. In this study, we focused on the following three aspects: (1) revealing the assembly mechanism of microbial communities in different host compartments; (2) mining genetic information of dominant species and key groups in microbial communities; and (3) exploring the potential microbial species in submerged plants that drive the biogeochemical cycle of lakes.

2. Materials and Methods

2.1. Sample Collection and Surface Sterilization

Mature V. natans specimens were collected from the experimental field of Anhui University and randomly divided into three parallel groups (30 plants per group). Non-rhizosphere sediments covered with V. natans vegetation were collected as sediment samples. The collected samples included sediment, rhizosphere soil, rhizoplane soil, roots, and leaves from this sampling site. All soil samples had the salt crusts and litter layer removed. They were placed in aseptic bags and transported back to the lab. All tissue samples were washed with 10 mM PBS buffer and surface sterilized with 70% ethanol for 45 s, then washed 5 times with sterile water. Finally, they were soaked in 0.1% mercuric chloride for 5 min and washed with sterile water 5 times. Sterile filter paper was used to absorb water. Next, 100 μL of sterile water from the final wash was applied to an LB agar plate, which was inverted and incubated at 30 °C for 72 h to assess the efficacy of surface sterilization of plant material. If contamination was observed on the surface of the agar medium, the plant material was discarded.

2.2. Microbial Genomic DNA Extraction, PCR, and Sequencing

The sediment and rhizosphere DNA were extracted from 250 mg soil using the FastDNA® Spin Kit for Soil (MP Biomedicals, Santa Ana, CA, USA). For rhizoplane DNA extraction, microbial cells were dislodged and collected from 5 g of roots and subjected to DNA extraction using the FastDNA® Spin Kit for Soil. All endophytic tissue samples were extracted by the improved CTAB method. The chimeric sequences such as host mitochondria and chloroplasts were removed, and the V5–V7 hypervariable regions of the 16S rDNA of the microorganism were amplified by using nested primers 799 F (5′-AACMGGATTAGATACCCKG-3′) and 1392R (5′-ACGGGCGGTGTGTRC-3′) and 799 F (5′-AACMGGATTAGATACCCKG-3′) and 1193R (5′-ACGTCATCCCCACCT-TCC-3′). The purified PCR products were sequenced on an Illumina Miseq PE300 platform.

2.3. Bioinformatics Analysis Process

Firstly, the raw data were subjected to quality control and filtering using Trimmomatic software. Next, the effective data were assembled using FLASH v.1.2.11 software (https://ccb.jhu.edu/software/FLASH/index.shtml, accessed on 14 April 2022). Operational taxonomic units (OTUs) were clustered based on 97% similarity using UPARSE v.11 software (http://www.drive5.com/uparse/, accessed on 14 April 2022). Chimeras were removed using UCHIME v.11 software (https://drive5.com/uchime, accessed on 14 April 2022). Finally, the RDP classifier was employed to annotate species classification for each sequence by comparing it with the Silva database.

2.4. Statistical Analysis

A bipartite network analysis of the bacterial community associated with the V. natans was performed using the QIIME script and visualized in Gephi v.0.9.2. Nodes in the network correspond to sample and bacterial OTUs, and the link indicates the presence of OTUs in the sample. For beta diversity, the Bray–Curtis dissimilarity distance matrices were used to perform a principal coordinates analysis (PCoA). Alpha diversity indices were calculated using the mothur v.1.30.2 software (https://www.moth-ur.org/wiki/Download_mothur, accessed on 14 April 2022).
A co-occurrence network analysis was performed for each bacterial community associated with the different plant compartments of V. natans to explore the significant relations among the OTUs, using the routine CoNet in Cytoscape v3.9.1. To build the network, we filtered out the OTUs with frequencies less than 0.05 then combined an ensemble of the Pearson and Spearman correlation coefficients and the Bray–Curtis (BC) and Kullback–Leibler (KLD) dissimilarity indices. Based on the topological position of the microbial community within the ecological network, several potential keystones were identified. The topological type of each node was determined using measures of within-model connectivity (Zi) and among-module connectivity (Pi). Based on the values of Zi and Pi, four categories were established: peripherals (Zi < 2.5, Pi < 0.62), connectors (Zi < 2.5, Pi > 0.62), module hubs (Zi > 2.5, Pi < 0.62), and network hubs (Zi > 2.5, Pi > 0.62). The latter three categories played a crucial role in network topology and were deemed key taxa. Network analysis was depicted using Gephi v.0.9.2. The biogeochemical cycle functionality of microorganisms was analyzed using FAPROTAX v2.2 (https://pages.uoregon.edu/slo-uca/LoucaLab/archive/FAPROTAX/lib/php/index.php, accessed on 14 April 2022), a database that was manually constructed to map prokaryotic taxa, such as genera or species, to ecologically related functions. These functions include nitrification, denitrification, or fermentation, based on literature from culture representatives. Using a Python script, the sample OTU table was converted into the following six tables: sample menu, process report, OTUs for each functional group, table with each functional group overlapping, annotation information used in the analysis process, and sub-table of the input sample OTU table (only OTUs related to a specific function are listed).

3. Results and Discussion

3.1. OTUs Distribution and Diversity

Bacterial diversity was assessed using Illumina MiSeq sequencing of the 16S rRNA gene in five distinct plant compartments, including sediment, rhizosphere, rhizoplane, root endosphere, and leaf endosphere. A total of 833 unique operational taxonomic units (OTUs, cut-off level 97%) were identified at an even sequencing depth of 24,100 sequences per sample. The highest number of different OTUs was found in sediment (774 OTUs), followed by rhizoplane (747 OTUs), rhizosphere (730 OTUs), root endosphere (535 OTUs), and leaf endosphere (325 OTUs). A bipartite network analysis showed that the sediment harbors the highest number of OTUs, which are shared, in part, with the root-associated compartments (Figure 1A). Among the sediment OTUs (774), 56.33% were shared between the rhizosphere and rhizoplane, while only 32.69% were in common with the root and leaf tissues. PCoA based on the Bray–Curtis distance revealed that PCoA1 explained 49.92% and PCoA2 19.98% of the total variance among the different communities (PERMANOVA: R2 = 0.878, p < 0.001; Figure 1B). The results suggest that the plant compartments exerted a significant influence on the bacterial community, with distinct samples exhibiting significant clustering. To investigate the effect of plant compartments on bacterial communities, α-diversity indices of bacterial communities in different samples were calculated (Figure 1C–E; Table 1). In all samples, the value of the coverage index was greater than 0.99, indicating that the sequencing results covered almost all bacteria in the samples (Figure 1C). In Figure 1D, bacterial richness (Chao1 index) was significantly greater in sediment (755.950 ± 7.417) than in other plant compartments (rhizosphere, 722.446 ± 28.159; rhizoplane, 709.794 ± 54.305; root endosphere, 515.799 ± 44.392; leaf endosphere, 311.312 ± 51.223). The Shannon index in the sediment samples had the highest value (5.455 ± 0.028), whereas the lowest values (2.546 ± 0.444) were in the leaf endosphere samples (Figure 1E). Contrary to the Chao1 index, the Shannon index increased gradually from rhizosphere (3.882 ± 0.121) to rhizoplane (3.969 ± 0.124) and then to root endosphere (4.050 ± 0.350). The alpha-diversity results showed that the bacterial communities were significantly different in different plant compartments.
To investigate the impact of plant compartments on the bacterial communities, we analyzed the distribution of OTUs and diversity indexes across different plant compartments. We examined the effect of V. natans on sediment by comparing the sediment with the rhizosphere. The OTUs number and α-diversity indexes were found to be higher in the sediment, consistent with previous studies that reported significantly lower bacterial community diversity associated with the rhizosphere of submerged macrophytes than that of the sediment [23]. Yan et al. used V. natans for rhizodegradation to repair PAH-contaminated sediments, which showed that the Shannon indices gradually decreased from the sediment bacterial community to the endosphere and then to the rhizoplane [24]. These results indicate a significant influence of V. natans on sediment. Moreover, it has been demonstrated that sediment is the primary source of the plant-associated microbiome, and the microbiome in all plant compartments is tightly regulated by a rigorous process [25]. This rigorous process is mainly caused by the plant immune defense mechanism; microorganisms must overcome it to survive and colonize in different plant compartments.
Most of the observed endosphere bacterial community originates from the rhizosphere, which is consistent with the two-step selection model proposed for the filtration of Phragmites australis root endophytic microbiome, where microorganisms colonize the rhizosphere and then select the roots [26]. This filtration effect between different root compartments is due to niche characteristic differences and the selection of the microbial community [14]. Compared to other plant compartments, there are few studies on bacterial communities in submerged plant leaves, and whether the filtration mechanism of bacterial communities in root compartments can also be applied to the filtration mechanism of bacterial communities in leaves has yet to be determined [18]. Since most of the endophytic bacteria originated from the sediment and shared a certain similarity with the sediment microbial community, the diversity of root endophytic bacteria was higher than that of leaf tissue, which might be attributed to the fact that the root system was the main site for the interaction between plants and the sediment [26]. Most bacteria found in leaves were also found in the root endosphere, indicating that the root endosphere may be the major source of bacteria in leaves. This is likely because endophyte bacteria can spread systemically inside the plant and colonize the leaves [27]. However, this may not be the only pathway for the leaf endosphere, as submerged plants are completely submerged in water, which contains a large number of microorganisms that can enter the plant directly via the leaf stomata. The results of this and previous studies support the idea that submerged macrophytes recruit and select bacteria instead of randomly obtaining bacteria from surrounding environments.

3.2. Taxonomic Composition Analysis at Different Levels

A total of 26 phyla and 71 classes were identified, including some unclassified bacteria (Figure 2). In Figure 2A, the predominant phyla included Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Desulfobacterota, and Chloroflexi, which could be detected in all samples with varied values. The most dominant phylum in the six samples was Proteobacteria, which accounted for 33.12~86.80% (mean 61.07%) of the relative abundance value, followed by Bacteroidetes (1.65~45.10%, mean 17.46%), Actinobacteria (1.06~23.99%, mean 6.92%), Firmicutes (1.14~9.59%, mean 5.84%), Desulfobacterota (0.02~9.37%, mean 2.25%), and Chloroflexi (0.02~8.20%, mean 2.13%). We observed significant enrichment (p < 0.05) of Proteobacteria in plant-associated communities (mean 68.06%), compared with sediment (33.12%). Bacteroidetes were significantly enriched (p < 0.05) in the leaf endosphere (45.10%) compared with other samples (mean 10.56%). In addition, Actinobacteriota, Desulfobacterota, and Chloroflexi were present at higher levels in sediment than in other plant compartments. The relative abundance of Firmicutes in endophytic tissues was lower than that in sediment, rhizosphere, and rhizoplane (p < 0.05).
At the class level, there was a significant difference in the relative abundance of dominant classes across plant compartments, as depicted in Figure 2B (p < 0.05). Within Proteobacteria, Gammaproteobacteria (22.43~54.15%; mean 36.54%) and Alphaproteobacteria (10.69~45.99%; mean 24.53%) exhibited absolute dominance in all samples, with their combined relative abundance exceeding 60%. However, Bacteroidia (45.10%) were the most dominant class in the leaf endosphere, whereas their abundance was less than 1% in sediment. Desulfuromonadia (5.96%) and Thermoleophilia (5.04%) were also the predominant dominant classes in sediment, whereas their abundance was less than 1% in other samples.
The hierarchical clustering heat map (Figure 3) drawn according to the top 30 genera of total abundance in five samples exhibited three groups. The sediment samples were separately aggregated, while the rhizosphere and rhizoplane samples, as well as the endophytic bacteria samples from the roots and leaves, were grouped. The results indicate that the endophytic bacterial samples (from roots and leaves) and the rhizosphere and rhizoplane samples have similar bacterial community compositions. Species distributions differed greatly across samples. Intrasporangium, Gaiella, Ramlibacter, Ellin6067, Bacillus, and Methylocystis were mainly distributed in sediment, while the relative abundance of these genera in other samples was low. Species belonging to Rhizobium, Flavobacterium, Pseudomonas, Massilia, and Polaromonas were particularly enriched in the rhizosphere and rhizoplane, as well as in root and leaf endospheres, when compared to sediment samples. More globally, the relative abundance of Limnohabitans was high in all samples. The relative abundance of unclassified_f_Comamonadaceae was higher in the root endosphere samples than in the other samples.
Although the results suggest that plant compartments have varying effects on the distribution of operational taxonomic units (OTUs) and α-diversity indices of bacterial communities, a comprehensive understanding of the relationship between plant compartments and bacterial communities requires an examination of the impact of plant compartments on the composition of bacterial communities. Our investigation revealed that each sample contained a unique set of bacteria. Proteobacteria play a crucial role in organic matter decomposition and cycling and are closely related to plant growth. They were dominant in all compartments related to V. natans in a similar way to other submerged plants such as Potamogeton crispus [28], Hydrilla verticillatea [29], and Potamogeton lucens [30]. However, we observed that the relative abundance of Alphaproteobacteria and Gammaproteobacteria in root-related compartments was significantly higher than that in sediment and leaf endophytes. Alphaproteobacteria and Gammaproteobacteria, as r-strategy organisms, can utilize a wide range of carbon substrates from roots, which explains their dominance in roots [31]. Several genera belonging to the phylum Proteobacteria were involved in various biogeochemical processes that control nutrient fluxes in aquatic ecosystems. For instance, Limnohabitans, Rhizobium, Pseudomonas, and Massilia, which are affiliated with Proteobacteria, were found to be positively related to improving metabolic capacity, promoting plant growth, degrading pollutants, and enhancing plant disease resistance [32,33,34] and were widely distributed in the plant compartments, except in sediment, of V. natans. The relative abundance of Bacteroidetes in the leaf endosphere was higher than that in other compartments. These bacteria can decompose soluble organic matter composed of macromolecules, which are widely involved in the biogeochemical cycle and nitrogen removal of freshwater ecosystems [35]. Flavobacterium, the most dominant genus of the phylum Bacteroidetes, is a highly abundant aerobic denitrifying microorganism that can obtain energy for its growth by degrading organic matter from algae such as cyanobacteria in the presence of sufficient light [36]. The relative abundance of Actinobacteria and Chloroflexi phyla was significantly higher in the sediment. These bacteria are key to nitrogen removal in the management of wastewater [37] and are involved in the degradation of plant-derived compounds such as cellulose [38]. Strikingly, the dominant genus Bacillus in the phylum Firmicutes affects plant growth through biological nitrogen fixation and the production of indole-3-acetic acid (IAA) and gibberellins (GA) [39]. However, Bacillus was not detected in root and leaf endosphere samples, suggesting that it fails to enter the internal plant tissues of its host after rhizosphere colonization. Together, these findings confirmed that plant-associated bacterial communities were enriched via tissue selection based on the character of the niche and the adaptation of bacteria to habitat conditions.

3.3. Specific and Shared Bacterial Assemblages

The Venn diagram showed specific OTUs for each plant compartment and a cluster of shared OTUs (Figure 4). The sediment samples contained 28 specific OTUs, accounting for 3.36% of the total OTUs, while other samples contained no specific OTUs. A total of 249 shared OTUs were present in the five samples; 15 of these OTUs were relatively abundant (>1%) (Figure 4A). At the genus level, OTU915 was annotated as Pseudomonas, accounting for 11.95% (Figure 4B); OTU548, OTU845, and OTU1436 were annotated as Flavobacterium, accounting for 20.50%; OTU975 and OTU543 were annotated as Rhizobium, accounting for 16.41%; OTU1144 and OTU1593 were annotated as Limnohabitans, accounting for 6.14%; OTU1903 was annotated as Massilia, accounting for 4.32%; OTU2741 was annotated as Intrasporangium, accounting for 2.72%; OTU1604 was annotated as Polaromonas, accounting for 2.24%; and OTU1073 was annotated as Duganella, accounting for 1.63%. OTU868 and OTU1001 could not be assigned to a specific genus and were only annotated as Oxalobacteraceae and Comamonadaceae, accounting for 1.26% and 1.01%, respectively. In brief, the shared OTUs (249) were mainly affiliated with Gammaproteobacteria, Alphaproteobacteria, Bacteroidetes, and Actinobacteria, accounting for 66.44% of the relative abundance of all samples sequenced. There were 682 shared OTUs for sediment and rhizosphere samples and 285 shared OTUs for plant compartments. All but 10 of the OTUs in the root and leaf endospheres were also present in the rhizosphere and rhizoplane. By contrast, 227 OTUs were found in the root endosphere but not in leaves; only 17 OTUs were present in leaves but not in the root endosphere. These results indicated that the plant compartments have a significant effect on the distribution of bacterial OTUs.
By defining shared OTUs from different plant compartments in overlapping areas of Venn diagrams, we found that the dominant phyla/classes were mainly Gammaproteobacteria, Alphaproteobacteria, Bacteroidetes, and Actinobacteria. These common groups were defined as the key groups shared among the microbial communities in V. natans. The predominance of these phyla or classes was consistent with their ability to use a variety of root carbon matrices from plants, promote plant growth, and degrade pollutants [31,35,37]. These dominant bacteria were also found in other studies, which assessed the composition of rhizosphere bacterial communities of two emergent plants (Phragmites australis and Triarrhena lutarioriparia) [25]. These dominant groups may have significantly contributed to many major ecosystem processes, including playing a role in organic matter flux and biomass production [40].

3.4. Co-Occurrence Network of Bacterial Communities in Different Plant Compartments

An analysis of the co-occurrence bacterial networks within different samples of V. natans showed different connectivity patterns directly and differently influenced by the plant compartments (Figure 5A). According to the total number of nodes and edges, the network complexity gradually decreased from sediment to rhizosphere, rhizoplane, root endosphere, and finally leaf endosphere (Table 2). We recorded a significantly higher number of co-occurrence interactions than mutual exclusions in all networks. The classification of the most connected nodes in each network is significantly different (Table 3). Actinobacteria mainly shaped the topology of the bacterial network in the sediment (up to 24.93% of the total degree of connection), followed by Gammaproteobacteria (up to 22.74%) and Alphaproteobacteria (up to 10.40%). The connectivity between rhizosphere and rhizoplane bacterial communities was mainly driven by Gammaproteobacteria (56.14% and 43.95%, respectively). In the root endosphere, connectivity among the bacterial community was primarily driven by Alphaproteobacteria (46.98%) and Gammaproteobacteria (42.19%), without Desulfobacterota. By contrast, Bacteroidetes mainly shaped the topology of the bacterial network in the leaf endosphere (46.66%), followed by Alphaproteobacteria (26.70%) and Gammaproteobacteria (22.88%), but without Desulfobacterot and Chloroflexi.
Subsequently, we identified the potential keystone taxa in the networks (peripheral nodes, connectors, module hubs, and network hubs) (Figure 5B). In this study, the majority of nodes in each network are peripheral nodes, without network hubs. Nodes belonging to module hubs and connectors are defined as keystone species. A total of 14 key groups (OTUs) were detected in the networks, including 13 connectors and 1 module hub. Specifically, four connectors were found in the sediment network, whereas only one module hub was detected in the rhizosphere network. The maximum number of keystones was found in the rhizoplane with six connectors. The root and leaf endospheres contain two and one connectors, respectively. Most of these module hubs and connectors were affiliated with four phyla or classes: Gammaproteobacteria, Alphaproteobacteria, Firmicutes, and Actinobacteria. These results suggested that plant compartments significantly affect the network structure and topological properties of individual OTUs and key microbial populations.
Co-occurrence network analysis can provide a better understanding of interactions within functional bacterial communities and enable the assessment of the topological roles of each taxon in the network [41]. Our study’s network analysis revealed distinct structural characteristics in various plant compartments, and the differentiation of these compartments’ niches may account for the differences in topological traits. The endophytic bacterial community’s niche differentiation is low due to the host plant’s space and nutrient limitations. Similarly, sediment microbial communities exhibit niche differentiation, which may result from root deposition effects, leading to the selection of specific bacteria combinations [42]. The networks of V. natans’ various samples were positively correlated, indicating that most bacterial colonies had extensive cooperative and symbiotic potential in their respective compartments. Extensive cooperation among bacterial groups can significantly impact ecological processes and host adaptation [43].
Research has demonstrated that hubs or connectors can act as keystone groups in stabilizing ecosystems [44]. The five networks of V. natans were found to have varying numbers of module hubs and connectors. The rhizosphere network had only one module hub, indicating a higher level of orderliness compared to other networks. It has been established that more module hubs can maintain and stabilize microbial community structure [45]. The rhizoplane network had the highest number of connectors, allowing it to organize a series of modules into a complete community. This enhances the efficiency of energy metabolism, nutrient cycling, and material transformation in the environment [46]. Alphaproteobacteria, Gammaproteobacteria, Actinobacteria, and Firmicutes were identified as key OTUs in the five networks. As mentioned earlier, these key phyla or classes play crucial ecological roles within bacterial communities. These findings highlight the significant influence of keystone species on bacterial community structure and functioning [8].

3.5. Functional Prediction of Bacterial Communities in Different Plant Compartments

Many microorganisms are involved in crucial biogeochemical processes and interspecies interactions. The putative functions of FAPROTAX are mainly used to analyze the functions of biogeochemical cycles of microorganisms, such as the circulatory functions of sulfur, carbon, hydrogen, and nitrogen. According to the classification annotation of 16S rRNA sequences, 833 OTUs were classified into 28 functional groups (Figure 6). Although the most common functions of bacterial communities in different plant compartments were aerobic chemoheterotrophy, chemoheterotrophy, and nitrogen fixation, there was no significant difference in relative abundance between the three functional types (p > 0.05; Figure 7). Chemoheterotrophy and aerobic chemoheterotrophy were mainly contributed by the abundant Proteobacteria. A comparison of the different compartments shows that the functional populations related to human and animal diseases (such as animal parasites or symbionts, human pathogens (all), and human pathogens (pneumonia)), ureolysis, and nitrate reduction were more abundant in the rhizosphere, while methylotrophy, methanol oxidation, methanotrophy, and hydrocarbon degradation were higher in the endosphere. These results indicate that the bacterial community in the V. natans samples mainly participates in biogeochemical functions, where the carbon cycle, nitrogen cycle, and fermentation are rich in autotrophic functions (phototrophy, photoautotrophy, oxygenic photoautotrophy), as well as in animal parasites or symbionts, human pathogens (all), methylotrophy, methanol oxidation, methanotrophy, and hydrocarbon degradation.
To further explore the relationship between plant compartments and microbial communities, we compared the functional potential of the bacterial communities in different plant compartments. FAPROTAX can effectively predict the variations of bacterial ecological functional profiles [47]. The dominant functions were chemoheterotrophy and aerobic chemoheterotrophy involved in C cycling, which is consistent with the research on the main potential functions of other aquatic plants [48]. Heterotrophic bacteria are often used as decomposers and are responsible for in situ pollution repair and degradation of organic matter in ecosystems [49]. This observation agrees with the idea that most bacteria have evolved flexible and diverse strategies to use available carbon sources [50]. Nitrogen is one of the key indicators of water quality in lake ecosystems, and its excessive emissions will lead to global problems, such as eutrophication, water deterioration, and even harm to human health [51]. Nitrogen fixation is the primary reaction in the biogeochemical cycling of nitrogen [52], and it is also one of the most important functions of the bacterial community related to V. natans. Similar studies reported that the sediment and epiphytic bacterial community of V. natans had important functions in nitrogen fixation [53]. Therefore, the functions of V. natans-related bacterial communities were mainly concentrated within the lake biochemical cycle, especially the cycle of carbon and nitrogen.

4. Conclusions

In summary, we observed significant effects of plant compartments on plant-related microorganisms. These effects extend to all aspects of the bacterial community, including the OTU distribution, diversity, community composition, co-occurrence networks, and the mutualism between the potential function of the microbiome and different plant compartments of the V. natans. There were significant differences between sediment and other plant compartments, although two endosphere samples (root and leaf endospheres) had some similarities. It was proven that different plant compartments of V. natans greatly influenced its microbial community in sediment. Our network and functional prediction analyses revealed different co-occurrence patterns of bacterial communities in different plant compartments and the relative abundance of functional groups. Although the application value of V. natans in water environment improvement and ecological restoration is gradually being explored, the technology based on V. natans for aquatic plant restoration inevitably has certain limitations. This study provides theoretical and technical data for the future work of isolating beneficial bacteria from plants.

Author Contributions

Conceptualization: L.W. and Z.Z.; data curation: L.W.; investigation: L.W., S.C. and S.Z.; methodology: Z.Z.; writing—original draft: L.W.; writing—review and editing: L.W. and Z.Z.; funding acquisition: Z.Z. All authors have read and agreed to the published version of the manuscript.

Funding

Major Science and Technology Project of Anhui Province in 2020: Research and application of key technologies for ecological restoration of heavily polluted water bodies based on micro-nano aeration, efficient microorganisms, and aquatic plants (No. 202003a07020014).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author Zhongze Zhou. E-mail: [email protected].

Acknowledgments

We would like to thank Shuyi Chen for her assistance with data collection. We would also like to thank Zhongze Zhou and Xiaoxin Ye for their advice on the paper. We would like to thank Marci Baun from the University of California, Los Angeles, for editing the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Plett, J.M.; Martin, F.M. Know your enemy, embrace your friend: Using omics to understand how plants respond differently to pathogenic and mutualistic microorganisms. Plant J. 2018, 93, 729–746. [Google Scholar] [CrossRef]
  2. Gao, T.; Shi, X.Y. Taxonomic structure and function of seed-inhabiting bacterial microbiota from common reed (Phragmites australis) and narrowleaf cattail (Typha angustifolia L.). Arch. Microbiol. 2018, 200, 869–876. [Google Scholar] [CrossRef]
  3. Srivastava, J.K.; Chandra, H.; Kalra, S.J.S.; Mishra, P.; Khan, H.; Yadav, P. Plant–microbe interaction in aquatic system and their role in the management of water quality: A review. Appl. Water Sci. 2016, 7, 1079–1090. [Google Scholar] [CrossRef]
  4. Lyu, Y.; Huang, R.; Zeng, J.; Wu, Q. Aquatic Macrophytes and Local Factors Drive Bacterial Community Distribution and Interactions in a Riparian Zone of Lake Taihu. Water 2020, 12, 432. [Google Scholar] [CrossRef]
  5. Chen, P.; Zhao, M.; Tang, F.; Hu, Y.; Peng, X.; Shen, S. The effect of plant compartments on the Broussonetia papyrifera-associated fungal and bacterial communities. Appl. Microbiol. Biotechnol. 2020, 104, 3627–3641. [Google Scholar] [CrossRef]
  6. Muller, D.B.; Vogel, C.; Bai, Y.; Vorholt, J.A. The Plant Microbiota: Systems-Level Insights and Perspectives. Annu. Rev. Genet. 2016, 50, 211–234. [Google Scholar] [CrossRef]
  7. Edwards, J.; Johnson, C.; Santos-Medellin, C.; Lurie, E.; Podishetty, N.K.; Bhatnagar, S.; Eisen, J.A.; Sundaresan, V. Structure, variation, and assembly of the root-associated microbiomes of rice. Proc. Natl. Acad. Sci. USA 2015, 112, E911–E920. [Google Scholar] [CrossRef]
  8. Hu, S.; He, R.; Wang, W.; Zhao, D.; Zeng, J.; Huang, R.; Duan, M.; Yu, Z. Composition and co-occurrence patterns of Phragmites australis rhizosphere bacterial community. Aquat. Ecol. 2021, 55, 695–710. [Google Scholar] [CrossRef]
  9. Wang, B.; Zheng, X.; Zhang, H.; Yu, X.; Lian, Y.; Yang, X.; Yu, H.; Hu, R.; He, Z.; Xiao, F.; et al. Metagenomic insights into the effects of submerged plants on functional potential of microbial communities in wetland sediments. Mar. Life Sci. Technol. 2021, 3, 405–415. [Google Scholar] [CrossRef]
  10. Singh, A.; Vyas, D.; Malaviya, P. Two-stage phyto-microremediation of tannery effluent by Spirodela polyrrhiza (L.) Schleid. and chromium resistant bacteria. Bioresour. Technol. 2016, 216, 883–893. [Google Scholar] [CrossRef]
  11. Emilson, E.J.S.; Carson, M.A.; Yakimovich, K.M.; Osterholz, H.; Dittmar, T.; Gunn, J.M.; Mykytczuk, N.C.S.; Basiliko, N.; Tanentzap, A.J. Climate-driven shifts in sediment chemistry enhance methane production in northern lakes. Nat. Commun. 2018, 9, 1801. [Google Scholar] [CrossRef] [PubMed]
  12. Cregger, M.A.; Veach, A.M.; Yang, Z.K.; Crouch, M.J.; Vilgalys, R.; Tuskan, G.A.; Schadt, C.W. The Populus holobiont: Dissecting the effects of plant niches and genotype on the microbiome. Microbiome 2018, 6, 31. [Google Scholar] [CrossRef] [PubMed]
  13. Mangeot-Peter, L.; Tschaplinski, T.J.; Engle, N.L.; Veneault-Fourrey, C.; Martin, F.; Deveau, A. Impacts of Soil Microbiome Variations on Root Colonization by Fungi and Bacteria and on the Metabolome of Populus tremula × alba. Phytobiomes J. 2020, 4, 142–155. [Google Scholar] [CrossRef]
  14. Xiao, X.; Chen, W.; Zong, L.; Yang, J.; Jiao, S.; Lin, Y.; Wang, E.; Wei, G. Two cultivated legume plants reveal the enrichment process of the microbiome in the rhizocompartments. Mol. Ecol. 2017, 26, 1641–1651. [Google Scholar] [CrossRef] [PubMed]
  15. Song, Z.; Wang, H.; Strong, P.J.; Li, Z.; Jiang, P. Plant impact on the coupled terrestrial biogeochemical cycles of silicon and carbon: Implications for biogeochemical carbon sequestration. Earth-Sci. Rev. 2012, 115, 319–331. [Google Scholar] [CrossRef]
  16. Andrade, G.V.S.; Rodrigues, F.A.; Nadal, M.C.; Dambroz, C.M.S.; Martins, A.D.; Rodrigues, V.A.; Ferreira, G.M.R.; Pasqual, M.; Buttros, V.H.; D’oria, J. Plant-endophytic bacteria interactions associated with root and leaf microbiomes of Cattleya walkeriana and their effect on plant growth. Sci. Hortic. 2023, 309, 111656. [Google Scholar] [CrossRef]
  17. Orozco-Mosqueda, M.D.C.; Santoyo, G. Plant-microbial endophytes interactions: Scrutinizing their beneficial mechanisms from genomic explorations. Curr. Plant Biol. 2021, 25, 100189. [Google Scholar] [CrossRef]
  18. Afzal, I.; Shinwari, Z.K.; Sikandar, S.; Shahzad, S. Plant beneficial endophytic bacteria: Mechanisms, diversity, host range and genetic determinants. Microbiol. Res. 2019, 221, 36–49. [Google Scholar] [CrossRef]
  19. Compant, S.; Duffy, B.; Nowak, J.; Cle´ment, C.; Barka, E.A. Use of plant growth-promoting bacteria for biocontrol of plant diseases: Principles, mechanisms of action, and future prospects. Appl. Environ. Microbiol. 2005, 71, 4951–4959. [Google Scholar] [CrossRef]
  20. Peng, H.; Ge, D.; Yuan, G.; Zou, D.; Fu, H.; Jeppesen, E. Effect of clonal fragmentation on the growth of Vallisneria natans (Lour.) Hara at contrasting nutrient and light conditions. Hydrobiologia 2021, 848, 903–912. [Google Scholar] [CrossRef]
  21. Huang, S.; Song, Q.; Li, Q.; Zhang, H.; Luo, X.; Zheng, Z. Damage of heavy metals to Vallisneria natans (V. natans) and characterization of microbial community in biofilm. Aquat. Toxicol. 2020, 225, 105515. [Google Scholar] [CrossRef] [PubMed]
  22. Li, J.; Yang, L.; Wu, Z. Toxicity of chlortetracycline and oxytetracycline on Vallisneria natans (Lour.) Hare. Environ. Sci. Pollut. Res. Int. 2021, 28, 62549–62561. [Google Scholar] [CrossRef] [PubMed]
  23. Zhao, D.; Wang, S.; Huang, R.; Zeng, J.; Huang, F.; Yu, Z. Diversity and composition of bacterial community in the rhizosphere sediments of submerged macrophytes revealed by 454 pyrosequencing. Ann. Microbiol. 2017, 67, 313–319. [Google Scholar] [CrossRef]
  24. Yan, H.; Yan, Z.; Wang, L.; Hao, Z.; Huang, J. Toward understanding submersed macrophyte Vallisneria natans-microbe partnerships to improve remediation potential for PAH-contaminated sediment. J. Hazard. Mater. 2022, 425, 127767. [Google Scholar] [CrossRef]
  25. Zhang, X.; He, R.; Su, R.; Zeng, J.; Zhou, Q.; Huang, R.; Zhao, D.; Guo, L.; He, F.; Yu, Z. Composition and co-occurrence network of the rhizosphere bacterial community of two emergent macrophytes and implications for phytoremediation. Mar. Freshw. Res. 2021, 72, 1053–1064. [Google Scholar] [CrossRef]
  26. He, R.; Zeng, J.; Zhao, D.; Huang, R.; Yu, Z.; Wu, Q. Contrasting Patterns in Diversity and Community Assembly of Phragmites australis Root-Associated Bacterial Communities from Different Seasons. Appl. Environ. Microbiol. 2020, 86, e00379-20. [Google Scholar] [CrossRef]
  27. Compant, S.; Clément, C.; Sessitsch, A. Plant growth-promoting bacteria in the rhizo- and endosphere of plants: Their role, colonization, mechanisms involved and prospects for utilization. Soil Biol. Biochem. 2010, 42, 669–678. [Google Scholar] [CrossRef]
  28. Chen, W.M.; Tang, Y.Q.; Mori, K.; Wu, X.L. Distribution of culturable endophytic bacteria in aquatic plants and their potential for bioremediation in polluted waters. Aquat. Biol. 2012, 15, 99–110. [Google Scholar] [CrossRef]
  29. Feng, J.; Cui, B.; Yuan, B.; Zhang, L.; Zhang, J.; Zhang, A.; Han, X.; Pan, L. Purification Mechanism of Low-Pollution Water in Three Submerged Plants and Analysis of Bacterial Community Structure in Plant Rhizospheres. Environ. Eng. Sci. 2020, 37, 560–571. [Google Scholar] [CrossRef]
  30. Shi, L.; Xia, P.; Lin, T.; Li, G.; Wang, T.; Du, X. Temporal Succession of Bacterial Community Structure, Co-occurrence Patterns, and Community Assembly Process in Epiphytic Biofilms of Submerged Plants in a Plateau Lake. Microb. Ecol. 2023, 85, 87–99. [Google Scholar] [CrossRef]
  31. Philippot, L.; Raaijmakers, J.M.; Lemanceau, P.; Putten, W.H.V.D. Going back to the roots: The microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 2013, 11, 789–799. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, N.; Wu, M.; Che, Y.; Kong, Y.; Wang, Q.; Sha, W.; Gong, Z.; Zhou, J. Effects of shining pondweed (Potamogeton lucens) on bacterial communities in water and rhizosphere sediments in Nansi Lake, China. Environ. Sci. Pollut. Res. Int. 2022, 29, 51665–51673. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, C.; Liu, S.; Zhang, Y.; Liu, B.; Zeng, L.; He, F.; Zhou, Q.; Wu, Z. Effects of Planted Versus Naturally Growing Vallisneria natans on the Sediment Microbial Community in West Lake, China. Microb. Ecol. 2017, 74, 278–288. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, H.; Hao, B.; Cai, Y.; Liu, G.; Xing, W. Effects of submerged vegetation on sediment nitrogen-cycling bacterial communities in Honghu Lake (China). Sci. Total Environ. 2021, 755, 142541. [Google Scholar] [CrossRef]
  35. Yang, Y.; Cheng, K.; Li, K.; Jin, Y.; He, X. Deciphering the diversity patterns and community assembly of rare and abundant bacterial communities in a wetland system. Sci. Total Environ. 2022, 838, 156334. [Google Scholar] [CrossRef]
  36. Chang, W.; Sun, J.; Pang, Y.; Zhang, S.; Gong, L.; Lu, J.; Feng, B.; Xu, R. Effects of different habitats on the bacterial community composition in the water and sediments of Lake Taihu, China. Environ. Sci. Pollut. Res. Int. 2020, 27, 44983–44994. [Google Scholar] [CrossRef]
  37. Wang, Q.; Xie, H.; Ngo, H.H.; Guo, W.; Zhang, J.; Liu, C.; Liang, S.; Hu, Z.; Yang, Z.; Zhao, C. Microbial abundance and community in subsurface flow constructed wetland microcosms: Role of plant presence. Environ. Sci. Pollut. Res. Int. 2016, 23, 4036–4045. [Google Scholar] [CrossRef]
  38. Wang, C.; Liu, S.; Zhang, Y.; Liu, B.; He, F.; Zhou, Q.; Wu, Z. Bacterial Communities and Their Predicted Functions Explain the Sediment Nitrogen Changes Along with Submerged Macrophyte Restoration. Microb. Ecol. 2018, 76, 625–636. [Google Scholar] [CrossRef]
  39. Liu, Y.; Zou, Y.; Kong, L.; Bai, G.; Liu, Z.; Wang, C.; Ding, Z.; He, F.; Wu, Z.; Zhang, Y. Effects of bentonite on the growth process of submerged macrophytes and sediment microenvironment. J. Environ. Manag. 2021, 287, 112308. [Google Scholar] [CrossRef]
  40. Xue, Y.; Chen, H.; Yang, J.R.; Liu, M.; Huang, B.; Yang, J. Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom. ISME J. 2018, 12, 2263–2277. [Google Scholar] [CrossRef]
  41. Jiao, C.; Zhao, D.; Zeng, J.; Guo, L.; Yu, Z. Disentangling the seasonal co-occurrence patterns and ecological stochasticity of planktonic and benthic bacterial communities within multiple lakes. Sci. Total Environ. 2020, 740, 140010. [Google Scholar] [CrossRef] [PubMed]
  42. Peng, M.; Wang, C.; Wang, Z.; Huang, X.; Zhou, F.; Yan, S.; Liu, X. Differences between the effects of plant species and compartments on microbiome composition in two halophyte Suaeda species. Bioengineered 2022, 13, 12475–12488. [Google Scholar] [CrossRef] [PubMed]
  43. Layeghifard, M.; Hwang, D.M.; Guttman, D.S. Disentangling Interactions in the Microbiome: A Network Perspective. Trends Microbiol. 2017, 25, 217–228. [Google Scholar] [CrossRef]
  44. Shi, S.; Nuccio, E.E.; Shi, Z.J.; He, Z.; Zhou, J.; Firestone, M.K. The interconnected rhizosphere: High network complexity dominates rhizosphere assemblages. Ecol. Lett. 2016, 19, 926–936. [Google Scholar] [CrossRef] [PubMed]
  45. Qi, G.; Ma, G.; Chen, S.; Zhao, X. Microbial Network and Soil Properties Are Changed in Bacterial Wilt-Susceptible Soil. Appl. Environ. Microbiol. 2019, 85, e00162-19. [Google Scholar] [CrossRef]
  46. Wan, X.; Gao, Q.; Zhao, J.; Feng, J.; Nostrand, J.D.V.; Yang, Y.; Zhou, J. Biogeographic patterns of microbial association networks in paddy soil within Eastern China. Soil Biol. Biochem. 2020, 142, 107696. [Google Scholar] [CrossRef]
  47. Zheng, W.; Wen, X. How exogenous influent communities and environmental conditions affect activated sludge communities in the membrane bioreactor of a wastewater treatment plant. Sci. Total Environ. 2019, 692, 622–630. [Google Scholar] [CrossRef]
  48. Liu, Y.; Guo, Z.; Zhang, P.; Du, J.; Gao, P.; Zhang, Z. Diversity and Structure of Vegetation Rhizosphere Bacterial Community in Various Habitats of Liaohekou Coastal Wetlands. Sustainability 2022, 14, 16396. [Google Scholar] [CrossRef]
  49. Wei, Z.; Liu, Y.; Feng, K.; Li, S.; Wang, S.; Jin, D.; Zhang, Y.; Chen, H.; Yin, H.; Xu, M.; et al. The divergence between fungal and bacterial communities in seasonal and spatial variations of wastewater treatment plants. Sci. Total Environ. 2018, 628–629, 969–978. [Google Scholar] [CrossRef]
  50. Gorke, B.; Stulke, J. Carbon catabolite repression in bacteria: Many ways to make the most out of nutrients. Nat. Rev. Microbiol. 2008, 6, 613–624. [Google Scholar] [CrossRef]
  51. Lucchetti, R.; Onotri, L.; Clarizia, L.; Natale, F.D.; Somma, I.D.; Andreozzi, R.; Marotta, R. Removal of nitrate and simultaneous hydrogen generation through photocatalytic reforming of glycerol over “in situ” prepared zero-valent nano copper/P25. Appl. Catal. B Environ. 2017, 202, 539–549. [Google Scholar] [CrossRef]
  52. Linz, A.M.; He, S.; Stevens, S.L.R.; Anantharaman, K.; Rohwer, R.R.; Malmstrom, R.R.; Bertilsson, S.; McMahon, K.D. Freshwater carbon and nutrient cycles revealed through reconstructed population genomes. PeerJ 2018, 6, e6075. [Google Scholar] [CrossRef] [PubMed]
  53. Zhu, H.Z.; Jiang, M.Z.; Zhou, N.; Jiang, C.Y.; Liu, S.J. Submerged macrophytes recruit unique microbial communities and drive functional zonation in an aquatic system. Appl. Microbiol. Biotechnol. 2021, 105, 7517–7528. [Google Scholar] [CrossRef] [PubMed]
Figure 1. OTU distribution and diversity of bacterial communities in samples from different compartments (sediment, rhizosphere, rhizoplane, root endosphere, and leaf endosphere). (A) Bipartite network analysis represented sample/OTU interactions. Different circle colors indicate different sample types of V. natans; the OTU node is silver, with the edges connecting sample nodes to OTU. Nodes are colored by sample type. (B) Based on the Bray–Curtis distance, principal coordinates analysis (PCoA) was conducted for bacterial communities in different plant compartments. (CE) α–diversity index of the bacterial community in different samples.
Figure 1. OTU distribution and diversity of bacterial communities in samples from different compartments (sediment, rhizosphere, rhizoplane, root endosphere, and leaf endosphere). (A) Bipartite network analysis represented sample/OTU interactions. Different circle colors indicate different sample types of V. natans; the OTU node is silver, with the edges connecting sample nodes to OTU. Nodes are colored by sample type. (B) Based on the Bray–Curtis distance, principal coordinates analysis (PCoA) was conducted for bacterial communities in different plant compartments. (CE) α–diversity index of the bacterial community in different samples.
Diversity 15 00676 g001
Figure 2. At different taxonomy levels, community composition and relative abundance of dominant bacterial groups in different samples of V. natans. (A,B) represent the community composition of dominant bacterial groups at phylum and class levels, respectively. (C,D) represent the relative abundance of dominant bacterial groups at phylum and class levels, respectively. (* p < 0.05, ** p < 0.01).
Figure 2. At different taxonomy levels, community composition and relative abundance of dominant bacterial groups in different samples of V. natans. (A,B) represent the community composition of dominant bacterial groups at phylum and class levels, respectively. (C,D) represent the relative abundance of dominant bacterial groups at phylum and class levels, respectively. (* p < 0.05, ** p < 0.01).
Diversity 15 00676 g002
Figure 3. Heat map graph of hierarchy cluster for the top 30 genera. The relative abundance values for the bacteria genera are depicted by the color intensity with the legend indicated at the right of the figure. Brick red and blue, respectively, represent high or poor enrichment of a genus.
Figure 3. Heat map graph of hierarchy cluster for the top 30 genera. The relative abundance values for the bacteria genera are depicted by the color intensity with the legend indicated at the right of the figure. Brick red and blue, respectively, represent high or poor enrichment of a genus.
Diversity 15 00676 g003
Figure 4. Microbial community shared by different plant compartments of Vallisneria natans. (A) Venn diagrams showing the shared and specified bacterial OTUs of different plant compartments. (B) The phylogenetic trees and distribution of the 14 most abundant (>1%) shared OTUs in all selected different plant compartments of Vallisneria natans.
Figure 4. Microbial community shared by different plant compartments of Vallisneria natans. (A) Venn diagrams showing the shared and specified bacterial OTUs of different plant compartments. (B) The phylogenetic trees and distribution of the 14 most abundant (>1%) shared OTUs in all selected different plant compartments of Vallisneria natans.
Diversity 15 00676 g004
Figure 5. The correlation networks of bacteria in different plant compartments. (A) co-occurrence network of bacterial communities associated with different plant compartments. The colors and sizes of network nodes correspond to bacterial phylum/class and the degree of connectivity of each node. The colors of the network links represent the correlations between nodes (green: positive correlation; red: negative correlation). (B) Within-module connectivity (Zi) and among-module connectivity (Pi) plots show the distribution of OTUs based on their topological roles. Each node represents an OTU in the bacterial community. The topological role of each OTU depends on the scatter plot of Zi and Pi.
Figure 5. The correlation networks of bacteria in different plant compartments. (A) co-occurrence network of bacterial communities associated with different plant compartments. The colors and sizes of network nodes correspond to bacterial phylum/class and the degree of connectivity of each node. The colors of the network links represent the correlations between nodes (green: positive correlation; red: negative correlation). (B) Within-module connectivity (Zi) and among-module connectivity (Pi) plots show the distribution of OTUs based on their topological roles. Each node represents an OTU in the bacterial community. The topological role of each OTU depends on the scatter plot of Zi and Pi.
Diversity 15 00676 g005
Figure 6. Based on the FAPROTAX database, the potential functions of bacterial communities in different plant compartments were displayed by heatmap.
Figure 6. Based on the FAPROTAX database, the potential functions of bacterial communities in different plant compartments were displayed by heatmap.
Diversity 15 00676 g006
Figure 7. The top ten functional groups in relative abundance in different samples. (* p < 0.05).
Figure 7. The top ten functional groups in relative abundance in different samples. (* p < 0.05).
Diversity 15 00676 g007
Table 1. Mean bacterial α-diversity indexes in the different compartments.
Table 1. Mean bacterial α-diversity indexes in the different compartments.
CompartmentOTUs NumberCoverageChao1Shannon
Sediment7740.998755.950 ± 7.4175.455 ± 0.028
Rhizosphere7300.994722.446 ± 28.1593.882 ± 0.121
Rhizoplane7470.994709.794 ± 54.3053.969 ± 0.124
Root endosphere5350.996515.799 ± 44.3924.050 ± 0.350
Leaf endosphere3250.997311.312 ± 51.2232.546 ± 0.444
An OTU was defined as reads with 97% sequence similarity. Values represent an average number of α-diversity indexes (±standard error).
Table 2. Bacterial co-occurrence network characteristics in each compartment niche.
Table 2. Bacterial co-occurrence network characteristics in each compartment niche.
NicheNodesEdgesAverage DegreeModularity aAverage Clustering Coefficient bAverage Path Distance c
Sediment516343713.3220.9030.77310.901
Rhizosphere197301230.5790.6690.8064.673
Rhizoplane185295833.0590.6430.7953.443
Root endosphere154288337.4420.5510.8042.628
Leaf endosphere8678118.1630.5700.7682.773
a: Degree of nodes tending to differentiate into different network modules. b: Degree of nodes tending to cluster together. c: Network path distance is the length of the shortest path between two nodes within the network.
Table 3. In that histogram, the relative abundance of the node with connectivity was reported at the phylum/class level for different samples.
Table 3. In that histogram, the relative abundance of the node with connectivity was reported at the phylum/class level for different samples.
Phylum/ClassSedimentRhizosphereRhizoplaneRoot EndosphereLeaf Endosphere
Alphaproteobacteria10.40%13.08%26.05%46.98%26.70%
Gammaproteobacteria22.74%56.14%43.95%42.19%22.88%
Bacteroidetes1.50%17.91%15.90%8.12%46.66%
Actinobacteria24.93%3.26%3.63%1.23%0.79%
Firmicutes6.99%8.01%9.36%0.60%2.97%
Desulfobacterota9.75%0.61%0.09%00
Chloroflexi6.95%0.40%0.61%0.32%0
Others16.73%0.59%0.41%0.56%0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wan, L.; Zhang, S.; Zhou, Z.; Chen, S. Plant Compartments Shape the Assembly and Network of Vallisneria natans-Associated Microorganisms. Diversity 2023, 15, 676. https://doi.org/10.3390/d15050676

AMA Style

Wan L, Zhang S, Zhou Z, Chen S. Plant Compartments Shape the Assembly and Network of Vallisneria natans-Associated Microorganisms. Diversity. 2023; 15(5):676. https://doi.org/10.3390/d15050676

Chicago/Turabian Style

Wan, Linqiang, Siyong Zhang, Zhongze Zhou, and Shuyi Chen. 2023. "Plant Compartments Shape the Assembly and Network of Vallisneria natans-Associated Microorganisms" Diversity 15, no. 5: 676. https://doi.org/10.3390/d15050676

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop