Abstract

Stripe rust caused by Puccinia striiformis f. sp. tritici (Pst) is a devastating wheat disease worldwide. Potential application of near-infrared spectroscopy (NIRS) in detection of pathogen amounts in latently Pst-infected wheat leaves was investigated for disease prediction and control. A total of 300 near-infrared spectra were acquired from the Pst-infected leaf samples in an incubation period, and relative contents of Pst DNA in the samples were obtained using duplex TaqMan real-time PCR arrays. Determination models of the relative contents of Pst DNA in the samples were built using quantitative partial least squares (QPLS), support vector regression (SVR), and a method integrated with QPLS and SVR. The results showed that the kQPLS-SVR model built with a ratio of training set to testing set equal to 3 : 1 based on the original spectra, when the number of the randomly selected wavelength points was 700, the number of principal components was 8, and the number of the built QPLS models was 5, was the best. The results indicated that quantitative detection of Pst DNA in leaves in the incubation period could be implemented using NIRS. A novel method for determination of latent infection levels of Pst and early detection of stripe rust was provided.

1. Introduction

Wheat stripe rust, caused by the biotrophic pathogen Puccinia striiformis f. sp. tritici (Pst), is an important wheat disease worldwide [15]. This disease can cause severe yield losses. Once an epidemic of this disease occurs, at least 10%–30% of wheat yield can be reduced [1]. In China, wheat stripe rust occurs in almost all wheat-growing regions, and it is the most devastating wheat disease [1, 4, 5]. This disease is always a great potential threat to the national wheat production in China. Since 1950, many severe epidemics of wheat stripe rust have occurred in China [1, 5], among which the four most destructive epidemics occurred in 1950, 1964, 1990, and 2002, resulting in yield losses of 6.0, 3.2, 1.8, and 1.3 million tons, respectively [4, 6].

Generally, an infection process of wheat stripe rust in the field can be divided into four stages, that is, the contact period, the penetration period, the incubation period, and the diseased period [1]. Under favorable environmental conditions, Pst urediospores landing on the surface of wheat leaves germinate, produce germ tubes and appressoria, and then penetrate into wheat leaf tissues [1, 5]. During the incubation period, a large quantity of Pst hyphae accumulate gradually in the infected leaf tissues [1, 5], but it is very difficult to observe disease symptoms on the surface of the infected wheat leaves with naked eyes, causing lots of difficulties for early monitoring and prediction of this disease. Once uredinia appear on the surface of wheat leaves, a large number of urediospores may be liberated after rupture of the uredinia and then be spread by air, which may result in secondary infections. So it is of great significance to conduct the early detection of Pst infections and the quantitative detection of pathogen amounts in the infected plants for early precise site-specific control of the disease, disease prediction, and control strategy making. In particular, it is very important to realize the early detection of pathogen amounts in the infected plants in the overwintering regions and the oversummering regions of Pst for the macro control of wheat stripe rust. Traditionally, the monitoring of wheat stripe rust is carried out via field investigation. This method is time-consuming and laborious. Moreover, using this conventional method, only the diseased wheat fields with disease symptoms can be surveyed, and the infected wheat leaves without symptom appearance in the incubation period cannot be identified accurately and rapidly. At present, some techniques and methods, including molecular biology techniques [710], hyperspectral remote sensing technology [1114], thermal infrared imaging technology [15], and near-infrared spectroscopy (NIRS) [16, 17], have been applied to the monitoring and early detection of wheat stripe rust.

There are many reports on early qualitative detection of wheat stripe rust. The Pst infection of a wheat leaf in the incubation period of stripe rust can be detected accurately and qualitatively by detecting the presence of wheat stripe rust pathogen using molecular biology techniques such as polymerase chain reaction (PCR) assay and loop-mediated isothermal amplification (LAMP) assay [79]. Hyperspectral remote sensing technology has been used to identify healthy wheat plants and Pst-infected wheat plants without symptoms [18]. It was reported that Pst latent infections in wheat leaves could be detected using thermal infrared imaging technology [15] and NIRS technology [16, 17] before symptoms appear. A study on the detection of Pst latent infections in wheat leaves incubated in an artificial climate chamber using NIRS technology was conducted by Li et al. [17], and the results obtained by analyzing the spectral data demonstrated that the latently infected leaves could be distinguished from healthy wheat leaves as early as one day after artificial inoculation with the suspension of Pst urediospores. However, there are relatively few studies on the early quantitative detection of Pst in the incubation period of wheat stripe rust. The quantitative detection of Pst in the latently infected wheat leaves was implemented mainly relying on a real-time PCR method [10, 19].

The detection process using a molecular biological method is complex, and the method has high requirements on technologies and instruments [20, 21]. With the development of science and technology, portable fluorescence quantitative PCR instruments have been applied gradually in recent years and the instruments for quantitative analysis in the field have been available [20]. However, it is still necessary for users to master the relevant molecular biology skills, and a large number of samples to be tested are required to be obtained. Only using molecular biology technology, it is very difficult to quickly get enough data in a large scale to provide information for disease forecasting and disease warning. The instruments used in hyperspectral remote sensing are very expensive. Thermal infrared imaging technology has high requirements on imaging resolution and temperature resolution of the used instrument. At the present time, detection methods of latent infection of pathogens based on molecular biology technology, hyperspectral remote sensing technology, and thermal infrared imaging technology are very difficult to be popularized in practice. Therefore, it is critical to explore a simple, convenient, accurate, and rapid method for early detection of wheat stripe rust.

As a kind of nondestructive analytical technology, NIRS can be used to carry out both qualitative analysis and quantitative analysis of samples [22, 23]. The analysis process of NIRS is simple and rapid, and it is a fast detection technology with low cost and is suitable for online analysis [22, 23]. It has been widely used in many fields such as agriculture, food industry, chemical industry, pharmaceutical industry, and petroleum industry [2127]. NIRS has been applied to implement early qualitative identification of wheat stripe rust and wheat leaf rust (caused by Puccinia triticina) before disease symptoms appeared [16] and assess disease severity of wheat stripe rust [28]. It has also been used to qualitatively identify urediospores of Pst and P. triticina and quantitatively determine the content of each pathogen in a mixture of two kinds of pathogens including Pst and P. triticina [29]. Moreover, using NIRS, wheat leaves infected with Pst could be identified before symptoms appear [17]. However, to the best of our knowledge, there are no reports on quantitative detection of Pst in wheat leaves in the incubation period of wheat stripe rust using NIRS.

In this study, a method based on NIRS was explored for quantitative determination of the relative content of Pst DNA in wheat leaves in the incubation period. Based on the acquired near-infrared spectral data of the latently Pst-infected wheat leaf samples and the data on the relative content of Pst DNA in the corresponding samples obtained using the duplex TaqMan real-time PCR arrays, dynamic changes of the relative contents of Pst DNA in the latently Pst-infected wheat leaves during the incubation period were investigated and quantitative determination models were built to realize the quantitative and rapid detection of the pathogen amounts in the latently Pst-infected wheat leaves. The aim of this study was to provide a method for rapid, nondestructive, and quantitative determination of latent infection levels of Pst in wheat leaves and early detection of wheat stripe rust. Furthermore, support information can be provided for prediction and control of this disease. Some methodological references can also be provided for the early quantitative and nondestructive detection of other diseases.

2. Materials and Methods

2.1. Materials

CYR 33, a dominant physiological race of Pst in China, was used in this study. A highly susceptible wheat cultivar Mingxian 169 was used to multiply the pathogen in an artificial climate chamber in the Laboratory of Plant Disease Epidemiology, Department of Plant Pathology, China Agricultural University. Artificial inoculation by spraying urediospores of CYR 33 onto the surface of seedling leaves of Mingxian 169 was conducted to obtain the latently Pst-infected wheat leaves.

2.2. Multiplication of Wheat Stripe Rust Pathogen

Urediospores of the Pst physiological race CYR33 were multiplied using a similar method as described by Cheng et al. [30]. After seeds of Mingxian 169 were soaked for 24 h in sterile water, the plump and well-germinated seeds were selected and sown in pots (10 cm in diameter) with approximately 25 seeds per pot. Then the pots were incubated in an artificial climate chamber at 11–13°C with 60%–70% relative humidity (RH) and 14 h of light per day (10,000 lux) and were watered when necessary. When the first leaves of wheat seedlings fully expanded, urediospores of CYR 33 stored in a liquid nitrogen container were taken out, reactivated in warm water of 40–45°C for 5 min, and then hydrated at 4°C for 12 h. A spore suspension was made with 0.2% Tween 80 solution for artificial spray inoculation. Immediately after inoculation, each pot with the inoculated wheat seedlings was covered with a clear glass cylinder that was covered with two layers of sterile cotton gauze on the top, and all pots with the inoculated wheat seedlings were transferred into a moist chamber under dark conditions at 11–13°C for 24 h. Subsequently, the inoculated wheat seedlings were placed into the artificial climate chamber under the conditions described above and incubated until a large number of urediospores were produced on the surface of the seedling leaves. The fresh urediospores were harvested for the following experiments or stored in the liquid nitrogen container for later use.

2.3. Collection of Latently Pst-Infected Wheat Leaves

Using the artificial spray inoculation method described above, healthy wheat seedlings with the fully expanded first leaves were inoculated with a 0.15 mg/mL spore suspension prepared with 3 mg harvested fresh urediospores of CYR33 and 0.2% Tween 80 solution. A total of 50 pots of wheat seedlings were inoculated for further experiments. Sixty wheat leaves with uniform size and same growth vigor were collected every 24 h after inoculation until disease symptoms appeared and the incubation period of wheat stripe rust ended. The incubation period was defined as the number of days from inoculation to rupture of the first uredinium [31], and it was 10 days in this study. Two leaves were treated as a sample for acquisition of near-infrared spectral data, and 30 samples per day were collected. A total of 300 samples were collected in this study.

2.4. Acquisition of Near-Infrared Spectral Data

The near-infrared spectra of the latently Pst-infected wheat leaf samples were acquired by using FT-NIR MPA spectrometer (Bruker, Germany). Wheat leaves of each sample were cut into small square fragments and then were placed into a sample cup (20 mm in diameter) for spectral measurement. Using integrating sphere diffuse reflectance method, the spectra in a range of 4000–12,000 cm−1 were measured with the spectral resolution of 8 cm−1 and the number of scan processes of 32. Each acquired spectrum contained 2100 wavelength points. After spectral acquisition, each sample was put into a 2 mL grinding tube. Each tube was numbered and immediately kept at −80°C for subsequent DNA extraction. On each day during the incubation period, the near-infrared spectra of 30 collected samples were acquired. In this study, a total of 300 spectra were obtained as shown in Figure 1.

2.5. DNA Extraction and Quantitative Measurements of DNA in Wheat Leaf Samples

Healthy wheat leaves, 3 mg of urediospores of Pst, and the latently Pst-infected wheat leaves after acquisition of near-infrared spectral data were used as samples for DNA extraction. The extracted wheat DNA from the healthy leaves and the extracted Pst DNA from the Pst urediospores were used to generate standard curves. DNA was extracted using a modification of the procedure described by Justesen et al. [32]. Each sample in a 2 mL grinding tube was mixed with 0.4 g of quartz sands, one glass bead, and 600 μL of 2% cetyl trimethyl ammonium bromide (CTAB) buffer. The CTAB buffer was prepared with 2 g of CTAB and 1 g of polyvinyl pyrrolidone by adding 10 μL of 1 M Tris-HCl (pH = 8.0), 2 μL of 0.5 M ethylene diamine tetraacetic acid (pH = 8.0), and 28 mL of 5 M NaCl into a final volume of 100 mL adjusted with deionized water, then was sterilized under high temperature and high pressure, and finally was added with 100 μL β-mercapto ethanol after cooling. The tube was shaken in a FastPrep-24 instrument (MP Biomedicals, Santa Ana, CA, USA) at 6.0 m/s for two periods of 40 s, with 5 min cooling on ice between them. After incubation in water bath at 65°C for 1 h with gentle shakes every 15 min, the tube was added with 600 μL of chloroform/isoamyl alcohol (v : v = 24 : 1) with mixing and then was centrifuged at the speed of 12,000 rpm for 10 min. The supernatant was transferred to a new clean 1.5 mL centrifuge tube, and 0.6 volumes of cold isopropanol at −20°C was added. The mixture in the tube was shaken gently and then was incubated at −20°C for 1 h. The tube was centrifuged at the speed of 12,000 rpm for 10 min, and then the liquid supernatant was abandoned. The resultant was added with 500 μL of 70% ethanol for rinse, shaken gently, and then centrifuged at 12000 rpm for 10 min. After discarding the supernatant and drying the resultant in the air, the DNA was dissolved in 50 μL of sterile double distilled water and kept at −20°C for later use.

Primers and probes used for the duplex TaqMan real-time PCR assays in this study were listed in Table 1. When the quantity of wheat DNA was determined using duplex TaqMan real-time PCR, TAG2315F and TAG2473R reported by Sandberg et al. [33] that were designed based on the DNA sequence of a prolamin gene of wheat were used as the primers, and the probe TAG-Pr1 was designed based on this DNA sequence. When the quantity of Pst DNA was determined using duplex TaqMan real-time PCR, the primers (Pst-F and Pst-R) and the probe Pst-P reported by Li et al. [34] that were designed based on the internal transcribed spacer (ITS) region sequence of Pst were used.

Each real-time PCR assay was performed in a volume of 20 μL containing 3.20 μL of MgCl2 (25 μM), 2.00 μL of dNTP (2500 μM), 2.00 μL of Taq buffer (10x), 0.60 μL of Taq (5 U/μL), 0.30 μL of each primer (Pst-F, Pst-R, TAG2315F, and TAG2473R, 10 μM each), 0.25 μL of each probe (Pst-P and TAG-Pr1, 10 μM each), 2.00 μL of template DNA, and double distilled water added to a final volume of 20 μL.

The real-time PCR amplification conditions were as follows: 1 cycle of initial denaturation at 95°C for 3 min and 40 cycles consisting of denaturation at 94°C for 20 s, annealing at 56°C for 30 s and extension at 72°C for 30 s (fluorescence signal detection was conducted under this condition).

To generate standard curves, tenfold serial dilutions of Pst DNA (10, 1, 10−1, 10−2, 10−3, and 10−4 ng/μL) and wheat DNA (100, 10, 1, 10−1, and 10−2 ng/μL) were made. The real-time PCR amplifications were conducted as described above, and the corresponding cycle threshold (Ct) values were recorded. Then two standard curves to quantify Pst DNA and wheat DNA from the samples consisting of the latently Pst-infected wheat leaves after acquisition of near-infrared spectral data were generated.

Using the extracted DNA from the latently Pst-infected samples as template DNA, the real-time PCR amplifications were conducted as described above, and the corresponding Ct value of each sample was obtained. According to the standard curves, the contents of Pst DNA and wheat DNA were quantitatively determined. And the relative content of Pst DNA could be calculated using the following formula: , where RCP is the relative content of Pst DNA, CP is the content of Pst DNA (ng), and CW is the content of wheat DNA (ng).

2.6. Establishment of Determination Models for Quantification of the Relative Content of Pst DNA in Wheat Leaves in the Incubation Period

The obtained near-infrared spectra of the latently Pst-infected wheat leaf samples were preprocessed by using three methods including multiplication scatter correction (MSC), standard normalized variate (SNV), and vector normalization (VN). The 300 obtained spectra of the samples were divided into training set and testing set based on a ratio of the training set to the testing set equal to 3 : 1. By the combined use of quantitative partial least squares (QPLS) and support vector regression (SVR), the kQPLS-SVR models to quantify the relative content of Pst DNA in wheat leaves in the incubation period were built based on the original near-infrared spectra and the spectral data obtained by using the three preprocessing methods.

As shown in Figure 2, to build a kQPLS-SVR model, firstly, m features were randomly selected from the spectral features (spectral attributes) of 2100 wavelength points, the number of principal components was set as n, and a QPLS model was built. Using this method, k QPLS models were built; secondly, a SVR model was built by using the predicted values of the k QPLS models as variable. Thus, a kQPLS-SVR model was obtained. The SVR model was built with radial basis function as the kernel function. Using the grid search algorithm, the penalty parameter C and the kernel function parameter g for the SVR model were optimized in a range of 2−8–28 with the searching step of 0.8. As the minimum mean squared error of the training set was achieved at a point within the grid, the corresponding values of C and g were regarded as the optimal parameters for building the SVR model. In this study, the number of the randomly selected spectral features (m) was set as 700 or 1400, the number of principal components (n) while building a QPLS model was set as 4, 8, or 12, and the number of the built QPLS models (k) was set as 5, 10, or 15. All the calculations and modeling processes described above were implemented using the software MATLAB 7.8.0 (R2009a) (MathWorks, Natick, MA, USA).

The coefficient of determination (R2), standard error of calibration (SEC), average absolute relative deviation (AARD) and relative prediction deviation (RPD) of the training set and R2, standard error of prediction (SEP), AARD, and RPD of the testing set were used to evaluate the kQPLS-SVR models built for quantification of the relative content of Pst DNA in wheat leaves in the incubation period. A value of R2 more than 0.5 denotes that the corresponding model can be used for rough screening and actual application [35]. The closer the value of R2 is to 1, the higher the accuracy of the model is. The accuracy of a model is also related to the value of RPD, and a higher value of RPD denotes that the model has greater prediction ability [35]. The less value of SEC, SEP, or AARD denotes that higher accuracy can be obtained using the model and that the model has greater prediction ability. According to these evaluation indicators described above, the selection of the optimal kQPLS-SVR model was conducted.

Determination models of the relative content of Pst DNA in wheat leaves in the incubation period were also built using individual methods including QPLS and SVR based on the same training set and testing set as used for building the optimal kQPLS-SVR model. A comparison of the effects of the three models was conducted according to R2, SEC, AARD, and RPD of the training set and R2, SEP, AARD, and RPD of the testing set. The SVR model was built and optimized as described above. The number of principal components used during building the QPLS model was determined by evaluating the prediction residual error sum of square (PRESS) [36]. Generally, when the minimum PRESS value is obtained, the value of the corresponding number of principal components may be the optimum, but it may lead to overfitting in this case. The value of PRESS can be calculated using the following formula [36]: where is the number of the samples in the training set, is the number of principal components used during modeling, is the predicted value of the ith sample of the training set, and is the actual value of the ith sample of the training set. In this study, the optimal number of principal components was determined using statistical method described as the following formula: , where is the number of principal components corresponding to the minimum PRESS value; the optimal number of principal components (f) is less than , and it should be as small as possible and satisfy the following condition: F(f) < Fα,m,m (α = 0.25, m is the number of degrees of freedom) [36].

3. Results

3.1. Quantitative Measurements of Pst DNA in Wheat Leaves in the Incubation Period Using Duplex TaqMan Real-Time PCR

After the real-time PCR amplifications were conducted with tenfold serial dilutions of Pst DNA, the standard curve for quantification of Pst DNA (Figure 3) was generated using common logarithmic value of the concentration of Pst DNA (lg(Pst DNA concentration)) as the ordinate and Ct value as the abscissa, and the corresponding linear regression equation was obtained. As shown in Figure 3, the equation of the standard curve for quantification of Pst DNA was as follows: y = −0.3019x+ 7.4752 (), where y is lg(Pst DNA concentration) and x is the Ct value. After the real-time PCR amplifications were conducted with tenfold serial dilutions of wheat DNA from healthy leaves, the standard curve for quantification of wheat DNA (Figure 4) was generated using common logarithmic value of the concentration of wheat DNA from leaves (lg(wheat leaf DNA concentration)) as the ordinate and Ct value as the abscissa, and the corresponding linear regression equation was also obtained. As shown in Figure 4, the equation of the standard curve for quantification of wheat DNA was as follows: y = −0.3008x+ 9.2660 (), where y is lg(wheat leaf DNA concentration) and x is the Ct value. Moreover, satisfactory amplification efficiencies of the real-time PCR arrays were obtained for both Pst DNA and wheat DNA, and they were 99.7% and 100.1%, respectively.

After the real-time PCR arrays were conducted with the extracted DNA from the latently Pst-infected wheat leaf samples, the contents of Pst DNA and wheat DNA in each sample were figured out based on the linear regression equations of the generated standard curves. The results showed that the relative content of Pst DNA increased exponentially with time (Figure 5). An equation to fit the daily changes of the relative contents of Pst DNA in wheat leaves after inoculation was built: y = 0.0096e0.9287x (), where y is the relative content of Pst DNA (%) in a wheat leaf sample after inoculation and x is the days post inoculation. In this study, attempts were made to build a linear regression model by using partial least squares method. The relative contents of Pst DNA in wheat leaves after inoculation changed in a range of 0.00385%–90.09%. Especially, all relative contents of Pst DNA in the latently Pst-infected wheat leaves in the first six days were no more than 0.2%. To make the data follow a normal distribution, a logarithmic transformation of the relative contents of Pst DNA in wheat leaves after inoculation was conducted in this study. Meanwhile, to ensure that each value after transformation was positive, each value of the relative contents of Pst DNA in wheat leaves after inoculation was multiplied by 105 before the logarithmic transformation. The linear regression model (y = 0.4036x+ 0.9861, where y is lg(relative content Pst DNA × 105) and x is the days post inoculation) built based on the common logarithmic values was shown in Figure 6. As shown in Figures 5 and 6, the coefficients of determination (R2) of the equations built with the days post inoculation as independent variable before and after logarithmic transformation were 0.8784 and 0.8750, respectively. It was demonstrated that there was no great difference between the values of R2 of the two equations. Therefore, in combination with the near-infrared spectral data, lg(CP × 100%/(CP + CW) × 105) in which CP was the content of Pst DNA in ng and CW was the content of wheat DNA in ng, that is, lg(relative content Pst DNA × 105) was treated as chemical value for further modeling to develop a determination model for quantification of the relative content of Pst DNA in wheat leaves in the incubation period.

3.2. Results of the kQPLS-SVR Models for Quantification of the Relative Content of Pst DNA in Wheat Leaves in the Incubation Period

To quantitatively determine relative content of Pst DNA in wheat leaves in the incubation period using NIRS technology, the results of the kQPLS-SVR models built with original near-infrared spectra were shown in Table 2. The results demonstrated that the best effects were achieved for the kQPLS-SVR model built when the number of the spectral attributes of the randomly selected wavelength points was 700, the number of principal components was 8, and the number of the built QPLS models was 5. For this model, the R2 value, SEC, AARD, and RPD of the training set were 0.9027, 0.3852, 0.1257, and 3.2057, respectively, and the R2 value, SEP, AARD, and RPD of the testing set were 0.8820, 0.4349, 0.1274, and 2.9107, respectively.

The results of the kQPLS-SVR models built based on the data obtained after preprocessing of original near-infrared spectra by using MSC were shown in Table 3. The model built when the number of the spectral attributes randomly selected was 700, the number of principal components was 8, and the number of the built QPLS models was 15 was better than others among the kQPLS-SVR models as shown in Table 3. For this model, the R2 value, SEC, AARD, and RPD of the training set were 0.9227, 0.3434, 0.1127, and 3.5961, respectively, and the R2 value, SEP, AARD, and RPD of the testing set were 0.8730, 0.4510, 0.1361, and 2.8062, respectively.

When the original near-infrared spectra were preprocessed by using the method SNV, the results of the built kQPLS-SVR models were shown in Table 4. Among the kQPLS-SVR models as shown in Table 4, the better effects were obtained for the model built when the number of the spectral attributes randomly selected was 1400, the number of principal components was 4, and the number of the built QPLS models was 15. For this model, R2, SEC, AARD, and RPD of the training set were 0.8963, 0.3976, 0.1117, and 3.1057, respectively, and R2, SEP, AARD, and RPD of the testing set were 0.8870, 0.4255, 0.1319, and 2.9748, respectively.

The results of the built kQPLS-SVR models based on the data obtained when VN was used as the preprocessing method were shown in Table 5. As demonstrated in Table 5, the effects of the model built when the number of the spectral attributes randomly selected was 700, the number of principal components was 8, and the number of the built QPLS models was 15 was better than others. For this kQPLS-SVR model, the value of R2, SEC, AARD, and RPD of the training set were 0.8232, 0.5193, 0.1814, and 2.3782, respectively, and the value of R2, SEP, AARD, and RPD of the testing set were 0.7964, 0.5711, 0.1833, and 2.2164, respectively.

The results shown in Tables 25 indicated that satisfactory effects could be obtained using the kQPLS-SVR models built based on the original near-infrared spectra and the spectral data obtained by using MSC and SNV. In contrast, for the built kQPLS-SVR models when the original near-infrared spectra were preprocessed by using the method VN, the values of both R2 and RPD of the training set and the testing set were less. For the optimal kQPLS-SVR model built based on the data obtained when SNV was used as the preprocessing method, the number of principal components was 4. It was less than the number of principal components used for building the optimal kQPLS-SVR model based on the original near-infrared spectra, the spectral data obtained by using MSC, or the spectral data obtained by using VN. In this case, there may be underfitting problem resulting in the reduction of the prediction ability of the optimal kQPLS-SVR model built when SNV was used as the preprocessing method. For the optimal kQPLS-SVR model built based on the data obtained when MSC was used as the preprocessing method, the value of R2 of the testing set was 0.8730. And it was less than the value of R2 of the testing set resulting from the optimal kQPLS-SVR model built based on the original near-infrared spectra. Therefore, the optimal kQPLS-SVR model built based on the original near-infrared spectra (the number of the spectral attributes of the randomly selected wavelength points was 700, the number of principal components was 8, and the number of the built QPLS models was 5) was regarded as the optimal kQPLS-SVR model to quantitatively determine the relative content of Pst DNA in wheat leaves in the incubation period.

3.3. Results of the QPLS Models and the SVR Model for Quantification of the Relative Content of Pst DNA in Wheat Leaves in the Incubation Period

Since the selected optimal kQPLS-SVR model to quantitatively determine the relative content of Pst DNA in wheat leaves in the incubation period was built based on the original near-infrared spectra, the QPLS models and the SVR model were built with the original near-infrared spectra based on the same training set and testing set as used for building the optimal kQPLS-SVR model. When the determination models of the relative content of Pst DNA in wheat leaves in the incubation period were built using QPLS, the number of principal components corresponding to the minimum PRESS value was calculated, and it was 8. The value of Fα,m,m was calculated using the software IBM SPSS Statistics 19.0 (IBM Corporation, Somers, NY, USA), and it was 1.09. When the number of principal components was 6, the value of F(f) was 1.07. And the value of F(f) was 1.02 when the number of principal components was 7. The two values of F(f) were less than the value of Fα,m,m. So the optimal number of principal components was set as 6 for building the QPLS model to quantitatively determine the relative content of Pst DNA in wheat leaves in the incubation period. As shown in Table 6, for the optimal QPLS model, the value of R2, SEC, AARD, and RPD of the training set were 0.8534, 0.4728, 0.1689, and 2.6117, respectively, and the value of R2, SEP, AARD, and RPD of the testing set were 0.8684, 0.4592, 0.1467, and 2.7565, respectively. When C = 256 and g = 0.1895, the optimal SVR model for quantification of the relative content of Pst DNA in wheat leaves in the incubation period was obtained. As shown in Table 6, for this SVR model, the value of R2, SEC, AARD, and RPD of the training set were 0.8887, 0.4119, 0.1277, and 2.9981, respectively, and the value of R2, SEP, AARD, and RPD of the testing set were 0.7971, 0.5702, 0.2110, and 2.2200, respectively.

A comparison of the effects of the three models, including the optimal kQPLS-SVR model, the optimal QPLS model, and the optimal SVR model, was conducted according to R2, SEC, AARD, and RPD of the training set and R2, SEP, AARD, and RPD of the testing set. The results showed that the effects of the optimal kQPLS-SVR model built based on the original near-infrared spectra (the number of the spectral attributes of the randomly selected wavelength points was 700, the number of principal components was 8, and the number of the built QPLS models was 5) were the best. Therefore, this optimal kQPLS-SVR model was selected as the optimal model to quantitatively determine the relative content of Pst DNA in wheat leaves in the incubation period.

4. Discussion

For the kQPLS-SVR models, the optimal QPLS model, and the optimal SVR model as shown in Tables 26, all the values of R2 of the training set and the testing set were more than 0.5. The results indicated that the correlation between the features of near-infrared spectra (absorbances) and the relative content of Pst DNA in wheat leaves in the incubation period was relatively high and that the built models could be used for rough screening. It was indicated that quantitative determination of the relative content of Pst DNA in wheat leaves in the incubation period using NIRS technology is feasible. After a comprehensive comparison of the determination results of the models as shown in Tables 26, the kQPLS-SVR model built with a ratio of the training set to the testing set equal to 3 : 1 based on the original near-infrared spectra in a range of 4000–12,000 cm−1 when the number of the spectral attributes of the randomly selected wavelength points was 700, the number of principal components was 8, and the number of the built QPLS models was 5 was regarded as the optimal model to quantitatively determine the relative content of Pst DNA in wheat leaves in the incubation period. The satisfactory effects were obtained using this optimal model. For this model, the value of R2, SEC, AARD, and RPD of the training set were 0.9027, 0.3852, 0.1257, and 3.2057, respectively, and the value of R2, SEP, AARD, and RPD of the testing set were 0.8820, 0.4349, 0.1274, and 2.9107, respectively. The results indicated that the proposed method based on NIRS could be used as a new method for rapid, nondestructive, and quantitative detection of Pst DNA content in wheat leaves. Meanwhile, a reference was provided for nondestructive determination of the DNA contents of other kinds of pathogens in plant hosts.

Changes of the near-infrared spectra of wheat leaves infected with Pst during the incubation period may be induced by many factors. The infection of Pst and subsequent expansion in plant hosts are dynamic change processes [1, 5]. After germination of a Pst urediospore, a growing germ tube emerges. Then an appressorium is formed at the tip of the germ tube. An infection peg grows from the appressorium and then enters the wheat leaf. During the incubation period, with the continuous spread of hyphae in the infected leaf tissues, changes of the photosynthesis, respiration, and transpiration of wheat leaves occur, affecting the synthesis and decomposition of organic substances [1, 3739]. Thus, the near-infrared spectra of the infected wheat leaf are influenced. Especially in the middle and later stage of the incubation period, a large quantity of Pst hyphae accumulate in the infected leaf tissues, the latent lesions on the infected leaves become more and more obvious, and the effects on various physiological indexes of wheat are increasing [1], more directly affecting the near-infrared spectra of the infected wheat leaves. In this study, according to the changes of the near-infrared spectra of the infected wheat leaves during the incubation period, the detection of the relative content of Pst DNA in the latently Pst-infected wheat leaves was performed using NIRS technology in combination with the duplex TaqMan real-time PCR arrays, providing a basis for modeling the quantitative relationship between the relative content of Pst DNA in the latently Pst-infected wheat leaves and the corresponding near-infrared spectral data.

The quantitative detection results obtained by using the duplex TaqMan real-time PCR arrays in this study showed that the quantity of Pst DNA exponentially increased as Pst extended in wheat leaf during the incubation period of wheat stripe rust, and it was consistent with the growth trend of Pst DNA reported by Pan et al. [10]. The results in this study showed that the relative contents of Pst DNA in wheat leaves after inoculation continuously increased. The relative contents of Pst DNA in the latently Pst-infected wheat leaves in the incubation period were more than 20% on 7 days post inoculation (dpi) and reached to approximately 90% on 10 dpi. It was demonstrated that there was a particularly obvious increase of the relative contents of Pst DNA in the latently Pst-infected wheat leaves from the seventh day to the tenth day during the incubation period. This indicated that it would have great significance to perform early detection and control of wheat stripe rust.

Wheat stripe rust is a kind of air-borne disease. As disease symptoms appear on wheat leaves, a large number of urediospores will be produced and serve as inocula for the further spread and epidemic of the disease [1, 35]. Accurate quantitative determination of the latent infection levels of wheat leaves by Pst can provide important information for estimating potential inoculum level. Early detection of Pst infection, early disease warning, and early disease control measure-making can effectively reduce the amount of inoculum sources and the dosage of pesticides. The results of this study demonstrated that latent infection of wheat seedlings caused by Pst could be detected as early as 24 h after inoculation based on NIRS. For the first time, dynamic changes of the amounts of wheat stripe rust pathogen in the infected leaves during the incubation period were investigated by using NIRS technology in combination with a real-time PCR method. A modeling method integrated with QPLS and SVR was used to establish the dynamic quantitative detection model of wheat stripe rust pathogen during the incubation period. A reference was provided for prediction of disease epidemic trend and early control of wheat stripe rust in this study.

5. Conclusions

It is critical to implement early quantitative detection of Pst and disease prediction accurately for the prevention and control of wheat stripe rust. In this study, a method based on NIRS to implement the quantitative determination of the relative content of Pst DNA in wheat leaves in the incubation period was investigated. Based on the data on the relative content of Pst DNA in wheat leaves in the incubation period obtained using the duplex TaqMan real-time PCR arrays and the corresponding near-infrared spectral data, the optimal model integrated with QPLS and SVR was obtained to quantify the relative content of Pst DNA in wheat leaves in the incubation period. Using this model, satisfactory results were achieved. The results indicated that rapid, nondestructive, and quantitative detection of the amount of Pst in the infected leaves during the incubation period could be realized using NIRS technology. A novel method for quantitatively estimating the latent infection levels of Pst in wheat leaves and early detection of wheat stripe rust was provided in this study. Furthermore, this proposed method can be used as a reference to establish models for determining the latent amounts of other kinds of pathogens in plant hosts and implement early nondestructive detection of the diseases caused by these agents. It is very helpful in getting the information about disease prevalence and inoculum amount in the field as early as possible, inhibiting pathogen accumulation, reducing the quantity of propagules, performing early precise site-specific control, disease forecasting, and making macro strategy for disease management. And it is also conducive to reducing the dosage of pesticides and increasing efficiency of control measures in plant disease management. Moreover, in this study, some basis was provided for the development of portable near-infrared spectrometers and sensors for quantitative detection of plant diseases.

Conflicts of Interest

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contributions

Yaqiong Zhao, Yilin Gu, and Feng Qin contributed equally to this paper.

Acknowledgments

This study was supported by the National Key Basic Research Program of China (2013CB127700), National Natural Science Foundation of China (31471726), and National Key Technologies Research and Development Program of China (2012BAD19BA04).