Use of RGB Vegetation Indexes in Assessing Early Effects of Verticillium Wilt of Olive in Asymptomatic Plants in High and Low Fertility Scenarios

: Verticillium Wilt of Olive, a disease caused by the hemibiotrophic vascular fungus Verticillium dahliae Kleb. presents one of the most important constraints to olive production in the world, with an especially notable impact in Mediterranean agriculture. This study evaluates the use of RGB vegetation indexes in assessing the effects of this disease during the biotrophic phase of host-pathogen interaction, in which symptoms of wilt are not yet evident. While no differences were detected by measuring stomatal conductance and chlorophyll ﬂuorescence, results obtained from RGB indexes showed signiﬁcant differences between control and inoculated plants for indexes Saturation, a*, b*, green Area (GA), normalized green-red difference index (NGRDI) and triangular greenness index (TGI), presenting a reduction in plant growth as well as in green and yellow color components as an effect of inoculation. These results were contrasted across two scenarios of mineral fertilization in soil and soil amended with two different olive mill waste composts, presenting a clear interaction between the host-pathogen relationship and plant nutrition and suggesting the effect of V. dahliae infection during the biotrophic phase was not related to plant water status.


Introduction
Verticillium Wilt of Olive (VWO), a vascular disease caused by the pathogenic soil-borne fungus Verticillium dahliae Kleb, is one of the most severe threats to olive production today [1]. Taxonomic classifications place this fungus in Division Ascomycota, Class Sordariomycetes, Subclass Hypocreomycetidae, Order incertae sedis, Family Plectosphaerellaceae [2]. Other classifications place it within Deuteromycetes [3] along with other fungi without known sexual reproduction, since V. dahliae is a haploid species for which no teleomorph has been identified to date.
Verticillium dahliae (Vd) can survive short periods of time in the form of mycelium or conidia in infected perennial plants, but a specific form of melanized mycelium known as microsclerotia (MS) are the main survival structures and primary infectious propagules. MS are incorporated into soil from degrading host tissue, remaining viable for up to 14 years [4] establishing the non-parasitic phase of the biological cycle. The parasitic phase starts when favorable conditions stimulate MS germination, author's knowledge, although other remote sensing technologies have been used in assessing VWO, this was the first time that RGB vegetation indexes were used for this purpose.
Information captured digitally by taking pictures of plants is contained in the red, green and blue digital numbers for each pixel. When these images are processed, RGB values are converted to values in IHS (Intensity-Hue-Saturation) color space, (https://micro.magnet.fsu.edu/primer/digitalimaging/ russ/colorspaces.html). Intensity represents luminance information, Hue refers to the color itself and represents the average wavelength of the color spectrum and Saturation defines the color purity from white to the corresponding primary color. These IHS values are used to obtain vegetation indexes, some of which represent proportions within an image. Indexes fitting this criteria are the green area index (GA) [32], which calculates percentage of green pixels in a given image by defining green as 60 • < Hue < 180 • , or the greener area index (GGA) which has a stricter definition with green as 80 • < Hue < 180 • , thus excluding pixels with a yellowish-green hue generally associated with senescent plant matter. The difference between these two is used to determine the crop senescence index (CSI) which estimates the percent of senescent plant matter relative to total canopy [30].
Other indexes represent different qualities of color by providing a numerical value that describes the mean obtained from the different pixels in a given image; such is the case for properties of color Intensity, Hue, Saturation and Lightness (obtained from IHS color space), as well as values obtained from color spaces CIELAB (L*a*b*) and CIELUV (L*u*v*), proposed by the CIE (Comission Internationale de l'Éclairage) and which attempt to produce an intuitive representation of human color perception. In each of these color spaces L represents lightness, with a maximum value of 100 corresponding to a perfect reflecting diffuser and a minimum value of 0 representing black. Hue is represented on a two-dimensional plane orthogonal to L* and is defined by two axes or color opponent coordinates, named a* and b* in CIELAB (Figure 1), with a* running from green (−a*) to red (+a*) and b* from blue (−b*) to yellow (+b*) on a scale from 100 to −100. Axes u* and v* in CIELUV space are analogous to a* and b* respectively, although calculated using a different formula [33]. Figure 1 details colors that correspond to the 'green quadrant' (−a*, +b*) of CIELAB color space.
Remote Sens. 2019, 11, x FOR PEER REVIEW 3 of 19 kernels [29] or assessing grain yield losses and resistance to yellow rust in wheat [30,31]. To the author's knowledge, although other remote sensing technologies have been used in assessing VWO, this was the first time that RGB vegetation indexes were used for this purpose. Information captured digitally by taking pictures of plants is contained in the red, green and blue digital numbers for each pixel. When these images are processed, RGB values are converted to values in IHS (Intensity-Hue-Saturation) color space, (https://micro.magnet.fsu.edu/primer/digitalimaging/russ/colorspaces.html). Intensity represents luminance information, Hue refers to the color itself and represents the average wavelength of the color spectrum and Saturation defines the color purity from white to the corresponding primary color. These IHS values are used to obtain vegetation indexes, some of which represent proportions within an image. Indexes fitting this criteria are the green area index (GA) [32], which calculates percentage of green pixels in a given image by defining green as 60° < Hue < 180°, or the greener area index (GGA) which has a stricter definition with green as 80° < Hue < 180°, thus excluding pixels with a yellowish-green hue generally associated with senescent plant matter. The difference between these two is used to determine the crop senescence index (CSI) which estimates the percent of senescent plant matter relative to total canopy [30].
Other indexes represent different qualities of color by providing a numerical value that describes the mean obtained from the different pixels in a given image; such is the case for properties of color Intensity, Hue, Saturation and Lightness (obtained from IHS color space), as well as values obtained from color spaces CIELAB (L*a*b*) and CIELUV (L*u*v*), proposed by the CIE (Comission Internationale de l'Éclairage) and which attempt to produce an intuitive representation of human color perception. In each of these color spaces L represents lightness, with a maximum value of 100 corresponding to a perfect reflecting diffuser and a minimum value of 0 representing black. Hue is represented on a two-dimensional plane orthogonal to L* and is defined by two axes or color opponent coordinates, named a* and b* in CIELAB (Figure 1), with a* running from green (-a*) to red (+a*) and b* from blue (-b*) to yellow (+b*) on a scale from 100 to -100. Axes u* and v* in CIELUV space are analogous to a* and b* respectively, although calculated using a different formula [33]. Figure 1 details colors that correspond to the 'green quadrant' (-a*, +b*) of CIELAB color space. Other indexes such as the normalized green-red difference index (NGRDI) and the triangular greenness index (TGI) may be derived from conventional images as well as from multispectral information, also contributing information relating to plant physiological status. NGRDI approximates biomass and N status, while TGI is a solid indicator of chlorophyll content and is more accurate at a canopy scale compared to NGRDI [34,35,36]. Other indexes such as the normalized green-red difference index (NGRDI) and the triangular greenness index (TGI) may be derived from conventional images as well as from multispectral information, also contributing information relating to plant physiological status. NGRDI approximates biomass and N status, while TGI is a solid indicator of chlorophyll content and is more accurate at a canopy scale compared to NGRDI [34][35][36].
The objectives of this study were (i) to study the role of olive mill waste compost in Verticillium wilt suppressiveness, (ii) to elucidate the interaction between Verticillium infection and nutrient availability by contrasting different mixes of soil and compost in two scenarios, one with and one without mineral Remote Sens. 2019, 11, 607 4 of 17 fertilization, iii. to test RGB vegetation indexes in detection and assessment of Verticillium wilt and contrast obtained results with classical physiological measurements.

Experimental Design
In order to evaluate the effects of Verticillium Wilt of Olive in planta, a greenhouse experiment was conducted at the University of Barcelona Torribera Campus (41 • 27 48.4" N 2 • 12 53.0" E) employing 240 young olive tree clones (Olea europaea) of the cultivar Picual, which has been characterized as highly susceptible to Verticillium wilt [37]. Each cutting was transplanted into a 10L black plastic pot during November 2016 employing 6 distinct substrate mixes, based on a soil obtained from an olive orchard in Abrera, Catalonia, Spain in which VWO has not been detected.
For each of the soils or soil-compost mixes, a control and a V. dahliae-inoculated (I) group were established, yielding a total of 6 experimental conditions in each mineral fertilization scenario. For this purpose, a V. dahliae defoliating pathotype isolate provided by Dr. Manuel Avilés (Universidad de Sevilla), obtained from 6 Petri dishes in which the fungus had grown in Czapek-Dox Agar media for 4 weeks in the dark at 28 • C, was transferred to a sterile 40 L fermenter and incubated for 5 days at an aeration rate of 6 L min −1 (pO 2 was adjusted at 100 ±5%), 1000 rpm, 25 • C, and at a non-buffered pH in Czapek-Dox Broth, yielding a conidial suspension. Plants were inoculated with this suspension at 10 6 CFU/mL soil at the beginning of the experiment and were re-inoculated once every 6 months coinciding with the start of spring and fall. Control and inoculated plants were maintained separate within the greenhouse in order to prevent cross-contamination, while plants of the different treatments were otherwise randomized. Environmental conditions were semi-controlled, since fluctuations in temperature depended on the weather, but were adjusted to some extent by controlling the opening of the greenhouse roof. Watering was carried out regularly according to plant demand and weeds were controlled manually as needed. Of the total 240 plants, 48 were randomly selected for sampling, establishing 4 replicates for each of the 6 experimental conditions in each mineral fertilization scenario (4 replicates × 3(C) × 2(I) × 2(MF)). The current study was carried out during the springtime of the second year since the start of the experiment, with 5 samplings carried out at ±3 hours around solar noon during the last week of each month spanning from February to June 2018. Each sampling time point (T) was labeled: T1, T2, T3, T4 and T5.

Physiological Measurements
Physiological measurements were carried out separately in young and mature leaves, with young leaves selected from the third pair of new growth counting downwards from the apical meristem, and mature leaves selected based on stem width and a lower position on the plant, always maintaining a necessary minimum leaf area required for the fluorometer and porometer leaf clips.
Leaf stomatal conductance (G s ) was determined by measuring leaf abaxial epidermises using a Decagon SC-1 Leaf Porometer (Decagon Devices Inc., Pullman, Washington, USA). Chlorophyll fluorescence was measured using a MINI-PAM fluorimeter (Walz, Effeltrich, Germany) determining photosystem II quantum efficiency (φPSII) by exposing leaf adaxial epidermises to a 0.8s saturating flash at 6000 µmol photons m 2 ·s 1 in vivo in daylight conditions within the greenhouse, as well as in severed leaves that were dark-adapted for 1 hour in order to determine photosystem II maximum quantum efficiency (Fv/Fm), corresponding with all chlorophyll reaction centers being open.
Soil temperature at the time of sampling leaf stomatal conductance and chlorophyll fluorescence was measured with a Decagon Pro-Check (Decagon Devices Inc., Pullman, Washington, USA).

RGB Vegetation Indexes
RGB images were obtained using a Panasonic Lumix DMC-GX7 16-megapixel resolution camera (Panasonic Corporation, Osaka, Japan) with a G1:1.7/20 ASPH Φ46 0.2m/0.66ft-∞ lens, with aperture and ISO programed in automatic mode, and a Smartphone as remote shutter-release button controller through Panasonic Image App. The camera was mounted on a tripod and four images were obtained for each plant, all taken from a distance of 2.3 m from the camera to the edge of the pot: two from a zenithal plane and two from a lateral plane, rotating the plants 90 • for maximum stereoscopy. The images were taken ±1 hour around solar noon under conditions of natural light inside the greenhouse, in a chamber that had translucid walls and an opaque ceiling, thus providing considerable light diffusion. Cloudy days were avoided for sampling, since passing clouds would cause undesirable immediate changes in light conditions. Images were saved in a 4592 × 3448 pixel JPG format. Gamma-correction was not applied. The mean calculated from the results obtained from each of the 4 different images taken per plant was treated as the final result.
Each image was processed using open-source software FIJI, a modified version of Image-J (https: //github.com/fiji), first cropping images to desired size, and subsequently eliminating all pixels not corresponding to the plant (Figure 2), thus reducing background noise considerably. This was achieved by taking a picture in which the plant was removed from the frame immediately after every 4 plants were photographed, thus keeping changes in ambient light minimal. To all effects, for each image of each plant, another virtually identical picture was obtained without the plant. The algorithm "Difference" within FIJI "Image Calculator" (Process>Image Calculator) was then employed to subtract one image from the other creating a new image based on the premise that closely identical pixels in both images corresponded to the background and pixels that were different between images to the plant. This resulting image was converted from RGB to 8-bit grayscale (Image>Type>8-bit) in order to eliminate color, and was further treated (Image>Adjust>Threshold) via algorithm "Otsu" setting a threshold that resulted in all pixels being either black or white, effectively creating a binary mask in which white pixels represented the plant and black pixels the background. This mask was superposed (Process>Image Calculator) via the algorithm "Add" on the original image substituting 'background' pixels with black pixels but leaving all 'plant' pixels intact, effectively eliminating any information from color not corresponding to the plant in order for it to not interfere with final results. Each image was processed using open-source software FIJI, a modified version of Image-J (https://github.com/fiji), first cropping images to desired size, and subsequently eliminating all pixels not corresponding to the plant ( Figure 2), thus reducing background noise considerably. This was achieved by taking a picture in which the plant was removed from the frame immediately after every 4 plants were photographed, thus keeping changes in ambient light minimal. To all effects, for each image of each plant, another virtually identical picture was obtained without the plant. The algorithm "Difference" within FIJI "Image Calculator" (Process>Image Calculator) was then employed to subtract one image from the other creating a new image based on the premise that closely identical pixels in both images corresponded to the background and pixels that were different between images to the plant. This resulting image was converted from RGB to 8-bit grayscale (Image>Type>8-bit) in order to eliminate color, and was further treated (Image>Adjust>Threshold) via algorithm "Otsu" setting a threshold that resulted in all pixels being either black or white, effectively creating a binary mask in which white pixels represented the plant and black pixels the background. This mask was superposed (Process>Image Calculator) via the algorithm "Add" on the original image substituting 'background' pixels with black pixels but leaving all 'plant' pixels intact, effectively eliminating any information from color not corresponding to the plant in order for it to not interfere with final results.
FIJI plugin Breedpix was employed to determine different vegetation indexes based on image pixel color properties, resulting in high-resolution information calculated from pure pixels of either vegetation or background, as opposed to airborne images with aggregated pixels that combine FIJI plugin Breedpix was employed to determine different vegetation indexes based on image pixel color properties, resulting in high-resolution information calculated from pure pixels of either vegetation or background, as opposed to airborne images with aggregated pixels that combine information from vegetation and background. The CSI was determined from the formula CSI = 100 · (GA − GGA) / GA. NGRDI was obtained from the formula: NGRDI = (Green DN − Red DN)/(Green DN + Red DN), with DN standing for digital number, measuring the difference between green and red band reflectances with a possible range of −1 to 1. TGI is a calculation of the area within a triangle with (λr, Rr), (λg, Rg), and (λb, Rb) as vertices, where λ is the wavelength (nm) and R is the reflectance for red (r), green (g), and blue (b) bands. Images used to determine CIELAB and CIELUV color coordinates a*, b*, u* and v*, as well as color properties Intensity, Hue, Saturation, and Lightness were further processed to eliminate black background pixels that would otherwise be factored into calculated means, therefore making it necessary to filter out this information by converting these images using Matlab (Mathworks, Natick, Massachusetts, USA) into new images lacking all black or near-black pixels by using RGB value ≤ 20, 20, 20 as a threshold. For the purpose simplifying analysis of results, values obtained for index a* (which were in all cases negative since green hues are represented by a* values between 0 and −100) were converted to their absolute value, and are therefore expressed as positive.

Leaf Pigment Extractions
In order to better understand changes in color detected by RGB vegetation indexes, starting at time point two (T2), one young and one mature leaf were severed from each plant and stored at −20 • C. Chlorophyll a, chlorophyll b, and total carotenoids were extracted and quantified by the method of Hiscox and Israelstam [39] with slight modifications. 50 mg of leaf fresh weight were obtained from each leaf and cut into thin strips (~1 mm wide) and placed in vials containing 3.5 mL dimethyl-sulphoxide (DMSO) and incubated at 65 • C for 4 h. Liquid extracts were transferred to a graduated tube and made up to 5 mL. OD values at 665.1, 649.1 and 480 nm were read using a spectrophotometer against a DMSO blank. Equations from Sumanta et al. [40] were used to determine concentrations (µg/mL) of chlorophyll a (Chl a), chlorophyll b (Chl b) and total carotenoids (Car): Chl a = 12.47·Abs665.1 − 3.62·Abs649.1, Chl b = 25.06·Abs649.1 − 6.5·Abs665.1, Car = (1000·Abs480 − 1.29·Chl a − 53.78·Chl b)/220. Results were later adjusted to µg pigment/mg FW and the mean value obtained from the young and mature leaf sampled from each plant was used as the final result per replicate.

Data Analysis
Data was processed using Centurion Statgraphics XVII (Statpoint Technologies Inc., Warrenton, Virginia, USA) to test for statistically significant differences between treatments using a Repeated Measures Analysis of Variance (rANOVA) with significance at p ≤ 0.05. Fisher's LSD test at p ≤0.05 was used for multiple comparisons and to check for factor significance at each time point. The Kolmogorov-Smirnov test was used to check for normality of residues, and non-normal distributions were transformed before analysis. Mauchly's sphericity test was used and in cases in which this assumption was violated the Huynh-Feldt correction was applied when epilson > 0.75 and the Greenhouse-Geisser correction when epsilon < 0.75. Data is provided in the Supplementary Materials.

Results
During the course of the experiment inoculated plants remained asymptomatic, presenting no evidence of wilt. We were however, able to positively re-isolate and identify Verticillium dahliae from branch segments obtained from inoculated plants and grown out in Petri dishes in Czapek-Dox agar medium [41]. Soil temperatures recorded at the time of sampling physiological parameters ( Figure 3

Physiological Measurements
Results obtained from physiological measurements carried out to determine ɸPSII, Fv/Fm and Gs did not present significant differences between control and inoculated treatments ( Table 2). Chlorophyll fluorescence measurements (ɸPSII, Fv/Fm) did however prove sensitive to compost amendments in mature leaves in both MF scenarios, showing significantly lower values for Soil-C1 plants when compared with Soil and Soil-C2 plants (Table 3). Values obtained for Gs were only significantly affected by the sampling time point (Table 2).

Physiological Measurements
Results obtained from physiological measurements carried out to determine φPSII, Fv/Fm and G s did not present significant differences between control and inoculated treatments ( Table 2). Chlorophyll fluorescence measurements (φPSII, Fv/Fm) did however prove sensitive to compost amendments in mature leaves in both MF scenarios, showing significantly lower values for Soil-C1 plants when compared with Soil and Soil-C2 plants (Table 3). Values obtained for Gs were only significantly affected by the sampling time point (Table 2). Table 2. p-values obtained from rANOVA for the experimental factors without mineral fertilization (MF-) and with mineral fertilization (MF+).

RGB Vegetation Indexes
RGB vegetation indexes Saturation, a*, b*, GA, NGRDI and TGI showed significant differences between control and inoculated plants in the MF-scenario (Table 2), while inoculation was not significant in MF+. All RGB indexes except lightness, v* and CSI were sensitive to compost amendments in the MF-scenario, and all except b* and v* in MF+. As for the effect of compost (Table 3), results obtained from RGB indexes measuring greenness (a*, u*, GA, GGA, NGRDI, TGI) in all cases showed higher values in correspondence with the addition of compost in MF-.
Index a* showed significantly higher values in Soil-C2 compared with Soil-C1 and Soil, GA and GGA showed higher values in Soil-C2 and Soil C1 compared with Soil, and values for u*, NGRDI and TGI increased significantly from those obtained for Soil when C1 was added, and once again when comparing Soil-C2 to Soil-C1. Indexes measuring yellow hues proved less sensitive, since only b* showed a significant increase in obtained values with the addition of compost. A significant increase in Hue was observed in compost amended treatments compared with Soil, while in the case of Saturation an increase was significant in Soil-C2 compared with Soil-C1 and Soil. Results obtained in the MF+ scenario showed a similar trend, with compost amended plants showing higher values for most parameters when compared with non-amended plants, although differences between the three soil-compost mixes were in certain cases reduced. Overall values obtained in MF+ increased in comparison with MF-for Hue and indexes measuring greenness, and remained similar or showed a slight reduction for Saturation and indexes measuring yellow hues. Table 3. Means obtained from soil and soil-compost mixes without mineral fertilization (MF-) and with mineral fertilization (MF+). Upon representing results obtained for RGB indexes GA, NGRDI and TGI, such that each combination of soil mix and inoculation is shown at each time point (Figure 4), it is evident in MFthat results obtained from Soil treatments were the overall lowest, while the addition of C1 showed an increase in obtained values, and the addition of C2 a further increase. Values obtained from inoculated plants were lower than their respective controls, although differences due to inoculation were reduced in both compost treatments for GA and in Soil-C2 for NGRDI and TGI, a trend that is accentuated at T4 and T5, resulting in the inoculation factor overall only showing statistical significance for GA between T1 and T3, at T2 for NGRDI and T2 and T3 for TGI. In comparison, values obtained from MF+ show that values from Soil treatments were higher than in MF-, with mean values closer to those obtained from compost treatments. Values in the MF+ scenario did not show significant differences due to inoculation, and differences between soil mixes were reduced, although in general Soil-C2 presented the highest values.  Values obtained from CIELAB color space represented as a function of a* and b* color-opponent coordinates graphically represent mean pixel color on a two dimensional-plane ( Figure 5). Differences between the means obtained for control or inoculated treatments depict the inoculation effect vector at the different sampled time points. In both MF-and MF+, inoculated treatments showed a reduction of a* and a reduction of b* when compared with controls in all cases except at T5 in MF-, where a negligible increment in a* could be observed. Significant differences in vector magnitude between controls and inoculated plants for both a* and b* color opponent coordinates were found in MF-at time points T1, T2 and T3, coinciding with the early season samplings, while no significant differences were observed in vector direction in either MF-or MF+. Values obtained from CIELAB color space represented as a function of a* and b* color-opponent coordinates graphically represent mean pixel color on a two dimensional-plane ( Figure 5). Differences between the means obtained for control or inoculated treatments depict the inoculation effect vector at the different sampled time points. In both MF-and MF+, inoculated treatments showed a reduction of a* and a reduction of b* when compared with controls in all cases except at T5 in MF-, where a negligible increment in a* could be observed. Significant differences in vector magnitude between controls and inoculated plants for both a* and b* color opponent coordinates were found in MF-at time points T1, T2 and T3, coinciding with the early season samplings, while no significant differences were observed in vector direction in either MF-or MF+.

Leaf Pigment Extractions
Values obtained from pigment extractions detected significant differences due to inoculation in the MF-scenario (except for the case of the Car/Chl ratio), while values obtained from MF+ were not sensitive to inoculation (Table 2). Compost treatments (Table 3) also showed significantly higher pigment concentrations compared to soil in MF-, while in MF+ the only result that proved sensitive to compost amendments was carotenoid content, which was significantly lower in Soil when compared with Soil-C1, with Soil-C2 showing intermediate values.
Results obtained for chlorophyll a, chlorophyll b, and total carotenoid concentrations from the different treatments across sampling time points T2-T5 ( Figure 6) showed less consistent trends when compared with RGB indexes. Significant differences due to inoculation were observed in MF-, and although overall, values were higher in inoculated treatments, T3 shows a reversal of this tendency for Soil-C1 plants in all three measurements as well as for the Soil treatment for carotenoid concentrations, while T4 shows lower values of Chl a for inoculated Soil-C1 plants. Differences due to compost were also significant for all three pigment measurements in MF-, presenting lower values in Soil treatments compared with Soil-C1, and in turn with Soil-C2, changes in this trend at certain time points resulted in the compost and time interaction being significant for all three pigment measurements. Similarly to what was observed for RGB indexes, values obtained from Soil

Leaf Pigment Extractions
Values obtained from pigment extractions detected significant differences due to inoculation in the MF-scenario (except for the case of the Car/Chl ratio), while values obtained from MF+ were not sensitive to inoculation (Table 2). Compost treatments (Table 3) also showed significantly higher pigment concentrations compared to soil in MF-, while in MF+ the only result that proved sensitive to compost amendments was carotenoid content, which was significantly lower in Soil when compared with Soil-C1, with Soil-C2 showing intermediate values.
Results obtained for chlorophyll a, chlorophyll b, and total carotenoid concentrations from the different treatments across sampling time points T2-T5 ( Figure 6) showed less consistent trends when compared with RGB indexes. Significant differences due to inoculation were observed in MF-, and although overall, values were higher in inoculated treatments, T3 shows a reversal of this tendency for Soil-C1 plants in all three measurements as well as for the Soil treatment for carotenoid concentrations, while T4 shows lower values of Chl a for inoculated Soil-C1 plants. Differences due to compost were also significant for all three pigment measurements in MF-, presenting lower values in Soil treatments compared with Soil-C1, and in turn with Soil-C2, changes in this trend at certain time points resulted in the compost and time interaction being significant for all three pigment measurements. Similarly to what was observed for RGB indexes, values obtained from Soil treatments in the MF+ scenario were overall higher when compared to MF-, and differences due to inoculation or compost treatments were reduced, such that the only significant factor in MF+ was time.
treatments in the MF+ scenario were overall higher when compared to MF-, and differences due to inoculation or compost treatments were reduced, such that the only significant factor in MF+ was time. Correlations between values obtained from pigment extractions and RGB indexes (Table 4), expressed as Pearson's r-value, generally range between 0.2 and 0.6 (positive or negative) with overall stronger correlations found in the MF-scenario. Chl b and the Car/Chl ratio showed the strongest correlations with RGB indexes in MF-, while in MF+ the strongest correlations were with the Car/Chl ratio and carotenoid concentrations, followed by Chl b. Positive correlations were found between chlorophyll contents and RGB indexes measuring green, as well as the Car/Chl ratio and indexes measuring yellow, while negative correlations were found between chlorophyll contents and indexes measuring yellow or the Car/Chl ratio and indexes measuring green. In the case of carotenoid contents, correlations did not follow a consistent trend.

Discussion
V. dahliae infection has been characterized as consistent with water deficit stress as a consequence of xylem blockages caused by the fungus as well as host defense mechanisms [9,42,43]. Results obtained for stomatal conductance and chlorophyll fluorescence in this study (Table 2) did not show differences due to inoculation, and since both these parameters may be used as indicators of water status these results suggest effects that mimic water stress were not present, notwithstanding other physiological changes detected. This lack of symptoms and an apparently slow disease progress may be explained by the putative disease suppressiveness of the soil-compost mixes employed as well as environmental conditions, given that the olive cultivar Picual is highly susceptible to Verticillium wilt [37,44] and the fact the defoliating pathotype V. dahliae employed in this study has been characterized as more virulent than the non-defoliating pathotype [8].
Of the environmental conditions affecting our plant pathosystem, temperature may have played a defining role in slowing disease development, since values recorded at sampling times T3-T5 (Figure 3) were in a range that has been characterized as detrimental to fungal growth for Verticillium [18], purportedly allowing the biotrophic host-pathogen interaction to stagnate and impede its progress to the necrotrophic phase. In fact, in agricultural settings, expression of wilt symptoms has been characterized to occur in the springtime and fall coinciding with favorable ambient temperatures, while yearly variations in temperature also cause a much slower disease progress when compared to that observed in growth chamber experiments with a constantly favorable temperature for fungal growth [11]. The results obtained from RGB indexes (Figures 4 and 5) generally show significant differences due to inoculation from T1 to T3, when temperatures were under 27 • C. A lack of differences detected for stomatal conductance suggests xylem blockages causing symptoms similar to water deficit stress may not occur in the biotrophic phase and possibly be exclusive of necrotrophic host-Vd interactions.
As detailed in Figure 1, a reduction in both green and yellow color components in CIELAB space tend toward gray, in effect reducing color saturation, and therefore coinciding with significantly lower values obtained for Saturation in inoculated plants. In this light, differences in the inoculation effect vector observed in Figure 5 can be described as (i) changes in vector magnitude, therefore a CIELAB equivalent to the Saturation result obtained from IHS space, or (ii) changes in vector direction, defined by the a*/b* ratio and equivalent in CIELAB space to the Hue result calculated from IHS, further supported by the highly similar correlations with pigment concentrations obtained for Hue and the a*/b* ratio ( Table 4).
The fact that the lower values obtained from indexes GA, NGRDI and TGI in inoculated plants were more pronounced in absence of compost amendments in the MF-scenario, coupled with the fact that values obtained for Soil treatments from indexes GA, NGRDI and TGI ( Figure 4) were higher in the MF+ scenario compared with MF-, as well as higher in general in the presence of compost amendments (Table 3) suggests host-pathogen interactions in inoculated plants had an effect that acted in an opposite direction to plant nutrition. This hypothesis is further substantiated by generally higher values found in Soil-C2 compared with Soil-C1, considering Soil-C2 was more fertile (Table 1), results obtained for GA ( Figure 4) which highly correlate with plant biomass [45] inferring increased plant growth as a consequence of mineral fertilization and compost amendments along with decreased plant growth due to inoculation, as well as a reduction in the inoculation effect vector magnitude when comparing MF+ to MF-( Figure 5), with no important changes in vector direction observed between mineral fertilization scenarios or across sampling time points. These findings, showing a clear interaction between host infection and nutrition are in accordance with results obtained in an assay carried out during the previous year employing the same experimental design [41], which showed reductions in leaf growth as well as reduced leaf N, K and S inputs in inoculated plants. Reduced differences observed between control and inoculated plants in results obtained from RGB indexes can thus be attributed to improved nutrition conferred by both mineral and organic fertility, although suppressiveness as a consequence of microbial activity from compost amendments cannot be ruled out and calls for further study.
The relationship between the observed changes in plant color inferred from RGB indexes and the results obtained from pigment extractions was, however, not clear-cut. Chlorophyll and carotenoid concentrations (Table 3) were higher in the MF+ scenario compared with MF-, and concentrations in the MF-scenario were higher in compost amended treatments when compared with Soil, coinciding with increases in green due to both compost and MF, measured by a*, u*, GA, GGA, NGRDI, and TGI, and decreases in yellow color, measured by b* for compost and MF and v* and CSI for MF.
In the case of inoculation however, results showed lower values of GA, NGRDI and TGI in inoculated plants in MF- (Figure 4), as well as a reduction in both green and yellow color components corresponding to a* and b* ( Figure 5), while chlorophyll and carotenoid concentrations were slightly higher in inoculated plants (Figure 6), suggesting symptoms of chlorosis characterized in the literature as a consequence of Verticillium wilt [46,47] may correspond to the necrotrophic phase. These findings seem to indicate that the relationship between pigments and plant color inferred from RGB indexes is complex and may depend on other factors not contemplated by our statistical model, such as leaf thickness or leaf angle in relation to the plant, since abaxial and adaxial leaf surfaces present clear differences of color in olive plants.
Results obtained for pigment concentrations and RGB indexes (Table 4) showed significant correlations, and excepting those obtained for carotenoids, in all cases Pearson's r-values pairing indexes measuring green with chlorophyll contents, or yellow with the Car/Chl ratio were positive, while an opposite pairing resulted in negative values. It is likely the effect of improved nutrition as a consequence of mineral fertilization and compost was the main factor contributing to positively to these correlations, with the opposite effect due to inoculation presumably subtracting from them. In fact, the lack of a consistent trend observed for correlations between RGB indexes and carotenoid contents is not remarkable, since estimation of carotenoid contents from spectral reflectance has proven difficult due to an overlap in chlorophyll and carotenoid absorption peaks, with estimation of Car/Chl ratio proving more successful [48]. The fact that the Car/Chl ratio showed high correlations overall in both MF-and MF+ scenarios can be thus interpreted as not coincidental, since perceived color is logically a result of the overlap of pigments of different colors. Since this ratio did not show, by itself, significant differences between control and inoculated plants, the changes in color measured by RGB indexes may not depend solely on leaf pigmentation, since the effect of an increase in both pigments' concentration due to inoculation that does not substantially modify their ratio may be offset by other factors, not entirely explaining nor contradicting results measured from RGB image analysis method.

Conclusions
This study presents evidence suggesting the effects of V. dahliae infection during the biotrophic phase have little or no relation to water deficit stress and a strong interaction with plant nutrition, showing an opposite effect to both organic and mineral fertilization, since in cases in which values obtained for the studied parameters in inoculated plants were lower than controls, organic or mineral fertilization resulted in higher values.
Compost amendments and mineral fertilization reduced the effect of the pathogen as measured from image analysis, suggesting improved plant nutrition may contribute positively towards plant disease tolerance in Verticillium Wilt of Olive.
RGB vegetation indexes, widely used in plant phenotyping and in assessing abiotic stress, showed considerable resolution in detecting changes in plant color that could be attributed to the inoculation factor, especially notable in a context of a lack of wilting symptoms. Although remote sensing methods have already proven effective in assessing VWO symptomatic plants showing considerable wilt, this study presents data showing promise in early detection of VWO as well as other biotic stress detection and assessment, warranting further study in order to develop reliable methods of detection that may be used with images obtained from aerial platforms as a cost-effective method for scanning large areas in field scenarios. The main challenge to be faced in future applications is differentiating between the effects of biotic and abiotic stress, due to complex interactions present between crop physiological status and pathogen infection.