Background: A large variation in seed coat colors and seed phenolic metabolites is present in common bean (Phaseolus vulgaris L.). The study of the relationships between seed coat color phenotype and the phenolic profile is an important step in the elucidation of the gene network involved in the phenylpropanoid biosynthetic pathway. However, this relationship is still poorly understood in this species. Results: A genome-wide association study (GWAS) was used to investigate the genomic regions associated with the synthesis of 10 flavonoids (5 anthocyanins and 5 flavonols) and with 10 seed coat color traits using a set of 308 common bean lines of the Spanish Diversity Panel (SDP) which have been genotyped with 11,763 SNP markers.. A total of 31 significant SNP-trait associations (QTNs) were identified, grouped in 20 chromosome regions: 6 for phenolic metabolites on chromosomes Pv01, Pv02, Pv04, Pv08, and Pv09, 13 for seed coat color on chromosomes Pv01, Pv02, Pv06, Pv07, and Pv10, and 1 including both types of traits located on chromosome Pv08. In all, 58 candidate genes underlying these regions have been proposed, 31 of them previously described in the phenylpropanoid pathway in common bean, and 27 of them newly proposed in this work based on the association study and their homology with Arabidopsis anthocyanin genes. Conclusions: Chromosome Pv08 was identified as the main chromosome involved in the phenylpropanoid pathway and in consequence in the common bean seed pigmentation, with three independent chromosome regions identified, Phe/C_Pv08(2.7) (expanding from 2.71 to 4.04 Mbp), C_Pv08(5.8) (5.89–6.59 Mbp), and Phe_Pv08(62.5) (62.58 to 63.28 Mbp). Candidate genes previously proposed by other authors for the color genes V and P were validated in this GWAS. Candidate genes have been tentatively proposed from this study for color genes B and Rk on Pv02, Asp on Pv07, and complex C on Pv08. These results help to clarify the complex network of genes involved in the genetic control of phenolic compounds and seed color in common bean and provide the opportunity for future validation studies.

Snap beans are a group of bean cultivars grown for their edible immature pods. The objective of this work was to characterize the diversity of pod phenotypes in a snap bean panel (SBP), comprising 311 lines collected in Europe, and establish a core set (Core-SBP) with the maximum diversity of pod phenotypes. Phenotyping of the SBP was carried out over two seasons based on 14 quantitative pod dimension traits along with three qualitative traits: pod color, seed coat color, and growth habit. Phenotypes were grouped into 54 classes using a hierarchical method, and a Core-SBP with one line per phenotype class was established. A further field-based evaluation of the Core-SBP revealed higher diversity index values than those obtained for the SBP. The Core-SBP was also genotyped using 24 breeder-friendly DNA markers tagging 21 genomic regions previously associated with pod trait control. Significant marker-trait associations were found for 11 of the 21 analyzed regions as well as the locus fin. The established Core-SBP was a first attempt to classify snap bean cultivars based on pod morphology and constituted a valuable source of characteristics for future breeding programs and genetic analysis.

Key message: QTL mapping, association analysis, and colocation study with previously reported QTL revealed three main regions controlling pod morphological traits and two loci for edible pod characteristics on the common bean chromosomes Pv01 and Pv06. Bean pod phenotype is a complex characteristic defined by the combination of different traits that determine the potential use of a genotype as a snap bean. In this study, the TUM RIL population derived from a cross between ‘TU’ (dry) and ‘Musica’ (snap) was used to investigate the genetic control of pod phenotype. The character was dissected into pod morphological traits (PMTs) and edible pod characteristics (EPC). The results revealed 35 QTL for PMTs located on seven chromosomes, suggesting a strong QTL colocation on chromosomes Pv01 and Pv06. Some QTL were colocated with previously reported QTL, leading to the mapping of 15 consensus regions associated with bean PMTs. Analysis of EPC of cooked beans revealed that two major loci with epistatic effect, located on chromosomes Pv01 and Pv06, are involved in the genetic control of this trait. An association study using a subset of the Spanish Diversity Panel (snap vs. non-snap) detected 23 genomic regions, with three regions being mapped at a position similar to those of two loci identified in the TUM population. The results demonstrated the relevant roles of Pv01 and Pv06 in the modulation of bean pod phenotype. Gene ontology enrichment analysis revealed a significant overrepresentation of genes regulating the phenylpropanoid metabolic process and auxin response in regions associated with PMTs and EPC, respectively. Both biological functions converged in the lignin biosynthetic pathway, suggesting the key role of the pathway in the genetic control of bean pod phenotype.

Key message: Association analysis, colocation study with previously reported QTL, and differential expression analyses allowed the identification of the consistent QTLs and main candidate genes controlling seed traits. Common beans show wide seed variations in shape, size, water uptake, and coat proportion. This study aimed to identify consistent genomic regions and candidate genes involved in the genetic control of seed traits by combining association and differential expression analyses. In total, 298 lines from the Spanish Diversity Panel were genotyped with 4,658 SNP and phenotyped for seven seed traits in three seasons. Thirty-eight significant SNP-trait associations were detected, which were grouped into 23 QTL genomic regions with 1,605 predicted genes. The positions of the five QTL regions associated with seed weight were consistent with previously reported QTL. HCPC analysis using the SNP that tagged these five QTL regions revealed three main clusters with significantly different seed weights. This analysis also separated groups that corresponded well with the two gene pools described: Andean and Mesoamerican. Expression analysis was performed on the seeds of the cultivar ‘Xana’ in three seed development stages, and 1,992 differentially expressed genes (DEGs) were detected, mainly when comparing the early and late seed development stages (1,934 DEGs). Overall, 91 DEGs related to cell growth, signaling pathways, and transcriptomic factors underlying these 23 QTL were identified. Twenty-two DEGs were located in the five QTL regions associated with seed weight, suggesting that they are the main set of candidate genes controlling this character. The results confirmed that seed weight is the sum of the effects of a complex network of loci, and contributed to the understanding of seed phenotype control.

The Fabada market class within the dry beans has a well‐differentiated seed phenotype with very large white seeds. This work investigated the genetic diversity maintained in the seed collections within this market class and possible genetic erosion over the last 30 years. A panel with 100 accessions was maintained in seed collections for 30 years, 57 accessions collected from farmers in 2021, six cultivars developed in SERIDA, and 16 reference cultivars were gathered and genotyped with 108,585 SNPs using the genotyping‐by‐sequencing method. Filtering based on genotypic and phenotypic data was carried out in a staggered way to investigate the genetic diversity among populations. The dendrogram generated from genotyping revealed 90 lines forming 16 groups with identical SNP profiles (redundant lines) from 159 lines classified as market‐class Fabada according to their passport data. Seed phenotyping indicated that 19 lines were mistakenly classified as Fabada (homonymies), which was confirmed in the dendrogram built without redundant lines. Moreover, this study provides evidence of genetic erosion between the population preserved for 30 years and the currently cultivated population. The conserved population contains 54.6% segregation sites and 41 different SNP profiles, whereas the cultivated population has 19.6% segregation sites and 26 SNP profiles. The loss of genetic variability cannot be attributed to the diffusion of modern cultivars, which increase genetic diversity (six new SNP profiles). The results allow for the more efficient preservation of plant genetic resources in genebanks, minimizing redundant accessions and incorporating new variations based on genotypic and phenotypic data.

Fabada is a traditional market class of common bean (Phaseolus vulgaris L.) cultivated in northern Spain, recognized for its distinctive seed phenotype with very large, white, and oblong seeds. High-throughput genotyping of this market class revealed that its genome is predominantly of Andean origin, with approximately 30 % introgression from the Mesoamerican gene pool. The de novo genome assembly of the A25 bean genotype, derived from the local cultivar ‘Andecha’ which belongs to the Fabada market class, is described. Two complementary technologies were used: PacBio (Sequel II) for long reads and Illumina (NovaSeq PE150) for short reads. The depth of coverage achieved was 112x for Illumina and 47x for PacBio. Genome assembly resulted in 468,9 Mbp (98.5 % BUSCO completeness), organized into 1363 scaffolds plus the mitochondrial and chloroplast genomes. Based on the reference genome (NCBI accession number GCF000499845.2), these scaffolds were organized into 11 chromosomes and 169 unplaced scaffolds. The mitochondrial genome was assembled based on that of P. vulgaris (NCBI accession number NC_045135). Two mitochondrial scaffolds were obtained, one of 371,437 bp, and the other of 11,183 bp. The chloroplast genome was assembled based on that of P. vulgaris (NCBI accession NC_009259.1), resulting in 161,310 bp. To the best of our knowledge, this is the first available genome of a common bean accession exhibiting recombination between the two major gene pools, Andean and Mesoamerican, and the second assembled genome of a European common bean.