Next Article in Journal
Can Policy Promote Agricultural Service Outsourcing? Quasi-Natural Experimental Evidence from China
Next Article in Special Issue
Legume Grains as an Alternative to Soybean Meal in the Diet of Intensively Reared Dairy Ewes
Previous Article in Journal
Persistence of Hg-Contaminated Soil Stabilization in Typical Areas of Dehua County, Fujian Province, China
Previous Article in Special Issue
Animals or Humans: What Do Greek Consumers Care More about When Buying Feta Cheese?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of the Genetic Diversity of Faba Bean (Vicia faba L.)

by
Eleni Avramidou
1,2,
Ioannis Ganopoulos
3,
Photini Mylona
3,
Eleni M. Abraham
2,
Irini Nianiou-Obeidat
4,
Maslin Osathanunkul
5,6 and
Panagiotis Madesis
1,7,*
1
Institute of Applied Biosciences (INAB), CERTH, 6th Km Charilaou-Thermi Road, 57001 Thessaloniki, Greece
2
Laboratory of Range Science, Faculty of Forestry and Natural Environment, Aristotle University, 55135 Thessaloniki, Greece
3
Hellenic Agricultural Organization Dimitra, Institute of Plant Breeding and Genetic Resources, 57001 Thermi, Greece
4
Laboratory of Genetics and Plant Breeding, School of Agriculture, Forestry and Natural Environment, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
5
Department of Biology, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
6
Research Center in Bioresources for Agriculture, Industry and Medicine, Chiang Mai University, Chiang Mai 50200, Thailand
7
Department of Agriculture Crop Production and Rural Environment, University of Thessaly, Fytokou St., 38446 Nea Ionia, Greece
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(2), 1016; https://0-doi-org.brum.beds.ac.uk/10.3390/su15021016
Submission received: 24 November 2022 / Revised: 27 December 2022 / Accepted: 29 December 2022 / Published: 5 January 2023
(This article belongs to the Collection Sustainable Livestock Production and Management)

Abstract

:
Faba bean (Vicia faba L.) is an important grain legume with high protein content (approximately 25–30%) and high nutritional value. It is broadly cultivated in temperate areas both for human consumption and as animal feed. According to FAOSTAT (2020), the total cultivated area of faba bean reached approximately 2.5 million ha, yielding more than 4.5 million tons. The characterization of the genetic diversity in faba bean is an important parameter for genetic and biodiversity studies, germplasm characterization, and for introducing genetic variability in plant breeding. The present study aims to assess the genetic diversity among 53 Greek, varied faba bean populations provided by the Hellenic Agricultural Organization “DEMETER” seed bank. To determine the genetic diversity of the studied populations, six SCoT DNA markers were used. A total of 114 loci were obtained with 37.95% being polymorphic and 62.05% monomorphic within or between populations. SCoT markers are a useful tool for the detection of genetic diversity among faba bean populations and encourage targeted crossing strategies. The present study is the first step towards the development of an efficient breeding program.

1. Introduction

Vicia faba, commonly named faba bean, is one of the oldest crops that is cultivated worldwide, and it is native to the Mediterranean–West Asia region [1,2]. Although the exact origin of faba bean is unknown, it is believed that it was one of the earliest food legumes to be domesticated since the Neolithic period [3]. Faba bean is one of the first domesticated food legumes and has a long history of cultivation [4]. The oldest known faba bean was first identified 14,000 years ago in the southern Levant [5].
Due to the high nutritional content of its seeds, faba bean is commonly utilized as a source of protein in the Mediterranean region for both human and animal nutrition [6]. Faba beans have a protein content of approximately 29% of the dry weight [7], which is higher compared to other common food legumes and makes it one of the main sources of protein for people in the Middle East, Latin America and Africa and for livestock feed in many developed countries [8]. According to [9], faba bean includes 13–14% of the cell wall slightly lignified, 3% saturated sugars, 6.5% oligosaccharides, and about 27–34% proteins, while starch makes up about 45% of the total dry mass. Additionally, apart from being a crucial food crop, it also contributes to the provision of animal feed and fodder and has a good impact on soil productivity for the cultivation of cereal crops [10] due to its ability to fix free nitrogen and to grow in different climatic zones [11]. In cooler climates, faba bean has advantages compared to other legumes such as soybean because it is better adapted to growth under low temperatures [12].
The use of legume crops in traditional farming systems forms a symbiosis with nodule-forming bacteria that have nitrogen-fixing ability, which provides major benefits to cropping systems and the environment and contributes to agricultural sustainability through soil improvement [4]. Compared to other grain legumes, faba bean is considered an excellent protein crop due to its ability to provide nitrogen inputs into temperate agricultural systems because of its wide adaptation [13,14].
According to FAOSTAT (2018), faba bean is the fourth most widely grown cool season legume after pea (Pisum sativum), chickpea (Cicer arietinum) and lentil (Lens culinaris Medik.) [15,16]. Only in 2016, 2.4 million ha of faba bean was harvested globally, and 4.5 million tons of dry grains was produced overall [17]. Faba bean production ranked seventh among all legumes crops worldwide [15].
Vicia faba is mainly a self-pollinating species. The authors of [18] stated that faba bean has one of the highest reported genome sizes among crop legumes, about 13,000 Mb, is a diploid with 2n = 2x = 12 chromosomes [19] and has a partial crosspollinated rate between 4 and 84% [20]. In comparison to the genome of the legume Medicago truncatula, this is 25 times larger [13].
Nevertheless, faba bean cultivars have antinutritional compounds such as tannins, primarily located in the seed coat, which limit their use for both human and animal consumption [21]. According to Zanollo et al. (2020), the presence of tannins decreases the digestibility and nutritional availability of protein, energy and starch in monogastric animals [21]. Faba bean tannins decrease protein digestibility in poultry due to the formation of tannin protein complexes. Tannins are also reported to reduce energy [22] and the digestibility of starch [23], but the effects were not always significant.
Analyzing genetic variation is a crucial component of genetic studies, biodiversity research, germplasm characterization and the creation of genetic variety in plant breeding. Moreover, it is important to estimate genetic variation in order to choose favorable genotypes. New techniques need to be created in order to estimate and utilize genetic variation in favor of a breeding program once it is realized that significant levels of genetic variation are not expressed in the phenotype. Recent technological developments in molecular genetics have made it possible to quantify genetic diversity at the DNA level by creating different molecular markers [2]. The major advantage of molecular markers is that they are not affected by environmental factors or by plant developmental stages [24].
Different types of markers, including isozymes, 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), have already been employed to assess the genetic variability of Vicia species and V. faba L. populations [2]. Intersimple sequence repeat (ISSR) markers were also utilized in order to explore the genetic diversity of V. faba accessions [25,26].
Recently, ref. [27] outlined a quick and innovative DNA marker technology called start codon targeted (SCoT) polymorphism. This marker was developed based on the short-conserved region flanking the Adenine-Thymine-Guanine (ATG) start codon in plant genes. SCoT markers are generally reproducible, and it is suggested that primer length and annealing temperature are not the sole factors determining reproducibility [28]. SCoT markers are similar to random amplified polymorphic DNA (RAPD) and intersimple sequence repeat (ISSR) because a single primer is used as forward and reverse. Moreover, as a PCR-based gene target technique, SCoT analysis has low cost and is effective to use [28]. SCoT primers were also used for the detection of the genetic diversity of Vicia sativa [29], chickpea (Cicer arietinum) [30], Dendrobium nobile L [28] and grape varieties (Vitis vinifera L.) [31].
The aim of the current study was to assess the genetic diversity for the first time among 53 Greek varied faba bean populations, provided by the Hellenic Agricultural Organization “DEMETER” seed bank, by using SCoT molecular markers. Detecting the genetic diversity among Greek faba bean populations is a crucial step in order to encourage targeted crossing strategies and to develop an efficient breeding program.

2. Materials and Methods

2.1. Plant Material

Twenty seeds per population, provided by the Hellenic Agricultural Organization “DEMETER” seed bank, were sown to pots in a controlled growth chamber, with a photoperiod of 14/10 h light/dark. The temperature was between 20 and 27 °C, with a mean temperature of 23 °C. Seeds were collected from different locations all over Greece (Figure 1, Table 1).

2.2. DNA Isolation and Marker Analysis

Three-week-old faba bean leaves from 53 selected genotypes were collected, and one gram of fresh young leaves of 53 V. faba populations was ground with liquid nitrogen and saved at −20 °C (five samples per population). Total genomic DNA isolation was carried out using the modified CTAB protocol (Doyle and Doyle 1991). The protocol was modified by the technicians of the lab, and to purify DNA, only chloroform was used. Additionally, centrifuge was carried out for 10 min at 16,000× g. The amount of DNA was quantified by a UV-Vis spectrophotometer (Q 5000), and then the samples were diluted to a 20 ng/μL working concentration.
PCR for SCoT analysis was performed in a total volume of 20 μL containing 20 ng total genomic DNA, 100 mM of each dNTP, 1.5 mM of Mg, 10 μΜ of primers and 5 u/μL of Kapa taq (KappaBiosystems, Bath, UK). To study the intra- and interpopulation diversity, 6 SCoT primers (SCoT13, SCoT14, SCoT15, SCoT61, SCoT66 and ScoT33) were used for PCR amplification (Table 2). After the review of the relevant literature, specific SCoT markers were chosen because they were found to be highly polymorphic for legume species and were also used in a variety of different plant species. PCR amplifications took place in SureCycler 8800 (Agilent Technologies) as follows: initial denaturation for 5 min at 94 °C followed by 40 cycles of 30 s at 94 °C, annealing at 50 °C for 90 s, and extension at 72 °C for 90 s. A 5 min step at 72 °C was programmed as a final extension.
Amplification products were separated by electrophoresis on 1.5% agarose gel and stained with ethidium bromide. Gel images were placed in a UVItec transilluminator, and a 100 bp or 1 Kb DNA ladder (Invitrogen, Carlsbad, CA, USA) was used as a size marker.
Binary data points denote the presence/absence of each distinguishable band across all samples for the same primer in both replicate sets of amplifications.

2.3. Data Analysis

SCoT markers are dominant markers, each band representing the phenotype at a single biallelic locus. Only bands that could be unambiguously scored were used in the analysis. SCoT amplified bands were scored for band presence (1) or absence (0), and a binary qualitative data matrix was formed [32].
The percentage of polymorphic loci (P), effective numbers of alleles (NE), gene diversity (expected heterozygosity, HE), Shannon’s diversity index (I) and unbiased genetic distances according to [33] were calculated using GenALEx ver.6.51b2 [32]. The hierarchical distribution of genetic variation among and within populations was also characterized by analysis of molecular variance (AMOVA) [3,34] using the GenALEx ver. 6.51b2 software [32], with variation being examined among and within populations. The tests were implemented using estimates of ΦST based on distances calculated from allelic data. Principal coordinate analysis (PCoA) was used to examine the genetic structure of V. faba populations using GenAlEx v. 6.51b2 [32] based on the standardized covariance of genetic distance for dominant markers. Additionally, Mantel tests [35] were conducted for the comparison of genetic and geographic distances using PopTools version 3.2.5 [36] with 99 iterations.

3. Results

A total of 114 loci were obtained from the six selected SCoT primers. Of all loci analyzed, 37.95% were polymorphic and 62.05% were monomorphic within or between populations. Gene diversity (GD) within the studied populations of faba bean ranged from 0.222 to 0.087, with a mean of 0.149 among the populations and the Shannon index (I) from 0.127 to 0.328 with a mean of 0.220 among the populations (Table 3). Among the populations, VF37 exhibited the highest level of polymorphism (P = 57.02%) and VF78 the lowest (P = 21.05%). In addition, VF37 and VF78 showed higher and smaller levels of genetic variation for Shannon’s information index (0.328 and 0.127, respectively) compared to other populations (Table 3).
The AMOVA analysis indicates that 56% of the total genetic variation was attributed to differences within populations, and the rest (44%) was attributed to differences among populations (Figure 2, Table 4).
Genetic differentiation (Nei’s genetic distance) [33] between faba bean populations ranged from 0.049 (VF3 to VF53 populations) to 0.655 (VF72 to VF39 populations) (Table 5).
The UPGMA dendrogram, based on the genetic distances between populations, depicted two main clusters, which are further subdivided. Cluster I and II grouped the highest number of accessions (eight), but cluster II was divided into two subclusters. The first subcluster contained eleven accessions (VF2, VF63, VF60, VF6, VF4, VF3, VF53, VF69, VF61, VF70 and VF37). Those subclusters are also divided into two clusters. The first one encompassed sixteen accessions (VF28, VF30, VF11, VF9, VF44, VF57, VF72, VF35, VF74, VF77, VF78, VF80, VF79, VF65, VF58 and VF62). The second on encompassed thirteen accessions (VF29, VF33, VF1, VF25, VF42, VF43, VF47, VF48, VF67, VF76, VF66, VF55 and VF56) (Figure 3). However, it is not able to discriminate faba bean accessions on the basis of their geographical origin. The only obvious correlation is among the populations VF77, VF78, VF79 and VF80, which are cultivars varieties. PCoA was used to examine associations among 53 accessions. The first principal component (PC1) and the second principal component (PC2) specified the localization of individuals. PC1 and PC2 in this analysis explained 20.39% and 15.30% of the overall variability, respectively. A different perspective on the genetic distances of faba bean accessions is provided by the PCoA (Figure 4). The plant species formed clusters (Figure 3), which is in accordance with the UPGMA dendrogram. PCoA was also unable to distinguish faba bean accessions according to their geographical origin and also confirmed the results observed in the UPGMA dendrogram.
Finally, the Mantel test showed a relatively high correlation between the genetic and the geographic distances (R = 0.878, p < 0.05, 99 iterations).

4. Discussion

The knowledge of genetic variation within collections of Greek faba bean genetic resources is crucial for the effective conservation and utilization of these resources in breeding programs and could be dramatically enhanced by using molecular markers [2,37].
Restriction fragment length polymorphism (RFLP) markers [38], random amplified polymorphic DNA (RAPD) markers [37], amplified fragment length polymorphism (AFLP) markers [39] and intersimple sequence repeats (ISSR) markers [25] were successfully used to characterize the genetic variability among faba bean genotypes [2,40]. This is the first attempt at exploring the genetic diversity of Greek faba bean populations using SCoT molecular markers.
Several types of molecular markers were utilized in order to assess the genetic diversity of legume crops and especially faba bean. In the present study, we used six SCoT molecular markers to characterize genetic diversity in 53 Greek faba bean populations. In this study, 37.95% of the alleles were polymorphic, which is much lower compared to the results that [25] found studying the genetic diversity of 20 Greek local faba bean populations by using ISSR markers (98.9%). Moreover, ref. [26] studying the genetic diversity of Ethiopian faba bean varieties by using ISSR markers, reported that the percentage of polymorphism was 85.5%. Additionally, the percentage of polymorphism of the present study is much lower compared to SRAP markers (96% and 98.2%, respectively) [37,41], SSR markers (99% and 87.8%, respectively) [18,42] and RAPD markers (75.96%) [43]. It is obvious that different kinds of molecular markers affect the manner of polymorphism [44]. The authors of [45] stated that genetic diversity within populations is also influenced by other factors apart from the selection of molecular markers, such as long periods of low population density or gene flow rates. Additionally, low genetic differentiation within populations can result in using either pollen or seeds or long-distance gene transmission [45]. This transmission can be explained by grazing from wild or domesticated animals. Faba bean is an outcrossing crop, with the ratio of outcrossing differing between environments and genotypes. Considering that bee species are the major vectors for pollen transfer, the rate of crosspollination varies across environments, resulting in the loss of genetic diversity [46].
Furthermore, regarding the Shannon diversity index, SCoT analysis revealed an average of 0.220, while in other studies, the average was 0.63 for SSR markers [42], 0.41 for ISSR markers [26] and 0.29 for RAPD markers [43]. The results of our study are really close to the results reported by [47]. They studied the genetic diversity of some species of genus Vicia using ISSR and ITS molecular techniques. Shannon diversity was estimated to be 0.24. The Shannon diversity index was also estimated to be 0.262 by [48] in a study that took place in Vicia species using SSAP markers from 56 accessions from Spain and Syria.
AMOVA showed that 56% of the total genetic variation was attributed to differences within populations and that the rest (44%) was attributed to differences among populations. SCoT markers are not such polymorphisms as ISSRs, SSRs or RAPDs markers.
In contrary to the results of this study, all other studies exhibited a high level of genetic variability within populations using ISSRs (75.4%) [25] for local Greek faba bean populations, SSRs (90%) for faba bean accessions from Turkey, Sweden, Australia, Finland and Egypt [18] and RAPDs (94%) for Tunisian populations [43]. Additionally, ref. [49] reported that after evaluating the genetic diversity of 65 accessions from 12 Vicia species from the National Grass Germplasm resource bank, using SSR molecular markers, genetic variability within populations was 89% and among populations was 11%.
The same results were also found by [50], who studied the genetic diversity among 20 pea varieties and 57 accessions from wild Pisum species, and AMOVA revealed the intergroup component of variance to be 29%, while the intragroup component was 71%.
The results of the present study are not in accordance with the prementioned studies, suggesting a level of heterogeneity that is higher than expected for faba bean since it is mainly self-pollinating. Yet, crosspollination and hybrid formation has been reported to result in heterogeneous mixtures of inbreeds and hybrids [37]. In addition, the possible explanation of the observed geographical ‘mixed’ accessions in these populations (Figure 3 and Figure 4) could be their moving from one place to another by human activities or by insect pollinators [51]. Grazing by wild or domesticated animals is the main human activity for the transaction of seeds and, consequently, the transaction of gene flow and exchange of genotypes in one or more cases. This is probably the reason that populations do not group according to geographical origin. This is in line with [39], who could not identify geographic partitioning of diversity in their analysis because of the low number of V. faba lines (n = 20) studied. However, the knowledge of the genetic diversity of the population under study is invaluable for the Greek Seed Bank collection as well as for faba breeders. The results obtained here provide important information on the genetic background of the Greek faba germplasm in the collection and offer at the same time the opportunity for the development of new collection and breeding strategies.
The authors of [52] studied the genetic diversity in 28 pea (Pisum sativum L.) genotypes using SSR markers from the Department of Genetics and Plant Breeding Institute of Agricultural Sciences in India. The 28 pea genotypes had a coefficient of genetic similarity that ranged from 0.11 to 0.73, demonstrating a high degree of genetic variety. In contrast, [53] found little variation (0.69–0.88) between cultivars of P. sativum ssp sativum and P. sativum ssp arvense. In comparison to results obtained with AFLP markers, ref. [54] showed a substantially greater similarity range (0.80–0.94) with RAPD markers in pea cultivars (0.85–0.94). Additionally, they reported that the estimated genetic diversity (0.05–0.82) among pea accessions in the study was higher than the one reported by [55] (0.0–0.66) and [56] (0.05–0.48) but comparable to that published by [57] (0.24–0.84). The greater estimated genetic distance could be attributed to variations between accessions as a result of genotype diversity in their ancestry.
According to UPGMA cluster analysis and the PCoA, the studied 53 faba bean populations were genetically evenly distributed regardless of their geographic origin, which probably suggests the mix of the genetic germplasm at some point, probably due to travel or commerce. Genetically distinguishable were the populations VF77, VF78, VF79 and VF80, which are cultivars, and they formed a group. Similar results were described in previous studies [58,59,60]. According to Malek et al. (2021) who studied 14 Algerian faba bean accessions, the UPGMA dendrogram and PCoA failed to classify the accessions according to their place of origin [58]. Additionally, the authors of [60] who studied 16 Tunisia faba bean local populations also did not reveal any population and geographical area interaction. These results show that there was significant gene flow among the communities that were under study. As a result, in addition to outcrossing, human activities may affect how genetic diversity and structure are distributed in faba bean germplasm. A partial sharing of their ancestral genetic polymorphism could probably account for the overall inability of the analyzed faba bean accessions to cluster according to their place of origin. However, the Mantel test showed a relatively high correlation between the genetic and the geographic distances (R = 0.878, p < 0.05, 99 iterations), which is not depicted in PCoA analysis and by the UPGMA dendrogram.
Management programs of faba bean should be put in place in order to maintain its genetic diversity at the best possible level. Thus, in future programs, there should be a sampling plant set up for as many individuals as possible, with emphasis on diverse and isolated local populations. We should also stress the usefulness of and need for molecular tools and especially molecular markers as a means to monitor small populations for which pedigrees are not always available. Moreover, Vicia faba is a traditional legume cultivated in varied environmental conditions. Thus, it becomes even more important for V. faba to identify stable and high-yielding genotypes from the collection expeditions in extremely different environments, having in mind the relatively low genetic diversity of the Greek genotypes. Furthermore, it is expected that the incorporation of the genomic studies of V. faba along with other relative species and the discovery of SNPS, the relative high-density SNP map and the synteny among them will offer new tools of high power for the genetic studies of V. fava.

5. Conclusions

Our results with SCoT markers suggest that the 53 Greek faba bean populations utilized in this study are genetically diverse and that SCoT molecular markers are reliable and powerful tools to evaluate genetic polymorphisms and relationships among faba bean genotypes in order to develop faba bean breeding programs. However, more research should be performed with the design of an expedition for the collection of more genetic resources from more extreme environments in Greece, and furthermore, more research should study the genetic makeup using diverse germplasm, also using RILs around the world with even more markers in order to fully understand the genetic makeup of the Greek varieties and develop an efficient breeding program towards higher yield with less inputs. Yet, the results presented here give the relevant scientific community the means to already start designing and implementing an efficient breeding program.
SCoT molecular markers are innovative compared with other molecular markers since they have many advantages: they have simple operation, they do not cost much, they have an abundant rate of polymorphism, and they are more conducive to operation for molecular-assisted breeding. This study was the first attempt to explore the genetic diversity of Greek faba bean populations using SCoT markers and develop an efficient breeding program. SCoT markers revealed a relatively high percentage of heterozygosity, which facilitates the development of an efficient breeding program. Such variation has profound implications in practical breeding and should be further exploited in order to develop synthetic varieties for improving faba bean productivity.

Author Contributions

Conceptualization, P.M. (Panagiotis Madesis), E.A. and E.M.A.; methodology, P.M. (Panagiotis Madesis), M.O., E.A., I.G., I.N.-O. and P.M. (Photini Mylona); validation, E.A., I.N.-O., I.G., E.M.A. and P.M. (Photini Mylona); resources, P.M. (Panagiotis Madesis), P.M. (Photini Mylona) and E.A.; writing—original draft preparation, E.A. and P.M. (Panagiotis Madesis); writing—review and editing, P.M. (Panagiotis Madesis), M.O., E.A., I.G., I.N.-O. and P.M. (Photini Mylona); visualization, E.A.; supervision, P.M. (Panagiotis Madesis); project administration, P.M. (Panagiotis Madesis); funding acquisition, P.M. (Panagiotis Madesis) and E.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially founded by Chiang Mai University and cofinanced by the European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation under the call RESEARCH-CREATE-INNOVATE (project code: Τ1EDK-04448).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We would like to thank Chiang Mai University for the support and NAGREF—Greek Seed Bank—for the seeds provided.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tanno, K.I.; Willcox, G. The Origins of Cultivation of Cicer arietinum L. and Vicia faba L.: Early Finds from Tell El-Kerkh, North-West Syria, Late 10th Millennium B.P. Veg. Hist. Archaeobotany 2006, 15, 197–204. [Google Scholar] [CrossRef]
  2. Abid, G.; Mingeot, D.; Udupa, S.M.; Muhovski, Y.; Watillon, B.; Sassi, K.; M’hamdi, M.; Souissi, F.; Mannai, K.; Barhoumi, F.; et al. Genetic Relationship and Diversity Analysis of Faba Bean (Vicia faba L. Var. Minor) Genetic Resources Using Morphological and Microsatellite Molecular Markers. Plant Mol. Biol. Rep. 2015, 33, 1755–1767. [Google Scholar] [CrossRef]
  3. Torres, A.M.; Avila, C.M.; Gutierrez, N.; Palomino, C.; Moreno, M.T.; Cubero, J.I. Marker-Assisted Selection in Faba Bean (Vicia faba L.). Field Crop. Res. 2010, 115, 243–252. [Google Scholar] [CrossRef]
  4. Khazaei, H.; O’Sullivan, D.M.; Stoddard, F.L.; Adhikari, K.N.; Paull, J.G.; Schulman, A.H.; Andersen, S.U.; Vandenberg, A. Recent Advances in Faba Bean Genetic and Genomic Tools for Crop Improvement. Legum. Sci. 2021, 3, e75. [Google Scholar] [CrossRef] [PubMed]
  5. Caracuta, V.; Weinstein-Evron, M.; Kaufman, D.; Yeshurun, R.; Silvent, J.; Boaretto, E. 4000-Year-Old Seeds Indicate the Levantine Origin of the Lost Progenitor of Faba Bean OPEN. Sci. Rep. 2016, 6, 37399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Mejri, S.; Mabrouk, Y.; Voisin, M.; Delavault, P.; Simier, P.; Saidi, M.; Belhadj, O. Variation in Quantitative Characters of Faba Bean after Seed Irradiation and Associated Molecular Changes. Afr. J. Biotechnol. 2012, 11, 8383–8390. [Google Scholar] [CrossRef]
  7. Warsame, A.O.; O’Sullivan, D.M.; Tosi, P. Seed storage proteins of faba bean (Vicia faba L.): Current status and prospects for genetic improvement. J. Agric. Food Chem. 2018, 66, 12617–12626. [Google Scholar] [CrossRef]
  8. Burstin, J.; Gallardo, K.; Mir, R.R.; Varshney, R.K.; Duc, G. Improving Protein Content and Nutrition Quality. Biol. Breed. Food Legum. 2011, 314–328. [Google Scholar] [CrossRef] [Green Version]
  9. Ayed, M.H. Effect of Partially Replacing Soybean by Broad Beans “Vicia faba Var Minor” on Digestibility, Growth Performances and Carcass Characteristics in Broiler Chickens. Adv. Environ. Biol. 2011, 5, 2332–2335. [Google Scholar]
  10. Duc, G. Faba Bean (Vicia faba L.). Field Crop. Res. 1997, 53, 99–109. [Google Scholar] [CrossRef]
  11. Gutierrez, N.; Torres, A.M. Characterization and Diagnostic Marker for TTG1 Regulating Tannin and Anthocyanin Biosynthesis in Faba Bean. Sci. Rep. 2019, 9, 16174. [Google Scholar] [CrossRef]
  12. Cooper, J.W.; Wilson, M.H.; Derks, M.F.L.; Smit, S.; Kunert, K.J.; Cullis, C.; Foyer, C.H. Enhancing Faba Bean (Vicia faba L.) Genome Resources. J. Exp. Bot. 2017, 68, 1941–1953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rispail, N.; Kaló, P.; Kiss, G.B.; Ellis, T.H.N.; Gallardo, K.; Thompson, R.D.; Prats, E.; Larrainzar, E.; Ladrera, R.; González, E.M.; et al. Model Legumes Contribute to Faba Bean Breeding. Field Crop. Res. 2010, 115, 253–269. [Google Scholar] [CrossRef] [Green Version]
  14. Cernay, C.; Ben-Ari, T.; Pelzer, E.; Meynard, J.-M.; Makowski, D. Estimating Variability in Grain Legume Yields across Europe and the Americas. Sci. Rep. 2015, 5, 11171. [Google Scholar] [CrossRef] [Green Version]
  15. Maalouf, F.; Hu, J.; O’Sullivan, D.M.; Zong, X.; Hamwieh, A.; Kumar, S.; Baum, M. Breeding and Genomics Status in Faba Bean (Vicia faba). Plant Breed. 2019, 138, 465–473. [Google Scholar] [CrossRef] [Green Version]
  16. Papastylianou, P.; Vlachostergios, D.N.; Dordas, C.; Tigka, E.; Papakaloudis, P.; Kargiotidou, A.; Pratsinakis, E.; Koskosidis, A.; Pankou, C.; Kousta, A.; et al. Genotype × Environment Interaction Analysis of Faba Bean (Vicia faba L.) for Biomass and Seed Yield across Different Environments. Sustainability 2021, 13, 2586. [Google Scholar] [CrossRef]
  17. FAOSTAT. Production Statistics of the Food Agriculture Organization of The United States. Available online: http://www.fao.org/faostat/en/#data/QA (accessed on 3 June 2021).
  18. Göl, Ş.; Doğanlar, S.; Frary, A. Relationship between Geographical Origin, Seed Size and Genetic Diversity in Faba Bean (Vicia faba L.) as Revealed by SSR Markers. Mol. Genet. Genom. 2017, 292, 991–999. [Google Scholar] [CrossRef] [PubMed]
  19. Johnston, J.S.; Bennett, M.D.; Rayburn, A.L.; Galbraith, D.W.; Price, H.J. Reference Standards for Determination of DNA Content of Plant Nuclei. Am. J. Bot. 1999, 86, 609–613. [Google Scholar] [CrossRef]
  20. Ellwood, S.R.; Phan, H.T.T.; Jordan, M.; Hane, J.; Torres, A.M.; Avila, C.M.; Cruz-Izquierdo, S.; Oliver, R.P. Construction of a Comparative Genetic Map in Faba Bean (Vicia faba L.); Conservation of Genome Structure with Lens Culinaris. BMC Genom. 2008, 9, 380. [Google Scholar] [CrossRef] [Green Version]
  21. Zanotto, S.; Khazaei, H.; Elessawy, F.M.; Vandenberg, A.; Purves, R.W. Do Faba Bean Genotypes Carrying Different Zero-Tannin Genes (Zt1 and Zt2) Differ in Phenolic Profiles? J. Agric. Food Chem. 2020, 68, 7530–7540. [Google Scholar] [CrossRef]
  22. Vilariño, M.; Métayer, J.P.; Crépon, K.; Duc, G. Effects of Varying Vicine, Convicine and Tannin Contents of Faba Bean Seeds (Vicia faba L.) on Nutritional Values for Broiler Chicken. Anim. Feed Sci. Technol. 2009, 150, 114–121. [Google Scholar] [CrossRef]
  23. Lacassagne, L.; Francesch, M.; Carré, B.; Melcion, J.P. Utilization of Tannin-Containing and Tannin-Free Faba Beans (Vicia faba) by Young Chicks: Effects of Pelleting Feeds on Energy, Protein and Starch Digestibility. Anim. Feed Sci. Technol. 1988, 20, 59–68. [Google Scholar] [CrossRef] [Green Version]
  24. Abdel-Razzak, H.S.; Alfrmawy, A.M.; Ibrahim, H.M.; El-Hanafy, A.A. Genetic Diversity in Faba Bean (Vicia faba L.) Using Inter-Simple Sequence Repeat (Issr) Markers and Protein Analysis. Life Sci. J. 2012, 9, 497–503. [Google Scholar]
  25. Terzopoulos, P.J.; Bebeli, P.J. Genetic Diversity Analysis of Mediterranean Faba Bean (Vicia faba L.) with ISSR Markers. Field Crop. Res. 2008, 108, 39–44. [Google Scholar] [CrossRef]
  26. Asfaw, B.M.; Dagne, K.; Wakayo, G.K.; Kemal, S.A.; Muleta, K.T. Genetic Diversity Study of Ethiopian Faba Bean (Vicia faba L.) Varieties Based on Phenotypic Traits and Inter Simple Sequence Repeat (ISSR) Markers. Afr. J. Biotechnol. 2018, 17, 433–446. [Google Scholar] [CrossRef]
  27. Collard, B.C.Y.; Mackill, D.J. Start Codon Targeted (SCoT) Polymorphism: A Simple, Novel DNA Marker Technique for Generating Gene-Targeted Markers in Plants. Plant Mol. Biol. Rep. 2009, 27, 86–93. [Google Scholar] [CrossRef]
  28. Bhattacharyya, P.; Kumaria, S.; Kumar, S.; Tandon, P. Start Codon Targeted (SCoT) Marker Reveals Genetic Diversity of Dendrobium Nobile Lindl., an Endangered Medicinal Orchid Species. Gene 2013, 529, 21–26. [Google Scholar] [CrossRef] [PubMed]
  29. Chai, X.; Dong, R.; Liu, W.; Wang, Y.; Liu, Z. Optimizing Sample Size to Assess the Genetic Diversity in Common Vetch (Vicia sativa L.) Populations Using Start Codon Targeted (SCoT) Markers. Molecules 2017, 22, 567. [Google Scholar] [CrossRef] [Green Version]
  30. Pakseresht, F.; Talebi, R.; Karami, E. Comparative Assessment of ISSR, DAMD and SCoT Markers for Evaluation of Genetic Diversity and Conservation of Landrace Chickpea (Cicer arietinum L.) Genotypes Collected from North-West of Iran. Physiol. Mol. Biol. Plants 2013, 19, 563–574. [Google Scholar] [CrossRef] [Green Version]
  31. Guo, D.L.; Zhang, J.Y.; Liu, C.H. Genetic Diversity in Some Grape Varieties Revealed by SCoT Analyses. Mol. Biol. Rep. 2012, 39, 5307–5313. [Google Scholar] [CrossRef]
  32. Peakall, R.O.D.; Smouse, P.E. GENALEX 6: Genetic Analysis in Excel. Population Genetic Software for Teaching and Research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  33. Ne, M. Estimation of Average Heterozygosity and genetic distance from a small number of individuals. Genetics 1978, 89, 583–590. [Google Scholar] [CrossRef] [PubMed]
  34. Excoffier, L.; Smouse, P.E.; Quattro, J.M. Analysis of Molecular Variance Inferred from Metric Distances among DNA Haplotypes: Application to Human Mitochondrial DNA Restriction Data. Genetics 1992, 131, 479–491. [Google Scholar] [CrossRef] [PubMed]
  35. Jensen, J.L.; Bohonak, A.J.; Kelley, S.T. Isolation by Distance, Web Service. BMC Genet. 2005, 6, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hood, G.M. PopTools, Version 3.2.5. Pest Animal Control CRC, Canberra. 2010. Available online: http://www.poptools.org (accessed on 20 July 2022).
  37. Alghamdi, S.S.; Al-Faifi, S.A.; Migdadi, H.M.; Khan, M.A.; El-Harty, E.H.; Ammar, M.H. Molecular Diversity Assessment Using Sequence Related Amplified Polymorphism (SRAP) Markers in Vicia faba L. Int. J. Mol. Sci. 2012, 13, 16457–16471. [Google Scholar] [CrossRef] [PubMed]
  38. Van De Ven, M.; Powell, W.; Ramsay, G.; Waugh, R. Restriction Fragment Length Polymorphisms as Genetic Markers in Vicia. Heredity 1990, 65, 329–342. [Google Scholar] [CrossRef] [Green Version]
  39. Zeid, M.; Schön, C.-C.; Link, W. Genetic Diversity in Recent Elite Faba Bean Lines Using AFLP Markers. Theor. Appl. Genet. 2003, 107, 1304–1314. [Google Scholar] [CrossRef] [PubMed]
  40. Duc, G.; Bao, S.; Baum, M.; Redden, B.; Sadiki, M.; Suso, M.J.; Vishniakova, M.; Zong, X. Diversity Maintenance and Use of Vicia faba L. Genetic Resources. Field Crop. Res. 2010, 115, 270–278. [Google Scholar] [CrossRef] [Green Version]
  41. Elshafei, A.A.M.; Amer, M.A.E.; Elenany, M.A.M.; Helal, A.G.A.E. Evaluation of the Genetic Variability of Faba Bean (Vicia faba L.) Genotypes Using Agronomic Traits and Molecular Markers. Bull. Natl. Res. Cent. 2019, 43, 106. [Google Scholar] [CrossRef] [Green Version]
  42. Tufan, H.; Erdoğan, C. Genetic Diversity in Some Faba Bean (Vicia faba L.) Genotypes Assessed by Simple Sequence Repeats. Biotechnol. Biotechnol. Equip. 2017, 31, 29–35. [Google Scholar] [CrossRef] [Green Version]
  43. Backouchi, I.Z.; Aouida, M.; Khemiri, N.; Jebara, M. Genetic Diversity in Tunisian Populations of Faba Bean (Vicia faba L.) Based on Morphological Traits and Molecular Markers. Genet. Mol. Res. 2015, 14, 7587–7596. [Google Scholar] [CrossRef] [PubMed]
  44. Stathi, E.; Kougioumoutzis, K.; Abraham, E.M.; Trigas, P.; Ganopoulos, I.; Avramidou, E.V.; Tani, E. Population Genetic Variability and Distribution of the Endangered Greek Endemic Cicer Graecum under Climate Change Scenarios. AoB Plants 2021, 12, plaa007. [Google Scholar] [CrossRef] [PubMed]
  45. Zhao, X.; Ma, Y.; Sun, W.; Wen, X.; Milne, R. High Genetic Diversity and Low Differentiation of Michelia Coriacea (Magnoliaceae), a Critically Endangered Endemic in Southeast Yunnan, China. Int. J. Mol. Sci. 2012, 13, 4396–4411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Gnanasambandam, A.; Paull, J.; Torres, A.; Kaur, S.; Leonforte, T.; Li, H.; Zong, X.; Yang, T.; Materne, M. Impact of Molecular Technologies on Faba Bean (Vicia faba L.) Breeding Strategies. Agronomy 2012, 2, 132–166. [Google Scholar] [CrossRef] [Green Version]
  47. El-Badan, G.E.; Amin, A.W.; Ashour, F.M.; El-Sadek, L.M. Characterization of the Genetic Diversity of Some Species of Genus Vicia Using ISSR and ITS Molecular Techniques. Egypt. J. Bot. 2022, 62, 475–492. [Google Scholar] [CrossRef]
  48. Martín Sanz, A.; Gilsanz Gonzalez, S.; Syed, N.H.; Suso, M.J.; Saldaña, G.C.; Flavell, A.J. Genetic Diversity Analysis in Vicia Species Using Retrotransposon-Based SSAP Markers. Mol. Genet. Genom. 2007, 278, 433–441. [Google Scholar] [CrossRef]
  49. Sun, W.-H.; Wu, F.-F.; Cong, L.; Jin, M.-Y.; Wang, X.-G. Assessment of Genetic Diversity and Population Structure of the Genus Vicia (Vicia L.) Using Simple Sequence Repeat Markers. Grassl. Sci. 2022, 68, 205–213. [Google Scholar] [CrossRef]
  50. Nasiri, J.; Haghnazari, A.; Saba, J. Genetic Diversity among Varieties and Wild Species Accessions of Pea (Pisum sativum L.) Based on SSR Markers. Afr. J. Biotechnol. 2009, 8, 3405–3417. [Google Scholar]
  51. Smýkal, P.; Horáček, J.; Dostálová, R.; Hýbl, M. Variety Discrimination in Pea (Pisum sativum L.) by Molecular, Biochemical and Morphological Markers. J. Appl. Genet. 2008, 49, 155–166. [Google Scholar] [CrossRef]
  52. Kumari, P.; Basal, N.; Singh, A.K.; Rai, V.P.; Srivastava, C.P.; Singh, P.K. Genetic Diversity Studies in Pea (Pisum sativum L.) Using Simple Sequence Repeat Markers. Genet. Mol. Res. 2013, 12, 3540–3550. [Google Scholar] [CrossRef]
  53. Samec, P.; Naginec, V. The Use of RAPD Technique for the Identification and Classification of Pisum sativum L. Genotypes. Euphytica 1996, 89, 229–234. [Google Scholar] [CrossRef]
  54. Simioniuc, D.; Uptmoor, R.; Friedt, W.; Rd On, F.O. Genetic Diversity and Relationships among Pea Cultivars Revealed by RAPDs and AFLPs. Plant Breed. 2008, 121, 429–435. [Google Scholar] [CrossRef]
  55. Tar’an, B.; Zhang, C.; Warkentin, T.; Tullu, A.; Vandenberg, A. Genetic diversity among varieties and wild species accessions of pea (Pisum sativum L.) based on molecular markers, and morphological and physiological characters. Genome 2005, 48, 257–272. [Google Scholar] [CrossRef]
  56. Ford, R.; Le Roux, K.; Itman, C.; Brouwer, J.B.; Taylor, P.W.J. Diversity Analysis and Genotyping in Pisum with Sequence Tagged Microsatellite Site (STMS) Primers. Euphytica 2002, 124, 397–405. [Google Scholar] [CrossRef]
  57. Cupic, T.; Tucak, M.; Popovic, S.; Bolaric, S.; Grljusic, S.; Kozumplik, V. Genetic Diversity of Pea (Pisum sativum L.) Genotypes Assessed by Pedigree, Morphological and Molecular Data. J. Food Agric. Environ. 2009, 7, 343–348. [Google Scholar]
  58. Malek, N.; Aci, M.M.; Khamassi, K.; Lupini, A.; Rouaissi, M.; Hanifi-Mekliche, L. Agro-Morphological and Molecular Variability among Algerian Faba Bean (Vicia faba L.) Accessions. Agronomy 2021, 11, 1456. [Google Scholar] [CrossRef]
  59. Rebaa, F.; Abid, G.; Aouida, M.; Abdelkarim, S.; Aroua, I.; Muhovski, Y.; Baudoin, J.P.; M’hamdi, M.; Sassi, K.; Jebara, M. Genetic Variability in Tunisian Populations of Faba Bean (Vicia faba L. Var. Major) Assessed by Morphological and SSR Markers. Physiol. Mol. Biol. Plants 2017, 23, 397–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Yahia, Y.; Hannachi, H.; Monforte, A.J.; Cockram, J.; Loumerem, M.; Zarouri, B.; Ferchichi, A. Genetic Diversity in Vicia faba L. Populations Cultivated in Tunisia Revealed by Simple Sequence Repeat Analysis. Plant Genet. Resour. Characterisation Util. 2014, 12, 278–285. [Google Scholar] [CrossRef]
Figure 1. Geographical distribution of faba bean accessions from Greece.
Figure 1. Geographical distribution of faba bean accessions from Greece.
Sustainability 15 01016 g001
Figure 2. Distribution of percentages of molecular variance among fifty-three faba bean populations.
Figure 2. Distribution of percentages of molecular variance among fifty-three faba bean populations.
Sustainability 15 01016 g002
Figure 3. UPGMA dendrogram based on Nei’s genetic distance of 53 faba bean populations.
Figure 3. UPGMA dendrogram based on Nei’s genetic distance of 53 faba bean populations.
Sustainability 15 01016 g003
Figure 4. Principal coordinate analysis (PCoA) of 53 faba bean accessions using 6 SCoT markers.
Figure 4. Principal coordinate analysis (PCoA) of 53 faba bean accessions using 6 SCoT markers.
Sustainability 15 01016 g004
Table 1. Data of faba bean accessions provided by the Hellenic Agricultural Organization “DEMETER” seed bank.
Table 1. Data of faba bean accessions provided by the Hellenic Agricultural Organization “DEMETER” seed bank.
SampleCodeSpeciesExpeditionCollection StateCollection Site
VF1IS-153/07V. fabaIKARIA-SAMOSSAMOS ISL.VOURLIOTES (37°47′6″ Ν, 26°50′52.8″ Ε)
VF2KThI-052/07V. fabaKERKIRA-THESPROTIA-IOANNINATHESPROTIAFROSINI (39°27′09.13″ Ν, 20°35′51.67″ Ε)
VF3IS-212/07V. fabaIKARIA-SAMOSSAMOS ISL.AMMOUDIES (37°37′42.08″ Ν, 26°47′36.79″ Ε)
VF4AO-060/07V. fabaAG. OROSAG. OROSKELI AG. MITROFANOUS (40°15.26″ Ν, 24°14′42″ Ε)
VF6KThI-051/07V. fabaKERKIRA-THESPROTIA-IOANNINATHESPROTIATSAGGARI (39°25.04″ Ν, 20°36′39″ Ε)
VF9AO-033/07V. fabaAG. OROSAG. OROSKELI AG. SKEPIS (40°15.26″ Ν, 24°14′42″ Ε)
VF11HL-158/07V. fabaHERAKLEION-LASITHILASITHIPINAKIANO (35°11′50.4″ Ν, 25°27′52.2″ Ε)
VF15HL-203/07V. fabaHERAKLEION-LASITHILASITHISXINOKAPSALA (35°3′10.8″ Ν, 25°52′58.8″ Ε)
VF17HL-077/07V. spHERAKLEION-LASITHIIRAKLIONAG. VARVARA (35°9′0″ Ν, 25°3′0″ Ε)
VF25ANP-158/07V. fabaAMORGOS-NAXOS-PAROSNAXOS ISL.MIRISI (37°9′19.9″ Ν, 25°30′35.5″ Ε)
VF27ANP-219/07V. fabaAMORGOS-NAXOS-PAROSPAROS ISLLEUKES (37°03′21″ Ν, 25°12′31″ Ε)
VF28ANP-086/07V. fabaAMORGOS-NAXOS-PAROSAMORGOS ISL.THOLARIA (37°47′6″ Ν, 26°50′52.8″ Ε)
VF29ANP-121/07V. fabaAMORGOS-NAXOS-PAROSNAXOS ISL.AG. ARSENIOS (37°3′39″ Ν, 25°23′27″ Ε)
VF30ANP-066/07V. fabaAMORGOS-NAXOS-PAROSAMORGOS ISL.KOLOFANA (36°50′00″ Ν, 25°54′00″ Ε)
VF31ANP-169/07V. fabaAMORGOS-NAXOS-PAROSNAXOS ISL.APIRANTHOS (37°04′21″ Ν, 25°31′19″ Ε)
VF33ANP-171/07V. fabaAMORGOS-NAXOS-PAROSNAXOS ISL.APIRANTHOS (37°04′21″ Ν, 25°31′19″ Ε)
VF35ROX-031/07V. fabaRODOPI-XANTHIXANTHIMYRODATO (40°58′33.5″ Ν, 24°55′36.09″ Ε)
VF36ROX-145/07V. fabaRODOPI-XANTHIRODOPIP.KROVILI (40°57′33″ Ν, 25°33′30″ Ε)
VF37ROX-019/07V. fabaRODOPI-XANTHIXANTHIMELISSA (41°01′01″ Ν, 24°53′42″ Ε)
VF38ROX-003/07V. spRODOPI-XANTHIXANTHIAG. ATHANASIOS (41°02′24.00″ Ν, 24°46′58.08″ Ε)
VF39MFS-022/07V. fabaMILOS-FOLEGANDROS-SIKINOSMILOS ISL.ZEFIRIA (36°42′3″ Ν, 24°29′26″ Ε)
VF40MFS-060/07V. fabaMILOS-FOLEGANDROS-SIKINOSMILOS ISL.DASIFNOS (36°44′38″ Ν, 24°25′21″ Ε)
VF41MFS-100/07V. fabaMILOS-FOLEGANDROS-SIKINOSFOLEGANDROS ISL.ANO MERIA (36°38′37″ Ν, 24°53′00″ Ε)
VF42SAS-118/07V. spSKIATHOS-ALONISOS-SKOPELOSSKIATHOS ISL.KATAVOTHRA (39°10′00″ Ν, 23°27′00″ Ε)
VF43RK-104/07V. spRODOS-KASTELORIZO-RORODOS ISLASKLIPIO (36°4′19″ Ν, 27°55′46″ Ε)
VF44XKA-082/07V. fabaCHALKIDIKICHALKIDIKISYKIA (40°2′20″ Ν, 23°56′27″ Ε)
VF46XKA-044/07V. fabaCHALKIDIKICHALKIDIKIMARATHOUSA (38°57′08″ Ν, 23°36′30″ Ε)
VF47RK-054/07V. fabaRODOS-KASTELORIZO-RORODOS ISLAG. ISIDOROS (36°09′55″ Ν, 27°50′59″ Ε)
VF48XKA-093/07V. fabaCHALKIDIKICHALKIDIKIAG. PARASKEVI (38°57′08″ Ν, 23°36′30″ Ε)
VF50SAS-009/07V. fabaSKIATHOS-ALONISOS-SKOPELOSALONNISOS ISL.PATITIRI (39°8′41″ Ν, 23°51′49″ Ε)
VF51XKA-091/07V. spCHALKIDIKICHALKIDIKIPALIOURI (39°56′52″ Ν, 23°39′53″ Ε)
VF53T-096/06V. spTRIKALA-KARDITSAKARDITSANEO IKONIO (39°16′28″ Ν, 22°13′10″ Ε)
VF55T-458/06V. spTRIKALA-KARDITSATRIKALAGAVROS (39°48′02.5″ Ν, 21°35′54.6″ Ε)
VF56P-078/06V. fabaPELOPONESE (ARKADIA-LAKONIA-ILIA-MESSINIA)MESSINIAΚAΚAΝA (37°18′0″ Ν, 21°44′27″ Ε)
VF57IK-095/06V. fabaITHAKA-KEFALONIAITHAKI ISL.PERACHORI (38°20′00″ Ν, 20°43′00″ Ε)
VF58P-091/06V. fabaPELOPONESE (ARKADIA-LAKONIA-ILIA-MESSINIA)MESSINIAVASILIKO (37°15′52″ Ν, 21°53′48″ Ε)
VF60T-336/06V. fabaTRIKALA-KARDITSATRIKALALIGARIA (39°30′38″ Ν, 21°42′3″ Ε)
VF61T-450/06V. fabaTRIKALA-KARDITSATRIKALASKEPARI (39°47′43″ Ν, 21°37′5″ Ε)
VF62IK-124/06V. fabaITHAKA-KEFALONIAKEFALONIA ISL.LOURDATA (38°6′58″ Ν, 20°38′6″ Ε)
VF63X-017/06V. fabaCHIOS-LEMNOSCHIOS ISL.XALKEION (38°20′0.06″ Ν, 26°5′54.42″ Ε)
VF65SK-071/06V. fabaSERRES-KILKISSERRESKATO POTAMIA KILKIS (40°57′24″ Ν, 22°56′04″ Ε)
VF66X-015/06V. fabaCHIOS-LEMNOSCHIOS ISL.XALKEION (38°20′0.06″ Ν, 26°5′54.42″ Ε)
VF67K-044/06V. fabaKOZANI-GREVENAKOZANISPILIA (40°35′26″ Ν, 21°46′36″ Ε)
VF69K-034/06V. fabaKOZANI-GREVENAKOZANIPTOLEMAIDA (40°30′52.99″ Ν, 21°40′43.00″ Ε)
VF70M-122/06V. fabaMITILINI (LESVOS)LESVOS ISL.AG. PARASKEVI (39°14′53″ Ν, 26°16′17″ Ε)
VF71P-215/06V. fabaPELOPONESE (ARKADIA-LAKONIA-ILIA-MESSINIA)ARKADIAARTEMISIO (37°40′34″ Ν, 22°22′40″ Ε)
VF72K-150/06V. fabaKOZANI-GREVENAGREBENAPONTINH (40°4′16″ Ν, 21° 40′36″ Ε)
VF74K-230/06V. fabaKOZANI-GREVENAGREBENAPALIOURIA (39°57′07″ Ν, 21°43′12″ Ε)
VF76IK-026/06V. fabaITHAKA-KEFALONIALEFKADA ISL.KARUA (38°45′0″ Ν, 20°39′0″ Ε)
VF77IRENAV. fabaCultivar
VF78DIVAV. fabaCultivar
VF79MELODIEV. fabaCultivar
VF80FAVELV. fabaCultivar
Table 2. Name and sequence of the SCoT molecular markers used in genetic diversity of faba bean accessions.
Table 2. Name and sequence of the SCoT molecular markers used in genetic diversity of faba bean accessions.
PrimersSequence 5′→3′Annealing TemperatureSize Range (bp)
SCoT 13ACGACATGGCGACCATCG502000–200
SCoT 14ACGACATGGCGACCACGC503000–100
SCoT 15ACGACATGGCGACCGCGA503000–100
SCoT 61CAACAATGGCTACCACCG503000–300
SCoT 66ACCATGGTACCAGCGAG504000–300
ScoT 33CCATGGCTACCACCGCAG503000–200
Table 3. Genetic diversity of fifty-three V. faba populations.
Table 3. Genetic diversity of fifty-three V. faba populations.
PopulationN 1NPB 2No. of Private BandsP (%) 3Shannon Index (I) 4GD 5
VF9555029.820.1700.115
VF11564034.210.1970.133
VF28554023.680.1380.094
VF30551029.820.1810.125
VF57554127.190.1590.108
VF72551022.810.1350.093
VF35550028.070.1590.107
VF44554034.210.1970.133
VF74545023.680.1370.093
VF77553028.070.1660.114
VF78544021.050.1270.087
VF80554036.840.2070.139
VF58554027.190.1540.104
VF62558039.470.2350.161
VF65561037.720.2110.142
VF79549028.950.1660.112
VF17570042.110.2500.171
VF31568045.610.2600.175
VF39563035.090.1980.133
VF41561035.090.2090.143
VF38566042.980.2540.174
VF27559038.600.2160.145
VF71562042.980.2560.175
VF51555041.230.2310.154
VF36568247.370.2660.178
VF46567038.600.2250.153
VF15559039.470.2320.159
VF40560043.860.2540.173
VF50563048.250.2810.191
VF1560036.840.2120.143
VF25563043.860.2650.182
VF29563036.840.2040.136
VF33567038.600.2290.157
VF42566043.860.2590.177
VF43560034.210.1940.131
VF47567045.610.2750.189
VF48563143.860.2410.160
VF55557034.210.2060.142
VF56561034.210.2000.136
VF66559042.980.2380.159
VF67558038.600.2250.153
VF76562046.490.2780.191
VF2564039.470.2380.164
VF6564041.230.2400.163
VF37567057.020.3280.222
VF70576048.250.2730.184
VF69575044.740.2570.174
VF61575036.840.2060.138
VF60561035.960.2030.136
VF63572046.490.2780.191
VF4562035.090.2040.139
VF3568046.490.2750.188
VF53562045.610.2650.180
MO5 37.950.2200.149 1
1 N: number of individuals; 2 NPB: number of polymorphic bands; 3 P (%): percentage of polymorphic bands; 4 GD: gene diversity; I: 5 Shannon’s information index.
Table 4. Analysis of molecular variance results for fifty-three faba bean populations based on SCoT molecular markers.
Table 4. Analysis of molecular variance results for fifty-three faba bean populations based on SCoT molecular markers.
SourcedfSSMSEst. Var.%
Among Pops522722,90652,364834444%
Within Pops2122256,80010,64510,64556%
Total2644979,706 18,989100%
StatValueP(rand ≥ data)
PhiPT0.4390.010
Table 5. Pairwise genetic distances (below diagonal) between the accessions of faba bean.
Table 5. Pairwise genetic distances (below diagonal) between the accessions of faba bean.
PopulationVF9VF11VF28VF30-VF63VF4VF3VF53
VF90.000
VF110.1390.000
VF280.2210.1380.000
VF300.2010.1170.1610.000
-----0.000----
VF630.4160.3920.4190.370-0.000
VF40.5670.5210.6070.545-0.3340.000
VF30.4610.4140.4810.411-0.2180.2180.000
VF530.4180.4080.4580.411-0.1970.3260.0490.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Avramidou, E.; Ganopoulos, I.; Mylona, P.; Abraham, E.M.; Nianiou-Obeidat, I.; Osathanunkul, M.; Madesis, P. Comparative Analysis of the Genetic Diversity of Faba Bean (Vicia faba L.). Sustainability 2023, 15, 1016. https://0-doi-org.brum.beds.ac.uk/10.3390/su15021016

AMA Style

Avramidou E, Ganopoulos I, Mylona P, Abraham EM, Nianiou-Obeidat I, Osathanunkul M, Madesis P. Comparative Analysis of the Genetic Diversity of Faba Bean (Vicia faba L.). Sustainability. 2023; 15(2):1016. https://0-doi-org.brum.beds.ac.uk/10.3390/su15021016

Chicago/Turabian Style

Avramidou, Eleni, Ioannis Ganopoulos, Photini Mylona, Eleni M. Abraham, Irini Nianiou-Obeidat, Maslin Osathanunkul, and Panagiotis Madesis. 2023. "Comparative Analysis of the Genetic Diversity of Faba Bean (Vicia faba L.)" Sustainability 15, no. 2: 1016. https://0-doi-org.brum.beds.ac.uk/10.3390/su15021016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop