The Potential Genetic Effect for Yield and Foliar Disease Resistance in Faba Bean (Vicia faba L.) Assessed via Morphological and SCoT Markers

Faba bean is considered one of the most prominent grain legumes, with high protein content for human food consumption and livestock feed. The present study evaluated the nature of gene action and determined the genetic diversity among different populations of three crosses for resistance to foliar diseases at the molecular level. Analysis of variance exposed significant differences among the generations for all measured traits. Both dominance and additive gene effects were essential, but dominance genes, for the most part, exhibited greater effects than additive ones. This indicates an essential role for dominant genes alongside the additives one in inheriting such traits. The third cross (Marina × Giza 40) gave desired significant and positive (additive × additive) values for the number of pods/plant, seeds/plant, and seed yield/plant, in addition to desirable negative values for chocolate spot and rust characteristics. Furthermore, assessing the lines under study using seven SCoT primers disclosed three bands with recorded molecular weights of 260, 207, and 178 bp, generated by SCoT-1, SCoT-4, and SCoT-7 primers, respectively. These bands exist in the resistant parent (Marina), which could be attributed to the high-disease-resistance phenotypes, and they are absent in the sensitive parent (Giza 40) and other putative sensitive lines. Based on the molecular profiles and the genetic similarity between parents and the selected lines, the highest similarity value (0.91) was detected between Marina genotype and BC1, revealing a high foliar disease resistance. Meanwhile, Giza 40 (susceptible to foliar diseases) exhibited the maximum value (0.93) with F2. Additionally, cluster analysis based on genetic relationships was performed, and a high level of correlation between the results of PCR-based SCoT analysis and the foliar disease reactions was observed in the field. Consequently, this study concluded that SCoT markers created reliable banding profiles for evaluating genetic polymorphism among faba bean lines, which could be a foundation for developing an efficient breeding program.


Introduction
Faba bean (Vicia faba L.) is one of the most prominent cool-season grain legumes cultivated worldwide as a vital source of protein for human food and animal feed uses.Faba bean plays a role in contributing to nitrogen fixation and soil biodiversity in cereals cropping systems [1,2].It is a diploid species with 2n = 12 that belongs to the Fabaceae family and can be divided into three subspecies based on seed size including large-seeded type (V.faba subsp major), medium-seeded type (V.faba subsp equina), and small-seeded type (V.faba subsp minor).The cross-pollination in faba beans has been observed to range from 8% to 84%, with an average of around 30% to 60% depending on the climatic factors and genotypes [3][4][5].
The Faba bean are classified as the fourth crop among legume crops after pea, chickpea, and lentil and are considered the seventh most produced legume [6,7], with fully mature (dry seeds) production reaching 4.84 million tons in 2017 [8].In Egypt, its productivity recorded 170,236 metric tons of total production from the total harvested area of 120,116 faddan (1 faddan = 0.42 ha) with a seed yield average of 1.42 tons/faddan (Economic Affairs Sector, Ministry of Agriculture 2021).However, the yield is still low due to the partial adaptability of available cultivars to wide-ranging environments and disease susceptibility [9,10].
Foliar and soil-borne diseases are major limiting factors in legume production since chocolate spot and rust are the major diseases affecting Faba bean productivity and quality.Two different pathogens cause chocolate spot diseases.Botrytis fabae sard is the most severe and the more specific disease in Egypt, especially in humid areas.Moreover, this disease has been reported to cause a dense reduction in yields, reaching over 60% of yield by destroying the plant leaves and limiting photosynthetic activity [11][12][13][14].The second disease is the Faba rust caused by Uromyces viciae-fabae (Pers.), a widespread disease with the highest identified virulence rate in the Egyptian populations.The virulence tendency varies from moderate to severe, with yield losses ranging 22-26% from the yield [15][16][17].
In faba bean germplasm, only a few sources of genetic resistance to chocolate spots have been known [18,19].Different tactics may apply for controlling chocolate spot, and more effort is needed in order to find resistant materials adapted to different environments [10,20].Furthermore, improving resistance genotypes for foliar diseases is the most vital objective for a breeding program.Breeders use several biometrical methods to formulate the most efficient breeding systems for estimating the genetic effects of genes controlling quantitative traits to improve the breeding program.Thus, generation means analysis, gene action, inbreeding depression, genetic parameters, and anticipated genetic gain from character selection provide the essential evidence for plant breeders to predict the effective breeding program.These parameters can be used in the first generations to obtain maximum improvement and to enhance yield potentials amongst a large selection [21,22].
The six-parameter model analysis includes parents (P 1 and P 2 ), first and second generations (F 1 and F 2 ), and the first two backcrosses (BC 1 and BC 2 ) is considered the best analysis method.This method was suggested by Gamble [23] and supplies breeders with the early genetic information of selected genotypes.
Heterosis and inbreeding depression are essential for explaining the genetic parameters in faba bean.Heterosis for seed yield is estimated to indicate heterotic effects in yield components.Consequently, heterosis resulting from the combined action and interaction between allelic and inter-allelic genes is an effective means by which faba bean and hybrid combinations will improve yield [22,24].The heterotic effects in faba bean ranged from significantly positive to significantly negative estimates, contributing to the improvement in yield and its components [25,26].Inbreeding depression reduces auto fertility and yield in the absence of pollinator, and yield can be minimized by 11% through the loss of heterosis, whereas high inbreeding depression recorded lower values in the F 2 s, in general, than the F 1 s for most traits in faba bean [26][27][28].
Estimating heritability is essential for predicting the expected genetic advance through selection in segregating populations.Genetic advances help us understand gene action's nature, which expresses many polygenic characteristics.High values of genetic advance indicate additive gene action, whereas low values point to non-additive gene action; subsequently, the heritability estimation will be dependable if escorted by a high genetic advance [29].
The success of any breeding program depends on the presence of sufficient genetic variability among used genotypes to permit effective selection.Mean analysis was developed to estimate the variance of genetic components and the acceptable knowledge obtained about the mode of inheritance and the nature of gene effects.The gene effects are essential for developing high-yielding varieties in faba bean [12].The yield potential remains unrealized in faba bean mainly due to the lack of success in hybrid breeding for the exploitation of heterosis.Therefore, improvements in seed yield and yield stability addressed through the component traits, namely, resistance against different pathogens and pests, are considered major breeding objectives in faba bean enhancement [25].
Assessing genetic variation is crucial in plant breeding as it enables the selection of genotypes that exhibit favorable traits [30].Various molecular markers have been used to demonstrate plant and other organisms' genetic variation [31,32].Molecular markers are fast, reliable, repeatable, unaffected by environmental conditions, and used to demonstrate genetic variation and select important plant agricultural characteristics [33].Many molecular techniques, such as random amplification of polymorphic DNAs (RAPDs), restriction fragment length polymorphisms (RFLPs), target region amplification polymorphisms (TRAPs), sequence-specific amplification polymorphisms (SSAPs), and amplified fragment length polymorphisms (AFLPs), were implied to detect the genetic variability of Vicia species and V. faba L. populations [30].The start codon targeted (SCoT) marker system is a relatively new molecular tool used in plant research that targets conserved regions flanking the ATG start codon in plant genes [34].This technique has many advantages, including the fact that it does not require prior sequence information, allowing it to be used with a wide range of plant species; it is highly reproducible and cost-effective; and it is linked to functional genes or regions surrounding these genes and their corresponding traits [35][36][37].Remarkably, based on the SCoT marker simplicity (it only requires a single primer for amplification), it produces a high percentage of polymorphism and an abundance of available genetic information [38].
The present study aimed at (i) determining the behavior of gene action and heritability of yield and its components in three crosses of faba bean using six populations; (ii) exploring the inbred depression in yield and its components in the same crosses and populations; (iii) assessing the genetic variation among selected lines at the molecular level; (iv) developing some molecular genetic markers (SCoT) associated with resistance and susceptibility to foliar diseases.The previous measurements are essential for improving foliar disease tolerance in faba bean varieties.

The Performance of Crosses
The mean performance and variance of all the studied traits in the three crosses (Nubaria 5 × Sakha 1, Nubaria 5 × Giza 40 and Marina × Giza 40) for the six populations (P 1 , P 2 , F 1 , F 2 , BC 1 , and BC 2 ) were presented in Table 1 Results revealed a difference between every two parents in all studied traits of the three crosses.The parental genotypes could be classified into two groups, the first group included the genotypes of Nubaria 5 (P 1 ), Sakha 1 (P 2 ), and Marina (P 4 ), which are considered resistant to foliar diseases.On the contrary, the second group involved Giza 40 (P 3 ), showed high susceptibility to foliar diseases, accompanied by affected yield components and seed yield/plant.The heterosis, inbreeding depression, and potency ratio of the three crosses for the studied traits are mentioned in Table 2.In most cases, all traits revealed significant and/or highly significant and positive heterotic effects regarding heterosis over better and midparent (BP and MP).Results showed positive and significant (p ≤ 0.05 and p ≤ 0.01) heterosis over mid-parent in all the measured traits through the three crosses except a chocolate spot in the second and third crosses and a rust trait in the first cross, which showed a significant and negative effect.Meanwhile, heterosis over better parent recorded the same trend except for some traits such as plant height in the third cross, chocolate spot in the first and second crosses, and rust in the first cross.In this light of that, both types of heterosis disclosed the desired best values of positive heterotic effects regarding the number of seeds/plant, number of branches/plant, seed yield/plant, and pods/plant in all crosses.However, the desired negative values of heterotic effects were detected with chocolate spot and rust traits in the second and first crosses, respectively.Applying the inbreeding depression revealed positive values for all traits except for a chocolate spot in the first and second crosses and rust diseases in the first and third crosses.Inbreeding depression significantly or highly significantly (p ≤ 0.05 and p ≤ 0.01) and positively increased and produced the maximum number of branches/plant in the first and third crosses, the number of pods/plant in the first cross, and chocolate spot and rust in the third and second crosses.Regarding the potence ratio, values greater than unity (p > +1) reflect over-dominance towards the better parent for most cases.Over-dominance manifested for all evaluated traits under different crosses except for the chocolate spot in all crosses, plant height, number of pods/plant, and rust trait in the third cross with scored values less than unity, demonstrating a partial dominance effect.
Applying scaling tests, A, B, C, and D (Table S1) revealed significant and positive or negative differences for all the assessed traits through the three crosses, indicating the presence of non-allelic interactions and the inadequacy of the simple model in interpreting the differences between population means.Hence, generation means as genetic analysis produces genetic parameters such as (m) mean effect, (a) additive, (d) dominance, (aa) additive × additive, (ad) additive × dominance, and (dd) dominance × dominance, as previously conducted by Gamble [23].These assessed parameters in the three crosses are mentioned in Table 3.The mean effect (m) showed highly significant (p ≤ 0.01) differences in the three crosses for all the evaluated characteristics, which specified the inheritance of quantitative effects.The effects of the additive gene (a) presented a highly significant effect (p ≤ 0.01) in the positive or negative direction for plant height in all three crosses, plus number of branches/plant, number of pods/plant, and rust through the first cross, in addition to the rust trait in the second cross.The estimates of the dominance gene effect (d) manifested a positive and highly significant effect (p ≤ 0.01) for plant height and number of branches/plant in the first cross, whereas the first and second crosses exhibited the same behavior with respect to number of pods/plant, number of seeds/plant, and seed yield/plant.Meanwhile, dominance with a negatively significant effect (p ≤ 0.01) was expressed for chocolate spot and rust traits (negative values are eligible) in the second and third crosses.
When going forward and estimating the values of additive × additive (aa), a positively significant effect for plant height, number of branches/plant, and seed yield/plant in the first cross was detected, whereas traits such as number of pods/plant and number of seeds/plant recorded this positive effect in both first and third crosses.Meanwhile, a values were negative for plant height, number of pods/plant, and number of seeds/plant in the second cross, as well as for foliar diseases in the second and third crosses.Doing the same analysis for additive × dominance (ad) demonstrated a highly significant and positive behavior through all the assessed traits in the third cross except two traits (number of pods/plant and rust), whereas the second cross disclosed the same trend with positive effect in traits such as plant height and foliar diseases.When applying dominance × domi-nance (dd) assessment, the second cross exhibited a desirable and highly significant effect regarding traits such as plant height, number of pods/plant, and number of seeds/plant.
The three crosses assessed for all measured traits in relation to heritability in broad (h 2 b ) and narrow (h 2 n ) senses, in addition to genetic advance (GS), are presented in Table 4.The heritability, in a broad sense, recorded the best values in the third cross for number of seeds/plant (90.1), whereas the first cross recorded the uppermost value (92.2, 91.63) for foliar diseases (chocolate spot and rust), respectively.Meanwhile, the lowest value was recorded (57.28) for seed yield/plant in the second cross.The estimated value of heritability in a narrow sense using the data of F 2 and backcross were evaluated, and they presented from high to moderate values for most traits.Traits including number of pods/plant and number of seeds/plant in the second cross, in addition to chocolate spot and rust in the first cross increasing to the maximum values.On the other hand, the lowest values were assigned for number of branches in the third cross and seed yield/plant in the second cross.Regarding the genetic advance under selection (GS%), the number of pods/plant and the number of seeds/plant in the first cross had the best percentages (80.3 and 77.55), respectively.Meanwhile, the minimum values were detected in plant height trait in all crosses.

Genetic Parameters
Genetic parameters are critical for the effectiveness of a polymorphism-based marker technique used in discriminatory genotypes.Table 6 shows that several genetic parameters were estimated using the online Marker Efficiency Calculator (iMEC) to evaluate the informative and discriminatory power of subsequent faba bean generations.Heterozygosity represents the direct count of heterozygosity in the population and is estimated based on the allele frequency of individuals in the population according to the Hardy-Weinberg equilibrium.The heterozygosity index (H) was obtained and ranged from 0.416 (primer SCoT-1) to 0.491 (primer SCoT-6), with an average of 0.462.Moreover, the polymorphism index content (PIC) analysis was performed to determine the efficiency of each SCoT primer in expressing polymorphic loci in lines under evaluation.The calculated (PIC) values for primers ranged from 0.330 (SCoT-1 and SCoT-7) to 0.371 (SCoT-6), with an average of 0.355.Going forward, the effective multiplex ratio (E) was assessed and pointed to 1.424 with primer SCoT-3 to 3.515 with SCoT-7 primer, with an average of 2.260.The arithmetic means heterozygosity (H.av) ranged from 0.0024 for SCoT-5 primer to 0.0050 for SCoT-6 primer, with an average of 0.0034.The highest value of the marker index (MI) was detected with SCoT-1 and SCoT-7 primer (0.0089), while the minimum value was 0.0049 and revealed by SCoT-3 primer with an average of 0.0074.The discriminating power (D) of the applied primers ranged from 0.505 (SCoT-1) to 0.875 (SCoT-3), with an average of 0.701.Resolving power (RP) is used to describe the capacity of the marker combination to detect the differences among various lines.RP values of the primers varied from 0.485 (SCoT-1) to 3.152 (SCoT-5), with an average of 1.818.

Genetic Distance and Similarity
Table S3 summarizes the estimated genetic distance between lines according to the Dice similarity coefficient.The genetic distance reflects the genetic relationships and the direction of the genetic improvement process.The data showed that the genetic similarity between pairs of lines ranged from 0.42 to 0.96.The biggest genetic similarity value was observed between BC 1 (1) and both BC 1 (2) and BC 1 (3); these lines shared a high resistance to foliar diseases, whereas the lowest value of genetic similarity was detected among F 2 (2) and F 2 (6), with a contradictory ability towards resistance to foliar diseases.Concerning Marina's parents, genetic similarity values ranged from 0.63 to 0.91, since the uppermost genetic similarity value was found between Marina and BC 1 (5) with high foliar disease resistance.
In comparison, the lowest genetic similarity value manifested among Marina and BC 2 (7).The previously mentioned lines exhibited a contradictory ability toward foliar disease resistance.Concerning Giza 40 parent, genetic similarity values ranged from 0.60 to 0.93.Giza 40 and F 2 (10) manifested the highest genetic similarity value while exhibiting susceptibility to foliar diseases.Contrarily, the minimum value in genetic similarity was among Giza 40 and F 2 (6), with a contradictory ability towards foliar disease resistance.

Discussion
Faba bean (Vicia faba L.) is considered one of the most popular and important l

Discussion
Faba bean (Vicia faba L.) is considered one of the most popular and important legume crops in the Mediterranean and Central Asia (its native area of origin).It can be used as a protein source in human food consumption in addition to livestock feed based on its high protein content in seed and straw.Biotic stresses such as foliar diseases and Orobanche affect faba bean production and quality since early disease development during flowering or the season could cause large losses and severe damage [39,40].
Concerning chocolate spot and rust diseases, BC 1 and BC 2 revealed the lowest numbers of infected plants, especially in the third cross, and this may refer to to the first parent (Marina P 1, ) in these crosses, which was not adapted for the other environmental factors.Likewise, this genotype (Marina) is resistant to foliar diseases, mainly to chocolate spot, especially in F 1 , F 2 , BC 1, and BC 2 plants, in addition to rust in BC 1 and BC 2 when compared with the other two crosses.
It was observed that F 1 mean values were higher than other populations and mostly exceeded the high parent with respect to most evaluated traits in the three crosses, which may refer to over-dominance.Then, selection could effectively improve these traits in the next generation, particularly the number of seeds/plant and seed yield/plant.Results obtained in all measured traits established that the variance of the segregated generations (F 2 and BCs) was greater than the non-segregated generations (P 1 , P 2 , and F 1 ), pointing to the environmental effects on the gene expression and effectiveness of selection for these traits.Therefore, parents are precisely selected to find the desired recombination in the segregated generations, and these results are in consonance with the findings of prior published reports by Ashrei et al., Hendawy, Koumber and El-Gammaal, Akhshi et al. [22,[41][42][43].
The heterosis percentage is a deviation of mean performance for the F 1 average from the mid-parent or better parent with respect to the traits.Both types of heterosis exhibited the best values and positive heterotic effects in the number of seeds/plant, number of branches/plant, seed yield/plant, and pods/plant in all crosses, whereas the desired negative values of heterotic effects were observed in the second and first crosses with chocolate spot and rust traits, respectively.These results indicated that the dominance direction was toward the best respective parent, and the heterotic effect might be due to the dominance and/or dominance × dominance effects.These findings are in harmony with those obtained by Ashrei et al. [22].Hybrid vigor means that trait analysis recorded statistically positive and significant different results.This means that the F 1 superior to the mid-parent might reflect the genetic variability in crossed parents and vice versa [26,44].The value of inbreeding depression is the opposite of the heterosis value.In the case of crossing between inbreed lines, their progeny may achieve higher performance than them (hybrid vigor), as reported by Falconer [45].The reduction in the mean performance of F 2 indicates an increment in inbreeding depression and vice versa, which implies the low genetic diversity of the parents [22,26,46].In the potence ratio, values of measured characteristics, more than unity, reflect over-dominance effects.Most of the evaluated traits had an over-dominance effect, suggesting that the selection must be delayed to late generations [42].
In the current study, the dominance degree appeared higher than the additive gene effect for most traits in the three crosses, indicating the significance of additive and dominance gene effects playing a vital role in the inheritance of these traits.Moreover, we can select desirable characters in the early generations, but the last generations are influential.
The gene effect of (aa) was relatively more important than the additive effects, whereas the dominance effect was less important.Furthermore, the effects of (dd) are higher in magnitude than the epistatic gene effects (aa) for all traits under study in the three crosses except the number of branches and number of pods/plant in the first cross, and number of seeds/plant and seed yield/plant in the third cross.The previous findings suggest that selection plays a vital role if delayed till dominance, and its epistatic effects are reduced to the lowest value, slowing down the progress of selection.These results agreed with the research reported by Akhshi et al., Attia and Salem, and Salama and Mohamed [43,[47][48][49].
Heritability estimation of different traits is essential for crop improvement programs and for predicting the response of selection depending on the magnitudes of genetic variance components for additive and dominance.High heritability, in a broad sense, indicates the essential role for the inheritance of additive-and non-additive-effect genes related to these traits.Moreover, high heritability in a narrow sense, indicates that the genetic variance of additives played a vital role in the existence of variability.These results are similar to those reported previously [22,26,44,46].Johnson et al. [50] reported that to predict selection results for the best individuals, we need to estimate the heritability value with genetic advances, which are more valuable than the value of heritability alone.Otherwise, Dixit et al. [51] disclosed that a high value of heritability is not always associated with a high genetic advance value.In this manuscript, the high value of genetic advance was associated with high and moderate heritability, in a narrow sense, for number of pods/plant in the first and third crosses, number of seeds/plant in the first cross, and seed yield/plant in the first and third crosses.Notably, moderate and low heritability plus predicted genetic advance indicated that the selection of faba bean in subsequent generations would be relatively more effective than in the early generation of F 2 .Similar findings were depicted by Abd El-zaher, Ashrei et al., Obiadalla-Ali et al., Abou-Zaid, and Aziz and Osman [21,22,44,52,53].
The obtained results revealed a high level of polymorphism, reflecting high divergence in the studied faba bean lines; subsequently, the SCoT technique is efficient in differentiating the lines under study and can be used to evaluate genetic variation among closely related cultivars.Our findings are in consonance with those reported by Nosair [54], who produced one hundred and eighty-three bands (183), generating 93.99% polymorphism of various Leguminosae species via SCoT markers.Furthermore, Essa et al. [55] observed a high level of polymorphism (reaching 70.93%) using six SCoT-PCR primers screened across eight faba bean cultivars.The present results indicate that the resistance genes present in faba bean genotypes may differ.However, more genotyping is still required in order to validate the suitability of these markers for marker-assisted faba bean breeding, and these loci could be effectively used in breeding programs.Heiba et al. [56] demonstrated that molecular markers exhibited interest-specific loci relating to faba bean genotypes' resistance to chocolate spot, whereas Bosily et al. [57] detected specific markers for barley leaf rust disease using five Scot primers.
The polymorphism index content (PIC) is performed to determine the efficiency of each SCoT primer in expressing polymorphic loci into the selected faba bean lines [58].The PIC values are more than 0.5 for highly informative markers, between 0.25-0.5 for reasonably informative markers, and less than 0.25 for slightly informative markers [59].The discrimination power of a marker represents its ability to distinguish between closely related genotypes; additionally, it can aid in the efficient evaluation of various markers [60].These findings indicated good sources of diversity, which will aid breeders in assessing genetic diversity and the relationship between different genotypes.The results confirmed that SCoT markers created trustworthy banding profiles, which were able to evaluate the genetic polymorphism among selected lines.The obtained findings coincide with those of Albrifcany and Askander [61], who confirmed that one could rely on the SCoT technique to evaluate the genetic diversity among the faba bean cultivars.The assessment of genetic similarity and genetic distance among plant cultivars helps adjust breeding programs to facilitate the selection of desirable genotypes among segregating and/or backcrossing populations [62].The results suggested that the SCoT markers showed considerable potential for identifying and discriminating faba bean lines via their resistance to foliar diseases.Their use for evaluating chocolate spot disease resistance in faba bean genotypes has been depicted by Heiba et al. [56].The current results showed that the SCoT marker technique is reliable in dividing sorts based on genetic distance [63].Finally, differences in lines clustering patterns may be due to marker sampling error, polymorphism level, or the number of loci and their coverage across the genome [64].

Field Experiment
This investigation was conducted at the Food Legumes Department, Sakha Agricultural Research Station, Agriculture Research Center (ARC), Egypt, during three growing winter seasons: 2019/2020, 2020/2021, and 2021/2022.Four genotypes of faba bean (Nubaria 5, Sakha 1, Giza 40 (belonging to V. faba subsp equina), and Marina (belonging to V. faba subsp minor)) were crossed in six population designs.The pedigree, origin, and reaction to foliar diseases of these genotypes are presented in Table 7.In the 2019/2020 season, the four faba bean genotypes were crossed in six populations fashioned under insect-free cages to obtain the F 1 hybrid seeds of three crosses as follows: (Nubaria 5 × Sakha 1), (Nubaria 5 × Giza 40), and (Marina × Giza 40).In the 2020/2021 sowing season, the F 1 seeds and their parents were planted, and F 1 plants were selfpollinated to obtain F 2 and backcrossed to both parents to obtain BC 1 and BC 2 for each cross.Re-hybridization was performed on the three crosses to obtain enough F 1 hybrid seeds under the same insect-free cage.
In the last season (2021/2022), a field experiment was constructed to evaluate the six populations' seeds (P 1 , P 2 , F 1 , F 2 , BC 1 , and BC 2 ) for the three crosses in a randomized complete block design (RCBD) in three replicates.Two ridges represented each entry for parents and F 1 s, eight ridges for F 2 s, and three ridges for each BC.The ridge was 3 m long and 60 cm apart.Seeds were planted on one side of the ridge at 20 cm hill spacing with one seed per hill.Data were recorded as an average of 15, 20, 100, and 30 on individual guarded plants from each plot for each cross's parents (P1 and P2), F 1 , F 2 , and BC (BC 1 and BC 2 ) generations, respectively.All recommendations related to faba bean agronomy and production applied during the growing seasons.Data were recorded on plant height (cm), number of branches/plant, number of pods/plant, number of seeds/plant, seed yield/plant (g), and foliar diseases (chocolate spot and rust).The chocolate spot and rust scale was 1-9 (1: highly resistant and 9: highly susceptible).Foliar disease reactions were measured according to Bernier et al. [65].

Statistical Analysis
Analysis of variance and mean performance of the studied traits measured the differences among the evaluated generations [66].The mean analysis among generations was performed using the method of Mather and Jinks [67] after the procedure scaling test to detect the presence of epistasis (non-allelic interaction).

Types of Gene Effect (Gene Action)
Data were analyzed to estimate the mean effects (m), additive (a), dominance (d), additive × additive (aa), additive × dominance (ad), and dominance × dominance (dd) via genetic analysis of the mean between generation according to Gamble [23].
Heterosis and inbreeding depression are calculated according to the method of Mather and Jinks [67], whereas the potency ratio was calculated based on the method of Petr and Frey [38] as follows: PR = F 1 − MP/0.5 HP − LP , where: F 1 = mean performance of F 1 ; MP = value of mid parent; HP = high parent value; LP = low parent value.Heritability was estimated in broad (h 2 b ) and narrow (h 2 n ) senses as mentioned by Mather [68], whereas predicted genetic advance (GS) was calculated as previously stated by Miller et al. [69].

DNA Extraction and SCoT Amplification Marker
To identify some molecular genetic markers associated with resistance and susceptibility to foliar diseases, the total genomic DNA was extracted by the DNA easy Plant Mini Kit (QIAGEN, Hilden, Germany) from the 14-day leaves of 33 lines [the parent Marina (resistant, coded 1), Giza 40 (susceptible, coded 2), the promising F 1 hybrid (Marina X Giza 40, coded 3), 10 selected lines from each of F 2 (coded 4-8 resistant, 9-11 moderate, and 12-13 susceptible), BC1 (Marina X Giza 40 with Marina, coded 14-18 resistant, 19-21 moderate, and 22-23 susceptible), and BC2 (Marina X Giza 40 with Giza 40, coded 24-28 resistant, 29-31 moderate and 32-33 susceptible)).A NanoDrop spectrophotometer was used to estimate the concentration and quality of extracted DNA.Seven SCoT primers were implied to detect polymorphism, as shown in Table S2.The amplification reaction was carried out in a 20 µL reaction volume containing 10 µL Master Mix (sigma), 2 µL primer (10 pcmol), 2 µL template DNA (10 ng), and 6 µL dH 2 O, according to Ibrahim et al. [67].PCR amplification was accomplished in a Perkin-Elmer/GeneAmp ® PCR System 9700 (PE Applied Biosystems) programmed to fulfill 40 cycles after an initial denaturation cycle for 5 min at 94 • C. Each cycle consisted of a denaturation step at 94 • C for 45 s, an annealing step at 50 • C for 50 s, and an elongation step at 72 • C for 1 min.The primer extension segment was extended to 7 min at 72 • C in the final cycle.The PCR products were resolved by electrophoresis in a 1.5% agarose gel containing ethidium bromide (0.5 µg/mL) in 1× TBE buffer at 95 volts.Gels containing bands were visualized and photographed under UV light using a Gel Documentation System (BIO-RAD 2000).

PCR Scoring and Data Analysis
For SCoT analysis, only clear and unambiguous bands were visually scored as either present (1) or absent (0) for all samples using Gel Analyzer software [70], and the final datasets included both polymorphic and monomorphic bands.
The Online Marker Efficiency Calculator (iMEC software) [71] was used to compute seven basic measures of polymorphism indices for individual markers, such as iMEC calculates heterozygosity index (H), polymorphism information content (PIC), discriminating power (D), effective multiplex ratio (E), marker index (MI), arithmetic mean heterozygosity (Hav), and resolving power (R).
Then, a binary statistic matrix was constructed and Dice's similarity matrix coefficients were calculated between lines using the unweighted pair group method with arithmetic averages (UPGMA) [72].Using the PAST software Version 1.91 [73], this matrix was used to create a phylogenetic tree (dendrogram) based on Jaccard similarity coefficients [74].

Conclusions
The present study used six populations (P 1 , P 2 , F 1 , F 2 , BC 1 , and BC 2 ) of three faba bean crosses (Nubaria 5 × Sakha 1; Nubaria 5 × Giza 40; and Marina × Giza 40).The scaling test showed that all of the evaluated traits significantly differed, indicating non-allelic interactions.BC 1 and BC 2 , in the first cross, verified the highest seed yield/plant, and they displayed the least infected plants in the third cross with respect to chocolate spot and rust diseases (negative values are desirable).Both heterosis types disclosed the desired positive heterotic effects regarding the number of seeds/plant, number of branches/plant, seed yield/plant, and pods/plant in all crosses, whereas the desired negative values of heterotic effects were detected with respect to chocolate spot and rust traits in the second and first crosses, respectively.Moreover, broad-sense heritability recorded the best values in the third cross for number of seeds/plant, whereas the first cross revealed the uppermost values for foliar diseases (chocolate spot and rust), respectively.The molecular profile for the cross Marina × Giza 40 and some selected lines of F 1 , F 2 , BC 1 , and BC 2 against seven SCoT primers generated 522 loci (357 polymorphic bands and 165 monomorphic bands), with average polymorphism (%) values reaching 74%.Three bands of molecular weight 260, 207, and 178 bp, generated via SCoT-1, SCoT-4, and SCoT-7 primers, respectively, were distinguished as being associated with the resistant parent (Marina).Based on the molecular profiles, genetic similarity among selected lines was estimated and found to range between 0.24 ang 0.96.The highest value (0.91) observed between Marina and BC1 (5) revealed a high resistance to foliar disease.Consequently, the SCoT markers can distinguish between lines for resistance or sensitivity to foliar diseases in developing an efficient breeding program.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12203645/s1,Table S1.Scaling test of the three crosses for the studied traits.Table S2.List of SCoT primers and their sequences.Table S3.Similarity coefficient (Dice measurement) of the selected genotypes based on SCoT banding profile.Figure S1.SCoT-PCR amplification using the selected primers.

Figure 2 .
Figure 2. Dendrogram for the 33 selected lines constructed from a similarity matrix compu cording to the Dice coefficient.

Figure 2 .
Figure 2. Dendrogram for the 33 selected lines constructed from a similarity matrix computed according to the Dice coefficient.

Table 1 .
Mean performance and variance for the six populations of the three crosses for the studied traits.BC 1 , and BC 2 in the first cross verified the highest seed yield/plant (110.94,121.14, 115.61, and 112.83 gm), respectively, whereas the third cross produced the minimum values of seed yield/plant.However, BC 1 and BC 2 displayed the plants least infected with chocolate spot and rust diseases, especially in the third cross.In general, F 1 mean value revealed the tallest plants (149.00 cm) and the highest number of seeds/plant (168.60) in the second cross; meanwhile, it had the highest number of branches/plant (7.00), pods/plant (49.47), and seed yield/plant (121.14) in the first cross.

Table 2 .
Heterosis, inbreeding depression, and potence ratio of the three crosses for the studied traits.
* and ** significant at 0.05 and 0.01 levels of probability, respectively.

Table 3 .
Gene action (Gamble's parameters) of the three crosses for the studied traits.
* and ** significant at 0.05 and 0.01 levels of probability, respectively.

Table 4 .
Estimates of heritability in broad (h 2 b ) and narrow (h 2 n ) senses plus genetic advance (GS) of the three crosses for the studied traits.

Table 5 .
Total bands (monomorphic and polymorphic), polymorphism percentage, and detected SCoT markers generated for resistance and susceptibility to foliar diseases using the seven SCoT markers across the selected faba bean lines.

Table 6 .
Indicators of genetic variability and information content of SCoT markers based on the analysis of the selected lines.

Table 7 .
Pedigree, origin, and reaction to foliar diseases of four parental faba bean genotypes under study.