The names of the metabolites and their KEGG IDs were matched using MetaboAnalyst (http://www.metaboanalyst.ca/) . Pathway enrichment analysis was performed using MBRole . The relevant pathways were acquired by inputting metabolite KEGG ID and running a Hypergeometric Test. The metabolic pathway illustration was generated by Cytoscape and MetScape .