Plant Pathol J > Volume 36(1); 2020 > Article
Mahmoud and Abd El-Fatah: Genetic Diversity Studies and Identification of Molecular and Biochemical Markers Associated with Fusarium Wilt Resistance in Cultivated Faba Bean (Vicia faba)

Abstract

Faba bean (Vicia faba L.) is one of the most important legume crops in Egypt. However, production of faba bean is affected by several diseases including fungal diseases. Fusarium wilt incited by Fusarium oxysporum Schlecht. was shown to be the most common wilt disease of faba bean in Assiut Governorate. Evaluation of 16 faba bean genotypes for the resistance to Fusarium wilt was carried out under greenhouse conditions. Three molecular marker systems (inter-simple sequence repeat [ISSR], sequence related amplified polymorphism [SRAP], and simple sequence repeat [SSR]) and a biochemical marker (protein profiles) were used to study the genetic diversity and detect molecular and biochemical markers associated with Fusarium wilt resistance in the tested genotypes. The results showed that certain genotypes of faba bean were resistant to Fusarium wilt, while most of the genotypes were highly susceptible. The percentage of disease severity ranged from 32.83% in Assiut-215 to 64.17% in Misr-3. The genotypes Assiut-215, Roomy-3, Marut-2, and Giza-2 were the most resistant, and the genotypes Misr-3, Misr-1, Assiut-143, Giza-40, and Roomy-80 performed as highly susceptible. The genotypes Assiut-215 and Roomy-3 were considered as promising sources of the resistance to Fusarium wilt. SRAP markers showed higher polymorphism (82.53%) compared with SSR (76.85%), ISSR markers (62.24%), and protein profile (31.82%). Specific molecular and biochemical markers associated with Fusarium wilt resistance were identified. The dendrogram based on combined data of molecular and biochemical markers grouped the 16 faba bean genotypes into three clusters. Cluster I included resistant genotypes, cluster II comprised all moderate genotypes and cluster III contained highly susceptible genotypes.

Faba bean (Vicia faba L.) is the most common fast food item in the Egyptian diet. It is an excellent source of protein and complex carbohydrates, dietary fiber, choline, lecithin, minerals and secondary metabolites (Ali et al., 2019; Rabey et al., 1992). Faba bean is subjected to attack by numerous pathogenic fungi wherever the crop is grown. Plant pathogenic fungi decrease crop productivity and lower the quality of seeds; it attacks the roots and stem base causing serious losses in seed germination and plant stand. Wilt disease incited by Fusarium oxysporum Schlecht. is one of the most important diseases affecting faba bean in the major growing areas of Egypt and causing serious losses on faba bean production (Abdul Wahid et al., 1998; Abou-Zeid et al., 1997; Ahmed et al., 2002; Elwakil et al., 2009; Mahmoud, 2016a; Mazen et al., 2008). Fusarium wilt (FW) of faba bean caused by F. oxysporum was first reported to occur in Egypt by Abdel Rehim (1962). It was found to be the causal pathogen of wilt of faba bean plants in many countries worldwide (Dong et al., 2014).
Cultivation of the resistant genotypes is the most effective method for integrated management of the disease and remarkable reduced the incidence of the wilt disease on faba bean. In this investigation, we aimed to detect some biochemical and/or molecular markers related to resistant genes for FW disease in some genotypes of faba bean. Therefore, these markers can be used in the early selection of resistant genotypes to provide the effort and cost resulting from the use of traditional breeding methods. Genetic variation among faba bean genotypes is imperative for their efficient utilization in plant breeding programs and effective conservation. Although traditional methods of screening plants for FW resistance are effective, they are time-consuming, expensive and affected by environmental factors. In such state of affairs, indirect selection using biochemical and/or molecular markers associated with resistance gene(s) could be an efficient method in marker-assisted selection to conquer the inexactness in the field evaluation (Tanksley et al., 1989).
Development of molecular markers provided reliable information for evaluating and assessing the genetic diversity of faba bean genotypes; random amplified polymorphic DNA (RAPD) (Alghamdi, 2009b; Yahia et al., 2014), inter-simple sequence repeat (ISSR) (Alghamdi et al., 2011; Terzopoulos and Bebeli, 2008), amplified fragment length polymorphism (Zeid et al., 2003), simple sequence repeat (SSR) (Tekalign et al., 2019; Tufan and Erdoğan, 2017; Yahia et al., 2014; Yang et al., 2012), and sequence related amplified polymorphism (SRAP) (Alghamdi et al., 2012; Mahmoud, 2016b). Molecular markers are sufficient in numbers, not influenced by environmental factors or by development stages (Bebeli and Kaltsikes, 1993).
In the present study, three molecular marker systems (ISSR, SRAP, and SSR markers) have been used to study genetic diversity and identify markers associated with FW resistance in tested faba bean genotypes.
These markers differ in technical principle, type of inheritance, reproducibility, distribution in the plant genome and amount of polymorphism (Li and Quiros, 2001; Provan et al., 1999; Zietkiewicz et al., 1994).
Okiror (2002) reported that the genetic basis of FW resistance is still not clear, and numbers of genes postulated to be involved vary from a single dominant gene to two complementary genes and even involvement of multiple factors. So, advent of genomic tools especially molecular markers has facilitated breeding in many crops leading to development of several improved cultivars/varieties with enhanced resistance/tolerance to biotic or abiotic stresses (Varshney et al., 2006).
SSR marker technique based on the use of microsatellite sequences as primers for amplification in the polymerase chain reaction (PCR) can be identified polymorphism around the disease-resistant (Zietkiewicz et al., 1994). The RAPD markers have been used successfully to tag various resistant genes in many crops like downy mildew resistance gene in broccoli (Giovannelli et al., 2002), powdery mildew resistance gene in pea (Srivastava et al., 2012), and YVMV resistance gene in mung bean (Dhole and Reddy, 2013). The ISSR technique has been used in tagging disease resistance genes in a number of crops (Reddy et al., 2002; Souframanien and Gopalakrishna, 2006).
Electrophoretic techniques present an exceptional opportunity to evaluate and study the protein structure among different genotypes.
The protein profiling of germplasm and genetic markers have been widely and effectively used to determine the taxonomic and evolutionary aspects of several crops, Vigna mungo (Ghafoor et al., 2002), Pigeonpea (Nagaraja et al., 2016), and Vicia faba (Abdel-Razzak et al., 2012; Alghamdi, 2009a).
Isozymes, peroxidase, polyphenol oxidase, superoxide dismutase, and pathogenesis-related proteins, such as chitinase and β-1,3-glucanase, have very important roles in resistance to plant pathogenic fungi (Prasad et al., 2003). Natural resistance is genetically controlled, and resistance genes are expressed through the production of biological products such as secondary metabolites, including phenolics, lignin, suberin, phytoalexins, alkaloids, terpenes, glycosides, degrading enzymes, and pathogenesis-related proteins, all of which contribute to disease resistance (Gupta et al., 2010).
We propose that the selection and propagation of resistant cultivars of faba bean could be an effective means for the management of FW disease and minimize the damage to this important crop.
The aims of this investigation were to (1) evaluate the response of different genotypes of faba bean to the wilt disease incited by Fusarium oxysporum; (2) detect the resistant faba bean genotypes for FW; (3) differentiate between 16 faba bean genotypes by ISSR, SRAP, and protein profiles; and (4) detect ISSR, SRAP, SSR, and protein markers associated with the FW resistance.

Materials and Methods

Source of isolates

Highly pathogenic isolate of Fusarium oxysporum (No. 10, which showed 73.75% disease severity according to the results of the previous study of Mahmoud, 2016a) was utilized in this study. The pathogen was previously isolated from naturally diseased faba bean plants, showing wilt symptoms. The diseased samples were collected from farmer’s fields in Assiut Governorate, Egypt. The identification of the isolates was carried out based on the morphological features of mycelia and spores as well as cultural characteristics, the obtained isolates were identified as Fusarium oxysporum Schlecht. emend. Snyder & Hansen. according to Booth (1971), Gerlach and Nirenberg (1982), Kleczewski and Egel (2011) and Snyder and Hansen (1940).

Greenhouse experiment

Reactions of 16 faba bean genotypes were evaluated. The experiment was carried out in the greenhouse of Plant Pathology Department, Faculty of Agriculture, Assiut University, Egypt during two seasons in 2016-2017 and 2017-2018. Inoculums of Fusarium oxysporum was prepared by transferring agar pieces cut from the fungal culture into 500-ml Erlenmeyer flasks containing barley medium (100 g barley, supplemented with 2 g glucose + 1 g yeast extract + 100 ml distilled water) and incubated at 25 ± 2°C for 15 days. Sterilized pots (25 cm in diameter) were filled with sterilized sandy-loam soil which mixed thoroughly with equal amounts of F. oxysporum inoculums at the ratio of 2% soil weight, mixed well and thoroughly irrigated. Soil infestation was carried out one day before sowing seeds. Each pot was planted with 10 surface-sterilized seeds of faba bean. Seeds were surface sterilized by immersion in 0.5% solution of sodium hypochlorite for 2-3 min, then the seeds were removed, washed three times in distilled sterilized water, and air dried on the top layer of two sterilized filter paper sheets before sowing. Pots containing 2% non-infested barley medium were used as control. Four pots were used as replicates and randomly distributed on the greenhouse benches. The prevailing temperatures during the infection test were 15 ± 2°C (minimum) and 22 ± 2°C (maximum).

Disease assessment

Disease severity was estimated after 10 weeks from planting date, as a wilting percentage, on the basis of root discoloration or leaf yellowing. Plants with typical FW symptoms were assessed according to the type of symptoms that were observed using numerical grades ranging from 0 to 5 as follows (Mahmoud, 2016a): 0, no visible symptoms; 1, 1 to <20% of plant leaves are yellow and of the vascular systems are light brown (discoloration); 2, 20 to <40% of plant leaves are yellow and of the vascular systems are brown (discoloration); 3, 40 to <60% of plant leaves are yellow and of the vascular systems are dark brown (discoloration); 4, 60 to <80% of plant leaves are yellow and of the vascular systems are dark brown (discoloration); 5, 80 to ≤100% of plant leaves are yellow and of the vascular systems are dark brown (discoloration) or completely dead plants.
Disease severity (%)=Σ[(N×0)+(N×1)+(N×5)]5T×100
Where: N = the number of plants corresponding to the numerical grade, 0, 1, 2, 3, 4, and 5; 5T = the total number of plants (T) multiplied by maximum numerical grade (5).

Statistical analysis

The results were analyzed using ANOVA test and the means differences were regarded as significant using least significant difference test at 5% level of probability according to Gomez and Gomez (1984).

Molecular and biochemical analysis

DNA extraction

Total genomic DNA was extracted from young leaves of five plants each genotype following the CTAB method described by Murray and Thompson (1980) with minor modification of the extraction buffer (2% CTAB, 100 mM Tris-HCl pH 8.0, 20 mM EDTA, 1.4 M NaCl, and 0.06% β-mercaptoethanol added just before use). The quantity and quality of the DNA were measured using spectrophotometer at 260-280 nm and checked on 1% agarose gel.

ISSR, SRAP, and SSR analyses

Fifteen ISSR primers, 20 SRAP primer combinations, and 10 SSR primer pairs, obtained from Metabion International AG Company (metabion international AG/metabion GmbH, München, Germany), were tested in this investigation to amplify the template DNA (Supplementary Table 1).
Amplification reactions were carried out in a 25 μl volumes, containing 1× reaction buffer, 0.2 mM of each dNTP, 2.0 mM MgCl2, 0.2 μM for (ISSR primer) or 0.3 μM of each of forward and reverse primer for (SRAP and/or SSR), 1.0 units of Taq DNA polymerase and 40 ng of genomic DNA template. Amplification conditions were carried out in a Lab Cycler (SensoQuest GmbH, Göttingen, Germany) with the following specification:
  • For ISSR; initial denaturation for 5 min at 94°C, 40 cycles of 1 min at 92°C, 1 min at 44°C and 2 min at 72°C, and final extension step for 10 min at 72°C.

  • For SRAP; 5 min at 94°C followed by 10 cycles: 1 min at 94°C, 1 min at 35°C and 2 min at 72°C then 35 cycles with 1 min at 94°C, 1 min at 50°C and 2 min at 72°C then finally extension for 10 min at 72°C.

  • For SSR; 5 min at 94°C (initial denaturation), 40 cycles of 1 min at 92°C, 1 min at 44°C and 2 min at 72°C, and final extension step for 10 min at 72°C.

The PCR products were separated in an Ultra-Pure agarose gel with different concentrations, 2% ISSR and 2.5% SRAP and SSR, at 80 volts for 3-3.5 h. Gels were stained with ethidium bromide (0.5 μg/ml) and DNA fragments were visualized using GelDoc-It 2 Imager (LTF Labortechnik GmbH & Co. KG, Wasserburg, Germany).

Data scoring and statistical analysis

The ISSR, SRAP, and SSR bands were scored as present (1) or absent (0) across 16 genotypes for each primer or primer combination. The three parameters viz., polymorphic information content (PIC), marker index (MI), and resolving power (RP), were calculated according to Ghislain et al. (1999), Powell et al. (1996), and Prevost and Wilkinson (1999), respectively.
A dendrogram was constructed based on Dice similarity coefficients (GS) (Dice, 1945) using the unweighted pair group method with arithmetic average (UPGMA) with the SAHN module of NTSYS-pc 2.20 (Rohlf, 2000).

Assessment of gene expression

Total soluble proteins and enzymes were extracted from infected and non-infected plants to study the changes in gene expression in resistant and sensitive genotypes response to infection with FW. Protein was extracted by crushing 1.0 g of sample tissue in 1.0 ml extraction buffer (0.1 M Tris-HCl + 2 mM EDTA, pH 7.8).

Protein analyses

Electrophoresis of protein patterns was separated on sodium dodecyl sulfate polyacrylamide gels (12%) according to Laemmli (1970) method. Coomassie Blue R stain was used to detect protein bands. Then, gels were destained by repeated immersions in a mixture of methanol: acetic acid: water (1:1:8, by volume). The molecular weight of protein bands was determined against protein markers which consisted of 175, 80, 58.46, 30, 27, 17, and 7.0 kD using GS 365 electrophoresis data system program version 3.01 (Microsoft Windows version).

Isozymes analysis

The gels were stained for two enzyme systems. The staining protocol for peroxidase (E.C. 1.11.1.7) was conducted according to Guikema and Sherman (1981). While, chitinase (E.C. 3.2.1.14) was stained according to the protocol described by Shen et al. (1991). Stained gels were placed in a light box to determine their isozyme banding patterns (BP). The isozyme band numbers were recorded, and their relative mobilities (RF) were calculated according to the following formula:
RF=Distance traveled by the bandDistance traveled by the tracking dye
Finally, the correlation between the different molecular marker systems as well as between molecular markers and biochemical markers were calculated using Mantel test (Mantel, 1967).

Results

Reactions of faba bean genotypes for resistance to FW disease under greenhouse conditions

Faba bean production areas in Egypt adversely affected by FW incited by F. oxysporum. The adult plants become severely diseased at flowering to late pod-filling stage, resulting in dull green foliage color followed by wilting of whole plant or individual branches, and lower yields or death of affected plant. This study was determined to evaluate reaction of 16 faba bean genotypes against F. oxysporum, under greenhouse conditions during two seasons in 2016-2017 and 2017-2018. Results presented in Table 1 showed that all tested genotypes manifested visible symptoms of FW. Susceptible genotypes revealed remarkable significant differences than the resistant genotypes. This result indicated the importance to cultivate the resistant genotypes of faba bean that weak affected by F. oxysporum.
The percentage of wilt disease severity ranged from 32.83% for the resistant genotype Assiut-215 to 64.17% for the highly susceptible genotype Misr-3. None of the faba bean genotypes presented reaction similar to immunity to F. oxysporum and only four genotypes (Assiut-215, Roomy-3, Marut-2, and Giza-2) performed as resistant genotypes. Most of the tested genotypes were moderately resistant, while five genotypes (Misr-3, Misr-1, Assiut-143, Giza-40, and Roomy-80) were highly susceptible, with disease severity ranged from 60.43% to 64.17%. No symptoms were observed in non-inoculated plants (control).

Molecular analysis

Molecular polymorphism among 16 faba bean genotypes was assessed using three molecular marker systems including 45 primers. Out of 45 primers or primer pairs, 17 (37.78%) (6 ISSR, 6 SRAP, and 5 SSR) were polymorphic among all genotypes and were used for further analysis. In our study, the monomorphic primers or primer combinations were excluded from the analysis.

ISSR analysis

A total of 53 scored bands were generated using the ISSR primers previously screened across 16 faba bean genotypes (Fig. 1, Supplementary Table 2). The number of amplified bands per primer ranged from seven (HB13, HB06) to 12 (HB09) with an average of 8.83 bands per primer (Table 2). Of the 53 amplified bands, 33 (62.24%) were polymorphic, ranging from 3 (HB14, HB12) to 9 (HB09), with a mean of 5.5 polymorphic bands per primer. The percentage of polymorphic bands across the primers ranged from 33.33% (HB14, HB12) to 88.89% (HB10), with an average of 62.24% polymorphism. The primer HB10 exhibited the highest values (0.3, 2.42, and 3.88), while the primer HB14 showed the lowest values (0.1, 0.31, and 1.25) for the three primer parameters (PIC, MI, and RP, respectively) (Table 2).
GS coefficients obtained with UPGMA algorithm using Dice coefficient (Dice, 1945) ranged from 0.754 (Romy-3 and Misr-1) to 0.959 (Sakha-1 and Sakha-4). The dendrogram grouped the 16 genotypes into three main clusters (with a GS of 0.825) (Fig. 1).

SRAP analysis

A total of 76 bands were observed using SRAP primer combinations (Fig. 2, Supplementary Table 3), of which 62 were polymorphic (82.53%), ranging from 5 (SRAP-5) to 13 (SRAP-3, SRAP-4) (Table 2). Average number of bands and polymorphic bands per primer combinations were 12.67 and 10.33, respectively. Percentage of polymorphic bands (PPB) ranged from 60% (SRAP-1) to 100% (SRAP-4) with an average polymorphism of 82.53% across all the genotypes (Table 2). The two SRAP primer combinations, SRAP-3 and SRAP-4 showed the highest PIC and MI values (0.28, 3.66 and 0.28, 3.64, respectively), while the two primer combinations, SRAP-3 and SRAP-6 exhibited the highest values for RP parameter (5.25 and 5.5, respectively) (Table 2).
The dendrogram using UPGMA analysis based on Dice GS of the SRAP markers of 16 genotypes is shown in Fig. 2. The GS varied from 0.63 (Misr-3 and Sakha-4) to 0.923 (Roomy-3 and Mariut-2). The 16 faba bean genotypes were grouped into five clusters (with a GS of 0.687) (Fig. 2).

SSR analysis

A total of 48 bands were amplified using SSR primer pairs (Fig. 3, Supplementary Table 4), of which 37 were polymorphic (76.85%), ranging from 5 (SSR-4) to 11 (SSR-2) (Table 2). Average number of bands and polymorphic bands per primer pairs were 9.6 and 7.4, respectively. PPB ranged from 63.64% (SSR-5) to 91.67% (SSR-2) with an average polymorphism of 76.85% across all the genotypes (Table 2). The SSR-2 primer pairs showed the highest mean values for three genetic parameters, PIC, MI, and RP (0.25, 2.74, and 4.13, respectively) (Table 2).
The dendrogram using UPGMA analysis based on Dice GS of the SSR markers of 16 genotypes is shown in Fig. 3. The GS varied from 0.721 (Giza-40 and Assiut-215) to 0.947 (Misr-1 and Misr-3). The 16 faba bean genotypes were grouped into four clusters (with a GS of 0.797) (Fig. 3).

Combined ISSR, SRAP, and SSR data

To obtain more accurate genetic estimates, a combined analysis was carried out using all the ISSR, SRAP, and SSR data.
A total of 177 DNA fragments were obtained using ISSR, SRAP, and SSR markers used to screen 16 faba bean genotypes, with an average of 10.41 bands/primer. Out of 177 fragments, 132 (73.69%) showed polymorphism. The two genotypes Romy-3 and Giza-2 displayed the highest number of DNA fragments (126 bands) followed by Assiut-215 and Mariut-2 (125 and 122 bands, respectively), while the genotype Misr-3 revealed the lowest number of bands (93 bands).
The three molecular marker systems were successful in characterizing 13 out of the 16 faba bean genotypes by generating unique positive and/or negative markers (Table 3). The two genotypes Assiut-215 and Sakha-4 were distinguished by four unique bands for each, while the two genotype Nobaria-3 and Giza-843 had the highest positive unique bands (2 bands for each).
The Dice GS ranged from 0.725 (Misr-3 and Assiut-215) to 0.919 (Roomy-3 and Mariut-2) (Supplementary Table 5). The Dice GS matrix based on combined data was used to generate a dendrogram showing genetic relationships among the genotypes. The dendrogram grouped the 16 genotypes into three main clusters (Fig. 4). Cluster I included four genotypes (Assiut-215, Roomy-3, Mariut-2, and Giza-2) with GS 0.889, which showed resistance reaction for FW. Cluster II was the largest cluster consisting of seven genotypes and was further divided into three subclusters with a GS of 0.85. All genotypes in this cluster showed moderate resistance to the FW. Cluster III comprised of five genotypes (Misr-3, Misr-1, Assiut-143 Giza-40, and Roomy-80), with a GS of 0.81, which were highly susceptible to the FW.
The Mantel test performed between Dice similarity matrices obtained using ISSR, SRAP, and SSR markers showed significant correlations between ISSR + SRAP (r = 0.282), ISSR + SSR (r = 0.348), SSR + SRAP (r = 0.418), and SRAP + (ISSR + SSR) (r = 0.338), which indicated good agreement among ISSR, SRAP, and SSR markers.
Among the polymorphic bands, 10 DNA fragments, 457 bp (HB10), 209 bp (HB14), 312 bp (HB12), 219 bp (HB06), 984 bp (SRAP-1), 720 bp (SRAP-3), 401 bp (SRAP-4), 291 bp (SRAP-6), 419 bp (SSR-1), and 230 bp (SSR-3) were present in resistance faba bean genotypes, but absent in moderate and susceptible genotypes. The DNA fragment at molecular size 401 bp (SRAP-5) was present in moderate and susceptible genotypes, but absent in resistant genotypes. Eight DNA fragments, 498 bp (HB14), 289 bp (SRAP-2), 319 bp (SRAP-4), 561 bp (SRAP-5) and 656 bp, 593 bp and 146 bp (SRAP-6) and 498 bp (SSR-1) were present in all moderate and resistant genotypes, but absent in susceptible genotypes. These findings suggested that these fragments could be served as markers associated with the resistance to FW in faba bean.
These specific bands putatively associated with the FW resistance were subsequently screened against the highest resistant genotype (Assiut-215) and highly susceptible genotype (Misr-3). In general, 18 fragments, 457 bp (HB10), 209 bp (HB14), 312 bp (HB12), 219 bp (HB06), 984 bp (SRAP-1), 720 bp (SRAP-3), 401 bp (SRAP-4), 291bp (SRAP-6), 419 bp (SSR-1), 230 bp (SSR-3), 498 bp (HB14), 289 bp (SRAP-2), 319 bp (SRAP-4), 561 bp (SRAP-5) and 656 bp, 593 bp and 146 bp (SRAP-6) and 498 bp (SSR-1) were present only in resistant genotype, but absent in highly susceptible genotype. On the other hand, the DNA fragment at molecular size 401 bp (SRAP-5) was observed in susceptible genotypes but absent in resistant genotypes. To confirm the association of ISSR, SRAP, and SSR markers tested with the FW resistance, 30 individual plants of the resistant genotype (Assiut-215) and 30 individual plants of the highly susceptible genotype (Misr-3) were screened.
In general, the DNA fragments, 457 bp (HB10), 209 bp (HB14), 312 bp (HB12), 219 bp (HB06), 984 bp (SRAP-1), 720 bp (SRAP-3), 401 bp (SRAP-4), 291 bp (SRAP-6), 419 bp (SSR-1), and 230 bp (SSR-3) were found only in 29 (96.67%) of all individual resistant plants, but not observe in all susceptible plants. Eight DNA fragments, 498 bp (HB14), 289 bp (SRAP-2), 319 bp (SRAP-4), 561 bp (SRAP-5) and 656 bp, 593 bp and 146 bp (SRAP-6) and 498 bp (SSR-1) were generated in 28 (93.33%) of resistant plants while not amplified in susceptible plants. The DNA fragment at molecular size 401 bp (SRAP-5) was not observed in all resistant plants, while this fragment was generated in all susceptible plants. These findings confirmed that these bands could be considered as specific markers associated with the FW resistance in these genotypes.

Biochemical markers

Protein profiling

Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) of total soluble protein was carried out for the 16 faba bean genotypes uninfected (control) and infected with FW (Fig. 5A and B, Supplementary Table 6). For control; electrophoretic BPs produced a total of twenty-two total protein bands, ranged in sizes from 7.3 kD to 40.1 kD as shown in Fig. 5A. Out of 22 protein bands, seven bands were polymorphic (31.82%), while 15 (68.18%) were common bands for all the genotypes tested. The dendrogram using UPGMA analysis based on Dice GS of the protein profiles of 16 genotypes is shown in Fig. 5C. The GS varied from 0.842 to 1.0 (Supplementary Table 7). The 16 faba bean genotypes were grouped into three main clusters (with a GS of 0.687) (Fig. 5C). All resistant genotypes were grouped in cluster I, while all moderate and susceptible genotypes were separated in cluster II and cluster III.
Among polymorphic protein profiles, protein bands at molecular weights 30.1 kD and 10.2 kD were expressed only in all resistant genotypes, but not expressed in all moderate and susceptible genotypes. On the other hand, protein bands 21.4 kD was expressed in all moderate and resistant genotypes, but absent in highly susceptible genotypes.
For infected plants; a new protein band of molecular weight 26.2 kD induced in all faba bean genotypes infected with FW as compared with control plants (Fig. 5B). Also, new band of molecular weight 29.7 kD appeared only in all resistant genotypes infected with FW as compared with control plants. On the other hand, protein band of molecular weight 36.1 kD disappeared only in highly susceptible faba bean genotype infected with FW as compared with control plants. While protein band at 12.1 kD disappeared from all moderate and susceptible genotypes infected with FW as compared with control (Fig. 5B). Such protein bands, which appeared only in resistant genotypes infected and uninfected with FW, could be used as protein markers for FW resistant in faba bean genotypes.

Chitinase isozyme

Seven chitinase isozyme patterns were expressed across control and infected faba bean genotypes (Fig. 6A and B). The four chitinase isozyme bands at RF values of 0.46, 0.53, 0.67, and 0.72 were monomorphic among all control and infected plants, while the three bands at RF values of 0.42, 0.77, and 0.82 were polymorphic. A chitinase band at RF value of 0.42 was found in all control and infected genotypes except the five genotypes, Misr-3, Misr-1, Assiut-143, Giza-40, and Romy-80. A band at RF value of 0.77 expressed only in the infected genotypes, Assiut-215, Romy-3, Maruit-2, and Giza-2 and it was not present in control and other infected genotypes. The chitinase band at RF value 0.82 appeared only in infected genotypes, Assiut-215, Maruit-2, and Giza-2 and it was not present in control and other infected genotypes.

Peroxidase isozyme

Peroxidase patterns showed nine bands at different RF values varying from 0.14 to 0.86, with four bands being polymorphic, and the other five bands at RF values of 0.41, 0.63, 0.69, 0.74, and 0.76 were found to be monomorphic among all control and infected faba bean genotypes (Fig. 6C and D). The peroxidase band at RF value of 0.14 appeared in all control and infected genotypes except five infected faba bean genotypes, Misr-3, Misr-1, Assiut-143, Giza-40, and Romy-80. The peroxidase band at RF values of 0.34 was found in all control plants and four infected genotypes. A band at RF value of 0.47 only appeared in the infected plants and was not present in control plants. The peroxidase band at RF value 0.86 appeared only in infected genotypes, Assiut-215, Romy-3, Maruit-2, and Giza-2 and it was not present in control and other infected genotypes.

Combined analysis based on molecular and biochemical markers

Based on combined analysis of ISSR, SRAP, SSR, protein patterns, and isozymes profiles, similarity matrix was developed by NTSYS-pc 2.20 computer package system (Rohlf, 2000). The Dice GS ranged from 0.74 (Assiut-215 and Giza-40) to 0.928 (Roomy-3 and Mariut-2) (Supplementary Table 8). The Dice GS matrix based on combined data was used to generate a dendrogram showing genetic relationships among the genotypes.
The dendrogram grouped the 16 genotypes into three main clusters (Fig. 6E). Cluster I included four genotypes (Assiut-215, Roomy-3, Mariut-2, and Giza-2), which showed resistant reaction for FW. Cluster II was the biggest cluster consisting of seven genotypes and was further divided into three sub-clusters with a GS of 0.868. All genotypes in this cluster showed moderate resistant to FW. The cluster III comprised five genotypes Misr-3, Misr-1, Assiut-143 Giza-40, and Roomy-80 with GS 0.836. These five genotypes in this cluster were highly susceptible for FW.
The Mantel test values between Dice similarity matrices of molecular markers showed significant correlations between ISSR + biochemical markers (r = 0.364), SRAP + biochemical markers (r = 0.451), SSR + biochemical markers (r = 0.611), and (ISSR + SRAP + SSR) + biochemical markers (r = 0.578) which indicated good agreement between biochemical and molecular markers.

Discussion

Faba bean is the most important food legume in Egypt and is widely grown in many areas of the world. FW of faba bean caused by F. oxysporum is an important vascular disease of the crop and considered as one of the major handicaps of production. It causes the plants to wilt and become less productive or even die. The cultivation of the resistant genotypes provides the most suitable method for control the disease and remarkable reduced the incidence of the wilt disease on faba bean. To select materials with potential for use in faba bean breeding programs and/or disease integrated management, 16 genotypes were evaluated for disease resistance in this study. The use of resistant cultivars is widely recognized as the safest, most economical and most effective method for protecting crops from disease (Johnson and Jellis, 1992).
Assessment of the genetic variability within a cultivated crop has important consequences in plant breeding and the conservation of genetic resources, and molecular and biochemical markers offer approaches to unveil the genetic diversity among genotypes based on polypeptide chains and nucleic acid polymorphisms. In the present study, protein profiles and three molecular maker systems, ISSR, SRAP, and SSR markers were applied to assess the level and pattern of genetic diversity of 16 faba bean genotypes. The results showed a wide genetic variability among different faba bean genotypes (Table 2). The polymorphism rate obtained using ISSR, SRAP, and SSR markers in the present study was higher than those obtained in other plant species including alfalfa (Vandemark et al., 2006), Paeonia (Hao et al., 2008), and buffalo grass (Budak et al., 2004). The high level of polymorphism indicates the genetic divergence among the studied genotypes, suggesting that these marker systems were ideal for the study of the genetic diversity among closely related genotypes. The applicability of ISSR, SRAP, and SSR as genetic markers to characterize the genetic diversity of faba bean was compared. SRAP markers showed higher PPB (82.53%) than the ISSR and SSR markers (62.24% and 76.85%, respectively). Moreover, the PIC, MI, and the RP were higher for SRAP than those obtained with ISSR and SSR markers (Table 2).
These results are relatively close with those of Terzopoulosa and Bebeli (2008) and Khalaf et al. (2015) they declared that the percentage of polymorphism within faba bean genotypes was 67.48% using ISSR primers. Similar results were also reported by Ammar et al. (2015) who found high polymorphism among faba bean genotypes using SRAP primers. These results also, in accordance with Tekalign et al. (2019) who study the genetic diversity among faba bean genotypes and found that the percentage of polymorphism was 68.5%.
The present study showed that the clustering of genotypes based on ISSR, SRAP, and SSR markers was slightly similar. The correlation between Dice similarity coefficients resulted in a significant correlation among ISSR, SRAP, and SSR data, indicated a good agreement among three molecular marker systems. Accordingly, Velasco-Ramírez et al. (2014) recommended that the information that is generated by various markers had a valuable capacity for genetic studies involving both, diversity and relationships. By this criterion, the combining different marker system was critically better for the study of the diversity as recommended through numerous reports (Amar and El Wahab, 2013; Costa et al., 2016; Jhang et al., 2010; Velasco-Ramírez et al., 2014; Wang et al., 2012).
Our results indicated the presence of a considerable genetic variability among tested faba bean genotypes. Genetic variation among faba bean genotypes based on ISSR, SRAP, and SSR analysis could be useful to select parents to be crossed for generating appropriate populations intended for both genome mapping and breeding purposes. The difference in resolution of three marker systems is that the three marker techniques targeted different sequences of the genome. The ISSR markers scattered throughout the genome which revealed the diversity of the entire genome, SRAP markers only amplified target sequence of open reading frame, while SSR marker target tandem sequence repeats of short DNA motifs (1-6 nucleotides long) (Li and Quiros, 2001; Provan et al., 1999; Zietkiewicz et al., 1994). Since each DNA maker system has its own advantages, it is important to use more than one DNA marker system, such as ISSR, SRAP, and SSR, in the analysis of genetic diversity. All three marker types were able to generate unique bands in particular genotypes but not in other genotypes. These unique bands could be used as positive and/or negative markers for this genotype.
Different kinds of electrophoretic methods based on protein patterns have been used for the identification and the characterization of cultivars (Yüzbaşıoğlu et al., 2008). Some researchers proposed that protein profiles may be useful as an indicator of taxonomic relationships within some species (Lioi et al., 2005), but others are of the opinion that this method is insufficient for discrimination at the cultivar level (Yüzbaşıoğlu et al., 2008). Our results revealed that SDS-PAGE of protein patterns gave specific BPs that distinguishing resistant faba bean genotypes from susceptible genotypes; however, the differentiations were not sufficient in distinguishing the genotypes. The results are in partial agreement with the findings of Alghamdi (2009a), Balkaya and Yanmaz (2002), and Duc (1997).
The results of the cluster analysis based on ISSR, SRAP, SSR, protein and combined ISSR + SRAP + SSR + protein data showed that the species obviously clustered into three major groups according to their resistance to FW. These results are in agreement with those obtained by Abdel-Razzak et al. (2012), who clarified that ISSR markers and protein analysis were helpful to recognize genetic variation among faba bean genotypes. The obtained results indicated management options for FW using resistant genotypes, especially genotypes Assiut-215, Roomy-3, Marut-2, and Giza-2. Resistance in these genotypes may be due to their genetic variations (Colditz et al., 2004; Mahmoud et al., 2015; Nyamsuren et al., 2003). These results also, in agreement with Kamboj et al. (1990) and Mahmoud et al. (2018) they mentioned that the inheritance of resistance to vascular wilt is controlled by independently segregating genes based on the reaction of individual plants.
In this study out of 132 polymorphic bands amplified by 17 markers; eight fragments generated by HB14 (498 bp), SRAP-2 (289 bp), SRAP-4 (319 bp), SRAP-5 (561 bp), SRAP-6 (146 bp, 593 bp, and 656 bp), and SSR-1 (498 bp) were present in all moderate and resistant genotypes, but absent in all highly susceptible genotypes. Furthermore, 10 DNA fragments, 457 bp (HB10), 209 bp (HB14), 312 bp (HB12), 219 bp (HB06), 984 bp (SRAP-1), 720 bp (SRAP-3), 401 bp (SRAP-4), 291 bp (SRAP-6), 419 bp (SSR-1), and 230 bp (SSR-3) were present in resistant faba bean genotypes, but absent in moderate and susceptible genotypes suggesting that these fragments could be considered as markers associated with the resistance to the FW in these genotypes. Moreover, these results indicate that the resistance genes present in faba bean genotypes may be different from each other. However, further genotyping is still required to validate the suitability of these markers for marker-assisted faba bean breeding. Similar results regarding genes responsible for the resistance to Fusarium vascular wilt were found by Eujayl et al. (1998), Rashid and Bernier (1993), and Rubio et al. (2003).
Our results in agreement with several reports in different crops, Parmar et al. (2009) reported the primer set SSR-67, a molecular marker linked to the FW disease resistance in tomato which could be employed for the discrimination of susceptible and resistant entity amongst the tomato population. Singh et al. (2013) analyzed genetic diversity among forty pigeonpea genotypes using 12 polymorphic SSR markers. They found that most of the FW resistant genotypes were grouped together in one cluster, and detected SSR markers associated with FW resistant. Mace et al. (2006) identified an association of eight SSR markers with rust and late leaf spot (LLS). While, Mondal and Badigannavar (2010) found three and four SSR markers associated with rust and LLS resistance in cultivated groundnut (Arachis hypogaea L.) germplasm, respectively. Singh et al. (2011) detected RAPD markers linked to FW resistance in castor bean (Ricinus communis L.).
Infection by pathogens can cause substantial biochemical alterations leading to harmful effects on plant health. Diseases may reduce plant growth and yield (Mofunanya and Edu, 2015). Because enzymes control biochemical reactions, and their synthesis is under the control of specific gene(s), any change in the activity and/or the number of BPs would reflect the pattern of gene expression and corresponding metabolic events in the cell (Neog et al., 2004).
Our study compared the SDS-PAGE protein patterns and isozyme bands of control and infected faba bean plants. The SDS-PAGE protein profile of total soluble protein from infected faba bean plants showed different band patterns from those of their respective control plants. Fusarium oxysporum infection caused disappearance of protein bands in some faba bean genotypes, while these bands were present in the control plants.
The absence of these protein bands in the infected plants can be explained as a result of the inhibition of the trigger for resistant gene transcription, which produces pathogen-related proteins, as previously reported by El-Khallal (2007).
On the other hand, some new protein bands were induced in infected plants, such that these bands act as a protein marker for resistance mechanisms allowing plants to become more resistant to the pathogens. Plants have flexible detection systems and most likely employ several recognition and signal transduction pathways to activate their defenses. Different types of proteins have been found to play specific roles in the plant defense mechanism and resistance to plant pathogens (Kumari and Vengadaramana, 2017).
Our results show that the, protein bands at molecular weights 30.1 kD and 10.2 kD were expressed only in all resistant genotypes, but not expressed in all moderate and susceptible genotypes. On the other hand, protein band 21.4 kD was expressed in all moderate and resistant genotypes, but absent in highly susceptible genotypes. For infected plants; protein band of molecular weight 26.2 kD induced in all faba bean genotypes infected with FW as compared with control plants (Fig. 4B). Also, new band of molecular weight 29.7 kD appeared only in all resistant genotypes infected with FW as compared with control plants. These findings suggested that these bands could be used as biochemical markers associated with the resistance to Fusarium vascular wilt in tested faba bean genotypes. These results are similar to the results obtained by El-Sayed (2006) who found that protein electrophoresis of two genotypes of faba bean infested with bruchid detected distinct bands which differentiated the most susceptible from the most resistant.
The induction of and changes in isozyme expression play important roles in the cellular defense against oxidative stress caused by pathogen infection (Almagro et al., 2009). Our results show that the isozyme patterns of peroxidase and chitinase of control and infected plants showed variations in BPs. New isozymes induced in infected plants were not observed in control plants. Furthermore, some isozyme bands disappeared in all or some infected faba bean genotypes. The induction of new isozymes in response to pathogen infection indicates that the induced isozymes have definitive roles in suppressing disease development in plant tissue. Several studies have investigated alterations in isozyme patterns in plants infected with different pathogens. Prasad et al. (2003) analyzed isozyme variability among Fusarium udum resistant cultivars of pigeon pea (Cajanus cajan L.). They concluded that the polymorphisms observed among isozymes such as peroxidases, phosphatases and esterases may be helpful in the characterization of resistance and susceptibility. Sahoo et al. (2007) indicated the presence of two new peroxidase bands in resistant genotypes of taro (Colocasia esculenta L.) under induced blight conditions. Saikia et al. (2005) analyzed chitinase isozymes in induced chickpea seedlings. They observed two chitinase isozyme bands were detected in diseased plants but not in healthy plants.
In conclusion, the present study reports genetic diversity of 16 elite faba bean genotypes, resistant/susceptible to FW. These genotypes were broadly grouped into two groups comprising most of the FW moderate and resistant genotypes into one group and susceptible genotypes in other groups. The present study revealed feasibility of prescreened ISSR, SRAP, SSR, protein, and isozyme markers to be used in genetic diversity analysis and their potential association with FW resistance. Diverse faba bean FW resistant/susceptible genotypes reported in this study will be useful to generate mapping populations for FW mapping and work in breeding programs for developed and release new resistant cultivars for Fusarium wilt.

Supplemental Materials

Notes

Electronic Supplementary Material

Supplementary materials are available at The Plant Pathology Journal website (http://www.ppjonline.org/).

References

Abdel-Razzak, HS, Alfrmawy, AM, Ibrahim, HM and El-Hanafy, AA 2012. Genetic diversity in faba bean (Vicia faba L.) using inter-simple sequence repeat (ISSR) markers and protein analysis. Life Sci J. 9:497-503.
Abdel Rehim, MA 1962. Studies on the organisms causing root rot and wilt of horse beans (Vicia faba var. equina) in U. A. R. PhD thesis. Faculty of Agriculture, Alexandria University, Qism Bab Sharqi, Egypt.
Abdul Wahid, OA, Ibrahim, ME and Omar, MA 1998. Occurrence of soil suppressiveness to Fusarium wilt disease of broad bean in Ismailia Governorate. J Phytopathol. 146:431-435.
crossref
Abou-Zeid, NM, El-Morsy, GA, Hassanein, AM and Arafa, MK 1997. Major organisms causing root-rot, wilt and their relative importance on faba bean, lentil and chickpea. Egypt J Agric Res. 75:529-542.
Ahmed, HFS, El-Araby, MM and Omar, SA 2002. Differential effect of jasmonic acid on the defense of faba bean against Fusarium wilt: modulation of other phytohormones and simple phenols. Int J Agric Biol. 4:447-453.
Alghamdi, SS 2009a. Chemical composition of faba bean (Vicia faba L.) genotypes under various water regimes. Pak J Nutr. 8:477-482.
crossref
Alghamdi, SS 2009b. Varietal identification and genetic purity assessment of F1 hybrid seeds using RAPD markers in faba bean (Vicia faba L.). Acta Hortic. 829:269-274.
Alghamdi, SS, Al-Faifi, SA, Migdadi, HM, Ammar, MH and Siddique, KHM 2011. Inter-simple sequence repeat (ISSR)-based diversity assessment among faba bean genotypes. Crop Pasture Sci. 62:755-760.
crossref
Alghamdi, SS, Al-Faifi, SA, Migdadi, HM, Khan, MA, El-Harty, EH and Ammar, MH 2012. Molecular diversity assessment using sequence related amplified polymorphism (SRAP) markers in Vicia faba L. Int J Mol Sci. 13:16457-16471.
crossref pmid pmc
Ali, MB, Haridy, AG and Mahmoud, AF 2019. Evaluation of faba bean genotypes for yield and resistance to Fusarium root rot under greenhouse and field conditions. Int J Biosci. 14:374-385.
Almagro, L, Gómez Ros, LV, Belchi-Navarro, S, Bru, R, Ros Barceló, A and Pedreño, MA 2009. Class III peroxidases in plant defence reactions. J Exp Bot. 60:377-390.
crossref pmid
Amar, MH and El Wahab, MA 2013. Comparative genetic study among Origanum L. plants grown in Egypt. J Biodivers Environ Sci. 3:208-222.
Ammar, MH, Alghamdi, SS, Migdadi, HM, Khan, MA, El-Harty, EH and Al-Faifi, SA 2015. Assessment of genetic diversity among faba bean genotypes using agro-morphological and molecular markers. Saudi J Biol Sci. 22:340-350.
crossref pmid pmc
Balkaya, A and Yanmaz, R 2002. Morphological properties of cultivar nominate selected Black Sea Region bean populations and identification by protein markers. Ankara Univ J Agric Sci. 9:182-188.
Bebeli, P and Kaltsikes, P 1993. New developments in varietal identification. In: Seed science and technology, eds. by AJG van Gastel, MA Pagnotta and E Porceddu, pp. 161-172. ICARDA, Aleppo, Syria.
Booth, C 1971. The genus Fusarium. Commonwealth Agricultural Bureaux (for the) Commonwealth Mycological Institute, Kew, UK. 237.
Budak, H, Shearman, RC, Parmaksiz, I, Gaussoin, RE, Riordan, TP and Dweikat, I 2004. Molecular characterization of buffalograss germplasm using sequence-related amplified polymorphism markers. Theor Appl Genet. 108:328-334.
crossref pmid
Colditz, F, Nyamsuren, O, Niehaus, K, Eubel, H, Braun, HP and Krajinski, F 2004. Proteomic approach: identification of Medicago truncatula proteins induced in roots after infection with the pathogenic oomycete Aphanomyces euteiches. Plant Mol Biol. 55:109-120.
crossref pmid
Costa, R, Pereira, G, Garrido, I, Tavares-de-Sousa, MM and Espinosa, F 2016. Comparison of RAPD, ISSR, and AFLP molecular markers to reveal and classify Orchard grass (Dactylis glomerata L.) germplasm variations. PLoS ONE. 11:e0152972
crossref pmid pmc
Dhole, VJ and Reddy, KS 2013. Development of a SCAR marker linked with a MYMV resistance gene in mungbean (Vigna radiata L. Wilczek). Plant Breed. 132:127-132.
crossref
Dice, LR 1945. Measures of the amount of ecologic association between species. Ecology. 26:297-302.
crossref
Dong, Y, Dong, K, Zheng, Y, Tang, L and Yang, ZX 2014 Faba bean Fusarium wilt (Fusarium oxysporum) control and its mechanism in different wheat varieties and faba bean intercropping system. J Appl Ecol. 25:1979-1987. (in Chinese)..
Duc, G 1997. Faba bean (Vicia faba L). Field Crops Res. 53:99-109.
crossref
El-Khallal, SM 2007. Induction and modulation of resistance in tomato plants against Fusarium wilt disease by bioagent fungi (arbuscular mycorrhiza) and/or hormonal elicitors (Jasmonic acid & Salicylic acid): 2-changes in the antioxidant enzymes, phenolic compounds and pathogen related- proteins. Aust J Basic Appl Sci. 1:717-732.
El-Sayed, ZG 2006. Genetic markers for some important traits in faba bean (Vicia faba L.). PhD thesis. Department of Genetics, Faculty of Agriculture, Ain Shams University, Cairo, Egypt.
Elwakil, MA, El-Refai, IM, Awadallah, OA, El-Metwally, MA and Mohammed, MS 2009. Seed-borne pathogens of faba bean in Egypt: detection and pathogenicity. Plant Pathol J. 8:90-97.
Eujayl, I, Erskine, W, Bayaa, B, Baum, M and Pehu, E 1998. Fusarium vascular wilt of lentil: inheritance and identification of DNA markers for resistance. Plant Breed. 117:497-499.
crossref
Gerlach, W and Nirenberg, HI 1982. The genus Fusarium: a pictorial atlas. Mitt Biol Bundesanst Land-Forstwirtsch Berlin-Dahlem. 209:1-406.
crossref
Ghafoor, A, Ahmad, Z, Qureshi, AS and Bashir, M 2002. Genetic relationship in Vigna mungo (L.) Hepper and V. radiata (L.) R. Wilczek based on morphological traits and SDS-PAGE. Euphytica. 123:367-378.
Ghislain, MD, Zhang, D, Fajardo, D, Huamán, Z and Hijmans, RJ 1999. Marker-assisted sampling of the cultivated Andean potato Solanum phureja collection using RAPD markers. Genet Resour Crop Evol. 46:547-555.
Giovannelli, JL, Farnham, MW and Wang, M 2002. Development of sequence characterized amplified region markers linked to downy mildew resistance in broccoli. J Am Soc Hortic Sci. 127:597-601.
crossref
Gomez, KA and Gomez, AA 1984. Statistical procedures for agricultural research. 2nd ed. John Wiley and Sons, New York, NY, USA. 704.
Guikema, JA and Sherman, LA 1981. Electrophoretic profiles of cyanobacterial membrane polypeptides showing heme-dependent peroxidase activity. Biochim Biophys Acta. 637:189-201.
crossref
Gupta, VK, Misra, AK, Gaur, RK, Jain, PK, Gaur, D and Sharma, S 2010. Current status of Fusarium wilt disease of guava (Psidium guajava L.) in India. Biotechnology. 9:176-195.
crossref
Hao, Q, Liu, Z-A, Shu, QY, Zhang, R, De Rick, J and Wang, L-S 2008. Studies on Paeonia cultivars and hybrids identification based on SRAP analysis. Hereditas. 145:38-47.
crossref pmid
Jhang, T, Kaur, M, Kalia, P and Sharma, TR 2010. Efficiency of different marker systems for molecular characterization of subtropical carrot germplasm. J Agric Sci. 148:171-181.
crossref
Johnson, R and Jellis, GJ 1992. Development in plant pathology. Vol. 1. Breeding for disease resistance. Kluwer Academic, Dordrecht, Netherlands. 206.
Khalaf, AEA, Afiah, SA, Khalil, RMA and Ahmed, AIS 2015. Molecular markers for resistance of chocolate spot disease in faba bean (Vicia faba L.) using ISSR-PCR. Middle-East J Sci Res. 23:558-567.
Kamboj, RK, Pandey, MP and Chaube, HS 1990. Inheritance of resistance to Fusarium wilt in Indian lentil germplasm (Lens culinaris Medik.). Euphytica. 50:113-117.
crossref
Kleczewski, NM and Egel, DS 2011. A diagnostic guide for Fusarium wilt of watermelon. Plant Health Prog. 12:27
crossref
Kumari, YSMAI and Vengadaramana, A 2017. Stimulation of defense enzymes in tomato (Solanum lycopersicum L.) and chilli (Capsicum annuum L.) in response to exogenous application of different chemical elicitors. Univers J Plant Sci. 5:10-15.
crossref
Laemmli, UK 1970. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 227:680-685.
crossref pmid
Li, G and Quiros, CF 2001. Sequence-related amplified polymorphism (SRAP), a new marker system based on a sample PCR reaction: its application to mapping and gene tagging in Brassica. Theor Appl Genet. 103:455-461.
crossref
Lioi, L, Piergiovanni, AR, Pignone, D, Puglisi, S, Santantonio, M and Sonnante, G 2005. Genetic diversity of some surviving on-farm Italian common bean (Phaseolus vulgaris L.) landraces. Plant Breed. 124:576-581.
crossref
Mace, ES, Phong, DT, Upadhyaya, HD, Chandra, S and Crouch, JH 2006. SSR analysis of cultivated groundnut (Arachis hypogaea L.) germplasm resistant to rust and late leaf spot diseases. Euphytica. 152:317-330.
crossref
Mahmoud, AF 2016a. Evaluation of certain antagonistic fungal species for biological control of faba bean wilt disease incited by Fusarium oxysporum. J Phytopathol Pest Manag. 3:1-14.
Mahmoud, AF 2016b. Genetic variation and biological control of Fusarium graminearum isolated from wheat in Assiut-Egypt. Plant Pathol J. 32:145-156.
crossref
Mahmoud, AF, Abou-Elwafa, SF and Shehzad, T 2018. Identification of charcoal rot resistance QTLs in sorghum using association and in silico analyses. J Appl Genet. 59:243-251.
crossref pmid
Mahmoud, AF, Hassan, MI and Amein, KA 2015. Resistance potential of bread wheat genotypes against yellow rust disease under Egyptian climate. Plant Pathol J. 31:402-413.
crossref pmid pmc
Mantel, N 1967. The detection of disease clustering and a generalized regression approach. Cancer Res. 27:209-220.
pmid
Mazen, MM, El-Batanony, NH, Abd El-Monium, MM and Massoud, ON 2008. Cultural filtrate of Rhizobium spp. and arbuscular mycorrhiza are potential biological control agents against root rot fungal diseases of faba bean. Global J Biotechnol Biochem. 3:32-41.
Mofunanya, A and Edu, EA 2015. Physiological and biochemical changes in Cucurbita moschata Duch Ex Poir inoculated with a Nigerian strain of moroccan watermelon mosaic virus (MWMV): Lagenaria breviflora isolate. Int J Plant Pathol. 6:36-47.
crossref
Mondal, S and Badigannavar, AM 2010. Molecular diversity and association of SSR markers to rust and late leaf spot resistance in cultivated groundnut (Arachis hypogaea L.). Plant Breed. 129:68-71.
crossref
Murray, MG and Thompson, WF 1980. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8:4321-4325.
crossref pmid pmc
Nagaraja, NR, Kumar, MBA, Prasad, CTM, Rajkumar, , Sinha, P and Jain, J 2016. Genetic diversity studies and screening for Fusarium wilt (Fusarium udum Butler) resistance in wild pigeonpea accessions, Cajanus scarabaeoides (L.) Thouars. Indian J Plant Genet Resour. 29:121-129.
crossref
Neog, B, Yadav, RNS and Singh, ID 2004. Peroxidase, polyphenol oxidase and acid phosphatase activities in the stigma-style tissue of Camellia sinensis (L) O. Kuntze following compatible and incompatible pollination. J Indian Inst Sci. 84:47-52.
Nyamsuren, O, Colditz, F, Rosendahl, S, Tamasloukht, MB, Bekel, T, Meyer, F, Kuester, H, Franken, P and Krajinski, F 2003. Transcriptional profiling of Medicago truncatula roots after infection with Aphanomyces euteiches (oomycota) identifies novel genes upregulated during this pathogenic interaction. Physiol Mol Plant Pathol. 63:17-26.
crossref
Okiror, MA 2002. Genetics of resistance to Fusarium udum in pigeonpea (Cajanus cajan (L.) Millsp.). Indian J Genet Plant Breed. 62:218-220.
Parmar, PP, Bhatt, KN, Oza, VP, Patel, AD, Kathiria, KB and Subramanian, RB 2009. Microsatellite marker associated with fusarium wilt resistance in tomato. World J Agric Sci. 5:389-393.
Powell, W, Machray, GC and Provan, J 1996. Polymorphism revealed by simple sequence repeats. Trends Plant Sci. 1:215-222.
crossref
Prasad, P, Eswara Reddy, NP, Anandam, RJ and Lakshmikantha Reddy, G 2003. Isozymes variability among Fusarium udum resistant cultivars of pigeonpea (Cajanus cajan (L.) (Millsp). Acta Physiol Plant. 25:221-228.
crossref
Prevost, A and Wilkinson, MJ 1999. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theor Appl Genet. 98:107-112.
crossref
Provan, J, Thomas, WTB, Forster, BP and Powell, W 1999. Copia-SSR: a simple marker technique which can be used on total genomic DNA. Genome. 42:363-366.
crossref
Rabey, JM, Vered, Y, Shabtai, H, Graff, E and Korczyn, AD 1992. Improvement of parkinsonian features correlate with high plasma levodopa values after broad bean (Vicia faba) consumption. J Neurol Neurosurg Psychiatry. 55:725-727.
crossref pmid pmc
Rashid, KY and Bernier, CC 1993. Genetic diversity among isolates of Rhizoctonia solani and sources of resistance in Vicia faba. Can J Plant Pathol. 15:23-28.
crossref
Reddy, MP, Sarla, N and Siddiq, EA 2002. Inter simple sequence repeat (ISSR) polymorphism and its application in plant breeding. Euphytica. 128:9-17.
Rohlf, FJ 2000. Statistical power comparisons among alternative morphometric methods. Am J Phys Anthropol. 111:463-478.
crossref pmid
Rubio, J, Hajj-Moussa, E, Kharrat, M, Moreno, MT, Millan, T and Gill, J 2003. Two genes and linked RAPD markers involved in resistance to Fusarium oxysporum f. sp. ciceris race 0 in chickpea. Plant Breed. 122:188-191.
crossref
Sahoo, MR, DasGupta, M, Kole, PC, Bhat, JS and Mukherjee, A 2007. Antioxidative enzymes and isozymes analysis of taro genotypes and their implications in Phytophthora blight disease resistance. Mycopathologia. 163:241-248.
crossref pmid
Saikia, R, Singh, BP, Kumar, R and Arora, DK 2005. Detection of pathogenesis-related proteins-chitinase and β-1,3-glucanase in induced chickpea. Curr Sci. 89:659-663.
Shen, SH, Chrétien, P, Bastien, L and Slilaty, SN 1991. Primary sequence of the glucanase gene from Oerskovia xanthineolytica. Expression and purification of the enzyme from Escherichia coli. J Biol Chem. 266:1058-1063.
crossref pmid
Singh, AK, Rai, VP, Chand, R, Singh, RP and Singh, MN 2013. Genetic diversity studies and identification of SSR markers associated with Fusarium wilt (Fusarium udum) resistance in cultivated pigeonpea (Cajanus cajan). J Genet. 92:273-280.
crossref pmid
Singh, M, Chaudhuri, I, Mandal, SK and Chaudhuri, RK 2011. Development of RAPD markers linked to Fusarium wilt resistance gene in castor bean (Ricinus communis L.). Genet Eng Biotechnol J. 28:1-9.
Snyder, WC and Hansen, HN 1940. The species concept in Fusarium. Am J Bot. 27:64-67.
crossref
Souframanien, J and Gopalakrishna, T 2006. ISSR and SCAR markers linked to the mungbean yellow mosaic virus (MYMV) resistance gene in blackgram (Vigna mungo (L.) Hepper). Plant Breed. 125:619-622.
crossref
Srivastava, RK, Mishra, SK, Singh, AK and Mohapatra, T 2012. Development of a coupling-phase SCAR marker linked to the powdery mildew resistance gene ‘er1’ in pea (Pisum sativum L.). Euphytica. 186:856-866.
crossref
Tanksley, SD, Young, ND, Paterson, AH and Bonierbale, MW 1989. RFLP mapping in plant breeding: new tools for old science. Nat Biotechnol. 7:257-264.
crossref
Tekalign, A, Derera, J, Sibiya, J and Mumm, RH 2019. Molecular analysis for genetic diversity and population structure of Ethiopian faba bean (Vicia faba L.) accessions. Plant Biol Crop Res. 1:1010.
Terzopoulos, PJ and Bebeli, PJ 2008. Genetic diversity analysis of Mediterranean faba bean (Vicia faba L.) with ISSR markers. Field Crops Res. 108:39-44.
crossref
Tufan, H and Erdoğan, C 2017. Genetic diversity in some faba bean (Vicia faba L.) genotypes assessed by simple sequence repeats. Biotechnol Biotechnol Equip. 31:29-35.
crossref
Vandemark, GJ, Ariss, JJ, Bauchan, GA, Larsen, RC and Hughes, TJ 2006. Estimating genetic relationships among historical sources of alfalfa germplasm and selected cultivars with sequence related amplified polymorphisms. Euphytica. 152:9-16.
crossref
Varshney, RK, Hoisington, DA and Tyagi, AK 2006. Advances in cereal genomics and applications in crop breeding. Trends Biotechnol. 24:490-499.
crossref pmid
Velasco-Ramírez, AP, Torres-Morán, MI, Molina-Moret, S, Sánchez-González, JDJ and Santacruz-Ruvalcaba, F 2014. Efficiency of RAPD, ISSR, AFLP and ISTR markers for the detection of polymorphisms and genetic relationships in camote de cerro (Dioscorea spp.). Electron J Biotechnol. 17:65-71.
crossref
Wang, Z, Wang, J-E, Wang, X-M, Gao, H-W, Dzyubenko, NI and Chapurin, VF 2012. Assessment of genetic diversity in Galega officinalis L using ISSR and SRAP markers. Genet Resour Crop Evol. 59:865-873.
crossref
Yahia, Y, Hannachi, H and Ferchichi, A 2014. Genetic diversity of Vicia faba L. based on random amplified polymorphic DNA and simple sequence repeat markers. Acta Bot Gallica Bot Lett. 161:151-158.
crossref
Yang, T, Bao, S-Y, Ford, R, Jia, T-J, Guan, J-P, He, Y-H, Sun, X-L, Jiang, J-Y, Hao, J-J, Zhang, X-Y and Zong, X-X 2012. High-throughput novel microsatellite marker of faba bean via next generation sequencing. BMC Genomics. 13:602
crossref pmid pmc
Yüzbaşioğlu, E, Açik, L and Özcan, S 2008. Seed protein diversity among lentil cultivars. Biol Plant. 52:126-128.
crossref
Zeid, M, Schön, C-C and Link, W 2003. Genetic diversity in recent elite faba bean lines using AFLP markers. Theor Appl Genet. 107:1304-1314.
crossref pmid
Zietkiewicz, E, Rafalski, A and Labuda, D 1994. Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics. 20:176-183.
crossref pmid
Fig. 1
Gel electrophoresis inter-simple sequence repeat (ISSR) primers in 16 faba bean genotypes and dendrogram developed from six ISSR primers. 1, Misr-3; 2, Nobaria-1; 3, Nobaria-2; 4, Assiut-215; 5, Romy-3; 6, Misr-1; 7, Assiut-143; 8, Nobaria-3; 9, Mariut-2; 10, Giza-843; 11, Giza-40; 12, Giza-2; 13, Sakha-1; 14, Wadi-5; 15, Sakha-4; 16, Romy-80.
ppj-36-011f1.jpg
Fig. 2
Gel electrophoresis sequence related amplified polymorphism (SRAP) primer combinations in 16 faba bean genotypes and dendrogram developed from six SRAP primers. 1, Misr-3; 2, Nobaria-1; 3, Nobaria-2; 4, Assiut-215; 5, Romy-3; 6, Misr-1; 7, Assiut-143; 8, Nobaria-3; 9, Mariut-2; 10, Giza-843; 11, Giza-40; 12, Giza-2; 13, Sakha-1; 14, Wadi-5; 15, Sakha-4; 16, Romy-80.
ppj-36-011f2.jpg
Fig. 3
Gel electrophoresis simple sequence repeat (SSR) primer combinations in 16 faba bean genotypes and dendrogram developed from five SSR primers. 1, Misr-3; 2, Nobaria-1; 3, Nobaria-2; 4, Assiut-215; 5, Romy-3; 6, Misr-1; 7, Assiut-143; 8, Nobaria-3; 9, Mariut-2; 10, Giza-843; 11, Giza-40; 12, Giza-2; 13, Sakha-1; 14, Wadi-5; 15, Sakha-4; 16, Romy-80.
ppj-36-011f3.jpg
Fig. 4
Dendrogram of 16 faba bean genotypes developed from combined inter-simple sequence repeat, sequence related amplified polymorphism, and simple sequence repeat markers data.
ppj-36-011f4.jpg
Fig. 5
(A-C) Protein banding pattern of 16 faba bean genotypes. Control (A), infected with Fusarium oxysporum (B), and dendrogram (C) of 16 faba bean genotypes developed from protein profiles under control condition. 1, Misr-3; 2, Nobaria-1; 3, Nobaria-2; 4, Assiut-215; 5, Romy-3; 6, Misr-1; 7, Assiut-143; 8, Nobaria-3; 9, Mariut-2; 10, Giza-843; 11, Giza-40; 12, Giza-2; 13, Sakha-1; 14, Wadi-5; 15, Sakha-4; 16, Romy-80.
ppj-36-011f5.jpg
Fig. 6
(A, B) Chitinase (Chit) banding pattern of 16 faba bean genotypes. Control (A) and infected with Fusarium oxysporum (B). (C, D) Peroxidase (PRX) banding pattern of 16 faba bean genotypes. Control (C) and infected with F. oxysporum (D), and dendrogram (E) of 16 faba bean genotypes developed from combined molecular and biochemical markers data. 1, Misr-3; 2, Nobaria-1; 3, Nobaria-2; 4, Assiut-215; 5, Romy-3; 6, Misr-1; 7, Assiut-143; 8, Nobaria-3; 9, Mariut-2; 10, Giza-843; 11, Giza-40; 12, Giza-2; 13, Sakha-1; 14, Wadi-5; 15, Sakha-4; 16, Romy-80.
ppj-36-011f6.jpg
Table 1
Reactions of faba bean genotypes for resistance to Fusarium wilt disease under greenhouse conditions
Cultivar Fusarium wilt severity (%)a Disease reactions

2016-2017 2017-2018 Average
Assiut-143 62.55 C 60.87 C 61.71 HS
Assiut-215 33.00 P 32.67 O 32.83 R
Giza-2 39.75 M 37.40 L 38.57 R
Giza-40 60.75 E 60.12 D 60.43 HS
Giza-843 48.85 G 46.00 H 47.42 MR
Marut-2 37.75 N 35.95 N 36.85 R
Misr-1 65.75 B 61.55 B 63.65 HS
Misr-3 66.25 A 62.10 A 64.17 HS
Nobaria-1 46.00 J 45.25 I 45.62 MR
Nobaria-2 47.25 H 44.00 K 45.62 MR
Nobaria-3 44.50 L 43.25 L 43.87 MR
Roomy-3 35.00 O 36.87 M 35.93 R
Roomy-80 62.37 D 59.40 E 60.88 HS
Sakha-1 45.55 K 44.50 J 45.02 MR
Sakha-4 46.50 I 47.65 F 47.07 MR
Wadi-5 49.75 F 46.70 G 48.22 MR

R, resistant, wilt severity >30 to 40%; MR, moderate resistant, wilt severity >40 to 50%; HS, highly susceptible, wilt severity ≥60%.

a Means within the same column followed by different letters are significantly different (P ≤ 0.05) based on least significant difference.

Table 2
Genetic diversity parameters obtained using 6 ISSR and 6 SRAP markers
Marker TB NPB PPB PIC MI RP
ISSR markers
 HB10 9 8 88.89 0.30 2.42 3.88
 HB14 9 3 33.33 0.10 0.31 1.25
 HB09 12 9 75.00 0.21 1.88 3.63
 HB13 7 4 57.14 0.14 0.54 1.50
 HB12 9 3 33.33 0.15 0.45 2.13
 HB06 7 6 85.71 0.24 1.43 2.13
  TOT 53 33 - - - -
  AVR 8.83 5.5 62.24 0.19 1.17 2.42
SRAP markers
 SRAP-1 15 9 60.00 0.22 1.97 4.63
 SRAP-2 13 12 92.31 0.21 2.52 3.38
 SRAP-3 14 13 92.86 0.28 3.66 5.25
 SRAP-4 13 13 100.00 0.28 3.64 4.75
 SRAP-5 6 5 83.33 0.26 1.29 2.00
 SRAP-6 15 10 66.67 0.25 2.47 5.50
  TOT 76 62 - - - -
  AVR 12.67 10.33 82.53 0.25 2.59 4.25
SSR markers
 SSR-1 10 7 70.00 0.20 1.38 2.50
 SSR-2 12 11 91.67 0.25 2.74 4.13
 SSR-3 8 7 87.50 0.22 1.55 2.13
 SSR-4 7 5 71.43 0.18 0.92 1.63
 SSR-5 11 7 63.64 0.16 1.14 2.38
  TOT 48 37 - - - -
  AVR 9.6 7.4 76.85 0.20 1.55 2.55
  TOT 177 132 - - - -
  AVR 10.41 7.76 73.69 0.21 1.78 3.10

ISSR, inter-simple sequence repeat; SRAP, sequence related amplified polymorphism; TB, total number of amplified bands; NPB, number of polymorphic bands; PPB, percentage of polymorphic bands; PIC, polymorphism information content; MI, marker index; RP, resolving power; TOT, total; AVR, average; SSR, simple sequence repeat.

Table 3
Faba bean genotypes identified by unique presence or absence of specific markers
Genotype Unique bands (bp)

Positive Negative
Misr-3 - 527 (HB10), 335 (SSR-1)
Nobaria-1 - 404 (SRAP-3), 742 (HB13)
Nobaria-2 954 (HB13) -
Assiut-215 333 (HB14) 358 (SRAP-2), 519 (SRAP-4), 248 (SSR-2)
Misr-1 - 308 (HB10), 186 (HB06)
Assiut-143 748 (SRAP-1) -
Nobaria-3 631 (SRAP-2), 372 (SSR-2) -
Mariut-2 562 (SRAP-4) -
Giza-843 230 (HB09), 152 (SSR-5) -
Giza-2 493 (HB09) -
Sakha-1 189 (SRAP-3) 313 (SSR-2)
Wadi-5 419 (SSR-4) -
Sakha-4 - 407, 310 (SRAP-2), 577 (HB09), 237 (SSR-5)


ABOUT
BROWSE ARTICLES
EDITORIAL POLICY
FOR CONTRIBUTORS
Editorial Office
Rm,904 (New Bldg.) The Korean Science & Technology Center 22,
Teheran-ro 7-Gil, Gangnamgu, Seoul 06130, Korea
Tel: +82-2-557-9360    Fax: +82-2-557-9361    E-mail: paper@kspp.org                

Copyright © 2022 by Korean Society of Plant Pathology.

Developed in M2PI

Close layer
prev next