Effects of pineapple rotation and residue amendment on banana Fusarium wilt disease incidence and relative abundance of Fusarium
The effects of different concentrations of pineapple and banana residues (0.1%, 1.0% and 10.0%) on spore germination were evaluated (Fig. 2A). The infusion of all parts of banana had a promoting effect on spore germination of the pathogen, while the infusion of all parts of pineapple had an inhibiting effect on spore germination, with the fruit infusion of pineapple having a highly significant effect on spore germination. The effect of pineapple and banana root exudates on spore germination, compared to banana root exudates, pineapple root exudates were significantly less able to promote spore germination of the pathogen. The fungus significantly increased the number of pathogenic spores compared to the control (Fig. 2B). We observed an overall significantly reduced Fusarium wilt disease incidence in both the pineapple rotation and residue amendment systems (Duncan’s t-test, p < 0.001) (Pr and Pp), which was significantly lower than those in the fallow and monoculture treatments (p < 0.001) (Fig. 2C). Moreover, the abundance of rhizosphere Fusarium (rotation treatment) (p = 0.001), as indicated by the residue treatment bulk Fusarium (p = 0.118), and the relative abundance of rhizosphere Fusarium (p < 0.001) were both positively correlated with banana Fusarium wilt DI (Fig. 2D), demonstrating disease suppression ability after rotation and residue addition with pineapple in the banana orchard. In the rhizosphere soil, the relative abundance of Fusarium was significantly lower in the rotation than in the monoculture treatments (p < 0.05, Fig. 2C). The qPCR results showed that significantly fewer Fusarium oxysporum were detected in the rotation rhizosphere (Duncan’s t-test, p < 0.05, Fig. S1).
Effects of pineapple rotation and residue amendment on soil microbial community structure
The PCoA plots (Fig. 3A, B, C and D) showed significant differences in the bacterial and fungal community composition in the rotation and residue amendment systems (PERMANOVA, p < 0.001), in the rotation system, there was a significant difference in bacteria (p (bulk) < 0.01, p (rhizosphere) < 0.01), moreover, the different rhizosphere soil microbial communities changed significantly from the corresponding initial soil microbial communities, and the microbial communities of different plants grown under the same soil type also differed significantly (p < 0.001) (Fig. S2). The sequencing results were analysed and detailed, and the microbial community richness and diversity sequencing results are shown in the Supplementary Material (Table S3). PERMANOVA and VPA analyses revealed that soil rhizosphere microbial communities were significantly influenced by soil type (p < 0.001) and soil pH (p < 0.001). The relative importance of soil type and its soil pH on soil and rhizosphere microbial communities was greater than that of plant species (Fig. 3 E, F, G and H).
Effect of pineapple rotation and residue amendment on taxonomic composition
Volcano plot analysis of the sequence results revealed rotation and residue amendment bacterial and fungal community compositions with specific respective sets of OTUs (Fig. 4A and B). With red indicating bulk soil and black indicating rhizosphere soil, we found both that the pineapple rotation residue amendment treatment had a higher OTU number in the bulk soil, and selecting an OTU with a relative abundance greater than 0.1% for Venn analysis revealed that rotation and residue treatments had more bacteria, while fungi in the bulk soil had significantly higher OTU numbers than those in the rhizosphere soil (Fig. 4C and D).
At the phylum level, the pineapple rotations and residue amendment, which had the same composition, consisting mainly of the bacterial phyla Proteobacteria, Acidobacteria, Bacteroidetes, Firmicutes, Gemmatimonadetes, Actinobacteria, Verrucomicrobia, Chloroflexi, Planctomycetes, Cyanobacteria/Chloroplast, and Thaumarchaeota and the fungal phyla Ascomycota, and Basidiomycetes, had the most abundant phyla in all samples (Fig. 5A and B). The pathogens were significantly negatively correlated with the bacterial phyla Thaumarchaeota (-0.482**), Nitrospirae (-0.401*), Firmicutes (-0.515**), and Ascomycota (-0.608**) (Table S5), which were generally consistent with previous findings on disease-suppressing soils (Shen et al. 2018a).
Analysis of the variability of the top 30 genera of bacteria and fungi showed that the variability of genera significantly increased in the soil. In the rotation, we found that the Burkholderia, Bacillus, Rhizobium, Sphingosinicella, Pseudomonas and Talaromyces genera were significantly increased in rhizosphere soils and in bulk soil, along with Gemmatimonas, Gp6, Nitrososphaera, Gp7, Penicillium and Mortierella (Fig. 5C). Among them, in the pineapple rotation, there were significant negative correlations with the number of Fusarium, such as Burkholderia spp. (Pearson: -0.558**, Spearman: -0.443**), Pseudomonas spp. (Pearson: NS, Spearman: -0.363*), and Talaromyces spp. (Pearson: -0.687**, Spearman: -0.530**) (Tables S6 and S8).
In the residue amendment, we discovered that the Pseudomonas, Sphingobium, Azohydromonas, Bacillus, Georgfuchsia, Rhizobium, Pseudoduganella, Klebsiella, Sphingomonas, Elaphocordyceps and Penicillium genera were significantly increased in rhizosphere soils, in the bulk soil, the Gemmatimonas, Gp6, Nitrososphaera, Opitutus, Gp3, GpXIII, Aspergillus and Chaetomium genera were significantly increased (Fig. 5D). Among them, in the residue treatment, there were significant negative correlations with the number of Fusarium, such as Pseudomonas spp. (Pearson: -0.524**, Spearman: -0.440**), Aspergillus spp. (Pearson: -0.378*, Spearman: -0.354*), and Penicillium spp. (Pearson: -0.453**, Spearman: -0.540**) (Tables S7 and S8).
Bulk and rhizosphere network construction through effects on specific microbial taxa
In this study, we constructed co-occurrence networks using random matrix theory (RMT) to determine the differences in bacterial and fungal assemblages (OTU relative abundance > 0.1%) in bulk and rhizosphere soils of the different treatments. All values of the calculated modularity index were larger than 0.4 (Table S9), suggesting typical module structures (Chen et al., 2020). Overall, pineapple rotation and residue amendment showed marked effects on the soil microbial network: the average path distance (GD), the average clustering coefficient (avgCC) and the modularity of the empirical networks were higher than those of the corresponding, identically sized random networks (Table S9). Here, we found that residue assemblages (in Fig. 6 C and D) formed more connected and more complex networks with fewer nodes but more connections (edges) between nodes compared with the bulk soil. There were many keystone taxa in the microbial communities whose removal could cause a dramatic shift in microbiome structure and function. Keystone taxa in network analysis can be computationally identified as hubs with a high within-module degree Zi (Zi ≥ 0.5 indicates that the nodes are “well connected” to other nodes in the module). The PBr and PRr treatments had some keystone taxa, such as Burkholderia and Pseudomonas, and no hub was found in the bulk or rhizosphere soil (Cf, Cn, Bm and Bb) (Fig. S3).
Effects of the soil properties, number of F.oxysporum, the bulk and rhizosphere microbial communities and key microorganism on banana Fusarium wilt disease incidence
To investigate the potentially relative important suppression predictors of banana Fusarium wilt disease incidence, we then used Linear models (LM) to identify the potential positive or negative effects of the number of Fusarium (including the number of F.oxysporum and the relative abundance of Fusarium), bacterial and fungal communities (including bulk and rhizosphere communities), and key microorganism Burkholderia, Pseudomonas, Talaromyces and Pseudomonas, Penicillium, Aspergillus (bulk and rhizosphere soil) in the crop rotation and residue amendment system on banana Fusarium wilt disease incidence, respectively (Table 1).
For microbial linear model in the banana-pineapple rotation system, importantly, Fungal-pcoa1 (rhizosphere) (F = 23.41, p = 0.001, Relative Importance = 10.89%), Fusarium relative abundance (rhizosphere) (F = 201.74, p < 0.001, Relative Importance = 12.70%), Burkholderia (rhizosphere) (F = 3.56, p = 0.092, Relative Importance = 13.79%), Talaromyces (bulk) (F = 64.34, p < 0.001, Relative Importance = 23.61%) and Talaromyces (rhizosphere) (F = 1.41, p = 0.265, Relative Importance = 19.09%) constrained disease incidence the most (with a relative importance more than 10%). Besides, based on linear regression analyses between disease incidence and selected microbial indicators, Fungal-pcoa1 (rhizosphere) (p =0.001), Fusarium relative abundance (bulk) (p < 0.001), Fusarium relative abundance (rhizosphere) (p < 0.001) and Talaromyces (bulk) (p < 0.001) have significant negative relationship to disease incidence (Table 1).
For microbial linear model in the pineapple residue amendment system, importantly, Fusarium relative abundance (rhizosphere) (F = 363.96, p < 0.001, Relative Importance = 25.49%) and Aspergillus (rhizosphere) (F = 16.15, p = 0.003, Relative Importance = 40.83%) constrained disease incidence the most (with a relative importance more than 10%). Besides, based on linear regression analyses between disease incidence and selected microbial indicators, Penicillium (bulk) (p =0.051) and Aspergillus (rhizosphere) (p = 0.003) have significant negative relationship to disease incidence (Table 1).
For physicochemical linear model, the content of soil organic matter (F = 8.94, p = 0.011, Relative Importance = 22.22%), the available potassium (F = 19.45, p = 0.001, Relative Importance = 36.45%) and the available phosphorus (F = 89.40, p < 0.001, Relative Importance = 16.37%), the available potassium (F = 223.59, p < 0.001, Relative Importance = 28.88%) and the available nitrogen (F = 7.69, p = 0.016, Relative Importance = 45.10%) constrained disease incidence the most (with a relative importance more than 10%) in the rotation and residue amendment system, respectively (Table S10).
Based on the above results, a conceptual model illustrating potentially path with important suppression predictors in intercropping system was constructed (Fig. 7). The conceptual mode indicated that two ways in banana plantations reduced the relative abundance of Fusarium by the soil characteristics and microbial community structure regulation. Among all the suppression predictors, the key physicochemical factors AP (in rotation system) and OM (in residue system) contents are significantly affected by pineapple rotation and rseidue amendment, respectively. And they were significant leading to changes of fungal community and beneficial microorganism. And the significant increases in bacterial genera (Burkholderia) and fungal genera (Talaromyces and Aspergillus) can directly affect the relative abundance of Fusarium, thereby reduce the incidence of banana (Fig. 7).