Julie Quenouille

Estelle Paulhiac

Benoît Moury

Benoit Moury

Alain Palloix
email: alain.palloix@avignon.inra.fr

Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants

Keywords: resistance durability, quantitative trait locus, resistance breakdown, Capsicum spp., PVY, eukaryotic translation initiation factor 4E, potyviruses, resistance spectrum

Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants

INTRODUCTION

The use of resistant cultivars provides efficient, economical and environmentally friendly crop protection against plant pathogens, and many cultivars carrying major resistance genes successfully protected cultivated crops against their main pathogens. When available in genetic resources, major resistance genes are very attractive for breeding because of their simple inheritance that makes rapid and cheap their introgression through backcrosses to high-yielding but susceptible cultivars, and because they are expected to confer a nearly-complete resistance against the targeted pathogen. These resistances involve a gene-for-gene interaction and are usually race-specific. The widespread deployment of such resistance genes in many elite cultivars imposes a strong selection pressure on the pathogen population leading to the appearance and/or the increase in frequency of pathogen variants overcoming the resistance (or resistance breaking (RB) variants). This adaptation of the target pathogen makes the resistance ineffective, losing its economical interest. In order to preserve the durability of cultivar resistance, management strategies were developed to handle the adaptation of the pathogen populations and limit the appearance of RB variants. When several major resistance genes are available to control a given pathogen species, pyramiding of the distinct genes into a same cultivar, alternation of resistance genes across cultivation cycles or association of cultivars with distinct resistance genes during the same cultivation cycle were proposed [START_REF] Wolfe | The current status and prospects of multiline cultivars and variety mixtures for disease resistance[END_REF][START_REF] Mundt | Relevance of integrated disease management to resistance durability[END_REF][START_REF] Pink | Strategies using genes for non-durable disease resistance[END_REF][START_REF] Mcdonald | Pathogen population genetics, evolutionary potential, and durable resistance[END_REF]. More recently, another strategy was suggested from experimental data showing that the durability of a major resistance gene was dependent on the plant genetic background in which it was introgressed. In 3 distinct pathosystems (plant-virus/fungus/nematode) the emergence and/or frequency of the pathogen variants breaking down the major resistance gene remained low or did not occur in the cultivars which combined the major resistance gene in a quantitatively resistant genetic background, but increased when the resistance gene was introgressed into a plant genotype with a susceptible genetic background [START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF][START_REF] Brun | Quantitative resistance increases the durability of qualitative resistance to Leptosphaeria maculans in Brassica napus[END_REF][START_REF] Fournet | Selection of nematodes by resistant plants has implications for local adaptation and cross-virulence[END_REF]. These results raise the hypothesis that the selective value of allelic variation at a resistance locus depends on the genetic background. It addresses the question of how genetic background impacts the effect of mutations and evolution [START_REF] Chandler | Does your gene need a background check? How genetic background impacts the analysis of mutations, genes, and evolution[END_REF]. Considering plant breeding practices it suggests that combining qualitative and quantitative resistance increases the durability of plant protection. This also raises a paradox, since the increase in resistance efficiency is expected to increase the selection pressure on the pathogen population and thus the risk to select adapted pathogen genotypes [START_REF] Poland | Shades of gray: the world of quantitative disease resistance[END_REF]. However, none of the pathosystem studies cited above formally demonstrated that the increase in durability resulted from the action of quantitative resistance factors, since the durability of the major resistance was compared between pairs of genotypes with fully different genetic backgrounds. Consequently the decrease in resistance breakdown frequency might be due to any other genes from those genetic backgrounds. Further on, [START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF] showed that the breakdown frequency of the major resistance allele pvr2 3 controlling Potato virus Y (PVY) was correlated to the quantitative resistance level in a set of pepper genotypes carrying pvr2 3 and differing in their genetic background. This suggested that the breakdown frequency of the major allele and the quantitative resistance shared a common genetic basis. When very few major resistance genes are available, this would open the possibility to preserve their durability. Moreover, it would provide tools for breeders to construct resistant cultivars with enhanced durability. In our study, we (i) examined whether resistance breakdown is a heritable trait, (ii) described the genetic architecture underlying the genetic background which affect the breakdown frequency of a major resistance gene and (iii) discussed how breeders can exploit these results to preserve the durability of major resistance genes. Using the pepper (Capsicum annuum)/Potato virus Y (PVY) system, our strategy consisted in (i) detecting QTLs affecting the breakdown frequency of the pvr2 3 major resistance allele, since breakdown frequency is a key determinant of resistance durability and

MATERIALS AND METHODS

Mapping population

A segregating doubled-haploid (DH) population comprising 350 lines was obtained from the F 1 hybrid between two Capsicum annuum inbred lines: 'Perennial' carrying the PVY major resistance allele pvr2 3 in a partially resistant genetic background is an inbred line issued from an Indian landrace and 'Yolo Wonder' carrying the PVY susceptibility allele pvr2 + . All DH lines were genotyped with a tetra-primer ARMS-PCR (amplification refractory mutation system -polymerase chain reaction) marker which targets SNP (single nucleotide polymorphism) signatures differentiating pvr2 + from pvr2 3 [START_REF] Rubio | Functional markers for selection of potyvirus resistance alleles at the pvr2-eIF4E locus in pepper using tetra-primer ARMS-PCR[END_REF]. For phenotyping, we used the 153 DH lines carrying the pvr2 3 allele but segregating for the genetic background.

Virus material

For evaluation of RB frequencies, we used a PVY chimera between PVY clones SON41p and LYE84.2 [START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF]. This PVY chimera, named "CI chimera", differs from SON41p only by the CI cistron, which was replaced by that of LYE84.2. The rest of the genome and more particularly the region directly involved in pvr2 3 RB (the VPg cistron) are strictly identical between the two PVY clones. This "CI chimera" was chosen because of its higher ability to break the pvr2 3 resistance down compared to the parental clones [START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF]. However, the breakdown capacity of the "CI chimera" was shown to be highly correlated with that of SON41p when inoculated to a large range of pepper genotypes [START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF]. The "CI chimera" and SON41p were not infectious per se in plants carrying pvr2 3 , i.e. in these plants, the "CI chimera" showed very weak accumulation detected through RT-PCR but not ELISA and suggesting that the RB mutation occurred in the resistant plant [START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF] whereas mutants of the "CI chimera" possessing single non-synonymous substitutions in the VPg cistron could be detected at high concentration in systemically infected leaves. In addition, the CI chimera showed the same distributions of RB mutations when inoculated on plants carrying pvr2 3 when compared to SON41p. The level of quantitative resistance controlled by the plant genetic background was evaluated using a single-nucleotide mutant of the "CI chimera" carrying an amino acid substitution at position 119 (aspartic acid to asparagine) of the VPg (named "CI chimera VPg-N") which conferred the breakdown capacity towards the pvr2 3 resistance gene and was shown to be the most frequent in isolates from field as well as from experimental evolution experiments [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF]. This mutant enables the pathogen to overcome the resistance conferred by pvr2 3 and to reveal the quantitative resistance due to the genetic background.

Measure of resistance breakdown frequency

The breakdown frequency of pvr2 3 was tested after inoculation of the PVY "CI chimera" to the 153 DH lines carrying the pvr2 3 allele. Throughout the test, plants were grown in a climate-controlled room at 20-22 °C, 12h light /day. Inoculum was obtained after virus propagation in Nicotiana tabacum cv. Xanthi plants. Fifty grams of Xanthi leaves showing PVY symptoms were crushed in 200 mL of 0.03M phosphate buffer (pH 7.0) supplemented with 2% (w:v) diethyldithiocarbamate, 4.5 g of active charcoal and 4.5 g of Carborundum. Thirty pepper seedlings with two expanded cotyledons (2-3 weeks after sowing) per genotype were inoculated mechanically on their cotyledons. Thirty-eight days post-inoculation (dpi), plants were submitted to virus detection by DAS-ELISA as described by [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF]. In these conditions, in each plant infected systemically, the virus population was shown to be composed of one or several VPg mutants [START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF]. For each pepper line, the RB frequency of pvr2 3 was assessed by the ratio of the number of systemically infected plants over the total number of inoculated plants. Two independent tests of 30 plants per DH line were performed.

Measure of quantitative resistance

To evaluate the level of the quantitative resistance due to the genetic background, we used the "CI chimera VPg-N" PVY mutant which overcomes the resistance conferred by pvr2 3 (infects all the genotypes carrying pvr2 3 ) and reveals the quantitative resistance due to the genetic background. Two different traits were assessed : the area under the disease progress curve (AUDPC) and the virus accumulation. The PVY "CI chimera VPg-N", was mechanically inoculated on 20 plants of each DH lines carrying pvr2 3 . The plants were grown in a climate-controlled room as previously described. The AUDPC assessment combined both intensity of symptoms and latency period: symptom intensity was evaluated, on the 20 plants, in a 3-class scale (1= no symptom, 2 = weak mosaic and 3 = severe mosaic/ necrosis) at 14, 21, 28 and 35 dpi. AUDPC was calculated using the formula:

where "n" is the number of symptom assessments, "y" the symptom intensity and "t" the time in dpi. At 36 dpi, the virus accumulation was independently evaluated on 10 individual plants per DH line by semi-quantitative DAS-ELISA, using a dilution range of extracts of PVY infected plants, and expressed relatively to a common reference sample incorporated in each ELISA plate, as in [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF].

Genotyping of molecular markers and linkage map

The 350 DH progeny was an extension of the former 'PY' mapping population [START_REF] Lefebvre | Towards the saturation of the pepper linkage map by alignment of three intraspecific maps including known-function genes[END_REF]. It was genotyped with 236 molecular markers: 108 single-nucleotide polymorphisms (SNPs, [START_REF] Nicolaï | Discovery of a large set of SNP and SSR genetic markers by high-throughput sequencing of pepper (Capsicum annuum)[END_REF], 29 SNPs [START_REF] Jung | Discovery of single nucleotide polymorphism in Capsicum and SNP markers for cultivar identification[END_REF], 47 simple sequence repeats (SSRs, [START_REF] Lee | Characterization and molecular genetic mapping of microsatellite loci in pepper[END_REF][START_REF] Nagy | Development, characterization, and transferability to other Solanaceae of microsatellite markers in pepper (Capsicum annuum L.)[END_REF][START_REF] Minamiyama | An SSR-based linkage map of Capsicum annuum[END_REF], 33 amplified fragment-length polymorphisms (AFLPs, [START_REF] Lefebvre | Towards the saturation of the pepper linkage map by alignment of three intraspecific maps including known-function genes[END_REF], one single-strand conformation polymorphism (SSCP-PM5, [START_REF] Fazari | The root-knot nematode resistance N-gene co-localizes in the Me-genes cluster on the pepper (Capsicum annuum L.) P9 chromosome[END_REF], two known genes and 16 new markers (supplementary table I) corresponding to previously mapped restriction fragment length polymorphisms (RFLPs) and/or conserved ortholog set (COS II) sequences and designed using the Sol Genomics network (www.solgenomics.net). The two known genes were pvr2 (the major resistance locus to PVY) genotyped using the tetra-ARMS marker as in [START_REF] Rubio | Functional markers for selection of potyvirus resistance alleles at the pvr2-eIF4E locus in pepper using tetra-primer ARMS-PCR[END_REF], and pvr6 coding for the eIF(iso)4E (eukaryotic initiation factor iso4E) protein genotyped with a new marker targeting a 82 nucleotide deletion differentiating alleles pvr6 + and pvr6 [START_REF] Ruffel | Simultaneous mutations in translation initiation factors eIF4E and eIF(iso)4E are required to prevent pepper veinal mottle virus infection of pepper[END_REF] (supplementary data I). The DH progeny was also genotyped for two phenotypic markers, the presence/absence of pungency in the fruits (locus C) and the erected/pendant fruit position (locus up).

The genetic linkage map was constructed using the Mapmaker software version 3.0b [START_REF] Lander | MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations[END_REF] with a LOD score threshold of 3.0, a maximum recombination fraction of 0.3 and distances between markers were calculated using the Kosambi mapping function. For each linkage group, marker order was verified with the "ripple" command and markers were conserved only if the LOD value was greater than 2.0. Three out of the 236 markers remained unlinked and 33 were deleted because they had no unique placement at LOD 2.0. The resulting map contained 200 markers distributed in 13 linkage groups (LG), which were assigned to the 12 pepper chromosomes through alignment of markers with previous pepper maps [START_REF] Lefebvre | Towards the saturation of the pepper linkage map by alignment of three intraspecific maps including known-function genes[END_REF][START_REF] Wu | A COSII genetic map of the pepper genome provides a detailed picture of synteny with tomato and new insights into recent chromosome evolution in the genus Capsicum[END_REF] Comment citer ce document : Quenouille, J., Paulhiac, E., Moury ). The total length of the map was 2457.7 cM with an average length interval between markers of 12.3 ± 12.4 cM. The increased length of this map compared to former pepper maps was mainly attributed to the larger size of the progeny (350 individuals) that allowed to declare significant a recombination rate up to 0.45 cM (slightly more than 50 Kosambi cM), and also to the assignment of the two linkage groups to the single P7 chromosome, resulting in longer chromosomes (supplementary data II).

Statistical analyses

Descriptive statistics were performed with the R software (http://www.r-project.org/). For each phenotypic trait, narrow-sense heritability (h²) was estimated using the formula h²=σ² G /(σ² G +σ² E /n), where σ² G corresponds to the genotypic variance, σ² E to the environment (error) variance and n to the number of replicates.

For AUDPC, mean distributions were close to a normal distribution and we analyzed nontransformed data for QTL detection. Means of virus accumulation (VA) were non-normally distributed. Thus, for QTL detection, these data were log transformed [ln (VA) +1] to approximate a normal distribution. Resistance breakdown frequency was a non-normally distributed trait with an excess of low values, and no transformation was efficient to improve the distribution.

QTL analysis by regression interval mapping

For a non-normal distribution, QTL detection can be performed with a standard regression-based interval mapping approach, which is robust against deviations from normality when the threshold value is calculated by permutations (Rebaï, 1997). QTL analyses were performed using iterative composite interval mapping (iQTLm) methods implemented in the MCQTL v.5.2.4 [START_REF] Jourjon | MCQTL: multi-allelic QTL mapping in multi-cross design[END_REF] software. MCQTL performed QTL mapping using a linear regression model and the iQTLm scan method allows a multiple QTL model to be determined automatically [START_REF] Charcosset | Heterosis in maize investigated using connected RIL populations[END_REF]. The QTL significance was determined using the MCQTL test, which is equal to -log 10 of the Fisher test Pvalue as described in the MCQTL v.5 reference manual (http://carlit.toulouse.inra.fr/MCQTL/). The significance threshold of the MCQTL test was estimated after performing 1000 permutation tests for each trait to reach a genome-wide type-I error rate of 0.05 and 0.1 for QTLs and cofactors, respectively. This resulted in a threshold of 3.4 for [log (VA) +1] and of 3.5 for RB and AUDPC. The cofactor threshold was fixed at 3.2 for the three traits.

Cofactor selection was performed by forward selection of markers. The QTL confidence intervals (CI) were defined as a 2 unit drop-off of the MCQTL test around the maximum MCQTL score to obtain approximately a 95% confidence interval : CI = [max (MCQTL test) -2; max (MCQTL test) + 2]. Epistasis between QTLs and the genetic background was tested using the iterative method implemented in MCQTL (threshold of 3.5). The part of phenotypic variation explained by additive and/or epistasic effects of each QTL (R² coefficients) and the global R² were estimated using the full QTL model in MCQTL.

QTL analysis by maximum likelihood interval mapping

For AUDPC and virus accumulation (log (VA) +1), QTL analyses were also performed using multiple QTL mapping (MQM) methods implemented in the R/qtl package of the R [START_REF] Broman | R/qtl: QTL mapping in experimental crosses[END_REF] software. We did not use this method for the RB trait because data were non-normally distributed. MQM is a maximum likelihood-based interval mapping method composed of three steps. First, missing data are "augmented", i.e. all genotype possibilities are modeled and their probability calculated from neighboring marker information. In the second step, significant markers are detected with the "mqmscan" function without cofactor, which corresponds to a simple interval with a backward elimination. In the last step, QTLs are detected using markers detected in step two as cofactors, which corresponds to a composite interval mapping, and a backward elimination is applied to remove the non-significant cofactors and to establish the final QTL model. Significance LOD thresholds were estimated for each trait, with a genome-wide type-I error rate of 0.05, after performing 1000 permutation tests. The LOD threshold was estimated at 2.6 for the two traits. For each QTL, the CI was defined as a 2-LOD drop-off around the maximum LOD score to obtain approximately a 95% confidence interval. R² coefficients were calculated with the "fitqtl" function of R/qtl.

RESULTS

Characterization of the phenotypic traits

To evaluate the frequency of pvr2 3 resistance breakdown, the "CI Chimera", which is not infectious per se toward pvr2 3 resistant plants, was inoculated on the set of 153 DH carrying pvr2 3 . In these plants, only RB mutants of the "CI Chimera" possessing single non-synonymous substitutions in the factor (VPg) conferring virulence towards the pvr2 3 allele were detected in previous studies [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF][START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF]. In our tests, mosaic symptoms appeared in 100% of the susceptible (pvr2 + ) control Yolo Wonder 2 weeks after inoculation, but mosaic or necrotic symptoms appeared later (4 to 5 weeks after inoculation) in some of the plants of some DH lines. Thirty-eight days after inoculation (dpi), the ELISA revealed an important phenotypic variability between the 153 DH lines, with frequencies of systemic infection, corresponding to RB frequencies, varying from 0% to 93.2% with a mean (± standard deviation) of 14.7% (± 23.8). The correlation between the two independent tests was high (ρ Spearman =0.75; p<10 -7 , supplementary data 3), and a heritability of 0.87 was estimated. The distribution of the trait was strongly skewed towards low values, with 70 DH lines displaying no resistance breakdown and 83 HD lines displaying a RB frequency between 1.7% and 93.2% (figure 1). The quantitative resistance was assessed using a pvr2 3 -breaking mutant of the "CI Chimera", the "CI Chimera VPg-N". With this virus, the pvr2 3 resistant lines displayed symptoms as soon as 14 dpi and 100% of the plants were infected at 35 dpi as attested by DAS-ELISA. Highly significant variations for virus accumulation (VA) and AUDPC were observed between DH lines with a heritability of 0.64 for Ln(VA+1) and 0.98 for AUDPC. For statistical and QTL analyses, VA values were log transformed [ln(VA+1)] to approximate a normal distribution. The ln(VA+1) varied from 0 to 1.43 with a mean equal to 0.67±0.35 and AUDPC varied from 29.75 to 63 with a mean equal to 46.6±7.9 (figure 1). Significant correlations (p<0.05) were observed between the three traits with Spearman's coefficients equal to 0.40, 0.33 and 0.32 for RB/VA, RB/AUDPC and VA/AUDPC, respectively.

Mapping QTLs for the quantitative resistance traits

Distributions of the quantitative resistance traits AUDPC and log(VA+1) were close to normality. For these two traits, QTL detection was performed using iQTLm and MQM methods. Two significant QTLs for VA and three for AUDPC were detected by both iQTLm and MQM methods. These QTLs were named VA-3, VA-6, A-1, A-3 and A-9, according to the trait considered and the chromosome location. Hereafter, we will detail only the results from iQTLm because results from both methods are similar. The position, significance and effect of each QTL are detailed in table 1 and figure 2.

VA-3 and VA-6 explained 34.5% and 15.7% of the VA variation, respectively. Epistasis tests did not reveal any interaction between QTLs and the genetic background. The final model combining the effects of the two significant QTLs explained 43.9% of VA phenotypic variation corresponding to 69% of the trait heritability (h²=0.64). A-1, A-3 and A-9 explained from 14.69%, 15.5% and 16.1% of the AUDPC variation. Epistasis tests did not reveal any interaction between these QTLs and the genetic background. The part of the AUDPC phenotypic variation explained by A-1, A-3 and A-9 was equal to 34.1% corresponding to 35% of the trait heritability and a large part of the variability remained unexplained (h²=0.98).

Mapping QTLs for the frequency of pvr2 3 -resistance breakdown

Because of the skewed distribution of the RB frequency, QTL detection was performed with the standard regression-based interval mapping approach [START_REF] Haley | A simple regression method for mapping quantitative trait loci in line crosses using flanking markers[END_REF]Rebaï, 1997) using the iQTLm [START_REF] Charcosset | Heterosis in maize investigated using connected RIL populations[END_REF] method implemented in MCQTL software [START_REF] Jourjon | MCQTL: multi-allelic QTL mapping in multi-cross design[END_REF]. Three significant QTLs, explaining 12.8%, 39.8% and 8.6% of the pvr2 3 RB frequency variation, were detected on chromosomes 1, 3 and 5, respectively. QTLs were named RB-1, RB-3 2 detail position, significance and effect of each QTL. Epistasis tests revealed that RB-3 interacts with three distinct positions in the pepper genome, on chromosomes 1, 5 and 6. The positions on chromosomes 1 and 5 interacting with RB-3 are included in the CIs of RB-1 and RB-5, respectively, suggesting interaction effects between RB-1 and RB-3 and between RB-3 and RB-5 (figure 2). The position on chromosome 6 that interacts with RB-3 did not show any additive effect and will be further named RB-6. The final model combining additive and epistasic effects of the significant QTLs explained 68.9% of the variation of pvr2 3 RB frequency, which corresponds to a great part (79%) of the trait heritability (h²=0.87). Because of the skewed segregation of the RB trait, the QTL effects (as well as heritability) may have been overestimated. Indeed, the residues of the final model did not fit a normal distribution (Shapiro-Wilk normality test W=0.87, p<10 -7 ). The analysis was reproduced after removing the 70 lines with a RB value of 0. The same additive QTLs were detected and the final model including these 3 additive QTLs explained 48% of the phenotypic variance (instead of 52% with the whole population) with a normal distribution of residues (W=0.98, p=0.55). The interaction effects were not significant anymore. This very probably resulted from the fact that removing the 70 individuals with RB=0 also removed most of the individuals carrying the Yolo Wonder allele at RB-3, resulting in highly unbalanced allele frequencies at this locus (71/12), reduced population size (n=83) and decreased statistical power. Looking for co-localization between QTLs, the confidence interval of A-1, proved fully included in the confidence interval of RB-1, as was the confidence interval of RB-3 within those of A-3 and VA-3. The major QTL affecting the RB frequency of pvr2 3 (RB-3) co-localizes with the major QTLs VA-3 and A-3 affecting VA and AUDPC, respectively. The QTL RB-1 also co-localizes with QTL A-1. The RB-6 locus displaying interaction effect with RB-3 also located at the LOD peak of VA-6, indicating that RB-6 colocalizes with VA-6 (figure 2).

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury

Effect of parental alleles on the RB frequency

All genotypes carrying the Yolo Wonder allele at QTL RB-3 showed a very low resistance breakdown frequency with a mean equal to 0.31% whereas, genotype carrying the Perennial allele at QTL RB-3 showed a much higher resistance breakdown frequency with a mean equal to 29.63%. The phenotypic values of the different allele combination between QTL RB-3 and epistatic QTLs are presented in Table 2. For genotypes carrying the Yolo Wonder allele at QTL RB-3, the resistance breakdown frequency was independent of the alleles at the interacting QTLs (mean of resistance breakdown frequency varies between 0.17% and 0.37%), so that these QTL effects were not significant in these genotypes. Contrarily, in genotypes carrying the Perennial allele at QTL RB-3, the RB frequency strongly depended on the alleles present at the interacting QTLs. The Yolo Wonder alleles at QTLs RB-1, RB-5 or RB-6 increased the RB frequency when combined with the Perennial allele at QTL RB-3, whereas the Perennial alleles at these QTLs strongly decreased the RB frequency (Table 2). Altogether, the combination of the Perennial alleles at the three epistatic QTLs (RB-1, RB-5 and RB-6) with the Perennial allele at RB-3 decreased the RB frequency from 69.7% to 5.3% (table 2). The high effect of these QTLs in genotypes carrying the Perennial allele at RB3, contrarily to those carrying the Yolo Wonder allele at this QTL probably caused the significant interaction (scale effect).

DISCUSSION

Durable management of plant resistance to pathogen populations was mostly promoted through the spatio-temporal management of resistance genes, but can also be achieved through management of resistance genes at the genotype scale: the cultivar genotype construction. In this paper we firstly succeeded to identify plant QTLs for a trait directly related to the evolution of the viral population: the frequency of emergence of viral variants adapted to the resistance. The colocation of the plant genetic factors modulating this viral evolution with the QTLs decreasing the virus accumulation further indicated that the more efficient the more durable was the resistance, which was not obvious since higher resistance levels are also expected to increase the selection pressure. Consequently, the high durability of the resistance of the Perennial pepper landrace was shown to result from a complex genetic architecture, including alleles at a major gene and at QTLs from the genetic background with complementary but also contradictory effects, which address the question of the selection processes leading to this combination in this landrace.

Heritability of the breakdown frequency of a major resistance gene

In this study, we measured the resistance breakdown (RB) frequency of a major resistance gene (pvr2 3 ) in a set of pepper DH lines carrying pvr2 3 but segregating for the genetic background. The observed RB frequency directly resulted from the frequency of appearance and from the [START_REF] Fabre | Key determinants of resistance durability to plant viruses: Insights from a model linking within-and between-host dynamics[END_REF]. Previous studies of this pathosystem showed that the RB mutation most probably occur in the resistant pepper host which allowed for weak PVY multiplication and that the genetic background acts at several levels, including the mutational pathways of the virus to overcome resistance and the speed of selection of the RB variants [START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF][START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF]. The genetic analysis showed that the RB frequency of pvr2 3 was highly variable between genotypes (from 0% to 93.2%) and highly heritable (h²=0.87). Because of the non-normal distribution of data for RB frequency, the h² value is most probably overestimated. However, it fits the Spearman correlation value between the 2 repeats of the test (0.75), demonstrating the genetic control of this trait. This result is not trivial. Indeed, during the RB processes, several factors could have increased the heterogeneity of results between plants and experiments and consequently could have decreased the heritability of RB frequency. Among these factors are the stochastic nature of (i) the appearance of RB mutations in PVY genome, (ii) the genetic drift acting on virus during systemic invasion of plants (Garcia-Arenal et al., 2001;[START_REF] French | Evolution of wheat streak mosaic virus: dynamics of population growth within plants may explain limited variation[END_REF] and (iii) the fact that at least 9 different VPg mutations, with contrasted effects on PVY competitiveness, can be involved in pvr2 3 breakdown [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Montarry | Genetic background matters: a plant-virus gene-for-gene interaction is strongly influenced by genetic contexts[END_REF]. Our study demonstrated that measuring the resistance breakdown frequency with a reasonable number of plants (60 per genotype) compensates these factors of heterogeneity and allows the genetic analyses. The genetic background dependency of pvr2 3 durability could be attributed to 4 distinct genetic factors with quantitative effects: one major-QTL (RB-3) and three additional QTLs (RB-1, RB-5 and RB-6) acting additively and/or in interaction with RB-3. The discovery of such genetic factors shows that direct selection for alleles increasing the durability of a major resistance is feasible and opens new tools for durable resistance breeding. Nevertheless, the successful detection of such genetic factors affecting the RB frequency of a major resistance gene depends on the ability to generate RB strains at high frequencies through experimental evolution. For many virus diseases, the experimental conversion of an avirulent virus to an RB variant is rather easy and was commonly used to determine viral determinants of pathogenicity [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Margaria | Evidence that the nonstructural protein of Tomato spotted wilt virus is the avirulence determinant in the interaction with resistant pepper carrying the Tsw gene[END_REF][START_REF] Hajimorad | Experimental adaptation of an rna virus mimics natural evolution[END_REF][START_REF] Lacroix | Effect of passage of a Potato virus Y isolate on a line of tobacco containing the recessive resistance gene va 2 on the development of isolates capable of overcoming alleles 0 and 2[END_REF]. Moreover, the RB frequency of resistance genes is highly correlated with the durability of these genes in the field [START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Janzac | Phenotype and spectrum of action of the Pvr4 resistance in pepper against potyviruses, and selection for virulent variants[END_REF][START_REF] Lacroix | Effect of passage of a Potato virus Y isolate on a line of tobacco containing the recessive resistance gene va 2 on the development of isolates capable of overcoming alleles 0 and 2[END_REF]. In tobacco, field observations revealing that va 2 -RB isolates were more frequent than va 0 -RB isolates [START_REF] Verrier | The ''va'' resistance to PVYN in Nicotiana tabacum: an assessment of the frequency of ''va''breaking PVYN strains based on seven years of field survey on a worldwide basis[END_REF] were confirmed by results from experimental evolution [START_REF] Lacroix | Effect of passage of a Potato virus Y isolate on a line of tobacco containing the recessive resistance gene va 2 on the development of isolates capable of overcoming alleles 0 and 2[END_REF]. Moreover, the distributions of RB mutations in viral genomes are usually similar in field and laboratory observations [START_REF] Moury | Mutations in Potato virus Y Genome-linked protein determine virulence toward recessive resistances in Capsicum annuum and Lycopersicon hirsutum[END_REF][START_REF] Ayme | Different mutations in the genome-linked protein vpg of potato virus y confer virulence on the pvr2 3 resistance in pepper[END_REF][START_REF] Hajimorad | Experimental adaptation of an rna virus mimics natural evolution[END_REF]. However, for some resistance genes which proved very durable in the field, it may be difficult to generate RB variants by direct inoculation of avirulent strains. In such cases, alternative protocols imposing strong inoculum pressures like graft-inoculation or inoculation through agroinfiltration can generate sufficient inoculum pressure to increase RB frequency [START_REF] Bruun-Rasmussen | The same allele of translation initiation factor 4E mediates resistance against two Potyvirus spp. in Pisum sativum[END_REF][START_REF] Janzac | Phenotype and spectrum of action of the Pvr4 resistance in pepper against potyviruses, and selection for virulent variants[END_REF]. These experimental procedures make it possible to measure the breakdown frequency of major resistance genes in many plant-virus pathosystems, providing powerful tools for improving the durability of virus resistance in plants. Considering other pathogens, experimental evolution may not be so relevant, due to their longer generation time, smaller populations and smaller mutation rates. However, if confirmed as a generic result, the relationship between RB QTLs and QTLs controlling the level of quantitative resistance may provide alternative selection criteria.

Alleles increasing quantitative resistance are favourable to the durability of major resistance genes.

Using the same pepper DH progeny, we performed a QTL analysis of quantitative resistance. To get rid of the pvr2 3 major effect and to reveal the quantitative resistance conferred by the genetic background, we used a pvr2 3 -RB PVY variant (the "CI Chimera VPg-N") which differs from the PVY variant used for the measurement of pvr2 3 RB frequency by a single nucleotide substitution. This mutation restores the interaction between the PVY VPg and the pepper eukaryotic translation initiation factor 4E1 (eIF4E1) encoded by pvr2 3 , causing the breakdown of the pvr2 3 -mediated resistance [START_REF] Charron | Natural variation and functionnal analyes provide evidence for co-evolution between plant eIF4E and potyviral VPg[END_REF]. Using this PVY variant, four QTLs affecting virus accumulation and/or symptom expression were detected. These four QTLs were also previously detected for the AUDPC trait with two other PVY strains, with identical allele directed effects: resistancefavourable Perennial alleles on chromosomes P1 and P6 (LG3 and LG10 in Caranta et al., 1997) and favourable Yolo Wonder alleles on chromosomes P9 and P3 (LG9 and LG4 in Caranta et al., 1997). In this experiment, distinct natural PVY isolates with distinct pvr2 3 breaking down mutations in their VPg were used, indicating no differences between the QTLs detected using distinct VPg mutants. Our present results show that three of the four QTLs affecting the RB frequency of pvr2 3 colocate with QTLs affecting quantitative resistance (viral accumulation and/or AUDPC). These colocations may result from genetic linkage between the different traits or from a pleiotropic effect of quantitative resistance factors on the RB frequency of pvr2 3 . Looking at the effect of parental alleles that colocate on chromosomes P1, P3 and P6, all the alleles that increased the quantitative resistance also decreased the breakdown frequency. Moreover, the QTL on chromosome P3 displayed the major effect on both virus accumulation and RB frequency. This concordance between parental allele effects for the different traits strongly suggests that colocations are due to pleiotropic effects and that alleles increasing the resistance level also increase the resistance durability. This is also consistent with the initial observations, in 3 pathosystems, that the durability of major genes is enhanced in partially resistant cultivars [START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF][START_REF] Brun | Quantitative resistance increases the durability of qualitative resistance to Leptosphaeria maculans in Brassica napus[END_REF]Fournet et al., 2012). However, it may look contradictory since higher resistance levels are expected to increase the selection pressure on the pathogen and the risk of selection of adapted pathogen populations [START_REF] Poland | Shades of gray: the world of quantitative disease resistance[END_REF]. In a recent study of the different mechanisms that could account for the higher durability of polygenic resistance, [START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF] showed the selection of adapted RB mutants by the polygenic resistance of Perennial was slower than by the monogenic resistance (pvr2 3 in a susceptible genetic background). In this same report, the additional quantitative resistance carried by the genetic background explained a large part of the variation in RB frequency. Dissecting the genetic background into QTLs indicate that breeding for quantitative resistance alleles will enhance the durability of major genes. Such a strategy is achievable for most pathosystems where breaking down variants or experimental evolution tests are not available.

Resistance durability as the result of a combination of mechanisms

At the landscape scale, it was proposed previously that diversifying the methods for pathogen control is an efficient strategy to limit adaptation of pathogen populations [START_REF] Mundt | Relevance of integrated disease management to resistance durability[END_REF][START_REF] Pink | Strategies using genes for non-durable disease resistance[END_REF]. Diversifying resistance mechanisms at the genotype scale can also have a positive effect on durability. In accordance with [START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF], our results raise the main impact of viral accumulation on the RB frequency of pvr2 3 , since all QTLs involved in viral accumulation are also involved in the RB frequency of pvr2 3 . Viral accumulation was measured in apical leaves and a decrease of viral accumulation can result from decreases of the efficiency of virus replication, cell-to-cell and/or systemic movement. Combining a major resistance gene with another gene involved in viral accumulation has already proved to enhance resistance durability in the PVY/tobacco system [START_REF] Acosta-Leal | Complementary functions of two recessive R-genes determine resistance durability of tobacco `Virgin A Mutant' (VAM) to Potato virus Y[END_REF]. In their study, [START_REF] Acosta-Leal | Complementary functions of two recessive R-genes determine resistance durability of tobacco `Virgin A Mutant' (VAM) to Potato virus Y[END_REF] showed that the RB frequency of the va recessive gene which confers a resistance to PVY cell-tocell movement in tobacco was reduced when va was combined with a second recessive gene (va2) conferring resistance to within-cell PVY accumulation. Two additional QTLs affected the RB frequency but had no detectable effect on viral accumulation (RB-1 and RB-5), suggesting distinct modes of action. QTL RB-1 colocates with QTL A-1 controlling symptom intensity and dynamics. Here we consider tolerance in the sense of [START_REF] Schafer | Tolerance to plant disease[END_REF] in which tolerant plants endure high pathogen loads without severe losses in field. Thus, QTL A-1 can be considered as a tolerance QTL since it reduces the AUDPC but it does not significantly affect viral accumulation. To our knowledge, mechanisms involved in plant tolerance are still unknown but they could also have an effect on the pvr2 3 RB frequency. Finally, QTL RB-5 did not colocate with any quantitative resistance QTL, suggesting that a third mechanism, independent of quantitative resistance, affects the pvr2 3 RB frequency. [START_REF] Quenouille | Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants[END_REF] also showed that beside virus accumulation, the plant genetic background can control two additional modes of action: making the mutational pathways for viral resistance breakdown more complex and unlikely and modulating the selection pressure exerted on viral populations. QTLs RB-1 and RB-5 may be candidates for such modes of actions.

A genetic trade-off between broad-spectrum resistance and resistance durability that can be solved by breeding strategy Among the 4 additive and/or epistatic QTLs affecting the RB frequency, the alleles from Perennial were shown to decrease the risk of pvr2 3 breakdown at 3 of these QTLs (RB-1, RB-5 and RB-6), whereas the Perennial allele increased strongly this risk for the major QTL RB-3. This result looks surprising since the pvr2 3 allele was previously shown to be highly durable in the Perennial landrace [START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF]. It is noteworthy that the peak of QTL RB-3 was localized at the pvr6 locus which encodes the eIF(iso)4E, an isoform of eIF4E1 encoded by pvr2. It was previously shown that the Perennial allele at pvr6 included a deletion of 82 nucleotides followed by a premature stop codon, compared to the pvr6 + allele from Yolo Wonder. The eIF(iso)4E protein encoded by the Perennial pvr6 allele was therefore non-functional [START_REF] Ruffel | Simultaneous mutations in translation initiation factors eIF4E and eIF(iso)4E are required to prevent pepper veinal mottle virus infection of pepper[END_REF] and can be considered as a deleted allele of pvr6 + . When combined with different pvr2 alleles, the Perennial pvr6 allele was shown to enlarge the spectrum of resistance to additional potyviruses, including Chilli veinal mottle virus (ChiVMV) and Pepper veinal mottle virus (PepVMV) [START_REF] Caranta | Polygenic resistance of pepper to potyviruses consists of a combination of isolate-specific and broad-spectrum quantitative trait loci[END_REF][START_REF] Moury | Serological, molecular, and pathotype diversity of Pepper veinal mottle virus and Chili veinal mottle virus[END_REF][START_REF] Rubio | Allele mining in the pepper gene pool provided new complementation effects between pvr2-eIF4E and pvr6-eIF(iso)4E alleles for resistance to pepper veinal mottle virus[END_REF]. Hence, the Perennial allele at the pvr6/RB-3 locus contributes to a gain of resistance against PepVMV and ChiVMV but also to a decrease of the durability of the pvr2 3 resistance to PVY together with a decrease in quantitative resistance traits VA and AUDPC. The Perennial inbred line was selected from a landrace issued from farm selection in India, where several potyviruses including ChiVMV are prevalent in pepper crops and it is highly likely that pvr2 3 and pvr6 have been selected locally to enlarge the spectrum of resistance against these potyviruses. Whether pvr6 is the functional gene for the RB-3 QTL remains to be determinated through functional validation. However, the tight colocation between pvr6 and RB-3 already reveals a trade-off between a large resistance spectrum against potyviruses when combined with a pvr2 resistance allele, and the durability of the pvr2 3 PVY resistance in the Perennial genotype.

Such a trade-off looks unusual since resistance gene pyramiding and broad-spectrum resistance are often considered favourable to durability [START_REF] Mcdonald | Pathogen population genetics, evolutionary potential, and durable resistance[END_REF][START_REF] Kou | Broad-spectrum and durability: understanding of quantitative disease resistance[END_REF]. This statement was recently reinforced by Le [START_REF] Van | Differential selection pressures exerted by host resistance quantitative trait loci on a pathogen population: a case study in an apple 9 Venturia inaequalis pathosystem[END_REF] who showed that broad spectrum resistance QTLs in apple did not exert any differential selection pressures on a mixed Venturia inaequalis inoculum, whereas narrow spectrum QTLs did. In pepper, the combination of pvr2 resistance alleles (particularly pvr2 1 and pvr2 2 alleles) with pvr6 is already used in breeding programs to create resistant varieties against a large range of potyviruses [START_REF] Rubio | Allele mining in the pepper gene pool provided new complementation effects between pvr2-eIF4E and pvr6-eIF(iso)4E alleles for resistance to pepper veinal mottle virus[END_REF] and such cultivars may lead to a premature breakdown of pvr2-mediated resistance by PVY. More generally, translation initiation factor-mediated virus resistance is widespread [START_REF] Robaglia | Translation initiation factors: a weak link in plant RNA virus infection[END_REF] and combination between mutated or knock-out alleles at different genes coding for to provide a way to breed for broad spectrum resistance. In tomato, the narrow range resistance conferred by a knock-out allele of eIF4E1 was enlarged when combined with a knock-out of eIF4E2 [START_REF] Mazier | Knock-Down of Both eIF4E1 and eIF4E2 Genes Confers Broad-Spectrum Resistance against Potyviruses in Tomato[END_REF]. Our study indicates combining mutant alleles from genes belonging to the same multigenic family should be used carefully, since the trade-off observed in pepper between large-spectrum and durability of resistance, may also occur in other plants species.

Our previous observations indicate introgressing resistance alleles from the plant germplasm into new elite cultivars with susceptible genetic backgrounds can endanger the long-term use of these genes as well as provide an evolutionary springboard to the pathogen for resistance breakdown [START_REF] Palloix | Durability of plant major resistance genes to pathogens depends on the genetic background, experimental evidence and consequences for breeding strategies[END_REF]. Investigating the background-dependence of the durability of pvr2 3 led to an improved understanding of how selection acted on this gene/genetic background combination [START_REF] Chandler | Does your gene need a background check? How genetic background impacts the analysis of mutations, genes, and evolution[END_REF]. In mass selection of plant populations by farmers, only the seeds from the healthiest and most productive individuals participate to the next cultivated generation. In North-West India where Perennial originated, plants combining pvr2 3 and pvr6 certainly gain a selective advantage in presence of ChiVMV. However, the rapid adaptation of PVY and/or other potyviruses would have counterselected this allelic combination in a few cultivation seasons. Over a few plant generations, and in heterogeneous plant populations, this multi-pathogen selection pressure promoted the recombinant individuals carrying favorable alleles at additional QTLs (RB-1, RB-5 and RB-6) which, in interaction with RB-3, decreased the risk of pvr2 3 resistance breakdown. Such a polygenic combination of co-adapted alleles contributed to increase resistance in efficiency, spectrum and durability.

Figure 1 :

 1 Figure 1: Histograms represent distributions of DH lines into resistance breakdown frequency, virus accumulation (ratio concentration of the sample/concentration of the control) and AUDPC classes. Arrows indicate the position of each parent in the phenotypic distribution. Per : Perennial, YW : Yolo Wonder.

Figure 2 :

 2 Figure 2: Map positions on Pepper genetic map of the significant quantitative trait loci (QTLs) for pvr2 3 resistance breakdown (RB) frequency, the area under the disease progress curve (AUDPC) and the virus accumulation (VA). On the left, the scale indicates map distances calculated by the Kosambi mapping function. For each QTL, the length of the block represents the confidence interval (CI) of the QTL using a MCQTL test unit fall of 1 and the length of the vertical bar represents the CI of the QTL using a MCQTL test unit fall of 2. For each QTL, the closest marker is underlined. Lines indicate epistatic effects between loci.

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	(ii) comparative mapping between these QTLs and QTLs affecting the quantitative resistance to
	PVY.
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	3

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	mapping (SIM)
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	6

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	7

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.Manuscrit d'auteur / Author manuscript and RB-5 according to the trait and the chromosome location. Table1 and figure

	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	8

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	9

, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html. RB variants in the resistant plants, two steps of the virus evolution toward RB that are considered as major components of the resistance durability

	accumulation dynamics of
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	10

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	11

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

  Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.In plants, a decrease of symptom intensity can result from resistance and/or tolerance mechanism.

	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	12

a rational basis for sustainable resistance breeding in plants.

Heredity, 112 (6)

, 579-587. DOI : 10.1038/hdy.2013.138

  Comment citer ce document : Quenouille, J., Paulhiac, E.,Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138 Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

	translation initiation factors is expected
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript
	Manuscrit d'auteur / Author manuscript	13

  Comment citer ce document : Quenouille, J., Paulhiac, E.,Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138 Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html. Version définitive du manuscrit publié dans / Final version of the manuscript published in Heredity, DOI : 10.1038/hdy.2013.138. The original publication is available at http://www.nature.com/hdy/index.html.

Comment citer ce document : Quenouille, J., Paulhiac, E., Moury, B., Palloix, A. (Auteur de correspondance) (2014). Quantitative trait loci from the host genetic background modulate the durability of a resistance gene: a rational basis for sustainable resistance breeding in plants. Heredity, 112 (6), 579-587. DOI : 10.1038/hdy.2013.138

ACKNOWLEDGEMENTS

The authors thank B. Mangin (INRA-Toulouse) for help in the use of MCQTL software, G. Girardot, P. Mistral, G. Nemouchi, K. Nozeran, B. Savio and V. Simon for technical assistance, A.M. Sage-Palloix for providing the Capsicum genetic resources, J-P. Bouchet for bioinformatic assistance and K. Lamour (Tennessee Univ., USA) for critical reading of the manuscript. This work was financially supported by the Comité Technique Permanent de la Sélection (CTPS, French Ministry of Agriculture and Fisheries), by the Agence Nationale de la Recherche (Project VirAphid ANR-2010STRA-001-01), and by the Région Provence Alpes Côte d'Azur (PACA).

DATA ARCHIVING

Data available from the Dryad Digital Repository: doi:10.5061/dryad.s4n76

CONFLICT OF INTEREST

The authors declare no conflict of interest.