Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

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

Abstract

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].

High genetic variation exists in B. cinerea populations, and four transposon genotypes (transposa, flipper-only, boty-only, and vacuma) have been distinguished [411]. The frequency of these transposon genotypes within a B. cinerea population depends on the growth stage of the grapevine, the organs infected [12], and the vineyard location [9,13].

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 [1517]. 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 [1517,19,20]. Moreover, disease incidence is higher on wounded than on unwounded berries [19,2123], 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) [2429]; (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,3234]. 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 [3941].

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,4346]. 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.

Model development

Life cycle of Botrytis cinerea

This section briefly describes the key elements of the B. cinerea life cycle used in model development, with emphasis on inoculum sources, production of conidia, and the main infection pathways described by Elmer and Michailides [3].

B. cinerea is an inhabitant of vineyards and grows and sporulates under a wide range of environmental conditions [35]. The fungus saprophytically colonizes various organic substrates including grape debris (trash) and weeds covering the soil [3,35,36], and infects both grape and alternative hosts [3].

B. cinerea mainly survives in grape debris (tendrils, mummified berries, dead leaves, and canes) as over-wintering mycelium [2528,50] and sclerotia [3,36,5153]. Abundant conidia are produced on these inoculum sources [3,28,53] under a range of environmental conditions [30]. Mature conidia become airborne with a circadian periodicity (spore concentrations in the air are highest at about midday), which is positively correlated with changes in temperature and wind velocity and negatively correlated with changes in relative humidity and the presence of dew [54]. Airborne conidia can easily travel from field to field [52]. Conidia have been considered the most important form of B. cinerea inoculum [3,35,47].

Conidia cause infection of inflorescences (at any growth stage) and young berries [3,14,35] through infection path I (conidial infection of the style and ovules), IIa (conidial infection of the stamens or petals), and IIb (fruit infection via the fruit pedicel) [3]. These infections cause either inflorescence blight or they develop into latent infections of berries [29,5559]. After veraison, latent infections become visible as rotted berries [56,58,60,61].

During the flowering stage, the pathogen also colonizes bunch trash (aborted flowers and fruitless calyptras and stamens), which are retained within the developing bunches (path III: conidial infection and extensive colonisation of floral debris in grape bunches; Elmer and Michailides [3]). B. cinerea over-summers as saprophytic mycelium on bunch trash, which has been considered a major source of inoculum within developing bunches [29]. Under favourable conditions, the mycelium-colonized bunch trash produces conidia, resulting in inoculation of the external surfaces of the ripening berries (path IV: conidial accumulation within the developing bunch; Elmer and Michailides [3]).

After veraison, a classical pre-harvest polycyclic epidemic can develop under favourable weather conditions; rot develops, and a new crop of conidia are dispersed to new infection sites (pathway V: conidial infection of ripening fruit; Elmer and Michailides [3]). In addition to conidial infection, ripening berries can be infected through contact with the aerial mycelium produced on adjacent infected berries (berry-to-berry infection) [17,18,43,62,63]. Thus, pathway Vb (berry-to-berry infection), which was not described by Elmer and Michailides [3], is included in the new model, and the infection of ripening berries by conidia (former pathway V) is renamed as pathway Va.

Because the life cycle of B. cinerea includes multiple sources of inoculum, conidia are generally present in the vineyard but at different densities [3,64,65].

Grape is susceptible to B. cinerea infection during two main periods: during the early developmental stages [66] and from veraison until ripening. In the latter period, susceptibility progressively increases, and this increase has been associated with structural and biochemical changes during berry maturation [18,22,67,68].

Modeling approach

The model was elaborated according to the principles of “systems analysis” [69]. The B. cinerea life cycle was divided into different state variables, and changes from one state to the following 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.

thumbnail
Fig 1. Relational diagram of the model simulating the life cycle of Botrytis cinerea.

Legend: boxes are state variables; line arrows show fluxes and direction of changes; valves define rates regulating these fluxes; diamonds show switches (i.e., conditions that open or close a flux); circles crossed by a line show parameters and external variables; dotted arrows show fluxes and direction of information from external variables to rates; circles are intermediate variables; clouds indicate outgoing variables. See Table 1 for acronym explanations.

https://doi.org/10.1371/journal.pone.0140444.g001

thumbnail
Table 1. List of variables, rates, and parameters used in the model describing the life cycle of Botrytis cinerea on grapevines.

https://doi.org/10.1371/journal.pone.0140444.t001

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: (1) with (2) (3)

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 aw≥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 berry-to-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: (4) where Teq = temperature equivalent as described for eq (2), with Tmin = 0°C and Tmax = 35°C;

WD = wetness duration (in hours); and (5)

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].

Eq (4) was developed by Ciliberti et al. [14]. Eq (5) was developed by using data from Ciliberti et al. [14] as described in S1 Appendix.

Relative infection severity in the first infection window is then calculated as: (6)

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: (7) where Teq = temperature equivalent as described for eq (2) with Tmin = 0°C and Tmax = 35°C, and (8)

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].

Eq (7) was developed by Ciliberti et al. [17], and eq (8) was developed based on data from Deytieux-belleau et al. [18] as explained in S1 Appendix.

Relative infection severity for conidial infection is then calculated as follows: (9)

Infection rate for berry-to-berry infection during the second infection window is calculated as follows: (10) where Teq = temperature equivalent as described in eq (2) with Tmin = 0, and Tmax = 30°C, and (11)

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: (12)

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 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 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).

thumbnail
Table 2. Summary characteristics (location, year, cultivar, and training system) of the vineyards used for validation of the Botrytis cinerea model. Also indicated are the observed incidence and severity of Botrytis bunch at maturity and the classification of the epidemics.

https://doi.org/10.1371/journal.pone.0140444.t002

thumbnail
Fig 2. Weather data and model output at Conselice in 2010 (CO-10).

A: daily data of temperature (red line, T in °C), relative humidity (green line, RH in %), rain (blue bars, in mm), and wetness duration (light blue area, WD in hours); B: predicted rate of spore production (SPOR) by Botrytis cinerea; C: predicted relative infection severity (RIS) on inflorescences and young clusters by conidia (RIS1; green bars) during the first infection period (green area), and on ripening berries by conidia (RIS2; blue bars) and mycelium (berry-to-berry infection) (RIS3; red bars) during the second infection period (orange area). Lines indicate the accumulated values of RIS1 (SEV1; green line) and of RIS2+RIS3 (SEV2+SEV3; purple line). Numbers in C indicate vine growth stages of full flowering (stage 65 of the scale of Lorenz et al. [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).

https://doi.org/10.1371/journal.pone.0140444.g002

thumbnail
Fig 3. Weather data and model output at S. Michelle in 2012 (SM-12).

A: daily data of temperature (red line, T in °C), relative humidity (green line, RH in %), rain (blue bars, in mm), and wetness duration (light blue area, WD in hours); B: predicted rate of spore production (SPOR) by Botrytis cinerea; C: predicted relative infection severity (RIS) on inflorescences and young clusters by conidia (RIS1; green bars) during the first infection period (green area), and on ripening berries by conidia (RIS2; blue bars) and mycelium (berry-to-berry infection) (RIS3; red bars) during the second infection period (orange area). Lines indicate the accumulated value of RIS1 (SEV1; green line) and of RIS2+RIS3 (SEV2+SEV3; purple line). Numbers in C indicate vine growth stages of full flowering (stage 65 of the scale of Lorenz et al. [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).

https://doi.org/10.1371/journal.pone.0140444.g003

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 [7881]. The analysis created two linear combinations of discriminant functions in the form: Fn = an + bn,1·ln(SEV1+1) + bn,2·ln(SEV2+SEV3+1), where n = 1 or 2 for the two discriminant functions, respectively. Briefly, the procedure automatically selected a first function (F1) that was able to separate the groups as much as possible; a second function (F2) 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 F1F2 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 (F1 and F2) by using the non-standardized canonical coefficients calculated in the previous section. This allowed us to locate the cases in the F1F2 space and to calculate the standardized distance (Duj) 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 F1F2 space as vectors, and matrix operations were used to resolve the equations [81]. The Mahalanobis distance was calculated in the form: where xu is the vector of case u, σj is the vector of the centroid of population j, and 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|xu) was calculated as follows [81]:

Where f(xu|j) is the likelihood of xu = 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(xu|j) with a multivariate approach as described by Huberty and Olejnik [81]:

Where |∑j| is the determinant of the covariance matrix of population, and Duj 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 (Pmild), intermediate (Pint), or severe (Psev) 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].

Results of Model Validation

Of the 21 B. cinerea epidemics, four were mild (with average disease incidence = 0.9±0.54, and average severity = 1.2 ± 1.08), 10 were intermediate (with disease incidence and severity equal to 50.6 ± 4.7 and 6.6 ± 1.4, respectively), and seven were severe (with disease incidence and severity of 82.6 ± 2.2 and 19.3 ± 5.3, respectively) (Table 2).

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 intermediate epidemics (0.20 ± 0.06 and 0.37 ± 0.16, respectively). The highest values of SEV2 and SEV3 were calculated at C0-10, where Botrytis bunch rot incidence was 81.8%, and severity was 15.9% at harvest (Table 2).

thumbnail
Fig 4. Model output (i.e., SEV1, SEV2, and SEV3) for 21 Botrytis cinerea epidemics.

The 21 epidemics are coded as in 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.

https://doi.org/10.1371/journal.pone.0140444.g004

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).

thumbnail
Table 3. Grouping of 21 Botrytis cinerea epidemics based on disease assessment in vineyards and as predicted based on the discriminant function analysis (DFA) using model output as discriminant variables.

https://doi.org/10.1371/journal.pone.0140444.t003

thumbnail
Fig 5. Assignment of 21 Botrytis cinerea epidemics to three severity groups based on the discriminant function analysis.

The 21 epidemics are coded as in 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.

https://doi.org/10.1371/journal.pone.0140444.g005

In this separation of epidemics into three groups, F1 accounted for 82.2% of the variation and F2 accounted for the remaining 17.2% (Table 4); canonical correlation coefficients (Table 4) and the significance of Wilks’ lambda tests (Table 5) all indicated that F1 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 F1 and that SEV2+SEV3 had a major effect in distinguishing the groups in F2 (Table 6).

thumbnail
Table 4. Statistics of the discriminant function analysis used to classify 21 Botrytis cinerea epidemics based on model output.

https://doi.org/10.1371/journal.pone.0140444.t004

thumbnail
Table 5. Significance of the discriminant function analysis used to classify 21 Botrytis cinerea epidemics based on model output.

https://doi.org/10.1371/journal.pone.0140444.t005

thumbnail
Table 6. Coefficients for each discriminant variable in each canonical function (F1 or F2) obtained in the discriminant analysis used to classify 21 Botrytis cinerea epidemics based on model output.

https://doi.org/10.1371/journal.pone.0140444.t006

Fig 6 shows the daily changes in the assignment of the 21 B. cinerea epidemics to one of three categories (mild, intermediate, or severe) based on the DFA of the daily model outputs (i.e., daily values of SEV1, SEV2, and SEV3). Five of these epidemics (at Z-12, Z-13, SM-12, SM-13, and RAV-12) remaining mild all season. The other epidemics were assigned as 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.

thumbnail
Fig 6. Assignment of 21 Botrytis cinerea epidemics to mild, intermediate, and severe groups based on model output.

Discriminant function analysis (DFA) was calculated daily between “inflorescences clearly visible” (growth stage 53 of Lorenz et al. [31]) and “berries ripe for harvest” (stage 89). The colors of the 21 horizontal bars indicate the daily assignment to epidemic groups based on DFA. The 21 epidemics are also coded on the left of each horizontal bar based on the severity of Botrytis bunch rot as observed in the vineyard at harvest and as indicated in Table 2. For both codes and bars, green, orange, and red indicate mild, intermediate, and severe epidemics, respectively. For example, the epidemic at BAG-12 was severe at harvest based on observation and progressed from mild to intermediate to severe based on DFA. Numbers in bars indicate the critical growth stages for fungicide applications; end of flowering (stage 69); pre-bunch closure (stage 77); veraison (stage 83).

https://doi.org/10.1371/journal.pone.0140444.g006

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 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 [8486]. 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, 9298]. 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,4446], 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 [100102]. 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 site-specific 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].

Supporting Information

S1 Appendix. Development of SUS1, SUS2 and SUS3 equations.

https://doi.org/10.1371/journal.pone.0140444.s001

(DOCX)

Acknowledgments

The authors thank M. Fermaud (INRA, Bordeaux, France) and R. Bugiani (Servizio Fitosanitario Regionale, Emilia-Romagna, Italy) for providing data used to validate the model.

Author Contributions

Conceived and designed the experiments: TC NC VR. Performed the experiments: TC NC VR. Analyzed the data: EGD VR. Contributed reagents/materials/analysis tools: EGD TC NC VR. Wrote the paper: EGD TC NC VR.

References

  1. 1. Williamson B, Tudzynski B, Tudzynski P, van Kan J a L. Botrytis cinerea: the cause of grey mould disease. Mol Plant Pathol. 2007;8: 561–580. pmid:20507522
  2. 2. Elad Y, Williamson B, Tudzynski P, Delen N. Botrytis spp. and diseases they cause in agricultural systems-an introduction. In: Elad Y, Williamson B, Tudzynski P, Delen N, editors. Botrytis: Biology, Pathology and Control. Springer Netherlands; 2007. pp. 1–6.
  3. 3. Elmer PAG, Michailides TJ. Epidemiology of Botrytis cinerea in orchard and vine crops. In: Elad Y, Williamson B, Tudzynski P, Delen N, editors. Botrytis: Biology, Pathology and Control. Springer Netherlands; 2007. pp. 243–272.
  4. 4. Giraud T, Fortini D, Levis C, Lamarque C, Leroux P, Lobuglio K, et al. Two sibling species of the Botrytis cinerea complex, transposa and vacuma, are found in sympatry on numerous host plants. Phytopathology. 1999;89: 967–973. pmid:18944743
  5. 5. Fournier E, Fortini D, Leroux P, Giraud T, Brygoo Y, Levis C. Characterization of Bc-hch, the Botrytis cinerea homolog of the neurospora crassa het-c vegetative incompatibility locus, and its use as a population marker. Mycologia. 2003;95: 251–261. pmid:21156611
  6. 6. Walker A-S, Gautier AL, Confais J, Martinho D, Viaud M, Le P Cheur P, et al. Botrytis pseudocinerea, a new cryptic species causing gray mold in French vineyards in sympatry with Botrytis cinerea. Phytopathology. 2011;101: 1433–1445. pmid:21830954
  7. 7. Fekete É, Fekete E, Irinyi L, Karaffa L, Árnyasi M, Asadollahi M, et al. Genetic diversity of a Botrytis cinerea cryptic species complex in Hungary. Microbiol Res. 2012;167: 283–291. pmid:22130648
  8. 8. Albertini C, Thebaud G, Fournier E, Leroux P. Eburicol 14α-demethylase gene (CYP51) polymorphism and speciation in Botrytis cinerea. Mycol Res. 2002;106: 1171–1178.
  9. 9. Muñoz G, Hinrichsen P, Brygoo Y, Giraud T. Genetic characterization of Botrytis cinerea populations in Chile. Mycol Res. 2002;106: 594–601.
  10. 10. DM Angelini RM, Masiello M, Rotolo C, Faretra F, De Miccolis A, Pollastro S. Molecular characterisation and detection of resistance to succinate dehydrogenase inhibitor fungicides in Botryotinia fuckeliana (Botrytis cinerea). Pest Manag Sci. 2014;70: 1884–1893. pmid:24481672
  11. 11. De Miccolis A, Milicevic T, Natale P, Lepore A, De Guido M, Pollastro S, et al. Botryotinia fuckeliana isolates carrying different transposons show differential response to fungicides and localization on host plants. J Plant Pathol. 2004;110: 208–214.
  12. 12. Martinez F, Dubos B, Fermaud M. The role of saprotrophy and virulence in the population dynamics of Botrytis cinerea in vineyards. Phytopathology. 2005;95: 692–700. pmid:18943786
  13. 13. Váczy KZ, Sándor E, Karaffa L, Fekete E, Fekete E, Arnyasi M, et al. Sexual recombination in the Botrytis cinerea populations in Hungarian vineyards. Phytopathology. 2008;98: 1312–1319. pmid:19000006
  14. 14. Ciliberti N, Fermaud M, Languasco L, Rossi V. Influence of fungal strain, temperature, and wetness duration on infection of grapevine inflorescences and young berry clusters by Botrytis cinerea. Phytopathology. 2015;105: 325–333. pmid:25354016
  15. 15. Nair N, Allen R. Infection of grape flowers and berries by Botrytis cinerea as a function of time and temperature. Mycol Res. 1993;97: 1012–1014.
  16. 16. Latorre BA, Rioja M. Efecto de la temperatura y de la humedad relativa sobre la germinacion de conidias de Botrytis cinerea. Cienc e Investig Agrar. 2002;29: 67–72.
  17. 17. Ciliberti N, Fermaud M, Roudet J, Rossi V. Environmental conditions affect Botrytis cinerea infection of mature grape berries more than the strain or transposon genotype. Phytopathology. 2015;105: 1090–1096. pmid:26218433
  18. 18. Deytieux-belleau C, Geny L, Roudet J, Mayet V, Donèche B, Fermaud M. Grape berry skin features related to ontogenic resistance to Botrytis cinerea. Eur J Plant Pathol. 2009;125: 551–563.
  19. 19. Nair N, Emmet R, Parker PE. Some factors predisposing grape berries to infection by Botrytis cinerea. New Zeal J Exp Agric. 1988;16: 257–263.
  20. 20. Broome J, English J, Marois J, Latorre BA, Aviles J. Development of an infection model for botrytis bunch rot of grapes based on wetness duration and temperature. Phytopathology. 1995;85: 97–102.
  21. 21. Kretschmer M, Kassemeyer H-H, Hahn M. Age-dependent grey mould susceptibility and tissue-specific defence gene activation of grapevine berry skins after infection by Botrytis cinerea. J Phytopathol. 2007;155: 258–263.
  22. 22. Mundy DC, Beresford RM. Susceptibility of grapes to Botrytis cinerea in relation to berry nitrogen and sugar concentration. New Zeal Plant Prot. 2007;60: 123–127.
  23. 23. Nelson K. Factors influencing the infection of table grapes by Botrytis cinerea. Phytopathology. 1951;41: 319–326.
  24. 24. Thomas A, Kotzé J, Matthee F. Development of a technique for recovery of soilborne sclerotia of Botrytis cinerea. Phytopathology. 1983;73: 1374–1376.
  25. 25. Nair N, Guilbaud-Oulton S, Barchia I, Emmet R. Significance of carry over inoculum, flower infection and latency on the incidence of Botrytis cinerea in berries of grapevines at harvest in New South Wales. Aust J Exp Agric. 1995;35: 1177–1180.
  26. 26. Mundy DC, Agnew RH, Wood PN. Grape tendrils as an inoculum source of Botrytis cinerea in vineyards—a review. New Zeal Plant Prot. 2012;65: 218–227.
  27. 27. Mundy DC, Haycock SR, Mclachlan ARG, Wood PN, Raw V. Tendrils as a source of seasonal carryover of Botrytis cinerea in vineyards. New Zeal Plant Prot. 2012;65: 236–240.
  28. 28. Jaspers MV, Seyb M, Trought MCT, Balasubramaniam R. Overwintering grapevine debris as an important source of Botrytis cinerea inoculum. Plant Pathol. 2013;62: 130–138.
  29. 29. Nair N, Parker PE. Midseason bunch rot of grapes: an unusual disease phenomenon in the Hunter Valley. Australas Plant Pathol. 1985;34: 302–305.
  30. 30. Ciliberti N, Fermaud M, Roudet J, Languasco L, Rossi V. Environmental effects on the production of Botrytis cinerea conidia on different media, grape bunch trash, and mature berries. Aust. J. Grape Wine Res. 2015; Accepted.
  31. 31. Lorenz D, Eichhorn K, Bleiholder H, Klose R, Meier U, Weber E. Phänologische Entwicklungsstadien der Weinrebe (Vitis vinifera L. ssp. vinifera). Vitic Enol Sci. 1994;49: 66–70.
  32. 32. Bulit J, Lafon R, Guillier G. Périodes favorables a l’application de traitments pour lutter contre la pourriture grise de la vigne. Phytiatrie-Phytopharmacie. 1970;19: 159–165.
  33. 33. Baldacci E, Belli G, Fogliani G. Osservazioni sul ciclo vitale della Botrytis cinerea Pers. nella vite. Not Mal delle Piante. 1962;62: 29–43.
  34. 34. Molot B. Modelling of grapevine grey mould by the potential state of infection system. Prog Agric Vitic. 1987;104: 355–358.
  35. 35. Jarvis WR. Botryotinia and Botrytis species: taxonomy, physiology and pathogenicity. Monograph no. 15. Ontario, Canada: Canada Department of Agriculture; 1977.
  36. 36. Coley-Smith J. Sclerotia and other structures of survival. In: Coley-Smith J, Verhoeff K, Jarvis W, editors. The Biology of Botrytis. London, UK: Academic Press; 1980. pp. 85–114.
  37. 37. Alavanja MCR, Hoppin J, Kamel F. Health effects of chronic pesticide exposure: cancer and neurotoxicity. Annu Rev Public Health. 2004;25: 155–197. pmid:15015917
  38. 38. Rossi V, Caffi T, Salinari F. Helping farmers face the increasing complexity of decision-making for crop protection. Phytopathol Mediterr. 2012;51: 457–479.
  39. 39. De Miccolis RM, Rotolo C, Masiello M, Pollastro S, Ishii H, Faretra F. Genetic analysis and molecular characterisation of laboratory and field mutants of Botryotinia fuckeliana (Botrytis cinerea) resistant to QoI fungicides. Pest Manag Sci. 2012;68: 1231–1240. pmid:22488841
  40. 40. Fernández-ortuño D, Grabke A, Li X, Schnabel G. Independent emergence of resistance to seven chemical classes of fungicides in Botrytis cinerea. Phytopathology. 2014;105: 424–432.
  41. 41. Shtienberg D. Rational management of Botrytis-incited diseases: integration of control measures and use of warning systems. In: Elad Y, Williamson B, Tudzynski P, Delen N, editors. Botrytis: biology, pathology and control. Springer; 2007. pp. 335–346.
  42. 42. Lafon R, Verdu D, Bulit J. Mise au point sur le traitement de la pourriture grise dans le vignoble. Rev Zool Agric Pathol végétale. 1972;1: 31–43.
  43. 43. Strizyk S. Modéle d’état potential d’infection. Application au Botrytis cinerea de la vigne. Association de coordination technique agricole (ACTA). Paris, France; 1985.
  44. 44. Zerbetto F, Rho G, Vercesi A, Cortesi P, Bisiach M. Il modello previsionale EPI-Botrytis nella difesa di Riesling Italico. Not sulle Mal delle piante. 1988;109: 28–36.
  45. 45. Vercesi A, Zerbetto F, Cortesi P, Bisiach M. Applicabilità di EPI-Botrytis ad alcuni vitini a grappolo compatto. Not sulle Mal delle piante. 1990;111: 125–158.
  46. 46. Egger E, Marinelli M, D’Arcangelo M. Valutazione del modello EPI-Botrytis in provincia di Arezzo. Giornate Fitopatologiche. Montevisiano, Italy: ATTI; 1994. pp. 157–166.
  47. 47. Rodríguez-Rajo FJ, Jato V, Fernández-González M, Aira MJ. The use of aerobiological methods for forecasting Botrytis spore concentrations in a vineyard. Grana. 2010;49: 56–65.
  48. 48. Ellison P, Ash G, McDonald C. An expert system for the management of Botrytis cinerea in Australian vineyards. I. Development. Agric Syst. 1998;56: 185–207.
  49. 49. Ash G, Mcdonald C, Ellison P. An expert system for the management of Botrytis cinerea in Australian vineyards. II. Validation. Agric Syst. 1998;56: 209–224.
  50. 50. Seyb A, Gaunt R, Trought M, Frampon C, Balasubramaniam R, Jaspers M V. Relationship between debris within grape bunches and Botrytis infection of berries. Proceedings of the New Zealand plant protection conference. 2000. p. 451.
  51. 51. Thomas A, Matthee F, Kotzé J. Survival of Botrytis cinerea from table grapevines in South Africa. Phytophylactica. 1981;13: 157–160.
  52. 52. Holz G, Coertze S, Williamson B. The Ecology of Botrytis on Plant Surfaces. In: Elad Y, Williamson B, Tudzynski P, Delen N, editors. Botrytis: Biology, Pathology and Control. Dordrecht: Springer Netherlands; 2007. pp. 9–24.
  53. 53. Nair N, Nadtotchei A. Sclerotia of Botrytis as a source of primary inoculum for bunch rot of grapes in New South Wales, Australia. J Phytopathol. 1987;119: 42–51.
  54. 54. Chastagner G, Ogawa JM, Manji B. Dispersal of conidia of Botrytis cinerea in tomato fields. Phytopathology. 1978;68: 1172–1176.
  55. 55. McClellan W, Hewitt W. Early botrytis rot of grapes: time of infection and latency of Botrytis cinerea Pers. in Vitis vinifera L. Phytopathology. 1973;63: 1151–1157.
  56. 56. Pezet R, Pont V. Infection florale et latence de Botrytis cinerea dans les grappes de Vitis vinifera (var. Gamay). Reveu suisse Vitic Arboric Hortic. 1986;18: 317–322.
  57. 57. Holz G, Gütschow M, Coertze S, Calitz FJ. Occurrence of Botrytis cinerea and subsequent disease expression at different positions on leaves and bunches of grape. Plant Dis. 2003;87: 351–358.
  58. 58. Viret O, Keller M, Jaudzems VG, Cole FM. Botrytis cinerea infection of grape flowers: light and electron microscopical studies of infection sites. Phytopathology. 2004;94: 850–7. pmid:18943105
  59. 59. Holz G, Coertze S, Basson E. Latent infection of Botrytis cinerea in grape pedicels leads to post harvest decay. Phytopathology. 1997;87: S43.
  60. 60. Keller M, Viret O, Cole FM. Botrytis cinerea infection in grape flowers: defense reaction, latency, and disease expression. Phytopathology. 2003;93: 316–322. pmid:18944341
  61. 61. Pezet R, Viret O, Perret C, Tabacchi R. Latency of Botrytis cinerea Pers.: Fr. and biochemical studies during growth and ripening of two grape berry cultivars, respectively susceptible and resistant to grey mould. J Phytopathol. 2003;151: 208–214.
  62. 62. Thomas C, Marois J, English JT. The effects of wind speed, temperature, and relative humidity on development of aerial mycelium and conidia of Botrytis cinerea on grape. Phytopathology. 1988;78: 260–265.
  63. 63. Bélanger M-C, Roger J-M, Cartolaro P, Fermaud M. Autofluorescence of grape berries following Botrytis cinerea infection. Int J Remote Sens. 2011;32: 3835–3849.
  64. 64. Vercesi A, Bisiach M. Indagine sulla fluttuazione del potenziale d’inoculo di Botrytis cinerea pers. in vigneto. Riv di Patol Veg. 1982;18: 13–48.
  65. 65. Yunis H, Elad Y. Survival of dicarboximide-resistant strains of Botrytis cinerea in plant debris during summer in Israel. Phytoparasitica. 1989;17: 13–21.
  66. 66. Hill G, Stelwaag-Kittler F, Huth G, Schlosser E. Resistance of grapes in different developmental stages to Botrytis cinerea. Phytopathol Zettschrift. 1981;102: 328–338.
  67. 67. Padgett M, Morrison JC. Changes in grape berry exudates during fruit development and their effect on mycelial growth of Botrytis cinerea. J Am Soc Hortic Sci. 1990;115: 269–273.
  68. 68. Valdés-Gómez H, Fermaud M, Roudet J, Calonnec A, Gary C. Grey mould incidence is reduced on grapevines with lower vegetative and reproductive growth. Crop Prot. 2008;27: 1174–1186.
  69. 69. Leffelaar P, Ferrari T. Some elements of dynamic simulation. In: Rabbinge R, Ward S, van Laar H, editors. Simulation and system management in crop protection. Wageningen, The Netherlands: Pudoc Scientific Publishers; 1989. pp. 19–45.
  70. 70. Rossi V, Caffi T. Effect of water on germination of Plasmopara viticola oospores. Plant Pathol. 2007;56: 957–966.
  71. 71. Sosa-Alvarez M, Madden LV, Ellis M. Effects of temperature and wetness duration on sporulation of Botrytis cinerea on strawberry leaf residues. Plant Dis. 1995;79: 609–615.
  72. 72. Kottek M, Grieser J, Beck C, Rudolf B, Rubel F. World Map of the Köppen-Geiger climate classification updated. Meteorol Zeitschrift. 2006;15: 259–263.
  73. 73. Vail ME, Marois JJ. Grape cluster architecture and the susceptibility of berries to Botrytis cinerea. Phytopathology. 1991;81: 188–191.
  74. 74. Raw V, Mundy DC, McLachlan ARG, Clifford C, and Walter M. Botrytis cinerea control on Sauvignon blanc using a recycling sprayer. N. Z. Plant Prot. 2012; 65: 249–255.
  75. 75. Molitor D, Rothmeier M, Behr M, Fischer S, Hoffmann L and Evers D. Crop cultural and chemical methods to control grey mould on grapes. Vitis. 2011; 50: 81–87.
  76. 76. Jacometti MA,Wratten SD and Walter M. Review: Alternatives to synthetic fungicides for Botrytis cinerea. Aust. J. Grape Wine Res. 2010; 16: 154–172.
  77. 77. Gent DH, De Wolf E, Pethybridge SJ. Perceptions of risk, risk aversion, and barriers to adoption of decision support systems and integrated pest management: an introduction. Phytopathology. 2011;101: 640–643. pmid:21117876
  78. 78. Hau B, Kranz J. Mathematics and statistics for analysis in epidemiology. In: Kranz J, editor. Epidemics of plant diseases. New York, USA: Springer; 1990. pp. 12–52.
  79. 79. Mendes A, Soares E, Azevedo J. Efficiency measures in the agricultural sector: with applications. New York, USA: Springer; 2013.
  80. 80. Quinn G, Keough M. Experimental design and data analysis for biologists. Cambridge, MA, USA: MIT Press, 2011.
  81. 81. Huberty C, Olejnik S. Applied manova and discriminant analysis. Wiley 2006.
  82. 82. Sanogo S, Yang XB. New Applications of selected multivariate statistical methods and their use in phytopathological research. Phytopathology. 2004;94: 1004–1006.
  83. 83. Rossi V, Giosuè S, Caffi T. Modelling plant diseases for decision making in crop protection. In: Oerke E-C, Gerhards R, Menz G, Sikora RA, editors. Precision Crop Protection—the Challenge and Use of Heterogeneity. Dordrecht: Springer Netherlands; 2010. pp. 241–258.
  84. 84. Workneh F, Allen TW, Nash GH, Narasimhan B, Srinivasan R, Rush CM. Rainfall and temperature distinguish between Karnal bunt positive and negative years in wheat fields in Texas. Phytopathology. 2008;98: 95–100. pmid:18943243
  85. 85. Salinari F, Giosue S, Tubiello F, Rettori A, Rossi V, Spanna F, et al. Downy mildew (Plasmopara viticola) epidemics on grapevine under climate change. Glob Chang Biol. 2006;12: 1299–1307.
  86. 86. Jiménez-Díaz RM, Olivares-García C, Landa BB, Jiménez-Gasco MM, Navas-Cortés JA. Region-Wide analysis of genetic diversity in Verticillium dahliae populations infecting olive in southern Spain and agricultural factors influencing the distribution and prevalence of vegetative compatibility groups and pathotypes. Phytopathology. 2010;101: 304–315.
  87. 87. Fermaud M, Le Menn R. Association of Botrytis cinerea with grape berry moth larvae. Phytopathology. 1989;79: 651–656.
  88. 88. Fermaud M, Le Menn R. Transmission of Botrytis cinerea to grapes by grape berry moth larvae. Phytopathology. 1992;82: 1393–1398.
  89. 89. Fermaud M, Gaunt R, Elmer P. The influence of Thrips obscuratus on infection and contamination of kiwifruit by Botrytis cinerea. Plant Pathol. 1994;43: 953–960.
  90. 90. Latorre B, Lillo C, Rioja M. Effect of timing on the efficacy of fungicide treatments applied against Botrytis cinerea of grapevine. Cienc e Investig Agrar. 2001;28: 61–66.
  91. 91. Wilcox WF. Controlling Botrytis: a perspective from the eastern USA. New Zeal Grapegrow Winemak. 2002;10: 22–27.
  92. 92. Agulhon R. La pourriture grise des raisins. Le maladies des plantes. Paris, France: ACTA; 1969.
  93. 93. Agulhon R. Quelques aspects de la lutte contre la pourriture grise en 1972. Vigne vins. 1973;220: 5–11.
  94. 94. Bisiach M. Nuove acquisizioni nella protezione antibotritica in viticoltura. Giornate Fitopatologiche. Torino, Italy: ATTI; 1975. pp. 753–769.
  95. 95. Bisiach M, Minervini G, Zerbetto F. L’epidemia di Botrytis cinerea in due località dell’Italia settentrionale e la possibilità di interventi fitoiatrici. Inf Agrar. 1978;24.
  96. 96. Bisiach M, Zerbetto F, Minervini G. Possibilità di lotta secondo il metodo fenologico e climatico. Inf Agrar. 1978;24: 2287–2290.
  97. 97. Gartel W. Les problèmes posés par la pourriture grise sur la vigne, comment les résoudre. Vigne vins. 1977;260: 15.
  98. 98. Bolay A, Rochaix M. Heurs et malheurs des fongicides systémiques dans la lutte contre la paurriture grise des raisins. Reveu suisse Vitic Arboric Hortic. 1975;7: 137–147.
  99. 99. Martinez F, Blancard D, Lecomte P, Levis C, Dubos B, Fermaud M. Phenotypic differences between vacuma and transposa subpopulations of Botrytis cinerea. Eur J Plant Pathol. 2003;109: 479–488.
  100. 100. Caffi T, Rossi V, Bugiani R. Evaluation of a warning system for controlling primary infections of grapevine downy mildew. Plant Dis. 2010;94: 709–716.
  101. 101. Caffi T, Legler SE, Rossi V, Bugiani R. Evaluation of a warning system for early-season control of grapevine powdery mildew. Plant Dis. 2012;96: 104–110.
  102. 102. Gent DH, Mahaffee WF, McRoberts N, Pfender WF. The use and role of predictive systems in disease management. Annu Rev Phytopathol. 2013;51: 267–289. pmid:23682914
  103. 103. Rossi V, Salinari F, Poni S, Caffi T, Bettati T. Addressing the implementation problem in agricultural decision support systems: the example of vite.net®. Comput Electron Agric.; 2014;100: 88–99.
  104. 104. Shtienberg D. Will decision-support systems be widely used for the management of plant diseases? Annu Rev Phytopathol. 2013;51: 1–16. pmid:23767845
  105. 105. Gabler FM, Smilanick JL, Mansour M, Ramming DW, Mackey BE. Correlations of morphological, anatomical, and chemical features of grape berries with resistance to Botrytis cinerea. Phytopathology. 2003;93: 1263–1273. pmid:18944326