Original Research ARTICLE
The Spanish Core Collection of Common Beans (Phaseolus vulgaris L.): An Important Source of Variability for Breeding Chemical Composition
- 1Miquel Agustí Foundation, Barcelona, Spain
- 2Nestle Purina, St. Louis, MI, United States
- 3Department of Agri-Food Engineering and Biotechnology, BarcelonaTech, Universitat Politecnica de Catalunya, Barcelona, Spain
The Iberian Peninsula is considered as a secondary center of diversity for the common bean, and the Spanish National Plant Genetic Resources Centre’s germplasm bank holds more than 3,000 Spanish accessions of Phaseolus vulgaris L. from which a core collection of 202 landraces has been selected. In order to encourage the use of this abundant resource, this study aimed to characterize genetic diversity, by measuring chemical composition in these core collections (in both the seed coat and cotyledon) using previously developed near infrared spectroscopy models. Crucially, these landraces in question all originated under similar agroclimatic conditions, allowing these field trials to be conducted in a single location without significantly altering the agronomic behavior of individual accessions. Using previously reported data, we also explored the correlations between chemical composition and culinary/sensory traits, as well as possible associations between chemical composition and seed coat color or gene pool (Middle American or Andean). The general Mahalanobis distance was >3 in only 11 of 1,950 estimations, confirming the robustness of the regression models previously developed. Variability was greater in seed coat than in cotyledon compounds and ranges for all compounds were wide: ash 34–94 g/kg, Ca 5–31 g/kg, dietary fiber 554–911 g/kg, Mg 2–4.4 g/kg, uronic acid 95–155 g/kg, protein 192–304 g/kg, starch 339–446 g/kg, amylose 208–291 g/kg, amylopectin 333–482 g/kg, and apparent amylose 241–332 g/kg. Accessions with white seed coats tended to be richer in ash, dietary fiber, uronic acid, and Ca, and accessions of the Middle American gene pool had on average 65% more Ca than the Andean gene pool. Strong genetic correlations were not identified between chemical and culinary/sensory traits. This is particularly positive with regards to plant breeding, as it means that synchronic improvement of nutritional composition and sensory traits is possible. The genetic diversity of chemical composition described in the Spanish core collection of beans therefore represents a promising opportunity to develop cultivars with superior nutritional profiles.
Dry beans (Phaseolus vulgaris L.) are an extremely important aliment, not only representing the main source of dietary protein for humans in several world regions but also contributing greatly to diet with starch, fiber, vitamins, and minerals (Hayat et al., 2014). Annual global production is presently approximately 26.5 million tons, most of which is used for human consumption (FAOSTAT, 2018). Although in the most developed countries, the importance of dried bean consumption has diminished with the increase of meat consumption, producing animal protein is far more expensive, and is generally considered unsustainable (Pimentel and Pimentel, 2003; Gonzalez et al., 2011). Red meat in particular also has negative effects on health, for example being associated with the development of cardiovascular disease and cancer (Micha et al., 2010; Pan et al., 2013; Gonzales et al., 2014; Ekmekcioglu et al., 2018). International organizations and specialists therefore recommend increasing the consumption of beans and other legumes to fulfill our nutritional needs and decrease inputs of food production (Leterme, 2002; WHO, 2015; FAO, 2016; De Ron et al., 2017; McDermott and Wyatt, 2017). To promote this change in diet, apart from optimizing crop yields, we also need to identify palatable varieties of beans and legumes that are both tied to our gastronomic cultures and provide the maximum amount of nutrients (especially protein). Breeding programs for beans have focused mainly on maximizing yield (Kelly et al., 1998; Singh, 1999) and improving resistance to both biotic and abiotic stresses (Ishitani et al., 2004; Miklas et al., 2005; Singh and Schwartz, 2009; Araújo et al., 2015). More recently, however, nutritional and sensory traits have also been added to the ideotypes of these crops (Vaz Patto et al., 2015).
Natural genetic variation for chemical composition in the cultivated genepool and wild crop relatives of Phaseolus vulgaris is particularly important as it offers us the opportunity to develop cultivars with superior nutritional profiles. For example, reported percentages of protein content range from 18 to 31%, from 50 to 76% for carbohydrates, and from 0.05 to 0.31% for Calcium (Ca) (Sathe et al., 1984a; Islam et al., 2002; Pinheiro et al., 2010; Hayat et al., 2014). This genetic diversity can fortunately be found conserved in agrosystems located in the centers of diversity for the species (in situ conservation) and in genebanks (ex situ conservation). Of more than 7.4 million accessions stored in approximately 1,750 genebanks (FAO, 2010), 260,000 belong to Phaseolus vulgaris L. Although significant efforts have been devoted to characterizing such germplasm collections for simple traits (mostly related to botanical aspects) (Wang et al., 2017) and for the most important agronomic traits (yield and resistance to pests and diseases) (Tanksley and McCouch, 1997), less is known about other important characteristics such as those related to nutritional, culinary, or sensory attributes (Smýkal et al., 2015). This is mainly due to the fact that these traits are quantitative and multigenic, with low heritability and strong genotype by environment interactions (GxE), thus making it difficult to assign a genotypic value for each accession. Moreover, the phenotyping of these traits is expensive in terms of both direct monetary input and human labor. However, the lack of phenotyping data available, regarding the accessions conserved ex situ, hinders utilization of this genetic diversity both in breeding programs and directly by farmers (Ramanatha Rao and Hodgkin, 2002; Hodgkin et al., 2003). Novel “omics” platforms enabling the massive analysis of the genome, transcriptome, proteome, and metabolome (Tanksley and McCouch, 1997; Prada, 2009; Langridge and Fleury, 2011), and information science (Michael et al., 2018), represent promising tools to find desirable agronomic alleles in seed banks and unlock the stored genetic diversity. The development of core collections that are representative of the genetic spectrum in the whole collection (van Hintum et al., 2000) has been proposed as an alternative to enable cost-effective characterization of the plant genetic resources held in genebanks. Initially, complex traits are phenotyped in the core collection, and if desirable qualities are found, similar accessions can then be identified in the whole collection (Prada, 2009). An example of this strategy is the Spanish core collection of beans of the Spanish National Plant Genetic Resources Center (CRF), comprising 202 accessions representative of the variability present in the total set of accessions collected in Spain. This collection was compiled based on location data (province, township, and altitude) and seed phenotype (color, shape, and size) (De la Rosa et al., 2000) and is an important resource because the Iberian Peninsula is considered a secondary center of diversity for Phaseolus vulgaris L. (Santalla et al., 2002). To date, several traits have been measured in this collection: flower, pod, and seed traits; growth habit; type of phaseolin and molecular markers (Pérez-Vega et al., 2009); resistance to pests and diseases [including halo blight (Pseudomonas syringae), common bacterial blight (Xanthomonas campestris) (Asensio et al., 2010), anthracnose (Colletotrichum lindemuthianum) (Pérez-Vega et al., 2006), and white mold (Sclerotinia sclerotiorum) (Pascual et al., 2010)]; and sensory and culinary traits (Rivera et al., 2016). However, nothing was known about the variability in chemical composition and nutritional potential of this collection. The main factors generally limiting substantial chemical composition studies and comparison of materials conserved in genebanks are the difficulties involved in analyzing so many samples. Nowadays, however, instrumental methods such as near infrared spectroscopy (NIRS) or nuclear magnetic resonance spectroscopy (Hacisalihoglu et al., 2010; Parlak and Güzeler, 2016) make it possible to analyze large collections of samples. With this objective, our team has developed NIRS predictive models for chemical composition and sensory traits (Plans et al., 2012, 2013, 2014), which show a robust capacity for accurate prediction.
The present study extends our previous work about the culinary and sensory traits of the accessions in the CRF’s core collection of common beans (Rivera et al., 2016). Here, we examine the collection to (i) describe the variability of chemical composition in the seed coat and cotyledon, (ii) compare and correlate this variability and its magnitude with the variability in sensory and culinary traits described in the previous study, and (iii) analyze the possible associations between chemical composition and seed coat color and the Middle American or Andean gene pool.
Materials and Methods
Plant Material and Field Trials
The 202 accessions from the Spanish core collection of common beans form a rich mosaic of colors and shapes (Figure 1), encompass 51 market classes (Santalla et al., 2001), and represent all areas of Spain where beans are cultivated. All of the seeds were sown in Sabadell (Northeast Spain: 41° 32′ 50.7″ N, 2° 4′ 14.7″ E), a location with loam soil with abundant Ca (6.35 mg/kg), low phosphorus (6 mg/kg), and a mild Mediterranean climate, which allows both short and long-cycle accessions to develop to maturity. We used a randomized two block design with 21 plants per block and accession (total 42 plants per accession). The experiment was conducted at a low density (29,167 plants/ha) to facilitate the individualized recording of data. A vertical plastic net was used to trellis the accessions with indeterminate growth. All plants were cultivated with the traditional management practices of the area, including drip irrigation and fert-irrigation application (NPK 110 kg/ha). Due to the lack of significance of the block factor for agronomic traits (data not shown), seeds from the two blocks were pooled in order to ensure sufficient sample was available for the analyses. Seeds from a total of 195 accessions were harvested and processed for analyses (seven accessions did not produce enough seeds for analyses).
FIGURE 1. Photograph showing the variety of shapes and colors represented in the Spanish core collection of beans.
The seed coat accounts for less than 10% of the total dry matter of bean seeds and has a chemical composition distinctly different to that of the cotyledon (Singh et al., 1968; Moraghan et al., 2006). To ensure accurate information about the two fractions, seed coat, and cotyledons were analyzed separately. Seeds from each harvested accession (approximately 50 g) were soaked in 150 ml distilled water for 24 h, drained, dried with filter paper, and weighed. The seed coat was then manually separated from the cotyledon, and the two fractions were dried (60°C for 48 h). Finally, dried samples were ground to 0.4 mm in a laboratory mill (Perten 3100, Perten Instruments Inc., Springfeld, IL, United States). Ground samples were stored in polyethylene bags at 4°C until spectroscopic analysis was conducted.
NIRS Recording to Estimate Chemical Composition
Near infrared spectroscopy is a well-established technique for determining the components of foods (Nicolaï et al., 2007). Models developed by means of multivariate analysis correlating physicochemical properties, and the spectra obtained with NIRS technology enable prediction of the value of a sample. In our study, we used NIRS predictive models previously developed by Plans et al. (2012), (2013), to estimate the following chemical components: ashes, Ca, magnesium (Mg), dietary fiber, and uronic acids in the milled seed coat; and protein, starch, amylose, amylopectin, and apparent amylose in the milled cotyledon. These traits were selected because of their relationship with sensory attributes driving consumer preferences, as reported in previous studies (Casañas et al., 2002, 2006; Mkanda et al., 2007). All results were expressed as g/kg dry matter. The estimations were calculated using regression models from NIRS measurements. Infrared spectra from the ground samples were recorded by a spectrophotometer (model 5000, Foss NIRSystems, Silver Spring, MD, United States) at every 2 nm between 1,100 and 2,500 nm, with a mean of 32 scans performed. Each sample was analyzed per triplicate, and the mean spectrum reading was used for calculations. The absorbance of each wavelength was transformed into Log (1/R), where R represents reflectance, due to this variable correlating better with chemical components than raw reflectance. To record spectra and import data, VISION software was used (version 2.51, Foss NIRSystems, Silver Spring, MD, United States).
Twenty out of the accessions under study were randomly selected in order to perform chemical analysis. For each accession, three biological replicates were characterized for the chemical components. Ash, Ca, dietary fiber, Mg, and uronic acid were quantified by reference methods following the methodology described in Plans et al. (2012), while protein, starch, amylose, amylopectin, and apparent amylose were analyzed as described in Plans et al. (2013). These reference analyses were used to increase the domain of the chemical constituents of the NIRS model, thus improving its robustness as described by Barton et al. (2000).
Calculating Chemical Composition
Two different approaches can be used to correlate NIRS data from samples with the constituents of interest (in this case, chemical composition): global models based on partial least squares regression (PLSR) and principal components regression (PCR) or local models based on the similarity of the unknown spectrum to the known spectra in a database. Classical regression techniques such as PLSR and PCR assume that the relationship between spectral and reference constituents is linear (Shenk and Westerhaus, 1991). The accuracy of NIRS in predicting chemical constituents is related to the structure and distribution of spectra in the calibration set, which must cover most of the possible variability in spectra encountered during routine analysis. For global PLSR models, increasing the domain of the calibration equation to include new diverse geographical and climatic areas, will require a larger number of representative samples in the calibration set. This, in turn, increases the complexity of the spectral reference model, making it necessary to compute more parameters to explain the variation in the constituents (Martens and Naes, 1989). The accuracy of prediction using global PLSR models usually decreases when the domain of the constituent increases.
Local calibration, also called memory-based regression, aims to overcome these difficulties. Since, in theory, an optimal prediction should be computed based on a specific calibration equation, we can use a small calibration dataset tailored to the unknown sample from a large database of spectra. This approach combines the advantage of global calibration (to cover a large constituent domain) and the accuracy obtainable with specific calibrations based on spectral similarities of a smaller set of samples. Various approaches can be used to calculate local models, with the most common algorithms being the classic similarities based on Mahalanobis distances (H) (Naes et al., 1990) and correlation coefficients (Shenk et al., 1997). Several authors found that local regression models significantly improved (20–30%) the standard error of prediction (SEP) compared to global calibration (Naes et al., 1990; Aastveit and Marum, 1993; Sinnaeve et al., 1994) in studies using sample selection and calibration methods based on PCR and PLSR.
We used global PLSR-derived correlations between NIR spectra and reference constituents (cotyledon: protein, starch, amylopectin, amylose, and apparent amylose; and seed coat: ash, calcium, dietary fiber, magnesium, and uronic acid). These models were previously developed using a wide range of genetic and environmental variability (with cultivation in many locations) (Plans et al., 2012, 2013). Furthermore, we increased the domain of the chemical constituents by using a set of 20 new accessions from the Spanish core collection of common beans. To avoid the problems inherent in using PLSR with increased domains, we used local regression (Ramirez-Lopez and Stevens, 2016) to correlate the spectra with the chemical referents. Local models were not developed using a single PLSR or PCR equation; rather a new calibration was created when an unknown sample was presented against the spectral database. The ability to predict unknown samples was based on global Mahalanobis distances (GH). This distance is normalized when measured in standard deviation units from the center of the selected set of samples. Based on a normal distribution, samples with a GH of greater than three would be classified as outliers and predictions in these samples would be considered questionable (Shenk et al., 1997).
To estimate the global variability for each trait, and the relationships between traits, we calculated their mean, and variance. Previous results on culinary, sensory, and external appearance traits obtained from the same samples (Rivera et al., 2016) were used to calculate the correlation with the results obtained in this study. Sensory traits were evaluated using regression models from NIRS measurements in ground cooked samples. These models were previously developed by using evaluations from a trained panel (Plans et al., 2014). By contrast, culinary traits were quantified directly by using a standardized cooking process, as described in Romero del Castillo et al. (2012). Correlation analysis was performed using the Pearson coefficient with the Bonferroni correction. For further analysis of the variability, the accessions were classified according to their seed coat color and origin gene pool (Middle American or Andean). The data concerning seed coat color were reported by Rivera et al. (2016). Seeds were visually evaluated and divided into two categories: white and colored (including: yellow, cream, gray brownish to greenish, brown, vinous brown, black, ochre, purple, rosy, bicolor, and tricolor). The origin gene pool was previously described by PérezVega et al. (2009), following study of the phaseolin protein pattern and 11 molecular markers [two sequence characterized amplified regions (SCAR) and nine simple sequence repeats (SSR)] in the collection. For each type of classification, analysis of variance (ANOVA) was carried out based on the linear model yij = μ + αi + 𝜀ij, where yi was an individual level for a specific trait; μ was the grand mean; αi was the effect of i group based on the type of seed coat color, or gene pool origin; and 𝜀ij was the random error for i groups with j replications of the model following a N∼(0, σ2). All factors were considered to be fixed. Finally, for seed coat color and gene pool classification, we performed normalized principal component analyses (PCA) and calculated 95% confidence ellipses around each cluster of accessions with the same category. We used R software (R Core Team, 2017); agicolae, PCAmethods, and ellipse packages for all statistical analyses (Murdoch and Chow, 2007; Stacklies et al., 2007; Mendiburu, 2017).
Overall Variability in the Collection
The GH values obtained were extremely low and reflect the reliability of our estimations using the NIRs models. Only 11 from a total of 1950 estimations presented a GH value greater than 3. These values corresponded to amylopectin, and in each case were not considered in the subsequent computations. Furthermore, the mean GH values for individual traits all scored below GH < 1 (except amylopectin; GH < 1.3), thus highlighting the robustness of our data regarding the different chemical components in both the seed coat and cotyledon (Shenk et al., 1997).
The proportion of seed coat ranged from 4% to 10% (coefficient of variation, CV = 11.45%), and the proportion of cotyledon ranged from 90% to 96% (CV = 0.91%). For the chemical traits measured in the cotyledon, the CV ranged from 5.35% for starch to 10.97% for proteins. In general, CV from the seed coat components measured is higher than the CV from the cotyledon. In the seed coat, CV ranged from 8.31% for uronic acids to 39.07% for Ca (Table 1).
TABLE 1. Variability for the chemical compounds analyzed in the cotyledon and in the seed coat in the Spanish bean core collection of beans.
Mean phenotypic values were used to estimate genotypic correlations among chemical compounds (Table 2). In the cotyledon components, negative correlations were identified between protein and starch (-0.60, p < 0.001), amylose (-0.37, p < 0.001), and apparent amylose (-0.34, p < 0.01). Positive correlations were identified between many other traits, the most significant being between amylose and apparent amylose (0.69, p < 0.001), and between starch and amylose (0.5, p < 0.001). With regard to the seed coat composition, several significant correlations were found, including between Ca and ashes (0.91, p < 0.001), uronic acid and dietary fiber (0.42, p < 0.001), and uronic acid and ashes (0.35, p < 0.01). In relation to correlations between seed coat and cotyledon compounds, no correlations were identified. Considering the wide range of genetic diversity explored, we also performed a correlation analysis between the chemical results obtained in this work and both sensory and culinary traits reported previously in the same samples (Rivera et al., 2016). This analysis included 174 out of the original 195 accessions, because for some accessions there were not enough available seeds to perform the sensory and culinary analysis, for which a large sample is needed. Several significant correlations were identified (Table 2). With regard to sensory traits, the most significant correlations were found between mealiness and the seed coat measurements of ashes, Ca, dietary fiber, and uronic acid (-0.62, -0.57, -0.42, and -0.47, p < 0.001, respectively), and between seed coat brightness and ashes, Ca, dietary fiber, and uronic acid (0.52, 047, 0.43, and 0.46, p < 0.001, respectively). Fewer significant correlations were found with the culinary traits, for instance between rate of water absorption and ashes, Ca, dietary fiber, and uronic acid (0.39, 0.36, p < 0.001; 0.33 p < 0.01; and 0.3, p < 0.05, respectively), and between cooking time and protein (0.33, p < 0.01). Percentage of white surface color was positively correlated with all chemical components measured in the seed coat, while 100 seed weight correlated only with Ca (-0.37, p < 0.001) and ashes (-0.34, p < 0.01) (Table 2).
TABLE 2. Genotypic correlations (Pearson coefficient with Bonferroni correction) between chemical composition traits measured in the cotyledon and seed coat, and sensory and culinary traits measured in a previous study (Rivera et al., 2016) using the same samples.
Effect of Seed Color on Chemical Variability
Seed color description, as reported in Rivera et al. (2016), revealed that 38 out of the 174 accessions studied had white seeds, and 136 accessions presented colored seeds. An ANOVA comparing the seed color classification (white or colored) was performed for all of the chemical traits measured (Table 3). Significant differences were found between white and colored seeds for all of the traits measured in the seed coat (ash, Ca, dietary fiber, Mg, and uronic acid), while solely for starch content regarding the traits measured in the cotyledon. White seeds tended to have higher concentrations for all of the traits measured in the seed coat and had lower starch content. The PCA created, considering all of the chemical traits evaluated in the seed coat and cotyledon with the different representation of the white and colored accessions, showed that the distribution of the accessions along the first (PC1) and second (PC2) principal components together explain 49% of the total variation (Figure 2). PCA showed that ash, Ca, uronic acids, and dietary fiber, which were correlated with PC2, had the higher contribution for clustering between groups, while the other chemical compounds had a low influence. Nevertheless, the confidence intervals for the two color groups overlapped, highlighting that within the colored group, there are genotypes with high values for the chemical traits analyzed.
TABLE 3. Analysis of variance to compare the levels of the different chemical components measured in the cotyledon and seed coat according to seed coat color classification proposed in Rivera et al. (2016).
FIGURE 2. Principal component analysis based on chemical composition of the seed coat and cotyledon. Accessions with white seeds are represented with green symbols and accessions with colored seeds are represented with red symbols. The dashed lines represent the 95% confidence ellipses.
Effect of Andean or Middle American Gene Pool on Chemical Variability
Previously, PérezVega et al. (2009) had used phaseolin and molecular markers to classify the Spanish core collection according to their origin gene pool. Forty three out of the 174 accessions studied were classified within the Middle American gene pool, and 131 within the Andean gene pool. Using this classification, we performed an ANOVA, which revealed significant differences among gene pools for the following components: in the cotyledon, protein, starch (p < 0.001), amylose, and apparent amylose (p < 0.01); in the seed coat, ashes, and Ca (p < 0.001) (Table 4). Accessions classified in the Middle American gene pool tended to have higher ash, Ca, and protein contents, with lower levels of starch, amylose, and apparent amylose. The greatest difference was found for the content of Ca in the seed coat, with Middle American samples yielding a content 65% higher than Andean materials (Table 4). The same PCA as Figure 2 performed considering all of the chemical traits measured in the seed coat and cotyledon, but with the different representation of the accessions of Middle American and Andean gene pool showed that ash and Ca, which were correlated with PC2, had the higher contribution for clustering between groups (Figure 3). This multivariate analysis highlights that within the Andean gene pool there are genotypes with high values for uronic acids and dietary fiber contents (11 and 14 out of the 20 higher scores for uronic acids and dietary fiber, respectively, belong to accessions from Andean gene pool).
TABLE 4. Analysis of variance comparing the content of the different analytes studied according to the gene pool classification reported by PérezVega et al. (2009).
FIGURE 3. Principal component analysis based on chemical composition of the seed coat and cotyledon, considering the gene pool classification reported by PérezVega et al. (2009). Accessions with Middle American origin are represented by yellow symbols and accessions with Andean classification are represented by gray symbols. The dashed lines represent the 95% confidence ellipses.
Use of NIR in this study enabled significant, in depth evaluation of the chemical composition of the Spanish core collection of beans. Data obtained here have added further knowledge to previously reported findings, helping make this one of the most well studied and characterized bean and legume collections. Detailed information, including pictures of seeds and flowers for each accession can be found on the website of the Spanish National Plant Genetic Resources Center: http://www.crf.inia.es/crfesp/paginaprincipaljudia.asp (verified 20 June 2018).
The chemical analysis performed on 20 samples enabled the inclusion of new reference values to the NIR models (Plans et al., 2012, 2013), thus increasing their robustness, as described by Barton et al. (2000), who showed the effect of adding new samples (∼10) to the memory based models reduced ¼ the SEP of the predictions (Barton et al., 2000). Accordingly, the prediction of all the components with the improved NIR models showed good precision, as for all the traits under study GH was below 3, which has been described as the threshold for considering a prediction robust (Shenk et al., 1997).
All the accessions showed a good agronomic behavior, completing the crop cycle, and enabling the harvesting of seeds. This was probably because all the accessions are landraces originating in similar agroclimatic conditions to those of the experimental field, thus yielding a similar phenotype to that of their area of origin. We identified great variability among the accessions of the core collection in some important chemical traits analyzed in the seed coats and cotyledons. The coefficients of variation were low in some cases [e.g., starch, amylose, amylopectin, apparent amylose, and uronic acid all had CV below 10% (Table 1)]. However, even in these cases, the extreme values were found to be far apart (with the maximum value at least 30% higher than the minimum for all of these traits). The proportion of seed coat ranged from 4% to 10%, CV (11.45%). These results are in accordance with those reported by Singh et al. (1968) and PérezVega et al. (2009), and slightly different from the variability described by Moraghan et al. (2006).
The few studies that have analyzed the chemical composition of the seed coat and cotyledon separately have reported values within the ranges found in this study (Singh et al., 1968; Moraghan et al., 2006). We were unable to find other studies that differentiated between components in the seed coat and cotyledon, but as the seed coat represents only a small proportion of the total weight of the beans, we can assume that the amounts found for the analytes in the cotyledon are not very different from what would be found in an analysis of the entire seed. The genetic differences identified for proteins and carbohydrates are within the range reported in studies considering the whole species (Sathe et al., 1984a,b; Hayat et al., 2014), and somewhat larger than those reported in specific collections (Paredes et al., 2009). Thus, it seems that the Spanish core collection represents a rich source of genetic variation for chemical composition in beans. Moreover, this diversity is present in genotypes that show good adaptation to the agroclimatic conditions of the Iberian Peninsula. Considering that this collection has been characterized by several traits, specific accessions can be selected for their agronomic, chemical, and sensory profiles, and used as elite genotypes to perform breeding programs devoted to obtain new varieties with good agronomic behavior, high contents of nutritional compounds, and sensory profiles close to consumer demands. This strategy can promote the consumption of this legume and achieve the objective of increasing the proportion of vegetable protein in the human diet.
With the aim of determining whether correlations exist between the assessed parameters which could potentially limit the progress in breeding programs for quality in common beans, we conducted a correlation analysis considering all of the traits studied in the cotyledon and seed coat. In this study, we identified significant positive correlations among the analyzed compounds in the seed coat, such as Ca and ash (0.92, p < 0.001). Moreover, results showed a significant correlation between protein and carbohydrates. These results are in accordance with Vargas-Torres et al. (2004) and Casañas et al. (2013), and imply that an increase in the concentration of protein leads to a decrease in starch content. No correlations have been found between the chemical composition of the seed coat and cotyledon, signaling that the chemical composition of both parts seems to be under independent genetic control. In accordance with this hypothesis, Casañas et al. (2013) identified independent quantitative trait loci (QTL) controlling chemical composition of the seed coat and cotyledon.
With regards to the relationship between chemical composition and sensory characteristics, it seems that ashes and Ca are the components with the strongest influence. Relative to the culinary traits, some authors have reported a tendency for genotypes with shorter cooking times to retain higher nutrient concentrations than those with longer cooking times (Wiesinger et al., 2016). In this work, the only relationship with cooking time was observed for protein. The lack of data for phytate content in the collection is an important shortcoming to complete this analysis, as this compound shows a strong relation with the culinary or sensory traits in beans (Casañas et al., 2002; Coelho et al., 2007). Our results also showed a negative correlation between Ca and weight of 100 seeds (-0.37, p < 0.001). In a study of eight accessions, Moraghan and Grafton (2001) reported a much higher negative correlation between these two characteristics. However, the diversity in Ca concentration observed in some accessions with a similar 100 seed weight highlights that there are other factors involved in determining this characteristic, such as the color of the seeds (Table 2) or the seed coat area/cotyledon weight ratios (Moraghan and Grafton, 2001).
Although there were many significant genotypic correlations between chemical components and culinary or sensory traits, most were not strong enough to make indirect selection efficient (i.e., the selection for a desired trait using another trait that is genetically linked). These findings contrast with those of studies examining smaller collections, especially those done within a single variety or a few varieties, which found genetic correlations between chemical components and sensory traits (Casañas et al., 2006; Rivera et al., 2015). Therefore, in a large collection with greater variability, these effects are understandably diluted. Moreover, considering the important GxE effects on quantitative traits, and more specifically on traits related to chemical composition (Shellie and Hosfield, 1991; Moraghan and Grafton, 2001; Florez et al., 2009) and sensory profile (Romero del Castillo et al., 2008), significant correlations should be validated in further studies considering different agroclimatic conditions. For instance, in our case, we have performed the experiment in a field with high content of Ca and low content of P. Soil composition influences the chemical composition of common beans (Sameni et al., 1980; Moraghan and Grafton, 2001). Thus new experiments should be conducted in sites with highly different edaphologic characteristics to complement this study.
Seed coat color in common beans is controlled by several loci that act independently or in an epistatic manner to affect the color and pattern (Moghaddam et al., 2014). The P locus is known as the core factor for all seed coat color genotypes. The presence of its recessive allele results in white seeds and flowers due to its epistatic effect on the expression of the other color and pattern genes (Bassett, 2007). The background color and different patterns of marking on the seed coat are caused by the accumulation of anthocyanin and phenolic substances, which influence nutritional value and are coded by various gene systems (Beninger et al., 1999, 2000; Caldas and Blair, 2009). Our results show that there are significant differences in seed coat composition (ash, Ca, dietary fiber, Mg, and uronic acid) between white and colored seeds, but no differences for cotyledon composition (except for starch) (Table 3 and Figure 2). Although the negative correlation between starch and protein (r = -0.60, p < 0.001), which was confirmed within both groups (white seeds, r = -0.53, p < 0.01; colored seeds, r = -0.60, p < 0.001), the significant differences between color groups for the starch are not reflected in the protein content. Mean values for protein content in colored accessions was lower with respect to white seeded, but these differences were not statically supported, mainly because the variation within groups was very high for this trait. In contrast, other authors found greater amounts of protein and other nutrients in the whole seeds of black beans compared to lighter varieties (Silva et al., 2012; Hacisalihoglu and Settles, 2013). Our results concur with the results reported by Casañas et al. (2013), who reported five QTL associated with ash, Ca, dietary fiber, Mg, and uronic acid content, which were mapped in the region of the P gene.
Common beans were independently domesticated in the Andean and Middle American areas. Many authors have studied this phenomenon and attributed differentiated characteristics to beans originating from each area (Gepts and Debouck, 1991; Singh et al., 1991; Schmutz et al., 2014), but the differences for chemical composition have been scarcely studied. Research performed by PérezVega et al. (2009), comparing molecular markers and phaseolins in the Spanish core collection of common bean, enabled the classification of accessions into the Middle American or Andean gene pool. Following this classification, we performed an ANOVA that showed significant differences for all chemical components in the cotyledon, with the exception of amylopectin, and significant differences only for ash and Ca in the seed coat (Table 4 and Figure 3). Considering that phaseolins are the major seed storage protein (constitute over 50% of total protein in beans), and that the concentration and type of phaseolin present contribute to the classification gene pool (Gepts et al., 1986; Singh et al., 1991), it is to be expected that cotyledon components are those with the greatest differences depending on gene pool. Results obtained from our study indicate that seeds from the Spanish core collection identified as of the Middle American gene pool have a higher concentration of protein, ash, and Ca, and lower concentrations of carbohydrates. Some of these results are in accordance with those obtained by Islam et al. (2002). In their study on the CIAT core collection of beans, they looked at these and other chemical components in whole seeds from the two major gene pools (Middle American and Andean), as well as from the North Andean Group and a mixed group. Their study showed that accessions from the Middle American gene pool contained higher Ca concentrations than those from the North Andean and Andean gene pools. However, their study also reported a lower concentration of phaseolin protein in the Middle American gene pool.
In conclusion, results of this work show that the variability in the Iberian Peninsula, a second center of diversity for P. vulgaris, is very high, and can be used as an important source for breeding more nutritional and palatable varieties or directly by farmers. Moreover, our results point out that there are not strong correlations between the most important nutritional and sensory attributes, which is an important finding, signaling that synchronic improvement of both traits is feasible. Significant differences for nutritional composition have been identified between colored and non-colored seeds and between gene pools (Middle American and Mesoamerican), contributing to the knowledge about the diversification process of this species. NIRS models improved in this work can be an useful technology for mass phenotyping of other sources of genetic diversity within the species.
FC, AR, MP, and JSa planned the study. AR, AR, and JSa conducted the experiments. MP carried out the improvement of models and obtained the data. AR, FC, MP, and JC wrote the manuscript. JSi revised the article critically. All authors read and approved the final manuscript.
This work has been funded by the National Institute for Agronomic Research (INIA), projects RF2010-00005-C05 and RTA 2011-00076-C02-02.
Conflict of Interest Statement
Aastveit, A. H., and Marum, P. (1993). Near-infrared reflectance spectroscopy: different strategies for local calibrations in analyses of forage quality. Appl. Spectrosc. 47, 463–469. doi: 10.1366/0003702934334912
Araújo, S. S., Beebe, S., Crespi, M., Delbreil, B., González, E. M., Gruber, V., et al. (2015). Abiotic stress responses in legumes: strategies used to cope with environmental challenges. Crit. Rev. Plant Sci. 34, 237–280. doi: 10.1080/07352689.2014.898450
Asensio, C., Asensio-S-Manzanera, M., Ibeas, A., and De la Rosa, L. (2010). Resistance to halo blight common bacterial blight and bacterial brown spot in Spanish common bean core collection. Annu. Report. Bean Improv. Coop. 53, 101–111.
Barton, F. E., Shenk, J. S., Westerhaus, M. O., and Funk, D. B. (2000). The development of near infrared wheat quality models by locally weighted regressions. J. Near Infrared Spectrosc. 8, 201–208. doi: 10.1255/jnirs.280
Beninger, C. W., Hosfield, G. L., and Bassett, M. J. (1999). Flavonoid composition of three genotypes of dry bean (Phaseolus vulgaris L.) differing in seedcoat color. J. Am. Soc. Hort. Sci. 124, 514–518.
Beninger, C. W., Hosfield, G. L., Bassett, M. J., and Owens, S. (2000). Chemical and morphological expression of the B and Asp seedcoat genes in Phaseolus vulgaris L. J. Am. Soc. Hortic. Sci. 125, 52–58.
Caldas, G. V., and Blair, M. W. (2009). Inheritance of seed condensed tannins and their relationship with seed-coat color and pattern genes in common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 119, 131–142. doi: 10.1007/s00122-009-1023-4
Casañas, F., Pérez-Vega, E., Almirall, A., Plans, M., Sabaté, J., and Ferreira, J. J. (2013). Mapping of QTL associated with seed chemical content in a RIL population of common bean (Phaseolus vulgaris L.). Euphytica 192, 279–288. doi: 10.1007/s10681-013-0880-8
Casañas, F., Pujolà, M., Bosch, L., Sánchez, E., and Nuez, F. (2002). Chemical basis for low sensory perception of Ganxet bean (Phaseolus vulgaris L.) seed coat. J. Sci. Food Agric. 82, 1282–1286. doi: 10.1002/jsfa.1193
Casañas, F., Pujolà, M., Romero, del Castillo, R., Almiralll, A., Sanchez, E., et al. (2006). Variability in some texture characteristics and chemical composition of common beans (Phaseolus vulgaris L.). J. Sci. Food Agric. 86, 2445–2449. doi: 10.1002/jsfa.2638
Coelho, C. M. M., de Mattos Bellato, C., Santos, J. C. P., Ortega, E. M. M., and Tsai, S. M. (2007). Effect of phytate and storage conditions on the development of the ‘hard-to-cook’ phenomenon in common beans. J. Sci. Food Agric. 87, 1237–1243. doi: 10.1002/jsfa.2822
De la Rosa, L., Lázaro, A., and Varela, F. (2000). Racionalización de la colección española de Phaseolus vulgaris L. II Seminario de judía de la Península Ibérica. Villaviciosa: Asoc Esp Leguminosas, 55–62.
Ekmekcioglu, C., Wallner, P., Kundi, M., Weisz, U., Haas, W., and Hutter, H. P. (2018). Red meat, diseases, and healthy alternatives: a critical review. Crit. Rev. Food Sci. Nutr. 58, 247–261. doi: 10.1080/10408398.2016.1158148
FAO (2010). Second Report on the State of the World’s Plant Genetic Resources for Food and Agriculture. Available at: http://www.fao.org/docrep/013/i1500e/i1500e00.htm
FAO (2016). International Year of Pulses. Available at: http://www.fao.org/pulses-2016/faq/en/
Florez, A., Pujolà, M., Valero, J., Centelles, E., Almirall, A., and Casañas, F. (2009). Genetic and environmental effects on chemical composition related to sensory traits in common beans (Phaseolus vulgaris L.). Food Chem. 113, 950–956. doi: 10.1016/j.foodchem.2008.08.036
Gepts, P., and Debouck, D. (1991). “Origin, domestication, and evolution of the common bean (Phaseolus vulgaris L.),” in Common Beans: Research for Crop Improvement, eds A. van Schoonhoven and O. Voysest (Wallingford: CAB International-CIAT), 7–53.
Gepts, P., Osborn, T., Rashka, K., and Bliss, F. (1986). Phaseolin-protein variability in wild forms and landraces of the common bean (Phaseolus vulgaris): evidence for multiple centers of domestication. Econ. Bot. 40, 451–468. doi: 10.1007/BF02859659
Gonzales, J. F., Barnard, N. D., Jenkins, D. J. A., Lanou, A. J., Davis, B., Saxe, G., et al. (2014). Applying the precautionary principle to nutrition and cancer. J. Am. Coll. Nutr. 33, 239–246. doi: 10.1080/07315724.2013.866527
Gonzalez, A. D., Frostell, B., and Karlsson-Kanyama, A. (2011). Protein efficiency per unit energy and per unit greenhouse gas emissions: potential contribution of diet choices to climate change mitigation. Food Policy 36, 562–570. doi: 10.1016/j.foodpol.2011.07.003
Hacisalihoglu, G., Larbit, B., and Settles, A. M. (2010). Near-infrared reflectance spectroscopy predicts protein, starch and seed weight in intact seeds of common bean (Phaseolus vulgaris L.). J. Agric. Food Chem. 58, 502–506. doi: 10.1021/jf9019294
Hacisalihoglu, G., and Settles, A. M. (2013). Natural variation in seed composition of 91 common bean genotypes and their possible association with seed coat color. J. Plant Nutr. 36, 772–780. doi: 10.1080/01904167.2012.754041
Hayat, I., Ahmad, A., Masud, T., Ahmed, A., and Bashir, S. (2014). Nutritional and health perspectives of beans (Phaseolus vulgaris L.): an overview. Crit. Rev. Food Sci. Nutr. 54, 580–592. doi: 10.1080/10408398.2011.596639
Ishitani, M., Rao, I., Wenz, P., Beebe, S., and Tohme, J. (2004). Integration of genomics approach with traditional breeding towards improving abiotic stress adaptation: drought and aluminium toxicity as case studies. Field Crop Res. 90, 35–45. doi: 10.1016/j.fcr.2004.07.004
Islam, F. M. S., Basford, K. E., Jara, C., Redden, R. J., and Beebe, S. (2002). Seed compositional and disease resistance differences among gene pools in cultivated common bean. Genet. Resour. Crop. Evol. 49, 285–293. doi: 10.1023/A:1015510428026
Mendiburu, F. (2017). Agricolae: Statistical Procedures for Agricultural Research. R package version 1.2-8. Avilable at: https://CRAN.R-project.org/package=agricolae
Micha, R., Wallace, S. K., and Mozaffarian, D. (2010). Red and processed meat consumption and risk of incident coronary heart disease stroke, and diabetes mellitus: an updated review and meta-analysis. Circulation 121, 2271–2283. doi: 10.1161/CIRCULATIONAHA.109.924977
Michael, H., Tinashe, C., Ruaraidh, S. H., Brad, K., Emily, M., Eric, W., et al. (2018). Plant genetic resources for food and agriculture: opportunities and challenges emerging from the science and information technology revolution. New Phytol. 217, 1407–1419. doi: 10.1111/nph.14993
Miklas, P. N., Kelly, J. D., Beebe, S. E., and Blair, M. W. (2005). Common bean breeding for resistance against biotic and abiotic stresses: from classical to MAS breeding. Euphytica 147, 105–131. doi: 10.1007/s10681-006-4600-5
Mkanda, A. V., Minnaar, A., and de Kock, H. L. (2007). Relating consumer preferences to sensory and physicochemical properties of dry beans (Phaseolus vulgaris). J. Sci. Food Agric. 87, 2868–2879. doi: 10.1002/jsfa.3046
Moraghan, J. T., Etchevers, J. D., and Padilla, J. (2006). Contrasting accumulations of calcium and magnesium in seed coats and embryos of common bean and soybean. Food Chem. 95, 554–561. doi: 10.1016/j.foodchem.2004.10.060
Moraghan, J. T., and Grafton, K. (2001). Genetic diversity and mineral composition of common bean seed. 2001. J. Sci. Food Agric. 81, 404–408. doi: 10.1002/1097-0010(200103)81:4<404::AID-JSFA822>3.0.CO;2-H
Murdoch, D., and Chow, E. D. (2007). Ellipse: Functions for Drawing Ellipses and Ellipse like Confidence Regions. R package version 0.4.1. Avilable at: https://cran.r-project.org/package=ellipse
Nicolaï, B. M., Beullens, K., Bobelyn, E., Peirs, A., Sacys, W., Theoron, K. I., et al. (2007). Nondestructive measurement of fruit and vegetable quality by means of NIR spectroscopy: a review. Postharvest Biol. Technol. 46, 99–118. doi: 10.1016/j.postharvbio.2007.06.024
Pan, A., Sun, Q., Bernstein, A. M., Manson, J. A. E., Willett, W. C., and Hu, F. B. (2013). Changes in red meat consumption and subsequent risk of type 2 diabetes mellitus: three cohorts of US men and women. JAMA Intern. Med. 173, 1328–1335. doi: 10.1001/jamainternmed.2013.6633
Paredes, M. C., Becerra, V., and Tay, J. U. (2009). Inorganic nutritional composition of common bean (Phaseolus vulgaris L.) genotypes race Chile. Chil. J. Agric. Res. 69, 486–495. doi: 10.4067/S0718-58392009000400002
Pascual, A., Campa, A., Pérez-Vega, E., Giraldez, R., Miklas, P. N., and Ferreira, J. J. (2010). Screening common bean for resistance to four Sclerotinia sclerotiorum isolates collected in northern Spain. Plant Dis. 94, 885–890. doi: 10.1094/PDIS-94-7-0885
Pérez-Vega, E., Campa, A., de la Rosa, L., Giraldez, R., and Ferreira, J. J. (2006). Evaluation for resistance to anthracnose of a core collection established from the CRF-INIA common bean collection. Annu. Rep. Bean Improv. Coop. 49, 201–202.
Pérez-Vega, E., Campa, A., de la Rosa, L., Giraldez, R., and Ferreira, J. J. (2009). Genetic diversity in a core collection established from the main Spanish bean genbank. Crop Sci. 49, 1372–1386. doi: 10.2135/cropsci2008.07.0409
Pinheiro, C., Baeta, J. P., Pereira, A. M., Domingues, H., and Ricardo, C. P. (2010). Diversity of seed mineral composition of Phaseolus vulgaris L. germplasm. J. Food Comp. Anal. 23, 319–325. doi: 10.1016/j.jfca.2010.01.005
Plans, M., Simó, J., Casañas, F., Romero, del Castillo, R., Rodriguez-Saona, L. E., et al. (2014). Estimating sensory properties of common beans (Phaseolus vulgaris L.) by near infrared spectroscopy. Food Res. Int. 56, 55–62. doi: 10.1016/j.foodres.2013.12.003
Plans, M., Simó, J., Casañas, F., and Sabaté, J. (2012). Near-infrared spectroscopy analysis of seed coats of common beans (Phaseolus vulgaris L.): a potential tool for breeding and quality evaluation. J. Agric. Food Chem. 60, 706–712. doi: 10.1021/jf204110k
Plans, M., Simó, J., Casañas, F., Sabaté, J., and Rodriguez-Saona, J. L. (2013). Characterization of common beans (Phaseolus vulgaris L.) by infrared spectroscopy: comparison of MIR, FT-NIR and dispersive NIR using portable and benchtop instruments. Food Res. Int. 54, 1643–1651. doi: 10.1016/j.foodres.2013.09.003
Ramirez-Lopez, L., and Stevens, A. (2016). Resemble: Regression and Similarity Evaluation for Memory-Based Learning in Spectral Chemometrics. R package version 1.2.2. Available at: https://CRAN.R-project.org/package=resemble
Rivera, A., Casquero, P., Mayo, S., Almirall, A., Plans, M., Simó, J., et al. (2016). Culinary and sensory traits diversity in the Spanish Core Collection of common beans (Phaseolus vulgaris L.). Span. J. Agric. Res. 14:e0701. doi: 10.5424/sjar/2016141-7726
Rivera, A., Roselló, S., and Casañas, F. (2015). Seed curvature as a useful marker to transfer morphologic, agronomic, chemical and sensory traits from Ganxet common bean (Phaseolus vulgaris L.). Sci. Hortic. 197, 476–482. doi: 10.1016/j.scienta.2015.10.006
Romero del Castillo, R., Almirall, A., Valero, J., and Casañas, F. (2008). Protected designation of origin in beans (Phaseolus vulgaris L.): towards an objective approach based on sensory and agromorphological properties. J. Sci. Food Agric. 88, 1954–1962. doi: 10.1002/jsfa.3303
Romero del Castillo, R., Costell, E., Plans, M., Simó, J., and Casañas, F. (2012). A standardized method of preparing common beans (Phaseolus vulgaris L.) for sensory analysis. J. Sens. Stud. 27, 188–195. doi: 10.1111/j.1745-459X.2012.00381
Sameni, A. M., Maftoun, M., Bassiri, A., and Sepaskhah, A. R. (1980). Growth and chemical composition of dry beans as affected by soil salinity and N fertilization. Plant Soil 54, 217–222. doi: 10.1007/BF02181847
Santalla, M., de Ron, A. M., and Voysest, O. (2001). “European bean market classes,” in Catalogue of Bean Genetic Resources, eds M. Amurrio, M. Santalla, and A. M. de Ron (Corunha: Fundación Pedro Barrie de la Maza), 79–94.
Santalla, M., Rodiño, A., and de Ron, A. M. (2002). Allozyme evidence supporting southwestern Europe as a secondary center of genetic diversity for common bean. Theor. Appl. Genet. 104, 934–994. doi: 10.1007/s00122-001-0844-6
Sathe, S. K., Desphande, S. S., Salunkhe, D. K., and Rackis, J. J. (1984a). Dry beans of phaseolus. A review. Part 1. Chemical composition: proteins. Crit. Rev. Food Sci. 20, 1–46. doi: 10.1080/10408398409527382
Sathe, S. K., Desphande, S. S., Salunkhe, D. K., and Rackis, J. J. (1984b). Dry beans of phaseolus. A review. Part 2. Chemical composition: carbohydrates, fiber, minerals, vitamins and lipids. Crit. Rev. Food Sci. 21, 41–93. doi: 10.1080/10408398409527396
Schmutz, J., McClean, P. E., Mamidi, S., Wu, G. A., Cannon, S. B., Grimwood, J., et al. (2014). A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 46, 707–713. doi: 10.1038/ng.3008
Shenk, J. S., and Westerhaus, M. O. (1991). Population definition, sample selection and calibration procedures for near infrared reflectance spectroscopy. Crop Sci. 31, 469–474. doi: 10.2135/cropsci1991.0011183X003100020049x
Silva, C. A., Abreu, A. F. B., Ramalho, M. A. P., and Maia, L. G. S. (2012). Chemical composition as related to seed color of common bean. Crop Breed. Appl. Biotechnol. 12, 132–137. doi: 10.1590/S1984-70332012000200006
Singh, S. P. (1999). “Integrated genetic improvement,” in Common Bean Improvement in the Twenty-First Century Developments in Plant Breeding, Vol. 7, ed. S. P. Singh (Dordrecht: Springer), 133–165. doi: 10.1007/978-94-015-9211-6
Smýkal, P., Coyne, C. J., Ambrose, M. J., Maxted, N., Schaefer, H., Blair, M. W., et al. (2015). Legume crops phylogeny and genetic diversity for science and breeding. Crit. Rev. Plant Sci. 34, 43–104. doi: 10.1080/07352689.2014.897904
Stacklies, W., Redestig, H., Scholz, M., Walther, D., and Selbig, J. (2007). PcaMethods–a bioconductor package providing PCA methods for incomplete data. Bioinformatics 23, 1164–1167. doi: 10.1093/bioinformatics/btm069
Vargas-Torres, A., Osorio-Díaz, P., Tovar, J., Paredes-López, O., Ruales, J., and Bello-Pérez, L. A. (2004). Chemical composition, starch bioavailability and indigestible fraction of common beans (Phaseolus vulgaris L.). Starke 56, 74–78. doi: 10.1002/star.200300205
Vaz Patto, M. C., Amarowicz, R., Aryee, A. N. A., Boye, J. I., Chung, H. J., Martín-Cabrejas, M. A., et al. (2015). Achievements and challenges in improving the nutritional quality of Food legumes. Crit. Rev. Plant Sci. 34, 105–143. doi: 10.1080/07352689.2014.897907
WHO (2015). World Health Organization. Available at: http://www.who.int/mediacentre/factsheets/fs394/en/
Wiesinger, J. A., Cichy, K. A., Glahn, R. P., Grusak, M. A., Brick, M. A., Thompson, H. J., et al. (2016). Demonstrating a nutritional advantage to the fast-cooking dry bean (Phaseolus vulgaris L.). J. Agric. Food Chem. 64, 8592–8603. doi: 10.1021/acs.jafc.6b03100
Keywords: Phaseolus vulgaris, common bean, diversity, genebanks, gene pool, nutrient composition, protein concentration, seed color
Citation: Rivera A, Plans M, Sabaté J, Casañas F, Casals J, Rull A and Simó J (2018) The Spanish Core Collection of Common Beans (Phaseolus vulgaris L.): An Important Source of Variability for Breeding Chemical Composition. Front. Plant Sci. 9:1642. doi: 10.3389/fpls.2018.01642
Received: 07 June 2018; Accepted: 23 October 2018;
Published: 13 November 2018.
Edited by:Spyridon Alexandros Petropoulos, University of Thessaly, Greece
Reviewed by:Paul John Hunter, University of Warwick, United Kingdom
Caspar Christian Cedric Chater, University of Sheffield, United Kingdom
Copyright © 2018 Rivera, Plans, Sabaté, Casañas, Casals, Rull and Simó. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Rivera, firstname.lastname@example.org