Bubble plots of the top 20 enriched KEGG pathways from the differentially expressed metabolites between the AMW and Control groups on days 60. In E-G, Number indicates the number of DE Metas which were annotated to the corresponding pathway: a bigger value indicates more DE Metas enriched in this pathway; Ratio indicates the ratio of the number of DE Metas which were enriched in the corresponding pathway to the number of annotated metabolites in this pathway: a bigger value indicates a higher enrichment level of DE Metas in this pathway; and p -values indicates the significance calculated by hypergeometric inspection. The color of symbols indicates their corresponding p -value: a smaller value indicates greater inspection reliability. VIP, variable importance in the projection; DE Meta, differentially expressed metabolites