Chemometric Investigation and Antimicrobial Activity of Salvia rosmarinus Spenn Essential Oils

To ensure the better production and sustainable management of natural resources, a chemometric investigation was conducted to examine the effect of cooperative and harvesting periods on the crop yields and chemical compositions of Salvia rosmarinus Spenn essential oils in the Oriental region of Morocco. The samples were collected from three cooperatives over nine time periods from January 2018 to April 2019. The chemical composition of Salvia rosmarinus Spenn essential oils was analyzed by gas chromatography coupled with mass spectrometry. The data from this study were processed by multivariate analyses, including principal component analysis (PCA) and hierarchical cluster analysis (HCA). The disc diffusion technique and a determination of the minimal inhibitory concentration were performed to study the antibacterial properties of the oils. Statistical analysis showed that the cooperative and harvest period have a significant effect on yields. The highest yield of essential oil was recorded in April 2019 at cooperative C1. The PCA and the HCA results were divided into two groups: Group A for the summer season and group B for the winter season. The samples collected during summer were characterized by a high amount of 1,8-cineole component and a high yield of essential oil, whereas the samples collected during winter were qualified by a high amount of α-pinene component and a low yield of essential oil. The antibacterial activity of Salvia rosmarinus Spenn essential oils showed that Mycobacterium smegmatis ATCC23857 and Bacillus subtilis ATCC 23857 are the most susceptible strains, stopping growth at 1/500 (v/v). The least susceptible strain is Escherichia coli ATCC25922, with an MIC value corresponding to 1/250 (v/v). The findings of this study could have a positive economic impact on the exploitation of rosemary in the Oriental region, especially during the best harvest periods, as they indicate how to obtain the best yields of oils richest in 1,8-cineole and α-pinene chemotypes.

the hot weather because the biosynthesis of EOs is more prominent in warm climates. This was shown by Rehman et al. [28], who found that the majority of enzymatic activities for the production of volatile compounds are enhanced by temperature.
Accordingly, the yield of S. rosmarinus essential oil was maximized by using plants from cooperative C1 that were harvested in April. The variability of rosemary's essential oil yields could have resulted from the origin of the plant, the harvest time [29,30], the environmental [27] and agronomic conditions [31,32], the plant's phenological stage (Before and during full flowering stages) [33,34], and the extraction method [35].  Figure 1 demonstrates the significant differences between the cooperatives and harvest periods in terms of the average yield of S. rosmarinus essential oil. There was a significant difference (p < 0.05) in the average yield of S. rosmarinus EO between cooperatives C1 and C2 as well as between cooperatives C1 and C3, while no significant differences (p > 0.05) were noticed among cooperatives C2 and C3 in terms of the average yield of S. rosmarinus EO. The highest average yields of the rosemary EOs was found in plants from cooperative C1. Regarding harvest periods, no significant differences were observed among periods P1 (January), P8 (February), P7 (December), and P6 (October) in terms of the average yield of S. rosmarinus EO. In addition, no significant differences were noted among periods P2 (March), P5 (September), P3 (May), and P9 (April), whereas significant differences were noted between the periods P2 (March), P4 (July), and P3 (May) in terms of the average yield of S. rosmarinus EO. However, highly significant differences were identified among P1 (January) and P9 (April), as well as between P8 (February) and P3 (May). The samples from Period 9 (collected in April) and the samples from Period 3 (collected in May) showed the highest average yields of S. rosmarinus EOs (3.1 ± 0.1 and 3.0 ± 0.1, respectively). The lowest average yield of S. rosmarinus EO (1.1 ± 0.1) was found in the plant collected in period 1 (January). According to two-way ANOVA, the yield of EOs was significantly affected by the harvest period and the cooperative. This finding may be explained by the environmental conditions in Figuig (a Saharan arid climate), especially the hot weather because the biosynthesis of EOs is more prominent in warm climates. This was shown by Rehman et al. [28], who found that the majority of enzymatic activities for the production of volatile compounds are enhanced by temperature.
Accordingly, the yield of S. rosmarinus essential oil was maximized by using plants from cooperative C1 that were harvested in April. The variability of rosemary's essential oil yields could have resulted from the origin of the plant, the harvest time [29,30], the environmental [27] and agronomic conditions [31,32], the plant's phenological stage (Before and during full flowering stages) [33,34], and the extraction method [35]. Table 1 shows the results of the major components of the chemical composition of S. rosmarinus essential oils identified by gas chromatographic analysis coupled with mass spectrometry (GC/MS) for the three cooperatives during January 2018 and April 2019. Rosemary EOs contained different components according to the harvest period. Eight components were identified, representing a total of 83.7% for the three cooperatives. The major constituents were 1,8-cineole (28.6-51.1%), α-pinene (9.9-16.2%), camphor (5.3-16.8%), β-pinene (2.2-8.0%), camphene (2.3-7.7%), myrcene (0-4.5%), α-terpineol (0-3.8%), and limonene (0-3.3%) for all the samples during all periods. The essential oil from cooperative C1 showed the highest percentage of components. The 1,8-cineole component was present in a higher proportion than the other components, with a maximum value during period 3 (May, 2018); the α-pinene component reached its maximum value during period 6 (October 2018).

Principal Component Analysis
Principal component analysis (PCA) allowed us to identify the correlation among the major compounds and the yield of rosemary essential oil from different cooperatives and harvest periods, as well as to classify similar samples in accordance with the yield and the major compounds of rosemary essential oil. PCA also enabled the explicit representation of individuals and variables to highlight the different profiles between samples according to the harvesting period and the cooperative. The individuals are represented by the 27 essential oil samples taken from three cooperatives during nine different harvest periods. The variables are represented by the average yield in essential oil and the main components of the chemical composition of rosemary. To verify the feasibility of the PCA, the Kaiser-Meyer-Olkin (KMO) index and the Bartlett test were used. Calculation of the KMO index resulted in a value of 0.8, which is greater than 0.5; and Bartlett test has obtained a probability of less than 5%, meaning the test is highly significant. Therefore, the PCA is attainable.
To determine the number of principal components, the components with eigenvalues were greater than one had to be retained [48]. Table 2 demonstrates that the first three principal components had eigenvalues greater than one and that they explained most of the variables. The cumulative percentage of these three components explained 82.9% of the data variability. The first component explained 47.5% of the variables: α-pinene, 1,8-cineole, camphor, α-terpineol, and yield. The second component explained 22.5% of the variables: camphor, β-pinene, and limonene. The third component explained 12.9% of the variables: myrcene. Table 2. Quality of representation of the variables on the first three components, eigenvalues of the components, percentages of variability explained, and cumulative percentages. The loading plot or distribution of the variables according to the first principal component (PC1) and the second principal component (PC2) (Figure 2) shows that the variables α-pinene, 1,8-cineole, camphor, yield, and α-terpineol are well represented in axis 1 (PC1), and the variables camphor, β-pinene, and limonene are well represented in axis 2 (PC2). Thus, the loading plot also shows the presence of some correlation among the studied variables. A positive correlation was identified among the yield of essential oil, α-terpineol, and 1,8-cineole. α-pinene and camphor also showed a positive correlation. However, 1,8-cineole and camphor were negatively correlated. Another negative correlation was observed between α-pinene and the yield of essential oil. These results indicate that when the yield of essential oil increased, the amount of 1,8-cineole increased. However, when the quantity of 1,8-cineole contained in essential oil was higher, the quantity of α-pinene and camphor was lower. that when the yield of essential oil increased, the amount of 1,8-cineole increased. However, when the quantity of 1,8-cineole contained in essential oil was higher, the quantity of α-pinene and camphor was lower. The score plot (Figure 3) or the distribution of individuals (harvest period and cooperative) in accordance with PC1 and PC2 shows that individuals were separated into two main groups. Group A consists of P9 (April), P3 (May), P4 (July), and P5 (September), or the summer season. Group B is composed of P6 (October), P1 (January), P7 (December), and P8 (February), or the winter season. The intermediary between the two groups was P2 (March), which shared both of their characteristics. No significant difference was found between the essential oils of the three cooperatives because the distribution of cooperatives in the groups was random. This finding may also be due to the locations of the cooperatives being in the same province (Figuig). The score plot (Figure 3) or the distribution of individuals (harvest period and cooperative) in accordance with PC1 and PC2 shows that individuals were separated into two main groups. Group A consists of P9 (April), P3 (May), P4 (July), and P5 (September), or the summer season. Group B is composed of P6 (October), P1 (January), P7 (December), and P8 (February), or the winter season. The intermediary between the two groups was P2 (March), which shared both of their characteristics. No significant difference was found between the essential oils of the three cooperatives because the distribution of cooperatives in the groups was random. This finding may also be due to the locations of the cooperatives being in the same province (Figuig). To obtain the biplot, the individual factors were projected on the variables (Figure 4). The graph shows that group A is dominated by the highest amount of 1,8-cineole and αterpineol compounds and a high yield of essential oil. However, group B presented the opposite trend of group A, which is characterized by α-pinene and camphor compounds, whereas group B was characterized by a low yield of essential oil and a low amount of 1,8-cineole and α-terpineol compounds, especially in P9 (October) and P7 (December) individuals. Period 2 (March) was characterized by an average yield of essential oil and an To obtain the biplot, the individual factors were projected on the variables (Figure 4). The graph shows that group A is dominated by the highest amount of 1,8-cineole and α-terpineol compounds and a high yield of essential oil. However, group B presented the opposite trend of group A, which is characterized by α-pinene and camphor compounds, whereas group B was characterized by a low yield of essential oil and a low amount of 1,8-cineole and α-terpineol compounds, especially in P9 (October) and P7 (December) individuals. Period 2 (March) was characterized by an average yield of essential oil and an average amount of the constituents of both groups. Several studies have shown the effect of the harvest period and climatic conditions on the level of certain compounds in rosemary essential oil [34]. Melito et al. [49] found that summer had the highest yield of S. rosmarinus essential oils, whereas winter had the lowest yield. These authors confirmed that the chemical composition of S. rosmarinus was significantly influenced by the seasons, observing that γ-terpinene and terpinolene values were high in summer (0.8 and 0.7, respectively) and lowest in winter. They also reported that β-pinene and bornyl acetate components were the most elevated in the spring. Ismaili et al. [50] reported that the best harvest time was spring, with a yield of 2.0%, and they also showed that the chemical composition did not vary according to the season. Lemos et al. [51] noted that yields varied during the harvest period, with the highest yield of rosemary essential oil occurring in April (0.90%). The authors of [34,51] showed that rainfall and light intensity directly influence the production of essential oil.

Hierarchical Cluster Analysis
To better visualize the classification of the studied samples (harvest period) according to the main components of chemical composition and the yield of rosemary essential oil, an HCA was performed ( Figure 5). As a validation of PCA results (scores plot), the individuals were divided into two main groups: Group A represented summer (periods P2, P3, P4, P5, and P9) and Group B represented winter (periods P2, P6, P7, and P8). March (P2) was included in both groups. Thus, it represented the intermediary between the two groups. Group A (winter season) was characterized by a high amount of α-pinene and camphor, a low amount of 1,8-cineole and α-terpineol compounds, and a low yield of essential oil. However, group B (summer season) was characterized by a high amount of 1,8-cineole and α-terpineol compounds, a high yield of essential oil (red color), and a low amount of α-pinene and camphor (blue color). The intermediate group (

Hierarchical Cluster Analysis
To better visualize the classification of the studied samples (harvest period) according to the main components of chemical composition and the yield of rosemary essential oil, an HCA was performed ( Figure 5). As a validation of PCA results (scores plot), the individuals were divided into two main groups: Group A represented summer (periods P2, P3, P4, P5, and P9) and Group B represented winter (periods P2, P6, P7, and P8). March (P2) was included in both groups. Thus, it represented the intermediary between the two groups. Group A (winter season) was characterized by a high amount of αpinene and camphor, a low amount of 1,8-cineole and α-terpineol compounds, and a low yield of essential oil. However, group B (summer season) was characterized by a high amount of 1,8-cineole and α-terpineol compounds, a high yield of essential oil (red color), and a low amount of α-pinene and camphor (blue color). The intermediate group (P2: March) was characterized by an average yield of essential oil and an average amount of S. rosmarinus main compounds. The result of both PCA and HCA indicated that the yield and chemical composition of essential oil were significantly influenced by the harvest period and environmental conditions. The majority of studies indicate that the yield of essential oil and the chemical composition of S. rosmarinus are impacted by environmental conditions. Bajalan et al. [52] found a correlation between the major compounds of rosemary essential oil and environmental factors. Sarmoum et al. [53] indicated that the production of rosemary essential oil was directly affected by water stress; this explains why the yield of rosemary essential oil was higher in the summer season and lower in the winter season. Nevertheless, Sadeh et al. [36] showed that the main compounds of rosemary essential oil were not significantly affected by season.

Antimicrobial Activity
The essential oils tested were composed of the average of the samples collected during the third period (May 2018), which belongs to summer season (high yield of S. rosmarinus essential oil and the highest proportion of 1,8-cienole compound). The antimicrobial activity was first qualitatively tested by the disc diffusion technique with the aim of selecting the most active EOs among those tested. These samples were evaluated against one Gram-negative (Escherichia coli) and two Gram-positive (Bacillus subtilis) (Mycobacte-

Antimicrobial Activity
The essential oils tested were composed of the average of the samples collected during the third period (May 2018), which belongs to summer season (high yield of S. rosmarinus essential oil and the highest proportion of 1,8-cienole compound). The antimicrobial activity was first qualitatively tested by the disc diffusion technique with the aim of selecting the most active EOs among those tested. These samples were evaluated against one Gram-negative (Escherichia coli) and two Gram-positive (Bacillus subtilis) (Mycobacterium smegmatis) species. Figure 6 illustrates the significant differences (p < 0.05) in the diameter of inhibition zones of the rosemary essential oil from three cooperatives: C1, C2, and C3. Significant differences were observed among the essential oil from all cooperatives and antibiotics in terms of all bacteria strain inhibition zones. However, no significant difference (p > 0.05) was noted between the essential oil from cooperatives C1 and C3, whereas a significant difference was found among the essential oils from the cooperatives C1 and C2, as well as antibiotics in terms of E. coli inhibition zones. In terms of inhibition zones of B. subtilis, a significant difference was noticed between essential oil from the cooperatives C1, C2, and C3, whereas no significant difference was observed between essential oils from the cooperatives C2 and antibiotics. Regarding M. smegmatis, again, significant differences were identified among the essential oils from the cooperatives C1 and C2, as well as antibiotics, whereas no significant difference was identified between the essential oils from cooperatives C2 and C3. Two-way ANOVA showed that there was a statistically significant difference among the means of the inhibition diameters and three cooperative EOs for the three strains used. C3, whereas no significant difference was observed between essential oils from the cooperatives C2 and antibiotics. Regarding M. smegmatis, again, significant differences were identified among the essential oils from the cooperatives C1 and C2, as well as antibiotics, whereas no significant difference was identified between the essential oils from cooperatives C2 and C3. Two-way ANOVA showed that there was a statistically significant difference among the means of the inhibition diameters and three cooperative EOs for the three strains used.  letters (a, b, c, d, e, f, and g) represent a significant difference (p < 0.05) based on Tukey test. Cn EO: essential oil from different cooperatives. Figure 6 shows that the essential oil from cooperative C1 had the highest activity, with inhibition zones of 16.6 ± 0.5 mm, 31.3 ± 2.1 mm, and 37.1 ± 1.5 mm for E. coli ATCC25922, B. subtilis ATCC 23857, and M. smegmatis MC2 155, respectively. In the essential oil from cooperative C2, the greatest activity was noted against M. smegmatis, with an inhibition diameter of 26.2 ± 1.1 mm. The lowest activity was observed against E. coli, with inhibition diameters of 19.0 ± 1.0 mm. The same results were observed in the essential oil from cooperative C3. The inhibition diameters obtained for M. smegmatis (27.0 ± 1.1 mm) were higher than those obtained for B. subtilis (18.6 ± 2.1 mm) and E. coli (15.2 ± 1.1 mm). Many studies have demonstrated that rosemary essential oil provides effective antibacterial activity against various microorganisms [38,40]. Rosemary essential oil was slightly less sensitive against Gram-negative bacteria (E. coli) than they against Gram-positive bacteria [54]. Megzari et al. [55] observed that rosemary essential oil presented antibacterial activity against the same bacterial strain, with diameters of inhibition similar to  letters (a, b, c, d, e, f, and g) represent a significant difference (p < 0.05) based on Tukey test. Cn EO: essential oil from different cooperatives. Figure 6 shows that the essential oil from cooperative C1 had the highest activity, with inhibition zones of 16.6 ± 0.5 mm, 31.3 ± 2.1 mm, and 37.1 ± 1.5 mm for E. coli ATCC25922, B. subtilis ATCC 23857, and M. smegmatis MC2 155, respectively. In the essential oil from cooperative C2, the greatest activity was noted against M. smegmatis, with an inhibition diameter of 26.2 ± 1.1 mm. The lowest activity was observed against E. coli, with inhibition diameters of 19.0 ± 1.0 mm. The same results were observed in the essential oil from cooperative C3. The inhibition diameters obtained for M. smegmatis (27.0 ± 1.1 mm) were higher than those obtained for B. subtilis (18.6 ± 2.1 mm) and E. coli (15.2 ± 1.1 mm). Many studies have demonstrated that rosemary essential oil provides effective antibacterial activity against various microorganisms [38,40]. Rosemary essential oil was slightly less sensitive against Gram-negative bacteria (E. coli) than they against Gram-positive bacteria [54]. Megzari et al. [55] observed that rosemary essential oil presented antibacterial activity against the same bacterial strain, with diameters of inhibition similar to those observed in our study. However, Bajalan et al. [56] stated that the essential oil of rosemary showed antibacterial activity against the four bacteria strains, especially E. coli, with an inhibition zone of 18.5 mm, followed by Staphylococcus aureus (14.6 mm), Klebsiella pneumoniae (13.9 mm), and Streptococcus agalactiae (13.1 mm). Abdellaoui et al. [57] revealed that the highest antibacterial activity was against E. coli, which was the most sensitive of the tested strains with the largest inhibition zone (18.4 mm). Chahboun et al. [58] showed that rosemary essential oil in the Taza region of Morocco was active against all strains, except for Pseudomonas aeruginosa, with the other strains having diameters between 12.0 and 22.0 mm.
The results of the minimum inhibition concentration are represented in Table 3. The essential oils showed important inhibitor activity against the microorganisms studied. In fact, the bacterial strains as a whole were inhibited at a 1/250 v/v concentration. M. smegmatis and B. subtilis were the most susceptible strain for all Eos, with growth ending at 1/500 (v/v). E. coli was the least sensitive microorganism, an MIC value of with 1/250 (v/v). All EOs exerted strong activity against M. smegmatis. Fadil et al. [59] reported that rosemary essential oil exhibited antimicrobial activity against Salmonella typhimurium, with an MIC value of 2% (v/v).The MIC value of rosemary essential oil reported in [60] against three bacterial strains (S.aureus, E.coli, and Pseudomonas aeruginosa) was 3.33, 1.67, and 6.67 µL/mL, respectively. According to [61], rosemary essential oil showed antimicrobial activity against E. coli, Salmonella choleraesuis, Listeria monocytogenes, and S. aureus, with an MIC value of 0.65, 0.65, 0.91, and 0.91 mg/mL, respectively. Table 3. Minimum inhibition concentration of S. rosmarinus essential oils tested. It can be concluded that the rosemary essential oil from cooperative C1 exhibited the best antimicrobial activity of the three cooperatives. This can be explained by a high amount of 1,8-cineole compound in rosemary essential oil from cooperative C1 (51.1 ± 1.8%) compared to the cooperatives C2 and C3 (46.2 ± 0.3 and 47.3 ± 0.8%, respectively), but the contributions of α-pinene and camphor compounds should also be noted. In addition, the analysis of the S. rosmarinus essential oil from cooperative C1 showed a high rate of camphor (7.3 ± 0.5%) compared to cooperatives C2 (5.8 ± 0.2%) and C3 (5.3 ± 0.1%). The essential oil from cooperative C2 represented the highest amount of α-pinene compound (11.6 ± 0.5%), followed by C3 (10.7 ± 0.5%) and C1, with a value of 9.9 ± 0.2%. The effect of synergy could also be the origin of this activity. Several authors demonstrated that chemical composition has an effect on antimicrobial activity, such as 1,8-cineole and α-pinene [62,63]. According to [64], the antibacterial activity against M. smegmatis of rosemary was most likely to occur when high quantities of 1,8-cineole were present. Moreover, Kovač et al. [65] demonstrated that the α-pinene compound exhibited antimicrobial activity against 16 strains of Campylobacter jejuni. Ojeda-Sana et al. [66] indicated that rosemary essential oil has a high antibacterial activity against Gram-positive bacteria (S. aureus and Enterococcus faecalis) and against the Gram-negative bacteria (E. coli and K. pneumoniae). The same authors also demonstrated that the α-pinene compound presented a broad antibacterial spectrum against Gram-negative and Gram-positive bacteria and that the 1,8-cineole compound showed antibacterial activity against pathogenic bacteria by causing a rupture of the cell membrane. It is known by the scientific community that the main compounds of EOs have a key role in antibacterial activity. However, the importance of minor components or a combination thereof could explain the complete activity of EOs against clinical and foodborne pathogens [64,67].

Cooperative Strains
The present study indicates that the most important period to exploit the essential oil of S. rosmarinus is in summer for the three studied cooperatives in Figuig Province. This is due to the abundance of the main components that affect antibacterial activity, including 1,8-cineole, α-pinene, and camphor.

Plant Material
The aerial parts of Salvia rosmarinus Spenn (S. rosmarinus) were collected from three cooperatives ( Table 4)   The present study indicates that the most important period to exploit the essential oil of S. rosmarinus is in summer for the three studied cooperatives in Figuig Province. This is due to the abundance of the main components that affect antibacterial activity, including 1,8-cineole, α-pinene, and camphor.

Plant Material
The aerial parts of Salvia rosmarinus Spenn (S. rosmarinus) were collected from three cooperatives ( Table 4)

Essential Oil Extraction
The aerial parts of the plants were dried in shade at ambient temperatures and then hydro distilled using a Clevenger-type apparatus [69]. This method is recommended for essential oil extraction because it is simpler and faster than other extraction methods, and it allows for a good yield and recovery of the compounds of the essential oil [70]. Each

Essential Oil Extraction
The aerial parts of the plants were dried in shade at ambient temperatures and then hydro distilled using a Clevenger-type apparatus [69]. This method is recommended for essential oil extraction because it is simpler and faster than other extraction methods, and it allows for a good yield and recovery of the compounds of the essential oil [70]. Each specimen was put in a 2 L glass flask containing 100 g of dry plant material and 800 mL of distilled water [71]. The mixture was then brought to a boil for three hours. The essential oil obtained was placed in an opaque green flask, dried using anhydrous sodium sulfate, and stored at 4 • C [72]. The yields of the EOs from the different samples were calculated using formula (1) indicated by [73]. All the experiments were carried out in triplicate.

Gas Chromatography/Mass Spectrometry Analysis
Gas chromatography analyses for all the samples were undertaken on a Hewlett-Packard (HP 6890) gas chromatograph equipped with an HP-5 capillary column (30 m × 0.25 mm, film thickness of 0.25 µm), an FID detector, and an injector fixed at 275 • C. The oven temperature was set at 50 • C for five min and was then raised to 250 • C at a rate of 4 • C/min. The carrier gas applied was N2 at 1.8 mL/min, employing a split mode of ratio: 1/50, flow: 72.1 mL/min. The samples were diluted to 1/50 in methanol, and the injection was added manually with an injected volume of 1.2 µL.
The chemical composition was realized by gas chromatography coupled with mass spectrometry (GC/MS). This analysis was performed on a Hewlett-Packard gas chromatograph (HP 6890) coupled with a mass spectrometer (HP 5973). The column used was a capillary column equipped with an HP-5MS (cross-linked 5% PHME siloxane) (30 m × 0.25 mm, film thickness 0.25 µm). The column temperature was set at 50 • C and raised to 250 • C at a rate of 2 • C/min. The carrier gas was He at 1.5 mL/min, and a split mode ratio was employed: 1/74.7, flow: 112 mL/min. The MS identity of the components was confirmed using the NIST 98 spectra library. The parameters exploited in the mass spectra were ionization voltages of 70 eV, the ion source temperature was 230 • C, and the scan mass range was 35-450 m/z. The identification of the components was also verified by the comparison of the compounds' elution order with their relative retention indices indicated in the literature. Each experiment was performed in triplicate.

Disc Diffusion Method
The antimicrobial activity for each essential oil was determined using a disc diffusion method based on that used in [74] to identify the antibacterial power of the oils. The bacterial suspension was adjusted to 10 8 CFU/mL and was then spread on a plate containing Muller-Hinton agar (MHA). Whatman paper sterile discs (6 mm in diameter) were soaked with 10 µL of essential oil samples and put on the inoculated MHA. In addition, the negative control was prepared in the same manner as the experimental test, and sterile distilled water was added rather than essential oil. A commercial disc of amoxicillin (30 µg and 86.2% purity) was used as a positive reference standard to determine the sensitivity of the tested strain. The diameters of the inhibition zone around the disks were measured after 24 h incubation at 37 • C. All assays were realized in triplicate.

Minimal Inhibitory Concentration
The minimal inhibitory concentration (MIC) of essential oils was determined in a 96-well microplate using the microdilution method, according to [75], with a slight modification. The EOs were diluted sequentially in Mueller-Hinton broth (MHB) with agar added at 0.15% w/v as an emulsifier. The final concentration of the essential oil was between 1/10 v/v and 1/2000 v/v. The 12th was regarded as a growth control, as it only contained the culture medium and the strain. Afterwards, 50 mL of bacterial suspension was supplemented in each well at a final concentration of 10 6 CFU/mL. After incubation at 37 • C for 18 h, the MIC was presented as the lowest essential oil concentration that showed a negative bacterial growth translated by a non-change in resazurin color. A positive MIC is detected by a reduction of blue dye resazurin to pink resorufin. Each experiment was performed in triplicate.

Statistical Analysis
To examine the significant differences between mean values, an analysis of variance (two-way ANOVA), and a Tukey's post hoc test were utilized (a probability of p ≤ 0.05 was considered statistically significant) using Origin Pro 2021 software. The percentage composition of 8 major compounds (1,8-cineole, camphor, α-pinene, camphene, β-pinene, myrcene, α-terpineol, and limonene) and the average yield of rosemary essential oil were calculated to determine the variable variation and the correlation between them, as well as to identify the similarities between samples (harvest period and cooperatives) with principal component analysis (PCA). For better visualization of the similarity between harvest periods, a hierarchical cluster analysis (HCA) was performed. This classification tool divided the groups according to the average yield and the 8 major components of rosemary essential oil. The PCA and the HCA were realized using SPSS (version 25) and JMP Pro 14 software.

Conclusions
In this study, we focused on explaining the variations in rosemary essential oil yield, chemical composition, and antibacterial activity in relation to cooperatives and harvest periods. The best yield was obtained in essential oils for plants from cooperative C1 for period 9 (April 2019). The highest amount of the 1,8-cineole component appeared in period 3 (May 2018). PCA and HCA indicated the presence of two groups: group A (winter), which was characterized by camphor and α-pinene, and group B (summer), which was characterized by a high quantity of 1,8-cineole and high yield of rosemary essential oil. These statistical analyses indicated that summer is the best harvest period to exploit the essential oil of this plant in Figuig Province. The antibacterial activity showed an inhibitory potential against all strains, especially among samples from the cooperative C1. Our findings could have a positive economic impact on the exploitation of rosemary in the Oriental region, especially during the best harvest periods, as they show the best yields and the oils richest in 1,8-cineole and α-pinene chemotypes. The use of chemometric tools has the potential to identify the compounds responsible for antibacterial activities and will also make it possible to exploit these results to determine the most active oil fractions against pathogenic strains. In this context, this study could be extended to other bacterial strains of medical and food interest.