Bioactive Compounds, Antioxidant Activity, and Antinutritional Content of Legumes: A Comparison between Four Phaseolus Species

Beans (Phaseolus spp.) are one of the most important legumes for their nutritional value and health benefits in many world regions. In addition to Phaseolus vulgaris, there are four additional species that are cultivated in many regions of the world and are a source of food for human consumption: P. lunatus, P. coccineus, P. polyanthus, and P. acutifolius. In this work, phenolic compounds, antioxidant activity, and anti-nutritional compounds of 18 bean accessions, corresponding to four different species of the genus Phaseolus, were analyzed. In addition, their physical characteristics, proximate composition, and amino acid content were determined in order to compare their phytochemical composition and nutritional value. The species closest to each other in terms of essential amino acid content were P. polyanthus with P. vulgaris and P. lunatus with P. coccineus. Furthermore, there was a strong positive correlation between antioxidant activity and flavonoids, anthocyanins, and lectins with all the accessions collected. Significant differences in the content of phenolic compounds were found among the bean species studied. Therefore, in addition to P. vulgaris, other species such as P. coccineus and P. lunatus have high biological and antioxidant potential that could be beneficial to human health when consumed as nutraceutical foods.


Introduction
Beans (Phaseolus spp.) belong to the legume family (Fabaceae), along with peas (Pisum sativum L.), soybeans (Glycine max (L.) Merril), and lentils (Lens culinaris Medik). Beans are the most important legume and are consumed by approximately 300 million people. Moreover, they are considered the second source of protein in East and Southeast Africa and the fourth source in America [1].
There are about 150 known species of beans around the world. In Mexico alone, there are approximately 65 species, 52 of which belong to the Phaseolus genus, and 31 species are endemic [2]. The largest number of species are distributed in Western Mexico, with P. vulgaris being the most studied species. Thus, different studies have quantified the content of phytochemicals in some varieties of P. vulgaris, such as: Total phenols, anthocyanins, tannins, flavonoids, lectins, phytic acid, and oligosaccharides [3,4]. In addition to the protein, carbohydrate, fiber, mineral, and vitamin contents [5].
23.75 ± 0.71 b 28.47 ± 0.77 h 0.73 ± 0.01 b-d 12.9 ± 1.1 c-f 1.48 ± 0.01% to −1.48 ± 0.08 (Table 1). Luminosity (*L) is the radiant intensity per unit area and is associated with the color black to white. Chromaticity (c) is an index similar to saturation and indicates the degree of deviation from grey to pure chromatic color. The hue (h) angle is measured from pure red to orange, yellow, green, etc. [11]. P. vulgaris accessions exhibited the most diverse range of color from black, red, brown and white tones, as 13 accessions out of 18 correspond to this species; the accessions with purple tones were both P. coccineus samples (CH-07 and CH-16); P. polyanthus accessions presented yellow (CH-06) and red (CH-11) tones; and P. lunatus accessions presented a pink tone .
Furthermore, of the accessions collected, 22% were small, 33% were medium, and 45% were large. The P. polyanthus and P. coccineus species were the widest, and P. lunatus showed the greatest thickness, while P. coccineus was the longest. However, P. vulgaris accessions were distinguished as the most heterogeneous species in terms of size, since a higher number of collected accessions of this species were found in the Mayan region of Mexico (Table 1).

Proximate Composition of Bean Seeds
Carbohydrates were the predominant element in bean seeds of all accessions ( Figure 1). Furthermore, according to analysis of variance, there was a significant difference between the bean species in terms of their percent moisture, ash, fats, and protein content (p-value < 0.05). Moreover, P. coccineus accessions had a higher percent by weight of raw fiber and fat, but a lower amount of protein than P. vulgaris; whereas, P. lunatus had a high percent of carbohydrates and a low amount of ash, fats and raw fiber. However, the accessions from P. polyanthus presented a greater composition of moisture and ash. Finally, the P. vulgaris accessions presented a higher percent of protein and a lower amount of moisture and carbohydrates ( Figure 1).

Essential Amino Acids
One of the main traits among the bean species studied is the content of essential amino acids. It was observed that there is a closeness between certain species, which in turn show certain distance from the rest of the other species, clustering in two main clades. The closest species to each other in terms of essential amino acids were P. coccineus with P. lunatus and P. polyanthus with P. vulgaris ( Figure 2). The latter two species showed a higher concentration of essential amino acids. Arginine was Molecules 2020, 25,3528 5 of 21 the major amino acid, ranging from 661.25 to 115.67 mg/100 g bean flour, and in descending order were tryptophan (108.06-10.34 mg/100 g), histidine (61.7-6.35 mg/100 g), threonine (38.34-8.59 mg/100 g), valine (37.17-8.49 mg/100 g), lysine (28.09-7.99 mg/100 g), leucine (26.92-11.45 mg/100 g), phenylalanine (18.44-4.11 mg/100 g), isoleucine (12.81-2.64 mg/100 g), and methionine (6.98-0.81 mg/100 g).
One of the main traits among the bean species studied is the content of essential amino acids. It was observed that there is a closeness between certain species, which in turn show certain distance from the rest of the other species, clustering in two main clades. The closest species to each other in terms of essential amino acids were P. coccineus with P. lunatus and P. polyanthus with P. vulgaris ( Figure 2). The latter two species showed a higher concentration of essential amino acids. Arginine was the major amino acid, ranging from 661.25 to 115.67 mg/100 g bean flour, and in descending order were tryptophan (108.06-10.34 mg/100 g), histidine (61.7-6.35 mg/100 g), threonine (38.34-8.59 mg/100 g), valine (37.17-8.49 mg/100 g), lysine (28.09-7.99 mg/100 g), leucine (26.92-11.45 mg/100 g), phenylalanine (18.44-4.11 mg/100 g), isoleucine (12.81-2.64 mg/100 g), and methionine (6.98-0.81 mg/100 g). Histidine and isoleucine presented a significant difference between P. coccineus, P. polyanthus, and P. vulgaris (p-value < 0.05). P. polyanthus was distinguished from the other species by its high content of threonine, valine, and isoleucine. While, P. vulgaris had the highest content of histidine, arginine, and tryptophan. In general, P. coccineus contained the lowest proportion of amino acids ( Table 2).  Histidine and isoleucine presented a significant difference between P. coccineus, P. polyanthus, and P. vulgaris (p-value < 0.05). P. polyanthus was distinguished from the other species by its high content of threonine, valine, and isoleucine. While, P. vulgaris had the highest content of histidine, arginine, and tryptophan. In general, P. coccineus contained the lowest proportion of amino acids (Table 2).

Total Phenols
The phenolic compound concentrations ranged from 1364.54 ± 0.5 to 114.9 ± 0.7 mg GAE/g dw (Table 3), where sample CH-09 (P. vulgaris, red seed hull) presented the highest concentration. The species concentration of total phenols in descending order was P. vulgaris > P. lunatus > P. polyanthus > P. coccineus. However, significant differences in total phenol concentration were found between bean species P. vulgaris and P. coccineus (Figure 3). The phenolic compound concentrations ranged from 1364.54 ± 0.5 to 114.9 ± 0.7 mg GAE/g dw (Table 3), where sample CH-09 (P. vulgaris, red seed hull) presented the highest concentration. The species concentration of total phenols in descending order was P. vulgaris > P. lunatus > P. polyanthus > P. coccineus. However, significant differences in total phenol concentration were found between bean species P. vulgaris and P. coccineus ( Figure 3).

Flavonoids
The flavonoid concentration of the accessions ranged from 19.94 ± 0.47 to 1.09 ± 0.33 mg QE/g dw (Table 3), and CH-09 (P. vulgaris, red seed hull) showed the highest concentration. Hence, the'concentration of flavonoids in descending order was P. vulgaris > P. lunatus > P. coccineus > P. polyanthus. Furthermore, there were significant differences between the bean species and their flavonoid concentrations (p-value = 0.007). Thus, the accessions collected from P. coccineus, P. lunatus, and P. polyanthus species showed significant differences with P. vulgaris accessions ( Figure 3).

Condensed Tannins
The concentration of tannins ranged from 0.53 ± 0.05 to 5.57 ± 0.02 mg CAE/g dw (Table 3). CH-10 (P. vulgaris, black) showed the highest concentration of condensed tannins. Moreover, red and black accessions of P. vulgaris showed a high concentration of tannins, followed by both accessions of P. coccineus (purple) and P. lunatus (pink). Both P. polyanthus (yellow and red) and P. vulgaris (brown and white) accessions had the lowest tannin content. As a result, the species P. polyanthus and P. lunatus displayed no significant differences, but there were a significant difference (p < 0.05) with P. coccineus and P. vulgaris (Figure 3).

Anthocyanins
The concentration of anthocyanins in the bean seed coat oscillated between not being detecting and 9.42 ± 0.03 mg of cyanidin 3-glucoside per gram of the bean seed coat (Table 3). Similar to tannins, CH-10 (P. vulgaris, black) showed the highest concentration of anthocyanins. Moreover, the red and black accessions of P. vulgaris species showed a high concentration of anthocyanins, followed by both species of P. coccineus (purple) and P. vulgaris (brown); both species of P. polyanthus (yellow and red), P. lunatus (pink), and P. vulgaris (white) did not show any anthocyanin content (Table 3).
In addition, when a multiple range analysis of the means in terms of tannins and seed coat color was performed, it was observed that there were significant differences between the black genotype and the white, brown, and red genotypes from all the species collected, with the black seeds presenting the highest tannin contents. Moreover, there is a positive correlation between the anthocyanins and the tannin content (Table 4). In the present study, the content of cyanidin 3-glucoside, genistein, ferulic acid, quercetin 3-glucoside, and kaempferol 3-glucoside was determined. Statistical differences were observed between the analyzed compounds (p-value = 0.031) (Figure 4). Cyanidin and kaempferol were only detected in P. phaseolus; whereas, quercetin was quantified in three of the four species analyzed (P. coccineus, P.'lunatus, and P. vulgaris). Similarly, genistein was found in P. coccineus, P. polyanthus and P. vulgaris. Ferulic acid was observed in P. lunatus and P. vulgaris. In the species P. polyanthus only genistein was identified while in P. vulgaris all the phenolic compounds analyzed were detected ( Figure 4).

Antioxidant Activity
Accession CH-18 (96.9% ARSA; 98.2% DRSA) showed the highest antioxidant activity in both the ABTS and DPPH assays. However, the accessions with the lowest antioxidant activity were CH-04 from the ABTS assay (56.7% ARSA) and CH-16 from the DPPH assay (78.9% DRSA). However, all bean accessions showed antioxidant activity with both methodologies (Table 3), but no significant differences were observed between the phenolic compounds and the antioxidant activity with the ABTS technique.
Interspecies analysis indicated that P. coccineus had the lowest phytic acid and lectin content, an intermediate value in trypsin inhibitory activity and oligosaccharide. P. lunatus showed an intermediate value of phytic acid, together with P. polyanthus, the highest trypsin inhibitory activity,

Antioxidant Activity
Accession CH-18 (96.9% ARSA; 98.2% DRSA) showed the highest antioxidant activity in both the ABTS and DPPH assays. However, the accessions with the lowest antioxidant activity were CH-04 from the ABTS assay (56.7% ARSA) and CH-16 from the DPPH assay (78.9% DRSA). However, all bean accessions showed antioxidant activity with both methodologies (Table 3), but no significant differences were observed between the phenolic compounds and the antioxidant activity with the ABTS technique.
Interspecies analysis indicated that P. coccineus had the lowest phytic acid and lectin content, an intermediate value in trypsin inhibitory activity and oligosaccharide. P. lunatus showed an intermediate value of phytic acid, together with P. polyanthus, the highest trypsin inhibitory activity, and the highest oligosaccharides content of the four species analyzed. P. polyanthus was distinguished by the highest phytic acid content. While, P. vulgaris exhibited the lowest trypsin inhibitory activity, the highest lectin activity, and the lowest oligosaccharide content than the rest of the Phaseolus species ( Figure 5).
Molecules 2020, 25, x FOR PEER REVIEW 11 of 22 and the highest oligosaccharides content of the four species analyzed. P. polyanthus was distinguished by the highest phytic acid content. While, P. vulgaris exhibited the lowest trypsin inhibitory activity, the highest lectin activity, and the lowest oligosaccharide content than the rest of the Phaseolus species ( Figure 5).

Discriminant Analysis
Additionally, a discriminant analysis was performed with the phenolic compound concentration and antioxidant activity and showed that the accessions were correctly grouped 83.3% of the time. Only two accessions of P. vulgaris exhibited atypical behavior. These samples were CH-12, which was grouped into P. polyanthus, and CH-01, which was grouped close to P. lunatus ( Figure 6a).
Discriminant analysis from the antinutritional compounds showed that the samples were correctly grouped 88.9% of the time (Figure 6b), and only two atypical samples were observed: CH-01 (P. vulgaris, red), which was grouped with P. lunatus, and sample CH-14 (P. vulgaris, red), which was grouped with the species P. coccineus.
Both discriminant analyses between phenolic compounds with antioxidant activity and antinutritional components were performed separately, and atypical behavior was observed in some accessions. However, combining all the variables described above, the accessions were correctly grouped 100% of the time. The closest varieties to each other were P. coccineus and P. lunatus ( Figure  6c).

Discriminant Analysis
Additionally, a discriminant analysis was performed with the phenolic compound concentration and antioxidant activity and showed that the accessions were correctly grouped 83.3% of the time. Only two accessions of P. vulgaris exhibited atypical behavior. These samples were CH-12, which was grouped into P. polyanthus, and CH-01, which was grouped close to P. lunatus (Figure 6a).
Discriminant analysis from the antinutritional compounds showed that the samples were correctly grouped 88.9% of the time (Figure 6b), and only two atypical samples were observed: CH-01 (P. vulgaris, red), which was grouped with P. lunatus, and sample CH-14 (P. vulgaris, red), which was grouped with the species P. coccineus.
Both discriminant analyses between phenolic compounds with antioxidant activity and antinutritional components were performed separately, and atypical behavior was observed in some accessions. However, combining all the variables described above, the accessions were correctly grouped 100% of the time. The closest varieties to each other were P. coccineus and P. lunatus ( Figure 6c).

Discussion
In this work, four bean species (P. vulgaris, P. coccineus, P. polyanthus, and P. lunatus) were tested. The results are consistent with another research, where it was observed that dark-colored seeds present the lowest luminosity and chromaticity [3]. Therefore, the dimensions collected from the P. vulgaris accessions (Table 2) could be considered within the ranges established for the bean races that have been reported to be domesticated in Mesoamerica [12].
However, the P. lunatus accession collected is similar in weight to the varieties that belong to the Andean germplasm [13]; the above might be due to the proximity of the Chiapas state with the countries in Central America, where this species is one of the most cultivated. Moreover, for P. coccineus, both species collected, according to the morphological analysis, would be classified as medium according to their length [14], although in terms of their weight, they would be classified as large varieties (Table 2).
Along with the purple color of the seed coat, it could be assumed that the accessions collected were more similar to the Andean germplasm, such as P. lunatus, which might be since the Mayan region encompasses other Central American countries [12,13]. This indicated that the different species of beans consumed in Chiapas were cultivated and spread throughout the Mayan territory. In the case of P. polyanthus, the data coincide with previous research [15], and the accessions collected were within the range of large seeds.

Discussion
In this work, four bean species (P. vulgaris, P. coccineus, P. polyanthus, and P. lunatus) were tested. The results are consistent with another research, where it was observed that dark-colored seeds present the lowest luminosity and chromaticity [3]. Therefore, the dimensions collected from the P. vulgaris accessions ( Table 2) could be considered within the ranges established for the bean races that have been reported to be domesticated in Mesoamerica [12].
However, the P. lunatus accession collected is similar in weight to the varieties that belong to the Andean germplasm [13]; the above might be due to the proximity of the Chiapas state with the countries in Central America, where this species is one of the most cultivated. Moreover, for P. coccineus, both species collected, according to the morphological analysis, would be classified as medium according to their length [14], although in terms of their weight, they would be classified as large varieties (Table 2).
Along with the purple color of the seed coat, it could be assumed that the accessions collected were more similar to the Andean germplasm, such as P. lunatus, which might be since the Mayan region encompasses other Central American countries [12,13]. This indicated that the different species of beans consumed in Chiapas were cultivated and spread throughout the Mayan territory. In the case of P. polyanthus, the data coincide with previous research [15], and the accessions collected were within the range of large seeds.
The protein content in legumes ranges from 17 to 40% [16]. In the particular case of common beans (P. vulgaris), the protein content depends on the type of bean, and ranges from 14 to 33% are reported. Additionally, the protein is enriched with amino acids such as lysine, phenylalanine, and tyrosine [2]. In the present work, it was found that the CH-09 accession (P. vulgaris) showed the highest protein content (30.6%), while the lowest content (19.5%) was found in sample CH-16 (P. coccineus). The protein content of CH-09 is higher than that reported in other Mexican varieties of common beans, where the maximum protein content is 26% [2]. Meanwhile, the CH-16 accession has 2.4% less protein than another purple variety (21.9%) of P. coccineus [17]. In addition to the variation between species and between varieties of the same species, other factors such as fertilizer use, agricultural practices, and soil-climate characteristics of cultivation sites can influence the protein content of bean seeds [17,18].
Cluster analysis between the species (Figure 2), based on the number of free amino acids, revealed a closeness between P. vulgaris and P. polyanthus, as well as between P. coccineus and P. lunatus, with arginine being the predominant amino acid in all species studied. Additionally, arginine was the major amino acid in all the accessions collected. Therefore, the accessions of beans collected could help to prevent the risk of vascular and heart diseases, improve the immune response and renal function, and stimulate the release of pancreatic insulin [19]. Thus, the accessions collected showed potential to produce a nutraceutical agent.
Phenolic compounds, including simple molecules such as phenolic acids and flavonoids, which are compounds with intermediate molecular weights, and tannins, with high molecular weights, were extracted and quantified for their total phenolic compounds, flavonoids and condensed tannins with different assays. The solvent mix extraction used was adequate since the concentrations reached in each assay were equal to or even higher than those reported in the literature [20].
The highest concentration of phenolic compounds and antioxidant activity were found in P. vulgaris accessions with red and black seed coats ( Figure 3). These results were consistent with another study, who found a positive correlation between total polyphenols and anthocyanins and their antioxidant activity from different varieties of P. vulgaris relative to their color [8]. In addition, the accessions collected from P. coccineus, P. lunatus, and P. polyanthus could reach equal or even higher values of phenolic compounds and antioxidant activity than the accessions from P. vulgaris with light colors. Thus, these results showed potential antioxidant activity in all species collected [2].
Furthermore, significant differences were observed between bean species and anthocyanin concentrations. The species that showed significant differences were P. coccineus and P. polyanthus with P. lunatus, and with P. vulgaris ( Figure 3). Similar results were reported when P. coccineus accessions were compared with P. vulgaris accessions, where a higher concentration of polyphenols, flavonoids, anthocyanins, and even higher antioxidant activity was found in the accessions exhibiting purple, black, and brown colors [18,21].
In addition, some authors found that the red genotype of different P. vulgaris varieties showed significant differences in anthocyanin production among purple-and brown-colored genotypes [22]. Nevertheless, in the present study, the black genotype presented significant differences in the concentration of anthocyanins (p-value = 0.002) among the seed coat color of all species. These results were consistent with those reported by Capistrán-Carabarin et al. [21], who found that native beans with black seed coats had the highest content of anthocyanins, even if they were from different species. There is some evidence that has reported the relationship between the anthocyanin and the tannin content and the color of the beans [22]. Thus, this work confirmed the relationship between anthocyanin and tannins with seed color. The black genotype showed a high concentration of both tannins and anthocyanins.
Certain reports have shown that the main flavonoids in some varieties of P. vulgaris are quercetin and kaempferol, and both compounds have been shown to decrease the risk of lung cancer and cardiovascular disease [2,3]. In this study, it was observed that kaempferol was only detected in some accessions of P. vulgaris, while quercetin was found in almost all the collected accessions, except for the P. polyanthus varieties where the presence of this flavonoid was not detected (Figure 4). However, the presence of genistein was observed in both accessions of P. polyanthus, at even higher concentrations than those detected in the accessions of P. coccineus and P. vulgaris. This isoflavonoid is characterized as a phytoestrogen compound that inhibits carcinogenic cells, including breast and prostate cancer, and has until now only been reported in some varieties of P. vulgaris during the germination stage, as well as in soybeans [23]. This could suggest that bean accessions collected in the Mayan region of Mexico have high potential to be consumed and improve human health by reducing the risk of some types of cancer.
Regarding ferulic acid content, P. lunatus presented an average concentration higher than the accessions of P. vulgaris collected (Figure 4). Therefore, this accession may have high antioxidant potential. In the case of anthocyanins, the present work quantified cyanidin 3-glucoside, which was only observed in some varieties of P. vulgaris (CH-05, CH-10, and CH-13). This confirms that the composition of anthocyanins in bean seeds can vary, even if the varieties have the same color and that the anthocyanin content depends on both external (climatic and soil conditions) and internal (genetic variation) factors [8,24].
In this sense, some studies have indicated that different species of legumes have a wide variation in their capacity to absorb, translocate, and store flavonoids, polyphenols, and anthocyanins, and together with their variable interaction with cellulose and pectin, this greatly limits the bioavailability of phenolic compounds. These interactions cause great diversity of the compounds present in the bean seeds [21,25]. Finally, it was observed that the phenolic compounds detected have some hydroxyl groups. Hence, it could be assumed that the greater the number of hydroxyl groups, the greater the antioxidant activity [26].
According to correlation analysis (Table 4), the antioxidant activity with the DPPH assay showed a strong positive relationship with flavonoids (r = 0.95) and anthocyanins (r = 0.96). These correlations could be related to the extraction system used, since there is some evidence that acetone solvent extraction exhibits a higher yield of flavonoids and ethanol solvent extraction exhibits a higher yield of anthocyanins [20,27]. Furthermore, some reports point out that the interaction between the compounds with antioxidant potential and the DPPH* radical is affected by the medium where these compounds were extracted [20]. Thus, the antioxidant activity depends not only on the concentration of phenolic compounds but also on the polarity and chemical structure of each phenolic compound [2,28]. Therefore, it could be assumed that the mechanism of interaction between the antioxidant potential of collected beans and the DPPH* radical is more structurally related and responds faster, reducing the number of DPPH* molecules in the medium [29].
Among the antinutritional compounds, reported in beans, are phytic acid, trypsin inhibitory activity, lectins and oligosaccharides, where significant differences were found among the four species evaluated. However, recent studies have reported that trypsin inhibitors and lectins could confer some health benefits. For example, trypsin inhibitors could confer protection against rotavirus, inhibit some types of carcinogenesis, and could be used as chemopreventive agents [10], and lectins may decrease lymphoma growth and could be used as diagnostic markers for tumors, as well as help in the prevention of obesity [2].
The antinutritional compounds were determined in raw bean seeds, therefore it is important to mention that there are several methods and processing techniques to reduce the levels of these antinutritional factors. Treatments such as soaking, milling, boiling, microwave cooking, autoclaving, hydrothermal processing, and germination, either individually or in combination, can inhibit or reduce the tannin content, trypsin inhibitor activity, phytic acids and hemagglutinins [30].
In this work, the results showed a positive correlation (r = 0.97) between lectins and antioxidant activity (DPPH assay), where accessions CH-16 (P. coccineus) and CH-12 (P. vulgaris) presented the lowest and highest contents, respectively. Therefore, it could be assumed that there is a relationship between lectins and antioxidant activity that could help prevent some types of diseases related to inflammation [31]. Thus, it could be confirmed that the accessions collected have a high potential to be used as precursors for nutraceutical products due to their chemical composition and high nutritional content.
Taking into account all the results, the bean species collected in the Mayan region of Mexico exhibit high potential for biological activity that could be used to develop products due their nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deep-rooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 ( Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 ( Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 ( Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 ( Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33]. nutritional, phenolic content, and their antioxidant effect could prevent some chronic diseases related to the oxidative stress [32].

Plant Material
Different bean accessions were collected in various rural communities in the Mayan region of Chiapas State, Mexico, in March 2019 (Table 5). The accessions were classified according to their physical and morphological characteristics within the bean species P. vulgaris, P. coccineus, P. lunatus, or P. polyanthus. These species are located in the Mayan region, particularly in the State of Chiapas, which is characterized by great agrodiversity. In this region, rural areas with low economic resources, a lack of social assistance and information on health issues, low academic preparation, and the deeprooted predominance of their culture, among other factors, influence the nutritional status of these communities [33].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 °C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Physicochemical Analysis
Once the plant material was collected, the bean seeds were washed to remove any foreign material and dried in a convection oven for 4 h at 40 • C (Convection oven plus HCX II, SAN-SON, Mexico). Afterwards, the seeds were ground to fine powder and kept in hermetically sealed containers for further determination of phenolic compounds, antioxidant activity, and content of antinutritional factors.

Bean Size Determination
One hundred bean seeds of each accession were randomly measured (length, with, and thickness) and weighed to determine if the seeds were small (<25 g), medium (25-40 g) or large (>40 g) [13].

Bean Color Evaluation
These parameters were determined according to the CIE L*A*B* scale using a spectrophotometer (CM-5, Konica Minolta Sensing Americas, Ramsey, NJ, USA) and were calculated according to the following equations: Brightness (L) with a range of 0 (black) to 100 (white), chromaticity (c) and hue (h); these last two parameters were calculated according to the following equations: c = (a *2 + b *2 ) 1/2 (1) where a* = this variable being positive represents red tones and this value being negative represents green tones; b* = this variable being positive represents yellow tones and this variable being negative represents blue tones [11].

Proximate Composition Analysis of Bean Seeds
This analysis was performed according to the AOAC international official methods of analysis: protein (960.52), moisture (945.38), ash (930.03), lipids (920.30), raw fiber (985.28), and total carbohydrates by difference from the proximate composition [34].

Extraction of Phenolic Compounds
Bean flour of the 18 accessions was obtained by milling (A11 mill, IKA, Staufen, Germany) and passing through MESH 60. Later, maceration was carried out at room temperature from 16 h by weighing 15 g of bean flour (dry weight, dw) and 150 mL of acetone solvent mixture (acetone, water, and acetic acid, 70:29.5:05 v/v/v). Then, the extracts were centrifuged at 3488× g for 15 min (SL 40R, Thermo Scientific, Waltham, MA, USA), washed with the solvent mixture once, and the supernatant was retained. Finally, the extract was concentrated using a rotary evaporator at 45 • C and 450 mbar (R-210, BUCHI, Flawil, Switzerland). The extracts were frozen at −20 • C until analysis [13,20].

Determination of Total Phenols
Folin-Ciocalteu (F-C) spectrophotometric method was performed [16]. Quantification of the phenolic compounds was performed with gallic acid (mg GAE/g dw) as follows: 50 µL of each sample was mixed in 3 mL of distilled water, and then 250 µL of F-C reagent (2N) and 750 µL of 7% (w/v) sodium carbonate were added. Later, the mixture was allowed to incubate for 8 min at room temperature (RT). Finally, 950 µL of distilled water was added, followed by homogenization and incubation in the dark for 2 h at RT, and the absorbance was measured at 765 nm (Infinite M200 Pro, TECAN, Männedorf, Switzerland).

Determination of Flavonoids
Aluminum chloride (AlCl 3 ) test was carried out [16,18]. Quantification of the flavonoid content was expressed as quercetin equivalents per gram of bean flour (mg QE/g dw) by mixing 50 µL of each sample, 700 µL of deionized water, and 250 µL of AlCl 3 [133 mg of AlCl 3 and 400 mg of sodium acetate in 100 mL of methanol-water-acetic acid (140:50:10 v/v/v)]. Then, the solution was mixed and kept at RT for 30 min for completion of the reaction. Finally, the absorbance was measured at 410 nm.

Determination of Condensed Tannins
Vanillin test was used to quantify the condensed tannins based on the catechin equivalents per gram of bean flour (mg CAE/g dw). Briefly, 500 µL of each sample, 3 mL of 4% (w/v) vanillin methanol solution, and 1.5 mL of 50% (v/v) sulfuric acid solution were mixed, gently stirred and incubated at RT for 15 min in the dark [36,37]. Later, the absorbance was measured at 550 nm.

Determination of Anthocyanins
This analysis was carried out from the bean seed coat using AOAC method 2005.02 [6]. One hundred milligrams of bean seed coat was weighed and extracted with 2.4 mL of acidified ethanol (1N HCl, 85:15 v/v) for 15 min. Then, the volume of the samples was adjusted to 5 mL with acidified ethanol and centrifuged at 3488× g for 15 min, and the supernatant was recovered to carry out the corresponding analyses. The absorbance at 520 and 700 nm was measured with a spectrometer (Infinite M200 Pro, TECAN, Männedorf, Switzerland). Finally, the following equation was used to calculate the concentration of anthocyanins as mg of cyanidin 3-glucoside per liter.

Quantification of Phenolic Compounds by UHPLC
Bean flour was weighed at 0.5 g and extracted into 10 mL of ethanol with 1% trifluoroacetic acid (TFA) and refrigerated at 4 • C for 24 h. Then, the samples were centrifuged (3,000 rpm) for 5 min, and the supernatant was collected. Finally, the residue was washed once and centrifuged under the same conditions; both supernatants were mixed for analysis by liquid chromatography (UHPLC).
UHPLC analysis was performed with a PDA Detector (Acquity Arc, Waters, Milford, MA, USA), and the gradient conditions started with 90% solvent A, followed by 82% A at 10 min, 72% A at 18 min, 60% A at 19 min, and 90% A at 23 min. Solvent A was water with 0.1% TFA, and solvent B was acetonitrile with 0.1% TFA, with a flow rate of 0.7 mL/min. The column used was a C18, 2.7 µm, 4.6 × 150 mm (Cortects, Waters, Milford, MA, USA), the oven was set at 30 • C, and the volume injected was 20 µL. The compounds quantified in this study were detected at different wavelengths; ferulic acid was analyzed at 250 nm, kaempferol 3-glucoside, quercetin 3-glucoside, and genistein were analyzed at 300 nm, and cyanidin 3-glucoside was analyzed at 550 nm. All samples were injected in duplicate [24]. This analysis was performed by mixing 20 µL of each sample and 270 µL of the ABTS + radical. A solution of 7 mM ABTS and 2.5 mM potassium persulfate was used, and the reaction was allowed to stabilize for 12 to 16 h in the dark at RT. The ABTS + solution was diluted in ethanol until an absorbance of 0.70 ± 0.02 was obtained at 734 nm (Infinite M200 Pro, TECAN, Männedorf, Switzerland). Reaction blanks were prepared with an ethanol solution and the control sample was deionized water (20 µL) and the ABTS + radical solution (270 µL). The absorbance was measured at 734 nm, 6 min after the initial mixing until the reaction was complete, which was 24 min after starting the reaction [28].

DPPH Assay (2, 2 Diphenyl-1-picrylhydrazyl)
This test was carried out by mixing 100 µL of each sample and 100 µL of the DPPH* radical. The absorbance was measured at 515 nm every 5 min after the beginning of the reaction until 30 min had elapsed [28].
Both ABTS + (ARSA %) and DPPH* (DRSA %) radical scavenging activity was measured following this equation: where Am = sample absorbance (sample and ABTS + radical solution or sample and DPPH* radical solution), Ab = blank absorbance (sample and ethanol or sample and methanol), Ac = control absorbance (deionized water and ABTS+ radical solution or deionized water and DPPH* radical solution).

Phytic Acid Assay
This assay was performed enzymatically from bean flour, using the commercial phytic acid (phytate)/total phosphorus Megazyme kit (K-PHYT, Megazyme, Wicklow, Ireland), according to the instructions provided by the supplier.

Oligosaccharide Assay
This assay was carried out enzymatically from bean flour, using the raffinose/sucrose/glucose Megazyme kit (K-RAFGL, Megazyme, Wicklow, Ireland), following the instructions provided by the supplier.

Determination of Lectin Activity
The bean flour extract used to determine lectin activity (LA) was prepared following the methodology proposed by previous authors [38], with some modifications. Phosphate buffer (PBS, pH 7.4) was used as a buffer solution with the flour in a ratio of 1:13 (w/v). The solution was left under continuous stirring for 18 h. The mixture was centrifuged at 12,000× g for 20 min, obtaining the supernatant for the agglutination test.
The erythrocyte preparation was carried out with human blood type A + in a heparin tube, taking 1 mL of blood and adding 1 mL of PBS. The erythrocytes were washed twice with the same solution and centrifuged (847× g for 8 min). Once the erythrocytes were obtained, a 4% (v/v) solution was prepared.
The agglutination test was carried out in 96-well ELISA plates, in which 50 µL of PBS and 50 µL of the bean flour extract sample were placed, and two-fold dilutions were made. Afterward, 50 µL of erythrocyte solution was added to each well and the mixture was allowed to stand at controlled RT for 1 h.
Controls were a 0.2% Phaseolus vulgaris L. lectin solution, used as a positive control, and a PBS solution, used as a negative control. The lectin activity was expressed as the agglutination activity, which was the inverse of the maximum dilution at which agglutination was observed in mg of protein (HAU/mg protein).

Determination of Trypsin Inhibitory Activity
Preparation of the bean sample was performed according to Kakade et al. [39]. One gram of raw bean flour was weighed, and 50 mL of sodium hydroxide (NaOH) was added followed by stirring for 3 h to ensure maximum extraction. Then, the solution was adjusted to pH 8.2 and centrifuged at 3,000× g for 10 min, and the precipitate was discarded.
The trypsin activity methodology was performed using the commercial trypsin activity colorimetric kit (MAK290-1KT, Sigma-Aldrich, St. Louis, MO, USA), with some modifications. The reagents were prepared according to the instructions of the kit, where the standard curve of p-nitroaniline (p-Na) was used at five concentrations (2, 4, 6, 8, and 10 nM). The positive control with trypsin binds to the substrate to generate p-Na, and the color intensity is proportional to the content of p-Na, allowing for measurement of the trypsin activity. Bean samples containing trypsin inhibitors may inhibit the formation of p-Na. One microliter of trypsin solution was added to the samples. The control was the sample without the addition of any inhibitor (100% trypsin formation).
Trypsin activity was calculated, and trypsin activity formation was compared with the control to determine the TIU, which is the amount of inhibitor that is required to inhibit one microgram of pure trypsin. The results are reported as TIU per g of sample.

Statisdical Analysis
Statistics performed include the average, standard deviation, and analysis of variance (ANOVA) with a confidence level of α = 0.05. In addition, correlation studies were performed with Pearson's correlation coefficients and discriminant and cluster analyses to reveal the relationship and differences between the species of beans collected. To assess the difference between mean values, Tukey's test was performed with a p < 0.05 significance level. The statistical software Statgraphics Centurion XVI (Version 16.1.03) (Statgraphics Technologies, The Plains, VA, USA) was used.

Conclusions
In all species collected, a strong correlation was obtained between the antioxidant activity (DPPH) and the content of flavonoids, anthocyanins, and lectin. Similarly, quercetin 3-glucoside and genistein were the main phenolic compounds found in almost all the accessions analyzed.
Furthermore, P. vulgaris accessions showed the highest concentration of phenolic compounds and antioxidant activity. Two accessions of P. vulgaris (CH-01 and CH-12) showed a similar composition of phenolic compound content and antioxidant activity as the accessions of P. lunatus and P. polyanthus. Moreover, two accessions of P. vulgaris (CH-01 and CH-14) presented similar compositions of antinutritional content and could be grouped with the accessions of P. lunatus and P. coccineus. However, it was observed that there is a clear division between species (100%) based on their phenolic compounds, protein content, antioxidant activity, and antinutritional components. Therefore, it can be assumed that accessions of different bean species have a high biological potential for the prevention of diseases such as cancer, heart disease, the immune system, and obesity, among others. In addition, this study aimed to encourage and expand human consumption of bean species that are currently underutilized, as well as to promote their production and conservation.