Effect of hosts on competition among clones and evidence of differential selection between pathogenic and saprophytic phases in experimental populations of the wheat pathogen Phaeosphaeria nodorum
© Sommerhalder et al; licensee BioMed Central Ltd. 2011
Received: 10 April 2011
Accepted: 1 July 2011
Published: 1 July 2011
Monoculture, multi-cropping and wider use of highly resistant cultivars have been proposed as mechanisms to explain the elevated rate of evolution of plant pathogens in agricultural ecosystems. We used a mark-release-recapture experiment with the wheat pathogen Phaeosphaeria nodorum to evaluate the impact of two of these mechanisms on the evolution of a pathogen population. Nine P. nodorum isolates marked with ten microsatellite markers and one minisatellite were released onto five replicated host populations to initiate epidemics of Stagonospora nodorum leaf blotch. The experiment was carried out over two consecutive host growing seasons and two pathogen collections were made during each season.
A total of 637 pathogen isolates matching the marked inoculants were recovered from inoculated plots over two years. Genetic diversity in the host populations affected the evolution of the corresponding P. nodorum populations. In the cultivar mixture the relative frequencies of inoculants did not change over the course of the experiment and the pathogen exhibited a low variation in selection coefficients.
Our results support the hypothesis that increasing genetic heterogeneity in host populations may retard the rate of evolution in associated pathogen populations. Our experiment also provides indirect evidence of fitness costs associated with host specialization in P. nodorum as indicated by differential selection during the pathogenic and saprophytic phases.
The evolution of pathogens is widely believed to be one of the major challenges facing agriculture and medicine [1, 2]. Experimental studies focused on the evolution of pathogens, including the emergence of virulence and pathogen adaptation to changing agricultural and medical practices, can provide critical information for more effective management of infectious diseases. In medicine, infectious diseases are mitigated mainly through the application of antimicrobial substances such as antibiotics. While pesticides such as fungicides are widely utilized in agricultural ecosystems, host resistance imposes many fewer environmental costs and is a more cost efficient approach to control plant diseases. In both agriculture and medicine, the efficacy of host resistance and antimicrobials usually decays over time as a result of the continuous evolution and adaptation of pathogens.
Plant pathogens are thought to evolve faster in agricultural ecosystems than in natural ecosystems [3–6]. Wild relatives are the primary sources of host resistance bred into cultivated crops. The disease resistance genes carried by these wild relatives of modern crop plants have coexisted with their pathogens for many thousands or millions of years in natural ecosystems. However, when these resistance genes are introgressed into modern crops and deployed in agricultural ecosystems, their value in controlling infectious diseases usually does not last for more than 10 years . Temporal analysis of population dynamics also is consistent with the hypothesis of rapid pathogen evolution in agricultural ecosystems. For example US-1 was the predominant genotype in Phytophthora infestans populations around the world until the 1990s  but this genotype is rarely recovered since 2000. In the UK, a single P. infestans genotype called Blue-13-A2 was first detected in southern England in 2003 at a very low frequency. By 2007, this genotype was detected in all populations sampled across the UK and accounted for more than half of 1000 isolates assayed (J. Zhan & D. Cooke, unpublished data). These P. infestans examples illustrate how pathogen populations can experience rapid turnover as new genotypes with greater fitness emerge, spread, out-compete and replace earlier genotypes.
The evolution of pathogens can be influenced by the type of resistance, the amount of diversity found in host populations and the type of cropping system [9–12]. Modern agriculture is dominated by species monocultures grown at a high density. In these agricultural ecosystems, it is common for a single host cultivar or genotype carrying a major resistance gene to be grown over a large area. The limited genetic diversity in the host populations coupled with intensive use of major resistance genes can lead to rapid shifts in associated pathogen populations. A mutant with higher fitness that emerges in a pathogen population as a result of a single mutation event can quickly increase in frequency through strong directional selection and spread across entire fields or regions through natural or human-mediated migration.
Multi-cropping, where the same annual crop is grown in the same field more than once during the same year, is another common practice in modern agriculture, especially in countries experiencing a shortage of arable land. This practice may further accelerate the evolution of plant pathogens because locally adapted pathogen genotypes with a high parasitic fitness can steadily increase in frequency due to the year-around availability of the living host (i.e. a "green bridge" allows the parasitic phase of the pathogen life cycle to occur continuously).
It is hypothesized that the evolution of plant pathogens in agricultural ecosystems can be retarded by increasing genetic diversity of the host populations, by using partial resistance encoded by several genes and by avoiding multi-cropping systems. Increasing genetic diversity in host populations by mixing plants carrying different major resistance genes (e.g. cultivar mixtures or multilines) is thought to be an ecologically and evolutionarily sound approach to control plant diseases, particularly for airborne pathogens of cereals . Increasing host diversity by using cultivar mixtures will impose disruptive selection on pathogen populations, i.e. pathotypes that are favored on one host will have lower fitness on the other hosts in the mixture [14–16], impeding their ability to evolve towards higher virulence, here defined as the damage a pathogen causes to its host . On the other hand, because many fungal pathogens have large effective population sizes  and exhibit a mixture of sexual and asexual reproduction [19, 20], they can rapidly obtain new pathogenicity factors through mutation or new combinations of pathogenicity factors through recombination and then maintain the novel combinations of pathogenicity factors through asexual reproduction. Thus extensive use of cultivar mixtures could lead to the development of complex races [21, 22] that would be able to infect a large number of host genotypes carrying different major resistance genes.
Though less efficient, partial resistance is thought to offer a more durable method to control plant diseases than major-gene resistance because it works against all pathogen strains and selects equally against all pathotypes [23, 24]. Partial resistance mediated by multiple genes is generally inherited as a quantitative trait [16, 25, 26], where each gene makes a minor but additive contribution to the overall resistance . But selection can increase the frequencies of genes encoding higher virulence in pathogen populations infecting partially resistant hosts and reduce the effectiveness of quantitative resistance [9, 19, 28–31] though possibly at a slower pace compared to major resistance genes .
In contrast to multi-cropping, in single cropping systems an annual crop is grown for only 6-9 months of the year or different crops are rotated annually, forcing pathogens to undergo a saprophytic phase in their life cycle in which different strains not only compete with each other but also with other microbial species for nutrients and habitats. Pathogen genotypes that have a high parasitic fitness on living hosts may have a low saprophytic fitness on the dead host biomass. This trade-off could delay the emergence of highly parasitic pathogen strains in agricultural ecosystems characterized by single cropping and regular crop rotations.
Recovery of inoculants
A total of 637 isolates matching the multilocus haplotypes of the inoculants were recovered from the inoculated plots. In addition to the inoculants, a large number of isolates (550) sampled over the course of this experiment were novel isolates with multilocus haplotypes that did not match the nine inoculants. The frequency of these novel haplotypes increased steadily over the course of the experiment. The majority of the novel genotypes were detected only once and the most frequent one was detected five times from two adjacent plots in 2005B.
Variation in genotype frequencies among inoculants
Comparisons of genotype frequency among P. nodorumpopulations from different hosts
Pair-wise comparisons for differences in genotype frequencies between Phaeosphaeria nodorum collections made from different hosts in 2004 based on a contingency χ2 test.
Pair-wise comparisons for differences in genotype frequencies between Phaeosphaeria nodorum collections made from different hosts in 2005A based on a contingency χ2 test.
χ2 tests for difference in genotype frequencies among Phaeosphaeria nodorum populations sampled from different hosts at the same point in time.
55.48 (32) **
51.67 (28) **
Changes in genotype frequencies among P. nodorumpopulations over time
Contingency χ2 tests for differences in genotype frequencies between Phaeosphaeria nodorum collections made from the same host treatment at different sampling times.
2004A vs. 2004B
2004B vs. 2005A
2004A vs. 2005A
17.50 (7) **
32.37 (16) **
15.87 (8) *
30. 13 (16)*
19.54 (8) **
31.64 (16) **
Contingency χ2 tests for differences in genotype frequencies between Phaeosphaeria nodorum collections sampled at different times.
44.67 (8) ****
40.57 (8) ****
Average selection coefficients and their standard deviations (in parentheses) for five released pathogen strains on each host treatment.
0.51 (0.08) a
0.03 (0.03) a
0.34 (0.13) a
0.15 (0.14) c
0.42 (0.09) b
0.48 (0.10) a
0.59 (0.09) b
0.44 (0.15) a
0.03 (0.03) c
0.82 (0.07) a
0.27 (0.10) b
0.37 (0.11) c
0.43 (0.15) a
0.52 (0.05) b
0.32 (0.10) c
0.30 (0.11) b
0.80 (0.06) a
0.40 (0.16) a
0.71 (0.28) a
0.29 (0.09) c
0.29 (0.12) b
0.83 (0.05) a
0.43 (0.14) a
0.67 (0.06) a
0.05 (0.04) b
Analysis of variance for selection coefficients of the five most frequent inoculated strains of Phaeosphaeria nodorum.
Cultivar * Isolate
Rapid change in the composition of P. nodorumpopulations
Because the epidemics were initiated by artificially inoculating the five host treatments with the same P. nodorum population (i.e. the mixture of nine marked strains in equal proportions), our null hypotheses were that the frequencies of the nine released isolates would be nearly equal in different host populations and that the genetic composition of these populations would not change over time. Instead, we observed significant differences in the frequencies of the released isolates and the majority of P. nodorum populations sampled from the five host treatments changed significantly over time. The differences in genotype frequency among the inoculated strains within a host population and among host populations sampled from different points in time could be due to random genetic drift or natural selection, but we believe that the observed differences in this case should be attributed mainly to selection. We have two lines of evidence supporting this hypothesis. 1) If genetic drift was the main factor, we would expect random changes in genotype frequencies among the P. nodorum populations sampled from different hosts. Instead, we found that temporal dynamics of the P. nodorum populations was strongly affected by the corresponding host populations as indicated by significant changes in genotype frequencies both at local (pair-wise comparisons, Table 4 columns 2-4) and global (multiple population comparison, Table 4 last column) levels of comparison and selection coefficients were strongly affected by host genotypes (Table 6). 2) The differences in genetic composition among pathogen populations from different hosts increased over time. Greater differences in genotype frequencies were observed among P. nodorum populations sampled from different hosts at late stages of the experiment compared to early stages of the experiment (Tables 3 and 5).
Phaeosphaeria nodorum requires 2-3 weeks to complete a cycle of asexual reproduction and the discharge of its pycnidiospores requires rain [45, 52–54]. The time intervals between the first collection and second collection were 28 days in the 2003-2004 experiment and 60 days in the 2004-2005 experiment, respectively. Using the meteorological data provided by the local weather station, we estimate that only one generation of asexual reproduction occurred between the first collection and the second collection in 2004 while two asexual generations occurred between the first collection and the second collection in 2005. The significant changes in population composition observed in our experiments indicate strong competition among pathogen genotypes and rapid adaptation to particular hosts, consistent with the hypothesis of rapid pathogen evolution in agricultural ecosystems. Elevated rates of pathogen evolution in agriculture have also been supported empirically for other plant-pathogen interactions. In Mycosphaerella graminicola, pathogen populations were collected three times during a single growing season from a susceptible host and rapid directional increases/decreases in genotype frequency were observed for all marked isolates across all replicates . Sequence analyses of plant cell wall degrading enzymes in M. graminicola  provide further evidence that the evolution of pathogens is accelerated in agricultural ecosystems. Montarry et al.  also detected a rapid change of population composition in the potato pathogen Phytophthora infestans.
The only strain carrying ToxA, SN99CH3.20a, began at a relatively high frequency (17-27%) in each host treatment but was always present at a lower frequency (0-13%) by the final 2005B sample, with an average decrease in frequency of 17%. By comparison, the two other strains (SNCH3.08a and SNCH3.10a) that had an overall decrease in frequency on each host treatment between the first and the final collections showed a decrease averaging less than 2%. A fitness cost associated with ToxA was proposed earlier to explain the observed differences in frequencies of ToxA positive strains among geographical P. nodorum populations . However, this experiment was not designed to determine whether there is a fitness cost associated with carrying ToxA and it is not clear whether any of the Swiss wheat cultivars used in this experiment carry the corresponding toxin sensitivity allele Tsn1. Therefore, we cannot conclude that the decrease in frequency of SN99CH3.20a reflected a fitness cost associated with ToxA.
The effect of host diversity and partial resistance on the evolution of P. nodorum
Both theoretical and empirical studies have demonstrated that increasing genetic diversity in host populations through deployment of cultivar mixtures offers a promising approach to control plant diseases, with the advantages of lower input costs and a reduction in ecological damage compared to use of fungicides while also providing greater yield stability [12, 22, 56–58]. Investigations of the effect of cultivar mixtures on the evolution of pathogens have been mainly theoretical [12, 33, 34]. It was hypothesized that increasing genetic diversity in host populations would retard the rate at which pathogens evolve [14, 15, 21, 59] because heterogeneity in the host population would lead to divergent selection pressure on the pathogen population . Several theoretical investigations support the hypothesis that increasing genetic diversity in host populations by using cultivar mixtures will delay the emergence of virulence against major resistance genes [22, 33, 34]. For quantitative resistance, mixing two cultivars in any proportions may reduce the final virulence attained by the pathogen population and prolong the time needed to reach the equilibrium point of highest virulence .
We determined the effect of host diversity on the population dynamics of P. nodorum by directly monitoring changes in frequencies of marked strains in a replicated field experiment. Our results support the hypothesis that increasing genetic diversity in host populations through deployment of cultivar mixtures can slow down the rate of evolution in pathogen populations. Multi-population comparisons indicated that the genetic structure of P. nodorum populations sampled from the host mixture did not change significantly over two years (Table 4 last column) and displayed the lowest variation in selection coefficients (Table 6). Although this experiment included only one host mixture, a similar evolutionary pattern was observed in field experiments involving other plant pathogens including Mycosphaerella graminicola on wheat  and Rhynchosporium secalis and Blumeria hordei on barley [38, 60].
It was postulated that partial resistance would retard the evolution of pathogens and thus increase the durability of resistance [24, 61]. A theoretical study indicated that virulence of pathogens would evolve slowly in the presence of partially resistant hosts . We tested this hypothesis in our experiments by comparing the changes in frequency of marked pathogen strains competing on susceptible and partially resistant cultivars. In a similar experiment conducted using the wheat-Mycosphaerella graminicola pathosystem, we found that the pathogen populations sampled from a partially resistant cultivar exhibited less change in genetic structure over time and smaller selection coefficients than those from a susceptible cultivar , consistent with the theoretical expectation. Results from this experiment also support the hypothesis. We detected less significant changes in P. nodorum populations sampled from the partially resistant cultivars Tamaro and Runal than from the susceptible cultivar Tirone and resistant cultivar Levis (Table 4). Multi-population comparisons revealed that the changes in genetic structure of the P. nodorum populations sampled from the two partially resistant cultivars were significant at the 5% level over the two-years of the experiment, while the differences were significant at the 1% level on the resistant and susceptible cultivars. This result suggests that directional selection also occurs in pathogen populations infecting partially resistant cultivars albeit at a slower pace, leading to the erosion of resistance . The adaptation to partial resistance was also observed in field experiments with the barley scald pathogen Rhynchosporium secalis  and the potato late blight pathogen Phytophthora infestans . Some theoretical analyses of host-pathogen co-evolution suggest that hosts with partial resistance can select for increased virulence in pathogen populations . Our findings agree with these predictions.
Differential selection on P. nodorumstrains during parasitic and saprophytic phases
It is hypothesized that saprophytic and parasitic phases of pathogen life cycles may select for different pathogen traits. During the parasitic phase, pathogen strains with a high capacity to exploit the host may have a selective advantage if they are able to produce greater numbers of viable offspring compared to strains with a lower capacity. But the traits that are favored during the parasitic phase of the life cycle may be selected against during saprophytic phases of the life cycle when living hosts are not available [7, 23]. The trade-offs that occur between parasitic and saprophytic phases of the life cycle may prevent or delay the emergence of high levels of pathogen virulence.
We found some evidence for differential selection between parasitic and saprophytic phases of the life cycle in the P. nodorum populations sampled from Runal, suggesting a fitness cost associated with high virulence during the saprophytic phase [10, 62]. On Runal, a significant difference in P. nodorum genetic structure was found between 2004A and 2004B but not between 2004A and 2005A, suggesting that selection occurring during the parasitic phase might be offset by selection that occurred during the saprophytic phase. The frequency distribution of isolate SN99CH2.04 also suggested that differential selection might occur between parasitic and saprophytic phases of the pathogen life cycle. This isolate increased in frequency on four of the five treatments (the mixture was the exception) during the parasitic phase of the disease cycle but decreased in frequency during the saprophytic phase on all treatments except the cultivar mixture (Figure 3). This finding indicates that this strain may exhibit higher relative fitness during the parasitic phase on the majority of living host tissue but lower competitive ability during the saprophytic phase on the dead host tissue. Using the same experimental approach, Abang et al.  also found that some isolates increased in frequency during the parasitic phase but decreased in frequency during the saprophytic phase in the barley pathogen Rhynchosporium secalis.
The lack of evidence for differential selection between parasitic and saprophytic phases in other hosts and isolates may be partially attributed to our sampling strategy. The 2005A collection was made from infected plants several months after the application of the wheat stubble inoculum. Thus at least one cycle of parasitic competition had likely occurred among the pathogen strains before this collection was made. If there was differential selection between the parasitic and saprophytic phases, selection for traits involved in establishment and reproduction during the initiation of the 2005A epidemics may have partially offset selection for traits involved in saprophytic competition. Further experiments with an additional population sample drawn at the beginning of the second cycle of the parasitic phase will be necessary to confirm this hypothesis
Only isolates derived from asexual reproduction of the nine inoculants were used to calculate selection coefficients and determine the effects of host genotypes on the population genetic structure of P. nodorum. Isolates having genotypes different from the nine inoculated strains (called novel isolates) were excluded. Because the contribution of mutation to the formation of new genotypes is expected to be trivial within the time scale of this experiment, we believe these novel isolates originated either via immigration from outside of the experimental plots or by recombination among inoculants and/or immigrants within the experimental plots (details in 39). The current paper focuses on the influence of host genotypes and diversity on clonal competition and we believe that excluding these novel isolates did not affect our interpretations.
We used both pair-wise and multiple population comparisons to evaluate the effects of evolutionary time (different sampling points) and host genotypes on the population dynamics of P. nodorum. In the multiple population comparisons, pathogen populations from different sampling points or hosts were considered simultaneously in a single analysis. This approach is useful to determine the overall pattern of evolutionary change in pathogen populations over hosts (e.g. last column in Table 4) but cannot be used to determine the sequential change in population structure over time within a host, for example whether the population genetic structure between the first (2004A) and the second (2004B) collection differs more than that between the first and the last (2005B) collection. For the latter case, we adopted pair-wise comparisons that included a Bonferroni correction. We believe that combining these approaches was necessary to achieve a comprehensive analysis of the data.
Understanding the evolutionary response of pathogen populations to host diversity and environmental changes (such as over-wintering or over-summering between growing seasons) is important for disease management. Many studies on host-pathogen interactions have focused on the development of mathematical models [12, 33, 63] to predict pathogen evolution in response to different strategies of resistance gene deployment [13, 61, 64–66]. Very few empirical studies have been conducted to test these theoretical models in agricultural ecosystems. Here, we present empirical evidence that strong selection occurs during both parasitic and saprophytic phases of the disease cycle. Evolution during the parasitic phase occurred most slowly on the cultivar mixture. The same result was also reported in similar experiments conducted with the wheat pathogen Mycosphaerella graminicola  and the barley pathogen Rhynchosporium secalis , suggesting that the observed pattern of evolution may be applicable for other splash-dispersed pathogens on cereals.
A two-year mark-release-recapture experiment was conducted at the Agroscope Changins-Wädenswil research center in Changins, Switzerland during the 2003-2004 winter wheat growing season on field allotment-34-North and the 2004-2005 winter wheat growing season on field allotment-35-North. Both fields were grown with a permanent meadow for at least three years prior to the mark-release-recapture experiment. Four commercial Swiss wheat cultivars, namely Levis, Runal, Tamaro and Tirone, were used in these experiments. The varieties differed in quantitative resistance to P. nodorum leaf infection according to disease assessments conducted at Changins between 2001 and 2002 . Cultivars Levis, Tamaro and Runal are partially resistant to P. nodorum leaf blotch with their levels of resistance decreasing in that order. Cultivar Tirone is susceptible to leaf blotch. Cultivar Levis is partially resistant on the leaves but not on glumes. The four cultivars and a 1:1 mixture of cultivars Runal and Tamaro (5 host treatments in total) were planted in a randomized complete block design (RCBD) with three replications. The field plots were 1.5 m in width and 4.5 m in length. Each wheat plot was surrounded by four equal-sized plots planted with the highly resistant triticale variety Tridel. The experiment was planted on 5 October 2003 in the first year and on 17 October 2004 in the second year using commercial seeds treated with the fungicide Coral (2.38% difenoconazole and 2.38% fludioxonil, 2 ml/kg seeds).
Nine P. nodorum isolates collected in 1999 from naturally infected wheat fields near Bern, Switzerland were chosen as inoculants for the 2003-2004 experiment. Each of the isolates had distinct multi-locus haplotypes when assayed with seven single-locus RFLP markers , ten polymorphic EST-derived microsatellite markers and one minisatellite marker . Only one of the nine isolates (SN99CH3.20a) carried the ToxA gene  that encodes a host specific toxin. After completing the field experiments it was discovered that one isolate (SN99CH3.19a) had been replaced by a contaminant of unknown origin, hereafter called C1. Following removal from long-term storage at -80°C, the isolates were first grown on Yeast Maltose Agar (YMA, yeast 4 gl-1, maltose 4 gl-1, sucrose 4 gl-1, agar 10 gl-1) at 21°C for ten days and then transferred to 1000 ml flasks containing 300 g of sterilized wheat kernels (cultivar Arina) in a dark incubator at 4°C. Three months later the infected wheat kernels were harvested and ground to a powder using a gristmill. The powdered kernels were mixed with distilled water and the spore suspension was filtered through cheese cloth and glass-wool. The spore suspension from each isolate was adjusted to 106 spores per ml using a hemacytometer and mixed in equal proportions. A surfactant (Tween 20) was added to the spore suspension at the rate of one drop per 50 ml. The aqueous spore suspension was applied onto disease-free wheat seedlings at growth stage 31 on 11 May 2004. Each field plot was sprayed with 500 ml of the calibrated spore suspension. To maximize the humidity and increase the probability of infection, inoculations were carried out in the late afternoon on a cloudy day and the inoculated seedlings were covered with a plastic tarp for 24 hours.
The source of primary inoculum in the 2004-2005 experiment was the infected straw and other plant debris saved from the first year's experiment. After harvesting the grain at the end of July 2004, the straw and other plant debris in each plot were collected and stored separately in burlap potato sacks for 3 months. The sacks were stored in a dry, dark and cool room to allow the development of the saprophytic phase without the risk of excessive moulding. At the beginning of tillering (Zadoks stage 13 to 21, 69), the straw was applied onto the corresponding host plots.
A total of four fungal collections were made across the two growing seasons. The first collection, hereafter called 2004A, was made on 4 June 2004 from the third or fourth full leaf  at three weeks after the artificial inoculation. The second collection, hereafter called 2004B, was made on 2 July 2004 from flag leaves. The third collection, hereafter called 2005A, was made on 11 April 2005 from the second true leaf and the last collection, hereafter called 2005B, was made on 10 June 2005 from the third true leaf. For each collection, 30 to 40 leaves were collected from each inoculated plot at intervals of approximately 20 cm along transects within the inner rows of the field plots. In most cases, only one isolate was made from each infected leaf. Because many lesions did not contain pycnidia, the total number of isolations made was much lower than the number of wheat leaves collected. For some collections with very low levels of infection, two isolations were made from the same leaf. In these cases, isolations were made from clearly separated lesions to minimize the probability of sampling isolates from the same infection. Our earlier work showed that P. nodorum isolations made from different lesions within an infected leaf usually contain different genotypes, suggesting they originate from different infection events .
DNA extraction and microsatellite data collection
DNA was extracted from each isolate using the DNeasy Plant Mini DNA extraction kit (Qiagen GmbH, Germany) according to the specifications of the manufacturer. The genotype of each isolate was determined using the same ten EST-derived microsatellite markers (SN1, SN3, SN5, SN11, SN15, SN16, SN17, SN21, SN22, and SN23) and one minisatellite marker (SN8) used to mark the nine inoculants. Multiplexed polymerase chain reactions (PCR) were carried out with fluorescently labeled primers using the same conditions described previously . The PCR products were first cooled on ice for 2 min and then separated on an ABI PRISM 3100 sequencer according to the manufacturer's instructions (Applied Biosystems). Fragment sizes were estimated and alleles were assigned using the program GENESCAN 3.7 (Applied Biosystems).
The multilocus haplotype (MLHT) for each isolate was formed by joining the alleles at each of the 11 marker loci. Isolates with the same MLHT were considered to be clones, the products of asexual reproduction for a particular genotype. Because fungal collections were found to consist of both inoculated and novel isolates (novel isolates are defined as isolates having MLHTs different from any of the nine inoculated isolates, 39), only isolates derived from asexual reproduction of the nine inoculants were used to calculate selection coefficients and determine the effects of host genotypes on the population genetic structure of P. nodorum. The novel isolates, which originated either via immigration from outside of the experimental fields or by recombination among inoculated isolates and immigrants, were considered in a separate publication .
Pair-wise and multiple population comparisons in genotype frequencies were performed using contingency χ2 tests to detect differences in pathogen populations sampled from different hosts or sampling times . Pair-wise comparisons were corrected using a sequential Bonferroni procedure as described previously . Multiple population comparisons were conducted by using all populations sampled from the same host across different sampling points simultaneously, generating a c x l contingency table, where c is the number of sampling points and l is the number of genotypes detected. The association between sampling time and differences in genotype frequency among populations was evaluated using simple linear correlation. In this analysis, differences in genotype frequency among populations were measured by GST [73, 74]
This estimate of selection coefficient measures the overall fitness of a genotype relative to the most-fit isolate in a population during the entire life cycle of the pathogen, taking into account its ability to infect, colonize, reproduce and spread . For example, a genotype with a selection coefficient of 0.30 has 30% lower fitness than the most-fit genotype. Selection coefficient is different from selection intensity, a term used mainly in breeding and quantitative genetics to quantify a potential evolutionary gain after selecting a pool of parents from a variable population. To make more robust estimates, selection coefficients were estimated for each host treatment by pooling together data from different replications using only five of the nine inoculants for the 2004A and 2004B collections. Isolate SN99CH3.23a was not recaptured during the course of the experiment and the other three inoculants were recovered at frequencies too low to make meaningful estimates of selection coefficients. Selection coefficients were not estimated for the 2005 collections due to small sample sizes (9-18 isolates/host treatment) remaining after the novel isolates were excluded from 2005B. The means and standard deviations of the selection coefficients were generated based on 100 resamples of the original genotype frequencies using the Excel add-in PopTools 2.7 (CSIRO, Australia). Tukey's significant differences implemented in SYSTAT were used to compare selection coefficients among the five inoculants.
This research was supported by the Swiss Federal Institute of Technology Grant TH-49a/02-1t. SSR data were collected using facilities of the Genetic Diversity Center at ETH Zurich. We thank S. Kellenberger for the field work, D. Gobbin, M. Lutz and P. Zaffarano for help with the data analysis, and V. Martinez, C. Phan, S. Seeholzer, and M. Zala for technical assistance.
- Anderson JB: Evolution of antifungal-drug resistance: mechanisms and pathogen fitness. Nature Rev Microbio. 2005, 3: 547-556.View ArticleGoogle Scholar
- McDonald BA, Linde C: Pathogen population genetics, evolutionary potential and durable resistance. Annu Rev Phytopathol. 2002, 40: 359-379.View ArticleGoogle Scholar
- Stahl EA, Bishop JG: Plant-pathogen arms races at the molecular level. Curr Opinion Plant Biol. 2000, 3: 299-304.View ArticleGoogle Scholar
- Laine AL: Evolution of host resistance: looking for coevolutionary hotspots at small spatial scales. Proc R Soc B. 2006, 273: 267-273.View ArticlePubMedGoogle Scholar
- Stukenbrock EH, McDonald BA: The origins of plant pathogens in agro-ecosystems. Annu Rev Phytopathol. 2008, 46: 75-100.View ArticlePubMedGoogle Scholar
- Brunner PC, Keller N, McDonald BA: Wheat domestication accelerated evolution and triggered positive selection in the β-Xylosidase enzyme of Mycosphaerella graminicola. PLoS One. 2009, 4: e7884-View ArticlePubMedPubMed CentralGoogle Scholar
- Kiyosawa S: Genetic and epidemiological modeling of breakdown of plant disease resistance. Annu Rev Phytopathol. 1982, 20: 93-117.View ArticleGoogle Scholar
- Goodwin SB, Cohen BA, Fry WE: Panglobal distribution of a single clonal linage of the Irish potato famine fungus. Proc Natl Acad Sci USA. 1994, 91: 11591-11595.View ArticlePubMedPubMed CentralGoogle Scholar
- Gandon S, Michalakis Y: Evolution of parasite virulence against qualitative and quantitative host resistance. Proc R Soc Lond B. 2000, 267: 985-990.View ArticleGoogle Scholar
- Bergelson J, Dwyer G, Emerson JJ: Models and data on plant-enemy coevolution. Annu Rev Genet. 2001, 35: 469-499.View ArticlePubMedGoogle Scholar
- Montarry J, Corbiere R, Lesueur S, Glais I, Andrivon D: Does selection by resistant hosts trigger local adaptation in plant-pathogen systems?. J Evol Biol. 2006, 19: 522-531.View ArticlePubMedGoogle Scholar
- Marshall B, Newton AC, Zhan J: Quantitative evolution of aggressiveness of powdery mildew under two-cultivar barley mixtures. Plant Pathol. 2009, 58: 378-388.View ArticleGoogle Scholar
- Mundt CC: Use of multiline cultivars and cultivar mixtures for disease management. Annu Rev Phytopathol. 2002, 40: 381-410.View ArticlePubMedGoogle Scholar
- Burdon JJ: Disease and Plant Population Biology. 1987, Cambridge University Press, CambridgeGoogle Scholar
- Dileone JA, Mundt CC: Effect of wheat cultivar mixtures on populations of Puccinia striiformis races. Plant Pathol. 1994, 43: 917-930.View ArticleGoogle Scholar
- Higashi T: Genetic studies on field resistance of rice to blast disease. Bull Tohoku Natl Agri Exp Stn. 1995, 90: 19-75.Google Scholar
- Ebert D, Hamilton WD: Sex against virulence: the coevolution of parasitic diseases. Trends Ecol Evol. 1996, 11: 79-82.View ArticlePubMedGoogle Scholar
- Zhan J, Mundt CC, McDonald BA: Using RFLPs to assess temporal variation and estimate the number of ascospores that initiate epidemics in field populations of Mycosphaerella graminicola. Phytopathology. 2001, 91: 1011-1017.View ArticlePubMedGoogle Scholar
- Andrivon D, Pilet F, Montarry J, Hafidi M, Corbiere R, Achbani EH, Pelle R, Ellisseche D: Adaptation of Phytophthora infestans to partial resistance in potato: evidence from French and Moroccan populations. Phytopathology. 2007, 97: 338-343.View ArticlePubMedGoogle Scholar
- Zhan J, Mundt CC, McDonald BA: Sexual reproduction facilitates the adaptation of parasites to antagonistic host environment: evidence from field experiment with wheat- Mycosphaerella graminicola system. Intl J Parasitol. 2007, 37: 861-870.View ArticleGoogle Scholar
- Parlevliet JE: Stabilizing selection in crop pathosystems: an empty concept or reality?. Euphytica. 1981, 30: 259-269.View ArticleGoogle Scholar
- Lannou C, Mundt CC: Evolution of a pathogen population in host mixtures: simple race - complex race competition. Plant Pathol. 1996, 45: 440-453.View ArticleGoogle Scholar
- Vanderplank JE: Disease Resistance in Plants. 1968, Academic Press, New YorkGoogle Scholar
- Simons MD: Polygenic resistance to plant disease and its use in breeding resistant cultivars. J Environ Qual. 1972, 1: 232-240.View ArticleGoogle Scholar
- Robinson RA: Plant Pathosystems. 1976, Springer, Berlin, GermanyView ArticleGoogle Scholar
- Maruyama K, Kikuchi F, Yokoo M: Gene analysis of field resistance to rice blast (Pyricularia oryzae) in Rikuto Norin Mochi 4 and its use for breeding. Bull Natl Inst Agr Sci. 1983, 35: 1-31.Google Scholar
- Parlevliet JE: What is durable resistance: a general outline. Durability of Disease Resistance. Edited by: Jacobs T, Parlevliet JE. 1993, Kluwer Academic Publishers, LondonGoogle Scholar
- Clifford BC, Clothier RB: Physiologic specialization of Puccinia hordeion barley hosts with non-hypersensitive resistance. Trans Br Mycol Soc. 1974, 63: 421-430.View ArticleGoogle Scholar
- Kolmer JA, Leonard KJ: Genetic selection and adaptation of Cochliobolus heterostrophusto corn hosts with partial resistance. Phytopathology. 1986, 76: 774-777.View ArticleGoogle Scholar
- Pink DAC, Lot H, Johnson R: Novel pathotypes of lettuce mosaic virus-Breakdown of a durable resistance?. Euphytica. 1992, 63: 169-174.View ArticleGoogle Scholar
- Schouten HJ, Benires JE: Durability of resistance to Globodera pallida. I. Changes in pathogenicity, virulence, and aggressiveness during reproduction on partially resistant potato cultivars. Phytopathology. 1997, 87: 862-867.View ArticlePubMedGoogle Scholar
- Zhan J, Mundt CC, Hoffer ME, McDonald BA: Local adaptation and effect of host genotype on the rate of pathogen evolution: an experimental test in a plant pathosystem. J Evol Biol. 2002, 15: 634-647.View ArticleGoogle Scholar
- Lannou C: Intrapathotype diversity for aggressiveness and pathogen evolution in cultivar mixtures. Phytopathology. 2001, 91: 500-510.View ArticlePubMedGoogle Scholar
- Segarra J: Stable polymorphisms in a two-locus gene-for-gene system. Phytopathology. 2005, 95: 728-736.View ArticlePubMedGoogle Scholar
- Goss EM, Larsen M, Chastagner GA, Givens DR, Grunwald NJ: Population genetic analysis infers migration pathway of Phytophthora ramorumin US nurseries. PLoS Pathog. 2009, 5: e1000583-View ArticlePubMedPubMed CentralGoogle Scholar
- Zhan J, Mundt CC, McDonald BA: Measuring immigration and sexual reproduction in field populations of Mycosphaerella graminicola. Phytopathology. 1998, 88: 1330-1337.View ArticlePubMedGoogle Scholar
- Zhan J, Mundt CC, McDonald BA: Estimating rates of recombination and migration in populations of plant pathogens--a reply. Phytopathology. 2000, 90: 324-326.View ArticlePubMedGoogle Scholar
- Abang MM, Baum M, Ceccarelli S, Grando S, Linde C, Yahyaoui A, Zhan J, McDonald BA: Differential selection on Rhynchosporium secalisduring the parasitic and saprophytic phases in the barley scald disease cycle. Phytopathology. 2006, 96: 1214-1222.View ArticlePubMedGoogle Scholar
- Sommerhalder RJ, McDonald BA, Mascher F, Zhan J: Sexual recombinants make a significant contribution to epidemics caused by the wheat pathogen Phaeosphaeria nodorum. Phytopathology. 2010, 100: 855-862.View ArticlePubMedGoogle Scholar
- Bennett RS, Milgroom MG, Sainudiin R, Cunfer BM, Bergstrom GC: Relative contribution of seed-transmitted inoculum to foliar populations of Phaeosphaeria nodorum. Phytopathology. 2007, 97: 584-591.View ArticlePubMedGoogle Scholar
- Eyal Z: The Septoria tritici and Stagonospora nodorum blotch diseases of wheat. Eur J Plant Pathol. 1999, 105: 629-641.View ArticleGoogle Scholar
- Cowger C, Silva-Rojas HV: Frequency of Phaeosphaeria nodorum, the sexual stage of Stagonospora nodorum, on winter wheat in North Carolina. Phytopathology. 2006, 96: 860-866.View ArticlePubMedGoogle Scholar
- Sommerhalder RJ, McDonald BA, Zhan J: The frequencies and spatial distribution of mating types in Stagonospora nodorum are consistent with recurring sexual reproduction. Phytopathology. 2006, 96: 234-239.View ArticlePubMedGoogle Scholar
- Shaner G: Effect of environment on fungal leaf blights of small grains. Annu Rev Phytopathol. 1981, 19: 263-296.View ArticleGoogle Scholar
- Bathgate JA, Loughman R: Ascospores are a source of inoculum of Phaeosphaeria nodorum, P. avenaria f. sp. avenaria and Mycosphaerella graminicola in Western Australia. Aust Plant Path. 2001, 30: 317-322.View ArticleGoogle Scholar
- Griffiths E, Ao HC: Dispersal of Septoria nodorum spores and spread of glume blotch in the field. Trans Br Mycol Soc. 1976, 67: 413-418.View ArticleGoogle Scholar
- Brennan RM, Fitt BDL, Taylor GS, Colhoun J: Dispersal of Septoria nodorum pycnidiospores by simulated raindrops in still air. Phytopathology. 1985, 112: 281-290.View ArticleGoogle Scholar
- Keller SM, McDermott JM, Pettway RE, Wolfe MS, McDonald BA: Gene flow and sexual reproduction in the wheat glume blotch pathogen Phaeosphaeria nodorum, (anamorph Stagonospora nodorum). Phytopathology. 1997, 87: 353-358.View ArticlePubMedGoogle Scholar
- Arseniuk E, Góral T, Scharen AL: Seasonal patterns of spore dispersal of Phaeosphaeria spp. and Stagonospora spp. Plant Dis. 1998, 82: 187-194.View ArticleGoogle Scholar
- Weber GF: Septoria diseases of wheat. Phytopathology. 1922, 12: 537-585.Google Scholar
- von Wechmar MB: Investigation on the survival of Septoria nodorum Berk. on crop residues. S African J Agri Sci. 1966, 9: 93-100.Google Scholar
- Shearer BL, Zadoks JC: The latent period of Septoria nodorum in wheat. 1. The effect of temperature and moisture treatments under controlled conditions. Neth J Plant Pathol. 1972, 78: 231-241.View ArticleGoogle Scholar
- Shearer BL, Zadoks JC: The latent period of Septoria nodorum in wheat. 2. The effect of temperature and moisture treatments under field conditions. Neth J Plant Pathol. 1973, 80: 48-60.View ArticleGoogle Scholar
- Solomon PS, Lowe RGT, Ian KC, Waters ODC, Oliver RP: Stagonospora nodorum: cause of stagonospora nodorum blotch of wheat. Mol Plant Pathol. 2006, 7: 147-156.View ArticlePubMedGoogle Scholar
- Stukenbrock EH, McDonald BA: Geographic variation and positive diversifying selection in the host specific toxin SnToxA. Mol Plant Pathol. 2007, 8: 321-332.View ArticlePubMedGoogle Scholar
- Zhu YY, Chen HR, Fan JH, Wang YY, Li Y, Chen JB, Fan JX, Yang SS, Hu LP, Leung H, Mew TW, Teng PS, Wang ZH, Mundt CC: Genetic diversity and disease control in rice. Nature. 2000, 406: 718-722.View ArticlePubMedGoogle Scholar
- Woldeamlak A, Grando S, Maatougui M, Ceccarelli S: Hanfets, a barley and wheat mixture in Eritrea: yield, stability and farmer preferences. Field Crops Res. 2008, 109: 50-56.View ArticleGoogle Scholar
- Frankow-Lindberg BE, Halling M, Hoglind M, Forkman J: Yield and stability of yield of single- and multi-clover grass-clover swards in two contrasting temperate environments. Grass Forage Sci. 2009, 64: 236-245.View ArticleGoogle Scholar
- Huang R, Kranz J, Welz HG: Selection of pathotypes of Erysiphe graminis f.sp. hordei in pure and mixed stands of spring barley. Plant Pathol. 1994, 43: 458-470.View ArticleGoogle Scholar
- Chin KM, Wolfe MS: Selection on Erysiphe graminis in pure and mixed stands of barley. Plant Pathol. 1984, 33: 535-546.View ArticleGoogle Scholar
- Parlevliet JE: Durability of resistance against fungal, bacterial and viral pathogens; present situation. Euphytica. 2002, 124: 147-156.View ArticleGoogle Scholar
- Bahri B, Kaltz O, Leconte M, de Vallavieille-Pope C, Enjalbert J: Tracking costs of virulence in natural populations of the wheat pathogen, Puccinia striiformis f.sp.tritici. BMC Evol Biol. 2009, 9: 26-View ArticlePubMedPubMed CentralGoogle Scholar
- Barrett JA: Pathogen evolution in multilines and variety mixtures. Zeitschriften für Pflanzenkrankheiten und Pflanzenschutz. 1980, 87: 383-396.Google Scholar
- Paillard S, Goldringer I, Enjalbert J, Trottet M, David J, de Vallavieille-Pope C, Brabant P: Evolution of resistance against powdery mildew in winter wheat populations conducted under dynamic management. II. Adult plant resistance. Theor Appl Genet. 2000, 101: 457-462.View ArticleGoogle Scholar
- De Meaux J, Mitchell-Olds T: Evolution of plant resistance at the molecular level: ecological context of species interactions. Heredity. 2003, 91: 345-352.View ArticlePubMedGoogle Scholar
- Long J, Holland B, Munkvold GP, Jannink JL: Response to selection for partial resistance to crown rust in oat. Crop Sci. 2006, 46: 1260-1265.View ArticleGoogle Scholar
- Collaud JF, Schwärzel R, Bertossa M, Menzi M, Anders M: Variétés de céréales recommandées par l'interprofession pour la récolte 2003. Revue suisse d'agriculture. 2002, 34: (insert)Google Scholar
- Stukenbrock HE, Banke S, Zala M, McDonald BA, Oliver RP: Isolation and characterization of EST-derived microsatellite loci from the fungal wheat pathogen Stagonospora nodorum. Mol Ecol Notes. 2005, 5: 931-933.View ArticleGoogle Scholar
- Zadoks JC, Chang TT, Conzak CF: A decimal code for the growth stages of cereals. Weed Res. 1974, 14: 415-421.View ArticleGoogle Scholar
- McDonald BA, Miles J, Nelson LR, Pettway RE: Genetic variability in nuclear DNA in field populations of Stagonospora nodorum. Phytopathology. 1994, 84: 250-255.View ArticleGoogle Scholar
- Everitt BS: The Analysis of Contingency Tables. 1977, Wiley, New YorkView ArticleGoogle Scholar
- Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225.View ArticleGoogle Scholar
- Nei M: F-statistics and analysis of gene diversity in subdivided populations. Ann Hum Genet. 1977, 41: 225-233.View ArticlePubMedGoogle Scholar
- Selander RK, Caugant DA, Ochman H, Musser JM, Gilmour MN, Whittam TS: Methods of multilocus enzyme electrophoresis for bacterial population genetics and systematics. Appl Environ Microbiol. 1986, 51: 873-884.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.