Introduction

The improper and uncontrolled use of antibiotics and the high ability of bacteria to adapt to adverse environmental conditions have led to the emergence of multidrug-resistant strains in bacteria, which is a significant threat to public health1. In the same way, resistance to various fungicides in many pathogenic strains of fungi is a similar consequence of the uncontrolled use of chemical antimicrobials2. The promising results of herbs as natural sources in pharmaceutical and food industries, along with the increasing public interest in herbal medicines, have led many researchers in various scientific and applied fields to measure the composition of each plant and the applied fields of their metabolites3. Hence, ethnobotanical and phytochemical investigations are essential in detecting secondary medicinal plant components to determine their biological properties. Essential oils, as rich sources of plant secondary metabolites, have antioxidant, insecticidal, fungicidal, bactericidal, and even herbicidal properties. To this end, the application of these compounds has increased in various industries4.

Due to the hydrophobic nature of essential oils, these compounds cause disruption and destruction in the cell membrane of fungi, yeasts, and bacteria, followed by considerable cellular damage and, eventually, cell death in these microorganisms. However, the intensity of these properties largely depends on the chemical profile of each essential oil5,6. The family of Apiaceae Lindl. includes economically important vegetables, herbs, and spices, where it is considered one of the most highly frequent plant families with well-known antioxidant, antibacterial, and antifungal properties7. In different organs of plants (e.g., fruits, aerial organs, and roots) of this family, there are schizogenous secretory canals and ducts for essential oils that store these valuable metabolites with different qualities and quantities in each organ8. Iran is an important country in the world regarding the abundance of species of the Apiaceae family, and is ranked third after China and Turkey9. The properties of many endemic and wild plants of this family in Iran are unknown throughout the world. Therefore, in addition to the ecological and genetic values of these plants, the identification of compounds, medicinal, and antioxidant properties of these plants also assumes crucial importance.

As it has been mentioned, the essence of Apiaceae family plants contains a wide range of compounds with varying biological properties. K. odoratissima Mozaff. is an endemic and endangered medicinal plant available in the flora of Iran, which is used as a traditional medicinal and spice plant (particularly in distribution areas)10. A major part of the essence of K. odoratissima plant consists of (Z)-ligustilide (phthalide) and germacrene D (monocyclic sesquiterpene). In addition to its traditional uses as a fresh and processed vegetable in dairy products, this plant is used as a medicine for the treatment of rheumatic disorders, blood pressure, indigestion, inflammation, and antispasmodic pains as well as the treatment of digestive problems, due to its compounds11. On the other hand, many properties and applications have been discovered in some more well-known species of the family plant. These plants have different uses based on the compounds existing in the fruit essence. For example, cuminaldehyde, as the main compound of Cuminum cyminum L, has commercial applications in the perfumery and cosmetics industry, especially as an antifungal, antibacterial and antioxidant agent12. Foeniculum vulgare and Pimpinella anisum fruit oil (rich in anethole) are consumed in moisturizing cream formulas, skin creams, lotions, ointments, other cosmetics, and personal care products. Furthermore, it is used as a perfumery agent for oral care products13. Carvone is one of the other compounds found in the essence of Anethum graveolens plant. This compound has anti-cholesterol and antioxidant activity in medicine, insecticidal properties for pest control, and an inhibitory effect on the growth of potato tubers. It has been proven its usability for potatoes in the post-harvest storage and transportation stages14. F. gummosa Boiss is another medicinal plant of the Apiaceae family of the genus Ferula L. that grows in southwest Iran and the foothills of Zagros Mountains15. Oleo-gum-resin contains valuable monoterpene secondary metabolites and its beneficial applications have been proven in various medicinal items such as laxative, expectorant, anti-nociceptive, antiseptic, anti-hysteric, and Anti-catarrh items, besides its use in glue manufacturing for industrial and cosmetics consumption16. K. odoratissima and F. gummosa are generally less known, and there are few reports about their essential oil properties, while these species are at risk of extinction in Iran. Besides them, Heracleum persicum Desf. (genus Heracleum) is another plant of the Apiaceae family. Unfortunately, there are limited studies pertaining to its antioxidant and antimicrobial activity. H. persicum is an annual Iranian plant that grows in humid and mountainous areas with an altitude of more than 1500 m in the north, northwest, center, and northeast of Iran, especially in the provinces of Mazandaran, Azerbaijan, Khorasan, Tehran, Yazd and Kerman17,18. It is normally used as a flavouring and spicy in traditional culture, while it is also used in folk Iranian medicine as a treatment for diarrhoea, indigestion, flatulence, infection, and vomiting constipation. Some studies have reported anticonvulsant properties for this plant19,20,21.

Although studies on some well-known species of the Apiaceae family have begun years ago, there are still shortages or contradictions between these reports. For example, prior reports have focused on the essential oil of the aerial part of these plants. In comparison, the fruit of many of these plants and their essential oil have consumption for spice, folk medicine, alcoholic beverages, flavourings, and perfumes. On the other hand, the high level of contradiction in reports on the essential oil properties of the Apiaceae family does not provide clear insight into each plant's potential compared to other family plants. For instance, Oroojalian et al.22 reported that the range of minimum inhibitory concentration (MIC) was about 7.5 and 3 mg/mL in Cuminum cyminum essential oil for bacteria such as Staphylococcus aureus and Escherichia coli. Whereas, other reports about the essential oil of this plant showed that MIC was about 25 and 50 µg/mL in S. aureus and E. coli bacteria23. These mentioned distinctions have existed in other plants related to this family. For example, regarding Foeniculum vulgare, the MIC range of fruit essential oil against S. aureus and E. coli bacteria has been reported equal to 250 and 125 µg/mL24. However, in another research, MIC was > 10.0 and 0.25 mg/mL against the same bacteria25. Alternatively, about Anethum graveolens fruit essential oil, the MIC range against S. aureus and E. coli ranged from > 10.0 mg/mL to 5 µg/mL in different studies26,27. Although the compositions of essential oils and bacteria strains may be different in these studies, the striking differences between mg/mL and µg/mL in inhibitory concentrations make it very hard to select and classify plants according to their specific properties, especially in intra-familial classes. In addition, these distinctions could be found in many reports concerning antifungal and antioxidant properties. Hence, the absence of any clear vision of the wild and common species of these plants makes an unreliable condition and prevents the expansion in the practical and scientific fields. Therefore, the purpose of the present study is to investigate some biological properties of the nine species of the Apiaceae family, considering the results of prior studies on the studied species. In this way, a comprehensive study of these species can lead to a summary of the properties of each plant species based on the chemical profile of their essential oils. In the same way, in order to accurately compare the effect of the essential oils of the studied plants on the destruction levels of bacteria, one gram-positive and one gram-negative bacterium were selected and the effects of these essential oils on the two bacteria were examined by Scanning Electron Microscope.

Material and methods

Materials

In the present study, the materials included Mueller hinton agar (MHA) (Merck), Mueller hinton broth (MHB) (HIMEDIA), Sabouraud dextrose agar (SDA) (QUELAB), Sabouraud dextrose broth (SDB) from (QUELAB), Gentamicin (10 µg per disc) and Ciprofloxacin (5 µg per disc) antibiotic disc diffusion (Padtanteb, Iran), Voriconazole Antibiotic disk diffusion (1 µg per disk) (RIDACOM). These materials were all purchased; moreover, all chemicals except 2,2-diphenyl-1-picrylhydrazyl (DPPH) (Sigma Aldrich) were purchased from Merck.

Plant material

Wild plant material was collected with the consent of the Natural Resources and Watershed Management Organization. Since some species are widely distributed and cultivated throughout the country, no legal requirements or permission was required to gather the samples. The collection of the plant material complied with the relevant Natural Resources and Watershed Management Organization and national and international guidelines and legislation. F. gummosa fruits (31° 01′ 29″ N 51° 21′ 53″ E) were freshly collected from Sisakht, Kohgiluyeh and Boyer-Ahmad Province, K. odoratissima fruits from the Khuyeh region in Chaharmahal and Bakhtiari Province (32° 41′ 11″ N 49° 49′ 42″ E), and H. persicum from Khalkhal, Ardabil, Iran during the fruiting period in August 2021. The exact species of F. gummosa (herbarium number: PM1393), K. odoratissima (herbarium number: PM1396) and H. persicum (Herbarium number: PM1392) were determined in the Herbarium of Pharmacology Department, Shiraz University of Medical Sciences, Shiraz, Iran. The fresh fruit of F. vulgare (var. azoricum), Petroselinum crispum (var. Neapolitanum) and A. graveolens was prepared from Pakan Bazr Company, Isfahan. Similarly, Fresh fruits of C. cyminum and P. anisum were purchased from Razavi Seed and Seedling Institute (Razavi Khorasan province) and B. persicum was prepared from Agricultural and Natural Resources Research Center of the Kerman province. The fruits of all the plants were dried (in the shade at room temperature). Then, the essential oil and methanol extracts were extracted and the samples were stored at 4 °C in the refrigerator.

Methanolic and EOs extraction of fruits

The essential oil of each plant was extracted several times using distillation with water by the Clevenger apparatus. The essential oils of each plant were stored at 4 °C in a dark and closed container in a dark place until use. The methanolic extract was obtained from the fruit (9 g) by a Soxhlet extractor with methanol solvent. Finally, the extract was purified by a rotary evaporator. To make the mother methanolic extract based on the amount of fruits, all the pure extract was again diluted with methanol to 30 mL (300 mg fruit/mL) and stored at 4 °C in dark and sealed containers. At the time of measuring the antioxidant properties, different concentrations (100, 50, 25, 12.5, 6.25 and 3.125 mg fruit/mL) of the main methanolic extract were prepared based on the amount of fruit extracted from the final solution volume.

Total phenolic and flavonoid content

The test solution for estimating the total phenolic content of each fruit contained methanolic extract and Folin–Ciocalteu reagent to which sodium bicarbonate was added after placing in Bain-marie. The test solution for estimating total flavonoid content also contained aluminum chloride, methanolic extract, potassium acetate, and distilled water that was placed in Bain-marie at 25 °C. Finally, the total phenolic content was measured according to the absorbance of the test solution at a wavelength of 760 nm and based on the absorbance of standard concentrations of gallic acid equivalent (GAE)28. The total flavonoid content was also measured based on the absorbance of test solutions at 415 nm using a spectrophotometer (Lambda 210 EZ) according to the absorption of standard concentrations of quercetin (quercetin equivalent or QE)29. The results are reported based on the mean of 3 replications with ± standard error.

DPPH radical scavenging activity

The antioxidant power of extracts and essential oils was reported based on the concentrations of the extract or essential oil scavenging 50% of DPPH radicals (IC50). To this end, 1 mL of different concentrations of the sample was added to 0.25 mL of 0.2 mM DPPH solution. Accordingly, the absorbance of each sample was read at 517 nm using the spectrophotometer after 30 min of rest at room temperature and in the dark30. The results were reported based on means of 3 replications with ± standard error.

Reducing power

To measure the reducing power, 2.5 mL of K3Fe(CN)6 10%, 2.5 mL of 0.2 M potassium phosphate with pH 6.6, and 1 ml of different concentrations of the sample were mixed. In addition, 2.5 mL of Trichloroacetic acid (10%) was added to the resulting mixture after its placement in Bain-marie at 50 °C for 25 min and was then centrifuged. Finally, the test solution containing 2.5 ml of distilled water, 0.5 mL of 0.1 ferric chloride, and 2.5 mL of the centrifuged mixture was read by the spectrophotometer. The final report was based on a concentration of the samples and was equal to 0.5 at a wavelength absorption of 700 nm (EC50)31.

Analysis of the essential oil

Gas chromatography–mass spectrometry (GC–MS) analyses were performed with an Agilent 7890 A gas chromatograph equipped with a HP-5MS fused silica column (30 m × 250 μm coated with 5% phenyl methyl silicone, 95% dimethylpolysiloxane, 0.25 μm film thickness), interfaced with an Agilent mass selective detector 5975C inter MSD. Oven temperature was programmed to rise from 60 to 220 °C at a rate of 3 °C/min; and the transfer line temperature was 260 °C. The carrier gas was Helium with a flow rate of 30.6 cm/s and a split ratio of 1:100. Mass scan range from 40 to 300 a.m.u. The identification of oil components was assigned according to the retention index of each GC–MS compound. For measuring the content of each compound, the spectra were considered to be simultaneous if the similarity index was high (%95%). The Kovats index (KI) was also calculated according to compatible Kovats indices (KI) on databases (WILEY, NIST, and PHEROBASE).

Microbial strains

Antimicrobial properties were measured on 4 bacteria and 2 fungi (prepared from the Persian Type Culture Collection). The bacteria included Pseudomonas aeruginosa (ATCC 9027), Escherichia coli (ATCC 25922), Bacillus subtilis (ATCC 6633), and Staphylococcus aureus (ATCC 6538). The fungi included Aspergillus flavus (ATCC 15548) and Candida albicans. It is noteworthy that the bacterial strains were stored at − 85 °C (in 20% glycerol).

Agar diffusion test

This procedure was performed according to the instructions of document M51-A2 of the Clinical Laboratory Standard Institute32. Bacterial and fungal suspensions were cultured in Petri dishes containing sterile food media to evaluate the antimicrobial properties based on disk Agar diffusion. It should be noted that Mueller hinton agar (MHA) and Sabouraud dextrose agar (SDA) were used for bacteria and fungi, respectively. The suspensions of microorganisms were prepared freshly and on a daily basis for the tests. The sterile Mueller Hinton Broth (MHB) medium was used for bacteria and Sabouraud dextrose broth (SDB) for fungi. Incubation conditions for bacterial suspension included 37 °C for 24 h. Its concentration was based on 0.5 McFarland Turbidity Standard. In preparing the fungus suspension, it was first prepared in a petri dish containing 7-day SDA Spore Bank. The suspension was then prepared by sterile 0.85% saline solution from the spore bank (~ 105 spore. mL−1), and 15 µL of EOs were added to the blank disks in each petri dish. For the better diffusion of the biomolecules of the essential oil, each petri dish was first exposed to 4 °C for 4 h. The incubation was then performed on bacterial samples at 37 °C for 24 h, yeast samples at 28 °C for 48 h, and fungal samples at 28 °C for 72 h. The antimicrobial activity of each essential oil was determined by a caliper based on inhibition zone (mm) measurements. The results were reported based on means of 3 replications with ± standard error. The controls of ciprofloxacin (5 μg per disc) and genatamicin (10 μg per disc) of Padtanteb, Iran were used for bacteria; and the Voriconazole (1 μg per disc) of RIDACOM brand was used for fungi and yeast.

Minimum inhibitory concentration (MIC), minimum bactericidal and fungicidal concentration test (MBC/MFC)

Owing to the liquid nature and different essential oils, the density of each one was separately determined to calculate and report different concentrations in mg/mL. MIC and MBC/MFC were estimated according to the broth microdilution method33. The content of each well of the 96-well micro-plates included 100 µL of different dilutions of essential oils (EOs), 95 µL of liquid culture medium (MHB for bacteria and SDB for fungi), and 5 µL suspension of bacteria and fungi. Different amounts of essential oil concentrations were prepared in 1% dimethyl sulfoxide (DMSO). The incubation conditions were 24 h and 37 °C for bacteria, 48 h and 28 °C for yeast, and 72 h and 28 °C for fungi. The MIC was estimated based on the concentration of essential oil that prevented the growth of bacteria or fungi. In fact, MBC/MFC was the minimum concentration of the essential oil preventing the growth of bacteria or fungi in the re-culture of MIC-related microplate wells in solid culture media.

Scanning electron microscope (SEM)

Overnight cultures of a gram-positive bacterium, B. subtilis, and a gram-negative bacterium, E. coli, were adjusted to McFarland 1 standard and then treated with EO constituents at the determined MIC/4 values. After a 1 h incubation period (37 °C) for each bacterium, the cells were harvested by centrifugation, washed five times PBS, and resuspended in sterile distilled PBS. 20 μL of suspension was spread onto a microscope slide and fixed in 3% glutaraldehyde for 2 h. The samples were then diluted in 20% to 100% dilutions in ethanol alcohol solutions. Finally, the slides containing the sample were coated with gold and evaluated using an LEO electron microscope, model 1455VP (Germany).

Statistical analysis

Principal Component Analysis (PCA) and Agglomerative Hierarchical Clustering (AHC) were analyzed by XLSTAT (2021 version). Statistical analyses were conducted utilizing the SAS version. 9.1. LSD test was used to compare the means.

Ethical approval

The study was approved by the Institutional Ethics Review Board of the Yasouj university (2021-02-02).

Human and animal participants

This article does not contain any studies with human or animal subjects.

Results and discussion

Secondary metabolite and antioxidant activity of fruits

The amount of total phenol and flavonoids in the extracts of the studied species showed a significant difference (Table 1). The antioxidant properties of methanolic extracts were determined using DPPH and reducing power methods. In many studies, the amounts of antioxidant properties were reported based on the pure solid extract obtained by various solvents that were then purified. In the present study, the measurement of antioxidant properties was based on milligrams (mg) of the extracted fruit since the plants were from different species. Thus, different weights of the solid extract were obtained in the extraction of each species based on a certain amount of the fruit. Therefore, determining the level of antioxidant activity solely based on the weight of the pure extract (solid extract) may not be a reasonable basis for determining the level of the antioxidant activity of the fruit in a comparative study. This is so because a particular pure plant extract may have considerable antioxidant levels. However, the weight ratio of the pure extract per mg of fruit and its antioxidant activity based on fruit weight can be lower than those of other plants. Also, since fruits of all species are used in powder as spices and traditional medicines, the comparison of antioxidant activity in fruits compared to pure extract can be more functional in the present study.

Table 1 Variance analysis and mean comparisons of total phenolic compounds (TPC), total flavonoid compounds (TFC) and antioxidant activity of Ferula gummosa (FG), Petroselinum crispum (PC), Anethum graveolens (AG), Cuminum Cyminum (CC), Bunium persicum (BP), Heracleum persicum (HP), Pimpinella anisum (PA), Kelussia odoratissima (KO) and Foeniculum vulgare (FV) fruits methanolic extract.

The amount of TFC and TPC levels in P. crispum, A. graveolens, and K. odoratissima fruits were estimated to be about 4.23 and 3.33 mg/mL, 4.17 and 1.36 mg/mL, and 6.05 and 2.04 mg/mL, respectively (Table 1). The antioxidant power of A. graveolens fruit based on the reducing power method with EC50 of 12.22 mg/mL was intermediate among all species. The antioxidant activity of K. odoratissima fruit was significantly lower than those of P. anisum, H. persicum, C. cyminum, A. graveolens and, F. gummosa fruits. Generally, the fruit of F. gummosa had the highest total phenol content (TPC) (21.23 mg/g) and total flavonoid content (TFC) (6.18 mg/g) and B. persicum had the lowest amount of phenolic and flavonoid compounds among others. Since the highest concentration of B. persicum was consumed to inhibit 50% of DPPH (IC50) and reduce ferric ions (Fe3 +) to ferrous ions (Fe2 +) (reducing power EC50), the fruit of this plant had the lowest antioxidant activity. The antioxidant activity of methanolic extract of F. gummosa was higher than those of all species. In H. persicum, P. anisum, and C. cyminum, TPC were 9.58, 7.07, and 5.41 mg/g and TFC were 2.82, 4.73, and 3.71 mg/g respectively (Table 1). In general, the second highest amount of TPC and TFC was observed in H. persicum and P. anisum fruits, respectively. In the same way, H. persicum had the second highest reducing power (EC50: 5.06 mg/mL). After F. gummosa fruit, P. anisum had the best efficacy in scavenging DPPH radicals with IC50 of 10.84 mg/mL (Table 1). Previous studies have demonstrated that levels of phenolic and flavonoid compounds in a plant species are directly and strongly associated with antioxidant powers such as reducing power and DPPH radical scavenging34,35. Therefore, the higher levels of TPC and TFC in F. gummosa and P. anisum fruit compared to the other plants were the main reasons for its antioxidant power. On the other hand, C. cyminum contains a smaller amount of TPC and TFC than some studied species; however, the fruit of this plant had the highest antioxidant activity against DPPH radicals with IC50 of 9.48 mg/mL (Table 1). These results can be attributed to the differences in the types of flavonoid and phenolic compounds of the different plants (Table 1). Different phenolic and flavonoid compounds have different antioxidant powers36,37.

Chemical composition and density of essential oils

According to Tables 2 and 3, some differences were revealed in the chemical profiles of essential oils extracted from all the fruits of the plants. The highest range of compounds was identified in K. odoratissima (31 compounds), and the smallest one was related to P. anisum (8 compounds). On the other hand, Table 4 shows that the density of each essential oil is different due to its unique chemical profile. The chemical type of F. gummosa, P. crispum and A. graveolens essential oils were mostly monoterpene hydrocarbons and oxygenated monoterpenes (Table 2). Also, the density ratios of these plants were about 0.859, 0.826 and 0.873 mg/µL, respectively (Table 4).

Table 2 Chemical compounds of Ferula gummosa (FG), Petroselinum crispum (PC), Anethum graveolens (AG), Cuminum Cyminum (CC) and Bunium persicum (BP) fruit essential oil.
Table 3 Chemical compounds of Heracleum persicum (HP), Pimpinella anisum (PA), Kelussia odoratissima (KO) and Foeniculum vulgare (FV) fruit essential oil.

The main compounds of A. graveolens essential oil included limonene (monoterpene hydrocarbons), carvone (oxygenated monoterpenoids) and dill apiole (phenylpropanoids) in amounts of 33.12, 54.2 and 4.04%, respectively. α-pinene and β-pinene were the main terpene components identified in the essential oils of the two plants, and they constituted 22.83% and 48.36% of the total essential oil in F. gummosa, and 45.3% and 29.70% of the total essential oil in P. crispum, respectively (Table 2). δ-3-carene was also a member of the monoterpene hydrocarbon family and accounted for 10.69% of F. gummosa essential oil (Table 2). Phenylpropanoids constituted 11.22% of the essential oil of P. crispum. Indeed, myristicin was the most important phenylpropanoid in P. crispum (9.34%) (Table 2). In addition to high amounts of monoterpene hydrocarbons and oxygenated monoterpenoids, C. cyminum and B. persicum also had considerable amounts of aldehyde. The highest aldehyde content in these essential oils belonged to cuminaldehyde and γ-terpinen-7-al, which were about 13.92% and 15.34% in C. cyminum and 41.84% and 3.32% in B. persicum. Also, their main original monoterpene hydrocarbons were p-cymene, β-pinene, and γ-terpinene (Table 2). In comparison to C. cyminum (0.914 mg/µL), the density of B. persicum (0.876 mg/µL) essential oil was smaller (Table 2). sesquiterpene hydrocarbon is another fundamental family of the chemical profile of essential oils. They comprise about 67.2% K. odoratissima essential oil. The major chemical compounds of K. odoratissima essential oil encompassed δ-cadinene, selinene, germacrene D, and γ-elemene. limonene (monoterpene hydrocarbon), and Z-ligustilide (phthalides) were the two other major components constituting 8.74% and 15.18% of the essential oil of this plant, respectively (Table 3). trans-anethole is a phenylpropanoid whose amounts in the fruits of P. anisum and F. vulgare were 75.43 and 63.35% of essential oil, respectively (Table 3). Also, the essential oil of these plants had the highest mass-to-volume ratio (mg/μL) compared to other plants (Table 4). Other important compounds of these plants comprised limonene (10.13%), carvone (9.15%) in P. anisum, and fenchone (16.75%) in F. vulgare essential oils (Table 3). The classes of the compositions of H. persicum essential oil were totally different. In addition, this essential oil has the lowest density (0.766 mg/µL) among other studied species of Apiaceae. Most compounds of H. persicum essential oil belonged to aliphatic ester family. 1-octyl acetate, n-octanol, and octyl-2-methylbutyrate were the main components of H. persicum essential oil, which made up 63.44, 8.38, and 6.03% of the essential oil of the fruit of this plant, respectively (Table 3). In general, the main compounds and the range of essential oils amounts current study are in line with previous evaluating of F. gummosa38, P. crispum39,40, A. graveolens41,42, C. Cyminum18, B. persicum43, H. persicum44, P. anisum45, K. odoratissima4 and F. vulgare24 composition.

Table 4 Density and antioxidant activity of Ferula gummosa (FG), Petroselinum crispum (PC), Anethum graveolens (AG), Cuminum Cyminum (CC), Bunium persicum (BP), Heracleum persicum (HP), Pimpinella anisum (PA), Kelussia odoratissima (KO) and Foeniculum vulgare (FV) fruit essential oil.

Antioxidant activity of fruits essential oil

We presented essential oil antioxidant activity through reducing power and DPPH methods. Through DPPH methods (IC50), the antioxidant potential of F. gummosa and P. crispum essential oils was evaluated at about 40.7 and 22.13 mg/mL (Table 4). Also, the reducing powers (EC50) of F. gummosa and P. crispum essential oils were estimated at about 83.38 and 183.33 mg/mL. Indeed, F. gummosa had a weaker DPPH radical scavenging capacity and stronger reducing power than P. crispum, while both of them contained high amounts of β-pinene and α-pinene (Table 4). Usually, the biological activity of essential oils varies in different plant species. At the same time, the difference in enantiomers content of a compound with different biological activity, such as α-Pinene and β-pinene enantiomers, have distinct properties46. The lowest antioxidant activity of essential oils belonged to F. vulgare. Also, H. persicum essential oil had a weaker antioxidant power than the most studied plants. The lowest consumption of essential oils for inhibiting 50% of DPPH (IC50) is related to B. persicum (Table 4). Therefore, in terms of scavenging DPPH radicals, the B. persicum essential oil (IC50: 10.81 mg/mL) has the highest potential. Afterwards, C. cyminum and P. anisum had a stronger antioxidant activity than others based on the DPPH method. The ranges of IC50, such as 114.87 to 287.56 µg/mL45 and 0.384 to 0.809 mg/mL47 were presented for the DPPH free radical scavenging capacity of the P. anisum essential oil. These findings are inconsistent with that of the present study. In a report released by Das et al.48, the IC50 of P. anisum (Varanasi, India) essential oil was equal to 12.31 µL/mL (for DPPH free radical scavenging), which is consistent with the finding obtained in this study (IC50 = 19.56 mg/mL), due to the range of density of the P. anisum essential oil (0.932 g/mL or mg/µL). The slight difference between them was due to differences in the content of constituents of essential oils. In previous studies, IC50 was reported to range from 7.9 to 13.6 mg/mL for DPPH scavenging capacity of C. cyminum essential oil49 and is in the acceptable range of the current study. It has been reported that among cumin essential oil compounds, the DPPH radical inhibitory power of γ-terpinene (IC50 = 22.73 mg/mL) was higher than those of other compounds such as cuminaldehyde (IC50 = 105.36 mg/mL), β-pinene (IC50 = 135.55 mg/mL) and p-cymene (IC50 = 122.25 mg/mL), and the scavenging power of the C. cyminum essential oil based on IC50 was 35.94 mg/mL50. Therefore, the antioxidant properties of essential oils such as C. cyminum and B. persicum can differ from the main pure compounds due to the interaction between the components and the chemical structure of the essential oil profile. Based on the DPPH method, the antioxidant power of K. odoratissima essential oil was weaker than those of all species (except F. vulgare). However, this plant has the highest reducing power with EC50 of 43.22 mg/Ml, while both of these methods are based on measuring the electron-donating capacity of an antioxidant. Therefore, there should not be huge contradictions in determining the antioxidant potential of several chemical compounds (like essential oil) in these methods. In some cases, the results of the reducing power method have little in consonance with the DPPH method regarding the antioxidant properties of these plants essential oils. For instance, according to the DPPH results, the essential oils of C. cyminum, B. persicum, and P. anisum show a high antioxidant power. However, the results obtained from their reducing power indicate a lower antioxidant power than some plants like F. gummosa and K. odoratissima (Table 4). In the reducing power methods, all reaction solutions were water-soluble. Owing to the nature of essential oils as nonpolar and hydrophobic molecules, the essential oil concentration was prepared in a non-aqueous solvent like methanol. When the reaction solution and different essential oil concentrations were mixed, an emulsion was formed to evaluate the reducing power of essential oils. It has been proven that different types of chemical compositions of essential oils affect the quality and stability of emulsions formed51,52. Therefore, unlike the DPPH method, the reducing power of the essential oil, in addition to its antioxidant power, is influenced by the quality and stability of the emulsion in different plants. This approach makes some misinterpretations in the antioxidant power of essential oils in this method. For this reason, unlike the antioxidant report of the methanol extract of the fruit (Table 1), the antioxidant power results of the essential oil via DPPH and power reduction methods could not be matched. Furthermore, the evaluation of the reducing power of essential oils at concentrations above 200 mg/mL was inestimable due to the instability of the emulsions created in the reaction solution. On the other hand, in the essential oils of A. graveolens, P. anisum, F. vulgare, and H. persicum, the increasing trend of reducing power was insignificant with an increasing concentration; therefore, EC50 (absorption 0.5) was not calculated (Table 4). In the previous study, it has been reported that in a comparison between the essential oils of 6 medicinal plants and their 5 main compounds, exactly the same trends (i.e., stability in absorption along with increasing concentration) were observed in trans-anethole and F. vulgare. Thus, the antioxidant activity did not increase with the increase in the concentration of anethole and F. vulgare essential oil53.

Antimicrobial activity

In the current study, positive controls were Gentamicin and Ciprofloxacin for antibacterial properties and Voriconazole for antifungal properties. As shown in Table 5, P. aeruginosa is the most resistant bacterium to the used chemical antibiotics. Furthermore, P. aeruginosa has more resistance to essential oils than other bacteria. P. aeruginosa species has a high resistance to external factors due to having a lipopolysaccharide and phospholipid membrane54. Also, the lipopolysaccharide membrane plays an essential role in the resistance of this bacterium to hydrophobic antibacterial compounds (e.g., essential oils)55.

Table 5 Antimicrobial activity chemical antibiotics.

In general, there are no detailed reports on the antimicrobial properties of some wild plants, such as F. gummosa, H. persicum and K. odoratissima. However, in Table 6, a brief report of the range of antimicrobial properties of the essential oils of Apiaceae family plants was obtained based on MIC and MBC. What is evident in the presented table is that the range of estimated values is vastly contradictory. In fact, based on a general summary of previous research results, it is impossible to distinguish the antimicrobial power of each plant's essential oil, either separately or comparatively within the family. In previous studies, the range of MIC and MBC values exists as values with units of µg/mL and mg/mL in some species (Table 6). Although the results can be different due to the difference in the bacterial strains and the chemical profile of the plants, many of the observed differences are beyond these factors. Due to the liquid nature of essential oils, it is common to report MIC and MBC values as µL/mL. However, special attention should be paid to the density of essential oils to express MIC and MBC based on µg/mL and mg/mL. As presented in Table 4, each essential oil has a different density according to the its constituent compounds. For example, according to Table 7, the MIC value in the essential oil of H. Persicum and F. vulgare plants against A. flavus was estimated as 1 mg/ML. According to the density of the essential oil of H. Persicum (0.766 mg/µL) and F. vulgare (0.990 mg/µL), this value is equivalent to 1.30 and 1.01 µL/mL, respectively. The lack of attention to the density and the approach expressed regarding unit conversion are the main reasons for the difference in various reports. In addition, a numerical range regarding the antimicrobial properties of the main purified compounds related to the essential oils of this family, such as α-pinene, limonene, carvone, trans-anethole and ligustilide has also been presented in Table 6. This can provide an overview of the range of antimicrobial effects of the main essential oil compounds of this family in the form of MIC and MBC values. In this study, according to the results of a wide range of effective concentrations in previous studies, the range of concentrations from 0.005 mg/mL (equivalent to 5 µg/mL) to 96 mg/mL was tested to enable the conduct of accurate evaluation and measurements. According to Tables 7 and 8A. graveolens, K. odoratissima and H. persicum have acceptable antibacterial properties. Regarding the results of MIC, MBC, and inhibition zone, the antibacterial effect of K. odoratissima, particularly against B. subtilis and E. coli is substantial. Nevertheless, despite the positive antibacterial response of H. persicum and A. graveolens essential oil, their antibacterial properties compared to other studied plants have intermediate positions. Against similar studied bacteria, Ruangamnart et al.26 reported that the range of MIC of A. graveolens essential oil (containing 44.61% limonene, 28.02% carvone and 19.98% dillapiole) was higher than 10 mg/Ml. While the reported chemical profile was partly different from that of the current study, especially in terms of dillapiole content. limonene constitutes 33.12% of the essential oil of A. graveolens and 8.74% of the essential oil of K. odoratissima. However, the antibacterial properties of this compound through the MIC method were reported to be approximately 16 µL/mL against E. coli72. In another study, the reported MIC values of 26.25 and 52.5 mg/mL for E. coli and S. aureus were claimed, respectively73. Carvone was another main A. graveolens essential oil compound (54.3%) and its MIC against S. aureus (ATCC 6538) and B. subtilis (ATCC 6633) with the same strains of the current study were evaluated to be about 5 and 10 mg/mL22.

Table 6 Summary of the MIC and MBC range of some Apiaceae essential oils and their main compounds against some pathogenic microorganisms in previous studies.

The current result also demonstrates that the inhibitory and bactericidal power of A. graveolens essential oil against S. aureus is better than B. subtilis; and the range of MIC and MBC has consonance with the biological activity of its main compounds. In the case of limonene, in addition to the antibacterial effects of the compound itself, high synergistic effects were observed between the compound and chemical antibiotic compounds against several bacterial species such as E. coli and S. aureus74. The presence of limonene, along with other main components of K. odoratissima essential oil, may increase the synergetic effects and consequently increase the antibacterial properties. As it was mentioned above, there are no specific reports about the antibacterial properties of K. odoratissima and H. persicum; however, the biological properties of any essential oil are directly related to its components. For instance, K. odoratissima contains a high amount of Z-ligustilide and its MICs were evaluated to be 1, 1, 8 and 10 mg/mL against B. subtilis, S. aureus, E. coli and P. aeruginosa, respectively (Table 8). Previously, it has been reported that Z-ligustilide compound purified from Cnidium officinale essential oil had high antibacterial properties, and its MIC against S. aureus was 1 mg/mL70. The amount of MIC in H. persicum essential oil was estimated to be 8, 10, 8, and 24 mg/mL against B. subtilis, S. aureus, E. coli and P. aeruginosa, respectively (Table 8). Its essential oil consists of high amount of acetate composition and octanol. These compounds such as octyl acetate and octanol were factors of antibacterial properties in the Zosima absinthifolia plant75. In the past, the concentrations of MIC of Heracleum rigens essential oil (containing 51.2% bornyl acetate, 9.62% limonene, and 3.94% octyl acetate) for S. aureus, Bacillus subtillis, E. coli, and P. aeruginosa were estimated to be 4.5, 2.25, 36, and 36 mg/mL, respectively. Based on MIC and MBC results, P. crispum, P. anisum, and F. vulgare essential oil up to 96 mg/mL had no antibacterial power on the tested bacteria (Tables 8). Furthermore, these essential oils inhibited the lowest zone (Table 7). Even the essential oil of P. crispum and P. anisum against B. subtilis bacterium did not form any inhibition zone. Approximately 75 and 63% of the essential oil of P. anisum and F. vulgare were trans-anethole. In the past, it has been determined that trans-anethole isolated from the Croton zehntneri plant did not have antibacterial properties against S. aureus, P. aeruginosa, and E. coli bacteria76. Also, the previous evaluation of trans-anethole activity against S. aureus indicated that this compound had a weak antimicrobial potential68. Presumably, trans-anethole has a non-toxic effect on bacteria so it can be the reason for the lack of antibacterial responses in the P. anisum and F. vulgare essential oil. Similar to the results of the current study, Fitsiou et al.77 studied the essential oils of several medicinal plants and found that the P. anisum essential oil did not have any antibacterial activity. On the contrary, there were reports of antibacterial properties of the essential oil of this plant on the same bacteria, MIC in the range of 0.06 to 0.25 µg/mL78 and 2.0 to 3.0 µg/mL60 were reported, which were inconsistent with the results of this study.

Table 7 Antimicrobial activity (inhibition zone) of Ferula gummosa (FG), Petroselinum crispum (PC), Anethum graveolens (AG), Cuminum Cyminum (CC), Bunium persicum (BP), Heracleum persicum (HP), Pimpinella anisum (PA), Kelussia odoratissima (KO) and Foeniculum vulgare (FV) fruits essential oil.

According to the results of Table 2, the predominant compounds of F. gummosa and P. crispum essential oils are two common compounds of α-pinene and β-pinene. However, there are significant differences between antibacterial properties of these essential oils in terms of MIC, MBC and inhibition zone traits. The weak antibacterial properties of P. crispum essential oil have previously been revealed79. In contrast, the amount of inhibition zone of F. gummosa essential oils indicated strong antibacterial properties. Similarly, these estimations indicate its MIC and MBC value against B. subtilis, S. aureus, and E. coli was 2 mg/mL. The possibility of differences in the enantiomeric profile of the main components of the essential oils of these plants is a major reason for this difference. It was reported that differences in levels of enantiomeric compositions of α-pinene ((1R)-(+)-α-pinene and (1S)-(−)-α-Pinene) in Juniperus communis essential oil could increase its lethality effects by 4 to 16 times. It has also been reported that (−)-α-pinene, as an enantiomeric composition of α-Pinene, has had no antibacterial or antifungal properties; however, (+)-α-pinene has had strong antimicrobial properties80. Therefore, the amounts of pinene detected in the essential oils of the two different plants, including F. gummosa and P. crispum, might have different chemical reactions due to the differences in the enantiomeric compositions of Pinene. Also, the most enantiomeric composition of α-pinene in F. gummosa essential oil has been reported to be (1R)-(+)-α-pinene38.

Table 8 Antimicrobial activity (MIC and MBC) of Ferula gummosa (FG), Petroselinum crispum (PC), Anethum graveolens (AG), Cuminum Cyminum (CC), Bunium persicum (BP), Heracleum persicum (HP), Pimpinella anisum (PA), Kelussia odoratissima (KO) and Foeniculum vulgare (FV) fruits essential oil.

Interactions and combined effects of all the essential oil components could be another reason for the differences in the properties of F. gummosa and P. crispum essential oils. This is due to the fact that the final properties of each essential oil are affected by all essential oil compounds and their antagonist and synergetic effects. For example, it is reported that (−)-α-pinene alone does not have high antibacterial properties; however, its use in combination with chemical antibiotics increases their antibacterial power by several hundred times65. Therefore, the interaction of the chemical components of the essential oil leads to differences in their properties due to the different chemical profile of each species. The δ-3-Carene compound is an active compound against bacteria with a MIC of about 0.16 g/L (equivalent to 0.16 mg/mL) against E. coli, B. cereus, and S. aureus69. It has also been determined that the reason for the high antibacterial power of 3-Carene compound can be attributed to the multiple mechanisms of its destructive effects on the morphological, physiological and molecular levels of bacteria. At the morphological level, this compound causes changes in the natural shape and structure of bacteria, disruption and destruction of bacterial membrane integrity, and leakage of intracellular substances. At the physiological levels, it reduces the metabolic function by impairing the ATP breakdown, intracellular Tricarboxylic acid cycle, and Glycolysis-related enzymes. At the molecular level, this compound binds to DNA, disrupts gene expression, and as a result, increases the antibacterial effect81. Therefore, the presence of this compound in F. gummosa essential oil, along with other essential oil components, can be a significant reason for the higher potency of its antibacterial properties than P. crispum essential oil (Tables 7, 8). The Antibacterial properties of C. cyminum and B. persicum were higher than others in all bacterial species based on the inhibition zone, MIC, and MBC. The concentration ranges of MIC, MBC, and inhibition zone for C. cyminum essential oil were 0.5–4, 1–10 mg/mL, and 24.88–37.23 mm and for B. persicum were 0.5–4, 1–6 mg/mL, and 25.23–39.33 mm, respectively (Tables 7, 8). In relation to the antibacterial properties of C. cyminum and B. persicum essential oils in the current study, the range of MIC and MBC had consonance with previous research findings. However, the expected difference could be attributed to differences in essential oils compounds and differences in the resistance of bacterial strains. Among the studied plants, H. persicum, C. cyminum and B. persicum essential oils had high antifungal activity against both A. flavus and C. albicans. In contrast, the essential oils of A. graveolens, P. anisum and F. vulgare had the lowest antifungal activity against C. albicans (Tables 7, 8). The essential oils of F. gummosa and P. crispum showed inhibition against A. flavus in 4 mg/mL concentration, but it seems that they had the lowest fungicidal properties against this fungus. Generally, it has been founded that since plant essential oils have hydrophobic components, they degrade and disrupt the membrane fluidity, change ergosterol biosynthesis, increase leakage, and change pH in the cytosolic space. Thus, they cause damage to cell organs and lead to the death of fungal and yeast cells5. In addition, the results showed that different compounds of the essential oils of each plant have different effects on different microorganisms. For example, it has been reported that α-pinene has no strong fungicidal effect against A. flavus82. According to the results of one study on the fungicidal effect of α-Pinene, the estimated amount of the MFC for A. flavus for this compound was approximately 15 times higher than the MFC for C. albicans80. Considering the amount of pinene in the essential oils of P. crispum and F. gummosa, it can be a reason for the superiority of these essential oils in the inhibition and lethality power of C. albicans compared to A. flavus. According to another report, cuminaldehyde has a high potency in damaging fungal hyphae, reducing ergosterol biosynthesis, growth inhibition, and aflatoxin production in A. flavus, as the MIC of fungal growth in this compound is about 0.64 μL/mL (equal to 0.62 mg/mL considering a density of about 0.97 g/mL for this compound)83. Previously, Das et al.48 have indicated the power of P. anisum essential oil in inhibiting A. flavus. anethole has been the main component of a P. anisum and F. vulgare essential oil. It has also been previously reported that the antifungal properties of Illicium verum and F. vulgare essential oils can be attributed to the high amount of anethole84,85.

Related research findings have revealed that the membranes and walls of bacteria are destroyed by plant essential oils. However, other destructive effects such as cytoplasmic coagulation, cytoplasmic content leakage, and protein damage have also been reported against bacteria with essential oil application86. One of the reasons for the destruction of the bacterial membrane is the hydrophobic nature of essential oils, which disrupts the two phospholipid layers of the bacterial membrane, although the destructive effects of each essential oil on bacteria may be different depending on the type and composition of each essential oil6. All the SEM images show the effect of membrane degradation on both bacterial species with the application of H. persicum and C. cyminum essential oils (Figs. 1, 2). In SEM images, some membrane destruction can be detected in both species of bacteria with the use of H. persicum essential oil, but the SEM imaging of C. cyminum show that the power of C. cyminum essential oil in membrane destruction is much more intense. In addition, the effects of C. cyminum on the destruction of each of the two bacteria, especially B. subtili, can be regarded as a mass of damaged cell membranes. Some bacterial cells with disrupted membranes can be observed in A. graveolens treatments. But the effects of A. graveolens essential oil were lower than those of H. persicum and C. cyminum. In addition, no membrane degradation by P. crispum essential oil was not observed. Nevertheless, the essential oil of this plant caused morphological changes in bacteria, including changes in the size of the bacteria and shrinkage of the bacterial membrane, especially E. coli. (Figs. 1, 2).

Figure 1
figure 1

Scanning electron micrographs of the effect of Cuminum cyminum (a), Heracleum persicum (b), Anethum graveolens (c) and Petroselinum crispum (d) essential oils on B. subtili.

Figure 2
figure 2

Scanning electron micrographs of the effect of Cuminum Cyminum (a), Heracleum persicum (b), Anethum graveolens (c) and Petroselinum crispum (d) essential oils on E. coli.

Principal component analysis (PCA) and agglomerative hierarchical clustering (AHC) on essential oils’ profiles and properties

AHC analysis is a statistical approach to categorizing the tested factors based on their characteristics. In the current research, this analysis was done three times based on essential oil compounds data (Fig. 3a), antimicrobial properties data (Fig. 3b), and data related to essential oil compounds and antimicrobial properties together (Fig. 3c). As presented in Fig. 3a, the essential oils of KO and HP plants are more distinct from other plants in terms of chemical profile, and both of them form a cluster independently. FG and PC with a high level of pinene and BP and CC with a high level of Aldehyde compounds form the other two clusters. Finally, PA and FV with high similarity to each other and together with AG formed a cluster (Fig. 3a). In Fig. 3b, it has been shown that there are three classifications based on antimicrobial properties. The cluster has a group with high antimicrobial properties, including CC, BP, FG and KO; a group with moderate antimicrobial properties, including HP and AG; and a group with weak antimicrobial properties, including FV, PA and PC. Clearly, the best and most accurate classification was obtained based on the analysis of all data related to essential oil compounds and biological activities (Fig. 3c).

Figure 3
figure 3

(a) Analysis based on features related to the chemical profile, (b) Analysis based on features related to the antimicrobial properties, (c) Analysis based on features related to the antimicrobial properties and chemical profile. (d) Principal Component Analysis based on antimicrobial activity of studied plants and their main compound groups. In figures, HP: Heracleum persicum, KO: Kelussia odoratissima, AG: Anethum graveolens, FV: Foeniculum vulgare, PA: Pimpinella anisum, CC: Cuminum Cyminum, BP: Bunium persicum, FG: Ferula gummosa and PC: Petroselinum crispum. BMIC, BMBC and BIZ: MIC, MBC and inhibition zone of B. subtilis, EMIC, EMBC and EIZ: MIC, MBC and inhibition zone of E. coli, BMIC, BMBC and BIZ: MIC, MBC and inhibition zone of B. subtilis, SMIC, SMBC and SIZ: MIC, MBC and inhibition zone of S. aureus, PMIC, PMBC and PIZ: MIC, MBC and inhibition zone of P. aeruginosa, CMIC, CMFC and CIZ: MIC, MFC and inhibition zone of C. albicans and AMIC, AMFC and AIZ: MIC, MFC and inhibition zone of A. flavus.

As a result, the differences between the plants have been determined more precisely in Fig. 3c based on essential oil compounds and biological activities and, thereby, a more complete and accurate view has been presented. The identification of compounds or common tests to evaluate the biological properties of essential oils is a costly, time-consuming, and sometimes dangerous process. On this basis, there is the possibility that a mathematical prediction model in the future can predict the properties of essential oils based on their constituent compounds. Since there are valuable differences and similarities in the grouping of Fig. 3a, b with Fig. 3c. for example, in all three modes of analysis, CC and BP have been placed in one group, and PA and FV have been placed in another group (Fig. 3a–c). This shows that it is possible to determine the potential of antimicrobial properties using modeling in advance, particularly in intra-familial classes based on the chemical composition of these plants. Reyes-Jurado et al.87 clearly state that common laboratory tests regarding the antimicrobial properties of essential oils are influenced by many variables such as human errors, source of essential oil, extraction method, amount of chemical compounds, the volume of inoculation, growth phase of microorganism, the strain of microorganisms to be tested, the culture medium used, the incubation time, the temperature, and the pH of water. Finally, according to the impact of the mentioned cases, there are distinctions and differences in the reports, even about a plant species (Table 6). As a result, as in the present study, one of the most effective methods to reduce these differences is to examine several plant species simultaneously and use techniques such as SEM to confirm the results of tests that rely solely on visual observations (such as MIC, MBC and inhibition zone). In this regard, the factors causing differences during the tests become uniform and, thereby, accurate comparisons are made with less error. However, the use of a mathematical model seems to be a more appropriate solution. Today, the mathematical modeling of chemical antibiotics is done in very advanced fields to check the effective dose even in the human body88. At the same time, there is no report regarding the application of these methods to investigate plant compounds' properties, even in vitro. In the present study, advanced modeling was not applicable due to the small number of tested samples. Therefore, the extent and number of laboratory samples in the future should be considered. Also, the placement of KO and HP in separate groupings in different states in AHC analysis shows that plant species should be selected based on similarities in chemical profiles to design the prediction model in the future. From each chemical profile with similar compositions, there should be several plant species for the conduct of a more accurate modeling. Accordingly, very contradictory chemical profiles should be removed from the prediction model. Another apparent contradiction is the position of FG in these groupings. Indeed, in Fig. 3a, the high level of pinene placed this plant in the same group as PC plant. While in their grouping based on antibacterial properties, the position of these two plants differed (Fig. 3a). In Fig. 3c, the FG plant was also placed in a group alone. It was pointed out that the enantiomers of a compound such as α-pinene have different biological properties. This trend shows that the difference in biological properties of enantiomeric composition can cause errors in the model.

Therefore, for the accomplishment of more accurate modeling, the enantiomeric composition levels of the main components of the essential oil should be well estimated. Based on the PCA analysis (Fig. 3d) and the location of independent variables (active observation), a high level of similarity between the essential oils of CC and BC plants and FG and PA plants was also observed. In the PCA diagram, the more acute the angle between the two factors, the stronger and more positive the correlation between the two factors. In this regard, the aldehyde compounds in the Apiaceae family can be considered potent bactericidal agents. Additionally, according to the distance of the traits indicating antibacterial properties (active variables) in different plants (active observation), the essential oils of two plants CC and BP (containing a high level of aldehyde compounds), had high antibacterial properties. Furthermore, the ester and alcohol compounds of plant essential oils are highly effective in killing fungi. HP plant essential oil shows the most fungicidal effect on these microorganisms.

Conclusion

Overall, the results of measuring the methanol extract of fruits by DPPH and reducing power methods showed that the fruits of F. gummosa, C. cyminum and P. anisum plants had higher antioxidant properties than other studied plants. According to the DPPH method, the antioxidant properties of essential oil C. cyminum and B. persicum were stronger than others. Considering previous studies and a comparative approach to investigate the properties of 9 plants of the Apiaceae family, the present study determined the potential of these essential oils in terms of antimicrobial properties and effective doses. In this regard, the essential oils of 3 plants, K. odoratissima, F. gummosa, and H. persicum, with their specific chemical profiles, had a favorable potential in inhibiting the growth of various bacteria. The essential oil of H. persicum was more effective than other plants in inhibiting C. albicans and A. flavus. Thus, the essential oil of this plant can be a suitable option for medical applications to inhibit this human pathogen C. albicans and to inhibit fungal food contamination and storage loss of A. flavus. The essential oils of C. cyminum and B. persicum plants had very high antibacterial properties compared to others, hence, both plants are more emphatically recommended for industrial and medical applications. Also, the results of AHC and PCA analysis showed that the creation of mathematical models for predicting the biological properties of essential oils based on their constituents can be effective and it is suggested to pay more attention to this finding in the future.