Next Article in Journal
Modulation of Sex Pheromone Discrimination by a UDP-Glycosyltransferase in Drosophila melanogaster
Next Article in Special Issue
Molecular Analysis of the Official Algerian Olive Collection Highlighted a Hotspot of Biodiversity in the Central Mediterranean Basin
Previous Article in Journal
Comprehensive Sequence Analysis of IQD Gene Family and their Expression Profiling in Grapevine (Vitis vinifera)
Previous Article in Special Issue
Characteristics of Microsatellites Mined from Transcriptome Data and the Development of Novel Markers in Paeonia lactiflora
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Serendipitous In Situ Conservation of Faba Bean Landraces in Tunisia: A Case Study

1
National Gene Bank of Tunisia (BNG), Street Yesser Arafet, Charguia 1, Tunis 1080, Tunisia
2
Institut National de la Recherche Agronomique de Tunisie (INRAT), Agricultural Applied Biotechnology Laboratory (LR16INRAT06), University of Carthage, Rue Hédi Karray, Menzah 1, Tunis PC 1004, Tunisia
3
Institut National de la Recherche Agronomique de Tunisie (INRAT), Field Crop Laboratory (LR16INRAT02), University of Carthage, Rue Hédi Karray, Menzah 1, Tunis PC 1004, Tunisia
4
Institute of Biosciences and Bioresources of the National Research Council (IBBR-CNR), Via Amendola 165/A, 70126 Bari, Italy
5
Department of Soil, Plant and Food Sciences (DiSSPA), Sect. Genetics and Plant Breeding, University of Bari Aldo Moro, Via Amendola 165/A, 70126 Bari, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally.
Submission received: 21 January 2020 / Revised: 13 February 2020 / Accepted: 16 February 2020 / Published: 24 February 2020
(This article belongs to the Special Issue Genetic Diversity Assessment and Marker-Assisted Selection in Crops)

Abstract

:
Cultivation of faba bean (Vicia faba L.) in Tunisia is largely based on improved varieties of the crop. However, a few farmers continue to produce local cultivars or landraces. The National Gene Bank of Tunisia (NGBT) recently launched a collection project for faba bean landraces, with special focus on the regions of the North West, traditionally devoted to cultivating grain legumes, and where around 80% of the total national faba bean cultivation area is located. The seed phenotypic features of the collected samples were studied, and the genetic diversity and population structure analyzed using simple sequence repeat markers. The genetic constitution of the present samples was compared to that of faba bean samples collected by teams of the International Center for Agricultural Research in the Dry Areas (ICARDA) in the 1970s in the same region, and stored at the ICARDA gene bank. The results of the diversity analysis demonstrate that the recently collected samples and those stored at ICARDA largely overlap, thus demonstrating that over the past 50 years, little genetic change has occurred to the local faba bean populations examined. These findings suggest that farmers serendipitously applied international best practices for in situ conservation of agricultural crops.

1. Introduction

Faba bean (Vicia faba L., 2n = 2x = 12) is a facultative cross-pollinating species with outcrossing rates varying between 1 and 55% depending on its environment; it belongs to the Fabaceae family, Faboideae subfamily, tribe of Fabeae, and is not interfertile with any other Vicia species [1]. The wild progenitor of V. faba is unknown, but recent archaeological excavations have allowed, in the Mont-Carmel (Mediterranean Levantine), the discovery of fossilized seeds that are compatible with a wild progenitor of this crop, dating as back as 14,000 ybp [2]. Considering other archaeological evidences, such as those relating to findings in Tel el Kerkh [3], northwest Syria, it is possible to hypothesize that this species has been domesticated since the Neolithic era, and that the wild progenitor, possibly distributed in small habitats, was entirely domesticated and then became extinct [2,3,4,5]. According to Cole [6] and Cubero [7], the spread of faba bean from its center of origin to other countries could have involved five routes. In the Mediterranean, in particular, faba bean mainly spread through two routes: the first across Anatolia to Greece, the Illyric coast (possibly the Danubian regions), and then to Italy; the second, beginning at the Nile Delta, moving towards the West, along the North African Mediterranean coast, to the Maghreb and then to the Iberian Peninsula. It is worth mentioning that, in this regard, North Africa and Tunisia in particular constitute a center of primary and secondary diversification of several agricultural and wild species [8].
In Tunisia, faba bean covers more than 70% (59,583 ha) of the total area annually devoted to grain legume crops [9]. Approximately half of the area is cultivated with grain legume types meant for fresh pod consumption; the rest as forage (plant and seeds) is mainly based on small seeded types. The average productivity in Tunisia is 1.03 t/ha, 40% below the world average [10]. This is mainly due to the parasitic weed broomrape (Orobanche crenata Forssk. and O. foetida Poir.) and drought stress occurring in faba bean-growing areas [11]. Until the last century, most crops consisted of landraces, often named after the farmer who selected them or after areas where they were grown [12]. Some landraces (‘Batata’, ‘Malti’, ‘Chemchali’ or ‘Masri local’) are still grown by farmers and seeds can be bought from local informal markets.
In recent years, thanks to significant achievements in crop breeding, modern high-yielding varieties are widely used in cultivation and have almost completely replaced local populations and landraces [11]. The increase in yield was obtained mainly with breeding programs targeted at tolerance to abiotic (heat and drought) and biotic (foliar diseases and parasitic weeds) stresses. Moreover, recent breeding efforts are directed towards the development of new cultivars with low anti-nutritional compounds (vicine and convicine), to improve the quality and utilization efficiency of faba in human diet and for livestock feed [13]. Different molecular tools were used to investigate genetic diversity in grain legume species [14,15,16,17]. Some attempts to evaluate genetic variability have also been made for Tunisian faba bean germplasm. The analysis of isozyme [18] and sequence specific amplified polymorphism (SSAP) [19] markers analyzed in nine Tunisian Vicia faba collections has indicated a certain degree of genetic cohesiveness. Using simple sequence repeat (SSR) markers, 16 faba bean accessions, selected from 42 populations collected across eight southern oases arid agro-ecosystems, were analyzed, evidencing genetic cohesiveness among the studied samples [19,20,21], together with a low level of variability among accessions. In both reports, the authors stated that intense seed exchange among farmers had led to a leveled degree of genetic diversity among those populations.
In the 1960s, major concerns focused on the genetic erosion of biodiversity, eventually leading to fostering of ex situ conservation efforts and the creation of gene banks [22]. As a result, in the latest decades of the 20th century, several different research centers organized collection missions for crop diversity in order to secure the local germplasm before it was completely lost [23]. Within the frame of “emergency” collections, from the seventies until the early nineties, ICARDA, among others, carried on a series of faba bean collection missions in North Africa. During that period, several faba bean samples were retrieved from different regions, in particular from Tunisia. At the end of the 20th century, taking into account the evident loss of genetic diversity in the Mediterranean [24], specific attention was paid to the practice of in situ conservation of crops. This is the conservation of agricultural genetic resources on farms located in the same areas where local communities had developed them, with specific attention to neglected crops [25,26]. According to Duc et al. [23], in situ conservation of biodiversity may contribute to the development of the best-adapted materials for local agronomic practices and involve farmers in the selection process through participatory breeding [27]. Within this frame, NGBT started an ongoing program of genetic resource collection in different areas of Tunisia. The aim was to preserve Tunisian crop gene-pools from genetic erosion and characterize local germplasm, thanks to an integrated approach, including on-farm conservation of local germplasm and landraces.
In order to better plan an in situ conservation strategy for faba bean and to understand the possible loss of genetic diversity in Tunisian Vicia faba germplasm, the genetic structure of the samples collected in recent years was compared with those collected by ICARDA in the 1970s. This paper reports on the results of this comparison.

2. Materials and Methods

2.1. Plant Material

The plant material used in the present study consisted of a collection of 51 Tunisian local faba bean samples (Table 1). It included 29 samples collected during 2016–2018 by the NGBT and 22 faba bean accessions collected by ICARDA, starting from the seventies until the early nineties. The NGBT samples were collected in the governorates of Beja and Jendouba (Figure S1) characterized by annual average rainfall of 800 and 600 mm, respectively. The passport data and ethno-botanical information of the NGBT samples are available at NGBT. The samples of ICARDA (labelled ICAR) were derived from collection missions conducted from the seventies until the early nineties in North Tunisia (Beja, Bizerte, and Siliana).

2.2. Seed Phenotypic Traits

Five seeds of each sample were randomly selected to determine average seed size. The three axial dimensions of seed length (L), width (W), and thickness (T) were measured using a Vernier caliper (Gilson Tools, Japan) with accuracy of 0.05 mm. The geometric mean diameter (Dg) was calculated by using the equation reported by Mohsenin [28]:
Dg = (L × W × T)1/3
The sphericity (φ) of faba bean seeds was calculated using the following formula:
φ = [(L × W × T)1/3/L]*100

2.3. DNA Extraction and SSR Assays

Total genomic DNA was extracted from fresh young leaves—five plants per sample—using the cetyltrimethyl ammonium bromide (CTAB) method, as described by Fulton et al. [29]. DNA concentration was determined using a NanoDropTM ND-2000 (Thermo Scientific, MA, USA) and the quality was verified by separation on 0.8% agarose gel. Equal DNA quantities of the five plants of the sample were then pooled, and all DNA samples were diluted to a standard working concentration of 50 ng/µl by adding ultrapure water (Gibco, Invitrogen, USA).
A set of 11 simple sequence repeats (SSRs) markers, retrieved from the literature [30,31,32] were used for this study (Table S1). A preliminary assay was carried out in order to evaluate the robustness of PCR reaction and the reproducibility of the fragments. In particular, different faba samples randomly chosen were analysed considering technical replicates. The PCR conditions for each SSR markers were set up at best conditions, considering an annealing temperature ranging from 45 to 60 °C. The fragments produced were separated on 2.0% agarose gel containing Nancy-520 DNA Gel Stain (Sigma Genosys, St. Louis, MO, USA), and visualized under UV light. The amplification reactions were performed in a final volume of 20 µl, containing the template DNA (50 ng), the 5′end of each forward primer with the M13 (21 bp) tail, the reverse primer (M13) labeled with fluorescent dye (FAM, VIC, PET, or NED). PCR reactions were performed in a thermal cycler (Bio-Rad Laboratories, Hercules, CA, USA) as follows: an initial denaturing step at 95 °C for 3 min, followed by 36 cycles of 94 °C for 20 s, 56 °C for 50 s, 72°C for 1 min, and a final extension step at 72 °C for 7 min. The amplification products were detected by automatic capillary sequencer ABI PRISM 3100 Genetic Analyzer (Applied Biosystems, Waltham, MA, USA), and the fragments were analyzed with GeneMapper genotyping software version 5.0 (Thermo Fisher Scientific, Waltham, MA, USA). The internal molecular weight standard was GeneScanTM 600 LIZ dye Size Standard (Thermo Fisher Scientific, Waltham, MA, USA).

2.4. Data Analysis

Hierarchical ascending classification (HAC) clustering analysis based on dissimilarity matrix of morphometric seed data (L, W, T, Dg, and φ) was performed to evaluate the relationship among the faba samples using XLSTAT statistical software ver. 2016.2 (Addinsoft Inc, New York, USA).
The genetic indices, number of different alleles (Na), Shannon’s information index (I), observed Heterozygosity (Ho), expected Heterozygosity (He), Fixation Index (F), and private alleles were calculated using GenAlEx version 6.5 [33]. The allelic data were used to obtain a similarity matrix, from which a dendrogram was constructed using the UPGMA algorithm with MEGA ver. 4 [34]. The molecular data were processed using STRUCTURE ver. 2.3.4 [35]. The number of sub-populations (K) was estimated by 10 independent runs for each K (from 1 to 10), applying the admixture model, 500000 Markov Chain Monte Carlo (MCMC) repetitions, and a 100000 burn-in period. The means of the log-likelihood estimates for each K were calculated. The true K was determined with the Evanno test [36] using STRUCTURE HARVESTER [37]. Analysis of molecular variance (AMOVA) was used to partition the genetic variation into inter- and intra-gene pool diversities in faba using GenAlEx program ver. 6.5, with 1000 permutations.

3. Results

3.1. Variation in Seeds Phenotypic Traits

Morphometric seed traits (L, W, T, Dg, and φ) were measured in the samples belonging to both the NGBT and ICARDA faba bean collection. The average mean of the three principal axial dimensions (L, W, and T) and Dg of the NGBT and ICAR groups are shown in Figure 1. No statistically significant differences for these morphometric seed traits were detected between the two groups. The values of faba bean sphericity (φ) were calculated by using the geometric mean diameter and length data. No significant difference was found for these traits when comparing the NGBT and ICAR groups.
The values of morphometric seed traits of each samples were used to calculate a dissimilarity matrix based on Euclidean distances. A dendrogram was obtained starting from the matrix, in which three main clusters are evidenced, corresponding to three seed types—large, medium (called also equina), and small (Figure 2). Each cluster included both NGBT and ICAR samples with no reference to the area of origin.

3.2. Molecular Variation of Faba Bean Collection

In order to evaluate the genetic diversity of 51 faba bean samples, a set of 11 SSR markers was used (Table 1 and Table S1). A total number of 94 alleles were identified, ranging from 3 (loci M22 and M46) to 22 (locus VFG41), with an average of 10.2 alleles per locus (Table S2).
AMOVA did not show a molecular diversity among the 29 NGBT faba beans when these were grouped according to the collection sites (El Hamra, Oued Ghrib and Fouazia), suggesting that the NGBT faba bean samples belonged to one genetically cohesive population (Table S3).
According to these results, we investigated the genetic diversity among samples collected in recent years by NGBT and ICAR during the 1970s, in the same Tunisian areas. Similar values for the number of different alleles, Shannon’s information index, Heterozygosity expected, Heterozygosity observed, Diversity index, and Fixation index were found between the two groups, although the ICAR samples appeared slightly more fixed (Table 2).
AMOVA indicated that the NGBT and ICAR groups have no statistical difference, at the molecular level (Table 3), as most of the diversity clearly appeared within groups, with only a limited variation occurring among them.
To define the genetic relationships among the faba bean samples, the similarity matrix obtained was also used to produce a UPGMA dendrogram (Figure 3). The clustering showed that the NGBT samples were admixed with ICAR ones, thus suggesting that the two groups could be considered a unique meta-population. Total faba bean collection was also evaluated with Bayesian clustering modeling, performed using SSR allelic data generated according to 11 SSR markers. As the clustering model presumes the underlying existence of K clusters, an Evanno test [36] was performed that yielded K=3 as the highest log-likelihood (Figure S2). Nevertheless, each sample analyzed showed to belong to all three clusters identified, and none of them predominantly pertained to a specific cluster (Figure 4). These results seemingly indicate that the Tunisian faba collection was structured in three subpopulations that do not correspond to NGTB and ICAR samples’ distinction.

4. Discussion

Faba bean is an important crop for sustainable agriculture in marginal areas and advanced agricultural systems, as it plays an important role in soil fertility and nitrogen fixation, being able to grow in diverse climatic and soil conditions. Although faba bean is less consumed in western countries as human food, it is considered one of the main sources of cheap protein and energy for many people in Africa, parts of Asia, and Latin America, where many people cannot afford to buy meat [38]. Seed size and shape are characters of polygenic control [39,40] and have undergone strong selection measures by farmers during evolution of the crop. Seed size is considered a key trait in the study of the historical evolution of this crop based on archaeological remains and findings [2]. This selection pressure can still be found in the habits of farmers, who manually select seeds to be sown in the next season [41], a habit that was followed in many Mediterranean regions until modern times and lost only on introduction of improved varieties. This manual selection procedure resulted in the formation of peculiar landraces, such as “Larga di Leonforte” [42]. However, this study noted that Tunisian farmers did not practice seed selection, and have not done so for the last 50 years.
The analysis of morphometric seed traits and molecular markers carried out in this work did not distinguish patterns in the distribution of morphological and genetic variations. All the samples collected by the NGBT teams appeared to belong to a single, genetically cohesive population. Similar results were observed for the faba bean samples obtained by ICARDA and collected from the same areas. When all faba samples were analyzed together, there were no differences between the NGBT and ICARDA groups. In fact, our molecular analyses and seed morphometric variation study demonstrated that the two groups belonged to a single meta population.
Seed size is one of the most important morphological traits responsible for yield, and a major target for breeding. Several studies have led to the mapping of QTLs/genes for seed weight/size in soybean [43], chickpea [44], and lentil [45]. In faba bean, the genetic control of these traits is still unclear, although consensus linkage maps have been produced [43,46].
Our data confirmed that morphometric seed traits were not associated to the markers used. This implies that any selection of lines based on seed traits still retains a variable level of genetic diversity potentially associated to other traits, such as adaptation or resistance, an issue to be taken into account in any faba bean-breeding program based on phenotypic data.
The absence of genetic differentiation in different collection sites might depend on many factors. The absence of human selection pressure does not force crop adaptation in a specific direction. The outbreeding habit of faba bean is a second factor; in fact, while the advanced breeding varieties favor inbreeding in search of higher stability, as requested by UPOV standards, the landraces are generally quite allogamous [47,48]. Recent studies have also demonstrated that the level of allogamy might depend on the species of pollinating insects [49]. A further, effective mechanism could be derived by the spontaneous seed exchange practice occurring among farmers. In informal seed markets, seed exchange by farmers favors the establishment of a landrace in a given environment with uniform agro-climatic characteristics. The farmers cultivating these landraces are actually the relics of a once wider cultivation area. All these factors act synergically to produce the observed genetic differentiation patterns.
The Evanno test yielded three subpopulations within the faba bean collection. In addition, each sample had a coefficient membership lower than 0.50, thus denying that it predominantly belonged to one of these subpopulations. This further supports the fact that the two meta-populations, NGBT and ICARDA, are subsamples of a unique population collected in the same area at different times.
The biological significance of the three subpopulations detected with STRUCTURE is unclear, and might depend on the presence of genetic signals that do not parallel clear-cut characteristics. Nevertheless, several authors have reported that hierarchical analyses, such as those based on STRUCTURE or similar software, relay on strict assumptions that might not completely apply to the case studies, thus resulting in incorrect evaluation of the population structure [50].
In conclusion, unintended conservation of ancient faba bean germplasm in Tunisian farms is witnessed, because the farmers cultivating faba bean landraces do not follow seed selection, as well as owing to concurring factors such as the high level of cross-pollination, and the consistent presence of pollinators (due to the limited use or absence of insecticides). This beneficial situation is an empirical application of best practices recommended by research institutions for on-farm conservation of plant genetic resources. In their quest to feed their families and their precious animals, the Tunisian farmers in the small villages mentioned above serendipitously stored and protected their faba bean genetic resources. These findings and implications should be further discussed in the broadest context possible. Future research directions should also be highlighted.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/2/236/s1, Figure S1: NGBT sampling Figure S2: Delta K for differing numbers of subpopulations (K) estimated for faba bean collection with 11 markers; Table S1: List of SSR markers used for genetic analysis of faba bean collection; Table S2: Number of different alleles (Na), Shannon’s information index (I), Heterozygosity observed (Ho), Heterozygosity expected (He), and Fixation Index (F) of 11 SSR markers used in faba bean collection; Table S3: Analysis of molecular variance (AMOVA) of 29 NGTB faba bean split according to collection sites (El Hamra, Oued Ghrib, and Fouazia).

Author Contributions

Conceptualization, E.B., K.K., D.P., M.M.F.-S., and G.M.; Formal analysis, M.M.M., D.D., and G.M.; Funding acquisition, C.M. and D.P.; Investigation, E.B., K.K., W.S., and M.M.F.-S.; Methodology, D.P., M.M.F.-S., and G.M.; Resources, E.B. and K.K.; Supervision, D.P., M.M.F.-S., and G.M.; Writing—original draft, E.B., K.K., D.P., M.M.F.-S., and G.M.; Writing—review & editing, D.P., M.M.F.-S., and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out within the framework of the project “Ressources phytogénétiques tunisiennes mieux conservées et valorisées”, funded by the Italian Cooperation to the Tunisian Ministry of Local Affairs and Environment.

Acknowledgments

The authors wish to thank Amine Hmid, Anis Khlij, and Alberto Dragotta from the Mediterranean Agronomic Institute of Bari, Italy (CIHEAM-IAM); Olfa Saddoud and Mbarek Ben Naceur from the National Gene Bank of Tunisia for project management and administration. Special thanks are also given to Mr. Mostafa Khemiri and Mr. Ali Ayari, agronomic engineer-extensionists and head of CTV Fernana and CTV Amdoun, respectively, for facilitating contacts with farmers and help during the collection missions. We give special thanks to Nicoletta Rapanà for excellent technical assistance.

Conflicts of Interest

The authors declare that they have no competing interests for this research.

References

  1. Leht, M.; Jaaska, V. Cladistic and phenetic analysis of relationships in Vicia subgenus Vicia (Fabaceae) by morphology and isozymes. Plant Syst. Evol. 2002, 232, 237–260. [Google Scholar] [CrossRef]
  2. Caracuta, V.; Barzilai, O.; Khalaily, H.; Milevski, I.; Paz, Y.; Vardi, J.; Regev, L.; Boaretto, E. The onset of faba bean farming in the Southern Levant. Sci. Rep. 2015, 5, 14370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tanno, K.; Willcox, G. The origins of Cicer arietinum L. and Vicia faba L.: Early finds from Tell el Kerkh, north-west Syria, late 10th millennium B.P. Veg. Hist. Archaeobot. 2006, 15, 197–204. [Google Scholar]
  4. Caracuta, V.; Weinstein-Evron, M.; Kaufman, D.; Yeshurun, R.; Silvent, J.; Boaretto, E. 14,000-year-old seeds indicate the Levantine origin of the lost progenitor of faba bean. Sci. Rep. 2016, 6, 37399. [Google Scholar] [CrossRef] [Green Version]
  5. Kosterin, O.E. The lost ancestor of the broad bean (Vicia faba L.) and the origin of plant cultivation in the Near East. Vavilov J Genet. PL. BR. 2014, 18, 831–840. [Google Scholar] [CrossRef] [Green Version]
  6. Hole, F. A Reassessment of the Neolithic Revolution. Paléorien 1984, 10, 49–60. [Google Scholar] [CrossRef]
  7. Cubero, J.I. On the evolution of Vicia faba L. Theor. Appl. Genet. 1974, 45, 47–51. [Google Scholar] [CrossRef]
  8. Valderrábano, M.; Gil, T.; Heywood, V.; de Montmollin, B. Conserving Wild Plants in the South and East Mediterranean Region; IUCN: Gland, Switzerland; Málaga, Spain, 2018; p. 146. [Google Scholar]
  9. Anonymous. Direction Générale de la Production Agricole (DGPA); Ministry of Agriculture Fisheries and Hydraulic Resources of Tunisia: El Ain, Tunisia, 2018. [Google Scholar]
  10. 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] [Green Version]
  11. Kharrat, M.; Ouchari, H. Faba bean status and prospects in Tunisia. Grain Legumes 2011, 56, 11–12. [Google Scholar]
  12. Hamza, N. La valorization et l’utilisation des ressources génétiques Locals, moyens de leur conservation durable en Tunisie. Ann. de l’INRAT 2ème Numéro Spécial Centen. 2015, 88, 26–32. [Google Scholar]
  13. Singh, R.K.; Bohra, N.; Sharm, L. Valorizing faba bean for animal feed supplements via biotechnological approach: Opinion. Biocatal. Agric. Biotechnol. 2019, 17, 366–368. [Google Scholar] [CrossRef]
  14. De Giovanni, C.; Pavan, S.; Taranto, F.; Di Rienzo, V.; Miazzi, M.M.; Marcotrigiano, A.R.; Mangini, G.; Montemurro, C.; Ricciardi, L.; Lotti, C. Genetic variation of a global germplasm collection of chickpea (Cicer arietinum L.) including Italian accessions at risk of genetic erosion. Physiol. Mol. Biol. Plants 2016. [Google Scholar] [CrossRef]
  15. Oliveira, H.R.; Tomàs, D.; Silva, M.; Viegas, W.; Veloso, M.M. Genetic diversity and population structure in Vicia faba L. landraces and wild related species assessed by nuclear SSRr. PLoS ONE 2016, 11, e015480. [Google Scholar] [CrossRef]
  16. Liu, Z.; Li, H.; Wen, Z.; Fan, X.; Li, Y.; Guan, R.; Guo, Y.; Wang, S.; Wang, D.; Qiu, L. Comparison of genetic diversity between Chinese and American soybean (Glycine max (L.)) accessions revealed by high-density SNPs. Front. Plant Sci. 2017, 8, 2014. [Google Scholar] [CrossRef] [Green Version]
  17. Pavan, S.; Bardaro, N.; Fanelli, V.; Marcotrigiano, A.R.; Mangini, G.; Taranto, F.; Catalano, D.; Montemurro, C.; De Giovanni, C.; Lotti, C.; et al. Genotyping by sequencing of cultivated lentil (Lens culinaris Medik.) highlights population structure in the Mediterranean gene pool associated with geographic patterns and phenotypic variables. Front. Genet. 2019, 10, 872. [Google Scholar] [CrossRef]
  18. Ouji, A.; Suso, M.J.; Rouaissi, M.; Abdellaoui, R.; El Gazzah, M. Genetic diversity of nine faba bean (Vicia faba L.) populations revealed by isozyme markers. Genes Genom. 2011, 33, 31–38. [Google Scholar] [CrossRef] [Green Version]
  19. Ouji, A.; El Bok, S.; Syed, N.H.; Abdellaoui, R.; Rouaissi, M.; Flavell, A.J.; El Gazzah, M. Genetic diversity of faba bean (Vicia faba L.) populations revealed by sequence specific amplified polymorphism (SSAP) markers. Afr. J. Biotechnol. 2012, 11, 2162–2168. [Google Scholar]
  20. 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. Res. 2014, 278–285. [Google Scholar] [CrossRef]
  21. Yahia, Y.; Guetat, A.; Walid, E.; Ferchichi, A.; Yahia, H.; Loumerem, M. Analysis of agromorphological diversity of southern Tunisia faba bean (Vicia faba L.) germplasm. Afr. J. Biotechnol. 2012, 11, 11913–11924. [Google Scholar]
  22. Dempsey, G.J. In Situ Conservation of Crops and Their Relatives: A Review of Current Status and Prospects for Wheat and Maize; Paper 96-08; CIMMYT NRG: Veracruz, Mexico, 1996. [Google Scholar]
  23. Damania, A.B. History, achievements, and status of genetic resources conservation. Agron. J. 2008, 100, S-27–S-39. [Google Scholar] [CrossRef]
  24. Maxted, N.; Bennett, S.J. Conservation, diversity and use of Mediterranean legumes. In Plant Genetic Resources of Legumes in the Mediterranean; Current Plant Science and Biotechnology in Agriculture; Maxted, N., Bennett, S.J., Eds.; Springer: Dordrecht, The Netherlands, 2001; Volume 39. [Google Scholar]
  25. Altieri, M.A.; Merrick, L. In situ conservation of crop genetic resources through maintenance of traditional farming systems. Econ. Bot. 1987, 41, 41–86. [Google Scholar] [CrossRef]
  26. Khoury, C.; Laliberté, B.; Guarino, L. Trends in ex situ conservation of plant genetic resources: A review of global crop and regional conservation strategies. Genet. Resour. Crop Evol. 2010, 5, 57–625. [Google Scholar] [CrossRef]
  27. Enjalbert, J.; Dawson, J.C.; Paillard, S.; Rhoné, B.; Rousselle, Y.; Thomas, M.; Goldringer, I. Dynamic management of crop diversity: From an experimental approach to on-farm conservation. C. R. Biol. 2011, 334, 458–468. [Google Scholar] [CrossRef]
  28. Mohsenin, N.N. Physical Properties of Plant and Animal Materials; Gordon and Breach Science Publishers: New York, NY, USA, 1970. [Google Scholar]
  29. Fulton, T.; Chunwongse, J.; Tanksley, S. Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Mol. Biol. Rep. 1995, 13, 207–209. [Google Scholar] [CrossRef]
  30. Ma, Y.; Yang, T.; Guan, J.; Wang, S.; Wang, H.; Sun, X.; Zong, X. Development and characterization of 21 EST-derived microsatellite markers in Vicia faba (fava bean). Am. J. Bot. 2011, 98, e22–e24. [Google Scholar] [CrossRef]
  31. Gong, Y.M.; Xu, S.C.; Mao, W.H.; Hu, Q.Z.; Zhang, G.W.; Ding, J.; Li, Z.Y. Generation and characterization of 11 novel est derived microsatellites from Vicia faba (Fabaceae). Am. J. Bot. 2010, 97, 69–71. [Google Scholar] [CrossRef]
  32. Zeid, M.M.; Mitchell, S.; Link, W.; Carter, W.; Nawar, A.; Fulton, T.; Kresovich, S. Simple sequence repeats (SSRs) in faba bean: New loci from Orobanche-resistant cultivar ‘Giza 402’. Plant Breed. 2009, 128, 149–155. [Google Scholar] [CrossRef]
  33. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [Green Version]
  34. Tamura, K.; Dudley, J.; Nei, M.; Kumar, S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol. Biol. Evol. 2007, 24, 1596–1599. [Google Scholar] [CrossRef]
  35. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar]
  36. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Earl, D.A.; von Holdt, B.M. Structure harvester: A website and program for visualizing structure output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]
  38. Alghamdi, S.S.; Migdadi, H.M.; Ammar, M.H.; Paull, J.G.; Siddique, K.H.M. Faba bean genomics: Current status and future prospects. Euphytica 2012, 186, 609–624. [Google Scholar] [CrossRef]
  39. Satovic, Z.; Avila, C.M.; Cruz-Izquierdo, S.; Díaz-Ruíz, R.; García-Ruíz, G.M.; Palomino, C.; Gutiérrez, N.; Vitale, S.; Ocaña-Moral, S.; Gutiérrez, M.V.; et al. A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia faba L.). BMC Genom. 2013, 14, 932. [Google Scholar] [CrossRef] [Green Version]
  40. 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 crops Res. 2010, 115, 243–252. [Google Scholar] [CrossRef]
  41. Negri, V. Landraces in central Italy: Where and why they are conserved and perspectives for their on-farm conservation. Genet. Resour. Crop Evol. 2003, 50, 871–885. [Google Scholar] [CrossRef]
  42. Gresta, F.; Avola, G.; Albertini, E.; Raggi, L.; Abbate, V. A study of variability in the Sicilian faba bean landrace ‘Larga di Leonforte’. Genet. Resour. Crop Evol. 2010, 57, 523–531. [Google Scholar] [CrossRef]
  43. Sun, Y.N.; Pan, J.B.; Shi, X.L.; Du, X.Y.; Wu, Q.; Qi, Z.M.; Jiang, H.W.; Xin, D.W.; Liu, C.Y.; Hu, G.H.; et al. Multi-environment mapping and meta-analysis of 100-seed weight in soybean. Mol. Biol. Rep. 2012, 39, 9435–9443. [Google Scholar] [CrossRef]
  44. Verma, S.; Gupta, S.; Bandhiwal, N.; Kumar, T.; Bharadwaj, C.; Bhatia, S. High-density linkage map construction and mapping of seed trait QTLs in chickpea (Cicer arietinum L.) using Genotyping-by-Sequencing (GBS). Sci. Rep. 2015, 5, 17512. [Google Scholar] [CrossRef] [Green Version]
  45. Verma, P.; Goyal, R.; Chahota, R.K.; Sharma, T.R.; Abdin, M.Z.; Bhatia, S. Construction of a genetic linkage map and identification of QTLs for seed weight and seed size traits in lentil, Lens culinaris Medik. PLoS ONE 2015, 10, e0139666. [Google Scholar] [CrossRef] [Green Version]
  46. Webb, A.; Cottage, A.; Wood, T. A SNP-based consensus genetic map for synteny based trait targeting in faba bean (Vicia faba L.). Plant Biotechnol. J. 2016, 14, 177–185. [Google Scholar]
  47. Suso, M.J.; Pierre, J.; Moreno, M.T.; Esnault, R.; Le Guen, J. Variation in outcrossing levels in faba bean cultivars: Role of ecological factors. J. Agric. Sci. 2001, 136, 399–405. [Google Scholar] [CrossRef]
  48. Link, W.; Edered, W.; Metz, P.; Buiel, H.; Melchinger, A.E. Genotypic and environmental variation for degree of cross-fertilization in faba bean. Crop Sci. 1994, 34, 960–964. [Google Scholar] [CrossRef]
  49. Marzinzig, B.; Brünjes, L.; Biagioni, S.; Behling, H.; Link, W.; Westphal, C. Bee pollinators of faba bean (Vicia faba L.) differ in their foraging behaviour and pollination efficiency. Agric. Ecosyst. Environ. 2018, 264, 24–33. [Google Scholar] [CrossRef]
  50. Lawson, D.J.; Van Dorp, L.; Falush, D. A tutorial on how not to over-interpret STRUCTURE and ADMIXTURE bar plots. Nat. Commun. 2018, 9, 3258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Length (L), width (W), thickness (T), and geometric mean diameter (Dg) averages of samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR). Lines represent the standard deviation.
Figure 1. Length (L), width (W), thickness (T), and geometric mean diameter (Dg) averages of samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR). Lines represent the standard deviation.
Genes 11 00236 g001
Figure 2. Dendrogram based on dissimilarity matrix calculated from morphometric seed traits in the faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR). The colors orange, purple, and green were used to distinguish between the medium, large, and small seed type clusters, respectively.
Figure 2. Dendrogram based on dissimilarity matrix calculated from morphometric seed traits in the faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR). The colors orange, purple, and green were used to distinguish between the medium, large, and small seed type clusters, respectively.
Genes 11 00236 g002
Figure 3. Dendrogram of faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR) resulting from the UPGMA cluster analysis based on similarity matrix obtained from 11 SSR allelic data.
Figure 3. Dendrogram of faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICAR) resulting from the UPGMA cluster analysis based on similarity matrix obtained from 11 SSR allelic data.
Genes 11 00236 g003
Figure 4. Membership coefficient (Q) mean of the faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA). The different colors indicate the three subpopulations detected using a Bayesian approach (blue: subpopulation 1; red: subpopulation 2, and green: subpopulation 3).
Figure 4. Membership coefficient (Q) mean of the faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA). The different colors indicate the three subpopulations detected using a Bayesian approach (blue: subpopulation 1; red: subpopulation 2, and green: subpopulation 3).
Genes 11 00236 g004
Table 1. List of faba bean samples included in the study.
Table 1. List of faba bean samples included in the study.
Id
Name
Local
Name
GovernorateLocationLongitude
(E)
Latitude
(N)
Seed
Type
NGBT 1MaltiBejaEl Hamra9.01164136.52219Large
NGBT 3Malti9.01164136.52219Large
NGBT 4Malti9.01164136.52219Large
NGBT 5Chemchali9.01164136.52219Small
NGBT 8Chemchali9.01164136.52219Small
NGBT 9 9.01250236.52136Medium
NGBT 10-9.01250236.52136Small
NGBT 13-9.01250236.52136Small
NGBT 34-9.01250236.52136Small
NGBT 35-9.01250236.52136Small
NGBT 50Chemchali9.01164136.52219Small
NGBT 51Chemchali9.01164136.52219Small
NGBT 52Chemchali9.01164136.52219Small
NGBT 53Chemchali9.01164136.52219Small
NGBT 55-9.01250236.52136Medium
NGBT 66-9.01250236.52136Medium
NGBT 16-JendoubaOued Ghrib8.41281536.37286Small
NGBT 18-8.41281536.37286Large
NGBT 21Bachar8.41420336.37533Small
NGBT 56-8.41420336.37533Large
NGBT 57-8.41420336.37533Large
NGBT 60-8.41281536.37286Small
NGBT 61-8.41281536.37286Large
NGBT 63-8.41420336.37533Large
NGBT 36MaltiFouazia8.40481136.40103Large
NGBT 38Malti8.40481136.40103Large
NGBT 48Malti8.40481136.40103Large
NGBT 62Malti8.40481136.40103Large
NGBT 64Malti8.40481136.40103Large
ICAR 22Localn.a.n.a.n.a.n.a.Small
ICAR 23Sevillen.a.n.a.n.a.n.a.Large
ICAR 24Misri 32Silianan.a.9.61667036.35000Small
ICAR 25Localn.a.n.a.n.a.n.a.Large
ICAR 26-Bizerten.a.n.a.n.a.Large
ICAR 27-Bejan.a.n.a.n.a.Medium
ICAR 28-n.a.n.a.n.a.Large
ICAR 29-n.a.n.a.n.a.n.a.Large
ICAR 30LocalBizerten.a.n.a.n.a.Small
ICAR 31Localn.a.n.a.n.a.n.a.Large
ICAR 32Localn.a.n.a.n.a.n.a.Small
ICAR 33-n.a.n.a.n.a.n.a.Small
ICAR 39LocalBizerten.a.n.a.n.a.Medium
ICAR 40Malti 24n.a.9.66667037.05000Medium
ICAR 41Malti 25n.a.9.66668037.06000Small
ICAR 42Misri 39Bejan.a.9.58333036.66670Small
ICAR 43Misri 41n.a.9.21667036.71670Small
ICAR 44-n.a.n.a.n.a.Medium
ICAR 45-n.a.n.a.n.a.Small
ICAR 46LocalBizerten.a.n.a.n.a.Large
ICAR 47Localn.a.n.a.n.a.n.a.Small
ICAR 68Localn.a.n.a.n.a.n.a.Large
n.a. data not available.
Table 2. Number of different alleles (Na), Shannon’s information index (I), Heterozygosity observed (Ho), Heterozygosity expected (He), Fixation Index (F), and private alleles of faba bean collection split in samples recently collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA).
Table 2. Number of different alleles (Na), Shannon’s information index (I), Heterozygosity observed (Ho), Heterozygosity expected (He), Fixation Index (F), and private alleles of faba bean collection split in samples recently collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA).

Samples
NaIHoHeFPrivate
Alleles
NGBT2971.4280.6210.6580.02723
ICARDA2261.3920.5460.6570.14512
Whole Collection5171.4100.5840.6570.086
Table 3. Analysis of molecular variance (AMOVA) of faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA).
Table 3. Analysis of molecular variance (AMOVA) of faba bean collection split in samples collected by the National Gene Bank of Tunisia (NGBT) and Agricultural Research in the Dry Areas (ICARDA).
Source of VariationdfSSMSEst. Var.p-Values
Among groups13.6503.650.01n.s
Within groups49161.943.313.3
Total50165.59 3.31
df = degree of freedom, SS = sum of squares, MS = mean squares, Est. Var. = estimate of variance, p-value.

Share and Cite

MDPI and ACS Style

Babay, E.; Khamassi, K.; Sabetta, W.; Miazzi, M.M.; Montemurro, C.; Pignone, D.; Danzi, D.; Finetti-Sialer, M.M.; Mangini, G. Serendipitous In Situ Conservation of Faba Bean Landraces in Tunisia: A Case Study. Genes 2020, 11, 236. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11020236

AMA Style

Babay E, Khamassi K, Sabetta W, Miazzi MM, Montemurro C, Pignone D, Danzi D, Finetti-Sialer MM, Mangini G. Serendipitous In Situ Conservation of Faba Bean Landraces in Tunisia: A Case Study. Genes. 2020; 11(2):236. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11020236

Chicago/Turabian Style

Babay, Elyes, Khalil Khamassi, Wilma Sabetta, Monica Marilena Miazzi, Cinzia Montemurro, Domenico Pignone, Donatella Danzi, Mariella Matilde Finetti-Sialer, and Giacomo Mangini. 2020. "Serendipitous In Situ Conservation of Faba Bean Landraces in Tunisia: A Case Study" Genes 11, no. 2: 236. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11020236

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