The FID signals were processed in MestReNova software (Version 14.1.0; Mestrelab Research S.L., Santiago de Compostela, Spain) and the 1 H-NMR spectra were manually phased, baseline corrected using the Whitaker smoother algorithm and the TSP signal was referenced at 0.0 ppm. The pre-processed spectra were then exported as JCAMP files to Chenomx NMR Suite 8.5 software, Professional Edition (Chenomx Inc., Edmonton, AB, Canada), for annotation and quantification, identifying a total of 35 compounds (Supplementary Table ). Concentrations of compounds not identified in the spectra were replaced with 1/3 of the lowest recorded concentration among all the spectra.