Genetic Diversity, Population Structure, and Linkage Disequilibrium in a Spanish Common Bean Diversity Panel Revealed through Genotyping-by-Sequencing

Fatty Acids in Dry Beans (Phaseolus vulgaris L.): A Contribution to Their Analysis and the Characterization of a Diversity Panel

Dry bean (Phaseolus vulgaris L.) is a crop of high nutritional interest widespread throughout the world. This research had two objectives. On the one hand, the development and validation of an analytical method to quantify fatty acids in dry beans based on the extraction and derivatization in a single step and later quantification by gas chromatography. On the other, its application to characterize the fatty acid content in a diversity panel consisting of 172 lines. The method was successfully validated in terms of accuracy, precision and robustness. Among the 14 fatty acids that constitute the fatty acid profile of dry bean, the most quantitatively important were linolenic acid, the major fatty acid in all cases, with an average value of 6.7 mg/g, followed by linoleic acid (3.9 mg/g), palmitic acid (2.9 mg/g) and oleic acid (1.5 mg/g). The concentrations of fatty acids in dry bean were influenced by the gene pool, with the Mesoamerican gene pool showing a higher content of palmitic, stearic, linoleic and linolenic acids and the Andean gene pool a higher level of cis-vaccenic acid. Also, the expression of fatty acid content showed high heritability. The information generated constitutes a robust database of interest in food technology, nutrition and breeding programs.

Identification of consistent QTL and candidate genes associated with seed traits in common bean by combining GWAS and RNA-Seq

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.

Factors Influencing the Emergence of Heterogeneous Populations of Common Bean (Phaseolus vulgaris L.) and Their Potential for Intercropping

The common bean is an important legume valued for its protein-rich seeds and its ability to fix nitrogen, making it a key element of crop rotation. In conventional agriculture, the emphasis is on uniformity and genetic purity to optimize crop performance and maximize yields. This is due to both the legal obligations to register varieties and the challenges of implementing breeding programs to create genetically diverse varieties. This paper focuses on the factors that influence the occurrence of heterogeneous common bean populations. The main factors contributing to this diversity have been described, including local adaptations, variable weather conditions, different pollinator species, and intricate interactions between genes controlling seed coat colour. We also discuss the benefits of intercropping common beans for organic farming systems, highlighting the improvement in resistance to diseases, and adverse environmental conditions. This paper contributes to a better understanding of common bean seed heterogeneity and the legal obligation to use heterogeneous populations.

Phaseolus vulgaris mutants reveal variation in the nuclear genome

Phaseolus vulgaris L. (common bean) is an essential source of proteins in the human diet worldwide. Bean breeding programs to increase genetic diversity based on induced mutagenesis have a long tradition in Bulgaria. Common bean varieties with high productivity, wide environmental adaptability, good nutritional properties, and improved disease resistance have been successfully developed. In this study, we aimed to investigate selected nuclear genome features, such as the genome size, the number and chromosomal distribution of 5S and 35S rDNA loci by using the fluorescence in situ hybridization (FISH), as well as the level of DNA damage in some local Bulgarian accessions and mutants of P. vulgaris. Flow cytometry analyses revealed no significant differences in genome size between analyzed lines except for one of the analyzed mutants, M19. The value of genome size 2C DNA is about 1.37 pg2C -1 for all lines, whereas it is 1.42 pg2C-1 for M19. The chromosome number remains the same (2n=22) for all analyzed lines. Results of FISH analyses showed that the number of 5S rDNA was stable among accessions and mutant lines (four loci), while the number of 35S rDNA loci was shown as highly polymorphic, varying between ten and sixteen, and displaying differences in the size and location of 35S rDNA loci between analyzed genotypes. The cell cycle profile was different for the analyzed genotypes. The results revealed that wide variation in genome organization and size as well as DNA damage characterizes the analyzed genetic resources of the common bean.

Mapping yield and yield-related traits using diverse common bean germplasm

Common bean (bean) is one of the most important legume crops, and mapping genes for yield and yield-related traits is essential for its improvement. However, yield is a complex trait that is typically controlled by many loci in crop genomes. The objective of this research was to identify regions in the bean genome associated with yield and a number of yield-related traits using a collection of 121 diverse bean genotypes with different yields. The beans were evaluated in replicated trials at two locations, over two years. Significant variation among genotypes was identified for all traits analyzed in the four environments. The collection was genotyped with the BARCBean6K_3 chip (5,398 SNPs), two yield/antiyield gene-based markers, and seven markers previously associated with resistance to common bacterial blight (CBB), including a Niemann-Pick polymorphism (NPP) gene-based marker. Over 90% of the single-nucleotide polymorphisms (SNPs) were polymorphic and separated the panel into two main groups of small-seeded and large-seeded beans, reflecting their Mesoamerican and Andean origins. Thirty-nine significant marker-trait associations (MTAs) were identified between 31 SNPs and 15 analyzed traits on all 11 bean chromosomes. Some of these MTAs confirmed genome regions previously associated with the yield and yield-related traits in bean, but a number of associations were not reported previously, especially those with derived traits. Over 600 candidate genes with different functional annotations were identified for the analyzed traits in the 200-Kb region centered on significant SNPs. Fourteen SNPs were identified within the gene model sequences, and five additional SNPs significantly associated with five different traits were located at less than 0.6 Kb from the candidate genes. The work confirmed associations between two yield/antiyield gene-based markers (AYD1m and AYD2m) on chromosome Pv09 with yield and identified their association with a number of yield-related traits, including seed weight. The results also confirmed the usefulness of the NPP marker in screening for CBB resistance. Since disease resistance and yield measurements are environmentally dependent and labor-intensive, the three gene-based markers (CBB- and two yield-related) and quantitative trait loci (QTL) that were validated in this work may be useful tools for simplifying and accelerating the selection of high-yielding and CBB-resistant bean cultivars.

Genetic control of pod morphological traits and pod edibility in a common bean RIL population

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.

Genetic erosion within the Fabada dry bean market class revealed by high-throughput genotyping

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.

High-density DArTSeq SNP markers revealed wide genetic diversity and structured population in common bean (Phaseolus vulgaris L.) germplasm in Ethiopia

INTRODUCTION

Common bean is one of the widely consumed food security crop in Africa, Asia, and South America. Understanding genetic diversity and population structure is crucial for designing breeding strategies.

MATERIALS

Two hundred and eighty-nine germplasm were recently collected from different regions of Ethiopia and introduced from CIAT to estimate genetic diversity and population structure using 11,480 DArTSeq SNP markers.

RESULTS

The overall mean genetic diversity and polymorphic information content (PIC) were 0.38 and 0.30, respectively, suggested the presence of adequate genetic diversity among the genotypes. Among the geographical regions, landraces collected from Oromia showed the highest diversity (0.39) and PIC (0.30). The highest genetic distance was observed between genotypes collected from SNNPR and CIAT (0.49). In addition, genotypes from CIAT were genetically more related to improved varieties than the landraces which could be due to sharing of parents in the improvement process. The analysis of molecular variance revealed that the largest proportion of variation was due to within the population both in geographical region (63.67%) and breeding status (61.3%) based classification. Model-based structure analysis delineated the 289 common bean genotypes into six hypothetical ancestoral populations.

CONCLUSIONS

The genotypes were not clustered based on geographical regions and they were not the main drivers for the differentiation. This indicated that selection of the parental lines should be based on systematic assessment of the diversity rather than geographical distance. This article provides new insights into the genetic diversity and population structure of common bean for association studies, designing effective collection and conservation for efficient utilization for the improvement of the crop.

Genome-wide association analysis reveals the optimal genomic regions for pod size in bean

The snap bean is the most commonly grown vegetable legume worldwide, and its pod size is both an important yield and appearance quality trait. However, the improvement of pod size in snap beans grown in China has been largely hindered by a lack of information on the specific genes that determine pod size. In this study, we identified 88 snap bean accessions and evaluated their pod size traits. Through a genome-wide association study (GWAS), 57 single nucleotide polymorphisms (SNPs) significantly associated with pod size were detected. Candidate gene analysis showed that cytochrome P450 family genes, WRKY, and MYB transcription factors were the predominant candidate genes for pod development, and eight of these 26 candidate genes showed relatively higher expression patterns in flowers and young pods. A significant pod length (PL) SNP and a single pod weight (SPW) SNP were successfully converted into kompetitive allele-specific polymerase chain reaction (KASP) markers and validated in the panel. These results enhance our understanding of the genetic basis of pod size, and also provide genetic resources for the molecular breeding of pod size in snap beans.

Genome-wide association study for the extractable phenolic profile and coat color of common bean seeds (Phaseolus vulgaris L.)

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.

Novel candidate loci for morpho-agronomic and seed quality traits detected by targeted genotyping-by-sequencing in common bean

Phaseolus vulgaris L., known as common bean, is one of the most important grain legumes cultivated around the world for its immature pods and dry seeds, which are rich in protein and micronutrients. Common bean offers a cheap food and protein sources to ameliorate food shortage and malnutrition around the world. However, the genetic basis of most important traits in common bean remains unknown. This study aimed at identifying QTL and candidate gene models underlying twenty-six agronomically important traits in common bean. For this, we assembled and phenotyped a diversity panel of 200 P. vulgaris genotypes in the greenhouse, comprising determinate bushy, determinate climbing and indeterminate climbing beans. The panel included dry beans and snap beans from different breeding programmes, elite lines and landraces from around the world with a major focus on accessions of African, European and South American origin. The panel was genotyped using a cost-conscious targeted genotyping-by-sequencing (GBS) platform to take advantage of highly polymorphic SNPs detected in previous studies and in diverse germplasm. The detected single nucleotide polymorphisms (SNPs) were applied in marker-trait analysis and revealed sixty-two quantitative trait loci (QTL) significantly associated with sixteen traits. Gene model identification via a similarity-based approach implicated major candidate gene models underlying the QTL associated with ten traits including, flowering, yield, seed quality, pod and seed characteristics. Our study revealed six QTL for pod shattering including three new QTL potentially useful for breeding. However, the panel was evaluated in a single greenhouse environment and the findings should be corroborated by evaluations across different field environments. Some of the detected QTL and a number of candidate gene models only elucidate the understanding of the genetic nature of these traits and provide the basis for further studies. Finally, the study showed the possibility of using a limited number of SNPs in performing marker-trait association in common bean by applying a highly scalable targeted GBS approach. This targeted GBS approach is a cost-efficient strategy for assessment of the genetic basis of complex traits and can enable geneticists and breeders to identify novel loci and targets for marker-assisted breeding more efficiently.

Innovation in the Breeding of Common Bean Through a Combined Approach of in vitro Regeneration and Machine Learning Algorithms

Common bean is considered a recalcitrant crop for in vitro regeneration and needs a repeatable and efficient in vitro regeneration protocol for its improvement through biotechnological approaches. In this study, the establishment of efficient and reproducible in vitro regeneration followed by predicting and optimizing through machine learning (ML) models, such as artificial neural network algorithms, was performed. Mature embryos of common bean were pretreated with 5, 10, and 20 mg/L benzylaminopurine (BAP) for 20 days followed by isolation of plumular apice for in vitro regeneration and cultured on a post-treatment medium containing 0.25, 0.50, 1.0, and 1.50 mg/L BAP for 8 weeks. Plumular apice explants pretreated with 20 mg/L BAP exerted a negative impact and resulted in minimum shoot regeneration frequency and shoot count, but produced longer shoots. All output variables (shoot regeneration frequency, shoot counts, and shoot length) increased significantly with the enhancement of BAP concentration in the post-treatment medium. Interaction of the pretreatment × post-treatment medium revealed the need for a specific combination for inducing a high shoot regeneration frequency. Higher shoot count and shoot length were achieved from the interaction of 5 mg/L BAP × 1.00 mg/L BAP followed by 10 mg/L BAP × 1.50 mg/L BAP and 20 mg/L BAP × 1.50 mg/L BAP. The evaluation of data through ML models revealed that R2 values ranged from 0.32 to 0.58 (regeneration), 0.01 to 0.22 (shoot counts), and 0.18 to 0.48 (shoot length). On the other hand, the mean squared error values ranged from 0.0596 to 0.0965 for shoot regeneration, 0.0327 to 0.0412 for shoot count, and 0.0258 to 0.0404 for shoot length from all ML models. Among the utilized models, the multilayer perceptron model provided a better prediction and optimization for all output variables, compared to other models. The achieved results can be employed for the prediction and optimization of plant tissue culture protocols used for biotechnological approaches in a breeding program of common beans.

A Core Set of Snap Bean Genotypes Established by Phenotyping a Large Panel Collected in Europe

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.

Dissecting the genetic control of seed coat color in a RIL population of common bean (Phaseolus vulgaris L.)

Three genes associated with the seed coat color in a TU/Musica RIL population were located on a genetic map, and two candidate genes proposed to control black seed coat in the TU genotype were characterized. Seed coat color is an important characteristic of common bean (Phaseolus vulgaris L.) associated with the marketability of dry bean cultivars, quality and nutritional characteristics of seed, as well as response to pathogens. In this study, the genetic control of seed coat color in a recombinant inbred line population (175 lines) obtained from the cross 'TU' × 'Musica' was investigated. Phenotypic segregation fitted 1:1 for white vs. nonwhite, and 3:1 for brown versus black, indicating the involvement of three independent genes, one controlling white color and two (with epistatic interaction) controlling black color. Using a genetic map built with 842 SNPs, the gene responsible for the white seed coat was mapped on the linkage group Pv07, in the position previously described for the P gene. For the black seed coat phenotype, two genes were mapped to the beginning of chromosomes Pv06 and Pv08, in the positions estimated for the V gene and the complex C locus, respectively, by classical studies. The involvement of these two genomic regions was verified through two crosses between three selected RILs exhibiting complementary and dominant inheritance, in which the TU alleles for both genes resulted in a black phenotype. Two genes involved in the anthocyanin biosynthesis pathway were proposed as candidate genes: Phvul.006G018800 encoding a flavonoid 3'5'hydroxylase and Phvul.008G038400 encoding MYB113 transcription factor. These findings add knowledge to the complex network of genes controlling seed coat color in common bean as well as providing genetic markers to be used in future genetic analysis or plant breeding.

Genetic diversity and inter-gene pool introgression of Mesoamerican Diversity Panel in common beans

Brazil is among the largest producers and consumers of common bean (Phaseolus vulgaris L.) and can be considered a secondary center of diversity for the species. The aim of this study was to estimate the genetic diversity, population structure, and relationships among 288 common bean accessions in an American Diversity Panel (ADP) genotyped with 4,042 high-quality single nucleotide polymorphisms (SNPs). The results showed inter-gene pool hybridization (hybrids) between the two main gene pools (i.e., Mesoamerican and Andean), based on principal component analysis (PCA), discriminant analysis of principal components (DAPC), and STRUCTURE analysis. The genetic diversity parameters showed that the Mesoamerican group has higher values of diversity and allelic richness in comparison with the Andean group. Considering the optimal clusters (K), clustering was performed according to the type of grain (i.e., market group), the institution of origin, the period of release, and agronomic traits. A new subset was selected and named the Mesoamerican Diversity Panel (MDP), with 205 Mesoamerican accessions. Analysis of molecular variance (AMOVA) showed low genetic variance between the two panels (i.e., ADP and MDP) with the highest percentage of the limited variance among accessions in each group. The ADP showed occurrence of high genetic differentiation between populations (i.e., Mesoamerican and Andean) and introgression between gene pools in hybrids based on a set of diagnostic SNPs. The MDP showed better linkage disequilibrium (LD) decay. The availability of genetic variation from inter-gene pool hybridizations presents a potential opportunity for breeders towards the development of superior common bean cultivars.

Genome-Wide Association Study and Genomic Prediction for Soybean Cyst Nematode Resistance in USDA Common Bean (Phaseolus vulgaris) Core Collection

Soybean cyst nematode (SCN, Heterodera glycines) has become the major yield-limiting biological factor in soybean production. Common bean is also a good host of SCN, and its production is challenged by this emerging pest in many regions such as the upper Midwest USA. The use of host genetic resistance has been the most effective and environmentally friendly method to manage SCN. The objectives of this study were to evaluate the SCN resistance in the USDA common bean core collection and conduct a genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers with SCN resistance. A total of 315 accessions of the USDA common bean core collection were evaluated for resistance to SCN HG Type 0 (race 6). The common bean core set was genotyped with the BARCBean6K_3 Infinium BeadChips, consisting of 4,654 SNPs. Results showed that 15 accessions were resistant to SCN with a Female Index (FI) at 4.8 to 9.4, and 62 accessions were moderately resistant (10 < FI < 30) to HG Type 0. The association study showed that 11 SNP markers, located on chromosomes Pv04, 07, 09, and 11, were strongly associated with resistance to HG Type 0. GWAS was also conducted for resistance to HG Type 2.5.7 and HG Type 1.2.3.5.6.7 based on the public dataset (N = 276), consisting of a diverse set of common bean accessions genotyped with the BARCBean6K_3 chip. Six SNPs associated with HG Type 2.5.7 resistance on Pv 01, 02, 03, and 07, and 12 SNPs with HG Type 1.2.3.5.6.7 resistance on Pv 01, 03, 06, 07, 09, 10, and 11 were detected. The accuracy of genomic prediction (GP) was 0.36 to 0.49 for resistance to the three SCN HG types, indicating that genomic selection (GS) of SCN resistance is feasible. This study provides basic information for developing SCN-resistant common bean cultivars, using the USDA core germ plasm accessions. The SNP markers can be used in molecular breeding in common beans through marker-assisted selection (MAS) and GS.

Gene Flow in Phaseolus Beans and Its Role as a Plausible Driver of Ecological Fitness and Expansion of Cultigens

The genus Phaseolus, native to the Americas, is composed of more than eighty wild species, five of which were domesticated in pre-Columbian times. Since the beginning of domestication events in this genus, ample opportunities for gene flow with wild relatives have existed. The present work reviews the extent of gene flow in the genus Phaseolus in primary and secondary areas of domestication with the aim of illustrating how this evolutionary force may have conditioned ecological fitness and the widespread adoption of cultigens. We focus on the biological bases of gene flow in the genus Phaseolus from a spatial and time perspective, the dynamics of wild-weedy-crop complexes in the common bean and the Lima bean, the two most important domesticated species of the genus, and the usefulness of genomic tools to detect inter and intraspecific introgression events. In this review we discuss the reproductive strategies of several Phaseolus species, the factors that may favor outcrossing rates and evidence suggesting that interspecific gene flow may increase ecological fitness of wild populations. We also show that wild-weedy-crop complexes generate genetic diversity over which farmers are able to select and expand their cultigens outside primary areas of domestication. Ultimately, we argue that more studies are needed on the reproductive biology of the genus Phaseolus since for most species breeding systems are largely unknown. We also argue that there is an urgent need to preserve wild-weedy-crop complexes and characterize the genetic diversity generated by them, in particular the genome-wide effects of introgressions and their value for breeding programs. Recent technological advances in genomics, coupled with agronomic characterizations, may make a large contribution.

GWAS of pod morphological and color characters in common bean

BACKGROUND

Common bean (Phaseolus vulgaris L.) is an important legume species which can be consumed as immature pods and dry seeds after re-hydration and cooking. Many genes and QTL, and epistatic interactions among them, condition pod morphological traits. However, not all them have been mapped or validated nor candidate genes proposed. We sought to investigate the genomic regions conditioning pod morphological and color characters through GWAS.

RESULTS

Single and multi-locus genome wide association analysis was used to investigate pod traits for a set of 301 bean lines of the Spanish Diversity Panel (SDP). The SDP was genotyped with 32,812 SNPs obtained from Genotyping by Sequencing. The panel was grown in two seasons and phenotypic data were recorded for 17 fresh pods traits grouped in four pod characters: pod length, pod cross-section, pod color, and number of seeds per pod. In all, 23 QTL for pod length, 6 for cross-section, 18 for pod color, 6 for number of seeds per pod and 9 associated to two or more pod characters were detected. Most QTL were located in the telomeric region of chromosomes Pv01, Pv02, Pv04, Pv08, Pv09 and Pv10. Eighteen detected QTL co-localized with 28 previously reported QTL. Twenty-one potential candidate genes involving developmental processes were detected underlying 11 QTL for pod morphological characters, four of them homologous to A. thaliana genes FIS2, SPL10, TTG2 and AML4 affecting silique size. Eight potential candidate genes involved in pigment synthesis, were found underlying five QTL for pod color.

CONCLUSIONS

GWAS for pod morphological and color characters in the bean Spanish Diversity Panel revealed 62 QTL, 18 co-localized with previously reported QTL, and 16 QTL were underlain by 25 candidate genes. Overall 44 new QTL identified and 18 existing QTL contribute to a better understanding of the complex inheritance of pod size and color traits in common bean and open the opportunity for future validation works.

Phenolic Content and Antioxidant Activity in Seeds of Common Bean (Phaseolus vulgaris L.)

Dry bean (Phaseolus vulgaris L.) is one of the most important pulses consumed in the world. Total phenolic content, total flavonoid content, total monomeric anthocyanin content and antioxidant capacity were determined, using ferric reducing antioxidant power and free radical scavenging activity, in 255 lines grown under the same environmental conditions. For all parameters analysed, there was a wide range of variability, with differences always above one order of magnitude. Phenolic compounds in beans with coloured coats were found to be more efficient antioxidants than those with completely white coats, and samples with more strongly coloured coats (red, cream, black, pink and brown) showed the highest antioxidant capacities. Based on the strong correlation detected between the variables, total phenolic content can be considered an appropriate indicator of antioxidant activity. The results provide a robust database for selecting those lines of greater functional and nutritional interest in terms of cultivation for direct consumption, for inclusions in food formulations or for use in future breeding programs.

Assessment of the Origin and Diversity of Croatian Common Bean Germplasm Using Phaseolin Type, SSR and SNP Markers and Morphological Traits

Landraces represent valuable genetic resources for breeding programmes to produce high-yielding varieties adapted to stressful environmental conditions. Although the common bean (Phaseolus vulgaris L.) is an economically important food legume for direct human consumption worldwide, common bean production in Croatia is based almost exclusively on landraces and there is no common bean breeding program. Information on phaseolin type and results of population structure and genetic diversity obtained by analysis of SSR and SNP markers, in combination with the morphological characterization of 174 accessions of 10 common bean landraces (morphotypes), enabled thorough classification of accessions. The accessions were classified into phaseolin type H1 (“S”) of Mesoamerican origin and phaseolin types H2 (“H” or “C”) and H3 (“T”) of Andean origin. By applying distance- and model-based clustering methods to SSR markers, the accessions were classified into two clusters at K = 2 separating the accessions according to the centres of origin, while at K = 3, the accessions of Andean origin were further classified into two clusters of accessions that differed in phaseolin type (H2 and H3). Using SNP markers, model-based analysis of population structure was performed, the results of which were consistent with those of SSR markers. In addition, 122 accessions were assigned to 14 newly formed true-type morphogenetic groups derived from three different domestication events: (1) Mesoamerican (H1A) (“Biser”, “Kukuruzar”, “Tetovac”, “Trešnjevac”), (2) Andean—indeterminate type (H2B1) (“Dan noć”, “Sivi”, “Puter”, ”Sivi prošarani”, “Trešnjevac”) and (3) Andean—determinate type (H3B2) (“Bijeli”, “Dan noć”, “Puter”, “Trešnjevac”, “Zelenčec”). The rest of the accessions could represent putative hybrids between morphogenetic groups. The differences between the true-type groups of accessions were further analysed based on nine quantitative traits, and the subsets of traits that best distinguish among centres of origin (A: Mesoamerican, B: Andean) and genetic groups (H1A, H2B1, H3B2) were proposed.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean

White mold (WM) is a devastating fungal disease affecting common bean (Phaseolus vulgaris L.). In this research, a genome-wide association study (GWAS) for WM resistance was conducted using 294 lines of the Spanish diversity panel. One single-locus method and six multi-locus methods were used in the GWAS. Response to this fungus showed a continuous distribution, and 28 lines were identified as potential resistance sources, including lines of Andean and Mesoamerican origin, as well as intermediate lines between the two gene pools. Twenty-two significant associations were identified, which were organized into 15 quantitative trait intervals (QTIs) located on chromosomes Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09. Seven of these QTIs were identified for the first time, whereas eight corresponded to chromosome regions previously identified in the WM resistance. In all, 468 genes were annotated in these regions, 61 of which were proposed potential candidate genes for WM resistance, based on their function related to the three main defense stages on the host: recognition (22), signal transduction (8), and defense response (31). Results obtained from this work will contribute to a better understanding of the complex quantitative resistance to WM in common bean and reveal information of significance for future breeding programs.

Genetic Diversity, Population Structure, and Andean Introgression in Brazilian Common Bean Cultivars after Half a Century of Genetic Breeding

Brazil is the largest consumer and third highest producer of common beans (Phaseolus vulgaris L.) worldwide. Since the 1980s, the commercial Carioca variety has been the most consumed in Brazil, followed by Black and Special beans. The present study evaluates genetic diversity and population structure of 185 Brazilian common bean cultivars using 2827 high-quality single-nucleotide polymorphisms (SNPs). The Andean allelic introgression in the Mesoamerican accessions was investigated, and a Carioca panel was tested using an association mapping approach. The results distinguish the Mesoamerican from the Andean accessions, with a prevalence of Mesoamerican accessions (94.6%). When considering the commercial classes, low levels of genetic differentiation were seen, and the Carioca group showed the lowest genetic diversity. However, gain in gene diversity and allelic richness was seen for the modern Carioca cultivars. A set of 1060 'diagnostic SNPs' that show alternative alleles between the pure Mesoamerican and Andean accessions were identified, which allowed the identification of Andean allelic introgression events and shows that there are putative introgression segments in regions enriched with resistance genes. Finally, genome-wide association studies revealed SNPs significantly associated with flowering time, pod maturation, and growth habit, showing that the Carioca Association Panel represents a powerful tool for crop improvements.

Characterization of extractable phenolic profile of common bean seeds (Phaseolus vulgaris L.) in a Spanish diversity panel

Phenolic compounds are important bioactive compounds in common bean (Phaseolus vulgaris L.). The aim of this work was the characterization of extractable phenolic profile (corresponding to 12 hydroxycinnamic acids and derivatives, 13 anthocyanins and 15 flavonols) in a bean diversity panel constituted by 220 lines, all grown under the same environmental conditions. Hydroxycinnamic derivatives were detected in all samples, while anthocyanins and flavonols were not detected in samples with completely white seed coats. In general, lines with black seeds showed higher contents of anthocyanins, followed by some red-seeded lines, while notable levels of flavonols were detected in market classes, including those with yellow, pink, and cream seed coats. However, a clear relationship between phenolic composition and seed phenotype could not be established, indicating the great influence of the genotype. This wide variability in the phenolic profiles analyzed is of particular interest for further breeding trials and the selection of varieties on the basis of this group of compounds.

Genome-wide association study of Striga resistance in early maturing white tropical maize inbred lines

BACKGROUND

Striga hermonthica (Benth.) parasitism militates against increased maize production and productivity in savannas of sub-Saharan Africa (SSA). Identification of Striga resistance genes is important in developing genotypes with durable resistance. So far, there is only one report on the existence of QTL for Striga resistance on chromosome 6 of maize. The objective of this study was to identify genomic regions significantly associated with grain yield and other agronomic traits under artificial Striga field infestation. A panel of 132 early-maturing maize inbreds were phenotyped for key agronomic traits under Striga-infested and Striga-free conditions. The inbred lines were also genotyped using 47,440 DArTseq markers from which 7224 markers were retained for population structure analysis and genome-wide association study (GWAS).

RESULTS

The inbred lines were grouped into two major clusters based on structure analysis as well as the neighbor-joining hierarchical clustering. A total of 24 SNPs significantly associated with grain yield, Striga damage at 8 and 10 weeks after planting (WAP), ears per plant and ear aspect under Striga infestation were detected. Under Striga-free conditions, 11 SNPs significantly associated with grain yield, number of ears per plant and ear aspect were identified. Three markers physically located close to the putative genes GRMZM2G164743 (bin 10.05), GRMZM2G060216 (bin 3.06) and GRMZM2G103085 (bin 5.07) were detected, linked to grain yield, Striga damage at 8 and 10 WAP and number of ears per plant under Striga infestation, explaining 9 to 42% of the phenotypic variance. Furthermore, the S9_154,978,426 locus on chromosome 9 was found at 2.61 Mb close to the ZmCCD1 gene known to be associated with the reduction of strigolactone production in the maize roots.

CONCLUSIONS

Presented in this study is the first report of the identification of significant loci on chromosomes 9 and 10 of maize that are closely linked to ZmCCD1 and amt5 genes, respectively and may be related to plant defense mechanisms against Striga parasitism. After validation, the identified loci could be targets for breeders for marker-assisted selection (MAS) to accelerate genetic enhancement of maize for Striga resistance in the tropics, particularly in SSA, where the parasitic weed is endemic.

Uncovering Phenotypic Diversity and DArTseq Marker Loci Associated with Antioxidant Activity in Common Bean

Antioxidants play an important role in animal and plant life owing to their involvement in complex metabolic and signaling mechanisms, hence uncovering the genetic basis associated with antioxidant activity is very important for the development of improved varieties. Here, a total of 182 common bean (Phaseolus vulgaris) landraces and six commercial cultivars collected from 19 provinces of Turkey were evaluated for seed antioxidant activity under four environments and two locations. Antioxidant activity was measured using ABTS radical scavenging capacity and mean antioxidant activity in common bean landraces was 20.03 µmol TE/g. Analysis of variance reflected that genotype by environment interaction was statistically non-significant and heritability analysis showed higher heritability of antioxidant activity. Variations in seed color were observed, and a higher antioxidant activity was present in seeds having colored seed as compared to those having white seeds. A negative correlation was found between white-colored seeds and antioxidant activity. A total of 7900 DArTseq markers were used to explore the population structure that grouped the studied germplasm into two sub-populations on the basis of their geographical origins and trolox equivalent antioxidant capacity contents. Mean linkage disequilibrium (LD) was 54%, and mean LD decay was 1.15 Mb. Mixed linear model i.e., the Q + K model demonstrated that four DArTseq markers had significant association (p < 0.01) for antioxidant activity. Three of these markers were present on chromosome Pv07, while the fourth marker was located on chromosome Pv03. Among the identified markers, DArT-3369938 marker showed maximum (14.61%) variation. A total of four putative candidate genes were predicted from sequences reflecting homology to identified DArTseq markers. This is a pioneering study involving the identification of association for antioxidant activity in common bean seeds. We envisage that this study will be very helpful for global common bean breeding community in order to develop cultivars with higher antioxidant activity.

Diversification and genetic structure of the western-to-eastern progression of European Phaseolus vulgaris L. germplasm

BACKGROUND

Common bean (Phaseolus vulgaris L.) is the most important food legume for direct human consumption around the world, as it represents a valuable source of components with nutritional and health benefits.

RESULTS

We conducted a study to define and explain the genetic relatedness and diversification level of common bean (Phaseolus vulgaris L.) germplasm from Portugal to Ukraine, along a western-to-eastern line of southern European countries, including Poland. This was based on the P. vulgaris genetic structure, and was designed to better describe its distribution and domestication pathways in Europe. Using the multi-crop passport descriptors that include geographic origin and different phaseolin types (corresponding to the Mesoamerican and Andean gene pools), 782 accessions were obtained from nine gene banks and 12 geographic origins. We selected 33 genome/ gene-related/ gene-pool-related nuclear simple sequence repeat markers that covered the genetic diversity across the P. vulgaris genome. The overall polymorphic information content was 0.800. Without specifying geographic origin, global structure cluster analysis generated 10 genetic clusters. Among the PvSHP1 markers, the most informative for gene pool assignment of the European P. vulgaris germplasm was PvSHP1-B. Results of AMOVA show that 89% of the molecular variability is shared within the 782 accessions, with 4% molecular variability among the different geographic origins along this western-to-eastern line of southern Europe (including Poland).

CONCLUSIONS

This study shows that the diversification line of the European P. vulgaris germplasm followed from the western areas of southern Europe (Portugal, Spain, Italy, Slovenia) to the more eastern areas of southern Europe. This progression defines three geographically separated subgroups, as the northern (Poland, Ukraine, Romania), southern (Albania, Bulgaria), and central (Bosnia and Herzegovina, Serbia, Hungary) areas of eastern Europe.

Genotyping-by-Sequencing Reveals Molecular Genetic Diversity in Italian Common Bean Landraces

The common bean (Phaseolus vulgaris L.) is one of the main legumes worldwide and represents a valuable source of nutrients. Independent domestication events in the Americas led to the formation of two cultivated genepools, namely Mesoamerican and Andean, to which European material has been brought back. In this study, Italian common bean landraces were analyzed for their genetic diversity and structure, using single nucleotide polymorphism (SNP) markers derived from genotyping-by-sequencing (GBS) technology. After filtering, 11,866 SNPs were obtained and 798 markers, pruned for linkage disequilibrium, were used for structure analysis. The most probable number of subpopulations (K) was two, consistent with the presence of the two genepools, identified through the phaseolin diagnostic marker. Some landraces were admixed, suggesting probable hybridization events between Mesoamerican and Andean material. When increasing the number of possible Ks, the Andean germplasm appeared to be structured in two or three subgroups. The subdivision within the Andean material was also observed in a principal coordinate analysis (PCoA) plot and a dendrogram based on genetic distances. The Mesoamerican landraces showed a higher level of genetic diversity compared to the Andean landraces. Calculation of the fixation index (FST) at individual SNPs between the Mesoamerican and Andean genepools and within the Andean genepool evidenced clusters of highly divergent loci in specific chromosomal regions. This work may help to preserve landraces of the common bean from genetic erosion, and could represent a starting point for the identification of interesting traits that determine plant adaptation.

Evidence for Strong Kinship Influence on the Extent of Linkage Disequilibrium in Cultivated Common Beans

Phaseolus vulgaris is an important grain legume for human consumption. Recently, association mapping studies have been performed for the species aiming to identify loci underlying quantitative variation of traits. It is now imperative to know whether the linkage disequilibrium (LD) reflects the true association between a marker and causative loci. The aim of this study was to estimate and analyze LD on a diversity panel of common beans using ordinary r² and r2 extensions which correct bias due to population structure (rS²), kinship (rV²), and both (rVS²). A total of 10,362 single nucleotide polymorphisms (SNPs) were identified by genotyping by sequencing (GBS), and polymorphisms were found to be widely distributed along the 11 chromosomes. In terms of r2, high values of LD (over 0.8) were identified between SNPs located at opposite chromosomal ends. Estimates for rV² were lower than those for rS². Results for rV² and rVS² were similar, suggesting that kinship may also include information on population structure. Over genetic distance, LD decayed to 0.1 at a distance of 1 Mb for rVS². Inter-chromosomal LD was also evidenced. This study showed that LD estimates decay dramatically according to the population structure, and especially the degree of kinship. Importantly, the LD estimates reported herein may influence our ability to perform association mapping studies on P. vulgaris.