# Assessing the durability and efficiency of landscape-based strategies to deploy plant resistance to pathogens

Loup Rimbaud, Julien Papaix, Jean-François Rey, Luke G. Barrett, Peter H. Thrall

<figure><img src="image_1.png" /></figure>

## To cite this version:

Loup Rimbaud, Julien Papaix, Jean-François Rey, Luke G. Barrett, Peter H. Thrall. Assessing the durability and efficiency of landscape-based strategies to deploy plant resistance to pathogens. PLoS Computational Biology, 2018, 14 (4), pp.e1006067. 10.1371/journal.pcbi.1006067. hal-02629098

HAL Id: hal-02629098

https://hal.inrae.fr/hal-02629098

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Assessing the durability and efficiency of landscape-based strategies to deploy plant resistance to pathogens

Loup Rimbaud¹*, Julien Papaïx²*, Jean-François Rey², Luke G. Barrett¹, Peter H. Thrall¹

1 CSIRO Agriculture and Food, Canberra, ACT, Australia, 2 BioSP, INRA, Avignon, France

© These authors contributed equally to this work.

* loup.rimbaud06@aptalumni.org

<figure><img src="image_4.png" /></figure>

## Abstract

Genetically-controlled plant resistance can reduce the damage caused by pathogens. However, pathogens have the ability to evolve and overcome such resistance. This often occurs quickly after resistance is deployed, resulting in significant crop losses and a continuing need to develop new resistant cultivars. To tackle this issue, several strategies have been proposed to constrain the evolution of pathogen populations and thus increase genetic resistance durability. These strategies mainly rely on varying different combinations of resistance sources across time (crop rotations) and space. The spatial scale of deployment can vary from multiple resistance sources occurring in a single cultivar (pyramiding), in different cultivars within the same field (cultivar mixtures) or in different fields (mosaics). However, experimental comparison of the efficiency (i.e. ability to reduce disease impact) and durability (i.e. ability to limit pathogen evolution and delay resistance breakdown) of landscape-scale deployment strategies presents major logistical challenges. Therefore, we developed a spatially explicit stochastic model able to assess the epidemiological and evolutionary outcomes of the four major deployment options described above, including both qualitative resistance (i.e. major genes) and quantitative resistance traits against several components of pathogen aggressiveness: infection rate, latent period duration, propagule production rate, and infectious period duration. This model, implemented in the R package *landsepi*, provides a new and useful tool to assess the performance of a wide range of deployment options, and helps investigate the effect of landscape, epidemiological and evolutionary parameters. This article describes the model and its parameterisation for rust diseases of cereal crops, caused by fungi of the genus *Puccinia*. To illustrate the model, we use it to assess the epidemiological and evolutionary potential of the combination of a major gene and different traits of quantitative resistance. The comparison of the four major deployment strategies described above will be the objective of future studies.

<figure><img src="image_5.png" /></figure>

OPEN ACCESS

**Citation:** Rimbaud L, Papaïx J, Rey J-F, Barrett LG, Thrall PH (2018) Assessing the durability and efficiency of landscape-based strategies to deploy plant resistance to pathogens. PLoS Comput Biol 14(4): e1006067. https://doi.org/10.1371/journal.pcbi.1006067

**Editor:** Samuel Alizon, CNRS, FRANCE

**Received:** September 4, 2017

**Accepted:** February 27, 2018

**Published:** April 12, 2018

**Copyright:** © 2018 Rimbaud et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files. The code of the model is implemented in the R package *landsepi*. Landscape Epidemiology and Evolution (version 0.0.2, https://cran.r-project.org/package=landsepi).

**Funding:** LR, LGB, and PHT were funded by grant CSP00192 from the Grain Research & Development Corporation (GRDC, https://grdc.com.au/). All authors were supported by the INRA/CSIRO linkage program. The funders had no role in
study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

## Author summary

There are many recent examples which demonstrate the evolutionary potential of plant pathogens to overcome the resistances deployed in agricultural landscapes to protect our crops. Increasingly, it is recognised that how resistance is deployed spatially and temporally can impact on rates of pathogen evolution and resistance breakdown. Such deployment strategies are mainly based on the combination of several sources of resistance at different spatiotemporal scales. However, comparison of these strategies in a predictive sense is not an easy task, owing to the logistical difficulties associated with experiments involving the spread of a pathogen at large spatio-temporal scales. Moreover, both the durability of a strategy and the epidemiological protection it provides to crops must be assessed since these evaluation criteria are not necessarily correlated. Surprisingly, no current simulation model allows a thorough comparison of the different options. Here we describe a spatio-temporal model able to simulate a wide range of deployment strategies and resistance sources. This model, implemented in the R package *landsepi*, facilitates assessment of both epidemiological and evolutionary outcomes across simulated scenarios. In this work, the model is used to investigate the combination of different sources of resistance against fungal diseases such as rusts of cereal crops.

## Introduction

The deployment of resistant cultivars in agricultural landscapes aims to reduce the ability of plant pathogens to cause disease on crops. However, the durability of plant resistance has often been limited by evolutionary changes in pathogen populations [1]. Typically, there are two main types of resistance. Although exceptions exist, qualitative (or 'major gene') resistance is usually monogenic and complete, i.e. conferring full immunity [2-4]. In contrast, quantitative resistance is mostly polygenic and partial, i.e. infection is still possible but pathogen development is reduced to a greater or lesser extent. Consequently, quantitative resistance is often described as affecting one or more components of pathogen aggressiveness (defined as the quantitative ability to colonise and cause damage to the host): lower rate of infection, longer latent period, reduced propagule production, shorter infectious period or lower toxin production [5-8].

Regardless of the source of resistance, pathogens have the potential to evolve infectivity (defined as the qualitative ability to infect the host) and aggressiveness in response to the selection posed by plant resistance [4, 9]. This adaptation of pathogens at the population level is generally driven by mutations from non-adapted pathogen genotypes, changes in frequencies of pre-existing adapted genotypes, or introductions via migration from distant areas [9-11]. Emergence of novel pathotypes can result in the breakdown of qualitative resistance or the erosion of quantitative resistance, and consequently restoration of the infectivity or aggressiveness of pathogen population. The economic costs induced by pathogen adaptation may be huge, due firstly to the yield losses directly associated with an epidemic, and secondly because of the significant research and breeding efforts required to identify new resistance sources and develop new resistant cultivars [12]. Several deployment strategies have been proposed to improve the cost-effectiveness of plant resistance and prevent the frequently documented breakdown of major genes after their uniform deployment over large areas [3, 4, 9, 13, 14]. These strategies rely on the use of quantitative resistance, which is thought to be more durable because it poses smaller selection pressure on pathogen populations [1], or on the management of host genetic diversity [15-17]. This diversity can be introduced in time through crop
rotations (e.g. recurring succession of different crops in the same field [18]). In space, different crops can be combined in the same field in cultivar mixtures [19, 20] or in different fields of the landscape as mosaics [3, 17]. Finally, several resistance sources can be stacked in the same cultivar through pyramiding [21, 22].

Few empirical studies have directly compared the performance of these different strategies, probably owing to the difficulty of implementing landscape scale experiments in practice (but see [23] for a comparison of all these strategies, except mosaics, using plastic tunnels). Models are not constrained by this difficulty, but surprisingly there are currently no published models enabling a global comparison. Indeed, most models focus on a single strategy (e.g. crop rotation [24, 25]; mixtures [26–28]; pyramiding [29, 30]; mosaics [31–42]) or a combination of strategies (e.g. mixture and pyramiding: [43]; mosaic and pyramiding: [44, 45]). Only a few studies explicitly compare two types of strategies [46–50] and only two studies evaluated more than two strategies [51, 52]. As a result, a comprehensive evaluation of different deployment schemes is complicated, and currently only feasible via pairwise comparisons [53, 54]. The situation for quantitative resistance is similar, since often only one [28, 34, 41, 42, 55], two [36, 37, 56], or a combination [26, 44, 49] of pathogen aggressiveness components are targeted, although quantitative resistance can affect several life-history traits of the pathogen. As articulated above, this current gap in our ability to predict which strategy will maximise our ability to control disease epidemics as well as pathogen evolutionary potential (or indeed whether these goals are compatible) emphasises the need for models that can compare different deployment schemes within the same framework, using standardised assumptions.

Comparison of different resistance deployment strategies requires the use of relevant evaluation criteria, which may vary depending on stakeholder objectives, and thus have an impact on the optimal strategy [39, 45]. Most published models focused only on one criterion, like resistance durability (i.e. the duration from initial deployment to the time when resistance is considered to have been overcome) (e.g. [30, 37, 47, 51, 52]), or epidemiological protection (i.e. reduction in pathogen population size and consequently in the proportion of diseased plants) (e.g. [26, 28, 32, 34–36, 43, 49, 50]). However, no correlation seems to exist between durability and epidemiological protection [45, 57], and those objectives can even be incompatible in severe epidemic contexts [33]. It is therefore essential to develop methods to assess deployment strategies against multiple evolutionary and epidemiological criteria [3, 17, 58].

The present study describes a demo-genetic model (i.e. it includes pathogen population demographic dynamics and its genetic evolution) which simulates the spread of a pathogen in a spatially explicit agricultural landscape. The model is flexible and can simulate mosaics, mixtures, rotation, and pyramiding of different major genes and up to four traits for quantitative resistance, acting on different components of pathogen aggressiveness (infection rate, latent period, infectious period, reproduction rate). Performance of resistance deployment is evaluated using several criteria from evolutionary (e.g. durability) or epidemiological (e.g. disease severity) perspectives. Thus, the model enables direct comparison of a range of spatio-temporal deployment strategies, but also enables investigation of the effects of landscape, epidemiological and evolutionary parameters on the ability of a given strategy to control disease. Although the main purpose of this paper is to provide a comprehensive description of the simulation model, we take advantage of the generality of the model and address three specific questions of interest to the scientific community. These questions aim to assess the potential of the combination of qualitative and quantitative resistance, given that the former prevents disease spread but is prone to breakdown, whereas the latter allows some disease development but to a lesser extent:

1. What is the durability of a major gene alone or combined with another source of resistance?
2. What is the level and speed of erosion of quantitative resistance alone or combined with a major gene?

3. What is the severity of epidemics in the landscape when qualitative and quantitative resistances are deployed alone or in combination?

The model can be parameterised to simulate various pathogen life histories. Here, we investigate the above questions in the context of rust diseases of cereal crops, caused by fungi of the genus *Puccinia* which can dramatically affect crop yields worldwide [59, 60]. Over the past decades, breeders have been engaged in an arms race against these pathogens, which exhibit high evolutionary potential in terms of their ability to overcome resistant crop cultivars following deployment [45, 60–63].

# Models

## Overview of the model

The model is stochastic, spatially explicit (the basic spatial unit is an individual field), based on a SEIR ('susceptible-exposed-infectious-removed') structure with a discrete time step. It simulates the spread (through clonal reproduction and dispersal) and evolution (via mutation) of a pathogen in an agricultural landscape, across cropping seasons split by host harvests which represent potential bottlenecks to the pathogen. The model is based on the model described in previous articles [40, 45], but has been considerably modified and extended to enable simulation of a wide array of deployment strategies: mosaics, mixtures, rotations and pyramiding of multiple major resistance genes which affect pathogen infectivity, and up to four quantitative resistance traits. These traits target different aggressiveness components of the pathogen, i.e. the infection rate, the duration of the latent period and the infectious period, and the propagule production rate. Initially, the pathogen is not adapted to any source of resistance and is only present on susceptible hosts. However, through mutation, it can evolve and may acquire infectivity genes (which leads to breakdown of major resistance genes) or increase aggressiveness (which leads to the erosion of the relevant quantitative resistance traits). Evolution of a pathogen toward infectivity or increased aggressiveness on a resistant host is often penalised by a fitness cost on susceptible hosts [2, 64–66]. Consequently, in the present model, pathogens carrying infectivity genes may have lower infection rates (cost of infectivity) on susceptible hosts relative to pathogens that do not carry these genes. Similarly, a gain in pathogen aggressiveness on quantitatively resistant hosts is penalised by decreased aggressiveness on susceptible hosts, leading to a trade-off.

The evolutionary outcome of a deployment strategy is assessed by measuring the time until the pathogen achieves the three steps to adapt to plant resistance: (d1) first appearance of adapted mutants, (d2) initial migration to resistant hosts and infection (also referred as 'arrival' or 'introduction' in invasion biology), and (d3) broader establishment in the resistant host population (i.e. the point at which extinction becomes unlikely). Epidemiological outcomes are evaluated using the Green Leaf Area (GLA) as a proxy for yield, and the area under the disease progress curve (AUDPC) to measure disease severity.

## Landscape and resistance deployment strategies

In this study, a cropping landscape is represented by both its physical structure (defined as the spatial arrangement of field boundaries) and its genetic composition (defined by the allocation of crop cultivars within and among individual fields).

**Landscape structure.** Both real and simulated landscape structures can be used as input to the model. In this study, the landscape structure is simulated using a T-tessellation
<figure><img src="image_6.png" /><figcaption>**Fig 1. Allocation of two cultivars in the landscape.** Firstly, a landscape structure is generated using T-tessellations. Secondly, the susceptible cultivar (SC, white) and resistant cultivar (RC, grey) are allocated to fields. Both the proportions of the surface coverage (horizontal axis: 50% in A and B, 80% in C and D) and level of spatial aggregation (vertical axis: high in A and C, low in B and D) for each cultivar are controlled via model input parameters.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g001

algorithm [67] (see S1 Fig) in order to control specific features such as number, area and shape of the fields, as described in previous studies [40, 45].

**Allocation of crop cultivars.** We used an algorithm based on latent Gaussian fields to allocate two different crop cultivars across the simulated landscapes (e.g. a susceptible and a resistant cultivar, denoted as SC and RC, respectively). This algorithm allows control of the proportions of each cultivar in terms of surface coverage, and their level of spatial aggregation (Fig 1). Briefly, a random vector of values is drawn from a multivariate normal distribution with expectation 0 and a variance-covariance matrix which depends on the pairwise distances between the centroids of the fields. Next, the crop cultivars are allocated to different fields
depending on whether each value drawn from the multivariate normal distribution is above or
below a threshold. The proportion of each cultivar in the landscape is controlled by the value
of this threshold. The sequential use of this algorithm allows the allocation of more than two
crop cultivars (see S2 Fig). Therefore, deployment strategies involving two sources of resis-
tance can be simulated by: (1) running the allocation algorithm once to segregate the fields
where the susceptible cultivar is grown, and (2) applying one of the following deployment
strategies to the remaining candidate fields:

1. Mosaics: two resistant cultivars (RC₁ and RC₂, carrying the first and the second resistance sources, respectively) are assigned to candidate fields by re-running the allocation algorithm;

2. Mixtures: both RC₁ and RC₂ are allocated to all candidate fields;

3. Rotations: RC₁ and RC₂ are alternatively cultivated in candidate fields, depending on the number of cropping seasons over which a given cultivar is grown before being rotated;

4. Pyramiding: all candidate fields are populated with RC₁₂, a resistant cultivar carrying both resistance sources.

Using this approach, a wide range of strategies can be simulated, noting that RC₁ and RC₂
are not necessarily deployed in balanced proportions or at a specific level of aggregation (i.e.
these can also be varied to explore their effects). The same approach can be used to deploy
more than two resistance sources. Alternatively, any crop allocation model (e.g. [68]) could be
used to provide a more realistic landscape composition.

**Host-pathogen genetic interaction**

All model parameters and values used in the simulations are listed in Table 1 (see S1 Text for details on model parameterisation to rust pathogens).

**Host genotype.** A host genotype (indexed by v) is represented by a set of binary variables, indicating the resistances it carries, denoted as $mg_g$ (g = 1, . . . ,G) for major genes and $qr_w$ (w = e,γ,r,Y, see also below) for quantitative resistance traits against the main components of pathogen aggressiveness.

**Pathogen genotype.** A pathogen genotype (indexed by p) is represented by a set of binary variables indicating whether it carries infectivity genes (variables $ig_g$) and four discrete variables giving the level of aggressiveness for each life-history trait (variables $ag_w$). All of these infectivity genes and aggressiveness components are assumed to evolve independently from each other.

**Infectivity matrix.** The interaction between potential major host resistance genes and associated pathogen infectivity genes is represented by a multiplicative factor of the infection rate of the pathogen. Thus, for major gene g, the possible interactions are summarised in an infectivity matrix (denoted as $INF^g$, Table 2).

In the infectivity matrix, $\rho_g$ is the efficiency of major resistance gene g on a non-adapted pathogen ($\rho_g = 1$ for a typical qualitative resistance gene conferring immunity [69], but $\rho_g < 1$ if a major gene is only partially expressed [3]), and $\theta_g$ is the cost of infectivity paid by an adapted pathogen on a susceptible host ($\theta_g = 0$ means absence of cost of infectivity, $\theta_g = 1$ means loss of infectivity on the susceptible cultivar).

**Aggressiveness matrix.** As for the infectivity matrices, the four potential traits for quantitative resistance (in the host) and the associated components for aggressiveness (in the pathogen) lead to four similar aggressiveness matrices (denoted as $AGG^w$ for component w,
<table><caption>Table 1. Summary of the parameters used in the model and values for rust pathogens.</caption><thead><tr><th>Notation</th><th>Parameter</th><th>Value(s) for rust pathogens used in the simulations <sup>a</sup></th></tr></thead><tbody><tr><td colspan="3"><strong>Simulation parameters</strong></td></tr><tr><td>Y</td><td>Number of simulated years</td><td>50 years</td></tr><tr><td>T</td><td>Number of time steps in a cropping season</td><td>120 days.year<sup>-1</sup></td></tr><tr><td>J</td><td>Number of fields in the landscape</td><td>{155; 154; 152; 153; 156}<sup>b</sup></td></tr><tr><td>V</td><td>Number of host cultivars</td><td>2</td></tr><tr><td colspan="3"><strong>Initial conditions and seasonality</strong></td></tr><tr><td>C<sub>v</sub><sup>0</sup></td><td>Plantation host density of cultivar v</td><td>0.1 m<sup>-2</sup> <sup>c</sup></td></tr><tr><td>C<sub>v</sub><sup>max</sup></td><td>Maximal host density of cultivar v</td><td>2 m<sup>-2</sup> <sup>c</sup></td></tr><tr><td>δ<sub>v</sub></td><td>Host growth rate of cultivar v</td><td>0.1 day<sup>-1</sup> <sup>c</sup></td></tr><tr><td>φ</td><td>Initial probability of infection</td><td>5.10<sup>-4</sup></td></tr><tr><td>λ</td><td>Off-season survival probability</td><td>10<sup>-4</sup></td></tr><tr><td colspan="3"><strong>Pathogen aggressiveness components</strong></td></tr><tr><td>e<sub>max</sub></td><td>Maximal expected infection rate</td><td>0.40 spore<sup>-1</sup></td></tr><tr><td>γ<sub>min</sub></td><td>Minimal expected latent period duration</td><td>10 days</td></tr><tr><td>γ<sub>var</sub></td><td>Variance of the latent period duration</td><td>9 days</td></tr><tr><td>Y<sub>max</sub></td><td>Maximal expected infectious period duration</td><td>24 days</td></tr><tr><td>Y<sub>var</sub></td><td>Variance of the infectious period duration</td><td>105 days</td></tr><tr><td>r<sub>max</sub></td><td>Maximal expected propagule production rate</td><td>3.125 spores.day<sup>-1</sup></td></tr><tr><td colspan="3"><strong>Pathogen dispersal</strong></td></tr><tr><td>g(.)</td><td>Dispersal kernel</td><td>Power-law function <sup>d</sup></td></tr><tr><td>a</td><td>Scale parameter</td><td>40</td></tr><tr><td>b</td><td>Width of the tail</td><td>7</td></tr><tr><td>π(.)</td><td>Contamination function</td><td>Sigmoid curve</td></tr><tr><td>κ</td><td>Related to position of the inflexion point</td><td>5.33 <sup>e</sup></td></tr><tr><td>σ</td><td>Related to position of the inflexion point</td><td>3 <sup>e</sup></td></tr><tr><td colspan="3"><strong>Host-pathogen genetic interaction</strong></td></tr><tr><td>G</td><td>Total number of major genes</td><td>{1; 2}</td></tr><tr><td>τ<sub>g</sub></td><td>Mutation probability for infectivity gene g <sup>f</sup></td><td>10<sup>-4</sup> g</td></tr><tr><td>τ<sub>w</sub></td><td>Mutation probability for aggressiveness component w <sup>f</sup></td><td>10<sup>-4</sup> g</td></tr><tr><td>Q<sub>w</sub></td><td>Number of pathotypes relative to aggressiveness component w</td><td>6 <sup>g</sup></td></tr><tr><td>ρ<sub>g</sub></td><td>Efficiency of major gene g</td><td>1.0 <sup>g</sup></td></tr><tr><td>ρ<sub>w</sub></td><td>Efficiency of quantitative trait w</td><td>0.5 <sup>g</sup></td></tr><tr><td>θ<sub>g</sub></td><td>Cost of infectivity of infective gene g</td><td>0.5 <sup>g</sup></td></tr><tr><td>θ<sub>w</sub></td><td>Cost of aggressiveness for component w</td><td>0.5 <sup>g</sup></td></tr><tr><td>β<sub>w</sub></td><td>Trade-off strength for aggressiveness component w</td><td>1.0 <sup>g</sup></td></tr></tbody></table>

<footnote><sup>a</sup> model parameterisation is detailed in S1 Text.</footnote>

<footnote><sup>b</sup> values for the five landscape structures.</footnote>

<footnote><sup>c</sup> same value for all cultivars (no cost of resistance).</footnote>

<footnote><sup>d</sup> g(||z' - z||) = (b-2)(b-1) / (2πσ²) · (1 + ||z' - z||)⁻¹ with ||z' - z|| the Euclidian distance between locations z and z' in fields i and i', respectively; the mean dispersal distance is given by: 2a/(b-3) = 20 m.</footnote>

<footnote><sup>e</sup> the position of the inflexion point of the sigmoid curve is given by the relation x₀ = ((σ - 1) / κσ)¹/² ≈ 0.5</footnote>

<footnote><sup>f</sup> probability for a propagule to change its infectivity or its aggressiveness on a resistant cultivar carrying major gene g or quantitative resistance trait w.</footnote>

<footnote><sup>g</sup> same value for all major genes, quantitative resistance traits, infectivity genes and aggressiveness components.</footnote>
<table><caption>Table 2. Infectivity matrix.</caption><thead><tr><th colspan="2">Infectivity matrix for major gene g (INF<sup>g</sup>)</th><th colspan="2">Host genotype v</th></tr><tr><th rowspan="2">Pathogen genotype p</th><th colspan="2">Susceptible (SC) mg<sub>g</sub>(v) = 0</th><th colspan="2">Resistant (RC) mg<sub>g</sub>(v) = 1</th></tr><tr><th>Non-infective ig<sub>g</sub>(p) = 0</th><th>1</th><th>1 - ρ<sub>g</sub></th><th>1</th></tr></thead><tbody><tr><th></th><th> Infective ig<sub>g</sub>(p) = 1</th><th>1 - θ<sub>g</sub></th><th></th><th>1</th></tr></tbody></table>

https://doi.org/10.1371/journal.pcbi.1006067.t002

Table 3). These matrices summarise the multiplicative effect of this interaction on the infection rate of the pathogen (w = e), the duration of the infectious period (w = Y) and the propagule production rate (w = r), or the divisor effect on the duration of the latent period (w = γ).

Here, $\rho_w$ is the efficiency of quantitative resistance with respect to aggressiveness component w for a non-adapted pathogen, $\theta_w$ the aggressiveness cost paid by adapted pathogens on susceptible hosts, and $Q_w$ the number of possible pathotypes (thus, quantitative resistance trait w is completely eroded in $Q_w$-1 steps). For non-adapted pathogens, every gain in aggressiveness on resistant hosts is penalised by a cost on susceptible hosts, according to the relationship [70]:

$$cost = 1 - \left(1 - gain^{\frac{1}{\beta_w}}\right)^{\beta_w} \quad (1)$$

with $\beta_w$ the strength of the trade-off for aggressiveness component w (Fig 2). $\beta_w<1$ means a weak trade-off (gain higher than cost), $\beta_w>1$ means a strong trade-off (gain lower than cost) and $\beta_w = 1$ represents a linear trade-off (gain equals cost). Thus, the aggressiveness associated with a pathogen genotype can be described by: $(1 - \rho_w) + gain \times \rho_w$ on resistant hosts; and by: $1 - cost \times \theta_w$ on susceptible hosts.

Note, there are $Q_w$ different pathotypes relative to aggressiveness component w, whereas there are only two pathotypes relative to each of the major infectivity genes. Thus the total number of pathogen genotypes is $P = 2^G \times \prod_{w = e, \gamma, Y, r} Q_w$. It is also noteworthy that the infectivity matrix (Table 2) shown here is a specific case of the more general aggressiveness matrix (Table 3 with $Q_w = 2$).

## Host and pathogen demo-genetic dynamics

The demo-genetic dynamics of the host-pathogen interaction are based on a SEIR structure. However, to avoid any confusion with the 'susceptible' cultivar, we have labelled this structure HLIR for 'healthy-latent-infectious-removed'. Thus, in the following, $H_{i,v,t}$, $L_{i,v,p,t}$, $I_{i,v,p,t}$, $R_{i,v,p,t}$ and $Pr_{i,p,t}$ respectively denote the number of healthy, latent, infectious and removed individuals (in this model, an 'individual' is a given amount of plant tissue, and is referred to as a 'host' hereafter for simplicity), and pathogen propagules, in field i (i = 1, ..., J), for cultivar v (v = 1, ..., V), pathogen genotype p (p = 1, ..., P) at time step t (t = 1, ..., T, X, Y). T is the number of time steps in a cropping season and Y the number of simulated years (i.e. cropping seasons). Since the host is cultivated, we assume there is no host reproduction, dispersal or natural

<table><caption>Table 3. Aggressiveness matrix.</caption><thead><tr><th colspan="2">Aggressiveness matrix for component w (AGG<sup>w</sup>)</th><th colspan="2">Host genotype v</th></tr><tr><th rowspan="2">Pathogen genotype p</th><th colspan="2">Susceptible (SC) qr<sub>w</sub>(v) = 0</th><th colspan="2">Resistant (RC) qr<sub>w</sub>(v) = 1</th></tr><tr><th>Non-aggressive ag<sub>w</sub>(p) = 1</th><th>1</th><th>1 - ρ<sub>w</sub></th><th>1</th></tr></thead><tbody><tr><th></th><th>...</th><th>...</th><th>...</th><th>...</th></tr><tr><th></th><th>...</th><th>...</th><th>...</th><th>...</th></tr><tr><th></th><th>Fully aggressive ag<sub>w</sub>(p) = Q<sub>w</sub></th><th>1 - θ<sub>w</sub></th><th></th><th>1</th></tr></tbody></table>

https://doi.org/10.1371/journal.pcbi.1006067.t003
<figure><img src="image_7.png" /><figcaption>**Fig 2. Trade-off relationship.** Levels of pathogen aggressiveness on resistant (RC) and susceptible (SC) cultivars are linked by a linear (solid curve, $\beta = 1$), a strong (dashed curve, $\beta = 1.5$ in this example), or a weak (dotted curve, $\beta = 0.5$) trade-off. The blue vertical line is related to resistance efficiency ($\rho_w = 0.9$) and the red horizontal line is related to the cost of aggressiveness ($\theta_w = 0.7$).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g002

mortality (leaf senescence near the end of the cropping season is considered as part of host harvest). Fig 3 gives a schematic representation of the model structure.

**Host growth.** Only healthy hosts (denoted as $H_{i,v,t}$) are assumed to contribute to growth of the crop (note, the model is flexible enough to relax this assumption; see Discussion). Thus, at each step t during a cropping season, the plant cover of cultivar v in field i increases as a

<figure><img src="image_8.png" /><figcaption>**Fig 3. Architecture of the simulation model.** Healthy hosts can be contaminated by propagules and may become infected. Following a latent period, infectious hosts start producing new propagules which may mutate and disperse across the landscape. At the end of the infectious period, infected hosts become epidemiologically inactive. Qualitative resistance usually prevents transition to the infected state, whereas quantitative resistance can affect several steps of the epidemic cycle but does not completely prevent infection. Green boxes indicate healthy hosts which contribute to crop yield and host growth, in contrast to diseased plants (i.e. symptomatic, red boxes) or those with latent infections (dark blue box).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g003
logistic function, and the new amount of healthy plant tissue is:

$$H_{i,v,t+1} = H_{i,v,t} \left[ 1 + \delta_v \times \left( 1 - \frac{N_{i,v,t}}{K_{i,v}} \right) \right] \quad (2)$$

with $\delta_v$ the growth rate of cultivar v (possibly, although not necessarily, lower for resistant cultivars due to fitness costs of resistance [71]); $N_{i,v,t} = H_{i,v,t} + \sum_{p=1}^{P} \{L_{i,v,p,t} + I_{i,v,p,t} + R_{i,v,p,t}\}$ the total number of hosts in field i for cultivar v and at time t; and $K_{i,v} = A_i \times C_v^{\max}$ the carrying capacity of cultivar v in field i, which depends on $A_i$, the area of the field, and $C_v^{\max}$, the maximal density for cultivar v. Note that when a mixture of several cultivars is present in the same field, decreased growth due to susceptible plants being diseased is not compensated by increased growth of resistant plants.

**Contamination of healthy hosts.** The healthy compartment (H) is composed of hosts which are free of pathogen propagules (H¹), as well as hosts contaminated (but not yet infected) by the arrival of such propagules (H²). At the beginning of each step, all healthy hosts are considered free of propagules (H¹). Then, as described in a previous modelling study [42], at time t in field i and for cultivar v, the number of contaminable hosts (i.e. accessible to pathogen propagules, denoted as $H_{i,v,t}^{\text{contaminable}}$) depends on the proportion of healthy hosts ($H_{i,v,t}^1$) in the host population ($N_{i,v,t}$):

$$H_{i,v,t}^{\text{contaminable}} \sim \text{Binomial} \left( H_{i,v,t}^1 ; \pi \left( \frac{H_{i,v,t}^1}{N_{i,v,t}} \right) \right) \quad (3)$$

with $\pi(x) = \frac{1-e^{-K_v \sigma}}{1-e^{-K}}$, a sigmoid function with $\pi(0) = 0$ and $\pi(1) = 1$, giving the probability for a healthy host to be contaminated. Here, we assume that healthy hosts are not equally likely to be contacted by propagules, for instance because of plant architecture. Moreover, as the local severity of disease increases, eventually the probability for a single propagule to contaminate a healthy host declines due to the decreased availability of host tissue.

Following the arrival of propagules of pathogen genotype p in field i at time t (denoted as $Pr_{i,p,t}^4$, see below), susceptible hosts become contaminated. The pathogen genotypes of these propagules are distributed among contaminable hosts according to their proportional representation in the total pool of propagules. Thus, for cultivar v, the vector describing the maximum number of contaminated hosts by each pathogen genotype (denoted as $[H_{i,v,t}^{\text{maxConta}}]_{p=1,...,p}$) is given by a multinomial draw:

$$[H_{i,v,t}^{\text{maxConta}}]_{p=1,...,p} \sim \text{Multinomial} \left( H_{i,v,t}^{\text{contaminable}} ; \left[ \frac{Pr_{i,p,t}^4}{\sum_{p=1}^{P} Pr_{i,p,t}^4} \right]_{p=1,...,p} \right) \quad (4)$$

However, the number of deposited propagules ($Pr_{i,p,t}^4$) may be smaller than the maximal number of contaminated hosts ($H_{i,v,p,t}^{\text{maxConta}}$). Thus, the true number of hosts of cultivar v, contaminated by pathogen genotype p in field i at t (denoted as $H_{i,v,p,t}^2$) is given by:

$$[H_{i,v,t}^1 \rightarrow H_{i,v,p,t}^2] = \min(H_{i,v,p,t}^{\text{maxConta}}; Pr_{i,p,t}^4) \quad (5)$$

**Infection.** Between t and t+1, in field i, contaminated hosts ($H_{i,v,p,t}^2$) become infected (state L) with probability $e_{v,p}$, which depends on the maximum expected infection efficiency, $e_{\text{max}}$,
and the interaction between host (v) and pathogen (p) genotypes:

$$[H_{i,v,p,t}^2 \rightarrow L_{i,v,p,t+1}] \sim \text{Binomial}(H_{i,v,p,t}^2; e_{v,p}) \qquad (6)$$

$$e_{v,p} = e_{max} \times AGG_{ag_e(p),qr_e(v)}^e \times \prod_{g=1}^{G} INF_{ig_g(p),mg_g(v)}^g \qquad (7)$$

**Latent period.** Infected hosts become infectious (state I) after a latent period (LI) drawn from a Gamma distribution (a flexible continuous distribution from which durations in the interval [0; +∞[ can be drawn) parameterised with an expected value, $\gamma_{exp}$, and variance, $\gamma_{var}$. The expected duration of the latent period depends on the minimal expected duration, $\gamma_{min}$, and the interaction between host (v) and pathogen (p) genotypes:

$$(LI)_{v,p} \sim \text{Gamma}(\gamma'_{expv,p}; \gamma'_{var}) \qquad (8)$$

$$\gamma'_{expv,p} = \frac{\gamma'_{min}}{AGG_{ag_r(p),qr_r(v)}^r} \qquad (9)$$

Note, the usual shape and scale parameters of a Gamma distribution, $\beta_1$ and $\beta_2$, can be calculated from the expectation and variance, exp and var, with: $\beta_1 = \frac{exp^2}{var}$ and $\beta_2 = \frac{var}{exp}$, respectively.

**Infectious period.** Finally, infectious hosts become epidemiologically inactive (i.e. they no longer produce propagules, thus are in state R, 'removed') after an infectious period (IR) drawn from a Gamma distribution parameterised with expected value, $Y_{exp}$, and variance, $Y_{var}$, similar to the latent period. The expected duration of the infectious period depends on the maximal expected duration, $Y_{max}$, and the interaction between host (v) and pathogen (p) genotypes:

$$(IR)_{v,p} \sim \text{Gamma}(Y_{expv,p}; Y_{var}) \qquad (10)$$

$$Y_{expv,p} = Y_{max} \times AGG_{ag_Y(p),qr_Y(v)}^Y \qquad (11)$$

**Pathogen reproduction.** In field i at time t, infectious hosts associated with pathogen genotype p produce a total number of propagules (denoted as $Pr_{i,p,t}^1$), drawn from a Poisson distribution whose expectation, $r_{exp}$, depends on the maximal expected number of propagules produced by a single infectious host per time step, $r_{max}$, and the interaction between host (v) and pathogen (p) genotypes:

$$Pr_{i,p,t}^1 \sim \text{Poisson}(\sum_{v=1}^{V} r_{expv,p}) \qquad (12)$$

$$r_{expv,p} = r_{max} \times AGG_{ag_r(p),qr_r(v)}^r \times I_{i,v,p,t} \qquad (13)$$

**Pathogen mutation.** The following algorithm is repeated independently for every potential infectivity gene g:

1. the pathotype (i.e. the level of adaptation with regard to major gene g, indexed by q; q = 1, ..., $Q_g$; with $Q_g = 2$ since the pathotype is either infective, or non-infective) of the pathogen propagules is retrieved from their genotype p;
2. propagules can mutate from pathotype q to pathotype q' with probability $m_{qq'}^g$ such as $m_{qq'}^g = \tau_g$ if $q' \neq q$ (hence $m_{qq}^g = 1 - \tau_g$ since $Q_g = 2$). Thus, in field i at time t, the vector of the number of propagules of each pathotype arising from pathotype q (denoted as $[M_{i,q,q',t}]_{q'=1,...,Q_g}$) is given by a multinomial draw:

$$[M_{i,q,q',t}^g]_{q'=1,...,Q_g} \sim \text{Multinomial}(Pr_{i,p,t}^1; [m_{qq'}^g]_{q'=1,...,Q_g}) \quad (14)$$

3. the total number of propagules belonging to pathotype q' and produced in field i at time t (denoted as $Pr_{i,q',t}^2$) is:

$$Pr_{i,q',t}^2 = \sum_{q'=1}^{Q_g} M_{i,q,q',t}^g \quad (15)$$

4. the new propagule genotype p' is retrieved from its new pathotype (q'), and the number of propagules is incremented using a variable denoted as $Pr_{i,p',t}^3$.

The same algorithm is also repeated independently for every potential aggressiveness component w. In this case, $\tau_g$ is replaced by $\tau_w$ and $Q_g$ by $Q_w$. In addition, pathogen evolution is assumed to be stepwise, i.e. in a single step, a given pathotype can only mutate to a closely related pathotype (with regard to aggressiveness component w). Thus, we have $m_{q(q-1)}^w = m_{q(q+1)}^w = \tau_w/2$ and $m_{qq}^w = 1 - \tau_w$. The exceptions are pathotypes fully adapted to either the susceptible host or the resistant host. These mutate towards less specialised pathotypes with a probability of $\tau_w$ to ensure their overall mutation probability is equivalent to that of other genotypes. In this model, it should be noted that the mutation probability $\tau_g$ (or $\tau_w$) is not the classic 'mutation rate' (i.e. the number of genetic mutations per generation per base pair), but the probability for a propagule to change its infectivity (or aggressiveness) on a resistant cultivar carrying major gene g (or quantitative resistance trait w). This probability depends on the classic mutation rate, the number and nature of the specific genetic mutations required to overcome major gene g (or improve aggressiveness component w), and the potential dependency between these mutations.

**Pathogen dispersal.** Propagules can migrate from field i (whose area is $A_i$) to field i' (whose area is $A_{i'}$) with probability $\mu_{ii'}$, computed from:

$$\mu_{ii'} = \frac{\int_{A_i} \int_{A_{i'}} g(\|z' - z\|) dz dz'}{A_i} \quad (16)$$

with $\|z' - z\|$ the Euclidian distance between locations z and z' in fields i and i', respectively, and g(.) the two-dimensional dispersal kernel of the propagules. The computation of this probability is performed using the CaliFloPP algorithm [72]. Thus, at time t, the vector of the number of propagules of genotype p migrating from field i to each field i' (denoted as $[D_{i',p,t}]_{t'=1,...,J}$) is:

$$[D_{i',p,t}]_{t'=1,...,J} \sim \text{Multinomial}(Pr_{i',p,t}^3; [\mu_{i'j}]_{j=1,...,J}) \quad (17)$$

and the total number of propagules arriving in field i' at time t (denoted as $Pr_{i',p,t}^4$) is:

$$Pr_{i',p,t}^4 = \sum_{j=1}^{J} D_{i',p,t} \quad (18)$$
We consider that propagules landing outside the boundaries of the simulated landscape are lost (absorbing boundary condition), and there are no propagule sources external to the simulated landscape.

**Seasonality.** Let $t^0(y)$ and $t'(y)$ denote the first and last days of cropping season y ($y = 1, \dots, Y$), respectively. The plant cover in field i for cultivar v at the beginning of cropping season y is set at $H_{i,v,t^0(y)} = A_i \times C_v^0 \times \mathbb{I}_{V_{(i)}=v}$, with $C_v^0$ the plantation density of cultivar v and $\mathbb{I}_{V_{(i)}}$ an indicative variable set at 1 when field i is cultivated with cultivar v and 0 otherwise. At the end of a cropping season, the host is harvested. We assume that the pathogen needs a green bridge to survive the off-season. This green bridge could, for example, be a wild reservoir or volunteer plants remaining in the field (e.g. owing to incomplete harvest or seedlings). The size of this reservoir imposes a bottleneck for the pathogen population. The number of remaining infected hosts in field i for cultivar v and pathogen genotype p (denoted by $I_{i,v,p,t'(y)}^*$) at the end of the off-season is given by:

$$I_{i,v,p,t'(y)}^* \sim \text{Binomial}(L_{i,v,p,t'(y)} + I_{i,v,p,t'(y)}; \lambda) \quad (19)$$

with $\lambda$ the survival probability of infected hosts. Depending on host (v) and pathogen (p) genotypes, the number of propagules produced by the remaining hosts during their whole infectious period (denoted by $Pr_{i,p,t'(y)}^*$) is drawn from a Poisson distribution:

$$Pr_{i,p,t'(y)}^* \sim \text{Poisson} \left( \sum_{v=1}^{V} \{ r_{max} \times AGG_{ag_r(p),qr_r(v)}^r \times Y_{max} \times AGG_{ag_Y(p),qr_Y(v)}^Y \times I_{i,v,p,t'(y)}^* \} \right) \quad (20)$$

These propagules can mutate and disperse exactly as happens during the cropping season, and constitute the initial inoculum for the next cropping season.

**Initial conditions.** At the beginning of a simulation, healthy hosts are planted in each field as previously described. The initial pathogen population is assumed to be totally non-adapted to the resistance host cultivars, and is only present in susceptible fields with probability $\phi$. Then the initial number of infectious hosts in these fields is: $I_{i,v} = 1, p = 1, t = 1 = \text{Binomial}(H_{i,v} = 1, t = 1; \phi)$.

The model utilises multiple criteria to enable the assessment of different deployment strategies with regard to both evolutionary and epidemiological outcomes.

## Evolutionary outputs

**Durability of qualitative resistance.** For a given major gene, several computations are performed: (d1) time to first appearance of a pathogen mutant; (d2) time to first true infection of a resistant host by such mutants; and (d3) time when the number of infections of resistant hosts by these mutants reaches a threshold above which mutant pathogens are unlikely to go extinct. These metrics characterise the three critical steps to the breakdown of qualitative resistance: mutation toward infectivity, immigration (or introduction) on resistant hosts through dispersal, and subsequent broader establishment in the resistant host population.

**Erosion of quantitative resistance.** Since pathogen adaptation to quantitative resistance is gradual, the three measures described above are computed for every step towards complete erosion of resistance (i.e. $Q_w-1$ levels). This allows us to characterise the point in time when quantitative resistance starts to erode, the final level of erosion, and the speed of erosion (i.e. the percentage of erosion per time unit from the time when quantitative resistance begins to erode).

**Durability of a deployment strategy.** A simulation run is divided into three periods: (1) the initial short-term period when all resistance sources are at their highest potential; (2) a transitory period during which a given deployment strategy is only partially effective; and (3) a
longer-term period when all the resistances have been overcome or completely eroded. To assess the end of the short-term period, the time to establishment (durability measure (d3)) is computed for every major gene, and every quantitative trait at the first level of erosion (ag<sub>w</sub>(p) = 2). The minimal value of these measures, denoted by D<sub>1</sub>, delimitates short-term and transitory periods. Similarly, the time to establishment is computed for every major gene, and for every quantitative trait at the highest level of erosion (ag<sub>w</sub>(p) = Q<sub>w</sub>). The maximal value of these measures, termed D<sub>2</sub>, delimits transitory and long-term periods. Fig 4 provides an example of a simulation run and the delimitation of these periods.

## Epidemiological outputs

The epidemiological impact of pathogen spread is evaluated by two different measures. Firstly, we use a measure termed Green Leaf Area (GLA), based on the Healthy Area Duration initially developed by Waggoner et al. [73]. The GLA represents the average number of productive hosts per time step and per surface unit, and is considered as a proxy for crop yield [74]. Secondly, we use the area under the disease progress curve (AUDPC), which is the average proportion of diseased hosts relative to the carrying capacity and represents disease severity [74]. We assume only healthy hosts (state H) contribute to crop production (Fig 3). Thus, for cultivar v during year y:

$$GLA_{v,y} = \frac{\sum_{t=t^0(y)}^{t^f(y)} \sum_{i=1}^{J} H_{i,v,t}}{T \times \sum_{i=1}^{J} A_i} \qquad (21)$$

$$AUDPC_{v,y} = \frac{\sum_{t=t^0(y)}^{t^f(y)} \sum_{i=1}^{J} \sum_{p=1}^{P} \{I_{i,v,p,t} + R_{i,v,p,t}\}}{\sum_{t=t^0(y)}^{t^f(y)} \sum_{i=1}^{J} K_{i,v}} \qquad (22)$$

<figure><img src="image_9.png" /><figcaption>Fig 4. Computation of output variables in a simulated example. Two major resistance genes are deployed as a mosaic composed of a susceptible cultivar (solid curve) and two resistant cultivars (dotted and dashed curves) carrying the two major genes. The dynamics of the proportion of healthy (A) or diseased (B) hosts is integrated every year into the Green Leaf Area (GLA) or the area under disease progress curve (AUDPC), respectively. The vertical blue lines mark the times to breakdown of the first (dotted line) and the second (dashed line) major genes. These time points delimit the short-term (green zone), transitory (grey) and long-term (red) phases of resistance breakdown.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g004
**Global epidemiological control.** The GLA and AUDPC of every cultivar ($GLA_v$ and $AUDPC_v$) as well as the whole landscape ($GLA_{TOT}$ and $AUDPC_{TOT}$) are averaged across the whole simulation run, to measure the global epidemiological performance of a deployment strategy.

**Short-term epidemiological control.** The average GLA and AUDPC of the susceptible cultivar is computed on whole cropping seasons from the beginning of the simulation until the end of the season preceding $D_1$ (see 'Evolutionary outputs' for the computation of $D_1$ and $D_2$), and denoted as $GLA_{ST}$ and $AUDPC_{ST}$, respectively (Fig 4, green zone). These outputs cannot be computed if a resistance is overcome or starts to erode before the end of the first cropping season.

**Epidemiological control during the transitory period.** The average GLA and AUDPC of the susceptible cultivar is computed on whole seasons from the beginning of the season following $D_1$ to the end of the season preceding $D_2$, and denoted as $GLA_{TP}$ and $AUDPC_{TP}$, respectively (Fig 4, grey zone). These outputs cannot be computed if there is not at least one complete season between $D_1$ and $D_2$.

**Long-term epidemiological control.** The average GLA and AUDPC of the whole landscape is computed on whole seasons from the beginning of the season following $D_2$ to the end of the simulation, and denoted as $GLA_{LT}$ and $AUDPC_{LT}$, respectively (Fig 4, red zone). These outputs cannot be computed if at least one of the resistances is not overcome or completely eroded by the end of the simulation.

## Application case: Durability and efficiency of combinations of different resistance sources

In this section, ten scenarios are simulated using the model to assess the evolutionary and epidemiological outcomes of deploying different pyramided combinations of qualitative and quantitative resistances. In this context, the model was parameterised to approximate biotrophic foliar fungal diseases as typified by rusts of cereal crops, caused by fungi of the genus *Puccinia*.

**Model parameterisation for rust diseases.** Within these pathosystems, propagules (called 'spores') are produced by sporulating lesions which develop on the leaves of infected hosts after a latent period of a few days, following which they are dispersed by wind. Thus, in this case, an 'individual host' in the model can be considered as a foliar site where a spore can land and potentially trigger the development of a localised infection. The parameter values associated with pathogen aggressiveness components (infection rate, latent period duration, sporulation duration, sporulation rate) were estimated using available data from the literature [29, 30, 33, 37, 39-42, 45, 50, 55, 56, 60, 65, 75-116]; other parameters were arbitrarily fixed (see Table 1 for parameter values and S1 Text for details). In addition, the time step was set at one day and simulations were run over 50 cropping seasons of 120 days per year. Given this parameterisation, the threshold for mutant pathogen establishment (and thus resistance breakdown) was set at 50,000 infections. Above this threshold, the probability of extinction of a pathogen genotype is below 1% (see S2 Text for details of how the threshold was calculated).

**Simulated scenarios.** A set of five landscape structures of about 150 fields within a total area of 2x2 km² was generated (see S1 Fig). Areas of the simulated fields ranged from 0.36 to 5.38 ha (mean: 2.60 ha). Of these fields, 80% were cultivated with a resistant cultivar, and the remaining 20% with a susceptible cultivar, with a low level of spatial aggregation (landscape typified by Fig 1D). In the first five simulated scenarios, the resistant cultivar carried a single resistance source: either a major gene (conferring complete immunity against non-infective pathogens, $\rho_g$ = 1, as typically described in gene-for-gene scenarios for many rusts [62]), or
one of the four traits for quantitative resistance (against infection rate, latent period duration, sporulation rate, or sporulation duration). The efficiency of quantitative resistance traits was set at $ρ_w = 0.5$, meaning that infection rate, sporulation rate or the duration of the sporulation period of non-aggressive pathogens on resistant hosts was reduced by 50%. Similarly, the latent period duration was increased in such a way that the number of epidemic cycles of a non-aggressive pathogen in a cropping season was reduced by 50% on resistant hosts. In the next five scenarios, the previous resistance sources (a major gene or one of the four quantitative resistance traits) were combined with a major gene to represent a single pyramided cultivar. In all scenarios, the mutation probabilities were set at $τ_g = τ_w = 10^{-4}$ in such way that a cultivar carrying a single major gene would be overcome in less than one year. Each scenario was simulated ten times on all five landscape structures, resulting in 50 replicates per scenario, and a total of 500 simulations.

The model is written using the C and R languages, and is available via the R package land-sepi (v0.0.2, [117]). Within the R software (v3.4.0, [118]) landscape structures were generated using the package RLiTe (http://kien-kieu.github.io/lite/rlite.html). The cultivar allocation algorithm used positive definite matrices generated with the *Exponential* function of the package fields (v8.10, [119]) and the *nearPD* function of the package Matrix (v1.2–11, [120]). The CaliFloPP algorithm was performed using the package RCALI (v0.2–18, [72]). Depending on the number of resistance sources, one simulation run takes 60 to 180 seconds on a standard desktop computer (2 cores, 2.30 GHz). Simulation results are available in supporting information (S1 File).

## Results

Several different scenarios were simulated to investigate three specific questions relative to the evolutionary and epidemiological outcomes of deploying combinations of qualitative and quantitative resistances against rust diseases of cereal crops. In this context, we varied landscape structure (i.e. spatial arrangement of field boundaries, see S1 Fig) but not its composition (i.e. proportion and spatial aggregation of the resistant cultivar). Landscape structure was never significant, thus in what follows we consider that every scenario was replicated 50 times.

## What is the durability of a major gene alone or combined with another source of resistance?

To assess the durability of a single major gene, the time period from the beginning of the simulation to appearance of mutants able to overcome the major gene (d1), first infection of the resistant cultivar (d2), and broader establishment on the resistant cultivar (d3) were used. Under the simulated conditions, mutants always appeared, dispersed across the landscape and established on a resistant cultivar carrying one major gene in less than one year (Fig 5A, ‘MG1’). When the resistant cultivar carried a second major gene (i.e. as a pyramid with two major genes), the first infection of a resistant host was delayed to 8.3 years (90% central range, $CR_{90}$: 0.4–23.4) on average, and the pathogen population was not established before, on average, 20.7 years ($CR_{90}$: 0.9–50.0) (Fig 5A, ‘MG2’).

In contrast, when a major gene was combined with a quantitative resistance with a 50% efficiency, the breakdown of the major gene was nearly as quick as if the major gene was alone, regardless of the pathogen life-history trait targeted by the quantitative resistance (Fig 5A, column 3 to 6, and inset). This conclusion is consistent with the results of an experimental study on pepper resistance against root-knot nematode [121], but differs from those of other studies carried out on different pathosystems, showing that a quantitative resistant background can significantly increase the durability of a cultivar carrying a major gene for resistance [122–
<figure><img src="image_10.png" /><figcaption>**Fig 5.** Evolutionary outcomes for *Puccinia rusts* after 50 years: Durability of the major gene. Time to appearance of mutants carrying the associated infectivity gene (white segment of bars), to first infection (light grey) or to establish (dark grey) on resistant hosts. The resistant cultivar carries a single major gene (MG1, efficiency $ρ_1$ = 100%), a combination of two major genes (MG2, efficiency $ρ_2$ = 100%) or the single major gene is combined with one of several quantitative resistance traits (columns 3 to 6, efficiency $ρ_w$ = 50% in A and $ρ_w$ = 90% in B). Inset: enlargement of scenarios showing short durability. Every scenario is replicated 50 times. Vertical bars represent the 90% central range (represented in grey for the time to first infection, for legibility).</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g005

124]. In order to test if this difference could be due to our assumption that quantitative resistance has a 50% efficiency, we replicated the numerical experiment with higher efficiencies ($ρ_w$ = 60%, 70%, 80%, 90%). The results of these new simulations indicated that when quantitative resistance efficiency is higher than 80%, durability of the major gene is improved in combination with quantitative resistance, especially if latent period is targeted (time to establishment delayed to 5.2 years, $CR_{90}$: 1.7–10.8, see S3A Fig). Furthermore, above 90% efficiency, quantitative resistance can increase the durability of the major gene compared to pyramiding of two major genes (average time to establishment between 11.0 and 30.9 years depending on the targeted trait, Fig 5B).

**What is the level and speed of erosion of quantitative resistance alone or combined with a major gene?**

When deployed alone, quantitative resistance was on average eroded by 26.4% ($CR_{90}$: 20–40), 20.8% ($CR_{90}$: 0–40), 22.0% ($CR_{90}$: 0–40), and 19.2% ($CR_{90}$: 0–40) by the end of the simulation, for resistances targeting respectively infection rate, duration of the latent period, sporulation rate and duration of the sporulation period (Fig 6). The associated average speeds of erosion were 2.16, 0.78, 1.44 and 1.02% per year, respectively. The targeted pathogen life-history trait appeared to have a significant effect on the final level (Kruskal Wallis $χ^2$ tests with 3 degrees of freedom, $p$ = $8.10^{-9}$) and speed ($p$ = $6.10^{-5}$) of resistance erosion. The combination with a major gene significantly affected the final level (Kruskal Wallis $χ^2$ test with 1 df, $p$ = $8.10^{-3}$), but not the speed ($p$ = 0.85) of erosion.

It is important to remember that the speed of erosion was computed from the time point when quantitative resistance started to erode (hence, following breakdown of the major gene) and not from the beginning of simulations. Assessing the average speed of erosion from the beginning of the simulations would greatly impact the results if the major gene was durable for several years. As previously described, this is not the case with a 50% efficiency of quantitative
<figure><img src="image_11.png" /><figcaption>**Fig 6.** Evolutionary outcomes for *Puccinia rusts* after 50 years: Final level of erosion of quantitative resistance traits. Quantitative resistance (efficiency $ρ_w$ = 50%) is deployed alone (top row) or in combination with a major resistance gene (bottom row). The red shading indicates the average speed of erosion from the time when quantitative erosion starts to erode to the time when the final level of erosion is reached, with darker shades representing faster erosion rates. Lower case letters indicate statistically different groups for final level (black) and speed (white) of erosion, according to pairwise comparison tests (Dunn's test with Holm correction). Every scenario is replicated 50 times. Vertical bars represent the 90% central range.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g006

resistance. However, with a 90% efficiency, the average speed of erosion of quantitative resistance from the beginning of the simulation was 2.5 to 3.0% per year when deployed alone, and only 0.7 to 1.7% per year when combined with a major gene. In this context, the combination of quantitative resistance with a major gene significantly delayed the start of quantitative resistance erosion (Kruskal Wallis $\chi^2$ test with 1 df, $p<10^{-15}$).

**What is the severity of epidemics when qualitative and quantitative resistances are deployed alone or in combination?**

The AUDPC of the susceptible cultivar (AUDPC<sub>SC</sub>), the resistant cultivar (AUDPC<sub>RC</sub>) and across the entire cropping landscape (AUDPC<sub>TOT</sub>) averaged across the whole simulation period were used as indicators of the severity of epidemics in the simulated scenarios (see S4 Fig for examples). As expected, regardless of the cultivar or the deployment scenario, disease severity in a landscape where a resistant cultivar is deployed was lower than in a fully susceptible landscape (dashed line in Fig 7).

When only a single resistance source was incorporated into the resistant cultivar, the most effective source was quantitative resistance against duration of the latent period (AUDPC<sub>RC</sub> = 0.25), followed by quantitative resistance against infection rate, sporulation rate and duration of the sporulation period (AUDPC<sub>RC</sub> = 0.32, 0.33 and 0.34, respectively), and finally major gene resistance (AUDPC<sub>RC</sub> = 0.37). When a major gene was combined with one of these
<figure><img src="image_12.png" /><figcaption>**Fig 7.** Epidemiological outcomes for *Puccinia rusts* after 50 years. Average area under disease progress curve (AUDPC) for a susceptible (white bars) and a resistant (dark grey bars, cropping ratio: 80%) cultivar, as well as for the whole landscape (light grey bars). The top row represents situations where the resistant cultivar carries a single major gene (MG, efficiency ρ₁ = 100%) or a single quantitative resistance trait (columns 2 to 5, efficiency ρ<sub>w</sub> = 50%). The bottom row shows results when a major resistance gene was added to each of the single resistance scenarios (e.g. the leftmost set of bars on the bottom row shows results for two major genes in combination; the rightmost set is when a major gene was combined with quantitative resistance against sporulation duration). Each scenario was replicated 50 times. The horizontal dashed line represents the average AUDPC in a fully susceptible landscape. Vertical bars represent the 90% central range.</figcaption></figure>

https://doi.org/10.1371/journal.pcbi.1006067.g007

resistance sources, the combination of two major genes became the most efficient strategy (AUDPC<sub>RC</sub> = 0.22) but also the most variable (CR<sub>90</sub>: 0.00-0.37). Below this, were combinations including quantitative resistance traits in the same order as above and with similar efficiencies (AUDPC<sub>RC</sub> = 0.24, 0.31, 0.32, and 0.33 against latent period, infection rate, sporulation rate and sporulation duration, respectively).

For the susceptible cultivar, disease severity was less variable between scenarios than for the resistant cultivar. The average AUDPC<sub>SC</sub> was 0.33, 0.37, 0.37, 0.36 and 0.36, respectively when a major gene or a quantitative resistance trait against infection rate, latent period, sporulation rate or duration of the sporulation period was deployed alone. However, when a major gene was combined with these sources, AUDPC<sub>SC</sub> was smaller (0.28, 0.30, 0.29, 0.29 and 0.30, respectively).

Globally, across the whole landscape, disease severity (AUDPC<sub>TOT</sub>), was very similar to levels of disease severity seen in the resistant cultivar (AUDPC<sub>RC</sub>), since this cultivar constituted a high proportion of the area being cropped (80%, see Fig 1D for an example). The global severity of disease could also be assessed using the GLA (focusing on healthy hosts) instead of the AUDPC (focusing on diseased hosts). In our context, GLA<sub>TOT</sub> was highly negatively correlated with AUDPC<sub>TOT</sub> (Pearson correlation coefficient of -0.92, p<10<sup>-15</sup>).
Discussion

This article describes the first spatiotemporal model with the ability to flexibly simulate the evolution of a plant pathogen in a cropping landscape, following the deployment of resistance within the four main categories of deployment strategies: mosaics, mixtures, rotations and pyr-amiding. In addition, the model includes every possible combination of different major genes (qualitative resistance conferring immunity to the plant) and traits for quantitative resistance against the main aggressiveness components of the pathogen (infection rate, latent period, infectious period, and reproduction rate). Because of its flexible parameterisation and ability to simulate an explicit landscape, the model offers the possibility to vary the number and type of deployed resistance sources, their relative proportion in surface coverage across the land-scape and their level of spatial (or temporal) aggregation, as well as the epidemiology and the evolutionary potential of the pathogen. In addition, although the present study used a simu-lated landscape in order to control particular features of the fields, the model also permits the use of real landscapes to match more specific applied contexts. This model is implemented in the package landsepi [117] for the statistical computing environment R. Consequently, it pro-vides a new and useful tool to assess the performance of a wide range of deployment options. In particular, it was used here to evaluate the combination of different sources of resistance, and it will be used in future studies to explore different spatio-temporal deployment strategies.

With respect to criteria used to evaluate the performance of different deployment strategies, these may vary depending on needs of different stakeholder groups (e.g. breeders and growers). Thus, the model generates a panel of outputs describing the epidemiological (i.e. the ability of different resistance sources to reduce the severity of epidemics and consequently their impact on crops) and the evolutionary (relative to the durability of these resistance sources) outcomes of a given deployment strategy.

Owing to similar epidemiological and evolutionary concepts, there are interesting parallels between the deployment of plant resistance to pathogens and the application of pesticides (for plants) or vaccines (for animals, including humans). Indeed, the proportion of resistant fields in our study is equivalent to the proportion of fields treated with a pesticide [29] or the proportion of vaccinated individuals in an animal population [125]. The efficiency of plant resistance is similar to pesticide dose [46], and the type of resistance may be compared to the type of pesticide or vaccine. For instance, quantitative resistance traits against pathogen infection rate, latent period and sporulation rate are analogous to imperfect vaccines with anti-infection, anti-growth or anti-transmission modes of action, respectively [126]. It is noteworthy that spatio-temporal deployment strategies for plant resistance have counterparts in the application of pesticides. Thus, mosaics of different resistant crops are analogous to mosaics of fields treated with different pesticides, rotations of resistant crops to periodic application of different pesticides, and crop mixtures and pyramids of resistance genes to combinations of pesticide molecules [53, 54].

The model is flexible enough to facilitate investigation of a wide range of host-pathogen sys-
tems. Although our application case focused on biotrophic rust fungi of cereal crops, other
pathogens transmitted by wind, rain or vector insects (e.g. many fungi, bacteria and viruses),
can be simulated by modifying in particular the dispersal kernel and the value of parameters
associated with pathogen aggressiveness components (see for example how to parameterise the
model to necrotrophic fungal pathogens in S1 Text). Nevertheless, some fixed features of the
model (e.g. clonal reproduction of the pathogen) may need to be modified if sexual reproduc-
tion (or recombination for viruses) significantly contribute to pathogen genetic diversity and
dynamics. In addition, modelling vector-borne bacteria and viruses transmitted in a persistent
manner may require incorporation of vector dynamics and behaviour into the model frame-
work. Future work on the model is planned to address some of these current limitations.
## A spatially explicit stochastic model of pathogen evolution

Three main properties distinguish models simulating the deployment of resistance and the interpretation of their results: whether they are demographic or demo-genetic, whether they are deterministic or stochastic, and whether they are spatial or not.

In contrast with purely demographic models, our model includes the genetic evolution of the pathogen in addition to pathogen population dynamics and host growth during a cropping season. This enables the explicit simulation of the appearance of new genotypes through mutation [29, 30, 46, 48], which represents the first step towards resistance breakdown or erosion (i.e. prior to migration to and broader establishment onto resistant hosts). The time required for the achievement of this first step can be the main determinant of the durability of some deployment strategies, like pyramiding [51]. This step is thus essential to comprehensively assess the relative performance of different deployment strategies.

We further note that the present model is stochastic, i.e. it relies on probabilistic computations of simulated biological processes. It is consequently well able to account for biologically realistic random events [127]. As illustrated by Lo Iacono et al. [35], who used a stochastic version of the model developed by van den Bosch et al. [39], the likelihood of extinction events in pathogen populations can considerably impact the performance of a deployment strategy.

Finally, because our model is able to simulate explicit landscapes, it accounts for spatial heterogeneity, which affects landscape connectivity and consequently the ability of the pathogen to disperse. Pathogen dispersal, in addition to being one of the unavoidable steps of adaptation to plant resistance, strongly shapes pathogen evolutionary dynamics [12, 41, 128]. Moreover, the spatial nature of this model enables a wide range of deployment options to be considered at different spatial scales, particularly given the possibility to vary deployment parameters like the proportion of the landscape across which resistant cultivars are planted [32, 51], as well as the level of spatial aggregation [28, 45, 48, 49, 129, 130].

## Computation of epidemiological and evolutionary outputs

Resistance durability (i.e. the period of time from initial deployment of a resistant cultivar to when resistance is considered to have been overcome) is a typical evolutionary target of resistance deployment strategies, but its computation in a model is not obvious. Proposed methods include targeting the point in time when the first adapted pathogens appear [29, 30, 39], when their prevalence [37, 48, 55, 131] or frequency in pathogen population [39, 51, 52] exceeds a threshold, or when productivity of the resistant cultivar drops below an arbitrary threshold [45].

Since different measures give different information, the present model includes several measures of durability: the time until first appearance of mutants, the time until first infection of a resistant host by such mutant, and the time until the prevalence of these mutants exceeds a threshold. The first measure assesses the ability of deployment strategies to reduce the probability of appearance of mutants by reducing pathogen population size. The difference between the first and the second measure provides information on the potential to hinder pathogen migration to resistant hosts, and the difference between the second and the last measure is related to the rate of establishment on resistant cultivars as a result of the balance between selection and genetic drift (occurring between seasons in this model). With respect to the establishment process, it is interesting to note that different stages of invasion can be targeted by changing the value of the threshold. Importantly, we considered that time to establishment was best computed using a prevalence threshold (i.e. total number of resistant hosts infected by adapted pathogens) as opposed to using a frequency-dependent threshold in the pathogen population (i.e. the proportion of adapted pathogens in the global pathogen population). This
is because the frequency of adapted pathogens may never (or always) exceed a threshold simply as a result of a very low (or very high) proportion of resistant hosts in the landscape relative to susceptible hosts.

As pointed out in the introduction, evolutionary and epidemiological outcomes are not necessarily correlated. In this context, several epidemiological outputs are used in this study to characterise the level of protection provided by a deployment strategy against the potential damage caused by an epidemic. Some previous studies focused on the final state of the simulations or on the point at which a stable evolutionary equilibrium is reached, thus proposing criteria related to the final proportion of healthy [43, 44] or infected [26, 49, 129] hosts. However, in addition to long-term measures related to stable equilibria, short-term and transitory periods measures are important because severe epidemics responsible for heavy losses may occur during these stages [125, 132]. These periods can be accounted for by averaging the number of healthy hosts over the whole simulation run using an analogy of our Green Leaf Area [35–37, 40, 42, 133] or the number of infected hosts using the area under disease progress curve (AUDPC) [28, 32, 33, 48, 50]. Interestingly, these measures offer the possibility to concentrate on different evolutionary phases, for example the short-term period following resistance deployment until resistance breakdown, and the long-term period once resistance is overcome [45].

In the present study, the proposed outputs are based on both GLA and AUDPC, computed not only for the whole simulation run to have a global snapshot of the epidemiological outcome, but also for different time periods: from initial resistance deployment until the first resistance is overcome, from the time when all resistances are overcome until the end of the simulation as well as during the transitory period. The performance of various deployment strategies can differ during these three periods, since they are associated with different epidemiological contexts. Thus, it is important to consider these different measures together. Moreover, the objectives of different stakeholders define whether such criteria should be computed from the AUDPC or from the GLA. Indeed, if the objective is to limit the amount of disease, the use of AUDPC-based variables may be more appropriate, since they target infected hosts. In contrast, GLA-based variables represent the amount of healthy hosts which generally represent the largest contribution to crop yield. Thus, this variable can be useful when considering the impact of a given deployment strategy on productivity, especially when resistance costs, different planting densities, or different host species are involved. It is interesting to note that here, only healthy hosts were assumed to participate in host growth and contribute to final yield, considering that the infection by a pathogen consumes host resources. For rust pathogens of cereal crops, this assumption seems reasonable since it has been shown that manual defoliation of wheat leaves decreases yield less than infection by *Puccinia striiformis* [73]. Nevertheless, this assumption can easily be changed by including hosts at different sanitary stages (e.g. latently infected but not yet diseased) in the logistic equation of host growth and in computation of the GLA.

## Combining different sources of resistances

In this study, we present an initial exploration of the model where we evaluate the potential of combining qualitative and quantitative resistance to control rust of cereal crops. We note that in this exploration some parameters (in particular mutation probabilities, see S1 Text) were arbitrarily fixed to study a simple and theoretical situation where a resistant cultivar carrying a single major gene would be rapidly overcome following deployment. This is because our intent in this case was to compare different resistance combinations rather than provide an absolute prediction of the durability and efficiency of a particular strategy. Regardless, altering the
mutation probabilities, as long as they are the same for all infectivity genes and aggressiveness components, changes our results quantitatively but not qualitatively (see S1 Text and the results of simulation performed with smaller mutation probabilities in S5 Fig).

**Durability of qualitative resistance.** We found greater durability of a major resistance gene when combined with another major gene than when combined with a quantitative resistance trait exhibiting moderate efficiency. The high durability of pyramided major genes is often explained by the low probability of the pathogen simultaneously acquiring the required mutations to infect the resistant cultivar, added to the severe fitness costs associated with these mutations [64]. In the scenario we evaluated, it seems that this probability was not low enough to prevent the rapid appearance of adapted mutants (Fig 5A, ‘MG2’, durability measure (d1)). On the other hand, the great difference between durability measures (d1) and (d2) indicates that a long time passed before these mutants were able to disperse to resistant fields and infect a resistant host. This delay can be explained by the low probability that a spore simultaneously mutates and disperses to a resistant field, added to the severe fitness costs imposed on mutant pathogens in susceptible fields (note, in our simulations the infection rate of a pathogen carrying two infectivity genes is reduced by 75% on susceptible hosts). Furthermore, the great difference between durability measures (d2) and (d3) suggests that many extinction events impeded mutant establishment on the resistant cultivar, likely due to bottlenecks between cropping seasons.

Conversely, the combination of a major resistance gene with a highly efficient source of quantitative resistance can considerably increase the durability of the major gene. It is noteworthy that the durability of the major gene becomes also more variable (compare the 90% central ranges in Fig 5 and S3A Fig). This is mainly attributed to an increasing proportion of simulations where the major gene has a very long durability or is still not overcome by the end of the 50-year simulated period (e.g. in 42% of the simulations performed with a major gene combined with a quantitative resistance targeting pathogen latent period; see Fig 5B ‘Latent period’).

**Durability of quantitative resistance.** Our results show that the presence of a major gene in plant genotype does not greatly affect pathogen adaptation to a moderately effective quantitative resistance. However, highly effective quantitative resistance delays the breakdown of the major gene which at the same time delays the start of quantitative resistance erosion. Once quantitative resistance begins to erode, one would expect complete erosion over time, as observed in experimental serial passages performed with a plant virus [134]. In our results, quantitative resistance was not completely eroded by the end of the simulations. This was not due to the number of simulated years, since a speed of erosion of about 2% per year should enable complete erosion within 50 years of simulation. Moreover, it is interesting to note that quantitative resistance was eroded by approximately 25% when the initial resistance efficiency was 50%, and eroded by about 50% when the initial efficiency was 90% (Fig 6 and S3C Fig). This suggests that, in our simulations, pathogen adaptation always converged towards the same level of aggressiveness, which is likely a fitness optimum, given the fitness costs associated with pathogen adaptation, the shape of the trade-off function, and the proportion of the resistant cultivar in the landscape.

**Epidemiological outcomes for combined qualitative and quantitative resistances.**

Complementing a moderately efficient quantitative resistance with a major gene seems to have a greater effect on disease severity on the susceptible cultivar than on the resistant cultivar. This is probably because under these conditions the major gene has a very short durability (see above and Fig 5A), and only slightly delays the erosion of quantitative resistance (see above and Fig 6). On the other hand, pathogen evolution towards increased infectivity and aggressiveness on the resistant cultivar (which is strongly favoured in our simulations since 80% of
the landscape is grown with a resistant cultivar) triggers a stronger fitness cost on the susceptible cultivar, than for situations when only the evolution of increased aggressiveness is necessary to adapt to host resistance. This results in decreased severity on the susceptible cultivar and consequently across the whole landscape. Our results also suggest that quantitative resistance against the duration of the latent period has the best potential to limit disease severity, especially when the resistance is highly efficient (see S6B Fig). This is in agreement with previous studies [8, 93, 135], which highlighted the key role of the number of epidemic cycles within a cropping season.

**General conclusions.** As expected, the combination of multiple resistance sources in a resistant cultivar results in improved evolutionary and epidemiological outcomes relative to the deployment of a single resistance. In addition to the greater evolutionary barrier this represents, such combinations are more likely to trigger high fitness costs for the pathogen, thus constraining its ability to be simultaneously well adapted to different cultivars. Our simulations indicate that the pyramiding of two major genes is highly durable and effective, because this completely blocks the infection of resistant hosts by non-adapted pathogens, and requires the simultaneous acquisition of two costly infectivity genes to be overcome. Nevertheless, despite its partial efficiency, quantitative resistance may provide good epidemiological control if it has a sufficiently high efficiency (see S6 Fig) and especially if it is complemented by a major gene, which can significantly delay the time when quantitative resistance starts to erode. In turn, quantitative resistance helps reduce disease severity once the major gene is overcome.

Using a demographic model, Pietravalle et al. [37] found that quantitative resistance against the infection rate or the propagule production rate of the pathogen are equally efficient. However, in a demo-genetic version of a similar model, Lo Iacono et al. [36] found a greater effect of quantitative resistance against pathogen infection rate. Our results are consistent with this second study. However, our model includes for the first time the four main pathogen aggressiveness components. According to our simulation results, the most promising trait for quantitative resistance is against the duration of the latent period, owing to its influence on the number of epidemic cycles the pathogen can complete during a growing season. Consequently, quantitative resistance targeting pathogen latent period could effectively supplement the use of major genes in breeding programs. Thus, latent periods of rust pathogens could be measured on different host cultivars, and quantitative trait loci (QTL) associated with longer latent periods could be identified through genome-wide association studies, and subsequently selected for in marker-assisted selection programs [136]. S1 Video illustrates the potential of the deployment of a cultivar combining a major gene and a quantitative trait against latent period.

It is important to remember that, based on our model assumptions and the values chosen for parameters associated with pathogen evolution, the scenario we simulated here was a relatively simple one. For example, we considered that the different quantitative traits evolve independently from each other by assuming that: (i) mutations are independent from each other; (ii) when a quantitative resistance trait is deployed, only the targeted aggressiveness component of the pathogen can evolve; and (iii) a fitness cost is paid on the susceptible cultivar on the same aggressiveness component. Nevertheless, there is considerable evidence that different pathogen aggressiveness components often vary in association with each other [8, 83], suggesting that they do not necessarily evolve independently. Several trade-offs between aggressiveness components have been described [6, 7], and such trade-offs may considerably impact the performance of a particular deployment strategy [56]. In addition, some aggressiveness component(s) not directly targeted by host resistance could counterbalance the affected component(s), like a compensation phenomenon (as simulated in some modelling studies [137, 138]), and the fitness cost on the susceptible cultivar could be paid on an aggressiveness component different from the evolving component. These aspects could be captured through the
use of statistical estimates of the genetic variance and covariance of pathogen life-history traits, issued from phenotype values (G-matrix [139]). In addition, even in the case of well documented pathosystems like rusts of cereal crops, calibration of some parameters can be challenging (e.g. effective sporulation rate and dispersal kernel, see S1 Text), either because data obtained in laboratory experiments may not be representative of what happens in the field, or simply because data are missing. To fill this gap and help calibrate simulation models, estimation models built with exactly the same parameters can be useful to estimate poorly known parameters from field data [140–142].

## Future research directions

In a recent review dealing with the combination of qualitative and quantitative resistances, Pilet-Nayel et al. [143] wrote “There is still a need to adequately choose resistance QTLs to create optimal combinations and limit QTL erosion”. The present work highlights that longer latent periods may be one promising target. Nevertheless, the general features of the model could be used to test if our conclusions hold in different contexts. Our simulations showed that the efficiency of quantitative resistance ($\rho_w$) has a strong impact on the durability and epidemiological control of plant resistance. Therefore, in future studies we plan to assess the influence of the efficiency of qualitative resistance by including partially effective major genes ($\rho_g<1$, [3]). We also envisage varying the mutation probabilities ($\tau_g$ and $\tau_w$), cost of infectivity ($\theta_g$), cost of aggressiveness ($\theta_w$), and number of steps to erode a trait for quantitative resistance ($Q_w$), to simulate various possible choices of major genes and traits for quantitative resistance. It would also be of interest to simulate and investigate the potential for combinations of several sources of quantitative resistance, whose performance on disease severity has been experimentally demonstrated [144]. Finally, we will also explore different spatiotemporal strategies to deploy plant resistance to pathogens. In particular, we plan to compare mosaics, mixtures, rotations and pyramids of genetic resistances with respect to their epidemiological and evolutionary outcomes. The challenge is to identify strategies for a given host-pathogen interaction that are not only durable and efficient, but also feasible and likely to be adopted. In this context, promising strategies identified with the simulation model could be experimentally tested in the field. In addition, broader consideration of the economic context will also help maximise the chances that identified strategies are considered for real deployment.

## Supporting information

* S1 Fig. Landscape structures used in this study. (PDF)

* S2 Fig. Spatial allocation of three cultivars across a cropping landscape. (PDF)

* S3 Fig. Evolutionary outcome when quantitative resistance is highly effective. (PDF)

* S4 Fig. Dynamics of diseased hosts (AUDPC) from different scenarios of resistance deployment. (PDF)

* S5 Fig. Evolutionary and epidemiological outputs when the mutation probabilities are smaller. (PDF)
S6 Fig. Epidemiological outcome when quantitative resistance is highly effective.

(PDF)

S1 Text. Model parameterisation to rust diseases.

(PDF)

S2 Text. Calculation of the threshold for pathogen establishment.

(PDF)

S1 Video. Dynamics of diseased hosts in the landscape following the deployment of a resistant cultivar (proportion: 80%; aggregation: low) carrying a major resistance gene (efficiency 100%) combined with quantitative resistance (efficiency 90%) against the latent period duration of the pathogen, in a 50-year simulation (one year per second). The left panel shows the area under disease progress curve (AUDPC). The vertical blue line marks the time to breakdown of the major gene; the vertical red line marks the point in time at which quantitative resistance begins to erode. The right panel shows the spatio-temporal dynamics of the epidemic in the simulated landscape.

(MP4)

S1 File. Compressed folder containing .txt files.

- *simul.txt* contains the simulation plan and the epidemiological outputs (AUDPC and GLA);

- *dur1mut.txt*, *dur1inf.txt* and *durEsta.txt* contain the evolutionary outputs computed for every of the 5 levels of resistance erosion and from the three measures of durability, respectively (number of time steps from the beginning of the simulation to: (d1) the first appearance of mutants; (d2) the first infection of a resistant host by a mutant; and (d3) the broader establishment of these mutants in the resistant host population).

(ZIP)

**Acknowledgments**

The authors thank Loïc Houde for computing assistance and Jeremy Burdon for stimulating discussions and for reviewing this manuscript. Simulations were performed using the Pearcey cluster within the CSIRO supercomputer platform.

**Author Contributions**

**Conceptualization:** Loup Rimbaud, Julien Papaïx, Luke G. Barrett, Peter H. Thrall.

**Funding acquisition:** Julien Papaïx, Luke G. Barrett, Peter H. Thrall.

**Investigation:** Loup Rimbaud, Julien Papaïx.

**Methodology:** Loup Rimbaud, Julien Papaïx.

**Software:** Jean-François Rey.

**Supervision:** Julien Papaïx, Luke G. Barrett, Peter H. Thrall.

**Writing – original draft:** Loup Rimbaud.

**Writing – review & editing:** Loup Rimbaud, Julien Papaïx, Luke G. Barrett, Peter H. Thrall.

**References**

1. Stuthman DD, Leonard KJ, Miller-Garvin J. Breeding crops for durable resistance to disease. Advances in Agronomy. Volume 95. San Diego, United States: Academic Press; 2007. p. 319–367.
2. Brown JKM. Durable resistance of crops to disease: A darwinian perspective. Annu. Rev. Phytopathol. 2015; 53(1):513–539. https://doi.org/10.1146/annurev-phyto-102313-045914. PMID: 26077539

3. Burdon JJ, Barrett LG, Rebetzke G, Thrall PH. Guiding deployment of resistance in cereals using evolutionary principles. Evol. Appl. 2014; 7(6):609–624. https://doi.org/10.1111/eva.12175 PMID: 25067946

4. Parlevliet JE. Durability of resistance against fungal, bacterial and viral pathogens; present situation. Euphytica. 2002; 124(2):147–156. https://doi.org/10.1023/a:1015601731446

5. Mundt CC. Durable resistance: A key to sustainable management of pathogens and pests. Infect. Genet. Evol. 2014; 27:446–455. https://doi.org/10.1016/j.meegid.2014.01.011 PMID: 24486735

6. Lannou C. Variation and selection of quantitative traits in plant pathogens. Annu. Rev. Phytopathol. 2012; 50(1):319–338. https://doi.org/10.1146/annurev-phyto-081211-173031 PMID: 22702351

7. Pariaud B, Ravigné V, Halkett F, Goyeau H, Carlier J, Lannou C. Aggressiveness and its role in the adaptation of plant pathogens. Plant Pathol. 2009; 58(3):409–424. https://doi.org/10.1111/j.1365-3059.2009.02039.x

8. Parlevliet JE. Components of resistance that reduce the rate of epidemic development. Annu. Rev. Phytopathol. 1979; 17(1):203–222. https://doi.org/10.1146/annurev.py.17.090179.001223

9. McDonald BA, Linde C. Pathogen population genetics, evolutionary potential, and durable resistance. Annu. Rev. Phytopathol. 2002; 40(1):349–379. https://doi.org/10.1146/annurev.phyto.40.120501.101443 PMID: 12147764

10. Lecoq H, Moury B, Desbiez C, Palloix A, Pitrat M. Durable virus resistance in plants through conventional approaches: a challenge. Virus Res. 2004; 100(1):31–39. https://doi.org/10.1016/j.virusres.2003.12.012 PMID: 15036833

11. Moury B, Fabre F, Montarry J, Janzac B, Ayme V, Palloix A. L'adaptation des virus de plantes aux résistances variétales. Virologie. 2010; 14(4):227–239. https://doi.org/10.1684/vir.2010.0311

12. Gilligan CA. Sustainable agriculture and plant diseases: an epidemiological perspective. Philos. Trans. R. Soc. B. 2008; 363(1492):741–759. https://doi.org/10.1098/rstb.2007.2181 PMID: 17827101

13. García-Arenal F, McDonald BA. An analysis of the durability of resistance to plant viruses. Phytopathology. 2003; 93(8):941–952. https://doi.org/10.1094/PHYTO.2003.93.8.941 PMID: 18943860

14. Johnson R. Genetic background of durable resistance. In: Lamberti F, Waller JM, Graaff NA, editors. Durable Resistance in Crops. Boston, MA: Springer New York; 1983. p. 5–26.

15. McDonald BA. How can we achieve durable disease resistance in agricultural ecosystems? New Phytol. 2010; 185(1):3–5. https://doi.org/10.1111/j.1469-8137.2009.03108.x PMID: 20088970

16. McDonald BA. Using dynamic diversity to achieve durable disease resistance in agricultural ecosystems. Trop. Plant Pathol. 2014; 39:191–196. https://doi.org/10.1590/S1982-56762014000300001

17. Zhan J, Thrall PH, Papaïx J, Xie L, Burdon JJ. Playing on a pathogen's weakness: using evolution to guide sustainable plant disease control strategies. Annu. Rev. Phytopathol. 2015; 53(1):19–43. https://doi.org/10.1146/annurev-phyto-080614-120040 PMID: 25938275

18. Curl EA. Control of plant diseases by crop rotation. The Botanical Review. 1963; 29(4):413–479. https://doi.org/10.1007/bf02860813

19. Mundt CC. Use of multiline cultivars and cultivar mixtures for disease management. Annu. Rev. Phytopathol. 2002; 40(1):381–410. https://doi.org/10.1146/annurev.phyto.40.011402.113723 PMID: 12147765

20. Wolfe MS. The current status and prospects of multiline cultivars and variety mixtures for disease resistance. Annu. Rev. Phytopathol. 1985; 23(1):251–273. https://doi.org/10.1146/annurev.py.23.090185.001343

21. Fuchs M. Pyramiding resistance-conferring gene sequences in crops. Cur. Opin. Virol. 2017; 26:36–42. https://doi.org/10.1016/j.coviro.2017.07.004 PMID: 28755651

22. Ellis JG, Lagudah ES, Spielmeyer W, Dodds PN. The past, present and future of breeding rust resistant wheat. Front. Plant Sci. 2014; 5:641. https://doi.org/10.3389/fpls.2014.00641 PMID: 25505474

23. Djian-Caporalino C, Palloix A, Fazari A, Marteu N, Barbary A, Abad P, et al. Pyramiding, alternating or mixing: comparative performances of deployment strategies of nematode resistance genes to promote plant resistance efficiency and durability. BMC Plant Biol. 2014; 14(1):1–13. https://doi.org/10.1186/1471-2229-14-53 PMID: 24559060

24. Crete R, Renton M, editors. A new model to investigate whether regional crop rotation strategies can protect crops from fungal pathogens. 21st International Congress on Modelling and Simulation; 2015; Australia.

25. Milne AE, Bell JR, Hutchison WD, van den Bosch F, Mitchell PD, Crowder D, et al. The effect of farmers? Decisions on pest control with *Bt* Crops: A billion dollar game of strategy. PLoS Comput. Biol. 2016; 11(12):e1004483. https://doi.org/10.1371/journal.pcbi.1004483 PMID: 26720851
26. Mikaberidze A, McDonald BA, Bonhoeffer S. Developing smarter host mixtures to control plant disease. Plant Pathol. 2015; 64(4):996-1004. https://doi.org/10.1111/ppa.12321

27. Skelsey P, Rossing WAH, Kessel GJT, Powell J, van der Werf W. Influence of host diversity on development of epidemics: An evaluation and elaboration of mixture theory. Phytopathology. 2005; 95(4):328-338. https://doi.org/10.1094/PHYTO-95-0328 PMID: 18943033

28. Xu XM, Ridout MS. Stochastic simulation of the spread of race-specific and race-nonspecific aerial fungal pathogens in cultivar mixtures. Plant Pathol. 2000; 49(2):207-218. https://doi.org/10.1046/j.1365-3059.2000_t01-1-00444.x

29. Bourget R, Chaumont L, Sapoukhina N. Timing of pathogen adaptation to a multicomponent treatment. PLoS ONE. 2013; 8(8):e71926. https://doi.org/10.1371/journal.pone.0071926 PMID: 23991006

30. Fabre F, Bruchou C, Palloix A, Moury B. Key determinants of resistance durability to plant viruses: Insights from a model linking within- and between-host dynamics. Virus Res. 2009; 141(2):140-149. https://doi.org/10.1016/j.virusres.2008.11.021 PMID: 19159653

31. Alstad DN, Andow DA. Managing the evolution of insect resistance to transgenic plants. Science. 1995; 268(5219):1894-1896. https://doi.org/10.1126/science.268.5219.1894 PMID: 17797533

32. Fabre F, Rousseau E, Mailleret L, Moury B. Durable strategies to deploy plant resistance in agricultural landscapes. New Phytol. 2012; 193(4):1064-1075. https://doi.org/10.1111/j.1469-8137.2011.04019.x PMID: 22260272

33. Fabre F, Rousseau E, Mailleret L, Moury B. Epidemiological and evolutionary management of plant resistance: optimizing the deployment of cultivar mixtures in time and space in agricultural landscapes. Evol. Appl. 2015; 8(10):919-932. https://doi.org/10.1111/eva.12304 PMID: 26640518

34. Holt J, Chancellor TCB. Modelling the spatio-temporal deployment of resistant varieties to reduce the incidence of rice tungro disease in a dynamic cropping system. Plant Pathol. 1999; 48(4):453-461. https://doi.org/10.1046/j.1365-3059.1999.00360.x

35. Lo Iacono G, van den Bosch F, Gilligan CA. Durable resistance to crop pathogens: An epidemiological framework to predict risk under uncertainty. PLoS Comput. Biol. 2013; 9(1):e1002870. https://doi.org/10.1371/journal.pcbi.1002870 PMID: 23341765

36. Lo Iacono G, van den Bosch F, Paveley N. The evolution of plant pathogens in response to host resistance: Factors affecting the gain from deployment of qualitative and quantitative resistance. J. Theor. Biol. 2012; 304:152-163. https://doi.org/10.1016/j.jtbi.2012.03.033 PMID: 22483999

37. Pietravalle S, Lemarié S, van den Bosch F. Durability of resistance and cost of virulence. Eur. J. Plant Pathol. 2006; 114(1):107-116. https://doi.org/10.1007/s10658-005-3479-7

38. Storer NP. A spatially explicit model simulating western corn rootworm (Coleoptera: Chrysomelidae) adaptation to insect-resistant maize. J. Econ. Entomol. 2003; 96(5):1530-1547. https://doi.org/10.1093/jee/96.5.1530 PMID: 14650529

39. van den Bosch F, Gilligan CA. Measures of durability of resistance. Phytopathology. 2003; 93(5):616-625. https://doi.org/10.1094/PHYTO.2003.93.5.616 PMID: 18942985

40. Papaïx J, Adamczyk-Chauvat K, Bouvier A, Kièu K, S. T, Lannou C, et al. Pathogen population dynamics in agricultural landscapes: The *Ddal* modelling framework. Infect. Genet. Evol. 2014; 27:509-520. https://doi.org/10.1016/j.meegid.2014.01.022 PMID: 24480053

41. Papaix J, Burdon JJ, Lannou C, Thrall PH. Evolution of pathogen specialisation in a host metapopulation: joint effects of host and pathogen dispersal. PLoS Comput. Biol. 2014; 10(5):e1003633. https://doi.org/10.1371/journal.pcbi.1003633 PMID: 24853675

42. Papaix J, Touzeau S, Monod H, Lannou C. Can epidemic control be achieved by altering landscape connectivity in agricultural systems? Ecol. Model. 2014; 284:35-47. https://doi.org/10.1016/j.ecolmodel.2014.04.014

43. Ohtsuki A, Sasaki A. Epidemiology and disease-control under gene-for-gene plant-pathogen interaction. J. Theor. Biol. 2006; 238(4):780-794. https://doi.org/10.1016/j.jtbi.2005.06.030 PMID: 16085107

44. Vyska M, Cunniffe N, Gilligan C. Trade-off between disease resistance and crop yield: a landscape-scale mathematical modelling perspective. J. R. Soc. Interface. 2016; 13:20160451. https://doi.org/10.1098/rsif.2016.0451 PMID: 27707906

45. Papaix J, Rimbaud L, Burdon JJ, Zhan J, Thrall PH. Differential impact of landscape-scale strategies for crop cultivar deployment on disease dynamics, resistance durability and long-term evolutionary control. Evol. Appl. early view. https://doi.org/10.1111/eva.12570

46. Hobbelen PHF, Paveley ND, van den Bosch F. The emergence of resistance to fungicides. PLoS ONE. 2014; 9(3):e91910. https://doi.org/10.1371/journal.pone.0091910 PMID: 24658678

47. Kiyosawa S. Theoretical comparison between mixture and rotation cultivations of disease-resistant varieties. Ann. Phytopathol. Soc. Jpn. 1972; 38(1):52-59. https://doi.org/10.3186/jjphytopath.38.52
48. Sapoukhina N, Durel C-E, Le Cam B. Spatial deployment of gene-for-gene resistance governs evolution and spread of pathogen populations. Theor. Ecol. 2009; 2(4):229-238. https://doi.org/10.1007/s12080-009-0045-5

49. Skelsey P, Rossing WAH, Kessel GJT, van der Werf W. Invasion of *Phytophthora infestans* at the landscape level: How do spatial scale and weather modulate the consequences of spatial heterogeneity in host resistance? Phytopathology. 2010; 100(11):1146-1161. https://doi.org/10.1094/PHYTO-06-09-0148 PMID: 20932163

50. Djidjou-Demasse R, Moury B, Fabre F. Mosaics often outperform pyramids: insights from a model comparing strategies for the deployment of plant resistance genes against viruses in agricultural landscapes. New Phytol. 2017; 216(1):239-253. https://doi.org/10.1111/nph.14701 PMID: 28776688

51. Lof ME, de Vallavieille-Pope C, van der Werf W. Achieving durable resistance against plant diseases: scenario analyses with a national-scale spatially explicit model for a wind-dispersed plant pathogen. Phytopathology. 2017; 107(5):580-589. https://doi.org/10.1094/PHYTO-05-16-0207-R PMID: 28095206

52. Lof ME, van der Werf W. Modelling the effect of gene deployment strategies on durability of plant resistance under selection. Crop Prot. 2017; 97:10-17. https://doi.org/10.1016/j.cropro.2016.11.031

53. Consortium REX. Heterogeneity of selection and the evolution of resistance. Trends Ecol. Evol. 2013; 28(2):110-118. https://doi.org/10.1016/j.tree.2012.09.001 PMID: 23040463

54. Consortium REX. Combining selective pressures to enhance the durability of disease resistance genes. Front. Plant Sci. 2016; 7:1916. https://doi.org/10.3389/fpls.2016.01916 PMID: 28066472

55. Papaïx J, Burdon JJ, Zhan J, Thrall PH. Crop pathogen emergence and evolution in agro-ecological landscapes. Evol. Appl. 2015; 8(4):385-402. https://doi.org/10.1111/eva.12251 PMID: 25926883

56. Bourget R, Chaumont L, Durel C-E, Sapoukhina N. Sustainable deployment of QTLs conferring quantitative resistance to crops: first lessons from a stochastic model. New Phytol. 2015; 206(3):1163-1171. https://doi.org/10.1111/nph.13295 PMID: 25623549

57. Johnson R. A critical analysis of durable resistance. Annu. Rev. Phytopathol. 1984; 22(1):309-330. https://doi.org/10.1146/annurev.py.22.090184.001521

58. Burdon JJ, Zhan J, Barrett LG, Papaïx J, Thrall PH. Addressing the challenges of pathogen evolution on the world's arable crops. Phytopathology. 2016; 106(10):1117-1127. https://doi.org/10.1094/PHYTO-01-16-0036-FI PMID: 27584868

59. Chen XM. Epidemiology and control of stripe rust *Puccinia striiformis* f. sp. tritici on wheat. Can. J. Plant Pathol. 2005; 27(3):314-337. https://doi.org/10.1080/07060660509507230

60. Park RF. Breeding cereals for rust resistance in Australia. Plant Pathol. 2008; 57(4):591-602. https://doi.org/10.1111/j.1365-3059.2008.01836.x

61. Boyd LA. Can robusig defeat an old enemy?—Yellow rust of wheat. J. Agric. Sci. 2005; 143:233-243. https://doi.org/10.1017/S0021859605005095

62. Thompson JN, Burdon JJ. Gene-for-gene coevolution between plants and parasites. Nature. 1992; 360(6400):121-125. https://doi.org/10.1038/360121a0

63. Papaïx J, Goyeau H, Du Cheyron P, Monod H, Lannou C. Influence of cultivated landscape composition on variety resistance: an assessment based on wheat leaf rust epidemics. New Phytol. 2011; 191 (4):1095-1107. https://doi.org/10.1111/j.1469-8137.2011.03764.x PMID: 21585391

64. Leach JE, Vera Cruz CM, Bai J, Leung H. Pathogen fitness penalty as a predictor of durability of disease resistance genes. Annu. Rev. Phytopathol. 2001; 39(1):187-224. https://doi.org/10.1146/annurev.phyto.39.1.187 PMID: 11701864

65. Laine AL, Barrès B. Epidemiological and evolutionary consequences of life-history trade-offs in pathogens. Plant Pathol. 2013; 62:96-105. https://doi.org/10.1111/ppa.12129

66. Thrall PH, Burdon JJ. Evolution of virulence in a plant host-pathogen metapopulation. Science. 2003; 299(5613):1735-1737. https://doi.org/10.1126/science.1080070 PMID: 12637745

67. Kiêu K, Adamczyk-Chauvat K, Monod H, Stoica RS. A completely random T-tessellation model and Gibbsian extensions. Spat. Stat. 2013; 6:118-138. https://doi.org/10.1016/j.spasta.2013.09.003

68. Castellazzi MS, Matthews J, Angevin F, Sausse C, Wood GA, Burgess PJ, et al. Simulation scenarios of spatio-temporal arrangement of crops at the landscape scale. Environ. Model. Softw. 2010; 25 (12):1881-1889. https://doi.org/10.1016/j.envsoft.2010.04.006

69. Flor HH. Current status of the gene-for-gene concept. Annu. Rev. Phytopathol. 1971; 9(1):275-296. https://doi.org/10.1146/annurev.py.09.090171.001423

70. Débarre F, Gandon S. Evolution of specialization in a spatially continuous environment. J. Evol. Biol. 2010; 23(5):1090-1099. https://doi.org/10.1111/j.1420-9101.2010.01966.x PMID: 20345821
71. Brown JKM, Rant JC. Fitness costs and trade-offs of disease resistance and their consequences for breeding arable crops. Plant Pathol. 2013; 62:83–95. https://doi.org/10.1111/ppa.12163

72. Bouvier A, Kiêu K, Adamczyk K, Monod H. Computation of the integrated flow of particles between polygons. Environ. Model. Softw. 2009; 24(7):843–849. https://doi.org/10.1016/j.envsoft.2008.11.006

73. Waggoner PE, Berger RD. Defoliation, disease, and growth. Phytopathology. 1987; 77(3):393–398. url: https://www.apsnet.org/publications/.../phyto77n11_1495.pdf.

74. Jeger MJ. Analysis of disease progress as a basis for evaluating disease management practices. Annu. Rev. Phytopathol. 2004; 42(1):61–82. https://doi.org/10.1146/annurev.phyto.42.040803.140427 PMID: 15283660

75. Mersha Z, Zhang SA, Hau B. Effects of temperature, wetness duration and leaf age on incubation and latent periods of black leaf mold (*Pseudocercospora fuligena*) on fresh market tomatoes. Eur. J. Plant Pathol. 2014; 138(1):39–49. https://doi.org/10.1007/s10658-013-0295-3

76. Sartwell PE. The distribution of incubation periods of infectious disease. Am. J. Hyg. 1950; 51(3):310–318. url: https://academic.oup.com/aje/article-pdf/51/3/310/.../51-3-310.pdf. PMID: 15413610

77. Rimbaud L, Dallot S, Delaunay A, Borron S, Soubeyrand S, Thébaud G, et al. Assessing the mismatch between incubation and latent periods for vector-borne diseases: The case of sharka. Phytopathology. 2015; 105(11):1408–1416. https://doi.org/10.1094/PHYTO-01-15-0014-R PMID: 26512749

78. Soubeyrand S, Sache I, Lannou C, Chadœuf J. A frailty model to assess plant disease spread from individual count data. J. Data Sci. 2007; 5:63–83. url: http://www.jds-online.com/v5-1.

79. Sackett KE, Mundt CC. The effects of dispersal gradient and pathogen life cycle components on epidemic velocity in computer simulations. Phytopathology. 2005; 95(9):992–1000. https://doi.org/10.1094/PHYTO-95-0992 PMID: 18943297

80. Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford, UK: Oxford University Press; 1992. 757 p.

81. Mikaberidze A, Mundt CC, Bonhoeffer S. Invasiveness of plant pathogens depends on the spatial scale of host distribution. Ecol. Appl. 2016; 26(4):1238–1248. https://doi.org/10.1890/15-0807 PMID: 27509761

82. Heagle AS, Moore MB. Some effects of moderate adult resistance to crown rust of oats. Phytopathology. 1970; 60:461–466. https://doi.org/10.1094/Phyto-60-461

83. Sache I, de Vallavieille-Pope C. Classification of airborne plant pathogens based on sporulation and infection characteristics. Can. J. Bot. 1995; 73(8):1186–1195. https://doi.org/10.1139/b95-128

84. Katsuya K, Green GJ. Reproductive potentials of races 15B and 56 of wheat stem rust. Can. J. Bot. 1967; 45(7):1077–1091. https://doi.org/10.1139/b67-113

85. Mortensen K, Green GJ. Assessment of receptivity and urediospore production as components of wheat stem rust resistance. Can. J. Bot. 1978; 56(16):1827–1839. https://doi.org/10.1139/b78-221

86. de Vallavieille-Pope C, Giosue S, Munk L, Newton A, Niks R, Østergård H, et al. Assessment of epidemiological parameters and their use in epidemiological and forecasting models of cereal airborne diseases. Agronomie. 2000; 20(7):715–727. https://doi.org/10.1051/agro:2000171

87. Teng PS, Close RC. Effect of temperature and uredinium density on urediniospore production, latent period, and infectious period of *Puccinia hordei* Otth. N. Z. J. Agric. Res. 1978; 21(2):287–296. https://doi.org/10.1080/00288233.1978.10427413

88. Mehta YR, Zadoks JC. Uredospore production and sporulation period of *Puccinia recondita* f. sp. tritici on primary leaves of wheat. Neth. J. Plant Pathol. 1970; 76(5):267–276. https://doi.org/10.1007/bf03041356

89. Sache I, de Vallavieille-Pope C. Comparison of the wheat brown and yellow rusts for monocyclic sporulation and infection processes, and their polycyclic consequences. J. Phytopathol. 1993; 138(1):55–65. https://doi.org/10.1111/j.1439-0434.1993.tb01360.x

90. Sache I. Effect of density and age of lesions on sporulation capacity and infection efficiency in wheat leaf rust (*Puccinia recondita* f.sp. tritici). Plant Pathol. 1997; 46(4):581–589. https://doi.org/10.1046/j.1365-3059.1997.d01-33.x

91. Tomerlin JR, Eversmeyer MG, Kramer CL, Browder LE. Temperature and host effects on latent and infectious periods and on uredidniospore production of *Puccinia recondita* f. sp. tritici. Phytopathology. 1983; 73:414–419. https://doi.org/10.1094/Phyto-73-414

92. Richardson KL, Vales MI, Kling JG, Mundt CC, Hayes PM. Pyramiding and dissecting disease resistance QTL to barley stripe rust. Theor. Appl. Genet. 2006; 113(3):485–495. https://doi.org/10.1007/s00122-006-0314-2 PMID: 16736138

93. Sandoval-Islas JS, Broers LHM, Mora-Aguilera G, Parlevliet JE, Osada-Kawasoe S, Vivar HE. Quantitative resistance and its components in 16 barley cultivars to yellow rust, *Puccinia striiformis* f. sp. hordei. Euphytica. 2007; 153(3):295–308. https://doi.org/10.1007/s10681-006-9236-y
94. Broers LHM. Components of quantitative resistance to yellow rust in ten spring bread wheat cultivars and their relations with field assessments. Euphytica. 1997; 96(2):215-223. https://doi.org/10.1023/a:1002916110347

95. Cromey MG. Adult plant resistance to stripe rust (*Puccinia striiformis*) in some New Zealand wheat cultivars. N. Z. J. Crop Hortic. Sci. 1992; 20(4):413-419. https://doi.org/10.1080/01140671.1992.10418058

96. Elahinia SA, Tewari JP. Assessment of two different sources of durable resistance and susceptible cultivar of wheat to stripe rust (*Puccinia striiformis* f. sp. *tritici*). Caspian J. Environ. Sci. 2005; 3 (2):117-122. url: http://cjes.guilan.ac.ir/article_948_5e0230f27bfd99b837bd4d79b5fd8591.pdf.

97. Ma H, Singh RP. Expression of adult plant resistance to stripe rust at different growth stages of wheat. Plant Dis. 1996; 80:375-379. https://doi.org/10.1094/PD-80-0375

98. McGregor AJ, Manners JG. The effect of temperature and light intensity on growth and sporulation of *Puccinia striiformis* on wheat. Plant Pathol. 1985; 34(2):263-271. https://doi.org/10.1111/j.1365-3059.1985.tb01358.x

99. Milus EA, Seyran E, McNew R. Aggressiveness of Puccinia striiformis f. sp. tritici isolates in the south-central United States. Plant Dis. 2006; 90(7):847-852. https://doi.org/10.1094/PD-90-0847

100. Milus EA, Kristensen K, Hovmøller MS. Evidence for increased aggressiveness in a recent widespread strain of *Puccinia striiformis* f. sp. *tritici* causing stripe rust of wheat. Phytopathology. 2009; 99(1):89-94. https://doi.org/10.1094/PHYTO-99-1-0089 PMID: 19055439

101. Park RF, Rees RG. Expression of adult plant resistance and its effect on the development of *Puccinia striiformis* f.sp. *tritici* in some Australian wheat cultivars. Plant Pathol. 1989; 38(2):200-208. https://doi.org/10.1111/j.1365-3059.1989.tb02134.x

102. Quan W, Hou G, Chen J, Du Z, Lin F, Guo Y, et al. Mapping of QTL lengthening the latent period of *Puccinia striiformis* in winter wheat at the tillering growth stage. Eur. J. Plant Pathol. 2013; 136(4):715-727. https://doi.org/10.1007/s10658-013-0201-z

103. Sørensen CK, Hovmøller MS, Leconte M, Dedryver F, de Vallavieille-Pope C. New races of *Puccinia striiformis* found in Europe reveal race specificity of long-term effective adult plant resistance in wheat. Phytopathology. 2014; 104(10):1042-1051. https://doi.org/10.1094/PHYTO-12-13-0337-R PMID: 24624957

104. Denissen CJM. Components of adult plant resistance to leaf rust in wheat. Euphytica. 1993; 70 (1):131-140. https://doi.org/10.1007/bf00029650

105. Eyal Z, Peterson JL. Uredospore production of five races of Puccinia recondita rob. ex desm. as affected by light and temperature. Can. J. Bot. 1967; 45(4):537-540. https://doi.org/10.1139/b67-058

106. Lehman JS, Shaner G. Genetic variation in latent period among isolates of *Puccinia recondita* f.sp. *tritici* on partially resistant wheat cultivars. Phytopathology. 1998; 86(5):633-641. https://doi.org/10.1094/Phyto-86-633

107. Rimé D, Robert C, Goyeau H, Lannou C. Effect of host genotype on leaf rust (*Puccinia triticina*) lesion development and urediniospore production in wheat seedlings. Plant Pathol. 2005; 54(3):287-298. https://doi.org/10.1111/j.1365-3059.2005.01174.x

108. Robert C, Bancal M-O, Lannou C. Wheat leaf rust uredospore production and carbon and nitrogen export in relation to lesion size and density. Phytopathology. 2002; 92(7):762-768. https://doi.org/10.1094/PHYTO.2002.92.7.762 PMID: 18943273

109. Frantzen J, van den Bosch F. Spread of organisms: can travelling and dispersive waves be distinguished? Basic Appl. Ecol. 2000; 1(1):83-92. https://doi.org/10.1078/1439-1791-00010

110. Sackett KE, Mundt CC. Primary disease gradients of wheat stripe rust in large field plots. Phytopathology. 2005; 95(9):983-991. https://doi.org/10.1094/PHYTO-95-0983 PMID: 18943296

111. Soubeyrand S, Enjalbert J, Sanchez A, Sache I. Anisotropy, in density and in distance, of the dispersal of yellow rust of wheat: Experiments in large field plots and estimation. Phytopathology. 2007; 97 (10):1315-1324. https://doi.org/10.1094/PHYTO-97-10-1315 PMID: 18943690

112. Farber DH, Medlock J, Mundt CC. Local dispersal of *Puccinia striiformis* f. sp. *tritici* from isolated source lesions. Plant Pathol. 2016; 66(1):28-37. https://doi.org/10.1111/ppa.12554 PMID: 28082751

113. Fitt BDL, Gregory PH, Todd AD, McCartney HA, Macdonald OC. Spore dispersal and plant disease gradients: a comparison between two empirical models. J. Phytopathol. 1987; 118(3):227-242. https://doi.org/10.1111/j.1439-0434.1987.tb00452.x

114. Xu XM, Ridout MS. Effects of initial epidemic conditions, sporulation rate, and spore dispersal gradient on the spatio-temporal dynamics of plant disease epidemics. Phytopathology. 1998; 88(10):1000-1012. https://doi.org/10.1094/PHYTO.1998.88.10.1000 PMID: 18944811

115. Van de Wouw AP, Marcroft SJ, Howlett BJ. Blackleg disease of canola in Australia. Crop and Pasture Science. 2016; 67(3-4):273-283. https://doi.org/10.1071/CP15221
116. Moreno A, Fereres A, Cambra M. Quantitative estimation of plum pox virus targets acquired and transmitted by a single *Myzus persicae*. Arch. Virol. 2009; 154(9):1391-1399. https://doi.org/10.1007/s00705-009-0450-5 PMID: 19597934

117. Rimbaud L, Papaïx J, Rey J-F. landsepi: Landscape Epidemiology and Evolution. R package version 0.0.2. 2018. url: https://cran.r-project.org/package=landsepi.

118. R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2012. 3551 p.

119. Douglas N, Furrer R, Paige J, Sain S. fields: Tools for spatial data. R package version 8.10. 2015. https://doi.org/10.5065/D6W957CT

120. Bates D, Maechler M. Matrix: Sparse and dense matrix classes and methods. R package version 1.2-11. 2017. url: https://CRAN.R-project.org/package=Matrix.

121. Sánchez-Solana F, Ros C, Guerrero MdM, Martínez V, Lacasa CM, Hernández A, et al. Effectiveness of quantitative resistance conferred by the genetic background of pepper in the control of root-knot nematodes and influence onto durability of Me1- and Me3-resistant genes in greenhouse conditions. Plant Breed. 2017; 136(5):759-766. https://doi.org/10.1111/pbr.12517

122. Quenouille J, Montarry J, Palloix A, Moury B. Farther, slower, stronger: how the plant genetic background protects a major resistance gene from breakdown. Molec. Plant Pathol. 2013; 14(2):109-118. https://doi.org/10.1111/j.1364-3703.2012.00834.x PMID: 23046402

123. Brun H, Chèvre A-M, Fitt BDL, Powers S, Besnard A-L, Ermel M, et al. Quantitative resistance increases the durability of qualitative resistance to *Leptosphaeria maculans* in *Brassica napus*. New Phytol. 2010; 185(1):285-299. https://doi.org/10.1111/j.1469-8137.2009.03049.x PMID: 19814776

124. Palloix A, Ayme V, Moury B. Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies. New Phytol. 2009; 183(1):190-199. https://doi.org/10.1111/j.1469-8137.2009.02827.x PMID: 19344475

125. Gandon S, Day T. The evolutionary epidemiology of vaccination. J. R. Soc. Interface. 2007; 4 (16):803-817. https://doi.org/10.1098/rsif.2006.0207 PMID: 17264055

126. Gandon S, Mackinnon MJ, Nee S, Read AF. Imperfect vaccines and the evolution of pathogen virulence. Nature. 2001; 414(6865):751-756. https://doi.org/10.1038/414751a PMID: 11742400

127. Restif O. Evolutionary epidemiology 20 years on: Challenges and prospects. Infect. Genet. Evol. 2009; 9(1):108-123. https://doi.org/10.1016/j.meegid.2008.09.007 PMID: 18977460

128. Wingen LU, Shaw MW, Brown JKM. Long-distance dispersal and its influence on adaptation to host resistance in a heterogeneous landscape. Plant Pathol. 2013; 62(1):9-20. https://doi.org/10.1111/j.1365-3059.2012.02621.x

129. Suzuki SU, Sasaki A. How does the resistance threshold in spatially explicit epidemic dynamics depend on the basic reproductive ratio and spatial correlation of crop genotypes? J. Theor. Biol. 2011; 276(1):117-125. https://doi.org/10.1016/j.jtbi.2011.02.002 PMID: 21300071

130. Lô-Pelzer E, Bousset L, Jeuffroy MH, Salam MU, Pinochet X, Boillot M, et al. SIPPOM-WOSR: A Simulator for Integrated Pathogen Population Management of phoma stem canker on Winter OilSeed Rape: I. Description of the model. Field Crops Res. 2010; 118(1):73-81. https://doi.org/10.1016/j.fcr.2010.04.007

131. Kiyosawa S. Genetics and epidemiological modeling of breakdown of plant disease resistance. Annu. Rev. Phytopathol. 1982; 20(1):93-117. https://doi.org/10.1146/annurev.py.20.090182.000521

132. Lion S, Gandon S. Spatial evolutionary epidemiology of spreading epidemics. Proc. R. Soc. B. 2016; 283:20161170. https://doi.org/10.1098/rspb.2016.1170 PMID: 27798295

133. van den Bosch F, Gilligan CA. Models of fungicide resistance dynamics. Annu. Rev. Phytopathol. 2008; 46(1):123-147. https://doi.org/10.1146/annurev.phyto.011108.135838 PMID: 18680425

134. Montarry J, Cartier E, Jacquemond M, Palloix A, Moury B. Virus adaptation to quantitative plant resistance: erosion or breakdown? J. Evol. Biol. 2012; 25(11):2242-2252. https://doi.org/10.1111/j.1420-9101.2012.02600.x PMID: 22947055

135. Leonard KJ, Mundt CC. Methods for estimating epidemiological effects of quantitative resistance to plant diseases. Theor. Appl. Genet. 1984; 67(2):219-230. https://doi.org/10.1007/bf00317041 PMID: 24258552

136. St Clair DA. Quantitative disease resistance and quantitative resistance loci in breeding. Annu. Rev. Phytopathol. 2010; 48(1):247-268. https://doi.org/10.1146/annurev-phyto-080508-081904 PMID: 19400646

137. Gandon S, Michalakis Y. Evolution of parasite virulence against qualitative or quantitative host resistance. Proc. R. Soc. B. 2000; 267(1447):985-990. https://doi.org/10.1098/rspb.2000.1100 PMID: 10874747
138. van den Berg F, Lannou C, Gilligan CA, van den Bosch F. Quantitative resistance can lead to evolutionary changes in traits not targeted by the resistance QTLs. Evol. Appl. 2014; 7(3):370–380. https://doi.org/10.1111/eva.12130 PMID: 24665339

139. Steppan SJ, Phillips PC, Houle D. Comparative quantitative genetics: evolution of the G matrix. Trends Ecol. Evol. 2002; 17(7):320–327. https://doi.org/10.1016/S0169-5347(02)02505-3

140. Cunniffe NJ, Laranjeira FF, Neri FM, DeSimone RE, Gilligan CA. Cost-effective control of plant disease when epidemiological knowledge is incomplete: modelling bahia bark scaling of citrus. PLoS Comput. Biol. 2014; 10(8):e1003753. https://doi.org/10.1371/journal.pcbi.1003753 PMID: 25102099

141. Parry M, Gibson GJ, Parnell S, Gottwald TR, Irey MS, Gast TC, et al. Bayesian inference for an emerging arboreal epidemic in the presence of control. Proc. Natl. Acad. Sci. USA. 2014; 111(17):6258–6262. https://doi.org/10.1073/pnas.1310997111 PMID: 24711393

142. Vinatier F, Gosme M, Valentin-Morison M. Explaining host–parasitoid interactions at the landscape scale: a new approach for calibration and sensitivity analysis of complex spatio-temporal models. Landscape Ecol. 2013; 28(2):217–231. https://doi.org/10.1007/s10980-012-9822-4

143. Pilet-Nayel M-L, Moury B, Caffier V, Montarry J, Kerlan M-C, Fournet S, et al. Quantitative resistance to plant pathogens in pyramiding strategies for durable crop protection. Front. Plant Sci. 2017; 8:1838. https://doi.org/10.3389/fpls.2017.01838 PMID: 29163575

144. Laloi G, Vergne E, Durel CE, Le Cam B, Caffier V. Efficiency of pyramiding of three quantitative resistance loci to apple scab. Plant Pathol. 2017; 66(3):412–422. https://doi.org/10.1111/ppa.12581