Genetic structure analysis of Fusarium oxysporum f. sp. phaseoli populations on common bean by using SSR markers

Document Type : Original Article

Authors

Department of Plant Protection, Faculty of Agriculture, University of Ilam, Ilam, Iran

Abstract

Determination of the genetic structure of the Fusarium oxysporum populations provides different levels of pathogen behavior and environmental compatibility in the management of root rot disease in bean farms. In this study, Short Sequence Repeat (SSR) markers were used to determine the genetic structure and estimate genetic diversity in sixty F. oxysporum isolates from five counties in Ilam province located in the west of Iran (Ilam, Ayvan, Malekshahi, Sirvan, Chardavol). A set of five microsatellite primer pairs revealed one allele in each locus across the populations. The average number of alleles (Na) and the mean number of effective alleles (Ne) observed among populations were 1.484, and 1.438 respectively. The number of Genetic diversity (H) and Shannon's coefficient (I) were also maximum in Chardavol (He = 0.442, I=0.333) but the minimum values were estimated in Ilam (He = 0.228, I = 0.333). The minimum genetic distance was found between Ayvan and Malekshahi (0.046) but the maximum genetic distance (0.151) was revealed between Ilam and Chardavol. The total gene diversity (Ht) and genetic variability between the subpopulations (Hs) were estimated to be 0.312 and 0.261, respectively. Gene diversity attributable to differentiation among the population (Gst) was 0.162, while gene flow (Nm) was 0.5717. Cluster analysis based on UPGMA showed the lowest genetic distance between Ayvan and Malekshahi. The dendrogram showed a clear break between the population from Sirvan and Chardavol and other remaining populations. Results from this study could be useful in breeding programs for developing root rot resistant cultivars.

Keywords


INTRODUCTION

Bean (Phaseolus vulgaris L.) is an important legume in the world which has high nutritional and economic value. Plant pathogens can cause important bean diseases among which wilt disease caused by Fusarium oxysporum f. sp. phaseoli (FOP) is one of the most important diseases of common bean (Karimian et al. 2010). Symptoms of Fusarium wilt initially appear as yellowing and wilting and in case of infection at younger ages, plants may also be stunted. F. oxysporum isolates exhibit great variability in morphology and aggressiveness (Abbas 1995, Belabid et al. 2004). Pathogenic isolates of F. oxysporum often show host specificity and based on the host may be subdivided into formae speciales (Armstrong et al. 1981). Control of F. oxysporum infection in the field is difficult because the pathogen can survive for a long time in the form of mycelium in infected plant debris or in the form of chlamydospores in soil (Haware et al. 1996). Fusarium wilt disease management can be conducted through the use of resistance cultivars (Jalali & Chand 1992).

The genetic diversity of plant pathogens plays a significant role in the understanding of plant pathogens evolution and their correlation with plant hosts, which helps us in applying any preventive measure against them (Karimian et al. 2010). Knowledge of genetic diversity is necessary for resistance development and shifts identification in races or population structure that might occur (McDonald 1997). Identification of diversity by morphological characteristics is highly variable in Fusarium isolates and these characteristics are influenced by cultural conditions. Many DNA based methods have been used to study diversity in pathogenic Fusarium populations (Kiprop et al. 2002, Sivaramakrishnan et al. 2002a, Belabid et al. 2004, Nourollahi and Aliaran 2017). Molecular markers have been used to characterize the genetic diversity of different isolates of Fusarium oxysporum (Edel et al. 1997, Kistler 2001, Alves-Santos et al. 1999) formae-speciales (Gherbawy 1999, Alves-Santos et al. 2002). RAPD is an effective method for discriminating of FOP isolates (Namvar Hamzanlouyi et al. 2006, Alves-Santos et al. 2002, Woo et al. 1996).

Short sequence repeats (SSRs) are powerful tools for taxonomic and population genetic studies (Britz et al. 2002). Alleles vary according to the number of present repeat units, although it has been shown that mutations could be responsible for allele length variations in SSR analysis (Burgess et al. 2001, Slippers et al. 2004). Because of the high resolution provided by SSRs (Enjalbert et al. 2005), they have also been used in different fungal species including Sclerotinia sclerotiorum (Sirjusingh & Kohn, 2001), Rhizoctonia solani (Mwang Òmbe et al. 2007) and Ascochyta rabiei (Nourollahi et al. 2011). Some researchers have already worked on a molecular variation of different F.oxysporum form specials (Stewart et al. 2006). Bogale et al. (2005) determined that polymorphism shown with eight SSR markers should be enough for study of genetic variability in F. oxysporum complex. Based on the studies, microsatellite markers could differentiate the four races of F. oxysporum f. sp. ciceris causing varied levels of wilting with differential host cultivars (Barve et al. 2001). The ISSR marker has been used for determining genetic variations between several populations of F. oxysporum f. sp. ciceris (Bayraktar & Dolar 2008).

Most studies of F. oxysporum have focused on plant pathogenic isolates (Mohammadi et al. 2004). It is difficult to control pathogenic fungi with a high level of genetic variation, as they adapt more quickly to any control strategy such as applying resistant cultivar. Therefore, knowledge of the genetic diversity of FOP has contributed to the development of disease control strategies (Kistler 2001). This study was conducted to analyze the diversity and genetic relationships of F. oxysporum isolates collected from common bean fields in Ilam province using SSR markers.

 

MATERIALS AND METHODS

 

Sampling and fungal isolation

 

Roots of bean plants with wilt symptoms and black lesions were randomly collected from five different regions including Ilam, Ayvan, Malekshahi, Sirvan, and Chardavol from 2015 to 2016, in the west of Iran. Isolates of each region were considered as a population (Table 1, Fig 1). Infected samples were cut into 2- to 5-mm-long pieces, surface sterilized by dipping into 5 % sodium hypochlorite, (5% NaOCl) for 2–3 min, washed three times with sterile distilled water, dried with sterile filter paper, and finally, were placed on potato dextrose agar (PDA) plates. All the samples were incubated for three days at 25 °C with a 12-h photoperiod to induce the conidia production. The fungus was isolated and purified using the hyphal tip and single spore methods (Hawker 1950). After the identification, characterized isolates were stored for a short time on SNA culture (Leslie et al. 2006) at 4 ºC, and for perennial time, spores of each isolate stored at 4 °C in micro tubs containing sand. F. oxyspourum isolates were initially identified according to the morphological and microscopic characteristics of the macroconidia, phialides, microconidia, chlamydospores, and colony growth traits described by Jens et al. 1991, Nelson et al. 1983, Barnett & Hunter 1972, Leslie et al. 2006, Maina et al. 2017.

 

Table 1. Geographical origin of Fusarium oxysporum f. sp. phaseoli populations from Ilam province.

Sampling region

Origin site

No.Isolates

Population number

Ilam

Ilam

14

1

Ayvan

Khoran

10

2

Malekshahi

Mehr

12

3

Sirvan

Lomar

9

4

Chardavol

Sarableh

15

5

 

Fig. 1 Geographical location of sampling sites of Fusarium oxysporum f. sp. phaseoli populations in Ilam province.

 

Pathogenicity test

Pathogenicity test was done in the greenhouse at Ilam University and one representative isolate from each population was randomly chosen for a pathogenicity test on local cultivar to evaluate the infective ability of isolates as described by Perveen et al. (2007). The standard root dip inoculation technique was used to analyze the pathogenicity of isolates on common bean plants (Alves-Santos et al. 1999). Briefly, roots of one-week-old seedlings were washed in running tap water; about 2 cm pieces was cut from each tip, and subsequently, the roots were dipped for three min in spore suspension with 1×106 spores/ml concentration. Inoculated seedlings were transplanted into plastic pots with sterilized vermiculite. Finally, the plants were incubated for 24 h in a controlled climatic chamber at 23-24 °C with 60-75 % relative humidity. Fusarium yellows severity was recorded at one-week intervals after inoculation by using the Centro International de Agricultura Tropical scale (Alves-Santos et al. 1999) score as 1 (no external symptoms) to 9 (dead or severely infected plants with 100 % of the foliage showing wilting, chlorosis, necrosis, and/or premature defoliation).

 

DNA extraction

DNA was extracted from 60 FOP isolates. To obtain the mycelia mass, liquid cultures were initiated by adding 2–4 mm2 pieces of filter papers containing the fungus to 250 mL Erlenmeyer flasks containing 100 mL PDB medium (potato dextrose broth). All the flasks were incubated at room temperature, approximately 25 °C on a rotary shaker for 6–8 days. Mycelium was collected by filtration through sterile filter paper with a vacuum funnel. Mycelia were harvested, and stored at −20 °C. DNA extraction was done according to the CTAB (Cetyltrimethyl-Ammonium Bromide) method (Doyle & Doyle 1987). Mycelia were ground in liquid nitrogen and suspended in a 2 % CTAB extraction buffer (1.4 M NaCl, 0.1 M Tris-HCl, pH 8.0, 20 mM EDTA, 0.2% β-mercaptoethanol). The samples were treated with five units RNAse at 37 °C for 30 min., and then extracted with chloroform-isoamyl alcohol 24:1 (v/v). DNA in the supernatant was precipitated with isopropanol, rinsed with ethanol, and adjusted to a final concentration of 20 ngμl-1 in TE (pH 7.4). The quality of the extracted DNA was visually checked on 0.8 % agarose gel.

A set of five locus-specific primer pairs for SSRs (Table 2) described by Bogale et al. (2005) were selected. For each primer pair, primer aliquots for each marker were prepared by mixing equimolar amounts of appropriate forward and reverse primers in 1X TE buffer (1 mM EDTA, 10 mM Tris–HCl, pH 8.0) and used for the amplification of individual microsatellite loci. PCR amplification was perfor-
med in a 25 μL volume reaction containing 2.5 μL of 10X PCR Buffer, 0.2 mM of dNTPs mix (100 mM of each dNTPs), 1.5 mM MgCl2, 1 μL of each forward and reverse primer, 0.5 U of Taq polymerase with 25 ng of template DNA. Amplification was performed using Biometra thermal cycler (USA). PCR conditions for SSR were as Follows: the PCR programmed had one initial denaturation step at
94 °C for 4 min Followed by 35 cycles of 94 °C for 30s, annealing temperature (59 -62°C) for 30s (appropriate annealing temperature were used for each primer set (Table 2)), and 72 °C for 5 min. The thermal cycles were terminated by a final extension of 10 min at 72 °C. The amplified products were resolved in 2.0 % agarose gel at 60 V cm-1, using Tris Boric Acid EDTA (1X TBE) buffer and strained with DNA Safe Stain and photographed under UV Trans Laminator with Gel documentation, Intas.

 

Molecular data analysis

Data analyses of populations were defined according to the geographic locations. The bands generated by SSR primers that were repeatable and visible with high intensity were scored manually for the presence (1) or absence (0) of bands in each isolate. The pair-wise distance among the isolates was calculated from the binary matrix using the simple mismatch coefficient (Sneath & Sokal 1973). Genetic similarity between pairs was estimated using Jaccard's similarity coefficient. Similarity coefficients were used for the construction of UPGMA dendrogram (Rolhf 1990). For each primer pair, the polymorphic information content (PIC) and marker index (MI) were calculated. The polymorphic information content (PIC) was calculated using PICi =2fi (1-fi), where i is the information of marker I, fi is the frequency of the amplified allele (presence of fragments), and (1 – fi) is the frequency of the null alleles (Roldan-Ruiz et al. 2000). The genetic variation was measured in terms of genetic diversity and was computed by averaging PIC estimates overall loci (Nei and Li 1979). The marker index (MI) was calculated by MI = PIC x EMR, where EMR is the “effective multiplex relationship” given by the product of the total number of fragments per primer (n) and the fraction of polymorphic fragments (β) (Varshney et al. 2007). Genotypic diversity (H) among isolates was estimated from allelic frequencies using the equation H = 1 - Σ Xi2, where, Xi is the frequency of the ith allele (Nei 1973). The coefficient of population subdivision (Gst) was computed as (HtHs)/Ht, where, Ht is the total genetic diversity and Hs is the average gene diversity overall subgroups (Nei 1973). The allele frequencies at polymorphic loci, the Nm values (effective migration rate), and the genetic identity among populations for characterizing genetic variation, observed number of alleles (Na), effective number of alleles (Ne), Nei's gene diversity He) and Shannon's information (index (I) was calculated in both origin sites and subspecies levels. Mean values of gene diversity in total populations (Ht), gene diversity between populations (Hs), the proportion of gene diversity attributable to differentiation among populations (Gst), and the estimate of gene flow from Gst (Nm) were obtained across loci (McDermot & MacDonald 1993). The relationships among populations were estimated using UPGMA dendrograms based on Nei’s genetic distances (1978), and also Analysis of molecular variance (AMOVA), which is the hierarchical distribution of genetic variation among and within populations, was assessed. Principal coordinate analysis (PCA) was performed to evaluate the genetic differences among isolates within populations. All above calculations were performed using POPGENE ver. 1.31 (Yu et al. 1999) and Gen Alex ver. 6.5 (Peakall & Smouse 2006).

 

Table 2. Characterization of SSR primers used in this study.

PIC (the polymorphic information content), EMR (effective multiplex relationship ratio) MI (marker index)

 

PIC

MI

EMR

Annealing

Repeat motif

Primer sequence (5′–3′)

Primer locus

0.367

109.2

300

60

 

(CA)9

 

F:TGGCTGGGATACTGTGTAATTG

R:TTAGCTTCAGAGCCCTTTGG

MB9

0.285

114

400

61

 

(AAC)6

 

F:TATCGAGTCCGGCTTCCAGAAC

R:TTGCAATTACCTCCGATACCAC

MB10

0.396

158.4

400

59

(CCA)5

 

F:CGTCTCTGAACCACCTTCATC

R:TTCCTCCGTCCATCCTGAC

MB14

0.405

162

400

60

 

(CA)21

F:ACTGATTCACCGATCCTTGG

R:GCTGGCCTGACTTGTTATCG

MB17

0.364

 

0.363

145.6

 

137.8

400

 

380

62

 

(CAACA)6

F:GGTAGGAAATGACGAAGCTGAC

R:TGAGCACTCTAGCACTCCAAAC

MB18

 

Average

 

 

 

 

 

 

 

 

 

 

 

 

RESULTS

 

Morphological features and pathogenicity test

A total of 60 isolates of FOP were isolated and purified from infected bean plants (Table 1). Based on morphological characteristics, the isolates were identified as FOP. The isolates showed significant variations in cultural characteristics, production of microconidia, macroconidia, and chlamydospores. Our morphological identifications were further confirmed by the molecular method using specific SSR primers used in this study for FOP isolates. Morphological identification was based on characteristics of the macroconidia, phialides, microconidia, chlamydospores, and colony growth traits (Maina et al. 2017). The colony diameter for the isolates at seven days of growth differed among the isolates incubated at 25-26 °C on PDA. Moreover, the isolates showed variations in their micro and macroconidia in terms of size and number of septa. The macroconidia in all isolates had more than one septum. They were slightly sickle-shaped with slightly foot-shaped basal cells and attenuated apical cells. The aseptate microconidia were produced abundantly from short monophialides. The chlamydospores were present at the terminal or intercalary positions, usually in singles and in pairs. The results of the pathogenicity test indicate that all FOP isolates were virulent and caused wilting and death of the seedlings. In the early stage, symptoms appeared as yellowing of the lower leaves, then symptoms progressed up toward upper leaves, at this stage, wilting plants were observed, Vascular discoloration was typically observed in infected plants. Typical symptoms of wilt disease were observed 20–40 days after inoculation.

 

Primers information and allelic frequency

 

Based on the microsatellite data, the polymorphic information content (PIC) varied from 0.258 (Primer MB10) to 0.405 (Primers MB17), which reflects the informative content of the used primers. EMR (effective multiplex relationship ratio) varied from three to four. The marker index (MI), which incorporates the informative content of the marker (PIC), the number of fragments per primer pair and the fraction of polymorphic fragments (EMR), varied from 109.2 (MB9) to 162 (MB17) (Table 2).

To investigate the genetic structure of F. oxysporum populations, the isolates were divided into five populations including Ilam, Ayvan, Malekshahi, Sirvan, and Chardavol. From the microsatellite loci used in this study, a total of 300 DNA bands (200-500 bp) were amplified. A total of 16 different alleles were produced by SSR primers (Table 3, Fig. 2).

 

The average number of bands per gene locus was one. The highest number of bands was 58 at MB14 and the lowest number at 55 at MB9. The average number of alleles (Na) observed among populations was 1.484 and the mean number of effective alleles (Ne) among populations (1.438), with the highest number of alleles (1.480) in Malekshahi population then that in Chardavol and Ayvan populations with 1.476 and 1.432, respectively and the lowest number of effective alleles was found in Sirvan population with 1.394 (Table 4).

 

Fig. 2. Banding pattern of MB17 (A) and MB14 (B) and MB10 (C) primers on Fusarium oxysporum f. sp. phaseoli. M indicates 1Kbp DNA ladder and F1, F2, …, F24 indicate the isolates code.

 

Table 3. Number of alleles at each locus across the five populations of Fusarium oxysporum f. sp. phaseoli populations from Ilam province.

SSR primer

Number of alleles per population

Total alleles

Chardavol

Sirvan

Malekshahi

Ayvan

Ilam

MB9

3

2

3

3

3

12

MB10

1

2

2

2

4

14

MB14

2

3

4

2

2

15

MB17

3

3

3

3

3

13

MB18

3

4

3

3

4

16

Average

2.4

2.8

3

2.6

3.2

14.2

 

 

Table 4. Genetic diversity estimates for Fusarium oxysporum f. sp. phaseoli populations from Ilam province based on five microsatellite loci.

Population

Na

Ne

I

H

Ilam

1.211

1.409

0.333

0.228

Ayvan

1.579

1.432

0.395

0.259

Malekshahi

1.579

1.480

0.431

0.288

Sirvan

1.368

1.394

0.363

0.240

Chardavol

1.684

1.476

0.442

0.292

Mean

1.484

1.438

0.393

0.261

 

Genetic diversity among populations

Analysis of obtained SSR data for polymorphism studies among sixty F. oxysporum isolates was performed. The total gene diversity (Ht) and gene diversities between subpopulations (Hs) were estimated to be 0.312 and 0.261, respectively. Gene diversity attributable to differentiation among populations (Gst) was 0.162, while gene flow (Nm) was 0.571. The average genetic distance was calculated among the five populations. Nei’s pairwise genetic distances between the populations varied from 0.046 to 0.151. The lowest genetic distance was found between Ayvan and Malekshahi, while the highest genetic distance was revealed between Ilam and Chardavol (Table 5). Cluster analysis (UPGMA) result revealed the genetic relationships between the populations based on the SSR data, the dendrogram revealed a distinction between the Chardavol population and the four remaining populations (Fig.3). PCA (Principal Component Analysis) using SSR data showed the genetic differences among isolates within populations and gene flow between different populations (Fig. 4) accordingly, the first and second principal coordinates was 27.42 % and 19.52 % of the variation, respectively. There was no clear separation among individuals from different populations; however, isolates in the same populations tended to gather together. This suggests the geographical regions of sampling play an important role in the formation of populations. The AMOVA of genetic variation in FOP populations revealed that 13 % of the variance occurred among populations and 87 % within populations (Table 6).

 

Fig. 3 Dendrogram of genetic relationships between five populations of Fusarium oxysporum f. sp. phaseoli from Ilam province reconstructed by UPGMA using Nei’s genetic distance.

Fig. 4 Principal Component Analysis based on SSR data from Fusarium oxysporum f. sp. phaseoli populations from Ilam province.

 

Table 5. Nei's original measure of genetic distance between pairs of Fusarium oxysporum f. sp. phaseoli populations from Ilam province.

 

Chardavol

Sirvan

Malekshahi

Ayvan

Ilam

Population

 

 

 

 

-

Ilam

 

 

 

-

0.075

Ayvan

 

 

-

0.046

0.054

Malekshahi

 

-

0.093

0.074

0.150

Sirvan

-

0.082

0.089

0.085

0.151

Chardavol

 

 

 

 

 

 

 

Table 6.  Analysis of molecular variance (AMOVA) within and between Fusarium oxysporum f. sp. phaseoli populations from Ilam province.

Source of variation

df.

Sum of squares (SS)

Mean of squares (MS)

Percentage of variation

P value

Among populations

4

29/590

7/398

11

0.126

within populations

55

149/376

2/716

89

 

Total

59

178/967

 

100

 

 

 

 

 

DISCUSSION

In this research, the population genetic variability of F. oxysporum from five bean growing regions in the west of Iran was analyzed. Essential information was generated and showed that SSR markers are powerful discriminative tools for studying diversity in F. oxysporum f. sp. phaseoli populations. The advantage of microsatellite markers is their high specificity, high polymorphism, good reproducibility and unambiguous score ability which provides a powerful means for taxonomic and population genetic studies (Tenzer et al. 1999, Britz et al. 2002). In microsatellite analysis, alleles modify according to the number of repeated units, and mutation has also been responsible for allelic difference (Burgess et al. 2001, Slippers et al. 2004, Nourollahi & Aliaran 2017). In this study, SSR marker grouped F. oxysporum f. sp. phaseoli isolates into three major groups based on their geographical regions. In the similar studies, 114 F. oxysporum f. sp. ciceris isolates and 101 F. oxysporum f. sp. lentis isolates from Iran were placed into two major groups (Nourollahi & Jalali 2016, Nourollahi & Aliaran 2017). In another study, based on ISSR and RAPD markers, 64 F. oxysporum f. sp. ciceris isolates from India were located into two major classes (Dubey & Singh 2008). In the most of worldwide studies on plant pathogenic fungi using molecular marker have highlighted the importance of molecular methods for illustrating genetic diversity within and between populations (Bentley et al. 1995, Sivaramakrishnan et al. 2002b, Belabid et al. 2004, Nourollahi et al. 2011,). In this study, multiple alleles were recorded in F. oxysporum f. sp. phaseoli populations, as reported in the similar studies on Fusarium species (Mohammadi et al. 2004, Nourollahi & Jalali 2016, Nourollahi & Aliaran 2017). The polymorphic features of SSRs often allow characterization of fungi at strain level (Barres et al. 2006). The variability of alleles per loci showed the level of polymorphism and directly used for showing of genetic evolution (Sanders 2002, Mwang Òmbe et al. 2007, Nourollahi & Jalali 2016). The results showed that there was a low level of genetic variation in the overall isolates of F. oxysporum f. sp. phaseoli in Ilam province. In these regions, low levels genetic diversities have been reported in among F. oxysporum isolates collected from chickpea and lentils farms (Nourollahi & Jalali 2016, Nourollahi & Aliaran 2017). In such studies, the low amount of variation within the populations was perhaps due to exchange of contaminated seeds between sampling regions, infected plant debris. Cultivation of these seeds leads to distribution of spores of the fungus in the farms (Nourollahi & Jalali 2016).

Analysis of molecular variance (AMOVA) showed that the genetic variance between populations is 13 % in contrast to 87 % for the within population component. Pair-wise genetic distances between the populations were also very similar and they showed a different level of diversity in comparison with the local and international isolates. Similar results were also reported in F. oxysporum f. sp. lentis (Nourollahi & Jalali 2016), F. oxysporum f. sp. phaseoli (Woo et al. 1996), F. oxysporum f. sp. ciceris isolates (Jimenez-Gasco et al. 2001, Nourollahi & Aliaran 2017) and in the Ethiopian F. oxysporum isolates using AFLP, SSR, and ITS (Bogale et al. 2006).

The low level of gene diversity (Gst = 0.162) was distinguished in these populations. The low Gst value directed little genetic variation among the populations and showed little indication for geographical subdivision (Bayraktar et al. 2010). The gene flow (Nm = 7.589) indices revealed high gene flow and low genetic separation in these regions. This study suggests that the gene flow between sampling area played key role in permanence of the fungus and the similarity of F. oxysporum populations. This subject confirmed by a study that showed F. oxysporum f. sp. ciceris is a monophyletic group (Jimenez-Gasco et al. 2003).

Gene flow is one of the evolutionary elements that can have a vital role in genetic variability of populations. In the absence of gene flow, genetic drift causes emerging different allele frequencies at neutral loci, leading to variation in populations (Keller et al. 1997). Knowledge of the occurrence severity, distribution, and genetic relatedness of F. oxysporum f. sp. phaseoli isolates is necessary for developing effective integrated disease management.

In this study, genetic structure of F. oxysporum populations was recognized and the result can provide essential information for plant breeders to introduce tolerant cultivars against Fusarium and developing integrated strategies for disease management.

 

ACKNOWLEDGEMENTS

The financial assistance provided by Ilam University, is gratefully acknowledged.

 

 

Abbas A. 1995. Variation in some cultural and physiological character sand host/pathogen interaction of Fusarium oxysporum f. sp. lentis, and inheritance of resistance to lentil wilt in Syria. Ph.D. Thesis, Faculty of Agriculture, University of Aleppo, Syria.
Alves-Santos FM, Benito EP, Eslava AP, Diaz-Minguez JM. 1999. Genetic diversity of Fusarium oxysporum from common bean fields in Spain. Applied Environmental Microbiology 65: 3335-3340.
Alves-Santos FM, Ramos B, Garcia-Sanchez MA, Eslava AP, Diaz-Minguez JM. 2002. A DNA based procedure for in planta detection of Fusarium oxysporum f. sp. phaseoli. Phytopathology 92: 237-244.
Armstrong GM, Armstrong JK. 1981. Formae speciales and races of Fusarium oxysporum causing wilt diseases, p. 391–399. In P. E. Nelson, T. A. Toussoun, and R. Cook (ed.), Fusarium: diseases, biology and taxonomy. The Pennsylvania State University Press, USA.
Barres B, Dutech C, Andrieux A, Caron H, Pinon J, Frey P. 2006. Isolation and characterization of 15 microsatellite loci in the poplar rust fungus, Melampsora larici-populina, and cross-amplification in related species. Molecular Ecology Note 6: 60–64.
Barve MP, Haware MP, Sainani MN, Ranjekar PK, Gupta VS. 2001. Potential of microsatellites to distinguish four races of Fusarium oxysporum f. sp. ciceris prevalent in India. Theoretical Applied Genetic 102: 138–147.
Bayraktar H, Dolar FS. 2008. Genetic diversity of wilt and root rot pathogens of chickpea. Turkey Journal of Agriculture 33: 1–10.
Bayraktar H, Türkkan M, Dolar FS. 2010. Characterization of Fusarium oxysporum f. sp. cepae from onion in Turkey based on vegetative compatibility and rDNA RFLP analysis. Phytopathology 10:1439–0434.
Belabid L, Baum M, Fortas Z, Bouzand Z, Eujal I. 2004. Pathogenic and genetic characterization of Algerian isolates if Fusarium oxysporum f. sp. lentis by RAPD and AFLP analysis. African Journal of Biotechnology 3: 25–31.
Bentley S, Pegg KG, Dale JL .1995. Genetic variation among a worldwide collection of isolates of Fusarium oxysporum f. sp. cubense analyzed by RAPD-PCR fingerprinting. Mycology Research 99: 1378–1384.
Bogale M, Wingfield BD, Wingfield MJ, Steenkamp ET.  2005. Simple sequence repeats markers for species in the Fusarium oxysporum complex. Molecular Ecology Notes, 5: 622-624.
Bogale M, Wingfield BD, Wingfield MJ, Steenkamp ET. 2006. Characterization of Fusarium oxysporum isolates from Ethiopia using AFLP, SSR and DNA sequence analyses. Fungal Diversity 23: 51–66.
Britz H, Coutinho T A, Wingfield BD, Wingfield MJ. 2002. Sequence characterized amplified polymorphic markers for the pitch canker pathogen Fusarium circinatum. Molecular Ecology Notes 3: 577–580.
Burgess T, Wing M, Wing B. 2001. Simple sequence repeat markers distinguish among morphotypes of Sphaeropsis sapinea. Applied Environment Microbiology 67: 354–362
Doyle JJ, Doyle JL. 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Photochemistry Bulletin 19: 11– 15.
Edel V, Steinberg C, Gautheron N, Alabouvette C. 1997. Populations of nonpathogenic Fusarium oxysporum associated with roots of four plant species compared to soil borne populations. Phytopathology 87: 693–697.
Enjalbert J, Duan X, Leconte M, Hovmøller MS, De Vallavielle-Pope C. 2005. Genetic evidence of local adaptation of wheat yellow rust (Puccinia striiformis f. sp. tritici) within France. Molecular Ecology 14: 2065- 2073.
Haware MP, Nene YL, Natarajan M. 1996. Survival of Fusarium oxysporum f. sp. ciceris in soil absence of chickpea. Phytopathologia Mediterranea 35: 9–12.
Hawker LE. 1950. Physiology of fungi. University of London Press Ltd, London, UK.
Jalali BL, Chand H. 1992. Disease of cereals and pulses, Chickpea wilt Plant disease of International importance Vol I. In: Singh US, Mukhtopadhyay AN, Kumar J, Chaube HS (eds) Prentice Hall Englewood, Cliff NJ pp 420–444.
Jimenez-Gasco MM, Perez-Artes E, Jimenez-Diaz RM. 2001. Identification of pathogenic races 0, 1B/C, 5 and 6 of Fusarium oxysporum f. sp. ciceris with random amplified polymorphic DNA (RAPD). European Journal of Plant Pathology 107: 237-248.
Jimenez-Gasco MM, Jimenez-Diaz RM. 2003. Development of a specific polymerase reaction-based assay for the identification of Fusarium oxysporum f. sp. ciceris and its pathogenic races 0, 1A, 5 and 6. Phytopathology 93: 200-209.
Karimian B, Javan-Nikkhah M, Abbasi M, Ghazanfari K. 2010. Genetic diversity of Fusarium oxysporum isolates from common bean and distribution of mating type alleles. Iranian Journal of Biotechnology 8. 2
Keller SM, McDermott JM, Pettway RE, Wolfe MS, McDonald BA. 1997. Gene flow and sexual reproduction in the wheat glume blotch pathogen Phaeosphaeria nodorum (anamorph Stagonospora nodorum). Phytopathology 87: 353-358.
Kiprop EK, Baudoin JP, Mwang’ombe AW, Kimani PM, Mergeai G, Maquet A. 2002. Characterization of Kenyan isolates of Fusarium udum from pigeon pea [Cajanus cajan (L.) Millsp.] By cultural characteristics, aggressiveness and AFLP analysis. Phytopathology 150: 517-525.
Kistler HC. 2001. Evolution in host specificity in Fusarium oxysporum. In: Summerell BA, Leslie JF, Backhouse D, Bryden WL, Burgess LW (eds) Fusarium. APS Press, USA.
Maina PK, Wachira PM, Okoth SA, Kimenju JW. 2017. Cultural Morphological and Pathogenic Variability among Fusarium oxysporum f. sp. phaseoli Causing Wilt in French bean (Phaseolus vulgaris L.). Journal of Advances in Microbiology 2: 1- 9.
McDonald BA. 1997. The population genetic of fungi: tools and techniques. Phytopathology 87: 448–453.
Mohammadi M, Aminipour M, Banhashemi Z. 2004. Isozyme analysis and soluble mycelial protein pattern in Iranian isolates of several formae speciales of Fusarium oxysporum. Phytopathology 152: 267–276.
Mwang`Ombe AW, Thiong G, Olubayo FM, Kiprop EK. 2007. DNA microsatellite analysis of Kenyan isolates of Rhizoctonia solani from common bean (Phaseolus vulgaris L.). Plant Pathology Journal 6: 66–71.
Namvar Hamzanlouyi H. 2006. Identification of vegetative compatibility groups (VCGs) of Fusarium oxysporum f. sp. phaseoli in northern and Razavi Khorasan provinces. Iranian Plant Protection Congress, Karaj, Iran. PP 128.
Nei M. 1973. Analysis of the genetic diversity in subdivided populations. Proceedings of the National Academy of Science USA 70: 3321–3323.
Nei M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583–590.
Nei M, Li WH. 1979. Mathematical model for studying genetic variation in terms of restriction endonuclease. Proceeding of National Academy Science USA 76: 5269–5273.
Nene YL, Haware MP. 1980. Screening chickpea for resistance to wilt. Plant Disease 64: 379 –380.
Nourollahi K, Javan–nikkhah M, Naghavi MR, Lichtenzveig J, Okhovat SM, Oliver RP, Ellwood SR. 2011. Genetic diversity and population structure of Ascochyta rabiei from the western Iranian Ilam and Kermanshah provinces using MAT and SSR markers. Mycological Progress 10: 1–7.
Nourollahi K, Madahjalali M. 2017. Analysis of population genetic structure of Iranian Fusarium oxysporum f. sp. lentis isolates using microsatellite markers. Australasian Plant Pathology 46: 35–42.
Nourollahi K, Aliaran A. 2017. Genetic structure analysis of Fusarium oxysporum f. sp. ciceris populations from chickpea in Ilam province, Iran. Mycologia Iranica 4: 93-102.
Peakall R, Smouse PE. 2006. GenAlEx 6: genetic analysis in excel population genetic software for teaching and research. Molecular Ecology Notes 6:288–295.
Perveen K, Haseeb A, Shukla PK. 2007. Management of Sclerotinia sclerotiorum on Mentha arvensis cv. Gomti. Journal of Mycology and Plant Pathology 37:33-36.
Roldan–Ruiz I, Dendauw JE, Van bockstaele E, Depicker A, Loose M. 2000. AFLP markers reveal high polymorphic rates in rye grasses (Lolium spp.). Molecular Breeding 6: 125–126.
Rolhf FJ. 1990. NTSYSPc Numerical taxonomy and multivariant analysis system Version 202 Applied Biostatics New York.
Sanders IR. 2002. Ecology and evolution of multigenomic arbuscular mycorrhizal fungi. Am Nat 160: 8– 41.
Sirjusingh C, Kohn LM. 2001. Characterization of microsatellite in the fungal plant pathogen Sclerotinia sclerotiorum.  Molecular Ecology Notes 1: 267–269.
Sivaramakrishnan S, Kannan S, Singh SD. 2002a. Detection of genetic variability in Fusarium udum using DNA markers. Indian Journal of Phytopathology 55: 258–263.
Sivaramakrishnan S, Kannan S, Singh SD. 2002b. Genetic variability of Fusarium wilt pathogen isolates of chickpea (Cicer arietinum L). Assessed by molecular markers. Mycopathology 155: 171–178.
Slippers B, Burgess T, Wingfield BD, Crous PW, Coutinho TA, Wingfield MJ. 2004. Development of simple sequence repeat markers for Botryosphaeria spp. with Fusicoccum anamorphs. Molecular Ecology 4: 675- 677
Sneath PHA, Sokal RR. 1973. Numerical taxonomy: The principles and practice of numerical classification. WH Freeman and Company San Francisco, USA.
Stewart JE, Kim M, James RL, Dumroese RK, Klopfenstein NB. 2006. Molecular characterization of Fusarium oxysporum and Fusarium commune isolates from a conifer nursery. Phytopathology 96:1124–1133.
Tenzer I, Ivanissevich SD, Morgante M, Gessler C. 1999. Identification Of microsatellitemarkers and their application to population genetics of Venturia inaequalis. Phytopathology 89: 748–753.
Varshney RK, Chabane K, Hendre PS, Aggarwal RK, Graner A. 2007. Comparative assessment of EST-SSR; EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild cultivated and elite barleys. Plant Science 173: 638–649.
Woo SL, Zoina A, Sorbo G, Lorito MD, Scala NBF, Noviello C. 1996. Characterization of Fusarium oxysporum f. sp. phaseoli by pathogenic races VCGs RFLPs and RAPD. Phytopathology 86:966–973.
Yu F, Yang R, Boyle T. 1999. POPGENE version 1.31 Microsoft windows-based software for population genetics analysis. University of Alberta and Centre for International Forestry Research Alberta.