Next Article in Journal
Remote Detection of Large-Area Crop Types: The Role of Plant Phenology and Topography
Previous Article in Journal
Assessment of Potential Land Suitability for Tea (Camellia sinensis (L.) O. Kuntze) in Sri Lanka Using a GIS-Based Multi-Criteria Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphological and Agronomic Characterization of Spanish Landraces of Phaseolus vulgaris L.

1
Institute for Plant Molecular and Cell Biology (IBMCP, UPV-CSIC), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
2
Laboratory for Plant Breeding and Conservation of Genetic Resources, Bioplant Center, University of Ciego de Avila, 69450 Ciego de Ávila, Cuba
3
Institute for Conservation and Improvement of Valencian Agrodiversity (COMAV, UPV), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
4
Mediterranean Agroforestry Institute (IAM, UPV), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
*
Author to whom correspondence should be addressed.
Submission received: 28 May 2019 / Revised: 26 June 2019 / Accepted: 3 July 2019 / Published: 9 July 2019

Abstract

:
Beans (Phaseolus vulgaris L.) originated on the American continent, specifically in the Mesoamerican zone, and their domestication took place independently in the Mesoamerican area and the Andean zone, giving rise to two well-differentiated genetic pools. It was also noted that the Andean wild populations originated from only a few thousand individuals from the Mesoamerican wild populations, which produced a great bottleneck in the formation of the Andean population. During centuries of cultivation in the Iberian Peninsula after its introduction in the 16th century, beans adapted to new environments, evolving numerous local landraces. Twenty-four local landraces of P. vulgaris from Spain were analyzed in the greenhouse during two consecutive seasons. From each genotype, five plants were grown and characterized for 17 quantitative and 15 qualitative traits using the International Board for Plant Genetic Resources (IBPGR) descriptors. Data were analyzed statistically by analysis of variance (ANOVA), principal component analysis (PCA), and cluster analysis. The results obtained indicate a high variability for most traits, especially those related to the yield and its components. The PCA and cluster analysis separated the landraces according to the color of the seed, the yield, and the pod and seed traits related to yield. Numerous traits exhibited interactions between the genotype and the environment. Most accessions reached higher yields in spring, in which solar radiation favors photosynthesis and, consequently, photoassimilation. The different response to the changing environment of the set of accessions studied in the present work is of great interest, and it can be exploited in breeding cultivars adapted to a broader range of environmental conditions.

1. Introduction

Common bean (Phaseolus vulgaris L.) is one of the most important legume crops. The world production of beans in 2017 was 57,496,465 tons, cultivated on 38,229,984 hectares. Asia is the continent with the highest production at 38,045,186 tons, and Europe is the continent with the lowest production at 2,305,846 tons. In the European Union, the countries with the highest production are France, Lithuania, and Spain with 370,050, 230,112, and 183,324 tons, respectively. Castilla y León is the community with the largest area and production of beans in Spain, followed by Galicia and Asturias [1].
According to archeological, botanic, historical, and linguistic data, Gepts and Debouck [2] deduced that the common bean originated in the area between northern Mexico and northeastern Argentina. Due to the wide distribution of its wild ancestor through Mesoamerica and South America, the species was domesticated in two separate areas, the Mesoamerican and the South Andean. This resulted in two different gene pools based on differences in the phaseolin seed protein [3], molecular markers [4,5,6,7], and morphological traits [2,8]. The Mesoamerican gene pool extends from Mexico through Central America to northern Colombia, whereas the Andean one covers Peru, Chile, Bolivia, and northern Argentina. However, more recently, Bitochi et al. [9] conducted a study to solve the ongoing debate on the origins of the common bean, by investigating the nucleotide diversity at five gene loci of a large sample that represents the entire geographical distribution of the wild forms of this species. Their results clearly indicate a Mesoamerican origin of the common bean and the occurrence of a bottleneck during the formation of the Andean gene pool.
The first varieties of P. vulgaris, corresponding to small-grain Mesoamerican genotypes, were introduced in Spain and Portugal around 1506 [10] and from South America around 1532 [11,12]. Later on, gradually, Andean genotypes of large grain were also introduced as new territories were explored in the New World [11]. This germplasm spread throughout the Spanish territory and adapted to new environments, in terms of temperature and cropping systems, especially in the north and northwest of the Iberian Peninsula, which is the area where the greatest genetic variability was found [13]. On the other hand, the cultivation of germplasm of both genetic stocks occurred in proximal areas, in small farms, which favored the genetic flow between them, contrary to the genetic isolation that existed in America. This led to recombination between both genetic stocks, resulting in new genetically intermediate forms with both Mesoamerican and Andean phenotypes. The Iberian Peninsula can be considered as a secondary center of diversification of this species since beans were cultivated here for centuries in different agrosystems [14]. Although commercial varieties largely displaced landraces, there are still farmers who grow them for self-consumption and sale in local markets, especially in the north and northwest of Spain, where there is a marked tradition of this crop.
A huge number of accessions of this crop were collected and conserved in germplasm banks, and its genetic diversity was assessed in many different studies. Some of the most important collections are stored at the Center for Tropical Agriculture (CIAT, Cali, Colombia), the Western Regional Plant Introduction Station, Pullman (Washington, USA), Institut für Pflanzengenetik und Kulturpflanzenforschung (IPK, Gatersleban, Germany), and the Vavilov Institute of Plant Industry (VIR, Russian Federation) [15]. In Spain, there are several institutions conserving accessions of the genus Phaseolus. The most important one is maintained in the Spanish Center for Genetic Resources, (CRF, Alcalá de Henares, Madrid), followed by the Biological Mission of Galicia (MBG, Pontevedra), the Agricultural Technology Institute of Castilla y León (ITACyL, Valladolid), the Agrifood Research and Technology Center of Aragon (CITA, Zaragoza), and the Institute of Conservation and Improvement of the Valencian Agrodiversity (COMAV, Polytechnic University of Valencia).
Average yields of common beans, as for practically all major crops, are much below 20–50% of their record yields, mostly due to abiotic stress factors, such as drought, soil salinization, or extreme temperatures, among others [16]. The forecasted effects of global warming are mostly longer and more frequent periods of drought, heat waves, and irregular meteorological patterns. They are already noticed in large areas of the world, including the Mediterranean [17]. These factors, in combination with the shortage of water available for irrigation, will probably increase the losses in agricultural production. Commercial varieties, generally selected on criteria related to their yields or tolerance to pests and parasites, are generally sensitive to abiotic stress. One of the challenges of modern agriculture is the use of genotypes more tolerant to drought and salinity; therefore, the screening of neglected varieties and local landraces stored in germplasm banks acquires a special relevance [18]. Phenotyping common beans for adaptation to abiotic stress, and analysis of its physiological mechanisms of response to stress were the object of numerous studies in the last decade [19,20,21], among many others. Similar comparative studies were performed by our group on different species of Phaseolus [22,23,24]. In this work, we present an exhaustive morphological characterization of 24 accessions of common bean local landraces stored at the seed bank of COMAV; they were chosen, maximizing their geographical origin and their morphological diversity, for studying their adaptation to drought and salinity (manuscript in preparation). The study of the phenotypic traits is essential to know the agronomic potential of these landraces and their use in plant breeding. The high variability found for morphological and agronomic traits in the different types of bean studied is remarkable and makes this set of accessions of great interest for their use in breeding.

2. Materials and Methods

2.1. Seed Material

Twenty-four local landraces of P. vulgaris from Spain were characterized for two consecutive growing seasons, in autumn and spring. Accessions were collected in the 1990s directly from growers. The origin of the analyzed accessions with their reference codes is indicated in Figure 1.
Six seeds of each genotype were germinated in individual pots of 1.6 L. Once plants formed the first three trifoliate true leaves and reached a minimum of 20 cm in height, they were transplanted to pots of 7 L on a substrate with peat–perlite–vermiculite (2:1:1) in the greenhouse. The experimental design included five individuals from each accession in one plot, with sufficient separation between them. Watering was carried out by drip irrigation with Hoagland solution provided by the greenhouse. The average temperature during the trial carried out in autumn was 20.3 °C (ranging from 15.9 to 24.1 °C), the relative humidity was 86.85%, and the solar radiation was 38.4 W∙m2. In spring, the average temperature was 20.4 °C (15.3–24.1 °C), relative humidity was 84.41%, and radiation 85.5 W∙m2 (Figure 2).

2.2. Descriptors Used for the Agronomic Characterization

During the two cycles of growth, 17 quantitative and 15 qualitative traits related to vegetative characteristics, flowers, pods, and seeds were recorded (Table 1). The descriptors developed by the International Board for Plant Genetic Resources (IBPGR) were used [25]. This descriptor list was prepared by IBPGR in a standard format following the advice of crop experts throughout the world. Five plants, and five pods and ten seeds from each plant were analyzed for each genotype in each cycle. Pods were collected when they were completely ripe. Productivity was calculated by multiplying the number of pods per plant by the number of seeds per pod and by the average weight of one seed.

2.3. Data Analysis

A two-way analysis of variance (ANOVA) was performed for all quantitative traits to check the interaction between the genotype and the environment. Quantitative traits were summarized as mean, standard deviation, coefficient of variation, and minimum and maximum values. The qualitative characteristics were expressed on scales, as shown in Table 1, and then graphically represented in the shape of stacked bar charts. A principal component analysis (PCA) was used to examine the association between the analyzed traits and the similarity among accessions. PCA was performed with all quantitative traits and the seed color, which was categorized using increasing values for darker colors. Data were analyzed using the software Statgraphics Centurion v.16 (Statpoint Technologies, Warrenton, VA, USA).
Additionally, for the quantitative traits, a Euclidean distance matrix based on standardized data was computed for clustering analysis by using the UPGMA (unweighted pair group method with arithmetic mean) method [28,29]. Goodness of fit of the resulting dendrogram was evaluated with the cophenetic correlation coefficient by using the Mantel test [30]. This analysis was made using the NTSYSpc2.0 software package (Applied Biostatistics Inc., Setauket, NY, USA).

3. Results

The seed morphology of the accessions was very variable in shape, size, and color (Figure 3). Five accessions, BGV1165, BGV3618, BGV3941, BGV4161, and BGV15856, exhibited a determinate growth habit; therefore, their plant height and yield-related traits differed significantly from those with indeterminate growth. For this reason, these accessions were excluded from the PCA.
A two-way ANOVA was run to check the possible effect of the environment and the interaction between the genotype and the environment. The results indicated the existence of a significant effect of the environment and the interaction between the genotype and the environment for the majority of traits, especially those related to yield and its components (Table 2). For this reason, data of the two seasons were not merged and were analyzed separately in the PCA and basic statistics calculations.

3.1. Principal Component Analysis and Cluster Analysis

To find out which quantitative traits explained most of the variation in the set of accessions studied, a principal components analysis was carried out. The first three components explained 50.6% of the variation contained in the data (21.7%, 16.7%, and 12.7% for the first, second, and third components, respectively). The traits more positively correlated with the first PC component were those related to the yield (Yield) and its components (PodL, Pplant, SeedWh, SeedL, SeedW, SeedH) (Table 3, Figure 4).
The color of the seed coat (Scoatdc) and the apparent seed veining (Aseed) were negatively correlated. Traits related to the size of the leaves (LeafL and Wcshe) were best correlated with the second principal component. Some of the components of yield, pods per plant (Pplant), and seeds per pod (SeedP) also had high positive correlations, while the seed coat color (Scoatdc), seed weight (SeedWh), and seed shape (Shseed) were negatively correlated (Table 3, Figure 4).
The accessions analyzed were dispersed onto the two main components on the scatter plot of the PCA scores, indicating a large variability (Figure 5). It was also patent a separation of the accessions according to the size and color of their seeds. The accessions with black seeds were located at the top of the graph, mainly due to the seed coat color and to the larger size of their leaves, the greater number of seeds per pod, and the narrower shape of the seeds. On the left side of the graph, mostly cultivars with whitish seeds, lower yield, and without apparent seed veining were grouped. Within this group, there was also a remarkable separation according to the second PC, mainly related to the size of leaves and some components of yield. The accessions with brown seeds were also very dispersed, but with a predominance on the right side of the graph, due to their higher yield and larger number of seeds per pod of most cultivars of this type. There were three accessions with different seed coat color; one of them, BGV143, had yellow to greenish-yellow seeds. This accession had very low yield and was located in the lower part of the graphic. Contrarily, BGV1169, with maroon seed coat, was highly productive and grouped with the brown seeds. The accession BGV1185, with gray, brownish to greenish seed color, grouped with the whitish ones, with a very low yield (Figure 5).
Regarding the cluster analysis, in the dendrogram obtained with the morphological traits by the UPGMA method, two main clusters can be observed (I and II, Figure 6). Cluster I includes the accessions encircled in yellow in Figure 5 and cluster II includes the ones encircled in green color. As expected, the accessions in the cluster analysis were grouped according to the traits related to yield (Yield) and its components (PodL, Pplant, SeedWh, SeedL, SeedW, SeedH). These traits were the most positively correlated with the first PC component in the PCA. In fact, the average values for these traits in the group of accessions included in cluster I (PodL = 15.44; Pplant = 48.91; SeedWh = 54.00; SeedL = 14.16; SeedW = 6.28; SeedH = 8.18; and Yield = 151.52) were consistently higher than the ones for the accessions included in cluster II (PodL = 13.71; Pplant = 22.83; SeedWh = 40.56; SeedL = 11.45; SeedW = 5.63; SeedH = 7.31; and Yield = 41.25). Also, as in the PCA analysis, the color of the seed was another relevant trait which contributed to the separation of the accessions. Thus, in cluster I, 75% of the accessions had seeds of brown color, while 55% of the seeds included in cluster II were white. The cophenetic correlation coefficient of the phenogram obtained from the 20 descriptors used in the analysis had a value of 0.86.

3.2. Description of the Variability of the Whole Set of Accessions and of the Groups Established

For the description of the variability found in the set of accessions, we grouped them according to the seed coat color, establishing three groups: black, brown, and whitish. The three accessions with different colors, BGV143, BGV1169, and BGV1185, are described separately, as well as the five accessions of determinate growth habit.

3.2.1. Variability of the Whole Set of Accessions

The whole set of accessions exhibited a high variability for most of the traits evaluated (Table 4 and Table 5). Especially relevant was the variation for yield and number of pods per plant, which ranged from 2.8 to 242.8 g per plant and from five to 72, respectively. The coefficient of variation for these traits was very high, 82.4 for yield and 59.4 for pods per plant. There was also a great variability for all traits related to inflorescence and some of the components of yield, with coefficient of variation (CV) values that ranged from 31.7 to 40.4. The traits that showed lower variation were those related to vegetative characteristics and the height and width of the seeds; in these cases, the CV values ranged from 15.5 to 29. This trend was similar in the two trials performed. Regarding the qualitative traits recorded, there was little variation for pod cross-section, which was predominately pear-shaped, with pod curvature being straight in general, pod beak position being downward in most cases, and apparent seed veining being absent in general. For the other traits, there was high variation (Figure 7 and Figure 8).

3.2.2. Description of the Three Accessions with Black-Colored Seeds

Accessions with black seeds were the most leafy, with a higher mean value of central leaflet length and width, long inflorescences and peduncles, and a high number of pods per plant (Pplant) (on average 43.47 in the first cycle and 44.40 in the second, compared with the mean values for the whole set of accessions that were 29.15 and 30.14, respectively). These cultivars had the smallest seeds and, therefore, more seeds per pod than those from other groups (Table 4 and Table 5). The average yield was high (92 in the first cycle, and 108 in the second), with a high coefficient of variation, as the yield ranged from 34 to 193 in the first cycle and from 33 to 239 in the second cycle. Qualitative traits (Figure 7 and Figure 8) showed only a relatively smaller variation in comparison to the other group of cultivars, without any seed coat patterns or apparent seed veining, with seeds either matte or medium brilliant, and only two shapes, either oval or cuboid. No variation was registered for the trait pod beak position (marginal in all cultivars), pod beak orientation (upward), or pod cross-section (pear-shaped).

3.2.3. Description of the Six Accessions with Brown-Colored Seeds

Accessions with brown seeds included the most robust plants, with plant height ranging from 1.75 to 2.95 m in the first cycle, and from 1.6 to 2.98 m in the second one, with a larger stem diameter (Table 4 and Table 5). The highest average yield was found for this group; the coefficient of variation for this characteristic was very high, 46% in both cycles, mainly due to the large variation in the number of pods per plant, ranging from five to 65 and from four to 65, respectively. The group had the largest seeds, most notable in the weight of seeds (100 seeds weighed 60 g in the first cycle and 54 g in the second one). Among the qualitative characteristics (Figure 7 and Figure 8), the seed coat pattern was outstanding and the most variable in comparison to the other groups. Furthermore, high variability was registered for the shape of seeds (including cultivars with ovoid, cuboid, and kidney-shaped seeds), the brilliance of seeds (medium and shiny), variable pod beak orientation (upward, straight, downward), and pod suture string (from stingless to very stringy). No variation was found for apparent seed veining, being absent in all cases.

3.2.4. Description of the Seven Accessions with White-Colored Seeds

Accessions with whitish seeds had generally smaller plants, with fewer pods per plant (26 in both cycles), smaller seeds, with a weight of 39 g for 100 seeds (in both cycles). The number of seeds per pod was also generally smaller and, therefore, the mean yield was much lower in comparison to the other accessions (54 g per plant in the first cycle and 51 in the second one), representing less than half of that registered for the group with brown seeds (Table 4 and Table 5). However, the coefficient of variation of this trait was high, as some cultivars were much more productive than the others (e.g., BGV 3176 with 132 g per plant). Among the qualitative traits, the seed shape was more variable than in other groups; most of the cultivars had kidney-shaped seeds, but there were also some with cuboid, round, and oval seeds (Figure 8). Seeds coat patterns were absent, but seed brilliance varied from matte to shiny, and an apparent seed veining was detected only in this group.

3.2.5. Description of the Accessions with Different Seed Colors and of the Accessions with a Determinate Growth Habit

There were three accessions with seed coat color different from the ones described above. The accession BGV1169 had maroon seeds, BGV1185 exhibited gray, brownish-greenish seeds, and BGV143 had yellow to greenish-yellow-colored seeds. The three accessions differed in many traits. For vegetative traits, BGV1169 developed bigger leaves, both in length and width, and a short period from flowering to the harvest of the pods (59 days). It was also the most productive, with 182 g per plant, with higher values for the components of yield. In contrast, the accession BGV1185 was the one with the lowest yield (20.8 g per plant) and with the most prolonged period from flowering to harvest (120 days). The accession BGV143 had intermediate values for all traits, as well as outstanding inflorescence length and pedicel length; it showed the shortest period from flowering to harvest (52 days). The behavior of the three accessions was similar in both cycles. The growth habit of BGV143, BGV1169, and BGV1185 was climbing indeterminate. The color of the pod was light green for the cultivars BGV143 and BGV1185 and normal green for BGV1169. The cross-section and curvature of the pod were pear-shaped and slightly curved for all cultivars. BGV1169 did not present ventral suture in the sheaths, while the rest were very fibrous. The pods of the cultivars BGV1169 and BGV1185 at physiological maturity were yellow and pale yellow for BGV143 with colored mottling or stripes. The position of the pods in the plant was combined for BGV143 and BGV1185, while for BGV1169 the pods were concentrated in the center. The position and orientation of the apex of the pods were marginal and downward (downward). The dry pod color was golden or deep yellow (BGV1169) and carmine red (BGV143 and BGV1185). The seed coat patterns were absent for the three cultivars, and the seed brilliance was medium for BGV1169 and BGV143 and shiny in BGV1185. The seed shape was kidney-shaped (BGV1169 and BGV1185) and kidney pronged (BGV143). In all three cultivars, no ribs were observed in the seeds.
Regarding the accessions with determinate growth habit, there were two with brown seeds and three with whitish seeds. All of them had small plants, less than 0.5 m in height with narrow leaves. Accessions with brown seeds were very early and showed a short period from flowering to harvest (46 days). All accessions exhibited a high number of flower buds per inflorescence and low numbers of pods per plant and seeds per pod, giving rise to a low yield. Accessions with brown seeds were more productive than those with whitish seeds, because of their higher pod length, number of seeds per pod, and length and width of the seeds.

4. Discussion

Analysis of the complete set of accessions highlighted the great variability for most of the 32 characteristics evaluated. A high diversity was also reported in previous works on bean collections of Spanish origin [31,32]. It is also evident the diversity of types in terms of size, color, and shape of seeds, as well as for other characteristics related to the main differences between the two gene pools, the Mesoamerican and the Andean; they include the size and shape of the terminal trifoliolate leaflet, flower bracteoles, the presence or absence of stripes at the outer base of the flower standard, and the pod beak position. Our findings support results obtained by other authors who demonstrated that, in the Iberian Peninsula, after centuries of joint cultivation, a mixture of the two pools of germplasm occurred. This was already revealed, based on morphological studies, phaseolin protein patterns, and molecular markers [32,33]. The existence of both pools represents an added value to the Spanish germplasm. This high variability of Spanish accessions is of great value for their adaptation to the environmental conditions of both areas of domestication of this crop. Additionally, they may be used as a bridge for the transfer of genes between accessions belonging to both gene pools, as already suggested by Rodiño et al. [32].
The cluster and PCA analyses were congruent, and they clearly separated the accessions for their yield and yield components, mainly those related to pods and seeds, as well as their seed color. These characteristics are much more important in explaining the diversity of bean types in Spain, in contrast to other traits related to plant growth and earliness. The lack of grouping of accessions of different origin reflects a lack of regional differentiation, probably caused by transport by humans and seed exchange between farmers.
The high number of traits exhibiting interactions between the genotype and the environment, when cultivated in spring and autumn, is in agreement with previous reports [34,35]. Also, Nienhuis et al. [36] and Sills et al. [37] reported environmental effects in the expression of quantitative traits in different types of beans. This means that the best genotype for one environment may not be the best for another one and, therefore, specific types of beans must be developed for cultivation in different environments [38]. In our trial, the main differential environmental factor between the two growing seasons, autumn and spring, was the solar radiation, as shown in Figure 2, while the temperature and relative humidity were controlled in the greenhouse. Approximately 50% of the accessions did not show a great variation in yield when comparing the two seasons. However, the other accessions showed variation, most reaching higher yields in spring, in which solar radiation favors photosynthesis and, consequently, photoassimilation, regardless of the geographic origin of the accessions. The different response to the changing environment of the set of accessions studied in the present work is of great interest, and it can be exploited in breeding cultivars adapted to a broader range of environmental conditions.
The standard deviation values for seed weight, seed length, and seed width were higher than those reported by Escribano et al. [39], who studied a set of 66 landraces collected in Galicia (Spain). The same occurred for traits related to vegetative characteristics, such as length (LeafL) and width of the leaflet (Wcshe), and other traits of the pods and seeds. The range of variation for these traits was also higher in the set of accessions studied in the present work than in the germplasm studied by Rodiño et al. [32], which consisted of a group of 388 cultivars collected in the major production regions of the Iberian Peninsula and a core collection of 52 accessions constructed from the set of 388 cultivars.
The coefficient of variation compares the relative amount of variability between crop traits, and it is an excellent indicator of the exploitable variability for breeders. It is remarkable that high values were found in this work for many traits, especially those related to the yield and its components, as well as for other traits related to flowering and earliness. This high source of variability can be employed to develop varieties with very different characteristics, adapted to different environments, and also tailored to market preferences, in both groups of cultivars, with whitish and brown seeds, which are the most economically interesting groups. The accessions studied in this work were tested for their adaptation to abiotic stresses. Some of them showed good behavior under stress conditions. The results are being analyzed, and a manuscript is in preparation. Those data, together with the morphological description reported here, will help breeders select the most promising accessions to be used in breeding programs. The germplasm employed in this work is available upon request in the vegetable gene bank of the COMAV of the Polytechnic University of Valencia.

Author Contributions

Conceptualization, M.J.D. and M.B.; Methodology, S.A. and L.Y.; Software, J.T.; Validation, M.J.D. and O.V.; Formal Analysis, E.S.; Investigation, S.A. and J.T.; Resources, M.J.D. and O.V.; Data Curation, E.M.; Writing-Original Draft Preparation, S.A.; Writing-Review & Editing, M.B., M.J.D. and O.V.; Visualization, J.T.; Supervision, M.B.; Project Administration, O.V.; Funding Acquisition, M.J.D. and O.V.

Funding

This research was funded by the Spanish “Ministerio de Economía, Industria y Competitividad”, grant number RFP2015-00013-00-00 and cofinanced with FEDER funds.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Food and Agriculture Organization of the United Nations (FAO). FAOSTAT. Available online: http://faostat3.fao.org/browse/Q/QC/S (accessed on 10 January 2019).
  2. Gepts, P.; Debouck, D.G. Origin, domestication and evolution of the common vean (Phaseolus vulgaris, L.). In Common Bean: Research for Crop Improvement; van Schoonhoven, A., Voysest, O., Eds.; CIAT: Wallingford, UK, 1991; pp. 7–43. [Google Scholar]
  3. Gepts, P. Biochemical evidence bearing on the domestication of Phaseolus (Fabaceae) beans. Econ. Bot. 1990, 44, 28–38. [Google Scholar] [CrossRef]
  4. Koening, R.; Gepts, P. Allozyme diversity in wild Phaseolus vulgaris: Further evidence for two major centers of genetic diversity. Theor. Appl. Genet. 1989, 78, 809–817. [Google Scholar] [CrossRef] [PubMed]
  5. Tohme, J.; Gonzalez, D.O.; Beebe, S.; Duque, M. AFLP analysis of gene pools of a wild bean core collection. Crop Sci. 1996, 36, 1375–1384. [Google Scholar] [CrossRef]
  6. Beebe, S.E.; Skroch, P.W.; Tohme, J.; Duque, E.; Pedraza, M.C.; Nienhuis, F. Structure of genetic diversity among common bean landraces of Middle American origin based on correspondence analysis of RAPD. Crop Sci. 2000, 40, 257–272. [Google Scholar] [CrossRef]
  7. Duran, L.A.; Blair, M.W.; Giraldo, M.C.; Macchiavelli, R.; Prophete, E.; Nin, J.C.; Beaver, J.S. Morphological and molecular characterization of common bean landraces and cultivars from the Caribbean. Crop Sci. 2005, 45, 1320–1328. [Google Scholar] [CrossRef]
  8. Singh, S.P.; Gepts, P.; Debouck, D.G. Races of common bean (Phaseolus vulgaris, Fabaceae). Econ. Bot. 1991, 45, 379–396. [Google Scholar] [CrossRef]
  9. Bitocchi, E.; Nanni, L.; Bellucci, E.; Rossi, M.; Giardini, A.; Zeuli, P.S.; Logozzo, G.; Stougaard, J.; McClean, P.; Attene, G.; et al. Mesoamerican origin of the common bean (Phaseolus vulgaris L.) is revealed by sequence data. Proc. Natl. Acad. Sci. USA 2012, 109, E788–E796. [Google Scholar] [CrossRef] [PubMed]
  10. Ortwin-Sauer, C. The Early Spanish Men; University of California Press: Berkeley, CA, USA, 1966. [Google Scholar]
  11. Brucher, O.B.; Brucher, H. The South American wild bean (Phaseolus aborigeneus Burk.) as ancestor of the common bean. Econ. Bot. 1976, 30, 257–272. [Google Scholar] [CrossRef]
  12. Debouck, D.G.; Smartt, J. Bean. In Evolution of Crop Plants, 2nd ed.; Smartt, J., Simmonds, N.W., Eds.; Longman Scientific and Technical: Harlow, UK, 1995; pp. 287–296. [Google Scholar]
  13. Santalla, M.; Rodiño, A.P.; De Ron, A.M. Allozyme evidence supporting southwestern Europe as a secondary center of genetic diversity for common bean. Theor. Appl. Genet. 2002, 104, 934–944. [Google Scholar] [CrossRef]
  14. Pinheiro, C.; Baeta, J.P.; Pereira, A.M.; Dominguez, H.; Ricardo, C. Mineral elements correlations in a Portugese germplasm collection of Phaseolus vulgaris. Integrating Legume Biology for Sustainable Agricultura. In Proceedings of the 6th European Conference on Grain Legumes, Lisbon, Portugal, 12–16 Noviembre 2007; pp. 125–126. [Google Scholar]
  15. Lioi, L.; Piergiovanni, A.R. European common bean. In Genetic and Genomic Resources of Grain Legume Improvement; Singh, M., Upadhyaya, H.D., Bisht, I.S., Eds.; Elsevier: Amsterdam, The Netherlands, 2013; pp. 11–40. [Google Scholar]
  16. Boyer, J. Plant productivity and environment. Science 1982, 218, 443–448. [Google Scholar] [CrossRef]
  17. Cramer, W.; Guiot, J.; Fader, M.; Garrabou, J.; Gatusso, J.P.; Iglesias, A.; Lange, M.A.; Lionello, P.; Llasat, M.C.; Paz, S.; et al. Climate change and interconnected risks to sustainable development in the Mediterranean. Nat. Clim. Chang. 2018, 8, 972–989. [Google Scholar] [CrossRef]
  18. Fita, A.; Rodriguez-Burruezo, A.; Boscaiu, M.; Prohens, J.; Vicente, O. Breeding and domesticating crops adapted to drought and salinity: A new paradigm for increasing food production. Front. Plant Sci. 2015, 6, 978. [Google Scholar] [CrossRef] [PubMed]
  19. Santos, M.; Ribeiro, R.M.E.; Pimentel, C. Photosynthetic parameters and leaf water potential of five common bean genotypes under mild water deficit. Biol. Plant. 2009, 53, 229–236. [Google Scholar] [CrossRef]
  20. Rosales, M.; Rodriguez-Valentin, R.; Olvera-Carrillo, Y.; Acosta-Gallegos, J.; Covarrubias, A. Physiological analysis of common bean (Phaseolus vulgarisL.) cultivars uncovers characteristics related to terminal drought resistance. Plant Physiol. Biochem. 2012, 56, 24–34. [Google Scholar] [CrossRef] [PubMed]
  21. Beebe, S.; Rao, I.; Blair, M.; Acosta-Gallegos, J. Phenotyping common beans for adaptation to drought. Front. Physiol. 2013, 4, 35. [Google Scholar] [CrossRef] [Green Version]
  22. Al Hassan, M.; Morosan, M.; Lopez-Gresa, M.; Prohens, J.; Vicente, O.; Boscaiu, M. Salinity-Induced variation in biochemical markers provides insight into the mechanisms of salt tolerance in common (Phaseolus vulgaris) and runner (P. coccineus) beans. Int. J. Mol. Sci. 2016, 17, 1582. [Google Scholar] [CrossRef]
  23. Morosan, M.; Al Hassan, M.; Naranjo, M.; Loprez-Gresa, M.B.M.; Vicente, O. Comparative analysis of drought responses in Phaseolus vulgaris (common bean) and P. coccineus (runner bean) cultivar. EuroBiotech J. 2017, 1, 247–252. [Google Scholar] [CrossRef]
  24. Arteaga, S.; Al Hassan, M.; Wijesinghe, C.; Yabor, L.; Llinares, J.; Boscaiu, M.; Vicente, O. Screening for Salt Tolerance in Four Local Varieties of Phaseolus lunatus from Spain. Agriculture 2018, 8, 201. [Google Scholar] [CrossRef]
  25. IBPGR. Phaseolus vulgaris Descriptors. In Plant Production and Protection Division; FAO: Rome, Italy, 1982; p. 32. [Google Scholar]
  26. Singh, S. A key for identification of different growth habits of Phaseolus vulgaris L. Annu. Rep. Bean Improv. Coop. 1982, 25, 92–94. [Google Scholar]
  27. Rodrigo-Miguez, A.P. Caracterización Morfoagronómica y Bioquímica de Judía Común (Phaseolus Vulgaris L.) de España. Ph.D. Thesis, Universidad de Santiago de Compostela, Santiago, Spain, 2000. [Google Scholar]
  28. Felsenstein, J. Numerical Taxonomy; W.H. Freeman & Company: San Francisco, CA, USA, 1973. [Google Scholar]
  29. Mohammadin, S.A.; Prasanna, B.M. Analysis of genetic diversity in crop plants—Salient statistical tools and considerations. Crop Sci. 2003, 43, 1235–1248. [Google Scholar] [CrossRef]
  30. Mantel, N. The detection of disease clustering and generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
  31. Gil, J.; De Ron, A. Variation in Phaseolus vulgaris in the Northwest of the Iberian Peninsula. Plant Breed. 1992, 109, 313–319. [Google Scholar] [CrossRef]
  32. Rodiño, P.; Santalla, M.; De Ron, A.; Singh, S. A core collection of common bean from the Iberian peninsula. Euphytica 2003, 131, 165–175. [Google Scholar] [CrossRef] [Green Version]
  33. Perez-Vega, E.; Campa, A.; De la Rosa, L.; Giraldez, R.; Ferreira, J. Genetic Diversity in a Core Collection Established from the Main Bean Genebank in Spain. Crop Sci. 2009, 49, 1377–1386. [Google Scholar] [CrossRef]
  34. Vaid, K.; Gupta, V.; Singh, R. Stability analysis in dry bean. Crop Improv. 1985, 12, 28–31. [Google Scholar]
  35. Escribano, M.; De Ron, A.; Amurrio, J. Diversity in agronomical traits in common bean populaltions from northwestern Spain. Euphytica 1994, 76, 1–6. [Google Scholar] [CrossRef]
  36. Nienhuis, J.; Singh, S. Combining ability analyses and relationships among yield, yield components and architectural traits in dry bean. Crop Sci. 1986, 26, 21–27. [Google Scholar] [CrossRef]
  37. Sills, G.; Nienhuis, J. Field plot technique affects snap bean yield evaluation. J. Am. Soc. Hortic. Sci. 1993, 188, 672–674. [Google Scholar] [CrossRef]
  38. Hosfield, G.; Uebersax, M.; Isleb, T. Seasonal and genotypic effects on yield and physico-chemical seed characteristics related to food quality n dry, edible beans. J. Am. Soc. Hort. Sci. 1984, 109, 182–189. [Google Scholar]
  39. Escribano, M.; Santalla, M.; Casquero, P.; De Ron, A. Patterns of genetic diversity in landraces of common bean (Phaseolus vulgaris L.) from Galicia. Plant Breed. 1998, 117, 49–56. [Google Scholar] [CrossRef]
Figure 1. Geographic origin of the 24 Phaseolus vulgaris cultivars.
Figure 1. Geographic origin of the 24 Phaseolus vulgaris cultivars.
Agriculture 09 00149 g001
Figure 2. Solar radiation registered in the greenhouse during the autumn (a) and spring trials (b).
Figure 2. Solar radiation registered in the greenhouse during the autumn (a) and spring trials (b).
Agriculture 09 00149 g002
Figure 3. Seed morphology of the 24 P. vulgaris accessions.
Figure 3. Seed morphology of the 24 P. vulgaris accessions.
Agriculture 09 00149 g003
Figure 4. Loading plot of the principal component analysis (PCA) conducted with the quantitative traits analyzed in 24 accessions of P. vulgaris. Abbreviations are as in Table 1.
Figure 4. Loading plot of the principal component analysis (PCA) conducted with the quantitative traits analyzed in 24 accessions of P. vulgaris. Abbreviations are as in Table 1.
Agriculture 09 00149 g004
Figure 5. Scatter plot of the PCA scores. In red, scores of spring trial; in blue, autumn trial. Yellow to greenish-yellow seeds (−), maroon seeds (◊), grey, brownish to greenish seeds (-), black seeds (□), brown, pale to dark seeds (Δ), whitish seeds (O). Accessions included in the yellow and green areas correspond to the ones included in the two main clusters in Figure 6.
Figure 5. Scatter plot of the PCA scores. In red, scores of spring trial; in blue, autumn trial. Yellow to greenish-yellow seeds (−), maroon seeds (◊), grey, brownish to greenish seeds (-), black seeds (□), brown, pale to dark seeds (Δ), whitish seeds (O). Accessions included in the yellow and green areas correspond to the ones included in the two main clusters in Figure 6.
Agriculture 09 00149 g005
Figure 6. UPGMA (unweighted pair group method with arithmetic mean) phenogram of the 19 accessions of Phaseolus cultivated in two consecutive seasons based on 17 quantitative descriptors. Phenetic relationships were based on Euclidean pairwise distances estimates. I and II represent the two main clusters obtained. In red, scores of spring season; in blue, autumn season.
Figure 6. UPGMA (unweighted pair group method with arithmetic mean) phenogram of the 19 accessions of Phaseolus cultivated in two consecutive seasons based on 17 quantitative descriptors. Phenetic relationships were based on Euclidean pairwise distances estimates. I and II represent the two main clusters obtained. In red, scores of spring season; in blue, autumn season.
Agriculture 09 00149 g006
Figure 7. Pod traits of the 24 P. vulgaris accessions evaluated.
Figure 7. Pod traits of the 24 P. vulgaris accessions evaluated.
Agriculture 09 00149 g007
Figure 8. Seed traits of the 24 P. vulgaris accessions evaluated.
Figure 8. Seed traits of the 24 P. vulgaris accessions evaluated.
Agriculture 09 00149 g008
Table 1. Quantitative and qualitative traits in the 24 Phaseolus vulgaris accessions. IBPGR—International Board for Plant Genetic Resources.
Table 1. Quantitative and qualitative traits in the 24 Phaseolus vulgaris accessions. IBPGR—International Board for Plant Genetic Resources.
Descriptor IBPGR [25]AcronymUnit/Scores
Quantitative
Vegetative
4.1.1. Leaflet length (LeafL)Measured on the terminal leaflet of third trifoliate leaf from pulvinus to leaf tip (cm)
Width of the central sheet ***(Wcshe)Measured on the terminal leaflet of the third trifoliate leaf (cm)
6.1.1. Hypocotyl length (HypL)Average length in cm, from 5 plants. Measured when primary leaves were fully expanded
6.1.9. Plant height (PlantH)Average height, in cm, at maturity, from 5 plants. Measured from cotyledon scar to tip of the plant
6.1.10. Stem diameter(StemD)Measured in mm at maturity for plants at crop density
Inflorescence and pod
4.2.2. Days to flowering(Dflow)Number of days from emergence to the stage where 50% of plants have set flowers
4.2.3. Flower buds per inflorescence(FlowBI)Average number of flower buds from 5 plants, examining one inflorescence from each plant
4.2.7. Pod length (PodL)Average length, in cm, of the largest fully expanded immature pod
6.2.9. Inflorescence length (InfloL)Average length, in cm, from 5 plants
6.2.10. Pedicel length (PedL)Average length, in cm, of older flower from 5 plants, examining one inflorescence from each plant
6.2.18. Pods per plant (PPlant)Average from 5 plants at crop density
Days to harvest ***(Dharv)Time period between seed sowing and mature pods sampling
Seed
6.3.1. Seeds per pod(SeedP)Average number of seeds from one pod, selected from 5 plants
6.3.3. Seed weight(SeedWh)Weight of 100 seeds, in grams
6.3.5.1. Length(SeedL)Measured parallel to the hilum (mm)
6.3.5.2. Width (SeedW)mm
6.3.5.3. Height (SeedH)Height, measured from hilum to opposite side (mm)
Qualitative
Vegetative
Plant type **(Planty)I. determinate bush, II. indeterminate straight, III. indeterminate crawling, IV. indeterminate climbing.
Inflorescence and Pod
4.2.6. Pod color (from the fully expanded immature pod) (Pcol)Dark purple (1), carmine red (2), purple stripe on green (3), carmine tripe on green (4), pale red stripe on green (5), dark pink (6), normal green (7), shiny green (8), dull green to silver gray (9), golden or yellow (10), pale yellow to white (11), other (12)
4.2.8. Pod cross-section(Pcross)Very flat (1), pear-shaped (2), round elliptic (3), figure of eight (4), other (5)
4.2.9. Pod curvature (of the fully expanded immature pod) (Pcurv)Straight (3), slightly curved (5), curved (7), recurving (9)
4.2.10. Pod suture string (Psuts)Stringless (0), few strings (3), moderately stringy (5), very stringy (7)
4.2.11. Pod color at physiological maturity (Pcmat)Dark purple (1), red (2), pink (3), yellow (4), pale yellow with colored mottling or stripes (5), persistent white (6)
6.2.12. Position of pods (PosP)Base (1), center (2), top (3), combination of 1, 2, and 3 (4), other (5)
6.2.15. Pod beak position (Pbeakp)Marginal (1), non-marginal (2)
6.2.15. Pod beak orientation (Pbeako)Upward (3), straight (5), downward (7)
6.2.17.2 Dry pod color(Dpcol)Dark purple (1), carmine red (2), purple stripe on green (3), carmine stripe on green (4), pale red stripe on green (5), dark pink (6), normal green (7), shiny green (8), dull green to silver gray (9), golden or deep yellow (10), pale yellow to white (11), other (12)
Seed
4.3.1. Seed coat patterns (Scoatp)Absent (0), constant mottled (1), striped (2), rhomboid spotted (3), speckled (4), circular mottling (5), marginal color pattern (6), broad-striped (7), bicolor (8), spotted bicolor (9), pattern around hilum (10), other (11)
4.3.2. Seed coat color (Scoatdc)Black (1), brown, pale to dark (2), maroon (3), gray, brownish to greenish (4), yellow to greenish yellow (5), pale-cream to buff (6), pure white (7), whitish (8), white, purple-tinged (9), chlorophyll green (10). green to olive (11), red (12), pink (13), purple (14), other (15)
4.3.4. Brilliance of seed (BSeed)Matte (3), medium (5), shiny (7)
4.3.5. Seed shape(SHseed)Round (1), oval (2), cuboid (3), kidney-shaped (4), kidney pronounced (5), fastigiate trunk (6)
6.3.2. Apparent seed veining(Aseed)Absent (0), Present (1)
** Singh [26], *** Rodrigo [27].
Table 2. Two-way analysis of variance (ANOVA) of cultivar and cycle effects and their interactions for the parameters considered. Numbers represent F-values at the 5% confidence level.
Table 2. Two-way analysis of variance (ANOVA) of cultivar and cycle effects and their interactions for the parameters considered. Numbers represent F-values at the 5% confidence level.
Dependent VariableIndependent Variable
CultivarCycleCultivar × Cycle
LeafL (Leaflet length)36.42 ***0.468 ns0.79 ns
Wcshe (Width of the central sheet)46.42 ***1.479 ns0.97 ns
HypL (Hypocotyl length)135.54 ***0.480 ns28.52 ***
PlantH (Plant height)492.86 ***3.24 ns0.28 ns
StemD (Stem diameter)40.29 ***0.055 ns9.85 ***
FlowBI (Flower buds per inflorescence)31.07 ***0.09 ns2.57 ***
PodL (Pedicel length)172.61 ***2.05 ns3.40 ***
InfloL (Inflorescence length)85.31 ***14.21 ***32.06 ***
PedL (Pedicel length)70.89 ***1.90 ns15.68 ***
Pplant (Pods per plant)222.09 ***3.92 *0.58 ns
SeedP (Seed per pod)31.24 ***7.09 **1.07 ns
SeedL (Length)237.020.423.76 ***
SeedW (Width)73.49 ***1.14 ns7.47 ***
SeedH (Height)71.75 ***6.35 *3.93 ***
Yield 343.70 ***2.27 ns4.22 ***
*p < 0.05; ** p < 0.01; *** p < 0.001; ns: not significant.
Table 3. Correlation coefficients between the first three principal components (PCs) and the morphological descriptors.
Table 3. Correlation coefficients between the first three principal components (PCs) and the morphological descriptors.
Component 1Component 2Component 3
LeafL0.090.39−0.08
Wcshe0.080.38−0.14
HypL0.090.040.17
PlantH−0.17−0.190.18
StemD0.090.20−0.22
Dflow0.110.050.12
FlowBI0.130.14−0.07
PodL0.220.070.26
InfloL0.13−0.21−0.41
PedL0.14−0.19−0.39
Pplant0.220.270.14
Dharv0.10−0.02−0.35
SeedP0.030.350.19
SeedWh0.36−0.250.04
SeedL0.33−0.180.28
SeedW0.33−0.15−0.13
SeedH0.35−0.220.02
Scoatdc−0.22−0.240.22
SHseed0.20−0.230.23
Aseed−0.29−0.120.15
Yield0.350.180.23
Table 4. Descriptive statistics for quantitative traits in each of the established groups and the whole set of accessions (number of accessions). Autumn growing season. Units of quantitative traits are as defined in Table 1. CV—coefficient of variation.
Table 4. Descriptive statistics for quantitative traits in each of the established groups and the whole set of accessions (number of accessions). Autumn growing season. Units of quantitative traits are as defined in Table 1. CV—coefficient of variation.
Accessions HypLPlantHStemDLeafLWcsheDflowDharvFlowBIInfloLPedLPplantPodLSeedPSeedLSeedWhSeedHSeedWYield
All accesionsMean10.531.705.1816.5111.4842.54110.254.177.653.7729.1514.265.0912.9647.147.775.9474.53
SD3.050.691.263.763.226.6523.081.422.421.5517.314.281.873.0717.291.201.1461.41
CV28.9840.3924.3522.7428.0215.6420.9434.0431.6941.1659.3830.0436.7823.6836.6815.4919.2482.39
Range2.800.322.759.805.4028.0071.001.002.711.105.005.601.007.0317.204.723.642.80
19.002.958.4026.5020.2055.00157.008.0014.608.1072.0028.509.0021.3981.4011.498.71242.75
Black (3)Mean10.731.855.5020.2615.8441.00117.334.058.294.2343.4712.666.1710.4533.897.405.7692.49
SD1.230.320.703.532.515.5729.691.431.882.0616.973.362.201.7415.671.011.1962.28
CV11.5017.4412.7717.4415.8313.5825.3035.3522.6748.6439.0326.5335.6516.6446.2513.5920.7167.33
Range9.10.994.21611369225.91.923817.519.25.624.3834.52
12.62.46.526.520.247150612.18.17218.5912.7750.399.778.71193.95
Brown (8)Mean10.672.285.8816.0510.8943.83114.173.986.973.6339.6613.945.6413.9160.097.955.96124.67
SD2.700.381.413.482.438.8622.991.391.461.4817.693.761.682.9517.821.471.3457.97
CV25.3016.8223.9921.6822.2820.2220.1434.8820.9840.8444.6226.9929.7521.2329.6518.4722.5646.50
Range4.51.752.99.87.8359413.81.155.629.4937.4553.6418.97
14.52.958.422.419.255143610.976521.2818.6981.410.748.16242.75
Whitish (10)Mean11.641.804.5216.3211.5545.29112.143.896.573.2026.3015.925.0212.5239.567.445.5254.07
SD2.930.200.842.922.782.6921.071.492.891.3210.515.251.892.8412.861.060.9238.14
CV25.2111.3418.4717.9024.055.9418.7938.2043.9541.3439.9632.9937.7122.7032.5014.2816.7170.54
Range7.81.22.7511.17.2418722.711.189.827.41194.723.8915.14
192.47.832519.850141814.67.84728.5918.0953.28.988.01146.83
Yellow to greenish yellow (1)Mean8.561.794.3314.049.444698.004.2511.785.2515.2014.225.0011.3453.208.697.7040.43
SD0.740.020.693.481.23--1.711.871.373.491.200.710.58-1.150.389.29
CV8.651.3416.0324.8013.08--40.1815.9126.1422.988.4714.145.16-13.234.9022.98
Range7.81.753.6710.17.34698293.21112.2410.1653.27.226.9929.26
9.51.815.32010.946986136.12015.2612.0653.211.498.1153.2
Gray, brownish to greenish (1)Mean6.181.806.0417.6813.5037.00157.004.337.674.0512.0013.954.5012.2338.607.905.6220.84
SD1.80-1.411.151.09--1.210.820.613.390.900.840.90-0.320.345.89
CV29.07-23.436.538.10--27.9510.6515.0428.266.4218.597.37-4.075.9928.26
Range4.81.83.621611.937157373.4813310.9238.67.395.313.89
9.21.87.131914.6371576951615.3513.7438.68.336.2627.79
Maroon (1)Mean13.662.264.8419.2914.6751.00110.003.806.482.3451.2017.425.6015.1063.708.567.46182.64
SD0.550.190.441.512.93--0.450.380.357.013.100.890.59-0.290.3925.02
CV4.008.259.187.8319.97--11.775.9214.9913.7017.8215.973.91-3.425.2313.70
Range12.82.084.33179.5511103624412.1414.1463.78.096.77156.95
14.22.55.3621.217.85111046.92.96219.8616.1663.78.998.14221.16
Determinate growth habit
Whithish (3)Mean6.970.544.7615.709.9841.00102.004.638.744.9713.3610.363.4112.0638.337.586.0019.44
SD3.020.130.514.772.593.614.581.452.781.584.072.921.083.5420.861.300.9413.69
CV43.2523.3310.7630.3625.908.794.4931.4631.7631.8830.4628.1631.8029.3554.4117.1915.6270.42
Range2.80.3354.29.85.4389824.232.866.527.0317.25.224.682.80
9.80.75.87261445107612.57.12215.1517.1558.99.567.4947.94
Brown (2)Mean12.220.425.4313.707.5830.5076.505.239.303.6613.2016.834.0716.8560.908.376.0132.50
SD1.910.082.121.900.763.547.781.171.510.734.211.491.141.662.400.810.6210.76
CV19.4916.3527.1612.779.5510.719.4919.4316.5816.5060.169.2428.5310.844.0610.1110.7163.35
Range90.3153.410.46.5287137.82.8714.5114.5159.27.114.6416.99
150.558.116.1933826124.82119.2521.3962.69.677.0450.97
Table 5. Descriptive statistics for quantitative traits in each of the established groups and the whole set of accessions (number of accessions). Spring growing season Units of quantitative traits are as defined in Table 1.
Table 5. Descriptive statistics for quantitative traits in each of the established groups and the whole set of accessions (number of accessions). Spring growing season Units of quantitative traits are as defined in Table 1.
Accessions HypLPlantHStemDLeafLWcsheDflowDharvFlowBIInfloLPedLPplantPodLSeedPSeedLSeedWhSeedHSeedWYield
All accessionsMean10.591.695.1317.0312.1142.25108.54.187.383.5730.1414.385.3412.9146.157.655.8976.06
SD3.030.661.063.873.268.4121.861.512.971.6218.023.981.762.9716.561.040.9961.71
CV28.6539.0220.7322.7526.8819.9220.1536.1540.245.3759.7827.6733.0223.0135.8913.616.7781.13
Range20.253.149.86.5277011.50.9846.816.5916.985.414.036.99
17.82.987.8726.520.555153925.5412.227727.8921.1776.2310.159.38232.82
Black (3)Mean12.151.885.0420.3115.8242.33121.334.527.473.4144.412.836.6710.5235.297.235.68108.08
SD2.950.180.973.441.686.4328.181.361.921.5218.714.311.832.0417.41.011.1375.21
CV24.289.3919.3116.9610.615.1923.2330.1625.744.6142.1333.6227.3919.4349.313.9119.8369.58
Range8.91.73.916.213.6359925.52.121847.2418.3745.764.133.56
17.82.47.8726.519.5471536127.87719.4913.3253.138.997.78232.82
Brown (8)Mean11.932.275.3416.9111.8444113.334.057.353.5641.7314.536.2313.5454.678.026.25123.87
SD1.30.411.193.92.789.0125.081.693.731.7818.743.421.552.4516.181.080.9157.58
CV10.9218.1522.2723.0423.520.4822.1341.8250.7450.0644.9123.5124.8518.1229.6113.4114.5346.48
Range9.81.63.359.87.8349123.81.946.82936.185.464.2718.688
14.52.987.8324.519.254147925.5412.226521.2817.776.239.659.38217.47
Whitish (10)Mean10.631.814.8815.9511.7943.86111.293.396.042.8726.0915.55.112.5639.397.435.4351.36
SD2.650.190.883.333.028.1320.851.272.031.0610.444.881.812.5311.380.980.9230.47
CV24.9210.6917.9520.8725.5818.5418.7437.5133.736.9640.0431.4635.5820.1428.8913.1916.9559.33
Range5.41.313.1411.17.5278513.20.9897.928.0119.5345.414.0317.17
15.62.57.825.220.553138614.364527.8817.449.19.367.9128.79
Yellow to greenish yellow Mean8.221.84.88149.5549984.411.484.2816.213.084.810.4956.17.857.2743.62
SD0.830.050.593.351.24--1.521.430.475.890.271.10.67-0.540.4515.86
CV10.082.5312.0723.9212.99--34.4712.4311.0136.362.0522.826.38-6.836.1536.36
Range7.21.763.9410.17.34998293.9912.839.6456.16.786.6924.23
9.11.885.452010.94998612.55.12313.5611.6356.18.437.9161.93
Gray, brownish to greenish Mean7.321.86.817.313.56401254.68.584.412.814.964.611.9337.056.865.2221.82
SD1.70.090.461.370.88--1.141.061.011.920.470.520.82-0.370.353.28
CV23.255.066.747.956.53--24.7912.422.915.033.1611.236.85-5.396.6115.03
Range 5.61.656.315.812.54012536.93.21014.5410.7137.056.084.4817.04
9.21.897.519.214.64012569.55.51515.6512.9137.057.395.6525.56
MaroonMean9.62.126.5519.515.1551104.336.42.655.217.685.6714.1859.018.156.82184.58
SD2.540.030.921.493.05--1.860.230.46.611.680.820.79-0.450.6222.11
CV26.431.36147.6620.21--42.973.5615.3811.989.4914.415.54-5.499.1411.98
Range 4.82.15.53179.55511036.22.14814.8412.7559.017.265.6160.50
11.52.157.3121.217.85511086.83.16319.1615.2759.018.587.68210.66
Determinate growth habit
Whitish (3)Mean7.280.544.8117.5911.033998.6758.114.9914.810.963.8912.4434.87.365.4821
SD3.640.150.784.591.813.617.091.073.662.352.932.830.893.6415.991.040.511.64
CV49.9328.1816.226.1116.399.257.1921.3845.1147.1119.8125.8722.9329.2745.9614.149.255.41
Range20.2513.82128359131.52.110726.5916.985.54.836.98
13.20.676.0526144210561510.12115.5517.247.918.96.7142.25
Brown (2)Mean12.610.414.8713.737.4227.5765.229.044.3114.315.264.2817.5970.698.486.5242.72
SD2.20.081.212.020.630.718.491.093.291.262.981.411.021.747.780.90.587.82
CV17.4719.2124.7714.738.542.5711.1620.9336.3629.1720.869.2423.799.891110.588.8318.31
Range70.2853.4910.46.5277035.12.3813.5114.8965.1897.155.5225.14
15.10.5317.1916.182882613.16.21917.2521.1776.1910.157.4954.49

Share and Cite

MDPI and ACS Style

Arteaga, S.; Yabor, L.; Torres, J.; Solbes, E.; Muñoz, E.; Díez, M.J.; Vicente, O.; Boscaiu, M. Morphological and Agronomic Characterization of Spanish Landraces of Phaseolus vulgaris L. Agriculture 2019, 9, 149. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9070149

AMA Style

Arteaga S, Yabor L, Torres J, Solbes E, Muñoz E, Díez MJ, Vicente O, Boscaiu M. Morphological and Agronomic Characterization of Spanish Landraces of Phaseolus vulgaris L. Agriculture. 2019; 9(7):149. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9070149

Chicago/Turabian Style

Arteaga, Sugenith, Lourdes Yabor, José Torres, Eva Solbes, Enrique Muñoz, María José Díez, Oscar Vicente, and Monica Boscaiu. 2019. "Morphological and Agronomic Characterization of Spanish Landraces of Phaseolus vulgaris L." Agriculture 9, no. 7: 149. https://0-doi-org.brum.beds.ac.uk/10.3390/agriculture9070149

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