Statistical analysis
Regression analyses were performed using models which best fit the
observed data considering the highest R2. Derivatives
of the regression equations were calculated to assess relativeg s change after exogenous application of
different compounds and maximum velocity (Vmax )
of g s responses. The g srate of change was plotted versus time and undergone clustering
classification for time-series comparison analysis using dtwclustpackage in R (Sardá-Espinosa 2019). The metabolomics data were analysed
by orthogonal partial least squares discriminant analysis (orthoPLS-DA)
and principal component analysis (PCA) using FactoMineR package (Lêet al. 2008). The heat maps and the hierarchical clustering
analysis (HCA) were carried out using the MetaboanalystR package on R
(Pang et al. 2020; R Core Team 2020).