Next Article in Journal
Seeding, Plating and Electrical Characterization of Gold Nanowires Formed on Self-Assembled DNA Nanotubes
Next Article in Special Issue
Identification of a Sesquiterpene Lactone from Arctium lappa Leaves with Antioxidant Activity in Primary Human Muscle Cells
Previous Article in Journal
Supercritical Extraction of Red Propolis: Operational Conditions and Chemical Characterization
Previous Article in Special Issue
Preparation, Characterization and Antioxidant Activities of Kelp Phlorotannin Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Enzymatic, Ultrasound, and Reflux Extraction Pretreatments on the Yield and Chemical Composition of Essential Oils

1
Faculty of Food Technology and Biotechnology, University of Zagreb, Pierottijeva 6, 10 000 Zagreb, Croatia
2
Faculty of Chemistry and Technology, University of Split, Ruđera Boškovića 35, 21 000 Split, Croatia
3
NMR Centre, Ruđer Bošković Institute, Bijenička cesta 54, 10 000 Zagreb, Croatia
4
Department of Natural and Health Sciences, Juraj Dobrila University of Pula, Zagrebačka 30, 52 100 Pula, Croatia
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(20), 4818; https://doi.org/10.3390/molecules25204818
Submission received: 8 September 2020 / Revised: 13 October 2020 / Accepted: 19 October 2020 / Published: 20 October 2020
(This article belongs to the Special Issue Antioxidants from Natural Sources: Separation and Characterization)

Abstract

:
The effect of different hydrodistillation pretreatments, namely, reflux extraction, reflux extraction with the addition of cell wall-degrading enzymes, and ultrasound, on the yield and chemical composition of essential oils of sage, bay laurel, and rosemary was examined. All pretreatments improved essential oil yield compared to no-pretreatment control (40–64% yield increase), while the oil quality remained mostly unchanged (as shown by statistical analysis of GC-MS results). However, enzyme-assisted reflux extraction pretreatment did not significantly outperform reflux extraction (no-enzyme control), suggesting that the observed yield increase was mostly a consequence of reflux extraction and enzymatic activity had only a minute effect. Thus, we show that ultrasound and reflux extraction pretreatments are beneficial in the production of essential oils of selected Mediterranean plants, but the application of enzymes has to be carefully re-evaluated.

1. Introduction

Bioactive compounds are secondary plant metabolites with potential for many applications in human use. Due to the high content of bioactive compounds in essential oils [1,2,3] and their ability to stop or delay the aerobic oxidation of organic matter, their application in the food industry, cosmetic industry, phytotherapy, medicine, and many other fields is constantly growing [4,5,6,7,8].
Essential oils are usually isolated by steam- or hydrodistillation. Although these methods have been successfully used for years, novel extraction procedures are being developed in recent decades, with the aim to increase the yield of essential oils and/or improve their chemical composition, while shortening the extraction time [9]. To get closer to that ideal goal, novel production procedures, like ultrasound or enzymatic degradation of the cell wall, are being applied as pretreatments to distillation. Benefits of ultrasound application in solid–liquid extraction are the intensification of mass transfer, improved solvent penetration into the plant tissue, and capillary effects. The collapse of cavitation bubbles near the cell walls is assumed to cause the cell disruption and at the same time a good penetration of the solvent into the cells through the ultrasonic jet [10]. Several studies have reported that the addition of ultrasonic pretreatment before hydrodistillation leads to a shortening of the extraction time [11,12,13,14], sometimes also accompanied by yield increase [12,13,15,16,17] and improvement in oil chemical composition [9,12,13,14].
Further, enzyme-assisted extraction has been intensively studied in the last decade, since the plant cell wall, as a structure composed mainly of resistant polymers like cellulose, xylan, lignin, and pectin, can reduce the extraction efficiency. Enzymes, like cellulase, xylanase, and pectinase, can degrade or disrupt cell wall components, thus enabling better release and more efficient extraction of bioactive compounds from plants and improving the bioactive content of essential oils and extracts [2,18,19,20,21,22]. However, while some authors reported that enzyme-assisted extraction enhanced essential oil yield and improved their chemical composition [2,19,21,22,23,24,25,26], others did not demonstrate such effect [1], indicating that further research of enzyme application in plant extractions is necessary.
We hypothesized that different hydrodistillation pretreatments, namely, ultrasound-assisted extraction, reflux extraction, and enzyme-assisted reflux extraction, can improve the yield of essential oils of selected Mediterranean wild plants: sage (Salvia officinalis L.), rosemary (Rosmarinus officinalis L.), and bay laurel (Laurus nobilis L.). Also, we hypothesized that these relatively mild pretreatments will not adversely affect the chemical composition of essential oils.

2. Results

2.1. Cellulase, Xylanase, and Pectinase Can Degrade Their Respective Substrates under the Reaction Conditions

Before applying cell wall-degrading enzymes (xylanase, pectinase, and cellulase) for the extraction of bioactive compounds from the plant material, their activity under the extraction conditions (1 h, 40 °C, in MiliQ water) was measured at a small scale by 3,5-dinitrosalicylic acid (DNSA) assay. The enzymes were mixed with their appropriate polymeric substrates in the MiliQ water and the reaction was incubated for 1 h at 40 °C. The results have shown that all enzymes were able to degrade their substrates: xylanase degraded xylan to 1.10 U/mL, pectinase pectin to 0.03 U/mL, and cellulase cellulose to 0.05 U/mL.

2.2. Optimization of Ultrasound Extraction Conditions

Optimal ultrasonic power and duration, chosen based on the overall maximum total phenol yield, were 30% of maximal ultrasonic power during 10 min (Table 1). In some cases, total phenol yield was slightly higher after 15 min ultrasound treatment (e.g. at 30% of maximal ultrasonic power for rosemary, and 90% for bay laurel). Still, since 15 min ultrasound treatment caused intensive solvent evaporation, the shorter time of 10 min was chosen as optimal.

2.3. Different Pretreatments Increased the Essential Oil Yield

All tested hydrodistillation pretreatments caused an increase in essential oil yield when compared to no-pretreatment control (HD): 56–64% for bay laurel, 40–60% for rosemary, and 45–55% for sage (Figure 1). However, only slight differences in essential oil yield were observed between different pretreatments for the same plant species. Importantly, none of the individual (HD-REX, HD-REC, HD-REP) or combined (HD-REPCX) enzymatic pretreatments significantly outperformed reflux extraction pretreatment (HD-RE). Since HD-RE served as a no-enzyme control, this indicates that enzymatic activity had only a minute effect on essential oil yield increase and that it should be mostly attributed to reflux extraction per se.
When comparing different plants, Clevenger hydrodistillation of bay laurel leaves consistently resulted in the highest oil yield (from 1.95 mL oil/100 g dry plant for HD-REP to 2.05 mL oil/100 g dry plant for HD-REX and HD-REPCX), followed by sage (from 1.45 mL oil/100 g dry plant for HD-RE to 1.55 mL oil/100 g dry plant for HD-REX and HD-REPCX) and rosemary (from 0.7 mL oil/100 g dry plant for HD-REC to 0.8 mL oil/100 g dry plant for HD-RE).

2.4. Overall Chemical Composition of Essential Oils Was Not Significantly Affected by Different Hydrodistillation Pretreatments

The chemical composition of obtained essential oils was analyzed by GC-MS (Tables S1–S3, Figure 2) and NMR (Supplementary file 1, Figures S1–S19). The pretreatments did not significantly affect the quality of essential oils (p < 0.05), that is, the oil composition was comparable to the no-pretreatment control, as confirmed by Spearman’s test performed using GC-MS data (Supplementary file 2, Figures S20–S22). GC-MS analysis showed that the essential oils of all plants were richest in monoterpenes: up to 74%, 71%, and 62% of total peak area in sage, bay laurel, and rosemary essential oils, respectively (Tables S1–S3, Figure 2). Sesquiterpenes were also well represented: up to 35%, 17%, and 56% of total peak area in sage, bay laurel, and rosemary essential oils, respectively. Dominant compounds (i.e., terpenes represented with more than 5% of total peak area) were confirmed by NMR (as described and represented in Supplementary file 1, Figures S1–S19 and Appendix 1): linalool, camphor, borneol, and berbenone in rosemary essential oil (Figure S19A), camphor, manool, α- and β-thujone, and veridiflorol in sage essential oil (Figure S19B), and α-terpenyl acetate, linalool, and methyleugenol in bay laurel essential oil (Figure S19C).
Seventy-one components were identified in the sage essential oils by GC-MS, representing up to 97% of total GC peak areas (Table S1). Sage essential oils were richest in oxygenated monoterpenes, such as α- and β–thujone (up to 23% and 10%, respectively), camphor (up to 17%), borneol (up to 8%), and 1,8-cineole (up to 8%). Additionally, some sesquiterpenes were present in significant quantities, such as manool and veridiflorol (up to 14% each). Among other compounds, phenylpropane derivates, eugenol (up to 3%), and methyleugenol (up to 2%) were detected with significant percentages, although much lower than monoterpene and sesquiterpene major compounds.
In the bay laurel essential oils, GC-MS allowed the identification of 84 components, covering up to 99% of the total GC profile (Table S2). Major detected monoterpenes were α-terpenyl acetate (up to 18%), 1,8-cineole (up to 27%), and linalool (up to 8%). Among sesquiterpenes, the most abundant were veridiflorol (up to 4%), trans-caryophyllene (up to 4%), bicyclogermacrene (up to 1.4%), and β-elemene (up to 1.6%). Other compounds detected in significant quantities were phenylpropane derivates, methyleugenol (up to 10%), and eugenol (up to 9%).
Finally, 60 components were detected in the rosemary essential oil, covering up to 95% of total GC peak areas (Table S3). Rosemary essential oil was richest in oxygenated monoterpenes, such as borneol (up to 24%), camphor (up to 18%), linalool (up to 6%), and 1,8–cineole (up to 9%). Major sesquiterpenes were berbenone (up to 22%) and t-muurolol (up to 6%). Other compounds, like eugenol (up to 6%) and methyleugenol (up to 4%), were also present in significant amounts.

3. Discussion

We have tested the effect of different pretreatments using selected Mediterranean plants: rosemary, sage, and bay laurel, and demonstrated that different hydrodistillation pretreatments can improve the yield of respective essential oils. At the same time, the pretreatments did not significantly affect the quality of oils and their composition was comparable to no-pretreatment controls. For the first time, we have demonstrated that reflux extraction (e.g., soaking of milled plant material at 40 °C for one hour) could significantly improve essential oil yield.
Rosemary and sage are aromatic, medicinal plants that have received particular attention in the Lamiaceae family due to their aromatic and chemical composition [2,27,28], while bay laurel is a valuable medicinal plant from the Lauraceae family, which has been widely used as a spice and flavoring agent [19]. Due to the numerous bioactive compounds, the essential oils of these plants show a wide range of biological activities such as antimicrobial, preservative, antioxidant, and antifungal, which makes them valuable in a range of applications, from medicinal to the food industry [29,30,31,32,33]. Considering these valuable properties, the improvement of current and the development of novel essential oil extraction procedures is a subject of many recent studies [28,34,35].
In this study, we have obtained similar volatile profiles of sage, rosemary, and bay laurel essential oils as reported elsewhere [2,3,19,35,36,37,38,39]. Also, the essential oil composition obtained after different pretreatments was comparable to the essential oils obtained by direct hydrodistillation. Other studies had also shown that enzymatic and ultrasonic pretreatments did not affect the overall composition of the oil [11,15,16,17,21,22,26], although some authors reported that the quantities of individual major components varied significantly in relation to the extraction technique used [1,2,12,13,14,19,24]. For instance, Boulila et al. [19] reported that enzyme pretreatment of bay laurel leaves resulted in increased concentration of oxygenated monoterpenes, and explained such effect by a possible increase in oxidation after cell wall disruption or the presence of oxidases in the enzyme preparation. Also, cellulase pretreatment reportedly increased the amount of cis-verbenol and camphor and decreased the amounts of 1,8-cineol, α-pinene, fenchene, and terpinen-4-ol in rosemary essential oil [2]. We have also observed that the ratios of some components varied among different pretreatments. For example, in sage essential oil, the amount of almost all monoterpene hydrocarbons increased after HD-REPCX. On the other hand, HD-US resulted in a higher amount of major oxygenated monoterpenes (1,8–cineole, α-thujone, camphor). In bay laurel essential oil, all pretreatments, except HD-REX, resulted in an increase of 1,8-cineole, which is in line with the findings of Boulila et al. [19]. In rosemary essential oil, an increase of monoterpene hydrocarbons was noticeable after HD-US, while HD-RE resulted in an increase of 1,8-cineole and linalool, as well as phenylpropane derivatives eugenol and methyleugenol. However, such slight changes in the quantity of individual components observed after different pretreatments were not significant and did not change the overall oil quality, as confirmed by the statistical analysis of GC-MS results.
Among the different pretreatments applied in this study, the simplest was reflux extraction (HD-RE). The incubation of finely milled plant material for 1 h at 40 °C in MiliQ water before hydrodistillation resulted in up to 60% essential oil yield increase. It was previously shown that such soaking of the plant material in water [40] or acidic medium [1] before distillation may increase the quantity of oil by enhancing the leaching of ingredients from the already disrupted cells. Also, swelling and hydration of plant material, which enlarge the pores in cell walls and increase the turgor pressure in the still intact plant cells, might lead to the enhanced diffusion of the oil ingredients into the soaking medium. For example, soaking of agarwood in lactic acid for 168 h [41] and soaking of thyme leaves in distilled water overnight at 50 °C improved the essential oil yield [40]. Thus, soaking of plant material before hydrodistillation, as applied here for selected Mediterranean plants, presents a simple and cost-effective treatment that results in a significant increase in essential oil yield.
In addition, we have tested whether the application of cell wall-degrading enzymes prior to hydrodistillation has a positive effect on the extraction efficiency. We have followed previously described protocols [2,19] and applied separate and combined pretreatments with cellulase, pectinase, and xylanase to plant material. Prior to the extractions, we used a small-scale enzyme assay and confirmed the activity of all enzymes under the reaction conditions. However, the application of enzymes did not result in a significant increase in the yield of the essential oil above the reflux extraction pretreatments, suggesting that, in our case, the enzymatic activity had only a minor effect on the yield. On the contrary, the positive impact of enzyme application on the yield of essential oils was reported in multiple recent studies [2,19,21,22,23,25,26]. For instance, enzyme-assisted extraction pretreatment was used prior to hydrodistillation of bay laurel [19], rosemary, and thyme leaves [2], and reportedly resulted in an increase in essential oil yield (up to 109% for thyme leaves). However, these studies were lacking a no-enzyme control, and the observed yield increase was calculated by comparison with no-pretreatment control. In our case, the extraction results were compared to no-pretreatment control (HD) as well as to reflux extraction pretreatment (HD-RE), that served as a no-enzyme control since it was performed in the same conditions as the enzyme-assisted extraction pretreatments (MiliQ water at 40 °C for 1 h), only without addition of enzyme(s). In comparison to no-pretreatment control, an increase in essential oil yield was significant for all enzyme pretreatments individually, as well as for their combination, and there were no significant differences between each individual pretreatment. Unexpectedly, reflux extraction resulted in approximately the same essential oil increase as the enzyme-assisted pretreatments, leading to the conclusion that essential oil yield increase can be more attributed to reflux extraction pretreatment per se (i.e., soaking of macerated plant material in warm water) than to enzymatic degradation of the cell wall. The positive effect of the soaking of plant material alone, as observed here, could not be properly assessed and compared with the effect of enzymes in the previous studies that were lacking no-enzyme controls [2,19,21,22,23,24,25,26]. Our results are congruent with the results of Dimaki and coworkers [1] who analyzed the effect of enzyme pretreatments on hydrodistillation and ultrasound-assisted extraction of Sideritis sp. essential oil. They have used both no-pretreatment and no-enzyme control (i.e., preincubation in acidic medium) and also found that enzymatic pretreatment was not superior to the mere soaking of plant material in an acidic buffer. Also, acidic and enzymatic pretreatments prior to the ultrasound-assisted extraction resulted in a similar increase in extraction yield in comparison to no-pretreatment control. In conclusion, our results and results of Dimaki et al. [1] point out that the application of enzymatic pretreatments with the aim to increase the yield of essential oils, although often reported, should be carefully re-evaluated.
Finally, we have demonstrated the beneficial effect of ultrasound pretreatment on the overall extraction efficiency. In comparison to the control experiment, ultrasound-assisted extraction pretreatment caused a 50% improvement in sage and rosemary and 60% in bay laurel essential oil yield, similar to the results of reflux-assisted extraction. This is comparable with the results of other studies [12,13,15,16,17]. For instance, 30 min-ultrasonic-maceration pretreatment prior to the steam distillation of milled Thymus vulgaris L., Mentha piperita L., and Origanum majorana L. leaves resulted in approximately 10% essential oil yield increase, compared to no-pretreatment control [15,16]. Ultrasound-assisted extraction as a pretreatment to hydrodistillation of grounded carrot seeds increased essential oil yield by approximately 33% [17]. Hydrodistillation of grounded Elettaria cardamomum L. seeds with ultrasonic pretreatment resulted in a 4.9% essential oil yield increase and extraction lasted less than an hour, compared to control, which lasted 6 h [13]. Even when there was no yield increase, the application of ultrasound resulted in a significant shortening of the essential oil isolation procedure [11,14]. In conclusion, we report a high increase in essential oil yield after ultrasound pretreatment when compared to other studies, accompanied by the shortening of the extraction time when compared to the reflux-assisted extraction pretreatment.
In conclusion, we have demonstrated the valuable effects of ultrasound and reflux extraction hydrodistillation pretreatments on the yield of essential oils of sage, rosemary, and bay laurel. The benefit of ultrasound is short extraction time, while reflux extraction is performed at mild temperature. Thus, both pretreatments are considered as gentle extraction procedures, resulting in unchanged overall oil quality, as demonstrated by essential oil composition analyzed here. Also, our results point out that when assessing the effect of enzymes as a hydrodistillation pretreatment, proper experiment controls are necessary. In our case, the soaking of plant material in the water at 40 °C was enough to disrupt the cells and additional enzyme activity was redundant.

4. Materials and Methods

4.1. Plant Material

Fresh leaves were collected from wild plants of rosemary (Rosmarinus officinalis L.), sage (Salvia officinalis L.), and bay laurel (Laurus nobilis L.) in the south Mediterranean region of Croatia in August 2018. The leaves were air-dried at room temperature (20 ± 2 °C) for one week. Dry leaves were packed in polyethylene bags and kept in a dark, dry, and cool place. Before being used for the extractions, the plant material was milled using a house blender (Tefal).

4.2. Chemicals

The following chemicals were used: 3,5-dinitrosalicylic acid (DNSA) (98%, Thermo Fisher Scientific, Maharahstra, India), sodium sulphite (Lach-ner, Brno, Chez Republic), sodium hydroxide (Lach-ner, Brno, Chez Republic), phenol (99 + %, Thermo Fisher Scientific, Maharahstra, India), potassium sodium tartrate (Sigma-Aldrich, Buchs, Switzerland), Folin–Ciocalteu reagent (Sigma-Aldrich, Buchs, Switzerland), anhydrous sodium carbonate (Lach-ner, Brno, Chez Republic), gallic acid (anhydrous) for synthesis (Merck, Darmstadt, Germany), cellulase (from Aspergillus niger) (Sigma-Aldrich, Tokyo, Japan), pectinase (from Aspergillus niger) (Sigma-Aldrich, Buchs, Switzerland), xylanase (from Theryomyces, expressed in Aspergillus oryzae) (Sigma-Aldrich, Søborg, Denmark), beechwood xylan (Biosynth, Berkshir, England, UK), pectin from citrus peel (74.0%, Sigma-Aldrich, Buchs, Switzerland), sodium carboxymethyl cellulose (Sigma-Aldrich, Buchs, Switzerland), xylose (99%, Sigma-Aldrich, Buchs, Switzerland, D-( + )-glucose (99.5%, Sigma-Aldrich, Buchs, Switzerland), D-( + )-galacturonic acid monohydrate (97.0%, Sigma Aldrich, Buchs, Switzerland), C9-C25 alkanes, deuterated chloroform for NMR spectroscopy (CDCl3-d with 0.03% v/v TMS, 99.80%, Eurisotop, Saint-Aubin, France).

4.3. 3,5-Dinitrosalicylic Acid (DNSA) Assay for the Determination of Enzyme Activity

The determination of enzyme activity was done according to a DNSA colorimetric method as described in Ghose [42] and Miller [43] with modifications. Enzymes (cellulase, pectinase, and xylanase) were added to MiliQ water to a final concentration of 0.02 mg/mL. Next, 10 mg of the appropriate substrate (cellulose, pectin, and xylan, respectively) was added to 1 mL of enzyme solution and the reaction was incubated at 40 °C for 60 min in a thermoshaker (Biosan, TS-100) at 900 rpm. The sample was then cleared by centrifugation (5 min, 1568× g). Next, sample aliquots (600 μL final volume) were mixed with 600 μL of DNSA reagent containing 10.0 g/L DNSA, 0.5 g/L sodium sulfite, 10 g/L sodium hydroxide, and 2 mL/L phenol. The mixtures were incubated for 15 min at 95 °C before adding 200 μL of 40 g/L potassium sodium tartrate solution (1.4 mL final volume). The samples were chilled on ice for 5 min, and then the absorbance at 575 nm was measured (Cary Series UV-Vis Spectrophotometer, Agilent Technologies). Product concentrations were calculated from calibration curves generated with the corresponding reducing sugars (glucose, galacturonic acid, xylose, respectively). If needed, the original samples were appropriately diluted to produce the absorption values within the range of the calibration curves, and dilutions were taken into account in the calculation of the enzymatic activity. One unit of enzymatic activity was defined as the amount of enzyme that releases 1 μmol of reducing sugar per minute under the specified assay conditions.

4.4. Determination of Total Phenols for Optimization of Ultrasonic Parameters

The total phenols of plant extracts were determined using the Folin–Ciocalteu method, according to the procedure of Singleton and Rossi [44]. Mass fraction of total phenols was expressed as mg of gallic acid equivalent (GAE) per g dry plant.

4.5. Extraction Procedures

Hydrodistillation (HD) was preceded by different pretreatments: hydrodistillation with reflux extraction pretreatment (HD-RE), hydrodistillation with reflux extraction pretreatment assisted with enzymes (pectinase, HD-REP; cellulase, HD-REC; xylanase, HD-REX; pectinase/cellulase/xylanase, HD-REPCX), and hydrodistillation with ultrasonic pretreatment (HD-U). Additionally, HD without pretreatment was performed as a negative control.
For each protocol, 20 g of milled plant material was mixed with 250 mL MiliQ water. For HD-RE, this mixture was subjected to reflux extraction at 40 °C for 1 h under stirring conditions. For enzyme-assisted pretreatments, we have followed previously described protocols [2,19] with slight modifications. Namely, the plant material in MiliQ water was subjected to reflux extraction at 40 °C for 1 h under stirring conditions in the presence of 10 mg pectinase, cellulase, or xylanase, that is, 0.5 mg of enzyme per g of dried plant material, or their combination (5 mg of each enzyme, i.e., 0.75 mg of total enzymes per g of dried plant material). For HD-U, the plant material in MiliQ water was treated with a 14 mm diameter ultrasonic probe (ultrasonic device UP200Ht, Hielscher, Teltow, Germany), at 30% of the maximal ultrasonic power for 10 min. Before the ultrasound extraction, optimal duration, and amplitude were determined (Table 1).
Next, the plant/water mixture was subjected to hydrodistillation using a Clevenger apparatus for 2.5 h. Every extraction was done once. The essential oil was stored in dark glass vials at 4 °C until further analyses.

4.6. Gas Chromatograph/Mass Spectrometer (GC-MS) Analysis

GC-MS analyses were carried out with an Agilent Technologies (Palo Alto, CA, USA) gas chromatograph model 7890 A equipped with a mass selective detector (MSD) model 5975 C (Agilent Technologies, Palo Alto, CA, USA) and an HP-5MS 5% phenyl-methylpolysiloxane capillary column (30 m × 0.25 mm, 0.25 μm film thickness, Agilent Technologies, Palo Alto, CA, USA). In brief, the injector and detector temperatures were 250 °C and 300 °C, respectively; the column temperature was held at 70 °C for 2 min, and was then increased from 70 to 200 °C at 3 °C/min, and was finally held at 200 °C for 18 min; 1.0 μL of the sample (10 µL of the oil dissolved in 1 mL of pentane) was injected using split mode (split ratio 1:50). Helium was used as carrier gas (1.0 mL/min). The MSD (EI mode) was operated at 70 eV, the ion source temperature was 230 °C, and the scan range was set to 30–350 amu. Identification of volatile constituents was based on the comparison of their retention indices (RIs), determined relative to the retention times of a homologous series of n-alkanes (C9–C25), with those reported in the literature and their mass spectra with those of authentic compounds available in our laboratories or those listed in the NIST 08 and Wiley 9 mass spectral libraries (matches higher than 90%). Relative concentrations of components were calculated by the area normalization method without considering response factors.

4.7. Statistical Analyses

All variables were log-transformed (using base 10 logs) to improve the data distribution and homogeneity of variances. The data were tested for normality by comparing histograms of the sample data to a normal probability curve, after which the null hypothesis (that data were normally distributed) was rejected. To compare and plot the effects of different pretreatments on the chemical composition of essential oils for each plant (as determined by GC-MS), we applied Spearman’s nonparametric measure of rank correlation using program R v. 3.2.0 [45]. The level of significance was set at p ≤ 0.05.

Supplementary Materials

The following are available online at, Figure S1: Sage 1H NMR spectrum (600 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 32 K data points; 256 scans; 0.37 Hz/point; 1 s delay), Figure S2: Sage 1H-1H COSY NMR spectrum (600 MHz, CDCl3-d, 25 °C), Figure S3: Sage 1H-1H TOCSY NMR spectrum (600 MHz, CDCl3-d, 25 °C), Figure S4: Sage 1H-13C HMQC NMR spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S5: Sage 1H-13C HMBC NMR spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S6: Sage 13C APT NMR spectrum (150 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 64 K data points; 44506 scans; 0.60 Hz/point; 1 s delay), Figure S7: Bay laurel 1H NMR spectrum (600 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 32 K data points; 256 scans; 0.37 Hz/point; 1 s delay), Figure S8: Bay laurel 1H-1H COSY NMR spectrum (600 MHz, CDCl3-d, 25 °C), Figure S9: Bay laurel 1H-1H TOCSY NMR spectrum (600 MHz, CDCl3-d, 25 °C), Figure S10: Bay laurel 1H-13C HMQC NMR spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S11: Bay laurel 1H-13C HMBC NMR spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S12: Bay laurel 13C APT NMR spectrum (150 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 64 K data points; ca. 30000 scans; 0.60 Hz/point; 1 s delay), Figure S13: Rosemary 1H NMR spectrum (600 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 32 K data points; 128 scans; 0.37 Hz/point; 1 s delay), Figure S14: Rosemary 1H-1H COSY spectrum (600 MHz, CDCl3-d, 25 °C), Figure S15: Rosemary 1H-1H TOCSY spectrum (600 MHz, CDCl3-d, 25 °C), Figure S16: Rosemary 1H-13C HMQC spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S17: Rosemary 1H-13C HMBC spectrum (CDCl3-d, 25 °C). The 600 MHz 1H NMR spectrum is shown at the top edge and a 150 MHz 13C NMR spectrum at the left-hand edge, Figure S18: Rosemary 13C APT NMR spectrum (150 MHz, 0.5 mL CDCl3; 5 mm sample tube; 25 °C; 64 K data points; ca. 34000 scans; 0.60 Hz/point; 1 s delay), Figure S19: (A) Rosemary, (B) sage and (C) bay laurel essential oils 1H NMR spectra at 600 MHz in CDCl3-d. The enumeration scheme used for the assignment of the NMR spectra is shown for every compound. Figure S20: Scatter plot showing correlations between different pretreatments (on the diagonal) regarding the chemical composition of sage essential oils. Significant p-values based on Spearman’s rank test are shown above the diagonal, while bivariate scatter plots are shown below the diagonal, Figure S21: Scatter plot showing correlations between different pretreatments (on the diagonal) regarding the chemical composition of bay laurel essential oils. Significant p-values based on Spearman’s rank test are shown above the diagonal, while bivariate scatter plots are shown below the diagonal, Figure S22: Scatter plot showing correlations between different pretreatments (on the diagonal) regarding the chemical composition of rosemary essential oils. Significant p-values based on Spearman’s rank test are shown above the diagonal, while bivariate scatter plots are shown below the diagonal. Table S1. Chemical composition of sage essential oils isolated by hydrodistillation with and without different pretreatments, Table S2. Chemical composition of bay laurel essential oils isolated by hydrodistillation with and without different pretreatments, Table S3. Chemical composition of rosemary essential oils isolated by hydrodistillation with and without different pretreatments.

Author Contributions

Conceptualization, A.B., and M.D.; methodology, M.D., Z.M., I.J., M.A., T.R., S.R., and D.V.-T.; validation, A.M., M.D., and A.B.; formal analysis, D.G.; investigation, A.M., Z.M., I.J., M.A., T.R., S.R., and D.V.-T.; resources, A.B..; writing—original draft preparation, A.M.; writing—review and editing, A.B. and M.D.; visualization, A.M., D.G., and S.R.; supervision, A.B., and M.D.; project administration, A.B. and D.G.; funding acquisition, A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the project “Interactions of freshwater pathogenic oomycetes and the environment” (InteractOomyc, UIP-2017-05-6267), Croatian Science Foundation.

Acknowledgments

The authors are thankful to Maja Macut and Gloria Hradec for technical assistance during the study and Professor Verica Dragović-Uzelac for her support and advice during this research.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Dimaki, V.D.; Iatrou, G.; Lamari, F.N. Effect of acidic and enzymatic pretreatment on the analysis of mountain tea (Sideritis spp.) volatiles via distillation and ultrasound-assisted extraction. J. Chromatogr. A 2017, 1524, 290–297. [Google Scholar] [CrossRef] [PubMed]
  2. Hosni, K.; Hassen, I.; Chaâbane, H.; Jemli, M.; Dallali, S.; Sebei, H.; Casabianca, H. Enzyme-assisted extraction of essential oils from thyme (Thymus capitatus L.) and rosemary (Rosmarinus officinalis L.): Impact on yield, chemical composition and antimicrobial activity. Ind. Crop. Prod. 2013, 47, 291–299. [Google Scholar] [CrossRef]
  3. Olmedo, R.H.; Asensio, C.M.; Grosso, N.R. Thermal stability and antioxidant activity of essential oils from aromatic plants farmed in Argentina. Ind. Crop. Prod. 2015, 69, 21–28. [Google Scholar] [CrossRef]
  4. Amorati, R.; Foti, M.C.; Valgimigli, L. Antioxidant activity of essential oils. J. Agric. Food Chem. 2013, 61, 10835–10847. [Google Scholar] [CrossRef]
  5. Lesgards, J.F.; Baldovini, N.; Vidal, N.; Pietri, S. Anticancer activities of essential oils constituents and synergy with conventional therapies: A review. Phyther. Res. 2014, 28, 1423–1446. [Google Scholar] [CrossRef]
  6. Martins, I.M.; Barreiro, M.F.; Coelho, M.; Rodrigues, A.E. Microencapsulation of essential oils with biodegradable polymeric carriers for cosmetic applications. Chem. Eng. J. 2014. [Google Scholar] [CrossRef] [Green Version]
  7. Ribeiro-Santos, R.; Carvalho-Costa, D.; Cavaleiro, C.; Costa, H.S.; Albuquerque, T.G.; Castilho, M.C.; Ramos, F.; Melo, N.R.; Sanches-Silva, A. A novel insight on an ancient aromatic plant: The rosemary (Rosmarinus officinalis L.). Trends Food Sci. Technol. 2015. [Google Scholar] [CrossRef]
  8. Trongtokit, Y.; Rongsriyam, Y.; Komalamisra, N.; Apiwathnasorn, C. Comparative repellency of 38 essential oils against mosquito bites. Phyther. Res. 2005, 19, 303–309. [Google Scholar] [CrossRef]
  9. Périno-Issartier, S.; Ginies, C.; Cravotto, G.; Chemat, F. A comparison of essential oils obtained from lavandin via different extraction processes: Ultrasound, microwave, turbohydrodistillation, steam and hydrodistillation. J. Chromatogr. A 2013, 1305, 41–47. [Google Scholar] [CrossRef]
  10. Toma, M.; Vinatoru, M.; Paniwnyk, L.; Mason, T.J. Investigation of the effects of ultrasound on vegetal tissues during solvent extraction. Ultrason. Sonochem. 2001, 8, 137–142. [Google Scholar] [CrossRef]
  11. Assami, K.; Pingret, D.; Chemat, S.; Meklati, B.Y.; Chemat, F. Ultrasound induced intensification and selective extraction of essential oil from Carum carvi L. seeds. Chem. Eng. Process. Process Intensif. 2012, 62, 99–105. [Google Scholar] [CrossRef]
  12. Lilia, C.; Abdelkader, A.; Karima, A.K.A.; Tarek, B. The Effect of Ultrasound Pre-treatment on the Yield, Chemical Composition and Antioxidant Activity of Essential Oil from Wild Lavandula stoechas L. J. Essent. Oil-Bear. Plants 2018, 21, 253–263. [Google Scholar] [CrossRef]
  13. Morsy, N.F.S. A short extraction time of high quality hydrodistilled cardamom (Elettaria cardamomum L. Maton) essential oil using ultrasound as a pretreatment. Ind. Crop. Prod. 2015, 62, 287–292. [Google Scholar] [CrossRef]
  14. Seidi Damyeh, M.; Niakousari, M.; Saharkhiz, M.J. Ultrasound pretreatment impact on Prangos ferulacea Lindl. and Satureja macrosiphonia Bornm. essential oil extraction and comparing their physicochemical and biological properties. Ind. Crop. Prod. 2016, 10, 22–28. [Google Scholar] [CrossRef]
  15. Kowalski, R.; Wawrzykowski, J. Effect of ultrasound-assisted maceration on the quality of oil from the leaves of thyme Thymus vulgaris L. Flavour Fragr. J. 2009, 24, 214–220. [Google Scholar] [CrossRef]
  16. Kowalski, R.; Kowalska, G.; Jamroz, J.; Nawrocka, A.; Metyk, D. Effect of the ultrasound-assisted preliminary maceration on the efficiency of the essential oil distillation from selected herbal raw materials. Ultrason. Sonochem. 2015, 24, 214–220. [Google Scholar] [CrossRef]
  17. Smigielski, K.B.; Majewska, M.; Kunicka-Styczyñska, A.; Gruska, R. The Effect of Ultrasound-Assisted Maceration on the Bioactivity, Chemical Composition and Yield of Essential Oil from Waste Carrot Seeds (Daucus carota). J. Essent. Oil-Bear. Plants 2014, 17, 1075–1086. [Google Scholar] [CrossRef]
  18. Balasubramaniam, V.G.; Ayyappan, P.; Sathvika, S.; Antony, U. Effect of enzyme pretreatment in the ultrasound assisted extraction of finger millet polyphenols. J. Food Sci. Technol. 2019, 74, 485–493. [Google Scholar] [CrossRef] [PubMed]
  19. Boulila, A.; Hassen, I.; Haouari, L.; Mejri, F.; Amor, I.B.; Casabianca, H.; Hosni, K. Enzyme-assisted extraction of bioactive compounds from bay leaves (Laurus nobilis L.). Ind. Crop. Prod. 2015, 74, 485–493. [Google Scholar] [CrossRef]
  20. Puri, M.; Sharma, D.; Barrow, C.J. Enzyme-assisted extraction of bioactives from plants. Trends Biotechnol. 2012, 30, 37–44. [Google Scholar] [CrossRef]
  21. Sowbhagya, H.B.; Srinivas, P.; Purnima, K.T.; Krishnamurthy, N. Enzyme-assisted extraction of volatiles from cumin (Cuminum cyminum L.) seeds. Food Chem. 2011, 127, 1856–1861. [Google Scholar] [CrossRef]
  22. Sowbhagya, H.B.; Srinivas, P.; Krishnamurthy, N. Effect of enzymes on extraction of volatiles from celery seeds. Food Chem. 2010, 120, 230–234. [Google Scholar] [CrossRef]
  23. Baby, K.C.; Ranganathan, T.V. Effect of enzyme pre-treatment on extraction yield and quality of cardamom (Elettaria cardamomum maton.) volatile oil. Ind. Crop. Prod. 2016, 89, 200–206. [Google Scholar] [CrossRef]
  24. Chandran, J.; Amma, K.P.P.; Menon, N.; Purushothaman, J.; Nisha, P. Effect of enzyme assisted extraction on quality and yield of volatile oil from black pepper and cardamom. Food Sci. Biotechnol. 2012, 21, 1611–1617. [Google Scholar] [CrossRef]
  25. Chávez-González, M.L.; López-López, L.I.; Rodríguez-Herrera, R.; Contreras-Esquivel, J.C.; Aguilar, C.N. Enzyme-assisted extraction of citrus essential oil. Chem. Pap. 2016, 70, 412–417. [Google Scholar] [CrossRef]
  26. Śmigielski, K.B.; Majewska, M.; Kunicka-Styczyńska, A.; Szczesna-Antczak, M.; Gruska, R.; Stańczyk, Ł. The effect of enzyme-assisted maceration on bioactivity, quality and yield of essential oil from waste carrot (Daucus carota) seeds. J. Food Qual. 2014, 37, 219–228. [Google Scholar] [CrossRef]
  27. Hamrouni-Sellami, I.; Rahali, F.Z.; Rebey, I.B.; Bourgou, S.; Limam, F.; Marzouk, B. Total Phenolics, Flavonoids, and Antioxidant Activity of Sage (Salvia officinalis L.) Plants as Affected by Different Drying Methods. Food Bioprocess Technol. 2013, 59, 257–268. [Google Scholar] [CrossRef]
  28. Ali, A.; Chua, B.L.; Chow, Y.H. An insight into the extraction and fractionation technologies of the essential oils and bioactive compounds in Rosmarinus officinalis L.: Past, present and future. TrAC Trends Anal. Chem. 2019, 118, 338–351. [Google Scholar] [CrossRef]
  29. Ali, M.B.; Mnafgui, K.; Feki, A.; Damak, M.; Allouche, N. In vitro Antidiabetic, anti-obesity and antioxidant proprities of Rosemary extracts. J. Adv. Chem. 2014, 10, 2305–2316. [Google Scholar] [CrossRef]
  30. Fidan, H.; Stefanova, G.; Kostova, I.; Stankov, S.; Damyanova, S.; Stoyanova, A.; Zheljazkov, V.D. Chemical Composition and Antimicrobial Activity of Laurus nobilis L. Essential oils from Bulgaria. Molecules 2019, 24, 804. [Google Scholar] [CrossRef] [Green Version]
  31. Soares, I.H.; Loreto, S.; Rossato, L.; Mario, D.N.; Venturini, T.P.; Baldissera, F.; Santurio, J.M.; Alves, S.H. In vitro activity of essential oils extracted from condiments against fluconazole-resistant and -sensitive Candida glabrata. J. Mycol. Med. 2015, 54, 950–956. [Google Scholar] [CrossRef] [PubMed]
  32. Fernández, N.J.; Damiani, N.; Podaza, E.A.; Martucci, J.F.; Fasce, D.; Quiroz, F.; Meretta, P.E.; Quintana, S.; Eguaras, M.J.; Gende, L.B. Laurus nobilis L. Extracts against Paenibacillus larvae: Antimicrobial activity, antioxidant capacity, hygienic behavior and colony strength. Saudi J. Biol. Sci. 2019. [Google Scholar] [CrossRef]
  33. Dammak, I.; Hamdi, Z.; Kammoun El Euch, S.; Zemni, H.; Mliki, A.; Hassouna, M.; Lasram, S. Evaluation of antifungal and anti-ochratoxigenic activities of Salvia officinalis, Lavandula dentata and Laurus nobilis essential oils and a major monoterpene constituent 1,8-cineole against Aspergillus carbonarius. Ind. Crop. Prod. 2019, 128, 85–93. [Google Scholar] [CrossRef]
  34. Ferreira, D.F.; Lucas, B.N.; Voss, M.; Santos, D.; Mello, P.A.; Wagner, R.; Cravotto, G.; Barin, J.S. Solvent-free simultaneous extraction of volatile and non-volatile antioxidants from rosemary (Rosmarinus officinalis L.) by microwave hydrodiffusion and gravity. Ind. Crop. Prod. 2020, 145, 112094. [Google Scholar] [CrossRef]
  35. Al Zuhairi, J.J.M.J.; Jookar Kashi, F.; Rahimi-Moghaddam, A.; Yazdani, M. Antioxidant, cytotoxic and antibacterial activity of Rosmarinus officinalis L. essential oil against bacteria isolated from urinary tract infection. Eur. J. Integr. Med. 2020, 38, 101192. [Google Scholar] [CrossRef]
  36. Vosoughi, N.; Gomarian, M.; Ghasemi Pirbalouti, A.; Khaghani, S.; Malekpoor, F. Essential oil composition and total phenolic, flavonoid contents, and antioxidant activity of sage (Salvia officinalis L.) extract under chitosan application and irrigation frequencies. Ind. Crop. Prod. 2018, 117, 366–374. [Google Scholar] [CrossRef]
  37. Khiva, Z.; Havani, M.; Garnar, A.; Kharchouf, S.; Amine, S.; Berrekhis, F.; Bouzoubae, A.; Zair, T.; Elhilali, F. Valorization of the Salvia officinalis L of the Morocco bioactive extracts: Phytochemistry, antioxidant and anticorrosive activities. J. King Saud Univ. Sci. 2019. [Google Scholar] [CrossRef]
  38. Hatipoglu, S.D.; Zorlu, N.; Dirmenci, T.; Goren, A.C.; Ozturk, T.; Topcu, G. Determination of volatile organic compounds in fourty five salvia species by thermal desorption-GC-MS technique. Rec. Nat. Prod. 2016, 10, 659–700. [Google Scholar]
  39. Russo, A.; Formisano, C.; Rigano, D.; Senatore, F.; Delfine, S.; Cardile, V.; Rosselli, S.; Bruno, M. Chemical composition and anticancer activity of essential oils of Mediterranean sage (Salvia officinalis L.) grown in different environmental conditions. Food Chem. Toxicol. 2013, 55, 42–47. [Google Scholar] [CrossRef]
  40. Awada, F.; Kobaissi, A.; Chokr, A.; Hamze, K.; Hayar, S.; Mortada, A. Factors affecting quantitative and qualitative variation of thyme (Origanum syriacum L.) essential oil in Lebanon. Adv. Environ. Biol. 2012, 6, 1509–1514. [Google Scholar]
  41. Nor Fazila, K.; Ku Halim, K.H. Effects of soaking on yield and quality of agarwood oil. J. Trop. For. Sci. 2012, 24, 557–564. [Google Scholar]
  42. Ghose, T.K. Measurement of cellulase activities. Pure Appl. Chem. 1987, 59, 257–268. [Google Scholar] [CrossRef]
  43. Miller, G.L. Use of Dinitrosalicylic Acid Reagent for Determination of Reducing Sugar. Anal. Chem. 1959, 31, 426–428. [Google Scholar] [CrossRef]
  44. Singleton, V.L.; Rossi, J.A. Colorimetry of total Phenolics with Phosphomolybdic-Phosphotungstic Acid reagents. Am. J. Enol. Vitic. 1965, 16, 144. [Google Scholar]
  45. R Development Core Team. R: A Language and Environment for Statistical Computing; R Development Core Team: Vienna, Austria, 2017. [Google Scholar]
Sample Availability: Samples of the compounds are not available from the authors.
Figure 1. Effect of different hydrodistillation pretreatments on the extraction yield of rosemary, sage, and bay laurel essential oil. Hydrodistillation without pretreatment (negative control)—HD; hydrodistillation with reflux extraction pretreatment—HD-RE; hydrodistillation with reflux extraction pretreatment assisted with enzymes: xylanase—HD-REX; cellulase—HD-REC; pectinase—HD-REP; pectinase + cellulase + xylanase—HD-REPCX; hydrodistillation with ultrasonic pretreatment—HD-US.
Figure 1. Effect of different hydrodistillation pretreatments on the extraction yield of rosemary, sage, and bay laurel essential oil. Hydrodistillation without pretreatment (negative control)—HD; hydrodistillation with reflux extraction pretreatment—HD-RE; hydrodistillation with reflux extraction pretreatment assisted with enzymes: xylanase—HD-REX; cellulase—HD-REC; pectinase—HD-REP; pectinase + cellulase + xylanase—HD-REPCX; hydrodistillation with ultrasonic pretreatment—HD-US.
Molecules 25 04818 g001
Figure 2. Overview of the chemical composition of sage (A), bay laurel (B), and rosemary (C) essential oils after different hydrodistillation pretreatments, as determined by GC-MS. Hydrodistillation without pretreatment (negative control)—HD; hydrodistillation with reflux extraction pretreatment—HD-RE; hydrodistillation with reflux extraction pretreatment assisted with enzymes: xylanase—HD-REX; cellulase—HD-REC; pectinase—HD-REP; pectinase + cellulase + xylanase—HD-REPCX; hydrodistillation with ultrasonic pretreatment—HD-US.
Figure 2. Overview of the chemical composition of sage (A), bay laurel (B), and rosemary (C) essential oils after different hydrodistillation pretreatments, as determined by GC-MS. Hydrodistillation without pretreatment (negative control)—HD; hydrodistillation with reflux extraction pretreatment—HD-RE; hydrodistillation with reflux extraction pretreatment assisted with enzymes: xylanase—HD-REX; cellulase—HD-REC; pectinase—HD-REP; pectinase + cellulase + xylanase—HD-REPCX; hydrodistillation with ultrasonic pretreatment—HD-US.
Molecules 25 04818 g002
Table 1. Optimization of ultrasound extraction parameters based on the total phenol yield.
Table 1. Optimization of ultrasound extraction parameters based on the total phenol yield.
TimeTotal Phenol Yield (mg/g)
30% of Max. Ultrasonic Power60% of Max. Ultrasonic Power90% of Max. Ultrasonic Power
Bay laurel5 min16.87 ± 0.188.86 ± 0.057.32 ± 0.14
10 min17.37 ± 0.538.27 ± 0.215.99 ± 0.00
15 min16.86 ± 0.188.30 ± 0.216.38 ± 0.15
Sage5 min29.46 ± 3.1432.95 ± 2.632.20 ± 2.9
10 min35.71 ± 5.634.24 ± 7.236.83 ± 3.9
15 min31.25 ± 5.431.76 ± 3.431.82 ± 0.32
Rosemary5 min93.98 ± 1.4146.35 ± 1.0650.16 ± 2.75
10 min97.44 ± 2.1245.27 ± 0.0751.22 ± 0.78
15 min103.44 ± 2.1243.08 ± 0.1443.12 ± 0.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Miljanović, A.; Bielen, A.; Grbin, D.; Marijanović, Z.; Andlar, M.; Rezić, T.; Roca, S.; Jerković, I.; Vikić-Topić, D.; Dent, M. Effect of Enzymatic, Ultrasound, and Reflux Extraction Pretreatments on the Yield and Chemical Composition of Essential Oils. Molecules 2020, 25, 4818. https://doi.org/10.3390/molecules25204818

AMA Style

Miljanović A, Bielen A, Grbin D, Marijanović Z, Andlar M, Rezić T, Roca S, Jerković I, Vikić-Topić D, Dent M. Effect of Enzymatic, Ultrasound, and Reflux Extraction Pretreatments on the Yield and Chemical Composition of Essential Oils. Molecules. 2020; 25(20):4818. https://doi.org/10.3390/molecules25204818

Chicago/Turabian Style

Miljanović, Anđela, Ana Bielen, Dorotea Grbin, Zvonimir Marijanović, Martina Andlar, Tonči Rezić, Sunčica Roca, Igor Jerković, Dražen Vikić-Topić, and Maja Dent. 2020. "Effect of Enzymatic, Ultrasound, and Reflux Extraction Pretreatments on the Yield and Chemical Composition of Essential Oils" Molecules 25, no. 20: 4818. https://doi.org/10.3390/molecules25204818

APA Style

Miljanović, A., Bielen, A., Grbin, D., Marijanović, Z., Andlar, M., Rezić, T., Roca, S., Jerković, I., Vikić-Topić, D., & Dent, M. (2020). Effect of Enzymatic, Ultrasound, and Reflux Extraction Pretreatments on the Yield and Chemical Composition of Essential Oils. Molecules, 25(20), 4818. https://doi.org/10.3390/molecules25204818

Article Metrics

Back to TopTop