POPULATION GENETIC STRUCTURE ANALYSIS OF sCleRotinia sCleRotioRuM ( LIB . ) DE BARY FROM DIFFERENT HOST PLANT SPECIES IN NORTHERN IRAN

The genetic structure of 65 Sclerotinia sclerotiorum isolates representing 52 field populations from three provinces in northern Iran were analyzed with Mycalial Compatibility Groupings (MCGs) and five polymorphic microsatellite loci. In total, 44 haplotypes were detected with 25 allele polymorphisms. A high level of genetic diversity was observed in about 67.6% (with clonal fraction = 0.077 to 0.20) of the regional populations studies, and the Shannon diversity index (Ho) for the whole region was found to be 0.86 (Htot). Partition of total diversity (Htot) showed that 64% corresponded to a variation in diversity within the S. sclerotiorum populations. By mycelial compatibility grouping (MCG) tests, the isolates were classified into 39 groups of which 26 MCGs were individual. Molecular and phenotypic analyses results of all of the isolates (except MCG4 and MCG23) were similar; however, the isolates in the MCG4 and MCG23 groups, with variable microsatellite haplotypes, were morphologically dissimilar. The results shown here were possibly due to high rates of outcrossing, as well as to the evolutionary potential within population of the pathogen in different locations in Iran. Nei’s genetic identity showed that populations from Golestan province and wild plants were the most diverse.


INTRODUCTION
Sclerotinia sclerotiorum (Lib.) de Bary is an important ascomycetous plant pathogenic fungus with wide geographic distribution and a diverse host range including many agronomic crops (Hartman, et al., 1999).Yield losses due to stem rot disease are variable and sometimes reach maximum level in susceptible plants (Purdy, 1979).This pathogen is the causal agent of sclerotinia stem rot in canola leading to serious losses in yield due to the lodging and premature shattering of seedpods (Gugel and Morrall, 1986).
Intraspecific variation in virulence (Marciano et al., 1983;Morral et al., 1972) and in such morphological characteristics as pigmentation of the mycelium, ascus, ascospore, sclerotial size and produc-tion (Boland and Smith, 1991;LeToumeau, 1979;Price et al., 1975;Purdy, 1955) have been reported in S. sclerotiorum.To identify the intraspecific variation of S. sclerotiorum, two presupposed unrelated criteria of mycelial compatibility groups (MCGs), and DNA fingerprinting have been employed (Kohn et al., 1990(Kohn et al., , 1991)).Therefore, knowledge of the population genetic structure and evolutionary potential of the pathogen will provide insight into the most suitable breeding strategy for durable resistance (Mc-Donald and Linde, 2002).The pathogens with higher evolutionary potential pose a greater risk of defeating resistance genes or counteracting other control methods such as applications of fungicides (McDonald and Linde, 2002).Two independent systems have been developed for the characterization of S. sclerotiorum genotypes; mycelial compatibility groups (MCGs) and DNA fingerprinting.

R E T R A C T E D
The mycelial compatibility-incompatibility grouping system is a useful method in studying the population dynamics of pathogenic and nonpathogenic isolates of fungi.This system is used for ascertaining the source of new races to a particular geographical area.The mycelial compatibility-incompatibility grouping system is an easy, quick, inexpensive, and macroscopic assay of the self/non-self recognition method and is determined using a side-by-side pairing system (Kohn et al., 1990).It has been suggested that MCGs represent genetically dissimilar individuals and each MCG is a particular genotype (Kohn et al., 1990).Infrequently, outbreeding in S. sclerotiorum, immigration of strains from other sites (Glass and Kuldau, 1992), genetic exchange, meiotic recombination (Carbone et al., 1999), mitotic recombination, transitory selection, selective neutrality, and diversifying selection (Kohli et al., 1992) are possible sources of the MCG diversity.DNA fingerprinting techniques can also be used to distinguish closely related fungal isolates.Southern hybridization of restriction-digested whole genomic DNA to a cloned probe containing a 4.5 kb repeated dispersed element of nuclear DNA from S. sclerotiorum was used by Kohn et al. (1991), and this method was used in several subsequent studies (Kohli et al., 1995;Cubeta et al. 1997).Microsatellites are widely dispersed and evenly distributed in the genome of eukaryotes and have been used to study intraspecific genetic variability within populations (Gaggiotti et al., 1999;Sirjusingh and Kohn, 2001).Microsatellite markers have high specificity, reproducibility, polymorphism and they are co-dominant markers (Vanderkoornhuyse et al., 2001).Recently, Sirjusingh and Kohn (2001) developed 23 microsatellite loci and 2 microsatellite-like polymorphic loci containing 2-10 alleles at each locus.
There are several reports concerning the population structure of S. sclerotiorum worldwide.Thirty-nine clones were identified among 66 isolates on canola (oilseed rape) in Canada from seven locations in Alberta, Saskatchewan and Manitoba (Kohli et al., 1992).Limited outcrossing among the S. sclerotiorum isolates was observed in North Carolina and California (Kohli and Kohn, 1998).Another study was focused on comparisons of the S. sclerotiorum populations collected from agricultural and wild populations in Norway (Kohn, 1995).It was found that there was genetic uniformity among populations on potato and canola; however, there was greater genetic diversity among wild populations.Sun et al. (2005) compared three S. sclerotiorum populations collected from across Europe, China and Canada.They found that genetic differentiation among and within populations was highly significant.
There is growing interest in oilseed rape and its production in Iran.Although S. sclerotiorum causes severe disease in many field and glasshouse grown crops all over the country, there is no report on analysis of the genetic variation in S. sclerotiorum in northern Iran.The aim of this research was to investigate the genetic structure within oilseed rape and the other hosts population of S. sclerotiorum in the Golestan, Mazandaran and Gillaun provinces of Iran based on molecular and morphological markers.

Isolates
Isolates of Sclerotinia sclerotiorum were collected from 52 rapeseed, lettuce, bean, tomato, and cucumber fields in Gilan, Mazandaran and Golestan, northern provinces of Iran, during the 2006-2007 growth season.Samples were collected from infected plants and the sclerotia were removed from the each plant sample.Single sclerotium was selected as an isolate.The sclerotia were surface sterilized for 1 min in 70% ethanol or 2 min in 2.5% sodium hypochlorite, rinsed in sterile distilled water, plated on potato dextrose agar medium (PDA) and then incubated at 22 ºC for two days.Each isolate was purified by transferring the single hyphal tip onto the fresh medium, and generated sclerotia were stored at -20ºC until use (Atallah et al., 2004;Cubeta et al., 1997;Willets et al., 1980).

Mycelial Compatibility Group Determination
For evaluation of inter-form MCG variability, five isolates from each field were paired together.0.5cmdiameter mycelial plugs were obtained from the edge of two-day-old colonies on PDA.Three mycelial plugs were paired in 6.5 cm-diameter Petri dishes containing PDA amended with 75 μl of Wilton's red food coloring per liter of culture medium (Kohn et al., 2006).Petri dishes were incubated in the dark at 22°C for 14 days.Pairings of the 65 isolates were performed in a pyramidal design, where groups of 10 isolates were paired in all-pairwise combinations.Three replications were made for each pairing.Only incompatible isolates were paired subsequently.Compatible isolates were distinguished by the fusion of mycelia, without an accumulation of red dye in the fusion zone.Incompatible reactions produced a barrage recognized by an obvious red line on the bottom side of Petri dishes or by the formation of aerial mycelia along the barrage line.Pairings that yielded questionable reactions were repeated to ensure accurate results.

Temperature treatments
The radial growth of the S. sclerotiorum isolates was assessed at five different temperatures, including the optimum temperature (22°C), as well as 10 and 5°C below and above (12, 17, 27 and 32°C).PDA plates were inoculated with a 5 mm-diameter plug of the colonized agar and incubated in the dark at the various temperatures.Colony diameter was measured daily and the values averaged until the colony neared the edge of the dish.In a completely randomized design, four Petri dish replicates of each isolate were used four times.Data were analyzed by ANOVA using MSTAT-C program and mean values were compared at the p≤0.05 level using Duncan's multiple range test.The color of each colony was noted 21 days after the initial growth on PDA.

DNA extraction
DNA was extracted using the rapid mini-preparation method (Liu Don et al., 2000) with some modifications.This procedure included the following steps: (i) A small lump of the mycelia was added to a 1.5 ml Eppendorf tube containing 500 μl of lysis buffer (400 mM Tris-HCl [pH 8.0], 60 mM EDTA [pH 8.0], 150 mM NaCl, 1% sodium dodecyl sulfate) using a sterile toothpick, then the tube was left at room temperature for 10 min.(ii) After adding 150 μl of potassium acetate (pH 4.8; which was made of 60 ml of 5 M potassium acetate, 11.5 ml glacial acetic acid and 28.5 ml distilled water), the tube was vortexed briefly and spun at 10,000g for 1 min.(iii) The supernatant was transferred to another 1.5 ml Eppendorf tube and centrifuged again as described above.After transferring the supernatant to a new 1.5-ml Eppendorf tube, an equal volume of isopropyl alcohol was added.The tube was mixed by inversion briefly.(iv) The tube was spun at 10,000 g for 2 min, the supernatant was discarded, and the pellet was air-dried.Finally, the pellet of purified DNA was diluted in 1X TE (10 mM Tris-HCl, pH 8.0, 0.1 mM EDTA) to a working concentration of 10-20 ng/μl and stored at 5°C (Hopwood et al., 1997).

Microsatellite primers, PCR conditions, separation of PCR products and data analysis
Five sets of microsatellite primers were used in this study, including (AGAT)14, (AAGC) 4 , (CATA) 25 , (CA) 9 , (GT) 10 and (TACA) 10 (Sirjusingh and Kohn, 2001).The PCR reaction mixture included 10-20 ng of the purified DNA and the reaction buffer (100 μM each of dATP, dCTP, dGTP and dTTP, 200 nM of microsatellite primer and 0.8 units of Taq polymerase, 10 mM Tris-HCl pH 8.3, 1.5 mM MgCl 2 , 50 mM KCl, 100 μg mL -1 gelatine, 0.05% tween 20 and 0.05% Nonidet P-40).The final reaction volume was adjusted to 50 μL with deionized H 2 O.All of the reagents were obtained from Fermentas Inc., USA.Amplification was carried out using initial denaturation at 95°C for 8 min, followed by 35 cycles of primer annealing at 59°C (for all microsatellite primers) and extension at 72°C for 60 s, with a 5 min extension at 72°C on the final cycle.The PCR products were separated on a denaturating agarose gel (2.6% w/v).Gels were stained with ethidium bromide, visualized under UV light and digitally documented with the gel documentation UVP-V system.The gel was run at 90 W for 90 min (Nicholson et al., 1997).All polymorphic alleles were identified from each microsatellite primer combination and bands representing alleles were scored as present (1) or absent (0).Nei's genetic distance matrix (Nei and Li, 1979) was prepared and bootstrap analysis with 2000 replications was performed to generate a dendrogram of unweighted pair-group mean analysis (UPGMA; Sokal and Michener, 1958) using the Treecon-1.3bprogram (Van de Peer and de Wachter, 1994).

Polymorphism Information Content (PIC)
The polymorphism information content determines the frequency of allele polymorphism in the gene locus to study population.PIC was calculated for each locus as follows: where pi was the relative frequency of ith alleles.

Determination of genotypic diversity
Isolates with the same allele combination for each of the five loci were located in a haplotype group.
Isolates present in a haplotype were obtained from clonal reproduction (asexual reproduction or sexual reproduction of the type homothalism).Then the clonal fraction was calculated for every three populations.The clonal fraction represents the frequency and prevalence of colonies in populations, calculated from N-G/ N, where N is the population size and G is the haplotype size in the population (Zhan et al., 2002), and the genotypic diversity is calculated from G/ N.

Shannon-Wiener's Index
Shannon-Wiener's index, H 0 , was calculated for each population as follows: where pi was the frequency of ith haplotype (Shannon and Wiener, 1949)

Morphological variability and mycelial compatibility groups among S. sclerotiorum isolates
The growth rate of the isolates differed significantly.
The isolates grew faster at 22°C, then at 27, 17 and 12°C, respectively.Isolates grew normally at 12, 17, 22 and 27°C, but weakly at 32°C.The colonies filled the 8-cm diameter Petri dishes in 3 days at 22°C.ANOVA for the growth rate data showed highly significant P-values among the isolates at all temperature treatments (Table 1).The color of the colonies on PDA growth medium varied among the isolates.There were three main colony colors: brown to dark brown pigmentation restricted to the region around the inoculation disc; beige pigmentation across the entire colony, and white (no pigmentation) (Fig. 1).Five isolates (R28, R30, R33, R36 and R51) exhibited brown pigmentation around the fungal discs and four isolates (R34, R38, R58 and R64) produced white colonies without any pigment.The colony color of the rest of the isolates was beige (Table 1).

Mycelial compatible group variability among S. sclerotiorum isolates
Assessments of compatibility were based on mycelia continuity between the interacting colonies without formation of either a strip of thin mycelium or aerial mycelium, and the uniform distribution of sclerotia in the plate.Evaluation of compatibility was based on the failure of the two colonies to fuse, which was reflected by the formation of a strip of this mycelium or aerial mycelial at the interaction zone (Kohli et al., 1992).Mycelia incompatibility can also be indicated by the formation of a dark line along the interaction zone associated with the red food dye (Kohli et al., 1992).Thirty-nine MCGs were determined among the 65 studied isolates: 26 isolates were established as independent MCGs.The isolate that belonged to the independent MCGs was compatible only with itself.MCG4 and MCG18 consisted of eight and three isolates, respectively.MCG1, MCG5 and MCG 23 included four isolates and eight MCG included only two isolates.The

R E T R A C T E D
remaining 26 isolates were compatible only with themselves (Table 1).

Genetic diversity among the S. sclerotiorum isolates
The microsatellite primers exhibited 25 clear polymorphic alleles from the 65 fungal samples.The number of polymorphic alleles per locus ranged from 2 to 8 (Table 2).
In order to determine the genetic relationships among populations of fungal isolates, a separate matrix was used for the data obtained from the 25 polymorphic alleles, and the fungal isolates were subsequently grouped by UPGMA cluster analysis.The dendrogram representing microsatellite marker polymorphisms among the fungal isolates is illustrated in Fig. 3.
The Shannon diversity index of haplotypes for the three regions (H tot ), ranged from 0.602 for Gillan province to 1.568 for Mazandaran province, with a mean of 1.083(H avr ).Overall, diversity was higher in the population from Mazandaran than in those from Gillan and Golestan.Partition of total diversity showed that 64% corresponded to a variation of diversity within S. sclerotiorum populations, while only 36% of diversity was responsible for variability among these populations (Table 3).
According to the dendrogram, nine major clusters can be defined among the 65 isolates.Although a high level of diversity was observed between the

R E T R A C T E D
clusters, some isolates within clusters were identical for all the microsatellite markers.For instance, isolates R28 and R33 had similar alleles for all the microsatellite markers.Overall, 44 distinct isolates were identified among the 65 isolates based upon microsatellite alleles.

DISCUSSION
The genetic diversity of 65 isolates of S. sclerotiorum on different hosts representing three geographic populations from Iran were studied using five microsatellite loci.In the present study, we found that microsatellite markers were very efficient in identifying genetic variation among the isolates, although we employed only five of the 25 marker sets described previously (Sirjusingh and Kohn, 2001) that had revealed polymorphisms among the Iranian isolates.When we compared our findings with those by Atallah et al. (2004), we found one more allele at (TACA) 10 , two more alleles at (GT) 10 , three more at (CATA) 25 and four more at (AGAT) 14 (AAGC) 4 loci.In contrast, Sexton and Howlett (2004) identified more alleles at five microsatellite loci than reported by Atallah et al. (2004).
We found 44 different clones (haplotypes) among the 65 isolates representing the population, indicating a high rate of variability in the region (Ĝ/N%=66.7%).On average this is similar to that found in Australian populations (Ĝ/N%=36 to 80% (Sexton and Howlett, 2004); Ĝ/N%=28 to 68% (Sexton et al., 2006)) by analyzing eight microsatellite loci, and Atallah et al. (2004) showed that 92% of the variability among 167 isolates was found within subpopulations in the Columbia basin of Washington State.The genotypic diversity obtained in Turkey (63%) (Mert-Türk et al., 2007) was also similar to the current study.
Although a high level of variation between clusters was observed, some of the isolates within a cluster were similar at all marker microsatellites.In 50% of the isolates, S. sclerotiorum divided into 9 large clusters, which indicates the diversity in these isolates is high.Based on this cluster, the population of the Golestan province completely separated from the other populations (Gilan and Mazandaran).In addition, the Sclerotia wild population isolated from wild mustard were placed in separate haplotype and mycelia compatibility groups.According to the semi-arid conditions of Golestan province compared to the temperate climates (Mazandaran and Gilan) this phenomenon was predictable and in accordance with the evolutionary history of population haplotypes (Carbone and Kohn, 2001), which indicated that isolates may

R E T R A C T E D
be grouped into subtropical, temperate, wild and two relatively recently evolved temperate subtropical populations.This geographic clustering infers that ecological conditions may contribute to adaptations associated with growth temperature range, light intensity, or sclerotial vernalization requirements.
Although we did not test for outcrossing, environmental conditions may facilitate outcrossing.The isolates were collected from fields in a geographical region of Iran where diverse crops are also cultivated.The environmental conditions in the region may favor sexual recombination within S. sclerotiorum.Oilseed rape is a new crop to this area and the high level of polymorphism may reflect the movement of S. sclerotiorum onto this crop from several wild plant hosts.
Clonal fractions of populations from Iran (ranging from 0.077 to 0.20) and from California (0.12) (Malvárez et al. 2007) are similar and confirm that most populations have a small clonal fraction.
Shannon's index of diversity was highest for the isolates from Mazandaran, which had more polymorphisms than those from the other regions.Therefore, the potential for the emergence of isolates resistant to fungicides and pathogenic on different canola cultivars is possible in Mazandaran.
Mycelial compatibility grouping, a phenotypic marker system controlled by multiple loci, was often associated with groups of identical or closely related microsatellite haplotypes (these results support the hypothesis that isolates within a single MCG and sharing the same microsatellite alleles may be clonal) except for MCG4 and MCG23 whose microsatellites showed polymorphisms between members of these two groups.For instance, MCG23 included 4 isolates (R50, R57, R58, and R59) that were all compatible when paired on PDA media.R50, R57 and

R E T R A C T E D
R59 shared identical microsatellite alleles; R58 differed from the three other isolates by four microsatellites and clustered closely with different isolates in the UPGMA dendrogram (Fig. 3).We confirmed the pairing of the isolates belonging to this MCG in two further replications, as the molecular data did not agree with the pairing. .
It may be indicated that isolates within a MCG tend to separate and form independent groups.The existence of unique MCGs in a sampling area suggests that new MCGs and thereby new genotypes are evolving.The emergence of new genotypes could indicate that MCGs could be adapted to specific field microclimates.In addition, the evolving of new genotypes may be associated with increased cropping of canola in such areas (Hambleton et al., 2002), leading to movement of S. sclerotiorum onto this crop from several other host plants, or due to an ascogenous system in sexual reproduction, indicating a sexual population in contrast to the clonal population structure as indicated by other studies (Kohli et al., 1992;Cubeta et al., 1997).
However, given the non-clonal nature of MCG4 and MCG23, it will be important to examine a large number of additional isolates to confirm this hypothesis.This could explain why members of one MCG change in one or more microsatellite genetic groups over time.The colony diameters of the isolates belonging to the same MCG were statistically significant in agreement with each other at all temperature points.In addition, the growth of all isolates belonging to MCG23 (R50, R57, and R59) was similar at all temperature points, except for R58 (Table 1).More variation in the mean growth rate was observed in isolates of MCG4.We presume that the variation in coloration is correlated with genetic differences within the population.The colony colors of the isolates within the same MCGs were similar, except for MCG4 and MCG23.MCG23 consisted of four isolates; the color of isolate R58 was dominantly white, while the rest were beige (Table 1).Carpogenic germination of S. sclerotiorum sclerotia is affected greatly by soil moisture and tem-perature.Hao et al. (2003) and Clarkson et al. (2004) reported that apothecium production by S. sclerotiorum is largely confined to a 10-15°C soil temperature and a -0.03 or -0.07 MPa soil moisture.Stem rot symptoms caused by S. sclerotiorum in northern Iran can be observed in mid-May but predominantly occur in early June.The average air temperature of the last 3 years was 14.2°C in April and 20.4°C in May.The rainfall for the same time period was 46.5 mm in April and 21.6 mm in May and relative humidity 82 and 79.3 percentage (data from the State Meteorological Station based in Gharakhail).It has been reported that 15-25°C air temperature and high humidity (optimum >95% RH; but also possible as low as 50-60%) are required for ascospore germination (Clarkson et al., 2004).Although we did not test for outcrossing, environmental conditions may facilitate outcrossing.

CONCLUSION
The isolates were collected from fields in a geographical region of Iran where diverse crops are cultivated.Given the high genetic diversity of this fungus, we suggest that the environmental conditions in the region may favor sexual recombination within S. sclerotiorum.Oilseed rape is a new crop to this area and the high level of polymorphism may reflect the movement of S. sclerotiorum onto this crop from several wild plant hosts.This is an original report on genetic variation within a population of S. sclerotiorum in Iran presenting molecular data and morphological characters demonstrating that the population is genetically diverse.

Fig. 2 .
Fig.2.Microsatellite analysis of Sclerotinia sclerotiorum isolates amplified with the (CATA) primer.PCR amplicons were separated on a 2.6% agarose gel in 0.5 X TBE.Bands were stained with ethidium bromide and visualized on a UV-Transeluminator.

Fig. 3 .
Fig. 3. Unweighted pair-group mean analysis dendrogram of genetic distance among the 65 fungal isolates based on Nei's coefficients.The numbers given above the lines indicate the bootstrap values of 2000 replicates.

Table 1 .
MCGs, colony color and growth rate of S. sclerotiorum isolates at different temperatures after 2 days of growth.
Notes:MCGs, mycelial compatibility groups.Colors were recorded 21 days after sub culturing on PDA, The numbers represent the colony diameters (inoculation disc subtracted) of the isolates on PDA medium 2 days post-transfer on 8 cm plates.Mean values within a column followed by the same letter are not significantly different at the P ≤ 0.05 level.

Table 2 .
Number of alleles recognized for each microsatellite marker set in the S. sclerotiorum population used in the present study compared with those from previous studies

Table 3 .
Shannon diversity index and Clonal fraction of haplotypes S. sclerotiorum populations from different locations.