Next Article in Journal
A Study on the Antimicrobial and Antibiofilm Peptide 1018-K6 as Potential Alternative to Antibiotics against Food-Pathogen Salmonella enterica
Next Article in Special Issue
Sustainability of the Olive Oil System
Previous Article in Journal
Dietary Curcumin Alleviated Acute Ileum Damage of Ducks (Anas platyrhynchos) Induced by AFB1 through Regulating Nrf2-ARE and NF-κB Signaling Pathways
Previous Article in Special Issue
Distinct Particle Films Impacts on Olive Leaf Optical Properties and Plant Physiology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioactive Potential of Minor Italian Olive Genotypes from Apulia, Sardinia and Abruzzo

1
Institute of Biosciences and BioResources, National Research Council (IBBR-CNR), Via Amendola 165/A, 70125 Bari, Italy
2
Spin off Sinagri s.r.l., University of Bari Aldo Moro, Via Amendola 165/A, 70125 Bari, Italy
3
Department of Soil, Plant and Food Sciences, University of Bari Aldo Moro, Via Amendola 165/A, 70125 Bari, Italy
4
Faculty of Bioscience and Technology for Agriculture, Food and Environment, University of Teramo, Via Renato Balzarini 1, 64100 Teramo, Italy
5
Dipartimento di Agraria, University of Sassari, Viale Italia 39/A, 07100 Sassari, Italy
6
Institute for Sustainable Plant Protection–Support Unit Bari, National Research Council (IPSP-CNR), Via Amendola 165/A, 70125 Bari, Italy
*
Author to whom correspondence should be addressed.
Foods 2021, 10(6), 1371; https://doi.org/10.3390/foods10061371
Submission received: 30 April 2021 / Revised: 5 June 2021 / Accepted: 7 June 2021 / Published: 14 June 2021
(This article belongs to the Special Issue Sustainability of Olive Oil System)

Abstract

:
This research focuses on the exploration, recovery and valorization of some minor Italian olive cultivars, about which little information is currently available. Autochthonous and unexplored germplasm has the potential to face unforeseen changes and thus to improve the sustainability of the whole olive system. A pattern of nine minor genotypes cultivated in three Italian regions has been molecularly fingerprinted with 12 nuclear microsatellites (SSRs), that were able to unequivocally identify all genotypes. Moreover, some of the principal phenolic compounds were determined and quantified in monovarietal oils and the expression levels of related genes were also investigated at different fruit developmental stages. Genotypes differed to the greatest extent in the content of oleacein (3,4-DHPEA-EDA) and total phenols. Thereby, minor local genotypes, characterized by stable production and resilience in a low-input agro-system, can provide a remarkable contribution to the improvement of the Italian olive production chain and can become very profitable from a socio-economic point of view.

1. Introduction

Olive (Olea europaea ssp. europaea) is considered among the historically and traditionally most important crops in the world and especially in the Mediterranean basin, where it has been cultivated for centuries not only for nutrition but also for cultural and religious reasons. In particular, the strategic geographical position of Italy in a temperate area has favored olive cultivation and therefore the enrichment of its germplasm over time, which is estimated to include about 800 cultivars [1,2]. The great diffusion of olive trees in the Italian territory highlights its importance for production, economy and local traditions, as proven and documented by innumerable historical catalogues and archives up to the XIII century B.C. [3].
As a consequence of the general awareness about the loss of plant genetic diversity and the drastic climate change currently underway, the attention of the scientific community has been recently put toward more sustainable agriculture, as opposed to the intensive, mono-cultivar farming systems, and to the safeguard of plant biodiversity as source of new interesting traits. Similarly to other cultivations, also for olive, a valid and intriguing opportunity to help in overcoming these issues is offered by the unexploited or still poorly characterized germplasm, including local or minor genotypes, i.e., autochthonous landraces generally spread at regional level and well adapted to specific pedoclimatic conditions in traditional groves with very low agronomic input [4]. These genotypes could represent an interesting reservoir of useful traits, such as nutraceutical and antioxidant compounds and/or resistances to environmental stresses, that makes their diffusion in the market and their utilization in breeding programs extremely promising. In this perspective, the first step for the valorization of these genotypes is represented by their fine characterization that provides the genetic identity of both plant material and derived products (drupes and oils) and that guarantees an unequivocal varietal recognition. Microsatellite markers (SSRs) are among the most simple, fast and economic molecular tools widely used for olive varietal identification [5,6,7,8] and for food tracking and tracing [9,10,11,12,13].
The extra virgin olive oil is considered a fundamental element of the Mediterranean diet [14], thus representing one of the pillars of the Italian economy. Being the most important fat source in the human diet, olive oil is characterized by a distinctive high content of mono-unsaturated fatty acids (MUFA) and by a certain level of secondary metabolites with nutraceutical properties, such as polyphenols [15]. The most abundant polyphenol class in olive oil is represented by the secoiridoids, i.e., complex molecules including oleuropein, ligstroside and their derivates (oleacein, oleocanthal, an isomer of oleuropein aglycone called 3,4-DHPEA-EA and the ligstroside aglycone also known as p-HPEA-EA), the latter compounds being produced during olive oil mechanical extraction. Polyphenols also include two important phenolic alcohols worthy to be mentioned for their antioxidant power, i.e., tyrosol and hydroxy-tyrosol [16].
A huge literature has been dedicated to the healthy value of polyphenols [17,18,19,20]. Numerous research has demonstrated their positive effects in contrasting cardiovascular diseases thanks to their anti-inflammatory and -oxidative actions [21,22,23,24]. In particular, polyphenols are able to bind low-density lipoprotein (LDL), avoiding its oxidation [25], modulate angiogenic responses [26], protect against endothelial dysfunction [27] and contribute to decrease the blood pressure [28]. Moreover, the positive effect of tyrosol and hydroxy-tyrosol against reactive oxygen species has been deeply investigated [29,30,31,32].
Furthermore, some phenolic compounds, especially oleuropein, ligstroside aglycone and derived molecules, are responsible for the pungency and bitter taste of olive oil [33,34,35,36]. For this reason, in recent decades, breeders and researchers focused on polyphenols as new trait of interest to satisfy the increasing request of high-quality olive oil for nutritional and organoleptic aspects [37,38,39].
The biosynthetic pathway of olive secoiridoids has been studied and proposed by several authors, even though it is not completely clarified yet. Transcriptomic and metabolomic analysis of developing olive fruits and leaves in different cultivars have shed light on the enzymes involved in polyphenol metabolism [40,41,42,43,44].
The main aim of this work was the valorization of some minor olive genotypes and cultivars, autochthonous of Central and Southern Italy, that are still poorly studied, but interesting for both their resilience and adaptation to low-input agriculture and for the quality of olive oil. The application of SSR markers has been necessary to obtain a clear molecular fingerprinting of these genotypes, as a recognized tool of primary importance for certifying plant material production, varietal tracing and authenticity testing. Moreover, an integrated approach of oil biochemical characterization, with particular focus on tyrosol, hydroxy-tyrosol and oleacein, coupled with a genetic expression study of key enzymes involved in the first steps of their biosynthesis, has contributed to clarify the molecular mechanisms underlying the olive polyphenol biosynthetic pathway and to further valorize the Italian olive germplasm.

2. Materials and Methods

2.1. Plant Material and Olive Sampling

In total, nine olive cultivars were selected for this work; they were grown in groves located in three different Italian regions: Sardinia, Abruzzo and Apulia.
The Sardinian cultivars “Sivigliana da olio”, Semidana and “Corsicana da olio” have been chosen for their suitability to be processed as table olives (Sivigliana da olio), for their resistance to some common diseases and high productivity (Corsicana da olio) and for the increasing interest of growers to introduce it in all Sardinian orchards (Semidana) [45,46]. Specialized non-irrigated orchards were located in the province of Sassari (Sardinia, Italy); olives of Sivigliana and Semidana cultivars were hand harvested in Ittiri (40°37′15.7″ NL 8°32′20.7″ EL), while samples of Corsicana were picked in Usini (40°39′39.0″ NL 8°31′32.7″ EL).
With regard to Abruzzo (Central Italy), the varieties Dritta, Tortiglione and Gentile dell’Aquila were selected as the most representative of the provinces of Pescara, Teramo and L’Aquila, respectively, also in consideration of their specific adaptation to the local pedoclimatic conditions correlated with peculiar agronomic and biochemical traits [47,48]. Dritta and Tortiglione were grown in orchards located in Notaresco, in the province of Teramo (42°38′54.6″ NL, 13°52′57.4″ EL), while Gentile dell’Aquila was grown in orchards from Vittorito, in the province of L’Aquila (42°08′16.6″ NL, 13°48′37.3″ EL).
In the Apulia region (Southern Italy), in addition to some well-known and notable cultivars for olive oil extraction, a great number of other varieties and/or landraces are present [3]. Among these, three have been considered in this study, specifically the varieties Bambina, Oliva Rossa and Cima di Melfi, for which few reports are currently available in the literature [49,50]. The Apulian genotype Bambina was cultivated in the orchard located in Gravina di Puglia (40°49′0″ NL, 16°25′0″ EL), while genotypes Cima di Melfi and Oliva Rossa were grown in the orchard of Putignano (40°51′0″ NL, 17°7′0″ EL), both in the province of Bari.
For the molecular characterization, young leaves of each cultivar were collected and stored at −20 °C before use. For gene expression analysis, drupes were harvested at different developmental stages of fruit, generally from mid-October to the beginning of November: T1, yellow-green olives; T2, turning olives; and T3, almost dark olives (Table 1). In order to minimize the effects of fruit asynchronous maturation within the same tree, drupes were harvested from the external parts of the canopy of trees. After harvesting, olives were immediately frozen in liquid nitrogen and stored at −80 °C until further processing. For two of the considered time-points (T2 and T3), one part of the collected drupes was immediately used for virgin olive oil extraction and biochemical characterization. An additional timepoint at 100% ripening (T4, fully dark olives, from the end of November to the beginning of December) was considered exclusively for oil production and biochemical characterization (Table 1). For each cultivar/genotype, three biological replicates were taken. Sampling was carried out for two consecutive years.

2.2. Genetic Characterization

Leaf samples were collected from three plants per cultivar, lyophilized and finely pulverized before use. Total genomic DNA was extracted from 200 mg of dry tissue according to [51], checked for quantity and quality and normalized to 50 ng µL−1. Olive genotypes were molecularly characterized by means of 12 nuclear SSRs (Table S1) [52,53,54], chosen on the basis of their suitability and reliability proven in several studies about olive variety identification [55,56,57]. Primer pairs were synthesized by Thermo Fisher Scientific (Waltham, MA, USA) and all forward primers were labeled with one of the following dyes: 6FAM™, NED™, VIC® and PET™. The amplification reactions were carried out in a final volume of 12.5 µL, using a T100 thermal cycler (Bio-Rad Laboratories, Segrate, MI, Italy) according to [9]. Two microliters of each PCR product was added to 0.5 µL of GeneScan™ 600 LIZ® Size Standard (Applied Biosystem, Foster City, CA, USA) and 9.5 µL of Hi-Di Formamide (Applied Biosystem, Foster City, CA, USA), and successively separated by capillary electrophoresis using an automatic sequencer ABI PRISM 3100 Avant Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Detection, sizing and data collection were carried out by means of the GeneMapper® genotyping software v.5.0 (Applied Biosystems, Foster City, CA, USA) as in [58]. To estimate the genetic distances among the considered genotypes, cluster analysis based on the Unweighted Neighbor Joining method was performed using DARWIN software v. 6.0.010 (http://darwin.cirad.fr, accessed on 26 April 2019), with 1000 bootstrap values for tree construction.

2.3. Oil Extraction and Characterization

Around 10 kg of drupes were collected from three trees per each sampling point and variety. Fruits were then divided into three aliquots representing the biological replicates. Oil extraction was performed within 5 h after harvesting at laboratory scale. The extraction system was made up of a semi-industrial scale hammer crusher (RETSCH GmbH 5657, Haan, Germany) working at 2850 rpm and a basket centrifuge with a bowl of 19 cm working at 2700 rpm (Marelli Motori S.p.A., Arzignano, VI, Italy) [59]. Briefly, about 1 kg of olives per sample and replicate was crushed and then the olive paste was transferred into the basket centrifuge for the oil recovery. Once extracted, the virgin olive oils (n = 3) were stored at 20 °C in 100 mL dark glass bottles until the analyses.
Two phenolic alcohols (i.e., hydroxyl-tyrosol and tyrosol) and the oleacein (3,4-DHPEA-EDA) were chosen as target phenolic compounds for HPLC analyses, being amongst the most abundant ones in olive oils and because of their powerful antioxidant activity and impact on oil sensorial feature [60]. These phenolic compounds were extracted by liquid–liquid extraction using a mixture of methanol/water (70/30 v/v) according to previous papers [61]. The extraction procedure was similar for both the Folin–Ciocalteu assay and for HPLC analysis, with the only difference being that in the latter case 250 μL of a 100 mg kg−1 solution of gallic acid as internal standard for quantification was added. Total phenolic compounds (TPC) were quantified spectrophotometrically [61] by means of a calibration curve of pure gallic acid and the results are expressed as gallic acid equivalent (GAE, mg kg−1). HPLC-DAD analysis was carried out as previously reported [61] using a UHPLC binary system (Dionex Ultimate 3000 RSLC, Waltham, MA, USA). The identification of hydroxy-tyrosol, tyrosol and 3,4-DHPEA-EDA was performed by comparing the peak retention times with those obtained by the injection of pure standards and/or with data in the literature [62]. The results are expressed as gallic acid equivalent (GAE, mg kg−1).

2.4. Polyphenolic Compound Gene Expression

The mesocarps of frozen olives were mechanically crushed by the use of a tissue lyser and total RNA was extracted from 100 mg of the obtained powder according to the manufacturer’s instructions of Spectrum Plant Total RNA Kit (Sigma-Aldrich, St. Louis, MO, USA). An additional step for on column genomic DNA digestion was added. RNA quantity and quality were checked by spectrophotometric measurement using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and by electrophoresis on 1.2% Certified Molecular Biology Agarose gel (Bio-Rad Laboratories, Segrate, MI, Italy) in 1X TBE buffer (1 M Trizma base, 1 M Boric Acid, 20 mM EDTA, pH 8.3).
Quantitative real-time polymerase chain reactions (qRT-PCRs) were carried out on three genes involved in the polyphenolic biosynthesis, named TYRD for tyrosine/dopa decarboxylase, CuAO for copper amine oxidase, and ALDH for alcohol dehydrogenase. The used primer pairs are listed in Supplementary Material (Table S2). The reverse transcription of 700 μg RNA samples was performed with the SuperScript™ VILO™ cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. qRT-PCR reactions were carried out using the SsoAdvanced Universal SYBER® Green Supermix (Bio-Rad Laboratories, Segrate, MI, Italy) and the CFX96 Touch Real-Time PCR Detection System (Bio-Rad Laboratories, Segrate, MI, Italy). Thermal cycling parameters were: initial denaturation at 95 °C for 3 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. The specificity of the amplification product per each primer pair was confirmed by evaluating the melting curve through an increase of 0.2 °C every 5 s from 65 to 95 °C. For qRT-PCR assay, each amplification reaction was run in triplicates of three biological replicates. The elongation factor 1α (EF1α, AM946404) was selected as a reference gene for normalization and the comparative Ct method (2−ΔΔCt method) was used to analyze the expression levels of the selected genes [63].

2.5. Statistical Analysis

All values of chemical and genetic analysis are means (n = 3) ± standard deviation (sd). The statistical analysis was performed, for samples of each region, by one-way analysis of variance (ANOVA) with sampling time as the group factor and post-hoc Fisher’s LSD test, using the software Statistica 10.0 for Windows. Differences were considered to be significant when p < 0.05.

3. Results

3.1. Genetic Characterization

The chosen nuclear microsatellites successfully amplified all the olive cultivars under analysis, with identical genetic profiles among the replicates of each cultivar. With the exception of missing data for only three loci related to the markers DCA09, DCA17 and GAPU101, unique alleles for all the used SSR were assigned to each cultivar, thus allowing an accurate varietal identification (Table 2 and Figure S1).
The Unweighted Neighbor Joining dendrogram allowed us to group the genotypes predominantly on the basis of their regional origin, except one. Indeed, as expected, three main clusters were identified (Figure 1). Cluster I included two of the three analyzed Abruzzo cultivars, i.e., Tortiglione and Dritta. The third Abruzzo cultivar, Gentile dell’Aquila, was included in Cluster II together with all the Sardinian cultivars (Semidana, Sivigliana and Corsicana). Finally, Cluster III encompassed cultivars from the Apulia region, i.e., Bambina, Cima di Melfi and Oliva Rossa.

3.2. Monovarietal Olive Oil Phenolic Content

This research focused on three of the main active constituents of olive polyphenols, i.e., tyrosol, hydroxy-tyrosol and oleacein (3,4-DHPEA-EDA), and on the TPC of monovarietal oils. In general, the concentration and the main trend of these compounds varied among the analyzed cultivars. For example, tyrosol content was higher in cultivar Tortiglione than in Dritta and Gentile dell’Aquila with no significant changes (p > 0.05) during ripening, while it significantly increased in Sardinian cultivars and in two of the Apulian genotypes, i.e., Bambina and Oliva Rossa. An opposite trend was observed in cultivar Cima di Melfi, that showed a significant drop in tyrosol content at complete maturity of drupes (Table 3; Figure 2A).
Hydroxy-tyrosol concentration was the lowest among the three studied compounds, generally showing its highest amount in two time-points, that is at T2 and T4. Among all the analyzed cultivars, two Sardinian cultivars (Semidana and Sivigliana) had the highest hydroxy-tyrosol content, with values of about 0.9 and 1.2 mg kg−1 olive oil, respectively. In particular, in Sivigliana cultivar, this compound first decreased when drupes were in advanced stage of ripening (T3) and then significantly increased at the end of fruit maturation. An exception to this behavior was recorded for cultivars Oliva Rossa and Corsicana, whose hydroxy-tyrosol contents were slightly higher at T3 (Table 3; Figure 2B).
The concentration of oleacein generally decreased with maturity stage in all the considered cultivars, with the only exception of Cima di Melfi and Gentile dell’Aquila, that showed a more evident and significant increase in this compound in fully ripe fruits (Table 3; Figure 2C). Cultivars Dritta and Bambina were the genotypes with the lowest amount of 3,4-DHPEA-EDA in their oils. Bambina VOO was already known for its particular phenolic pattern, which shows a significant contribution of flavonoids and a lesser amount of secoiridoid derivatives [49].
Finally, changes in the total phenol content (TPC) were generally recorded in most of the cultivars. Indeed, the concentration of total phenols was almost constant during drupe maturation or slightly increasing at the last stage of fruit ripening (T4), except for the Abruzzo cultivar Tortiglione and the Apulian cultivar Oliva Rossa, that showed a drastic and significant drop in the total phenol content in comparison with the first stage of fruit ripening (T4 vs. T2) (Table 3; Figure 2D). Sardinian cultivars showed, on the other hand, a slight but significant decrease in TPC at the T4 sampling. Anyway, the oils obtained from these cultivars, together with oils derived from most of the Apulian genotypes, were the ones with the highest and most stable amount of total phenols during all stages of fruit ripening.

3.3. Gene Expression Analysis

The expression profiles of the selected genes did not follow the same trend among the cultivars. In most genotypes, the mRNA levels of TYRD were quite low and remained almost constant as the fruit development proceeded (Figure 3A). An exception was represented by cultivars Sivigliana and Gentile dell’Aquila, in which this gene was highly expressed at the first drupe sampling (T1) and then drastically dropped at the second sampling, albeit with no significance. On the contrary, in cultivars Semidana and Tortiglione, TYRD mRNA levels remarkably increased from the second to the third sampling, thus reaching the highest values when drupes were almost dark and mature.
In some cultivars, the CuAO expression pattern did not notably change during fruit maturation, while in other genotypes, such as Semidana and Bambina, a slight increase was observed up to the third sampling time. Interestingly, the CuAO expression level was noteworthy in one of the Apulian cultivars, i.e., Cima di Melfi, not for its general trend but specifically for its abundance during the entire drupe maturation process, since it was extremely and significantly different (higher) in comparison to that detected in all other cultivars (Figure 3B).
In the majority of genotypes, ALDH was almost exclusively expressed during the T3 stage of fruit development, with a remarkable increase starting at the turning phase of drupes, in particular in the cases of Gentile dell’Aquila and Cima di Melfi, and with a less pronounced increment observed in the cases of Bambina, Semidana and Sivigliana. On the contrary, ALDH expression levels in the remaining cultivars were very low (Figure 3C).

4. Discussion

The identification of minor olive cultivars with favorable agronomic and nutraceutical features among the available and still uncharacterized Italian germplasm, and their fine characterization with different approaches, could effectively contribute to both the valorization and the spread of these cultivars worldwide. The possibility of their introduction in the olive-growing sector is reinforced by the availability of a molecular fingerprint, that can provide a crucial tool for the identification of these genotypes and the traceability of their oils. The use of microsatellite markers in this study has allowed not only a fine characterization of the considered genotypes, but also clarification of the phylogenetic relationships among them and assessment of an identification key. The Unweighted Neighbor Joining dendrogram mostly separated the genotypes according to their geographic origin as reported in many other studies about the evaluation of the genetic diversity of the Italian olive germplasm [5,64]. Surprisingly, only one of the analyzed genotypes has escaped this kind of clustering; indeed, the Gentile dell’Aquila cultivar, belonging to the autochthonous cultivars of Abruzzo region, grouped with all the Sardinian genotypes instead of with the other cultivars from Central Italy (Figure 1). The explanation of this result can be traced back to the complex migration history of the olive. In 2015, Dìez and colleagues [65] demonstrated the existence of three different gene pools (Q1, Q2 and Q3) in the Mediterranean area and the presence of a broad ‘mosaic’ group as a mixture of Q1 and Q2. The Italian cultivars considered in that research mainly belonged to the Q2 group. Anyway, the authors postulated that the wide occurrence of the ‘mosaic’ cultivars might be indicative of admixture events. Thereby, both the outcrossing nature of O. europaea ssp. europaea as well as the exchanges related to human migration over the centuries [64,66] could explain a certain level of genetic admixture between Gentile dell’Aquila and the Sardinian cultivars or their origin from a common wild ancestor.
Over the last decades, virgin olive oil (VOO) became the symbol of good nutritional habits. This role has been mostly linked to some features such as the high content of MUFA and of phenolic compounds, which allow us to recognize VOO as a functional food [67]. As proof of that, the European Food Safety Authority (EFSA) has permitted the declaration “Olive oil polyphenols contribute to the protection of blood lipids from oxidative stress” when olive oil contains a minimum of “5 mg hydroxy-tyrosol and its derivatives per 20 g of olive oil” [68,69]. However, VOO phenolic content not only plays an important nutritional role but also has a great impact on the organoleptic features and on the stability of the product [60]. Given the wide spectrum of features directly linked to olive phenolics, it is clear how important the understanding of their synthesis and regulation is.
Polyphenol accumulation has been reported to be a complex process considerably varying among genotypes, tissues, developmental stages and in response to different agronomic and environmental conditions [70,71,72]. In accordance with other studies, there was a component of the variability of the phenolic composition observed among the analyzed monovarietal virgin oils that mostly depended on the genotype. On average, the Apulian cultivars generally showed the highest values of the detected compounds at the complete maturation stage of drupes, followed by the Sardinian cultivars. This could also be related to the different cultivation areas and environmental temperatures of the considered groves, since the climates of Apulia and Sardinia are generally similar to each other and usually milder than that of Central Italy (Abruzzo). As reported by [73], the pathways related to olive phenol compounds could be differentially modulated on the basis of altitude and temperatures, phenols biosynthesis being prolonged in temperate areas.
Quantitative differences in phenolic compounds were highlighted both among the analyzed cultivars and among fruits of the same cultivar at different ripening stages. With regard to total phenols (TPC), a sharp decrease is generally and mainly observed between June and September, as result of the oleuropein hydrolysis during drupe maturation (high β-glucosidase activity) [74,75]. Indeed, the very early stages of drupe formation are dominated by secoiridoids (above all oleuropein), while in the last ripening stages, simple phenols and flavonoids become the major components. In our study, the drupe collecting period between October and December did not allow us to highlight such remarkable differences in total phenol content, this period being characterized by a higher stability of phenol metabolism (Table 3; Figure 2D). Thereby, with the only exception of cultivars Tortiglione and Oliva Rossa, that showed a sharp decrease from T2 to T3 stage, the variation of TPC was minimal in all cultivars during the fruit maturation period considered in this study. An intermediate behavior was registered in oils obtained by Sardinian cultivars, that showed a slight, but significant, decrease only at the overripe (T4) stage.
On the contrary, the content of single compounds was very variable among cultivars and sampling times. Secoiridoid derivatives, i.e., aglycon forms of the secoiridoid glucosides usually formed during oil extraction by the enzymatic hydrolysis of oleuropein, demethyloleuropein and ligstroside [76], were reported to be the most abundant phenols in monovarietal olive oils [72]. Among those, oleacein (3,4-DHPEA-EDA) is an important compound, structurally similar to oleocanthal and with comparable pharmacological properties [77]. In our study, the oleacein levels considerably differed among all genotypes, progressively decreasing as fruit maturation proceeded (Table 3; Figure 2C). Indeed, the higher values were found at the T2 time-point, although with remarkable differences among cultivars (in the range between 2 and 63 mg/kg). As a result of an opposite behavior to this general trend, oils obtained by fully ripe drupes (T4) of the Apulian cultivar Cima di Melfi and of the Abruzzo cultivar Gentile dell’Aquila were those with the highest oleacein concentration. As already mentioned, variation in the VOOs’ phenolic composition can be ascribable to several factors, including the technological steps of olive oil extraction, which is far from being merely a mechanical process. Indeed, starting from drupe crushing, several enzymes are set free whose activities are strongly influenced by the genotype, the maturity degree and the environmental conditions (time, temperature, atmosphere) [78,79]. Moreover, secoiridoid derivatives, such as 3,4-DHPEA-EDA, are among the most affected olive phenols during the extraction process. For example, in Tuscan olive oils, oleacein was found only a long time after milling, likely due to the activity of different esterases [80]. Thus, linearity between drupe genetic expression and VOO phenols could be difficult to see and the identification of clear markers of maturity is a complex task [81].
The content of minor phenols such as tyrosol and hydroxy-tyrosol was also evaluated in the considered genotypes, due to their relevance as bioactive molecules with beneficial properties on human health [16,82]. Hydroxy-tyrosol is a simple alcohol conjugated to form oleuropein derivatives, usually highly expressed in young olive fruits [42,43]. In accordance with other reports, the hydroxy-tyrosol amount decreased from semi-green to nearly ripe drupes (with the exception of Corsicana, Oliva Rossa and Gentile dell’Aquila), but then an interesting increase was again observed in the final developmental stage, i.e., in totally black fruits (Table 3; Figure 2B). This result is generally attributed to oleuropein catabolism rather than to new hydroxy-tyrosol biosynthesis, as a consequence of the recycling of all the valuable molecules of complex compounds and their conversion to simple phenols [42,44,83]. Moreover, the oleuropein catabolism in the mature fruit is necessary to reduce the bitter taste typically associated with this compound.
Changes in tyrosol content with ripening did not follow a clear pattern among the studied varieties. In some cases (such as Semidana and Cima di Melfi) an opposite trend with the oleacein content was found. These variations could be supposed to be linked to the genetic expression of the enzymes in the biosynthetic pathway (Figure 3). Indeed, tyrosol could be initially synthesized—and mostly found—in free form, while subsequently it is used as a precursor for oleacein production (Table 3; Figure 2A). However, as already reported [81], finding a straightforward relationship between gene expression and phenotype is very difficult, as the interaction with agronomic and pedoclimatic conditions should be also taken into account. Moreover, the content of phenolic compounds depends also on the activity of other enzymes such as the olive β-glucosidases (not considered in this study) which strongly affect their release in the oily phase. Different patterns of tyrosol accumulation with respect to the cultivar and the maturity degree have been already observed by other authors [84].
When this study started, the exact polyphenol biosynthetic pathway was still not clarified in detail; therefore, among genes annotated in phenol biosynthesis, for expression analysis our choice fell on three genes (TYRD, CuAO and ALDH) on the basis of the involvement of the corresponding encoded enzymes in the first steps of this biosynthetic pathway and in relation to the compounds detected in VOOs under characterization. Tyrosine is referred to as the precursor of the secoiridoid class [85]; in particular, the biosynthesis of tyrosol and hydroxy-tyrosol has been proposed to proceed through the formation of tyramine or, otherwise, it can follow the path including L-DOPA and dopamine as intermediates (Figure 4). Initially, an enzyme similar to a tyrosine/DOPA decarboxylase (TYRD) was recognized to be responsible for the conversion of both tyrosine in tyramine and L-DOPA in dopamine [40,86] (Figure 4A). Recently, these two conversions have been clarified to take place by the action of two different enzymes, i.e., a tyrosine decarboxylase (TDC) in the first way and a DOPA decarboxylase (DCC) in the parallel path (Figure 4B) [44]. Moreover, in the previous proposed version, the subsequent steps leading to the conversion of tyramine to tyrosol and dopamine to hydroxy-tyrosol required two enzymes called CuAO and ALDH [87]. Now, the presence of other intermediates in both paths and the action of an additional enzyme, that is, a phenylacetaldehyde reductase (PAR), have been recognized (Figure 4B) [44]. The reaction of the conversion of tyrosol to hydroxy-tyrosol and vice versa seems to involve a polyphenol oxidase (PPO) [42]. Since most of the genes responsible for secoiridoid formation have been identified in a restricted number of species, some steps of the biosynthesis of oleuropein are still unclear, even though the pathway from geranyl diphosphate to secologanin has been elucidated [42,44].
The relative expression of the selected genes was monitored during three stages of fruit maturation, at T1, T2 and T3 sampling times. The expression analysis showed the general tendency of these mRNA levels to decrease or remain more or less stable in semi-green fruits and then to increase, slightly or significantly in a few cases, in semi-dark drupes (Figure 3). These results are in line with the biochemical characterization of monovarietal oils at T2 and T3 stages, which proved that at these time-points the tyrosol and hydroxy-tyrosol contents were progressively diminishing, following the same pattern as the genetic expression data. Some exceptions to this trend have been surprisingly observed for few varieties and in particular for TYRD and ALDH genes. Indeed, these two genes were highly expressed, respectively, in the Tortiglione and Semidana genotypes, and in Cima di Melfi and Gentile dell’Aquila, with a remarkable increase in their transcript levels from T2 to T3 time-points. These unexpected but interesting results, not yet reported in other olive cultivars, to our knowledge, certainly require further investigation. In order to clarify the role of TYRD and ALDH genes in these cultivars, their expression profile should be investigated at more developmental stages of drupes during more consecutive years. Protein expression experiments as well as enzymatic activity assays could also help us to understand the involvement of these genes in polyphenol metabolism of these minor cultivars. At the moment, we could hypothesize that these genotypes are effectively characterized by a high level of gene expression even towards the end of drupe maturation, contrarily to that reported in this study and, more in general, in the literature. It is also plausible that the corresponding enzymes are conveyed to other metabolisms or that some post-translational control takes place and makes these enzymes unfunctional. Exceptions apart, the general agreement of qRT-PCR results with the biochemical outcomes confirmed that hydroxy-tyrosol biosynthesis is under a transcriptional control, as also reported by other authors [40,42,43,44].

5. Conclusions

Unveiling the genetic potentialities of minor cultivars/genotypes could open the doors to their valorization. Besides the safeguarding of local varieties, the information provided in this work could be useful for breeding programs with selection of genotypes that, while being resilient and better adapted to specific environmental conditions, could also exploit their potentiality in terms of accumulation of bioactive compounds with important consequences on the nutritional, sensory and stability properties of VOO. With a view to an eco-sustainable and low-input agriculture, the introduction of minor cultivars may offer interesting features that, together with the extraction technology, could address the needs of a worldwide market with growing interest towards VOOs. Thus, the results and information obtained in this study are encouraging for the valorization of the still poorly explored Italian germplasm, with important both economic and scientific consequences.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/foods10061371/s1, Figure S1: Example of SSR electropherograms, Table S1: Microsatellite markers used in this study, Table S2: List of primers utilized for qRT-PCR analysis.

Author Contributions

Conceptualization, W.S., C.M. and F.C.; investigation, S.G., I.M., G.S., F.F., P.C. and A.P.; writing—original draft preparation, W.S.; writing—review and editing, W.S., G.S., F.F., C.D.D.M., P.C., A.P., C.M. and F.C.; supervision, W.S., C.M. and F.C.; funding acquisition, F.C., W.S. and C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was mainly funded by the Italian AGER 2 Project, grant no. 2016-0105, and it was in part supported by the CNR project NUTR-AGE (FOE-2019, DSB.AD004.271) and the Apulia Region Project “Approcci di Next Generation Sequencing per l’analisi di variabilità e di espressione genica in genotipi di olivo autoctoni pugliesi” L.R. 28 December 2018 n. 67, art. 37.

Conflicts of Interest

The authors declare that they have no competing financial or personal interests that could influence the work reported in this paper.

References

  1. Muzzalupo, I. Olive Germplasm–Italian Catalogue of Olive Varieties; InTech: Rijeka, Croatia, 2012; p. 430. [Google Scholar] [CrossRef] [Green Version]
  2. Bartolini, G.; Cerreti, S.; Briccoli Bati, C.; Stefani, F.; Zelasco, S.; Perri, E.; Petruccelli, R. Oleadb: Worldwide database for the management of genetic resources of olive (Olea europaea L.). In Proceedings of the X National Congress on Biodiversity, Rome, Italy, 3–5 September 2014; pp. 18–23. [Google Scholar]
  3. Miazzi, M.M.; di Rienzo, V.; Mascio, I.; Montemurro, C.; Sion, S.; Sabetta, W.; Vivaldi, G.A.; Camposeo, S.; Caponio, F.; Squeo, G.; et al. REGEROP: An integrated project for the recovery of ancient and rare olive germplasm. Front. Plant Sci. 2020, 11, 73. [Google Scholar] [CrossRef] [Green Version]
  4. Rotondi, A.; Cultrera, N.G.; Mariotti, R.; Baldoni, L. Genotyping and evaluation of local olive varieties of a climatically disfavoured region through molecular, morphological and oil quality parameters. Sci. Hortic. 2011, 130, 562–569. [Google Scholar] [CrossRef]
  5. Muzzalupo, I.; Vendramin, G.G.; Chiappetta, A. Genetic biodiversity of Italian olives (Olea europaea L.) germplasm analyzed by SSR markers. Sci. World J. 2014, 2014. [Google Scholar] [CrossRef] [Green Version]
  6. Mousavi, S.; Mariotti, R.; Regni, L.; Nasini, L.; Bufacchi, M.; Pandolfi, S.; Baldoni, L.; Proietti, P. The first molecular identification of an olive collection applying standard simple sequence repeats and novel expressed sequence tag markers. Front. Plant Sci. 2017, 8, 1283. [Google Scholar] [CrossRef] [PubMed]
  7. Sebastiani, L.; Busconi, M. Recent developments in olive Olea europaea L.) genetics and genomics: Applications in taxonomy, varietal identification, traceability and breeding. Plant Cell Rep. 2017, 36, 1345–1360. [Google Scholar] [CrossRef] [PubMed]
  8. El Bakkali, A.; Essalouh, L.; Tollon, C.; Rivallan, R.; Mournet, P.; Moukhli, A.; Hayat Zaher, H.; Mekkaoui, A.; Hadidou, A.; Sikaoui, L.; et al. Characterization of worldwide olive germplasm banks of Marrakech (Morocco) and Córdoba (Spain): Towards management and use of olive germplasm in breeding programs. PLoS ONE 2019, 14, e0223716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Montemurro, C.; Miazzi, M.M.; Pasqualone, A.; Fanelli, V.; Sabetta, W.; di Rienzo, V. Traceability of PDO olive oil ‘‘terra di Bari’’ using high resolution melting. J. Chem. 2015, 2015. [Google Scholar] [CrossRef] [Green Version]
  10. Sabetta, W.; Miazzi, M.M.; di Rienzo, V.; Fanelli, V.; Pasqualone, A.; Montemurro, C. Development and application of protocols to certify the authenticity and the traceability of Apulian typical products in olive sector. Riv. Ital. Sost. Grasse 2017, 94, 37–43. [Google Scholar]
  11. Xanthopoulou, A.; Ganopoulos, I.; Bosmali, I.; Tsaftaris, A.; Madesis, P. DNA fingerprinting as a novel tool for olive and olive oil authentication, traceability, and detection of functional compounds. In Olives and Olive Oil as Functional Foods: Bioactivity, Chemistry and Processing; Shahidi, F., Kiritsakis, A., Eds.; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2017; pp. 587–601. [Google Scholar]
  12. Gomes, S.; Breia, R.; Carvalho, T.; Carnide, V.; Martins-Lopes, P. Microsatellite High-Resolution Melting (SSR-HRM) to track olive genotypes: From field to olive oil. J. Food Sci. 2018, 8, 2415–2423. [Google Scholar] [CrossRef]
  13. Piarulli, L.; Savoia, M.A.; Taranto, F.; D’Agostino, N.; Sardaro, R.; Girone, S.; Gadaleta, S.; Fucilli, V.; De Giovanni, C.; Montemurro, C.; et al. A robust DNA isolation protocol from filtered commercial olive oil for PCR-based fingerprinting. Foods 2019, 8, 462. [Google Scholar] [CrossRef] [Green Version]
  14. International Olive Council (IOC). 2019. Available online: https://www.internationaloliveoil.org (accessed on 5 April 2021).
  15. Uylaşer, V.; Yildiz, G. The historical development and nutritional importance of olive and olive oil constituted an important part of the Mediterranean diet. Crit. Rev. Food Sci. Nutr. 2014, 54, 1092–1101. [Google Scholar] [CrossRef]
  16. Tuck, K.L.; Hayball, P.J. Major phenolic compounds in olive oil: Metabolism and health effects. J. Nutr. Biochem. 2002, 13, 636–644. [Google Scholar] [CrossRef]
  17. Harwood, J.L.; Yaqoob, P. Nutritional and health aspects of olive oil. Eur. J. Lipid Sci. Technol. 2002, 104, 685–697. [Google Scholar] [CrossRef]
  18. Stark, A.H.; Madar, Z. Olive oil as a functional food: Epidemiology and nutritional approaches. Nutr. Rev. 2002, 60, 170–176. [Google Scholar] [CrossRef] [PubMed]
  19. Piroddi, M.; Albini, A.; Fabiani, R.; Giovannelli, L.; Luceri, C.; Natella, F.; Rosignoli, P.; Rossi, T.; Taticchi, A.; Servili, M.; et al. Nutrigenomics of extra-virgin olive oil: A review. BioFactors 2017, 43, 17–41. [Google Scholar] [CrossRef]
  20. Jimenez-Lopez, C.; Carpena, M.; Lourenço-Lopes, C.; Gallardo-Gomez, M.; Lorenzo, J.M.; Barba, F.J.; Prieto, M.A.; Simal-Gandara, J. Bioactive Compounds and Quality of Extra Virgin Olive Oil. Foods 2020, 9, 1014. [Google Scholar] [CrossRef] [PubMed]
  21. Beauchamp, G.K.; Keast, R.S.; Morel, D.; Lin, J.; Pika, J.; Han, Q.; Lee, C.; Smith, A.B.; Breslin, P.A.S. Ibuprofen-like activity in extra-virgin olive oil. Nature 2005, 437, 45–46. [Google Scholar] [CrossRef] [PubMed]
  22. Servili, M.; Esposto, S.; Fabiani, R.; Urbani, S.; Taticchi, A.; Mariucci, F.; Selvaggini, R.; Montedoro, G.F. Phenolic compounds in olive oil: Antioxidant, health and organoleptic activities according to their chemical structure. Inflammopharmacology 2009, 17, 1–9. [Google Scholar] [CrossRef]
  23. Lucas, L.; Russell, A.; Russell, K. Molecular mechanisms of inflammation. Anti-inflammatory benefits of virgin olive oil and the phenolic compound oleocanthal. Curr. Pharm. Des. 2011, 17, 754–768. [Google Scholar] [CrossRef]
  24. Yubero-Serrano, E.M.; Lopez-Moreno, J.; Gomez-Delgado, F.; Lopez-Miranda, J. Extra virgin olive oil: More than a healthy fat. Eur. J. Clin. Nutr. 2019, 72, 8–17. [Google Scholar] [CrossRef] [Green Version]
  25. De la Torre-Carbot, K.; Chávez-Servín, J.L.; Jaúregui, O.; Castellote, A.I.; Lamuela-Raventós, R.M.; Nurmi, T.; Poulsen, H.E.; Gaddi, A.V.; Kaikkonen, J.; Zunft, H.F.; et al. Elevated circulating LDL phenol levels in men who consumed virgin rather than refined olive oil are associated with less oxidation of plasma LDL. J. Nutr. 2010, 140, 501–508. [Google Scholar] [CrossRef]
  26. Calabriso, N.; Massaro, M.; Scoditti, E.; D’Amore, S.; Gnoni, A.; Pellegrino, M.; Storelli, C.; De Caterina, R.; Palasciano, G.; Carluccio, M.A. Extra virgin olive oil rich in polyphenols modulates VEGF-induced angiogenic responses by preventing NADPH oxidase activity and expression. J. Nutr. Biochem. 2016, 28, 19–29. [Google Scholar] [CrossRef]
  27. Storniolo, C.E.; Rosello-Catafau, J.; Pintó, X.; Mitjavila, M.T.; Moreno, J.J. Polyphenol fraction of extra virgin olive oil protects against endothelial dysfunction induced by high glucose and free fatty acids through modulation of nitric oxide and endothelin-1. Redox Biol. 2014, 2, 971–977. [Google Scholar] [CrossRef] [Green Version]
  28. Lockyer, S.; Rowland, I.; Spencer, J.P.E.; Yaqoob, P.; Stonehouse, W. Impact of phenolic-rich olive leaf extract on blood pressure, plasma lipids and inflammatory markers: A randomised controlled trial. Eur. J. Nutr. 2017, 56, 1421–1432. [Google Scholar] [CrossRef] [Green Version]
  29. Moreno, J.J. Effect of olive oil minor components on oxidative stress and arachidonic acid mobilization and metabolism by macrophages RAW 264.7. Free Radic. Biol. Med. 2003, 35, 1073–1081. [Google Scholar] [CrossRef]
  30. Martin, M.A.; Ramos, S.; Granado-Serrano, A.B.; Rodríguez-Ramiro, I.; Trujillo, M.; Bravo, L.; Goya, L. Hydroxy-tyrosol induces antioxidant/ detoxificant enzymes and Nrf2 translocation via extracellular regulated kinases and phosphatidylinositol-3-kinase/protein kinase B pathways in HepG2 cells. Mol. Nutr. Food Res. 2010, 54, 956–966. [Google Scholar] [CrossRef] [PubMed]
  31. Giordano, E.; Davalos, A.; Visioli, F. Chronic hydroxy-tyrosol feeding modulates glutathione-mediated oxido-reduction pathways in adipose tissue: A nutrigenomic study. Nutr. Metab. Cardiovasc. Dis. 2014, 24, 1144–1150. [Google Scholar] [CrossRef] [PubMed]
  32. Orak, H.H.; Karamac, M.; Amarowicz, R.; Orak, A.; Penkacik, K. Genotype-related differences in the phenolic compound profile and antioxidant activity of extracts from olive (Olea europaea L.) leaves. Molecules 2019, 24, 1130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Caponio, F.; Gomes, T.; Pasqualone, A. Phenolic compounds in virgin olive oils: Influence of the degree of olive ripeness on organ-oleptic characteristics and shelf-life. Eur. Food Res. Technol. 2001, 212, 329–333. [Google Scholar] [CrossRef]
  34. Siliani, S.; Mattei, A.; Innocenti, L.B.; Zanoni, B. Bitter taste and phenolic compounds in extra virgin olive oil: An empirical relationship. J. Food Qual. 2006, 29, 431–441. [Google Scholar] [CrossRef]
  35. Beltrán, G.; Ruano, M.T.; Jiménez, A.; Uceda, M.; Aguilera, M.P. Evaluation of virgin olive oil bitterness by total phenol content analysis. Eur. J. Lipid Sci. Technol. 2007, 109, 193–197. [Google Scholar] [CrossRef]
  36. Pedan, V.; Popp, M.; Rohn, S.; Nyfeler, M.; Bongartz, A. Characterization of phenolic compounds and their contribution to sensory properties of olive oil. Molecules 2019, 24, 2041. [Google Scholar] [CrossRef] [Green Version]
  37. Serrilli, A.M.; Padula, G.; Bianco, A.; Petrosino, L.; Ripa, V. Searching new olive (Olea europaea L.) cultivars: Analysis of phenolic fraction in olive selections derived from a breeding program. Adv. Hort. Sci. 2008, 22, 104–109. [Google Scholar]
  38. Riachy, E.M.; Priego-Capote, F.; Rallo, L.; Luque-de Castro, M.D.; León, L. Phenolic composition of virgin olive oils from cross breeding segregating populations. Eur. J. Lipid Sci. Tech. 2012, 114, 542–551. [Google Scholar] [CrossRef]
  39. Pérez Rubio, A.G.; León, L.; Sanz, C.; Rosa, R.D.L. Fruit Phenolic Profiling: A New Selection Criterion in Olive Breeding Programs. Front. Plant Sci. 2018, 9, 241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Alagna, F.; Mariotti, R.; Panara, F.; Caporali, S.; Urbani, S.; Veneziani, G.; Esposto, S.; Taticchi, A.; Rosati, A.; Rao, R.; et al. Olive phenolic compounds: Metabolic and transcriptional profiling during fruit development. BMC Plant Biol. 2012, 12, 162. [Google Scholar] [CrossRef] [Green Version]
  41. Alagna, F.; Geu-Flores, F.; Kries, H.; Panara, F.; Baldoni, L.; O’Connor, S.E.; Osbourn, A. Identification and characterization of the iridoid synthase involved in oleuropein biosynthesis in olive (Olea europaea L.) fruits. J. Biol. Chem. 2015, 291, 5542–5554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Mougiou, N.; Trikka, F.; Trantas, E.; Ververidis, F.; Makris, A.; Argiriou, A.; Vlachonasios, K.E. Expression of hydroxy-tyrosol and oleuropein biosynthetic genes are correlated with metabolite accumulation during fruit development in olive, Olea europaea, cv Koroneiki. Plant Physiol. Biochem. 2018, 128, 41–49. [Google Scholar] [CrossRef] [PubMed]
  43. Guodong, R.; Jianguo, Z.; Xiaoxia, L.; Ying, L. Identification of putative genes for polyphenol biosynthesis in olive fruits and leaves using full-length transcriptome sequencing. Food Chem. 2019, 300, 125246. [Google Scholar] [CrossRef]
  44. Guodong, R.; Zhang, J.; Liu, X.; Li, X.; Wang, C. Combined Metabolome and Transcriptome Profiling Reveal Optimal Harvest Strategy Model Based on Different Production Purposes in Olive. Foods 2021, 10, 360. [Google Scholar] [CrossRef]
  45. Bandino, G.; Mulas, M.; Sedda, P.; Moro, C. Le varietà di Olivo Della Sardegna; Consorzio Interprovinciale per la Frutticoltura di Cagliari, Oristano e Nuoro: Cagliari, Italy, 2001; p. 256. [Google Scholar]
  46. Deiana, P.; Santona, M.; Dettori, S.; Molinu, M.G.; Dore, A.; Culeddu, N.; Azara, E.; Naziri, E.; Tsimidou, M.Z. Can all the Sardinian varieties support the PDO “Sardegna” virgin olive oil? Eur. J. Lipid Sci. Technol. 2019, 121, 1800135. [Google Scholar] [CrossRef]
  47. Lombardo, N.; Alessandrino, M.; Godino, G.; Madeo, A.; Ciliberti, A.; Pellegrino, M. Caratterizzazione della crescita, produttività, inoliazione e qualità dell’olio di 10 cultivar di olivo autoctone abruzzesi. In Proceedings of the National Congress “Maturazione e Raccolta Delle Olive: Strategie e Tecnologie Per Aumentare la Competitività in Olivicoltura”, ARSSA, Regione Abruzzo Alanno, Italy, Alanno (PE), Italy, 1 April 2006; p. 238. [Google Scholar]
  48. Albertini, E.; Torricelli, R.; Bitocchi, E.; Raggi, L.; Marconi, G.; Pollastri, L.; Di Minco, G.; Battistini, A.; Papa, R.; Veronesi, F. Structure of genetic diversity in Olea europaea L. cultivars from central Italy. Mol. Breed. 2011, 27, 533–547. [Google Scholar] [CrossRef]
  49. Squeo, G.; Difonzo, G.; Silletti, R.; Paradiso, V.M.; Summo, C.; Pasqualone, A.; Caponio, F.C. Bambina, una varietà minore pugliese: Profilo di maturazione, composizione delle drupe e caratterizzazione chimica dell’olio vergine. Riv. Ital. Sost. Grasse 2019, 94, 143–149. [Google Scholar]
  50. Squeo, G.; Silletti, R.; Mangini, G.; Summo, C.; Caponio, F. The potential of Apulian olive biodiversity: The case of Oliva Rossa Virgin Olive Oil. Foods 2021, 10, 369. [Google Scholar] [CrossRef] [PubMed]
  51. Spadoni, A.; Sion, S.; Gadaleta, S.; Savoia, M.A.; Piarulli, L.; Fanelli, V.; Di Rienzo, V.; Taranto, F.; Miazzi, M.M.; Montemurro, C.; et al. A simple and rapid method for genomic DNA extraction and microsatellite analysis in tree plants. J. Agric. Sci. Technol. 2019, 21, 1215–1226. [Google Scholar]
  52. Sefc, K.M.; Lopes, M.S.; Mendonça, D.; Rodrigues Dos Santos, M.; Laimer Da Câmara Machado, M.; Da Câmara Machado, A. Identification of microsatellite loci in olive (Olea europaea L.) and their characterization in Italian and Iberian olive trees. Mol. Ecol. 2000, 9, 1171–1193. [Google Scholar] [CrossRef]
  53. Carriero, F.; Fontanazza, G.; Cellini, F.; Giorio, G. Identification of simple sequence repeats (SSRs) in olive (Olea europaea L.). Theor. Appl. Genet. 2002, 104, 301–307. [Google Scholar] [CrossRef] [PubMed]
  54. De La Rosa, R.; James, C.M.; Tobutt, K.R. Using microsatellites for paternity testing in olive progenies. Hort. Sci. 2004, 39, 351–354. [Google Scholar] [CrossRef] [Green Version]
  55. Pasqualone, A.; Montemurro, C.; Summo, C.; Sabetta, W.; Caponio, F.; Blanco, A. Effectiveness of microsatellites DNA markers in checking the identity of PDO extra virgin olive oil. J. Agric. Food Chem. 2007, 55, 3857–3862. [Google Scholar] [CrossRef] [PubMed]
  56. Baldoni, L.; Cultrera, N.G.; Mariotti, R.; Ricciolini, C.; Arcioni, S.; Vendramin, G.G.; Buonamici, A.; Porceddu, A.; Sarri, V.; Ojeda, M.A.; et al. A consensus list of microsatellite markers for olive genotyping. Mol. Breed. 2009, 24, 213–231. [Google Scholar] [CrossRef]
  57. Boucheffa, S.; Miazzi, M.M.; di Rienzo, V.; Mangini, G.; Fanelli, V.; Tamendjari, A.; Pignone, D.; Montemurro, C. The coexistence of oleaster and traditional varieties affects genetic diversity and population structure in Algerian olive (Olea europaea L.) germplasm. Gen. Res. Crop Evol. 2017, 64, 379–390. [Google Scholar] [CrossRef]
  58. Di Rienzo, V.; Sion, S.; Taranto, F.; D’Agostino, N.; Montemurro, C.; Fanelli, V.; Sabetta, W.; Boucheffa, S.; Tamendjari, A.; Pasqualone, A.; et al. Genetic flow among olive populations within the Mediterranean basin. PeerJ 2018, 6, e5260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Caponio, F.; Catalano, P. Hammer crushers vs disk crushers: The influence of working temperature on the quality and preservation of virgin olive oil. Eur. Food Res. Technol. 2001, 213, 219–224. [Google Scholar] [CrossRef]
  60. Bendini, A.; Cerretani, L.; Carrasco-Pancorbo, A.; Gómez-Caravaca, A.M.; Segura-Carretero, A.; Fernández-Gutiérrez, A.; Lercker, G. Phenolic molecules in virgin olive oils: A survey of their sensory properties, health effects, antioxidant activity and analytical methods. An overview of the last decade Alessandra. Molecules 2007, 12, 1679–1719. [Google Scholar] [CrossRef] [PubMed]
  61. Zago, L.; Squeo, G.; Bertoncini, E.I.; Difonzo, G.; Caponio, F. Chemical and sensory characterization of Brazilian virgin olive oils. Food Res. Int. 2019, 126, 108588. [Google Scholar] [CrossRef] [PubMed]
  62. International Olive Council (IOC). Determination of Biophenols in Olive Oils by HPLC; COI/ T.20/Doc. No 29/Rev.1; International Olive Council: Madrid, Spain, 2017. [Google Scholar]
  63. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  64. Baldoni, L.; Tosti, N.; Ricciolini, C.; Belaj, A.; Arcioni, S.; Pannelli, G.; Germanà, M.A.; Mulas, M.; Porceddu, A. Genetic structure of wild and cultivated olives in the central Mediterranean basin. Ann. Bot. 2006, 98, 935–942. [Google Scholar] [CrossRef]
  65. Dίez, C.M.; Trujillo, I.; Martinez-Uriroz, N.; Barranco, D.; Rallo, L.; Marfil, P.; Gaut, B.S. Olive domestication and diversification in the Mediterranean basin. New Phytol. 2015, 206, 436–447. [Google Scholar] [CrossRef]
  66. Erre, P.; Chessa, I.; Munoz-Diez, C.; Belaj, A.; Rallo, L.; Trujillo, I. Genetic diversity and relationships between wild and cultivated olives (Olea europaea L.) in Sardinia as assessed by SSR markers. Genet. Resour. Crop Evol. 2010, 57, 41–54. [Google Scholar] [CrossRef]
  67. Covas, M.I.; Ruiz-Gutiérrez, V.; De La Torre, R.; Kafatos, A.; Lamuela-Raventós, R.M.; Osada, J.; Visioli, F. Minor components of olive oil: Evidence to date of health benefits in humans. Nutr. Rev. 2006, 64 (Suppl. 4), S20–S30. [Google Scholar] [CrossRef]
  68. The European Parliament and The Council of the European Union. Commission Regulation (EU) No 1924/2006 of 20 December 2006 on Nutrition and Health Claims Made on Foods. Off. J. Eur. Union 2006, 1924. Available online: https://eur-lex.europa.eu/eli/reg/2006/1924/oj#document1 (accessed on 5 April 2021).
  69. Commission Regulation (EU) No 432/2012 of 16 May 2012 Establishing a List of Permitted Health Claims Made on Foods, Other than Those Referring to the Reduction of Disease Risk and to Children’s Development and Health. Available online: https://eur-lex.europa.eu/eli/reg/2012/432/oj (accessed on 5 April 2021).
  70. Malik, N.S.A.; Bradford, J.M. Recovery and stability of oleuropein and other phenolic compounds during extraction and processing of olive (Olea europaea L.) leaves. J. Food Agric. Environ. 2008, 6, 8–13. [Google Scholar]
  71. Baiano, A.; Terracone, C.; Viggiani, I.; Del Nobile, M.A. Effects of cultivars and location on quality, phenolic content and antioxidant activity of extra-virgin olive oils. J. Am. Oil Chem. Soc. 2013, 90, 103–111. [Google Scholar] [CrossRef]
  72. Miho, H.; Díez, C.M.; Mena-Bravo, A.; Sánchez de Medina, V.; Moral, J.; Melliou, E.; Magiatis, P.; Rallo, L.; Barranco, D.; Priego-Capote, F. Cultivar influence on variability in olive oil phenolic profiles determined through an extensive germplasm survey. Food Chem. 2018, 266, 192–199. [Google Scholar] [CrossRef]
  73. Bruno, L.; Picardi, E.; Pacenza, M.; Chiappetta, A.; Muto, A.; Gagliardi, O.; Muzzalupo, I.; Pesole, G.; Bitonti, M.B. Changes in gene expression and metabolic profile of drupes of Olea europaea L. cv Carolea in relation to maturation stage and cultivation area. BMC Plant Biol. 2019, 19, 428. [Google Scholar] [CrossRef]
  74. Talhaoui, N.; Gómez-Caravaca, A.M.; León, L.; De la Rosa, R.; Fernández-Gutiérrez, A.; Segura-Carretero, A. Pattern of Variation of Fruit Traits and Phenol Content in Olive Fruits from Six Different Cultivars. J. Agric. Food Chem. 2015, 63, 10466–10476. [Google Scholar] [CrossRef] [PubMed]
  75. Velázquez-Palmero, D.; Romero-Segura, C.; García-Rodríguez, R.; Hernández, M.L.; Vaistij, F.E.; Graham, I.A.; Pérez, A.G.; Martínez-Rivas, J.M. An oleuropein b-glucosidase from olive fruit is involved in determining the phenolic composition of Virgin Olive Oil. Front. Plant Sci. 2017, 8, 1902. [Google Scholar] [CrossRef] [Green Version]
  76. Servili, M.; Selvaggini, R.; Esposto, S.; Taticchi, A.; Montedoro, G.; Morozzi, G. Health and sensory properties of virgin olive oil hydrophilic phenols: Agronomic and technological aspects of production that affect their occurrence in the oil. J. Chromatogr. A 2004, 1054, 113–127. [Google Scholar] [CrossRef]
  77. Paiva-Martins, F.; Fernandes, J.; Rocha, S.; Nascimento, H.; Vitorino, R.; Amado, F.; Santos-Silva, A. Effects of olive oil polyphenols on erythrocyte oxidative damage. Mol. Nutr. Food Res. 2009, 53, 609–616. [Google Scholar] [CrossRef]
  78. Artajo, L.S.; Romero, M.P.; Motilva, M.J. Transfer of phenolic compounds during olive oil extraction in relation to ripening stage of the fruit. J. Sci. Food Agric. 2006, 86, 518–527. [Google Scholar] [CrossRef]
  79. Squeo, G.; Difonzo, G.; Summo, C.; Crecchio, C.; Caponio, F. Study of the influence of technological coadjuvants on enzyme activities and phenolic and volatile compounds in virgin olive oil by a response surface methodology approach. LWT 2020, 133, 109887. [Google Scholar] [CrossRef]
  80. Cecchi, L.; Migliorini, M.; Cherubini, C.; Giusti, M.; Zanoni, B.; Innocenti, M.; Mulinacci, N. Phenolic profiles, oil amount and sugar content during olive ripening of three typical Tuscan cultivars to detect the best harvesting time for oil production. Food Res. Int. 2013, 54, 1876–1884. [Google Scholar] [CrossRef]
  81. Kalua, C.M.; Allen, M.S.; Bedgood, D.R.; Bishop, A.G.; Prenzler, P.D. Discrimination of olive oils and fruits into cultivars and maturity stages based on phenolic and volatile compounds. J. Agric. Food Chem. 2005, 53, 8054–8062. [Google Scholar] [CrossRef] [PubMed]
  82. Cirilli, M.; Caruso, G.; Gennai, C.; Urbani, S.; Frioni, E.; Ruzzi, M.; Servili, M.; Gucci, R.; Poerio, E.; Muleo, R.M. The Role of Polyphenoloxidase, Peroxidase, and β-Glucosidase in Phenolics Accumulation in Olea europaea L. Fruits under Different Water Regimes. Front. Plant Sci. 2017, 8, 717. [Google Scholar] [CrossRef] [Green Version]
  83. Ortega-Garcia, F.; Blanco, S.; Peinado, M.A.; Peragon, J. Polyphenol oxidase and its relationship with oleuropein concentration in fruits and leaves of olive (Olea europaea) cv. “Picual” trees during fruit ripening. Tree Physiol. 2008, 28, 45–54. [Google Scholar] [CrossRef] [PubMed]
  84. Morelló, J.R.; Romero, M.P.; Ramo, T.; Motilva, M.J. Evaluation of L-phenylalanine ammonia-lyase activity and phenolic profile in olive drupe (Olea europaea L.) from fruit setting period to harvesting time. Plant Sci. 2005, 168, 65–72. [Google Scholar] [CrossRef]
  85. Ryan, D.; Antolovich, M.; Herlt, T.; Prenzler, P.D.; Lavee, S.; Robards, K. Identification of phenolic compounds in tissues of the novel olive cultivar hardy’s mammoth. J. Agric. Food Chem. 2002, 50, 6716–6724. [Google Scholar] [CrossRef]
  86. Facchini, P.J.; De Luca, V. Differential and tissue-specific expression of a gene family for tyrosine/dopa decarboxylase in opium poppy. J. Biol. Chem. 1994, 269, 26684–26690. [Google Scholar] [CrossRef]
  87. Saimaru, H.; Orihara, Y. Biosynthesis of acteoside in cultured cells of Olea europaea. J. Nat. Med. 2010, 64, 139–145. [Google Scholar] [CrossRef]
Figure 1. Unweighted Neighbor Joining dendrogram illustrating the clusterization of the olive cultivars under study, and a map of Italy with indication of the considered regions and the corresponding drupe sampling sites.
Figure 1. Unweighted Neighbor Joining dendrogram illustrating the clusterization of the olive cultivars under study, and a map of Italy with indication of the considered regions and the corresponding drupe sampling sites.
Foods 10 01371 g001
Figure 2. Trend of variation in some polyphenolic compound contents in VOOs of the analyzed genotypes, during three developmental stages of olive fruits (T2, T3, T4). (A) Tyrosol; (B) hydroxy-tyrosol; (C) oleacein (3,4-DHPEA-EDA); (D) total phenol content (TPC).
Figure 2. Trend of variation in some polyphenolic compound contents in VOOs of the analyzed genotypes, during three developmental stages of olive fruits (T2, T3, T4). (A) Tyrosol; (B) hydroxy-tyrosol; (C) oleacein (3,4-DHPEA-EDA); (D) total phenol content (TPC).
Foods 10 01371 g002
Figure 3. Expression profiles of the analyzed genes involved in the first step of the phenol biosynthetic pathway during fruit development (T1 = yellow-green drupes, T2 = turning drupes, T3 = almost dark drupes). The mRNA levels were determined by qRT-PCR and relatively expressed as ΔΔCt. (A) TYRD = tyrosine/dopa decarboxylase; (B) CuAO = copper amine oxidase; (C) ALDH = alcohol dehydrogenase. Data are means (n = 3) ± standard deviation. Different letters indicate significant differences according to one-way ANOVA followed by Fisher’s LSD post-hoc test (p = 0.05).
Figure 3. Expression profiles of the analyzed genes involved in the first step of the phenol biosynthetic pathway during fruit development (T1 = yellow-green drupes, T2 = turning drupes, T3 = almost dark drupes). The mRNA levels were determined by qRT-PCR and relatively expressed as ΔΔCt. (A) TYRD = tyrosine/dopa decarboxylase; (B) CuAO = copper amine oxidase; (C) ALDH = alcohol dehydrogenase. Data are means (n = 3) ± standard deviation. Different letters indicate significant differences according to one-way ANOVA followed by Fisher’s LSD post-hoc test (p = 0.05).
Foods 10 01371 g003
Figure 4. Simplified representation of the biosynthetic pathway of some important olive polyphenols. The previously identified enzymes and reactions (A) as well as the recently identified ones, with additional steps in the biosynthesis (B), are reported. Abbreviations: 3,4-DHPEA-EDA: elenolic acid linked to 3,4-dihydroxyphenyl ethanol (3,4-DHPEA); 4-HPA: 4-hydroxyphenylacetic acid; 3,4-DHPA: 3,4-dihydroxyphthalic acid; TYRD: tyrosine/DOPA decarboxylase; PPO: polyphenol oxidase; CuAO: copper amine oxidase; ALDH: alcohol dehydrogenase; DCC: DOPA decarboxylase; PAR: phenylacetaldehyde reductase. Adapted from [40,44].
Figure 4. Simplified representation of the biosynthetic pathway of some important olive polyphenols. The previously identified enzymes and reactions (A) as well as the recently identified ones, with additional steps in the biosynthesis (B), are reported. Abbreviations: 3,4-DHPEA-EDA: elenolic acid linked to 3,4-dihydroxyphenyl ethanol (3,4-DHPEA); 4-HPA: 4-hydroxyphenylacetic acid; 3,4-DHPA: 3,4-dihydroxyphthalic acid; TYRD: tyrosine/DOPA decarboxylase; PPO: polyphenol oxidase; CuAO: copper amine oxidase; ALDH: alcohol dehydrogenase; DCC: DOPA decarboxylase; PAR: phenylacetaldehyde reductase. Adapted from [40,44].
Foods 10 01371 g004
Table 1. Overview of the four drupe-collecting stages and type of analysis carried out in this study.
Table 1. Overview of the four drupe-collecting stages and type of analysis carried out in this study.
Sampling Time-PointDrupe StatusPerformed Analysis
T1yellow-greenGene expression
T2turningGene expression + oil biochemical characterization
T3almost darkGene expression + oil biochemical characterization
T4fully darkOil biochemical characterization
Table 2. Molecular fingerprinting of the analyzed cultivars by means of 12 microsatellite markers (SA, AB and AP, respectively, indicate Sardinian, Abruzzo and Apulian varieties). Allele sizes are reported in base pair (bp).
Table 2. Molecular fingerprinting of the analyzed cultivars by means of 12 microsatellite markers (SA, AB and AP, respectively, indicate Sardinian, Abruzzo and Apulian varieties). Allele sizes are reported in base pair (bp).
CultivarMicrosatellite Marker
DCA03DCA05DCA09DCA13DCA15DCA17DCA18GAPU45GAPU71bGAPU101EMO90EMOL
Corsicana (SA)232253198208162182120124246257109115177183181181127144192200186188190198
Semidana (SA)245253202206162174122140246257117117177179183185130144192218186188190198
Sivigliana (SA)243253206208162174120120246257113117179183181181127144198200188194190198
Tortiglione (AB)229245198206186186120120246246113113177185195195124127192192188198198198
Dritta (AB)229229208212186194120120246257113113177177181185127127190206186190192192
Gentile dell’Aquila (AB)243249206206162162120120246246--179181183183144144200206188194198198
Bambina (AP)243253198206--120120246266117117177177183185127144--188188198198
Cima di Melfi (AP)243253198206184184120120246246143143177181181185124144182182190190192198
Oliva Rossa (AP)239253206206184194120122246273143143177181185185127144182218188188200214
“-“ indicates missing data.
Table 3. Phenolic compounds and total phenolic content (TPC) detected in VOOs of all the analyzed cultivars from fruits at different ripening stages (T2 = turning drupes, T3 = almost dark drupes, T4 = fully dark drupes).
Table 3. Phenolic compounds and total phenolic content (TPC) detected in VOOs of all the analyzed cultivars from fruits at different ripening stages (T2 = turning drupes, T3 = almost dark drupes, T4 = fully dark drupes).
RegionGenotypeDrupe Collecting Time-PointHydroxy-Tyrosol (3,4-DHPEA)Tyrosol (p-HPEA)Oleacein (3,4-DHPEA-EDA)TPC
SardiniaCorsicanaT20.14 ± 0.12 d0.32 ± 0.02 f33.95 ± 0.21 a572 ± 3 a
CorsicanaT30.27 ± 0.01 d0.36 ± 0.02 f30.06 ± 0.17 b558 ± 1 a
CorsicanaT40.13 ± 0.01 d0.43 ± 0.01 ef9.20 ± 0.43 e511 ± 1 b
SemidanaT20.81 ± 0.28 bc1.18 ± 0.17 b13.32 ± 3.95 d479 ± 23 de
SemidanaT30.74 ± 0.19 c1.32 ± 0.08 b9.69 ± 0.61 e469 ± 9 ef
SemidanaT40.91 ± 0.15 abc1.55 ± 0.12 a7.78 ± 0.67 e453 ± 3 f
SiviglianaT21.07 ± 0.20 ab0.65 ± 0.10 d20.83 ± 1.46 c503 ± 10 bc
SiviglianaT30.70 ± 0.15 c0.55 ± 0.04 de7.88 ± 0.55 e487 ± 6 cd
SiviglianaT41.14 ± 0.29 a0.83 ± 0.04 c12.72 ± 1.46 d477 ± 10 de
ApuliaBambinaT20.40 ± 0.06 a1.51 ± 0.06 c4.85 ± 0.52 b358 ± 15 b
BambinaT30.36 ± 0.06 a1.86 ± 0.07 b10.75 ± 1.76 a406 ± 10 a
BambinaT40.32 ± 0.09 a2.35 ± 0.13 a3.12 ± 1.83 b392 ± 6 a
Cima di MelfiT20.75 ± 0.01 a3.01 ± 0.01 a17.49 ± 1.89 c386 ± 5 c
Cima di MelfiT30.39 ± 0.06 b1.36 ± 0.01 b36.79 ± 2.76 b502 ± 20 b
Cima di MelfiT40.41 ± 0.05 b1.00 ± 0.05 c65.52 ± 12.26 a730 ± 11 a
Oliva RossaT20.30 ± 0.05 ab0.48 ± 0.02 b55.34 ± 2.07 a677 ± 47 a
Oliva RossaT30.35 ± 0.05 a0.45 ± 0.01 b21.36 ± 3.46 b478 ± 5 b
Oliva RossaT40.22 ± 0.06 b0.59 ± 0.04 a15.46 ± 1.77 c416 ± 19 c
AbruzzoDrittaT20.12 ± 0.02 c0.25 ± 0.01 b3.71 ± 0.24 de129 ± 27 d
DrittaT30.03 ± 0.01 c0.23 ± 0.02 b0.50 ± 0.09 e39 ± 2 f
DrittaT40.09 ± 0.03 c0.26 ± 0.04 b1.41 ± 0.08 e71 ± 9 e
Gentile dell’AquilaT20.05 ± 0.01 c0.17 ± 0.03 d6.19 ± 1.37 d123 ± 38 d
Gentile dell’AquilaT30.07 ± 0.02 c0.13 ± 0.05 cd1.09 ± 0.51 e106 ± 26 de
Gentile dell’AquilaT40.09 ± 0.01 c0.16 ± 0.02 cd16.80 ± 1.35 c176 ± 15 c
TortiglioneT20.59 ± 0.12 a0.38 ± 0.07 a42.37 ± 3.78 a804 ± 68 a
TortiglioneT30.29 ± 0.08 b0.21 ± 0.02 b24.44 ± 3.72 b377 ± 30 b
TortiglioneT40.28 ± 0.06 b0.34 ± 0.01 a18.34 ± 1.61 c395 ± 6 b
Data are means (n = 3) ± standard deviation. Different letters for each parameter and region indicate significant differences according to one-way ANOVA followed by Fisher’s LSD post-hoc test (p = 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sabetta, W.; Mascio, I.; Squeo, G.; Gadaleta, S.; Flamminii, F.; Conte, P.; Di Mattia, C.D.; Piga, A.; Caponio, F.; Montemurro, C. Bioactive Potential of Minor Italian Olive Genotypes from Apulia, Sardinia and Abruzzo. Foods 2021, 10, 1371. https://doi.org/10.3390/foods10061371

AMA Style

Sabetta W, Mascio I, Squeo G, Gadaleta S, Flamminii F, Conte P, Di Mattia CD, Piga A, Caponio F, Montemurro C. Bioactive Potential of Minor Italian Olive Genotypes from Apulia, Sardinia and Abruzzo. Foods. 2021; 10(6):1371. https://doi.org/10.3390/foods10061371

Chicago/Turabian Style

Sabetta, Wilma, Isabella Mascio, Giacomo Squeo, Susanna Gadaleta, Federica Flamminii, Paola Conte, Carla Daniela Di Mattia, Antonio Piga, Francesco Caponio, and Cinzia Montemurro. 2021. "Bioactive Potential of Minor Italian Olive Genotypes from Apulia, Sardinia and Abruzzo" Foods 10, no. 6: 1371. https://doi.org/10.3390/foods10061371

APA Style

Sabetta, W., Mascio, I., Squeo, G., Gadaleta, S., Flamminii, F., Conte, P., Di Mattia, C. D., Piga, A., Caponio, F., & Montemurro, C. (2021). Bioactive Potential of Minor Italian Olive Genotypes from Apulia, Sardinia and Abruzzo. Foods, 10(6), 1371. https://doi.org/10.3390/foods10061371

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop