Next Article in Journal
RWLMod—Potential Model to Study Plant Tolerance in Drought Stress Conditions
Next Article in Special Issue
Sustainability Trait Modeling of Field-Grown Switchgrass (Panicum virgatum) Using UAV-Based Imagery
Previous Article in Journal
Proteomics and Interspecies Interaction Analysis Revealed Abscisic Acid Signalling to Be the Primary Driver for Oil Palm’s Response against Red Palm Weevil Infestation
Previous Article in Special Issue
The Impact of Different Types of Hydrocarbon Disturbance on the Resiliency of Native Desert Vegetation in a War-Affected Area: A Case Study from the State of Kuwait
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hyperspectral Reflectance Response of Wild Rocket (Diplotaxis tenuifolia) Baby-Leaf to Bio-Based Disease Resistance Inducers Using a Linear Mixed Effect Model

1
Council for Agricultural Research and Economics (CREA), Research Centre for Vegetable and Ornamental Crops, Via Cavalleggeri 25, 84098 Pontecagnano Faiano, Italy
2
Council for Agricultural Research and Economics (CREA), Research Centre for Vegetable and Ornamental Crops, Via Salaria 1, 63030 Monsampolo del Tronto, Italy
3
Council for Agricultural Research and Economics (CREA), Research Centre for Agriculture and Environment, Via Celso Ulpiani 5, 70125 Bari, Italy
4
Department of Engineering and Geology (InGeo), “Gabriele D’Annunzio” University of Chieti-Pescara, Via dei Vestini 31, 66013 Chieti, Italy
*
Author to whom correspondence should be addressed.
Plants 2021, 10(12), 2575; https://doi.org/10.3390/plants10122575
Submission received: 11 November 2021 / Revised: 23 November 2021 / Accepted: 24 November 2021 / Published: 25 November 2021

Abstract

:
Baby leaf wild rocket cropping systems feeding the high convenience salad chain are prone to a set of disease agents that require management measures compatible with the sustainability-own features of the ready-to-eat food segment. In this light, bio-based disease resistance inducers able to elicit the plant’s defense mechanism(s) against a wide-spectrum of pathogens are proposed as safe and effective remedies as alternatives to synthetic fungicides, to be, however, implemented under practical field applications. Hyperspectral-based proximal sensing was applied here to detect plant reflectance response to treatment of wild rocket beds with Trichoderma atroviride strain TA35, laminarin-based Vacciplant®, and Saccharomyces cerevisiae strain LAS117 cell wall extract-based Romeo®, compared to a local standard approach including synthetic fungicides (i.e., cyprodinil, fludioxonil, mandipropamid, and metalaxyl-m) and a not-treated control. Variability of the spectral information acquired in VIS–NIR–SWIR regions per treatment was explained by three principal components associated with foliar absorption of water, structural characteristics of the vegetation, and the ecophysiological plant status. Therefore, the following model-based statistical approach returned the interpretation of the inducers’ performances at field scale consistent with their putative biological effects. The study stated that compost and laminarin-based treatments were the highest crop impacting ones, resulting in enhanced water intake and in stress-related pigment adjustment, respectively. Whereas plants under the conventional chemical management proved to be in better vigor and health status than the untreated control.

1. Introduction

Wild rocket (Diplotaxis tenuifolia (L.) DC) is a cruciferous perennial herb, spontaneous in the Mediterranean Basin. In the last 25 years, Italy has become one of the major European producers of wild rocket using the species in intensified cultivation systems devoted to harvesting fresh-cut baby-leaf for the high convenience salad chain. It is sown with precision seed drills on 1.8–2.2 m width beds under polytunnels, fertigated with sprinklers and mechanically cut at complete foliar development. The harvested product quickly enters the cold chain, is minimally processed by washing, wiping, and bagging, and distributed through retail nets across several countries as ready-to-eat preparations. Packaged wild rocket meets the consumer preferences for the characteristic pungent–aromatic flavor associated with glucosinolates [1] and a few other nutraceutical properties (i.e., vitamins, antioxidants, fibers, and low calories) [2].
This makes the market very sensitive to the sustainability levels of the production process from field to shelf, conceived with a considerable reduction in the applied synthetic fungicides [3,4]. Nevertheless, this vegetable crop as well as all the other baby leaf species is susceptible to a plethora of both specific and non-specific pathogens that significantly reduce yields and impair their quality. As a consequence, the deployment of non-traditional effective control measures that include biological strategies is necessary [5].
Biologically-based disease resistance inducers may be biomimetic compounds or substances sourced from plants or microbes, or non-pathogenic microorganisms capable of eliciting the plant’s own defense mechanism(s) through microbe/pathogen-associated molecular patterns and/or the recognition of host-derived damage-associated molecular patterns to enhance their innate defense response against upcoming broad-spectrum diseases to varying degrees [6,7,8]. Therefore, innovative resistance activators can then be used as biopesticides in plant protection protocols as a safer alternative to synthetic chemicals to reduce the environmental disease management footprint and stimulate plant performances [9]. However, they still need to improve their efficacy under field conditions [10].
Non-invasive technologies may be helpful to optimize the field applications of these plant-targeted protectants [11]. Plants react to exogenous application of plant resistance inducers by possibly activating many metabolic pathways involved in biochemical and mechanical defense responses including shifts of cytosolic ion content, oxidative burst, synthesis of enzymes, proteins, and other secondary metabolites related to the defense, in addition to the activation of resistance-related hormones [12]. Changes in leaf composition and plant health may be detected in the reflected electromagnetic radiation once it is captured by optoelectronic sensors such as hyperspectral ones.
In optics, reflectance measures the ability of a given surface or material to send back part of the incident light on it. In particular, a hyperspectral sensor can record the part of the electromagnetic radiation of a natural (sun) or artificial light source that is reflected by a leaf at a very fine spectral resolution in the range of wavelengths between 350–2500 nm. The spectrum is divided into three regions called, in sequence: Visible (VIS), between 350 and 750 nm; Near Infrared (NIR) until 1400 nm; and the Short-Wave Infrared (SWIR) region until 2500 nm. Each region has been associated with many parameters describing the plant status [13]. Hyperspectral data analysis represents a very effective and sustainable tool for evaluating changes induced in the plant by abiotic and biotic stresses. Recently, hyperspectral data have been adopted on a large scale in the detection of biotic stresses on plants as in the case of the sudden spread of Xylella infection on entire olive cultivations in the Mediterranean area [14,15,16]. Some vegetation indices based on reflectance data in the VIS range have been used to assess changes occurring in the plant health status through the related effect on pigments such as carotenoids, anthocyanins, and chlorophyll. However, these variations might be due to either specific or non-specific infections [17,18,19,20,21], hence further laboratory phytopathogenic analysis on the leaves remains necessary.
The transition from low reflectance values in the red to high values in the infrared spectral range is very rapid: this portion of the spectrum, called Red Edge, is more indicative of the chlorophyll content than that of water [22,23,24,25]. Moreover, it is influenced by the cell structures of leaves that poorly absorb in the NIR because of the multiple scattering of radiation by the mesophyll. As for SWIR, the overlap between information on water content and organic compounds makes data interpretation more difficult. Statistical processing [26], mathematical regressions [16], and radiative transfer models [27,28,29] have been used for this purpose.
This leads to hyperspectral data being used as indicators of possible stress in the plant, although a direct relationship between the alteration observed in the spectra and its cause has not been discovered yet.
The analysis and interpretation of spectral data are further complicated by the way in which agronomic trials are generally conducted. The purpose of the agronomic experiments is to test whether the compared treatments (in this case the use of different plant resistance inducers) have any effect on the supposed response variable. If one treatment is taken as a control (e.g., zero treatment), the experiment will consist of testing whether every other treatment has an effect compared to the control treatment. The biggest challenge in an agronomic trial is to be able to separate the intrinsic variation of the response variable from that induced by the experimental treatments. In traditional agronomic trials, this is achieved by replicating each treatment according to a well-defined experimental design. Traditional statistical methods based on the design then allow for the determination of the probability that any measured difference between treatments is due to chance (null hypothesis).
Many times, when the experiment is conducted in a confined space such as a greenhouse, or on-farm, for purely practical reasons, there is a tendency to follow a more systematic pattern, with one treatment, for example, assigned to a particular part of the field. In addition, there may be too few plots per treatment (repetitions) to assess the underlying variability, and furthermore, such variability may be correlated [30].
These experiments very often fail to meet the fundamental assumptions required by classical statistical methods. It is therefore necessary to use more complex statistical methods [31,32] that are based on a model-based statistical approach [33]. This consists in describing both the variation and the correlations between the observations of the response variable using a statistical model. In this regard, the theory of linear mixed effects models (LMM) [34,35,36] allows for the total variance to be broken down into that which is attributable to fixed effects, corresponding to the treatments, and that which is attributable to random effects. The latter are linked to the intrinsic spatial variability of the agronomic system, which cannot be described by fixed effects and can be estimated by the covariance/correlation function of residuals.
Therefore, with a view to increasing sustainability in greenhouse cultivations of wild rocket, the use of naturally derived products can be proposed, but their actual interaction with the plants must be continuously monitored. For this purpose, among micro- and macroorganisms, certain active components described as plant defense promoters can be found and formulated in some successful experimental and commercial products to induce resistance in the plants. These may differ in the level of purity of the active ingredients, ranging from microorganisms to a single molecule, in the hormonal signaling pathways involved in the elicited plant reactions (i.e., oxidative burst, cell-wall fortification, etc.) and in the ranges of efficacy. Furthermore, proper organic soil management with compost amendments that promote plant growth and development by improving soil chemical, physical, and microbiological quality and fertility can synergize with inputs to the aerial part, leading to greater beneficial effects.
The specific objective of the work is to characterize the hyperspectral response of greenhouse-grown wild rocket to the application of three immunity-stimulating active ingredients (resistance inducers), split over plots amended and non-amended with green compost, using the mixed effect model theory to account for the actual experimental conditions.

2. Results

2.1. Principal Component Analysis (PCA)

Four principal components (PCs) with eigenvalues greater than 1 were extracted. Among these, however only three PCs explaining a cumulative variance proportion greater than 97% were retained (Table 1).
The most relevant bands for PC1 fell within the ranges 1405–1545 nm and 1855–2500 nm with a small plateau at 1435–1475 nm, a peak at 1885 nm, and two other plateaus at 2055–2105 and 2345–2385 nm. Since water peaks fall in SWIR regions (wavelengths centered at about 1450, 1940, and 2200 nm), PC1 appears to be related to foliar absorption of water [29,37]. However, these absorbance bands may also be related to the characteristics of chemical components such as cellulose, starch, and proteins. In particular, the peak at 1885 nm could be attributable to a stretch absorption mechanism of chemical groups such as OH and CO in cellulose and hemicellulose, and the plateau at 2345–2385 seems to be linked to a deformation mechanism of the CH bond [27]. The range 2055–2105 nm was highly correlated with N content [38], while the range 1435–1475 nm overlapped with the main absorbance peak of water at 1450.
PC1 could therefore be used as an indicator of the leaf water status as well as the content of cellulose, starch, and proteins, whose variations might be caused by mechanisms of action of the inducers through the biochemical pathways they influence (i.e., the closure of stomata) [39,40,41,42]. It is worth noting that since the data refer to reflectance, the relationship between PC1 and water, cellulose, starch, and protein contents is inverse.
Regarding PC2, the main loadings fell within the range 715–1375 nm. More specifically, the subrange up to 935 nm was more related to the structural characteristics (LAI) of the vegetation [22]. The secondary water absorption peak at 970 nm, often used to compute vegetation indices of water content [43], fell into the next subrange, as did the other peculiar water absorbance peak at 1200 nm [44]. PC2 can therefore be interpreted primarily as an indicator of the structural characteristics of the leaf cells and only secondarily as an indicator of the water content. Unlike PC1, the relationship between PC2 and LAI is direct as it is well-known that better leaf vigor produces higher reflectance in the NIR [22].
PC3 was considered, despite the loadings not exceeding the established limit value of 80, because of the additional information it can provide regarding the spectral response of wild rocket to the action of inducers. The path of its loadings (Figure 1) showed two peaks between 365–695 nm in VIS. The first peak had a maximum at 495 nm, a value related to the carotenoid content as an indicator of plant stress [45,46], while the second peak at 650 nm was related to chlorophyll content, as green plants have a main absorption peak in the red. In fact, the two leaf contents are closely related and their Car/Chla ratio can be used as an indicator of the ecophysiological state of leaves and plants [47]: a high ratio denotes severe stress.
Therefore PC3, despite the small portion of variance explained, can be taken as an indicator of plant health as high reflectivity in this VIS range might be associated with a low carotenoid content, intended as an indicator of stress, and efficient chlorophyll function.

2.2. Linear Mixed Model (LMM)

Before any further processing, the first three principal components, standardized to mean 0, were tested against the assumption of normal distribution. From the examination of Table 2, only PC2 showed significant deviations from normality, as can be verified by the results of the basic statistics and normality tests. It was therefore necessary to perform the relative rank transformation only for PC2, since the pronounced skewness of its distribution might affect the linearity of the mixed model to be estimated. From this point onward, all statistical processing is understood to be carried out on the rank-transformed PC2 values (rPC2) and the raw values of PC1 and PC3.
The results of the Levene’s test applied to the three variables are shown in Table 3. Assuming a reference probability level of 0.05, the fixed effect COMPOST is non-homogeneous only for PC1, while the effect TREATMENT is non-homogeneous for PC1 and PC3.
Both Moran’s I test and Geary’s c test were significant at p < 0.05 for all three variables (Table 4). However, these are overall tests: more information on the type and intensity of the spatial association can be obtained from residual-fitted variogram models.
Regarding the residual variogram models, a spherical model with nugget was fitted for PC1 and PC3 and a pure nugget effect model for rPC2 (Table 5). Although the parameters of the variogram models were significantly different from zero, it can be seen that the nugget effect for PC1 and PC3 was 87% and 96%, respectively, of the total sill (nugget + partial sill). These results show that although spatial autocorrelation occurred for PC1 and PC3, the structured component (partial sill) was only a small percentage (13% and 4%, respectively) of the total spatial variance because the spatially uncorrelated component (nugget effect) was highly prevalent. This can be explained based on the particular experimental scheme in which an intrinsic variability of the plant–soil system was superimposed on an artificial variability imposed by man through randomization in treatment assignment. However, these variogram models were only used to initialize the iterative process for jointly estimating the random effects parameters and fixed effects coefficients of the mixed effects models to be determined.
Table 6 shows the results of the mixed-effects model estimation for PC1, separated between the estimates of the stochastic effect parameters (Table 6a) and the statistical significance of the fixed effects (Table 6b). Regarding the variogram of the residuals, two separate models were applied for the two levels of the COMPOST effect due to the heteroscedasticity of variance. For the NO level. a pure nugget effect was fitted, revealing a lack of spatial structure. For the YES level, it was possible to fit a spherical model, which, when compared with the previous least squares model (Table 5), had a shorter interval but a higher partial sill. In any case, it represents a short-range correlation with a spatially uncorrelated component more than twice as large as the structured one.
Unlike for the TREATMENT effect, it was not possible to fit five separate variogram models due to insufficient degrees of freedom because of the small number of observations/replications.
Regarding the fixed effects, both main effects and their interaction were highly significant.
Only the most significant differences between the various levels of the two effects and their interaction at p < 0.05 are shown in Table 7.
For the COMPOST effect, the mean of the NO level for PC1 was larger than the mean of the YES level. Recalling the meaning attributed to PC1, it can be stated that plants treated with compost might have a higher water content.
For the effect of TREATMENT, CTR mean was significantly lower than that of the CER, TRI, and CHE levels and therefore might be associated with higher water content.
For rPC2, it was possible to fit a spherical model with nugget for the random effect (Table 6), in contrast to the least squares model (Table 5), which was the pure nugget effect. However, once again, this was short-range variability (<3 m) and with a spatially structured component that was only 14% of the spatially uncorrelated.
Regarding the fixed effects, only TREATMENT was significant, and in Table 7, only the differences between the TREATMENT levels that were significant are shown. Based on the interpretation of rPC2, which correlated positively with LAI, we can state that CHE treatment showed better leaf vigor than all remaining levels.
For PC3, there was sufficient agreement between the least squares estimates (Table 5) of the random effect and the RML estimates (Table 6a). However, in the latter case, there was a slight increase in the spatially structured component, which remained a small part (11%) of the spatially uncorrelated component. Both fixed effects and their interaction were highly significant (Table 6b).
Based on the interpretation of PC3, plants treated with compost appeared to be healthier than untreated plants and the same was true for plants under CHE treatment compared to those in the other levels (Table 7). These results are consistent with those obtained with rank transformed PC2.
In summary, the above analysis showed that the compost-treated plants had a higher water content than untreated plants. Similarly, the plants in level CTR of TREATMENT appeared wetter than those in the remaining levels.
In contrast, the plants in CHE appeared to be more luxuriant and healthier than those in the other levels. Statistical analysis showed that the apparently winning combination (YES-CHE) differed significantly for both PC1 and PC3 from the combinations YES-CTR and NO-CTR, resulting in lower water content but lower stress level (Table 7).
The results of the traditional analysis of variance, assuming that the residuals were spatially uncorrelated, are also shown for comparison in Table 8. Indeed, there were no significant differences and this was quite consistent with what was expected, given the small proportion of the structured component of the variance compared to the spatially uncorrelated component. However, the values of the F statistic tended to be higher than those of the mixed effects models, underlining the risk of committing type I errors if the spatial correlation is neglected [48].

3. Discussion

The statistical results obtained are consistent with a biological interpretation, which reinforces the idea that the spectral response of the plant can be used as an effective and reliable indicator of its health. PC1 summarizes information about N content and other biochemical compounds, which to date, has come from several studies regarding the SWIR region on potato and other mapped vegetation [49,50]. As far as PC2 is concerned, the main loadings fell in the ranges more related to the LAI as confirmed by recent studies on rice and maize in both proximal and remote sensing [51,52], and to the water absorption peaks that are used to calculate new vegetation indices associated with water content in different plant species [53]. Finally, PC3 can explain the ratio Car/Chla as an indicator of plant stress [54].
It is generally recognized that the incorporation of compost into the soil increases the water available to plants [55], delaying the possible wilting associated with drought [56] and thus protecting and/or enhancing photosynthetic activity [57]. On the other hand, the lack of statistical significance of the COMPOST effect for rPC2, which was associated with plant vigor, can be explained on the basis of the reduced nutrient supply capacity shown by green composts in the presence of a large fraction of non-labile carbon, whose degradation implies the net immobilization of N [58]. In addition, compost in combination with LAM and TRI had a positive effect on the water content (PC1) and in combination with LAM and CHE on the general health of the plant (PC3). In contrast, compost in combination with CER did not produce any positive effects in terms of either water content or LAI.
However, the resistance inducers used in this study, on the basis of their specific characteristics, might be implicated in the physiological processes underlying the interpretation of the PCs. Antagonistic fungi belonging to the genus Trichoderma have been reported to induce a resistance response into plants through multiple hormonal signaling pathways that modulate jasmonic acid, ethylene, and salicylic acid levels toward a wide-spectrum of pathogens [59]. Their biocontrol efficacy might result in the modulation of plant growth and yield improvement [60]. Kumar and Kumar [61] reported that root colonization of Trichoderma sp. can induce the production of stress enzymes such as peroxidase and glutathione reductase, which may be responsible for decreasing disease incidence in Brassica juncea. In a different way in cabbage, Trichoderma treatments increased the transcript levels of genes related to photosynthesis and sucrose transport, PR proteins, chitinases, and oxidases [62]. Yeast cell-wall extract, which carries polysaccharidic and peptidic polymers and oligomers of highly variable molecular mass, acts as MAMPs in inducing defense-related events through SA signaling [40,41,42,63]. However, there is no evidence in the literature that it has an impact on the reflectance of plants. On the other hand, concerning laminarin, which is a water-soluble glucan storage polysaccharide extracted from brown algae (i.e., Laminaria digitata Hudson, Lamouroux), it has been shown that it can elicit defense reactions in several plant species [64] via salicylic acid and reactive oxygen species pathways [65]. This is most likely due to the association with bound β-1,3–1,6 glucosyl residues [39]. It is also worth pointing out that in this study, the LAM effect on TRI was significantly higher in all pairwise COMPOST × TREATMENT interactions relative to PC3 associated with indicating stress occurrence. Consistently, laminarin has been reported as an unconventional elicitor of plant secondary metabolites [66]. In Arabidopsis, this molecule increased chloroplast stability by activating the antioxidant system under stress conditions [67]. Consequently, with regard to PC3, the current hyperspectral study indicated that LAM treatment associated with the compost effect is linked to the improved plant health status.

4. Materials and Methods

4.1. Experimental Design

The trial was conducted under a multi-tunnel greenhouse located at the CREA-Research Center for Vegetable and Ornamental Crops (Pontecagnano Faiano, Italy, 40°38′54.0″ N 14°53′21.4″ E). Wild rocket cv Yeti (Maraldi Sementi, Italy), suitable for the winter–spring cropping cycle, was precision sown on 19 November 2019 (2500 seeds m−2) and grown on 1.3 m wide beds under two polyethylene tunnels (20 × 7.2 m) with four beds each equally spaced by 0.5 m.
The experimental design was a two-way split-plot with the external effect (COMPOST) represented by the application or not of green compost as an amendant (two levels: Yes/No). The internal effect (TREATMENT) was represented by the application of inducers including five treatments: (1) two applications (at seeding and pre-emergence) of Trichoderma atroviride strain TA35 from CREA Collection [68] at a dose of 1000 L ha−1 (105 cells mL−1) (TRI); (2) four weekly spray applications of commercial laminarin-based product (Vacciplant®, Arysta Lifescience s.a.s., France) at dose of 1 L ha−1 (LAM); (3) four weekly spray applications of Saccharomyces cerevisiae strain LAS117 cell wall extract-based commercial product (Romeo®, Agrauxine, France) at a dose of 750 g ha−1 (CER); (4) local standard approach including one spray with Cyprodinil, Fludioxonil, Mandipropamide, and Metalaxil-M (CHE); and (5) no treatment (CTR) used as reference control. Main plots were not-amended or amended with green compost (10 t ha−1 dry matter) (COMPOST). The experimental design consisted of four blocks (each one including a single bed but split into two tunnels North/South); each block was split longitudinally into two halves randomly assigned to one level of COMPOST; each half of the block was randomly split into the five levels of the internal effect. The trial then consisted of 40 experimental units (plots) with a size of 3.0 × 1.3 m (Figure 2 and Figure 3).

4.2. Hyperspectral Reflectance Measurements

Leaf spectral measurements were performed in the spectral range 350–2500 nm, using a portable non-imaging spectroradiometer (FieldSpec® 4 Hi-Res, ASD Inc., Boulder, CO, USA) through a fiber-optic contact probe (ASD Plant Probe; ASD Inc., USA) with a spectral resolution of 3 nm in VIS-NIR and 8 nm in SWIR, a 10 mm field of view, and an integrated halogen reflector lamp. The instrument was warmed up for 90 min before measurements to increase the quality and homogeneity of acquired data. Calibration was obtained with the pre-calibrated white 99% spectral reference panel. Each sample scan represents a mean of 10 reflectance spectra.
Within each plot, 12 spots spaced 50 cm apart were selected based on the orientation of the leaves in sunlight, with preference given to the top and most exposed leaves; the reflectance measurement for each spot was obtained as an average of three replicates. The reflectance data were then averaged over 10-nm intervals, thus reducing the number of wavelengths from 2151 to 215, smoothing the spectra and keeping the risk of over-fitting low [69]. Pre-processing methods were applied to reduce the impact of multiplicative and additive effects of possible backscattering within the instrument [70]. The pre-processing first involved the splice correction, in order to minimize the inconsistency recorded at the spectral intervals among the three detectors of the spectroradiometer: the VIS-NIR range 350–1000 nm, which is characterized by a sampling interval of 1.4 nm and the two ranges in the NIR-SWIR, 1000–1800 nm and 1800–2500 nm, with a sampling interval of 2 nm. After that, two pretreatments were performed on the reflectance spectra: (1) multiple scattering correction (MSC); and (2) smoothing/denoising with Savitzky–Golay polynomials [71]. Multiple scattering correction (MSC) works mainly when the scattering effect is the dominant source of variability and removes additive and multiplicative components [72].
The first order Savitzky–Golay (SG) polynomial algorithm reduces the random noise of the measurements. The algorithm is based on a moving polynomial fit of any order and the filter size consists of (2n + 1) points, where n is the half-width of the smoothing window (w). The polynomial fit interpolates the points between the 2n’s. A window size (w) of 11 (w = 2n + 1) and the second polynomial order were applied here [73].
Splice correction was achieved using ViewSpecPro software (Analytical Spectral Devices Inc., Boulder, CO, USA). The other preprocessing methods were performed using ParLeS software [74].

4.3. Principal Component Analysis

Principal component analysis is a widely used dimensionality reduction technique that allows for the extraction of principal components (PCs), expressed as a linear combination of variables [75,76]. The PCs found in this way do not represent directly observable variables and must therefore be interpreted from a scientific-rational perspective. Mathematically, PCA works on the correlation matrix and extracts a few PCs equal to the number of measured variables, but only a part of them is useful.
The number of PCs was fixed by selecting the PCs with eigenvalues greater than or equal to 1 (Kaiser’s criterion) [77]. The eigenvalues refer to the share of variability “explained” by the PC and take on descending values as the first PC moves toward the last one. The result can then be subjected to rotation by various methods. Methods using orthogonal rotations preserve the independence of the PCs. The VARIMAX procedure was used in this study.
The most important parameters to evaluate were:
The amount of variance “explained” by both the set of the retained PCs (cumulative) and by each PC individually;
PC loadings, which describe the strength of the relationship between the PC and the variable being measured.
PCA was performed with the SAS/FACTOR procedure of the statistical software package SAS/STAT (release 9.4 SAS ANALYTICS U software).
For each PC, the loading values were multiplied by 100. The value of 80 was chosen as a limit above which to select the most relevant wavelengths.
The raw spectral data were preferred over the pretreated ones because with the same number of PCs, they explained a greater proportion of variance, probably due to the specific measurement modalities of the plant probe that reduced the random error.

4.4. Statistical Analysis

The assumption of univariate normality for each PC was checked with three non-parametric tests (Kolmogorov–Smirnov, Cramer-von Mises, and Anderson–Darling) [78].
If the variable showed large departures from normal distribution, to improve the linearization of the mixed effect model, PC was transformed to normal scores (yi) using Blom’s formula [79]:
y i = Φ 1 ( r i 3 / 8 ) ( n + 1 / 4 )
where Φ−1 is the inverse cumulative normal (PROBIT) function; ri is the rank of the ith observation; and n is the number of observations that have non-missing values for the ranking variable.
The spatial association of the residuals from the OLS models, obtained assuming independence of the residuals, was verified in two ways:
(1) By calculating the spatial autocorrelation statistics of Moran (1950) [80] and Geary [81] and comparing them to the null hypothesis statistics of completely random spatial model; and
(2) By fitting an authorized mathematical model to the experimental variogram of the OLS residuals, according to the least squares estimation (LSE) technique and testing the statistical significance of its parameters.

4.5. Linear Mixed Effect Model (LMM)

A LMM can be written in the form:
z = + ε
where z is a vector of length n corresponding to the measurements of the response variable; and M is the fixed effects design matrix of size n × t, where t is the total number of different experimental treatment levels. In the case under study, there are two levels for the COMPOST effect, five for the TREATMENT, and 10 for the interaction effect. β is the vector of length t of the coefficients of the fixed effects and ε the vector of length n of the random effects. The product represents the fixed effects. All entries of the matrix M (Mij) can assume the values 0 or 1. The entries of the first column (M1j) are all equal to 1 while the other Mij entries are 1 if observation i corresponds to the treatment level j, otherwise they are 0. Therefore, β1 corresponds to the overall mean and the other βj correspond to the adjustment of this mean concerning each treatment level j. The elements of the vector ε (ei), also known as residuals, are assumed to follow a normal distribution with mean 0 and covariance matrix Cov[ei, ej]. This matrix may differ from the identity matrix, and the random effects may also be spatially correlated.
Usually, the covariance is assumed to be a function of the distance between the locations xi and xj. If by h we denote the distance between xi and xj, the covariance model takes the general form:
Cov[ei, ej] = c(h) = σ [f(h)]
where σ is known as the partial sill and with the symbol f, one of the authorized mathematical functions used to represent a covariance function, is indicated. The models, for which f(h) is the same for all pairs of observations equally distant in a given direction, are called stationary models (of the second order if the average, assumed constant, is known). If f(h) does not depend on direction, then the covariance structure is said to be isotropic. Many isotropic covariance models exist, however, in this study, we used the spherical model:
c(h) = σ 0 + σ [1 − 1.5(h/ρ) + 0.5(h/ρ)3]  if h< ρ
c(h) = 0  if h≥ ρ
This model has three parameters that must be determined from the data: the nugget variance (σ 0), the partial sill (σ)m and the range ρ. These spatial models are drawn from geostatistics and refer to one of the numerous manuals [82,83,84] for an interpretation of these parameters. It is commonly preferred to use the variogram instead of the covariance function given by the formula:
γ (h) = σ 0 + σ − c(h)
Furthermore, complexity in determining the covariance model may arise from the heteroscedasticity of the variance relative to the individual PCs tested, as this may result in heterogeneity in the covariance function. The homogeneity of the variance was tested with Levene’s test [85].
The method of fitting the spherical model to the values of mean covariance for each class of distance, h, is based on an iterative procedure, aimed at maximizing the log likelihood of the residuals using the restricted maximum likelihood method (REML) [86]. Once the spatial covariance function is estimated, the fixed effects coefficients are obtained as generalized least squares estimates [86]. For comparison, traditional models (OLS), which assume residuals are normally distributed with zero mean but not correlated spatially, were also determined.
All statistical analyses were performed using statistical software package SAS/STAT (release 9.4 SAS ANALYTICS U software) and in particular, LMMs were estimated using the procedure MIXEDe.

5. Conclusions

This study demonstrated that hyperspectral proximal sensing of the wild rocket crop is suitable for tracking the field performance of canopy treatments with resistance inducers, based on their specific interactive mechanisms with the plant. Compost and laminarin-based treatments had the greatest impact on the crop, resulting in increased water uptake and in stress-related pigment regulation, respectively. Ultimately, the plant under the conventional chemical management proved to be more vigorous and healthier than the untreated control.
Furthermore, in order to produce reliable results in agronomic trials that are actually useful to farmers, knowledge of the spatial variability of the response variable can no longer be ignored. Otherwise, there is a risk of assigning to the experimental factor an effect that is instead attributable to the intrinsic variation of the crop–soil system. It is hoped that the theory of mixed effects models will become increasingly familiar in the agronomic community.

Author Contributions

Conceptualization, C.P. and A.C.; Methodology, A.C.; Software, A.G., C.R. and N.N.; Formal analysis, C.R. and A.C.; Investigation, C.P., A.G. and N.N.; Data curation, A.G., C.R. and A.C.; Writing—original draft preparation, C.P., C.R. and A.C.; Writing—review and editing, C.P. and A.C.; Supervision and project administration, C.P.; Funding acquisition, C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Agriculture, Food and Forestry Policies (MiPAAF), with the project AgriDigit, sub-project Tecnologie digitali integrate per il rafforzamento sostenibile di produzioni e trasformazioni agroalimentari (AgroFiliere) [DM 36503.7305.2018 of 20/12/2018].

Data Availability Statement

The data presented in this study are available in the article.

Acknowledgments

The field trial was conducted under a naturally ventilated polyethylene greenhouse patented CREA-OPUS ET VITA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bell, L.; Lignou, S.; Wagstaff, C. High glucosinolate content in rocket leaves (Diplotaxis tenuifolia and Eruca sativa) after multiple harvests is associated with increased bitterness, pungency, and reduced consumer liking. Foods 2020, 9, 1799. [Google Scholar] [CrossRef] [PubMed]
  2. Pennisi, G.; Orsini, F.; Castillejo, N.; Gómez, P.A.; Crepaldi, A.; Fernández, J.A.; Egea-Gilabert, C.; Artés-Hernández, F.; Gianquinto, G. Spectral composition from led lighting during storage affects nutraceuticals and safety attributes of fresh-cut red chard (Beta vulgaris) and rocket (Diplotaxis tenuifolia) leaves. Postharvest Biol. Technol. 2021, 175, 111500. [Google Scholar] [CrossRef]
  3. Minut, M.; Roșca, M.; Hlihor, R.M.; Cozma, P.; Gavrilescu, M. Modelling of health risk associated with the intake of pesticides from Romanian fruits and vegetables. Sustainability 2020, 12, 10035. [Google Scholar] [CrossRef]
  4. Bozdogan, A.M. Assessment of total risk on non-target organisms in fungicide application for agricultural sustainability. Sustainability 2014, 6, 1046–1058. [Google Scholar] [CrossRef] [Green Version]
  5. Gullino, M.L.; Gilardi, G.; Garibaldi, A. Ready-to-eat salad crops: A plant pathogen’s heaven. Plant Dis. 2019, 103, 2153–2170. [Google Scholar] [CrossRef] [Green Version]
  6. Jamiołkowska, A. Natural compounds as elicitors of plant resistance against diseases and new biocontrol strategies. Agronomy 2020, 10, 173. [Google Scholar] [CrossRef] [Green Version]
  7. Abbasi, S.; Sadeghi, A.; Omidvari, M.; Tahan, V. The stimulators and responsive genes to induce systemic resistance against pathogens: An exclusive focus on tomato as a model plant. Biocatal. Agric. Biotechnol. 2021, 33, 101993. [Google Scholar] [CrossRef]
  8. Torres-Rodriguez, J.A.; Reyes-Pérez, J.J.; Castellanos, T.; Angulo, C.; Quiñones-Aguilar, E.E. A vegetation research derived from aviris. In Proceedings of the Eighth Annual JPL Airborne Earth Science Workshop, Pasadena, CA, USA, 8–14 February 1999; pp. 8–14. [Google Scholar]
  9. Burketova, L.; Trda, L.; Ott, P.G.; Valentova, O. Bio-based resistance inducers for sustainable plant protection against pathogens. Biotechnol. Adv. 2015, 33, 994–1004. [Google Scholar] [CrossRef] [PubMed]
  10. Alexandersson, E.; Mulugeta, T.; Lankinen, Å.; Liljeroth, E.; Andreasson, E. Plant resistance inducers against pathogens in Solanaceae species—From molecular mechanisms to field application. Int. J. Mol. Sci. 2016, 17, 1673. [Google Scholar] [CrossRef] [PubMed]
  11. Deery, D.; Jimenez-Berni, J.; Jones, H.; Sirault, X.; Furbank, R. Proximal remote sensing buggies and potential applications for field-based phenotyping. Agronomy 2014, 4, 349–379. [Google Scholar] [CrossRef] [Green Version]
  12. Oliveira, M.D.M.; Varanda, C.M.R.; Félix, M.R.F. Induced resistance during the interaction pathogen x plant and the use of resistance inducers. Phytochem. Lett. 2016, 15, 152–158. [Google Scholar] [CrossRef] [Green Version]
  13. Silva-Perez, V.; Molero, G.; Serbin, S.P.; Condon, A.G.; Reynolds, M.P.; Furbank, R.T.; Evans, J.R. Hyperspectral reflectance as a tool to measure biochemical and physiological traits in wheat. J. Exp. Bot. 2018, 69, 483–496. [Google Scholar] [CrossRef] [Green Version]
  14. Almeida, R.P.P. Can Apulia’s olive trees be saved? Science. 2016, 353, 346–348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Schneider, K.; Van der Werf, W.; Cendoya, M.; Mourits, M.; Navas-Cortes, J.A.; Vicent, A.; Lansink, A.O. Impact of Xylella fastidiosa subspecies pauca in European olives. Proc. Natl. Acad. Sci. USA 2020, 117, 9250–9259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Riefolo, C.; Antelmi, I.; Castrignanò, A.; Ruggieri, S.; Galeone, C.; Belmonte, A.; Muolo, M.R.; Ranieri, N.A.; Labarile, R.; Gadaleta, G.; et al. Assessment of the hyperspectral data analysis as a tool to diagnose Xylella fastidiosa in the asymptomatic leaves of olive plants. Plants 2021, 10, 683. [Google Scholar] [CrossRef]
  17. Gitelson, A.; Merzlyak, M.N. Quantitative estimation of chlorophyll-a using reflectance spectra: Experiments with autumn chestnut and maple leaves. J. Photochem. Photobiol. B. Biol. 1994, 22, 247–252. [Google Scholar] [CrossRef]
  18. Gitelson, A.A.; Zur, Y.; Chivkunova, O.B.; Merzlyak, M.N. Assessing carotenoid content in plant leaves with reflectance spectroscopy. Photochem. Photobiol. 2002, 75, 272–281. [Google Scholar] [CrossRef]
  19. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef]
  20. Maimaitiyiming, M.; Ghulam, A.; Bozzolo, A.; Wilkins, J.L.; Kwasniewski, M.T. Early detection of plant physiological responses to different levels of water stress using reflectance spectroscopy. Remote Sens. 2017, 9, 745. [Google Scholar] [CrossRef] [Green Version]
  21. Zarco-Tejada, P.J.; Camino, C.; Beck, P.S.A.; Calderon, R.; Hornero, A.; Hernández-Clemente, R.; Kattenborn, T.; Montes-Borrego, M.; Susca, L.; Morelli, M.; et al. Previsual symptoms of Xylella fastidiosa infection revealed in spectral plant-trait alterations. Nat. Plants 2018, 4, 432–439. [Google Scholar] [CrossRef]
  22. Rouse, J.W.; Haas, R.; Schell, J.; Deering, D. Monitoring vegetation systems in the great plains with erts. In Proceedings of the Third Earth Resources Technology Satellite-1 Symposium, Washington, DC, USA, 10–14 December 1973; pp. 10–14. [Google Scholar]
  23. Guyot, G.; Baret, F.; Major, D. High spectral resolution: Determination of spectral shifts between the red and the near infrared. Int. Arch. Photogramm. Remote Sens. 1988, 11, 750–760. [Google Scholar]
  24. Filella, I.; Peñuelas, J. The red edge position and shape as indicators of plant chlorophyll content, biomass and hydric status. Int. J. Remote Sens. 1994, 15, 1459–1470. [Google Scholar] [CrossRef]
  25. Merton, R.; Huntington, J. Early simulation results of the ARIES-1 satellite sensor for multi-temporal vegetation research derived from AVIRIS. In Proceedings of the Eighth Annual JPL Airborne Earth Science Workshop, Pasadena, CA, USA, 8–14 February 1999; pp. 9–11. [Google Scholar]
  26. Huang, Z.; Turner, B.J.; Dury, S.J.; Wallis, I.R.; Foley, W.J. Estimating foliage nitrogen concentration from HYMAP data using continuum removal analysis. Remote Sens. Environ. 2004, 93, 18–29. [Google Scholar] [CrossRef]
  27. Fourty, T.; Baret, F.; Jacquemoud, S.; Schmuck, G.; Verdebout, J. Leaf optical properties with explicit description of its biochemical composition: Direct and inverse problems. Remote Sens. Environ. 1996, 56, 104–117. [Google Scholar] [CrossRef]
  28. Ceccato, P.; Flasse, S.; Tarantola, S.; Jacquemoud, S.; Grègoire, J.M. Detecting vegetation leaf water content using reflectance in the optical domain. Remote Sens. Environ. 2001, 77, 22–33. [Google Scholar] [CrossRef]
  29. Jacquemoud, S.; Ustin, S.L. Application of radiative transfer models to moisture content estimation and burned land mapping. In Proceedings of the 4th International Workshop on Remote Sensing and GIS Applications to Forest Fire Management, Ghent, Belgium, 5–7 June 2003; pp. 3–12. [Google Scholar]
  30. Marchantant, B.; Rudolph, S.; Roques, S.; Kindred, D.; Gillingham, V.; Welham, S.; Coleman, C.; Sylvester-Bradley, R. Establishing the precision and robustness of farmers’ crop experiments. Field Crops Res. 2019, 230, 31–45. [Google Scholar] [CrossRef]
  31. Brus, D.J.; De Gruijter, J.J. Random sampling or geostatistical modelling? Choosing between design-based and model-based sampling strategies for soil (with discussion). Geoderma 1997, 80, 1–44. [Google Scholar] [CrossRef]
  32. Rodrigues, M.S.; Corá, J.E.; Castrignanò, A.; Mueller, T.G.; Rienzi, E. A spatial and temporal prediction model of corn grain yield as a function of soil attributes. Agron. J. 2013, 105, 1878–1887. [Google Scholar] [CrossRef]
  33. Diggle, P.J.; Ribeiro, P.J. Model-Based Geostatistics; Springer: New York, NY, USA, 2007. [Google Scholar]
  34. Lark, R.M.; Cullis, B.R.; Welham, S.J. On spatial prediction of soil properties in the presence of a spatial trend: The empirical best linear unbiased predictor (E-BLUP) with REML. Eur. J. Soil Sci. 2006, 57, 787–799. [Google Scholar] [CrossRef]
  35. Cafarelli, B.; Castrignanò, A.; De Benedetto, D.; Palumbo, A.D.; Buttafuoco, G. A linear mixed effect (LME) model for soil water content estimation based on geophysical sensing: A comparison of a LME model and kriging external drift. Environ. Earth Sci. 2015, 73, 1951–1960. [Google Scholar] [CrossRef]
  36. Ferré, C.; Castrignanò, A.; Comolli, R. Comparison between spatial and non-spatial regression models for investigating tree–soil relationships in a polycyclic tree plantation of Northern Italy and implications for management. Agrofor. Syst. 2018, 93, 1–16. [Google Scholar] [CrossRef]
  37. Carter, G.A. Primary and secondary effects of water content on the spectral reflectance of leaves. Am. J. Bot. 1991, 78, 916–924. [Google Scholar] [CrossRef]
  38. Kokaly, R.F.; Clark, R.N. Spectroscopic determination of leaf biochemistry using band-depth analysis of absorption features and stepwise multiple linear regression. Remote Sens. Environ. 1999, 67, 267–287. [Google Scholar] [CrossRef]
  39. Klarzynski, O.; Plesse, B.; Joubert, J.M.; Yvin, J.C.; Kopp, M.; Kloareg, B.; Fritig, B. Linear β-1,3 Glucans are elicitors of defense responses in tobacco. Plant Physiol. 2000, 124, 1027–1038. [Google Scholar] [CrossRef] [Green Version]
  40. Boava, L.P.; Kuhn, O.J.; Pascholati, S.F.; Di Piero, R.M.; Furtado, E.L. Effect of acibenzolar-S-methyl and Saccharomyces cerevisiae on the activation of Eucalyptus defences against rust. Australas. Plant Pathol. 2009, 38, 594–602. [Google Scholar] [CrossRef]
  41. Guo, J.; Sun, K.; Zhang, Y.; Hu, K.; Zhao, X.; Liu, H.; Wu, S.; Hu, Y.; Zhang, Y.; Wang, Y. SlMAPK3, a key mitogen-activated protein kinase, regulates the resistance of cherry tomato fruit to Botrytis cinerea induced by yeast cell wall and β-glucan. Postharvest Biol. Technol. 2021, 171, 111350. [Google Scholar] [CrossRef]
  42. Lemaitre-Guillier, C.; Dufresne, C.; Chartier, A.; Cluzet, S.; Valls, J.; Jacquens, L.; Douillet, A.; Aveline, N.; Adrian, M.; Daire, X. VOCs are relevant biomarkers of elicitor-induced defences in grapevine. Molecules 2021, 26, 4258. [Google Scholar] [CrossRef] [PubMed]
  43. Peñuelas, J.; Piñol, J.; Ogaya, R.; Filella, I. Estimation of plant water concentration by the reflectance water index (R900/R970). Int. J. Remote Sens. 1997, 18, 2869–2875. [Google Scholar] [CrossRef]
  44. Gao, B.C. NDWI—a normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  45. Armstrong, G.A.; Hearst, J.E. Carotenoids 2: Genetics and molecular biology of carotenoid pigment biosynthesis. FASEB J. 1996, 10, 228–237. [Google Scholar] [CrossRef] [PubMed]
  46. Demmig-Adams, B.; Adams, W.W. The role of xanthophyll cycle carotenoids in the protection of photodynthesis. Trends Plant Sci. 2007, 1, 21–26. [Google Scholar] [CrossRef]
  47. Penuelas, J.; Filella, I.; Gamon, J.A. Assessment of photosynthetic radiation-use efficiency with spectral reflectance. New Phytol. 1995, 131, 291–296. [Google Scholar] [CrossRef]
  48. Schabenberger, O.; Gotway, C.A. Statistical Methods for Spatial Data Analysis; Taylor & Francis Inc.: Boca Raton, FL, USA, 2005. [Google Scholar] [CrossRef]
  49. Sun, H.; Liu, N.; Wu, L.; Chen, L.; Yang, L.; Li, M.; Zhang, Q. Water content detection of potato leaves based on hyper-spectral image. IFAC 2018, 51, 443–448. [Google Scholar] [CrossRef]
  50. Hennessy, A.; Clarke, K.; Lewis, M. Hyperspectral classification of plants: A review of waveband selection generalisability. Remote Sens. 2020, 12, 113. [Google Scholar] [CrossRef] [Green Version]
  51. Din, M.; Zheng, W.; Rashid, M.; Wang, S.; Shi, Z. Evaluating hyperspectral vegetation indices for leaf area in-dex estimation of Oryza sativa L. at diverse phenological stages. Front. Plant Sci. 2017, 8, 820. [Google Scholar] [CrossRef] [Green Version]
  52. Mananze, S.; Pôças, I.; Cunha, M. Retrieval of maize leaf area index using hyperspectral and multispectral data. Remote Sens. 2018, 10, 1942. [Google Scholar] [CrossRef] [Green Version]
  53. Li, H.; Yang, W.; Lei, J.; She, J.; Zhou, X. Estimation of leaf water content from hyperspectral data of different plant species by using three new spectral absorption indices. PLoS ONE 2021, 16, e0249351. [Google Scholar] [CrossRef]
  54. Feret, J.-B.; François, C.; Asner, G.P.; Gitelson, A.A.; Martin, R.E.; Bidel, L.P.R.; Ustin, S.L.; le Maire, G.; Jacquemoud, S. PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments. Remote Sens. Environ. 2008, 112, 3030–3043. [Google Scholar] [CrossRef]
  55. Curtis, M.J.; Claassen, V.P. Compost incorporation increases plant available water in a drastically disturbed serpentine soil. Soil Sci. 2005, 170, 939–953. [Google Scholar] [CrossRef]
  56. Nguyen, T.; Fuentes, S.; Marschner, P. Effects of compost on water availability and gas exchange in tomato during drought and recovery. Plant Soil Environ. 2012, 58, 495–502. [Google Scholar] [CrossRef] [Green Version]
  57. Qiu, Z.; Esan, E.O.; Ijenyo, M.; Gunupuru, L.R.; Asiedu, S.K.; Abbey, L. Photosynthetic activity and onion growth response to compost and Epsom salt. Int. J. Veg. Sci. 2020, 26, 535–546. [Google Scholar] [CrossRef]
  58. Hartz, T.K.; Costa, F.J.; Schrader, W.L. Suitability of composted green waste for horticultural uses. HortScience 1996, 31, 961–964. [Google Scholar] [CrossRef] [Green Version]
  59. Yuan, M.; Huang, Y.; Ge, W.; Jia, Z.; Song, S.; Zhang, L.; Huang, Y. Involvement of jasmonic acid, ethylene and salicylic acid signaling pathways behind the systemic resistance induced by Trichoderma longibrachiatum H9 in cucumber. BMC Genom. 2019, 20, 144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Alfiky, A.; Weisskopf, L. Deciphering Trichoderma–plant–pathogen interactions for better development of biocontrol applications. J. Fungi 2021, 7, 61. [Google Scholar] [CrossRef] [PubMed]
  61. Kumar, P.; Kumar, C. Molecular and enzymatic approach to study Trichoderma harzianum induced disease resistance in Brassica juncea against Albugo candida. J. Plant Dis. Prot. 2018, 125, 167–175. [Google Scholar] [CrossRef]
  62. Liu, S.Y.; Liao, C.K.; Lo, C.T.; Yang, H.H.; Lin, K.C.; Peng, K.C. Chrysophanol is involved in the biofertilization and biocontrol activities of Trichoderma. Physiol. Mol. Plant Pathol. 2016, 96, 1–7. [Google Scholar] [CrossRef]
  63. Narusaka, M.; Minami, T.; Iwabuchi, C.; Hamasaki, T.; Takasaki, S.; Kawamura, K.; Narusaka, Y. Yeast cell wall extract induces disease resistance against bacterial and fungal pathogens in Arabidopsis thaliana and Brassica crop. PLoS ONE 2015, 10, e0115864. [Google Scholar] [CrossRef]
  64. Aziz, A.; Poinssot, B.; Daire, X.; Adrian, M.; Bézier, A.; Lambert, B.; Joubert, J.M.; Pugin, A. Laminarin elicits defense responses in grapevine and induces protection against Botrytis cinerea and Plasmopara viticola. Mol. Plant Microbe Interact. 2003, 16, 1118–1128. [Google Scholar] [CrossRef] [Green Version]
  65. Gauthier, A.; Trouvelot, S.; Kelloniemi, J.; Frettinger, P.; Wendehenne, D.; Daire, X.; Joubert, J.M.; Ferrarini, A.; Delledonne, M.; Flors, V.; et al. Correction: The sulfated laminarin triggers a stress transcriptome before priming the SA- and ROS-dependent defenses during grapevine’s induced resistance against Plasmopara viticola. PLoS ONE 2018, 13, e0194327. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Gururaj, H.B.; Giridhar, P.; Ravishankar, G.A. Laminarin as a potential non-conventional elicitor for enhancement of capsaicinoid metabolites. Asian J. Plant Sci. Res. 2012, 2, 490–495. [Google Scholar]
  67. Wu, Y.R.; Lin, Y.C.; Chuang, H.W. Laminarin modulates the chloroplast antioxidant system to enhance abiotic stress tolerance partially through the regulation of the defensin-like gene expression. Plant Sci. 2016, 247, 83–92. [Google Scholar] [CrossRef] [PubMed]
  68. Manganiello, G.; Nicastro, N.; Caputo, M.; Zaccardelli, M.; Cardi, T.; Pane, C. Functional hyperspectral imaging by high-related vegetation indices to track the wide-spectrum Trichoderma biocontrol activity against soil-borne diseases of baby-leaf vegetables. Front. Plant Sci. 2021, 12, 1–21. [Google Scholar] [CrossRef] [PubMed]
  69. Shepherd, K.D.; Walsh, M.G. Development of reflectance spectral libraries for characterization of soil properties. Soil Sci. Soc. Am. J. 2002, 66, 988–998. [Google Scholar] [CrossRef]
  70. Martens, H.; Næs, T. Multivariate Calibration; John Wiley & Sons: New York, NY, USA, 1989. [Google Scholar]
  71. Savitzky, A.; Golay, M.J.E. Smoothing and differentiation of data by simplified least squares procedures. Anal. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  72. Næs, T.; Isakson, T.; Fearn, T.; Davies, T. A User-Friendly Guide to Multivariate Calibration and Classification; NIR Publications: Chichester, UK, 2004. [Google Scholar] [CrossRef]
  73. Colombo, C.; Palumbo, G.; Di Iorio, E.; Sellitto, V.M.; Comolli, R.; Stellacci, A.M.; Castrignanò, A. Soil organic carbon variation in alpine landscape (Northern Italy) as evaluated by diffuse reflectance spectroscopy. Soil Sci. Soc. Am. J. 2014, 78, 794–804. [Google Scholar] [CrossRef]
  74. Rossell, R.A.V. ParLeS: Software for chemometric analysis of spectroscopic data. Chemom. Intell. Lab. Syst. 2008, 90, 72–83. [Google Scholar] [CrossRef]
  75. Gorsuch, R.L. Factor Analysis; W. B. Saunders: Philadelphia, PA, USA, 1974. [Google Scholar]
  76. Cattell, R.B. The Scientific Use of Factor Analysis; Plenum: New York, NY, USA, 1978. [Google Scholar]
  77. Yong, A.G.; Pearce, S. A beginner’s guide to factor analysis: Focusing on exploratory factor analysis. Tutor. Quant. Methods Psychol. 2013, 9, 79–94. [Google Scholar] [CrossRef]
  78. Thode, H.C. Testing for Normality; Marcel Dekker: New York, NY, USA, 2002. [Google Scholar] [CrossRef]
  79. Blom, G. Statistical Estimates and Transformed Beta-Variables; John Wiley & Sons, Inc.: New York, NY, USA, 1958. [Google Scholar]
  80. Moran, P.A.P. Notes on continuous stochastic phenomena. Biometrika 1950, 37, 17–23. [Google Scholar] [CrossRef]
  81. Geary, R.C. Testing for normality. Biometrika 1947, 34, 209–242. [Google Scholar] [CrossRef] [Green Version]
  82. Goovaerts, P. Geostatistics for Natural Resources Evaluation; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
  83. Wackernagel, H. Multivariate Geostatistics: An Introduction with Applications, 3rd ed.; Springer: Berlin, Germany, 2003. [Google Scholar]
  84. Webster, R.; Oliver, M.A. Geostatistics for Environmental Scientists, 2nd ed.; Wiley: Chichester, UK, 2007. [Google Scholar]
  85. Levene, H. Robust tests for equality of variances. In Contributions to Probability and Statistics: Essays in Honor of Harold Hoteling; Olkin, I., Hotelling, H., Eds.; Stanford University Press: Redwood City, CA, USA, 1960; pp. 278–292. [Google Scholar]
  86. Littell, R.C.; Milliken, G.A.; Stroup, W.W.; Wolfinger, R.D.; Schabenberger, O. SAS System for Mixed Models; SAS Institute: Cary, NC, USA, 1996; Volume 633, pp. 130–134. [Google Scholar]
Figure 1. Graph of the loadings of the first PC1 (blue line), the second PC2 (red line), and the third PC3 (green line) >80. The dashed line represents the chosen threshold for loading values.
Figure 1. Graph of the loadings of the first PC1 (blue line), the second PC2 (red line), and the third PC3 (green line) >80. The dashed line represents the chosen threshold for loading values.
Plants 10 02575 g001
Figure 2. Experimental design of the two-way split-plot field trial carried out under a multi-tunnel greenhouse located at the CREA-Research Center for Vegetable and Ornamental Crops.
Figure 2. Experimental design of the two-way split-plot field trial carried out under a multi-tunnel greenhouse located at the CREA-Research Center for Vegetable and Ornamental Crops.
Plants 10 02575 g002
Figure 3. Experimental field trial consisting of a baby-leaf wild rocket cultivation carried out under greenhouse located at the CREA-Research Center for Vegetable and Ornamental Crops.
Figure 3. Experimental field trial consisting of a baby-leaf wild rocket cultivation carried out under greenhouse located at the CREA-Research Center for Vegetable and Ornamental Crops.
Plants 10 02575 g003
Table 1. Basic statistics of the PCs.
Table 1. Basic statistics of the PCs.
VariablesMeanMedianStd DeviationSkewnessKurtosis
PC1−0.0455−0.05940.6970.3580.230
PC2−0.06410.08230.769−1.1832.185
PC30.0110.00030.8050.2970.671
Table 2. Normality tests of the PCs.
Table 2. Normality tests of the PCs.
TestStatisticPC1PC2PC3p-ValuePC1PC2PC3
Kolmogorov-SmirnovD0.0300.1000.032p > D>0.1500<0.0100>0.1500
Cramer-von MisesW-Qu0.0951.5570.048p > W-Qu0.1352<0.0050>0.2500
Anderson-DarlingA-Qu0.7018.8870.393p > A-Qu0.0707<0.0050>0.2500
Table 3. Levene’s test of variance homogeneity for PC1, Gaussian transformed PC2 (rPC2) and PC3.
Table 3. Levene’s test of variance homogeneity for PC1, Gaussian transformed PC2 (rPC2) and PC3.
EffectsDFPC1rPC2PC3
F Valuep > FF Valuep > FF Valuep > F
Compost15.140.02380.470.49190.730.3931
Treatment44.070.00301.080.36481.140.3352
Table 4. Tests of spatial autocorrelation of the three PCs.
Table 4. Tests of spatial autocorrelation of the three PCs.
VariableCoefficientObservedExpectedStd DevZp > |Z|
PC1Moran’s I0.094−0.0020.0253.850.0001
Geary’s c0.9021.000.026−3.730.0002
rPC2Moran’s I0.065−0.0020.0252.690.0071
Geary’s c0.9321.000.026−2.620.0089
PC3Moran’s I0.050−0.0020.0252.080.0373
Geary’s c0.9411.000.026−2.230.0256
Table 5. LSE estimates of the variogram parameters for PC1, rPC2, and PC3.
Table 5. LSE estimates of the variogram parameters for PC1, rPC2, and PC3.
PC1rPC2PC3
ParameterValuet Valuep > |t|Valuet Valuep > |t|Valuet Valuep > |t|
Nugget0.37654.78<0.00010.917250.62<0.00010.55690.35<0.0001
Partial sill0.0537.12<0.0001---0.0222.830.0115
Range (m)6.4111.98<0.0001- -11.196.63<0.0001
Table 6. Results of the mixed effect model estimation for the PCs.
Table 6. Results of the mixed effect model estimation for the PCs.
a. Random Effect
VariableVariogram Model ParametersEffectRML EstimateStandard Errorz Valuep > Z
PC1NuggetCompost NO0.2590.0594.42<0.0001
Partial sillCompost YES0.1000.0422.400.0081
Range (m)Compost YES2.770.8403.300.0005
NuggetCompost YES0.2560.0347.62<0.0001
rPC2Partial sill 0.1140.0621.840.0332
Range (m) 7.253.032.390.0084
Nugget 0.8980.06414.09<0.0001
PC3Partial sill 0.0620.0421.480.0689
Range (m) 11.186.681.670.0470
Nugget 0.5570.03814.46<0.0001
b. Fixed Effect
PC1rPC2PC3
EffectF Valuep > FF Valuep > FF Valuep > F
Compost13.400.00150.040.85076.330.0185
Treatment3.170.01553.390.01064.920.0008
Compost × Treatment4.420.00210.990.41453.870.0050
Table 7. The most relevant significative LSE differences between the levels of each effect and the ones of interaction.
Table 7. The most relevant significative LSE differences between the levels of each effect and the ones of interaction.
VariablesEffectCompostTreatmentCompostTreatmentEstimatesStandard Errort ValuePr > |t|
PC1CompostNO YES 0.3050.0833.660.0015
Treatment CER CTR0.2210.1111.990.0486
Treatment TRI CTR0.2980.1102.710.0076
Treatment CHE CTR0.3760.1113.370.0010
Compost × TreatmentNOLAMYESTRI0.48030.16302.950.0041
Compost × TreatmentYESLAMNOTRI−0.64510.1630−3.960.0002
Compost × TreatmentNOCHEYESCTR0.3960.1632.430.0170
Compost × TreatmentYESCHENOCTR0.3560.1632.180.0318
Compost × TreatmentYESCHEYESCTR0.5480.1683.270.0019
rPC2Treatment CER CHE−0.6040.166−3.640.0003
Treatment LAM CHE−0.3570.160−2.240.0262
Treatment TRI CHE−0.3320.157−2.120.0356
Treatment CHE CTR0.3460.1602.170.0321
PC3CompostNO YES −0.2590.103−2.520.0185
Treatment CER CHE−0.3670.126−2.920.0040
Treatment LAM CHE−0.3510.121−2.900.0041
Treatment TRI CHE−0.5000.119−4.18<0.0001
Treatment TRI CTR0.2640.1222.170.0313
Compost × TreatmentYESLAMNOTRI0.55320.18423.000.0034
Compost × TreatmentYESLAMYESTRI0.52480.17602.980.0033
Compost × TreatmentYESCHENOCTR0.3650.1822.010.0478
Compost × TreatmentYESCHEYESCTR0.4230.1692.510.0131
Table 8. Results of the mixed effect model estimation for the PCs.
Table 8. Results of the mixed effect model estimation for the PCs.
PC1rPC2PC3
EffectF Valuep > FF Valuep > FF Valuep > F
Compost25.17<0.00010.530.468618.64<0.0001
Treatment4.910.00072.410.04845.180.0004
Compost × Treatment6.34<0.00011.640.16404.130.0027
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pane, C.; Galieni, A.; Riefolo, C.; Nicastro, N.; Castrignanò, A. Hyperspectral Reflectance Response of Wild Rocket (Diplotaxis tenuifolia) Baby-Leaf to Bio-Based Disease Resistance Inducers Using a Linear Mixed Effect Model. Plants 2021, 10, 2575. https://doi.org/10.3390/plants10122575

AMA Style

Pane C, Galieni A, Riefolo C, Nicastro N, Castrignanò A. Hyperspectral Reflectance Response of Wild Rocket (Diplotaxis tenuifolia) Baby-Leaf to Bio-Based Disease Resistance Inducers Using a Linear Mixed Effect Model. Plants. 2021; 10(12):2575. https://doi.org/10.3390/plants10122575

Chicago/Turabian Style

Pane, Catello, Angelica Galieni, Carmela Riefolo, Nicola Nicastro, and Annamaria Castrignanò. 2021. "Hyperspectral Reflectance Response of Wild Rocket (Diplotaxis tenuifolia) Baby-Leaf to Bio-Based Disease Resistance Inducers Using a Linear Mixed Effect Model" Plants 10, no. 12: 2575. https://doi.org/10.3390/plants10122575

APA Style

Pane, C., Galieni, A., Riefolo, C., Nicastro, N., & Castrignanò, A. (2021). Hyperspectral Reflectance Response of Wild Rocket (Diplotaxis tenuifolia) Baby-Leaf to Bio-Based Disease Resistance Inducers Using a Linear Mixed Effect Model. Plants, 10(12), 2575. https://doi.org/10.3390/plants10122575

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