“…For the HCA, normalized signal intensities of metabolites (unit variance scaling) were represented as a color spectrum [ 34 ]. Significantly regulated metabolites between groups were identified based on variable importance in projection (VIP) values ≥ 1 and absolute Log2FC ≥ 1 [ 34 ]. VIP values were extracted from the orthogonal partial least squares-discriminant analysis (OPLS-DA) result, which also included score plots and permutation plots, and was generated using the R package MetaboAnalystR [ 34 ].…”