Next Article in Journal
Research on Weed Reverse Detection Methods Based on Improved You Only Look Once (YOLO) v8: Preliminary Results
Next Article in Special Issue
Hidden Secrets of Mangrove Swamp Rice Stored Seeds in Guinea-Bissau: Assessment of Fungal Communities and Implications for Food Security
Previous Article in Journal
Integrated Ex-Situ Conservation and Ornamental Evaluation of the Vulnerable and Protected Greek Endemic Campanula laciniata L.: A Multifaceted Approach
Previous Article in Special Issue
The Role of Crop Wild Relatives and Landraces of Forage Legumes in Pre-Breeding as a Response to Climate Change
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diversity of Common Bean (Phaseolus vulgaris L.) and Runner Bean (Phaseolus coccineus L.) Landraces in Rural Communities in the Andes Highlands of Cotacachi—Ecuador

by
María José Romero-Astudillo
1,2,*,
César Tapia
3,
Joaquín Giménez de Azcárate
4 and
Diego Montalvo
5
1
Agrobiodiversity and Food Security Research Group—GIASSA, Agricultural and Environmental Science Faculty, Universidad Técnica del Norte, Ibarra 100105, Ecuador
2
Doctoral Programme in Agriculture and Environment for Development, Campus Terra, Universidad de Santiago de Compostela, 27002 Lugo, Spain
3
Santa Catalina Experimental Station, Instituto Nacional de Investigaciones Agropecuarias (INIAP), Quito 170201, Ecuador
4
Department of Botany, Campus Terra, Universidad de Santiago de Compostela, 27002 Lugo, Spain
5
Facultad Latinoamericana de Ciencias Sociales, FLACSO, Quito 170201, Ecuador
*
Author to whom correspondence should be addressed.
Submission received: 6 June 2024 / Revised: 28 June 2024 / Accepted: 29 June 2024 / Published: 29 July 2024
(This article belongs to the Special Issue Seeds for Future: Conservation and Utilization of Germplasm Resources)

Abstract

:
Phaseolus species are cultivated worldwide as a primary food source for human consumption. Common bean (Phaseolus vulgaris L.) and runner bean (Phaseolus coccineus L.) landraces are often cultivated together. The purpose was to document the traditional knowledge held by farmers regarding bean landraces, to describe the diversity through seed morphological descriptors (five quantitative and five qualitative traits), to understand the geographic distribution using Moran’s I statistic, and to analyze the diversity through the Shannon–Wiener Diversity–Equity (H′) index and the Shannon J’ homogeneity index at the community level to better understand the richness of this geographic area. A total of 361 common and runner bean accessions were collected in 10 Andean rural communities of Cotacachi, North Ecuador. We identified 47 landraces, predominantly kept by indigenous female farmers, with limited access to formal education, highlighting the role of this disadvantaged population in agrobiodiversity conservation. The cluster analysis revealed three groups (cophenetic correlation coefficient = 0.6). In the principal component (PC) analysis, 94% of the variation was explained by PC1 and PC2. A positive spatial autocorrelation (Moran’s I: 0.24; z-score: 2.20; p-value: 0.03) was identified, confirming a spatial structure. The Morochos community showed the highest diversity (H′ = 1.55). The information of the diversity and distribution of common and runner bean variability provided in this study is a contribution to further research focused on conservation.

1. Introduction

The Fabaceae family is an important contributor to food systems [1], and it provides 33% of the human dietary protein consumption [2,3]. It has about 770 genera and 19,500 species [4,5]. There are genera within the family that are important food sources for human feeding. The genera Phaseolus, Pisum, Vicia, Cicer, Glycine, Lens, Cajanus, Vigna, Arachis, and Lupinus are considered an alternative to animal protein and can improve the daily diet of families in most areas of the world [3,5,6,7]. The genus Phaseolus includes 75 infrageneric taxa [8]. There are species within the Phaseolus genus such as P. vulgaris, P. lunatus, P. coccineus, P. dumosus, and P. acutifolius that are cultivated around the world as a source of food for human consumption [9,10,11,12]. Common bean (Phaseolus vulgaris L.) and runner bean (Phaseolus coccineus L.) are usually cultivated together [13]. These two species are very closely related and are partially inter-crossable when P. vulgaris is the female parent. Common bean is an annual plant predominantly autogamous, while runner bean is mainly allogamous and considered a perennial plant [13,14].
Common bean is the most planted legume in the world [1,15]. It is estimated that the area harvested in 2022 was about 36,792,490 ha [16]. It involves wild and domesticated forms [17]. It has a significant nutritional value and health properties [14,18,19,20], as well as socioeconomic importance [21,22,23]. Furthermore, it is one of the most studied and ex situ conserved species worldwide [24], as well as one of the most cultivated species in traditional farming [25].
There have been very few studies into the origins of P. coccineus. As mentioned by Bitocchi et al. [14], two studies have addressed the domestication of this species. On the other hand, research on the history of P. vulgaris suggests that there was a divergence in wild bean populations approximately 165,000 years ago, originating from a Mesoamerican ancestral population [26,27,28]. Common bean was domesticated after this divergence [29] about 8000 years ago [27,30] independently in two main geographic locations: Mesoamerica and the Andes, deriving two major eco-geographically distinct domesticated gene pools [22,27,31,32,33,34,35,36]. Nevertheless, important gene exchange between these two gene pools and genetic material movement have occurred since domestication [15]. In general, Mesoamerican genotypes show small- or medium-sized seeds, while the Andean genotypes are typically larger [37]. Due to the process of domestication, a wide number of varieties have been generated [18,35].
Traditional varieties, known as landraces or farmers’ varieties, are an important element of agrobiodiversity [1]. Crop landraces are local varieties with agricultural, cultural, and historical value, usually genetically diverse [38,39], that have been adapted to a specific ecosystem through a selection process developed by farmers [15]. Landrace diversity is important and should be considered as a supply of genetic resources for crop breeding [1]. Farmers usually give these landraces different names based on criteria such as growth habit, phenology, seed morphology, and use [40].
In this context, traditional-based agriculture (small-scale production) is an important source of food from the perspective of cultural identity and values. Usually managed by households [41], it allows to conserve crop genetic diversity, as well as to support ecosystem functioning, resilience, and productivity [42,43]. The conservation of crop genetic resources is important for local and global food security, as well as to face population growth and global warming [44,45]. That is why landrace in situ conservation should be considered to improve the agroecosystem structure [43].
Unfortunately, the world is facing a considerable rate of agrobiodiversity loss [38,46,47]. There is an increasing reduction in genetic variability due to intensive anthropic interventions [8], the progressive substitution of landraces with high-yield modern varieties by farming communities [43], land use change, and habitat degradation [16,44], as well as cultural transformation and diet changes [23,48].
That is why understanding the dynamics of plant diversity within a specific area and its cultural environment is basic for the conservation of agrobiodiversity. By increasing the knowledge of these factors, we can better appreciate the importance of preserving the unique diversity and cultural heritage of a region. Consequently, inventories play an important role in obtaining information about farmers’ varieties, which can then be used to promote conservation and reduce the risk of genetic erosion [1,17]. Even when a genetic diversity assessment is a fundamental principle of crop breeding programs and plays an important role in germplasm conservation [49], by documenting these varieties, we can better understand their diversity and characteristics, leading to improve conservation efforts and sustainable agricultural practices. Furthermore, the collection and characterization of bean landraces are important to prevent the loss of unknown genetic resources [50].
In this context, we decided to work on this research as a contribution to better understand the bean diversity of an area known for its traditional-based agriculture. The area we selected is Cotacachi, in the north of Ecuador, which is considered an important center of agrobiodiversity [51,52]. Being crossed by the Andes Mountain range, it has several altitudinal belts, which favor climatic diversification, responsible for the development of extensive agricultural biodiversity. The agricultural calendar is regulated by the rainy seasons, which are characterized by higher levels of rainfall, occurring in March and April and October and November, with average rainfall of 111 mm/month and 96 mm/month, respectively [53]. In contrast, the dry season takes place during the months of July and August, with an average precipitation of 16 mm/month [53], accompanied by intense solar radiation and strong winds [54]. Throughout the remainder of the year, the area experiences moderate rainfall. The climate is high Andean equatorial, with constant temperatures throughout the year. At 2500 m.a.s.l., the average annual temperature is around 13.5 °C, and the average annual rainfall is approximately 820 mm, with a notable scarcity of water for agricultural use in the dry season [53,55].
The main source of food and economic income in the study area is family farming [56,57]. Farming is conducted starting at approximately 2300 m.a.s.l. in the Inter-Andean Valley and extended upwards along the Cotacachi Volcano slopes, until reaching an elevation of approximately 3300 m.a.s.l. [58]. Their farms have a surface area of less than 1 ha [59]. This agriculture system is locally called chakra (a Kichwa term, Kichwa words written in italic), and it is defined as an ancestral agricultural system of the Kichwa indigenous peoples in the Ecuadorian Andean region, as well as in other areas of Perú and Bolivia [60]. It is a small plot next to the house displaying a wide range of agrobiodiversity where a variety of traditional agricultural practices are applied [61]. The traditional practices of seed selection and on-farm savings are still very common in these indigenous communities [62]. The water supply is usually limited; thus, about 55% of families do not have access to water for irrigation purposes [63].
Therefore, the aims of this research were to record farmers’ knowledge on bean landraces (P. vulgaris and P. coccineus) and to describe its diversity through morphological characterization of the seeds to provide a view of the available genetic resources of the species in Cotacachi as a contribution to better understanding the agrobiodiversity richness in the area.

2. Materials and Methods

2.1. Study Area

This research was conducted in three administrative zones: Imantag, El Sagrario, and Quiroga (Figure 1) in the Andes Highlands (altitude >2500 m.a.s.l.) of Cotacachi Canton. It is located in the north of Ecuador, approximately 100 km from the capital, Quito.
Cotacachi has one of the highest concentrations of indigenous populations (41%) in the country [64,65]. The population is mainly bilingual (Spanish and Kichwa) [66]. Most of the communities in the high Andean area are part of Union of Peasant and Indigenous Organizations of Cotacachi (UNORCAC), which is a second-tier, non-profit organization created in 1977. UNORCAC is made up of 41 communities and various indigenous organizations. Their members have strong cultural and territorial identity and usually revalue and preserve ancestral knowledge and traditional practices [63].

2.2. Sampling

We randomly selected 87 smallholder farmers (Table 1), which represent 10% of the total farmers from 10 communities (out of 41 communities that are part of UNORCAC). These communities were chosen based on their focus on family farming and its role as traditional agriculture management representative examples in the area. Communities not centered in agriculture were not considered in this study.
We divided the research into three parts. The first one involved farmers’ information on the diversity of beans, the second part included seed morphological characterization conducted on 361 bean accessions collected in the chakras of the selected farmers, and the third part was related to the geographical distribution of the landraces and calculation of diversity indices at the community level.

2.3. Data Collection

2.3.1. Farmers’ Information

We conducted semi-structured interviews with 87 smallholder farmers to obtain information about their agricultural diversity. This method has progressively been used to document farmers’ knowledge related to their crop landraces and traditional agricultural practices [67,68]. Before the interview, informed consent was obtained from each farmer. The interview included information about the names of the landraces as farmers know them, eliminating synonymy (different names referring to the same landrace) and homonymy (one name applied to two or more different landraces). We used farmer-named varieties, because it has been stablished that it reflects farmers’ own understanding of their crop diversity [69]. At the same time, we focused on the perception of farmers about the landraces they believe are in risk of disappearing. The interview had three sections, as shown in Table 2.
The interviews took place in the chakras, and they were conducted mostly in Spanish. In a few cases, a Kichwa interpreter was required.

2.3.2. Morphological Characterization

Seed samples of bean landraces were collected in chakras from 10 communities (Table 1) between June 2019 and January 2020. Each sample had approximately 50 seeds, although, in some landraces, seeds were scarce (even less than 10 seeds/sample). The samples were sealed and labeled for further analysis. We collected a total of 361 accessions of different bean landraces.
In order to find intraspecific diversity, our study was focused on seed morphology through qualitative and quantitative descriptors, as other studies have shown that seed descriptors are the most discriminant for common bean [23,35,70,71,72]. Moreover, characterization has several purposes, including the provision of statistical data on the analyzed germplasm and the identification of any absent genotypes within a specific collection [73].

Qualitative Descriptors

We used five seed qualitative descriptors: primary color, secondary color, color around hilum, coat pattern, and shape.
Primary and secondary seed coat color was assessed in the dry grain by direct observation under natural daylight, establishing the primary and secondary colors using a modified scale from the International Board for Plant Genetic Resources [74]. Additionally, the hilum-surrounding color was defined, as it is considered an essential part of seed appearance [72]. To assess this trait, we used a modified scale [75].
The seed coat pattern was determined through the classification proposed by the International Board for Plant Genetic Resources [74]. To assess seed shape, we used a classification in three categories: oval-, cuboid-, and kidney-shaped, as applied in previous studies [71,76].

Quantitative Descriptors

We selected five seed quantitative descriptors previously used by Sinkovič et al. [35]. By using a digital caliper, we measured in millimeters (mm) 10 fully developed unharmed seeds in each sample in order to define the seed length (L) [highest parallel to hilum], seed width (W) [hilum to opposite side], seed thickness (T) [highest perpendicular to hilum], seed length/width ratio (L/W), and seed width/thickness ratio (W/T).

2.3.3. Geographic Distribution and Diversity Analysis

The geographic distribution was analyzed with software ArcGIS v. 10.8 [77], using Universal Transverse of Mercator (UTM) Projection, World Geodetic System (WGS) 1984 Datum, and 17 South Zone. Spatial autocorrelation Moran’s I statistic, with the Euclidean distance method, was performed to explain the spatial relation of the number of bean landraces among the chakras [78].
The Biodiversity Professional version 2 statistical analysis program [79] was employed to analyze the diversity of P. vulgaris and P. coccineus. Calculation of the agrobiodiversity indices was conducted at the community level using the Shannon–Wiener Diversity–Equity (H′) index, as well as the Shannon J’ homogeneity index.

2.4. Data Analysis

Statistical analyses were performed using the software InfoStat v. 2008 [80]. Statistical parameters for the quantitative descriptors included the mean, minimum and maximum values, standard deviation (SD), and coefficient of variation (CV) to assess the variability of the traits.
In addition, principal component and cluster analyses were executed [81]. Principal component analysis was, mainly, used to explore relationships among variables [82]. Cluster analysis was performed using the Ward method with Gower distance on standardized variables. It allowed to build a dendrogram, clustering the studied landraces into different groups according to their characteristics [69]. All the quantitative and qualitative traits of each accession were used. Analysis of variance (ANOVA) was performed, and the Fisher’s Least Significant Difference (LSD) method was used to compare means (α = 0.05).

3. Results

3.1. Farmers’ Information

The results show that 89% of the farmers identify themselves as indigenous, while the 11% remaining considered themselves as mestizo. Concerning the language, most farmers were bilingual (Kichwa and Spanish), although, in six cases, we needed an interpreter, because they only spoke Kichwa. Most of the farmers were female (63%), with an average age of 52 years old. In terms of education level, 41.4% of the farmers did not receive any formal education, 40.2% attended primary (elementary) school, 14.9% secondary education, and only 3.4% had access to higher education. Among the respondents, 4.6% rent the land they use to produce crops, while 8% use community land, and most of the farmers (87.4%) reported that they own their land. When considering the aspect of water accessibility, a significant majority of farmers (87.4%) mentioned that they do not have access to irrigation water and depend on weather conditions.
Regarding the perception of landraces at risk of disappearing, 26% of the farmers reported that Yana Vaca Poroto was the landrace they considered to be at risk of disappearing. However, we found this landrace in 10 of the 87 chakras. In contrast, we recorded two landraces (Poroto Campeon and Sucu Tomate Pintado) that were scarce in the study area, as they were found in only one chakra each.

3.2. Bean Landrace Diversity

Our study identified 47 landraces of bean species (Figure 2), from which 40 were common bean landraces (P. vulgaris) and 7 were runner bean landraces (P. coccineus). Common bean was usually planted in the middle of the chakras, while runner bean was planted along the borders. Locally, the generic name for common bean is poroto and, for runner bean, is popayan. As in previous research [40,83], we found that common names given to beans by local farmers in the Kichwa language are mainly based on seeds’ color and morphological characteristics. In our study, local names of the landraces seem to be related specially to the seed color, pattern, and shape, as described in other studies [84].
The results of this research show a wide variability in the morphological characteristics of bean seeds. This diversity is usually common in traditional varieties that have faced less selection processes than the commercially developed ones [20]. These results are consistent with the findings of other studies [18,71]. Freitas et al. [18] assessed 50 accessions of common bean, representing the diversity cultivated on the Island of Madeira, an area with a wide range of ecological conditions. The common bean materials evaluated in that study showed a broad intraspecific variability. It is important to mention that those materials had, presumably, their origin in the Andean gene pool, based on their results of morphological traits. On the other hand, Loko et al. [71] evaluated 57 common bean accessions from the Benin Republic based on seed morphological traits. As expected, they found a wide range of variability among landraces as well.
Landraces with larger populations were Canario (present in 31 out of 51 chakras) and Poroto Grande (28 out of 51 chakras). In contrast, Poroto Campeon and Sucu Tomate Pintado were found to have limited populations, being restricted to a single side each (Table 3). Farmers expressed their preference for cultivating Canario, especially motivated by commercial purposes and the consumer demand it generates. It is important to mention that Canario is not a landrace; rather, it is an improved variety.

3.3. Morphological Characterization

Concerning the primary color of the bean seed, we found a wide variability among landraces, which varied from light to dark colors. The most frequently observed were pale cream to buff (23%), white (16%), and orange (13% of the 361 samples collected), with smaller percentages found in black and gray colors. As mentioned by Abera et al. [1], seed color is the main morphological characteristic used to identify bean landraces. This attribute not only allows a better communication between researchers and farmers but also promotes improved mutual understanding.
Regarding the secondary seed color, 41% of the samples showed no color, 29% had black as the secondary color of the seed, and 7% of the samples showed brown and red colors. Brown was the most common color around the hilum (47% of the samples). A total of 41% of the samples had no seed coat pattern. Regarding the seed shape, the most frequent one was kidney-shaped (51%).
The mean, standard deviation, coefficient of variation, and minimum and maximum values of three quantitative variables can be found in the Supplementary Materials. The results reveal that seed thickness exhibits greater variability (CV = 15.21) compared to the other morphological traits, as evidenced by its higher coefficient of variation.
In the Supplementary Materials are presented summary statistics for the five quantitative bean traits. Our results show that the mean values for the attribute seed length (L) varied from 10.3 mm (Yura Allpa Poroto) to 19.5 mm (Paco Bolon). The seed width (W) values ranged from 4.2 mm (Sucu Allpa Poroto) to 13.6 mm (Paco Bolon), and seed thickness (T) values varied between 3.9 mm (Matambre Negro) and 9.5 mm (Canario). The lowest seed length/width ratio (L/W) was 1.1 mm (Canario), while the highest ratio registered was 3.3 mm (Sucu Allpa Poroto). Regarding the seed width/thickness ratio (W/T), the values differed from 0.8 mm (Sucu Allpa Poroto) to 2.2 mm (Paco Bolon).
The coefficient of variation for the variable seed length (L) was higher for the landrace Sucu Pintado. This landrace presented the greatest variability in the characteristic seed thickness (T) as well. Regarding the variable seed width (W), the Sucu Allpa Poroto landrace had the highest heterogeneity, which is considered a high intra-varietal variability.
The results of the cluster analysis revealed that the 47 landraces could be clustered into three groups (Figure 3). The cophenetic correlation coefficient was 0.6.
The first group (G1) represented by 22 landraces showed seed average dimensions as follows: L = 13.84 mm, W = 10.04 mm and, T = 7.41 mm, which were the highest within the three groups; the aspect ratio values, L/W and W/T, were 1.38 and 1.37, respectively; the accessions in the group displayed mainly white seed primary color, no secondary color, brown color around the hilum, no pattern, and an oval shape. The seeds of G1 were the largest and roundest of all the groups.
The second group (G2) contained six landraces with the following seed average dimensions: L = 13.51 mm and W = 8.27 mm, which were the lowest among the groups, and T = 6.20 mm, L/W = 1.61, and W/T = 1.34; the seeds in this group exhibited a primary color ranging from pale cream to buff, generally without a secondary color or pattern, brown color around the hilum was predominant, and most of the seeds were kidney-shaped.
The accessions within the third group (G3) included 19 landraces, characterized by the following seed dimensions: L = 13.69 mm, W = 8.29 mm, T = 6.06 mm—this last one was the lowest of all three groups—L/W = 1.67, and W/T = 1.39; the main seed primary color was pale cream to buff, black was the principal secondary color and color around the hilum, and most of the accessions showed a striped pattern and kidney shape.
Principal component analysis showed that 94% of the variation was explained by PC1 and PC2 (Figure 4 and Table 4). The first principal component (PC1) explained the highest percentage of variance (51%), and it was correlated with seed thickness (T) and the ratio L/W. The second principal component (PC2) justified 43% of the variation and was highly correlated with the seed length (L) and the ratio W/T. Most of the landraces of G1 were located on the right side of the plot (red circles). The P. coccineus landraces were concentrated in the upper right quadrant (blue circles), exhibiting a distinct boundary from the other groups. Landraces from G2 and G3 were placed on the left side of the plot (green circles). There were also outliers in the group, represented by landraces Poroto Campeon, Sucu Allpa Poroto, Yana Pintado, and Rayado Poroto.
Significant differences (p < 0.05) were detected on the traits of seed width (W), seed thickness (T), and the L/W ratio (Table 5). On the other hand, we did not find statistical differences on the seed length (L) trait, nor the W/T ratio.

3.4. Geographic Distribution

As shown in Figure 5, we found bean landraces distributed in the three administrative zones. The number of landraces in each chakra ranged from 1 to 6 (26%), 7 to 12 (24%), and 13 to 18 (8%).
As seen in Figure 6, we found a positive spatial autocorrelation (Moran’s I: 0.24; z-score: 2.20; p-value: 0.03). Since the index value is positive, we assume that the number of landraces is associated with close locations [85]. Given the z-score of 2.20, there is a less than 5% likelihood that this clustered pattern could be the result of random chance, which confirms the presence of a spatial structure.

3.5. Diversity Analysis

The findings show that the values of the Shannon index (H′) vary within a range of 0.45 to 1.55 (Table 6). The Morochos community exhibits the highest P. vulgaris and P. coccineus diversity among the communities.
In terms of homogeneity (Shannon J’), the Chilcapamba community exhibits the highest value (0.99), suggesting that the P. vulgaris and P. coccineus landraces in this community are more similar to each other. This implies a reduced diversity within the community.

4. Discussion

We found a wide variability of beans (Figure 2) distributed in the study area (Figure 5) (46 landraces and a commercial variety) regularly linked to traditional reasons. This diversity was predominantly kept by indigenous female farmers with limited access to formal education. Their agriculture depends on weather conditions, as most of them do not have access to irrigation water. Considering that bean diversity in the study area is conserved by a disadvantaged population, this study highlights the need for public policies that improve their living conditions and, consequently, indirectly favor agrobiodiversity conservation.
The Canario variety was the most abundant within the studied chakras (Table 3). It was geographically distributed in the three administrative zones. The relatively high abundance of this variety might be explained by its commercial purpose. Most of the interviewees stated that they chose to cultivate this particular landrace due to its high demand among consumers. On the other hand, the second-most abundant landrace was Poroto Grande. Similar to Canario, this landrace was distributed throughout the study area. Coincidentally, both share a similar color tone of their seeds. Canario was classified as light yellow, while Poroto Grande showed a medium yellow brown color. This suggests that yellow tones represent a desirable trait among farmers and consumers. Other authors have reported a preference for yellow tone common bean as well [86,87,88]. The interest in yellow beans can be attributed to their better digestibility characteristics, as well as the less intestinal distress they cause [87,89]. As a complement, it is worth mentioning that, although probably not linked to consumers’ preference, research has revealed that yellow bean varieties have a higher protein content compared to red seeds [90]. Being a commercial variety, numerous farmers cultivate it for trade due to their preferences and its superior yield. Consequently, this could represent a high risk of genetic erosion, as landraces and other low-yielding cultivars might eventually be displaced. This finding leads to exploring strategies for diversifying the production and promoting the cultivation of traditional varieties.
Consistent with previous studies [91,92], our findings also highlight the significance of seed shape and color combination in influencing the preferences of both farmers and consumers. The significant diversity registered in seed color suggests that farmers’ preferences probably influenced the distribution and conservation of these landraces based on seed colors.
Our findings suggest that Poroto Campeon and Sucu Tomate Pintado landraces are the ones that are, in fact, facing a substantial risk of disappearing. This result also implies that farmers’ perception may not always be entirely accurate, and further information should be considered to contrast the information given by farmers.
The cophenetic correlation coefficient of 0.6 in the cluster analysis indicates reasonably well-defined clusters based on the morphological seed traits within the three groups (Figure 3). The landraces in G1 were the most abundant (22), and their seeds showed the largest dimensions of all the groups, exhibiting an oval shape and mainly white seed primary color with no pattern. The seeds in G2 represented the smallest group (six landraces) and showed the shortest length and width dimensions, pale cream to buff primary color generally without a pattern, and most of the seeds were kidney-shaped. The seeds in G3 included 19 landraces and presented the lowest thickness of all the groups, the main seed primary color was pale cream to buff, and most of the accessions showed a striped pattern (with black as the principal secondary color) and kidney shape. These results could lead to selecting landraces for future breeding programs to obtain desirable characteristics based on the consumer preferences.
The principal component analysis revealed that most of the landraces of G1 were located on the right side of the plot (Figure 4). The P. coccineus landraces were concentrated on the upper right side, clearly separated from the other groups. In contrast, the landraces from G2 and G3 were located on the left side of the plot. The first principal component (Eigenvalue = 2.53) explained the highest percentage of variance (51%) and was correlated with the seed thickness (T) and the L/W ratio (Table 5). This result should be considered for further study due to its contribution to the total diversity.
The positive Moran’s index indicates that chakras with higher diversity tend to be distributed closer to each other, as do those with lower diversity (Figure 6). This clustering tendency suggests that there are factors that might be influencing this spatial distribution [93]. In this study, the factors potentially influencing this distribution may likely involve seed exchanges among neighbors.
The highest Shannon diversity index (H′) in this study was found in the Morochos community (H′ = 1.55) (Table 6). It was higher compared to the reports of [69] on common bean (H′ = 0.74). The high diversity found in this community makes it a potential site for both sourcing genetic resources and implementing conservation programs.
This diversity could be attributed to the traditional practices of indigenous communities in the study area, where landraces are usually kept in the chakras and seeds are frequently passed down from one generation to the next.

5. Conclusions

This study describes the morphological seed characteristics of 361 georeferenced samples of bean landraces cultivated in three administrative zones: Imantag, El Sagrario, and Quiroga in the Andes Highlands of Cotacachi Canton in Ecuador. Based on the seed morphological descriptors, as well as the information provided by farmers, we identified 46 bean landraces and 1 commercial variety. The landraces with a wider distribution and population throughout the study area were Canario and Poroto Grande, probably due to the preference by farmers and consumers for yellow common bean. Contrary to farmers’ perception, our findings revealed that landraces Poroto Campeon and Sucu Tomate Pintado had limited distribution and scarce populations in the study area. This observation raises awareness about the potential vulnerability of these landraces, representing a possible risk of their disappearance. The high diversity found in the Morochos community (H′ = 1.55) makes it a promising geographic area for sourcing germplasm for further study. The diversity in the studied common and runner bean landraces, as well as their distribution, was confirmed by morphological traits, diversity indices, spatial autocorrelation, and cluster and principal component analyses, and the results represent a contribution to further research focused on the conservation of bean diversity in the Andean area of Cotacachi.

Supplementary Materials

The following Supporting Information can be downloaded at https://www.mdpi.com/article/10.3390/agronomy14081666/s1: Table S1. Variations in the quantitative morphological traits: seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio.

Author Contributions

Conceptualization, M.J.R.-A., C.T. and J.G.d.A.; Methodology, M.J.R.-A., C.T., J.G.d.A. and D.M.; Investigation, D.M. and M.J.R.-A.; Validation, M.J.R.-A., C.T. and J.G.d.A.; Formal analysis, M.J.R.-A., C.T. and J.G.d.A.; Data curation, M.J.R.-A. and D.M.; Writing—original draft preparation, M.J.R.-A.; Writing—review and editing, M.J.R.-A., C.T., J.G.d.A. and D.M.; Visualization, M.J.R.-A., C.T., J.G.d.A. and D.M.; Supervision, C.T. and J.G.d.A.; Project administration, M.J.R.-A. and C.T.; Funding acquisition, M.J.R.-A. and C.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the International Treaty on Plant Genetic Resources for Food and Agriculture, grant number PR-268-Ecuador. Publication funding was provided by the Universidad Técnica del Norte.

Data Availability Statement

Acknowledgments

The authors acknowledge the farmers who kindly provided the information, as well as the seed samples, for this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Abera, B.; Berhane, M.; Nebiyu, A.; Ruelle, M.L.; McAlvay, A.; Asfaw, Z.; Tesfaye, A.; Woldu, Z. Diversity, use and production of farmers’ varieties of common bean (Phaseolus vulgaris L.; Fabaceae) in southwestern and northeastern Ethiopia. Genet. Resour. Crop Evol. 2020, 67, 339–356. [Google Scholar] [CrossRef]
  2. De la Rosa, C.; Lozano, L.; Castillo-Ramírez, S.; Covarrubias, A.A.; Reyes, J.L. Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes. Genome Biol. Evol. 2020, 12, 2355–2369. [Google Scholar] [CrossRef] [PubMed]
  3. Smýkal, P.; Coyne, C.J.; Ambrose, M.J.; Maxted, N.; Schaefer, H.; Blair, M.W.; Berger, J.; Greene, S.L.; Nelson, M.N.; Besharat, N.; et al. Legume Crops Phylogeny and Genetic Diversity for Science and Breeding. Crit. Rev. Plant Sci. 2015, 34, 43–104. [Google Scholar] [CrossRef]
  4. Catarino, S.; Brilhante, M.; Essoh, A.P.; Charrua, A.B.; Rangel, J.; Roxo, G.; Varela, E.; Moldão, M.; Ribeiro-Barros, A.; Bandeira, S.; et al. Exploring physicochemical and cytogenomic diversity of African cowpea and common bean. Sci. Rep. 2021, 11, 12838. [Google Scholar] [CrossRef] [PubMed]
  5. Shavanov, M.V. The role of food crops within the Poaceae and Fabaceae families as nutritional plants. IOP Conf.Ser. Earth Environ. Sci. 2021, 624, 012111. [Google Scholar] [CrossRef]
  6. Brezeanu, C.; Brezeanu, P.M.; Calara, M.; Ambarus, S.; Cristea, T.O. The Role of Food Legume Species in the Context of Sustainable Agriculture, Food Security, Agrobiodiversity, Conservation and Human Health. Sci. Papers Ser. B Hortic. 2021, 65, 382–391. [Google Scholar]
  7. Molares, S.; Ladio, A. The Usefulness of Edible and Medicinal Fabaceae in Argentine and Chilean Patagonia: Environmental Availability and Other Sources of Supply. J. Evid. Based Complement. Altern. Med. 2012, 2012, 901918. [Google Scholar] [CrossRef] [PubMed]
  8. Acosta-Diaz, E.; Hernandez-Torres, I.; Amador-Ramirez, M.D.; Padilla-Ramirez, J.S.; Zavala-García, F.; Baez-Gonzalez, A.D. Collection and characterization of wild species of Phaseolus (Fabaceae) in northeastern Mexico for ex situ conservation. Plant Ecol. Evol. 2015, 148, 119–127. [Google Scholar] [CrossRef]
  9. Alcázar-Valle, M.; Lugo-Cervantes, E.; Mojica, L.; Morales-Hernández, N.; Reyes-Ramírez, H.; Enríquez-Vara, J.N.; García-Morales, S. Bioactive Compounds, Antioxidant Activity, and Antinutritional Content of Legumes: A Comparison between Four Phaseolus Species. Molecules 2020, 25, 3528. [Google Scholar] [CrossRef]
  10. Bisht, Y.; Jugran, A.K.; Mishra, P.K.; Pareek, N.; Sharma, H.; Singh, A.V. Assessment of genetic diversity and population structure of Indian common bean accessions using microsatellite markers. Plant Genet. Resour. C 2023, 21, 182–193. [Google Scholar] [CrossRef]
  11. Myers, J.R.; Formiga, A.K.; Janick, J. Iconography of Beans and Related Legumes Following the Columbian Exchange. Front. Plant Sci. 2022, 13, 851029. [Google Scholar] [CrossRef]
  12. Neto, W.V.D.; Medeiros, A.M.; Carvalho, L.C.B.; Ferreira, C.D.; Lopes, A.C.D.; Gomes, R.L.F. Selection of Landraces ff Lima Bean for Family Agriculture. Rev. Caatinga 2022, 35, 137–147. [Google Scholar] [CrossRef]
  13. Sicard, D.; Nanni, L.; Porfiri, O.; Bulfon, D.; Papa, R. Genetic diversity of Phaseolus vulgaris L. and P. coccineus L. landraces in central Italy. Plant Breed. 2005, 124, 464–472. [Google Scholar] [CrossRef]
  14. Bitocchi, E.; Rau, D.; Bellucci, E.; Rodriguez, M.; Murgia, M.L.; Gioia, T.; Santo, D.; Nanni, L.; Attene, G.; Papa, R. Beans (Phaseolus ssp.) as a Model for Understanding Crop Evolution [Review]. Front. Plant Sci. 2017, 8, 722. [Google Scholar] [CrossRef]
  15. Ramirez-Villegas, J.; Khoury, C.K.; Achicanoy, H.A.; Mendez, A.C.; Diaz, M.V.; Sosa, C.C.; Debouck, D.G.; Kehel, Z.; Guarino, L. A gap analysis modelling framework to prioritize collecting for ex situ conservation of crop landraces. Divers. Distrib. 2020, 26, 730–742. [Google Scholar] [CrossRef]
  16. Food and Agriculture Organization of the United Nations. FAOSTAT Statistical Database; FAO: Rome, Italy, 1997. [Google Scholar]
  17. Kouam, E.B.; Ndomou, M.; Gouado, I.; Pasquet, R.S. Assessment of the Genetic Diversity of Cultivated Common Beans (Phaseolus vulgaris L.) from Cameroon and Kenya Using Allozymes Markers. J. Exp. Biol. Agric. Sci. 2017, 5, 87–97. [Google Scholar] [CrossRef]
  18. Freitas, G.; Ganança, J.F.T.; Nóbrega, H.; Nunes, É; Costa, G.; Slaski, J.J.; de Carvalho, M.A.P. Morphological evaluation of common bean diversity on the Island of Madeira. Genet. Resour. Crop Evol. 2011, 58, 861–874. [Google Scholar] [CrossRef]
  19. Kumar, S.; Singh, M.; Malhotra, N.; Blair, M.W.; Sharma, J.P.; Gupta, R. Introgression of anthracnose resistance into the background of locally adapted common bean landraces. Euphytica 2021, 217, 52. [Google Scholar] [CrossRef]
  20. Lázaro, A.; Villar, B.; Aceituno-Mata, L.; Tardío, J.; De la Rosa, L. The Sierra Norte of Madrid: An agrobiodiversity refuge for common bean landraces. Genet. Resour. Crop Evol. 2013, 60, 1641–1654. [Google Scholar] [CrossRef]
  21. De Lima, M.S.; Carneiro, J.E.D.; Carneiro, P.C.S.; Pereira, C.S.; Vieira, R.F.; Cecon, P.R. Characterization of genetic variability among common bean genotypes by morphological descriptors. Crop Breed. Appl. Biotechnol. 2012, 12, 76–84. [Google Scholar] [CrossRef]
  22. Delfini, J.; Moda-Cirino, V.; Ruas, C.d.F.; Neto, J.d.S.; Ruas, P.M.; Buratto, J.S.; Ruas, E.A.; Gonçalves, L.S.A. Distinctness of Brazilian common bean cultivars with carioca and black grain by means of morphoagronomic and molecular descriptors. PLoS ONE 2017, 12, e0188798. [Google Scholar] [CrossRef]
  23. Scarano, D.; Rubio, F.; Ruiz, J.J.; Rao, R.; Corrado, G. Morphological and genetic diversity among and within common bean (Phaseolus vulgaris L.) landraces from the Campania region (Southern Italy). Sci. Hortic. 2014, 180, 72–78. [Google Scholar] [CrossRef]
  24. Pipan, B.; Meglič, V. Diversification and genetic structure of the western-to-eastern progression of European Phaseolus vulgaris L. germplasm. BMC Plant Biol. 2019, 19, 442. [Google Scholar] [CrossRef]
  25. Carvalho, M.S.; de Oliveira Moulin Carias, C.M.; Silva, M.A.; da Silva Ferreira, M.F.; Posse, S.C.P.; Guilhen, J.H.S.; Ferreira, A. Microsatellites and agronomic approaches reveal the diversity of beans (Phaseolus vulgaris L.) cultivated in Espírito Santo—Brazil, by family farms. Genet. Resour. Crop Evol. 2023, 70, 731–747. [Google Scholar] [CrossRef]
  26. Bitocchi, E.; Nanni, L.; Bellucci, E.; Rossi, M.; Giardini, A.; Zeuli, P.S.; Logozzo, G.; Stougaard, J.; McClean, P.; Attene, G.; et al. Mesoamerican origin of the common bean (Phaseolus vulgaris L.) is revealed by sequence data. Proc. Natl. Acad. Sci. USA 2012, 109, E788–E796. [Google Scholar] [CrossRef]
  27. Gaut, B.S. The complex domestication history of the common bean. Nat. Genet. 2014, 46, 663–664. [Google Scholar] [CrossRef]
  28. MacQueen, A.H.; Khoury, C.K.; Miklas, P.; McClean, P.E.; Osorno, J.M.; Runck, B.C.; White, J.W.; Kantar, M.B.; Ewing, P.M. Local to continental-scale variation in fitness and heritability in common bean. Crop Sci. 2022, 62, 767–779. [Google Scholar] [CrossRef]
  29. Bitocchi, E.; Bellucci, E.; Giardini, A.; Rau, D.; Rodriguez, M.; Biagetti, E.; Santilocchi, R.; Zeuli, P.S.; Gioia, T.; Logozzo, G.; et al. Molecular analysis of the parallel domestication of the common bean (Phaseolus vulgaris) in Mesoamerica and the Andes. New Phytol. 2013, 197, 300–313. [Google Scholar] [CrossRef]
  30. Rendón-Anaya, M.; Montero-Vargas, J.M.; Saburido-Álvarez, S.; Vlasova, A.; Capella-Gutierrez, S.; Ordaz-Ortiz, J.J.; Aguilar, O.M.; Vianello-Brondani, R.P.; Santalla, M.; Delaye, L.; et al. Genomic history of the origin and domestication of common bean unveils its closest sister species. Genome Biol. 2017, 18, 60. [Google Scholar] [CrossRef]
  31. Ron, A.M.; Gonzalez, A.M.; Rodino, A.P.; Santalla, M.; Godoy, L.; Papa, R. History of the Common Bean Crop: Its Evolution Beyond its Areas of Origin and Domestication. Arbor Cienc. Pensam. Cult. 2016, 192, a317. [Google Scholar] [CrossRef]
  32. Meza, N.; Rosas, J.C.; Martín, J.P.; Ortiz, J.M. Biodiversity of common bean (Phaseolus vulgaris L.) in Honduras, evidenced by morphological characterization. Genet. Resour. Crop Evol. 2013, 60, 1329–1336. [Google Scholar] [CrossRef]
  33. Palčić, I.; Karažija, T.; Petek, M.; Lazarević, B.; Ćustić, M.H.; Gunjača, J.; Liber, Z.; Carović-Stanko, K. Relationship between origin and nutrient content of Croatian common bean landraces. J. Cent. Eur. Agric. 2018, 19, 490–502. [Google Scholar] [CrossRef]
  34. Rodriguez, M.; Rau, D.; Bitocchi, E.; Bellucci, E.; Biagetti, E.; Carboni, A.; Gepts, P.; Nanni, L.; Papa, R.; Attene, G. Landscape genetics, adaptive diversity and population structure in Phaseolus vulgaris. New Phytol. 2016, 209, 1781–1794. [Google Scholar] [CrossRef]
  35. Sinkovič, L.; Pipan, B.; Sinkovič, E.; Meglič, V. Morphological Seed Characterization of Common (Phaseolus vulgaris L.) and Runner (Phaseolus coccineus L.) Bean Germplasm: A Slovenian Gene Bank Example. Biomed Res. Int. 2019, 2019, 6376948. [Google Scholar] [CrossRef]
  36. Wu, D.; Hought, J.; Baseggio, M.; Hart, J.P.; Gore, M.A.; Ilut, D.C. Genomic characterization of the Native Seeds/SEARCH common bean (Phaseolus vulgaris L.) collection and its seed coat patterns. Genet. Resour. Crop Evol. 2019, 66, 1469–1482. [Google Scholar] [CrossRef]
  37. Zeffa, D.M.; Nogueira, A.F.; Buratto, J.S.; Oliveira, R.B.; Neto, J.D.; Moda-Cirino, V. Genetic Variability of Mineral Content in Different Grain Structures of Bean Cultivars from Mesoamerican and Andean Gene Pools. Plants 2021, 10, 1246. [Google Scholar] [CrossRef]
  38. Blanch-Ramirez, J.; Calvet-Mir, L.; Aceituno-Mata, L.; Benyei, P. Climate change in the Catalan Pyrenees intersects with socioeconomic factors to shape crop diversity and management. Agron. Sustain. Dev. 2022, 42, 91. [Google Scholar] [CrossRef]
  39. Pedrali, D.; Proserpio, C.; Borgonovi, S.M.; Zuccolo, M.; Leoni, V.; Borgonovo, G.; Bernardi, A.M.; Scarafoni, A.; Pagliarini, E.; Giorgi, A.; et al. Nutritional Characterization and Novel Use of “Copafam” Bean (Phaseolus coccineus L.) for the Sustainable Development of Mountains Areas. Sustainability 2022, 14, 3409. [Google Scholar] [CrossRef]
  40. Loko, L.; Toffa, J.; Adjatin, A.; Akpo, A.J.; Orobiyi, A.; Dansi, A. Folk taxonomy and traditional uses of common bean (Phaseolus vulgaris L.) landraces by the sociolinguistic groups in the central region of the Republic of Benin. J. Ethnobiol. Ethnomed. 2018, 14, 52. [Google Scholar] [CrossRef]
  41. Soleri, D.; Cuevas, F.A.; Garcia, H.C.; Cleveland, D.A.; Smith, S.E. The Household Context of In Situ Conservation in a Center of Crop Diversity: Self-Reported Practices and Perceptions of Maize and Phaseolus Bean Farmers in Oaxaca, Mexico. Sustainability 2022, 14, 7148. [Google Scholar] [CrossRef]
  42. Brouwer, B.; Winkler, L.; Atterberry, K.; Jones, S.; Miles, C. Exploring the role of local heirloom germplasm in expanding western Washington dry bean production. Agroecol. Sustain. Food Syst. 2016, 40, 319–332. [Google Scholar] [CrossRef]
  43. Raggi, L.; Pacicco, L.C.; Caproni, L.; Álvarez-Muñiz, C.; Annamaa, K.; Barata, A.M.; Negri, V. Analysis of landrace cultivation in Europe: A means to support in situ conservation of crop diversity. Biol. Conserv. 2022, 267, 109460. [Google Scholar] [CrossRef]
  44. Hirata Sanches, V.; Adams, C.; Ferreira, F.F. An integrated model to study varietal diversity in traditional agroecosystems. PLoS ONE 2022, 17, e0263064. [Google Scholar] [CrossRef] [PubMed]
  45. Mina-Vargas, A.M.; McKeown, P.C.; Flanagan, N.S.; Debouck, D.G.; Kilian, A.; Hodkinson, T.R.; Spillane, C. Origin of year-long bean (Phaseolus dumosus Macfady, Fabaceae) from reticulated hybridization events between multiple Phaseolus species. Ann. Bot. 2016, 118, 957–969. [Google Scholar] [CrossRef]
  46. Contreras-Toledo, A.R.; Cortés-Cruz, M.; Costich, D.E.; Rico-Arce, M.d.L.; Magos Brehm, J.; Maxted, N. Diversity and conservation priorities of crop wild relatives in Mexico. Plant Genet. Resour. 2019, 17, 140–150. [Google Scholar] [CrossRef]
  47. FAO. The State of the World’s Biodiversity for Food and Agriculture; FAO Commission on Genetic Resources for Food and Agriculture Assessments: Rome, Italy, 2019; p. 572. [Google Scholar]
  48. Torres, L.J.E.; Montes, I.G.; Van Etten, J.; Aguero, F.R.; Sibelet, N. Why is there agrobiodiversity lost? case of the chagra inga in the colombian Amazon. Mundo Amaz. 2020, 11, 11–38. [Google Scholar] [CrossRef]
  49. Narendra, D.; Dahiya, G.S.; Janghel, D.K.; Dohare, S. Diversity Analysis for Seed Yield and its Component Traits among Faba Bean (Vicia faba L.) Germplasm Lines. Legume Res. 2022, 45, 689–694. [Google Scholar] [CrossRef]
  50. Perazzini, R.; Leonardi, D.; Ruggeri, S.; Alesiani, D.; D’Arcangelo, G.; Canini, A. Characterization of Phaseolus vulgaris L. Landraces Cultivated in Central Italy. Plant Foods Hum. Nutr. 2008, 63, 211–218. [Google Scholar] [CrossRef]
  51. Deaconu, A.; Mercille, G.; Batal, M. The Agroecological Farmer’s Pathways from Agriculture to Nutrition: A Practice-Based Case from Ecuador’s Highlands. Ecol. Food Nutr. 2019, 58, 142–165. [Google Scholar] [CrossRef]
  52. Intriago, R.; Gortaire Amézcua, R.; Bravo, E.; O’Connell, C. Agroecology in Ecuador: Historical processes, achievements, and challenges. Agroecol. Sustain. Food Syst. 2017, 41, 311–328. [Google Scholar] [CrossRef]
  53. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  54. Camacho, J. Good to Eat, Good to Think: Food, Culture and Biodiversity in Cotacachi. In Development with Identity: Community, Culture and Sustainability in the Andes; Robert, R., Ed.; CABI Publishing: Cambridge, MA, USA, 2006; pp. 156–170. [Google Scholar] [CrossRef]
  55. Zehetner, F.; Miller, W.P. Shaping an Andean Landscape: Processes Affecting Topography, Soils and Hydrology in Cotacachi. In Development with Identity: Community, Culture and Sustainability in the Andes; Robert, R., Ed.; CABI Publishing: Cambridge, MA, USA, 2006; pp. 49–58. [Google Scholar] [CrossRef]
  56. Gascon, J. Residential tourism and depeasantisation in the Ecuadorian Andes. J. Peasant. Stud. 2016, 43, 868–885. [Google Scholar] [CrossRef]
  57. Moates, A.S.; Campbell, B.C. Incursion, Fragmentation and Tradition: Historical Ecology of Andean Cotacachi. In Development with Identity: Community, Culture and Sustainability in the Andes; Robert, R., Ed.; CABI Publishing: Cambridge, MA, USA, 2006; pp. 59–82. [Google Scholar] [CrossRef]
  58. Skarbø, K. The Cooked is the Kept: Factors Shaping the Maintenance of Agro-biodiversity in the Andes. Hum. Ecol. 2014, 42, 711–726. [Google Scholar] [CrossRef]
  59. Melby, C.L.; Orozco, F.; Averett, J.; Munoz, F.; Romero, M.J.; Barahona, A. Agricultural Food Production Diversity and Dietary Diversity among Female Small Holder Farmers in a Region of the Ecuadorian Andes Experiencing Nutrition Transition. Nutrients 2020, 12, 2454. [Google Scholar] [CrossRef]
  60. Caulfield, M.E.; Hammond, J.; Fonte, S.J.; Florido, M.A.; Fuentes, W.; Meza, K.; Navarette, I.; Vanek, S.J.; van Wijk, M. Unpicking the Inter-relationships Between Off-Farm Livelihood Diversification, Household Characteristics, and Farm Management in the Rural Andes. Front. Sustain. Food Syst. 2021, 5, 724492. [Google Scholar] [CrossRef]
  61. FAO. Globally Important Agricultural Heritage Systems. Andean Chakra: An Ancestral Agricultural System of Kichwas Cotacachi. 2023. Available online: https://www.fao.org/giahs/giahsaroundtheworld/designated-sites/latin-america-and-the-caribbean/andean-chakra/en/ (accessed on 12 January 2023).
  62. Skarbø, K. Multiple trends in interspecific crop diversity: A longitudinal case study from the Ecuadorian Andes. Genet. Resour. Crop Evol. 2016, 63, 1319–1343. [Google Scholar] [CrossRef]
  63. Fernández, Á. Suministro de agua Segura a Comunidades Indígenas de Cotacachi—Ecuador. Master’s Thesis, Cooperación y Desarrollo Universidad de Castilla, La Mancha, Spain, October 2013. Available online: https://fundaciongeneraluclm.es/wp-content/uploads/2018/05/24-%C3%81NGEL-FERN%C3%81NDEZ.pdf (accessed on 8 December 2023).
  64. Chassagne, N.; Everingham, P. Buen Vivir: Degrowing extractivism and growing wellbeing through tourism. J. Sustain. Tour. 2019, 27, 1909–1925. [Google Scholar] [CrossRef]
  65. Hair, C. A View from Above: Alternative Perspectives on Smallholder Livehoods and Agrobiodiversity Conservation in Northern Ecuador. (Publication Number 1781). Ph.D. Thesis, University of Southern Mississippi, Hattiesburg, MS, USA, 2020. Available online: https://aquila.usm.edu/dissertations/1781/ (accessed on 23 September 2023).
  66. Skarbø, K. Reconfiguration of Andean Fields: Culture, Climate and Agrobiodiversity. Ph.D. Thesis, University of Georgia, Athenes, GA, USA, 2012. Available online: https://getd.libs.uga.edu/pdfs/skarbo_kristine_201205_phd.pdf (accessed on 19 December 2023).
  67. Calvo-Iglesias, M.S.; Crecente-Maseda, R.; Fra-Paleo, U. Exploring farmer’s knowledge as a source of information on past and present cultural landscapes: A case study from NW Spain. Landsc. Urban Plan. 2006, 78, 334–343. [Google Scholar] [CrossRef]
  68. Segnon, A.C.; Achigan-Dako, E.G.; Gaoue, O.G.; Ahanchédé, A. Farmer’s Knowledge and Perception of Diversified Farming Systems in Sub-Humid and Semi-Arid Areas in Benin. Sustainability 2015, 7, 6573–6592. [Google Scholar] [CrossRef]
  69. Jarvis, D.I.; Brown, A.H.D.; Cuong, P.H.; Collado-Panduro, L.; Latournerie-Moreno, L.; Gyawali, S.; Tanto, T.; Sawadogo, M.; Mar, I.; Sadiki, M.; et al. A global perspective of the richness and evenness of traditional crop-variety diversity maintained by farming communities. Proc. Natl. Acad. Sci. USA 2008, 105, 5326–5331. [Google Scholar] [CrossRef]
  70. Kouam, E.B.; Kamga-Fotso, A.M.-A.; Anoumaa, M. Exploring agro-morphological profiles of Phaseolus vulgaris germplasm shows manifest diversity and opportunities for genetic improvement. J. Agric. Food Res. 2023, 14, 100772. [Google Scholar] [CrossRef]
  71. Loko, L.; Orobiyi, A.; Adjatin, A.; Akpo, J.; Toffa, J.; Djedatin, G.; Dansi, A. Morphological characterization of common bean (Phaseolus vulgaris L.) landraces of Central region of Benin Republic. J. Plant Breed. Crop Sci. 2018, 10, 304–318. [Google Scholar] [CrossRef]
  72. Sadohara, R.; Long, Y.; Izquierdo, P.; Urrea, C.A.; Morris, D.; Cichy, K. Seed coat color genetics and genotype × environment effects in yellow beans via machine-learning and genome-wide association. Plant Genome 2022, 15, e20173. [Google Scholar] [CrossRef]
  73. Voysest, O. Variedades de Fríjol en América Latina y su Origen; Centro Internacional de Agricultura Tropical (CIAT): Cali, Colombia, 1983; 87p. [Google Scholar]
  74. International Board for Plant Genetic Resources. Phaseolus Vulgaris Descriptors. Rome, Italy. 1982, p. 15. Available online: https://cgspace.cgiar.org/server/api/core/bitstreams/f9a56aa4-2514-4aad-b961-cc04650ffcb5/content (accessed on 28 May 2019).
  75. Muñoz Arguedas, G.; Giraldo Avila, G.; Fernández de Soto, J. Descriptores Varietales: Arroz, Frijol, Maíz, Sorgo. 1993. Available online: https://hdl.handle.net/10568/54651 (accessed on 28 May 2019).
  76. Kanwar, R.; Mehta, D.K. Survey, collection and seed morphometric characterization of French bean (Phaseolus vulgaris L.) landraces of Himachal Pradesh. Legume Res. 2018, 41, 333–341. [Google Scholar] [CrossRef]
  77. ESRI. Environmental Systems Research Institute; Software ArcGIS v. 10.8; ESRI: Redlands, CA, USA, 2022. [Google Scholar]
  78. Carrijo, T.B.; da Silva, A.R. Modified Moran’s I for Small Samples. Geogr. Anal. 2017, 49, 451–467. [Google Scholar] [CrossRef]
  79. McAleece, N.; Gage, J.D.G.; Lambshead, P.J.D.; Paterson, G.L.J. BioDiversity Professional Statistics Analysis Software; Jointly Developed by the Scottish Association for Marine Science and the Natural History Museum: London, UK, 1997. [Google Scholar]
  80. Di Rienzo, J.; Casanoves, F.; Balzarini, M.; Gonzalez, L.; Tablada, E.; Robledo, C. InfoStat, Versión 2008; (Version v. 2008); Grupo InfoStat, FCA, Universidad Nacional de Córdoba: Cordoba, Argentina, 2008. [Google Scholar]
  81. Long, J.; Zhang, J.; Zhang, X.; Wu, J.; Chen, H.; Wang, P.; Wang, Q.; Du, C. Genetic Diversity of Common Bean (Phaseolus vulgaris L.) Germplasm Resources in Chongqing, Evidenced by Morphological Characterization. Front. Genet. 2020, 11, 697. [Google Scholar] [CrossRef]
  82. Alabi, T.; Adebola, P.; Asfaw, A.; De Koeyer, D.; Lopez, A.; Asiedu, R. Spatial Multivariate Cluster Analysis for Defining Target Population of Environments in West Africa for Yam Breeding. Int. J. Appl. Geospat. Res. 2019, 10, 1–30. [Google Scholar] [CrossRef]
  83. Raggi, L.; Tiranti, B.; Negri, V. Italian common bean landraces: Diversity and population structure. Genet. Resour. Crop Evol. 2013, 60, 1515–1530. [Google Scholar] [CrossRef]
  84. Soleri, D.; Worthington, M.; Aragón-Cuevas, F.; Smith, S.E.; Gepts, P. Farmers’ Varietal Identification in a Reference Sample of Local Phaseolus Species in the Sierra Juárez, Oaxaca, Mexico. Econ. Bot. 2013, 67, 283–298. [Google Scholar] [CrossRef]
  85. Siabato, W.; Guzmán, J. La autocorrelación espacial y el desarrollo de la geografía cuantitativa. Cuad. Geogr. Rev. Colomb. Geogr. 2019, 28, 1–22. [Google Scholar] [CrossRef]
  86. Bruno, A.; Katungi, E.; Stanley, N.T.; Clare, M.; Maxwell, M.G.; Paul, G.; Patrick, R.; Richard, E. Participatory farmers’ selection of common bean varieties (Phaseolus vulgaris L.) under different production constraints. Plant Breed. 2018, 137, 283–289. [Google Scholar] [CrossRef]
  87. Kelly, J.D.; Awale, H.E.; Wiersma, A.T.; Cichy, K.A.; Wright, E.M. Registration of ‘Yellowstone’ yellow bean. J. Plant Regist. 2021, 15, 265–270. [Google Scholar] [CrossRef]
  88. Xiong, M.; Zhao, M.; Lu, Z.-X.; Balasubramanian, P. Genotypic variation for phenolic compounds in developing and whole seeds, and storage conditions influence visual seed quality of yellow dry bean genotypes. Can. J. Plant Sci. 2019, 100, 284–295. [Google Scholar] [CrossRef]
  89. Hooper, S.D.; Bassett, A.; Sadohara, R.; Cichy, K.A. Elucidation of the low resistant starch phenotype in Phaseolus vulgaris exhibited in the yellow bean Cebo Cela. J. Food Sci. 2021, 86, 3975–3986. [Google Scholar] [CrossRef]
  90. Bento, J.A.C.; Bassinello, P.Z.; Carvalho, R.N.; Souza Neto, M.A.D.; Caliari, M.; Soares Júnior, M.S. Functional and pasting properties of colorful bean (Phaseolus vulgaris L) flours: Influence of the cooking method. J. Food Process. Preserv. 2021, 45, e15899. [Google Scholar] [CrossRef]
  91. Asfaw, A.; Almekinders, C.J.M.; Blair, M.W.; Struik, P.C. Participatory approach in common bean (Phaseolus vulgaris L.) breeding for drought tolerance for southern Ethiopia. Plant Breed. 2012, 131, 125–134. [Google Scholar] [CrossRef]
  92. Canci, H.; Yeken, M.Z.; Kantar, F.; Bozkurt, M.; Ciftci, V.; Ozer, G. Assessment of variation in seed morphological traits in Phaseolus sp. landraces from western Anatolia. Banat’s J. Biotechnol. 2019, 10, 75–88. [Google Scholar] [CrossRef]
  93. Aral, N.; Bakır, H. Spatiotemporal pattern of COVID-19 outbreak in Turkey. GeoJournal 2023, 88, 1305–1316. [Google Scholar] [CrossRef]
Figure 1. Map of the province of Imbabura showing the location of bean sampling sites (smallholder farms) in the study area (three administrative zones within Canton Cotacachi).
Figure 1. Map of the province of Imbabura showing the location of bean sampling sites (smallholder farms) in the study area (three administrative zones within Canton Cotacachi).
Agronomy 14 01666 g001
Figure 2. Bean landraces found in the study area. P. vulgaris landraces: 1. Allpa Poroto; 2. Bolon Morado; 3. Bolon Pintado; 4. Bolon Rojo; 5. Caca de Conejo; 6. Café Pintado; 7. Canario; 8. Cargabello; 9. Frejol Duro; 10. Frejol Tomate; 11. Josico; 12. Killu Allpa Poroto; 13. Killu Canario; 14. Matambre Negro; 15. Poroto Campeon; 16. Poroto Conejo; 17. Poroto Grande; 18. Puca Pintado; 19. Poroto Pintado; 20. Puca Poroto; 21. Puca Vaca Pintado; 22. Racu Pintado Poroto; 23. Rayado Poroto; 24. Sara Poroto; 25. Sucu Allpa Poroto; 26. Sucu Pintado; 27. Sucu Poroto; 28. Sucu Poroto Pintado; 29. Sucu Rayado; 30. Sucu Tomate Pintado; 31. Sucu Vaca Poroto; 32. Tomate Vaca Poroto; 33. Yana Pintado; 34. Yana Poroto; 35. Yana Sucu Poroto; 36. Yana Vaca Poroto; 37. Yura Allpa Poroto; 38. Yura Pintado; 39. Yura Pintado Poroto; 40. Yura Poroto and P. coccineus landraces: 41. Paco Bolon; 42. Paco Bolon Rayado; 43. Popayan Morado; 44. Popayan Morado Rayado; 45. Popayan Pintado; 46. Paco Popayan Pintado; 47. Yura Popayan.
Figure 2. Bean landraces found in the study area. P. vulgaris landraces: 1. Allpa Poroto; 2. Bolon Morado; 3. Bolon Pintado; 4. Bolon Rojo; 5. Caca de Conejo; 6. Café Pintado; 7. Canario; 8. Cargabello; 9. Frejol Duro; 10. Frejol Tomate; 11. Josico; 12. Killu Allpa Poroto; 13. Killu Canario; 14. Matambre Negro; 15. Poroto Campeon; 16. Poroto Conejo; 17. Poroto Grande; 18. Puca Pintado; 19. Poroto Pintado; 20. Puca Poroto; 21. Puca Vaca Pintado; 22. Racu Pintado Poroto; 23. Rayado Poroto; 24. Sara Poroto; 25. Sucu Allpa Poroto; 26. Sucu Pintado; 27. Sucu Poroto; 28. Sucu Poroto Pintado; 29. Sucu Rayado; 30. Sucu Tomate Pintado; 31. Sucu Vaca Poroto; 32. Tomate Vaca Poroto; 33. Yana Pintado; 34. Yana Poroto; 35. Yana Sucu Poroto; 36. Yana Vaca Poroto; 37. Yura Allpa Poroto; 38. Yura Pintado; 39. Yura Pintado Poroto; 40. Yura Poroto and P. coccineus landraces: 41. Paco Bolon; 42. Paco Bolon Rayado; 43. Popayan Morado; 44. Popayan Morado Rayado; 45. Popayan Pintado; 46. Paco Popayan Pintado; 47. Yura Popayan.
Agronomy 14 01666 g002
Figure 3. Cluster analysis based on bean seed morphological traits (primary color, secondary color, color around hilum, coat pattern, shape, seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio).
Figure 3. Cluster analysis based on bean seed morphological traits (primary color, secondary color, color around hilum, coat pattern, shape, seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio).
Agronomy 14 01666 g003
Figure 4. Principal component analysis based on the seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio.
Figure 4. Principal component analysis based on the seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio.
Agronomy 14 01666 g004
Figure 5. Geographic distribution of bean landraces in the three administrative zones of Canton Cotacachi, classified by the number of landraces in the smallholder farms.
Figure 5. Geographic distribution of bean landraces in the three administrative zones of Canton Cotacachi, classified by the number of landraces in the smallholder farms.
Agronomy 14 01666 g005
Figure 6. Spatial autocorrelation Moran’s I report for the number of bean landraces distributed in ten communities of the three administrative zones in Canton Cotacachi.
Figure 6. Spatial autocorrelation Moran’s I report for the number of bean landraces distributed in ten communities of the three administrative zones in Canton Cotacachi.
Agronomy 14 01666 g006
Table 1. Number of smallholder farmers selected for the study in each of the ten communities in three administrative zones of Canton Cotacachi.
Table 1. Number of smallholder farmers selected for the study in each of the ten communities in three administrative zones of Canton Cotacachi.
CommunityTotal Number of Smallholder FarmersSelected Smallholder Farmers
El Cercado14614
San Pedro889
Topo Grande10310
Iltaqui364
Chilcapamba727
Morochos12513
San Antonio del Punje404
Cumbas Conde11911
Arrayanes354
Morlán10811
Total87287
Table 2. Sections in the semi-structured interview conducted with 87 smallholder farmers of ten communities in three administrative zones of Canton Cotacachi.
Table 2. Sections in the semi-structured interview conducted with 87 smallholder farmers of ten communities in three administrative zones of Canton Cotacachi.
SectionVariableCategories/Description
1. General informationFarmer’s name
Age
Gendermale, female
Ethnic groupIndigenous, mestizo
LanguageKichwa, Spanish, Kichwa and Spanish
Education levelno formal education, primary (elementary), secondary, higher education
Access to irrigation wateryes, no
Land ownershipown, leased, community, borrowed
Geographic coordinateslatitude, longitude, altitude
2. Bean diversity Identification of landracesCould you name the bean landraces do you grow in your chakra?
Can we get samples of the varieties from your chakra?
Can you identify the name of the landraces you have?
Seed’s sourceInherited, market, public institution, seed exchange
Production’s purposeSelf-consumption, market, self-consumption and market, seed production
3. Perception of riskLandraces at risk of disappearing What landraces do you think that are disappearing?
What do you think that are the main reason?
Table 3. Bean landrace names given by smallholder farmers and their respective translations from Kichwa and name meaning.
Table 3. Bean landrace names given by smallholder farmers and their respective translations from Kichwa and name meaning.
Landrace Translation, Name MeaningPhoto 1AF 2Landrace Translation, Name MeaningPhoto 1AF 2
(Farmers’ Name) (Farmers’ Name)
Canario 3 “Canary” (seed color)731 Puca Poroto “Red bean” (seed color)206
Poroto Grande “Big bean” (seed size)1728 Yura Pintado “White spotted” (seed color)386
Popayan Morado “Purple bean” (seed color)4314 Bolon Pintado “Round mottled”
(seed shape, color)
35
Josico Unknown meaning1113 Bolon Rojo “Red round” (seed color, shape)45
Yana Poroto “Black bean” (seed color)3413 Paco Bolon Rayado “Buff, round, striped” (seed color, shape, pattern)425
Puca Pintado “Red mottled” (seed color)1812 Poroto Conejo “Rabbit bean” (seed appearance)165
Sara Poroto “Corn” (crop association)2412 Racu Pintado Poroto “Thick speckled bean”
(seed size, pattern)
225
Yura Allpa Poroto “Ground white bean” (habit, seed color)3712 Rayado Poroto “Stripped bean” (seed pattern)235
Allpa Poroto “Ground bean” (habit)111 Yana Pintado “Black mottled” (seed color, pattern)335
Killu Canario “Yellow canary” (seed color)1311 Sucu Allpa Poroto “Ashen ground bean” (seed color, habit)254
Cargabello Juxtaposed word meaning: high yield810 Sucu Pintado “Ashen spotted” (seed color, pattern)264
Frejol Tomate “Orange bean” (seed color)1010 Sucu Rayado “Ashen striped” (seed color, pattern)294
Popayan Morado Rayado “Purple stripped bean” (seed color, pattern)4410 Yura Popayan “White bean” (seed color)474
Sucu Poroto “Ashen bean” (seed color)2710 Caca de Conejo “Rabbit excrement” (seed appearance)53
Yana Vaca Poroto “Black spottled bean” (seed color, pattern)3610 Frejol Duro “Hard bean” (seed consistency)93
Popayan Pintado “Speckled bean” (seed pattern)459 Tomate Vaca Poroto “Orange cow bean” (seed color, pattern)323
Sucu Poroto Pintado “Ashen spotted bean” (seed color, pattern)229 Yura Pintado Poroto “White spotted bean” (seed color, pattern)393
Café Pintado “Brown speckled” (seed color, pattern)68 Bolon Morado “Purple round” (seed color, shape)22
Matambre Negro “Black kill hunger” (seed color, juxtaposed word meaning: kill hunger)148 Killu Allpa Poroto “Yellow ground bean” (seed color, habit)122
Yana Sucu Poroto “Orange black bean” (seed color)358 Puca Vaca Pintado “Red cow spotted” (seed color, pattern)212
Paco Bolon “Light yellow brown, round” (seed color, shape)417 Sucu Vaca Poroto “Ashen cow bean” (seed color, pattern)312
Poroto Pintado “Mottled bean” (seed pattern)197 Poroto Campeon “Champion bean” (yield)151
Yura Poroto “White bean” (seed color)407 Sucu Tomate Pintado Red orange speckled (seed color, pattern)301
Paco Popayan Pintado “White speckled bean” (seed color, pattern)466
1 See Figure 3; 2 AF = absolute frequency; 3 Canario is not considered a landrace; instead, it is an improved variety.
Table 4. Eigenvalues and the contribution of each principal component based on the seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio.
Table 4. Eigenvalues and the contribution of each principal component based on the seed length (L), seed width (W), seed thickness (T), L/W ratio, and W/T ratio.
ComponentEigenvalue Contribution Rate (%)Accumulative Contribution Rate (%)
PC12.5351%51%
PC22.1743%94%
PC30.296%100%
Table 5. ANOVA for the five bean seed quantitative traits.
Table 5. ANOVA for the five bean seed quantitative traits.
Traitp-Value
Seed Lenght (L)0.3473
Seed Width (W)<0.0001
Seed Thickness (T)<0.0001
L/W<0.0001
W/T0.3439
Table 6. Shannon–Wiener Diversity–Equity (H′) index and Shannon J’ homogeneity index for P. vulgaris and P. coccineus landraces at the community level.
Table 6. Shannon–Wiener Diversity–Equity (H′) index and Shannon J’ homogeneity index for P. vulgaris and P. coccineus landraces at the community level.
IndexArrayanesChilcapambaCumbas CondeEl CercadoItalquíMorlánMorochosSan Antonio PungeSan PedroTopo Grande
Shannon H′ Log Base 101.3220.9871.261.4550.9411.3151.5490.4521.2461.002
Shannon J’0.9710.9870.9530.9430.9410.9650.960.9460.9580.962
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Romero-Astudillo, M.J.; Tapia, C.; Giménez de Azcárate, J.; Montalvo, D. Diversity of Common Bean (Phaseolus vulgaris L.) and Runner Bean (Phaseolus coccineus L.) Landraces in Rural Communities in the Andes Highlands of Cotacachi—Ecuador. Agronomy 2024, 14, 1666. https://doi.org/10.3390/agronomy14081666

AMA Style

Romero-Astudillo MJ, Tapia C, Giménez de Azcárate J, Montalvo D. Diversity of Common Bean (Phaseolus vulgaris L.) and Runner Bean (Phaseolus coccineus L.) Landraces in Rural Communities in the Andes Highlands of Cotacachi—Ecuador. Agronomy. 2024; 14(8):1666. https://doi.org/10.3390/agronomy14081666

Chicago/Turabian Style

Romero-Astudillo, María José, César Tapia, Joaquín Giménez de Azcárate, and Diego Montalvo. 2024. "Diversity of Common Bean (Phaseolus vulgaris L.) and Runner Bean (Phaseolus coccineus L.) Landraces in Rural Communities in the Andes Highlands of Cotacachi—Ecuador" Agronomy 14, no. 8: 1666. https://doi.org/10.3390/agronomy14081666

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop