Impact Factor 5.753 | CiteScore 8.2
More on impact ›

ORIGINAL RESEARCH article

Front. Plant Sci., 17 August 2018 | https://doi.org/10.3389/fpls.2018.01195

Hyperspectral Canopy Sensing of Wheat Septoria Tritici Blotch Disease

Kang Yu1*, Jonas Anderegg1, Alexey Mikaberidze2, Petteri Karisto2, Fabio Mascher3, Bruce A. McDonald2, Achim Walter1 and Andreas Hund1
  • 1Crop Science Group, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
  • 2Plant Pathology Group, Institute of Integrative Biology, ETH Zurich, Zurich, Switzerland
  • 3Plant Breeding and Genetic Resources, Strategic Research Division Plant Breeding, Agroscope, Nyon, Switzerland

Producing quantitative and reliable measures of crop disease is essential for resistance breeding, but is challenging and time consuming using traditional phenotyping methods. Hyperspectral remote sensing has shown potential for the detection of plant diseases, but its utility for phenotyping large and diverse populations of plants under field conditions requires further evaluation. In this study, we collected canopy hyperspectral data from 335 wheat varieties using a spectroradiometer, and we investigated the use of canopy reflectance for detecting the Septoria tritici blotch (STB) disease and for quantifying the severity of infection. Canopy- and leaf-level infection metrics of STB based on traditional visual assessments and automated analyses of leaf images were used as ground truth data. Results showed (i) that canopy reflectance and the selected spectral indices show promise for quantifying STB infections, and (ii) that the normalized difference water index (NDWI) showed the best performance in detecting STB compared to other spectral indices. Moreover, partial least squares (PLS) regression models allowed for an improvement in the prediction of STB metrics. The PLS discriminant analysis (PLSDA) model calibrated based on the spectral data of four reference varieties was able to discriminate between the diseased and healthy canopies among the 335 varieties with an accuracy of 93% (Kappa = 0.60). Finally, the PLSDA model predictions allowed for the identification of wheat genotypes that are potentially more susceptible to STB, which was confirmed by the STB visual assessment. This study demonstrates the great potential of using canopy hyperspectral remote sensing to improve foliar disease assessment and to facilitate plant breeding for disease resistance.

Introduction

Septoria tritici blotch (STB), caused by a fungal pathogen Zymoseptoria tritici, is one of the most damaging foliar diseases in wheat, causing crop yield losses and resulting in higher production costs due to fungicide applications (Fones and Gurr, 2015; Torriani et al., 2015). It will be difficult to meet the growing global demands for wheat whilst reducing the environmental costs associated with wheat production without better control of STB based on breeding varieties with improved STB resistance (O'Driscoll et al., 2014). Remote sensing and mapping techniques have been suggested as new tools to improve control of crop diseases (Pinter et al., 2003; Oerke et al., 2014). This idea is based on approaches of precision farming, where site-specific management of fertilizer and fungicide applications are performed based on input from remote sensing and geographic information systems (Pinter et al., 2003; Steddom et al., 2005; Mattupalli et al., 2015).

Advanced remote sensing techniques enhance the possibility for detecting crop diseases in a rapid and non-invasive way (Mahlein, 2016; Wakie et al., 2016). Yet, there remains a lack of investigations and/or applications of these technologies to advance breeding by improving the efficiency of selection in large, genetically diverse plant populations under field conditions (Deery et al., 2014; Walter et al., 2015). Advances in breeding for disease resistance are often limited by a lack of available resistant varieties and can also be constrained by deficiencies associated with conventional phenotyping methods—for instance, the precision, reproducibility, and throughput of conventional visual assessments (Steddom et al., 2005; O'Driscoll et al., 2014; Mahlein, 2016). To better meet the specific needs of crop breeders, remote, and proximal sensing techniques will need to generate accurate plant phenotypes in genetically diverse populations under field conditions (Furbank and Tester, 2011; Kirchgessner et al., 2017; Yang et al., 2017).

Hyperspectral remote sensing has been shown to be a promising tool for monitoring foliar diseases over scales ranging from individual leaves to fields to entire regions (Mahlein et al., 2013; Ashourloo et al., 2014b; Wakie et al., 2016; Yuan et al., 2016). For instance, hyperspectral reflectance measurements of wheat leaves were successfully used to detect leaf rust and estimate disease severity (Ashourloo et al., 2014a,b). Hyperspectral imaging was successfully used for the detection of Fusarium head blight in wheat grains (Bauriegel and Herppich, 2014). For field crops, foliar diseases often occur simultaneously at different leaf layers and sometimes infections begin in the lower canopy. Canopy-level hyperspectral signatures responding to certain disease symptoms can be complicated by crop morphology and canopy structure (Zarco-Tejada et al., 2001; Yu et al., 2014; Mahlein, 2016). Therefore, selection of spectral features that can qualitatively and quantitatively measure a disease irrespective of the different canopy structures likely to exist across a large and diverse population is critical in order for hyperspectral remote sensing to be useful for a plant breeding program.

Some leaf diseases can cause losses of chlorophyll and other diseases may trigger leaf water deficits (Delalieux et al., 2009; Yu et al., 2014; Oerke et al., 2016). As a result, infected plants may present different characteristics at the chlorophyll and/or water absorption bands of their spectral signatures compared to healthy plants. Canopy-level sensing methods and spectral features should be compared for their ability to detect these differing signatures of infection under field conditions. It will be particularly useful to better understand the capabilities of different spectral features to detect and quantify common crop diseases in large and genetically diverse breeding populations that possess diverse canopy structures.

Disease detection can be based on vegetation indices that use approximation methods to quantify leaf biochemical changes, e.g., changes in leaf pigments and water content (Bürling et al., 2011; Ashourloo et al., 2014a; Yu et al., 2015). Hyperspectral narrow vegetation indices enhance the power to detect plant diseases by using high resolution narrow bands that enable characterization of subtle differences in plants and canopies (Delalieux et al., 2009; Mahlein et al., 2013). For instance, hyperspectral disease severity indices have been developed specifically for the detection of wheat leaf rust (Ashourloo et al., 2014a). Different plant diseases can also produce disease-specific spectral signatures that allow for the discrimination of different diseases (Mahlein et al., 2013). For disease resistance breeding, an accurate measure of disease infection across a large and genetically diverse population is critical. Therefore, it is important to know which spectral indices have the best capability to provide an accurate measure of a crop disease. Comparisons of spectral features and spectral indices for crop diseases will also advance sensor development for field-based plant phenotyping and precision disease management.

The objective of this study was to investigate the potential of canopy-level hyperspectral measurements to detect and accurately measure STB infection in a genetically diverse population composed of 335 elite winter wheat varieties, and to determine suitable hyperspectral vegetation indices and the most sensitive spectral features for STB detection and quantification under field conditions.

Materials and Methods

Plant Materials and Experimental Design

Main Experiment

A field experiment was conducted during the 2015/2016 winter wheat growing season. The experiment was conducted on a 1 ha field site of the ETH Field Experiment Station Lindau-Eschikon, Switzerland (Kirchgessner et al., 2017). A total of 335 European wheat varieties selected from the GABI-wheat panel (Kollers et al., 2013) were grown in two blocks (replications). Within each block, each variety was grown in a randomly located plot (1.2 × 1.7 m), except for the check plots (grown as controls for field variability) that were planted in 21 evenly distributed locations at which the wheat variety CH CLARO was planted. The sowing and harvest dates were 13 October 2015 and 28 July 2016, respectively.

Fertilization and fungicide applications in the main experiment followed the local recommendations. Fungicide applications were made three times during the growing season with applications on 6 April (BBCH 31; Input® containing prothioconazol and spiroxamine, Bayer CropScience Ltd, Germany), 25 May (BBCH 51; Aviator Xpro containing bixafen and prothioconazole, Bayer CropScience Ltd, Germany) and 6 June (BBCH 65, Osiris® containing epoxiconazol and metconazol, BASF, Germany). Details of fertilizer and fungicide applications were reported in Karisto et al. (2018). In this study, for the purpose of binary classification, the main experimental plots treated with fungicides were considered healthy (Z. tritici un-treated canopies and referred to as the treatment T0, see next section for T1), though some of these plots were later shown to develop significant levels of STB (Karisto et al., 2018). All STB infections in the T0 plots were natural, a result of the ascospores that immigrated into the experimental plots from wheat fields in the same region. The precipitation between April and end of July was 573 mm, 36% above the reference period 1981–2010 (MeteoSchweiz, 2016). This high precipitation was conducive to development of a strong STB epidemic.

Septoria Tritici Blotch (STB) Reference Experiment

The STB reference experiment was located next to the blocks of the main experiment. Four varieties (CH CLARO, DRIFTER, RUNAL, and TITLIS), included also in the set of 335 varieties of T0, were used for the disease treatment (Z. tritici treated, hereafter referred to as T1). As shown in Karisto et al. (2018), DRIFTER is very susceptible, CH CLARO and RUNAL have intermediate resistance, and TITLIS is resistant to STB. Specifically, DRIFTER ranked 2/335 (the 2nd most susceptible variety), CH CLARO ranked 108/335 (moderately susceptible/resistant), RUNAL ranked 151/335 (moderately susceptible/resistant), and TITLIS ranked 305/335 (resistant). The four STB reference varieties were grown in four 2.45 × 3 m plots with two replications, in the absence of fungicide protection. Four Swiss strains of Z. tritici were inoculated into these plots on 10 May, 2016 in central areas (0.5 × 0.5 m) of each plot by spraying the whole canopy until runoff using 60 ml of a spore suspension containing 106 spores/ml of each strain on each plot. In addition to the STB infection caused by the four inoculated strains, we observed significant STB infection that developed by natural infection from neighboring wheat fields, but we found practically no other diseases (Karisto et al., 2018).

Canopy Spectral Measurements

Hyperspectral reflectance measurement was performed across different growth stages (Table 1) over wheat plots using an ASD FieldSpec® 4 spectroradiometor (ASD Inc.) that is equipped with an optical fiber that uses a field of view of 25°. The fiber was mounted on a tripod arm that allows for measuring from a nadir view at a height of 1.5 m above the ground. Canopy reflectance was collected at 1–2 locations for the main experiment plots (T0) and at 3–5 locations for the larger STB disease reference plots (T1), and for each location the spectrum was determined based on the average of 15–25 spectral records. A Spectralon® white reference panel was used for calibration before measuring canopy reflectance, and the calibration was repeated every 10–15 min or immediately following significant changes in light conditions.

TABLE 1
www.frontiersin.org

Table 1. Measurements for canopy spectra, visual assessment (VA), and leaf sample collections (C) for STB metrics.

STB Disease Assessment

Leaf Level Imagery Assessment

Leaves exhibiting obvious STB lesions were sampled from each cultivar on 20 May 2016 (Collection 1, see Table 1) and 4 July 2016 (Collection 2) (Karisto et al., 2018). In each plot of the main experiment, 16 infected leaves were collected, transported to the laboratory in cooled boxes, imaged with flatbed scanners, and analyzed using the method described in Stewart et al. (2016) and Karisto et al. (2018). In this study, the percentage of leaf area covered by lesions (PLACL), the density of pycnidia per unit leaf area (ρ-leaf) and the density of pycnidia per unit lesion area (ρ-lesion) were used as quantitative measures of STB leaf infections.

Canopy Level Visual Assessment

Visual assessment of STB symptoms was performed on 20 May, 21 June and 29 June 2016 (Table 1). STB severity in each plot was scored from 1 to 9 (1 and 9 denote no disease and complete infection, respectively), following a logistic (non-linear) progression based on the symptoms of the top three leaf layers (Michel, 2001). This assessment is the standard procedure for identifying STB-resistant germplasm to use for resistance breeding. In our study, the scores ranged only between 1 and 5 in the main experiment as a result of the fungicide treatments; meaning that only mild to moderate STB infections were observed. Based on visual assessments (VA) of STB at three time points (Table 1), the area under disease progress curve (AUDPC) was calculated for each variety (Karisto et al., 2018).

Data Analysis

Canopy reflectance values for individual plots were averaged and served as input data for the calculation of hyperspectral narrow-band indices based on the formulae listed in Table 2. Binary logistic regression (BLR) was employed for the evaluation of each spectral index for the detection of plant diseases at different time points of the STB epidemic. The c-statistic of spectral indices based BLR models can be used as a measure to assess the discriminatory performance of spectral indices (Delalieux et al., 2009; Yu et al., 2014). We calculated the c-statistic for each spectral index in detecting the event of interest (disease) based on the four reference varieties. The c-statistic value is equivalent to the area under the receiver-operating-characteristic (ROC) curve. We used a general rule that considers: 0.7 ≤ c < 0.8 as acceptable discrimination; 0.8 ≤ c < 0.9 as excellent discrimination; and c ≥ 0.9 as outstanding discrimination (Hosmer and Lemeshow, 2005).

TABLE 2
www.frontiersin.org

Table 2. Hyperspectral vegetation indices (dimensionless) used in this study for wheat STB disease detection.

Spearman's correlation analysis was applied to study the relationships between the spectral indices and STB measurements to assess their monotonic relationships (whether linear or not). We used partial least squares (PLS) regression (PLSR) to characterize the STB infection severity in the full spectrum of 350–2,500 nm, as well as partial least squares discriminant analysis (PLSDA) for the differentiation between diseased and healthy plots. PLS is a multivariate method which relates two matrices, X and Y, i.e., explanatory and response matrices, by extracting a number of components (also known as latent variables) to model the variations of both matrices (Wold et al., 2001). The variable importance in projections (VIP) of PLS models were computed to evaluate the importance of individual bands in modeling STB metrics, and VIP ≥ 0.8 was considered as the threshold for important contributions (Wold et al., 2001; Eriksson et al., 2013). The optimal number of components was determined by finding the lowest root mean squared error of prediction (RMSEP) and predicted residual sum of squares PRESS (van der Voet, 1994).

Spectral data for the four STB disease reference varieties (CH CLARO, DRIFTER, RUNAL, and TITLIS) were used to calibrate the BLR and PLSDA models. The remaining 331 varieties were used as independent validation data to test the classification performance. The sampling periods in which STB metrics were best correlated (Week 3, Table 1) with canopy spectra were used for the calibration and comparison of classification models. The validation dataset also includes data of the four reference varieties in Week 2, which is used to test the capability of early detection of STB infections.

Classification performance is illustrated in a confusion matrix with the following measures: Sensitivity=TP(TP+FN) (true positive rate, TPR), Specificity=TN(TN+FP) (true negative rate, TNR), Accuracy=(TP+TN)(TP+TN+FP+FN), where T, F, P, and N denote true, false, positive and negative, respectively, and the Kappa statistic. The Kappa statistic measures the agreement (1 represents perfect agreement) between the actual (the assigned treatments) and model predicted classes, and it accounts for chance agreement (McHugh, 2012). Hence, it is a balanced measure for unbalanced classification problems as in this study where we have a much smaller sample size for the diseased treatment.

Kappa=pa-pe1-pe    (1)
pe=(TN+FP)×(TN+FN)+(FN+TP)×(FP+TP)(TP+TN+FP+FN)2    (2)

where pa is the actual agreement (identical to accuracy), whereas pe is the expected agreement by chance (random accuracy) that can be calculated as Equation (2). Letters T and F used in the equation terms denote true and false, respectively, and P and N denote positive and negative, respectively.

All analyses were performed in the R environment (R Core Team, 2016). PLSR was implemented using the r package “pls” (Mevik and Wehrens, 2007). PLSDA was implemented using the package “mixOmics” (Rohart et al., 2017).

Results and Discussion

Canopy Spectral Characteristics of STB Infection

At the early stage of the STB infection, only marginal differences were visible between the diseased and healthy canopies (Figures 1A,I) and there were no distinct differences in the canopy reflectance (Figure 1E). As the epidemic progressed, STB symptoms became obvious in the diseased canopies (Figures 1BD,JL) and differences in canopy reflectance between diseased and healthy plots of the same variety became more pronounced (Figures 1FH).

FIGURE 1
www.frontiersin.org

Figure 1. Comparisons of wheat canopy RGB images and hyperspectral reflectance curves between the Zymoseptoria tritici un-treated (T0, images I–L) and treated (T1, images A–D) plots at four periods. (A,E,I) Week 0 was at the end of May 2016; (B,F,J) Week 1 was at mid-June 2016; (C,G,K) Week 2 was at the beginning of July; and (D,H,L) Week 3 was at mid-July 2016 (see also Table 1 for a full list of sampling periods and corresponding growth stages).

By the late stage of the STB epidemic, canopy reflectance in the green region was much lower for the diseased plots than for the healthy plots (Figure 1H), whereas reflectance in the red region was higher in diseased compared to healthy plots (Figures 1G,H, 2A). Reflectance in the green and red regions is associated mainly with green canopy fraction and leaf pigments (Ustin et al., 2009). Plant diseases often cause chlorotic and necrotic leaf surfaces that lack chlorophyll and thus affect leaf and canopy reflectance and transmittance in the red and green regions (Oerke et al., 2016). Therefore, the reflectance changes in the visible region observed in the diseased canopies suggest that STB infections cause reductions in leaf chlorophyll and corresponding reductions in photosynthetic activity (Mahlein, 2016; Oerke et al., 2016).

FIGURE 2
www.frontiersin.org

Figure 2. (A) Canopy hyperspectral reflectance and reflectance means of the Zymoseptoria tritici treated wheat plots (T1) in comparison with the Z. tritici untreated plots (T0) of the four reference varieties (CH CLARO, DRIFTER, RUNAL and TITLIS) at the stage of significant infection (Week 3). (B) Variable importance in projections (VIP, dimensionless) of the first three components (comp 1, 2, and 3) of PLSDA model in characterizing the difference of the two treatments (VIP ≥ 0.8 is considered as the threshold of influential contribution). Gray and blue bars indicate the spectral bands of NDVI and NDWI respectively (see formulae in Table 2).

Reflectance in the NIR (750–1,300 nm) region was distinctly lower in the diseased canopies than in the healthy canopies at week 3 (Figure 2A). Reflectance in the NIR is governed mainly by the internal leaf structure as well as by the structure of the whole canopy. It appears that STB infection caused damage to the internal and surface structures of leaves that increased water loss, resulting in a significant decrease in NIR reflectance in the diseased canopies (Figure 2A). The STB-induced canopy spectral signature also showed typical vegetation stress characteristics in the shortwave infrared (SWIR, 1,300–2,500 nm) regions, where reflectance in the diseased plots was higher than in the healthy plots (Figures 1G,H, 2A). Leaf reflectance in the SWIR region is associated with water absorption as well as proteins and carbohydrates (Kumar et al., 2001; Thenkabail et al., 2011), suggesting that STB infections decreased the water status of the diseased canopy leaves leading to low absorption and high SWIR reflectance. However, the observed canopy spectral responses in the chlorophyll and water absorption regions might be similar to those caused by other diseases or stresses. The changing spectral characteristics over time associated with the development of particular diseases or stresses may make it possible to identify the particular causes of stress by measuring canopy spectra across several time points.

VIP scores provide a measure of the importance of variables used in PLS modeling, and therefore the VIP scores over the full spectrum indicate the contribution of spectral bands in characterizing the response variables, i.e., the STB infection metrics used in this study. Figure 2B highlights the predominant contribution (VIP ≥ 0.8 is considered as influential contribution) of the NIR region and water absorption bands in differentiating between the healthy and diseased canopies. Healthy canopies yielded low reflectance and high absorption in the water absorption bands in the NIR and SWIR regions. This is attributed to the fact that healthy leaves have more compact leaf mesophyll layers and fewer air-water interfaces than diseased leaves, resulting in decreased reflectance in the NIR region (Kumar et al., 2001).

Spectral Indices for STB Detection

Correlation Between Canopy STB Metrics and Spectral Indices

Spearman's correlation was measured between canopy STB metrics (STB scores and AUDPC) and spectral indices (Table 3, Figure S1 and Figure 3). STB scores correlated positively with PSRI and SIPI, whereas they were negatively related to other spectral indices. The correlation coefficients for the third assessment (VA3) were larger than for the two earlier assessments, consistent with the continuous increase in STB infections. As expected, AUDPC yielded higher correlation with the spectral indices than the STB scores for individual assessments. STB scores and AUDPC both resulted in positive correlations with PSRI and negative correlations with NDWI (Figure S1 and Figure 3), suggesting associations between the amount of STB infection and the degree of leaf senescence and water status. Other studies reported that diseased leaves are prone to suffer severe water deficits as a result of stomatal opening disorders (Oerke et al., 2016), and water deficit has been shown to play a key role in accelerating leaf senescence (Rivero et al., 2007). We hypothesize that this explains the observed positive correlations between PSRI and the STB visual assessments.

TABLE 3
www.frontiersin.org

Table 3. Correlation coefficients of the Spearman's correlations between the spectral indices and the two canopy measures (visually assessed STB scores and corresponding AUDPC values) at three visual assessments (VA1, VA2, and VA3) and the three leaf STB metrics (PLACL, ρ-leaf, and ρ-lesion) at two collections (C1 and C2) of leaf samples.

FIGURE 3
www.frontiersin.org

Figure 3. Area under disease progress curve (AUDPC) of visual scores plotted as a function of spectral indices. Correlation coefficients (r) are for the Spearman's correlations between the spectral indices and AUDPC of STB scores (***p < 0.001). Loess fit lines are added to the scatter plot to show the trend of the relationship. Results are based on the data of the whole population (335 varieties) used in the main experiment.

STB scores and AUDPC were negatively correlated to the greenness- or chlorophyll-related spectral indices such as NDVI, ND705, mSR705, and CAI, in agreement with studies that reported negative correlations between plant disease severity and NDVI and visible band spectral indices (Steddom et al., 2005; Ashourloo et al., 2014a; Jansen et al., 2014). In contrast, SIPI was positively related to the STB scores and AUDPC despite using the visible bands related to plant pigments. SIPI estimates the pigment ratio between carotenoid and chlorophyll a (Car/Chla) and its increase is often related to stresses and subsequent changes in physiological and phenological status (Peñuelas et al., 1995b). Similarly, positive correlations have also been reported between the SIPI and visual rating scores of soybean diseases (Bajwa et al., 2017). Therefore, SIPI shows the potential of being a surrogate of plant disease severity under field conditions.

Correlation Between Leaf STB Metrics and Spectral Indices

Table 3 shows the correlation coefficients for the relationships between the three leaf-based STB metrics (PLACL, ρ-leaf, and ρ-lesion) and the spectral indices for the two collections. The correlation coefficients for the second collection were higher than for the first collection, which is in agreement with the results of the visual assessments and confirms the continuous increase in STB infections during the epidemic. Compared to the first collection, PLACL and ρ-leaf doubled their correlation coefficients in the second collection. In contrast, the correlation for ρ-lesion increased only marginally compared to the first collection. These results indicate that the canopy level spectra might be of limited value for detecting subtle differences associated with the number of pycnidia occurring within lesions.

The relationships between the spectral indices and leaf STB metrics for the second collection are illustrated in Figures 4, 5 and Figure S2. The spectral indices correlating negatively with PLACL are mainly related to vegetation fraction (e.g., NDVI), water content (NDWI), and chlorophyll or photosynthetic capacity (e.g., ND705, PRI). In contrast, PSRI and SIPI were positively related to the PLACL (Figure 4). Similar to the canopy STB metrics, SIPI was positively related to PLACL, ρ-leaf, and ρ-lesion, suggesting that the physiological and/or phenological changes in diseased leaves can be captured in hyperspectral signatures at both leaf and canopy levels (Peñuelas et al., 1995b; Yang et al., 2007; Bajwa et al., 2017). Here, increases in SIPI and PSRI were associated with increased degrees of lesion development on leaves, which resulted in decreased photosynthetic activities and also enhanced the senescence-related signatures that were observed from the canopy.

FIGURE 4
www.frontiersin.org

Figure 4. STB percentage of leaf area covered by lesions (PLACL) plotted as a function of spectral indices for the second collection. Correlation coefficients (r) are for the Spearman's correlations between the spectral indices and PLACL (***p < 0.001). Loess fit lines are added to the scatter plot to show the trend of the relationship. Results are based on the data of the whole population (335 varieties) used in the main experiment.

FIGURE 5
www.frontiersin.org

Figure 5. STB density of pycnidia per unit leaf area (ρleaf) plotted as a function of spectral indices for the second collection. Correlation coefficients (r) are for the Spearman's correlations between the spectral indices and ρleaf (***p < 0.001). Loess fit lines are added to the scatter plot to show the trend of the relationship. Results are based on the data of the whole population (335 varieties) used in the main experiment.

For PLACL and ρ-leaf, PSRI and SIPI yielded the same magnitude of correlations (Figures 4, 5). However, the correlations were significantly weaker for the relationships with the ρ-lesion (Figure S2). PLACL is the best indicator of STB damage to the host plant because it directly measures the loss in photosynthetic area due to infection, whereas ρ-lesion is the best indicator of pathogen reproduction on the host plant because it is a direct measure of the ability of the pathogen to convert infected plant tissue into pathogen reproductive structures. The combined measure ρ-leaf provides the best overall measurement of relative resistance to STB (Karisto et al., 2018). Our results suggest that canopy-level hyperspectral indices may be less useful for identifying resistance that affects pathogen reproduction (i.e., ρ-lesion), but could be very useful for identifying resistance that affects host damage (i.e., PLACL). Further improvements may be possible by applying high spatial resolution hyperspectral imaging sensors (Mahlein, 2016; Oerke et al., 2016).

Overall, consistent trends (either positive or negative relationships) were observed in the correlations between canopy spectral indices and both canopy- and leaf-level STB metrics. Spectral indices related to chlorophyll and water status were always negatively related to STB metrics, whereas the senescence-related indices were always positively related. Our results suggest that canopy hyperspectral measurements can provide a quantitative measure of STB damage on wheat and provide a useful complement to traditional visual disease ratings.

Discriminant Analysis Using Spectral Indices

Spectral indices are also used to discriminate stressed or diseased plants from healthy ones. Thus, understanding which spectral indices are suitable for identifying STB resistance is critical. BLR discriminant analysis shows the performances of spectral indices in distinguishing the diseased canopies from the healthy canopies. C-statistic values indicate the capability of spectral indices in the early (weeks 1–3, heading to fruiting stages) and late (weeks 4–6, ripening) periods (Table 4). Water band indices NDWI and WI yielded the best performance for discrimination for the early period, with c-statistics of 0.92 and 0.84, respectively, followed by other indices developed for chlorophyll estimation. Compared to healthy canopies, diseased canopies yielded high NDWI and low WI values at each time point throughout the entire observation period (Figures 6A,B). We also found that spectral indices performed differently in early (weeks 1–3, heading to fruiting stages) and late periods (weeks 4–6, ripening) of the STB epidemic (Table 4 and Figures 6C,D). Most of the spectral indices performed better in the early periods than in the late periods, except for the NPQI and DSI that showed better performances in the late periods. The presence of the spikes might contribute to the different performances of spectral indices at different sampling periods, since different indices and use of different spectral bands can have varied sensitivity to canopy structure changes due to spike presence.

TABLE 4
www.frontiersin.org

Table 4. The c-statistic of logistic regression model showing the performance of spectral indices in discriminating between the healthy (T0) and disease (T1) treatments of the four reference varieties [bold font highlights c-values ≥ 0.8 for the early (Weeks 1–3, heading to fruiting stages) and later stages (Weeks 4–6, ripening stage)].

FIGURE 6
www.frontiersin.org

Figure 6. Boxplots show the capability of (A) NDWI and (B) WI in discriminating between the Zymoseptoria tritici un-treated (T0) and treated canopies (T1) of the four reference varieties (CH CLARO, DRIFTER, RUNAL, and TITLIS) in different sampling periods, and ROC plots show the comparison of NDWI and WI at the early (C, Weeks 1–3) and late (D, Weeks 4–6) periods.

NDWI was the best index for detecting STB disease (Table 4 and Figures 6C,D), suggesting new insights for developing new crop sensors for early detection of diseases at the field or regional scales. Compared to the VIP scores of NDVI bands, both spectral bands of NDWI yielded significant contributions (VIP >1) in characterizing the differences between healthy and STB diseased canopies (Figure 2B). Until now, most of the hyperspectral cameras carried on near ground and low-altitude platforms (e.g., UAVs) for agricultural applications use only visible and NIR bands (e.g., 400–900 nm) (Garcia-Ruiz et al., 2013; Bareth et al., 2015), which constrains the potential of hyperspectral remote sensing for disease detection and disease mapping due to the absence of water absorption bands.

BLR models based on the selected spectral indices were applied to the independent validation dataset containing 331 varieties in week 3 (fruiting to early-ripening stages) as well as the same varieties in week 2 (early-fruiting stage). Our results indicate that NDWI enabled the best classification when considering sensitivity and specificity together, with an overall accuracy of 95% (kappa = 0.56) (Table 5). Model test sensitivity (i.e., TPR) ranged from 16.7 to 91.7% and test specificity (i.e., TNR) ranged from 37.2 to 97.2% (Table 5). NPQI and HI showed the best ability to correctly identify the diseased canopies, with sensitivities of 91.7 and 70.8%, respectively, whereas NDWI and CAI showed the best ability to correctly identify the healthy canopies, yielding specificities of 97.2 and 95.5%, respectively.

TABLE 5
www.frontiersin.org

Table 5. Confusion matrix of the disease classification using BLR models based on spectral indices and the PLSDA model when applied to the validation dataset (n = 758).

Full-Spectrum Analysis for STB Detection

Estimation of STB Metrics Using PLSR Models

To further investigate to what extent the full spectrum is able to account for the variations in STB metrics, PLSR models were constructed for the prediction of the STB scores and AUDPC, respectively. The optimal number of components used in the PLSR models for the STB scores and AUDPC was 4 and 3, respectively (Figure S3). The most influential bands (VIP >2) for STB scores and AUDPC were at the red-edge to NIR (700–800 and 1,100–1,150 nm) (Figures 7A,B). Red edge is associated with chlorophyll content and physiological changes due to stresses (Kumar et al., 2001; Sims and Gamon, 2002; Gitelson et al., 2003), which explains the high VIP scores. The VIP method allowed for selection of the spectral bands accounting for the variations in leaf pigment, water status, and internal structural changes. For the sake of comparison with the selected spectral indices, Figure 7 shows the relationships of the predicted STB scores and AUDPC against their measured values, as well as their Spearman's correlations. The PLSR model improved the prediction for STB scores, yielding a correlation coefficient of 0.41 (Figure 7C). Similarly, the PLSR model improved the prediction for AUDPC compared to the best index (SIPI, r = 0.42) and yielded an increased correlation coefficient (r = 0.50, Figure 7D).

FIGURE 7
www.frontiersin.org

Figure 7. VIP scores of the PLSR models for (A) the STB scores and (B) AUDPC. The three value breaks show different ranges of VIP values, which indicates the most influential wavelengths. PLSR model predicted values of (C) STB scores and (D) AUDPC plotted against their measured values with linear fit. Pearson and Spearman correlations produced the same correlation coefficients (r).

For the three leaf-based STB metrics, we constructed a multivariate PLSR model, i.e., with multiple response variables (PLACL, ρ-leaf, and ρ-lesion). To determine the suitable number of components, the PLSR model was first trained with 10 components, and the optimal number of components was determined based on the RMSEP and PRESS statistics (van der Voet, 1994). Results of cross validation showed that the use of four–six components allowed for low error rates and a small number of components (Figure S4). Therefore, the PLSR model for the leaf STB metrics was calibrated by using six components. The PLSR model improved the predictions for leaf STB metrics compared to spectral indices, yielding correlation coefficients of 0.60, 0.58, and 0.38 for PLACL, ρ-leaf, and ρ-lesion respectively (Figure 8).

FIGURE 8
www.frontiersin.org

Figure 8. The validation for the correlations between the observed and predicted values of PLACL, ρ-leaf and ρ-lesion using the PLSR model employing the number of components (comps) of 4, 5, and 6. Correlation coefficients (r) are for the Spearman's correlations.

In contrast to the selected spectral indices that use 2–3 bands, PLSR models make use of the full spectrum and select a subset of important bands that account for the more complete variations in infected canopies. Spectral bands yielding the highest VIP scores (≥ 4) were located in the red edge (730 nm) and NIR (e.g., peaks at 930, 1150 nm) regions (Figure S5), which again suggests the strong association of these spectral bands with STB infections and subsequent physiological changes of plants. Although the PLSR model accounted for smaller variations in the ρ-lesion than in the PLACL and ρ-leaf, the extent of improvement relative to spectral indices was largest for the ρ-lesion (PLSR vs. NDWI).

VI-based modeling would need different indices for individual canopy- or leaf-level STB metrics, suggesting that several individual models may be needed to achieve a certain accuracy (Mahlein et al., 2013; Jansen et al., 2014). In contrast, as observed here, the PLSR model allowed for improvement in linking canopy reflectance signatures to the three leaf STB metrics, PLACL, ρ-leaf, and ρ-lesion, by calibrating one model for multiple response variables (Figure 8). Further investigations of its potential for quantifying multi-infections and diseases are required, especially for the purpose of multiple diseases (Mahlein et al., 2013). Nevertheless, full-spectrum PLSR modeling shows promise for quantifying many traits associated with plant disease infections, e.g., damage to host plants and pathogen reproduction.

Classification of Healthy and Diseased Canopies Using PLSDA Model

Similar to the calibration of BLR models, canopy spectral data of the four reference varieties were used to calibrate the PLSDA model. The PLSDA model was first trained with ten components to determine the optimal number of components for the model, by employing a 5-fold cross-validation to compute the estimation errors. Results showed that the use of four components allowed for low prediction errors whilst maintaining high efficiency of the model (Figure S6), and thus it was decided to use four components for calibrating the model.

Selection of highly informative bands could improve the efficiency and accuracy of PLSDA models for estimation, as some bands have low signal-to-noise ratio and some bands might be redundant in term of similar responses to certain characteristics of plants and canopy (Thenkabail et al., 2011). In this study, band selection for the PLSDA model was performed based on the VIP statistics (Wold et al., 2001), and a total of 1,086 spectral bands that yielded VIP scores higher than 0.8 were selected for the PLSDA model. Results showed that the PLSDA model allowed for complete separation between the diseased and healthy canopies (Figure 9A), and spectral bands at wavelength 520–580, 700–1,300, 1,500–1,580, 1,760–1,800 nm and 2100–2400 were effective for STB disease classification (Figure 9B, Table S1). Other studies have also shown that NIR and SWIR spectral bands are very efficient for the detection of plant diseases (Delalieux et al., 2009; Bauriegel and Herppich, 2014).

FIGURE 9
www.frontiersin.org

Figure 9. (A) Scatterplot shows the observations that are projected onto the first two components of the PLSDA model, suggesting the complete separation of the two groups (0 = healthy, 1 = diseased). (B) Variable importance in projections (VIP) scores at different wavelengths projected in the first three components (comp 1, 2, and 3). (C) Scatter plot shows the relationship between NDWI and AUDPC, and red circles highlight the observations of Zymoseptoria tritici un-treated plots (T0) that were classified as diseased canopies by the PLSDA model. (D) Boxplot shows the difference in NDWI between the correctly classified and misclassified T0 plots. (E) Boxplot shows the difference in AUDPC between the correctly classified and misclassified T0 plots.

The trained PLSDA model was applied to the independent dataset [710 healthy (T0) and 48 diseased (T1) plots/canopies] to test the capability of classification. Results showed that the overall classification accuracy was 93% (Table 5). Model test sensitivity (i.e., TPR) was 95.8%, and test specificity (i.e., TNR) was 92.7%, i.e., only two diseased (T1) plots were misclassified as healthy plots (Table 5). 52 plots were classified as belonging to the diseased canopies, accounting for only 7% of the 710 healthy plots. This finding led us to consider whether the misclassified plots were planted to cultivars that were highly susceptible to STB. PLACL of individual leaves in the 52 plots ranged from 21 to 82 with a mean value of 50 (data not shown), placing these plots mainly in the upper range of the overall STB rating score (≥ 3, (Figure S7). In 22 of these plots, PLACL was higher than 50, indicating that these plots had significant STB infections despite the intensive use of fungicides. Hence, the classification of these plots as diseased by the model was reasonable (an expected misclassification), providing additional evidence for the power of this method. Moreover, the NDVI and NDWI values for these plots were relatively low compared to the correctly classified healthy plots (Figure 9C and Figure S7), in agreement with the negative correlations between NDVI or NDWI, on the one hand, and STB scores and AUDPC, on the other hand (Figure S1 and Figure 3). Incorrect classification of the 30 plots with PLACL values lower than 50 might be due to special phenotypic characteristics—for instance, a group of wheat genotypes having lighter green leaves. Further investigations will be needed to determine the cause of these misclassifications.

Calculated based on three independent measures of disease over time, AUDPC is often used in resistance breeding because it reduces the rating bias inherent in individual assessments (Li et al., 2006; Karisto et al., 2018). Here, the AUDPC values for the 52 misclassified plots spanned a very wide range compared to NDWI, which is in line with the capabilities of separation between diseased and healthy canopies (Figures 9D,E). NDVI indicates overall greenness of the canopy, and dark green canopies (high NDVI) are more likely to be visually rated at low STB scores. Results showed that there was a certain degree of consistency between the visual assessment and spectral indices using the visible region (e.g., NDVI, Figure S7), which also suggests that visual assessment might underestimate the infection when mild infection (e.g., scores ranged from 0 to 5 in this study) occurs in canopies of high NDVI. In contrast, NDWI improved separation as compared to NDVI and AUDPC, implying the possibility of hyperspectral canopy sensing of mild infections or early occurrence of STB infections using NIR to SWIR bands (Figure 9C). However, NDVI is also indicative of disease severity when the infection is severe (Jansen et al., 2014). Compared to visual assessments that depend largely on individual skill and training, canopy spectral metrics tend to be more reproducible and less prone to human bias (Steddom et al., 2005).

The PLSDA training dataset contains only four varieties, which might have constrained the power of the PLSDA model. However, the validation on an independent dataset of 331 varieties shows the promise of combining canopy hyperspectral measurements and PLSDA modeling for crop disease detection. Compared to spectral indices, the PLSDA model allowed for the best-balanced ability, yielding model test sensitivity and specificity of 95.8 and 92.7%, respectively and the overall accuracy of 93% (kappa = 0.60) (Table 5). Having the highest sensitivity, the PLSDA model suggests the great potential of sensing early infections before significant visible symptoms develop in the canopy.

Despite the promise of hyperspectral remote sensing for plant disease phenotyping, further investigations including a wider spectrum of plant diseases, an even larger number of crop varieties, and a broader range of environments will be required to improve crop disease detection and differentiation, and determine whether the desired precision, reproducibility and throughput needed for resistance breeding can be achieved (Furbank and Tester, 2011; O'Driscoll et al., 2014; Mahlein, 2016). Moreover, further investigations will be needed to determine whether canopy reflectance and discriminant analysis methods can differentiate STB from other symptoms such as nutrient or water stresses (Bürling et al., 2011). Future work should also investigate the possibility that fungicide applications might affect the comparison of canopy hyperspectral responses. Progress in addressing these questions will advance plant disease phenotyping and enable improved mapping for precise applications of fungicides in the field, especially from low-altitude sensor carrying platforms such as UAVs (Garcia-Ruiz et al., 2013; Yuan et al., 2016; Yang et al., 2017).

Conclusions

In this study we investigated the feasibility of using canopy hyperspectral data to detect wheat STB infection in a population of 335 elite European wheat varieties. Our results demonstrated that hyperspectral indices used for the estimation of leaf chlorophyll and water status were negatively related to the canopy and leaf STB metrics, whereas the senescence related spectral indices were positively related. NDWI is more indicative of early infections compared to NDVI. These findings lead us to hypothesize that STB infections induced leaf biochemical changes and accelerated plant senescence, resulting in early changes in the spectral signature associated with water absorption. Full-spectrum analysis through the partial least squares (PLS) method allowed for the prediction of canopy- and leaf-level STB severity metrics, suggesting the promise of using canopy spectra as new tools to quantify STB damage across a wide array of wheat genotypes. PLS discriminant analysis empowered the canopy hyperspectral measurements to discriminate STB diseased canopies from healthy canopies with high accuracy, and showed the potential of using the NIR to SWIR region to discriminate diseased canopies in this study. Additionally, hyperspectral narrow bands in the red edge, NIR to SWIR regions were identified as highly efficient bands for disease discrimination. Our findings suggest that canopy hyperspectral measurements have great potential to facilitate high throughput plant phenotyping for disease resistance breeding. Further investigations that lead to a better understanding of the interaction between disease progress, canopy structure, and physiological senescence will further improve our ability to identify and measure diseases using remote sensing technologies.

Author Contributions

KY, AM, and AH designed the experiment. KY and JA performed the canopy spectral measurement. AM and PK performed the leaf sampling and leaf STB imagery assessment. FM performed the visual assessment of STB infection. KY analyzed the data and wrote the draft. BM, AW, and AH supervised the project and designed the research. All authors read, revised, and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the Swiss Federal Office for Agriculture FOAG, Switzerland. AM and PK were supported by the SNSF Ambizione (Grant No. PZ00P3_161453). We thank Mr. Hansueli Zellweger, Ms. Maria Ines Falconi, and Mr. Williams Okrah for their assistances in performing spectral measurements. We would also like to thank IPK Gatersleben and Delley Seeds and Plants Ltd for the supply of wheat varieties.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01195/full#supplementary-material

References

Ashourloo, D., Mobasheri, M. R., and Huete, A. (2014a). Developing two spectral disease indices for detection of wheat leaf rust (Pucciniatriticina). Remote Sens. 6, 4723–4740. doi: 10.3390/rs6064723

CrossRef Full Text | Google Scholar

Ashourloo, D., Mobasheri, M. R., and Huete, A. (2014b). Evaluating the effect of different wheat rust disease symptoms on vegetation indices using hyperspectral measurements. Remote Sens. 6, 5107–5123. doi: 10.3390/rs6065107

CrossRef Full Text | Google Scholar

Bajwa, S. G., Rupe, J. C., and Mason, J. (2017). Soybean disease monitoring with leaf reflectance. Remote Sens. 9:127. doi: 10.3390/rs9020127

CrossRef Full Text | Google Scholar

Bareth, G., Aasen, H., Bendig, J., Gnyp, M. L., Bolten, A., Jung, A., et al. (2015). Low-weight and UAV-based hyperspectral full-frame cameras for monitoring crops: spectral comparison with portable spectroradiometer measurements. Photogramm. Fernerkund. Geoinf. 2015, 69–79. doi: 10.1127/pfg/2015/0256

CrossRef Full Text | Google Scholar

Bauriegel, E., and Herppich, W. B. (2014). Hyperspectral and chlorophyll fluorescence imaging for early detection of plant diseases, with special reference to Fusarium spec. Infect. Wheat Agric. 4, 32–57. doi: 10.3390/agriculture4010032

CrossRef Full Text | Google Scholar

Bürling, K., Hunsche, M., and Noga, G. (2011). Use of blue–green and chlorophyll fluorescence measurements for differentiation between nitrogen deficiency and pathogen infection in winter wheat. J. Plant Physiol. 168, 1641–1648. doi: 10.1016/j.jplph.2011.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Deery, D., Jimenez-Berni, J., Jones, H., Sirault, X., and Furbank, R. (2014). Proximal remote sensing buggies and potential applications for field-based phenotyping. Agronomy 4, 349–379. doi: 10.3390/agronomy4030349

CrossRef Full Text | Google Scholar

Delalieux, S., Auwerkerken, A., Verstraeten, W. W., Somers, B., Valcke, R., Lhermitte, S., et al. (2009). Hyperspectral reflectance and fluorescence imaging to detect scab induced stress in apple leaves. Remote Sens. 1, 858–874. doi: 10.3390/rs1040858

CrossRef Full Text | Google Scholar

Eriksson, L., Byrne, T., Johansson, E., Trygg, J., and Vikström, C. (2013). Multi- and Megavariate Data Analysis Basic Principles and Applications. Malmö: Umetrics Academy.

Fones, H., and Gurr, S. (2015). The impact of Septoria tritici blotch disease on wheat: an EU perspective. Fungal Genet. Biol. 79, 3–7. doi: 10.1016/j.fgb.2015.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Furbank, R. T., and Tester, M. (2011). Phenomics – technologies to relieve the phenotyping bottleneck. Trends Plant Sci. 16, 635–644. doi: 10.1016/j.tplants.2011.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamon, J. A., Serrano, L., and Surfus, J. S. (1997). The photochemical reflectance index: an optical indicator of photosynthetic radiation use efficiency across species, functional types, and nutrient levels. Oecologia 112, 492–501. doi: 10.1007/s004420050337

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, B. (1996). NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 58, 257–266. doi: 10.1016/S0034-4257(96)00067-3

CrossRef Full Text | Google Scholar

Garcia-Ruiz, F., Sankaran, S., Maja, J. M., Lee, W. S., Rasmussen, J., and Ehsani, R. (2013). Comparison of two aerial imaging platforms for identification of Huanglongbing-infected citrus trees. Comput. Electron. Agric. 91, 106–115. doi: 10.1016/j.compag.2012.12.002

CrossRef Full Text | Google Scholar

Gitelson, A. A., Gritz, Y., and Merzlyak, M. N. (2003). Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 160, 271–282. doi: 10.1078/0176-1617-00887

PubMed Abstract | CrossRef Full Text | Google Scholar

Hosmer, D. W., and Lemeshow, S. (2005). Applied Logistic Regression. 2nd Edn. New York, NY: John Wiley & Sons.

Jansen, M., Bergsträsser, S., Schmittgen, S., Müller-Linow, M., and Rascher, U. (2014). Non-invasive spectral phenotyping methods can improve and accelerate cercospora disease scoring in sugar beet breeding. Agriculture 4, 147–158. doi: 10.3390/agriculture4020147

CrossRef Full Text | Google Scholar

Karisto, P., Hund, A., Yu, K., Anderegg, J., Walter, A., Mascher, F., et al. (2018). Ranking quantitative resistance to Septoria tritici blotch in elite wheat cultivars using automated image analysis. Phytopathology 108, 568–581. doi: 10.1094/PHYTO-04-17-0163-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirchgessner, N., Liebisch, F., Yu, K., Pfeifer, J., Friedli, M., Hund, A., et al. (2017). The ETH field phenotyping platform FIP: a cable-suspended multi-sensor system. Funct. Plant Biol. 44, 154–168. doi: 10.1071/FP16165

CrossRef Full Text | Google Scholar

Kollers, S., Rodemann, B., Ling, J., Korzun, V., Ebmeyer, E., Argillier, O., et al. (2013). Whole genome association mapping of Fusarium head blight resistance in european winter wheat (Triticum aestivum L.). PLOS ONE 8:e57500. doi: 10.1371/journal.pone.0057500

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, L., Schmidt, K., Dury, S., and Skidmore, A. (2001). “Imaging spectrometry and vegetation science,” in Imaging Spectrometry - Basic Principles and Prospective Applications Remote Sensing and Digital Image Processing, eds F. D. van der Meer and S. M. De Jong (Dordrecht: Springer), 111–155.

Google Scholar

Laudien, R., Bareth, G., and Doluschitz, R. (2003). “Analysis of hyperspectral field data for detection of sugar beet diseases,” in Proceedings of the EFITA Conference (Debrecen), 375–381.

Google Scholar

Li, Z. F., Xia, X. C., Zhou, X. C., Niu, Y. C., He, Z. H., Zhang, Y., et al. (2006). Seedling and slow rusting resistance to stripe rust in chinese common wheats. Plant Dis. 90, 1302–1312. doi: 10.1094/PD-90-1302

CrossRef Full Text | Google Scholar

Mahlein, A.-K. (2016). Plant disease detection by imaging sensors – parallels and specific demands for precision agriculture and plant phenotyping. Plant Dis. 100, 241–251. doi: 10.1094/PDIS-03-15-0340-FE

CrossRef Full Text | Google Scholar

Mahlein, A.-K., Rumpf, T., Welke, P., Dehne, H.-W., Plümer, L., Steiner, U., et al. (2013). Development of spectral indices for detecting and identifying plant diseases. Remote Sens. Environ. 128, 21–30. doi: 10.1016/j.rse.2012.09.019

CrossRef Full Text | Google Scholar

Mattupalli, C., Komp, M. R., and Young, C. A. (2015). Integrating geospatial technologies and unmanned aircraft systems into the grower's disease management toolbox. APS Features. Available online at: http://www.apsnet.org/publications/apsnetfeatures/Pages/Toolbox.aspx

McHugh, M. L. (2012). Interrater reliability: the kappa statistic. Biochem. Med. 22, 276–282.

PubMed Abstract | Google Scholar

Merzlyak, M. N., Gitelson, A. A., Chivkunova, O. B., and Rakitin, V. Y. (1999). Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening. Physiol. Plant. 106, 135–141. doi: 10.1034/j.1399-3054.1999.106119.x

CrossRef Full Text | Google Scholar

MeteoSchweiz (2016). Climate Normals Zürich / Affoltern. Reference period 1981–2010. Federal Office of Meteorology and Climatology MeteoSwiss. Available online at: http://www.meteoswiss.admin.ch

Mevik, B.-H., and Wehrens, R. (2007). The pls package: principal component and partial least squares regression in R. J. Stat. Softw. 18, 1–23. doi: 10.18637/jss.v018.i02

CrossRef Full Text | Google Scholar

O'Driscoll, A., Kildea, S., Doohan, F., Spink, J., and Mullins, E. (2014). The wheat–Septoria conflict: a new front opening up? Trends Plant Sci. 19, 602–610. doi: 10.1016/j.tplants.2014.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Oerke, E.-C., Mahlein, A.-K., and Steiner, U. (2014). “Proximal sensing of plant diseases,” in Detection and Diagnostics of Plant Pathogens Plant Pathology in the 21st Century, eds M. L. Gullino and P. J. M. Bonants (Dordrecht: Springer), 55–68.

Google Scholar

Oerke, E. C., Herzog, K., and Toepfer, R. (2016). Hyperspectral phenotyping of the reaction of grapevine genotypes to Plasmopara viticola. J. Exp. Bot. 67, 5529–5543. doi: 10.1093/jxb/erw318

PubMed Abstract | CrossRef Full Text | Google Scholar

Oppelt, N. (2002). Monitoring of Plant Chlorophyll and Nitrogen Status Using the Airborne Imaging Spectrometer AVIS. Available online at: https://edoc.ub.uni-muenchen.de/354/ (Accessed July 12, 2017).

Peñuelas, J., Filella, I., Lloret, P., Muñoz, F., and Vilajeliu, M. (1995a). Reflectance assessment of mite effects on apple trees. Int. J. Remote Sens. 16, 2727–2733. doi: 10.1080/01431169508954588

CrossRef Full Text | Google Scholar

Peñuelas, J., Frederic, B., and Filella, I. (1995b). Semi-empirical indices to assess carotenoids/chlorophyll a ratio from leaf spectral reflectance. Photosynthetica 31, 221–230.

Google Scholar

Peñuelas, J., and Inoue, Y. (1999). Reflectance indices indicative of changes in water and pigment contents of peanut and wheat leaves. Photosynthetica 36, 355–360. doi: 10.1023/A:1007033503276

CrossRef Full Text | Google Scholar

Pinter, J., Paul, J., Hatfield, J. L., Schepers, J. S., Barnes, E. M., Moran, M. S., et al. (2003). Remote sensing for crop management. Photogramm. Eng. Remote Sens. 69, 647–664. doi: 10.14358/PERS.69.6.647

CrossRef Full Text | Google Scholar

R Core Team (2016). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing Available online at: https://www.R-project.org/

Rivero, R. M., Kojima, M., Gepstein, A., Sakakibara, H., Mittler, R., Gepstein, S., et al. (2007). Delayed leaf senescence induces extreme drought tolerance in a flowering plant. Proc. Natl. Acad. Sci. U.S.A. 104, 19631–19636. doi: 10.1073/pnas.0709453104

PubMed Abstract | CrossRef Full Text | Google Scholar

Rohart, F., Gautier, B., Singh, A., and Lê Cao, K. A. (2017). mixOmics: an R package for'omics feature selection and multiple data integration. PLoS Comput. Biol. 13:e1005752. doi: 10.1371/journal.pcbi.1005752

PubMed Abstract | CrossRef Full Text | Google Scholar

Sims, D. A., and Gamon, J. A. (2002). Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 81, 337–354. doi: 10.1016/S0034-4257(02)00010-X

CrossRef Full Text | Google Scholar

Steddom, K., Bredehoeft, M. W., Khan, M., and Rush, C. M. (2005). Comparison of visual and multispectral radiometric disease evaluations of cercospora leaf spot of sugar beet. Plant Dis. 89, 153–158. doi: 10.1094/PD-89-0153

CrossRef Full Text | Google Scholar

Stewart, E. L., Hagerty, C. H., Mikaberidze, A., Mundt, C. C., Zhong, Z., and McDonald, B. A. (2016). An improved method for measuring quantitative resistance to the wheat pathogen Zymoseptoria tritici using high-throughput automated image analysis. Phytopathology 106, 782–788. doi: 10.1094/PHYTO-01-16-0018-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Thenkabail, P. S., Lyon, J. G., and Huete, A. (2011). Hyperspectral Remote Sensing of Vegetation. Boca Raton, FL: CRC Press.

Torriani, S. F., Melichar, J. P., Mills, C., Pain, N., Sierotzki, H., and Courbot, M. (2015). Zymoseptoria tritici: a major threat to wheat production, integrated approaches to control. Fungal Genet. Biol. 79, 8–12. doi: 10.1016/j.fgb.2015.04.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Ustin, S. L., Gitelson, A. A., Jacquemoud, S., Schaepman, M., Asner, G. P., Gamon, J. A., et al. (2009). Retrieval of foliar information about plant pigment systems from high resolution spectroscopy. Remote Sens. Environ. 113 (Suppl. 1), S67–S77. doi: 10.1016/j.rse.2008.10.019

CrossRef Full Text | Google Scholar

van der Voet, H. (1994). Comparing the predictive accuracy of models using a simple randomization test. Chemom. Intell. Lab. Syst. 25, 313–323. doi: 10.1016/0169-7439(94)85050-X

CrossRef Full Text | Google Scholar

Wakie, T. T., Kumar, S., Senay, G. B., Takele, A., and Lencho, A. (2016). Spatial prediction of wheat septoria leaf blotch (Septoria tritici) disease severity in Central Ethiopia. Ecol. Inform. 36, 15–30. doi: 10.1016/j.ecoinf.2016.09.003

CrossRef Full Text | Google Scholar

Walter, A., Liebisch, F., and Hund, A. (2015). Plant phenotyping: from bean weighing to image analysis. Plant Methods 11:14. doi: 10.1186/s13007-015-0056-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wold, S., Sjöström, M., and Eriksson, L. (2001). PLS-regression: a basic tool of chemometrics. Chemom. Intell. Lab. Syst. 58, 109–130. doi: 10.1016/S0169-7439(01)00155-1

CrossRef Full Text | Google Scholar

Yang, C.-M., Cheng, C.-H., and Chen, R.-K. (2007). Changes in spectral characteristics of rice canopy infested with brown planthopper and leaffolder. Crop Sci. 47, 329–335. doi: 10.2135/cropsci2006.05.0335

CrossRef Full Text | Google Scholar

Yang, G., Liu, J., Zhao, C., Li, Z., Huang, Y., Yu, H., et al. (2017). Unmanned aerial vehicle remote sensing for field-based crop phenotyping: current status and perspectives. Front. Plant Sci. 8:1111. doi: 10.3389/fpls.2017.01111

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, K., Gnyp, M. L., Gao, L., Miao, Y., Chen, X., and Bareth, G. (2015). Estimate leaf chlorophyll of rice using reflectance indices and partial least squares. Photogramm. Fernerkund. Geoinf. 2015, 45–54. doi: 10.1127/pfg/2015/0253

CrossRef Full Text | Google Scholar

Yu, K., Leufen, G., Hunsche, M., Noga, G., Chen, X., and Bareth, G. (2014). Investigation of leaf diseases and estimation of chlorophyll concentration in seven barley varieties using fluorescence and hyperspectral indices. Remote Sens. 6, 64–86. doi: 10.3390/rs6010064

CrossRef Full Text | Google Scholar

Yuan, L., Pu, R., Zhang, J., Wang, J., and Yang, H. (2016). Using high spatial resolution satellite imagery for mapping powdery mildew at a regional scale. Precis. Agric. 17, 332–348. doi: 10.1007/s11119-015-9421-x

CrossRef Full Text | Google Scholar

Zarco-Tejada, P. J., Miller, J. R., Noland, T. L., Mohammed, G. H., and Sampson, P. H. (2001). Scaling-up and model inversion methods with narrowband optical indices for chlorophyll content estimation in closed forest canopies with hyperspectral data. IEEE Trans. Geosci. Remote Sens. 39, 1491–1507. doi: 10.1109/36.934080

CrossRef Full Text | Google Scholar

Keywords: field-based phenotyping, hyperspectral remote sensing, partial least squares discriminant analysis, plant phenotyping, plant spectroscopy, resistance breeding, Septoria tritici blotch, wheat disease

Citation: Yu K, Anderegg J, Mikaberidze A, Karisto P, Mascher F, McDonald BA, Walter A and Hund A (2018) Hyperspectral Canopy Sensing of Wheat Septoria Tritici Blotch Disease. Front. Plant Sci. 9:1195. doi: 10.3389/fpls.2018.01195

Received: 31 October 2017; Accepted: 26 July 2018;
Published: 17 August 2018.

Edited by:

Chunjiang Zhao, Beijing Academy of Agricultural and Forestry Sciences, China

Reviewed by:

Sébastien Saint-Jean, AgroParisTech Institut des Sciences et Industries du Vivant et de L'environnement, France
Conxita Royo, Institut de Recerca i Tecnologia Agroalimentàries (IRTA), Spain

Copyright © 2018 Yu, Anderegg, Mikaberidze, Karisto, Mascher, McDonald, Walter and Hund. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Kang Yu, kang.yu@alumni.ethz.ch

Present Address: Kang Yu, KU Leuven, Remote Sensing & Terrestrial Ecology, Department of Earth and Environmental Sciences, Leuven, Belgium