Physicochemical indexes and variation during composting
Total phosphorus (TP) content was not different between WJP and CK (Fig. 1). From days 0 to 21, the TP contents of CaP and CaP + WJP increased. The content of resin phosphorus increased in all treatments. The content of the resin phosphorus of CaP treated was higher than that of WJP and CaP + WJP after 42 days. The total phosphorus content of sodium bicarbonate treated with CaP increased the fastest, while the inorganic phosphorus content of sodium bicarbonate treated with WJP increased the fastest from days 21 to 42. The inorganic phosphorus content of sodium hydroxide increased throughout the treatment cycle, and the increase of the three composting treatments was greater than that of CK. The content of total sodium hydroxide phosphorus was similar to that of sodium hydroxide inorganic phosphorus, which generally increased (Fig. 1). Compared to control, EC of all treatments decreased at first and then increased on the 14th day, and the CaP treatment was the highest. On the 21th and 35th days, both WJP and CaP increased at first and then decreased. After 42 days, EC of CK and CaP + WJP treatments increased, while the WJP and CaP showed the opposite trend. The pH range of the three compost treatments and the control ranged from 6.12–7.91 (Fig. S2).
Effect of WJP and CaP on bacterial and fungal communities during composting
From the compost samples of different composting stages (days 0, 21, and 42) and different treatments (CK, CaP, WJP, and CaP + WJP), 320,945 bacterial sequences were obtained, which were divided into 15759 ASVs (Table S1). A total of 320,945 fungal sequences were obtained, which were divided into 2457 ASVs (Table S2). According to the rarefaction curve and rank abundance curve, the changes in the bacterial communities in the compost were effectively reflected in the high-throughput sequencing result (Fig. S3). At the family level, composting significantly changed bacterial and fungal communities (Fig. 2A and Fig. 2B). For the CK compost group, the levels of dominant bacteria at on day 0 were Enterobacteriaceae (6.72%), Burkholderiaceae (5.70%), and Pseudomonadaceae (3.8%). The dominant bacteria on day 21 were Enterobacteriaceae (23.98%), Sphingobacteriaceae (13.88%), and Rhizobiaceae (11.13%). The dominant bacteria on day 42 were Sphingobacteriaceae (15.27%), Rhizobiaceae (11.93%), and Enterobacteriaceae (9.39%) (Fig. 2A). For the WJP-treated compost group, the dominant bacteria on day 21 were Enterobacteriaceae (28.93%), Actinomycetaceae (10.17%), and Sphingobacteriaceae (8.21%). However, the proportions of these three bacteria were 5.09%, 19.55%, and 5.42%, respectively on day 42. Rhizobiaceae increased from 6.63–10.25%, and Flavobactereae increased from 4.33–10.72%. For the CaP-treated compost group, Enterobacteriaceae decreased from 36.87% on day 21 to 10.59% on day 42, and Xanthomonadaceae decreased from 11.02–4.65%. Actinomycetaceae increased from 7.02% on day 21 to 13.46% on day 42, and Flavobacteriaceae increased from 1.65–14.23%. For the CaP + WJP-treated compost group, Enterobacteriaceae decreased from 12.92–6.48%. For fungi, Dipodascaceae and Pichiaceae were the main dominant fungi in the composting process (Fig. 2B). The bacterial Chao1 value of compost samples with different treatments increased, while the Shannon and Simpson indices decreased (Fig. 2C). The results showed an increase in the number of bacterial but a decrease in their diversity. Decline Chao1, Shannon, and Simpson indices for fungi indicate that fungi abundance and diversity reduced (Fig. 2D).
The results of principal component analysis (PCA) show that the samples treated are closely gathered (Fig. 3A and 3B). Taxa for different treatment compost groups were identified through random forest classification (Fig. 3C and 3D). It was found that WJP treatment increased abundance of Pedobacter, Gluconacetobacter and Anaerosporobacter. CaP treatment increased the abundance of Pedobacter, Actinomycetaceae, Anaerosporobacter, Delftia, and Gluconobacter. The abundance of Brevundimonas, Chryseobium, Azotobacter, Ketogulonicigenium, Weissella, Acetobacter, Sphingobacterium and Procabacter was reduced (Fig. 3C). For fungi, WJP treatment increased the abundance of had an increase in Pichia abundance, Dipodascus abundance increased when treated with CaP.
Differentially abundant core and specific taxa under WJP or CaP treatment
We also used DESeq software to screen differential abundant bacteria in the composting process based on a fold change of > 2 or < 0.5 and an adjusted P < 0.05(Fig. 4). In the CK21D vs WJP21D comparison group, 100 differentials bacterial ASVs, including 31 enriched and 69 depleted ASVs, were identified. These differential bacteria were distributed in three phyla: Bacteroidetes (22ASVs), Firmicutes (11 ASVs), and Proteobacteria (67ASVs). A total of 109 differential bacterial ASVs, including 43 enriched and 66 depleted ASVs, were identified in the CK21D vs. CaP21D comparison group. These differential bacteria were distributed in three phyla: Bacteroidetes (22 ASVs), Firmicutes (14 ASVs), and Proteobacteria (73ASVs). In the CK42D vs WJP42D comparison group, 121 differentials bacterial ASVs, including 58 enriched and 63 depleted ASVs, were identified. These differential bacteria were distributed in four phyla: Actinobacteria (5 ASVs), Bacteroidetes (31 ASVs), Firmicutes (28 ASVs), and Proteobacteria (57 ASVs). A total of 108 differential bacterial ASVs, including 43 enriched and 65 depleted ASVs, were identified in the CK42D vs. CaP42D comparison group. These differential bacteria were distributed in four phyla: Actinobacteria (1 ASV), Bacteroidetes (29 ASVs), Firmicutes (28 ASVs), and Proteobacteria (50 ASVs).
The species composition of the bacterial communities changed significantly as a result of a co-occurrence network analysis built at the ASV level (Fig. 5 and Table 1). In CK0D, a total of 134 nodes and 1002 edges, including 60.98% and 39.02% positive correlations of the ecological network, were obtained. In composting 21D, ecological networks of 145, 130, and 141 nodes with 859, 643, and 852 edges were obtained for the CK, WJP, and CaP treatment groups, respectively. In composting 42D, networks of 144, 108, and 116 nodes with 792, 516, and 561 edges were obtained for the CK, WJP, and CaP treatment groups, respectively.
Effects of WJP and CaP on microbiome functions in the composting process
PICRUSt2 software was used to determine the functional difference of microbiota in the composting process of WJP or CAP treatments. In Fig. 6, the top 10 KEGG pathways of all nine comparison groups are shown. In the CK0D vs. CK21D comparison group, aerobactin biosynthesis, coenzyme M biosynthesis I, and the superpathway of methylglyoxal degradation increased, while the superpathway of bacteriochlorophyll a biosynthesis, chlorophyllide a biosynthesis I (aerobic, light-dependent), and factor 420 biosynthesis were depleted. In the CK0D vs. CK42D comparison group, aerobactin biosynthesis, pyrimidine deoxyribonucleotides de novo biosynthesis IV, and pyrimidine deoxyribonucleotides biosynthesis from CTP increased, while chlorophyllide a biosynthesis I (aerobic, light-dependent), factor 420 biosynthesis, and vitamin E biosynthesis (tocopherols) were depleted. In the CK21D vs. WJP21D comparison group, the superpathway of demethylmenaquinol-6 biosynthesis II, chondroitin sulfate degradation I (bacterial), and the superpathway of bacteriochlorophyll a biosynthesis increased, while mycolyl-arabinogalactan-peptidoglycan complex biosynthesis, isoprene biosynthesis II (engineered), and coenzyme B biosynthesis were depleted. In the CK21D vs. CaP21D comparison group, the superpathway of bacteriochlorophyll a biosynthesis, chondroitin sulfate degradation I (bacterial), and D-cycloserine biosynthesis increased, while mycolyl-arabinogalactan-peptidoglycan complex biosynthesis, isoprene biosynthesis II (engineered), and coenzyme B biosynthesis were depleted. In the CK42D vs. WJP42D comparison group, p-cumate degradation, p-cymene degradation, and reductive acetyl coenzyme A pathway increased, while the superpathway of bacteriochlorophyll a biosynthesis, L-valine degradation I, and S-methyl-5-thio-α-D-ribose 1-phosphate degradation were depleted. In the CK42D vs. CaP42D comparison group, D-cycloserine biosynthesis, adenosine nucleotides degradation IV, and the reductive acetyl coenzyme A pathway increased, while sucrose degradation II (sucrose synthase), mycolyl-arabinogalactan-peptidoglycan complex biosynthesis, and pyrimidine deoxyribonucleotides de novo biosynthesis IV was depleted.
Effects of WJP and CaP on metabolites in the composting process
We used liquid chromatography-mass spectrometry (LC-MS) to identify and quantitatively analyze the compost metabolites. There were 4312 metabolites identified in total. Based on principal component analysis (PCA), the treatment and control groups were significantly different, and metabolomes in the same group we closely clustered (Fig. S4). We further identified the significantly different accumulated metabolites (DAMs) between different treatments (Fig. 7). A total of 90, 159, 139, and 181 DAMs were identified in the CK21D vs. WJP21D, CK21D vs. CaP21D, CK42D vs. WJP42D, and CK21D vs. CaP21D comparison groups, respectively (Fig. 7A). These DAMs were classified into 12 main classes (Fig. 7B).
Combined analysis of microbiomes and metabolomes
In order to elucidate the changes in specific microorganisms and metabolites during the composting process, the data from microbiomes and metabolomes were correlated. Firstly, the correlations between CK21D, WJP21D, CaP21D, CK42D, WJP42D, and CaP42D were determined, and we calculated the top 10 highest abundance bacterial genera in six group samples, as well as the 10 VIP metabolites with the highest difference (Fig. 8A). Actinomycetaceae, Flavobacterium and Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium were negatively correlated with metabolites. Sphingobacterium, Stenotrophomonas, Klebsiella, Lactobacillus, Comamonas, and Pseudomonas were positively correlated with metabolites. Second, we calculated the top 20 most relevant microbial ASVs and differential metabolites (Fig. 8B). In the CK21D vs. WJP21D comparison group, ASV_36821 (Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium) and compound_3970 (Styrene), compound_2346 (Hexanoic acid), and compound_2561(L-Proline) were positively correlated. ASV_11077 (Weissella) and compound_0206 ((5E,8Z)-4,7-dihydroxy-2-methyl-2,3,4,7-tetrahydrooxecin-10-one), compound_ 1907 (Diacetoxyscirpenol), and compound_0208 ((5E)-7-methylidene-10-oxo-4-(propan-2-yl)undec-5-enoic acid) were positively correlated. In the CK21D vs. CaP21D comparison group, ASV_21585 (Sphingobacterium), ASV_20551(Acetobacter), and ASV_26983 (Sphingobacterium) were positively correlated with most DAMs. In the CK42D vs. WJP42D comparison group, ASV_12156 (unclassified_Rhodobacteraceae), ASV_63724 (Azotobacter), ASV_55811 (Sphingobacterium), and ASV_ 59630 (Dysgonomonas) were positively correlated with compound_2726 (Mannitol), compound_1655 (Choline), compound_0650 (2',4'-Dihydroxy-3,4,6'-trimethoxydihydrochalcone), and compound_1654 (Choline O-Sulfate). In the CK42D vs. CaP42D comparison group, ASV_244 (Sphingobacterium), ASV_13381 (Dysgonomonas), and ASV_50383 (unclassified_Enterobacteriaceae) were positively correlated with most DAMs.
Third, redundancy analysis (RDA) was used to analyze the relationship between physical and chemical indexes and the microorganisms of compost samples (Fig. 9). Total N and P component had a positive correlation with most DAMs in Fig. 9. EC had a positive correlation with compound_3740 (Quebrachitol) and compound_1656(Choline), and a negative correlation with pH. Compound_0206 ((5E,8Z)-4,7-dihydroxy-2-methyl-2,3,4,7-tetrahydrooxecin-10-one), compound_2346 (Hexanoic acid), compound_2375 (Hydrocinnamic acid), and compound_2726 (Mannitol) were positively correlated with various phosphorus components.