Abstract

Faba bean (Vicia faba L.) is a legume crop cultivated for its nutritious seeds that are an important worldwide source of human food and feed. Seeds characterization is a prerequisite step for faba bean quality improvement. The morphophysical characterization of the seeds of twenty-four local faba bean accessions following the UPOV descriptors and the AOAC International standards was carried out and assessed with an approach based on Euclidean statistical model. “205 Bulk” was the unique accession harboring white hilum color seed which is linked to low convicine grain content. Irregular seed shape was the most observed character among the studied accessions except “Badii” displaying an elliptical seed shape; therefore, seed shape did not allow discrimination within our Tunisian germplasm. Interestingly, the physical characters of the seeds showed significant diversity between the accessions for all the measured parameters. A highly significant variability was observed for axial, length, and width dimensions of seeds, with “Memdouh” being the longest and largest seed accession, whereas “01-02” was the shortest and narrowest. Classification of the studied faba bean germplasm accessions based on morphophysical characters using clustering by Euclidean distance revealed three different groups. Moreover, multivariate PCA analysis further classified the faba bean accessions into four main clusters. Correlation study performed by using Spearman’s test established positive correlations within physical parameters of seeds such as between mean length and mean width of seeds. Therefore, using morphophysical parameters screening, valuable phenotypes have been selected for deeper physiological characterization and further breeding programs.

1. Introduction

Faba bean, Vicia faba L., a diploid (2n = 12) Leguminosae species, also known as broad bean, is the third most important food grain legume after soybean and pea [1, 2]. It is mainly cultivated all over the world for its rich protein content, and therefore it serves for both food and feed use [3]. Moreover, introduced in intercropping system with cereals, it contributes to nitrogen fixation from soil, soil fertility improvement, and thus crop yield increase [1]. Faba bean is cultivated in almost seventy countries, on 2.58 million ha with a yield of 5.43 million tons [4]. In Tunisia, faba bean with its two seed sizes (large = broad bean or small = field bean) is the most cultivated Fabaceae with 60 thousand hectares representing 70% of the total grown grain legume area [5]. Furthermore, in 2018, the Tunisian Ministry of Agriculture accorded a budget to buy 150.000 tons of field bean from farmers to re-wholesale them later to fodder manufactures. This is a first step of the five-year plan national strategy to promote faba bean and increase the acreage up to double. Indeed, this approach is focusing on encouraging farmers to introduce faba bean into the cropping system as this crop plays a key role in agriculture sustainability. However, this strategy needs other alternatives to develop a value chain for this crop, which could create a demand market for highlighting the flour of faba bean in food industry. Consequently, a focus on physical properties and seed quality traits is the key for that purpose. Moreover, diversity study through the analysis of the physical properties of seeds is a preliminary step for faba bean germplasm characterization and improvement. Several attempts to analyze the morphological diversity and to characterize the physical properties of faba bean seeds have been realized, whether concerning varieties, accessions, cultivars, line, or even mutants [613]. However, no one paid great attention to the Tunisian germplasm that is rich in autochthonous accessions scattered among farmers. Therefore, in our present study, we initiated an analysis of the genetic variability of Tunisian faba bean germplasm based on seeds physical parameters. Our work aims ultimately to highlight the technological value of faba bean grains of Tunisian origin and to valorize their nutraceutical potentialities. The results of our study were further validated through powerful statistical assessment using Euclidean clustering, PCA multivariate analysis, and Spearman’s correlation to verify the significance and meaningfulness of each seed character.

2. Materials and Methods

2.1. Plant Material

Twenty-four Tunisian accessions of Vicia faba L. were used in this study. Seeds were obtained from three different origins: local populations and old landraces (16) provided by farmers and conserved “ex situ” in Tunisian National gene bank; Tunisian germplasm (4) collected during the seventies and preserved in the German Gene Bank (Institute of Plant Genetics and Crop Research); and improved varieties (4) as reference.

2.1.1. Local Populations

Part of the genetic material used in this work was collected from farmers’ own seed stock from different localities in northwestern of Tunisia. Indeed, only two farmers (Farm 1 and Farm 2) who are still growing their old own local landraces were found in “Amdoun” region from “Béja” governorate with 400 m altitude and average annual rainfall of 600 mm. Seeds were collected from Farm 1 and Farm 2. Therefore, a third farmer (Farm 3) was found in “Fernana” region from Jendouba governorate (altitude 220 m and 700 mm of average rainfall). The collection took place in spring 2018 from a raw seed sample of faba bean from the previous cropping year (2017-2018). It concerned only the farmers who cultivate the local varieties or landraces if the selection was not been done yet on that material. The size of the collected sample was sufficiently large (1 kg/local population/farmer) to represent the existing variability. The seed sample was divided into subgroups (homogeneous lots) distinguished according to seed morphology or seed type (field bean “small type” or broad bean “large seed type”). Table 1 shows details of Vicia faba seeds collections, subpopulations and origins.

2.1.2. Collection from IPK Genebank

The Genebank of the Leibniz Institute of Plant Genetics and Crop Plant Research at the Leibniz-Institut für Pflanzengenetik und Kulturpflanzenforschung (IPK, Germany) provided the other part of seeds material (Table 1). This material corresponds to accessions collected in Tunisia in the 1970s, kept and conserved by ex situ conservation in the IPK Genebank.

2.1.3. Improved Local Varieties

In addition to these local populations, four Tunisian improved cultivars were considered as references and for comparison study (Table 1).

2.2. Phenotypic Variability of Vicia faba L. Based on Seeds Traits

The phenotypic characterization of the seeds was carried out according to several morphological and physical parameters related to the technological quality of the beans following UPOV [15] and the agri-food industry AOAC international standards [16].

2.2.1. Morphological Parameters

The following morphological parameters were considered:(i)Shape of the seed: The characterization of the seed shape is based on the standards described in the universal descriptors according to UPOV [15] (Figure 1(a)).(ii)Hilum pigmentation: The pigmentation of the hilum is observed on the dry seed and scored as follows: black (1), white, or pigment-less (0) according to Khamassi et al. [14].(iii)Seeds dimension: The three axial dimensions of the seed (length, width, and thickness) were measured using a Vernier caliper (Gilson Tools, Japan) with accuracy of 0.05 mm (Figure 1(b)). Then, 10 randomly selected seeds (blindly taken from their bags of conservation) for each variety, landrace or local population and subgroup within each of them were used to estimate seeds’ dimensions.

2.2.2. Physical Properties of Seeds

The arithmetic mean of the seed diameter (Ad, mm) was calculated according to the formula given by Joshi et al. [17, 18]. This formula is based on three dimensions: length (L, mm), width (W, mm) and thickness (T, mm) as follows:

The geometric mean of the seed diameter (Gd, mm) was calculated using the formula reported by Rotimi [19]:

The degree of seed sphericity (Φ, %) is an indicator of seed shape. According to Rotimi [19] and Mohsenin [20], the degree of sphericity was calculated based on the mean values in mm of major diameter (a), minor diameter (b), and intermediate diameter (c) of the seed, according to the following formula:

The elongation (E) and the degree of flattening (A) of the seeds (dimensionless) are also a seed shape indicators. These parameters were determined as follows:

The seed aspect ratio (Ra, %) parameter, expressed as percentage, was determined according to Hara et al. [21] and Omobuwajo et al. [22]:

The seed volume (V, mm3) was calculated according to the below formula [20]:

The surface area (Sa, mm2) was determined by the following formula [20]:

In order to determine the one hundred seed weight (OHSW, g), a sample consisting of 100 seeds of each accession was weighed using a precision scale with three random replicates.

2.2.3. Statistical Analyses

To test the existence of significant differences between the local accessions studied, a single-factor analysis of variance (ANOVA) was applied using the SAS software. The comparison of the average morphological characteristics of the seed was made by means of the LSD test at the alpha threshold of 5%. Genetic diversity among cultivars based on biological and physicochemical data was examined using the clustering analysis with Euclidean distance matrix and the principal components analysis (PCA). Diversity among the 18 faba bean genotypes was assessed by estimating Euclidean distance as the formula suggested by Shifriss and Sacks [23]:where Xik = performance of the ith parent for kth character, Xjk = Performance of the jth parent for kth character, Sk = standard deviation of the kth character, and Ed: Euclidean distance.

The analysis of the relationship between all studied parameters was performed using several statistical techniques. Principal component analysis (PCA) of 18 faba bean seeds’ physical properties data was performed with the R language. An absolute value of 0.50 was used in the loading matrices to select the characteristics in a particular principal component (PC). Correlations between variables were calculated with Spearman’s correlation coefficients. All tests were considered as statistically significant when is less than 0.05. All analyses were carried out by using the computer R programming software.

3. Results and Discussion

3.1. Morphological Seeds Features

The black pigmentation of the hilum (Figure 2(b)) is present in all accessions evaluated except the seeds of “205 Bulk” accession, which are characterized by a white (pale or colorless) hilum (Figure 2(a)). This white pigmentation of the hilum is a recessive trait with a maternal effect [14, 24]. According to Khamassi et al. [14], this trait of interest would serve as morphological marker for the selection of lines carrying the zero vicine and convicine (VC-). The variety carrying the VC- gene is of interest for monogastric feeding and as food for human deficient G6PD enzyme allowing the avoidance of metabolic disease called favism [25]. Although faba bean seeds harbor several nutritional benefits, their consumption is restricted and limited by the vicine and convicine compounds that are responsible for favism disease [2628]. Therefore, discovering a morphological trait associated with low VC phenotype would be an easy and simple marker for farmers, but also a considerable gain of time and energy in breeding programs for detecting low/free CV faba bean phenotypes [2628]. Indeed, in German winter faba beans, the hilum color has been used as morphological marker in breeding works [29]. Recently, a linkage of 5 cM has been found between Hc gene (responsible for colorless hilum) and VC- gene only in the genitor line that carries VC- gene [2628] and this is applicable only to lines where one of the parent is carrying both VC- in linkage with Hc gene. Nevertheless, not all colorless hilum faba bean genotypes are VC-, it should not be considered as pleiotropic linkage [14]. In any case, the discovery of genetic sources of low VC are an important result from our study. Actually, low VC grain content cultivars are being registered under several catalogues and databases, and their cultivation is growing in many countries [28]. Irregular seed form (Figures 2(c) and 2(d)) was the main morphology observed in all accessions, except the variety “Badii” whose seeds harbor elliptic form (Figure 2(e)). The seed shape in this case did not allow discrimination between the studied subpopulations.

3.2. Physical Properties of Seeds

The axial dimensions of seeds analyzed by accession are shown in Table 2. Analysis of variance for all axial dimensions showed a highly significant difference between the evaluated accessions (). The lengths of the studied Vicia faba seeds vary between 9 and 22.80 mm; the average widths vary from 6.5 to 17 mm and the average thicknesses from 4.20 to 6.20 mm. Those results are quite similar to those observed in previous studies that dealt with local either small or large Vicia faba germplasm collections. Indeed, seeds length varied between 11.80 and 15 mm; width ranged between 7.5 and 10 mm in an indigenous Omani faba bean collection [9]. The study carried out on Turkish faba bean varieties showed bigger values of length, width, and thickness, ranging respectively from 18 to 22 mm, 13 to 16 mm, and 13 to 16 mm [30]. The seeds of the “Memdouh,” “01_05,” and “02_05” accessions are the longest (22.45 mm) while the seeds of the accessions “202_lot 2,” “01_02,” and “FAB 6334” with a length around 9 mm are the shortest (Figure 3). However, the average length of “Chourouk” seeds is statistically similar to those of “03_04” and “02_06.” Also, the accessions “03_05” and “202_lot 3” presented an average grain length similar to that of the Najeh reference (10.70 mm). The “FAB 6434” and “Badii” accessions produced seeds ranging in length from 10 to 10.20 mm. The width of the seeds varies depending on the accession of Vicia faba (Table 2). “Memdouh” and “FAB 7111” accessions produced the largest seeds (16.7 mm), while “Badii,” “01_02,” and “205 bulk” accessions gave the narrowest seeds (6.6 mm). In addition, “02_06” and “202_lot3” accessions had an average seed width similar to those of the “Chourouk” and “Najeh” references. Statistical analysis for the arithmetic (Ad) and geometric (Gd) diameters of the seeds showed highly significant differences () between Vicia faba accessions. Table 2 shows the average values of the arithmetic and geometric diameter recorded for the different accessions. It shows that the average arithmetic diameter and the average geometrical diameter vary respectively from 6.87 to 14.93 mm and from 6.67 to 12.64 mm. The seeds of accession “Memdouh” are characterized by the highest arithmetic and geometric diameters, while the seeds of accession “01_02” are characterized by smaller arithmetic and geometric diameters compared to other accessions. The axial dimensions (length, width, and thickness) allowed deducing the shape of the seeds, characterized by their degree of flattening, elongation, average index of sphericity, and aspect ratio. Analysis of variance for the various form parameters (E, A, Φ, and Ra) showed highly significant differences between the evaluated accessions (). Table 2 shows the average values of Vicia faba seed shape. The average values of the degree of flattening are between 1.22 and 3.40. “Memdouh” accession is the most flattened with a degree of flattening of 3.40. On the other hand, “205_Bulk-2” accession has the most convex seed which is characterized by a relatively weak flattening index (1.22) (Table 2). Average values of seed elongation range from 1.16 to 1.63. The seeds of accession “01_05” are the most elongated with an average elongation of the order of 1.63, whereas the seeds of “FAB 7111” are the least elongated (E = 1.16). The “03_05,” “205_Bulk,” and “03_04” accessions have an average elongation of 1.53, which is statistically similar to that of “Badii” (1.57) (Table 2). The seed aspect ratio connects the width to the length of the seeds. The average values of this aspect ratio ranged from 0.61 to 0.87 for “01_05” and “FAB 7111” accessions, respectively (Table 2). As a result, the seeds of “01_05” have a strong tendency towards the oblong form. These results confirm the elongated shape of the seeds of accession “01_05” characterized by an average elongation of 1.63. Mean spherical indices ranged from 54.40 to 80.64% (Table 2). This gives information on the strong tendency of Vicia faba seed form to a spherical shape. The seeds of the “202_lot 2” accession are the most spherical with a mean sphericity index of 80.64%, whereas the seeds of the “Memdouh,” “01_05,” and “02_04” accessions are less spherical (54, 40%) than the other accessions. The average sphericity index (67.86%) of the “03_06,” “FAB 6641,” and “79_5” accessions are statistically similar to that of the “Chourouk” and “Badii” references. Some Indian faba bean varieties reached 80–83% of sphericity index [31]. Indeed, the study of sphericity remains necessary to deduce the shape of the seed for the design of transport and storage equipment. Ixtaina et al. [32] and Mpotokwane et al. [33] showed that a high sphericity of seeds allows them to roll rather than slide during transport. The average area of the Vicia faba seed in the studied accessions varies between 140.10 and 488.60 mm2 (Table 2). The “Memdouh,” “01_05,” and “02_05” accessions are characterized by the largest area of seeds. On the other hand, the accessions “01_02” and “FAB 6434” are characterized by a weak Sa (140.10 mm2), while the other varieties recorded intermediate values. The average volume of the Vicia faba seeds for the studied accessions is between 156 and 1020 mm3 (Table 2). The accessions “01_05,” “02_05” and the reference variety “Memdouh” display the highest average volume. On the other hand, a relatively low average volume characterized the “01_02” and “FAB 6434” accessions compared to the other accessions that recorded intermediate values. By using an artificial neural network-based method, Aboukarima et al. [13] have shown that seeds volumes of Egyptian faba bean varieties ranged from 700 to 930 mm3. The weight of one hundred seeds is a very important parameter for the estimation of yield in Vicia faba. The weight of one hundred seeds of the accessions studied varies from 47 to 248.5 g (Table 2). The obtained results showed that accession “01_02” has the lightest seeds (OHSW = 47 g), whereas “03-02” had the heaviest seeds (248.5 g) which is almost twice the one obtained by Hacıseferoǧulları et al. [30] for Turkish varieties. A Palestinian landrace of faba bean could reach 239 g in terms of OHSW [34]. Upadhyaya et al. [35] have argued that the OHSW parameter is highly heritable and often used for the distinction between genotypes. Several authors have reported that OHSW is necessary for the design of aerodynamic cleaning equipment [32]. In addition, this parameter is very useful for calculating the sowing dose. There was variability in the measured dimensions between the accessions. This parameter makes it possible to have an idea about the varietal genetic diversity and could be a determining parameter in the selection. Moreover, Koura et al. [36] have shown that physical properties vary according to variety. These characteristics are very useful in the different separation, sorting, and conditioning processes [37].

3.3. Classification of the Studied Faba Bean Germplasm Accessions Based on Morphological Characters Using Clustering by Euclidean Distance and Principal Component Analysis
3.3.1. Hierarchical Classification Using Matrix of Euclidean Distances

The genetic similarity dendrogram is an effective means for evaluating the genetic relationship between the different studied accessions. Euclidean distances were calculated between all pairs of individual trees based on standardized values of eleven quantitative seed traits and the matrix of Euclidean distance was used for cluster analysis in R language, as shown in Figure 4. The graphical representation of the matrix of Euclidean distances, which lies between the studied twenty-four faba bean genotypes based on morphological characteristics, revealed three different groups. Group 1 included ten faba accessions (“03_03,” “03_07,” “0.3_02,” “03_01,” “FAB 5042,” “FAB 7111,” “Memdouh,” “02_04,” “02_05,” and “01_05”). Group 2 combined the two faba landraces (“202 lot 2” and “202 lot 3”). Group 3 contained twelve faba accessions (“205 Bulk,” “Najeh,” “01_02,” “02_06,” “03_04,” “Chourouk,” “03_06,” “FAB 6641,” “FAB 6434,” “Badii,” “79-5,” and “03_05”).

3.3.2. Multivariate Analysis by PCA Method

The interaction between eleven quantitative seed traits were studied by the PCA which serves to appreciate which combination of seed traits could leads to high quality of the landraces, varieties, and local populations of faba bean germplasm in terms of their significance and value in marketing or for adoption by farmer according to seed size or suitability for industrial use. The multivariate PCA analysis was performed by using R language as illustrated in Figure 5. The first two components of the PCA explained 87.38% of the total variation among a subset of 24 of faba beans germplasm included in the analysis. The first component (axis 1) explained 67.36% of the variation, followed by 20.02% for the second components (axis 2). Axis 1 was positively correlated with the arithmetic diameter (0.997), seed length (0.989), geometric diameter (0.988), grain area (0.984), seed width (0.977), and seed volume (0.975) and negatively correlated with sphericity index with a value of −0.847. On the other hand, axis 2 was positively correlated with the percentage of aspect ratio (0.992) and negatively correlated with the seed elongation with a value of −0.992. Moreover, we showed that the PCA method classified the studied faba accessions into four main clusters (Figure 4), where cluster1 was composed by “202 lot 2,” “202 lot 3,” “Najeh,” “01_02,” “03_06,” and “FAB 6641,” whereas cluster 2 grouped “FAB 7111,” “FAB 5042,” “03_03,” “03_06,” “0.3_01,” “03_02,” and “Memdouh.” Cluster 3 combined “FAB 6434,” “205 Bulk,” “79-5,” “03_04,” “03_05,” “02_06,” “Badii,” and “Chourouk.” Finally, cluster 4 grouped “02_04,” “02_05,” and “01_05.” The regrouping of the latter accessions in the same cluster 4 by PCA analysis was essentially related to a remarkable similarity of their large seed size, aspect ratio, and seed elongation, while the eight faba accessions composing the cluster 3 had similarity of the two quantitative seed traits (sphericity index and seed elongation). On the other hand, cluster 1 grouped six faba accessions that had resemblance in numerous seed trait parameters such as seed elongation, arithmetic diameter, seed length, geometric diameter, grain area, seed width, and seed volume. Finally, the seven faba accessions included in cluster 2 had high similarity in the aspect ratio, arithmetic diameter, seed length, geometric diameter, grain area, seed width, and seed volume.

3.4. Correlation Study between the Physical Characteristics of the Studied Seeds

The association between the different measured parameters has been established by Spearman’s correlation test. The correlation matrix showed that there was a highly positive significant correlation (r = 0.948, ) between mean length and mean width. The accessions of Vicia faba with the longest seeds are also the largest. In addition, the average length was high positively related to other physical parameters such as arithmetic diameter (r = 0.986, ), geometric diameter (r = 0.973, ), surface area (r = 0.973, ), degree of flattening (r = 0.878, ), and volume (r = 0.969, ), and high negatively influenced by the sphericity index (r = − 0.903, ). On the other hand, the average width was high positively correlated with arithmetic diameter (r = 0.973, ), geometric diameter (r = 0.970, ), surface area (r = 0.970, ), degree of flattening (r = 0.908, ), and volume (r = 0.970, ), and high negatively influenced by the sphericity index (r = −0.788, ). Therefore, this showed that the axial dimensions had a positive effect on certain physical properties and that the longest and widest seeds have a low sphericity index. Furthermore, it was noted that there was not a significant connection of the average thickness with all the studied parameters like the average width (r = −0.01, ) and the average length of the seeds (r = −0.32, ). Moreover, the aspect ratio parameter was high negatively correlated with the elongation (r = −0.990, ). Additionally, the average volume of the seeds was positively associated with the parameters length (r = 0.969, ), width (r = 0.970, ), arithmetic diameter (r = 0.994, ), geometric diameter (r = 0.998, ), flattening (r = 0.853, ), and surface area (r = 0.998, ), and negatively correlated with sphericity (r = −0.809, ). However, there was no correlation between thickness and the flattening index (r = −0.368, ). It has been also found that the OHSW is positively correlated with certain physical parameters such as length, width, arithmetic diameter, geometric diameter, surface, and volume. On the other hand, this parameter is negatively correlated with the sphericity index (r = −0.7637, ). It follows from all the foregoing that the seeds of the accessions “Memdouh,” “01_05,” “02_05,” “02_04,” “03_01,” “FAB 7111,” “03_02,” “FAB 5042,” “03_03,” and “03_07” are correlated to lengths. The highest averages are the widest, the most flattened with the largest arithmetic and geometric diameters and also they possess the largest surface area and the highest volume. Nevertheless, they are characterized by the lowest sphericity index.

4. Conclusions

Even the morphological characters (shape of seeds and hilum color) of seeds did not allow significant discrimination among the Tunisian faba bean germplasm, our results highlighted the presence in our studied germplasm of one accession (“205 bulk”) with low VC grain content. Indeed, the research of low VC content genotypes started from the 1970s by screening hundreds of lines, and the discovery of the first low VC genotype was in 1989 by Duc et al. [38]. Particularly, the “205 bulk” accession that could a valuable candidate genitor for potential breeding programs aiming to improve nutritional value and food quality of faba bean cultivars by crossing with low VC genotype. Moreover, with almost 250 g of OHSW, the “03_02” accession is suitable for horticultural practices and very interesting for consumers. Overall, the physical properties of seeds appeared to be interesting markers for Tunisian faba bean germplasm discrimination. Moreover, several seed characters were positively correlated one to each other, helping in screening for better faba bean quality phenotype. Our work led to the identification of some Tunisian faba bean accessions with desirable traits, which constitutes a valuable source for breeding material with the aim of improving agronomic and nutritional traits as well as adaptative characters. From all the foregoing, this work has provided useful data for the technological treatment of seeds, food quality improvement, elaboration of strategies for the conservation and sustainable management of the genetic heritage, and better use of germplasm and for global genetic improvement of Vicia faba.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors wish to thank the Ministry of Agriculture, Hydraulic Resources and Maritime Fisheries and the Ministry of Higher Education and Scientific Research for their financial support to this work. They are also thankful to the IPK Genebank for free donation of the Tunisian faba beans conserved in their genebank.