Antioxidant and Neuroprotective Activity of Extra Virgin Olive Oil Extracts Obtained from Quercetano Cultivar Trees Grown in Different Areas of the Tuscany Region (Italy)

Neurodegenerative diseases are driven by several mechanisms such as inflammation, abnormal protein aggregation, excitotoxicity, mitochondrial dysfunction and oxidative stress. So far, no therapeutic strategies are available for neurodegenerative diseases and in recent years the research is focusing on bioactive molecules present in food. In particular, extra-virgin olive oil (EVOO) phenols have been associated to neuroprotection. In this study, we investigated the potential antioxidant and neuroprotective activity of two different EVOO extracts obtained from Quercetano cultivar trees grown in two different areas (plain and hill) of the Tuscany region (Italy). The different geographical origin of the orchards influenced phenol composition. Plain extract presented a higher content of phenyl ethyl alcohols, cinnammic acids, oleacein, oleocanthal and flavones; meanwhile, hill extract was richer in lignans. Hill extract was more effective in protecting differentiated SH-SY5Y cells from peroxide stress thanks to a marked upregulation of the antioxidant enzymes heme oxygenase 1, NADPH quinone oxidoreductase 1, thioredoxin Reductase 1 and glutathione reductase. Proteomic analysis revealed that hill extract plays a role in the regulation of proteins involved in neuronal plasticity and activation of neurotrophic factors such as BDNF. In conclusion, these data demonstrate that EVOOs can have important neuroprotective activities, but these effects are strictly related to their specific phenol composition.


Introduction
In the "era" of the economic prosperity, the non-communicable diseases such as cardiovascular diseases, neurodegenerative diseases, cancer, obesity and diabetes, are collectively accountable for almost 70% of all deaths worldwide [1]. The complexity of these diseases makes it difficult to fight them with single-target molecules. Indeed, the research is shifting its strengths on multitarget compounds [2,3] and in recent years it is focusing on bioactive molecules naturally present in food as potential "weapons" against chronic degenerative diseases. On the base of these premises, nutrition could represent a key factor of lifestyle for counteracting the development and progression of such pathologies. In particular, the Mediterranean diet, rich in vegetables, cereals, fruits, legumes and monounsaturated fatty acids mainly derived from olive oil, has been associated with a reduced incidence of neurodegenerative diseases and enhanced cognitive performance [4,5].
Neurodegenerative diseases represent an increasingly public health problem, especially in the aging population. These pathologies are multifactorial non communicable diseases driven by several but linked mechanisms such as inflammation, abnormal protein aggregation, excitotoxicity, oxidative stress and mitochondrial dysfunction [6][7][8][9][10]. In addition, the accumulation of proteins and lipids modified by the excessive production of reactive oxygen species (ROS) leads to a self-feeding loop [11]. Oxidative stress consists in an imbalance condition where the production of reactive species (mainly ROS) exceeds their detoxification, leading to an abnormal accumulation of radical species and oxidative damages. Furthermore, oxidative stress triggers and is triggered by mitochondrial dysfunction [12]. Stated the key role played by mitochondria in energy metabolism and modulation of redox homeostasis, their dysfunction might contribute to the pathogenesis of neurodegeneration. Due to its high oxygen requirements, relative poorness of antioxidant defenses, elevated levels of unsaturated fatty acids and abundance in iron content, the brain is particularly sensible to oxidative stress [13,14]. Consequently, neurons are prone to oxidative injury and the physiological process of aging together with the aging-associated diseases share the promotion of a redox imbalance [15]. Of note, neurons are cells characterized by limited self-renewal, implying that damages can accumulate over time [13]. Regardless of how oxidative stress is caused, when encountered, cells counteract the damaging effect of this condition by activating the endogenous antioxidant defense system, which is unfortunately compromised in the context of neurodegeneration. Extensive studies during the past decade have proven the notion that the transcription factor NF-E2-related factor 2 (Nrf2) is an essential element for up-regulating the antioxidant defense system [16][17][18].
In addition, neuronal growth and survival are ensured by the neurotrophic signaling pathway and alteration of specific neurotrophic factors leads to brain degeneration [19]. In particular, the brain-derived neurotrophic factor (BDNF) is essential for the survival and normal functioning of mature neurons and its level are markedly depressed in neurodegenerative diseases [20][21][22]. These relationships suggest the BDNF and the Nrf2 signaling pathway are potential targets for supporting neuronal survival and regeneration of impaired neuronal structures and synaptic connectivity.
Unfortunately, to date, no therapeutic strategies are available for neurodegenerative diseases.
It has been suggested that the positive role of the Mediterranean diet is probably related to the presence of a huge quantity and different phenols [23,24]. One of the main source of phenols from the Mediterranean diet is the consumption of extra-virgin olive oil (EVOO) in amounts going from 30 to 50 g/day [25], representing the principle fat consumed in the Mediterranean diet [26]. In vitro and in vivo studies demonstrated preventive and protective properties of EVOO against neurodegeneration. The major polar phenolic compounds present in EVOO such as oleacein (dialdehydic form of decarboxymethyl elenolic acid linked to hydroxytyrosol, 3,4-DHPEA-EDA), oleuropein aglycone isomer (3,4-DHPEA-EA), oleocanthal (the dialdehydic form of decarbox-ymethylelenolic acid linked to tyrosol, p-HPEA-EDA), ligstroside aglycon (p-HPEA-EA), hydroxytyrosol (3,, tyrosol (p-HPEA), have been demonstrated to possess anti-inflammatory [27], anti-microbial and anti-viral activities [28], to counteract oxidative stress and modulate survival signaling pathways [29]. Many factors can influence EVOO phenolic composition Antioxidants 2021, 10, 421 3 of 24 and concentration including cultivars, climate conditions, ripening stage of the olives but also olive oil production processes and storage [30,31]. Therefore, the definition of an exact phenolic content is hard and referring to different studies the amount could range from 200 to 1000 mg/kg [28,30,[32][33][34][35]. Furthermore, the specific phenolic pattern of different EVOOs could influence their biological activity such as neuroprotection.
In this study, we investigated the potential antioxidant and neuroprotective effects of two different EVOO extracts obtained from Quercetano cultivar trees grown in two different areas of the Tuscany region (Italy). In particular, we characterized the extracts by HPLC-DAD/MS and evaluated their neuroprotective activity against H 2 O 2 -induced oxidative stress in differentiated SH-SY5Y cells. Moreover, proteomics combined with pathways analyses, was used to investigate the molecular events related to the protective effects and to gain insight into the underlying mechanisms of neuroprotection.

EVOO Sample Extraction
EVOO samples were obtained from the fruits of the olive cultivar Quercetano (Olea europaea L.) which were picked at the same stage of maturity during the same month season 2019 in two olive groves located at different altitude (hill, plain) and with some different soil type (clay, silt, sandy) in Tuscany region (Italy). The phytoextracts of hill and plain EVOOs were prepared essentially as previously described [36,37]. Briefly, EVOO (3 g) was mixed with n-hexane (12 mL) and acetonitrile (15 mL). The resulting mixture was homogenized using a vortex mixer and centrifuged at 4000 rpm for 5 min, at 25 • C. Then, the acetonitrile phase was collected and evaporated under reduced pressure, to obtain the phytoextracts.

Cell Culture and Treatments
The SH-SY5Y human neuroblastoma cell line was obtained from Sigma-Aldrich (Milan, Italy). Cells were grown in DMEM supplemented with 10% (v/v) FBS, 2 mM glutamine, 50 U/mL of penicillin and 50 µg/mL of streptomycin and maintained at 37 • C in a humidified incubator with 5% CO 2 as reported in [40]. Cell differentiation was induced reducing serum levels of the medium to 1% with 10 µM of RA for seven days prior to treatments. Differentiated SH-SY5Y were treated with plain and hill extracts for 24 h. Extracts were dissolved in DMSO. The control group was treated with an equivalent volume of the vehicle alone. Oxidative stress was induced exposing cells to 700 µM H 2 O 2 for different times depending on the assay.

MTT and Lactate Dehydrogenase Activity Assays
Cells were treated with different concentrations of plain and hill extracts (1 µg/mL-500 µg/mL) for 24 h. The induction of oxidative stress was achieved with 700 µM of H 2 O 2 for 1 h. Cell viability was evaluated by measuring formazan formation as previously reported [41]. Cells were incubated with 0.5 mg/mL of MTT solution for 1 h at 37 • C. After incubation, MTT was removed and 100 µL of DMSO were added, the absorbance was recorded at λ = 595 nm using a microplate spectrophotometer (VICTOR3 V Multilabel Counter; PerkinElmer, Wellesley, MA, USA). Lactate dehydrogenase (LDH) activity was evaluated in the culture medium and the test was performed using the Lactate Dehydrogenase Activity Assay Kit according to the manufacturer's instructions.

Intracellular ROS Production Assay
The fluorescent probe DCFH-DA was used to monitor the production of intracellular reactive oxygen species [42]. Differentiated SH-SY5Y were treated with different concentrations of plain and hill extracts for 24 h and then incubated with 10 µg/mL DCFH-DA in DMEM w/o FBS for 30 min. After DCFH-DA removal, cells were incubated with 200 µM H 2 O 2 for 15 min. Cell fluorescence was measured using 485 nm excitation and 535 nm emission with a microplate spectrofluorometer (VICTOR3 V Multilabel Counter).

Intracellular GSH Levels Assay
The levels of reduced glutathione (GSH) were evaluated by a fluorometric assay, using the fluorescent probe MCB [43]. Briefly, after treatments, differentiated SH-SY5Y were incubated with 50 µM MCB in serum-free medium for 30 min at 37 • C. After incubation, fluorescence was measured at 355 nm (excitation) and 460 nm (emission) with a microplate spectrofluorometer (VICTOR3 V Multilabel Counter; PerkinElmer).

RNA Extraction
Total RNA was extracted using RNeasy Mini Kit (QIAGEN GmbH, Hilden, Germany), following the manufacturer's protocol. The yield and purity of the RNA were measured using NanoVue Spectrophotometer (GE Healthcare, Milano, Italy).

Analysis of mRNA Levels by Reverse Transcriptase Polymerase Chain Reaction
cDNA was obtained by reverse transcribing mRNA starting from 1 µg of total RNA using iScript cDNA Synthesis Kit (BIO-RAD, Hercules, CA, USA), following the manufacturer's protocol. The subsequent polymerase chain reaction (PCR) was performed in a total volume of 10 µL containing 2.5 µL (12.5 ng) of cDNA, 5 µL SsoAdvanced Universal SYBR Green Supermix (BIO-RAD) and 0.5 µL (500 nM) of each primer (Sigma-Aldrich). The primers used are reported in Table 1, 18S rRNA was used as reference gene.

Proteomic Analysis
For proteomic analysis, differentiated SH-SY5Y cells were treated with plain and hill extracts (10 µg/mL) for 24 h as described above. At the end of treatments, cells were collected and washed with PBS. After centrifugation (1000× g for 5 min), the resulting pellets were immediately frozen and stored at −80 • C until use. For proteomic studies, each condition was performed in triplicate.
Cell pellets were resuspended in rehydration solution [44] and protein contents of resulting protein extracts were measured with an RC-DC Protein Assay from Bio-Rad.
The 2DE was carried out as previously described [45]. Briefly, 200 µg of proteins were filled up to 450 µL in rehydration solution. Immobiline Dry-Strips (GE Health Care Europe; Uppsala, Sweden); 18 cm, nonlinear gradient (pH 3-10) were rehydrated overnight in the sample and then transferred to the Ettan IPGphor Cup Loading Manifold (GE Healthcare) for isoelectrofocusing (IEF). The second dimension (Sodium Dodecyl Sulphate-Polyacrylamide Gel Electrophoresis; SDS-PAGE) was carried out by transferring the proteins to 12% polyacrylamide, running at 16 mA per gel and 10 • C for about 16 h, using the Protean ® Plus Dodeca Cell (Bio-Rad). The gels were stained with Ruthenium II tris (bathophenanthroline disulfonate) tetrasodium salt (Cyanagen Srl, Bologna, Italy) (RuBP). ImageQuant LAS4010 (GE Health Care) was used for the acquisition of images. The analysis of images was performed using Same Spot (v4.1, TotalLab; Newcastle Upon Tyne, UK) software. The spot volume ratios among the three different conditions (Control, hill, plain) were calculated using the average spot normalized volume of the three biological replicates performed in duplicate. The software included statistical analysis calculations.

Spot Digestion and Protein Identification
The gel pieces were trypsin digested and analyzed by LC-MS as previously described [46]. Each digested spot sample was analyzed by LC-MS/MS using a Proxeon EASY-nLCII (Thermo Fisher Scientific, Milan, Italy) chromatographic system coupled to a Maxis HD UHR-TOF (Bruker Daltonics GmbH, Bremen, Germany) mass spectrometer. The raw data were processed using PEAKS Studio v7.5 software (Bioinformatic Solutions Inc, Waterloo, ON, Canada) using the function "correct precursor only". The mass lists were searched against the nextprot database including isoforms (version as of June 2017; 42,151 entries) using 10 ppm and 0.05 Da as the highest error tolerances for parent and fragment ions, respectively. Carbamidomethylation of cysteines was selected as fixed modification and oxidation of methionines and deamidation of asparagine and glutamine as variable modifications allowing two missed cleavages.

Bioinformatic Analysis
For gene ontology analysis, including differential molecular function and biological processes, PANTHER software (Protein Analysis Through Evolutionary Relationships; http://www.pantherdb.org/genes/batchIdSearch.jsp (accessed 3 January 2021), was used to classify genes into distinct categories of molecular functions and biological processes. Proteins differentially expressed obtained both from hill vs Control and plain vs Control comparisons were functionally analyzed using QIAGEN's Ingenuity Pathway Analysis (IPA, QIAGEN Redwood City, USA, www.qiagen.com/ingenuity, Build version: 321501M Content version: 21249400, accessed 15 January 2021) with the aim to determine the predominant canonical pathways and interaction network involved. Swiss-Prot accession numbers and official gene symbols were inserted into the software along with corresponding comparison ratios and p values. Based on known associations in the literature, canonical pathways associated with differentially expressed proteins were defined. A comparison of the different analyses was created and the upstream regulators whose activity appears to change in a significant manner according to the activation z-score value were shown. Finally, the impact of activated or inhibited regulators on downstream functions and diseases were investigated. The use of an algorithm allowed to merge upstream and downstream results from the upstream regulator, through one or more iterations. The networks were merged only if the overlap of protein targets was possible and of statistical significance (Fisher's Exact Test). Higher scoring hypotheses were those with more consistent causal paths represented by a high Consistency Score.

Western Blot Analysis
Western blot (WB) was performed in order to validate peculiar increase of expression of both DPYSL2 and BDNF found in hill with respect to plain and Control samples with 2-DE and RT-PCR, respectively. Aliquots of protein samples (5 µg for DPYSL2 and 50 µg for BDNF) were mixed with Laemmli solution, run in 4-15% polyacrylamide gels (Mini-PROTEAN ® Precast Gels, Biorad, Hercules, CA, USA) using a mini-Protean Tetracell (Biorad, Hercules, CA, USA) and transferred onto nitrocellulose membranes (0.2 µm) using a Trans-Blot Turbo transfer system (Biorad) as previously described [45]. Anti-DPYSL2 (Cell Signaling Technology, Beverly, MA, USA) and anti-BDNF (Genetex, Irvine, CA, USA) antibodies were used at 1:1000 dilution. Moreover, an anti-β-actin (Merck KGaA, Darmstadt, Germany), was used for internal normalization. HRP-goat anti-rabbit secondary antibody was used at 1:10,000 dilution. Immunoblots were developed using the enhanced chemiluminescence detection system (ECL). The chemiluminescent images were acquired using LAS4010 (GE Health Care Europe, Upsala, Sweden). The immunoreactive specific bands were quantified using Image Quant-L software.

Statistical Analysis.
All analyses were performed at least in triplicate and values were expressed as mean ± standard error. In experiments with SH-SY5Y cell cultures, one-way ANOVA was used to compare differences among groups followed by Dunnett's or Bonferroni's test (Prism 5; GraphPad Software, San Diego, CA, USA). Differences at the level p < 0.05 were considered statistically significant. In 2DE experiments, a comparison among the different treatments was performed. The significance of the differences of normalized volume for each spot was calculated by the software Same Spot including the analysis of variance (ANOVA test). The protein spots significantly differentially expressed were cut out from the gel and identified by LC-MS analysis.

Phenolic Content
The soil characteristic of the hill and plain orchards are reported in Table 2. Interestingly, the two soils were very different in terms of silt, sandy, pH and mineral composition.
The analysis of phenolic components of EVOOs extracts from the hill and plain orchards shows a similar qualitative composition, but a different quantitative composition. In particular, the total phenol content is higher in the plain than in the hill extract (Table 3). Moreover, phenyl ethyl alcohols (tyrosol, hydroxytyrosol), oleacein (3,4-DHPEA-EDA), oleocanthal (p-HPEA-EDA), cinnammic acids (p-coumaric acid, ferulic acid), flavones (luteolin, apigenin) are higher in the plain extract, vice versa lignans (pinoresinol, ace-  a Asterisks after the phenolic substance's name indicate the level of significance (One way ANOVA, Tukey's test for pairwise comparison) in the difference between Hill and Plain extracts: * p < 0.05; ** p < 0.01; *** p < 0.001; Analyses were carried out by the "Regional laboratory for soil analysis and plant production" http://www.agriligurianet.it (accessed 21 January 2021).

EVOO Extracts-Mediated Neuroprotection against Oxidative Damage
To assess the potential cytotoxicity of the EVOO extracts, differentiated SH-SY5Y cells were exposed to increasing concentrations of the extracts for 24 h (Figure 1). Treatment with the hill extract led to a significant decrease of cell viability in respect to control cells at concentrations higher than 50 µg/mL, meanwhile treatment with the plain extract significantly reduced cell viability at concentrations higher than 100 µg/mL. Interestingly, both extracts significantly increased cell viability at the lowest concentrations. Figure 1. Viability of differentiated SH-SY5Y treated with EVOO extracts. Cells were treated with increasing concentrations (1 to 500 μg/mL) of EVOO extracts for 24 h and the cellular viability was measured by MTT assay as described in Materials and Methods. Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Dunnett's test. * p < 0.05 with respect to CTRL.
To assess the neuroprotective activity of the extracts, SH-SY5Y cell were pretreated with 10 μg/mL of the extracts and after 24 h were exposed to H2O2 to induce oxidative stress ( Figure 2). This extract concentration was selected as it is one order of magnitude lower than the toxic concentration and can be considered safe. Peroxide induced a strong and significant reduction of cell viability in respect to control cells, meanwhile both hill and plain extracts were able to significantly increase cell viability in respect to cells exposed to H2O2 (Figure 2A). Of note, the hill extract significantly increased cell viability in respect to the plain extract. Concerning LDH release, a marker of no-specific cell damage, H2O2 triggered a strong and significant increase of LDH release in respect to control cells ( Figure 2B). Only the hill extract was able to significantly reduce the activity of LDH in the culture medium in respect to H2O2 exposed cells, confirming the higher antioxidant activity of the hill in respect to the plain extract. Neuroprotective activity of EVOO extracts against H2O2 in SH-SY5Y cells. Cells were treated with 10 μg/mL of the EVOO extracts for 24 h and then exposed to H2O2. (A) Cell viability was measured by MTT assay and (B) LDH was measured as LDH activity in in the culture medium as described in Materials and Methods. Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Bonferroni's test. * p < 0.05 with respect to CTRL; ° p < 0.05 with respect to H2O2; § p < 0.05 with respect to plain. Figure 1. Viability of differentiated SH-SY5Y treated with EVOO extracts. Cells were treated with increasing concentrations (1 to 500 µg/mL) of EVOO extracts for 24 h and the cellular viability was measured by MTT assay as described in Materials and Methods. Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Dunnett's test. * p < 0.05 with respect to CTRL.

Antioxidant Properties of Hill and Plain Extracts
To assess the neuroprotective activity of the extracts, SH-SY5Y cell were pretreated with 10 µg/mL of the extracts and after 24 h were exposed to H 2 O 2 to induce oxidative stress ( Figure 2). This extract concentration was selected as it is one order of magnitude lower than the toxic concentration and can be considered safe. Peroxide induced a strong and significant reduction of cell viability in respect to control cells, meanwhile both hill and plain extracts were able to significantly increase cell viability in respect to cells exposed to H 2 O 2 ( Figure 2A). Of note, the hill extract significantly increased cell viability in respect to the plain extract. Concerning LDH release, a marker of no-specific cell damage, H 2 O 2 triggered a strong and significant increase of LDH release in respect to control cells ( Figure 2B). Only the hill extract was able to significantly reduce the activity of LDH in the culture medium in respect to H 2 O 2 exposed cells, confirming the higher antioxidant activity of the hill in respect to the plain extract. To assess the neuroprotective activity of the extracts, SH-SY5Y cell were pretreated with 10 μg/mL of the extracts and after 24 h were exposed to H2O2 to induce oxidative stress ( Figure 2). This extract concentration was selected as it is one order of magnitude lower than the toxic concentration and can be considered safe. Peroxide induced a strong and significant reduction of cell viability in respect to control cells, meanwhile both hill and plain extracts were able to significantly increase cell viability in respect to cells exposed to H2O2 (Figure 2A). Of note, the hill extract significantly increased cell viability in respect to the plain extract. Concerning LDH release, a marker of no-specific cell damage, H2O2 triggered a strong and significant increase of LDH release in respect to control cells ( Figure 2B). Only the hill extract was able to significantly reduce the activity of LDH in the culture medium in respect to H2O2 exposed cells, confirming the higher antioxidant activity of the hill in respect to the plain extract. Neuroprotective activity of EVOO extracts against H2O2 in SH-SY5Y cells. Cells were treated with 10 μg/mL of the EVOO extracts for 24 h and then exposed to H2O2. (A) Cell viability was measured by MTT assay and (B) LDH was measured as LDH activity in in the culture medium as described in Materials and Methods. Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Bonferroni's test. * p < 0.05 with respect to CTRL; ° p < 0.05 with respect to H2O2; § p < 0.05 with respect to plain. Neuroprotective activity of EVOO extracts against H 2 O 2 in SH-SY5Y cells. Cells were treated with 10 µg/mL of the EVOO extracts for 24 h and then exposed to H 2 O 2 . (A) Cell viability was measured by MTT assay and (B) LDH was measured as LDH activity in in the culture medium as described in Materials and Methods. Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Bonferroni's test. * p < 0.05 with respect to CTRL; • p < 0.05 with respect to H 2 O 2 ; § p < 0.05 with respect to plain.

Antioxidant Properties of Hill and Plain Extracts
The potential antioxidant activity of hill and plain extracts was investigated through the assessment of intracellular ROS and reduced glutathione (GSH) levels, the main endogenous intracellular antioxidant [47] (Figure 3).
The potential antioxidant activity of hill and plain extracts was investigated through the assessment of intracellular ROS and reduced glutathione (GSH) levels, the main endogenous intracellular antioxidant [47] (Figure 3). . Each bar represents means ± SEM of at least three independent experiments. Data were analyzed by one-way ANOVA followed by Dunnett's test or Bonferroni's test. * p < 0.05 with respect to CTRL; ° p < 0.05 with respect to H2O2.
The ability of plain and hill extracts to counteract H2O2-induced ROS production was evaluated by the DCFH-DA assay. Differentiated SH-SY5Y cells were treated with 10 μg/mL of hill and plain extracts for 24 h and then exposed to H2O2. As reported in Figure  3A, both extracts significantly reduced ROS levels in respect to H2O2 treated cells. Then, SH-SY5Y cells were treated with EVOO extracts for 24 h and GSH basal levels were measured by the MCB assay ( Figure 3B). In agreement with the previous data, both extracts were able to significantly increase GSH levels compared to control cells.

Modulation of Antioxidant and Pro-Survival Genes by Plain and Hill Extracts
Since the hill extract demonstrated a higher neuroprotective activity against oxidative stress compared to the plain extract, we next drilled down the mechanism underlying such protection.
To this purpose, we measured the expression of Nrf2-driven antioxidant genes, heme oxygenase 1 (HMOX1), thioredoxin reductase 1 (TXNRD1), NADPH quinone oxidoreductase 1 (NQO1) and glutathione reductase (GSR) in SH-SY5Y cells after treatment with the extracts (Figure 4). The plain extract was able to significantly increase the expression of HMOX1 and NQO1 genes in respect to control cells. On the other hand, the hill extract not only significantly upregulated all four of the genes in respect to control, but it also significantly increased their expression in respect to the plain extract, suggesting a higher ability of the hill extract to enhance the antioxidant defense system.
Since neurotrophins, such as BDNF, play a fundamental role in neuronal survival [19], we investigated the effect of the treatment with the extracts on BDNF gene expression ( Figure 5). Once again, the hill extract demonstrated a higher bioactivity significantly upregulating BDNF gene in respect both to control and plain extract treated cells. The plain extract did not influence the expression of this neurotrophin. The ability of plain and hill extracts to counteract H 2 O 2 -induced ROS production was evaluated by the DCFH-DA assay. Differentiated SH-SY5Y cells were treated with 10 µg/mL of hill and plain extracts for 24 h and then exposed to H 2 O 2 . As reported in Figure 3A, both extracts significantly reduced ROS levels in respect to H 2 O 2 treated cells. Then, SH-SY5Y cells were treated with EVOO extracts for 24 h and GSH basal levels were measured by the MCB assay ( Figure 3B). In agreement with the previous data, both extracts were able to significantly increase GSH levels compared to control cells.

Modulation of Antioxidant and Pro-Survival Genes by Plain and Hill Extracts
Since the hill extract demonstrated a higher neuroprotective activity against oxidative stress compared to the plain extract, we next drilled down the mechanism underlying such protection.
To this purpose, we measured the expression of Nrf2-driven antioxidant genes, heme oxygenase 1 (HMOX1), thioredoxin reductase 1 (TXNRD1), NADPH quinone oxidoreductase 1 (NQO1) and glutathione reductase (GSR) in SH-SY5Y cells after treatment with the extracts (Figure 4). The plain extract was able to significantly increase the expression of HMOX1 and NQO1 genes in respect to control cells. On the other hand, the hill extract not only significantly upregulated all four of the genes in respect to control, but it also significantly increased their expression in respect to the plain extract, suggesting a higher ability of the hill extract to enhance the antioxidant defense system.
Since neurotrophins, such as BDNF, play a fundamental role in neuronal survival [19], we investigated the effect of the treatment with the extracts on BDNF gene expression ( Figure 5). Once again, the hill extract demonstrated a higher bioactivity significantly upregulating BDNF gene in respect both to control and plain extract treated cells. The plain extract did not influence the expression of this neurotrophin.   Figure 5. Expression of BDNF in SH-SY5Y cells treated with the extracts. Cells were treated with hill and plain extracts (10 µg/mL) for 24 h. Real time-PCR was performed to detect BDNF mRNA levels. Data are expressed as relative abundance compared to untreated cells. Each bar represents mean ± SEM of three independent experiments. Data were analyzed with one-way ANOVA followed by the Bonferroni's test. * p < 0.05 vs. CTRL; § p < 0.05 vs. plain. Figure 6A illustrates a representative two-dimensional electrophoresis (2DE) image of differentiated SH-SY5Y cell proteins. An average of 1220 ± 120 spots was found within a nonlinear pH range from 3 to 10. Comparative analysis was performed between hill or plain extract treated and control cell samples. Venn diagram ( Figure 6B) shows the number of common and exclusive proteins between the two comparisons. One-hundred fifteen protein spots were found differentially expressed in hill extract treated cells with comparison to control cells, of which 26 spots were also differentially expressed in plain extract treated cells compared to control cells. Spots, which showed a significant fold change of expression ≥ 1.2 were subsequently subjected to nano-LC-ESI-MS/MS analysis and identified. Table 4 reports the list of identified proteins, which are exclusive of hill extract treated cells together with their MW, pI, peptides and coverage values of MS/MS, ratio and p values. A list of common proteins both to hill and plain extract treated cells is shown in Table 5. More than one identification was reported for some spots when MW and pI were not distinguishable. Figure 6A illustrates a representative two-dimensional electrophoresis (2DE) image of differentiated SH-SY5Y cell proteins. An average of 1220 ± 120 spots was found within a nonlinear pH range from 3 to 10. Comparative analysis was performed between hill or plain extract treated and control cell samples. Venn diagram ( Figure 6B) shows the number of common and exclusive proteins between the two comparisons. One-hundred fifteen protein spots were found differentially expressed in hill extract treated cells with comparison to control cells, of which 26 spots were also differentially expressed in plain extract treated cells compared to control cells. Spots, which showed a significant fold change of expression ≥ 1.2 were subsequently subjected to nano-LC-ESI-MS/MS analysis and identified. Table 4 reports the list of identified proteins, which are exclusive of hill extract treated cells together with their MW, pI, peptides and coverage values of MS/MS, ratio and p values. A list of common proteins both to hill and plain extract treated cells is shown in Table 5. More than one identification was reported for some spots when MW and pI were not distinguishable.

Bioinformatic Analysis
Using both PANTHER gene classification and Ingenuity Pathway Analysis (IPA) software, we found that the main biological processes, in which the deregulated proteins of hill and plain extract treated SH-SY5Y cells are involved, included "biogenesis," "metabolic process", and "cellular process." All proteins found differentially expressed in hill extract treated cells were analyzed by IPA to discover the most enriched canonical pathways, possible upstream regulators and downstream effects. The software generated two main networks, "Gene expression, RNA damage and repair, RNA post-transcriptional modifications" and "Cellular development, cellular growth and proliferation, nervous system development and function" with 46 and 32 score values, respectively.
The downstream analysis unveiled key biological and cellular functions such as organismal survival, nervous system development, cell death and survival, gene expression and RNA damage and repair. Moreover, all proteins found differentially expressed concurred in an upstream analysis to predict activation or inhibition of potential transcription factors or molecules. Table 6 shows a list of top upstream regulators based on p-value. Moreover, an algorithm connection among upstream regulators, data set molecules and downstream functions or diseases, generated four main regulators with high consistency scores. The following regulators, fibroblast growth factor 2 (FGF2), interleukin-5 (IL5) and transcription factor Sp1 (SP1) revealed that the functional categories most impacted by the presence of the hill EVOO extracts were apoptosis, organismal death, cell viability and migration with consistency score of 13568 (Figure 7). This analysis presumed SP1, IL5 and FGF2 activate 5 to 7 proteins and 2 functions (cell viability and migration) whereas they inhibit organismal death and apoptosis. In addition, the proto-oncogene tyrosine-protein kinase receptor (RET) regulator leads to activation of neurite growth (consistency score = 1732) through increased expression of dihydropyrimidinase-related protein 2 (DPYSL2), heat shock protein A8 (HSPA8) and neurosecretory protein VGF (VGF) (Figure 7). RET regulator also activates RET itself.

Bioinformatic Analysis
Using both PANTHER gene classification and Ingenuity Pathway Analysis (IPA) software, we found that the main biological processes, in which the deregulated proteins of hill and plain extract treated SH-SY5Y cells are involved, included "biogenesis," "metabolic process", and "cellular process." All proteins found differentially expressed in hill extract treated cells were analyzed by IPA to discover the most enriched canonical pathways, possible upstream regulators and downstream effects. The software generated two main networks, "Gene expression, RNA damage and repair, RNA post-transcriptional modifications" and "Cellular development, cellular growth and proliferation, nervous system development and function" with 46 and 32 score values, respectively.
The downstream analysis unveiled key biological and cellular functions such as organismal survival, nervous system development, cell death and survival, gene expression and RNA damage and repair. Moreover, all proteins found differentially expressed concurred in an upstream analysis to predict activation or inhibition of potential transcription factors or molecules. Table 6 shows a list of top upstream regulators based on p-value. Moreover, an algorithm connection among upstream regulators, data set molecules and downstream functions or diseases, generated four main regulators with high consistency scores. The following regulators, fibroblast growth factor 2 (FGF2), interleukin-5 (IL5) and transcription factor Sp1 (SP1) revealed that the functional categories most impacted by the presence of the hill EVOO extracts were apoptosis, organismal death, cell viability and migration with consistency score of 13568 (Figure 7). This analysis presumed SP1, IL5 and FGF2 activate 5 to 7 proteins and 2 functions (cell viability and migration) whereas they inhibit organismal death and apoptosis. In addition, the proto-oncogene tyrosine-protein kinase receptor (RET) regulator leads to activation of neurite growth (consistency score = 1732) through increased expression of dihydropyrimidinase-related protein 2 (DPYSL2), heat shock protein A8 (HSPA8) and neurosecretory protein VGF (VGF) (Figure 7). RET regulator also activates RET itself

Validation of BDNF and DPYSL2 Expression in SH-SY5Y Cells Treated with Hill and Plain EVOO Extracts
Western blot analysis was used to validate the different expression level of the BDNF protein in hill and plain extract treated cells as suggested by RT-PCR data. Moreover, the different expression level of DPYSL2 was assayed to validate 2DE data. For each tested protein, the optical density of the immunoreactive band was normalized by the optical density of the β-actin immunoreactive band. Representative immunoblots and bar graphs of normalized density values are shown in Figure 8. A single immunoreactive band with apparent molecular weight of approximately 28 kDa and 71 kDa was obtained for BDNF and DPYSL2, respectively. In particular, the protein recognized by the anti-BDNF antibody likely corresponds to unglycosilated pre-pro-BDNF. BDNF and DPYSL2 expression was affected by cell treatment with EVOO extracts, but only the hill extract caused a significant increased expression of BDNF compared both to control and plain extract. On the other hand, both EVOO extracts were able to significantly increase DPYSL2 expression (p < 0.001) compared to control.

Validation of BDNF and DPYSL2 Expression in SH-SY5Y Cells Treated with Hill and Plain EVOO Extracts
Western blot analysis was used to validate the different expression level of the BDNF protein in hill and plain extract treated cells as suggested by RT-PCR data. Moreover, the different expression level of DPYSL2 was assayed to validate 2DE data. For each tested protein, the optical density of the immunoreactive band was normalized by the optical density of the β-actin immunoreactive band. Representative immunoblots and bar graphs of normalized density values are shown in Figure 8. A single immunoreactive band with apparent molecular weight of approximately 28 kDa and 71 kDa was obtained for BDNF and DPYSL2, respectively. In particular, the protein recognized by the anti-BDNF antibody likely corresponds to unglycosilated pre-pro-BDNF. BDNF and DPYSL2 expression was affected by cell treatment with EVOO extracts, but only the hill extract caused a significant increased expression of BDNF compared both to control and plain extract. On the other hand, both EVOO extracts were able to significantly increase DPYSL2 expression (p < 0.001) compared to control.

Discussion
Quercetano is an endemic cultivar typical of a specific area of Tuscany named "Piana Versiliese". Quercetano is a millennial tree characterized by small fruits resistant to olive fly thanks to their slow maturation. Due to the characteristics of "Piana Versiliese" Quercetano olive trees grow both in the hill and plain and are influenced by changes in microclimate, soil composition and water availability. Here, we report the characterization of two Quercetano EVOO extracts, hill and plain and their different neuroprotective activity in SH-SY5Y cells.
The different geographical origin of the orchards influenced both the total phenol content and the quantity of single phenols. These results are in agreement with other studies that evidenced how environmental factors such as temperature, soil moisture, water status, light and nitrogen, affect the biosynthesis of phenols and their bioaccumulation [48][49][50][51].
In particular, the plain extract presented a higher content of phenyl ethyl alcohols, cinnamic acids, oleacein (3,4-DHPEA-EDA) and oleocanthal (p-HPEA-EDA) and flavones, meanwhile hill extract was richer in lignans and oleuropein and ligstroside aglycons derivatives (3,4-DHPEA-EA and p-HPEA-EA, respectively). These data are in agreement with the results of Agiomyrgianaki et al. [52], that investigated the composition of EVOO extracts obtained from olives of the same cultivar grown in different geographical areas of Greece. They observed that the total hydroxytyrosol and lignans, pinoresinol and

Discussion
Quercetano is an endemic cultivar typical of a specific area of Tuscany named "Piana Versiliese". Quercetano is a millennial tree characterized by small fruits resistant to olive fly thanks to their slow maturation. Due to the characteristics of "Piana Versiliese" Quercetano olive trees grow both in the hill and plain and are influenced by changes in microclimate, soil composition and water availability. Here, we report the characterization of two Quercetano EVOO extracts, hill and plain and their different neuroprotective activity in SH-SY5Y cells.
The different geographical origin of the orchards influenced both the total phenol content and the quantity of single phenols. These results are in agreement with other studies that evidenced how environmental factors such as temperature, soil moisture, water status, light and nitrogen, affect the biosynthesis of phenols and their bioaccumulation [48][49][50][51].
In particular, the plain extract presented a higher content of phenyl ethyl alcohols, cinnamic acids, oleacein (3,4-DHPEA-EDA) and oleocanthal (p-HPEA-EDA) and flavones, meanwhile hill extract was richer in lignans and oleuropein and ligstroside aglycons derivatives (3,4-DHPEA-EA and p-HPEA-EA, respectively). These data are in agreement with the results of Agiomyrgianaki et al. [52], that investigated the composition of EVOO extracts obtained from olives of the same cultivar grown in different geographical areas of Greece. They observed that the total hydroxytyrosol and lignans, pinoresinol and 1acetoxytyrosol, were the most influential variables in discriminating the Greek olive oil samples according to the different geographical divisions.
The potential antioxidant activity of the extracts has been investigated using the neuroblastoma cell line SH-SY5Y differentiated with retinoic acid. These cells are considered as neuronal precursors and differentiate into more mature neuronal phenotypes under selected growth conditions such as retinoic acid [53,54]. For this reason, they are widely utilized in in vitro studies to dissect out pathogenetic mechanisms of neurodegenerative disorders [55,56].
Of note, the treatment with both extracts increased cell viability of differentiated SH-SY5Y cells. This effect is probably due to an enhancement of the mitochondrial respiration. In fact, MTT assay measures cell viability in terms of reductive activity as enzymatic conversion of the tetrazolium compound to water insoluble formazan crystals by dehydrogenases occurring in the mitochondria of living cells [57]. This hypothesis is supported by the data of Grewal et al. [58], that demonstrated the treatment with ligstroside and, at a lower extent, with 3,4-DHPEA-EDA enhanced mitochondrial respiration in SH-SY5Y-APP695 cells, a cellular model of early AD.
Both extracts showed a good activity in counteracting H 2 O 2 induced oxidative stress, even if the hill extract was significantly more effective than the plain extract. In particular, the hill extract showed a higher ability to increase cell viability in respect to cells exposed to peroxide and to decrease LDH release. On the other hand, the two extracts had comparable effect in reducing intracellular ROS levels and increasing GSH levels.
To better characterize the antioxidant mechanism underpinned by this higher protective activity of the hill extract, we evaluated the expression of four fundamental antioxidant enzymes, namely HMOX1, NQO1, TXNRD1 and GR. HMOX1 catalyzes the degradation of heme to biliverdin, which is subsequently converted to bilirubin by biliverdin reductase and both biliverdin and biliverdin reductase have antioxidant and anti-inflammatory activity [59,60]. NQO1 catalyzes the two-electron reduction of quinones and protects cells from the dangerous effects of semiquinones produced by the one-electron reduction of quinones catalyzed by cytochrome P450 reductase. The semiquinones are known to enter in the redox cycling to generate oxygen free radicals [61]. TXNRD1 is a member of the pyridine nucleotide disulfide oxidoreductase family and catalyzes the reduction of oxidized thioredoxin, using NADPH as the electron donor [62]. Glutathione reductase catalyzes the NADPH-dependent reduction of oxidized glutathione to reduced glutathione and plays a key role in providing adequate levels of reduced GSH [63]. Interestingly, the expression of these antioxidant enzymes is regulated by nuclear factor erythroid 2-related factor 2 (Nrf2) through the binding to the antioxidant response elements (ARE), a characteristic sequence present in the promoter region of antioxidants and phase II enzymes including HMOX1, NQO1, TXNRD1 and GSR [64][65][66]. The hill extract demonstrated a markedly higher ability in boosting up the endogenous antioxidant defense system in respect to the plain extract. Indeed, the hill extract upregulated all four enzymes meanwhile the plain extract just slightly increased the expression of HMOX1 and NQO1. We hypothesize that the different effects of the extracts on the expression of these enzymes could be related to their different phenol composition. The hill extract is characterized by a higher content of the lignans, pinoresinol and acetoxypinoresinol, oleuropein aglycone and ligstroside aglycon. Li et al. [67] showed pinoresinol activates the transcription factor Nrf2 and, in agreement with our data, increases the expression of one of its downstream targets, NQO1. Ligstroside derivatives belong to the secoiridoid class of biophenols, which are molecules exclusively present in plants belonging to the Olearaceae group including Olea europaea [68]. In SH-SY5Y-APP695 cells, Grewal et al. [58] demonstrated ligstroside increases the activity of glutathione peroxidase 1, another downstream target of Nrf2 suggesting the higher expression of antioxidant enzymes induced by the hill extract could be associated to a higher presence of compounds able to activate the Nrf2 pathway. Other plant-derived compounds including curcumin, sulforaphane, gamma oryzanol and resveratrol exert neuroprotective and antioxidant activity via Nrf2 [69][70][71][72]. Gamma oryzanol, in particular, could be considered for future studies in comparison to EVOOs since found mainly in rice bran oil.
Moreover, the hill extract, but not the plain extract, significantly increased the expression of BDNF. The neurotrophin BDNF is one of the most studied and well characterized neurotrophic factor in the CNS. BDNF interacts with tropomyosin-related kinase B (TrkB) and p75 cellular receptors and supports maintenance of normal brain function, neurite outgrowth and synaptic plasticity [73,74]. Several lines of evidence indicate BDNF levels play an important role in the pathogenesis of several neurodegenerative diseases, including Alzheimer's disease (AD) [75]. In the brain of AD patients or aged subjects, BDNF expression levels are lower compared with healthy subjects [76,77]. In vitro, BDNF protects neurons against Aβ 1-42 and Aβ [25][26][27][28][29][30][31][32][33][34][35] induced toxicity [78] and promotes the dephosphorylation of tau protein [79].
Given the higher neuroprotective activity shown by the hill in respect to the plain extract in SH-SY5Y cells, we performed a proteomic analysis to identify further molecular pathways modulated by these extracts. Proteomic data suggested a beneficial effect of both extracts. In fact, both hill and plain extracts induced differential expression of proteins belonging to categories involved in metabolic and cellular processes, biogenesis and cellular component organization and response to stimulus. Particularly, an increase of proteins involved in redox balance, antioxidant defenses (CS, NDUFAF7, UQCRC2) and mitochondrial morphology (IMMT) was observed while cancer promotion proteins (DIABLO, TCPT) decreased. These expression changes agree with our functional results and documented roles of EVOO polyphenols in protecting neuronal cells [80]. Indeed, these polyphenols possess the ability to restore the redox balance and hence, the optimal neuronal function, not only as antioxidants but also as mild pro-oxidants, with ensuing upregulation of cell antioxidant defenses [80]. In consideration of the major proteome changes induced by treating SH-SY5Y cells with the hill extract, we focused our attention on related findings. In cells treated with the hill extract, a large part of the eighty-nine exclusive deregulated proteins belongs to classes of RNA binding, translational and cytoskeletal proteins. Of note, the highest expression increase was observed for proteins, which are involved in cytoplasmic RNA trafficking, pre-mRNA splicing (HNRNPA3, NONO) and neuronal differentiation and maintenance (ELAVL3). IPA enrichment analysis unveiled cellular processes and components affected by the hill EVOO extract. In particular, IPA revealed inhibition of such signaling cascades as those triggered by eNOS, VEGF and ILK suggested the hill extract components were able to induce a deregulation of neuronal cell proteins able to counteract oxidative stress, proliferative and angiogenetic effects potentially promoted by their activation. On the other hand, IPA analysis also revealed activation of Unfolded Protein Response (UPR) and BAG2 pathways, resulting from enhanced transcription of ER chaperones (HSPA1A/HSPA1B, HSPA4, HSPA5, HSPA8), folding enzymes (UBQLN2) and other components of the protein degradation machinery (PSMA2, PSMD2) therefore, suggesting a potential role of the hill extract in reducing the accumulation of misfolded proteins. In addition, IPA upstream regulator analysis pointed out activation of four main regulators (FGF2, IL5, SP1 and RET) with high confidence score. Specifically, FGF2, IL5 and SP1, through twelve target proteins, promote cell viability and migration and, on the other side, inhibit apoptosis and organismal death confirming the protective role of hill extract. In addition, a regulatory role of SP1 for neuronal development and synaptogenesis has been suggested in mice astrocytes [81].
Further to this point, RET regulator activation suggested by overexpression of dihydropyrimidinase-related protein 2 (DPYSL2), neurosecretory protein VGF (non-acronymic) and heat shock cognate 71 kDa protein (HSPA8) induced activation of neurite growth function, also called neuritogenesis. Neuritogenesis is a complex task regulated by an interplay of various neurotrophic factors. Recently, Moutinho et al. [82] have demonstrated that a human VGF-derived antidepressant neuropeptide promotes neurite outgrowth in SH-SY5Y cells in association with BDNF suggesting the important role of these peptides in neuroplasticity linked with learning, memory, depression and chronic pain. Moreover, overexpression of DPYSL2 can induce neural stem cells to differentiate into neurons. In fact, DPYSL2 is considered a novel regulator for neural stem cell differentiation in rats playing a pivotal role in neuronal development and polarity [83]. Similar effects have been also described for the above mentioned ELAVL3. Indeed, Ogawa et al. [84] observed a Elavl3 -/mouse model develops a deficit in axonal transport, abnormalities in neuronal polarity of Purkinje cells and a slowly progressive axonal degeneration. Therefore, the increase of VGF, DPYSL2 and ELAVL3 observed in our cells suggests the hill extract may favor a reprogramming of neuronal cells to activate the neurogenesis machinery which, in turn, may promote a protection from neuronal injury. Concerning neurogenesis, a positive synergic effect of FGF2 and BDNF has been described in dopaminergic motoneurons [85]. In animal models of Parkinson's disease (PD), BDNF enhances the survival of dopaminergic neurons, improves dopaminergic neurotransmission and motor performance [86]. Moreover, as previously underlined, neurotrophins prevent cell death and support neuronal proliferation and maturation enhancing both growth and function of affected neurons in Alzheimer's disease (AD) and PD [75,86,87]. Despite the under threshold of BDNF z-score value derived by IPA, both transcript and protein expression analysis by PCR and WB confirmed a significant increase of the neurotrophic factor expression in SH-SY5Y cells treated with the hill extract.
Overall, these results suggest a peculiar effect of the hill compared to the plain extract, particularly for triggering the over-expression of key antioxidant enzymes, regulating the expression of proteins involved in neuronal plasticity and potentially activating neurotrophic factors. Both quality and particularity of EVOO depend on several factors [88][89][90] and among them the cultivar is the main. However, we demonstrated that the environmental conditions and soil can deeply affect even the composition of EVOOs obtained from the same cultivar. In respect to the polar polyphenolic composition, the hill extract was characterized by a higher percentage of lignanes and of two secoiridos, namely 3,4-DHPEA-EA and p-HPEA-EA, even if the plain extract contains a higher percentage of total secoiridoids.
Recent evidence suggests a correlation between specific components of extracts and neuronal plasticity. Yu et al. [91] observed that lignan pinoresinol improves memory impairment in a mouse model of dementia influencing the regulatory mechanisms involved in synaptic plasticity, thus making it a promising candidate agent to treat AD. In addition, a neuroprotective and anti-inflammatory activity of pure lignans or lignans extracted from Viburnum erosum have been described [92]. Moreover, different studies [58,[93][94][95] have described the comparable effects of secoiridoids. In fact, ligstroside-fed mice show improved spatial working memory and enhanced cognitive function [58]. In addition, both in aged and early AD mice, ligstroside expands the lifespan with outstanding performance on mitochondrial bioenergetics [58]. Finally, a role of oleuropein aglycone on α-synuclein aggregation has been reported [94]. Indeed, a recent study has highlighted oleuropein aglycone stabilizes the monomeric α-synuclein and favors the growth of non-toxic aggregates [94]. The evidence reported in literature, which assign lignans and secoiridoids among promising candidates to treat brain disorders, are in line with molecular and function changes observed in SH-SY5Y cells treated with the hill EVOO extract.

Conclusions
Our data highlight the impact of climatic areas on EVOO phenolic composition and, consequently, how this influences the antioxidant and neuroprotective effects of EVOO extracts. In particular, the hill EVOO extract was more effective in counteracting the oxidative stress and boosting up the endogenous antioxidant system in respect to the plain EVOO extract. Moreover, proteomic analysis revealed that the hill extract also modulates important molecular pathways related to neuronal survival and neurogenesis. Future studies are needed to characterize the role of individual compounds whose content is significantly higher in hill extract, in respect to the plain extract, in the biological activity observed.

Data Availability Statement:
The data used to support the findings of this study are available upon request to the authors.