1. Introduction
Phlomoides rotata (Benth. ex Hook. f.) Mathiesen (PR), also known as
Lamiophlomis rotata (Benth.) Kudô, is a medicinal herb called "Duyiwei" (
Lamiophlomis herba) in Chinese. It belongs to the
Phlomoides Moench of Lamiaceae and is found in high altitudes in China [
1,
2,
3,
4]. Traditionally, the root, rhizome, or the entire herb were used for medicinal purposes [
5,
6]. However, nowadays, only the aerial parts are used [
4], and the digging of the root is prohibited due to PR being listed as a first-class endangered Tibetan medicine [
6,
7]. The underground parts of PR are known for promoting blood circulation, eliminating stasis, reducing swelling, and providing analgesic effects. The above-ground parts are used for treating grasserie, fractures, injuries from falls, osteomyelitis, gunshot injuries, and edema pain [
2,
5,
6,
7].
Lamiophlomis herba is yellowish-brown or sallow, with a bitter and flat taste [
4,
5,
7]. It was first documented in the classical masterpiece of Tibetan Medicine, Somaratsa [
7], and has been used for more than 2000 years in traditional Tibetan medicine, known as "Daba" and "Dabuba," to treat traumatic injuries, rheumatic arthritis (RA), and grasserie [
4,
5,
6,
7]. In clinical practice, PR is typically used directly without any prior processing, mainly for pain relief [
6]. Additionally, PR is commonly used as a key ingredient in combination with other Chinese herbs such as
Curcuma longa,
Salvia miltiorrhiza, and
Pyrrosia lingua [
6,
7]. PR is also utilized as an ingredient in various health products, including health drinks, soap, wine, mouth rinses, and biological toothpastes [
6].
Due to the low content of volatile chemicals, the focus of research on PR has primarily been on its involatile compounds. So far, at least 223 chemical constituents have been isolated from PR, including iridoids, flavonoids, phenylethanoid glycosides, polysaccharides, organic acids, volatile oils, and others [
6,
7]. The main compounds found are iridoid glycosides, which are responsible for the analgesic effect [
6,
7] and serve as characteristic metabolites (CMs) to assess the quality of Duyiwei [
7]. However, C13-norisoprenoids have been largely overlooked, with only one previously reported compound, 5
β, 6
α-dihydroxy-3
β-(
β-D-glucopyranosyloxy)-7-megastigmen-9-one [
8]. Additionally, fourteen organic acids, including palmitic acid (PA), have been isolated and identified [
6]. As of now, there is only one published study that has reported the chemical composition of the essential oils (Es) extracted from PR using steam distillation, and another study has reported the lipophilic composition in the CH
2Cl
2 extracted part. Additional details and data can be found in
supplemental Table S1 [
9,
10]. The main compounds found in terms of content are fatty acids (FAs), especially long-chain FAs such as PA, myristic acid (MA), oleic acid (OA), and linoleic acid (LA). However, some identifications are still debatable, such as OA, linoleic acid ethyl ester, and cyclohexenylacetic acid, taking into consideration the values of linear retention indices (LRIs) and the mass spectra [
9,
10]. Only three compounds, MA, PA, and LA, were detected in both studies [
9,
10]. The identification of CMs from PR EOs is of great importance as these CMs play a crucial role in distinguishing these EOs from others rapidly [
11,
12,
13].
In recent years, there have been limited studies on the in-depth biological effects of volatile oils from PR [
7]. However, the petroleum ether extracted part has been reported to possess anti-tumor activities, suggesting that the volatile compounds may also have such activities [
14]. Currently, there have been no tests conducted on the antioxidant activities (AOAs) of volatile chemicals extracted from PR. Oxidative stress has been associated with various diseases, including RA, cancer, and diabetes [
15,
16]. Additionally, supplemental FAs play an important role in maintaining the balance between oxidation and antioxidation in cells [
17,
18,
19,
20,
21,
22,
23,
24]. The effects of FAs on oxidant injury seem to be linked to the degree of unsaturation [
18]. PA, a type of saturated FAs (SFAs), can increase oxidative stress in cells in a concentration-dependent manner [
23,
24]. This is because it can induce overexpression of the pro-oxidant protein p66Shc [
23] or react with cells to generate reactive oxygen species (ROS), reduce the content of NO, and make cells more susceptible to oxidative stress [
24]. On the other hand, stearic acid, another type of SFAs, has been reported to protect pulmonary artery endothelial cells from oxidant injury [
18]. Generally, polyunsaturated FAs (PUFAs) can reduce oxidant injury [
19,
20,
21,
22], although contradictory results have also been reported [
18].
To the best of our knowledge, no single extract or compound from PR has been clinically applied for disease treatment thus far. Therefore, it is necessary to conduct studies and develop potentially therapeutic extracts or compounds [
7]. Building upon previous research [
25], our focus in this study is on the volatile chemicals and their AOAs. We have selected three samples of PR from Tibet to evaluate the CMs present in the Es and their AOAs. The process involved isolating crystals (Cs) from the Es using cryoprecipitation, which resulted in crystal-free Es (CFs). Subsequently, a comprehensive chemical profiling study was conducted on the Es, Cs, and CFs using GC-MS and GC-FID techniques. Additionally,
in vitro assays have been conducted to assess the AOAs of the Es, separated parts, and four CMs, namely PA, MA, LA, and OA. The findings of this study will contribute to establishing a theoretical foundation for the utilization of PR Es.
2. Results
2.1. Extraction and Separation
From PR with Voucher No. L8, L9, and L10, a total of 0.29 g, 0.26 g, and 0.19 g of light yellow Es were obtained, corresponding to 418 μL, 405 μL, and 238 μL, respectively. The yields of the Es were calculated as 0.13%, 0.13%, and 0.08%, respectively, based on the ratio of the volume of the Es (in mL) to the weight of the PR (in g) (mL/g). The densities of the Es were measured as 0.69, 0.64, and 0.8, respectively. These Es have a fresh and elegant smell. Crystals were separated from the EOs at temperatures of 4 °C or -4 °C. The EOs were then obtained by removing the crystals.
2.2. Chemical Composition of the Es from PR
A total of 107 compounds have been both qualified and quantified. Among them, 31 compounds have been previously reported [
9,
10], while 76 compounds are being reported for the first time from the EOs of PR. These compounds belong to various classes, including FAs, esters of FAs, C13-norisoprenoids, AMs, C18-norditerpenoids, ADs, Als, Others, and Introduced. The relative contents of these compounds are listed in
Figure 2.
Due to the low concentrations of samples, only eight compounds were detected by GC-MS using DB-5. These compounds include hexanal, 1-octen-3-ol, limonene, linalool, α-terpineol, hexahydrofarnesyl acetone, methyl palmitate, and PA. Among them, the contents of limonene, α-terpineol, and PA are relatively high.
The 9-hexadecenoic acid reported previously [
9] is most likely corresponding to the 9
E-hexadecenoic acid detected in this study, based on their LRIs values.
The EOs, crystals, and EOs free of crystals primarily consist of FAs, particularly LCFAs. PA stands out as the most prominent FA, which aligns with the reported findings [
9,
10]. Additionally, MA, OA, and LA are also significant, as previously reported [
9,
10]. The content of PA is relatively higher in crystals but relatively lower in EOs free of crystals compared to the corresponding EOs.
Regarding the esters of FAs, the major compounds are methyl palmitate, methyl oleate, methyl linoleate, and methyl linolenate.
These seven types of C13-norisoprenoids are being reported for the first time from the EOs of PR. Their content is relatively low, and the content in the EOs is relatively higher compared to that in the corresponding crystals or EOs free of crystals.
Among the AMs, linalool is relatively prominent. Within the C18-norditerpenoid group, hexahydrofarnesyl acetone stands out. Phytol is a noteworthy ADs. Tricosane and pentacosane are two notable Als. In addition, compounds such as hexanal, 1-octen-3-ol, and
β-caryophyllene should also be observed. A total of eight compounds, including limonene,
α-terpineol,
β-pinene,
β-myrcene,
p-cymene,
γ-terpinene, isoterpinolene, and terpinen-4-ol, are likely introduced from the EOs extracted from the peels of Nanfengmiju (
Citrus kinokuni Hort. ex Tanaka), a variety of
C. reticulata, which were studied simultaneously [
11].
2.3. Principal Component Analysis (PCA) and Partial Least Squares Discriminant Analysis (PLS-DA) of the Metabolites in Es of PR
Figure 3.
A total of 50 compounds were screened based on the data obtained by GC-FID and GC-MS with a FFAP column. The screening standard required a VIP (Variable Importance in Projection) value of not less than 1.
Figure 3.
A total of 50 compounds were screened based on the data obtained by GC-FID and GC-MS with a FFAP column. The screening standard required a VIP (Variable Importance in Projection) value of not less than 1.
Figure 4.
PCA and PLS-DA plots. In the sample names, the "_1" and "_2" indicate that the data was obtained from GC-FID and GC-MS with a FFAP column, respectively. The compounds represented in the loadings of the PCA plot are as follows: EF10 - Methyl linolenate, F13 - PA, O28 - 1,8(2H,5H)-Naphthalenedione, hexahydro-8a-methyl-, cis-, EF1 - Nonanoic acid, 9-oxo-, methyl ester, A5 - Hexadecane, F8 - Tridecanoic acid, O6 - Hexanal, F4 - Nonanoic acid, O23 - Cedrol, I3 - p-Cymene, A2 - Farnesane, I5 - γ-Terpinene.
Figure 4.
PCA and PLS-DA plots. In the sample names, the "_1" and "_2" indicate that the data was obtained from GC-FID and GC-MS with a FFAP column, respectively. The compounds represented in the loadings of the PCA plot are as follows: EF10 - Methyl linolenate, F13 - PA, O28 - 1,8(2H,5H)-Naphthalenedione, hexahydro-8a-methyl-, cis-, EF1 - Nonanoic acid, 9-oxo-, methyl ester, A5 - Hexadecane, F8 - Tridecanoic acid, O6 - Hexanal, F4 - Nonanoic acid, O23 - Cedrol, I3 - p-Cymene, A2 - Farnesane, I5 - γ-Terpinene.
It is evident that within each group, the six samples are relatively closer to each other, while between the groups, they are relatively less close.
2.4. Metabolic Pathway Analysis
As shown in
Figure 5, the metabolite pathway listed in sequence are fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, cutin, suberine and was biosynthesis, fatty acid elongation, fatty acid degradation.
2.5. AOAs of Es, Cs, RCs, and Four CMs
In the DPPH and ABTS assays, only the IC50 of vitamin C (Vc) was directly detected, while the IC50 of other samples was deduced from the results obtained (
Supplemental Table S3).
In the DPPH assay (
Figure 6), among the nine samples, RC10 exhibits the highest radical scavenging activity (RSA) at 14.7% in a concentration of 110 μg·mL
-1. PA demonstrates some AOAs at a concentration of 100 μmol·L
-1, but its RSA values become negative as the concentration increases. The RSA values of OA and LA in the high concentration group show a positive correlation with concentration, although the initial value of OA is negative. At 1200 mmol·L
-1, the RSA value of OA is 6.39%, whereas that of LA is 35%. Vc demonstrates the lowest IC50 value, indicating the strongest antioxidant effect.
The ABTS assay (
Figure 7) reveals that the RSA values of RCs are generally higher than those of the corresponding Es or Cs, indicating that the Cs may contain compounds with antagonistic AOA. Interestingly, the IC50 value of Vc is higher in the ABTS assay compared to the DPPH method, while the IC50 value of QE is lower in the ABTS assay. It is important to note that most of the samples exhibit improved RSA values when compared to the results obtained from the DPPH assay. This improvement can be attributed to the higher reactivity of ABTS radical cations [
27].
In the FRAP assay, the AOAs are measured based on the concentration of Fe2+ (mmol·L-1). At a concentration of 300 mmol·L-1, the FRAP values of OA and LA are significantly higher than that of Vc at 15 μmol·L-1, and similar to that of QE at 200 and 100 μmol·L-1, respectively. However, as the concentration increases to 600 and 1200 mmol·L-1, the FRAP values of OA and LA gradually decrease. It is noteworthy that the FRAP values of nine samples are similar to that of Vc at 10 μmol·L-1.
Overall, OA and LA exhibit some pro-oxidation activities (POAs) in the lower concentration group, but demonstrate good AOAs in the higher concentration group.
4. Materials and Methods
4.1. Plant Materials, Reagents and Chemicals
The aboveground portion of PR was collected from three populations: L8, L9, and L10. These populations corresponded to the same sample numbers used in a previous research study [
25]. The authenticity of the collected populations was confirmed by Professor Yi Zhang, from the Chengdu University of Traditional Chinese Medicine (CUTCM), Chengdu, China, and using internal transcribed spacer 2 (ITS2) DNA barcodes, as described in the previous study [
25]. Voucher samples of L8, L9, and L10 were deposited in the College of Ethnic Medicine at CUTCM, Chengdu, China, and the Chongqing Academy of Chinese Materia Medica, Chongqing, China.
Table 1.
The origins of PR and the corresponding GenBank accession numbers of the ITS2 sequences, as mentioned in the previous study [
25].
Table 1.
The origins of PR and the corresponding GenBank accession numbers of the ITS2 sequences, as mentioned in the previous study [
25].
Voucher No. |
Sources |
GPS Coordinates |
GenBank Accession Number |
L8 |
BianBa, LeiWuQi and NaQu counties of Tibet |
E: 93° W: 31° |
KP699743/45-4750-51/54 |
L9 |
L10 |
The following reagents and chemicals were used in this study: n-hexane (high-performance liquid chromatography grade), linalool (98%+), p-cymene (99%+), α-terpineol, and nonane (98%) (98%+), which were produced by Adamas Reagent Company Ltd.; d-limonene (96%) produced by Acros Organics, USA; γ-terpinene (97%) produced by Wako Pure Chemical Industries, Ltd., Japan; PA produced by CATO; n-alkanes standard solution of C10-C25 produced by Dr. Ehrenstorfer Inc, Germany; n-octacosane (99%) produced by Aldrich; DPPH, ABTS powder, potassium persulfate (K2S2O8), and Vc. All these reagents and chemicals were supplied by Shanghai Titan Scientific Co.,Ltd., China.
4.2. Extraction and Separation
Weighed powders of L8, L9, and L10 (315 g each) were placed in separate round-bottomed flasks. Pure water (3150 mL, 10 volumes) was added to each flask. The mixtures were soaked for 0.5 hours at 40 °C. The Es were extracted three times from each powder using hydrodistillation with a Clevenger-type apparatus. Each extraction lasted for 5 hours. n-Hexane was employed as the collecting solvent. The collected Es, which appeared as light yellow, were treated with anhydrous Na2SO4 to eliminate any residual water.
To evaluate crystallization, the Es of L8, L9, and L10 were stored at different temperatures: 4 °C, -4 °C, and -80 °C. Cs were obtained either at 4 °C or -4 °C. However, at -80 °C, the RCs of all three samples were in a solid state. Therefore, there were three sets of samples for L8, L9, and L10: E, C, and RC. Specifically, there were nine samples: E8, E9, and E10, C8, C9, and C10, RC8, RC9, and RC10. Each sample was stored in separate screw-capped vials at 4 °C.
4.3. The Identification and Quantitation of Chemicals in the Es, Cs, and RCs
4.3.1. Sample Preparation
The samples of E8, E9, E10, C8, C9, C10, RC8, RC9, and RC10 were prepared for analysis as follows: For GC-FID and GC-MS detection using a DB-5 column (30 m×0.25 mm i.d., 0.25 μm film thickness), the samples were diluted in the ratio of Vsample: Vn-hexane (HPLC) 1:1000 (0.1%). For GC-MS detection using a FFAP column (30 m×0.32 mm×0.5 μm), the samples were diluted in the ratio Vsample: Vn-hexane (HPLC) 1:250 (0.4%).
4.3.2. GC Analyses
GC analyses were performed using the prepared samples. During the analyses, the retention times and peak areas of the compounds were recorded. The obtained chromatograms and mass spectra were analyzed to identify and quantify the target compounds present in the samples.
GC-FID analyses were conducted using a GC-2010 instrument from Shimadzu, Japan, equipped with a DB-5 column. The following parameters were used: The oven temperature was programmed from 60 °C with a 3-minute hold, and then ramped up to 250 °C at a rate of 2.5 °C per minute. The final temperature was held for 2 minutes. Nitrogen was used as the carrier gas, with a constant flow rate of 1.7 mL per minute. Both the injector and the detector were maintained at 250 °C. The sample was split in a ratio of 5:1. Each sample was injected once with a volume of 1 μL.
GC-MS analyses were performed using a GCMS-TQ8040 instrument from Shimadzu, Japan. Separately, either a DB-5 column or a FFAP column was used for the GC-MS analyses. The following parameters were used.
For DB-5 column analysis: The oven temperature was programmed from 60 °C with a 3-minute hold, and then ramped up to 280 °C at a rate of 2.5 °C per minute. The final temperature was held for 2 minutes. Helium was used as the carrier gas, with a constant flow rate of 1 mL per minute. The sample was split in a ratio of 100:1. A solvent delay of 3.0 minutes was implemented. The injector was maintained at 250 °C, while the ion-source and interface were maintained at 200 °C and 250 °C, respectively. Mass spectra were acquired at 70 eV with a scan rate of 3.9 scans per second, covering the mass range from m/z 25 to 450 amu. Each sample was injected once with a volume of 1 μL.
For FFAP column analysis: The oven temperature was programmed from 60 °C with a 3-minute hold, and then ramped up to 230 °C at a rate of 2.5 °C per minute. The final temperature was held for 2 minutes. The rest of the parameters (carrier gas, splitting ratio, solvent delay, injector, ion-source, interface temperature, mass spectra, injection times, and injection volume) remained the same as for the DB-5 column analysis.
4.3.3. Identification and Quantitation
4.3.3.1. Identification
For compound identification, either a NIST 14 or a NIST 17 MS database was utilized. Initially, the peaks in the TICs were identified using probability-based matching. However, due to the presence of overlapped and embedded peaks in the TICs, the identification results may sometimes be inaccurate. In such cases, characteristic ion peaks were selected and compared with the NIST 14 or 17 database, as well as the mass spectra of known standards. By employing a combination of probability-based matching and comparing characteristic ion peaks, the compound identification process becomes more reliable, particularly when dealing with complex samples with overlapping peaks. This approach enhances the overall accuracy and confidence in the identification results.
The LRI value was calculated using the equation (1) proposed by Van Den Dool and Kratz [
54]. In equation (1), “t
n” and “t
n+1” represent the retention time (RT) of the
n-alkanes (C10-C25, C28, and the detected C26-C27, C29) with the corresponding number (“n”) of carbons. "t
x" represents the RT of the detected compound x (t
x), where “t
n ≤ t
x ≤ t
n+1”. Furthermore, the RT of
n-alkane C30 (t
30) was deduced by analyzing the TIC obtained through GC-MS using the FFAP column.
The calculated LRIsb1, LRIsb2, and LRIsd, which correspond to the RTs detected in GC-FID using a DB-5 column, GC-MS using a DB-5 column, and GC-MS using a FFAP column, respectively, were compared with the LRIa and LRIc recorded in the NIST 17 MS database for the corresponding identified chemicals.
4.3.3.2. Quantification
The relative area percentage of each compound was calculated using peak area normalization.
4.4. PCA and PLS-DA
In R version 4.3.1, the plots of VIP, PCA, PLS-DA for the nine samples and chemicals were generated. Variables with a VIP score greater than 1 are considered significant. The selection of compounds was based on the criteria that the VIP score > 1 and a significance level of P ≤ 0.05.
4.5. Metabolic Pathway Analysis
The metabolic pathway analysis of the CMs was performed using MetaboAnalyst 6.0, a web-based tool (
https://www.metaboanalyst.ca/), as well as the KEGG system (
https://www.kegg.jp/kegg/). Specifically, the results of the metabolic pathway analysis were visualized in “
Figure 5”, which was generated using MetaboAnalyst 6.0.
4.6. AOAs of Es, Cs, RCs, and Four CMs
The AOAs of nine samples, E8, E9, E10, C8, C9, C10, CF8, CF9, and CF10, were tested by DPPH, ABTS, and FRAP assays, respectively.. In addition, four CMs including PA, MA, OA, and LA, as well as two reference compounds, Vc and QE, were also included in the analysis.
4.6.1. Sample Preparation
All the subjects were diluted in methanol (MeOH). Since there was a limited amount of volatile oils available, only three low concentrations, namely 50, 80, and 110 μg·mL-1, were set up for the nine samples of Es, Cs, and RCs. The four CMs (PA, MA, OA, LA) and the reference compound QE were diluted to specific concentrations of 100, 200, and 400 μmol·L-1, respectively. In the case of OA and LA, an additional higher concentration group was established, with concentrations of 300, 600, and 1200 mmol·L-1. Furthermore, the reference compound Vc was diluted to concentrations of 5, 10, and 15 μmol·L-1 for the analysis.
4.6.2. DPPH Assay
A slight modification was made to the experimental procedure [
55]. A volume of 100 μL of the sample at various concentrations was added to individual wells of a 96-well microplate. Subsequently, 100 μL of DPPH solution (100 μmol·L
-1) also diluted with MeOH was added to each well. The microplate was then incubated in darkness at room temperature for 30 minutes. After the incubation period, the absorbance of the reaction mixture was measured at 517 nm using a microplate reader. Each sample was analyzed in triplicate, and MeOH was used as the blank control. The RSA was calculated using the following equation:
Where ABlank is the absorbance of the blank control (MeOH) and ASample is the absorbance of the reaction mixture. The RSA value indicates the percentage of DPPH radical scavenged by the sample, with higher values indicating stronger antioxidant activity.
4.6.3. ABTS Assay
A slight modification was made to the method described in a previous study [
56]. To prepare the ABTS radical cation (ABTS
•+) solution, 5 mL of a 7 mM aqueous ABTS solution was reacted with 88 µL of a 140 mM K
2S
2O
8 aqueous solution (resulting in a final concentration of 2.45 mM for K
2S
2O
8). The reaction mixture was kept in darkness at room temperature for 16 hours. Subsequently, the radical cation solution was diluted with MeOH, typically around 30-50 times, until its absorbance reached a value of 0.7±0.02 at 734 nm. For each sample, 100 μL was added to 100 μL of the ABTS radical solution. The mixture was thoroughly mixed at room temperature for 6 minutes. Following the incubation period, the absorbance at 734 nm was measured using a microplate reader. The calculation method for RSA was consistent with that used in the DPPH assay.
4.6.4. FRAP Assay
A slight modification was made to the method described previously [
53]. For each sample, 100 μL was added to 100 μL of the FRAP working solution. The FRAP working solution was composed of acetic acid buffer (0.3 mol·L
-1), TPTZ (2, 4, 6-Tris (2-pyridyl)-1, 3, 5-triazine) solution (10 mM), and FeCl
3 (20 mM) solution, with a volume ratio of 10:1:1. The sample and the FRAP working solution were thoroughly mixed and then placed in darkness at 37 °C for a 30-minute incubation period. After incubation, the absorbance of the mixture at 593 nm was immediately measured using a microplate reader. The increase in absorbance at this wavelength indicates the reduction capacity of the samples and provides information on their antioxidant potential.
To establish a calibration curve, 0.1 mL of Fe(II) aqueous solutions with concentrations ranging from 0.01 to 0.2 mM were mixed with 0.1 mL of the FRAP reagent. In this measuring system, the total antioxidant capacity was determined in terms of Fe(II) equivalents. The concentration of FeSO4 (in mmol·L-1) was calculated based on the absorbance value obtained from the standard curve after the reaction, which was referred to as the FRAP value. A higher FRAP value indicates a stronger antioxidant activity.
5. Conclusions
A comprehensive analysis of the Es of PR has revealed a total of 107 identified and quantified compounds. Among these compounds, 76 have been reported for the first time from the Es of PR. Based on their content, seven compounds, including PA, MA, LA, OA, methyl palmitate, hexahydrofarnesyl acetone, and phytol, have been selected as CMs. Additionally, trans-β-damascenone has been chosen as a CM based on its distinct flavor profile. In summary, these eight compounds and hexanal have been identified as the CMs in the Es of PR. Metabolic pathways related to these CMs include fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, cutin, suberine and wax biosynthesis, fatty acid elongation, and biosynthesis of terpenoids. The interplay and relationships between these metabolic pathways have also been depicted.
Interestingly, RCs have demonstrated stronger AOAs, while Cs have shown weaker AOAs compared to Es in general. This disparity is likely influenced by the varying content of PA in these samples. Specifically, the percentage of PA, the most abundant compound, has been found to be higher in Cs but lower in RCs compared to Es.
The accuracy of chemical identification and quantification was somewhat affected due to the inappropriate dilution of samples. Furthermore, no synergistic AOAs were detected in the screened CMs. In future research, we should exercise caution during sample dilution and to place a stronger emphasis on investigating the synergistic effects of CMs, particularly focusing on samples that exhibit stronger AOAs.
This study significantly contributes to our understanding of the composition of the Es of PR, shedding light on the potential applications and utilization of these volatile oils, which are rich in LCFAs.
Author Contributions
Conceptualization, Z.P. and J.W.; methodology, C.X. and J.W.; software, C.X., J.L., X.Y., A.H. and J.W.; validation, Z.P., C.X., X.Y., Y.S., A.H. and J.W.; investigation, C.X. and J.W.; resources, Z.P.; data curation, C.X., X.Y. and J.W.; writing—original draft preparation, C.X., J.L., X.Y. and J.W.; writing—review and editing, Z.P., C.X., X.Y., Y.S., A.H., N.M.K., S.A. and J.W.; supervision, Z.P. and J.W.; project administration, Z.P. and J.W.; funding acquisition, Z.P. and J.W. All authors have read and agreed to the published version of the manuscript.