Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2006, p. 185-191, Vol. 72, No. 1
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.1.185-191.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Phytopathologie et Méthodologies de la Détection, INRA, F-78026 Versailles, France,1 Protection des Plantes, INA-PG, F-78850 Thiverval-Grignon, France,2 Rothamsted Research, Harpenden AL5 2JQ, United Kingdom3
Received 7 April 2005/ Accepted 3 October 2005
|
|
|---|
|
|
|---|
The epidemiology of the stem canker disease has been extensively investigated (19, 39) since the first descriptions of the causal fungus by Tode in 1751 and Desmazières in 1849 (34). In Europe, OSR is grown as a winter crop and L. maculans produces pseudothecia in the autumn following sexual outcrosses on infected debris from the previous growing season. Airborne ascospores discharged from these pseudothecia infect leaves of the recently emerged crop, resulting in the development of leaf lesions (phoma leaf spots) that contain the asexual fruiting bodies, pycnidia. Asexual pycnidiospores can cause secondary infections, which are thought to be important in Australia but rare in Canada and Europe, where the disease is thought to be monocyclic (39). From the initial necrotrophic phase when leaf lesions occur, the fungus enters a symptomless biotrophic phase when, during winter and early spring, the fungus grows systemically in intercellular spaces down the petiole toward the crown at the stem base. During late spring and early summer in Europe, L. maculans again becomes necrotrophic and causes necrosis of the crown tissues that may result in the plant lodging prior to harvest (20, 21).
Strategies for stem canker disease management include cultural practices such as crop rotation, isolation of the crops from infected stubble of the previous growing season crops, and stubble management. Some control is achieved through the use of fungicides but, at present, disease control relies mainly on the use of disease-resistant cultivars. There are two types of disease resistance in B. napus: polygenic quantitative resistance (29, 30) or specific resistance, which involves major resistance gene(s) (19, 39). Nine avirulence genes (AvrLm genes), namely, AvrLm1 to AvrLm9, have been identified in L. maculans (2, 6-8), and the corresponding nine major resistance genes (Rlm1 to Rlm9) have been described in Brassica sp. (13), as expected in a host-pathogen system with gene-for-gene interactions (17). The Rlm genes effectively control the disease as long as the corresponding avirulent allele (AvrLm) dominates the pathogen population (3, 9, 31, 35). However, large-scale cropping of cultivars with a single effective Rlm gene results in the rapid breakdown of major gene resistance within a few years after commercial release of the cultivars due to the adaptation of the L. maculans field populations. Such breakdown of resistance was documented in France for Rlm1 (35) and in Australia for the breakdown of the "Surpass 400" resistance (23).
Our knowledge of genetic variation in L. maculans populations and the role of the evolutionary processes such as mutation, migration, genetic drift, selection, and/or recombination in these rapid shifts in the field population is still very limited. Previous studies of genetic variation in L. maculans populations from Canada and Australia found high levels of genetic variation in field populations of L. maculans (10, 11, 24), which is consistent with the hypothesis that sexual reproduction is an important part of the life cycle of L. maculans. However, reports of the distribution of genetic variation among and within field populations have been inconsistent (10, 11, 24). To date, there have been no studies of the genetic structure of L. maculans field populations in France and other countries of Western Europe. In these regions, previous work has focused mainly on field population race structure and has used virulence markers to assess genetic diversity in the fungus (5, 9, 22, 35).
Our objective here was to assess genetic variation in French populations of L. maculans with minisatellite loci (15), mating type, and virulence markers (8). We tested the hypotheses that the field populations were randomly mating, that populations from different fields were genetically separable, and that epidemics are initiated by ascospores. The present study provides the first estimates of genetic structure of field populations of L. maculans from France and suggests that genetic structure of populations from Western Europe is similar to that in Australia.
|
|
|---|
Pathogenicity tests.
Races of the L. maculans isolates were determined by inoculating the 401 isolates onto a B. napus differential set, comprising fixed lines or commercial cultivars with as few Rlm alleles as possible (8). Cotyledons of the fully susceptible cv. Westar (no known Rlm allele), the cvs. Columbus (Rlm1 and Rlm3), Bristol (Rlm2 and Rlm9), 22-1-1 (Rlm3), Jet Neuf (Rlm4), 23-2-1 (Rlm7), Falcon-MX (Rlm4 and Rlm6), and Samourai-MX (Rlm1 and Rlm6) or Darmor-MX (Rlm6) were inoculated according to established protocols (2, 7). At 14 to 27 days after inoculation, symptoms were scored from 10 to 12 plants by using the IMASCORE rating scale comprising six infection classes (IC1 to IC6), where IC1 to IC3 corresponded to avirulent isolates (AvrLm) and IC4 to IC6 corresponded to virulent isolates (avrLm) (7).
Mating type determination.
A multiplex PCR, developed to rapidly characterize the mating type of L. maculans isolates (12), was used to determine the distribution of the Mat1-1 and Mat1-2 alleles within the 401 field isolates collected. The distribution of mating type alleles was analyzed with a
2 test.
Minisatellite analyses.
We evaluated six polymorphic minisatellite loci (MinLm1 to MinLm6), all on different chromosomes in L. maculans (4, 15). New specific primers based on the published sequences (accession numbers AJ621802 and AJ621805, respectively) were designed for MinLm3 and MinLm5 by using Oligo v.5.0 software (Molecular Biology Insights, Inc., Cascade, Iowa). Primers MinLm3-ULG (5'-GGCTCGGTCGGTTAGTTA-3') and MinLm3-LLG (5'-AATGATGTACAGGACGGGATTT-3') or primers MinLm5-ULG (5'-GCCGCCCGCCGCCTTACC-3') and MinLm5-LLG (5'-GAGCTCCTGCGCCACAGTG-3') hybridized to the flanks of the minisatellite locus and allow amplification of the repeated motifs only. MinLm2, MinLm4, and MinLm6 were amplified by using primers that had been previously described (15).
For PCR, genomic DNA was extracted from pycnidiospores by using the DNeasy 96 plant kit (QIAGEN S.A., Courtaboeuf, France) as previously described (4). PCR amplification was done in a total volume of 15 µl containing 0.2 µM of each deoxynucleotide triphosphate, 1.2 µM of each primer, 0.6 U of Taq DNA polymerase (Qbiogen, Illkirch, France), 1.5 µl of a 10x reaction buffer supplied with the enzyme, and 10 to 30 ng of genomic DNA. PCR amplifications were done in an Eppendorf Mastercycler EP Gradient thermocycler (Eppendorf, Le Pecq, France), with 30 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 80 s, with a final extension at 72°C for 5 min. The cycling conditions were the same for the amplification of all minisatellite loci. The PCR products were separated by electrophoresis on a 2.5% SeaKem LE agarose gel (FMC, Rockland, Maine) and visualized by staining with ethidium bromide and UV illumination. Allele sizes were determined by using Quantity One 1-D Analysis software (Bio-Rad, Marnes-la-Coquette, France) and comparing the bands for the alleles with those of the 1-kb ladder (Invitrogen, Cergy Pontoise, France) and with PCR products obtained after amplification of alleles of known size for each minisatellite. Each allele identified at the six minisatellite loci was previously sequenced from at least two isolates to determine the exact allele size and to confirm that the size polymorphism observed after electrophoresis was due to variability in the number of tandemly repeated core motifs (4, 15). Sequencing of the same allele amplified from at least three different individuals also confirmed that sequences of the same size are homologous and size homoplasy did not occur (4, 15). The data were scored as the number of repeat units for each minisatellite.
Statistical analysis.
For statistical analyses, the 401 isolates were divided into a "leaf-scale group" comprising 50 isolates collected from one leaf in Oucques and 43 isolates collected from five leaves in Versailles and a "field-scale group" comprising 308 isolates. The field scale group was structured into four field populations, each of which was further subdivided into four subpopulations corresponding to four 2-m2 plots.
Genetic diversity, which was studied separately for each field population and for all fields sampled, was estimated from genotype and gene diversity. Genotypes of isolates in both the field- and leaf-scale groups were determined by combining information from the alleles at minisatellite loci and at the mating type locus. Genotype diversity (G) was then estimated as the percentage of unique genotypes obtained in each sample unit, in each field population, and in the overall sample units. A discriminatory index (DI), based on Simpson's index of diversity, was used to estimate the discriminatory power of the six minisatellites. DI measures the probability that two individuals randomly selected from a small sample will have the same haplotype. DI is determined by the number and relative frequencies of the different haplotypes and is calculated as follows:
![]() |
POPGENE version 1.31 software (40) was used to compute allele frequency, Nei's gene diversity (h) (27) for each population and Nei's measures of genetic identity (I) between populations over all loci (28). The mean gene diversity (H) was the mean of h over loci within each population.
Gametic linkage disequilibria were computed by using GENEPOP version 3.3 software (32) for each pair of loci within each of the four field populations and in the overall field samples. GENEPOP performs the Fisher exact test using a Markov chain to test for gametic disequilibrium between two loci. The null hypothesis H0 was: "genotypes at one locus are independent from genotypes at the other locus." The significance of departure from linkage equilibrium was assessed by using Bonferroni adjusted P values (33). Linkage disequilibrium also was estimated by pairwise comparison of haplotypes using the index of multilocus gametic disequilibrium
d (1), which tests to what extent individuals that are identical (or different) at one locus are more likely to be identical (or different) at the other loci. To determine whether an observation deviates significantly from the null hypothesis of random mating (no linkage disequilibrium among loci), the observed value was compared to the results of 500 randomized data sets. The program MULTILOCUS version 1.3 (1) was used for both estimation and randomization.
Heterogeneity of allele frequencies among field populations was investigated for each locus with contingency tables analysis and the Fisher exact test as implemented in the GENEPOP program. The null hypothesis tested was H0: "allele frequencies are homogeneous across populations." Genetic structure was analyzed with Nei's coefficient of population differentiation (GST), which indicates the proportion of the total genetic variation attributable to population differentiation. Analyses of molecular variance (AMOVA) were computed with Arlequin version 2.000 software (37) and used to estimate how genetic variability is partitioned among fields, among plots within fields, and within plots.
statistics associated with these components of variance (
CT,
SC, and
ST, respectively) were estimated and their significance tested by using a nonparametric permutation approach (16).
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Allelic frequencies and gene diversities at two avirulence loci of Leptosphaeria maculans in French field populations
|
|
View this table: [in a new window] |
TABLE 2. Frequencies of the most common minisatellite alleles, genic differentiation, and gene diversityd for minisatellite loci across four field populations of Leptosphaeria maculans
|
|
View this table: [in a new window] |
TABLE 3. Number of alleles at the six minisatellite loci in field populations of Leptosphaeria maculans
|
|
View this table: [in a new window] |
TABLE 4. Genotype diversity at three spatial scales in field populations of Leptosphaeria maculans in France
|
The Grignon and Oucques populations contained the most alleles over the six loci (Table 3). These two populations also had alleles at the MinLm2 and MinLm3 loci that occurred at low frequency and that were not present in the other field populations (Table 3). Two or three alleles at MinLm2, MinLm3, MinLm5, and MinLm6 were more frequent and were present at similar frequencies in all of the analyzed populations (Table 2). MinLm1 is physically close to the AvrLm1 avirulence gene, and the allele frequencies at these two loci were significantly different across the four populations (Tables 1 and 2). MinLm1 was excluded from analysis of gene diversity and distribution of genetic diversity because this locus was in linkage disequilibrium with AvrLm1 (4, 18). Nei's measures of gene diversity averaged over the four loci MinLm2, MinLm3, MinLm5, and MinLm6 were similar among the four field populations and ranged from 0.68 in the Grignon population to 0.75 in the Oucques population (Table 2).
Mating type distribution and linkage disequilibria in field populations.
PCR products corresponding to one of the two mating type idiomorphs could be amplified from all 401 isolates. Both mating types were present in the populations at all of the tested spatial scales, including different lesions on the same leaf, different plots within a field, and different fields within the country. No significant departures from the 1:1 ratio between the two mating types were observed at the field scale (Table 5). Linkage disequilibria were assessed with MinLm1, MinLm2, MinLm3, MinLm5, and MinLm6 and using two different approaches. Within field populations, only two pairwise tests of linkage disequilibria were significant in the Le Rheu and Versailles populations, but none of them were significant after Bonferroni's corrections for multiple comparisons were made. When the complete field-scale data set was analyzed as a single group, none of the pairwise comparisons were significant. The multilocus linkage disequilibrium values (
d) obtained for the four field populations were close to zero and not significant (P > 0.05). No significant deviation from random association of alleles was detected within the total population (
d < 0.01, P = 0.40), suggesting that recombination regularly generates new genotypes at the field level. Linkage disequilibria between the minisatellite loci and the polymorphic avirulence loci (AvrLm1 and AvrLm4) also were assessed in field populations. Overall, only one pairwise test of linkage disequilibria, which corresponded to the MinLm1 and AvrLm1 loci, was significant after Bonferroni's corrections. These two loci were in significant linkage disequilibria only within the Le Rheu and Oucques populations.
|
View this table: [in a new window] |
TABLE 5. Distribution of the mating type alleles of Leptosphaeria maculans in field populations in France
|
ST = 0.0007, P = 0.14; data not shown). The results of the AMOVA, which partitions the genetic variation hierarchically from field plot level to an among-field level, revealed that most of the allelic variability (>99.5%) was found in the 2-m2 field plots and that the genetic variation among field populations accounted for only 0.11% of the total variability (Table 7). Nei's measure of genetic identity was high for each pairwise comparison of field populations (Table 6). The highest estimate of genetic identity was obtained from a pairwise comparison of the Versailles and Grignon populations, which also were the geographically closest locations (Table 6). |
View this table: [in a new window] |
TABLE 6. Pairwise comparisons of population differentiation (GST) and Nei's unbiased measure of gene identity (I) in Leptosphaeria maculans field populations from France
|
|
View this table: [in a new window] |
TABLE 7. Analysis of molecular variance of the multilocus haplotypes for four Leptosphaeria maculans populations isolated from OSR fields in France
|
|
|
|---|
The 1:1 distribution of mating type alleles we observed also has been observed in field populations of L. maculans in Australia (10) and is expected in randomly mating populations (26). Milgroom (26) proposed several tests for random mating in natural populations of fungal pathogens. Among these, the frequent occurrence of fungal sexual structures, a high level of genotypic diversity, and gametic linkage equilibrium, all provide evidence that is consistent with sexual reproduction. In L. maculans, the frequent occurrence of pseudothecia has been documented in many instances (19, 39). The present study further supports the importance of sexual reproduction by showing that minisatellite loci were in linkage equilibrium and that genotypic diversity was very high even at small spatial scales in field populations of L. maculans in France.
Linkage disequilibrium was detected between the MinLm1 and AvrLm1 loci, which are physically close and genetically linked. The genetic distance between these loci ranges from 0 to >30 cM depending on the L. maculans genetic map considered (L. Gout, unpublished data). The suppression of recombination observed in some genetic maps may be related to nonhomologous sequences present in this region in the avirulent and virulent parental isolates used (L. Gout, unpublished data). The nonhomologous genome organization of the AvrLm1 region in avirulent and virulent isolates has been conserved in field populations from Western Europe based on the lack of apparent recombination between markers spanning this region (5). In the Le Rheu and Oucques populations, the presence of avirulent and virulent isolates could have reduced recombination near AvrLm1 and have led to the linkage disequilibrium between MinLm1 and AvrLm1 via "hitchhiking" to AvrLm1, which would be selected for when cultivars carrying Rlm1 were planted. MinLm1 and AvrLm1 were in linkage equilibrium in the Grignon and Versailles populations, but >95% of the isolates in these field populations were virulent. The lack of disequilibrium in this genomic region in these populations could result because the recombination rate in the region is not the same when populations are mainly virulent as it is when both avirulent and virulent isolates coexist in the population.
L. maculans has a mixed reproductive system (39). Ascospores are thought to be the major source of inoculum for the annual epidemics. Our results are consistent with this hypothesis and similar to those of studies of field populations from Canada (24) or Australia (10). Pycnidia are produced throughout the growing season, but the role of the asexual pycnidiospores in the epidemiology of the disease remains unclear, and secondary infections via pycnidiospores are thought to be rare in Western Europe (39). Our collections were mainly made at a time corresponding to the initiation of epidemics, and the variation we observed is consistent with the hypothesis that ascospores are the primary inoculum. Additional samples collected throughout the growing season are needed to determine the role of asexual spores in the epidemiology of L. maculans in France.
Most of the total genetic variability in these populations of L. maculans (>99.5%) was distributed at a small spatial scale, i.e., within 2-m2 field plots. These populations also had a very high level of gene (0.68
H
0.75) and genotypic (88%
G
100%) diversity even over small spatial scales. The four field populations, despite being collected in two different years and from sites more than 500 km apart, may be coevolving parts of a large panmictic population. Ascospores are airborne, can potentially be dispersed over distances of several kilometers, can survive up to 6 weeks (19, 39), and could enable this genetic homogenization.
In contrast to the minisatellite loci, the level of diversity at the virulence loci was low, with only two of the six characterized avirulence loci (AvrLm1 and AvrLm4) being polymorphic. Polymorphism was similar across the four field populations for AvrLm4 but not for AvrLm1. The two corresponding resistance genes, Rlm1 and Rlm4, have been widely deployed in France, with cultivars carrying Rlm4 in use since the 1970s (36) and cultivars carrying Rlm1 in use since the early 1990s (36). Thus, subpopulations of L. maculans may be differentiated based on selectable markers, e.g., the AvrLm1 locus, as a result of spatial heterogeneity in the use of Rlm1 cultivars (35). This genic differentiation of subpopulations may be transitory, and the genetic exchange occurring in French populations of L. maculans should homogenize the frequencies of these virulence alleles. Increased use of cultivars carrying Rlm gene(s) across all OSR growing regions could give the same result, but this scenario is less likely since the commercial success of cultivars with a given Rlm gene should be inversely proportional to the frequency of the corresponding virulence allele in the pathogen populations. Indeed, the frequencies of the avrLm4 virulence allele are very similar across all of the field populations we studied, and all of the isolates characterized in the present study carried the avrLm2 virulence allele. The corresponding resistance gene, Rlm2, was used in France prior to Rlm4 and is still present in many cultivars (36).
A low level of differentiation between field populations of L. maculans also was observed in Australia, where the high level of genetic variability revealed by amplified fragment length polymorphism analyses was distributed mainly within field populations (10). In contrast, AMOVA done on RAPD [random(ly) amplified polymorphic DNA] data attributed only 55% of the total variability to differences within field populations in Canada. Simultaneous analysis of field populations from these continents together with some European field populations is now needed (i) to evaluate the genetic structure of worldwide populations, (ii) to identify potential subpopulations, and (iii) to estimate the level of intercontinental genetic exchange. The amount of intercontinental genetic exchange is of particular interest since gene flow, due to natural long distance dispersal of propagules or linked to human activities, may lead to genetic exchange between populations located on different continents. Thus, the durability of a resistance gene on one continent could be altered by the durability of the same resistance gene on another. Resistance gene durability is important because only a few Rlm genes are known (13) and maintaining the effectiveness of each resistance gene is critical for the sustainable control of stem canker. The minisatellite markers we used in the present study were useful in determining the amount of genetic variability within and between field populations in France and may be useful for a global population study.
Our study is the first to evaluate the genetic variability of L. maculans field populations in France. Based on the level of gametic linkage equilibria observed and the frequency-dependent selection of mating type alleles, we showed that sexual recombination is an important part of this fungus's life cycle and occurs regularly in L. maculans field populations in France. Indeed, we showed that ascospores could be the primary source of inoculum and could be responsible for the leaf-spotting phase of the disease in both autumn and winter. High levels of gene flow also were inferred from the low levels of population differentiation. McDonald and Linde (25) recently proposed a framework to predict the evolutionary potential of pathogen populations based on analysis of their genetic structure. According to their model, pathogens such as L. maculans with high gene flow and mixed reproduction system have moderate to high evolutionary potential. However, the recent breakdown of the Rlm1 resistance occurred in a time frame (only 3 years) that is similar to that predicted for pathogens with high evolutionary potential, which pose the greatest "risk" of breaking down novel specific resistance genes. Fungal pathogens such as rusts or mildews assigned to this highest risk category (25) often have mixed reproductive systems and high levels of migration (genotype flow) through long-distance dispersal. Thus, the rapid breakdown of Rlm1 now questions the level of migration among L. maculans field populations.
We thank L. Bousset (INRA Bio3P, Le Rheu, France), E. Fournier (INRA, Versailles, France), and J. Enjalbert (INRA, Grignon, France) for helpful discussions during the preparation of the manuscript; anonymous reviewers for comments and suggestions; K. Louvard (INRA PMDV, Versailles, France) for isolate collections and DNA extractions; M. Chabirand and L. Coudard (INRA PMDV, Versailles, France) for plant management; J. P. Narcy and J. Roux (INRA-PMDV) for technical assistance; and C. Cogoluènhes and F. Le Coz, from Institut National Agronomique Paris-Grignon, for their contributions to this study.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»