Significant differences in garlic growth performance between adjacent fields
Analysis of garlic plant height revealed that the mean plant height at site H was 1.44-fold greater than that at site L (Wilcoxon rank-sum test, P < 0.001) (Fig. 1B–C). To investigate the underlying factors contributing to this discrepancy, we characterized the physicochemical properties and microbial community compositions of rhizosphere soils from both sites.
The AP content was significantly greater at site H than at site L
To determine whether rhizosphere soil physicochemical properties differed between sites H and L, we measured key parameters in garlic rhizosphere soil samples. Among the analyzed properties, the available phosphorus (AP) content at site H was significantly greater than that at site L (Wilcoxon rank-sum test, P < 0.05), whereas no significant differences were detected in pH, EC, TN, TP, NH₄⁺-N, NO₃⁻-N, SOC, or the C/N ratio (all *P > 0.05) (Fig. 2). These results highlight AP content as a key factor contributing to the observed differences in garlic growth between the two sites.
Significant differences in rhizosphere microbial community diversity
To evaluate bacterial and fungal diversity between sites H and L, α diversity indices were calculated for bacteria, core bacterial taxa (defined as ASVs with a mean relative abundance > 0.1% across samples), fungi, and core fungi. The Shannon index of the core bacteria at site H was significantly greater than that at site L (Wilcoxon rank-sum test, P < 0.05) (Fig. 3C). No significant differences were detected in bacterial, fungal, or core fungal α diversity (all P > 0.05) (Fig. 3A, B, D).
PCoA based on Bray–Curtis dissimilarities revealed distinct clustering patterns for bacterial and core bacterial communities between sites (PERMANOVA, P < 0.001) (Fig. 3E). In contrast, the fungal and core fungal communities presented no significant structural differences (PERMANOVA, P > 0.05) (Fig. 3F). Additionally, the exogenous nutrient levels were consistent between the two fields. These results suggest that bacterial community diversity and composition are critical factors contributing to the observed garlic growth disparities between sites H and L.
Higher stability of bacterial communities at location H
To evaluate rhizosphere microbial network complexity and stability, ecological networks were constructed for bacteria and fungi at sites H and L. Network analysis revealed that bacterial and fungal networks at site H presented greater complexity and connectivity than did those at site L, characterized by more nodes, a higher average degree, a larger network diameter, greater centralization, and greater closeness centralization (Fig. 4A–B; Table 1). Specifically, the bacterial network at site L presented an increased average path length, indicating decreased responsiveness to external perturbations (Table 1), whereas the fungal network at site L presented a decreased average path length, reflecting faster adaptive capacity (Table 1). Site L also harbored more keystone nodes, suggesting that garlic actively recruits critical microbes under stress. Community stability, measured by the AVD index, was significantly greater at site H for both bacteria and fungi (Fig. 4D). These results underscore bacterial community stability as a key determinant of garlic growth disparities between sites.
Table 1
Network properties of bacteria and fungi
| |
Bacteria_H
|
Bacteria_L
|
Fungi_H
|
Fungi_L
|
|
Nodes
|
218
|
219
|
88
|
88
|
|
Edges
|
2224
|
2048
|
370
|
367
|
|
Connectance
|
0.094026128
|
0.085794479
|
0.096656217
|
0.095872518
|
|
Average Degree
|
20.40366972
|
18.70319635
|
8.409090909
|
8.340909091
|
|
Average Path Length
|
1.879037574
|
1.985647338
|
2.298343637
|
2.251869175
|
|
Diameter
|
4.023809524
|
3.85206652
|
6.111667761
|
5.307278838
|
|
Centralization Degree
|
0.104130554
|
0.083930292
|
0.13322884
|
0.122518286
|
|
Centralization Closeness
|
0.141528417
|
0.113956041
|
0.195748564
|
0.163661951
|
|
The Number of Keystone Nodes
|
36
|
41
|
7
|
11
|
The AP content is significantly correlated with rhizosphere bacterial community composition
RDA was used to explore the influence of abiotic factors, specifically soil physicochemical properties, on bacterial and fungal communities. The first two axes of the RDA explained 30.63% of the total variation in the bacterial communities. Among these factors, AP had the most significant impact on the bacterial community structure (R2 = 0.7638, P = 0.009) (Fig. 5A, Table 2). The RDA results revealed that differences in AP content affected the distribution of bacterial communities, whereas the distribution of fungal communities was not influenced by these abiotic factors (Fig. 5A, B, Table 2). These findings suggest that AP and bacterial communities are crucial factors contributing to the significant differences in garlic growth between sites H and L.
Table 2
Results of RDA permutation test for bacteria, fungi communities and soil physical and chemical properties
|
Name
|
Bacteria
|
|
Fungi
|
|
R2
|
Padj
|
Individual (%)
|
|
R2
|
Padj
|
Individual (%)
|
|
PH
|
0.179336
|
0.5076
|
1.04
|
|
0.010451
|
0.901
|
0
|
|
EC
|
0.002077
|
0.992
|
0
|
|
0.04049
|
0.858857
|
0
|
|
NH4+-N
|
0.057199
|
0.892286
|
0
|
|
0.343983
|
0.711
|
3.36
|
|
NO3−-N
|
0.321534
|
0.306
|
1.95
|
|
0.106433
|
0.858857
|
1.09
|
|
AP
|
0.763771
|
0.009
|
3.8
|
|
0.069747
|
0.858857
|
0
|
|
TN
|
0.139292
|
0.582
|
0
|
|
0.100856
|
0.858857
|
0
|
|
TP
|
0.180848
|
0.5076
|
0
|
|
0.101308
|
0.858857
|
0
|
|
C/N
|
0.20139
|
0.5076
|
0.15
|
|
0.140689
|
0.858857
|
1.08
|
|
SOC
|
0.007786
|
0.992
|
0
|
|
0.012459
|
0.901
|
0.16
|
Location H is enriched with numerous microorganisms associated with P transformation
To quantify the contributions of soil physicochemical properties and microbial communities to garlic plant height at sites H and L, we performed hierarchical partitioning analysis. The results revealed that the bacterial community had the greatest contribution to garlic plant height, followed by AP (Fig. 6E). To further clarify the role of bacterial communities in determining garlic plant height, we analyzed the species composition of rhizosphere microbial communities at both sites. On the basis of the ASV classification results, Sphingomonas was the most abundant genus among the rhizosphere bacteria at both sites. However, the second most abundant genera differed: Flavobacterium at site H and Bacillus at site L. The most abundant fungus at site H was Mortierella, whereas at site L, it was Fusarium (Fig. 4A, B), a common soil-borne plant pathogen. To identify the bacterial and fungal genera that differed between the microbial communities at sites H and L, we conducted linear discriminant analysis effect size (LEfSe) analysis. Significant differences were found in the enriched fungi and bacteria between the two sites (Fig. 6A, B). At site H, 36 bacterial taxa met the score criteria, whereas 22 bacterial taxa met the score criteria at site L (Fig. 6C). We also found that the main enriched genera at site H, Arthrobacter and Thauera, are typically involved in the transformation of minerals such as P and sulfur, which promotes plant nutrient absorption and inhibits plant pathogens (Fig. 6C). Among the fungal taxa, 5 groups at site L and only 1 group at site H met the linear discriminant analysis (LDA) score criteria (Fig. 6D). Additionally, the enriched Penicillium at site L can cause stem-base rot and fruit rot in vegetables and fruits (Fig. 6D). Overall, these results suggest that the differences in garlic growth at sites H and L are influenced mainly by differences in the composition of rhizosphere microorganisms.