| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2007, p. 5162-5172, Vol. 73, No. 16
0099-2240/07/$08.00+0 doi:10.1128/AEM.00507-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Institut National de la Recherche Agronomique, UMR Santé Végétale (INRA-ENITAB), Institut des Sciences de la Vigne et du Vin, BP 81, 33883 Villenave d'Ornon Cedex, France,1 Institut National de la Recherche Agronomique, UMR Santé de la Vigne et Qualité du Vin (INRA-ULP), BP 20507, 68021 Colmar Cedex, France2
Received 6 March 2007/ Accepted 9 June 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Numerous experimental and phylogenetic studies in extensive geographical contexts have already examined the pattern of the evolution of genes coding for antimicrobial drugs or pesticide targets in animals, plants, and bacteria. For instance, the question of parallel genetic adaptation to drugs and pesticides has already been addressed for various insecticide resistances (2-4, 14, 42, 55), for herbicide resistance (11, 16), for resistance to a nematocide (17), and for resistance to antibiotics (8, 32). However, to our knowledge, no study has ever addressed the question of the number of origins of fungicide resistance alleles in plant pathogen populations. To study the evolutionary process that leads to the onset of fungicide resistance alleles, we investigated the evolution of the target site for quinone outside inhibiting (QoI) fungicides in the plant pathogen Plasmopara viticola.
Launched in the late 1990s, strobilurin and the other fungicides belonging to the group of QoI fungicides are an important class of chemicals for the management of a broad range of fungal diseases in agricultural systems. These fungicides are effective against three major groups of plant pathogens: oomycetes, basidiomycetes, and ascomycetes (9, 34). The QoI fungicides act as mitochondrial respiration inhibitors by binding to the ubiquinol oxidation site (Qo) formed by the cytochrome b and iron-sulfur protein domains within the cytochrome bc1 complex (10, 25). By compromising ATP production in mitochondria, QoI fungicides affect the energy-demanding stages of pathogen development, such as spore germination.
The first QoI fungicide (azoxystrobin) was launched for use in agricultural pest management in 1996. One year later, the development of field-reduced sensitivity to azoxystrobin in populations of the plant pathogen Mycosphaerella fijiensis was first observed (49). As for other single-site inhibitors, the longevity of QoI fungicides was thus quickly challenged in the field, and resistant isolates have been described for numerous plant pathogen species, such as Alternaria solani, Blumeria graminis, Mycosphaerella graminicola, Plasmopara viticola, Podosphaera fusca, Puccinia species, Magnaporthe oryzae, and Venturia inaequalis (5, 21, 22, 30, 33, 35, 38, 44, 49, 50, 52, 54, 56). In most cases, resistance to QoI fungicides is conferred by a point mutation in the mitochondrial cytochrome b gene, giving rise to a substitution from glycine to alanine at position 143 (G143A) of the amino acid sequence (25, 29, 35). In addition, another mutation at amino acid sequence position 129 of cytochrome b resulting in a substitution from phenylalanine to leucine (F129L) has also been reported to confer resistance to QoI fungicides in species such as Pythium aphanidermatum, Magnaporthe oryzae, and Alternaria solani, although the level of resistance is lower than that conferred by the G143A substitution (25, 35, 44).
P. viticola (Berkl. and Curt.) Berl. and de Toni., the causal agent of grapevine downy mildew disease, is a native species from North America that was introduced several times into Europe in the late 1870s (28). The disease spread rapidly through most of continental Europe to become one of the most important vineyard diseases (23). Chemical control is the most effective measure currently used to protect grapes from downy mildew disease. QoI fungicides were introduced into France for pest management in viticulture in 1998. Only 2 years later, QoI fungicide resistance was detected in P. viticola populations of most French vineyards (39), and the nonsynonymous G143A substitution was detected using a quantitative PCR method (48, 51). Moreover, the results of a recent comparative survey of the cytochrome b gene structure in various plant pathogen species suggest that the F129L substitution of the cytochrome b gene could also be involved in resistance to QoI fungicides in P. viticola (29).
The recent selection of well-characterized mutations associated with QoI fungicide resistance in natural populations of P. viticola provides a unique opportunity to understand the mechanisms leading to the evolution of fungicide resistance in a plant-pathogen system. Because the mitochondrial genome does not recombine, QoI fungicide resistance is a good model to investigate the number of origins of resistance alleles. In this respect, we address three fundamental questions concerning the evolution of QoI fungicide resistance haplotypes. (i) How many mutations are involved in resistance to QoI fungicides in grapevine downy mildew populations? (ii) How often do these mutations arise in natural populations of grapevine downy mildew populations, i.e., is resistance the result of a single mutation in an allele that spreads subsequently to different areas or rather the result of multiple convergent acquisitions of the same mutation in the gene coding for the cytochrome b target? (iii) What is the genetic structure of QoI fungicide resistance alleles in French P. viticola populations after 5 years of use of this group of fungicides?
Based on the complete sequence of the cytochrome b gene, nucleotide polymorphism variations were analyzed in a broad range of sensitive and resistant isolates of P. viticola collected in European and American vineyards. Patterns of nucleotidic polymorphisms confirmed the occurrence of a unique resistance mutation in the cytochrome b gene, while phylogenetic analyses of the complete fragment provided evidence for at least two independent origins of QoI fungicide-resistant haplotypes in natural populations of P. viticola. In addition, the genetic structure of QoI fungicide-resistant mitochondrial haplotypes in 17 French populations is discussed in relation to local fungicide selective pressures.
| MATERIALS AND METHODS |
|---|
|
|
|---|
A total of 819 downy mildew lesions on Vitis vinifera were collected. Each sample consisted of 1 cm2 of a fresh leaf colonized by sporulating downy mildew (also called oil spots). Since previous studies have shown that a single lesion rarely includes more than one multilocus genotype (27), we treated each lesion (oil spot) as a single strain.
Samples came from 51 localities in France (n = 799); 20 localities in Italy, Switzerland, Germany, Portugal, Spain, Greece, and Romania (n = 51); and 1 locality in the United States (n = 14) (Table 1). In France, sampling was carried out during June and July of the grapevine growing season in 2003 (Table 1). For samples collected abroad, the years of collection ranged from 1997 to 2004.
|
DNA extraction.
Oil spots were freeze-dried overnight, and DNA was extracted from each of them (containing infected grapevine leaf tissue and sporangia) according to the standard cetyltrimethylammonium bromide-phenol-chloroform method described previously (58) with a slight modification (phenol-chloroform was replaced with isoamyl alcohol and chloroform at a dilution ratio of 24:1). After extraction, DNA was precipitated with isopropanol, and the pellet was washed with 75% ethanol and resuspended in 50 µl of sterile water.
Fungicide in vitro sensitivity tests.
To characterize P. viticola isolate sensitivity to the QoI fungicides, 27 single lesions were multiplied onto leaves of Vitis vinifera cv. Cabernet-Sauvignon. Sporangia were collected and suspended in sterile water at 4°C. The inoculation was done by depositing fifteen 10-µl droplets onto the adaxial face of each leaf. The leaves were incubated in petri dishes under controlled conditions (21°C) with a 16 h photoperiod. After 7 days, freshly produced sporangia were harvested to inoculate the test units.
Six discriminatory doses (10, 100, 300, 600, 1,000, 1,500 mg liter–1) of an aqueous solution of famoxadone fungicide (96%) were applied in the bioassay. For sensitive isolates that did not grow at the fungicide concentration of 10 mg/liter, eight other fungicide doses were applied to determine the MIC: 0.1, 0.2, 0.5, 0.8, 1, 2, 5, and 10 mg liter–1. For each concentration (including a control), 3 ml of fungicide was sprayed using a hand-held sprayer (Ecospray) onto the lower side of 10 detached grapevine leaf disks placed in a petri dish. The upper sides of the leaves were placed in contact with filter paper impregnated with 4 ml of water. After the surfaces were dried at room temperature, the leaf disks were inoculated with 3 droplets of 10 µl per disk of an inoculum of 20,000 sporangia per ml. After 7 days of inoculation, development of mildew on each disk was visually estimated as the proportion of leaf area with sporulation (24). For each fungicide concentration, 10 leaf disks were tested. The average score for each fungicide concentration was converted to a percentage of inhibition by comparison with the untreated disks. The fungicide concentration inhibiting 50% of growth (IC50) and the MIC of growth (MIC) discriminated two classes of isolates according to their sensitivity to famoxadone: famoxadone-sensitive isolates (IC50 and MIC <10 mg liter–1) and highly resistant isolates (IC50 and MIC >100 mg liter–1). For all isolates, DNA was extracted for subsequent sequencing to characterize molecular variations in the cytochrome b gene sequence.
Isolation of cytochrome b gene and flanking regions.
Initial amplification of the mitochondrial genome fragment was obtained by performing a PCR with degenerate primers (P1) designed from two oomycete sequences (Phytophthora infestans, NC002387, and Phytophthora megasperma, L16863) and a reverse primer (P2) specific to P. viticola cytochrome b (partial sequence, AX577570) (Table 2; Fig. 1). Each reaction was carried out in a PCR of 30 µl containing 0.2 µM of each primer, 0.33 mM of each deoxynucleoside triphosphate (dNTP), 1x reaction buffer, 1.5 mM of MgCl2, 0.33 mM of each dNTP, 0.5 U of Taq polymerase, and 1 µl of genomic DNA. PCR cycles consisted of 95°C for 4 min, 35 cycles of 95°C for 45 s, 57°C for 45 s, and 72°C for 50 s, and a final extension step of 72°C for 7 min. The PCR products were purified and used for sequencing.
|
|
Mitochondrial sequence diversity.
A mitochondrial fragment of 2,281 bp including the cytochrome b gene was amplified and sequenced for 56 European and 14 American samples (Table 1). In addition, a partial sequence was generated on the 5' end, the middle part, and the 3' end of the complete fragment (length ranging from 700 bp to 1,500 bp) for 119 additional samples. This resulted in a data set of 182 mitochondrial DNA sequences that allowed the detection of 23 polymorphic sites (single nucleotide polymorphisms [SNPs]). Primers used for the sequencing are presented in Fig. 1 and Table 2.
SNP typing.
SNP typing was performed by using a primer extension genotyping assay with fluorescence polarization (FP) detection (12) implemented in the AcycloPrimeTM-FP SNP detection kit (PerkinElmer Life Science, Inc.) from a preamplified 630-bp SNP-containing fragment synthesized using primers P7 and P9 (Table 2; Fig. 1). The PCR was performed in a 30-µl volume containing 0.13 µM of each primer, 0.2 mM of each dNTP, 1x reaction buffer, 1.5 mM of MgCl2, 0.33 mM of each dNTP, 0.5 U of Taq polymerase (Silverstar; Eurogentec Inc.), and 1 µl of genomic DNA. Subsequent procedures were conducted as follows: (i) the AcycloPol polymerase was used to extend by one base a SNP primer (P15 to 17) (Table 2; Fig. 1) that ends immediately upstream of a targeted SNP position using one of two fluorescent dye-labeled terminators, depending on the SNP profile previously determined, and (ii) the identity of the terminal base added was determined by the increased FP of its linked dye (scanned and analyzed by Vector3TM, PerkinElmer Life Science, Inc.). PCR cycles for the primer extension reaction consisted of 95°C for 6 min and 35 cycles of 95°C for 15 s and 55°C for 30 s. We used a lower Tm temperature (51°C) for the assay using P16 so as to optimize the result. The reaction was prepared on a 96-well dark plate and subsequently visualized on a Vector3TM scanner.
Intraspecific phylogenetic analyses.
Two graphical representations obtained by the following methods were used to illustrate the evolutionary relationships among haplotypes derived from 70 mitochondrial sequences of P. viticola samples from Europe and North America. First, a phylogram was generated from phylogenetic analysis as implemented in PAUP version 4.0b10 (53) using the criterion of maximum parsimony (MP) under the exhaustive search option. The resultant MP tree was rooted using the midpoint method as implemented in PAUP. Statistical support of nodes was evaluated using MP bootstraps (19) with 1,000 replicates. Second, a haplotypic network was constructed using the median-joining network approach (7) as implemented in the program NETWORK, version 4.111 (available at http://www.fluxus-engineering.com/sharenet.htm).
Population genetic analysis.
The genetic structure was examined by analysis of molecular variance (AMOVA) using the Arlequin 2.001 software package (18). This method was used to partition the genetic variance within populations, among populations within groups, and among groups. The populations were grouped into regions according to their geographical origin. Levels of significance were determined through 1,000 random permutation replicates. Fisher exact tests for population differentiation were performed using Genepop version 3.1b (47).
Nucleotide sequence accession numbers.
The nucleotide sequences detected in this study have been deposited in GenBank under the accession numbers AY696297 and DQ459459 to DQ459469.
| RESULTS |
|---|
|
|
|---|
Mitochondrial sequence variability.
A set of 70 complete mitochondrial sequences with 2,281 nucleotides and 119 partial sequences (between bp 700 and 1500) was examined representing 71 different localities in European countries and 14 samples from the United States. Twenty-three polymorphic sites were detected among all the sequences (Table 3). Except for one nonsynonymous mutation in cytochrome b (position 1256), the 22 polymorphic sites were either synonymous substitutions or mutations in noncoding regions. Ten polymorphic sites were parsimony informative, and 11 different mitochondrial haplotypes were identified (Table 3).
|
|
|
Table 4 demonstrates that the mutation at amino acid position 143 is associated with the high level of famoxadone resistance in P. viticola. Moreover, the analysis of 189 complete or partial mitochondrial sequences from European and American samples failed to detect any mutation at amino acid position 129 of the cytochrome b gene. This result was also confirmed by SNP typing assays on 87 additional samples (data not shown) using the extended SNP primer P14 designed to detect this putative mutation.
Evolutionary relationships between mitochondrial haplotypes.
MP and median-joining network analyses were performed to infer the evolutionary relationships among 11 mitochondrial haplotypes identified from 70 complete sequences of the mitochondrial fragment including the cytochrome b gene (Fig. 2; Table 1). Only one MP tree with a tree length of 25 was found (Fig. 2A). The tree splits into two major clades (I and II) consisting of European haplotypes. These clades are supported by robust MP bootstrap values of 91 and 92, respectively (Fig. 2A). Both phylogenetic analyses yielded genealogy patterns revealing an unambiguous differentiation between European and American populations. The distance from either European clade to the most central missing intermediate among all haplotypes was found to be four mutational steps, whereas the distance between two European clades was six mutational steps (Fig. 2B). In contrast to the looser connection between American haplotypes, three haplotypes within each European clade were connected by a single mutational step (Fig. 2B).
At least two origins of QoI fungicide-resistant alleles.
Since the European QoI fungicide-sensitive (S) and -resistant (R) isolates were distributed between two well-supported groups separated by six mutational steps (clade I and clade II) (Fig. 2), we concluded that there are four major European haplotypes: IS and IIS for sensitive haplotypes and IR and IIR for resistant haplotypes (Fig. 1 and 2). No difference in famoxadone sensitivity was found between IS and IIS or between IR and IIR in bioassays (Table 4). According to the phylogenetic analyses, the two resistance haplotypes (IR and IIR) derived by one mutational step from the two QoI fungicide-sensitive haplotypes IS and IIS, respectively. At least two origins of QoI fungicide resistance haplotypes were thus revealed by the presence of several resistant alleles that differed from their potential susceptible progenitors only by the presence or absence of the resistance-associated mutation itself.
Mitochondrial haplotype distribution in French vineyards.
An SNP typing assay was used to survey the geographical distribution of various mitochondrial haplotypes from a panel of 676 P. viticola isolates collected in 17 French localities (Table 5). Depending on the downy mildew epidemics, the number of samples collected and genotyped ranged from 11 to 86, with a mean of 40 (±19) samples per population. Primers P7 and P9 (Fig. 1; Table 2) were used for DNA amplification and yielded a product of 630 bp. European clades I or II were established with extended SNP primers P16 and P17, and the extended primer P15 was used to detect the G143A mutation (Table 2). The four major European haplotypes, IS, IR, IIS, and IIR, were typed as GCA, CCA, GTG, and CTG, respectively, at nucleotidic positions 1256, 1281, and 1608 (Table 3). The same assay also differentiated between American and European haplotypes, since all American isolates shared the same specific haplotypic pattern (GTA), while no European isolate showed the American haplotypic pattern.
|
Concerning QoI fungicide resistance, at least one of the haplotypes conferring resistance was found in all populations sampled. Overall, the percentage of resistant haplotypes (IR and IIR) was 24.3% (±15.6%) in French populations but with considerable local differences, as indicated by the high standard deviation (Fig. 3). For instance, Loire and Rhone Valley populations are still apparently less penetrated by QoI fungicide resistance alleles (<10%), whereas in the southwestern and Champagne vineyards, the frequency of resistance haplotypes is rather high (32.1% and 21.3%, respectively). The effect of QoI fungicide treatments on the mitochondrial genetic structure was tested by comparing the latter in treated and untreated populations in Bordeaux vineyards. The Fisher exact test indicated that populations treated with QoI fungicides in 2003 presented a significantly higher number of resistance haplotypes than those not treated (P = 0.022).
|
|
| DISCUSSION |
|---|
|
|
|---|
Confirming previous findings (25, 29, 48, 51), we found that one major mutation at amino acid position 143 in cytochrome b provides a high level of resistance to QoI fungicide in P. viticola. However, contrary to Grasso et al. (29), we did not find any supplementary mutation at amino acid position 129 of the cytochrome b gene in 276 P. viticola isolates screened in 17 French vineyards. The lack of detection of the F129L substitution indicates that the frequency of this allele is below 0.36% in our sample, an estimate in agreement with monitoring data obtained by Sierotzki et al. (48), who evaluated the F129L mutation frequency to be 0.4% in French P. viticola populations in 2003.
Origins and spread of QoI fungicide resistance haplotypes.
Accumulation studies on QoI fungicide resistance have shown that a single replacement in a single gene occurred across a wide range of plant pathogens, clearly indicating parallel evolution of the resistance phenomena in different species (29). However, genetic data on the rate at which adaptive substitutions arise at the infraspecific level are still lacking. For the first time in a plant-pathogen system, we used phylogenetic arguments to demonstrate that a selected substitution (G143A) in the cytochrome b gene has occurred at least twice during the course of evolution of QoI fungicide resistance alleles in European populations of P. viticola. This finding confirms that parallel evolution of resistance to fungicide is acting not only between but also within species. This result is in agreement with expectations from the mutation-selection balance theory. The equilibrium allele frequency under mutation-selection balance for a haploid is expressed by µ/s, where µ is the mutation rate and s is the selection coefficient. Assuming complete selection (s = 1) against the resistant haplotype, the equilibrium frequency of a mutation will be that of the neutral substitution rate per site for mtDNA, which has been estimated to be approximately 2.10–8. Given that the population size of grapevine downy mildew must readily exceed 108, this calculation, while extremely rough, shows that it would not be surprising to find multiple independent origins of QoI fungicide resistance haplotypes in populations of P. viticola.
The double origins of resistance haplotypes estimated from our data set could be a low estimate of the number of recurrent substitutions that might have occurred at the selected site of the cytochrome b gene. Indeed, mtDNA genetic variability of P. viticola was extremely low in Europe. Two hypotheses could explain this. First, given the absence of normal recombination in the mitochondrial genome, any advantageous mutation such as the G143A will lead to the fixation of all other polymorphisms in the population by a process known as hitchhiking (40). This hypothesis is unlikely, since we did not find more haplotypes in European populations that never received QoI fungicide applications (data not shown). The alternative explanation is that the bottlenecking of populations resulting from the recent and multiple introduction of P. viticola in Europe (i.e., since 1878) (28) has greatly reduced mitochondrial variability in European populations. This lack of variability may have reduced our ability to detect phylogenetic events resulting in the emergence of the different resistance haplotypes. To improve our estimation of the rate of these recurrent mutations (G143A), it would be useful to analyze the mitochondrial variability of native populations of P. viticola (in the United States) that present a higher level of nuclear and mitochondrial genetic diversity, as illustrated by our results. Analysis of the mitochondrial polymorphism and the derived phylogenies of sensitive and resistant haplotypes from American isolates would certainly allow more precise evaluation of the number of independent replacements of the amino acid conferring resistance to QoI fungicides in P. viticola. The relevance of this approach is enhanced by the fact that grapevine downy mildew populations have recently been described to develop resistance to QoI fungicides in the United States (A. Baudoin, personal communication). Although recurrent mutation is probably an important initial force driving the evolution of QoI fungicide resistance-associated point mutations in this pathogen species, the role of migration in the spread of resistance haplotypes across vineyards should not be underestimated. Indeed, when combined with directional selection, long-distance migration of sporangia (by passive migration or by human activities) could be a primary determinant of the spread of these kinds of adaptive mutations over broadly distributed vineyards. The impact of gene flow between distant populations on the final distribution of resistance alleles has been well illustrated in Culex pipiens: Raymond et al. (46) showed that frequent migration events, including passive transportation, have resulted from the worldwide spread of resistance alleles to organophosphorous insecticides.
Patterns of genetic variation of mtDNA.
Mitochondrial genetic structure estimated by SNP typing of 17 populations was high (FST = 0.11, P < 10–3); most of the genetic variability was found within populations, and genetic variability was loosely structured by geography at the regional scale. The high mitochondrial genetic structure observed contrasts with previous estimates of population differentiation of P. viticola obtained using microsatellite markers (15, 26). Indeed, population differentiation observed on a scale similar to that found in our study has provided evidence for a low genetic structure (FST = 0.018 among French and Swiss populations, P < 10–3). The difference between these two FST estimates is not surprising, since it has been demonstrated that the mitochondrial genome is often more structured than the nuclear genome, owing to its smaller effective size (45) and to small fitness differences between mtDNA genomes (6). In our case, the differences observed between nuclear and mitochondrial markers likely resulted from differences in selective regimen. Indeed, population differentiation by adaptive genetic variations results from selection and genetic drift, while population differentiation by neutral polymorphism is only driven by genetic drift. Given that mtDNA variations in P. viticola influence pathogen fitness, populations of the pathogen might be under strong directional selection for local adaptation to fungicide pressure. In agreement with this, we found a significant increase in resistance haplotype (IR and IIR) frequency in populations treated with QoI fungicides during the sampling year. Similarly, removing treated populations from the AMOVA resulted in a drastic reduction of genetic differentiation among populations within regions. These results demonstrate that, in the same year, the mitochondrial genetic structure responds to the spatial heterogeneity resulting from fungicide treatments. The ability of P. viticola to quickly adapt to fungicide pressure is likely enhanced by its ability to undergo clonal amplification of resistant genotypes by several asexual generations during the growing season.
Even though we detected a significant global effect of fungicide treatments on mitochondrial population differentiation, it is still rather difficult to establish a solid relation between resistance haplotype frequencies in populations and the selective pressure exerted by fungicides. In fact, selection pressure might not be uniform among vineyards because many factors contribute to fungicide treatment efficacy, such as disease pressure, local climatic conditions, timing, intervals, and rates of spraying programs. Thus, selection by fungicides leads to a pattern of local adaptation of the pathogen and, as evidenced by the high genetic structure, this pattern is not homogenized by gene flow among populations during the grapevine growing season (26).
Another significant result of this study is that all populations either treated or not treated with fungicides presented a large proportion of QoI fungicide-resistant isolates. This large proportion of isolates could result from local migration, because control plots should be considered as small, untreated islands surrounded by treated vineyards in which selection for resistance is operating. Moreover, the untreated plots in this study likely received fungicide applications during the years preceding the sampling year. Thus, resistant genotypes may have remained in control plots from one year to another. The latter hypothesis raises the question of the factors determining the maintenance of resistance haplotypes in the field, in particular when chemical control is suspended. Fungicide resistance alleles are frequently costly, and their cost is variable (13). When resistance alleles carry a fitness cost, the rise in resistance is followed by a fall when chemical control is suspended. Maintenance of resistance alleles in untreated plots could indicate an absence of cost for QoI fungicides. However, another hypothesis could be that the initial fitness cost gradually decreases, owing to subsequent mutations at secondary loci that compensate for the cost (37). Unfortunately, reliable data on the existence of a cost associated with the QoI fungicide resistance mutations in the absence of fungicide pressure are lacking for P. viticola. Clearly, more studies are now needed to understand the complex interaction between mutations, migration, and competition between these QoI fungicide resistance haplotypes.
Implications for the evolution of resistance.
Many interacting factors contribute to the appearance and spread of fungicide resistance alleles in pathogen populations. The degree of risk of resistance appearance and development is dependent on the mode of action of the fungicide, the selective pressure resulting from its application, the local conditions for disease development, and the evolutionary potential of the pathogen target (migration rate, mutation rate, and effective population size) (41). In this study, we show that the large population size of P. viticola combined with strong selective pressure has resulted in at least two independent origins of QoI fungicide resistance alleles. We contend that the recurrent appearance of single resistant-associated mutations may be an important way in which pathogen or pest species respond to the challenge imposed by chemical selection. This process is probably linked with the kind of genomic modification involved in the resistance, with SNPs associated with target site insensitivity being more readily subject to recurrent mutations than more complex genome rearrangements (gene amplification or gene regulation). In accordance with this hypothesis, multiple origins of resistance alleles have already been described for numerous resistances to pesticides and drugs in plants, animals, protozoa, and bacteria. A parallel evolution has been evidenced for resistance to cyclodienes, pyrethroids, and dichlorodiphenyltrichloroethane insecticides in several insect species (1-4, 20, 42, 55), for resistance to herbicides inhibiting acetyl-CoA carboxylase in plants (16), for resistance to chloroquine in malaria (57), and for resistance to beta-lactamate in bacteria (8, 32). There is now growing evidence for independently acquired adaptive substitutions during the evolution of resistances to pesticides and drugs. This finding will certainly help to define the scale at which resistance evolves and will provide valuable insights into our models in order to devise appropriate resistance management strategies for microbial populations.
| ACKNOWLEDGMENTS |
|---|
We thank the following colleagues and institutes for help in sample collection: L. Delbac from INRA, Marc Raynal from the Institut Technique de la Vigne et du Vin (ITV) in Bordeaux, Marie-Laure Panon from the Comité Interprofessionel des Vins de Champagne (CIVC), D. Gobbin and C. Gessler from ETH Zurich, and M. Kennelly and D. Gadoury from Cornell University. We are grateful to R. Caroli, S. Menaouli, P. Cosson, V. Leger, J. Jolivet, and G. Taris for their laboratory assistance and J. L. Genet from Dupont Protection for the gift of famoxadone. We thank M. Milgroom and R. Cooke for reviewing a previous version of the manuscript.
| FOOTNOTES |
|---|
Published ahead of print on 22 June 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||