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).