Abstract

Rainfed olive groves have been traditionally tilled in order to reduce the competition for water and nutrients. In sloping Mediterranean olive groves, this practice leads to high erosion rates, resulting in a reduction in soil fertility. Cover crops have been employed as a sustainable olive grove management strategy, but previous studies found differences in their effect on fruit load and there is scarce information on their influence on the virgin olive oil (VOO) quality. The aim of this study is to evaluate the effect of different soil management strategies on olive and oil yield and VOO physicochemical and sensory characteristics in a rainfed olive grove (238 trees·ha−1) of Cornicabra cultivar, the main in Central Spain. No effect of soil management was found in olive or oil yield along three cropping seasons. VOO quality was mainly influenced by the year, but slight differences were found in the driest year between the studied treatments. Small differences in fruit ripening, pigments, and several fatty acids as well as in sensory parameters were found.

1. Introduction

Olive tree (Olea europaea L.) is one of the most drought-tolerant tree crops in the world [1]. They are able to tolerate low availability of soil water by means of morphological and physiological adaptation [2, 3]. Three periods of high sensitivity to water stress can be considered [4]: flowering; fruit growth, 6 to 10 weeks after blooming; and oil accumulation, 18 to 22 weeks after blooming. Water stress in the first period reduces the number of fruits, while water stress in the other two periods reduces fruit fresh weight. A positive relationship between water stress of olive trees and phenolic compounds [5, 6], volatile profile [7, 8], monounsaturated fatty acid levels, and sensory properties [9] of olive oil has been described.

In the Mediterranean basin, rainfall and underground water resources are the only supplies for most of the olive trees [10]. In this area, traditional olive groves were designed with low plant density according to total rainfall and paying attention to canopy size [3]. Actually around 50% out of 2,500,000 ha of olive groves in Spain is tilled [11] to avoid weed competition for water and nutrients, in order to increase olive tree yield. This conventional practice gives rise to a large area of bare soil that is prone to erosion processes, one of the most important land degradation driving processes in Mediterranean areas [12, 13].

Cover crops in olive groves have proved to be an effective tool in reducing soil erosion [1416] and improving water storage and physical properties [17], but its effect on olive tree yield is controversial. A negative effect on fruit yield due to cover crops was reported by Gucci et al. [18] in a high-density olive grove and by Caruso et al. [19] in a young and intensive olive grove. Nevertheless, other authors [17, 2022] did not find a significant reduction in olive yield in mature olive groves.

The influence of cover crops on food quality has been studied in vineyards [2327], apple orchard [28], eggplants [29], and sweet cherries [30]. In olive groves, the influence of soil management on some oil quality parameters has been seldom studied [18, 19, 21, 31] and to our knowledge there are no previous studies on the effect of cover crops on sensory characteristics of virgin olive oil.

Cornicabra cultivar is the most common in Central Spain, being the second cultivar in extension in Iberian Peninsula [32] covering over 270,000 ha [33]. Although Cornicabra cv. is considered one of the most drought-tolerant cultivars in Spain [33], small amounts of water increase olive and oil yield [34, 35]. The oil of Cornicabra is valued for its high stability and food sensory characteristics [36], with medium bitterness and pungency mainly due to the content in total phenols [37] and high oleic acid and low linoleic acid content [38]. Phenolic compounds, oleic acid, and monounsaturated fatty acids (MUFA) of VOO have beneficial effects on human health [39], and thus a higher content in total phenols and a high ratio of MUFA to polyunsaturated fatty acids (PUFA) are desirable in VOOs.

In order to recommend a more sustainable practice to rainfed olive growers of Central Spain, we performed a comparative trial in an intensive olive grove of Cornicabra cv. to evaluate the effect of three different cover crops and minimum tillage on olive yield and physicochemical and sensory virgin olive oil parameters.

2. Materials and Methods

2.1. Field Conditions, Climatic Data, and Plant Material

This study was performed in an experimental olive grove (Olea europaea L.) located in Central Spain, in southern Madrid (40°4′21′′N, 3°31′11′′W). The average elevation is 540 m.a.s.l., and the slope ranges from 9 to 12%. The soil is classified as Haplic Gypsisol [40], with a xeric moisture regime. The field capacity (FC) was 0.21 m3·m−3 and permanent wilting point (PWP) was 0.08 m3·m−3. This soil has a moderate silt content (29%) and it is low in organic matter (1.1%).

The climate is semiarid Mediterranean, with long hot summers (quite often above 35°C in July and August) and cold winters (2.7°C on average in December). The average annual temperature is 13.6°C and reference evapotranspiration (ET0) is 1112 mm. The annual precipitation is approximately 390 mm with high inter- and intra-annual variability [41, 42]. During the last 17 years, the total precipitation has strongly decreased (40%) in comparison with a previous period (1933–1969) in the study area [43]. Since it is a rainfed olive grove, as most of the region, the amount and distribution of rainfall and ET0 are very important. Thus, there is an automatic weather station in the olive grove to record the temperature, rainfall amount, relative humidity, solar radiation, wind speed and direction, and atmospheric pressure every 10 minutes. With the previous data, Penman-Monteith ET0 was calculated according to FAO [44].

The olive plantation was established in 2004 with trees in 6 × 7 m spacings (238 trees·ha−1), with an area of approximately 3 ha.

2.2. Soil Management and Treatments

Before the start of the trial, in November 2010, the whole area was tilled with a chisel at 0.30 m depth following the traditional management. The tillage was performed in 6 m wide interrows. Each treatment was performed in 3 consecutive interrows. The treatments consisted of the following: (1) barley annual cover (Hordeum vulgare L.) that was seeded each autumn at 70 kg·ha−1; (2) sainfoin (Onobrychis viciifolia Scop.), a legume that was seeded each year (42 kg·ha−1); (3) purple false brome (Brachypodium distachyon L. P. Beauv.), henceforth Brachypodium, a permanent grass cover seeded once (the first year) during the study period at 40 kg·ha−1; and (4) control, consisting of one pass per year with a chisel at 0.15–0.20 m deep in mid-November. All of the treatments were mechanically mowed once in the spring (during the first fortnight of May), except in 2013, when the vegetation was mowed twice (the second cut was at the end of May) because of high vegetation growth resulting from the abundant rains of that spring. Plant debris was left on the surface.

2.3. Soil Moisture Measurements

Three EC-5 soil moisture sensors (Decagon Devices Inc.) per treatment were placed in the interrow at 30 cm depth plugged into a data logger, which recorded data each 15 minutes. Data were fortnightly downloaded to a laptop.

2.4. Olive Harvesting, Yield and Oil Elaboration, and Oil Content Determination

Fruit olives were harvested in three consecutive cropping seasons: 2011/2012, 2012/2013, and 2013/2014. Ten healthy olive trees per treatment were randomly selected each year and fruits were hand-picked, when the maturity index (MI) was around 3.5. This MI was selected in order to obtain an optimal oil content and quality [45]. Yield per tree was measured in field. Olives collected of the same treatment were mixed up and 3 subsamples of 5 kg (3 subsamples × 4 treatments) were carried to the laboratory. MI was determined at the laboratory for each subsample of fruits following the method of Beltrán et al. [45] based on skin and pulp fruit colors. Olive oils from each subsample were extracted in the last two years using the Abencor system (MC2 Ingenierias y Sistemas), determining the industrial oil yield [46]. Oil yield per hectare was obtained multiplying fruit yield per tree by industrial oil yield by the number of trees per hectare.

2.5. Analysis of Virgin Olive Oil

In 2012 and 2013, different VOO parameters (2 years × 4 treatments × 3 repetitions) were analyzed.

2.5.1. Physicochemical Quality Parameters

Free acidity is expressed as percentage of oleic acid; peroxide value (PV) is expressed as milliequivalents of active oxygen per kilogram of oil (meq O2/kg); and UV spectrophotometric indices (, , and extinction coefficients) were determined following the analytical methods described in the European Commission Regulation 2568/91 and later amendments [47].

2.5.2. Total Phenolic Content

Total phenolic compounds were determined after methanol extraction and subsequent reaction with Folin-Ciocalteu reagent and measured at a wavelength of 725 nm [48]. Results are expressed as mg of caffeic acid equivalents per kg of oil.

2.5.3. Pigments (Chlorophylls and Carotenoids) and Chromatic Coordinates

Carotenoids and chlorophylls were determined at a wavelength of 470 and 670 nm, respectively, in cyclohexane, according to the method of Minguez-Mosquera et al. [49]. The concentrations of chlorophyll and carotenoids were expressed as mg of pheophytin and lutein per kg, respectively. Chromatic coordinates were measured by the software CINTRAL to obtain the color according to the CIEL    method [50].

2.5.4. Fatty Acid Composition

The fatty acid composition of the oils was determined by gas chromatography as fatty acid methyl esters (FAMEs). FAMEs were prepared by saponification/methylation with sodium methylate according to European Regulations (EEC 2568/91) [47]. A chromatographic analysis was performed in an Agilent Technologies 6850 series II Network GC system gas chromatograph equipped with a 60 m × 0.25 mm × 0.20 μm film thickness fused capillary column Supelco 24111 (Agilent Technologies) coupled to a flame ionization detector. Samples were introduced into the column at 170°C during a period of 30 min; after this time, the temperature was increased by 5°C/min to 200°C and maintained for 12 min. The flow rate of He, used as carrier gas, was 0.5 mL/min. Injector and flame ionization detector temperatures were 230°C and 250°C, respectively. FAMEs were identified by comparing their retention times with those of standard compounds.

Fatty acid composition was calculated as the percentage of total fatty acids, after their conversion to methyl esters, according to regulation EEC 2568/91 [47].

2.5.5. Sensory Analysis

Sensory analysis was performed by the “Panel de Catadores de Aceite de Oliva Virgen de la Comunidad de Madrid,” according to the method described in the European Commission Regulation (EC) 640/2008 [51]. This method allows for the classification of VOOs based on the detection of negative attributes (fusty/muddy, musty, winey, rancid, wet wood, and others) as well as the measurement of the intensity of three positive attributes (green or ripe fruitiness, bitterness, and pungency). The panel was constituted at least by 8 trained tasters that scored the descriptors on a normalized sheet (from 0 to 10).

Chemical determinations were performed twice per sample, obtaining the mean for statistical analysis.

2.6. Statistical Analysis

Statistical analyses were performed with the software package SPSS 19 (SPSS Inc.) for Windows. When needed, data were log-transformed before parametric testing. Repeated-measures ANOVA was performed for yield data and MAN(C)OVA analyses were performed on the VOO data to assess the differences between treatments and years, employing MI as covariate. LSD test () was used to establish significant differences between groups.

After MANOVA a Discriminant Function Analysis was performed to assess how the dependent variables discriminate groups.

3. Results and Discussion

3.1. Meteorology and Soil Water Content

Year 2012 registered the lowest total rainfall amount in the studied period (237 mm) while 2011 and 2013 were quite similar, 292 and 300 mm, respectively. Distribution of the rainfall varied among years: autumn was the rainiest season in 2011 and 2012, while in 2013 spring was the rainiest (Figure 1). Reference evapotranspiration (ET0) values were similar among the three years: 1209, 1217, and 1218 in 2011, 2012, and 2013, respectively. Summers have been characterized by high ET0 and scarce rainfalls.

Control was the treatment with the lower soil moisture at 30 cm depth along the three years. During the rainfalls of spring in 2011 and 2013, these differences decreased; the same happened with the autumn rainfall in 2012 (Figure 2). In summer, all the treatments endured several days below the PWP (0.08 m3·m−3), mainly in 2011, which implied water stress to the olive trees in the phase of oil accumulation [4].

3.2. Fruit and Oil Yield

Olive harvest dates were quite similar, 1 December 2011, 11 December 2012, and 2 December 2013. MI was around 3.5 in 2011 and 4 in 2012. In 2013, we had to harvest with a low MI, about 3, because of frost alert at the beginning of December. Despite the low soil moisture measured during all the studied period in all the treatments, thanks to the olive tree high capacity to cope with water stress [52], especially in Cornicabra cv. [33], olive trees under semiarid climate without irrigation have low fruit production on a year basis. Mean and standard deviation of fruit yield per tree and oil yield per hectare are shown in Figures 3(a) and 3(b). 2013 was the year with the highest yield in both parameters (9.1 kg·tree−1 of fruit and 476 kg·ha−1 of oil on mean), whereas 2012 was the year with the lowest yield (1.6 kg·tree−1 of fruit and 90 kg·ha−1 of oil on mean) being considered as “off” year (low or no yield). Fruit and oil yields were affected by year but not by treatment; no interaction was detected (Table 1). There were statistically significant differences between the three crop seasons in both fruit and oil yield (2012 < 2011 < 2013) according to LSD test. This alternate bearing is typical of olive tree, mainly in Cornicabra cultivar [32], with high yields in the “on” years (heavy yield) followed by “off” years of low yields.

This absence of differences in fruit and oil yield among different soil management systems was also observed by other researchers between tilled and nontilled [17, 20, 22, 53] and rainfed and irrigated olive groves in Cornicabra cv. [35]. By contrast, Moriana et al. [54] found a reduction in Picual cv. under water stress and Caruso et al. [19] reported over 40% olive yield reduction in young trees of Frantoio cv. managed with cover crops regarding tilling management. Corleto and Cazzato [21] found a reduction in fruit yield with perennial species as cover crops but not with annual legume species regarding conventional tillage on Coratina cv. Therefore, proper management of cover crops is crucial to avoid fruit yield reductions, mainly species selection and time of mowing.

3.3. VOO: Physicochemical Parameters

All the VOO samples had their analytical parameters within the ranges established for the highest quality category “Extra Virgin Olive Oil” in the Regulation 1348/2013 [55]. These thresholds are 0.8% for acidity, 20 meq O2·kg−1 for peroxides, 2.50 for , 0.22 for , and 0.01 for .

MANOVA results are shown in Table 2 for the physicochemical determinations in the VOOs of 2012 and 2013. MI was included as a covariate due the important effect of this index over several VOO parameters [5, 35, 56] but in our study it was not statistically significant on either of the measured variables (data not shown), so it was excluded as a covariate and included as a dependent variable in the analysis. Treatment and year influenced significantly most of the measured parameters (Table 2). There was no statistically significant effect of treatment or year in peroxide index or , while and were only affected by the year.

Table 3 shows VOO parameters as a function of the treatments and years. In relation to the effect of the different treatments on the fruit, legume coverage led to earlier olive maturation (3.7) with barley resulting in the opposite effect (3.3). Control and Brachypodium VOOs had slightly higher values of free acidity (0.22 and 0.21%, resp.), but these differences were not due to the maturity of the VOOs considering that both treatments had intermediate maturity (3.5). Regarding the pigments, legume VOOs had lower values for carotenoid, chlorophyll, and total pigments because of a higher MI similar to observations by other authors [57, 58]; the highest chlorophyll to carotenoid ratio was found in barley (1.3) while control had the lowest (1.1). As control had the same MI as Brachypodium, the differences in pigment content should be explained through changes in the plant-soil relationship produced by the soil management. VOOs obtained from legume coverage had a significantly higher value of color lightness () and lower value of greenish color () than the other treatments, which is due to the lower content in carotenoid and chlorophyll pigments [59].

Comparing the results by year, MI was lower in 2013 than in 2012 despite an earlier ripening date in 2013 due to a higher fruit load of the trees, similar to the results of Dag et al. [60]. In 2012, VOOs had slightly higher acidity (0.26%) than in 2013 (0.15%) probably due to the higher MI. UV absorbance at 270 and 232 nm indicates the presence of oxidized compounds; some authors found that and increased with water stress [10, 35] as in the present study. The maximum level of total phenols accumulation is reached at veraison [61] that is closely related to the MI of VOOs in 2013, although the higher content in total phenols was in 2012. This contradiction could be related to the scarce precipitation of that year, as a linear relationship has been described between total phenols and water stress by some authors [9, 52]. Nevertheless, Palese et al. [62] stated that the differences in total phenols, , and other quality parameters between rainfed and irrigated olive trees were due to the different plant crop load and fruit ripening pattern instead of water availability. VOOs of 2012 had a significantly higher value of color lightness and lower greenish color, with less content in carotenoids, chlorophyll, and total pigments and lower chlorophyll to carotenoid ratio than in 2013, which is related to a higher MI [57, 58].

Figure 4 shows the mean and standard deviation per treatment and year for those parameters with significant interaction. Free acidity (Figure 4(a)) was higher in control in 2012 regarding the other treatments, while in 2013 this place corresponded to Brachypodium. These changes were not correlated with MI but they could be related to soil water content (Figure 2) and other changes in plant-soil relationship. Barley had a different behavior regarding total phenols (Figure 4(b)) in both years, the highest content in 2012 and the lowest in 2013. In 2013, no effect of treatment was observed in chlorophyll to carotenoid ratio (Figure 4(c)) and (Figure 4(e)) parameters, while important differences between treatments appeared in 2012. These differences in total phenols and pigments may be responding to a higher water stress of barley on the driest year over the olive trees as was pointed out by Hernández et al. [63]. Regarding parameter (Figure 4(d)), control and barley were more similar to legume in 2012 and to Brachypodium in 2013. These differences among treatments should be explained through the effect of cover crops on soil moisture dynamic along the crop season and available nutrients for olive trees, as was described in vineyards in which legume cover crops increased soil N availability and vine N uptake [24] affected must quality.

3.4. VOO: Fatty Acid Composition

Most of the fatty acids were mainly influenced by year (Table 4) and not by soil management, as is the case of oleic acid, for instance, which remained stable across treatments. Similar results were also reported by other authors in relation to water stress [35, 37] but oleic acid increases throughout fruit ripening involving a reduction in palmitic acid [35]. On mean oleic acid percentage (79.5%) was similar to Cornicabra olive oils described by other authors [36, 37], but with a lower percentage in linoleic acid (3.1%) similar to observations by Gómez-Rico et al. [7]. Margaric, stearic, linoleic, and linolenic acids were influenced by treatment besides the year. MUFA to PUFA ratio and oleic to linoleic acid ratio were influenced by treatment and year and by their interaction. Although these statistically significant differences were observed (Table 4), their magnitude was very low (Table 5).

Brachypodium and barley treatments give rise to oils with higher content of linoleic (3.18%) and linolenic (0.6%) acids, while control had higher content of margaric and stearic acids, 0.05 and 3.09%, respectively (Table 5). The highest rates in MUFA to PUFA and oleic to linoleic acid corresponded to legume, 23 and 27%, respectively, while Brachypodium had the lowest, 21 and 25%, respectively. Linoleic acid increases during ripening while stearic [64] and palmitic acids and MUFA to PUFA and oleic to linoleic acid ratios decrease [57, 65]. Nevertheless, these results do not match the MI of the different treatments: barley VOOs had lower MI, followed by control and Brachypodium and the most mature legume VOOs. Therefore, barley should have the highest stearic acid percentage and in fact it had the lowest value (2.85%), although this difference with the other VOOs was very low. In the same way, legume VOO should have lower MUFA to PUFA and oleic to linoleic acid ratios, but it was the treatment with the highest value. Gómez-Rico et al. [35] found that MUFA to PUFA ratio was higher in oils obtained under rainfed conditions. Thus, legume could be exerting higher water stress on the olive trees, although it kept more soil water content than control (Figure 2), or legume cover could have a positive effect on the nutrients of the soil [24, 29, 66] influencing VOO fatty acids.

The most mature VOOs of 2012 had higher content of linoleic acid [64] and lower content of palmitic acid [35]. Stearic acid does not accumulate during the ripening process [64], which is the reason why the content in stearic acid in 2012, the year with the higher MI, was significantly lower than in 2013. MUFA to PUFA ratio was expected to be higher in water stressed olive trees [35] but this fact was not observed between the driest 2012 and 2013 probably due to the predominant effect of maturity in olive fruits, as was pointed out by Dag et al. [60].

Figure 5 shows mean and standard deviation for those parameters with significant interaction. Brachypodium VOOs had the highest palmitoleic acid content in 2012 (Figure 5(a)), and control VOOs of 2013 had the highest percentage of stearic acid (Figure 5(b)). The higher value of oleic to linoleic acid (Figure 5(c)) and MUFA to PUFA (Figure 5(d)) ratios was reached in 2013 by legume. The differences in parameters that were not explained by MI or soil water content could be related to changes in plant-soil relationships.

3.5. VOO: Sensory Analysis

As in previous parameters, year had a higher influence on bitterness and pungency than soil management (Table 6). Barley treatment produced slightly less bitter and pungent VOOs, over 1 point below the other VOOs. Fruitiness was not affected by treatments. In relation to year effect, in 2013, VOOs had around 3 points less in bitterness and pungency, but with similar fruitiness despite a lower MI in this year. These “sweeter” 2013 VOOs can be explained by the higher rainfall amount that year which produced a reduction in total phenols content (Table 3) as has been observed by other authors [5, 62].

MANOVA results are shown in Table 7. Bitterness and pungency were influenced by treatment, by year, and by their interaction (Table 7). The interaction of treatment and year also was statistically significant for fruitiness.

Figure 6 shows mean and standard deviation of fruitiness (Figure 6(a)), bitterness (Figure 6(b)), and pungency (Figure 6(c)) of the oils for the different treatments and years. Control treatment almost did not change fruitiness, bitterness, or pungency between years, opposite to cover crops VOOs that diminished 2 points in fruitiness in 2013 regarding the previous year. Brachypodium reduced less than 1 point of the bitterness and pungency from 2012 to 2013, while legume slightly increased both parameters. In 2012, Brachypodium VOOs stood out among the other VOOs for higher bitterness and pungency, while in 2013 barley stood out for less bitter and pungent VOOs. Despite these statistically significant differences, their magnitude was not high enough to state that the VOOs were sensorially different.

3.6. Global Treatment Effect

The MANOVA was followed up with discriminant analysis employing all the variables. As previously mentioned, the main influence of the year on the VOO quality was due to a rainfall amount and crop load higher in 2013/2014 crop season. Due to this important effect of the year on the VOOs quality, discriminant analysis was conducted separately in each year (Figure 7). The most useful variables from all the measured ones to discriminate treatments were MI, free acidity, total phenols, carotenoid and chlorophyll pigments, and color, MUFA to PUFA ratio, and pungency.

The differences in VOOs between treatments were greater in the driest year, 2012, as can be seen in the longer scale in both functions of this year (±40 in -axis and ±15 in -axis) regarding 2013 (±6 in -axis and ±6 in -axis). In 2012 (Figure 7(a)), the first discriminant functions explained 93.3% of the variance, canonical , with a cumulative variance of 99.9% with the second, canonical . The first axis separates control from barley VOOs. Control VOOs were characterized by high MI and free acidity, while barley VOOs were characterized by low MI and higher content of total phenols and carotenoid and chlorophyll pigments, but a low MUFA to PUFA ratio. The second axis differentiates Brachypodium with the greenish VOOs (low ) from legume VOOs with higher MI and less content of carotenoid and chlorophyll pigments.

In 2013 (Figure 7(b)), the first discriminant functions explained 57.2% of the variance, canonical , with a cumulative variance of 88.4% with the second, canonical . Legume and control VOOs were quite similar, characterized by a higher content in total phenols and MUFA to PUFA ratio. On the other hand, Brachypodium and barley VOOs had higher content of carotenoid and chlorophyll pigments, but Brachypodium had more free acidity, greenish color, and lower MUFA to PUFA ratio; and barley VOOs were mainly characterized by low content of total phenols and “sweeter” oils.

The influence of water stress, fruit ripeness, crop load, and nutrient availability and their interactions could explain changes in VOO parameters between soil management treatments [7, 10, 60, 62]. High contents of phenolic compounds, MUFA, and oleic acid have beneficial properties on health [39], but, despite the statistically significant differences found in VOOs between treatments or years, they were too low to get involved in changes in the beneficial effect on health.

4. Conclusions

Our results support the adoption of cover crops as sustainable management of rainfed olive groves in Central Spain in the drought-tolerant cultivar Cornicabra. Cover crops have not reduced fruit or oil yield either in heavy or in low yield years. Thus, the economic profits of land farmers are maintained, while improving environmental benefits with reducing soil erosion and increasing soil organic carbon.

Virgin olive oil quality was mainly influenced by the year, but slight differences appeared in some physicochemical and sensory characteristics assigned to soil management. In the year with higher rainfall amount, with more soil water content in autumn and high fruit load, there were fewer differences between the treatments of cover crops and minimum tillage. Therefore, slight differences were found in VOOs in the driest year, which could be caused by differences in soil water availability and nutrient uptake by olive tree. More research is required to investigate the influence of cover crops on water and nutrients uptake in Mediterranean semiarid olive groves and their influence on the final product. To overcome the year-by-year differences due to the alternate bearing cycle of olive trees and annual variations of Mediterranean climate, more cropping seasons are needed to find out whether the trends seen in this study will consolidate in the future.

Competing Interests

The authors declare that they have no competing interests regarding the publication of this paper.

Acknowledgments

This work was supported by the Madrid Institute for Rural, Agricultural and Food Research and Development (IMIDRA) under Grant FP12-CVO and S2013/ABI-2717 AGRISOST Projects.