Skip to main content

Metagenomic analysis reveal the phytoremediation effects of monocropping and intercropping of halophytes Halogeton glomeratus and Suaeda glauca in saline soil of Northwestern China

Abstract

Aims

Planting halophytes is a widely used method of phytoremediation for saline soils. The succulent halophytes Halogeton glomeratus and Suaeda glauca are widely used for remediation of saline soil in the arid region of Northwestern China. However, whether intercropping of H. glomeratus and S. glauca can increase the improvement effect for saline soil is yet to be proved.

Materials and methods

Therefore, this study analyzed three phytoremediation planting modes: monocropping of H. glomeratus (Hg), monocropping of S. glauca (Sg), and H. glomeratus and S. glauca intercropping (Hg||Sg). These were applied in field experiments, with biomass and soil physicochemical properties measured for each treatment, and the mechanism was analyzed using macrogenomics.

Results

After harvesting the halophytes after one season, the Hg treatment had the highest dry biomass and soil total dissolved salt content was reduced; correspondingly, soil pH were decreased and soil organic matter content were increased. The results showed that Actinobacteria, Acidobacteria and Proteobacteria were the dominant phylum under the four treatments. This suggests that Hg treatment was more capable of producing microorganisms favorable to saline soil remediation.

Conclusions

Thus, H. glomeratus monocropping is a more effective phytoremediation strategy for saline soil in the dry zone of Northwestern China.

Graphical Abstract

Peer Review reports

Introduction

Soil is an important resource for human survival, but soil salinization is a growing problem globally, leading to soil degradation and seriously affecting agricultural production [2]. The total area of saline soil globally is about 1 billion hectares, of which China has about 100 million hectares, accounting for more than one-tenth of the world’s total area, making it the third largest country in the world in terms of saline land distribution [4, 51, 61]. In Northwestern China, salinity and drought usually occur together, and have become one of the main factors restricting agricultural development in the arid regions of Northwestern China [56]. Due to its special soil conditions, the utilization and improvement of saline soils has been an important topic in agricultural science research [64]. Earlier studies have shown that reducing soil salinity by planting halophytes to improve saline soil is one of the most effective ways to solve the problem of soil salinization and significantly improve quality of saline soils [33, 62]. In recent years, saline soil improvement methods centered on phytoremediation have gradually gained attention and have progressed [1, 16, 21, 25, 26, 37, 38, 40]. The impact of halophytes on soil microorganisms, however, is not clear.

Macrogenomics is a technique to study microbial diversity and community structure and function, which reveals the microbial ecological characteristics of an environment by sequencing and analyzing all microbial genes in that environment [35, 42, 48, 55, 57]. A deep understanding of the soil physicochemical properties and microbial diversity after planting halophytes, as well as the isolation and screening of salt-alkali-tolerant microbial communities, is of significance for improving saline soil. As the most active part of soil, microorganisms exhibit high response characteristics to soil material components, physicochemical properties, and the microenvironment [46]. In the process of vegetation succession and restoration, soil microorganisms may play a series of roles in the plant rhizosphere [47]. Rhizosphere microorganisms, as regulators of nutrient transformation and transport between soil and roots, are the soil microbial community most affected by plants [6]. Their role in plant adaptation has been recognized, and many microorganisms beneficial for plants to cope with salt stress have been isolated and reapplied to saline soil to improve crop yield [45]. Therefore, clarifying the dominant species and functions of bacterial diversity in the inter-root soil of saline plants can be used to improve the quality of saline soil microenvironments by utilizing salt-tolerant microorganisms in the inter-root soil of saline plants [15, 56].

In Northwestern China, Halogeton glomeratus and Suaeda glauca are the dominant plants used for salt removal. Halogeton glomeratus is a dicotyledonous herbaceous plant in the Chenopodiaceae family and has strong drought and salt tolerance. Its most significant feature is that the plant stems and leaves are highly fleshy, gradually changing salinized “waste fields” into tillable “treasure fields,” and is a characteristic plant resource in the arid region of Northwestern China. Suaeda glauca, a member of the Amaranthaceae family, is a pioneer species of dicotyledonous halophyte used for ecological phytoremediation of saline and alkaline soils [10, 12]. Planting halophytes to improve saline soil is a proven practical and feasible approach. Studies have shown that cotton–halophyte intercropping can reduce soil salinity and improve crop production [33]. The content of salt and sodium ions (Na+) was decreased by intercropping with S. salsa and Zea mays L [53]. In moderately saline soils, intercropping of Arthrocaulon macrostachyum L. and tomatoes or planting tomatoes after harvest of a previous crop of A. macrostachyum L. reduced soil salinity and increased the yield of tomatoes [20]. Both H. glomeratus and S. glauca have been widely used to improve saline soil as they are highly tolerant to both salt and alkali. It is generally believed that intercropping of H. glomeratus and S. glauca can improve the coverage of saline soil and improve the soil, but whether this method can improve saline soil more effectively has not been reported.

In this experiment, we took H. glomeratus and S. glauca planted on saline and alkaline land in Northwestern China as research objects; analyzed the mechanism of saline soil improvement using macrogenomic methods; compared the differences in plant growth ability, soil physicochemical properties, and the effect of saline soil improvement of the two plants under different cultivation modes (monocropping and intercropping); and verified the results using potting experiments. The study aimed to (1) investigate the impact of three remediation methods on the physicochemical properties and microbial communities of the rhizosphere in saline soil; (2) compare monocropping and intercropping of H. glomeratus and S. glauca in saline soil improvement; and (3) analyze the relationship between the metagenome and environmental factors. The aim of this research is to enhance the development and utilization of saline plant resources by comparing the different effects of monocropping and intercropping on saline soil improvement, and to provide a scientific basis for saline soil improvement.

Materials and methods

Experimental design and sampling collection

The field trial was performed on saline soils in Northwestern China (39°13′N, 100°17′E; Fig. 1A) [32]. This region has a temperate continental desert climate with low precipitation, receiving an average annual rainfall of 117.2 mm. In November 2021, seeds of H. glomeratus and S. glauca were collected from the saline flat of Hetan Town, Huining County, Baiyin City, Gansu Province (36°12′N, 104°97′E), and cleaned to remove impurities and stored in a low temperature seed cupboard at 4℃. The H. glomeratus and S. glauca seeds were sown in May 2022 in saline soils. In November 2022, the saline soils were sampled. A total of 10 rows of monocropping of H. glomeratus (Hg), monocropping of S. glauca (Sg), and H. glomeratus and S. glauca intercropping (Hg||Sg) were planted in plots of 5 m × 10 m (Fig. 1B). Seeds were planted at 10-cm intervals in each plot. Plots without plants were used as controls (CK) [28]. Clinical trial number: not applicable.

Fig. 1
figure 1

Geographic location of saline land in Northwestern China (A). Planting pattern diagram (B): (a) schematic diagram (left) and sample diagram (right) of a H. glomeratus monocropping system; (b) schematic diagram (left) and sample diagram (right) of a S. glauca monocropping system; and (c) schematic diagram (left) and sample diagram (right) of H. glomeratus and S. glauca intercropping system. Seedling growth of Hg, Sg, and Hg||Sg for salt concentration of 0, 50, 100, 150, and 200 mmol·L− 1 (C). Abbreviations: Hg, H. glomeratus monocropping; Sg, S. glauca monocropping; Hg||Sg, H. glomeratus and S. glauca intercropping

After harvesting the Hg and Sg, four treatment areas were selected to collect topsoil samples from saline soils. In each plot, sampling was repeated three times. Each composite sample consisted of 5-cm-diameter soil cores taken to a depth of 40 cm, giving a total of 12 soil samples and three replicates for each of the four treatments. After debris removal, soil samples were filtered through a 2-mm sieve, then each soil sample was divided into two parts, and soil samples from the same sampling area were crushed and mixed evenly into sterile aseptic sealed bags, which were stored in a foam cryopreservation box [8, 60]. The first sample was air-dried to determine the physical and chemical properties of the soil and the second sample was stored in liquid nitrogen at − 80 °C for the extraction of DNA from soil microbes. After planting on saline-alkali soil for one season, 3 samples of 5 m×5 m were randomly selected from each repeat of halophyte plants. After cutting, 3 samples were cleaned, separated according to different tissues and organs, placed in the oven at 105℃ for 30 min, and then dried to constant weight at 80℃ and crushed. Weigh an appropriate amount of the above-ground sample, decompose it, cool it to constant weight, dissolve it with distilled water, and transfer it to a 100 m| volumetric bottle at constant volume to prepare the liquid to be measured. Na+, K+, were determined by flame spectrophotometry. Ca2+ and Mg2+ were determined by EDTA titration. Each indicator is repeated 3 times to take the average value. In the pot experiment, we used glass bottles of diameter of 6.5 cm and Hoagland nutrient solution to prepare medium with pH 5.8. The NaCl environment was simulated using five different salt concentrations of 0, 50, 100, 150, and 200 mmol·L− 1 [44]. The seeds of Hg and Sg were evenly planted on the medium, with 16 seeds for each medium. The glass bottles were placed flat in the incubator and the biomass of seedlings in each bottle was determined after 30 days of culture.

Determination of dry biomass, soil pH, and nutrient contents

The harvested materials were dried at 60 °C for 3 days and referred to as dry biomass. After planting on saline-alkali soil for one season, 3 samples of 5 m×5 m were randomly selected from each repeat of halophyte plants. After cutting, 3 samples were cleaned, separated according to different tissues and organs, placed in the oven at 105℃ for 30 min, and then dried to constant weight at 80℃ and crushed. Determination of Na, Ca, Mg, P [3], total water-soluble salt (TDS), soil organic matter (OM) [29], pH, total nitrogen (TN), alkaline dissolved nitrogen (AN), total phosphorus (TP), available phosphorus (AP), total potassium (TK), available potassium (AK) [63].

Determination of soil microbial community

Genomic DNA was extracted using HiPure Bacterial DNA Kits (Magen, Guangzhou, China) in accordance with the manufacturer’s instructions. Genome sequencing was performed on an Illumina Novaseq 6000 sequencer. Pair-end technology (PE 150) was used [52]. Genes were predicted using MetaGeneMark (version 3.38) on the basis of the final assembly contigs (> 500 bp). Bowtie (version 2.2.5) was used to count the number of reads, and the reads were re-aligned to the predicted gene. Non-redundant genes with gene read counts > 2 were included in the final gene catalog [54]. The R package ggplot2 was used to plot the box plot showing the number of genes in each group [34].

Chao1, ACE, Shannon, and Simpson indices were calculated using the Python scikit-bio package (version 0.5.6). Comparison of α-index between groups was calculated using Kruskal–Wallis H-test in the vegan package of R. Bray–Curtis distance matrix based on gene/taxon/function abundance was generated using the vegan package. Multivariate statistical techniques including RDA (redundancy analysis) and NMDS (non-metric multidimensional scaling) of Bray–Curtis distances were calculated using the vegan package and plotted using ggplot2. Permanova and Anosim tests were calculated using the vegan package. Heatmap graphs were plotted using the R package pheatmap. Venn analysis between groups was performed using the R project Venn Diagram package and upset plot using the UpSetR package to identify unique and common species or functions. Comparison of species/functions between groups was calculated using Welch’s t-test in the vegan package. Comparison of species/functions between groups was calculated using analysis of variance in the vegan package.

We used several complementary approaches to annotate the assembled sequences. The unigenes were annotated using DIAMOND (version 0.9.24) by alignment with those deposited in various protein databases, including the National Center for Biotechnology Information (NCBI) non-redundant protein (Nr) database and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Circular layout representations of functional gene abundance were plotted using Circos (version 0.69-3). Species comparisons between groups were calculated using Metastats (version 20090414). Differentially enriched KEGG pathways were identified according to their reporter score from the z-scores of individual Kos (KEGG orthologs). An absolute value of the reporter score = 1.96 or higher (95% confidence according to a normal distribution) was used as a detection threshold for pathways that differed significantly in abundance. Biomarker features of species and functions in each group were screened using LEfSe software (version 1.0).

Statistical analysis

Statistical analyses were performed using Microsoft Excel 2016 and SPSS25.0 (IBM, Armonk, New York, USA) with a two-way ANOVA. Evenness and Shannon diversity indices were assessed for each sample, and these indices were compared for four treatment soil samples using a Kruskal–Wallis test [39]. The relationship between soil pH and α-diversity of the soil bacterial community was analyzed using Pearson’s correlation coefficients (r). Statistical analyses were performed using IBM SPSS Statistics (version 25.0; IBM Corp., Armonk, NY, USA). The significance level was set at P ≤ 0.05 [7]. Values represent the mean ± standard error (SE) for the three repeated treatments.

Results

Changes of aboveground biomass of halophytes after planting

In the pot experiment, the Hg biomass under 0, 50, 100, 150, and 200 mmol·L− 1 NaCl treatment was the largest, consistent with the field experiment. The growth of Hg was stronger than that of Sg and Hg||Sg (Fig. 1C).

The aboveground biomass of H. glomeratus and S. glauca under three improved modes was measured (Table 1). The biomass of Hg was 12.35% and 5.68% higher than for Hg||Sg and Sg, respectively. Aboveground biomass of the different treatments differed significantly: Hg > Hg||Sg > Sg.

Table 1 Biomass of different cropping of hg, Sg, and Hg||Sg for field and pot experiments

Soil properties affected by planting different saline-tolerant plants

After planting halophytes, Na content in soil decreased significantly, including 0.732 mg/kg under Hg treatment, 0.289 mg/kg under Sg treatment, and 0.530 mg/kg under Hg||Sg treatment. Ca content in soil decreased significantly, including 10.279 mg/kg under Hg treatment, 5.097 mg/kg under Sg treatment, and 8.781 mg/kg under Hg||Sg treatment. Soil Mg content decreased significantly, including 2.125 mg/kg under Hg treatment, 1.310 mg/kg under Sg treatment, and 1.282 mg/kg under Hg||Sg treatment. However, P content in soil increased significantly, with an increase of 0.905 mg/kg under Hg treatment, 0.464 mg/kg under Sg treatment, and 0.662 mg/kg under Hg||Sg treatment. It is worth noting that there was a significant correlation between halophytes before and after seeding except for bare soil (P < 0.05).

The impact of halophytes on soil pH was related to cropping type (P < 0.05, Table 2). There were significant differences in pH between the Hg and Hg||Sg sampled areas during seedling and after harvesting (Table 2). Before and after planting halophytes, the soil pH of sample areas of Hg, Sg, and Hg||Sg decreased to 0.52, 0.19, and 0.34, respectively. For Hg, Sg, and Hg||Sg treatments, the TDS content in the soil decreased by 1.138, 0.744, and 0.923 g/kg, respectively, with the greatest decrease for Hg. Plant type significantly affected soil OM content. For Hg, Sg, and Hg||Sg treatments, the OM content in the soil of each treatment increased by 4.124, 1.113, and 2.563 g/kg, respectively. The order of OM content in the soil after different halophyte treatments follows: Hg > Hg||Sg > Sg.

The impact of halophytes on soil TN content is shown in Table 2, the highest increase in soil TN content was Hg (0.231 g/kg). Under different coverage of halophytes, the highest increase AN content in the rhizosphere soil of vegetation was for Hg (27.399 mg/kg). The order of TN and AN contents in rhizosphere soil for the halophyte treatments follows: Hg > Hg||Sg > Sg. Soil TP content is mainly influenced by soil-forming parent material and climate type, but there is also a relationship with vegetation activities. TP content only increase after Hg treatment (0.091 g/kg), because there are no statistical different among any other samples. The AP content also increased, with the greatest increase for Hg (5.297 mg/kg). After treatment with different halophytes, the TP and AP contents in the rhizosphere soil were as follows: Hg > Hg||Sg > Sg. Apart from CK, the highest reduction in soil TK content was for Hg (3.712 g/kg), while the lowest reduction was for Sg (1.196 g/kg). Compared to seedling stage, soil AK levels differed markedly at harvest stage (P < 0.05). Overall, after treatment with different halophytes, the order of TK and AK content in rhizosphere soil was Hg > Hg||Sg > Sg.

Table 2 Change of soil physicochemical properties content at sow halophytes for hg, sg and Hg||Sg in 2022

As can be seen from Table 3, the aboveground content of Na+ in halophytes was the highest under Hg, while the aboveground content in halophytes under Sg was the lowest. The aboveground part content of K+ in halophytes was the highest under Sg, and the lowest under Hg||Sg. The content of Ca2+ in the aboveground part of halophytes was the highest under Hg, and the lowest under Sg. The content of Mg2+ in the aboveground part of halophytes was the highest under Hg and the least under Sg.

Table 3 Comparison of salt ion content under three phytoremediation treatments

Changes in soil microbiota after planting different saline-tolerant plants

Metagenomic sequencing was used to analyze the characteristics and communities of soil microorganisms [61]. Metagenome sequencing produced 121,882.49 Mbp of raw data. After removing low-quality reads, 121,405.36 Mbp of clean data were obtained. After sequencing, we obtained a total of 3.2 × 105 to 4.7 × 105 contigs and 2.0 × 105 to 5.1 × 105 open reading frames (ORFs) (Tables S1 and S2). A total of six kingdom, 161 phylum, 254 class, 426 order, 754 family, 2257 genus, and 9100 species level classifications were obtained by metagenomic sequencing. According to the species annotation, the microbial composition of soil samples was mainly bacteria (over 84.11%) (Fig. 2A). Analysis of species differences at the phylum level in a Wayne diagram showed seven species unique to CK, two species unique to Hg, one species unique to Sg, and a total of 114 species in CK, Hg, Sg, and Hg||Sg, which had 160, 48, 54, and 54 unique genera, respectively; with 930 genera shared by the four treatments (Fig. 2B and C).

Fig. 2
figure 2

Kingdom-level relative abundance composition of soil microorganisms (A). Venn map of soil microbial species at phylum and genus levels in different covers (BC). Alpha diversity of soil bacterial communities (P > 0.05) (DG). The microbial community structure was shown by NMDS. Samples grouped according to different halophyte coverage methods are surrounded by 95% confidence ellipses (H). Phylum clustering map based on unweighted pair-group method (I). The abbreviations indicate the same as in Fig. 1

In this experiment, the effects of bare ground and Hg, Sg, and Hg||Sg treatments on the Chao, ACE, Shannon, and Simpson indices of microbial communities in remediated saline soils were investigated. The α-diversity index of soil microorganism communities showed significantly decreased species α-diversity (Chao 1 and ACE), and significantly increased species α-diversity (Shannon and Simpson) (Fig. 2D–G). The microbial community in the bare land (CK) significantly differed from those for the three halophyte cover treatments according to NMDS (Fig. 2H) [18]. The NMDS analysis showed that with a stress rate of 0.007, among the three phytoremediation treatments, there was a clear difference in theβ-diversity estimates of soil microorganisms proportions based on NMDS. It is noteworthy that the symbols were distributed in four quadrants each, representing the soil microbiota of CK, Hg, Sg, and Hg||Sg. The distribution variation between the CK, Hg, Sg, and Hg||Sg treatments showed a significant effect of the cover plants on soil microorganism communities.

Fig. 4The LEfSe method identifies the significantly different abundant taxa of soil bacterial communities (A). The taxa with significantly different abundances between treatments are represented by colored dots, and from the center outwards they represent the kingdom, phylum, class, order, family, genus, and species levels. Colored shadows represent trends of significantly different taxa. Only taxa with an LDA significance threshold > 3.5 are shown in this figure. Each colored point has an effect size LDA score as shown in Fig. S1. The RDA of bacterial community and environmental factors (B). Heatmap of the correlations between phyla and soil factors (C), note ** P < 0.01; *P < 0.05. The abbreviations indicate the same as in Fig. 1.

The sequencing results were annotated by species, and the UPGMA (unweighted pair-group method using arithmetic averages) of each soil sample was constructed with weighted UniFrac distance matrix at gate-level means, non-weighted group average method clustering tree (Fig. 2I). The bacterial groups with the highest relative abundance in the rhizosphere soil were Proteobacteria and Actinobacteria. Bacteroidetes, Gemmatimonadetes, and Chloroflexi were all in low abundance in the soil. The cluster diagram showed that the community structure composition of the four soil microorganisms at the phylum level was similar for Hg and Sg treatment. Although the dominant phyla in the soils of the four treatments were essentially the same, the relative abundance of individual dominant phyla differed significantly at different stages. For example, the abundances of Proteobacteria in CK, Hg, Sg, and Hg||Sg treatments were 29.63%, 25.71%, 22.21%, and 23.18%, respectively (Fig. 2I). Chloroflexi had a significant positive correlation with salt tolerance, whereas Proteobacteria had a significant negative correlation. Among the top 20 genera, the dominant phylum in CK, Hg, Sg, and Hg||Sg was Proteobacteria, with 32.86%, 22.78%, 22.76%, and 21.6%, respectively (Fig. 3A). Furthermore, the Hg, Sg, and Hg||Sg treatments had significantly higher relative abundances of Chloroflexi, Candidatus_Tectomicrobia, and Acidobacteria compared with CK.

Fig. 3
figure 3

Circos map of bacteria distribution in four soil samples at phyla level (A) and Circos map of differential KEGG functional categories (B). The abbreviations indicate the same as in Fig. 1.

The KEGG database was used to predict the gene function of non-redundant gene sets from saline soils covered by different halophytes. A total of 1,046,716 reads could correspond with the functional annotation results. There were 72,189 items in category Cellular Processes, 29,423 in category Human Diseases, and 40,508 in category Organismal Systems. There were 743,954 in Metabolism, 70,364 in Genetic Information Processing, and 90,278 in Environmental Information Processing. Among the four treatments, the most microbial functional genes belonged to Metabolism, while the least belonged to Human Diseases and Organismal Systems (Fig. 3B).

The LEfSe analyses (LDA values shown in Fig. S1) were used to check for significant differences in the abundance of soil bacterial taxa and to confirm the biological relevance of soil microbial species under different treatments (Fig. 4A). In particular, Gammaproteobacteria, Xanthomonadales, Proteobacteria, and Deltaproteobacteria were the most altered soil microorganism taxa for the CK; Betaproteobacteria, Acidobacteria, Burkholderiales, and Rhodospirillales were the most altered for Hg treatment; Sphingomonas and Gemmatimonadetes were the most altered for Sg treatment; and Actinobacteria, Chloroflexi, and Propionibacteriales were the most altered for Hg||Sg treatment.

Fig. 4
figure 4

The LEfSe method identifies the significantly different abundant taxa of soil bacterial communities (A). The taxa with significantly different abundances between treatments are represented by colored dots, and from the center outwards they represent the kingdom, phylum, class, order, family, genus, and species levels. Colored shadows represent trends of significantly different taxa. Only taxa with an LDA significance threshold > 3.5 are shown in this figure. Each colored point has an effect size LDA score as shown in Fig. S1. The RDA of bacterial community and environmental factors (B). Heatmap of the correlations between phyla and soil factors (C), note ** P < 0.01; *P < 0.05. The abbreviations indicate the same as in Fig.1

The KEGG heat map shows that for CK, the Two component system and Fatty acid metabolism are high, but decrease after planting Hg and Sg. After planting Sg, Valine, leucine, and isoleucine degradation, Propanoate metabolism, Citrate cycle (TCA cycle), Glyoxylate and dicarboxylate metabolism, Aminoacyl-tRNA biosynthesis, Purine metabolism, Carbon fixation pathways in prokaryrotes, ABC transporters, Butanoate metabolism, Alanine, aspartate, and glutamate metabolism, and Amino sugar and nucleoside sugar metabolism significantly increased. After planting Hg, there were no significant changes. After planting of Hg||Sg, Amino sugar and nucleotide sugar metabolism increased, while Fatty acid metabolism decreased. After Hg||Sg processing, Cysteine and methionine metabolism, Metabolic pathways, Oxidative physiology, Methane metabolism, Microbial metabolism in diverse environments, Quorum sensing, Carbon metabolisms, Pyruvate metabolism, Biosynthesis of amino acids, Glycolysis/gluconeogenesis, Biosynthesis of secondary metabolites, and Glycine, serine, and threonine metabolism significantly increased (Fig. S2).

Interactions between environmental factors and microbial communities

The RDA was used out to understand the effect of changes in key soil physical and chemical properties on the composition of soil microbial communities under four saline alkaline soil treatments. Two RDA axes explained 96.94% of the differences in soil microbial community structure (Fig. 4B), indicating that environmental factors significantly influenced bacterial community structure. The RDA plots used Mg, pH, TK, OM, TP, Na, AN, P, TN, TDS, and Ca as soil environmental variables (represented by blue arrows). Only Mg was present in the third quadrant, with other environmental factors concentrated in the first and fourth quadrants (Fig. 4B). A greater effect on the microbial community composition was indicated by the longer exposure to TN. A lesser effect on microbial community composition was indicated by the shorter exposure to Ca. The dominant phyla, including Proteobacteria, Acidobacteria, and Actinobacteria, were the main contributors to microbial community changes. The RDA revealed that soil microbial groups were distributed in different quadrants according to functional analysis of CK, Hg, Sg, and Hg||Sg patches. Soil OM, TP, Na, AN, P, TN, and AK were significantly correlated with Sg soil microbial functions. Soil TK and pH were significantly correlated with soil microbial functions in the Hg treatment.

The microorganisms whose components were in the top 20 with obvious differences among samples were selected, and their correlations with environmental factors analyzed. The TDS, Ca, Mg, pH, TN, AN, and AK were closely related to the abundance of dominant bacteria. The TDS content was positively correlated with Candidatus_Aenigmarchaeota (r = 0.719, P < 0.01). The Ca content was positively correlated with Candidatus_Diaphorotes (r = 0.587, P < 0.05), and negatively correlated with Candidatus_Lokiaarchaeota (r = − 0.835, P < 0.01) and with Candidatus_levy bacteria (r = − 0.596, P < 0.05). The Mg content was positively correlated with Crenarchaeota (r = 0.723, P < 0.01) and Euryarchaeota (r = 0.589, P < 0.05). Soil pH was positively correlated with Candidatus_korarchaeota (r = 0.583, P < 0.05). The TN was positively correlated with Candidatus_Heimdallarchaeota (r = − 0.741, P < 0.01), Bacteroidetes (r = − 0.821, P < 0.01), and Proteobacteria (r = − 0.889, P < 0.01). The TN was negatively correlated with Candidatus_levy bacteria (r = 0.699, P < 0.05). The AK was positively correlated with Candidatus_Aenigmarchaeota (r = 0.663, P < 0.05). The TN had the greatest influence on the change of microbial composition and quantity among the 13 indicators of soil. In the correlation heatmap, the main factors that significantly influenced the functions and relative abundances of the soil microbiome were soil Mg, pH, TK, OM, TP, Na, AN, P, TN, TDS, and Ca (Fig. 4C).

The four soil samples were grouped by screening the top 200 most abundant species at the genus level, after removing species without names; Pearson’s correlations were calculated to screen for relationships with Pearson’s correlation |r| > 0.8, P < 0.05. A co-occurrence network analysis was performed and a microbial co-occurrence network constructed for each treatment group of the four soil samples. In addition, the modular structure of each network was revealed by visualizing nine main modules within these networks (Fig. 5). Chloroflexi was the phylum with the highest number of nodes in the network in the four soil samples. In addition, Bacteroidetes and Deinococcus-Thermus had the highest number of nodes for the CK; Proteobacteria and Bacteroidetes had the highest number of nodes for the Hg treatment; Planctomycetes and Bacteroidetes nodes were the most numerous for the Sg treatment; and Bacteroidetes and Deinococcus-Thermus nodes were the most numerous for Hg||Sg treatment.

Fig. 5
figure 5

Network of co-occurring microbial genera on the basis of correlation analysis. A link represents a strong (Pearson’s correlation |r| > 0.8) and significant (P < 0.05) correlation. Node size is proportional to the relative abundance of genera; the thickness of the node connection (edge) is proportional to the value. The abbreviations indicate the same as in Fig. 1

Discussion

Effects of halophyte mulch on soil physicochemical properties

The changing of soil properties is an essential element in changing microbial communities. In general, plants suitable for desalting in salinized soils should have high salt tolerance and produce more biomass [43]. After planting cotton and alfalfa, the aboveground biomass of cotton and Suaeda salsa intercropping was higher than that of cotton and alfalfa intercropping during 2014–2016, by 20.5%, 26.3%, and 28.7%, respectively [33].In our study, Hg produced the most biomass, 9.6 t ha− 1, in salinized soils with pH 8.36 (Table 1). The field experiment and pot experiment results showed that the strength of the effect of different saline plants in improving saline soils was Hg > Hg||Sg > Sg (Table 1). Our findings are similar to those of other studies using halophytes for soil desalination.

Changes in plant growth conditions, soil pH, and nutrient levels are likely to be the main causes of changes in structure of the soil microbial community [22]. As shown in Table 2, the content of AN, AP, and AK in soil decreased when no halophytes were planted to repair saline-alkali soil (CK), because the release and conversion ability of microorganisms to these nutrients was generally reduced in saline soil [63]. On the contrary, the content of AN, AP and AK in soil increased after one season of halophyte harvest, and the effect was most obvious after halophyte treatment. This is because after phytoremediation, these nutrients are increased in the saline soil [23]. The release and conversion ability of microorganisms to these nutrients was generally reduced in saline soil. After Elaeagnus angustifolia was planted in the salt-alkali land of the Songnen Plain, soil pH and electrical conductivity decreased, while contents of TP, TN, TK, DOC, DOM, and AK increased [63]. After planting halophytes, soil pH in the Hg, Sg, and Hg||Sg sample areas decreased, and the TDS content in the halophyte rhizosphere decreased, with the greatest decrease for Hg. In particular, H. glomeratus monocropping significantly reduced soil pH and salinity compared with bare soil and the seedling stage (Table 2). Results showed that soil pH, electrical conductivity, OM, and TK were the main physical and chemical properties affecting microbial communities. Therefore, halophytes desalinate soil by absorbing salt from soil and accumulating dry matter in aboveground parts, with Hg working best.

Effects of planting halophytes on soil microbial communities

To further explore the impact of H. glomeratus and S. glauca monocropping with that of H. glomeratus and S. glauca intercropping on improving saline soil in Northwestern China, the soil microbial communities and their functions were studied using metagenomic sequencing technology. Many studies have shown that soil salinity significantly inhibits the number, composition, and diversity of most soil microorganisms [58, 65]. Results confirmed that phytoremediation treatment significantly increased the diversity and abundance of soil bacteria [57]. The α-diversity analysis revealed a decreasing trend in the ACE of the three improved Chao1 treatments, indicating that planting halophytes led to a reduction in bacterial abundance in saline-alkali soil. However, the Shannon and Simpson indices of the three improved treatments showed increasing trends, suggesting that halophyte planting increased saline microbial community diversity (Fig. 3). Feng et al. [11] showed that planting S. salsa can reduce salt stress in soil, prevent degradation of saline soils, and increase the diversity of soil bacterial communities, thus benefiting saline soils in coastal areas.

Microorganisms have different adaptability to different soil conditions [41]. Earlier studies found that Chloroflexi has desiccation tolerance and drought resistance, and the relative abundance of Chloroflexi was also high in this study – the bacteria in this phylum can produce energy through photosynthesis, reflecting their photoautotrophic properties [27]. In addition, Chloroflexi and Acidobacteria can promote nutrient uptake by the root system and maintain the balance between it and the microenvironment [30]. Du et al. [9] found that Acidobacteria and Proteobacteria were the dominant bacterial groups in the inter-root soils of nine plants in saline and alkaline soils in Daqing, when they studied the microbial community structure in inter-root soils. In this study, dominant microorganisms such as Proteobacteria, Actinobacteria, Acidobacteria, and Chloroflexi were present in the soils of the four treatments (Fig. 5A), and previous studies have shown that Actinobacteria and Acidobacteria as the dominant bacterial phyla are all salinophilic bacteria, and their occupancy in salinized soils is high. The soil in this experimental area is alkaline, which provides a good environment for the growth of Actinomycetes (Actinomycetes like a slightly alkaline environment) and makes them widely distributed in saline soil of Northwestern China [13, 59]. The relative abundances of Chloroflexi and Actinobacteria increased significantly with phytoremediation treatment in this study. After the modified treatment, pH in the saline soil tended to decrease, which was favorable for survival of Acidobacteria [19]. The dominant bacterial phyla in this study area were consistent with previous studies [14, 63].

There are more microbial species in bare soil at the phylum and genus levels 9 (Fig. 2B-C). This can be attributed to the presence of various special environmental conditions in these soils, resulting in the selection of microorganisms with extreme tolerance. After plant restoration, the soil environment is improved, which helps the growth of some new microorganisms, but may lead to the reduction of the original saline-tolerant microorganisms [17]. Therefore, although there are fewer microbial species in the repaired soil, there is usually a significant improvement in function and soil quality [31].

In addition, the LEFSe analysis showed that the three types of phytoremediation treatments caused significant changes in several groups of microorganisms (Fig. 4A) [7]: Gammaproteobacteria of the Proteobacteria was the most changed soil microorganism taxa under the CK; Betaproteobacteria was the most altered soil microorganism taxa under the Hg treatment – Betaproteobacteria is a class in the Proteobacteria and is most closely related to γ-Proteobacteria; under Sg treatment, Sphingomonas was the most altered soil microorganism taxa, and has high metabolic regulation mechanism and gene regulation ability, giving it great application potential in environmental remediation and plant growth promotion; and Actinobacteria was the most altered soil microorganism taxa under Hg||Sg treatment, is one of the most widespread bacteria in soil, and has attracted widespread attention for its excellent ability to degrade plant residues. The metagenomic analysis, which comprehensively reflects the variation in the functional constitution of soil microbes, showed that metabolism is the main function of soil microorganisms in saline soils to alleviate the adverse effects of salinized soils by comparing with the KEGG level A database (Fig. 3B).

Relationship between microorganisms and soil physicochemical properties

Soil microorganisms bridge the gap between soil and plant, and microorganisms can convert nutrients provided by soil into forms that can be absorbed and used by plants. At the same time, microorganisms can respond to the soil through altering community structure to maintain a dynamic balance of the soil microenvironment [36, 49, 50]. This study measured and analyzed the physicochemical properties of saline soils under four treatments and found that soil salinity decreased after saline plant phytoremediation treatment, with the decrease reaching 93.89% after Hg phytoremediation treatment, and soil OM content increased. The RDA and Pearson correlation analysis (Fig. 4B and C) showed that specific bacterial groups were associated with specific soil factors, and soil nutrient content was an important environmental factor influencing soil microbial community diversity and richness [63], while OM and pH also played significant roles in microbiota composition. For example, the dominant bacterial phylum Acidobacteria was positively correlated with soil pH and TK.

Many studies have also found that soil properties can directly affect microbial community structure [5, 24]. In this study, TN, OM, and Mg contents and pH were the main environmental factors causing differences in soil microbial community structure in saline soil of Northwestern China. There was clear separation of CK, Hg, Sg, and Hg||Sg treatments (Fig. 2H). Combined with the changes of network topological characteristics, the average clustering coefficients of Hg and Sg treatments significantly increased with the improvement, while the average clustering coefficients of Hg||Sg were basically unchanged, indicating that the aggregation and interaction of the microbial community co-occurrence network increased after Hg and Sg treatment. Thus, the Hg and Sg treatments not only improved the symbiotic relationship between microorganisms, but also made a closer connection between microbial communities.

Conclusion

This study comprehensively studied the structure composition and function of rhizosphere soil microbial communities of halophytes H. glomeratus and S. glauca in saline soil in Northwestern China by analyzing the biomass of halophytes, soil physicochemical properties, and metagenomic analysis (Graphical Abstract). The soil microbial community composition changed under different phytoremediation, and the ameliorative effect of various halophytes on saline soil showed that Hg, Sg, and Hg||Sg increased the diversity and function of the soil microbial community. Proteobacteria, Acidobacteria, Actinobacteria, and Chloroflexi were the main microorganisms utilized by the halophytes. In addition, microbial community changes were mainly caused by plant, TN, OM, Mg, and pH changes. The results showed that monocropping and intercropping could improve saline soil effectively, but interestingly, monocropping of H. glomeratus had the best effect on saline soil, as proven in the pot experiment. The results are important for guiding saline soil phytoremediationin the semi-arid area of Northwestern China using halophyte H.glomertaus.

Data availability

The datasets used and analysed during the current study are available from the corresponding author on reasonable request. Additionally, all relevant processed data, including tables, figures, and any other supplementary materials, are provided in the supplementary files accompanying this manuscript. For the plant materials used in this study, we have provided detailed information on the species, cultivar(s), and any specific growth conditions in the Methods section. Therefore, we encourage researchers to contact us for any further information or data requests that may not be immediately available in this manuscript or its supplementary materials.

Abbreviations

Abbreviation:

Full name

Hg:

Halogeton glomeratus monocropping

Sg:

Suaeda glauca monocropping

Hg||Sg:

Halogeton glomeratus and Suaeda glauca intercropping

TDS:

Total water-soluble salt

OM:

Soil organic matter

TN:

Total nitrogen

AN:

Alkaline dissolved nitrogen

TP:

Total phosphorus

AP:

Available phosphorus

TK:

Total potassium

AK:

Available potassium

RDA:

Redundancy analysis

NMDS:

Non-metric multidimensional scaling

NCBI:

National Center for Biotechnology Information

KEGG:

Non-redundant protein (Nr) database and the Kyoto Encyclopedia of Genes and Genomes

Kos:

KEGG orthologs

SE:

Standard error

References

  1. Ahmadi F, Mohammadkhani N, Servati M. Halophytes play important role in phytoremediation of salt-affected soils in the bed of Urmia Lake, Iran. Sci Rep. 2022;12:12223. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-022-16266-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Ajay S. (2015) Soil salinization and waterlogging: A threat to environment and agricultural sustainability. Ecological Indicators. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ecolind.2015.04.027

  3. Araya T, Mlahlwa AV, Elbasit MAMA, Newete SW. The impact of Tamarix invasion on the soil physicochemical properties. Sci Rep. 2022;12:5750. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41598-022-09797-3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Bai X, Zhang E, Wu J, Ma D, Zhang C, Zhang B, Liu Y, Zhang Z, Tian F, Zhao H, Wang B. Soil fungal community is more sensitive than bacterial community to modified materials application in saline–alkali land of Hetao Plain. Front Microbiol. 2024;15. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2024.1255536.

  5. Bissett A, Richardson AE, Baker G, Thrall PH. Long-term land use effects on soil microbial community structure and function. Appl Soil Ecol. 2011;51:66–78. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.apsoil.2011.08.010.

    Article  Google Scholar 

  6. Chang W, Chen W, Hu Y, Wang Z. Bacillus altitudinis LZP02 improves rice growth by reshaping the rhizosphere microbiome. Plant Soil. 2024;498:279–94. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-023-06435-3.

    Article  CAS  Google Scholar 

  7. Cheng H, Wu B, Wei M, Wang S, Rong X, Du D, Wang C. Changes in community structure and metabolic function of soil bacteria depending on the type restoration processing in the degraded alpine grassland ecosystems in Northern Tibet. Sci Total Environ. 2021;755:142619. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2020.142619.

    Article  CAS  PubMed  Google Scholar 

  8. Conforti M, Castrignanò A, Robustelli G, Scarciglia F, Stelluti M, Buttafuoco G. Laboratory-based Vis–NIR spectroscopy and partial least square regression with spatially correlated errors for predicting spatial variation of soil organic matter content. CATENA. 2015;124:60–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.catena.2014.09.004.

    Article  CAS  Google Scholar 

  9. Du YX, Xie BM, Cai HS, Tang L, Guo CH. Structural and functional diversity of rhizosphere microbial community of nine plant species in the Daqing saline-alkali soil region. Acta Ecol Sin. 2016;36:740–7. https://doiorg.publicaciones.saludcastillayleon.es/10.5846/stxb201404020621.

    Article  CAS  Google Scholar 

  10. Duan H, Ma Y, Liu R, Li Q, Yang Y, Song J. Effect of combined waterlogging and salinity stresses on euhalophyte Suaeda Glauca. Plant Physiol Biochem. 2018;127:231–7. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2018.03.030.

    Article  CAS  PubMed  Google Scholar 

  11. Feng Y, Xu X, Liu J, Han J, Lu H. Planting Suaeda salsa improved the soil properties and bacterial community diversity in a coastal mudflat. Land Degrad Dev. 2023;34:3262–71. https://doiorg.publicaciones.saludcastillayleon.es/10.1002/ldr.4681.

    Article  Google Scholar 

  12. Gerhardt KE, Gerwing PD, Greenberg BM. Opinion: taking phytoremediation from proven technology to accepted practice. Plant Sci. 2017;256:170–85. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plantsci.2016.11.016.

    Article  CAS  PubMed  Google Scholar 

  13. Gohel SD, Singh SP. Molecular phylogeny and diversity of the Salt-Tolerant Alkaliphilic Actinobacteria Inhabiting Coastal Gujarat, India. Geomicrobiol J. 2018;35:775–89. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/01490451.2018.1471107.

    Article  CAS  Google Scholar 

  14. Guan Y, Jiang N, Wu Y, Yang Z, Bello A, Yang W. Disentangling the role of salinity-sodicity in shaping soil microbiome along a natural saline-sodic gradient. Sci Total Environ. 2021;765:142738. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2020.142738.

    Article  CAS  PubMed  Google Scholar 

  15. Heng T, Liao RK, Wang ZH, Wu WY, Li WH, Zhang JZ. Effects of combined drip irrigation and sub-surface pipe drainage on water and salt transport of saline-alkali soil in Xinjiang, China. J ARID LAND. 2018;10:932–45. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s40333-018-0061-7.

    Article  Google Scholar 

  16. Hou J, Liu W, Wu L, Hu P, Ma T, Luo Y, Christie P. Modulation of the efficiency of trace metal phytoremediation by Sedum plumbizincicola by microbial community structure and function. Plant Soil. 2017;421:285–99. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-017-3466-8.

    Article  CAS  Google Scholar 

  17. Iannucci A, Canfora L, Nigro F, De Vita P, Beleggia R. Relationships between root morphology, root exudate compounds and rhizosphere microbial community in durum wheat. Appl Soil Ecol. 2021;158:103781. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.apsoil.2020.103781.

    Article  Google Scholar 

  18. Jiang S-Q, Yu Y-N, Gao R-W, Wang H, Zhang J, Li R, Long X-H, Shen Q-R, Chen W, Cai F. High-throughput absolute quantification sequencing reveals the effect of different fertilizer applications on bacterial community in a tomato cultivated coastal saline soil. Sci Total Environ. 2019;687:601–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2019.06.105.

    Article  CAS  PubMed  Google Scholar 

  19. Jiang Z, Zhang P, Wu Y, Wu X, Ni H, Lu Q, Zang S. Long-term surface composts application enhances saline-alkali soil carbon sequestration and increases bacterial community stability and complexity. Environ Res. 2024;240:117425. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.envres.2023.117425.

    Article  CAS  PubMed  Google Scholar 

  20. Jurado C, Díaz-Vivancos P, Gregorio B-E, Acosta-Motos JR, Hernández JA. Effect of halophyte-based management in physiological and biochemical responses of tomato plants under moderately saline greenhouse conditions. Plant Physiol Biochem. 2024;206:108228. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.plaphy.2023.108228.

    Article  CAS  PubMed  Google Scholar 

  21. Kearl J, McNary C, Lowman JS, Mei C, Aanderud ZT, Smith ST, West J, Colton E, Hamson M, Nielsen BL. Salt-tolerant halophyte rhizosphere Bacteria stimulate growth of Alfalfa in Salty Soil. Front Microbiol. 2019;10. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2019.01849.

  22. Khan A, Kong W, Khan S, Nawab J. Community succession and drivers of CO2-fixing microbes in recently deglaciated soils on the Tibetan Plateau. J Soils Sediments. 2023;23:1901–12. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11368-023-03446-6.

    Article  CAS  Google Scholar 

  23. Kumar A, Singh S, Mukherjee A, Rastogi RP, Verma JP. Salt-tolerant plant growth-promoting Bacillus pumilus strain JPVS11 to enhance plant growth attributes of rice and improve soil health under salinity stress. Microbiol Res. 2021;242:126616. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.micres.2020.126616.

    Article  CAS  PubMed  Google Scholar 

  24. Lammel DR, Nüsslein K, Tsai SM, Cerri CC. Land use, soil and litter chemistry drive bacterial community structures in samples of the rainforest and cerrado (Brazilian Savannah) biomes in Southern Amazonia. Eur J Soil Biol. 2015;66:32–9. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.ejsobi.2014.11.001.

    Article  CAS  Google Scholar 

  25. Lastiri-Hernández MA, Álvarez-Bernal D, Moncayo-Estrada R, Cruz-Cárdenas G, Silva García JT. Adoption of phytodesalination as a sustainable agricultural practice for improving the productivity of saline soils. Environ Dev Sustain. 2021;23:8798–814. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s10668-020-00995-5.

    Article  Google Scholar 

  26. Leewis MC, Kasanke C, Uhlik O, Leigh MB. Long-term legacy of phytoremediation on plant succession and soil microbial communities in petroleum-contaminated sub-arctic soils. SOIL. 2024;10:551–66. https://doiorg.publicaciones.saludcastillayleon.es/10.5194/soil-10-551-2024.

    Article  Google Scholar 

  27. Leite HMF, Calonego JC, Rosolem CA, Mendes LW, de Moraes LN, Grotto RMT, Araujo FF, Pereira APA, Melo VMM, Araujo ASF. Cover crops shape the soil bacterial community in a tropical soil under no-till. Appl Soil Ecol. 2021;168:104166. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.apsoil.2021.104166.

    Article  Google Scholar 

  28. Li B, Wang J, Yao L, Meng Y, Ma X, Si E, Ren P, Yang K, Shang X, Wang H. Halophyte Halogeton glomeratus, a promising candidate for phytoremediation of heavy metal-contaminated saline soils. Plant Soil. 2019;442:323–31. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-019-04152-4.

    Article  CAS  Google Scholar 

  29. Li F, Guo Y, Wang Z, Mu Y. Influence of different phytoremediation on soil microbial diversity and community composition in saline-alkaline land. Int J Phytoremediation. 2022;24:507–17. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/15226514.2021.1955240.

    Article  CAS  PubMed  Google Scholar 

  30. Li H, Xu Z, Yang S, Li X, Top EM, Wang R, Zhang Y, Cai J, Yao F, Han X. Responses of Soil Bacterial communities to Nitrogen Deposition and Precipitation Increment are closely linked with Aboveground Community Variation. Microb Ecol. 2016. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s00248-016-0730-z.

    Article  PubMed  Google Scholar 

  31. Li X, Chen D, Li B, Yang Y. Cd accumulation characteristics of Salvia tiliifolia and changes of rhizospheric soil enzyme activities and bacterial communities under a cd concentration gradient. Plant Soil. 2021;463:225–47. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-021-04905-0.

    Article  CAS  Google Scholar 

  32. Li Z, Deng X, Jin G, Mohmmed A, Arowolo AO. Tradeoffs between agricultural production and ecosystem services: a case study in Zhangye, Northwest China. Sci Total Environ. 2020;707:136032. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2019.136032.

    Article  CAS  PubMed  Google Scholar 

  33. Liang J, Shi W. Cotton/halophytes intercropping decreases salt accumulation and improves soil properties and crophysicochemical p productivity in saline-alkali soils under mulched drip irrigation: a three-year field experiment. Field Crops Res. 2021;262. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.fcr.2020.108027.

  34. Liao F, Qian J, Yang R, Gu W, Li R, Yang T, Fu X, Yuan B, Zhang Y. Metagenomics of gut microbiome for migratory seagulls in Kunming city revealed the potential public risk to human health. BMC Genomics. 2023;24:269. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12864-023-09379-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/su142416396.

    Article  CAS  Google Scholar 

  36. Liu Z, Fu B, Zheng X, Liu G. Plant biomass, soil water content and soil N:P ratio regulating soil microbial functional diversity in a temperate steppe: a regional scale study. Soil Biol Biochem. 2010;42:445–50. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.soilbio.2009.11.027.

    Article  CAS  Google Scholar 

  37. Ma Y, Rajkumar M, Oliveira RS, Zhang C, Freitas H. Potential of plant beneficial bacteria and arbuscular mycorrhizal fungi in phytoremediation of metal-contaminated saline soils. J Hazard Mater. 2019;379:120813. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jhazmat.2019.120813.

    Article  CAS  PubMed  Google Scholar 

  38. Manousaki E, Kalogerakis N. Halophytes—An Emerging Trend in Phytoremediation. Int J Phytoremediation. 2011;13:959–69. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/15226514.2010.532241.

    Article  CAS  PubMed  Google Scholar 

  39. Molefe RR, Amoo AE, Babalola OO. Metagenomic insights into the bacterial community structure and functional potentials in the rhizosphere soil of maize plants. J Plant Interact. 2021;16:258–69. https://doiorg.publicaciones.saludcastillayleon.es/10.1080/17429145.2021.1936228.

    Article  CAS  Google Scholar 

  40. Nouri H, Chavoshi Borujeni S, Nirola R, Hassanli A, Beecham S, Alaghmand S, Saint C, Mulcahy D. Application of green remediation on soil salinity treatment: a review on halophytoremediation. Process Saf Environ Prot. 2017;107:94–107. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.psep.2017.01.021.

    Article  CAS  Google Scholar 

  41. Paungfoo-Lonhienne C, Yeoh YK, Kasinadhuni NRP, Lonhienne TGA, Robinson N, Hugenholtz P, Ragan MA, Schmidt S. Nitrogen fertilizer dose alters fungal communities in sugarcane soil and rhizosphere. Sci Rep. 2015;5:8678. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/srep08678.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Rat A, Naranjo HD, Krigas N, Grigoriadou K, Maloupa E, Alonso AV, Schneider C, Papageorgiou VP, Assimopoulou AN, Tsafantakis N, Fokialakis N, Willems A. Endophytic Bacteria from the roots of the Medicinal Plant Alkanna tinctoria Tausch (Boraginaceae): exploration of Plant Growth promoting properties and potential role in the production of Plant secondary metabolites. Front Microbiol. 2021;12. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2021.633488.

  43. Rengasamy P. Soil salinization. Oxford University Press; 2016.

  44. Roe RAL, MacFarlane GR. The potential of saltmarsh halophytes for phytoremediation of metals and persistent organic pollutants: an Australian perspective. Mar Pollut Bull. 2022;180:113811. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.marpolbul.2022.113811.

    Article  CAS  PubMed  Google Scholar 

  45. Salwan R, Sharma A, Sharma V. Microbes mediated plant stress tolerance in saline agricultural ecosystem. Plant Soil. 2019;442:1–22. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-019-04202-x.

    Article  CAS  Google Scholar 

  46. Santonja M, Rancon A, Fromin N, Baldy V, Hättenschwiler S, Fernandez C, Montès N, Mirleau P. Plant litter diversity increases microbial abundance, fungal diversity, and carbon and nitrogen cycling in a Mediterranean shrubland. Soil Biol Biochem. 2017;111:124–34. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.soilbio.2017.04.006.

    Article  CAS  Google Scholar 

  47. Schmid MW, van Moorsel SJ, Hahl T, De Luca E, De Deyn GB, Wagg C, Niklaus PA, Schmid B. Effects of plant community history, soil legacy and plant diversity on soil microbial communities. J Ecol. 2021;109:3007–23. https://doiorg.publicaciones.saludcastillayleon.es/10.1111/1365-2745.13714.

    Article  Google Scholar 

  48. Shan S, Wei Z, Cheng W, Du D, Zheng D, Ma G. Biofertilizer based on halotolerant microorganisms promotes the growth of rice plants and alleviates the effects of saline stress. Front Microbiol. 2023;14. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2023.1165631.

  49. Shi W, Li M, Wei G, Tian R, Li C, Wang B, Lin R, Shi C, Chi X, Zhou B, Gao Z. The occurrence of potato common scab correlates with the community composition and function of the geocaulosphere soil microbiome. Microbiome. 2019;7:14. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40168-019-0629-2.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Steenwerth KL, Drenovsky RE, Lambert JJ, Kluepfel DA, Scow KM, Smart DR. Soil morphology, depth and grapevine root frequency influence microbial communities in a Pinot noir vineyard. Soil Biol Biochem. 2008;40:1330–40. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.soilbio.2007.04.031.

    Article  CAS  Google Scholar 

  51. Tang L, Zhan L, Han Y, Wang Z, Dong L, Zhang Z. Microbial community assembly and functional profiles along the soil-root continuum of salt-tolerant Suaeda glauca and Suaeda salsa. Front Plant Sci. 2023;14. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fpls.2023.1301117.

  52. Visagie CM, Magistà D, Ferrara M, Balocchi F, Duong TA, Eichmeier A, Gramaje D, Aylward J, Baker SE, Barnes I, Calhoun S, De Angelis M, Frisvad JC, Hakalova E, Hayes RD, Houbraken J, Grigoriev IV, LaButti K, Leal C, Lipzen A, Ng V, Pangilinan J, Pecenka J, Perrone G, Piso A, Savage E, Spetik M, Wingfield MJ, Zhang Y, Wingfield BD. IMA genome-F18. IMA Fungus. 2023;14:21. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s43008-023-00121-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Wang B, Peng Q, Wang R, Yu S, Li Q, Huang C. Efficient Microcystis removal and sulfonamide-resistance gene propagation mitigation by constructed wetlands and functional genes analysis. Chemosphere. 2022;292:133481. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.chemosphere.2021.133481.

    Article  CAS  PubMed  Google Scholar 

  54. Wang H, Zhang M, Qu C, Fei Y, Liang J, Bai W, Zhao W, Xiao G, Liu G. Characterization and correlation of Microbiota and higher alcohols based on metagenomic and metabolite profiling during Rice-Flavor Baijiu Fermentation. Foods. 2023;12:2720. https://doiorg.publicaciones.saludcastillayleon.es/10.3390/foods12142720.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Wang L, Zhang X, Xu S. Is salinity the main ecological factor that influences foliar nutrient resorption of desert plants in a hyper-arid environment? BMC Plant Biol. 2020a;20:461. https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-020-02680-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Wang Y, Liu L, Luo Y, Awasthi MK, Yang J, Duan Y, Li H, Zhao Z. Mulching practices alter the bacterial-fungal community and network in favor of soil quality in a semiarid orchard system. Sci Total Environ. 2020b;725:138527. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2020.138527.

    Article  CAS  PubMed  Google Scholar 

  57. Wang Y, Liu L, Yang J, Duan Y, Luo Y, Taherzadeh MJ, Li Y, Li H, Awasthi MK, Zhao Z. The diversity of microbial community and function varied in response to different agricultural residues composting. Sci Total Environ. 2020c;715:136983. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2020.136983.

    Article  CAS  PubMed  Google Scholar 

  58. Yamamoto K, Shiwa Y, Ishige T, Sakamoto H, Tanaka K, Uchino M, Tanaka N, Oguri S, Saitoh H, Tsushima S. Bacterial diversity Associated with the Rhizosphere and Endosphere of two Halophytes: Glaux maritima and Salicornia europaea. Front Microbiol. 2018;9. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fmicb.2018.02878.

  59. Yandi S, Tao M, Huakun Z, Honglin L, Zhonghua Z, Li M, Ruimin Q, Hongye S, Tao C, Jingjing W, Xue H. Characterization of the plant–soil feedback index in alpine meadow degradation and recovery: a field experiment. Front Environ Sci. 2023;10. https://doiorg.publicaciones.saludcastillayleon.es/10.3389/fenvs.2022.1097030.

  60. Yin F, Zhang F. Reclamation of abandoned saline-alkali soil increased soil microbial diversity and degradation potential. Plant Soil. 2022;477:521–38. https://doiorg.publicaciones.saludcastillayleon.es/10.1007/s11104-022-05451-z.

    Article  CAS  Google Scholar 

  61. Yue Y, Shao T, Long X, He T, Gao X, Zhou Z, Liu Z, Rengel Z. Microbiome structure and function in rhizosphere of Jerusalem artichoke grown in saline land. Sci Total Environ. 2020;724:138259. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.scitotenv.2020.138259.

    Article  CAS  PubMed  Google Scholar 

  62. Zhang B, Wang N. Study on the harm of saline Alkali Land and its Improvement Technology in China. IOP Conf Series: Earth Environ Sci. 2021;692:042053. https://doiorg.publicaciones.saludcastillayleon.es/10.1088/1755-1315/692/4/042053.

    Article  Google Scholar 

  63. Zhang C, Zhou X, Wang X, Ge J, Cai B. Elaeagnus angustifolia can improve salt-alkali soil and the health level of emphasizingsizing the driving role of core microbial communities. J Environ Manage. 2022;305:114401. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jenvman.2021.114401.

    Article  CAS  PubMed  Google Scholar 

  64. Zhang G, Bai J, Zhai Y, Jia J, Zhao Q, Wang W, Hu X. Microbial diversity and functions in saline soils: a review from a biogeochemical perspective. J Adv Res. 2023. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.jare.2023.06.015.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Zheng Q, Hu Y, Zhang S, Noll L, Wanek W. Soil multifunctionality is affected by the soil environment and by microbial community composition and diversity. Soil Biol Biochem. 2019;136:107521. https://doiorg.publicaciones.saludcastillayleon.es/10.1016/j.soilbio.2019.107521.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank the professionals of BioMed Proofreading LLC for English corrections and copyediting this manuscript.

Funding

This work was supported by National Natural Science Foundation of China (32001514, 31960072); Industrial Support Project of Colleges and Universities in Gansu Province (2021CYZC-12); Key Projects of Natural Science Foundation of Gansu Province (20JR10RA507; 22JR5RA880); China Agriculture Research System (Grant CARS-05-04B-2); Fuxi Talent Project of Gansu Agricultural University (Ganfx-03Y06; GAUfx-04Y011).

Author information

Authors and Affiliations

Authors

Contributions

JW and MS conceived the study and designed the experimental approach. They also supervised data collection and contributed significantly to the interpretation of the results.YL, PL, ES and BL performed the experiments, and was responsible for the initial data analysis. They also assisted in writing the Methods section of the manuscript.HW and XS contributed to the literature review, helping to frame the research question and position the study within the broader context of the field. They also revised the Introduction and Discussion sections.YM, XM, KY and HZ prepared all figures and tables, ensuring they accurately reflect the data and comply with journal guidelines. They also provided valuable input on the overall structure and flow of the manuscript.All authors were involved in critically reviewing the manuscript and providing feedback throughout the writing process. The final version of the manuscript was approved by all authors.

Corresponding author

Correspondence to Huajun Wang.

Ethics declarations

Ethics approval and consent to participate

The present study was conducted in accordance with the ethical principles outlined in the relevant guidelines and regulations. The study objectives, methodology, potential risks and benefits, and data protection measures, and determined that the study complied with all ethical standards and did not pose undue harm to any participants or the environment. For Plant Collection and Studies: The collection of plant samples was conducted in accordance with local laws and regulations, and where necessary, permission was obtained from the relevant authorities or land owners. Special care was taken to minimize any potential impact on the environment and to ensure the sustainability of the plant populations studied. Specifically, for the collection of soil samples and plant material, we obtained written permission from the private landowner(s) on whose property the sampling took place. This permission was granted after explaining the purpose of the study and ensuring that the landowner(s) understood the nature and scope of the sampling activities. We also ensured that our sampling activities did not cause any harm to the environment or the landowner’s property. In summary, we adhered to the highest ethical standards in conducting this research, and obtained all necessary permissions and consents prior to the commencement of fieldwork.

Consent for publication

We, the undersigned authors of this manuscript titled " Metagenomic analysis reveal the phytoremediation effects of monocropping and intercropping of halophytes Halogeton glomeratus and Suaeda glauca in saline soil of Northwestern China”, hereby give our full consent for the publication of this work in BMC Plant Biology. We confirm that: (1) We have all reviewed the manuscript and its contents, and we agree that it represents our collaborative research efforts accurately and fairly. (2) Each author has made a significant intellectual contribution to the conception, design, execution, or interpretation of the reported study, and is fully accountable for the accuracy and integrity of the work. (3) We have no ethical, legal, or financial conflicts of interest that would prevent the publication of this manuscript. (4) We acknowledge that we have been offered the opportunity to revise the manuscript should the editors or reviewers suggest any changes to improve the quality or clarity of the work. (5) We grant BMC Plant Biology the exclusive right to publish this manuscript in all forms and media, including electronic and print formats, for the duration of copyright. (6) We understand that the manuscript, if accepted for publication, will be made freely accessible online under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. (7) We agree to be identified as authors of this work and to have our names appear in the published article as listed below.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, J., Song, M., Yao, L. et al. Metagenomic analysis reveal the phytoremediation effects of monocropping and intercropping of halophytes Halogeton glomeratus and Suaeda glauca in saline soil of Northwestern China. BMC Plant Biol 25, 213 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-025-06225-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12870-025-06225-2

Keywords