Natural diversity in phenolic components and antioxidant properties of oregano (Origanum vulgare L.) accessions, grown under the same conditions

Oregano (Origanum vulgare L.) is a rich source of biologically active components such as phenolic compounds. Here, seven pot grown O. vulgare accessions belonging to three subspecies (subsp. virens, subsp. vulgare and subsp. gracile) were investigated for their content in sixteen bioactive phenolic compounds as well as their antioxidant capacities (DPPH• and FRAP tests), total phenolic content (TPC) and total flavonoid content (TFC) in order to identify the most suitable ones on an industrial level. HPLC analyses showed that rosmarinic acid (659.6–1646.9 mg/100 g DW) was by far the most abundant constituent, followed by luteolin (46.5–345.4 mg/100 g DW), chicoric acid (36.3–212.5 mg/100 g DW), coumarin (65.7–193.9 mg/100 g DW) and quercetin (10.6–106.1 mg/100 g DW), with variability in concentration depending on the accession and subspecies. The highest level of rosmarinic acid and TPC was obtained from Ardabil accession (subsp. virens). There was a significant and positive correlation between rosmarinic acid and antioxidant activity (r = 0.46). TFC significantly correlated to TPC (r = 0.57) as well as to chicoric acid (r = 0.73). Cluster (CA) and principal component (PCA) analyses classified the investigated accessions in three different groups. Such natural variabilities in phenolics provide the possibility of using elite plants for nutraceutical and pharmaceutical industries and domestication of highly antioxidative accessions of oregano.

www.nature.com/scientificreports/ Rosmarinic acid, luteolin, quercetin, apigenin, scutellarein and their derivatives are the major phenolic acids and flavonoids that have been detected in oregano species 10 .
Through inhibiting the initiation of oxidizing chain reaction, antioxidants inhibit or delay the oxidation of molecules. Because of its high phenolic content, oregano is considered as one of the most favorite natural antioxidants 11 . The antioxidant capacity of O. vulgare extract is mainly related to its phenolic constituents 8,12 . The phenolic antioxidants take over various pharmacological properties, such as antidiabetic, antiulcer, antiviral, cytotoxic, antitumor and anti-inflammatory activities 13 , and are mainly responsible for the health effects of O. vulgare. Potential anticancer characteristics of phenolic acids and flavonoids have been previously reported by Shukla and Gupta (2010) 14 .
Our previous investigations showed that there is a great variability in terms of essential oil compounds among the different subspecies and/or accessions of oregano herb 5,15 . However, there is also a big challenge as to whether there is variation in phenolics and antioxidant capacity among the accessions and subspecies of oregano. A previous study has shown a high variability in rosmarinic acid content and oxygen radical absorbance capacity (ORAC) and total phenolic content (TPC) of European oregano accessions 11 . As a result, it can be presumed that such variability in rosmarinic acid and other phenolic compounds can also be seen in Iranian oregano accessions. Considering the great importance of phenolic compounds in human nutrition and their health effects, the current study aimed to investigate the variability in main phenolic compounds (including rosmarinic acid, chicoric acid, apigenin, luteolin and others), total phenolic content (TPC), total flavonoids contents (TFC) and antioxidant capacity among different accessions of Iranian oregano belonging to three subspecies (subsp. virens, subsp. vulgare and subsp. gracile). This research is a continuation of an ongoing oregano breeding program, pointing out the elite accessions of oregano in terms of phenolic components for domestication. In the meantime, the results generated from the current study will be useful to introduce the elite accessions with high antioxidant capacity into food and pharmaceutical industries.
Plant and soil materials. Seeds of the seven accessions of O. vulgare belonging to three subspecies (subsp. virens, subsp. vulgare and subsp. gracile) were obtained from the seed gene bank of the Research Institute of Forest and Rangeland in Tehran, Iran. The seed gene bank declared that seeds of all accessions were obtained under national and international guidelines and the seed were prepared under the supervision and permission of Maragheh University and all authors comply with all the local and national guidelines. The voucher specimens of the plants were deposited at the herbarium of Department of Horticultural Sciences, University of Tehran, Karaj, Iran. Geological characteristics of the seed, collection sites, subspecies and their voucher numbers are presented in Table 1.
Oregano seeds were sown in a plastic germination tray filled with coco peat:perlite mixture (70:30, w:w) and kept in a glass greenhouse at the University of Maragheh, Maragheh, Iran. On 12 April 2020, the seedlings with about 10 cm height were transferred into 7 L pots. The growth medium was composed of combined proportions of field soil, silt, leaf mold and perlite (45:25:20:10, v:v). The soil contained 0.09% N, 21.8 mg kg −1 available P, and 487 mg kg −1 available K; the medium pH was 7.2 and EC was adjusted to be 1.2 dS m −1 . During the growth seasons, the plants were irrigated regularly. The pots were placed in the greenhouse under natural daylight with a maximum and minimum day temperature of 33.5 and 17.5 °C, respectively. To warrant the well growth of the plants, the pots were fed with half Hoagland-based solution six times during the growth period (Jons, 2014). The harvest was performed at full flowering stage on 15 August 2020. The plants were cut from 5 cm above the soil and dried in oven at 40 °C, after which their dry drug weight (flowering aerial part) was measured (g/pot). For each accession, 12 individual plants were grown and each of the four harvested oregano plants was bulked together to obtain three replications (n = 3) for the extraction purpose.
Plant extraction procedure. For extraction, 200 mg of leaf and inflorescence powder from each dried sample was dissolved in 20 mL of 80% methanol (methanol-water mixture in 80:20 proportion) and shaken for www.nature.com/scientificreports/ 72 h at room temperature (25 °C). Then, the supernatant was filtered using a Whatman filter paper (No. 4) and the residue was re-extracted with the same method 16 .
Determination of total phenolic compounds. In order to determine the total phenolic content of oregano accessions, the Folin-Ciocalteu method was used (Spanos and Wrolstad 1990). For this purpose, 10 μL of plant extract were added to 500 μL of 10% Folin-Ciocalteu's reagent. Next, 500 μL of 1% saturated sodium carbonate solution were added to the mixture and incubated for 2 h at 25 °C. After incubation, the solution absorbance was read at 765 nm using a microplate reader (Spectromax-M5-USA). Gallic acid with different concentrations was used to make the calibration curve. The total phenolics content was reported as mg gallic acid equivalent (GAE) g −1 dry weight 17 .
Determination of total flavonoids. According to the outline described by Quettier-Deleu et al. (2000) 18 , total flavonoids content of oregano extracts was measured using AlCl 3 reagent. Briefly, 2 mL of 2% methanolic aluminum chloride solution was added to 2 mL of the extract and then stirred slightly. The procedure was continued by adding 6 mL of 5% potassium acetate solution. The mixture was incubated for 40 min at 25 °C and after that the absorbance of samples was read at 415 nm using a microplate reader (Spectromax-M5-USA). For quantification of absorbance, quercetin was used as a standard. The amount of flavonoids in each extract was expressed in terms of mg quercetin equivalent (QE) g −1 dry weight.

RP-HPLC analysis of phenolic compounds.
To identify the phenolic components, oregano extracts were analyzed using high-performance liquid chromatography (HPLC, Shimadzu, Japan) coupled with a photodiode array detector (SPD-M20A, Shimadzu, Japan). Phenolics were separated using a phenomenex column (Gemini-5 µm NX-C18 110 Å, LC Column 250 × 4.6 mm). The mobile phase was composed of solvent B (20 mM phosphoric acid) and solvent D (methanol). The program began with isocratic elution with 10% B (1.5 min), then a linear gradient was used until 40 min, increasing B to 95%, 5 min in 95% B, decreased to 10% in 30 s, 6 min in 10%. The flow rate was 1 mL min −1 . Ten μL from each sample were injected into instrument. Scanning was performed from wavelength range of 200 to 400 nm. The identification of phenolic components was performed by comparing their pure standard retention time and the UV spectra with those of standards. For quantification of the components, different concentrations of each standard (10-80 μg mL −1 ) were injected under the same operative conditions of extracts and the calibration curve of each compound was plotted. The results were exposed as mg of each compound per 100 g of dry weight (mg/100 g DW).
DPPH-radical scavenging activity. The  Ferric-reducing antioxidant power (FRAP) assay. Ferric-reducing antioxidant power (FRAP) of the extracts was determined using the method explained by Benzie and Strain (1996) 19 with minor modifications. Following the procedure, the FRAP reagent was prepared by mixing 100 mL of acetate buffer (300 mM; pH 3.6), 10 mL of TPTZ (10 mM in HCl), and 10 mL of FeCl 3 .6H 2 O (20 mM in water). The fresh mixture was incubated at 37 °C until use. Acetate buffer was prepared by mixing 0.31 g of sodium acetate trihydrate with 1.6 mL of glacial acetic acid and reaching the volume of 100 mL using doubly distilled water. The reaction was done by mixing 380 μL of FRAP reagent and 20 μL of plant extracts in the wells and incubating at 37 °C for 10 min. The absorbance of samples was read at 593 nm using Spectromax-M5 microplate reader (USA). Quantification was done when the calibration curve of Fe 2 SO 4 .7H 2 O (25-1000 µM) was plotted. The FRAP antioxidant activity was expressed in terms of μM Fe (II).

Statistical analysis.
Analysis of variance (ANOVA) and least significant difference (LSD) test at a 5% probability level was performed using SAS statistical software (version 9.1, SAS Inst., USA). Pearson correlation coefficients were computed using IBM SPSS (SPSS, version 22, USA). Cluster analysis and principal component analysis (PCA) were performed using Xlstat software, 2018.

Results
Drug yield. The results from ANOVA showed that there were significant differences (p < 0.01) among oregano accessions in terms of drug yield (flowering aerial parts). As depicted in Fig. 1, a high variability was observed among the oregano accessions in terms of drug yield, ranging from 32.3 to 50.5 g/pot. According to mean compression values, the highest drug yield (50.5 g/pot) was obtained for Gilan accession (subsp. virens), followed by Kaleybar (43.2 g/pot; subsp. vulgare) and Ardabil (37.5 g/pot; subsp. virens) accessions. However, Arasbaran (32.3 g/pot; subsp. vulgare), Baneh (33.7 g/pot; subsp. gracile) and Mazandaran (33.8 g/pot; subsp. virens) accessions yielded the lowest drug weights in comparison with others. www.nature.com/scientificreports/ Total polyphenol and flavonoid contents. As depicted in Fig. 1, based on the mean compression of TPC in oregano accessions, a high variability was observed. According to the results, total polyphenols contents varied from 20.3 to 35.5 mg GAE g −1 DW. The highest and lowest TPC were observed in Mazandaran (subsp. virens) and Baneh (subsp. gracile) accessions, respectively. This creates a variance of roughly 44% between the lowest and highest concentration levels of TPC. However, there was no significant difference among the Ardabil, Arasbaran, Gilan, Namin and Kaleybar accessions in terms of TPC (Fig. 1). According to the ANOVA results, TFC were largely variable among the investigated accessions (p < 0.01). Overall, Ardabil and Arasbaran accessions had the highest and lowest TFC (22.2 and 7.7 mg QE g −1 DW, respectively) ( Fig. 1). This creates about a three-fold variation among the mentioned accessions.
Antioxidant properties. Figure 1 shows the free radical-scavenging activities of the seven oregano accessions which were assessed using the DPPH method. According to the obtained results, there was not a high difference among the oregano accessions in terms of DPPH antioxidant activity. However, the maximum antioxidant capacity belonged to Baneh population (IC 50 = 91.3 µg mL −1 ). It is notable to mention that the IC 50 content (50% reduction in DPPH concentration) for ascorbic acid as the positive control was 9.95 µg mL −1 . On the other hand, no significant difference was found among the oregano accessions in terms of DPPH antioxidant capacity. According to the mean compression, the FRAP activity of oregano accessions ranged from 347.7 to 493.6 (µmol g −1 DW). Overall, Gilan accession (subsp. virens), with 493.6 µmol g −1 DW, showed the highest FRAP antioxidant activity; however, there were no statistically significant differences between other investigated accessions (Fig. 1). The FRAP activity of ascorbic acid as the positive control was 1946.6 µmol g −1 .
Multivariate analysis of phenolic antioxidants. Cluster analysis (CA) and principal components analysis (PCA) were performed using a data matrix composed of 133 data points (19 variables × 7 observations) to display the grouping among oregano accessions based on their phenolic components and antioxidant attributes (Fig. 2). The first two components explained 69.03% of the total variance. According to the CA and PCA analyses, the mentioned accessions were classified into three different groups (Fig. 2). The first cluster was composed of three accessions (namely, Mazandaran (subsp. virens), Gilan (subsp. virens) and Baneh (subsp. gracile) including the lowest rosmarinic acid and luteolin content, and DPPH radical scavenging activity, and higher amounts of apigenin. Furthermore, Kaleybar accession is in close angle with kaempferol and placed in a distinct group (Fig. 2). However, Arasbaran (subsp. vulgare), Ardabil (subsp. virens) and Namin (subsp. virens) accessions possessed the highest rosmarinic acid, luteolin, vanillic acid, chicoric acid, coumarin and DPPH radical scavenging activity and were placed in the same group. Simple correlation analysis showed a significant relationship between the TFC and TPC content (r = 0.58, p < 0.01), as well as to chicoric acid (r = 0.73, p < 0.01). Furthermore, correlation analysis revealed that rosmarinic acid content is moderately correlated with DPPH antioxidant activity (r = 0.46, p < 0.05). Moreover, there were negative correlations between vanillic acid and caffeic acid content (r = −0.92, p < 0.01), and syringic acid (r = −0.93, p < 0.01). Statistical analysis also revealed moderate and significant relationships between the chlorogenic acid and caffeic acid (r = 0.67), syringic acid and benzoic acid (r = 0.70), syringic acid and naringenin (r = 0.67) ( Table 3).

Discussion
O. vulgare is a popular medicinal herb, which is frequently used as antioxidant, antifungal, antimicrobial, expectorant, stimulant, carminative, anti-cancer and ant-aging agent among other applications in the food and beverage industries 7,9 . Natural antioxidants isolated from oregano extracts have been found to have many health effects such as mental well-being and decreasing postsurgical wounds infection 6,9 . The quality of oregano plants based on antioxidant efficacy is significantly influenced by subspecies and geographical origin. Such variation could be due to the plant genetics and environmental conditions, such as light intensity, soil conditions and water availability 20 . The results of the present work showed a high variation in drug yield of Iranian oregano accessions. Concerning the planting in the same controlled environmental and growing conditions, such differences in drug yield may be related to the genetic factors of oregano subspecies 15 . Genetic breeding in oregano subspecies can potentially aid in achieving superior agronomic properties, higher yields and desirable cultivars rich in phenolic compositions 5,11 . Therefore, accessions with higher biomass yield can be considered for higher production of secondary metabolites in breeding programs. O. vulgare is a rich source of phenolic acids and flavonoids with a strong antioxidant activity, which qualify the plant as a potential functional food and medicine 10,11,21 . Polyphenols can fluctuate broadly in their arrangement  www.nature.com/scientificreports/ and overall classification; however, all of them share the common structure of comprising a minimum of one aromatic ring and one or more hydroxyl groups 22 . As a matter of fact, phenolic compounds are considered as important ingredients of plants, fruits and vegetables 23 . Here, a high variation in TPC and TFC was observed among the investigated accessions belonging to three different subspecies. The highest TPC was found in Mazandaran (subsp. virens) accession, which was about 42% higher than that of Baneh (subsp. gracile) accession. Regarding flavonoids, the highest concentration was obtained from Ardabil (subsp. virens), while the lowest was seen in Arasbaran (subsp. vulgare) accession. Among these accessions, the highest TFC was about three times higher than the lowest one. In agreement with the present study, a comparative study on European oregano demonstrated that there is a high variation among the oregano accessions in terms of TPC and antioxidant activity 11 . Compared with our study, Dambolena et al. (2010) 24 found a smaller variation in the TPC of oregano subspecies; however, they reported a range of 17 to 19.4 mg GAE/g DW for TPC that was in agreement with our findings. Yan et al. (2016) 11 reported a higher amount of TPC (ranging from 93 to 135 mg GAE/g DW) in European oregano subspecies; however, the TPC in our subspecies is comparable to those reported by Pizzale et al. (2002) 25 and Capecka et al. (2005) 26 . These disparate results may imply that several factors, such as extraction methods, genetic controls, environmental parameters, growing and post-harvesting conditions as well as gene expression patterns may affect plant metabolism 27 .
Oregano contains different classes of polyphenolic compounds with multiple biological effects, such as antioxidant activity 28 . In the present study, a high antioxidant activity was observed among the oregano accessions and was confirmed by the abundance of their polyphenol compounds. The highest DPPH scavenging activities was detected in the Namin accession; however, the highest FRAP activity was obtained in Gilan accession. Interestingly, both accessions belonged to the virens subspecies. In accordance with our results, a high variation in antioxidant activity has previously been reported for oregano accessions 11 . Consistent with our results, the study of Lamien-Meda et al. (2010) 29 showed a similar range for DPPH antioxidant activity among the 19 sage (Salvia officinalis L.) accessions (ranging from 64 to 132 mg TE/g DW). However, the upper range in FRAP activity was observed in our study. The total and specific polyphenol compounds are the key accountable agents for the antioxidant activity of oregano species 30 . Previous studies on oregano showed a high antioxidant activity of dry leaves in olive oil as well as a high increase in oxidative stability of fried chips 9 .
Rosmarinic acid is one of the major phenolic acids which have been reported in the oregano plant 31 . This component is also found in other species of the Lamiaceae family, such as thyme (Thymus daenensis and Thymus vulgaris) 32 , basil (Ocimum spp.) 23 , mint (Mentha spp.) 33 and sage (Salvia officinalis) (Lamien-Meda et al., 2010). In the current study, a high variation in rosmarinic acid content was found among the oregano subspecies (ranged from 659.6 to 1646.9 mg/100 g DW). An approximately similar variation (ranging from 7.2 to 41.3 mg/g DW) was previously reported in oregano accessions by Yan et al. (2016) 11 . Interestingly, based on a previous study 34 , the average content of rosmarinic acid for subspecies vulgare was 18.2 mg/g DW which was in good agreement with our results. This variation in rosmarinic acid content provides valuable data for oregano breeders for achieving high quality plants with high antioxidant capacity for food and medicinal applications. Results also demonstrate that the main bioactive component of accession with the lowermost to the uppermost concentration varied among the oregano accessions, for instance about 6.6-fold for luteolin, tenfold for quercetin, threefold for coumarin and 5.8-fold for chicoric acid (Table 2). Our findings agree with those of others, specifying that oregano species were reported to be rich sources of phenolic acids such as rosmarinic acid 11 . In accordance with our results on phenolics, literature review illustrates the existence of diverse phenolic constituents in the Oregano species, such as rosmarinic acid, hydroxycinnamic acid, apigenin, luteolin, quercetin, scutellarein, apigenin-7-O-glucoside, luteolin-7-O-glucoside, and luteolin-7-O-glucuronide 11,18,35 . The composition of phenolic compounds in oregano depends on the accession, subspecies, genotype, and geographical and environmental conditions. It has been shown that the flavonoid and phenolic compounds of oregano can be used to differentiate between oregano chemotypes within the same species. For instance, the concentration of rosmarinic acid has been reported to vary between the chemotype within the species of O. vulgare subp. hirtum, O. vulgare subsp. vulgare and O. syriacum 34 .
Correlation analysis showed some relations among phenolic compounds and antioxidant attributes. However, there was not a significant correlation among the TFC and antioxidant activity. The results demonstrated a moderate correlation among DPPH antioxidant activity and rosmarinic acid, but no significant correlation was observed among FRAP activity and rosmarinic acid content. Similar to our results, Yan et al. (2016) 11 reported that there is no correlation between the antioxidant capacity measured by oxygen radical absorbance capacity and the concentration of rosmarinic acid. It can be specified that other phenolics and constituents can also be considered as antioxidant agents. In another study, the authors ascribed the antioxidant capacity of the similar fairly extract to the amount of rosmarinic acid and other bioactive polyphenol constituents 36 .

Conclusion
Oregano is a widely used medicinal herb due to its special phenolic compounds, aroma transfer and medicinal properties. In this study, the genetic variation in phenolic components and antioxidant activity of Iranian oregano accessions was investigated. A high variation was observed between the highest and lowest quantity of concentration of rosmarinic acid and other phenolics. These results might be utilized for breeding oregano plants containing high amount of phenolic components and exhibiting high antioxidant activity. Upgraded species and accessions with great polyphenolic compounds and small progressive duration is of high interest to plant producers. The extracts of herbs as well as other plant materials are gradually gaining attention in numerous industrial applications based on natural products. As a result, identification of different accessions which are richest in special components is very important for food and pharmaceutical industries. Finally, the survey of our study offers a broad-ranging variability among the Iranian oregano populations for utilizing elite accessions www.nature.com/scientificreports/ and subspecies (especially Ardabil and Namin accessions belonging to virens subspecies) based on phenolic components and antioxidant attributes for breeding programs and domestication purposes.