Introduction: Deep eutectic solvents (DES) and natural deep eutectic solvents (NADES) have emerged as sustainable alternatives to conventional organic solvents for the extraction of phenolic compounds from aromatic plants. Multivariate chemometric tools provide a powerful means of evaluating complex extraction datasets and identifying the factors governing phenolic selectivity across different solvent systems.
Methods: Multivariate chemometric analysis was applied to evaluate the relationships among phenolic compound profiles of ten Lamiaceae species extracted with five solvent systems — two hydrophobic DES, two NADES, and 70% ethanol — yielding a dataset of 50 samples and 12 quantified phenolic compounds. Prior to analysis, data were standardized by z-score normalization. Principal Component Analysis (PCA), Hierarchical Cluster Analysis (HCA), and heatmap visualization were performed to assess sample grouping and identify key discriminating variables.
Results and Discussion: PCA revealed that the first two principal components accounted for 41.8% and 23.7% of total variance, respectively, with a cumulative explained variance of 65.5%. The score plot demonstrated a clear separation of samples primarily according to solvent polarity rather than botanical origin, indicating that extraction solvent is the dominant factor governing phenolic profiles in Lamiaceae species. The menthol/methyl salicylate DES formed a well-separated cluster driven by high rosmarinic acid, chlorogenic acid, and quercetin loadings, while NADES and ethanol extracts showed partial overlap, reflecting comparable polarity ranges. HCA confirmed two primary clusters corresponding to lipophilic DES and hydrophilic extraction systems. Loadings analysis identified rosmarinic acid, chlorogenic acid, and quercetin as principal discriminating variables.
Conclusions: These results demonstrate the utility of chemometric tools in characterizing solvent-dependent phenolic fingerprints and support rational solvent selection for targeted extraction from complex plant matrices.
Acknowledgments: This work is supported by Provincial Secretariat for Higher Education and Scientific Research, Province of Vojvodina (Grant No. 003874431 2025 09418 003 000 000 001).