A Mechanistic Model of Botrytis cinerea on Grapevines That Includes Weather, Vine Growth Stage, and the Main Infection Pathways

A mechanistic model for Botrytis cinerea on grapevine was developed. The model, which accounts for conidia production on various inoculum sources and for multiple infection pathways, considers two infection periods. During the first period (“inflorescences clearly visible” to “berries groat-sized”), the model calculates: i) infection severity on inflorescences and young clusters caused by conidia (SEV1). During the second period (“majority of berries touching” to “berries ripe for harvest”), the model calculates: ii) infection severity of ripening berries by conidia (SEV2); and iii) severity of berry-to-berry infection caused by mycelium (SEV3). The model was validated in 21 epidemics (vineyard × year combinations) between 2009 and 2014 in Italy and France. A discriminant function analysis (DFA) was used to: i) evaluate the ability of the model to predict mild, intermediate, and severe epidemics; and ii) assess how SEV1, SEV2, and SEV3 contribute to epidemics. The model correctly classified the severity of 17 of 21 epidemics. Results from DFA were also used to calculate the daily probabilities that an ongoing epidemic would be mild, intermediate, or severe. SEV1 was the most influential variable in discriminating between mild and intermediate epidemics, whereas SEV2 and SEV3 were relevant for discriminating between intermediate and severe epidemics. The model represents an improvement of previous B. cinerea models in viticulture and could be useful for making decisions about Botrytis bunch rot control.


Introduction
Botrytis cinerea Pers.: Fr., the anamorph of Botryotinia fuckeliana (de Bary) Whetzel, is the causal agent of grey mould, which is a worldwide disease that causes serious losses in more than 200 host species [1,2]. Grey mould (also called Botrytis bunch rot) is one of the most important diseases of grapevines [3].
B. cinerea infection pathways differ for conidia vs. mycelium; conidia infect inflorescences, young clusters, and ripening berries while mycelium is responsible for berry-to-berry infection [3]. Grape inflorescences are more susceptible at flowering (beginning, full, and end of flowering) than at earlier growth stages or at fruit swelling or berries groat-sized stages [14]. On inflorescences and young grape clusters, infection severity increases with hours of wetness, and the optimal temperature for infection is about 20°C [15][16][17]. The susceptibility of berries increases from veraison to ripening [18], and at ripening B. cinerea can infect all berries in a cluster and cause heavy crop losses. Infection incidence in mature grape berries is higher at temperatures between 15 and 25°C than at other temperatures and increases with increasing hours of wetness or high relative humidity [15][16][17]19,20]. Moreover, disease incidence is higher on wounded than on unwounded berries [19,[21][22][23], especially when wetness is short or humidity is low [17].
Control of B. cinerea is difficult because: (i) the pathogen can produce large numbers of conidia on multiple inoculum sources (overwintered grape and herb debris, bunch and leaf trash, and rotted berries) [24][25][26][27][28][29]; (ii) grapevines are susceptible at multiple growth stages [14,30]; (iii) different infection pathways exist [3]; and (iv) infection can occur under a range of environmental conditions, and these conditions differ depending on the infection pathway [14,17,30]. This complexity has caused growers to rely heavily on routine application of fungicides at four specific grape growth stages: A, end of flowering (growth stage 69 of Lorenz et al. [31]; B, pre-bunch closure (growth stage 77); C, veraison (growth stage 83); and D, before harvest (before growth stage 89) [2,20,[32][33][34]. This fungicide scheduling may result in unnecessary sprays [35,36], which is unacceptable because of environmental and public health considerations [37,38] and because of the increased risk of fungicide resistance [39]. There is therefore the need to improve B. cinerea control in vineyards by applying fungicides only when necessary and by following suitable anti-resistance strategies [39][40][41].
To reduce the risk of unnecessary sprays, different weather-based methods have been developed. The "15-15 empirical rule" [32,42] and the French EPI-Botrytis model [43] were in some cases associated with a reduction in the number of fungicide applications, but experiments carried out in different countries and seasons were inconsistent [34,[43][44][45][46]. Other weather-based models have been developed by Broome et al. [20], Nair and Allen [15], and Rodríguez-Rajo et al. [47], but these models lack of robust field validation. An expert system was also developed in Australia, which estimates multiple infection risks based on the available knowledge about B. cinerea biology, economic thresholds, and fungicides [48]. The expert system was validated in 1990 to 1994 but it was unable to improve the grower's practice [49].
To date, none of the above models is extensively used by viticulturists. One reason could be the failure of these methods to account for the complexity of B. cinerea epidemiology, as briefly explained in the previous paragraph (i.e., high sporulation potential, susceptibility of the different grape growth stages, multiple infection pathways, infection occurring under a wide range of environmental conditions, and interactions between environmental conditions and infection pathway). It follows that a method for predicting infection risk by B. cinerea is still needed.
A new model predicting risk of grapevine infection by B. cinerea was developed in this research. This model is mechanistic and attempts to account for the full complexity of the B. cinerea life cycle and of Botrytis bunch rot epidemiology. The model was developed using the results of recent publications that investigated the effect of environmental conditions on the biology and epidemiology of B. cinerea isolates belonging to different transposon genotypes [14,17,30]. The model was validated against 21 independent Botrytis bunch rot epidemics in Italy and France. one were determined by rate variables depending on environmental conditions and host growth stages, as shown in the relational diagram of the model (Fig 1 and Table 1). The model calculations begin when grape inflorescences are clearly visible and ends when berries are ripe for harvesting, with a time step of 1 day.

Model description
The model considers that multiple inoculum sources are present in the vineyard all season. These inoculum sources include dead, infested plant material; affected inflorescences and young berries; infested floral debris; and affected ripening berries (Fig 1). The relative abundance of conidia on these sources (CISO) on any day i of the grape-growing season depends on the rate at which the mycelium grows and saprophytically colonises the source tissue (MYGR) and on the rate of spore production (SPOR), as follows: with MYGR n ¼ 3:78 Â Teq 0:9 n Â ð1 À Teq n Þ À Á 0:475 Â Mf n ð2Þ SPOR n ¼ 3:7 Â Teq 0:9 n Â ð1 À Teq n Þ À Á 10:49 Â À3:595 þ 0:097 Â RH n À 0:0005 Â RH 2 In these equations, Teq = temperature equivalent in the form Teq = (T-Tmin) / (Tmax-Tmin); with T = daily average temperature (°C); Tmin = minimum temperature for mycelial growth or sporulation (0°C); Tmax = maximum temperature for mycelial growth (40°C) or sporulation (35°C); and Mf = factor accounting for the effect of moisture of the medium (e.g., bunch trash) on mycelial growth. With respect to Mf, B. cinerea mycelium grows quickly when the water activity of the medium is a w !0.95 (unpublished data), which corresponds to residue moisture greater than 60 to 80% depending on temperature and sorption or desorption conditions [70]. We therefore assumed that residues contain sufficient humidity for mycelial growth in the hours of the day when rain, R!0.2 mm or wetness duration, WD!30 min or RH (daily average relative humidity)!90% (i.e., moist hours); therefore, Mf = number of moist hours/24.
Eqs (2) and (3) were derived from Ciliberti et al. [14] and Ciliberti et al. [30], respectively. These equations were averaged over a 7-day period because researchers previously reported that this is the duration for sporulation of a single lesion [30,71].
The model assumes that, on any day, there are favourable conditions for conidia to disperse and settle on host plant surfaces, so that DISRi = 1 and DEPRi = 1 (Fig 1). This assumption is consistent with our knowledge about B. cinerea epidemiology, as described in the section "Life cycle of Botrytis cinerea".
The model considers two main infection periods (or infection windows) corresponding to two grape-growing periods (Fig 1). The first infection window occurs between "inflorescences clearly visible" (stage 53 of the scale of Lorenz et al. [31] and "berries groat-sized, bunches begin to hang" (stage 73); in this period the model calculates infection by conidia on inflorescences and young clusters (pathway I and II). The second infection window occurs between "majority of berries touching" (stage 79) and "berries are ripe for harvest" (stage 89); in this period, the model calculates infection on ripening clusters by conidia (pathway Va) and berryto-berry infection by aerial mycelium (pathway Vb). Compared to the relational diagram of Fig  1, the model described in this paper does not consider latent infections and incubation periods (latent growth and INC in Fig 1) because there is insufficient quantitative information about the effect of weather and host conditions on latency and incubation.
In the first infection window (stages 53 to 73), the model calculates an infection rate on inflorescences and young clusters (INF1) as:  Table 1 for acronym explanations.
where Teq = temperature equivalent as described for eq (2), with Tmin = 0°C and Tmax = 35°C; WD = wetness duration (in hours); and In the latter equation, SUS1 = relative susceptibility of the inflorescences and young clusters, and GS = growth stage of the plant based on the stages of the scale of Lorenz et al. [31].
Relative infection severity in the first infection window is then calculated as: In the second infection window (stages 79 to 89), the model calculates two infection rates on ripening berries: one for conidial infection (INF2) and another for berry-to-berry infection (INF3) (Fig 1).
Infection rate for conidial infection is calculated as follows: where Teq = temperature equivalent as described for eq (2) with Tmin = 0°C and Tmax = 35°C, and In the latter equation, SUS2 = relative susceptibility of the ripening berries to B. cinerea conidia infection, and GS = growth stage of the plant based on the scale of Lorenz et al. [31].
Relative infection severity for conidial infection is then calculated as follows: Infection rate for berry-to-berry infection during the second infection window is calculated as follows: where Teq = temperature equivalent as described in eq (2) with Tmin = 0, and Tmax = 30°C, and In the latter equation, SUS3 = relative susceptibility of the ripening berries to B. cinerea mycelium infection, and GS = growth stage of the plant based on the scale of Lorenz et al. [31].
Eq (10) was developed by Ciliberti et al. [17], and eq (11) was developed based on one experiment from Ciliberti et al. [17] as described in S1 Appendix.
Relative infection severity for berry-to-berry infection is then calculated as follows: Daily values of relative infection severity (i.e., RIS1, RIS2, and RIS3) are finally accumulated over the time of the infection window that they refer to. These accumulated values produce new variables (named SEV1, SEV2, and SEV3, respectively), which provide a picture of the total risk of infection. Two examples of model output are shown in Figs 2 and 3 (epidemics at CO-09 and SM-12 in Table 2, respectively) as follows: weather data are shown in Figs 2A and 3A, sporulation on inoculum sources in Figs 2B and 3B, and the daily and accumulated values of relative infection severity in Figs 2C and 3C. The considered period (i.e., between grape growth stages 53 and 89) was 153 days long at CO-10 and 129 days long at SM-12. At both locations, > 350 mm of rain fell in this period. This rain was more regularly distributed at CO-10 than at SM-12, and the number of hours with wetness was higher at CO-10 than at SM-12 (1060 vs. 530 h, respectively), especially during the first infection window (619 h for CO-10 and 254 h for SM-12) (Figs 2A and 3A). As a consequence, sporulation was abundant and regular at CO-10 during the first infection window, and a constant increase of RIS1 was observed, with an accumulated  [31]), fruit set (stage 71), and berries developing (stage 83). The first infection window extends from inflorescences clearly visible" (stage 53) to "berries groat-sized, bunches begin to hang" (stage 73). The second infection window extends from "majority of berries touching" (stage 79) to "berries are ripe for harvest" (stage 89). value SEV1 = 2.64 at the end of this period (Fig 2B and 2C). In the same period at SM-12, the model predicted lower levels of sporulation and lower values of RIS1, with only two peaks at the end of April and the beginning of June, resulting in an accumulated value SEV1 = 1.10 at the end of the first infection period (Fig 3B and 3C). During the second infection period, the  [31]), fruit set (stage 71), and berries developing (stage 83). The first infection window extends from "inflorescences clearly visible" (stage 53) to "berries groat-sized, bunches begin to hang" (stage 73). The second infection window extends from "majority of berries touching" (stage 79) to "berries are ripe for harvest" (stage 89). Epidemiological Model of Botrytis cinerea on Grapevines wetness duration decreased in both vineyards, and lower values of sporulation were predicted (Figs 2 and 3). At CO-10, high values of RIS2 and RIS3 were predicted only in early August and in the last 2 weeks before harvest. At SM-12, the model predicted favourable conditions for B. cinerea 1 week before harvest. Consequently, the values of SEV2+SEV3 were higher for CO-10 (2.61) than for SM-12 (0.5).

Model validation Field assays
The model output was validated against real data obtained from 2009 to 2014 from 21 Botrytis bunch rot epidemics in untreated plots in 12 experimental vineyards. These data have not been previously used for model development. Ten of these vineyards were in Italy, and two were in France (Table 2), under Cfb and Cfc climate types, respectively (C = warm temperate; f = fully humid precipitation; b = warm summer; c = cold summer) [72]. The vineyards were cropped with grape varieties highly susceptible to B. cinerea [48,57,73] and were managed as usual for the viticultural area, with the exception that no fungicides were used to control B. cinerea. Standard, electronic weather stations were installed in the vineyard borders (with sensors at 1.5 m above the ground) to measure temperature, relative humidity, wetness duration, and rainfall on an hourly basis. Growth stages of vines were periodically assessed in the vineyards according of the scale of Lorenz et al. [31]. Weather and vine growth stages were used as model inputs to calculate model output (i.e., SEV1, SEV2, and SEV3). At full ripening (growth stage 89), disease was assessed on a minimum of 100 random bunches per plot (on at least 20 vines per plot) and in at least four replicate plots per vineyard. Disease incidence (DI) and severity (DS) were visually assessed as the percentage of bunches with Botrytis rot and as the percentage of the surface of the bunch affected by rot and showing typical B. cinerea sporulation, respectively. Based on disease assessment, the 21 epidemics were placed into three groups based on the following classification rules: severe, DI ! 75% (and DS ! 15%); intermediate, 74 ! DI ! 25% (and DS < 15%); and mild, DI < 24% (and DS < 5%). Threshold values for these groups were selected based on the literature and personal communications from grape growers, advisors, and wineries. Epidemic classified as mild (DS<5%) result in yield losses the growers are willing to accept; disease severities of 3% [74], 5% [75] or 6% [76] were considered thresholds for fungicide application. The need of applying fungicides with intermediate epidemics (DS<15%) depends on the risk perception by the grower [77]; risky growers prefer to avoid treatments as they consider the risk of yield losses as acceptable. Finally, no growers renounce to protect the crop when the expected DS is ! 15%.

Ethics statements
The owners of the vineyards used for model validation gave permission to conduct the study on these sites.

Data analysis
A discriminant function analysis (DFA) was performed to determine whether SEV1, SEV2, and SEV3 (discriminating variables) were effective in predicting whether epidemics fell into one of three groups: severe, intermediate, or mild [78][79][80][81]. The analysis created two linear combinations of discriminant functions in the form: F n = a n + b n,1 Áln(SEV1+1) + b n,2 Áln(SEV2+SEV3 +1), where n = 1 or 2 for the two discriminant functions, respectively. Briefly, the procedure automatically selected a first function (F 1 ) that was able to separate the groups as much as possible; a second function (F 2 ) was then selected that was both uncorrelated with the first function and that provided as much further separation as possible. For each epidemic group, the centroid (i.e., the mean value of each group) was calculated.
The two functions were used to: (i) locate the cases (epidemics) into the F 1 -F 2 space based on the coordinates obtained by solving the two discriminant functions; (ii) define the case distance from the centroid of each group (as the Mahalanobis distance; [81]); (iii) determine the probability that each case belonged to each of the three groups and determine group membership based on these probabilities; and (iv) obtain information on the effect of the discriminant variables on group membership. For (iii), prior probability of membership was set proportional to the number of epidemics in each group (i.e., 0.19 for mild, 0.48 for intermediate, and 0.33 for severe epidemics). For (iv), the canonical coefficients (CCs), standardized canonical coefficients (SCCs), and the correlation coefficients (COCs) of the two functions were calculated. The magnitude of CCs and SCCs is an indicator of the weight of each variable in each of the discriminant functions; COCs indicate the discriminant power of each variable in these functions.
The Wilks' lambda was calculated as a measure of how well each function separated cases into the three groups; it is equal to the proportion of the total variance in the discriminant scores not explained by differences among the groups. Lambda ranges from 0 to 1, with 0 indicating that group means differ (i.e., that the function strongly differentiates groups), and 1 indicating that all group means are the same. A chi-square statistic was used to test the hypothesis that the means of the functions listed were equal across groups (a small significance value for this test indicates that the discriminant function does better than chance at separating the groups) [80].
DFA was performed using the discriminant procedure of the statistical software SPSS (ver. 21.0, IBM SPSS Statistics, IBM Corp., New York, USA).

Evaluation of model predictions during the season
DFA enables prediction of new cases based on the previously established structure [82]. In this section, DFA was evaluated to anticipate the prediction of the final severity of the B. cinerea epidemic (i.e., mild, intermediate, or severe) by using the daily model output, i.e., the daily values of SEV1, SEV2, and SEV3.
The daily model outputs (the natural logarithms of daily values of SEV1 and SEV2+SEV3) were used to solve the two canonical discriminant functions (F 1 and F 2 ) by using the non-standardized canonical coefficients calculated in the previous section. This allowed us to locate the cases in the F 1 -F 2 space and to calculate the standardized distance (D uj ) of each case (u) to the centroid of each group (j) (i.e., mild, intermediate, or severe) by using the Mahalanobis distance [81]. Because the DFA considered three groups, a multivariate approach was used: the different cases and centroids were represented in the F 1 -F 2 space as vectors, and matrix operations were used to resolve the equations [81]. The Mahalanobis distance was calculated in the form: where x u is the vector of case u, σ j is the vector of the centroid of population j, and X À1 j is the inverse of the covariance matrix of population j.
To assign each case u to one of the j groups, the posterior probability of membership P(j|x u ) was calculated as follows [81]: Where f(x u |j) is the likelihood of x u = j, and α j is the prior probability of group membership (i.e., the proportion of the 21 vineyards used for model validation belonging to each epidemic group as described in the previous section). Assuming a normal distribution for the population of each j group, we defined the density value f(x u |j) with a multivariate approach as described by Huberty and Olejnik [81]: Where |∑ j | is the determinant of the covariance matrix of population, and D uj is the Mahalanobis distance of each case u to the centroid of each group j.
With this procedure, the probability that the B. cinerea epidemic would be mild (P mild ), intermediate (P int ), or severe (P sev ) was calculated for each day and for each of the 21 epidemics described in the "field assays" section. For each day between growth stages 53 and 89, each epidemic was assigned to one of the three severity groups by using the maximum-likelihood principle, i.e., by selecting the group with the highest probability of occurrence [81].
Model predictions for these epidemics are shown in Fig 4. For mild epidemics, the average values of SEV1, SEV2, and SEV3 were 0.94 ± 0.13, 0.15 ± 0.06, and 0.34 ± 0.06, respectively. The lowest values were calculated at Z-12 (Fig 4), where the Botrytis bunch rot incidence at harvest was 1.3% (Table 2). Average values of SEV1 were higher for intermediate (2.46 ± 0.17) and severe epidemics (1.81 ± 0.35) than for mild epidemics (0.94 ± 0.13). Values of SEV2 and SEV3 were higher for severe epidemics (0.28 ± 0.06 and 0.76 ± 0.19, respectively) than for  Table 2; their symbols (dots) are green, orange or red based on Botrytis bunch rot incidence and severity as indicated in Table 2 (green, orange, and red indicate mild, intermediate, and severe epidemics, respectively). Triangles represent the average of each group. SEV1 = relative severity of the infections caused by conidia and accumulated over the first infection period ("inflorescences clearly visible" to "berries groat-sized, bunches begin to hang"); SEV2 = relative severity of the infections caused by conidia and accumulated over the second infection period ("majority of berries touching" to "berries are ripe for harvest"); and SEV3 = relative severity of the infections caused by mycelium (berry-to-berry infection) and accumulated over the second infection period. and SEV3 were calculated at C0-10, where Botrytis bunch rot incidence was 81.8%, and severity was 15.9% at harvest ( Table 2).
DFA correctly classified 17 of the 21 epidemics (i.e., 81% of correct classification) ( Table 3). In three of the four misclassifications, the real epidemic was more severe than the predicted one: epidemics at RAV-12, VI-14, and CO-13 were all severe based on disease assessment but were classified as mild or intermediate (Fig 5). In contrast, the epidemic at RUS-10 was intermediate but was classified as severe by the DFA (Fig 5). Using cross-validation, 15 of the 21 epidemics (71.4%) were correctly classified (Table 3).
In this separation of epidemics into three groups, F 1 accounted for 82.2% of the variation and F 2 accounted for the remaining 17.2% (Table 4); canonical correlation coefficients   Table 2; the color of the code indicates the severity of Botrytis bunch rot as assessed in the vineyard at harvest (see Table 2). The color of the dot indicates group assignment based on DFA. For both codes and dots, green, orange, and red indicate mild, intermediate, and severe epidemics, respectively. Crosses represent the centroids of each group based on DFA. Coefficients of the canonical discriminant functions are shown in Table 6.
doi:10.1371/journal.pone.0140444.g005 (Table 4) and the significance of Wilks' lambda tests ( Table 5) all indicated that F 1 had greater effect on group separation than F2. Canonical coefficients (CC), standardized canonical coefficients (SCC), and correlation coefficients (COC) collectively indicated that SEV1 had a major effect in distinguishing the groups in F 1 and that SEV2+SEV3 had a major effect in distinguishing the groups in F 2 ( Table 6).   indicates that all group means are the same. b Chi-square statistic used to test the hypothesis that the means of the functions listed are equal across groups. c A small significance value indicates that the discriminant function does better than chance at separating the groups.
doi:10.1371/journal.pone.0140444.t005 Table 6. Coefficients for each discriminant variable in each canonical function (F 1 or F 2 ) obtained in the discriminant analysis used to classify 21 Botrytis cinerea epidemics based on model output. intermediate starting from the first week of May (at CO-13, VI-14, and COT-13) to the first of June (CO-11 and Z-11), which was between 47 and 13 days before vine growth stage 69 (i.e., end of flowering) (Fig 6). For the epidemic at Z-11, the assignment of intermediate severity occurred at stage 69 (Fig 6). The change from mild to intermediate epidemics always occurred during the first infection period of the model (i.e., growth stages 53 to 73) (Fig 6), when the model calculates SEV1. In the five cases in which the epidemic was assigned to the group of severe epidemics (C-09, C-10, C-12, BAG-12, and RUS-10), the change from intermediate to severe occurred after growth stage 83 (veraison), between 3 days (CO-09) and 8 days (CO-12) before growth stage 89 (berries ripe for harvest) (Fig 6). At CO-12, the classification jumped from mild to severe.

Discussion
In this work, a new, mechanistic model for B. cinerea infections was elaborated [83]. In creating this model, we divided the B. cinerea life cycle into different state variables, and changes from one state to the following state were determined by rate variables that depended on environmental conditions and that were represented in mathematical equations. Equations were Epidemiological Model of Botrytis cinerea on Grapevines previously published on the effects of environmental variables and grape phenology on B. cinerea sporulation and on B. cinerea infection of inflorescences, young berries clusters, and mature berries [14,17,30]. The ability of these equations to correctly reproduce the biological processes investigated was demonstrated by their goodness of fit [14,17,30]. Although these equations were not validated against independent data, they were considered robust enough to be used in this model [83]. Because these previous papers showed that, the response of different strains belonging to the main transposon genotypes found in vineyards to environmental conditions was uniform, our model considered the B. cinerea population of a vineyard to be uniform. The model describes three key stages of the B. cinerea life cycle in vineyards: i) production of conidia on overwintered sources; ii) infection of inflorescences and young clusters by conidia and production of conidia on bunch trash (pathway I and II of Elmer and Michailides [3]); and iii) infection of mature berries by conidia (pathway Va) and by mycelium (i.e., berry-to-berry infection; pathway Vb), as well as production of conidia on affected, ripening berries.
On a daily basis, the model predicts and then accumulates the relative infection severity over two infection windows corresponding to the two grape-growing periods relevant for B. cinerea infection: i) between "inflorescences clearly visible" and "berries groat-sized, bunches begin to hang" and ii) ripening berries [3].
The model was validated in 21combinations of vineyard × year in different grape-growing areas of Italy and in France. A DFA analysis was used to evaluate the ability of the model output (SEV1, SEV2, and SEV3) to predict whether Botrytis bunch rot epidemics were mild, intermediate, or severe. DFA has been previously used in botanical epidemiology to compare epidemics and to evaluate the importance of multiple influencing variables on epidemic development [84][85][86]. Classification of the epidemics into the three groups was based on a prior probability of group membership, which was the real proportion of the 21 epidemics that were mild, intermediate, or severe. Under this assumption, a vineyard has a lower probability to develop a mild B. cinerea epidemic than an intermediate or severe epidemic (i.e., 0.19 vs. 0.48 and 0.33, respectively). The prior probabilities found in the data set based on observations of the 21 epidemics were consistent with those of a larger data set resulting from the analysis of 111 papers published in the journal Giornate Fitopatologiche between 1962 and 2012 (http:// fitogest.imagelinenetwork.com/it/giornate-fitopatologiche/). In this larger data set, 14% of the Botrytis bunch rot epidemics in untreated plots were mild, 59% were intermediate, and 26% were severe.
The DFA analysis correctly classified more than 80% of the 21 epidemics, indicating that the model output provided a reliable picture of reality. Three of the four misclassified epidemics were underestimated, i.e., the model predicted a lower level of disease than was observed in the field. Although the model predicts the development of Botrytis bunch rot epidemics based on grape phenology and environmental conditions, the severity of Botrytis bunch rot is influenced by other factors that were not included in the model. These other factors include wounds caused by abiotic factors [17,19,22,23,35] and by the grape moth Lobesia botrana [19,87,88], thrips [89], birds, and powdery mildew (Erysiphe necator) [35]. The severity of Botrytis bunch rot may also be influenced by cultural practices such as excessive nitrogen fertilization or irrigation [22]. These factors may have contributed to the model underestimation.
The results of the DFA showed that the variable SEV1 had the most influence in determining the severity of B. cinerea rot on mature bunches. SEV1 concerns B. cinerea infections occurring between "inflorescences clearly visible" and "berries groat-sized, bunches begin to hang". These early season infections may result in inflorescence and young berry rot, in latent infections, and in bunch trash colonisation [3]. We assessed bunch rot at bunch maturity but did not evaluate the accuracy of SEV1 in predicting rot of inflorescences and young clusters. Therefore, the relevance of SEV1 observed in this work may be caused by an indirect estimation of: i) latent infections of young berries that may have occurred at flowering; ii) the amount of bunch trash represented by aborted, infected flowers and young berries; and iii) the incidence of bunch trash colonized by B. cinerea. These three effects all increase Botrytis bunch rot at harvesting. That latent infections contribute to final disease incidence has been established [25,55,57,60,61], but the extent to which the onset of latent infections contributed to final disease was inconsistent and ranged from 5 to 100%; this inconsistency is probably explained by differences in locations, grape varieties, and the methods used to assess latency.
DFA also showed that SEV2 and SEV3, which account for the infections of ripening berries, had a minor effect on determining the severity of Botrytis epidemic. SEV2 and SEV3 relate to pathways Va and Vb (both were considered as pathway V by Elmer and Michaillides [3]), which have been previously considered to be relevant in several wine-growing regions [17,18,47,60,62,63,90,91]. SEV2 estimates the risk for infection of ripening berries by conidia, whereas SEV3 estimates the risk for berry-to-berry infection by aerial mycelium. The presence of colonized bunch trash and rotted berries within the bunch (sources of aerial conidia and mycelium) are essential for infection pathways Va and Vb, respectively. Thus, a high value of SEV2 and/or SEV3 (indicating favourable conditions for conidial and berry-to-berry infection, respectively) may result in no or low disease if there are no or few affected berries within the bunch. Because the presence of colonised bunch trash and the initial presence of infected berries during ripening are both indirectly accounted for by SEV1, the effect of SEV2 and SEV3 on bunch rot may also depend on SEV1. For this reason, an in-vineyard assessment of trash colonisation by B. cinerea during pre-bunch closure and/or of Botrytis rot incidence at veraison could greatly improve the accuracy of the model in predicting the infection severity of ripening berries.
DFA can also be used to predict the final bunch rot severity in an ongoing epidemic. Among the 21 epidemics considered in this work, prediction of the final disease severity changed from mild to intermediate before the end of flowering, and changed from intermediate to severe after veraison. Thus, it seems that SEV1 was the most important variable for discriminating between mild and intermediate epidemics, whereas SEV2 and SEV3 determined the change between intermediate and severe epidemics.
Currently, fungicide treatments are mainly applied based on vine phenology, with four vine growth stages being considered critical: end of flowering (A), pre-bunch closure (B), veraison (C), and before harvest (D) [33,92,93]. This method is easy to follow and provides excellent protection against Botrytis bunch rot in different grape-growing areas [42,[92][93][94][95][96][97][98]. Unfortunately, these treatments do not take into account the real risk of bunch infection at each growth stage, and they lead to expensive, unnecessary fungicide sprays [42,96]. To identify high-risk periods for B. cinerea infection, researchers have proposed the use of empirical rules [32,33], weather-driven models [15,20,47], and an expert system [43]. However, few or no validation data exist for some of the weather-driven models [15,20,47], validation results were inconsistent for others [34,[44][45][46], and no clear benefits have been demonstrated by the use of these methods [49].
The model proposed here: i) considers the key aspects of the B. cinerea life cycle; ii) integrates recent knowledge on host susceptibility at different growth stages and on the influence of weather on sporulation and on the main infection pathways; iii) is based on experiments [14,17,30] that considered the genetic variability of the B. cinerea populations in vineyards [12,99]; and iv) has been validated in different locations and under different environmental conditions. The model can then be regarded as an improvement of the previous Botrytis models in viticulture. Utility of the model in scheduling fungicides to control Botrytis bunch rot, however, should be verified with experiments [100][101][102]. After that, this new model could be used instead of the 15-15 rule or the previous models for making decisions about fungicide application at specific vine growth stages.
In the future, the model could be incorporated into existing decision support systems for sustainable vineyard management [38,103,104]. The model could also be improved with sitespecific information affecting the susceptibility to B. cinerea infection. This site-specific information could include the susceptibility of specific grape varieties [48,73,105] and the factors that predispose berries to infection, such as wounds, insect attacks, incidence of other diseases, or cultural practices [17,19,22,23,35,89].