3.1.Species annotation and alpha diversity analysis
A total of 6 composite samples, including 3 soil fungi (G, H and Q groups) and 3 soil bacteria (L, K and J groups) samples (with three biological replicates per sample), were subjected to high-throughput sequencing using the Illumina NovaSeq system. A total of 1542728 raw reads and 1438245 effective tags, including 764899 raw reads and 724099 effective tags from cabbage-cultivated soil microorganism (average: 84989 reads per sample) and 777829 raw reads and 714146 effective tags from pepper-cultivated soil microorganism (average: 864265 reads per sample). Based on the least number of sequences detected among the assayed samples, we randomly selected and compared 80455 and 79350 effective reads of cabbage-cultivated soil microorganism and pepper-cultivated soil microorganism per sample. Results from Fig. 1AC divulged that 512, 655, 723 fungus total OTUs as soon as 5524, 5334 and 5301 bacteria total OTUs were identified in control (G), 0~20 cm (H) and 20~40 cm (Q), respectively. 698, 766, 757 fungus total OTUs were identified while 5373, 5671 and 5624 bacteria total OTUs of was observed in control (L), 0~20 cm (K) and 20~40 cm (J). This data suggests that cabbage and pepper increase the total fungus OTUs as comparison to the control, but the effect of the former on the fungus total OTUs from deep soil is more significant than the later (Table 1).
Table 1. Data preprocessing statistics and quality control statistics
|
Sample Name
|
Raw Tags
|
Clean Tags
|
Effective Tags
|
Effective%
|
G.1
|
93,862
|
93,262
|
88,860
|
93.04
|
G.2
|
80,537
|
80,014
|
77,791
|
94.95
|
G.3
|
83,209
|
82,719
|
77,140
|
90.82
|
H.1
|
87,112
|
86,664
|
83,051
|
93.95
|
H.2
|
80,365
|
79,916
|
76,825
|
93.46
|
H.3
|
83,718
|
83,258
|
79,329
|
93.77
|
Q.1
|
88,299
|
87,788
|
82,271
|
91.96
|
Q.2
|
82,140
|
81,630
|
76,204
|
90.55
|
Q.3
|
85,657
|
85,233
|
82,628
|
94.72
|
Subtotals
|
764,899
|
760,484
|
724,099
|
|
Average
|
84,989
|
84,498
|
80,455
|
|
L.1
|
87,589
|
86,784
|
68,828
|
78.19
|
L.2
|
81,029
|
80,641
|
76,255
|
92.59
|
L.3
|
87,058
|
86,679
|
82,701
|
94.42
|
J.1
|
91,492
|
91,053
|
81,935
|
89.3
|
J.2
|
90,893
|
90,430
|
85,243
|
92.37
|
J.3
|
87,872
|
87,373
|
82,288
|
93.15
|
K.1
|
83,667
|
83,257
|
77,891
|
91.27
|
K.2
|
79,009
|
78,653
|
75,016
|
94.01
|
K.3
|
89,220
|
88,730
|
83,989
|
92.98
|
Subtotals
|
777829
|
773600
|
714146
|
|
Average
|
86425
|
85956
|
79350
|
|
Totals
|
1542728
|
1534084
|
1438245
|
|
G, H, Q means these samples from the control, 0~20 cm soil layer and 20~40 cm layer in cab-bage-cultivated soil while L, K and J presented these samples from the control, 0~20 cm soil layer and 20~40 cm layer in pepper-cultivated soil, respectively. 1, 2 and 3 represent the repetition.
Cabbage decreased the bacteria total OTUs while pepper increased that. In addition, H and Q groups contained 158 and 226 fungus specific OTUs, and K and J groups contained 189 and 224 fungus specific OTUs. Similarly, H and Q groups contained 703 and 770 bacteria specific OTUs, while K and J groups contained 731 and 876 bacteria specific OTUs (Fig. 1AB), indicating that both cultivation of cabbage and pepper increase the specific OTUs in soil.
To determine whether the amount of sequencing data generated was reasonable we constructed dilution curves by plotting the amount of extracted data (y-axis) against the alpha diversity index (x-axis), which reflects the abundance and diversity of a given species in the sample (Fig. 1C). The alpha diversity index was determined the Shannon index, which is often used along with the observed_species to reflect the alpha diversity index; Similar, the greater the observed_species, the higher the community diversity (Fig. 2A). The larger the value of the Shannon index, the higher the community diversity (Fig. 2B). Our results showed the diversity of microorganisms was enhanced after cabbage and pepper were cultivated. Microorganism diversity of shallow soils is more abundant than that of deep soil in cabbage-cultivated soil, but the contrary results were observed in pepper-cultivated soil.
3.2. Beta diversity and phylogenetic information visualization
Analysis of weighted unifric and unweighted unifric is a non-parametric test used to test whether the differences between groups (two or more) are significantly greater than those within a group (Fig. 3AB). In addition, we also evaluated beta diversity at the level of OTUs, which were defined based on a similarity cut-of of 97%. To compare the composition of the identified microbial communities within different plant compartments, Bray-Curtis dissimilarity was performed at the OTU level. The smaller the P-value, the greater the significance of the differences between groups (Fig. 3CD). A statistically significant P-value (P<0.05) and change of weighted unifric and unweighted unifric after cabbage and pepper were planted indicates that the grouping is meaningful. Principal coordinate analysis (PCA) of measured samples showed that all groups were dispersed as comparison to the control (Fig. 3E).
Phylogenetic visualization analysis of the similarity among the fungus samples showed that excluding the unclassified OTUs, a total of 89 genera belonging to 8 families were identified, and 67 classes genera from Ascomycota family, and 14 classes is distributed into the Basidiomycota. Other classes only possessed a little percentage. Cladosporium from Ascomycota possessed the higher percentage than others in all groups, and percentage of Cladosporium was decreased significantly when cabbage and pepper were planted in measured groups (Fig. 4A). Phylogenetic visualization analysis of the similarity among the six samples showed that 100 genera belonging to 14 families were identified, and 38 classes is presented in Proteobacteria family. 16, 11 and 12 classes were distributed in Myxococcota, Synergistota and Actinobacteria families, respectively. Sphingomonas from Proteobacteria family and Pseudarthrobacter from Actinobacteria family possessed the higher percentage than other (Fig. 4B).
3.3. Correlation between the relative abundance of dominant species and collinearity
Fig. 5A shows the microbial species composition and relative abundance (mean relative abundance>1%) at the phylum level in fungus groups. Excluding unclassified and others, the six groups were mainly composed of fungi belonging to the phyla Ascomycota, Basidiomycota, Mortierellomycota and Chrytridiomycota, with the former three phyla representing the core fungal communities. In all measured groups, Ascomycota possessed the highest relative abundance of fungus than any other, suggesting Ascomycota is the primary fungus species in soil. However, Although Ascomycota was dominant at all groups, its relative abundance gradually decreased with the increased deep of soil. Contrarily, increased relative abundance in Mortierellomycota and others were observed. This revealed that pepper and cabbage weakened dominance of Ascomycota and promote the balance of fungal distribution with increased soil depth.
Fig. 5B shows the microbial species composition and relative abundance (mean relative abundance>1%) at the genus level in bacteria groups. The six groups were mainly composed of Proteobacteria, Actinobacteria, Bacteroidota, Myxococcota, Firmicutes, Acidobacteriota, Verrucomicrobiota, Synergistota, Fusobacteriota, Actinobacteriota, Gemmatimonadetes, Crenarchaeota, Chloroflexi and Gemmatimonadota. In all measured groups, Proteobacteria was the most dominant genus, followed by Firmicutes, Hygrocybe, Actinobacteriota, and Synergistota. Pseudarthrobacter from Actinobacteria was the most dominant genus, followed by Sphingomonas from Proteobacteria and Bacillus from Firmicutes. The relative abundance of these microbial species is decreased when cabbage and pepper were cultivated. This data suggests that pepper and cabbage weakened expression of some dominant bacteria and encouraged the development of other weak microbial communities.
3.4. Analysis of different flora and differentially expressed abundance
LEfSe is used to determine the genetic or functional features that best explain the differences between two or more groups of samples under different biological conditions or environments as well as the degree to which these features influence the differences between groups. Our results showed that Asmcomycota was the main fungal communities in cabbage-cultivated fungus groups (Fig. 6A); Mortierellomycota in pepper-cultivated fungus groups (Fig. 6B); Hymenobacteraceae in cabbage-cultivated bacteria groups (Fig. 6C); Firmicutes, and Acidobacteria in pepper-cultivated bacteria groups (Fig. 6D).
Expression abundance of 35 primary classes were further determined, and found that abundance of Cladosporium from Ascomycota is reduced among all species after when crops were planted. In addition, species of increased expression possessed more ratio than the control among measured fungus (Fig. 7A). Similar results also were observed in bacterial. No new species were observed in the measured microorganisms (Fig. 7B). The data indicated that as a whole cabbage and pepper promote the diversity to change the relative expression of microorganism but not to generate new species. In addition, 8 classes of fungus, including Mortierellomycota, Basidiomycota, Zoopagomycota, Ascomycota, Chytridiomycota, Rozellomycota, Basidiobolomycota and Aphelidiomycota, generate the interaction in measured soil. Ascomycota and Chytridiomycota showed the direct obvious interact relationship (Fig. 7C). The relationship also was observed between Ascomycota and Basidiobolomycota. 25 classes of bacteria showed the interaction relationship in soil. Except the unidentified species, Proteobacteria, Acidobacteria and Actinobacteria are the main species generating the interaction (Fig. 7D).