Next Article in Journal
Hydrogen Peroxide Is Involved in Methane-Alleviated Cadmium Toxicity in Alfalfa (Medicago sativa L.) Seedlings by Enhancing Cadmium Chelation onto Root Cell Walls
Previous Article in Journal
Metabolome Shift in Centella asiatica Leaves Induced by the Novel Plant Growth-Promoting Rhizobacterium, Priestia megaterium HyangYak-01
Previous Article in Special Issue
An Extended Application of the Fast Multi-Locus Ridge Regression Algorithm in Genome-Wide Association Studies of Categorical Phenotypes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions

by
Stephanie Mariel Alves
1,
Giselly Figueiredo Lacanallo
2,
Maria Celeste Gonçalves-Vidigal
1,*,
Mariana Vaz Bisneta
1,
Andressa Gonçalves Vidigal Rosenberg
1 and
Pedro Soares Vidigal Filho
1
1
Pós-Graduação em Genética e Melhoramento, Universidade Estadual de Maringá, Av. Colombo, 5790, Maringá 87020-900, Brazil
2
Departamento de Agronomia, Universidade Estadual de Maringá, Av. Colombo, 5790, Maringá 87020-900, Brazil
*
Author to whom correspondence should be addressed.
Plants 2024, 13(18), 2638; https://doi.org/10.3390/plants13182638
Submission received: 17 August 2024 / Revised: 4 September 2024 / Accepted: 18 September 2024 / Published: 21 September 2024

Abstract

:
Exploring genetic resources through genomic analyses has emerged as a powerful strategy to develop common bean (Phaseolus vulgaris L.) cultivars that are both productive and well-adapted to various environments. This study aimed to identify genomic regions linked to morpho-agronomic traits in Mesoamerican and Andean common bean accessions and to elucidate the proteins potentially involved in these traits. We evaluated 109 common bean accessions over three agricultural years, focusing on traits including the grain yield (YDSD), 100-seed weight (SW), number of seeds per pod (SDPD), number of pods per plant (PDPL), first pod insertion height (FPIH), plant height (PLHT), days to flowering (DF), and days to maturity (DPM). Using multilocus methods such as mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, we identified 36 significant SNPs across all chromosomes (Pv01 to Pv11). Validating these SNPs and candidate genes in segregating populations is crucial for developing more productive common bean cultivars through marker-assisted selection.

1. Introduction

The common bean (Phaseolus vulgaris L.) holds a crucial role in human nutrition, both in Brazil and globally [1]. This significance stems from its rich nutritional profile, which includes proteins, micronutrients, and essential minerals [2]. Furthermore, the common bean is widely cultivated across diverse regions worldwide, showcasing its adaptability to a variety of climatic conditions and soil types [3,4].
The genetic variability among common bean cultivars in Brazil serves as an invaluable resource for the development of more productive varieties that are suited to various environmental conditions [5]. Predominantly, the Brazilian germplasm is composed of landrace varieties from the Mesoamerican gene pool, which provides greater genetic diversity compared with the Andean gene pool [6,7,8,9]. The domestication of the common bean fostered genetic diversification, driving physiological and morphological changes. These adaptations, including variations in growth habit, vegetative cycle, photoperiod response, seed shape, color, and size, enabled cultivated plants to thrive in different climates, distinguishing them from their wild progenitors [10,11,12].
Increasing productivity is a central objective in plant breeding programs. Achieving this requires a thorough understanding of the genetic factors that influence productivity and its components. In common bean cultivation, productivity is tied to a variety of morphological, agronomic, and physiological traits. Among these, plant architecture, as well as flowering and maturation periods, play significant roles in the plant’s adaptability and overall productivity [13,14,15].
However, enhancing these traits through traditional breeding methods demands extensive fieldwork, involving evaluations across multiple environments and years. This process is not only costly, but also time-intensive. In this context, the use of molecular markers has emerged as a valuable tool. Molecular markers facilitate the identification of specific genomic regions associated with productivity traits. These regions can be targeted through marker-assisted selection, thereby streamlining the breeding process and making it more efficient and cost-effective [13,16].
Next Generation Sequencing technologies has catalyzed advancements in genomic research, leading to the development of more detailed and dense genetic maps with numerous markers [17]. This genetic information is increasingly being integrated into common bean breeding programs, enabling the application of advanced methodologies such as Genomic Selection (GS), Quantitative Trait Loci (QTL) mapping, and Genome-Wide Association Studies (GWAS) in breeding for complex traits [13,14,15,16,17].
Genetic mapping can be performed using two primary approaches: linkage mapping, which relies on genetic linkage, and association mapping, which is based on linkage disequilibrium (LD). The LD-based approach, commonly known as genome-wide association studies (GWAS), is a powerful tool for identifying genes and QTLs. In GWAS, phenotypic data for a trait of interest are collected from a large group of individuals, along with their genetic information. These data are then analyzed to pinpoint regions of the genome associated with the observed phenotypes. These regions may harbor genes or regulatory elements that are crucial to trait expression, indicating their involvement in the manifestation of the studied trait [18,19].
Genome-wide association studies (GWAS) have proven to be an efficient method for the genetic mapping of quantitative traits and are widely used in the analysis of agronomic characteristics [16,20,21,22,23,24,25]. Furthermore, other studies have focused on identifying genomic associations related to traits such as grain yield and plant phenology in common beans [26,27,28,29]. Most research on genetic variability has focused on germplasm banks, which are vital repositories of genetic diversity. The Germplasm Bank of the Núcleo de Pesquisa Aplicada à Agricultura (Nupagri) at Universidade Estadual de Maringá (UEM) holds traditional cultivars (accessions) from various regions of the country and abroad, alongside elite lines developed for research.
Genome association studies of accessions adapted to specific environments represent an effective approach to investigating quantitative traits. Additionally, multilocus methods have been relatively underexplored in common bean breeding. Thus, this study aims to identify genomic regions associated with morphoagronomic traits in Mesoamerican and Andean common bean accessions from the Nupagri Germplasm Bank.

2. Results

2.1. Phenotypic Data

The average temperatures were suitable for normal plant development during all experimental periods; the average daily temperature in 2019 was 25.39 °C; in 2020, it was 24.32 °C; and in 2021, it was 23.73 °C. However, the precipitation was insufficient for the optimal development of common bean plants. In 2019, the accumulated rainfall observed was only 346.8 mm, 315.4 mm in 2020, and 470.2 mm in 2021. Thus, in all evaluated years irrigation water management was implemented. In fact, in 2019 and 2021 rainfall was observed during harvesting leading to losses in bean yield (Figure 1).
All traits evaluated in the experiments conducted in 2019, 2020, and 2021 are available in the supplementary file (Table S1). Statistically significant differences and genotype by environment (GxE) interactions were significant across all traits analyzed (Table 1). The average grain yields of the common bean accessions were 990.44 kg ha⁻1 in 2019, 1736.26 kg ha⁻1 in 2020, and 1093.57 kg ha⁻1 in 2021. Notably, the grain yield was highest in 2020, indicating particularly favorable growing conditions during that year, likely due to rainfall at the beginning of the vegetative phase (Figure 1).
The 100-seed weight of the accessions averaged 28.00 g, 27.92 g, and 28.10 g for the years 2019, 2020, and 2021, respectively. Although these variations are minor, they can significantly impact seed quality and commercial value. For the number of seeds per pod and the number of pods per plant, the averages were 3.21, 3.60, and 2.73, and 10.04, 9.67, and 10.42, respectively, across the same years (Table 1). These variations suggest that the plants are adapting to the specific conditions of each environment, presenting opportunities for productivity improvement.
Both the first pod insertion height and plant height showed variability, indicating adaptations in plant architecture to varying environmental conditions. Notably, higher averages were observed in 2020, likely due to the more favorable rainfall conditions. Similarly, the variations in days to flowering and maturity reflected differing responses to the growth cycle among accessions across the years. Although the averages between years showed minimal differences, the significant interaction between genotype and years suggested that certain genotypes may have responded more strongly to specific conditions in a particular year (Table 1). These findings underscore the importance of considering the growth cycle when selecting for more adapted varieties. Overall, the results underscore the phenotypic diversity among the accessions while clearly demonstrating the significant influence of environmental conditions on the traits evaluated. These findings emphasize the critical importance of ongoing research to improve crop performance across diverse environments.

2.2. Genome-Wide Association Study (GWAS)

The five ML-GWAS methods identified a total of 36 SNPs (single-nucleotide polymorphisms) associated with the eight evaluated morphoagronomic traits—grain yield (YDSD), 100-seed weight (SW), number of seeds per pod (SDPD), number of pods per plant (PDPL), first pod insertion height (FPIH), plant height (PLHT), number of days to flowering (DF), and number of days to maturity (DPM)—across all chromosomes: Pv01 through Pv11 (Table 2). The largest number of significant SNPs was found on chromosome Pv04, with a total of six SNPs. Chromosomes Pv01, Pv06, and Pv07 each harbored five significant SNPs (Table 2).

2.2.1. SNPs for Yield and Its Components

A total of 15 significant SNPs were identified for yield and its components, including YDSD, SW, SDPD, and PDPL. Among these, five SNPs associated with YDSD were found on chromosomes Pv02, Pv04, Pv05, Pv07, and Pv08 across the three years analyzed (Table 2). Notably, SNP ss715647997 was observed in the combined analysis, as well as in the year 2019 (Figure 2a). This SNP is located at 48,720,860 bp on chromosome Pv02 within the gene model Phvul.002G328800 that encodes a serine–threonine protein kinase. It explained 13.38 to 18.32% of the phenotypic variation.
For 100-seed weight (SW), SNPs ss715639786, ss715645689, and ss715645602 were detected on chromosomes Pv06, Pv07, and Pv09, respectively, showing significant associations in the combined analysis and accounting for 4.83% to 57.29% of the phenotypic variation (Figure 2b). SNP ss715645689 was also identified in 2019 and 2020, while SNP ss715645602 was noted in 2021. These SNPs were mapped to genomic positions 22,116,966 bp, 819,750 bp, and 33,444,042 bp, respectively, in the reference genome of common Andean beans [11,30], wherein genes Phvul.006G105000, Phvul.007G011800, and Phvul.009G226000 were annotated, encoding for ankyrin repeat-containing protein, methionyl-tRNA synthetase, and aspartyl protease family protein, respectively.
The SDPD trait was significantly associated in 2020 with SNP ss715645673 on Pv06, explaining 4.92% to 8.24% of the phenotypic variation. This SNP is located at position 27,760,992 bp and is within the genomic region containing the gene model Phvul.006G166700, which encodes an abc transporter c family member 11-related protein.
For the PDPL trait, four significant SNPs were identified on chromosomes Pv01, Pv02, Pv10, and Pv11 across all three years analyzed (Table 2). In the combined analysis, SNPs ss715645301 and ss715647997 on chromosomes Pv01 and Pv02, respectively, explained between 8.61% and 43.10% of the phenotypic variation. In the reference genome, close to the ss715645301 marker, the model gene Phvul.001G264600 is found, which encodes a C2H2-type zinc finger protein. Interestingly, the other SNP ss715647997 is the same marker found to be associated with YDSD, therefore showing a pleiotropic effect.

2.2.2. SNPs for Plant Architecture

Twelve significant SNPs were identified for plant architecture traits, including the first pod insertion height (FPIH) and plant height (PLHT). Among these, six SNPs associated with FPIH were detected on chromosomes Pv03, Pv04, Pv05, Pv06, and Pv09 across the three years analyzed (Table 2), three of them were found in combined analysis. Notably, two SNPs on chromosome Pv04, ss715646910 and ss715649971, were particularly significant. SNP ss715646910 is located in the genomic region containing the gene model Phvul.004G014500, which encodes a protein farnesyltransferase subunit beta (FNTB), and explains 9.25% to 10.92% of the phenotypic variation. SNP ss715649971 is located within the model gene Phvul.004G016000, which encodes leucine-rich repeat-containing protein, and accounts for 7.02% to 12.51% of the phenotypic variation (Table 2). It is noteworthy that SNP ss715649971 was also observed in the years 2019 and 2021 (Figure 2d). On chromosome Pv06, SNP ss715645752 was found to be associated with FPIH, explaining 6.28% to 21.51% of the phenotypic variation. This SNP is located at genomic position 25,791,849 bp, close to the gene model Phvul.006G144200, which encodes serine-threonine protein kinase (Table 2; Figure 2d). This SNP was also observed in 2019 and 2020.
For the plant height (PLHT) trait, six significant SNPs were found to be associated with this trait on chromosomes Pv01, Pv02, Pv04, and Pv07, wherein four were consistently found in combined analysis. On chromosome Pv01, two notable SNPs were identified: ss715647368 and ss715645852. SNP ss715647368 explained 18.12% to 38.65% of the phenotypic variation and was observed in all three years analyzed. It showed pleiotropy, affecting plant height (PLHT), days to flowering (DF), and days to maturity (DPM) (Table 2; Figure 2e). SNP ss715647368 is located at the position 45,894,030 bp close to the gene model Phvul.001G193200 that encodes a protease s28 pro-x carboxypeptidase-related protein. SNP ss715645852, also on chromosome Pv01, explained 7.10% to 10.16% of the phenotypic variation and was also identified in 2019. This SNP is located at position 49,657,760 bp close to the gene model Phvul.001G236100 that encodes a poly(a) RNA polymerase gld2 protein.

2.2.3. SNPs for Plant Phenology

A total of eleven SNPs significantly associated with plant phenology traits, specifically days to flowering (DF) and days to maturity (DPM), were identified across chromosomes Pv01, Pv03, Pv04, Pv06, Pv07, and Pv10 (Table 2). For DF, four SNPs were significant in the combined analysis (Figure 3a). On chromosome Pv01, two notable SNPs were found: ss715647368 (previously described with pleiotropy effect) and ss715639536 at position 46,027,600 bp. SNP ss715647368 explained 12.71% to 27.53% of the phenotypic variation, and SNP ss715639536 accounted for 24.72% to 24.97%. SNP ss715647368 was observed in 2019 and 2020, while ss715639536 was identified in all years analyzed. SNP ss715639536 is positioned near the gene model Phvul.001G194100, encoding peroxiredoxin (alkyl hydroperoxide reductase subunit C).
On chromosomes Pv07 and Pv10, SNPs ss715646355 and ss715640116 were identified, explaining phenotypic variations ranging from 1.48% to 5.72% and 14.42% to 56.84%, respectively. Both SNPs were detected across all analyzed years, as well as in the combined analysis (Table 2). SNP ss715646355 is located at the beginning of chromosome 3,770,008 bp in the reference genome of Andean common beans (v1.1) close to the model gene Phvul.007G046900 that encodes 16S rRNA (uracil(1498)-N(3))-methyltransferase/M(3)U(1498)-specific methyltransferase protein. SNP ss715640116 on Pv10 is positioned at 7,867,881 bp close to Phvul.010G049800 that encodes a phospholipase A1 (DAD1).
For the trait days to maturity (DPM), nine significant SNPs were identified (Figure 3b); four of them were significant in the combined analysis. The SNP ss715639271 on Pv01 was significantly associated with DPM in the combined analyses and the year 2019, explaining 21.53% to 45.08% of the phenotypic variation. This SNP is located at the end of the chromosome at the position 45,746,595 close to the model gene Phvul.001G191600, which encodes beta-fructofuranosidase (E3.2.1.26, sacA). The SNP ss715647636 on chromosome Pv03 was observed in the combined analysis and the year 2019, accounting for 2.44% to 12.81% of the phenotypic variation. This SNP located at 3,963,582 close to Phvul.003G038400 encodes act domain-containing protein 5.
On chromosome Pv06, two SNPs associated with the days to maturity (DPM) trait were identified in the combined analysis: ss715648492 and ss715645677. SNP ss715648492, located at position 25,619,371 bp, is close to the gene model Phvul.006G142100, which encodes armadillo/beta-catenin-like repeat-containing protein. This SNP explained between 5.51% and 16.31% of the phenotypic variation and was detected in the overall analysis, 2019, and 2021 (Figure 3b). SNP ss715645677, located near the gene model Phvul.006G162500, which encodes PPR repeat (PPR)//PPR repeat family (PPR_2)//Pentatricopeptide repeat domain (PPR_3), explained phenotypic variation ranging from 3.70% to 5.22%. This SNP was observed in both the combined analysis and specifically in 2019.

3. Discussion

3.1. Genotype × Environment Interaction

Germplasm banks play a crucial role in conserving and exploring the genetic variability in different bean accessions. The 109 accessions from the Nupagri Germplasm Bank encompass a diverse range of grain types, including Carioca, Black, Jalo, Red, Rosinha, Roxinho, White, and Bolinha. This genetic diversity represents a valuable reservoir of agronomic traits, which are vital for developing more productive cultivars that are resistant to diseases and adapted to various environmental conditions.
Genotype × environment (G × E) interaction studies are crucial for understanding the phenotypic and genotypic variation in common bean accessions of Mesoamerican and Andean origins. These studies provide insights into how genetic diversity expresses itself under different environmental conditions and help in identifying genotypes that are better adapted to specific environments. Observations of G × E interactions revealed that all evaluated traits exhibited unique behaviors, emphasizing that such interactions can produce varying phenotypic expressions, even for the same genotypes across different environments and years. The presence of G × E interactions is commonly observed in linkage mapping and association mapping studies, where it influences the interactions between QTL × environment and/or SNP × environment [15,27,31].

3.2. Common Bean Genomic Regions Associated with Morphological and Agronomic Traits

GWAS has proven to be a valuable tool for exploring genetic diversity, particularly in common bean cultivation. Studies focusing on accessions adapted to specific environments have demonstrated the effectiveness of this approach in investigating quantitative traits, as genomic regions significantly associated with traits across multiple environmental conditions have been identified. Additionally, markers specific to particular environmental conditions have also been discovered [15,27]. Multi-locus methods, which have been relatively underutilized in common bean research, are now gaining prominence for their ability to offer greater statistical power and precision in detecting SNP markers in GWAS [16,25]. This approach is particularly valuable for complex traits that are strongly influenced by environmental factors [32,33]. Multi-locus analyses of accessions from the Nupagri Germplasm Bank have led to significant advances in understanding the genetic control of the evaluated morpho-agronomic traits.
In this study, significant variations were observed among 109 accessions from the Nupagri Germplasm Bank. GWAS revealed genomic regions associated with SNPs for morphoagronomic traits across all chromosomes: Pv01 through Pv11. Similarly, Delfini et al. [15] identified markers related to productivity and its components across all chromosomes in their study of 178 Mesoamerican common bean accessions using ML-GWAS. Different genomic regions were found associated with grain yield (YDSD) and seed weight (SW) across different years. This variability can be attributed to the quantitative nature of these traits, which makes them highly responsive to environmental conditions. Depending on climatic and precipitation patterns, specific regions on the bean chromosomes exhibited varying impacts on these traits. For plant height (PLHT), days to flowering (DF), and days to maturity (DPM), distinct associations were also observed for each trait across years. Nevertheless, the Pv01 regions consistently showed a major effect on controlling these traits throughout the years. In the combined analysis of years; most associations mirrored those identified in one or more individual years, depending on the trait. This approach enriches the findings by revealing consistent patterns and strengthening the overall associations observed.
The GWAS identified significant regions on chromosomes Pv06, Pv07, and Pv09 associated with 100-seed weight (SW). While initial studies mapped 100-seed weight to chromosome Pv07 [34], further research has shown it to be a polygenic quantitative trait distributed across all chromosomes [11,27,34,35,36,37,38]. Many of the identified alleles are linked to domestication processes and ecogeographic factors [11,27,35,36], which aligns with our findings of SW being a quantitative and complex trait among Mesoamerican and Andean accessions.
In the functional annotation of the seed weight (SW) trait, a notable protein identified was an aspartyl protease family protein with similarity to cellulose synthase, which belongs to the CESA gene superfamily. This superfamily, found in Arabidopsis thaliana and other seed plants, encodes enzymes essential for cellulose production. The CESA genes are divided into groups based on their function and structure, such as CESA1, CESA2, and CESA3, found in dicotyledonous plants [39,40].
For productivity traits such as grain yield (YDSD) and number of pods per plant (PDPL), a pleiotropic SNP was identified on chromosome Pv02. Similar pleiotropic effects were observed in studies by Delfini et al. [15] and Nkhata et al. [41], who reported comparable findings on chromosome Pv11. These pleiotropic effects suggest that these markers may influence multiple traits in common beans. The genomic region containing this SNP is functionally annotated to a serine/threonine protein kinase, a key player in various cellular processes, including cell cycle regulation, cell signaling, and stress response [42].
For traits related to plant architecture and phenology, including the first pod insertion height (FPIH), plant height (PLHT), days to flowering (DF), and days to maturity (DPM), SNPs were identified on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, Pv07, Pv09, and Pv10. These traits were also mapped to different chromosomes in Mesoamerican and Andean accessions, which exhibit various growth habits, particularly on chromosomes Pv01, Pv04, Pv06, Pv08, Pv09, Pv10, and Pv11 [14,27,28,43]. The results of this study further evidence the importance of chromosomes Pv01 and Pv10 in controlling bean architecture, as previously described in the literature. Notably, a SNP named ss715647368, located at position 45,894,030 bp on chromosome Pv01, was identified as pleiotropic, indicating a possible relationship between plant architecture and phenology traits (Table 2, Figure 2 and Figure 3). Furthermore, SNP ss715647368 was functionally annotated to a PROTEASE S28 PRO-X carboxypeptidase on chromosome Pv01, which is homologous to the S28 serine carboxypeptidase family protein (AT5G65760) in Arabidopsis thaliana. These proteins play crucial roles in various biological processes, including growth, development, and the stress response. S28 family carboxypeptidases are particularly vital for regulating protein and peptide maturation, modulating hormonal activity, and enabling plants to adapt to challenging conditions, such as nutrient availability fluctuations and environmental stress [44].

4. Materials and Methods

4.1. Plant Material

This study utilized a total of 109 common bean accessions (55 Mesoamerican and 54 Andean) from the Nupagri Germplasm Bank at the Núcleo de Pesquisa Aplicada à Agricultura, Universidade Estadual de Maringá, Paraná, Brazil. These accessions were collected from farmers in the Brazilian states of Goiás, Mato Grosso do Sul, Minas Gerais, Paraná, and Santa Catarina (Figure 4). Descriptions of the evaluated accessions, such as name, origin, gene pool, and market class, are provided in Table S2.

4.2. Phenotyping Evaluation

The experiment was conducted at the Centro de Treinamento em Irrigação (CTI) of the Universidade Estadual de Maringá (UEM) in Maringá, Paraná, Brazil (latitude 23°25′ S and longitude 51°57′ W), during the period from July to November in the agricultural years of 2019, 2020, and 2021. A simple 11 × 11 lattice design was used for the evaluations, covering a total of 109 accessions (Figure 4), along with 12 additional accessions—BRS Pérola, BAT 447, IPR Uirapuru, SEA 5, BRS Esteio, BGF 04, BGF 09, BGF 59, BGF 92, BGF 155, BGF 203, and BGF 204—used to optimize the experimental design.
The experiment included three replications, with plots consisting of 2 m rows spaced 0.5 m apart. The planting density was 12 plants per linear meter, with a border of 4 common bean rows. The plants were fertilized and irrigated according to the specific needs of the crop.
The following morpho-agronomic traits were measured: (a) days to flowering (DF)—the number of days from emergence until 50% of the plants in the plot showed at least one flower; (b) days to physiological maturity (DPM)—the number of days from emergence until 90% of the pods lost their green coloration and began to dry. In the post-harvest phase, the following traits were evaluated: (c) first pod insertion height (FPIH)—the average height, in centimeters, from the soil surface to the first pod insertion; (d) final plant height (PLHT)—the average height, in centimeters, measured from the soil to the tip of the uppermost leaves; (e) number of pods per plant (PDPL)—the total number of pods in the plot divided by the number of plants in the plot; (f) seeds per pod (SDPD)—the average number of seeds in 10 randomly selected pods from each plot; (g) 100-seed weight (SW)—the weight, in grams, of 100 seeds sampled from each plot; and (h) grain yield (YDSD)—measured per plot at 13% moisture and expressed in kg ha−1. The crop ontology for agricultural data for the common bean was used as a reference to unify the name of traits (https://cropontology.org/term/CO_335:ROOT (accessed on 4 December 2023)).

4.3. Phenotypic Data Analysis

The phenotypic data from each experiment were analyzed using SELEGEN REML/BLUP software (https://www.conectagem.com/softwares (accesses on 7 February 2022)) [45]. Genotypic values were predicted through the best linear unbiased predictor (BLUP) approach, with variance components for random factors estimated via the restricted maximum likelihood (REML) method, following model 52:
y = X r + Z g + W b + T i + e
where y = data vector, r = vector of repetition effects (assumed as fixed) added to the general average, g = vector of genotypic effects (assumed as random), b = vector of block effects (assumed as random), i = vector of genotype × environment interaction effects (random), and e = vector of errors or residues (random). The uppercase letters represent the incidence matrices for the aforementioned effects.
The significance of the model effects was assessed using the Likelihood Ratio Test (LRT), which provided the basis for deviance analysis (ANADEV). Model fit was evaluated by comparing analyses with and without the effects of genotype and genotype × environment interaction. The deviance of the model excluding the specific effect was subtracted from the deviance of the full model, and the resulting value was compared against the chi-square (χ²) distribution with 1 degree of freedom at the 1% and 5% probability levels.
Mathematically, this process is represented as follows:
L R T = 2 l n M L   o f   r e d u c e d   m o d e l M L V   o f   c o m p l e t e   m o d e l
where ln denotes the natural logarithm, while ML stands for maximum likelihood estimation. The analyses evaluated the effect of each line on the phenotype, considering fixed effects and random effects to obtain the best linear unbiased estimators (BLUEs) and best linear unbiased predictors (BLUPs), respectively.

4.4. Genotyping of the Common Bean Accessions Using BARCBean6K_3 BeadChip

Genomic DNA was extracted from trifoliolate leaves of each accession using the DNeasy Plant Mini Kit (Qiagen, CA, USA) according to the manufacturer’s instructions. The DNA concentration was measured with a Qubit fluorometer, following the manufacturer’s guidelines. The DNA samples were genotyped using 5398 SNP markers on the BARCBean6K_3 Illumina BeadChip, adhering to the Infinium HD Assay Ultra Protocol (Illumina Inc., San Diego, CA, USA). The BeadChip was imaged with the Illumina BeadArray Reader to assess fluorescence intensity [46]. Automatic allele calling for each locus was conducted using Genome Studio Genotyping Module v1.8.4 software (Illumina, San Diego, CA, USA), with all allele calls subsequently verified visually. Data were filtered by excluding SNPs with more than 25% missing calls and minor allele frequencies less than 0.05, resulting in 4633 high-quality SNPs for further analysis.

4.5. Genome-Wide Association Study

The GWAS was conducted using adjusted genotypic values for the evaluated traits. These values were derived from REML/BLUP analyses for the years 2019, 2020, and 2021, as well as from the harmonic mean of relative performance of genotypic values (HMRPVG) for the combined analysis of these agricultural years. The population structure with K = 2 for the GWAS was determined using a Bayesian clustering approach with Structure v2.3.4 [47]. The analysis was configured with a burn-in period of 10,000 and 50,000 Markov Chain Monte Carlo (MCMC) iterations, exploring cluster numbers (K) from 2 to 11 across 20 independent runs. Analysis was performed using R software (version 4.0.3) [48] and the multilocus mixed models tool mrMLM.GUI (version 4.0.2) [49], employing five methods: mrMLM [50], FASTmrMLM [51], FASTmrEMMA [52], ISIS EM-BLASSO [53], and pLARmEB [21] for SNP identification. A significance threshold of LOD ≥ 3 was applied across all methods. To reduce false positive associations and enhance statistical power, both population structure and kinship matrix were incorporated into the models. The kinship matrix was calculated using mrMLM.GUI v4.0.2. Model adequacy was evaluated through Q-Q plots, which displayed the observed quantile distribution of –log10 (p) on the y-axis versus the expected quantile distribution of –log10 (p) on the x-axis.

4.6. Candidate Gene Identification

Candidate genes within the genomic regions associated with the traits of interest were identified by aligning these regions with the G19833 reference genome v1.0 [11] through the National Center for Biotechnology Information database available at https://www.ncbi.nlm.nih.gov (accessed on 4 April 2023). Functional annotations of each model gene were confirmed using the Phytozome database available at https://phytozome-next.jgi.doe.gov (accessed on 5 April 2023) [30] to predict their potential roles in morphoagronomic traits. This process involved linking the genes to specific biological functions and pathways, thereby aiding in the understanding of their contribution to the traits under investigation.

5. Conclusions

This study successfully identified 36 SNPs associated with morphoagronomic traits across all chromosomes (Pv01 to Pv11) in common bean accessions of both Mesoamerican and Andean origins. These SNPs offer valuable insights for improving the productivity and quality of common bean cultivars. Moreover, candidate genes were identified for the most significant SNPs, providing a strong foundation for future research. Ultimately, these findings have the potential to contribute to the development of more productive and resilient common bean varieties, thereby advancing breeding programs and enhancing overall crop performance.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13182638/s1, Table S1. Phenotypic data from common beans accessions from Nupagri-UEM germplasm bank evaluated during the years of 2019, 2020, and 2021. Table S2 Accessions from Nupagri-UEM germoplasm bank, gene pool, origin, and market class evaluated in the present study.

Author Contributions

Design and supervision of the research: S.M.A., G.F.L., A.G.V.R. and P.S.V.F.; conduct of the experiments: S.M.A.; analyzed the data: S.M.A., G.F.L., M.V.B. and M.C.G.-V.; original draft of the manuscript: S.M.A. and G.F.L.; revision and editing of the manuscript: M.V.B., A.G.V.R., M.C.G.-V. and P.S.V.F. All authors have read and agreed to the published version of the manuscript.

Funding

Federal Funding Institution National Council for Scientific and Technological Development (CNPq), Grant/Award Numbers: 317493/2021-3 and 403505/2023-2; and Coordination for the Improvement of Higher Education Personnel (Capes), Award number 88887.470337/2019-00 and 88882.448853/2019-01.

Data Availability Statement

All data are presented within the article.

Acknowledgments

S.M.A. is grateful for the grant from CAPES (Brasilia, Brazil) and for the support and assistance provided by professors and students of the Núcleo de Pesquisa Aplicada à Agricultura (Nupagri) and Centro de Treinamento em de Irrigação (CTI) of the Universidade Estadual de Maringá (UEM).

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Food and Agriculture Organization of the United Nations (FAO). Statistical Database, Food and Agriculture. Crops: Beans, Dry. 2022. Available online: http://www.fao.org/faostat/en/#data/QC/visualize (accessed on 30 May 2024).
  2. Thompson, H.J.; McGinley, J.N.; Neil, E.S.; Brick, M.A. Beneficial Effects of Common Bean on Adiposity and Lipid Metabolism. Nutrients 2017, 9, 998. [Google Scholar] [CrossRef] [PubMed]
  3. Jensen, E.S.; Peoples, M.B.; Boddey, R.M.; Gresshoff, P.M.; Hauggaard-Nielsen, H.; Alves, B.J.R.; Morrison, M.J. Legumes for mitigation of climate change and the provision of feedstock for biofuels and biorefineries. A review. Agron. Sustain. Dev. 2012, 32, 329–364. [Google Scholar] [CrossRef]
  4. Jeuffroy, M.H.; Baranger, E.; Carrouée, B.; De Chezelles, E.; Gosme, M.; Hénault, C.; Schneider, A.; Cellier, P. Nitrous oxide emissions from crop rotations including wheat, oilseed rape and dry peas. Biogeosciences 2013, 10, 1787–1797. [Google Scholar] [CrossRef]
  5. Burle, M.L.; Fonseca, J.R.; Kami, J.A.; Gepts, P. Microsatellite diversity and genetic structure among common bean (Phaseolus vulgaris L.) landraces in Brazil, a secondary center of diversity. Theor. Appl. Genet. 2010, 121, 801–813. [Google Scholar] [CrossRef]
  6. Blair, M.W.; Brondani, R.V.P.; Diaz, L.M.; Del Peloso, M.J. Diversity and population structure of common bean from Brazil. Crop Sci. 2013, 53, 1983–1993. [Google Scholar] [CrossRef]
  7. Valdisser, P.A.M.R.; Pereira, W.J.; Almeida Filho, J.E.; Müller, B.S.F.; Coelho, G.R.C.; De Menezes, I.P.P.; Vianna, J.P.G.; Zucchi, M.I.; Lanna, A.C.; Coelho, A.S.G.; et al. In-depth genome characterization of a Brazilian common bean core collection using DArTseq high-density SNP genotyping. BMC Genom. 2017, 18, 423. [Google Scholar] [CrossRef]
  8. Valentini, G.; Gonçalves-Vidigal, M.C.; Elias, J.C.F.; Moiana, L.D.; Mindo, N.N.A. Population Structure and Genetic Diversity of Common Bean Accessions from Brazil. Plant Mol. Biol. Rep. 2018, 36, 897–906. [Google Scholar] [CrossRef]
  9. de Carvalho Paulino, J.F.; de Almeida, C.P.; Song, J.; Carbonell, S.A.M.; Chioroto, A.F.; Benchimol-Reis, L.L. Genetic diversity and inter-gene pool introgression of Mesoamerican Diversity Panel in common beans. J. Appl. Genet. 2021, 62, 585–600. [Google Scholar] [CrossRef]
  10. Gepts, P.; Papa, R. Evolution during Domestication. Encycl. Sci. 2003. [Google Scholar] [CrossRef]
  11. Schmutz, J.; Mcclean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Jenkins, J.; Shu, S.; Song, Q.; Chavarro, C.; et al. A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef]
  12. Bitocchi, E.; Rau, D.; Bellucci, E.; Rodriguez, M.; Murgia, M.L.; Gioia, T.; Santo, D.; Nanni, L.; Attene, G.; Papa, R. Beans (Phaseolus spp.) as a model for understanding crop evolution. Front. Plant Sci. 2017, 8, 722. [Google Scholar] [CrossRef] [PubMed]
  13. Kamfwa, K.; Cichy, K.A.; Kelly, J.D. Genome-wide association study of agronomic traits in common bean. Plant Genome 2015, 8, plantgenome2014-09. [Google Scholar] [CrossRef] [PubMed]
  14. Resende, R.T.; De Resende, M.D.V.; Azevedo, C.F.; Fonseca e Silva, F.; Melo, L.C.; Pereira, H.S.; Souza, T.L.P.O.; Valdisser, P.A.M.R.; Brondani, C.; Vianello, R.P. Genome-wide association and regional heritability mapping of plant architecture, lodging and productivity in Phaseolus vulgaris. G3 Genes|Genomes|Genet. 2018, 8, 2841–2854. [Google Scholar] [CrossRef]
  15. Delfini, J.; Moda-Cirino, V.; Dos Santos Neto, J.; Zeffa, D.M.; Nogueira, A.F.; Ribeiro, L.A.B.; Ruas, P.M.; Gepts, P.; Gonçalves, L.S.A. Genome-Wide Association Study Identifies Genomic Regions for Important Morpho-Agronomic Traits in Mesoamerican Common Bean. Front. Plant Sci. 2021, 12, 748829. [Google Scholar] [CrossRef] [PubMed]
  16. Ikram, M.; Han, X.; Zuo, J.-F.; Song, J.; Han, C.-Y.; Zhang, Y.-W.; Zhang, Y.-M. Identification of QTNs and their candidate genes for 100-seed weight in soybean (Glycine max L.) using multi-locus genome-wide association studies. Genes 2020, 11, 714. [Google Scholar] [CrossRef]
  17. Pérez-de-Castro, A.M.; Vilanova, S.; Cañizares, J.; Pascual, L.; Blanca, J.M.; Diez, M.J.; Prohens, J.; Pico, B. Application of genomic tools in plant breeding. Curr. Genom. 2012, 13, 179–195. [Google Scholar] [CrossRef]
  18. Korte, A.; Farlow, A. The advantages and limitations of trait analysis with GWAS: A review. Plant Methods 2013, 9, 29. [Google Scholar] [CrossRef]
  19. Burghardt, L.T.; Young, N.D.; Tiffin, P. A guide to genome-wide association mapping in plants. Curr. Protoc. Plant Biol. 2017, 2, 22–38. [Google Scholar] [CrossRef]
  20. Wang, M.; Jiang, N.; Jia, T.; Leach, L.; Cockram, J.; Waugh, R.; Ramsay, L.; Thomas, B.; Luo, Z. Genome-wide association mapping of agronomic and morphologic traits in highly structured populations of barley cultivars. Theor. Appl. Genet. 2012, 124, 233–246. [Google Scholar] [CrossRef]
  21. Zhang, J.; Feng, J.Y.; Ni, Y.L.; Wen, Y.J.; Niu, Y.; Tamba, C.L.; Yue, C.; Song, Q.; Zhang, Y.-M. pLARmEB: Integration of least angle regression with empirical Bayes for multilocus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef]
  22. Cui, Y.; Zhang, F.; Zhou, Y. The application of multi-locus GWAS for the detection of salt-tolerance loci in rice. Front. Plant Sci. 2018, 9, 1464. [Google Scholar] [CrossRef] [PubMed]
  23. Su, J.; Ma, Q.; Mei, L.; Hao, F.; Wang, C. Multi-Locus Genome-Wide Association Studies of Fiber-Quality Related Traits in Chinese Early-Maturity Upland Cotton. Front. Plant Sci. 2018, 9, 1169. [Google Scholar] [CrossRef] [PubMed]
  24. Ward, B.P.; Brown-Guedira, G.; Kolb, F.L.; Van Sanford, D.A.; Tyagi, P.; Sneller, C.H.; Griffey, C.A. Genome-wide association studies for yield-related traits in soft red winter wheat grown in Virginia. PLoS ONE 2019, 14, e0208217. [Google Scholar] [CrossRef]
  25. Han, X.; Xu, Z.-R.; Zhou, L.; Han, C.-Y.; Zhang, Y.-M. Identification of QTNs and their candidate genes for flowering time and plant height in soybean using multi-locus genome-wide association studies. Mol. Breed. 2021, 41, 39. [Google Scholar] [CrossRef]
  26. Nemli, S.; Asciogul, T.K.; Kaya, H.B.; Kahraman, A.; Eşiyok, D.; Tanyolac, B. Association mapping for five agronomic traits in the common bean (Phaseolus vulgaris L.). J. Sci. Food Agric. 2014, 94, 3141–3151. [Google Scholar] [CrossRef]
  27. Moghaddam, S.M.; Mamidi, S.; Osorno, J.M.; Lee, R.; Brick, M.; Kelly, J.; Miklas, P.; Urrea, C.; Song, Q.; Cregan, P.; et al. Genome-Wide Association Study Identifies Candidate Loci Underlying Agronomic Traits in A Middle American Diversity Panel of Common Bean. Plant Genome 2016, 9, plantgenome2016-02. [Google Scholar] [CrossRef]
  28. Almeida, C.P.D.E.; Paulino, J.F.D.E.C.; Carbonell, S.A.M.; Chiorato, A.F.; Song, Q.; Di Vittori, V.; Rodriguez, M.; Papa, R.; Benchimol-Reis, L.L. Genetic Diversity, Population Structure, and Andean Introgression in Brazilian Common Bean Cultivars after Half a Century of Genetic Breeding. Genes 2020, 11, 1298. [Google Scholar] [CrossRef] [PubMed]
  29. Mutari, B.; Sibiya, J.; Shayanowako, A.; Chidzanga, C.; Matova, P.M.; Gasura, E. Identification of genomic regions of dry bean (Phaseolus vulgaris L.) associated with agronomic and physiological traits under drought stressed and well-watered conditions using genome-wide association study. BioRxiv 2022, 20, 2022-11. [Google Scholar] [CrossRef]
  30. Goodstein, D.M.; Shu, S.; Howson, R.; Neupane, R.; Hayes, R.D.; Fazo, J.; Mitros, T.; Dirks, W.; Hellsten, U.; Putnam, N.; et al. Phytozome: A comparative platform for green plant genomics. Nucleic. Acids Res. 2012, 40, D1178–D1186. [Google Scholar] [CrossRef]
  31. Elias, J.C.F.; Gonçalves-Vidigal, M.C.; Vaz Bisneta, M.; Valentini, G.; Vidigal Filho, P.S.; Gilio, T.A.S.; Moda-Cirino, V.; Song, Q. Genetic Mapping for Agronomic Traits in IAPAR 81/LP97-28 Population of Common Bean (Phaseolus vulgaris L.) under Drought Conditions. Plants 2021, 10, 1568. [Google Scholar] [CrossRef]
  32. Zhang, Y.-M.; Jia, Z.; Dunwell, J.M. The applications of new multi-locus GWAS methodologies in the genetic dissection of complex traits. Front. Plant Sci. 2019, 10, 100. [Google Scholar] [CrossRef]
  33. Yang, Y.; Chai, Y.; Zhang, X.; Lu, S.; Zhao, Z.; Wei, D.; Chen, L.; Hu, Y.-G. Multi-locus GWAS of quality traits in bread wheat: Mining more candidate genes and possible regulatory network. Front. Plant Sci. 2020, 11, 1091. [Google Scholar] [CrossRef]
  34. Koinange, E.M.K.; Singh, S.P.; Gepts, P. Genetic control of the domestication syndrome in common bean. Crop Sci. 1996, 36, 1037–1045. [Google Scholar] [CrossRef]
  35. Blair, M.W.; Iriarte, G.; Beebe, S. QTL analysis of yield traits in an advanced backcross population derived from a cultivated Andean x wild common bean (Phaseolus vulgaris L.) cross. Theor. Appl. Genet. 2006, 112, 1149–1163. [Google Scholar] [CrossRef]
  36. Blair, M.W.; Díaz, L.M.; Buendía, H.F.; Duque, M.C. Genetic Diversity, seed size associations and population structure of a core collection of common beans (Phaseolus vulgaris L.). Theor. Appl. Genet. 2009, 119, 955–972. [Google Scholar] [CrossRef] [PubMed]
  37. Pérez-Veja, E.; Pañeda, A.; Rodríguez-Suárez, C.; Campa, A.; Giraldez, R.; Ferreira, J. Mapping of QTLs for morpho-agronomic and seed quality traits in a RIL population of common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 2010, 120, 1367–1380. [Google Scholar] [CrossRef]
  38. González, A.M.; De La Fuente, M.; De Ron, A.M.; Santalla, M. Protein markers and seed size variation in common bean segregating populations. Mol. Breed. 2010, 25, 723–740. [Google Scholar] [CrossRef]
  39. Roberts, A.W.; Bushoven, J.T. The cellulose synthase (CESA) gene superfamily of the moss Physcomitrella patens. BMC Plant Biol. 2007, 63, 207–2019. [Google Scholar] [CrossRef]
  40. Zhang, B.; Gao, Y.; Zhang, L.; Zhou, Y. The plant cell wall: Biosynthesis, construction, and functions. J. Int. Plant Biol. 2021, 63, 251–272. [Google Scholar] [CrossRef]
  41. Nkhata, W.; Shimelis, H.; Melis, R.; Chirwa, R.; Mzengeza, T.; Mathew, I.; Shayanowako, A. Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean. PLoS ONE 2021, 16, e0250729. [Google Scholar] [CrossRef]
  42. Wang, H.; Wang, Y.; Sang, T.; Lin, Z.; Li, R.; Ren, W.; Shen, X.; Zhao, B.; Wang, X.; Zhang, X.; et al. Cell type-specific proteomics uncovers a RAF15-SnRK2.6/OST1 Kinase cascade in guard cells. J. Int. Plant Biol. 2023, 65, 2122–2137. [Google Scholar] [CrossRef] [PubMed]
  43. MacQueen, A.H.; White, J.W.; Lee, R.; Osorno, J.M.; Schmutz, J.; Miklas, P.N.; Myers, J.I.M.; McClean, P.; Thomas, E. Genetic associations in four decades of multi-environment trials reveal agronomic trait evolution in common bean. Genetics 2020, 215, 267–284. [Google Scholar] [CrossRef] [PubMed]
  44. Havé, M.; Marmagne, A.; Chardon, F.; Masclaux-Daubresse, C. Nitrogen remobilization during leaf senescence: Lessons from Arabidopsis to crops. J. Exp. Bot. 2017, 68, 2513–2529. [Google Scholar] [CrossRef] [PubMed]
  45. Resende, M.D.V. Software Selegen-REML/BLUP: A useful tool for plant breeding. Crop Breed. Appl. Biotechnol. 2016, 16, 330–339. [Google Scholar] [CrossRef]
  46. Song, Q.; Jia, G.; Hyten, D.L.; Jenkins, J.; Hwang, E.Y.; Schroeder, S.G.; Osorno, J.M.; Schumtz, J.; Jackson, S.A.; Mcclean, P.E.; et al. SNP Assay Development for Linkage Map Construction, Anchoring Whole Genome Sequence and Other Genetic and Genomic Applications in Common Bean. G3 Genes|Genomes|Genet. 2015, 5, 2285–2290. [Google Scholar] [CrossRef]
  47. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef] [PubMed]
  48. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistics Computing: Vienna, Austria, 2020; Available online: https://www.r-project.org/ (accessed on 30 May 2022).
  49. Zhang, Y.-W.; Tamba, C.L.; Wen, Y.-J.; Pei, L.; Ren, W.-L.; Ni, Y.-L.; Gao, J.; Zhang, Y.-M. mrMLM v4.0.2: An R platform for multi-locus genome-wide association studies. Genom. Proteom. Bioinform. 2020, 18, 481–487. [Google Scholar] [CrossRef]
  50. Wang, S.B.; Feng, J.Y.; Ren, W.L.; Huang, B.; Zhou, L.; Wen, Y.J.; Zhang, J.; Dunwell, J.M.; Xu, S.; Zhang, Y.-M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 2016, 6, 19444. [Google Scholar] [CrossRef]
  51. Tamba, C.L.; Zhang, Y.M. A FAST mrMLM algorithm for multi-locus genome-wide association studies. BioRxiv 2018, 1–34. [Google Scholar] [CrossRef]
  52. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wang, S.-B.; Dunwell, J.M.; Zhang, Y.-M.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief Bioinform. 2018, 19, 700–712. [Google Scholar] [CrossRef]
  53. Tamba, C.L.; Ni, Y.L.; Zhang, Y.M. Iterative sure Independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Precipitation (mm) and minimum and maximum temperature (degrees Celsius) in the experiments conducted at the Technical Irrigation Center (CTI) in Maringá, PR, Brazil, in 2019, 2020, and 2021.
Figure 1. Precipitation (mm) and minimum and maximum temperature (degrees Celsius) in the experiments conducted at the Technical Irrigation Center (CTI) in Maringá, PR, Brazil, in 2019, 2020, and 2021.
Plants 13 02638 g001
Figure 2. Manhattan and QQ plots obtained in GWAS for the yield and its components and plant architecture: (a) YDSD, grain yield (kg ha−1); (b) SW, 100-seed weight (g); (c) PDPL, number of pods per plant; (d) FPIH, first pod insertion height (cm); (e) PLHT, plant height (cm). The gray threshold line was considered for LOD score ≥ 3. The pink dots indicate the SNPs significantly associated with traits.
Figure 2. Manhattan and QQ plots obtained in GWAS for the yield and its components and plant architecture: (a) YDSD, grain yield (kg ha−1); (b) SW, 100-seed weight (g); (c) PDPL, number of pods per plant; (d) FPIH, first pod insertion height (cm); (e) PLHT, plant height (cm). The gray threshold line was considered for LOD score ≥ 3. The pink dots indicate the SNPs significantly associated with traits.
Plants 13 02638 g002
Figure 3. Manhattan and QQ plots obtained in GWAS for phenology: (a) DF, number of days for flowering; (b) DPM, number of days for maturity. The gray threshold line was considered for LOD score ≥ 3. The pink dots indicate the SNPs significantly associated with traits.
Figure 3. Manhattan and QQ plots obtained in GWAS for phenology: (a) DF, number of days for flowering; (b) DPM, number of days for maturity. The gray threshold line was considered for LOD score ≥ 3. The pink dots indicate the SNPs significantly associated with traits.
Plants 13 02638 g003
Figure 4. Accessions of common beans from the Nupagri Germplasm Bank (photos taken by UEM/ASC).
Figure 4. Accessions of common beans from the Nupagri Germplasm Bank (photos taken by UEM/ASC).
Plants 13 02638 g004
Table 1. Traits analyzed in common bean accessions belonging to the Nupagri Germplasm Bank evaluated over three agricultural years, Maringá, PR, Brazil.
Table 1. Traits analyzed in common bean accessions belonging to the Nupagri Germplasm Bank evaluated over three agricultural years, Maringá, PR, Brazil.
TraitsYearAccessions§PEV
AverageIntervalLRT (x2)
Yield and its Components
Grain yield
(kg ha−1)
2019990.44191.63–1798.61215.24 ****10,221.53
20201736.26335.15–3723.92193.41 **47,059.12
20211093.57292.21–2469.66149.20 **30,490.65
100-seed weight (g)201928.0015.22–49.73611.54 ****0.99
202027.9216.21–54.63565.88 **1.27
202128.1014.15–44.50421.76 **2.32
Number of seeds per pod20193.212.12–4.54135.70 ****0.05
20203.601.91–5.99154.82 **0.09
20212.732.05–3.5238.85 **0.05
Number of pods per plant201910.044.34–15.14119.59 ****1.37
20209.673.50–17.35119.90 **1.32
202110.425.45–21.0060.21 **3.81
Architecture
First pod insertion height (cm)201913.3411.32–18.14150.65 ****0.30
202013.5311.21–21.71152.08 **0.56
202113.1411.19–16.0151.08 **0.44
Plant height (cm)201949.0037.56–58.63309.44 ****1.17
202050.1236.17–60.80492.35 **0.68
202147.8837.77–58.08134.13 **3.33
Phenology
Number of days for flowering201942.7030.99–60.991208.73 ****0.06
202042.2631.73–60.96925.23 **0.22
202142.7130.00–60.001281.12 **0.0009
Number of days for maturity201981.4673.65–89.59360.75 ****1.03
202081.4973.57–91.30292.75 **1.32
202180.8973.28–87.69263.45 **2.35
** = significant by the chi-square test with one degree of freedom: at 1% (6.63); LRT (x2) = likelihood-ratio test (chi-square); §  = interaction between genotype treatments and years; PEV = variance in the error of prediction for genotypic values.
Table 2. SNPs associated with yield and its components, plant architecture, and phenology identified using different methods in common bean accessions from the Nupagri Germplasm Bank, Maringá, PR, Brazil.
Table 2. SNPs associated with yield and its components, plant architecture, and phenology identified using different methods in common bean accessions from the Nupagri Germplasm Bank, Maringá, PR, Brazil.
§ TraitsYearChrPosition (pb)SNP# LOD Score−Log10 P¥ SNP Effect R2 (%)Methods
YDSD2019248,720,860ss7156479974.555.33−147.00~−124.6412.43~17.291, 5
2020425,245,150ss7156509184.725.50−513.05~−361.059.81~16.551, 2
739,295,293ss7156404875.35~6.106.16~6.94−303.83~−235.6811.69~17.701, 2, 4, 5
855,278,106ss7156393593.994.74148.49~238.745.73~12.371, 2
202155,399,253ss7156486724.635.41−299.91~−150.7312.70~12.963, 4
C *248,720,860ss7156479974.855.64−203.83~−175.2513.38~18.321, 4
SW20197819,750ss7156456893.21~9.193.92~10.11−7.55~−2.9913.48~21.512, 3, 4
1024,598,176ss7156415433.46~4.954.18~5.74−6.90~−3.7910.13~21.713, 4
202023,118,774ss7156395023.41~9.104.13~10.01−6.48~−3.382.08~29.971, 3, 4
540,128,137ss7156466975.01~5.905.81~6.732.75~6.4711.28~15.642, 3
7819,750ss7156456895.75~10.156.57~11.09−7.30~−3.4517.83~25.232, 3, 4, 5
2021933,444,042ss7156456026.41~7.897.26~8.78−3.03~−2.035.30~11.751, 4, 5
C *622,116,966ss7156397864.49~30.635.26~31.54−6.21~−3.5218.43~57.294, 5
7819,750ss7156456893.10~3.383.80~4.10−9.49~−2.348.15~33.452, 3, 4
933,444,042ss7156456025.17~7.425.98~8.29−2.13~−1.964.83~5.744, 5
SDPD2020627,760,992ss7156456733.77~4.004.51~4.76−0.35~−0.324.92~8.241, 5
PDPL2019151,819,821ss7156453013.22~15.143.93~16.160.80~1.7710.11~43.391, 2, 3, 4, 5
248,720,860ss7156479974.80~8.155.59~9.04−0.97~−0.888.83~10.741, 4, 5
1039,797,018ss7156463303.18~4.093.89~4.85−1.03~−0.612.92~8.131, 2
20211039,797,018ss7156463303.13~5.643.83~6.46−2.38~−0.955.07~15.241, 2, 3, 4, 5
1148,780,038ss7156507486.69~6.867.54~7.72−1.59~−1.129.46~19.151, 4
PDPLC *151,819,821ss7156453013,54~14,894.27~15.910.95~1.9212.22~43.101, 2, 3, 4, 5
248,720,860ss7156479974.56~6.895.34~7.76−1.04~−0.998.61~9.361, 4, 5
FPIH201934,083,079ss7156479654.86~7.705.65~8.59−0.46~−0.304.32~12.854, 5
41,627,690ss7156499714.03~9.594.78~10.520.36~0.637.00~21.681, 2, 4, 5
625,791,849ss7156457523.09~6.383.79~7.23−0.83~−0.569.98~21.741, 2
931,832,898ss7156476264.23~5.394.99~6.210.42~0.507.12~10.122, 4, 5
2020540,141,917ss7156466993.53~4.094.25~4.85−0.95~−0.887.20~8.471, 2
625,791,849ss7156457524.45~5.825.22~6.65−2.31~−0.9414.47~19.871, 2, 3, 5
202141,627,690ss7156499713.33~4.374.05~5.140.29~0.738.98~23.461, 2, 3, 4
C *41,503,482ss7156469103.93~4.414.67~5.18−1.05~−0.599.25~10.923, 4
41,627,690ss7156499714.02~6.354.77~7.200.39~0.527.02~12.511, 2
625,791,849ss7156457523.04~6.303.74~7.15−0.89~−0.596.28~21.511, 2, 4
PLHT2019145,894,030ss7156473685.24~15.566.05~16.59−6.00~−2.0418.04~43.651, 2, 3, 4, 5
149,657,760ss7156458523.90~4.064.65~4.82−1.25~−1.587.43~10.311, 2, 5
41,224,240ss7156468966.177.00−1.69~−2.2514.20~21.931, 2
41,982,297ss7156478194.955.751.53~2.029.91~15.001, 2
2020145,894,030ss7156473684.49~12.035.27~13.01−6.18~−3.0929.53~33.982, 3, 4, 5
2021145,894,030ss7156473684.58~8.205.35~9.10−5.17~−1.9420.00~30.531, 2, 3, 4, 5
245,319,921ss7156459644.73~5.025.51~5.82−1.82~−1.287.57~15.231, 2, 4
73,648,568ss7156463533.253.96−1.26~−0.823.45~8.131, 2
C *145,894,030ss7156473685.29~15.616.10~16.64−6.16~−2.1118.12~38.651, 2, 3, 4, 5
149,657,760ss7156458523.81~4.014.55~4.76−1.62~−1.267.10~10.161, 2, 5
41,224,240ss7156468965.11~6.265.91~7.10−2.33~−1.378.59~22.181, 2, 4
41,982,297ss7156478193.74~4.994.48~5.781.14~2.085.04~15.051, 2, 4
DF2019145,894,030ss7156473683.834.57−3.75~−2.5612.80~27.401, 2
146,027,600ss71563953611.42~11.9412.38~12.91−3.68~−3.5324.52~26.654, 5
73,770,008ss7156463553.02~3.983.71~4.73−2.19~−1.131.49~5.601, 2, 5
107,867,881ss71564011610.39~12.4811.35~13.46−3.35~−3.9823.79~33.551, 2, 5
2020145,894,030ss7156473683.514.24−3.90~−2.5310.82~25.801, 2
146,027,600ss7156395367.80~11.418.69~12.37−3.41~−3.2217.46~19.894, 5
73,770,008ss7156463554.134.88−2.37~−1.743.07~5.721, 2
107,867,881ss71564011611.0412.003.69~4.3925.10~35.501, 2
2021146,027,600ss71563953610.57~16.5611.52~17.60−3.93~−3.3023.26~33.104, 5
73,770,008ss7156463553.48~3.694.20~4.42−1.33~−1.281.93~2.271, 2
748,520,131ss7156485804.10~5.644.85~6.472.50~3.323.01~5.304, 5
107,867,881ss7156401167.94~9.198.83~10.112.50~5.1714.42~56.841, 2, 5
C *145,894,030ss7156473683.894.63−3.79~−2.5712.71~27.531, 2
146,027,600ss71563953611.14~12.1512.11~13.13−3.59~−3.5724.72~24.974, 5
73,770,008ss7156463553.02~3.993.72~4.75−2.19~−1.131.48~5.331, 2, 5
107,867,881ss71564011610.54~12.6511.48~13.643.57~4.0223.78~33.671, 2, 5
DPM2019145,746,595ss7156392717.39~7.708.27~8.58−2.80~−2.2221.53~22.621, 2
146,027,600ss7156395363.42~7.594.14~8.47−2.31~−2.2422.66~24.254, 5
DPM201933,963,582ss7156476363.03~4.843.73~5.630.83~1.392.44~6.834, 5
44,497,307ss7156495284.93~7.785.73~8.672.52~2.7929.33~32.581, 2
625,619,371ss7156484925.99~6.076.82~6.90−1.99~−1.6412.34~16.311, 2
627,370,471ss7156456773.27~5.773.98~5.601.05~2.155.01~5.223, 4
2020146,027,600ss7156395363.75~7.154.49~8.01−2.46~−1.5813.04~22.071, 2, 5
34,083,079ss7156479654.33~5.075.09~5.87−2.65~−1.8415.28~24.251, 2
1039,577,266ss7156463184.35~5.165.12~5.962.12~2.8120.14~26.891, 2
2021145,894,030ss7156473684.48~6.055.25~6.88−3.47~−1.748.44~10.711, 2, 3, 4
44,497,307ss7156495284.90~16.145.69~17.182.97~5.9530.35~51.061, 2, 3, 4
625,619,371ss7156484923.27~3.803.99~4.54−1.83~−1.315.50~9.451, 2
C *145,746,595ss7156392718.22~12.349.11~13.33−3.18~−2.2622.92~45.081, 2
33,963,582ss7156476363.31~4.334.02~5.201.18~1.954.69~12.811, 5
625,619,371ss7156484923.51~6.284.23~7.13−1.63~−1.125.51~11.861, 2
627,370,471ss7156456773.22~3.703.92~4.430.93~2.193.70~5.153, 4
§ Traits: YDSD, grain yield (kg ha−1); SW, 100-seed weight (g); SDPD, number of seeds per pod; PDPL, number of pod per plant; FPIH, first pod insertion height (cm); PLHT, plant height (cm); DF, number of days for flowering; DPM, number of days for maturity; C *, combined analysis of years using the harmonic mean of relative performance of genotypic values (HMRPVG); Chr, chromosome; # LOD value ≥ 3; ¥ SNP effect, Single-nucleotide polymorphism effect on evaluated traits found through GWAS using multilocus mixed models; indicates the percentage of phenotypic variation explained by each SNP; methods: 1, mrMLM; 2, FASTmrMLM; 3, FASTmrEMMA; 4, pLARmEB; 5, ISIS EM-BLASSO.
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

Alves, S.M.; Lacanallo, G.F.; Gonçalves-Vidigal, M.C.; Vaz Bisneta, M.; Vidigal Rosenberg, A.G.; Vidigal Filho, P.S. Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions. Plants 2024, 13, 2638. https://doi.org/10.3390/plants13182638

AMA Style

Alves SM, Lacanallo GF, Gonçalves-Vidigal MC, Vaz Bisneta M, Vidigal Rosenberg AG, Vidigal Filho PS. Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions. Plants. 2024; 13(18):2638. https://doi.org/10.3390/plants13182638

Chicago/Turabian Style

Alves, Stephanie Mariel, Giselly Figueiredo Lacanallo, Maria Celeste Gonçalves-Vidigal, Mariana Vaz Bisneta, Andressa Gonçalves Vidigal Rosenberg, and Pedro Soares Vidigal Filho. 2024. "Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions" Plants 13, no. 18: 2638. https://doi.org/10.3390/plants13182638

APA Style

Alves, S. M., Lacanallo, G. F., Gonçalves-Vidigal, M. C., Vaz Bisneta, M., Vidigal Rosenberg, A. G., & Vidigal Filho, P. S. (2024). Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions. Plants, 13(18), 2638. https://doi.org/10.3390/plants13182638

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