Essential Oils Obtained from Aerial Eugenia punicifolia Parts: Chemical Composition and Antiproliferative Potential Evidenced through Cell Cycle Arrest

Essential oils (EOs) of the leaves of three Eugenia punicifolia specimens from two different Reservation Parks, namely Parque Nacional das Nascentes do Rio Parnaíba (EpNRP-I and EpNRP-II) and Parque Nacional da Chapada das Mesas (EpCM), in the state of Maranhão, Brazil, were extracted by hydrodistillation and investigated by gas chromatography coupled to mass spectrometry. Principal component and hierarchical cluster analyses indicated differences between the samples. Antiproliferative EOs activity was determined for U-251 (glioblastoma), MCF-7 (breast adenocarcinoma), NCI/ADR-RES (multidrug-resistant ovarian adenocarcinoma), OVCAR-3 (ovarian adenocarcinoma), HT-29 (colorectal adenocarcinoma), and HaCaT (non-tumor keratinocyte) cell lines applying the colorimetric method using 3-(4,5-dimethyl2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT) to determine the GI50 (50% growth inhibition) concentration. The extraction yields of the analyzed EOs were 0.58, 1.42 and 0.84%. The main constituents identified in two samples were α-pinene (49.75%), 1,8-cineole (13.77%) and α-terpineol (7.32%), and in the third sample, germacrene B (16.25%), (E)-caryophyllene (13.21%) and β-pinene (12.81%). The main GI50 results for sample EpNRP-I were noted for the U-251 (2.13 μg mL) and MCF-7 (6.72 μg mL) tumor lines. For the non-tumoral line HaCaT, the calculated GI50 was higher than the positive control comprising doxorubicin hydrochloride (13.35 μg mL). In addition, a flow cytometry analysis revealed that this same sample arrests the cell cycle of the MCF-7 line in the second interphase stage.


Introduction
Myrtaceae has more than 5600 species distributed between 130 to 140 genera. 1,2 It is native to the Cerrado and other biomes in the tropical and subtropical areas of South America, Australia, and Asia. In Brazil, it is considered one of the ten richest families of angiosperms and also one of the ten most expressive in the Cerrado region, where at least 14 of the 23 genera occurring in the country could be found. 1,3 It is also known for producing high content of terpenes, the main constituents of its essential oils. 4 Several research groups seek to study in a non-predatory way the potential of Brazilian biodiversity, in particular, the aromatic flora and its essential oils from the North and Northeast of Brazil. Studies are being carried out and adding knowledge about this rich flora, either to develop new technologies from natural products or to learn about the phytochemical characteristics of these species. This work is part of the effort that a group of researchers from federal universities in the North and Northeast of Brazil has been making, whose main objective is to investigate the potential of aromatic plants in the Amazon biome.
A database that lists several aromatic species collected in the nine states of the Brazilian Amazon has been published, 5 confirming not only the richness of the Amazonian flora but also its phytochemical potential. This study has been underway for more than 30 years, which reflects a significant example of how much remains to be discovered in the Brazilian biodiversity. The Cerrado biome in Brazil has been rapidly deforested, with about 2.2% of its area each year. Before, it represented about 2 million km 2 , almost a quarter of the total area of the country. Today, less than 1.0% of its original area is part of protected reserves in Brazil. Despite this, Cerrado is still responsible for about 12.0% of the known plant species. 6 Eugenia L., belonging to the Myrtaceae family, presents approximately 1000 species, with occurrence in Central and South America, and few in the African continent. They are shrubs or small trees with twigs glabrous or pubescent when young. 7 Essential oils (EOs) obtained from Eugenia species are noteworthy for displaying significant economic and pharmacological potential, such as cytotoxic, 1 antibacterial, 8 antifungal, 9 and insecticide activities. 10 The compounds present in the EOs of this genus, such as the monoterpenes αand β-pinene, limonene, and γ-terpinene, the sesquiterpenes germacrene B, α-humulene, and (E)-caryophyllene, as well the phenylpropanoid eugenol, have been described as displaying cytotoxic, and anticancer activities. 2,5,11,12 The antiproliferative activity of the EOs of nine Eugenia species has recently been described 13 highlighting promising results for E. cuspidifolia and E. tapacumensis in inhibiting cell migration and cell proliferation of HCT-116 colorectal cell line colonies. In addition, the EOs are not genotoxic and, therefore, do not cause DNA damage.
Based on the studies reported 2,3,14 for the Eugenia genus, it was observed that some common constituents existing in the genus essential oils have already shown antineoplastic activities, either isolated or in synergy with each other. One Eugenia species that has been the focus on work by our group is Eugenia punicifolia (Kunth) DC. (syn. E. ambigua O.Berg, E. benthamii O.Berg, E. fruticulosa DC., E. insipida Cambess., E. sancta DC., Myrtus punicifolia Kunth, among many others). This is an endemic species with a wide distribution in Brazil, occurring in the biomes of the Amazon and Atlantic Forest and the Cerrado and Caatinga of the Northeast. It is a shrub of up to 3 m and known as "cereja-da-praia", "cereja-do-campo", "murta", "pitanga-do-campo", and others with antinociceptive, anti-inflammatory and antibacterial properties described in literature, 12,15,16 but no one related to anticancer activity.
Therefore, studies on the chemical composition of EOs obtained from the aerial parts of E. punicifolia, allied to antiproliferative screening, may contribute to phytochemical and biological knowledge on Eugenia species. In this context, the aims of the present study are to evaluate the chemical composition and the potential technological application of Eugenia punicifolia EOs, aiming at the promising use of their chemical and biological properties in favor of the Brazilian society.

Essential oil extraction and yield calculation
Leaf samples (50 g) were submitted to hydrodistillation using a Clevenger-type apparatus (3 h) (LaborQuimi, São Paulo, Brazil). 2,5 The oils were dried over anhydrous sodium sulfate (Merck-Millipore, São Paulo, Brazil), and the percentage contents were calculated on basis of the dry weight of plant material. The moisture content of the samples was calculated using an Infrared Moisture Balance (Genaka, São Paulo, Brazil) for water loss measurement. The procedure of the moisture content was performed in duplicate.

Oil-composition analysis
The analysis of the EO samples was performed on a gas chromatograph coupled to mass spectrometer (GCMS)-QP2010 Ultra system (Shimadzu Corporation, Tokyo, Japan), equipped with the GCMS-Solution software containing the NIST 11, 17 FFNSC 2, 14 and Adams 18 libraries. An Rxi-5ms (30 m × 0.25 mm; 0.25 μm film thickness) silica capillary column (Restek Corporation, Bellefonte, PA, USA) was used. The conditions of analysis were: injector temperature of 250 °C; oven temperature programming of 60-240 °C (3°C min −1 ); the helium as carrier gas, adjusted to a linear velocity of 36.5 cm s −1 (1.0 mL min −1 ); split mode injection for 1.0 μL of the sample (oil 3.0 μL:hexane 500 μL); split ratio 1:20; ionization by electronic impact at 70 eV; ionization source and transfer line temperatures of 200 and 250 °C, respectively. The mass spectra were obtained by automatic scanning every 0.3 s, with mass fragments in the range of m/z 35-400. The retention index was calculated for all volatile components using a homologous series of C8-C40 n-alkanes (Sigma-Aldrich, Saint Louis, USA), according to the linear equation of Van Den Dool and Kratz. 19 The quantitative data regarding the volatile constituents were obtained using a GC 2010 series, coupled to flame ionization detector (FID), operated under similar conditions of the GC-MS system. The components were identified by comparing their retention indexes and mass spectra (molecular mass and fragmentation pattern) with those existing in the GCMS-solution system libraries.

Antiproliferative activity assay
The colorimetric assay 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT), Sigma-Aldrich, St. Louis, Missouri, USA), which indirectly evaluates cell viability through the activity of mitochondrial reductase enzyme, 20 was used to evaluate cell proliferation as described by Franco et al. 21 Briefly, 5 tumor cell lines: U-251 (glioblastoma), MCF-7 (adenocarcinoma), NCI/ ADR-RES (multi-drug resistant ovarian adenocarcinoma), PC-3 (prostate adenocarcinoma), OVCAR-3 (ovarian adenocarcinoma), HT-29 (colorectal adenocarcinoma) and one non-tumoral cell line HaCaT (keratinocyte), were seeded into 96-well plates (5 × 10 3 cells well −1 ) and allowed to grow overnight to adhere. The EO samples were dissolved in dimethyl sulfoxide (DMSO, Merck, Darmstadt, Germany) and diluted with Roswell Park Memorial Institute (RPMI) medium to final concentrations ranging from 1.6 to 100 μg mL −1 . The final concentration of DMSO in the culture medium was kept constant, below 0.1% (v/v). A control plate (T0 plate) was used to detect the basal number of cells. As a positive control, doxorubicin was used at the same concentrations. After 48 h, 100 μL MTT reagent (0.25 mg mL −1 per well) was added followed by incubation for 4 h at 37 °C. The absorbance (value) at 470 nm wavelength was measured using a microplate reader (Epoch BIOTEK ® , Winooski, Vermont, USA). The GI 50 (concentration required to inhibit 50% cell growth) values were determined by nonlinear regression analysis using software Origin 8.0. 22 The selectivity index (SI) was also calculated. SI is an analysis that allows measuring how more selective the compound is to tumor cells than non-tumor cells. This analysis can be obtained by dividing the GI 50 value of a non-tumor cell line by a tumor cell line. In this assay, the SI was calculated by the following equation 1: (1) MCF-7 was the cell line chosen to continue the studies, according to the sensitivity of the sample on the antiproliferative activity assay.

Soft agar clonogenic assay
A mix containing MCF-7 cells (1 × 10 5 ) agar, and RPMI medium supplemented twice with 20% of FBS (fetal bovine serum) + 1% P/S (penicillin/streptomycin) per well was seeded in a 6-well plate and incubated at 37 °C in a 5% CO 2 incubator. Every three days the samples were treated with DMSO (1%) and EpNRP-I (7 µg mL −1 ). After 21 days in culture, the medium was removed, the wells were washed with distilled water and the colonies were then fixed with formaldehyde (4%) and stained with crystal violet (0.005%). The wells were photographed, and the images analyzed in the ImageJ software, 23 for quantification of the colonies.

Cell cycle analysis
This experiment was done using Guava ® Cell Cycle reagent (Austin, TX, USA), as previously described. 24 Briefly, MCF-7 cells (5 × 10 4 cells well −1 ) in 12-well plates with the complete medium were incubated at 37 °C with 5% CO 2 for 24 h. Then, the complete medium was replaced by RPMI medium without fetal bovine serum (FBS) for cell cycle synchronization and cells were incubated for another 24 h. After that, the cells were treated with DMSO (negative control), etoposide (positive control, 5.89 μg mL −1 or 10 mM) and EpNRP-I (7 μg mL −1 ) during 24 h. The cells were harvested, collected, centrifuged (5 min, 2500 rpm), and the supernatant discarded. After fixation with 70% cold ethanol (24 h, 4 °C), each cell suspension was centrifuged and washed with phosphate-buffered saline (PBS), the supernatant was discarded, and then the Guava Cell Cycle reagent (200 µL cell −1 suspension) was added. After 20 min at room temperature in the dark, each cell suspension was analyzed (5000 events per replicate) by flow cytometry. Using the Guava Cell cycle ® software, 25 the subpopulations at G1, S, and G2/M phases of the cell cycle were quantified in percentage. The analyses were done in biological triplicate of two experiments.

Statistical analysis
Results were expressed as mean ± standard deviation (SD) of three independent experiments conducted in duplicate. Statistical analyses were performed with the GraphPad Prism 5 software. 26 t-Test and analysis of variance (ANOVA) followed by Bonferroni post hoc test were used, and p-values less than 0.05 were considered significant. For the cytotoxic activity assay, the linear regression of the curves was obtained using the mean growth percentage and, calculated with Origin ® software. 22 The oils components with a percentage higher than 1% of the total oil were subjected to hierarchical cluster analysis (HCA) and principal component analysis (PCA) using SPSS v22.0 software. 27 In the case of HCA, the dendrogram was produced using Ward's method of hierarchical clustering with squared Euclidean distance between oil samples.
It is important to take into account that geographic and environmental factors influence plant growth in different locations, while soil water and nutrient availability, sunlight, circadian and seasonal effects, and herbivory, among others, may alter EO composition. 29 The EpNRP-I and EpNRP-II oils samples showed chemical characteristics quite different from the EpCM sample. However, most primary mono-and sesquiterpenes constituents of the oils from E. punicifolia were the same, as α-pinene, germacrene B, 1,8-cineole, (E)-caryophyllene, β-pinene, (E)-β-ocimene, α-terpineol, bicyclogermacrene, limonene, and germacrene D, all of them arranged by their biosynthetic pathways and in an increasing elution order as can be seen in Table 1, and with the chemical structures displayed in Figure 1. In the mevalonate pathway, acyclic monoterpenes, such as (E)-β-ocimene, result from the conversion of geraniol or linalool, leading to the other monoterpene constituents identified in E. punicifolia oils. 30 (E)-β-Ocimene was found only in EpCM oil, suggesting that the plant sampled in Chapada das Mesas was younger than the samples of EpNRP-I and EpNRP-II plants collected in the Rio Parnaíba region, since it is a precursor molecule in the biosynthesis of other monoterpenes identified in the oils of E. punicifolia.
The oils of two E. punicifolia specimens collected at Manaus-Caracaraí Road and Manaus-ltacoatiara Road, Amazonas State, Brazil, displayed three sesquiterpene hydrocarbons, (E)-caryophyllene, α-humulene, and δ-cadinene, as primary constituents. 31 In the oils of E. punicifolia sampled at the localities of Serra Negra and Madre de Deus, Pernambuco State, Brazil, predominated linalool, (E)-caryophyllene and α-terpineol. 28 Also, the oil of an E. punicifolia specimen that is growing wild in sandbanks of Maracanã municipality, Pará State, Brazil, showed (E)-caryophyllene, bicyclogermacrene, germacrene D, and (E)-β-ocimene as the main volatile components. 28 Additionally, the oil of another E. punicifolia specimen occurring at Restinga area of Macaé City, Rio de Janeiro, Brazil, presented α-cadinol, 10-epi-eudesmol, paradisiol, 7-epi-α-selinene, (E)-caryophyllene, and β-selinene as major compounds. 32 The application of an HCA resulted in the dendrogram presented in Figure 2a, which indicates the composition of elements representing similarities between EpNRP-1 and EpNRP-1, while indicating that the EpCM sample was different from both. The PCA main components (PC1 and PC2) accounted for 100% of the total data variability, where PC1 explained 82.87% and component PC2, 17.13% (Figure 2b). Regarding the EpNRP-I and II samples, the monoterpene hydrocarbon α-pinene was the major constituent, present at approximately 50% of both samples. The EpCM sample, on the other hand, presented the sesquiterpene hydrocarbon germacrene B as the constituent present at the highest concentration in the EO, at approximately 16%. These data explain the PCA and HCA results.

Antiproliferative activity of EOs samples
The EpNRP-II oil sample was not subjected to the antiproliferative analysis, due to a very similar volatile composition when compared to EpNRP-I and reduced availability. α-Pinene, the monoterpene hydrocarbon present in the EpNRP-I oil sample as its principal constituent (49.8%) and present in the EpCM oil sample at low concentrations (5.2%), displayed significant antiproliferative activity in some tumor cells, as previously reported. 3 Based on this information and the significant α-pinene content in E. punicifolia oils, complementary antiproliferative bioassays for the EpNRP-I and EpCM oil samples were carried out, including a third sample consisting of a mixture of EpCM oil and a commercial α-pinene standard (99%), in order to represent a similar percentage value to that of α-pinene in the EpNRP-I oil sample. The third sample doped with the α-pinene standard was named EpCM-d. Additionally, the α-pinene standard was also separately subjected to bioassays to verify if it was directly responsible for antiproliferative activity. The antiproliferative activities of the samples were plotted on graphs representing cell growth percentages according to sample concentration. When the growth parameter remained at zero or close to zero in the graphs, a cytostatic sample, where no cell development occurs, was characterized. On the other hand, a cytocidal sample displays negative cell growth. The growth inhibition (GI 50 ) results of the E. punicifolia EpNRP-I, EpCM, and EpCM-d oil samples, as well as for the α-pinene standard and doxorubicin hydrochloride used as the positive control, are displayed in Figures 3 and 4.
EpNRP-I oil sample showed the lowest GI 50 values for all tumor cell lines evaluated (Table 2), with glioblastoma (U-251, GI 50 2.13 μg mL −1 ) and breast adenocarcinoma (MCF-7, GI 50 6.72 μg mL −1 ) cell lines being the most sensitive to the treatment, with SI greater than 2.0 ( Table 3). According to the National Cancer Institute (NCI), significant SI values are considered to be greater than or equal to 2.0. 33 Doxorubicin presents SI values greater than 2.0 only for MCF-7 tumor cells. In addition, EpNRP-I showed cytotoxic effects for the non-tumoral keratinocytes line (HaCaT), however, its effect was 43 times smaller when compared to the positive control of the assay, the chemotherapy drug doxorubicin hydrochloride ( Table 2).
Although conventional chemotherapy has been successful to some extent, the main drawbacks of chemotherapy are its poor bioavailability, high-dose requirements and adverse side effects, proving the need and importance of researching new molecules with more effective anticancer activity and with no or minimal side effects.
Considering the low GI 50 value for MCF-7, the SI = 2.0, together with the fact that this tumor type is the one that most affects Brazilian women and the existing treatments confer so many adverse effects to patients, this line was selected to continue the investigation of the antiproliferative activity of EpNRP-I. Although the focus of the study was the investigation of the antiproliferative activity of essential oil in breast tumor lineage, the result obtained for the glioblastoma lineage is quite interesting. Due to the presence of the blood-brain barrier, the implementation of new approaches for the treatment of glioblastoma is much more difficult than in the case of other tumors. The selectivity of this barrier reduces the bioavailability of hydrophilic drugs in the central nervous system. The volatile constituents of essential oils are believed to easily cross the blood-brain barrier due to their small molecular size and lipophilic nature. 34 Therefore, we encourage the continuity of studies with the U-251 cell line.
The EpCM oil sample doped with the α-pinene standard and the α-pinene standard sample alone did    not present low GI 50 , in contrast to that observed for EpNRP-I. Although α-pinene has been reported 3,35 as an important monoterpenoid with activity against various human cancer cell lines through apoptosis stimulation by mitochondrial function disruption, reactive oxygen species (ROS) formation, caspase-3 activation, heterochromatin aggregation, DNA disintegration, and cell surface phosphatidylserine exposure, our results suggest that the active compound responsible for the antiproliferative action may not be α-pinene, or at least not only α-pinene. Therefore, the hypothesis of a synergistic effect, resulting from diverse oil constituent action, is much more likely. 36 Another important result obtained in this assay comprised GI 50 values for HaCaT (keratinocytes). EpNRP-I was 43-fold less toxic to this non-tumoral cell line compared to doxorubicin, a chemotherapeutic drug widely applied in the medical clinic context. This result is very interesting, as it may predict a selective effect for the assessed tumoral cell line. However, despite advances in chemotherapy treatments, many patients do not undergo this type of intervention, mainly due to the effects caused by non-tumor cells, proving the need and importance of researching new molecules with more effective anticancer activity and minor effects. The antiproliferative activity of the extract was further confirmed by a soft agar clonogenic assay. This assay measures the potential of cells to expand into colonies, which reflects the ability of cells to proliferate. 37 As presented in Figure 5, there was a reduction of 63% in the number of cell colonies after treatment with EpNRP-I. This result corroborates the findings obtained through the MTT assay and highlights EpNRP-I potential antiproliferative activity.

Cell cycle analysis
It was possible to observe an increase from 9.3 to 20.0% in the number of MCF-7 (breast adenocarcinoma) cells in the second phase of interphase (S phase) when the EpNRP-I oil sample of E. punicifolia was applied. Also, the cells treated with etoposide, the positive control, increased from 9.3 to 26.7%, considering the same treatment time. 38 Consequently, in the analysis of the EpNRP-I oil sample, there was a reduction from 87.5 to 78.0% in the number of cells in the first phase of the interphase (G1 interval), and in the positive control, this reduction was from 87.5 to 71.8%. The results were compared to negative control DMSO (see Figure 6).

Conclusions
The findings reported herein indicate differences in the chemical composition of the EO samples obtained from aerial Eugenia punicifolia parts from two different locations. In addition, an expressive antiproliferative potential was noted for the EpNRP-I oil sample. However, no relationship between this activity and a single constituent was observed. A flow cytometry analysis indicates that the sample induces cell cycle arrest in the genetic material duplication phase in cells belonging to MCF-7 cell line. The results of the present study showed the potential chemical composition and the technological application of Eugenia punicifolia essential oil, in favor of Brazilian society.