Variation in Flavonoid Compounds, Volatiles and Yield Related Traits in Different Iranian Rosa damascena Mill. Cultivars Based on SPME Arrow and LC-MS/MS

Damask rose (Rosa damascena Mill.) is an aromatic industrial plant with different applications. Selection of cultivars with high-value metabolites such as flavonoids—with acceptable yields—can lead to elite cultivars for mass propagation in various industries. A field experiment was carried out in a randomized complete block design (RCBD) to evaluate metabolites and some yield-related morphological data. In the present investigation, for the first time 13 flavonoid components of nine Iranian damask rose cultivars were compared using LC-MS/MS. As a result, 13 flavonoids were identified, most of which were reported for the first time in rose petals. Phloridzin (72.59–375.92 mg/100 g dw), diosmetin (82.48–153.16 mg/100 g dw) and biochanin A (0–1066.89 mg/100 g dw) were the most abundant, followed by trans-chalcone (0–106.29 mg/100 g dw) and diosmin (41.55–84.57 mg/100 g dw). Levels of naringenin also ranged from 3.77 in B111 to 54.70 mg/100 g dw in C294, while luteolin varied from 4.37 in B111 to 28.87 mg/100 g dw in C294. The SPME Arrow technique also was applied to determine the real aroma of the studied cultivars. Phenethyl alcohol was the most abundant compound, in the range of 69.28 to 77.58%. The highest citronellol/geraniol (C/G) was observed in D234 (4.52%) and D237 (4.30%), while the lowest amount belonged to A104 (1.28%). Rose oxide, as the most crucial factor for odor, ranged from 0.06% in D237 to 0.15% in D211. Based on cluster and principal component analysis (PCA), D234 cultivar can be suggested as a promising cultivar with high yield, high C/G content and high rose oxide, while D234 and C294 were the most valuable cultivars in terms of flavonoids with high yield. Finally, these cultivars can be introduced for further breeding programs and industrial cultivation.


Introduction
Damask rose (Rosa damascena Mill.) is an aromatic plant belonging to the Rosaceae family, with approximately 200 species and 18,000 cultivars.Roses are mostly scattered in the northern hemisphere [1,2].R. damascena has been considered as the most promising species for industrial purposes [3].This species has been cultivated in many parts of the world, especially in Turkey and Iran [2,4].Different rose products such as rose water, oil, concrete, and absolute have been produced from R. damascena oil.Based on the preliminary production, distilled oil from damask rose was obtained in Persia [5,6].A Persian physician (Avicenna) distilled rose flowers for medical applications.Different activities were determined for rose petals, including anti-anxiety [7], antibacterial and anti-Trichomonas activity [8] and anti-inflammatory and antioxidant activity [9].The presence of volatiles and non-volatiles in R. damascena provides insights for further genetic improvement of the cultivars.R. damascena has been used for both ornamental and industrial applications [10] and this kind of diversity leads to different rose chemotypes in different countries.The essential oil yield and components of rose oil are attributed to different environmental and genetic factors.Most of the previous research has been reported on the essential oil variation of R. damascena in different countries including India [11], Bulgaria [12,13], Iran [14], and Pakistan [15].In most of these studies, the oils were analyzed by GC/MS instruments.The citronellol/geraniol ratio (C/G), is one of the main factors determining the quality of oil in rose petals [2].However, since the aroma of R. damascena is one of the most important factors for its quality, solid-phase microextraction (SPME)-based techniques can provide more insightful results [16].There is limited research regarding the SPME of rose petals, including in Egyptian [17], Turkish [18], and Chinese roses [19].Moreover, another neglected issue in the previous researches is the analysis of non-volatiles, especially phenolic acids and flavonoids in rose petals using the LC/MS technique.However, for this important aspect, there is also limited research literature on rose petals harvested in Turkey [20] and Morocco [21] using the HPLC method.Furthermore, most of the previous reports on Iranian rose damask cultivars focus on essential oil variation based on GC/MS analysis [2,14,22].Furthermore, there is limited research regarding the determination of one or two components, e.g., kaempferol and quercetin based on HPLC analysis of rose petals [23], and there is no report on the comprehensive polyphenolic compounds of Iranian rose petals based on LC/MS methodology in different cultivars.So, this report for the first time provides new information on new polyphenolic compounds of Iranian rose petals, as well as real aroma rose petals, using the highly validated method SPME Arrow.Finally, the combination of volatiles, non-volatiles, and morphological data provides new insights for further selection and introduction of new cultivars for further industrial or breeding purposes in rose petals.
The objectives of this study were (1) to analyze nine Iranian rose damask cultivars in terms of volatiles using SPME Arrow coupled to GC/MS, (2) to determine the phenolic and flavonoid compounds of the cultivars studied using the LC-MS/MS technique, (3) to compare the traits related to morphology and yield in studied plants and (4) to use multivariate analyses to find the correlation among metabolite and morphological traits.

Plant Materials
The cultivars are the result of a long-term project aiming to collect, compare, and select cultivars and accessions with superior traits, such as yield, resistance to environmental stresses, blooming period, essential oil content and composition, and flower characteristics from all over Iran.Table 1 presents the characteristics of the cultivars.The botanical identification was performed by Prof. Mozaffarian using Flora Iranica [24] and finally, the cultivars are selected and provided by the Research Institute of Forests and Rangelands of Iran.The studied cultivars are illustrated in Figure 1.

Field Experimentand Sampling
Nine cuttings of R. damascena cultivars were cultivated in a randomized block design (RCBD) on 20 March 2019 with three replicates and the traits were evaluated in three-yearold rose damask cultivars.Information on the cultivars studied is illustrated in Table 1.The study was carried out at the research farm of the Isfahan University of Technology (32 • 38 ′ N, 51 • 39 ′ E) and an altitude of 1620 m above sea level.Each experimental unit consisted of two rows of 2 m in length.The soil texture was clay loam with pH = 7.The LC/MS analysis was carried out according to the protocol explained by Rahimmalek et al., 2023 [25].The drying of the samples was performed before methanolic extraction.For this purpose, the rose petals were subjected to shade drying in 42% humidity, 25 • C for 10 days.After the drying process, 1 g of dried powdered rose petals was placed in a plastic falcon (50 mL of volume).Then, 15 mL of pure analysis-grade methanol (Sigma-Aldrich, Steinheim, Germany) was added.The mixtures prepared for extraction were placed for 1 h on the rotary shaker (120 rpm), then the samples were centrifuged (13,000 rpm, 10 min) and the supernatant was kept in a separate container.Finally, the samples were extracted two times with the same procedure.Then, 45 mL of the extract was transferred to a round-bottomed flask (100 mL) and partially evaporated by a rotary evaporator (Heidolph, Germany), then transferred to smaller round-bottomed flask (10 mL) and finally the solvent was completely evaporated.The extract was then resolved in 5 mL of pure hypergrade for LC-MS LiChrosolv ® methanol (Sigma-Aldrich, Steinheim, Germany).Finally, the extracts were centrifuged at 13,000 rpm for 10 min and diluted 10 times with methanol.All extracts were prepared in three replicates.
Identification was performed based on the optimised MRM mode (Table 2) according to polyphenol standards as a mixture (MetaSci, Toronto, ON, Canada) while quantification was carried out based on the external standard method.

HS-SPME-GC-MS 2.5.1. Sample Preparation
The rose petals were dried at 25 • C for 6 days in the shade.Subsequently, 50 mg of powdered samples along with 2-undecanone as the internal standard were used for analysis.

HS-SPME-GC-MS
The extraction was carried out with 1.10 mm DVB/C-WR/PDMS SPME Arrow fibre (Shimadzu, Kyoto, Japan).Before extraction, the samples were conditioned at 45 • C for 5 min and then the volatile components were extracted for 30 min at the same temperature.The desorption of the analytes was performed in the apparatus injector.
For this purpose, Shimadzu GCMS QP 2020 Plus (Shimadzu, Kyoto, Japan) equipped with a Zebron ZB-5 MSi capillary column ( The identification of analytes was based on the experimentally obtained mass spectra and linear retention indexes.As a reference the library of Flavors and Fragrances of Natural and Synthetic Compounds 3.0 (FFNSC 3.0) was applied.The analyte quantity calculation was based on peak area normalization.The analyses were carried out in three replicates.

Statistical Analysis
Analysis of variance (ANOVA) was performed using the SAS statistical programme ver.9.4, (SAS Inc., Cary, NC, USA) and GLM method.Fisher's least significant difference was applied for mean comparisons at p ≤ 0.05.Cluster and principal component analysis (PCA) analyses were designed by STATGRAPHICS ver.18.2.04software and Pearson correlation coefficient calculation were performed using SAS ver.9.4.

Chemical Composition of Rose Petals Based on SPME Arrow
The SPME Arrow results revealed high variation among the studied rose cultivars.The main components of the rose petal oil were phenethyl alcohol, citronellol, and geraniol (Table 3).Phenethyl alcohol ranged from 69.28% in B215 to 77.58% in A104, while in previous research, it ranged from 44.5% to 47.2% according to the distillation method [2].Phenyl ethyl alcohol is known for a fragrant rose-like odor with lilac honey note [17].Citronellol was also in the range of 4.48% in A104 to 9.61% in A105.Geraniol also ranged from 1.36% in D234 to 3.49% in A104.The highest amount of citronellol/geraniol (C/G) was observed in D234 (4.52%) and D237 (4.30%) followed by A105 (4.26%), while the lowest amount belonged to A104 (1.28%).Yaghoobi et al. [2] also reported higher ranges for C/G (3.3-6.9%)compared to our studied cultivars based on distilled oil analysis.These researchers highlighted that damasks that grow in semi-arid, cold condition possessed a higher C/G ratio, while those that grow in the same condition with low altitude had a lower C/G ratio [2].Similarly, in the present investigation, the highest C/G ratio based on the SPME method was obtained in D234 from Bardsir in high altitude, cold temperature, and semi-arid condition (Table 3).Rose oxide is one of the most important factors that determines the odor of rose petals.In the present investigation, it ranged from 0.062% in D237 from Maymend, Fars with the lowest altitude to 0.15% in D211.It might be suggested that the rose oxide odor increases in high altitudes and arid conditions.This amount is much higher than those reported by Bulagrian rose petals [12], but less than Indian [10] and Chinese [19].Mohsen et al. [17] also reported a similar range for rose oxide in Egyptian rose petals.These researchers also reported a relatively similar range for geraniol in Egyptian rose petals, suggesting that SPME provides the real aroma of the rose petals.However, there are no reports regarding the SMPE of Iranian rose petals for comparison.Aroma is one of the most important factors determining the quality of rose petals, especially for the production of perfumes.Among the compounds, rose oxide is considered the unique scent of damask rose that exerts a sweet floral aroma [3,17].However, citronellol/geraniol (C/G) is also crucial for the aroma and quality of rose petals.These compounds can be highly affected by harvesting time, plant growth stage, temperature, and drying methods [22].In the present investigation the samples were harvested at full bloom stage, in the morning at 25 • C .In most cases, harvesting in the early stage of petals leads to a decrease in citronellol and geraniol content [17].Means with different letters are statistically significant by LSD at 5% level probability.

Flavonoid Compounds of Rose Petals Based on LC/MS
In the present investigation, the polyphenolic compounds of Iranian rose petals were determined using LC/MS analysis.Table 4 shows the qualitative and quantitative analysis of the rose petals results.The MRM details of the flavonoids are illustrated in Table 2. Consequently, 13 flavonoids were identified, most of which were reported for the first time in rose petals.In terms of the total flavonoid content, the highest was found for the C294 sample followed by the D234 sample, while the lowest content was observed for the A104 sample.The total amounts of flavonoids found in the present study were significantly lower than those found by Alizadeh and Fattahi [10] in Iranian rose petals.Such differences can occur due to different aspects, such as weather conditions or soil fertility; moreover, there was significant difference in sample preparation, since Alizadeh and Fattahi [10] extracted flavonoids from fresh rose petals with 80% methanol, assisted by ultrasound, while we have used dried rose petals extracted with pure methanol without ultrasound assistance.Biochanin A (0-1066.89 mg/100 g dw), phloridzin (72.59-375.92mg/100 g dw) and diosmetin (82.48-153.16mg/100 g dw) were the most abundant in nine studied cultivars followed by trans-chalcone (0-106.29 mg/100 g dw) and diosmin (41.55-84.57mg/ 100 g dw).More commonly identified in roses and other Rosaceae plants [26,27], in previous studies naringenin ranged from 3.77 in B111 to 54.70 mg/100 g dw in C294, luteolin ranged from 4.37 in B111 to 28.87 mg/100 g dw in C294, while kaempferol was in the range of 8.56 in B111 to 67.81 mg/100 g dw in C294.
Previous reports have mostly determined the flavonoids of rose petals [17,28].Kumar et al. [29] determined kaempferol, quercetin, as important flavonoids in Indian rose petals, which is consistent with our results.However, the present study introduced new flavonoids with their amount in nine studied Iranian cultivars, providing new insights for producing of health benefit products from elite cultivars.

Analysis of Variance and Multivariate Analyses
Based on analysis, some cultivars and some traits were significantly different (Tables 5 and 6).In terms of the studied cultivars, no significant differences were observed (Table 7), except for the length of the bud, the diameter of the bud and dry weight of semi-open buds (Table 7).According to ANOVA analysis, some of the traits studied were non-significant (Table 6).Therefore, multivariate analyses were applied for a better interpretation of the data.The principal component analysis (PCA) of the morphological and essential oil components revealed that their main components interpreted 66.67% of total variation.Based on the results, PC1 revealed a positive correlation of benzyl alcohol, rose oxide, citronellol and geraniol with the dry weight of opened and semi-opened buds and number of petals and negative correlation with other morphological traits.PC2 showed 17.91% of total variation and revealed the positive correlation of geraniol, benzyl alcohol with bud diameter, petal length, dry weight of semi-opened bud, and negative correlation with others, such as flower fresh weight.PC3 also revealed a positive correlation of benzyl alcohol, rose oxide, and citronellol with some morphological traits such as bud diameter, dry and fresh weight (Figure 2).
Previous reports have mostly determined the flavonoids of rose petals [17,28].Kumar et al. [29] determined kaempferol, quercetin, as important flavonoids in Indian rose petals, which is consistent with our results.However, the present study introduced new flavonoids with their amount in nine studied Iranian cultivars, providing new insights for producing of health benefit products from elite cultivars.

Analysis of Variance and Multivariate Analyses
Based on analysis, some cultivars and some traits were significantly different (Tables 5 and 6).In terms of the studied cultivars, no significant differences were observed (Table 7), except for the length of the bud, the diameter of the bud and dry weight of semi-open buds (Table 7).According to ANOVA analysis, some of the traits studied were non-significant (Table 6).Therefore, multivariate analyses were applied for a better interpretation of the data.The principal component analysis (PCA) of the morphological and essential oil components revealed that their main components interpreted 66.67% of total variation.Based on the results, PC1 revealed a positive correlation of benzyl alcohol, rose oxide, citronellol and geraniol with the dry weight of opened and semi-opened buds and number of petals and negative correlation with other morphological traits.PC2 showed 17.91% of total variation and revealed the positive correlation of geraniol, benzyl alcohol with bud diameter, petal length, dry weight of semi-opened bud, and negative correlation with others, such as flower fresh weight.PC3 also revealed a positive correlation of benzyl alcohol, rose oxide, and citronellol with some morphological traits such as bud diameter, dry and fresh weight (Figure 2).Cluster analysis based on essential oil components also confirmed the results of PCA analysis in most cases.However, based on this analysis, D234, D237 and D231 were categorized in separate groups that suggests their different phytochemical patterns (Figure 3).Group 1 included A104 and A105, while group 2 consisted of B215, B111, C193, and C231 with moderate C/G content.
Group 1 included A104 and A105, while group 2 consisted of B215, B111, C193, and C231 with moderate C/G content.Based on yield-related traits obtained in the previous assays in the studied cultivars (Table 1) and new metabolite analyses, the D234 cultivar can be suggested as a promising cultivar with high yield, high C/G content and high rose oxide and can be introduced for further industrial cultivation.The second cultivar, A105, has a relatively high yield and high aroma quality with a high C/G ratio and rose oxide.
Multivariate analyses were also performed for flavonoids in rose petals.Based on cluster analysis, two major groups were obtained in which D234 and C294 were classified in group 2 with high amounts of biochanin A and diosmetin, and the other cultivars were classified in group 1 (Figure 3).Group 1, divided into two subgroups, included B111, A105, A104, and B215 with low amount of chalcone.Subgroup 2 consisted of C193, C251, and D237 with a high amount of chalcone and a moderate amount of diosmetin.Interestingly, D234 and C294 were the most valuable cultivars in terms of flavonoids.For this point of view, both cultivars can be introduced for further breeding programs and industrial cultivation as they also have high yields based on previous results (Table 1).

Correlation Analysis of Metabolites and Environmental Factors
Rose oxide and citronellol/geraniol showed a high positive correlation (r = 0.76) (Table 8).Biochanin A and chalcone in the flavonoid class also revealed a high positive correlation (r = 0.84).Among flavonoids, diosmin showed negative correlation with average temperature (r = −0.78).It suggests that diosmin increases mainly at lower temperature.Previous reports also highlighted that high temperature can decrease the amount of some flavonoids by degrading their structure [30].The stability and biological properties of flavonoids can be highly affected by temperature [31].Most of the flavonoids are sensitive to high temperatures; however, glycosylated flavonoids are more resistant to high temperatures than aglycon ones.The degradation of flavonoids depends on their structural solidity, and degradation can lead to a decrease in their antioxidant capacity.So, for selecting high-diosmin rose petals, the low temperature condition might be beneficial [31].Elevation had a positive correlation with citronellol (r = 0.49), while a negative correlation with geraniol (r = −0.74).Ghavam et al. [32] obtained similar results in damask rose for citronellol amount as they suggested that higher altitudes can lead to increased biosynthesis of some components such as citronellol, due to higher exposure to UV radiation and Based on yield-related traits obtained in the previous assays in the studied cultivars (Table 1) and new metabolite analyses, the D234 cultivar can be suggested as a promising cultivar with high yield, high C/G content and high rose oxide and can be introduced for further industrial cultivation.The second cultivar, A105, has a relatively high yield and high aroma quality with a high C/G ratio and rose oxide.
Multivariate analyses were also performed for flavonoids in rose petals.Based on cluster analysis, two major groups were obtained in which D234 and C294 were classified in group 2 with high amounts of biochanin A and diosmetin, and the other cultivars were classified in group 1 (Figure 3).Group 1, divided into two subgroups, included B111, A105, A104, and B215 with low amount of chalcone.Subgroup 2 consisted of C193, C251, and D237 with a high amount of chalcone and a moderate amount of diosmetin.Interestingly, D234 and C294 were the most valuable cultivars in terms of flavonoids.For this point of view, both cultivars can be introduced for further breeding programs and industrial cultivation as they also have high yields based on previous results (Table 1).

Correlation Analysis of Metabolites and Environmental Factors
Rose oxide and citronellol/geraniol showed a high positive correlation (r = 0.76) (Table 8).Biochanin A and chalcone in the flavonoid class also revealed a high positive correlation (r = 0.84).Among flavonoids, diosmin showed negative correlation with average temperature (r = −0.78).It suggests that diosmin increases mainly at lower temperature.Previous reports also highlighted that high temperature can decrease the amount of some flavonoids by degrading their structure [30].The stability and biological properties of flavonoids can be highly affected by temperature [31].Most of the flavonoids are sensitive to high temperatures; however, glycosylated flavonoids are more resistant to high temperatures than aglycon ones.The degradation of flavonoids depends on their structural solidity, and degradation can lead to a decrease in their antioxidant capacity.So, for selecting highdiosmin rose petals, the low temperature condition might be beneficial [31].Elevation had a positive correlation with citronellol (r = 0.49), while a negative correlation with geraniol (r = −0.74).Ghavam et al. [32] obtained similar results in damask rose for citronellol amount as they suggested that higher altitudes can lead to increased biosynthesis of some components such as citronellol, due to higher exposure to UV radiation and reduced temperature [22].However, the average annual rainfall revealed a different trend compared to elevation in terms of these compounds.However, a previous report on Citrus hystrix revealed that citronellol has a negative correlation with rainfall intensity [33].Although a high variation was reported within individual flavonoid groups in response to elevational changes, it can be suggested that long days with cool night temperatures are the most crucial one for the biosynthesis of flavonoids [34].

Conclusions
Comprehensive research was performed to evaluate the variation of flavonoids, aroma and yield related traits in Iranian damask rose cultivars to introduce the potent ones for further breeding and industrial programs.This research provided new information on most flavonoids in Iranian rose petals.Some new non-volatiles were introduced with strong health benefit compounds.Furthermore, the amounts of citronellol/geraniol and rose oxide were screened in the studied cultivars and some cultivars were introduced based on multivariate analyses and correlation between metabolites and yield-related traits.Consequently, D234 along with C294 were introduced as the most valuable cultivars with high yield, aroma, and major flavonoid components.Finally, the present data provided new insights for further breeding programs in introduced cultivars to improve the volatiles or non-volatiles in damask rose.

Foods
5. Bud stage to full flowering takes approximately 8-11 days.For metabolite analyses, the petals were harvested from semi-opened flowers on 22-25 April 2022 in 18.3-20 • C and 34-35.6%relative humidity at 5 am in the morning.

Figure 2 .
Figure 2. PCA analysis based on flavonoid compounds in nine damask rose cultivars.Figure 2. PCA analysis based on flavonoid compounds in nine damask rose cultivars.

Figure 2 .
Figure 2. PCA analysis based on flavonoid compounds in nine damask rose cultivars.Figure 2. PCA analysis based on flavonoid compounds in nine damask rose cultivars.

Figure 3 .
Figure 3. Dendrogram generated from cluster analysis of nine rose damask cultivars based on flavonoid compounds using WARD based on the squared Euclidean dissimilarity calculated.

Figure 3 .
Figure 3. Dendrogram generated from cluster analysis of nine rose damask cultivars based on flavonoid compounds using WARD based on the squared Euclidean dissimilarity calculated.

Table 1 .
The characteristics of studied damask rose cultivars.

Table 2 .
Flavonoid MRM details in studied rose damask cultivars.

Table 3 .
Chemical composition of the studied cultivars based on SPME Arrow.

Table 4 .
The mean of flavonoid compounds based on LC/MS.In each row means followed by the same letter are not significantly different at p = 0.05.
Note: means followed by same letter are not significantly different according to LSD's test at 5% level.

Table 5 .
Analysis of variance of block and cultivar effect on morphological characteristics of nine rose cultivars.
*, ** and ns are respectively significant at the 5%, 1% probability level and the absence of significant difference.

Table 7 .
Comparison of the average effect of genotype on the morphological characteristics of 9 rose genotypes.

Table 8 .
Pearson correlation coefficients between major flavonoids, essential oils compounds and environmental conditions.