Skip to main content

Effects of microbial organic fertilizers on Astragalus membranaceus growth and rhizosphere microbial community



The application of excessive chemical fertilizers during the cultivation of Astragalus membranaceus leads to a decline in the quality of this medicinal plant as well as the soil’s sustainable productivity. In this study, we developed a special microbial organic fertilizer for A. membranaceus and investigated its effects on plant growth and rhizosphere microbial communities.


The root biomass and main active components of A. membranaceus in different growth stages were measured to assess the impacts of microbial organic manure on plant growth. Meanwhile, 16S rRNA and ITS1 amplicons were amplified and high-throughput sequencing was performed to detect the dynamic impacts of microbial organic manure on rhizosphere microbial communities.


The results demonstrated that microbial organic manure significantly increased wet and dry weights of A. membranaceus seedlings and the accumulation of two effective components (flavonoids and saponin) in bacterial fertilizer treatment groups are significant higher than the control group. Research on rhizosphere microbial flora shows that the number and polymorphism of bacteria and fungi were decreased after the application of special fertilizer during the rapid growth period of plant and then gradually increased with seedling growth. The community structure of bacteria was regulated after the application of special fertilizer, and the beneficial bacteria for plant growth are enriched. Functional profiles prediction showed that significant shifts in metabolic functions impacting KEGG pathways of the microbial fertilizer treatment groups are related to metabolism and biosynthesis.


The results indicate that the microbial organic manure can improve A. membranaceus growth by providing appropriate nutrients and regulating the rhizosphere microbial community which has good potential in ecological cultivation of A. membranaceus.


Soil microorganisms contribute to nutrient recycling, system stability, and anti-jamming capabilities and sustainable productivity of the soil’s ecosystem (Heijden et al. 2010; Richard et al. 2017; Singh 2015). The application of natural microbial fertilizer, which is produced by plant-growth-promoting bacteria (PGPB), is a potential alternative for a more sustainable and environmentally friendly agricultural development (Ahemad and Kibret 2014; Ferreira et al. 2019). Many previous studies have shown the positive effects of different microbial fertilizers on major crops, which have been proven to effectively increase the nutrient utilization of crops, thereby increasing plant production in a cost-effective and environmentally friendly way (Tabassum et al. 2017). They can also improve soil porosity, structural stability, nutrient availability, and biological activity (Francis et al. 2010; Wang et al. 2011).

Astragalus membranaceus, a member of the Fabaceae family, is an important herb that has been used in traditional Chinese medicine (Liu et al. 2011). Astragalus membranaceus has significant effects in improving the immune system and cardiovascular system of the body, and it also has significant therapeutic effects in anti-aging, anti-virus, and anti-tumor (James and Ko 2007; Yang et al. 2013; Zhou et al. 2018). With over-harvesting of the wild A. membranaceus and high demand in the market, the available wild resources are becoming scarce. Thus, the cultivation of A. membranaceus continues to expand and the amount of chemical fertilizers and pesticides is large during the planning process. Unfortunately, these chemicals impart various negative effects on the agricultural ecosystem, including soil compaction and degradation, reduction of soil microbial diversity, ground-water resources pollution, and atmospheric contamination (Savci 2012; Sharma and Paliyal 2014; Yang et al. 2012).

As a type of Leguminosae plant, the roots of A. membranaceus usually interact with rhizobia in the soil, resulting in root nodule symbiosis which provides nutrients for plant growth and promotes the accumulation of active ingredients. Azotobacter fertilizer, which is mainly composed of native screened rhizobia, is an environmentally friendly and long-acting fertilizer for A. membranaceus. In the prophase experiment, we isolated two strains of rhizobia from the root nodule of A. membranaceus, Rhizobium sp. t16, and Sinorhizobium sp. t21, which demonstrate high nitrogen fixation activity (Zhiquan et al. 2016). Laboratory pot-experiment results showed that the two strains could significantly promote the growth of A. membranaceus seedlings (Zhiquan et al. 2016). However, future studies are necessary to determine the effects of organic microbial fertilizer in the field. The selection of inoculation method and carrier material of PGPB should be according to the viability of microorganisms, availability of plant, and feasible application. Solid carrier materials composed of nutrients and organic matter are considered the most suitable materials for field applications. (Tabassum et al. 2017). Meanwhile, returning stems and leaves to the field can affect soil structure and gas exchange and can also improve soil fertility. However, the effects of organic microbial fertilizer with stems and leaves powder as the carrier material for field plant growth and rhizosphere microbial communities remain unclear.

Therefore, we performed a field experiment to access the effects of special manure composed of stems, leaves powder, and two native strains of rhizobia for field plant growth and microbial community diversity. The purpose of this research was (i) to examine the effects of microbial organic fertilizers on A. membranaceus growth and accumulation of the main medicinal components (ii) to assess the effects of microbial organic fertilizers on rhizosphere microbial community. We hypothesized that microbial organic fertilizers could promote A. membranaceus growth and affect the function and diversity of rhizosphere microbial communities.

Materials and methods

Site description

The field experimental site was located in Jingle County, Shanxi Province of China (38° 36′ 51.47″ N, 111° 94′ 54.32″ E). The farm area is located in a typical monsoon climate at an altitude of 2032 m, with a mean annual rainfall of 447.6 mm, humidity of 61%, and temperature of 7.3 °C. The soil basic physicochemical properties of the experimental site are as follows: total organic matter 4.57 g/kg, pH 8.88, total N 0.42 g/kg, total P 0.47 g/kg, total K 21.81 g/kg, available N 27.63 mg/kg, available P 7.21 mg/kg, and available K 81.24 mg/kg.

Preparation of the microbial organic manure

Two strains of nitrogen-fixing bacteria, Rhizobium sp. t16 and Sinorhizobium sp. t21, were isolated from the root nodule of A. membranaceus in Hunyuan County, Shanxi Province, and deposited in China General Microbiological Culture Collection Center, with the access number of CGMCC No. 9676 and 9677 (Zhiquan et al. 2016). Before the experiments, the bacteria were inoculated into Luria-Bertani liquid medium for fermentation at the temperature of 28 °C, and the concentration of bacteria was regulated at 108 bacteria per mL.

The 3-year-old plants of A. membranaceus in the Jingle planning base were selected for the stems and leaves collection. The stems and leaves were harvested in May and August with 100 g for each plant. The plant tissue were washed with distilled water and dried at 70 °C and then crushed into powder via 60 mesh sieves for subsequent steps. Bacteria liquid was added at a powder to liquid ratio of 7:3. Different types of special manure containing the wet powder and bacteria solution were prepared for use.

Experimental design

The field experiment was conducted at May 2018 and three treatments using different combinations of bacterial solution and stem and leaf powder were performed based on randomized block design. The treatments were designed as follows: no fertilizer (control group), strain t16 bacterial solution + powder (group 4), strain t21 bacterial solution + powder (group 7), and t16 + t21 bacterial solution + powder (group 9). Five hundred milligrams of microbial fertilizer was applied to each group of experimental plot before sowing (1 m × 3 m) and each treatment was conducted with three replications. The soil in the control group was only plowed before sowing. The row spacing and plant distance for sowing is 25 cm and 15 cm in each plot when sowing. This study includes 4 groups of experimental plots according to the experimental designs (CK, group 4, group 7, and group 9).


Collected samples included rhizosphere soil (soil directly associated with the roots), soil for physicochemical parameters and roots. Soil samples before sowing were collected at 20 cm deep below the surface in triplicate (CK0). Rhizosphere soil samples were collected in triplicate at 30 days after sowing (CK1, 4.1, 7.1, 9.1) and 90 days after sowing (CK2, 4.2, 7.2, 9.2). The definitions of different samples for subsequent analysis were as follows: the soil control sample, CK0, was collected at day 0 and corresponded to the sowing date of A. membranaceus. The soil samples collected at 30 days (CK1, 4.1, 7.1, 9.1), corresponding to the seedling stage (also the rapid growth period), including the germination stage (0-15 days). The soil samples collected at 90 days (CK2, 4.2, 7.2, 9.2) correspond to the stable growth period and the plant flowering period.

Rhizosphere soil samples were gathered by brushing off the remaining soil adhering to the roots with sterile brushes, mixed together to reduce errors, placed in 10 ml sterile centrifuge tubes and then immediately stored at −80 °C for DNA extraction. Five seedlings of A. membranaceus were randomly collected from each group (30 days, 60 days, 90 days, and 120 days), placed in clean plastic, stored at icebox and taken back to the laboratory for growth parameters analysis (fresh weight, dry weight, root diameter, and effective components). Roots were cleaned with distilled water, absorbed moisture, measured the fresh weight, root diameter, dried at 70 °C to constant quality and measured the dry weight.

The dried roots of Astragalus (120 days) were crushed and passed through a 0.25-mm sieve, and the main medicinal components of Astragalus were extracted by flash extraction method. With rutin as the control standard curve, the total flavonoids content was determined by ultraviolet spectrophotometry. A standard curve was prepared with astragaloside IV as a control, and the total saponin content was determined by vanillin-concentrated sulfuric acid colorimetry. Anhydrous glucose was used as a control standard curve, and the polysaccharide content was determined by the phenol-sulfuric acid colorimetric method (Chun-Hong et al. 2011; Xi et al. 2010; Zheng et al. 2011). The measured values were statistically analyzed using ANOVA (p < 0.05) with SPSS to detect the significant difference between groups.

DNA extraction

DNA of soil samples collected at day 0, day 30, and day 90 was extracted with an Omega E.Z.N.A soil DNA kit purchased from Feiyang BIOTECH Co., Ltd. (Guangzhou, China) according to the manufacturer’s instructions (Dineen et al. 2010). The DNA was gathered with 0.8% agarose gel electrophoresis and then stored at −80 °C for further analysis.

Illumina MiSeq sequencing

The V4 region of 16S rRNA amplicons for bacterial community analysis were amplified by universal primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′). For fungi amplicons, we amplified the internal transcribed spacer 1 (ITS 1) region with universal primers 1743F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and 2043R (5′-GCTGCGTTCTTCATCGATGC-3′). Amplicons were generated using ExTaq HS (TaKaRa Bio Inc., Shiga, Japan) by polymerase chain reaction (PCR) under the following conditions: 95 °C for 2 min, followed by 35 cycles at 95°C for 30 s, annealing at 55 °C for 1 min, extension at 72 °C for 1 min, and a final extension at 72 °C for 10 min. All amplicons of bacteria and fungi were then sequenced by an Illumina MiSeq platform and MiSeq Reagent Kit v1 (Illumina, Inc., Santiago, CA, USA) at the Beijing Genomics Institute (Shenzhen, China).

Raw data generating from Illumina MiSeq platform were cleaned using python scripts as follows: (1) confirmation of barcode and adaptor completeness; (2) delete sequences which containing more than one ambiguous base (N); and (3) removal of sequences shorter than 100 bp. Cleaned pair-end reads were then connected into tags using the COPE software (Connecting Overlapped Pair-End, V 1.2.1) (Liu et al. 2012). Sequencing noise was removed using the pre.cluster tool in the MOTHUR software package v. 1.35.1 (Schloss et al. 2009). Chimera detection and removal were conformed using UCHIME v. 4.2 (Edgar et al. 2011). MOTHUR was also used to determine the Operational Taxonomy Unit (OTU) with a sequence identity threshold of 97% (Schloss et al. 2009). OTUs annotation was performed by local NCBI-BLAST v.2.9.0 and Usearch v.11.0 (Edgar 2018) and 16S rRNA training set 9 in the RDP database and ITS1 sequences in UNITE database were used as annotation references.

Statistical analysis

OTUs sequences and reads data were used for Alpha diversity analysis to assess the richness and diversity within the microbial communities. The Alpha diversity index was calculated by the MOTHUR package, including Simpson diversity, Shannon diversity, Chao 1 indices, as well as observed species and phylogenetic diversity (Kemp and Aller 2004). Principal coordinate analysis (PCoA) was performed with QIIME (Caporaso et al. 2010) to measure dissimilarity at phylogenetic distances based on UniFrac analysis, and the results were visualized using KING (Crawford et al. 2009). Statistically significant differences between experimental and control groups were determined using the non-parametric Mann-Whitney U test and multivariate analysis of variance in R v.3.6.0 platform. Metagenome functional profiling of bacterial community was predicted with 16S rRNA reads clustered at a 97% identity threshold using PICRUSt (Langille et al. 2013) and the reference Greengenes (version 13.5) database (DeSantis et al. 2006). Significantly enriched functional categories in different treatment groups were accessed using LEfSE (Segata et al. 2011) with a p value < 0.05 and a linear discriminant analysis (LDA) score >2.


Effects of microbial organic fertilizer on the growth of A.membranaceus

The seedlings of A. membranaceus were harvested at 30 days, 60 days, 90 days, and 120 days after sowing. The fresh weight, taproot diameter, and dry weight of the roots were measured (Fig. 1a-c). The results showed no difference between the samples obtained 30 days after sowing. On the 60th and 90th days, there were higher values in the experimental group compared to the control group, but there was no statistical difference. On the 120th day, the fresh weight and dry weight of groups 4 and 9 were statistically higher than that of group CK and 7 (p < 0.05) (Fig. 1a). The results suggest that the first 30 days were necessary for the seeds to germinate into the main root; the root was short with low weight, and there were no obvious interactions with soil microorganisms. Lateral roots had grown by the 60th or 90th day of growth but were not developed enough to make a difference between the groups. On the 120th day, the number of lateral roots was large enough to support plant growth, indicating that the application of bacterial fertilizer to plants with well-developed roots might be more beneficial than application to other plants.

Fig. 1

The effect of different fertilizers on the growth of Astragalus. Note: 1a, fresh weight of Astragalus; 1b, dry weight of Astragalus; 1c, root diameter of Astragalus. CK, without microbial fertilizer; 4, strain t16 + powder; 7, strain 21 + powder; 9, t16 + t21 + powder. Different lower case letters refer to significant differences between groups based on one-way ANOVA (p < 0.05)

Determination of main medicinal components

The main medicinal components of Astragalus roots (sampled at 120 days after seedling) were measured, including the total flavonoids content, the total saponin content, and the polysaccharide content (Fig. 2). The results showed that the accumulation of flavonoids and saponin in bacterial fertilizer treatment group 7 and group 9 are higher than in the group CK, but the polysaccharides showed no difference between groups. Interestingly, the accumulation of flavonoids and saponin in treatment group 7 (p < 0.05) and group9 (p < 0.01) is both significantly higher than in group CK. However treatment group 4 showed no significant difference compared with the group CK (Fig. 2a-b).

Fig. 2

The effect of different fertilizers on the main medicinal components (flavonoids, saponin, and polysaccharide) of A.membranaceus at 120 days after sowing. LSD test was used to detect the significant difference between groups and at a p value <0.05 marked as * and a p value <0.01 marked as **. Note: CK, without microbial fertilizer; group 4, strain t16 + powder; group 7, strain 21 + powder; group 9, strain t16 + strain t21 + powder.

Abundance and diversity of the bacterial and fungal microbiota

Illumina MiSeq sequencing was utilized to measure the microbial abundance and community diversity during the growth of A. membranaceus. A total of 103,396 high-quality 16S rRNA reads were generated from 9 sample sets for bacteria, with an average of 11,488 reads per sample (range: 10,332-12,214). Meanwhile, 370,055 high-quality ITS reads were obtained for fungi throughout the sequencing process, with an average of 41,117 reads per sample (range: 40,121-42,286) (Table 1). High-quality pair-end reads were then connected to tags and clustered into 982 OTUs for bacteria and 564 OTUs for fungi (both at 97% similarity). Finally, we obtained an average of 925 OTUs for bacteria (range: 880-949) and 334 for fungi (range: 311-365) after removing singletons. These sequences were subjected to further analysis.

Table 1 Sample information, microbial diversity and sequence abundance

Rarefaction curves of observed species and Simpson diversity indices were plotted for bacteria and fungi communities respectively (Fig. 3a-d). Although the rarefaction curve was not parallel to the x-axis, the Simpson diversity index reached saturation, indicating that most of the microbial diversity was captured in samples during the growth of A. membranaceus. Shannon index showed that bacterial diversity in treatment group 9 was higher than that in other groups both in the rapid growth period (6.07) and stable growth period (6.2) of A. membranaceus, and treatment group 7 had the lowest bacterial diversity among all the experimental groups. However, diversity indices between experimental and control groups showed no significant differences in bacterial and fungal diversity (p > 0.05) (Table 1).

Fig. 3

Rarefaction curves and Shannon index curves of different groups. Note: a, rarefaction curves of bacteria; b, rarefaction curves of fungi; c, Shannon curves of bacteria; d, Shannon curves of fungi

Characterization of bacterial and fungal communities at phylum and genus levels

To examine the changes in the bacterial communities alongside plant growth, the relative abundance of bacterial taxa was compared at both phylum and genus levels between the experimental and control groups. There were eight bacterial phyla identified in the special manure group, including Proteobacteria, Acidobacteria, Actinobacteria, Gemmatimonadetes, Chloroflexi, Nitrospirae, Bacteroidetes, and Planctomycetes. Relative abundance analysis of bacterial phyla throughout the entire A. membranaceus growth process revealed that Proteobacteria, Acidobacteria, and Actinobacteria were the predominant bacterial phyla in both experimental and control groups (Fig. 4a). The proportion of Proteobacteria remained prominent throughout the growth period at an average level of 34% once it reached dominance. Chloroflexi increased during plant growth, from the relative abundance of 3.16% before sowing (CK0) to 9.77% at the stable growth period (CK2). The abundance of Actinobacteria decreased from 28.52% before sowing (CK0) to 9.82% by day 30 (CK1) and then increased to 21.53% at day 90 (CK2). The observed decrease in the stable growth period of Actinobacteria might be due to interspecific competition with Proteobacteria and Chloroflexi, which experienced a significant increase in the same period. The inhibitory effect on Actinobacteria is more obvious in the treatment groups, especially in group 7.

Fig. 4

Relative abundance of bacteria at phyla (a) and genus (b) levels of different groups. X-axis shows the samples collected before sowing (CK0), 30 days after sowing (CK1, 4.1, 7.1, 9.1), and 90 days after sowing (CK2, 4.2, 7.2, 9.2)

At the genus level, uncultured bacteria were the dominant genera in all groups, ranged from 20.46 to 26.79% (Fig. 4b). During growth, the abundance of Arthrobacter decreased from 9.31% before sowing (CK0) to 2.33% at the stable growth period (CK2). In the treatment group 7, the abundance of Arthrobacter even dropped to 0.47% at the stable growth period (7.2). In contrast, the abundance of Nitrospira has a tendency to expand during the growth of A. membranaceus.

Both bacteria strain Rhizobium sp.t16 and Sinorhizobium sp.t21 were approximately estimated in the sequencing results. However, the reads which are identical with the two strains were low in all of the day 30 and day 90 samples, ranging from 102 reads/10700 reads to 508 reads/11723 reads in treatment groups and 27 reads/12174 reads to 94 reads/12174 reads for control group. But, the reads number of the two strains were still more in treatment groups than in control group.

Compared with bacteria, the fungal microbiota in A. membranaceus growth was less diverse. Only five fungal phyla were detected in experimental groups, compared with three fungal phyla in control group. Ascomycota and Zygomycota were the dominant phyla which widely existed in all groups, accounting for 52.99% and 41.17% of the average relative abundance of fungi at the phylum level (Fig. 5a). Ascomycota and Zygomycota contributed 34.62% and 62.83% in the rapid growth period, and 71.36% and 19.51% in the later stage, respectively.

Fig. 5

Relative abundance of fungi at phyla (a) and genus (b) levels in different groups. X axis shows the samples collected before sowing (CK0), 30 days after sowing (CK1, 4.1, 7.1, 9.1), and 90 days after sowing (CK2, 4.2, 7.2, 9.2)

At the genus level, Mortierella dominated the early period in experimental groups, accounting from 32.29 to 77.09%, and then decreased sharply at the stable growth period of plant (Fig. 5b). The abundance of Gibberella in experimental group 4.2 and 9.2 is significant less than in control group CK2. Paraphoma decreased sharply from 21.38% in group CK0 to an average of 1.63% in other groups with the growth of plant, which might be the result of inhibition of environmental fungi by root exudates (Badri and Vivanco 2009). Chaetomium was the second most abundant genus in the growth period, during which abundance increased from 3.75% on day 30 to 14.68% on day 90. The abundance of Preussia in experimental groups 4.2, 7.2, and 9.2 are all higher than that in the control group CK2 (1.36%), especially in group 4.2 (13.11%).

Effects of organic microbial fertilizer on the regulation of bacterial and fungal community structure

To understand the bacterial community structure, a heatmap was generated to reveal differences in genera between experimental and control groups (Fig. 6). At the rapid growth period of A. membranaceus, the genera structure in treatment group 7.1 and control group CK1 has the most significant difference, while in the stable growth period it is between the group 9.2 and the control group CK2. The dominant strains are different in each group based on the plant growth periods, and there are also differences in each treatment groups. Meanwhile, PCoA and cluster analysis of the bacterial microbiota were performed to evaluate the similarities between bacterial communities in all groups. Both weighted (PC1 variance = 37.12%, PC2 variance = 18.07%; Fig. 7a) and unweighted (PC1 variance = 50.89%, PC2 variance = 17.37%; Fig. 7b) PCoA was performed. The PCoA result based on weighted UniFrac values revealed an overlap of bacterial communities between rapid and stable growth periods of A. membranaceus, but cluster analysis indicated a clear difference between the two periods based on the weighted UniFrac metric (Fig. 8). The result indicated that the bacterial communities in treatment group 7.1 and control group CK1 had the most significant difference in rapid growth period of A. membranaceus. In the stable growth period, it was treatment group 9.2 and control group CK2 which had the significant difference in bacterial communities.

Fig. 6

Heatmap of abundant bacterial genera present in the microbial community of different groups. Note: Group A, samples of 0 day (CK0); group B, samples of 30 day (CK1, 4.1, 7.1, 9.1); group C, samples of 90 day (CK2, 4.2, 7.2, 9.2)

Fig. 7

Principal coordinate analysis of bacterial communities based on (a) weighted and (b) unweighted UniFrac metrics of samples in different groups. Note: Group A, samples of 0 day (CK0); group B, samples of 30 day (CK1, 4.1, 7.1, 9.1); group C, samples of 90 day (CK2, 4.2, 7.2, 9.2)

Fig. 8

Cluster analysis of the bacterial microbiota in different groups based on weighted UniFrac distances. Note: Group A, samples of 0 day (CK0); group B, samples of 30 day (CK1, 4.1, 7.1, 9.1); group C, samples of 90 day (CK2, 4.2, 7.2, 9.2)

Multivariate analysis was performed to estimate the similarities of fungal microbiota communities during the different growth period (PC1 variance = 67.56%, PC2 variance = 10.97%; Fig. 9a). PCoA result of fungal microbiota communities based on unweighted UniFrac index indicated that samples collected at the same time point were gathered together rather than between different treatment groups (Fig. 9b).

Fig. 9

Principal coordinate analysis of fungal communities based on (a) weighted and (b) unweighted UniFrac metrics of samples in different groups. Note: The values on the coordinate axis indicate the difference of the component in different groups. Group A, samples of 0 day (CK0); group B, samples of 30 day (CK1, 4.1, 7.1, 9.1); group C, samples of 90 day (CK2, 4.2, 7.2, 9.2)

Metagenome functional content prediction of bacterial community

In order to compare the potential functions of the bacterial community between different treatments, functional content was predicted from amplicon data using PICRUSt and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologs. When assessed with LEfSE at LDA > 2 and p value < 0.05, the microbiome of treatment groups were significantly enriched in eight functional categories compared to the control group (Table 2). These enriched functional categories were related to amino acid metabolism (e.g., alanine, aspartate and glutamate metabolism), biosynthesis of other secondary metabolites (e.g., tropane, piperidine, and pyridine alkaloid biosynthesis), metabolism of cofactors and vitamins (e.g., ubiquinone and other terpenoid-quinone biosynthesis), lipid metabolism (e.g., fatty acid biosynthesis), translation (e.g., translation factors), and carbohydrate metabolism (e.g., pentose phosphate pathway).

Table 2 Functional profile of the microbial community in different groups


Effects on the plant growth and accumulation of active ingredients

In the present study, we first developed the special manure of bacterial fertilizer and utilized it for cultivation of A. membranaceus in the field. The results suggest that the fertilizer can promote plant growth in the field which is consistent with numerous studies of nitrogen-fixing bacteria (Chi et al. 2010; Gage 2004; Gopalakrishnan et al. 2015; Rincon-Rosales et al. 2009; Zhiquan et al. 2016) as well as recent results from our laboratory using another type of bacterial fertilizer (Liang et al. 2019).

The main effective components of the A. membranaceus roots, including flavonoids, saponin, and polysaccharide were measured to evaluate the effect of bacterial manure on the accumulation of active ingredients in host. The result revealed that rhizobial bacteria of the fertilizer (strain Rhizobium sp. t16 and Sinorhizobium sp. t21) mainly interact with the plant root to increase the biosynthesis of secondary metabolites and the simultaneous inoculation of these two rhizobia will make this effect more obvious which indicates a special interaction between the two nitrogen bacteria. Previous studies have shown that inoculation with plant growth-promoting rhizobium could increase secondary metabolite content in marigold and Mentha piperita (Cappellari et al. 2017; Cappellari et al. 2013; Santoro et al. 2016), and our experiments have obtained consistent results. The molecular mechanism of the two nitrogen fixing bacteria in promoting the accumulation of secondary metabolites in host still needs further study.

Effects on the microbial abundance and community structure

Furthermore, the rhizosphere microbial abundance and community diversity during the plant growth were measured by Illumina MiSeq sequencing data. Diversity indices demonstrated an increase in bacterial community diversity with organic microbial fertilizer application, consistent with a previous report (Tabassum et al. 2017). However, this increase was not significant in this experiment. With the growth of plants, the abundance of indigenous bacteria in soil was inhibited, and some new strains have gained the dominant position. The dominant strains of the treatment group were also different compared with the control group. These results suggest that bacteria were suppressed during germination then promoted when seedlings unearthed, presumably due to root development, soil porosity, and sufficient oxygenation. Temperature, humidity, and biological interactions between microorganisms or between microbial and plant could inhibit microorganisms which have less tolerance to the environmental challenges, unlike adaptable microbes that thrive during A.membranaceus growth process (Tabassum et al. 2017). The result indicated that interspecific competition could contribute to the decline of the microbial community diversity, allowing adaptation of the microorganisms to changing environment for survive and eventually thrive (O’Sullivan et al. 2015). In the present study, the community diversity of fungi remained stable during both experimental and control groups, indicated that the fungal species could adapt to both environments and tolerate the transition between different periods.

In addition, the results of the relative abundance of rhizosphere bacteria throughout the entire A. membranaceus growth process shows three predominant phyla, Proteobacteria, Actinobacteria, and Acidobacteria, which is consistent with previous reports (Liu et al. 2015; Marsh et al. 2014). At genus level, the results showed that Pseudomonas, Lactobacillus, Corynebacterium, Bacillus, and Acinetobacter were the dominant bacteria during the host growth. Acinetobacter and Pseudomonas genera are able to degrade biosurfactants and hydrocarbons (Hamme et al. 2003). Bacillus species are confirmed an important role during the stable growth period of plant by secreting microbial inhibitors and plasminogen. Bacillus subtilis can produce antibacterial and bacteriolytic substances, such as subtilisin, organic acids, which can inhibit the growth and reproduction of pathogenic bacteria (Wang et al. 2019). Bacillus subtilis can also produce proteases, lipases, amylases, cellulases that promote the biochemical process during the A. membranaceus growth, which contribute to the accumulation of secondary metabolites, such as amino acids, organic acids, and aglycones (Yuan et al. 2019). Thus, these bacteria may play an important role in inhibiting the growth of pathogenic bacteria and promoting plant growth and accumulation of secondary metabolites through the secreted metabolites of these organisms during the stable growth period of A. membranaceus.

The abundance of the bacteria strain Rhizobium sp.t16 and Sinorhizobium sp.t21 were detected. Interestingly, their abundance in rhizosphere soil did not increase significantly after application of bacterial fertilizer. The results revealed that both rhizobial bacteria of the fertilizer mainly interact with the plant root to form special nodules to function in nitrogen fixation; thus, the distribution in the soil is relatively low and does not change significantly following application of fertilizer.

In this study, we also observed time point effects on the differentiation of bacterial community structure along with the plant growth, suggesting an interaction between plant development and the dynamics of rhizosphere bacteria community. The bacterial communities were relatively stable during the rapid growth period, but fluctuated considerably during the stable growth. Combining this with the physiological indices of A. membranaceus seedlings, we hypothesized that the nutrients consumed by the seedlings during the rapid growth period were mainly stored in seeds, and only a small part was applied to the soil, so there was no significant difference in the physiological characteristics of the seedlings of each group. Contrarily, during the stable growth period, the nutrient requirement rapidly increased. On one hand, special fertilizer could provide nutrients in time because it contains nutrients suitable for the growth of A. membranaceus (Liang et al. 2019). On the other hand, bacteria interact with root tissues by forming nodules on lateral roots, playing a role in nitrogen fixation and providing nutrients for plants. Meanwhile, plants secrete organic components through the root tips to regulate the microbial community structure of soil’s rhizospheres, and play a role in promoting growth (Badri and Vivanco 2009).

The functional content prediction of rhizosphere bacterial community was also conformed in our study. The functional approaches showed that significant shifts of metabolic functions impacting KEGG pathways in the organic microbial fertilizer treatment groups are related to metabolism and biosynthesis. These results provide a new perspective for further research on the relationship between microbiome and plants. More importantly, these functions of the bacterial community may play an important role in promoting the main medicinal components accumulation of A. membranaceus.

In addition to the bacterial community, the present study was also aimed to assess the fungal community structure. Shannon index results indicated that regulation of fungal community structure by the special fertilizer is also related to the growth stage of the plant. Microbial diversity of fungi in rhizosphere soil was inhibited during germination stage. During stable growth, there was a gradual development of the root system, with secretion of diverse compounds such as malic acid, sugars, amino acids, vitamins, various secondary metabolites, proteins, and mucilage to the soil to recruit beneficial microorganisms, which resulted in increased diversity of fungi in the soil’s rhizospheres.

Therefore, our study suggested that the special fertilizer developed in this research can significantly promote A. membranaceus growth and medical active components accumulation by regulating the community structure of rhizospheres in soil in different growth stages. This special fertilizer shows potential value in application for many crops, especially traditional Chinese medicinal herbs, due to its low cost and easy use. Moreover, our study implies a potential effect of the special fertilizer in restoring and rehabilitating degraded cropland soils by altering rhizosphere microbial communities.


In this study, a special bacterial fertilizer composed of stems and leaves powder of A. membranaceus and two native strains of rhizobia was applied in the planting of A. membranaceus, and the bacterial and fungal diversity was studied using high-throughput sequencing. The special bacterial fertilizer promoted plant growth and the main medicinal components accumulation by providing appropriate nutrients, promoting nitrogen-fixing activity, and regulating the rhizosphere microbial community. The special manure developed in this study has a lot of potential value for fertilizing crops as well as rehabilitating degraded cropland soils by altering the rhizosphere microbial communities.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. Ahemad M, Kibret M (2014) Mechanisms and applications of plant growth promoting rhizobacteria: current perspective. J King Abdulaziz Univ Ence 26:1–20

    Google Scholar 

  2. Badri DV, Vivanco JM (2009) Regulation and function of root exudates. Plant Cell Environ 32:666–681

    CAS  PubMed  Article  Google Scholar 

  3. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R (2010) Qiime allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. Cappellari LDR, Chiappero J, Santoro MV, Giordano W, Banchio E (2017) Inducing phenolic production and volatile organic compounds emission by inoculating mentha piperita with plant growth-promoting rhizobacteria. Entia Horticulturae 220:193–198

    CAS  Article  Google Scholar 

  5. Cappellari LDR, Santoro MV, Nievas F, Giordano W, Banchio E (2013) Increase of secondary metabolite content in marigold by inoculation with plant growth-promoting rhizobacteria. Appl Soil Ecol 70:16–22

    Article  Google Scholar 

  6. Chi F, Yang P, Han F, Jing Y, Shen S (2010) Proteomic analysis of rice seedlings infected by sinorhizobium meliloti 1021. Proteomics. 10:1861–1874

    CAS  PubMed  Article  Google Scholar 

  7. Chun-Hong LI, Ji T, Bing HE, Shun-Han X (2011) Determination of total saponins in astragalus membranaceus by ultraviolet spectroscopy. J Chongqing Med Univ 36:954–955

    Google Scholar 

  8. Crawford P, Crowley J, Sambandam N, Muegge B, Costello E, Hamady M (2009) Regulation of myocardial ketone body metabolism by the gut microbiota during nutrient deprivation. Proc Natl Acad Sci 106:11276–11281

    CAS  PubMed  Article  Google Scholar 

  9. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL (2006) Greengenes, a chimera-checked 16s rrna gene database and workbench compatible with arb. Appl Environ Microbiol 72:5069–5072

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. Dineen S, Aranda R, Anders D, Robertson J (2010) An evaluation of commercial DNA extraction kits for the isolation of bacterial spore DNA from soil. J Appl Microbiol 109:1886–1896

    CAS  PubMed  Article  Google Scholar 

  11. Edgar RC (2018) Accuracy of taxonomy prediction for 16s rRNA and fungal its sequences. Peerj. 6:e4652

    PubMed  PubMed Central  Article  Google Scholar 

  12. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R (2011) Uchime improves sensitivity and speed of chimera detection. Bioinformatics. 27:2194–2200

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. Ferreira CMH, Soares HMVM, Soares EV (2019) Promising bacterial genera for agricultural practices: an insight on plant growth-promoting properties and microbial safety aspects. Sci Total Environ 682:779–799

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  14. Francis I, Holsters M, Vereecke D (2010) The gram-positive side of plant–microbe interactions. Environ Microbiol 12:1–12

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. Gage DJ (2004) Infection and invasion of roots by symbiotic, nitrogen-fixing rhizobia during nodulation of temperate legumes. Microbiol Mol Biol Rev 68:280–300

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. Gopalakrishnan S, Sathya A, Vijayabharathi R, Varshney RK, Gowda CLL, Krishnamurthy L (2015) Plant growth promoting rhizobia: challenges and opportunities. 3 Biotech 5:355–377

    PubMed  Article  PubMed Central  Google Scholar 

  17. Hamme JV, Singh A, Ward O (2003) Recent advances in petroleum microbiology. Microbiol Mol Biol Rev 67:503–549

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. Heijden MGAVD, Bardgett RD, Straalen NMV (2010) The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol Lett 11:296–310

    Article  Google Scholar 

  19. James AE, Ko JKS (2007) Astragalus saponins induce growth inhibition and apoptosis in human colon cancer cells and tumor xenograft. Carcinogenesis. 28:1347–1355

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  20. Kemp P, Aller J (2004) Bacterial diversity in aquatic and other environments: what 16s rDNA libraries can tell us. FEMS Microbiol Ecol 47:161–177

    CAS  PubMed  Article  Google Scholar 

  21. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C (2013) Predictive functional profiling of microbial communities using 16s rRNA marker gene sequences. Nat Biotechnol 31:814–821

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. Liang H, Chen C, Chen J, Wang M, Zhi C, Haijun G, Jianping L (2019) Effect of nitrogen-fixing bacteria fertilizer on quality and yield of astragalus daodi in Shanxi province. J Shanxi Agric Sci 47:841–844

    Google Scholar 

  23. Liu B, Yuan J, Yiu S, Li Z, Xie Y, Chen Y (2012) Cope: An accurate k-mer-based pair-end reads connection tool to facilitate genome assembly. Bioinformatic. 28:2870–2874

    CAS  Article  Google Scholar 

  24. Liu J, Zhao Z, Chen H (2011) Review of astragali radix. Chinese Herb Med 3:90–105

    CAS  Google Scholar 

  25. Liu W, Xi X, Sudu Q, Kwok L, Guo Z, Hou Q (2015) High-throughput sequencing reveals microbial community diversity of Tibetan naturally fermented yak milk. Ann Microbiol 65:1741–1751

    CAS  Article  Google Scholar 

  26. Marsh AJ, O’Sullivan O, Hill C, Ross RP, Cotter PD (2014) Sequence-based analysis of the bacterial and fungal compositions of multiple kombucha (tea fungus) samples. Food Microbiol 38:171–178

    CAS  PubMed  Article  Google Scholar 

  27. O’Sullivan D, Fallico V, O’Sullivan O, McSweeney P, Sheehan J, Cotter P (2015) High-throughput DNA sequencing to survey bacterial histidine and tyrosine decarboxylases in raw milk cheeses. BMC Microbiol 15:1–12

    Article  CAS  Google Scholar 

  28. Richard J, Manuela P, Antonella S, Anna K, Stanislav K (2017) The role of soil microorganisms in plant mineral nutrition—current knowledge and future directions. Front Plant Sci 8:1617

    Article  Google Scholar 

  29. Rincon-Rosales R, Lloret L, Ponce E, Martínez-Romero EM (2009) Rhizobia with di¡erent symbiotic e/ciencies nodulate acaciella angustissima in Mexico, including sinorhizobium chiapanecum sp. Nov. Which has common symbiotic genes with Sinorhizobium mexicanum. FEMS Microbiol Ecol 67:103–117

    CAS  PubMed  Article  Google Scholar 

  30. Santoro MV, Bogino PC, Natalia N, Rosario CLD, Giordano WF, Erika B (2016) Analysis of plant growth-promoting effects of fluorescent pseudomonas strains isolated from mentha piperita rhizosphere and effects of their volatile organic compounds on essential oil composition. Front Microbiol 7:1085

    PubMed  PubMed Central  Article  Google Scholar 

  31. Savci S (2012) Investigation of effect of chemical fertilizers on environment. Apcbee Proc 1:287–292

    CAS  Article  Google Scholar 

  32. Schloss P, Westcott S, Ryabin T, Hall J, Hartmann M, Hollister E (2009) Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537–7541

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C (2011) Metagenomic biomarker discovery and explanation. Genome Biol 12:R60

    PubMed  PubMed Central  Article  Google Scholar 

  34. Sharma U, Paliyal SS (2014) Effect of long-term use of fertilizers and organics on phosphorus dynamics in rainfed maize-wheat cropping system. J Indian Soc Soil ence 62:216–223

    Google Scholar 

  35. Singh JS (2015) Plant–microbe interactions: a viable tool for agricultural sustainability. Appl Soil Ecol 92:45–46

    Article  Google Scholar 

  36. Tabassum B, Khan A, Tariq M, Ramzan M, Aaliya K (2017) Bottlenecks in commercialisation and future prospects of pgpr. Appl Soil Ecol 121:102–117

    Article  Google Scholar 

  37. Wang W, Niu J, Zhou X, Wang Y (2011) Long-term change in land management from subtropical wetland to paddy field shifts soil microbial community structure as determined by plfa and t-rflp. Pol J Ecol 59:37–44

    CAS  Google Scholar 

  38. Wang Z, Li Y, Zhuang L, Yu Y, Liu J, Zhang L, Gao Z, Wu Y, Gao W, Ding G-C, Wang Q (2019) A rhizosphere-derived consortium of bacillus subtilis and trichoderma harzianum suppresses common scab of potato and increases yield. Comput Struct Biotechnol J 17:645–653

    PubMed  PubMed Central  Article  Google Scholar 

  39. Xi X, Wei X, Wang Y, Chu Q, Xiao J (2010) Determination of tea polysaccharides in camellia sinensis by a modified phenol-sulfuric acid method. Arch Biol Sci 62:669–676

    Article  Google Scholar 

  40. Yang B, Xiao B, Sun T (2013) Antitumor and immunomodulatory activity of astragalus membranaceus polysaccharides in h22 tumor-bearing mice. Int J Biol Macromol 62:287–290

    CAS  PubMed  Article  Google Scholar 

  41. Yang ZC, Zhou LD, Lv YZ, Li H (2012) Long-term effects of crop residual and inorganic fertilizers on yield and soil organic matter for a winter wheat-maize system in north china plain. Adv Mater Res 356-360:2523–2530

    CAS  Article  Google Scholar 

  42. Yuan K, Song P, Li S, Gao S, Wen J, Huang H (2019) Combining metabolic flux analysis and adaptive evolution to enhance lipase production in bacillus subtilis. J Ind Microbiol Biotechnol 46:1091–1101

    CAS  PubMed  Article  Google Scholar 

  43. Zheng Y, Li C, Feng S, Huang X, Liu Y, Chen X, Di D (2011) Study on the content determination of total flavonoids in Olea europaea l. Leaves. Guang Pu Xue Yu Guang Pu Fen XI 31:547–550

    CAS  PubMed  Google Scholar 

  44. Zhiquan X, Zhongwei T, Hao L, Ran Z, Jianping L (2016) Isolation and application of nitrogen-fixing bacteria in rhizosphere of astragalus membranaceus bunge in Shanxi. J Shanxi Agric Univ (Nat Sci Ed) 36:483–488

    Google Scholar 

  45. Zhou R, Chen H, Chen J, Chen X, Wen Y, Xu L (2018) Extract from astragalus membranaceus inhibit breast cancer cells proliferation via pi3k/akt/mtor signaling pathway. BMC Complement Altern Med 18:83

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references


Not applicable


This study was supported by the Key Technologies Research and Development Program (2019YFC1710800), Shanxi Astragalus Resource Industrialization and Industrial Internationalization Collaborative Innovation Center Project (HQXTCXZX2016-004), Shanxi Provincial Key Research and Development Project (CN) (201903D321204) and Science and Technology Innovation Funds of Shanxi Agricultural University (2018YJ31).

Author information




JL and ZS designed the study and wrote the manuscript. ZX, ZY and JN performed the experiments. ZX and ZS analyzed the data. ZC revised the manuscript. The authors read and approved of the final manuscript.

Corresponding author

Correspondence to Zhi-Yong Shi.

Ethics declarations

Ethics approval and consent to participate

This article does not contain any studies with human participants or laboratory animals.

Consent for publication

Not applicable

Competing interests

The authors declare that there are no conflicts of interests.

Additional information

Publisher’s Note

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

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liang, JP., Xue, ZQ., Yang, ZY. et al. Effects of microbial organic fertilizers on Astragalus membranaceus growth and rhizosphere microbial community. Ann Microbiol 71, 11 (2021).

Download citation


  • Astragalus membranaceus
  • Microbial organic fertilizer
  • High-throughput sequencing
  • Rhizosphere microbial community
  • Diversity