Comparative Evaluation of Apoptosis Induction Using Needles, Bark, and Pollen Extracts and Essential Oils of Pinus eldarica in Lung Cancer Cells

: Lung cancer is one of the leading causes of cancer-related mortality worldwide. Although effective clinical drugs for treating advanced stages are available, interest in alternative herbal medicines has gained momentum. Herbal extracts are potent antioxidants that reportedly inhibit the growth of various cancer cell lines. In the present study, we investigated the effects of essential oils and hexane, methanolic, and aqueous extracts, obtained from various parts (bark, needles, and pollen) of Pinus eldarica against human lung cancer (A549) cells. First, the DPPH radical scavenging activities of P. eldarica extracts and essential oils were examined, which revealed that methanolic extracts presented higher antioxidant activity than the other extracts and essential oils. Next, A549 cells were exposed to various concentrations of the extracts and essential oils for 48 h. P. eldarica extracts/essential oil-treated lung cancer cells demonstrated a signiﬁcant decrease in cell proliferation, along with an induction of apoptotic cell death, particularly, the pollen hexane extract, bark essential oil, and methanolic needle extract showed superior results, with IC 50 values of 31.7, 17.9, and 0.3 µ g/mL, respectively. In the cell cycle analysis, treatment of A549 cells with the methanolic needle and pollen hexane extracts led to apoptosis and accumulation of cells in the sub-G1 phase. Further, exposure to the bark essential oil and methanolic needle extract decreased the cell population in the G2/M phase. Notably, treatment with the pollen hexane extract, bark essential oil, and methanolic needle extract resulted in caspase-3 activation, poly (ADP-ribose) polymerase cleavage, Bcl-2 downregulation, and Bax and p53 regulation in A549 cells. Furthermore, these extracts and essential oils decreased the migration, and colony formation of A549 cells. These ﬁndings provide experimental evidence for a new therapeutic effect of P. eldarica against human lung cancer and suggest P. eldarica as a potential chemopreventive natural resource for developing novel cancer therapeutics.


Introduction
Cancer is a term that encompasses a large family of diseases in which abnormal cells undergo uncontrolled division and can spread to other tissues through the blood and lymphatic systems. Cancer is the second leading cause of death worldwide and was reportedly responsible for 9.6 million deaths in 2018 [1]. Among these, lung cancer is the most commonly diagnosed cancer (11.6% of all cases), ranking first in terms of reported cancer deaths (18.4% of all cancer deaths) [1]. Throughout the history of human life, natural products derived from living organisms such as plants have been utilized as a primary source of medicinal agents [2]. Plants provide an extensive reservoir of therapeutic agents, representing a significant structural variety and remarkable chemical entities in modern medicine. Numerous studies have also reported that natural products derived from plants possess anticancer properties [3].
Gymnosperms are seed-bearing vascular plants, including cycads, ginkgo, yews, and conifers, in which the ovules or seeds are not enclosed within the ovary. Pinostrobin, pycnogenol, enzogenol, and pinocembrin are important anticancer agents of gymnosperm origin. Paclitaxel (taxol), isolated from the bark of yew trees (Taxus sp.), is reportedly cytotoxic to human prostate, ovarian, breast, and lung cancer cells [4,5]. The Pinus genus is the largest in the family Pinaceae [6]. These coniferous, evergreen, and typically monoecious trees are cultivated in most areas of the Northern Hemisphere [7]. Pinus eldarica, also named Tehran pine, is native to the South Caucasus between Europe and Asia and has been extensively grown in Iran. In traditional medicine, different parts of this plant (bark, needles, and nuts) have been used for the treatment of several conditions, including skin wounds, bronchial asthma, dermatitis, and skin irritation, and have reportedly shown antibacterial, antimutagenic, antitumor, anti-inflammatory, and glucose-lowering activities [8][9][10][11]. Essential oils and extracts from pine needles and bark are rich in terpenoids such as caryophyllene, limonene, α-pinene, and phenolic acids [12][13][14]. Pine extracts act as antioxidants by reacting with free radicals, thus protecting DNA against oxidative damage and apoptosis induced by hydroxyl radicals [15][16][17]. Phytochemical evidence has revealed that pine extracts contain catechin and taxifolin, two polyphenolic compounds that afford a potent anticancer product, Pycnogenol ® [8,18]. Apoptosis plays an important role in modulating immunological competency, development, and homeostasis through changes in cell morphology, chromatin condensation, DNA fragmentation, caspase-3 activation, Bax and p53 elevation, Bcl-2 reduction, and poly (ADP-ribose) polymerase (PARP) cleavage [19]. For instance, a pinecone extract of Pinus koraiensis has revealed a beneficial effect in cancer prevention by inducing apoptosis in lung cancer cells through caspase-3 activation [20]. Therefore, induction of apoptosis in cancer cells with plant-derived products can be employed as an effective strategy against tumor development.
Several studies have reported the pharmacological properties and chemical diversity of P. eldarica extracts; however, the antineoplastic effects of P. eldarica against human lung cancer cells and their underlying mechanisms of action are yet to be reported. In the present study, we aimed to evaluate the antioxidant, cytotoxic, and proapoptotic effects of hexane, methanolic, and aqueous extracts, and essential oils, of various parts (bark, needles, and pollen) of P. eldarica on lung cancer (A549) cells, with the possible elucidation of cellular and molecular mechanisms. The results of this study could help pave the way for developing novel cancer therapeutic agents, opening new horizons for the application of natural products in lung cancer therapy.

DPPH Radical Scavenging Assay
The 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical scavenging assay was performed to determine the antioxidant activity of extracts and essential oils obtained from the various parts (bark, needle, and pollen) of P. eldarica (Figure 1). In the DPPH assay, a higher DPPH radical scavenging value indicates a higher antioxidant activity. The DPPH radical scavenging activity of methanolic extracts in all groups were significantly higher than that of other extracts and essential oils (p < 0.01). Among the active methanolic extracts, the methanolic bark extract (mean value 93.9%) demonstrated significantly higher antioxidant activity (p < 0.01) than the methanolic pollen extract (mean value 77.57%) and no significant difference with methanolic needle extract (mean value 83.90%).
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 18 of other extracts and essential oils (p < 0.01). Among the active methanolic extracts, the methanolic bark extract (mean value 93.9%) demonstrated significantly higher antioxidant activity (p < 0.01) than the methanolic pollen extract (mean value 77.57%) and no significant difference with methanolic needle extract (mean value 83.90%).

Figure 1.
DPPH radical scavenging (%) of the four different extracts (hexane, methanolic, aqueous extracts, and essential oils (each 1 mg/mL)) prepared from the Pinus eldarica bark, needles, and pollen. Data are expressed as mean ± standard deviation (SD), representative of three independent experiments.

In Vitro Evaluation of Cytotoxic Effects of P. eldarica Extracts and Essential Oils on A549 Lung Cancer Cells
Cytotoxicity of P. eldarica extracts and essential oils on A549 lung cancer cells was measured using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. We observed >80% cell death in A549 cells treated with 100 μg/mL of all hexane extracts, the methanolic needle extract, and needle essential oil ( Figure 2). The half-maximal inhibitory concentration (IC50) of the methanolic needle extract at 0.3 μg/mL revealed its high toxicity, comparable with that of doxorubicin (DOX; a positive control). Among pollen fractions, the hexane extract with an IC50 value of 31.7 μg/mL was the most toxic; among bark fractions, the essential oil was the most toxic, with an IC50 value of 17.9 μg/mL (Table 1). Based on these findings, the methanolic needle extract, hexane pollen extract, and bark essential oil were selected for further investigation. DPPH radical scavenging (%) of the four different extracts (hexane, methanolic, aqueous extracts, and essential oils (each 1 mg/mL)) prepared from the Pinus eldarica bark, needles, and pollen. Data are expressed as mean ± standard deviation (SD), representative of three independent experiments.

In Vitro Evaluation of Cytotoxic Effects of P. eldarica Extracts and Essential Oils on A549 Lung Cancer Cells
Cytotoxicity of P. eldarica extracts and essential oils on A549 lung cancer cells was measured using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. We observed >80% cell death in A549 cells treated with 100 µg/mL of all hexane extracts, the methanolic needle extract, and needle essential oil ( Figure 2). The half-maximal inhibitory concentration (IC 50 ) of the methanolic needle extract at 0.3 µg/mL revealed its high toxicity, comparable with that of doxorubicin (DOX; a positive control). Among pollen fractions, the hexane extract with an IC 50 value of 31.7 µg/mL was the most toxic; among bark fractions, the essential oil was the most toxic, with an IC 50 value of 17.9 µg/mL (Table 1). Based on these findings, the methanolic needle extract, hexane pollen extract, and bark essential oil were selected for further investigation.

Figure 2.
Cytotoxic activity of the four different fractions (hexane, methanolic, aqueous extracts, and essential oils) prepared from Pinus eldarica (A: Bark, B: Needles, and C: Pollen) against the A549 cell line after 48 h of treatment with various concentrations of essential oil, and methanolic, aqueous, and hexane extracts. Data are expressed as mean ± standard deviation (SD), representative of three independent experiments. * p < 0.05 and ** p < 0.01 vs. control.

Evaluating Apoptosis Induction by DAPI Staining
4',6-Diamidino-2-phenylindole (DAPI) staining was performed to evaluate apoptosis induction in the A549 cell line treated with P. eldarica extracts, essential oils, and DOX (as a positive control), based on IC50 values. The degree of nuclear fragmentation, as an indicator of apoptotic cells, is shown in Figure 3. Compared with untreated control cells, the treated A549 cells revealed increased nuclear fragmentation based on the morphology assessment. A large number of apoptotic cells with high DNA fragmentation were observed following treatment with the methanolic needle extract (Figure 3d). To show the apoptotic cell population, apoptosis analysis was performed by flow cytometry in the next step. . Cytotoxic activity of the four different fractions (hexane, methanolic, aqueous extracts, and essential oils) prepared from Pinus eldarica (A: Bark, B: Needles, and C: Pollen) against the A549 cell line after 48 h of treatment with various concentrations of essential oil, and methanolic, aqueous, and hexane extracts. Data are expressed as mean ± standard deviation (SD), representative of three independent experiments. * p < 0.05 and ** p < 0.01 vs. control.

Evaluating Apoptosis Induction by DAPI Staining
4 ,6-Diamidino-2-phenylindole (DAPI) staining was performed to evaluate apoptosis induction in the A549 cell line treated with P. eldarica extracts, essential oils, and DOX (as a positive control), based on IC 50 values. The degree of nuclear fragmentation, as an indicator of apoptotic cells, is shown in Figure 3. Compared with untreated control cells, the treated A549 cells revealed increased nuclear fragmentation based on the morphology assessment. A large number of apoptotic cells with high DNA fragmentation were observed following treatment with the methanolic needle extract (Figure 3d). To show the apoptotic cell population, apoptosis analysis was performed by flow cytometry in the next step.

Analysis of Apoptosis using Flow Cytometry
A549 cells were exposed to IC50 concentrations of the hexane pollen extract, bark essential oil, and methanolic needle extract for 48 h. Then, to investigate the rate of apoptosis and necrosis, both treated and untreated control cells were stained with Annexin V/propidium iodide (PI) and analyzed through flow cytometry ( Figure 4). Quantitative results of Annexin V/PI revealed increased amount of late apoptotic cells 23.4% and early apoptotic cells 43.2% in the methanolic needle extract-treated group compared to control (p < 0.01). The hexane pollen extract-treated group showed considerably more viable cells (70.1%) than the methanolic needle extract-and bark essential oil-treated groups, where viable cells were reduced to 33.3% and 43.4%, respectively (p < 0.05).

Analysis of Apoptosis Using Flow Cytometry
A549 cells were exposed to IC 50 concentrations of the hexane pollen extract, bark essential oil, and methanolic needle extract for 48 h. Then, to investigate the rate of apoptosis and necrosis, both treated and untreated control cells were stained with Annexin V/propidium iodide (PI) and analyzed through flow cytometry ( Figure 4). Quantitative results of Annexin V/PI revealed increased amount of late apoptotic cells 23.4% and early apoptotic cells 43.2% in the methanolic needle extract-treated group compared to control (p < 0.01). The hexane pollen extract-treated group showed considerably more viable cells (70.1%) than the methanolic needle extract-and bark essential oil-treated groups, where viable cells were reduced to 33.3% and 43.4%, respectively (p < 0.05).

Scratch Migration Assay
A scratch migration assay was performed to observe the anticancer effects of P. eldarica extracts and essential oil on the migration of A549 cells. Wound areas in A549 cells after 24 and 48 h of treatment with the hexane pollen extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), methanolic needle extract (0.3 µg/mL), and the control group (no treatment) were measured ( Figure 5A). The cells incubated with 0.3 µg/mL of methanolic needle extract inhibited cell migration across the initial scratch zone. As shown in the migration rate chart ( Figure 5B), after 24 and 48 h, the wound area in the pollen hexane extractand bark essential oil-treated cells was greater than that of cells incubated with the control medium, indicating that these extracts decreased cell migration when compared with the control group (p < 0.05).

Scratch Migration Assay
A scratch migration assay was performed to observe the anticancer effects of P. eldarica extracts and essential oil on the migration of A549 cells. Wound areas in A549 cells after 24 and 48 h of treatment with the hexane pollen extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), methanolic needle extract (0.3 µg/mL), and the control group (no treatment) were measured ( Figure 5A). The cells incubated with 0.3 µg/mL of methanolic needle extract inhibited cell migration across the initial scratch zone. As shown in the migration rate chart ( Figure 5B), after 24 and 48 h, the wound area in the pollen hexane extract-and bark essential oil-treated cells was greater than that of cells incubated with the control medium, indicating that these extracts decreased cell migration when compared with the control group (p < 0.05).

Cell Cycle Analysis and Apoptosis Evaluation
Treatment of A549 cells for 48 h with the methanolic needle extract and pollen hexane extract resulted in apoptosis induction and cell accumulation in the sub-G1 phase, which were comparable with those observed in the control group ( Figure 6A,B). Exposure of cancer cells to the bark essential oil and methanolic needle extract decreased the cell population in the G2/M phase.

Cell Cycle Analysis and Apoptosis Evaluation
Treatment of A549 cells for 48 h with the methanolic needle extract and pollen hexane extract resulted in apoptosis induction and cell accumulation in the sub-G1 phase, which were comparable with those observed in the control group ( Figure 6A,B). Exposure of cancer cells to the bark essential oil and methanolic needle extract decreased the cell population in the G2/M phase.

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on Bax/Bcl 2 Ratio
Next, to investigate the effects of P. eldarica extracts and essential oil on the Bcl-2 family in A549 cells, the expression of Bcl-2 (antiapoptotic) and Bax (proapoptotic) proteins was detected via Western blot analysis, and the Bax/Bcl-2 ratio was determined (Figure 7). A549 cells were treated with the IC50 concentration of the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), methanolic needle extract (0.3 µg/mL), and control group without treatment for 48 h, revealing a decrease in Bcl-2 and increase in Bax protein expression levels ( Figure 7A). The Bax/Bcl-2 ratio ( Figure 7B) increased by 6.8-and 27.2-times in the pollen hexane extract-and methanolic needle extract-treated cells, respectively, when compared with that in the control cells.

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on Caspase-3 Activation and PARP Cleavage
Treatment of A549 cells with IC50 concentrations of P. eldarica extracts and essential oil for 48 h significantly increased caspase-3 cleavage (p < 0.05) (Figure 8). PARP is the most crucial protein cleaved after caspase-3 activation. The treatment of cells with 0.3

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on Bax/Bcl 2 Ratio
Next, to investigate the effects of P. eldarica extracts and essential oil on the Bcl-2 family in A549 cells, the expression of Bcl-2 (antiapoptotic) and Bax (proapoptotic) proteins was detected via Western blot analysis, and the Bax/Bcl-2 ratio was determined ( Figure 7). A549 cells were treated with the IC 50 concentration of the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), methanolic needle extract (0.3 µg/mL), and control group without treatment for 48 h, revealing a decrease in Bcl-2 and increase in Bax protein expression levels ( Figure 7A). The Bax/Bcl-2 ratio ( Figure 7B) increased by 6.8and 27.2-times in the pollen hexane extract-and methanolic needle extract-treated cells, respectively, when compared with that in the control cells.

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on Bax/Bcl 2 Ratio
Next, to investigate the effects of P. eldarica extracts and essential oil on the Bcl-2 family in A549 cells, the expression of Bcl-2 (antiapoptotic) and Bax (proapoptotic) proteins was detected via Western blot analysis, and the Bax/Bcl-2 ratio was determined (Figure 7). A549 cells were treated with the IC50 concentration of the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), methanolic needle extract (0.3 µg/mL), and control group without treatment for 48 h, revealing a decrease in Bcl-2 and increase in Bax protein expression levels ( Figure 7A). The Bax/Bcl-2 ratio ( Figure 7B) increased by 6.8-and 27.2-times in the pollen hexane extract-and methanolic needle extract-treated cells, respectively, when compared with that in the control cells.

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on Caspase-3 Activation and PARP Cleavage
Treatment of A549 cells with IC50 concentrations of P. eldarica extracts and essential oil for 48 h significantly increased caspase-3 cleavage (p < 0.05) (Figure 8). PARP is the most crucial protein cleaved after caspase-3 activation. The treatment of cells with 0.3 µg/mL of the methanolic needle extract and 31.7 µg/mL of the pollen hexane extract induced PARP cleavage, compared with those observed in the control group (p < 0.05).

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on p53
Western blot analysis (Figure 9) revealed that protein expression levels of p53 were significantly higher in groups treated with the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), and methanolic needle extract (0.3 µg/mL) than in the control group.

Colony Formation Assay
We investigated the colony formation ability of A549 cells in the presence or absence of P. eldarica extracts and essential oil for 14 days. The colony numbers were significantly reduced after treatment with IC50 concentrations of the pollen hexane extract, bark essential oil, and methanolic needle extract ( Figure 10A). Compared with those in the control group, the colony numbers in the methanolic needle extract-treated group were reduced by over 80% (Figure 10B).

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on p53
Western blot analysis (Figure 9) revealed that protein expression levels of p53 were significantly higher in groups treated with the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), and methanolic needle extract (0.3 µg/mL) than in the control group.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 18 µg/mL of the methanolic needle extract and 31.7 µg/mL of the pollen hexane extract induced PARP cleavage, compared with those observed in the control group (p < 0.05).

Western Blot Analysis: Effects of P. eldarica Extracts and Essential Oil on p53
Western blot analysis ( Figure 9) revealed that protein expression levels of p53 were significantly higher in groups treated with the pollen hexane extract (31.7 µg/mL), bark essential oil (17.9 µg/mL), and methanolic needle extract (0.3 µg/mL) than in the control group.

Colony Formation Assay
We investigated the colony formation ability of A549 cells in the presence or absence of P. eldarica extracts and essential oil for 14 days. The colony numbers were significantly reduced after treatment with IC50 concentrations of the pollen hexane extract, bark essential oil, and methanolic needle extract ( Figure 10A). Compared with those in the control group, the colony numbers in the methanolic needle extract-treated group were reduced by over 80% (Figure 10B).

Colony Formation Assay
We investigated the colony formation ability of A549 cells in the presence or absence of P. eldarica extracts and essential oil for 14 days. The colony numbers were significantly reduced after treatment with IC 50 concentrations of the pollen hexane extract, bark essential oil, and methanolic needle extract ( Figure 10A). Compared with those in the control group, the colony numbers in the methanolic needle extract-treated group were reduced by over 80% ( Figure 10B).

Discussion
P. eldarica needles and bark are rich in polyphenols, tannins, and terpenoids. Essential oils and extracts of pine needles and bark have various physiological and pharmacological activities. Extracts and essential oils of P. eldarica are known to possess antioxidant, antiinflammatory, anticancer, and apoptotic properties, correlating with their effects on cyclooxygenase activity, nitric oxide synthesis, and apoptotic protein regulation [16,21,22]. Previous studies show that needles essential oil of P. eldarica was rich in D-germacrene and main component in the bark was limonene. P. eldarica pollen essential oil was rich in α-pinene [12,23]. Iravani et al. demonstrated that bark extract of P. eldarica had high amounts of polyphenolic compounds such as catechin, ferulic acid, and taxifolin [18].
In the study of Sarvmeili et al., MTT-based cytotoxic effects of essential oils and hydroalcoholic extracts of P. eldarica leaf and bark were evaluated on MCF-7 and HeLa cells [24]. Since this plant is native to Iran, it was necessary to conduct this more extensive, confirm the initial data, and identify molecular pathways. Therefore, we sequentially extracted P. eldarica from hydrophobic to hydrophilic compounds by hexane, methanol, and water extracts and essential oils from three parts of P. eldarica. Then, the most common lung cancer cell line, A549 was chosen and cellular and molecular mechanisms of extracts and essential oils were evaluated.
The results showed that the scavenging effects of methanolic extracts were higher than those of the other extracts and essential oils. These results indicate that methanolic extracts possess considerably higher antioxidant activities, and among them, the methanolic bark and needle extracts demonstrated high DPPH radical scavenging activity (93.9% and 83.90%, respectively). The hydroalcoholic extract of P. massoniana bark had the DPPH radical scavenging activity in the study of Yu et al. [25]. According to Dudonne et al., these extracts are potential natural antioxidant sources [26].
Patient-derived cancer cell lines could be used for the study of the effects of drugs on cancer cells. Scientists use different methods, such as monolayer in vitro assays, 3D cultures, bioreactors, or in vivo evaluations with immunocompromised animals [27,28]. In the present study, we used a monolayer A549 cell line to evaluate the biological effects of essential oils and extracts. For this purpose, the cytotoxic effects of essential oils and hexane, methanolic, and aqueous extracts of P. eldarica bark, needles, and pollen were evaluated using the A549 lung cancer cell line. Based on the MTT assay results, the cell mortality rate in all groups increased in a dose-dependent manner, and the methanolic needle extract, with an IC50 value of 0.3 µg/mL, presented the highest cytotoxicity when compared

Discussion
P. eldarica needles and bark are rich in polyphenols, tannins, and terpenoids. Essential oils and extracts of pine needles and bark have various physiological and pharmacological activities. Extracts and essential oils of P. eldarica are known to possess antioxidant, anti-inflammatory, anticancer, and apoptotic properties, correlating with their effects on cyclooxygenase activity, nitric oxide synthesis, and apoptotic protein regulation [16,21,22]. Previous studies show that needles essential oil of P. eldarica was rich in D-germacrene and main component in the bark was limonene. P. eldarica pollen essential oil was rich in α-pinene [12,23]. Iravani et al. demonstrated that bark extract of P. eldarica had high amounts of polyphenolic compounds such as catechin, ferulic acid, and taxifolin [18].
In the study of Sarvmeili et al., MTT-based cytotoxic effects of essential oils and hydroalcoholic extracts of P. eldarica leaf and bark were evaluated on MCF-7 and HeLa cells [24]. Since this plant is native to Iran, it was necessary to conduct this more extensive, confirm the initial data, and identify molecular pathways. Therefore, we sequentially extracted P. eldarica from hydrophobic to hydrophilic compounds by hexane, methanol, and water extracts and essential oils from three parts of P. eldarica. Then, the most common lung cancer cell line, A549 was chosen and cellular and molecular mechanisms of extracts and essential oils were evaluated.
The results showed that the scavenging effects of methanolic extracts were higher than those of the other extracts and essential oils. These results indicate that methanolic extracts possess considerably higher antioxidant activities, and among them, the methanolic bark and needle extracts demonstrated high DPPH radical scavenging activity (93.9% and 83.90%, respectively). The hydroalcoholic extract of P. massoniana bark had the DPPH radical scavenging activity in the study of Yu et al. [25]. According to Dudonne et al., these extracts are potential natural antioxidant sources [26].
Patient-derived cancer cell lines could be used for the study of the effects of drugs on cancer cells. Scientists use different methods, such as monolayer in vitro assays, 3D cultures, bioreactors, or in vivo evaluations with immunocompromised animals [27,28]. In the present study, we used a monolayer A549 cell line to evaluate the biological effects of essential oils and extracts. For this purpose, the cytotoxic effects of essential oils and hexane, methanolic, and aqueous extracts of P. eldarica bark, needles, and pollen were evaluated using the A549 lung cancer cell line. Based on the MTT assay results, the cell mortality rate in all groups increased in a dose-dependent manner, and the methanolic needle extract, with an IC 50 value of 0.3 µg/mL, presented the highest cytotoxicity when compared with other groups. After 48 h of treatment with the pollen essential oil and extracts group, the pollen hexane extract (IC 50 = 31.7 µg/mL) was found to be more cytotoxic than the other employed treatments. Our findings, with those observed in previous reports on the effects of pine essential oils and extracts on cancer cell lines, revealing time-and dosedependent effects. Reportedly, the cytotoxic effect of P. massoniana bark hydroalcoholic extract on a human breast cancer cell line was demonstrated at a higher IC 50 value than 32.5 µg/mL [25]. Compared to our results, the IC 50 values of P. eldarica bark essential oil and hexane, aqueous, and methanolic extracts were 17.9, 18.0, 25.1, and 209 µg/mL, respectively. P. eldarica essential oil and extracts have been shown to prevent the growth of several cancer cell lines, with the essential oil selectively inhibiting the growth of human breast cancer MCF-7 and HeLa cells, presenting IC 50 values of 0.032 and 0.038 µg/mL, respectively [24].
Herein, we observed that the pollen hexane extract, bark essential oil, and methanolic needle extract of P. eldarica significantly suppressed the migration of lung cancer cells at 24 and 48 h. The migration rates of experimental groups were calculated, with the methanolic needle extract revealing the lowest rate (1.87 µm/h in 48 h). Recent studies have shown that pine extracts can inhibit lung cancer cell migration by reducing the expression of matrix metalloproteinases [29,30]. Matrix metalloproteinases are proteolytic enzymes that are secreted from tumor cells, involved in tumor metastasis by degrading the extracellular matrix. The result of the colony formation assay showed that colony numbers in treated groups decreased significantly compared to the control group (p < 0.01). Some studies suggest that natural products can inhibit colony formation by the change in tumor suppressor genes expression [31].
We explored the anticancer mechanisms of P. eldarica extracts and essential oils on A549 lung cancer cells. Apoptosis is a crucial pathway for reducing mutated or growing cancer cells [32]. Reportedly, some natural products have been shown to inhibit the growth of cancer cell lines by inducing apoptosis [33,34]. To evaluate the effects of P. eldarica extracts and essential oils on apoptosis, several experiments, including DAPI staining, Annexin V/PI staining, and cell cycle analysis, was performed. Morphological changes during apoptosis were observed by DAPI staining under a fluorescence microscope. This occurs by nucleic acid denaturation in the cytoplasm. DAPI binds to DNA and forms a complex that emits bluish-white fluorescence by ultraviolet light exciting, so small amounts of DNA can be visualized. Condensed and fragmented chromatin, as specific properties of apoptotic cells, were detected in all extract-and essential oil-treated groups, with a large number of apoptotic cells observed in the methanolic needle extract-treated group containing 0.3 µg/mL of extract; these findings were in agreement with our results obtained following the Annexin V/PI test. In the Annexin V/PI assay, the apoptotic cell population of A549 cells was increased from 0.5% in the control group to 66.6%, 56.5%, and 29.7% in the methanolic needle extract-, bark essential oil-, and pollen hexane extract-treated groups, respectively. The methanolic needle extract-treated cells showed the highest population of apoptotic cells when compared with the other two groups. Thus, we conclude that this group showed suitable efficacy as a potential therapeutic option.
In the cell cycle analysis, A549 cells treated with the methanolic needle extract (0.3 µg/mL) and pollen hexane extract (31.7 µg/mL) revealed cell accumulation in the sub-G1 phase, and the percentage of cells in the S phase was reduced in the methanolic needle extract-treated group. Ma et al. have investigated the cytotoxic effects of P. massoniana extract on HeLa cells, demonstrating that cell cycle arrest occurs in the sub-G1 phase [35]; this was in agreement with our results. Sub-G1 cells were non-apoptotic and viable. Some cells frequently returned and then proceeded to mitosis (hypodiploid and hyperdiploid cells), before arrested in the subsequent cell-cycle interphase. Arresting cells in sub-G1 inhibit cells entering mitosis [36]. Therefore, methanolic needle extract and pollen hexane extract caused post-mitotic arrest in proceeding cell cycles. This process was p53-dependent because when cells divided, p53 arrested the cell cycle in G1 [36]. Furthermore, we observed increased p53 protein expression levels in experimental groups when compared with those in the control group. p53 is a known tumor suppressor gene [37]. In several types of cancers, p53 is damaged owing to mutations, and p53 protein levels are decreased, leading to malignant tumors [38,39].
The cell cycle analysis showed that the accumulation of A549 cells in the S phase was observed after 48 h treatment with the bark essential oil (17.9 µg/mL). Cell cycle arrest at the S phase may result from different mechanisms of action, for example, interfering with cell cycle progression and apoptosis induction, probably by reducing Bcl-2 protein expression. Mitochondria play a critical role in apoptosis, releasing proapoptotic proteins from their membranes into the cytosol. Bcl-2 family proteins are essential mitochondrial regulators of apoptosis [40]. The primary proapoptotic member of these proteins is Bax, which can initiate the release of cytochrome c and subsequent activation of caspase-3 and caspase-9. Bcl-2 acts as an antiapoptotic protein. Bcl-2 inhibits cytochrome c translocation and inhibits caspase-3 activation and apoptosis initiation [19]. High Bax protein expression levels correlate with low Bcl-2 expression levels, which would stimulate apoptosis. The ratio of Bax/Bcl2 expression levels is a diagnostic indicator of cancer [41]. In the present study, Western blot analysis was performed to analyze changes in protein expression of Bcl-2 and Bax in the pollen hexane extract, bark essential oil, methanolic needle extract, and control groups. Based on our Western blot analysis results, the ratio of Bax/Bcl-2 expression levels, as a predictor of apoptosis, increased 6.8-and 27.2-times in the pollen hexane extract-(p < 0.05) and methanolic needle extract-treated groups (p < 0.01) when compared with the control, respectively. In recent studies, pine needle essential oil reportedly downregulated Bcl-2 protein expression and induced apoptosis in oral squamous cell carcinoma [42]. This was in agreement with the results of the present study, which showed that the methanolic needle extract induced apoptosis at a concentration of 0.3 µg/mL by altering Bcl-2 and Bax protein expression in A549 cells more efficiently than other evaluated extracts.
In normal cells, caspases exist in the procaspase form with specific molecular weights. However, upon stimulation of cells with an apoptogenic agent, caspases are activated via proteolytic cleavage. After activation, caspases can cleave other proteins, such as PARP, leading to apoptosis [43]. A previous report has reported the ability of caspases to induce apoptosis in cancer cells following treatment with natural plant-derived extracts [20]. For instance, pine cone water extract reportedly activates caspase-3 and induces apoptosis in a lung cancer cell line [20]. We observed elevated levels of cleaved caspase-3 and cleaved PARP in A549 cells treated with IC 50 concentrations of the pollen hexane extract, bark essential oil, and methanolic needle extract. These results suggest that activation of the caspase pathway is essential for inducing apoptosis in A549 cells following treatment with P. eldarica extracts.
Our results showed that P. eldarica extracts mainly, methanolic needle extract, can prevent and treat cancer by various mechanisms. Cell cycle arrest, inhibition of cancer cell migration, and proliferation were some of the anticancer effects. Induction of apoptosis by downregulating of antiapoptotic gene products, such as Bcl-2, and stimulation of apoptotic protein expression such as Bax, caspases, and p53 lead to apoptosis, which was confirmed with the Annexin V/PI test.

Plant Material
P. eldarica needles, pollen, and bark were collected in May and June 2018 from Tabriz district in Iran and identified (No.4036) by the Herbarium of Pharmacy Faculty, Tabriz University of Medical Sciences. All air-dried specimens were sequentially extracted in the order of hexane, methanol, and water by using the Soxhlet apparatus. Then, the extracts were concentrated using a rotary evaporator (Heidolph, Germany), dried using a freeze dryer (Christ Alpha1-4, Germany), and stored at −20 • C until use. For essential oil preparation, 100 g of P. eldarica needles, bark, and pollen powders were hydro-distilled in a Clevenger device separately. The resulting essential oil was dried using anhydrous sodium sulfate and stored in a sealed vial at 4 • C until further analysis. Different concentrations of the extracts and essential oil (0.001, 0.01, 0.1, 1, 10, and 100 µg/mL) were prepared. DMSO was used as the solvent to prepare the stock solution (max 1% v/v).

DPPH Radical Scavenging Assay
The DPPH radical scavenging activity was determined by adding 1 mL of the DPPH ethanolic solution (0.1 mM) to 1 mL of each extract (1 mg/mL), shaken, and incubated at 37 • C for 30 min in the dark. Then, the absorbance of each extract was measured at 517 nm by ultraviolet-visible (UV-Vis) spectrophotometry (Ultrospec 2000, Pharmacia Biotech, Amersham, UK) [44]. DPPH ethanolic solution was used as a control, and 70% ethanol was used as a blank. The radical scavenging activity as a percentage of DPPH discoloration was calculated using the following equation: DPPH activity (%) = UV Absorbance of extract UV Absorbance of control × 100.

Cell Culture
The A549 cell line was obtained from the Pasteur Institute (Tehran, Iran). Cells were grown in RPMI-1640 supplemented with 10% FBS and 1% penicillin/streptomycin. The media was sterilized using 0.22-µm filters and stored at 4 • C before use. Cells were grown in RPMI-1640 medium at 37 • C in a 5% CO 2 incubator.

MTT-Based Cytotoxicity Assay
The cytotoxicity of the essential oil and different extracts of P. eldarica against the A549 cell line was determined by employing a rapid colorimetric assay using MTT [45][46][47][48][49]. In this assay, soluble MTT is metabolized by the mitochondrial enzyme activity of viable cells into an insoluble formazan dye product. Subsequently, formazan was dissolved in DMSO and measured spectrophotometrically at 540 nm. Briefly, 180 µL of a cell suspension (5 × 10 4 cells/mL of medium) was seeded in 96-well microplates and incubated for 24 h (37 • C, 5% CO 2 humidified air). Then, prepared concentrations of essential oil and extracts ranging from 0 to 100 µg/mL were added to each well. The plates were incubated for another 48 h under the same conditions. DOX and 1% DMSO were used as positive and negative controls, respectively. To assess cell survival, 20 µL of MTT solution (5 mg/mL in PBS) was added to each well, and the plates were incubated at 37 • C for 3 h. Then, the old media containing MTT was removed and 150 µL DMSO was added to each well and pipetted to dissolve the formed formazan crystals. Absorbance was measured using an ELISA plate reader (Awareness Technology, Inc., Palm City, FL, USA) at 540 nm. Each extract concentration was repeated three times, and cell survival in the negative control was assumed to be 100%. The IC 50 value was calculated using GraphPad Prism software version 6.0.

DAPI Staining
DAPI staining was performed to evaluate nuclei fragmentation. A549 cells (2 × 10 6 cells/well) were seeded into 6-well plates. After incubation, cultured lung cancer cells at 70% confluence were treated with IC 50 concentrations of extracts and essential oil for 48 h. Next, the cells were washed with PBS, fixed with paraformaldehyde (4%) for 20 min and permeabilized with 0.15% Triton X-100 for 15 min. Subsequently, the cells were stained in the dark with 0.1% DAPI for 15 min, washed with PBS, and observed by fluorescence microscopy (Olympus, Seoul, Korea) [50].

Annexin V/PI Flow Cytometry
A549 cells (2 × 10 6 cells/well) were seeded in 6-well plates and incubated in RPMI-1640 medium for 48 h. Subsequently, the IC 50 concentration of the extracts and essential oil was added to the wells and incubated for 48 h. Then, the cells were trypsinized and washed with cold PBS, supernatants were removed, and 100 µL binding buffer was added at icecold temperature. Finally, lung cancer cells were stained with Pl and Annexin V-fluorescein isothiocyanate (FITC) according to the manufacturer's instructions. A FACSCalibur flow cytometer (BD Biosciences, San Jose, CA, USA) was used to analyze apoptosis rate. The data were analyzed using the FlowJo software package (Treestar, Inc., San Carlos, CA, USA) [51].

Scratch Migration Assay
To examine the migration of A549 cells in two-dimensional monolayers, cells were seeded in 6-well plates and incubated in RPMI-1640 medium. After the cells reach to confluent monolayer (at least 70% confluency), a thin linear scratch "wound" was created using a 200 µL pipette tip. The monolayers were then washed with PBS to remove detached cells and treated with IC 50 concentrations of extracts and essential oil; images were captured at 0, 24, and 48 h after treatment with a digital camera connected to a bright-field microscope. Wound width was calculated as the distance between two edges of the scratch using CorelDRAW software [52]. The migration rate was calculated using the following equation: Migration rate (µm/h) = initial wound width (µm) − final wound width(µm) duration of migration (h)

Cell Cycle Analysis
A549 cells (3 × 10 5 cells/mL) were seeded in 6-well plates. After 24 h, cells were treated with IC 50 concentrations of extracts and essential oil. Cells were trypsinized after 48 h and fixed in ice-cold 70% ethanol at 4 • C for 2 h. A549 cells were centrifuged and washed with PBS, treated with 1 mL PI master mix containing 950 µL PBS, 40 µL PI, and 10 µL RNase A to stain cellular DNA, incubated at room temperature for 30 min, and analyzed using a FACSCalibur flow cytometer (BD Biosciences, San Jose, CA, USA) [50].

Colony Formation Assay
The colony formation assay was performed as previously described [54]. A549 cells were seeded (1000 cells/well) in 6-well plates in triplicate in RPMI-1640 medium. After 24 h, the medium was removed and replaced with the IC 50 concentration of respective extracts and essential oil. The group with no added treatment was considered as the control. After 14 days, the colonies were stained with 0.5% crystal violet solution containing 25% methanol for 15 min at 37 • C and washed with excess dye. The colony numbers were counted using the ImageJ software.

Statistical Analysis
The data were processed using GraphPad Prism version 6.0 (GraphPad Software, Inc., La Jolla, CA, USA) for averages and standard deviation (SD) by performing one-way and two-way analyses of variance. Data are presented as mean values ± SD, and p-values less than 0.05 were considered significant.

Conclusions
In the present study, extracts and essential oils of various parts of P. eldarica revealed robust antioxidant activity. Especially, the pollen hexane extract, bark essential oil, and methanolic needle extract of P. eldarica inhibited the proliferation and cell cycle progression in A549 lung cancer cells, mainly through Bcl-2 downregulation, stimulation of p53 expression, and induction of apoptosis in a caspase-dependent manner. Furthermore, the pollen hexane extract, bark essential oil, and methanolic needle extract suppressed the migration, and colony formation of lung cancer cells. Therefore, our findings suggest that P. eldarica extracts and essential oils, especially the methanolic extract of P. eldarica needles (the most active extract), can be potentially developed as chemopreventive natural resources for cancer prevention and treatment, with further investigations warranted to determine the main components responsible for anticancer activities.