Community structures of mangrove endophytic and rhizosphere bacteria in Zhangjiangkou National Mangrove Nature Reserve

Bacterial communities play an important role in mangrove ecosystems. In order to gain information on the bacterial communities in mangrove species and rhizospheres grown in Zhangjiangkou National Mangrove Nature Reserve, this study collected root, branch, and leaf samples from five mangrove species as well as rhizosphere and non-rhizosphere samples and analyzed the community structure of endophytic bacteria and bacteria in rhizosphere and non-rhizosphere using Illumina high-throughput sequencing technique. Bacteria in 52 phyla, 64 classes, 152 orders, 295 families, and 794 genera were identified, which mainly belonged to Proteobacteria, Cyanobacteria, Actinobacteria, Firmicutes, Bacteroidetes, Fusobacteria, and Nitrospirota. At each taxonomic level, the community structure of the rhizosphere bacteria varied slightly with mangrove species, but endophytic bacteria differed greatly with plant species. The diversity indices of endophytic bacteria in branch and leaf samples of Acanthus ilicifolius were significantly lower, and endophytic bacteria in the plant tissues had higher abundance in the replication/repair and translation Clusters of Orthologous Genes functional categories but lower abundance in the carbohydrate metabolism category. This study helps to understand the community structure and diversity characteristics of endophytic and rhizosphere bacteria in different mangrove plants. Provide a theoretical basis for in-depth research on the functions of mangrove ecosystems.

www.nature.com/scientificreports/ the microbial communities in mangrove plants consist of bacteria [12][13][14] , which play an important role in nutrient transformation in mangrove ecosystems 3 .
The objectives of this study were to analyze the composition and diversity of endophytic bacteria in five representative mangrove species as well as their rhizospheres bacterial communities by means of Illumina high throughput sequencing 15 and to test if endophytic and rhizosphere bacterial community composition differences were associated with plant species in this Reserve.

Sequencing results
The V5-V7 regions of bacterial 16S rRNA gene were amplified using the DNA prepared from the roots, branches, and leaves of the trees of five mangrove species as well as rhizosphere and non-rhizosphere samples.Following high-throughput sequencing, a total of 4,956,926 valid sequences were obtained, of which the average length was 383 bp, and the library coverage was 0.942-0.996.The data shown in Fig. 1 revealed that as the number of sequences increased, the rarefaction curves for different samples gradually reached a plateau, indicating that in this study the number of sequences was close to saturation and the sequencing data were reasonable and capable of reflecting the structure of most of the bacterial communities in the samples.

Composition of bacterial communities
A total of 52 phyla, 64 classes, 152 orders, 295 families, and 794 genera of bacteria were identified in 63 samples collected from five mangrove species and related soils.Figure 2 shows the number of common OTUs shared among different samples.The roots (Fig. 2A), branches (Fig. 2B), leaves (Fig. 2C), and rhizospheres (Fig. 2D), related to five mangrove species shared 211, 138, 131, and 1212 common OTUs, respectively.Among them, soil samples had the largest number of shared OTUs, followed by root samples.Branch and leaf samples had a similar number of common OTUs.The different numbers of common OTUs from samples reflected the differences among them, specifically there were small differences between the soil samples but distinct differences between the branch and leaf samples for all five mangrove species.
At the phylum level (Fig. 3A), the bacteria in root, branch, leaf, and rhizosphere of the five mangrove species as well as non-rhizosphere mainly included the Proteobacteria, Cyanobacteria, Actinobacteria, Firmicutes, Bacteroidetes, Fusobacterium, Nitrospirae, Germmatimonadetes, Acidobacteria, and Chloroflexi phyla, of which the Proteobacteria was the dominant phylum.The second most dominant bacterial phylum in root sample ZQ.A of K. candel was Firmicutes (20.82%), and the second most dominant phylum in the roots of the other four species of mangrove plants was Actinobacteria.The second most dominant phylum in branch sample ZQ.B of K. candel and    gymnorrhiza was Burkholderiaceae (9.56%).For all the soil samples regardless of rhizosphere or non-rhizosphere, the identified bacteria at the family level had low abundance and low difference in community compositions.
At the genus level (Fig. 3E), the dominant bacterial genera and proportions in the tissue samples of mangrove plants were as follows: 16 For the leaf and soil samples collected from the forests of different mangrove species, the bacteria identified at the genus level had low abundance and there were no distinct differences in bacterial abundance among the samples.

Alpha diversity of bacterial communities
The alpha diversity indices calculated based on the number of OTUs of the root, branch, leaf, and soil samples collected from the five mangrove species are presented in Table 1, which were highly variable with plant species and plant organs, but less variable among rhizospheres.The inter-group differences of the alpha diversity indices among the samples are shown in Table 2.The difference between root sample ZB.A of A. marina and root sample ZL.A of A. ilicifolius was significant in the number of observed species, Chao1 index, and phylogenetic diversity (PD) index, and highly significant in the Shannon index (Table 2).There was a significant difference in Shannon Boxplots of Alpha-diversity indices are present in Fig. 4. Observed species (Fig. 4A) and Chao1 index (Fig. 4B) reflected the OUT abundance in samples.Shannon (Fig. 4C) and PD (Fig. 4D) indices represented the diversity of OUT in the samples.The greater the Chao1, the higher the expected species richness of microbiota.The greater the Shannon, the higher the diversity of microbiota.

Beta diversity of bacterial communities
PCoA analysis was performed based on the number of OTUs of the root, branch, leaf, and soil samples related to five mangrove species (Fig. 5).Clustering was performed to construct an unweighted pair group method with arithmetic mean (UPGMA) cluster tree of the samples (Fig. 6).The results indicated that the distance between the soil samples collected from the forests of five mangrove species was relatively small, and the UPGMA cluster tree also showed that there was no distinct difference between the soil samples, suggesting that all the soil samples had similar bacterial compositions and that the mangrove species had little effect on the soil bacterial community.

Prediction of bacterial functions
The functions of the bacteria in the root, branch, leaf, and soil samples collected from the five mangrove species were predicted.The top 10 functional categories with the highest abundances at Clusters of Orthologous Genes (COG) level 2 were used to generate Fig. 7 showing the relative abundances of COG functional categories in the samples.The top 10 COG functional categories, namely, membrane transport, amino acid metabolism, carbohydrate metabolism, replication and repair, energy metabolism, translation, metabolism of cofactors and vitamins, cellular processes and signaling, and lipid metabolism, were roughly the same for the bacteria in all samples.There was little difference in the predicted functional categories between the samples collected from K. candel, A. marina, A. corniculatum, and B. gymnorrhiza, while the relative abundances of the replication/repair and translation categories of the bacteria in the root, branch, and leaf samples of A. ilicifolius were higher, and the relative abundance of the carbohydrate metabolism category in the A. ilicifolius tissue samples was lower compared to the other four mangrove species.
The data in Fig. 8 show that the difference between mangrove rhizosphere bacteria and endophytic bacteria was distinct in the categories of metabolism, cellular processes, and organismal systems.The endophytic bacteria in the tissue samples of A. ilicifolius had higher abundance in organismal systems, genetic information processing, and human diseases, while the endophytic bacteria of K. candel, A. marina, A. corniculatum, and B. gymnorrhiza had higher abundance in environmental information processing, metabolism, and cellular processes.The results indicate that there may be specific bacteria that affect the growth and development of the mangrove plants by a species-specific means.Rhizosphere and non-rhizosphere were analyzed.Organic matter, pH, EC, available and total nitrogen, phosphorus, and potassium among the samples appeared to be largely similar (Table 3).

Mental test statistical analysis at the genus level
To further analyze the correlation between the main mineral elements of soil and the correlation with plant rhizosphere bacteria.We performed Mantel test analysis.As shown in Fig. 9, organic matter has a very significant positive correlation with available potassium and total nitrogen, and a very significant negative correlation with total phosphorus.Available phosphorus and available potassium showed extremely significant negative correlations with the corresponding total potassium and total phosphorus respectively.In this study, the top ten bacterial genera with relative abundance were selected.The analysis results showed that the bacterial genera with significant correlation with soil chemical properties were Methylobacterium, Staphylococcus, Vibrio, Alcaligenes and Pseudomonas.The analysis showed that Methylobacterium showed a significant positive correlation with changes in pH and conductivity.There was a significant positive correlation between the changes in the abundance of Staphylococcus and the changes in available potassium content.The relative abundance of this bacterial genus was relatively high in the rhizosphere soil samples of Kandelia candel and Bruguiera gymnor-rhiza.Vibrio has the highest relative abundance (8.26%) in the rhizosphere soil samples of Bruguiera gymnor-rhiza.There was a strong correlation between the abundance changes of Alcaligenes and the changes in available phosphorus and hydrolyzable nitrogen contents.The relative abundance of this bacterial genus was the highest in the rhizosphere  soil samples of Acanthus ilicifolius.There is a significant correlation between Pseudomonas and the total potassium content, and it is mostly distributed in the rhizosphere soil samples of Avicennia marina.

Discussion
The Illumina high-throughput sequencing technology was used in this study to examine the community structure of rhizosphere bacteria and endophytic bacteria associated with five mangrove species.Endophytic bacteria were highly diverse, and the bacterial community compositions had a degree of correlation with mangrove species  www.nature.com/scientificreports/ that was indicated by the distribution of the number of OTUs among the samples (Fig. 2).These results suggest that mangrove forests are rich in microbial resources, contributing to the productivity of mangrove ecosystems.
In general, the endophytic bacterial community of mangrove plants was dominated by Alphaproteobacteria and Gammaproteobacteria.This was like the Kandelia candel habitat in the Jiulong River Estuary in my country, which also has the highest abundance of Proteobacteria.The coexisting endophytic bacteria of the halophyte Aster tripolium L. growing in the salty meadows of central Poland are dominated by Bacillus species of the Firmicutes 16,17 .The bacteria in Alphaproteobacteria have oligotrophic properties that enable them to live in low-nutrient environments.Although mangrove ecosystems are rich in organic matter, they are deficient in nutrients overall 18 .The enrichment of Alphaproteobacteria, which includes most phototrophic bacteria, such as Rhizobacteria and Rhodospirillum, helps increase the productivity of mangrove ecosystems.Consistent with the results of multiple previous studies 4,19 , the present study found that the bacterial community in the rhizosphere of mangroves was dominated by Gammaproteobacteria and Deltaproteobacteria. Deltaproteobacteria contains most of the known sulfate-reducing bacteria, which play a key role in the deposition and cycling of nitrogen, carbon, and sulfur 20 .The bacteria in Deltaproteobacteria are anaerobic, so it could be speculated that the anaerobic conditions in mangrove soils might prompt the selection of specific microbial groups 21 .
Studies have shown that mangrove ecosystems have abundant functional bacteria resources that are related to the decomposition of organic matter and energy sources, and those bacteria have a degree of association with mangrove species.Pseudomonadales was the dominant order of endophytic bacteria in B. gymnorrhiza plants in the present study.Pseudomonadales contains nitrogen-fixing bacteria 22 that can decompose glucose and are also involved in nitrate reduction 23 .Rhodospirillales and Salinisphaerales were the dominant orders of endophytic bacteria in A. marina plants.The bacteria in the Rhodospirillales order have anaerobic properties, can  utilize light to perform anoxygenic photosynthesis for growth 19 and may be able to increase the productivity of mangrove ecosystems.The bacteria in Salinisphaerales are considered to be halophilic and tolerant to high levels of salinity.It is speculated that the bacteria in Salinisphaerales can help improve the adaptability of mangroves to salt-rich environments 24 .Enterobacterales and Pseudonocardiales were found to be the dominant orders of endophytic bacteria in A. corniculatum trees.Most of the bacteria in Enterobacterales are involved in nitratereducing processes and can catalyze the oxidation of nitrite.Pseudonocardiales is a kind of Actinomycetes that may be related to disease prevention in mangrove plants.Bacillales, which can promote plant growth, was the dominant order of endophytic bacteria in K. candel plants 23 .
An important finding of this study is that there were small differences in microbial compositions between nonrhizosphere and rhizosphere and there were no specific associations between plants species and root microbes (Fig. 3).The analysis of PCoA (Fig. 5) and UPGMA cluster tree (Fig. 6) also showed that the composition of bacterial communities in the soil samples related to five mangrove species was relatively similar.This finding to some extent disapproved our hypothesis that microbial community compositions in soils should be associated with plant species.It has been well documented that in terrestrial soils plant species affect rhizosphere microbial communities 4,6,[25][26][27][28] .Conversely, based on the plant-soil feedback theory, soil microbes can significantly affect plant communities by altering soil physicochemical and trophic properties or regulating plant coexistence [29][30][31][32] .Our finding may suggest that mangrove plant species differ from those grown in terrestrial soils and have a limited impact on soil microbial communities.This is consistent with the research results of Ping et al., who believe that mangrove bacteria are more susceptible to abiotic factors such as temperature, rather than cross-domain biotic factors 33 .This could be attributed to several factors: (1) soil properties in the mangrove ecosystem could be significantly affected by water properties, and the dynamics of water movement may alter soil physical and chemical properties and substantially reduce plant rhizosphere effects.Similar results were also reported by Gong et al. 32 in which environmental conditions and historical events played an important role in shaping the bacterial communities in mangrove swamps along the coast of Beibu Gulf in Guangxi, China.(2) Five species of mangrove plants in the present study grew in similar localities in the Reserve, soil physicochemical properties and nutrient status may be similar (Table 3).A long-time adaptation to similar soil conditions may diminish the effects of plant species to soil microbes.However, there was a difference in the composition of the endophytic bacterial communities among the five species of mangrove plants, indicating that the mangrove species had a direct effect on the bacterial community structure.Thus, it is tempting to speculate that soil properties in the area where samples were collected in mangrove could play more important role than plant species in influencing bacterial compositions in the rhizosphere.Further studies in this and other mangrove nature reserve are warranted to testify this proposition.

Conclusions
This study compared the dominant bacteria related to five species of mangrove plants and their rhizospheres in Zhangjiangkou National Mangrove Nature Reserve and analyzed the functions of the dominant bacteria.Our study showed that composition of bacterial communities in soil samples of the five mangrove species was relatively similar, which could be an indication that environmental conditions where mangrove plant growth may have substantial effects on soil microbial communities.On the other hand, there were significant differences in endophytic bacterial compositions among plant species, suggesting that plant communities can directly or indirectly influence microbial communities by altering the quantity and quality of litter 34,35 .These results provide some new information on microbial communities in the Zhangjiangkou National Mangrove Nature Reserve.As the Reserve has more than 2000 hectares, more research on microbial community compositions either as endophytes or in their rhizospheres are needed to confirm these findings.

Sample collection
The sampling took place in aforementioned Zhangjiangkou National Mangrove Nature Reserve (Fig. 10), which features a subtropical marine monsoon climate.The annual average temperature is 21.2 °C and annual precipitation is 1714 mm.In November 2020, the roots, branches, and leaves of five mangrove species (Kandelia candel, Aegicera corniculatum, Avicennia marina, Acanthus ilicifolius, and Bruguiera gymnorrhiza) as well as their rhizospheres were collected.The abbreviations for K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza were ZQ, ZT, ZB, ZL, and ZM.Root samples were labeled as ZQ.A, ZT.A, ZB.A, ZL.A, and ZM.A.Samples from branches were labeled as ZQ.B, ZT.B, ZB.B, ZL.B, and ZM.B. Leaf samples were stamped as ZQ.C, ZT.C, ZB.C, ZL.C, and ZM.C.Their rhizospheres samples were marked as ZQ.D, ZT.D, ZB.D, ZL.D, and ZM.D. The soil collected from the same area with without any roots as the bulk soil, which was used as control and labeled as ZCK.D.
For sample collection, healthy and robust trees were randomly selected at the sampling sites, and their branches and leaves were collected.The roots in the soil layer depth of 20-40 cm and within 1 m distance from the trunk of each tree were collected.The soil particles less than 1 cm in size on the roots were collected as rhizosphere samples after removing impurities, such as plant residues, roots, and small stones.Three individual trees per species were randomly selected for sampling, Thus, there were 60 root, branch, leaf, and soil samples with additional three non-rhizosphere samples collected away from any trees.The collected samples were packed and sealed in sterilized bags immediately and stored at − 20 °C within 24 h after collection.
The collected soil samples, both rhizosphere and non-rhizosphere ones, were analyzed according to the methods described by Sun et al. 36 for organic matter contents, pH and EC levels, and concentrations of available and total nitrogen, phosphorus, and potassium.Genomic DNAs from root, branch, root, and soil samples were extracted using the Plant Genomic DNA Extraction Kit produced by Tiangen Biotech (Beijing, China) according to the manufacturer's instructions.The concentration and purity of extracted genomic DNA were detected using agarose gel electrophoresis, and the qualified DNA samples were stored at − 20 °C for later use.Polymerase chain reaction (PCR) primers were designed based on the conserved regions of the bacterial 16S rRNA gene, and sequencing adapters were added to the ends of the primers.PCR amplification was performed, and the PCR products were purified, quantified, and homogenized to form a sequencing library.The constructed libraries were subjected to quality control.Qualified libraries were used to construct small fragment libraries for sequencing based on the Illumina NovaSeq sequencing platform and paired-end sequencing method.The sequencing was performed by Novogene (Beijing, China).

Data analysis
The raw reads obtained through sequencing were processed to obtain valid sequences that were then subjected to operational taxonomic unit (OTU) clustering and the analysis of species at a 97% similarity level 37 .The rarefaction curve and hierarchical clustering of the samples were plotted using R software.The number of observed species, Shannon index, Ace index, library coverage, and phylogenetic diversity index were computed using Qiime software 38 .R software was used in the plotting of principal coordinates analysis (PCoA) results.PICRUSt analysis was used to predict the metabolic functions of the bacterial populations based on the OTU tree in the Greengene database as well as the gene information included in the OTU tree 39 .

Determination and analysis of soil physical and chemical indicators
Determination of soil pH using potentiometric method.Soil organic matter was measured using potassium dichromate oxidation-external heating method 40 .Total N was determined using the Kjeldahl method.Hydrolyzable N was determined using the alkaline hydrolysis-diffusion method.Total P was determined using the alkali fusion-molybdenum antimony colorimetric method.The available P was measured using hydrochloric acid-sulfuric acid solution extraction method 40 .Total K was determined using alkali fusion-flame photometry.The available K was determined using ammonium acetate extraction-flame photometry.The soil conductivity was measured using the electrode method 40 .Mental test is analyzed in the "vegan" package of R (version 4.1.2) 41.

Figure 1 .
Figure 1.Rarefaction curves for plant samples of five mangrove species and their rhizospheres and nonrhizosphere at the 97% similarity level, where ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Figure 2 .
Figure 2. Distribution of the number of operational taxonomic units (OTUs) among the plant tissue samples of five mangrove species and soil samples (A: root sample; B: branch sample; C: leaf sample; and D: soil sample), where ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively.

Figure 3 .
Figure 3. Relative abundances of the bacterial populations in the plant tissue samples of five mangrove species and soil samples.For the group name, ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample (A: phylum level; B: class level; C: order level; D: family level; E: genus level).

Figure 4 .
Figure 4. Boxplots of the diversity indices of the bacteria in the plant tissue samples of five mangrove species.(A) number of observed species; (B) Chao1 index; (C) Shannon index; (D) phylogenetic diversity index.For the label of group, ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Figure 5 .
Figure 5. Principal coordinates analysis (PCoA) of the plant tissue samples of five mangrove species and their rhizospheres samples based on the unweighted UniFrac distance.Where ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Figure 6 .
Figure 6.Unweighted pair group method with arithmetic mean (UPGMA) cluster tree of the plant tissue samples of five mangrove species and their rhizospheres based on the unweighted UniFrac distance.

Figure 7 .
Figure 7. Histogram of the relative abundances of Clusters of Orthologous Genes (COG) functional categories of the dominant bacteria in the plant tissue samples of five mangrove species and their rhizospheres.Where ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Figure 8 .
Figure 8.Heat map of bacterial gene functional prediction using the PICRUSt algorithm.Each row refers to a KEGG pathway, and each column is a sample where ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Figure 10 .
Figure 10.The sites where plant, rhizosphere, and soil samples were collected at Zhangjiangkou National Mangrove Nature Reserve.Roots, branches, and leaves of five mangrove species (Kandelia candel, Aegicera corniculatum, Avicennia marina, Acanthus ilicifolius, and Bruguiera gymnorrhiza) as well as their rhizospheres were collected (the satellite imagery was used map software: Gugu Map Downloader (official version V1.05); URL: http:// www.gggis.com).Three plants were randomly selected per species.The abbreviations for K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza were ZQ, ZT, ZB, ZL, and ZM.Root samples were labeled as ZQ.A, ZT.A, ZB.A, ZL.A, and ZM.A.Samples from branches were labeled as ZQ.B, ZT.B, ZB.B, ZL.B, and ZM.B. Leaf samples were stamped as ZQ.C, ZT.C, ZB.C, ZL.C, and ZM.C.Their rhizosphere samples were marked as ZQ.D, ZT.D, ZB.D, ZL.D, and ZM.D. The soil collected from the same area with without any roots as the bulk soil, which was used as control and labeled as ZCK.D.
At the order level, the differences in bacterial community composition in the soil samples from different mangrove species again were small, and the majority of the bacteria were in the unidentified Gammaproteobacteria class and Rhizobiales.Vibrionales accounted for 8.86% of the bacteria in the soil samples related to B. gymnorrhiza, which was the highest proportion compared with other four mangrove species.For bacteria at the family level (Fig.3D), the dominant bacterial families were Burkholderiaceae (19.54%) and Planococcaceae(13.49%)in root sample ZQ.A of K. candel.The prevailing bacterial families in root sample ZT.A of A. corniculatum were Beijerinckiaceae (53.45%) and Pseudonocardiaceae (14.46%).The abundant bacterial family in root sample ZB.A of A. marina was Vibrionaceae (17.67%).The dominant bacterial families in root sample ZM.A of B. gymnorrhiza were Burkholderiaceae (11.70%) and Vibrionaceae (23.10%).Beijerinckiaceae was the most abundant in branch sample ZQ.B of K. candel, accounting for 47.75% of bacterial families.The dominant bacterial family in branch sample ZT.B of A. corniculatum was Beijerinckiaceae (18.29%).The prevailing bacterial families in branch sample ZT.B of A. marina were Beijerinckiaceae (11.42%) and Salinisphaeraceae (23.73%).The dominant bacterial family in branch sample ZM.B of B. gymnorrhiza was Staphylococcaceae (28.58%).Beijerinckiaceae (26.27%) was the most abundant bacterial family in leaf sample ZQ.C of K. candel.The dominant bacterial families in leaf sample ZT.C of A. corniculatum were Enterobacteriaceae (38.86%) and Staphylococcaceae (11.86%).The prevailing bacterial families in leaf sample ZB.C of A. marina were Salinisphaeraceae (18.06%) and Rhodospirillaceae (24.48%).The dominant bacterial family in leaf sample ZM.C of B.

Table 1 .
Alpha diversity index of five mangrove plant and their rhizospheres samples.For the samples, ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.Vol:.(1234567890)Scientific Reports | (2023) 13:17127 | https://doi.org/10.1038/s41598-023-44447-2www.nature.com/scientificreports/index between root sample ZB.A of A. marina and the root sample ZT.A of A. corniculatum.The branch samples between ZB.B of A. marina and ZL.B of A. ilicifolius had significant differences in the number of observed species.Significant differences occurred between ZB.B of A. marina and ZQ,.B of K. candel in the number of observed species and PD index.ZL.B of A. ilicifolius and ZM.B of B. gymnorrhiza showed significantly high in the number of observed species and the phylogenetic diversity index, and highly significant in the Shannon index.The difference between branch samples of ZL.B (A. ilicifolius) and ZT.B (A. corniculatum) was highly significant in the Shannon index.Branch samples between ZQ.B of K. candel and ZM.B of B. gymnorrhiza as well as between ZQ.B of K. candel and ZB.B of A. marina had significant differences in the number of observed species and the phylogenetic diversity index.The difference between leaf sample ZL.C of A. ilicifolius and the leaf samples of the other four mangrove species was significant in the number of observed species, Chao1 index, Shannon index, and phylogenetic diversity index.In addition, the significant difference in soil samples occurred between ZB.D of A. marina and ZM.D of B. gymnorrhiza only in the number of observed species and PD index.

Table 2 .
Inter-group differences of Alpha diversity index of five mangrove species and their rhizospheres samples.ZQ, ZT, ZB, ZL, and ZM represent K. candel, A. corniculatum, A. marina, A. ilicifolius, and B. gymnorrhiza, respectively and A, B, C, and D after the species symbol indicate root, branch, leaf, and soil samples, and ZCK.D is non-rhizosphere sample.

Table 3 .
pH, electrical conductivity (EC), organic matter, available and total nitrogen, phosphorus, and potassium in rhizosphere and non-rhizosphere sampled from Zhangjiangkou National Mangrove Nature Reserve.