Genome-Wide Association Study Identifies Candidate Loci Underlying Agronomic Traits in a Middle American Diversity Panel of Common Bean

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.

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.

Identification of loci associated with water use efficiency and symbiotic nitrogen fixation in soybean

Soybean (Glycine max) production is greatly affected by persistent and/or intermittent droughts in rainfed soybean-growing regions worldwide. Symbiotic N2 fixation (SNF) in soybean can also be significantly hampered even under moderate drought stress. The objective of this study was to identify genomic regions associated with shoot carbon isotope ratio (δ13C) as a surrogate measure for water use efficiency (WUE), nitrogen isotope ratio (δ15N) to assess relative SNF, N concentration ([N]), and carbon/nitrogen ratio (C/N). Genome-wide association mapping was performed with 105 genotypes and approximately 4 million single-nucleotide polymorphism markers derived from whole-genome resequencing information. A total of 11, 21, 22, and 22 genomic loci associated with δ13C, δ15N, [N], and C/N, respectively, were identified in two environments. Nine of these 76 loci were stable across environments, as they were detected in both environments. In addition to the 62 novel loci identified, 14 loci aligned with previously reported quantitative trait loci for different C and N traits related to drought, WUE, and N2 fixation in soybean. A total of 58 Glyma gene models encoding for different genes related to the four traits were identified in the vicinity of the genomic loci.

Genetic variation in a tepary bean (Phaseolus acutifolius A. Gray) diversity panel reveals loci associated with biotic stress resistance

Tepary bean (Phaseolus acutifolius A. Gray), indigenous to the arid climates of northern Mexico and the Southwest United States, diverged from common bean (Phaseolus vulgaris L.), approximately 2 million years ago and exhibits a wide range of resistance to biotic stressors. The tepary genome is highly syntenic to the common bean genome providing a foundation for discovery and breeding of agronomic traits between these two crop species. Although a limited number of adaptive traits from tepary bean have been introgressed into common bean, hybridization barriers between these two species required the development of bridging lines to alleviate this barrier. Thus, to fully utilize the extant tepary bean germplasm as both a crop and as a donor of adaptive traits, we developed a diversity panel of 422 cultivated, weedy, and wild tepary bean accessions which were then genotyped and phenotyped to enable population genetic analyses and genome-wide association studies for their response to a range of biotic stressors. Population structure analyses of the panel revealed eight subpopulations and the differentiation of botanical varieties within P. acutifolius. Genome-wide association studies revealed loci and candidate genes underlying biotic stress resistance including quantitative trait loci for resistance to weevils, common bacterial blight, Fusarium wilt, and bean common mosaic necrosis virus that can be harnessed not only for tepary bean but also common bean improvement.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

A robust SNP-haplotype assay for Bct gene region conferring resistance to beet curly top virus in common bean (Phaseolus vulgaris L.)

Beet curly top virus (BCTV), which is synonymous with curly top virus (CTV), causes significant yield loss in common bean (snap and dry beans) cultivars and several other important crops. Common bean cultivars have been found to be resistant to CTV, but screening for resistance is challenging due to the cyclical nature of epidemics and spotty feeding by the leafhopper that vectors the virus. We used an SNP dataset for the Snap Bean Association Panel (SnAP) agro-inoculated with CTV-Logan (CA/Logan) strain to locate the Bct gene region to a 1.7-Mb interval on chromosome Pv07 using genome-wide association study (GWAS) analysis. Recombinant lines from the SnAP were used to further narrow the Bct region to a 58.0-kb interval. A missense SNP (S07_2970381) in candidate gene Phvul.007G036300 Exonuclease V (EXO5) was identified as the most likely causal mutation, and it was the most significant SNP detected by GWAS in a dry bean population (DBP) naturally infected by the CTV-Worland (Wor) strain. Tm-shift assay markers developed for SNP S07_2970381 and two linked SNPs, S07_2970276 and S07_2966197, were useful for tracking different origins of the Bct EXO5 candidate gene resistance to CTV in common bean. The three SNPs identified four haplotypes, with haplotype 3-1 (Haplo3-1) of Middle American origin associated with the highest levels of CTV resistance. This SNP-haplotype assay will enable breeders to track resistance sources and to develop cultivars with better CTV resistance.

Combination of meta-analysis of QTL and GWAS to uncover the genetic architecture of seed yield and seed yield components in common bean

Increasing seed yield in common bean could help to improve food security and reduce malnutrition globally due to the high nutritional quality of this crop. However, the complex genetic architecture and prevalent genotype by environment interactions for seed yield makes increasing genetic gains challenging. The aim of this study was to identify the most consistent genomic regions related with seed yield components and phenology reported in the last 20 years in common bean. A meta-analysis of quantitative trait locus (QTL) for seed yield components and phenology (MQTL-YC) was performed for 394 QTL reported in 21 independent studies under sufficient water and drought conditions. In total, 58 MQTL-YC over different genetic backgrounds and environments were identified, reducing threefold on average the confidence interval (CI) compared with the CI for the initial QTL. Furthermore, 40 MQTL-YC identified were co-located with 210 SNP peak positions reported via genome-wide association (GWAS), guiding the identification of candidate genes. Comparative genomics among these MQTL-YC with MQTL-YC reported in soybean and pea allowed the identification of 14 orthologous MQTL-YC shared across species. The integration of MQTL-YC, GWAS, and comparative genomics used in this study is useful to uncover and refine the most consistent genomic regions related with seed yield components for their use in plant breeding.

Genome-wide association mapping for component traits of drought tolerance in dry beans (Phaseolus vulgaris L.)

Understanding the genetic basis of traits of economic importance under drought stressed and well-watered conditions is important in enhancing genetic gains in dry beans (Phaseolus vulgaris L.). This research aims to: (i) identify markers associated with agronomic and physiological traits for drought tolerance and (ii) identify drought-related putative candidate genes within the mapped genomic regions. An andean and middle-american diversity panel (AMDP) comprising of 185 genotypes was screened in the field under drought stressed and well-watered conditions for two successive seasons. Agronomic and physiological traits, viz., days to 50% flowering (DFW), plant height (PH), days to physiological maturity (DPM), grain yield (GYD), 100-seed weight (SW), leaf temperature (LT), leaf chlorophyll content (LCC) and stomatal conductance (SC) were phenotyped. Principal component and association analysis were conducted using the filtered 9370 Diversity Arrays Technology sequencing (DArTseq) markers. The mean PH, GYD, SW, DPM, LCC and SC of the panel was reduced by 12.1, 29.6, 10.3, 12.6, 28.5 and 62.0%, respectively under drought stressed conditions. Population structure analysis revealed two sub-populations, which corresponded to the andean and middle-american gene pools. Markers explained 0.08-0.10, 0.22-0.23, 0.29-0.32, 0.43-0.44, 0.65-0.66 and 0.69-0.70 of the total phenotypic variability (R2) for SC, LT, PH, GYD, SW and DFW, respectively under drought stressed conditions. For well-watered conditions, R2 varied from 0.08 (LT) to 0.70 (DPM). Overall, 68 significant (p < 10-03) marker-trait associations (MTAs) and 22 putative candidate genes were identified across drought stressed and well-watered conditions. Most of the identified genes had known biological functions related to regulating the response to drought stress. The findings provide new insights into the genetic architecture of drought stress tolerance in common bean. The findings also provide potential candidate SNPs and putative genes that can be utilized in gene discovery and marker-assisted breeding for drought tolerance after validation.

Meta-QTL Analysis for Yield Components in Common Bean (Phaseolus vulgaris L.)

Common bean is one of the most important legumes produced and consumed worldwide because it is a highly valuable food for the human diet. However, its production is mainly carried out by small farmers, who obtain average grain yields below the potential yield of the species. In this sense, numerous mapping studies have been conducted to identify quantitative trait loci (QTL) associated with yield components in common bean. Meta-QTL (MQTL) analysis is a useful approach to combine data sets and for creating consensus positions for the QTL detected in independent studies. Consequently, the objective of this study was to perform a MQTL analysis to identify the most reliable and stable genomic regions associated with yield-related traits of common bean. A total of 667 QTL associated with yield-related traits reported in 21 different studies were collected. A total of 42 MQTL associated with yield-related traits were identified, in which the average confidence interval (CI) of the MQTL was 3.41 times lower than the CIs of the original QTL. Most of the MQTL (28) identified in this study contain QTL associated with yield and phenological traits; therefore, these MQTL can be useful in common bean breeding programs. Finally, a total of 18 candidate genes were identified and associated with grain yield within these MQTL, with functions related to ubiquitin ligase complex, response to auxin, and translation elongation factor activity.

Modified screening method of middle american dry bean genotypes reveals new genomic regions on Pv10 associated with anthracnose resistance

Anthracnose, caused by the fungal pathogen Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.-Scrib., is one of the most devastating diseases in dry bean (Phaseolus vulgaris L.) with seed yield losses up to 100%. Most anthracnose resistance genes thus far identified behave in a dominant manner and were identified by seedling screening. The Middle American Diversity Panel (MDP; n=266) was screened with a modified greenhouse screening method to evaluate the response to anthracnose race 73. Thirty MDP genotypes exhibited resistance to the race of which 16 genotypes were not known to contain anthracnose resistance genes to race 73. GWAS with ~93,000 SNP markers identified four genomic regions, two each on Pv01 and Pv10, associated race 73 resistance. A likelihood-ratio-based R2 analysis indicated the peak four SNP markers are responsible for 26% of the observed phenotypic variation, where one SNP, S10_072250, explains 23% of the total variation. SNP S10_072250 is associated with a new region of anthracnose resistance and is in an intron of a ZPR1-like gene. Further greenhouse testing of the 16 resistant lines without previously known resistance to race 73 revealed various levels of resistance under various levels of disease pressure. Disease resistance was further characterized in the field using four representative genotypes. GTS-900 and Remington exhibited field resistance while Merlot and Maverick were susceptible. Field testing with two different fungicide regimes revealed the resistant genotypes had no significant disease differences. The results suggest resistance to anthracnose may differ at various growth stages and that breeders have been selecting for major genes at early seedling stages while ignoring the effect of alternative genes that may be active at later stages. The newly identified resistant lines may be related to Age Related Resistance (ARR) and could be exploited as parental sources of anthracnose resistance in addition to already known major genes. The physical localization of the multiple regions of resistance confirms the presence of two clusters of disease resistance genes on Pv01 and identifies two new regions of anthracnose resistance on Pv10 possibly associated with ARR. Future research should look at the mode of inheritance of this resistance and its effect when combined with other anthracnose resistance loci.

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.

Characterization of fungal pathogens and germplasm screening for disease resistance in the main production area of the common bean in Argentina

The common bean (Phaseolus vulgaris L.) is the most important grain legume in the human diet, mainly in Africa and Latin America. Argentina is one of the five major producers of the common bean in the world, and the main cultivation areas are concentrated in the northwestern provinces of this country. Crop production of the common bean is often affected by biotic factors like some endemic fungal diseases, which exert a major economic impact on the region. The most important fungal diseases affecting the common bean in Argentina are white mold caused by Sclerotinia sclerotiorum, angular leaf spot caused by Pseudocercospora griseola, web blight and root rot caused by Rhizoctonia solani, which can cause production losses of up to 100% in the region. At the present, the most effective strategy for controlling these diseases is the use of genetic resistance. In this sense, population study and characterization of fungal pathogens are essential for developing cultivars with durable resistance. In this review we report diversity studies carried out on these three fungal pathogens affecting the common bean in northwestern Argentina, analyzing more than 200 isolates by means of molecular, morphological and pathogenic approaches. Also, the screening of physiological resistance in several common bean commercial lines and wild native germplasm is reviewed. This review contributes to the development of sustainable management strategies and cultural practices in bean production aimed to minimize yield losses due to fungal diseases in the common bean.

Genome-Wide Association Mapping of Virulence Genes in Wheat Karnal Bunt Fungus Tilletia indica Using Double Digest Restriction-Site Associated DNA-Genotyping by Sequencing Approach

Tilletia indica is a quarantine fungal pathogen that poses a serious biosecurity threat to wheat-exporting countries. Acquiring genetic data for the pathogenicity characters of T. indica is still a challenge for wheat breeders and geneticists. In the current study, double digest restriction-site associated-DNA genotyping by sequencing was carried out for 39 T. indica isolates collected from different locations in India. The generated libraries upon sequencing were with 3,346,759 raw reads on average, and 151 x 2 nucleotides read length. The obtained bases per read ranged from 87 Mb in Ti 25 to 1,708 Mb in Ti 39, with 505 Mb on average per read. Trait association mapping was performed using 41,473 SNPs, infection phenotyping data, population structure, and Kinship matrix, to find single nucleotide polymorphisms (SNPs) linked to virulence genes. Population structure analysis divided the T. indica population in India into three subpopulations with genetic mixing in each subpopulation. However, the division was not in accordance with the degree of virulence. Trait association mapping revealed the presence of 13 SNPs associated with virulence. Using sequences analysis tools, one gene (g4132) near a significant SNP was predicted to be an effector, and its relative expression was assessed and found upregulated upon infection.

Genome-Wide Association Study Reveals That PvGUX1_1 Is Associated with Pod Stringlessness in Snap Bean (Phaseolus vulgaris L.)

Simple Summary

Using 138 snap bean accessions as plant materials, we investigated their suture strings across two years. With the goal of identifying the gene(s) responsible for the formation of suture strings, we conducted a genome-wide association study. A strong association signal was found in a 266.19 kb region on Chr02. Within the region, 23 candidate genes were identified. Importantly, the sequence and gene expression of PvGUX1_1 differed significantly between sutured pods and non-sutured pods. In addition, PvGUX1_1 was also a domesticated locus that diverged from PvGUX1_2 during an early stage. The results obtained in this study can provide important information for the improvement of pod quality in snap beans.

Abstract

Suture strings are a particularly important pod trait that determine the quality and texture of snap beans (Phaseolus vulgaris L.). The St locus on chromosome 2 has been described as a major locus associated with suture strings. However, the gene and genetic basis underlying this locus remain unknown. Here, we investigated the suture strings of 138 snap bean accessions across two years. A total of 3.66 million single-nucleotide polymorphisms (SNPs) were obtained by deep resequencing. Based on these SNPs, we identified a strong association signal on Chr02 and a promising candidate gene, PvGUX1_1. Further analysis revealed that the 2 bp deletion in the exon of PvGUX1_1 was significantly associated with stringlessness. Comparative mapping indicated that PvGUX1_1 was a domesticated locus and diverged from PvGUX1_2 during an early stage. Our study provides important insights into the genetic mechanism of suture string formation and useful information for snap bean improvement.

Doing Genetic and Genomic Biology Using the Legume Information System and Associated Resources

In this chapter, we introduce the main components of the Legume Information System ( https://legumeinfo.org ) and several associated resources. Additionally, we provide an example of their use by exploring a biological question: is there a common molecular basis, across legume species, that underlies the photoperiod-mediated transition from vegetative to reproductive development, that is, days to flowering? The Legume Information System (LIS) holds genetic and genomic data for a large number of crop and model legumes and provides a set of online bioinformatic tools designed to help biologists address questions and tasks related to legume biology. Such tasks include identifying the molecular basis of agronomic traits; identifying orthologs/syntelogs for known genes; determining gene expression patterns; accessing genomic datasets; identifying markers for breeding work; and identifying genetic similarities and differences among selected accessions. LIS integrates with other legume-focused informatics resources such as SoyBase ( https://soybase.org ), PeanutBase ( https://peanutbase.org ), and projects of the Legume Federation ( https://legumefederation.org ).

Genome-Wide Association Study Identifies Genomic Regions for Important Morpho-Agronomic Traits in Mesoamerican Common Bean

The population growth trend in recent decades has resulted in continuing efforts to guarantee food security in which leguminous plants, such as the common bean (Phaseolus vulgaris L.), play a particularly important role as they are relatively cheap and have high nutritional value. To meet this demand for food, the main target for genetic improvement programs is to increase productivity, which is a complex quantitative trait influenced by many component traits. This research aims to identify Quantitative Trait Nucleotides (QTNs) associated with productivity and its components using multi-locus genome-wide association studies. Ten morpho-agronomic traits [plant height (PH), first pod insertion height (FPIH), number of nodules (NN), pod length (PL), total number of pods per plant (NPP), number of locules per pod (LP), number of seeds per pod (SP), total seed weight per plant (TSW), 100-seed weight (W100), and grain yield (YLD)] were evaluated in four environments for 178 Mesoamerican common bean domesticated accessions belonging to the Brazilian Diversity Panel. In order to identify stable QTNs, only those identified by multiple methods (mrMLM, FASTmrMLM, pLARmEB, and ISIS EM-BLASSO) or in multiple environments were selected. Among the identified QTNs, 64 were detected at least thrice by different methods or in different environments, and 39 showed significant phenotypic differences between their corresponding alleles. The alleles that positively increased the corresponding traits, except PH (for which lower values are desired), were considered favorable alleles. The most influenced trait by the accumulation of favorable alleles was PH, showing a 51.7% reduction, while NN, TSW, YLD, FPIH, and NPP increased between 18 and 34%. Identifying QTNs in several environments (four environments and overall adjusted mean) and by multiple methods reinforces the reliability of the associations obtained and the importance of conducting these studies in multiple environments. Using these QTNs through molecular techniques for genetic improvement, such as marker-assisted selection or genomic selection, can be a strategy to increase common bean production.

Genome-Wide Association Study of Seed Folate Content in Common Bean

Plant-derived folates (Vitamin B₉) are essential components of the human diet. They provide one-carbon units that are required for the synthesis of nucleic acids and proteins, and folate deficiency is associated with numerous adverse health conditions. The development of high-folate cultivars of common bean (Phaseolus vulgaris L.) and other staple crops is an important tool to combat folate deficiency. A population of 96 P. vulgaris accessions, representing major North American market classes, was grown in 2 years in Ontario, Canada. The population was genotyped for 5,361 molecular markers with an Illumina Infinium platform. Total folate was extracted from mature seeds using the tri-enzyme extraction method and quantified based on a microbiological assay with Lactobacillus rhamnosus. Significant genetic diversity for folate content was observed among the population in both years of study, and folate content had a range 113-222 μg per 100 g of seeds. Quantitative trait loci (QTL) for seed folate content were identified based on a genome-wide association study (GWAS). Six QTL were identified on Chr. 4, 6, 8, and 11, with three in each year of field trials. Both QTL on Chr. 11 occurred in genomic regions that were syntenic to seed folate QTL detected in previous work with P. vulgaris, Z. mays, and O. sativa. Candidate genes were identified for these QTL that might be targets for the development of molecular markers for selecting P. vulgaris cultivars with improved seed folate content. This work reports the largest survey of genetic diversity for seed folate content in P. vulgaris and identified several genotypes, including SCN4, Bat 93, OAC Redstar, and Pompadour 1014, that would be useful for breeding beans with higher than average folate levels.

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.

Genetic dissection of days to flowering via genome-wide association studies in Turkish common bean germplasm

Common bean is a nutrient-dense legume crop serving as a source of food for millions of people. Characterization of unexplored common bean germplasm to unlock the phenotypic and genetic variations is still needed to explore the breeding potential of this crop. The current study aimed to dissect the genetic basis having association for days to flowering (DF). A total of 188 common bean accessions collected from 19 provinces of Turkey were used as plant material under five environments and two locations. Analysis of variance (ANOVA) revealed that genotypes and genotype by environment interaction have significant effects on DF. A total of 10 most stable accessions were evaluated from stability analysis. Overall maximum (75) and minimum (54) DF were observed for Hakkari-51 and Mus-46 accessions, respectively. The implemented constellation plot divided studied germplasm according to their DF and growth habit. A total of 7900 DArTseq markers were used for association analysis. Mixed linear model using the Q + K Model resulted a total of 18 DArTseq markers from five environments. DArT-8668385 marker identified in Bolu during 2016 was also associated with DF in Sivas during 2017. Combined data of five years resulted a total of four markers (DArT-22346534, DArT-3369768, DArT-3374613, and DArT-3370801) having significant association ( p  <  0.01 ) for DF. DArT-22346534 present on Pv 08 accounted a maximum of 9.89% variation to the studied trait. A total of four putative candidate genes were predicted from sequences reflecting homology to identified four DArTseq markers. We envisage that exploitation of identified DArTseq markers will hopefully beneficial for the development of new common bean varieties having better adaptation ability to changing climatic conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01029-8.

Genome Wide Association Mapping of Root Traits in the Andean Genepool of Common Bean (Phaseolus vulgaris L.) Grown With and Without Aluminum Toxicity

Common bean is one of the most important grain legumes for human diets but is produced on marginal lands with unfavorable soil conditions; among which Aluminum (Al) toxicity is a serious and widespread problem. Under low pH, stable forms of Al dissolve into the soil solution and as phytotoxic ions inhibit the growth and function of roots through injury to the root apex. This results in a smaller root system that detrimentally effects yield. The goal of this study was to evaluate 227 genotypes from an Andean diversity panel (ADP) of common bean and determine the level of Al toxicity tolerance and candidate genes for this abiotic stress tolerance through root trait analysis and marker association studies. Plants were grown as seedlings in hydroponic tanks at a pH of 4.5 with a treatment of high Al concentration (50 μM) compared to a control (0 μM). The roots were harvested and scanned to determine average root diameter, root volume, root surface area, number of root links, number of root tips, and total root length. Percent reduction or increase was calculated for each trait by comparing treatments. Genome wide association study (GWAS) was conducted by testing phenotypic data against single nucleotide polymorphism (SNP) marker genotyping data for the panel. Principal components and a kinship matrix were included in the mixed linear model to correct for population structure. Analyses of variance indicated the presence of significant difference between genotypes. The heritability of traits ranged from 0.67 to 0.92 in Al-treated and reached similar values in non-treated plants. GWAS revealed significant associations between root traits and genetic markers on chromosomes Pv01, Pv04, Pv05, Pv06, and Pv11 with some SNPs contributing to more than one trait. Candidate genes near these loci were analyzed to explain the detected association and included an Al activated malate transporter gene and a multidrug and toxic compound extrusion gene. This study showed that polygenic inheritance was critical to aluminum toxicity tolerance in common beans roots. Candidate genes found suggested that exudation of malate and citrate as organic acids would be important for Al tolerance. Possible cross-talk between mechanisms of aluminum tolerance and resistance to other abiotic stresses are discussed.

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.

Sources of Resistance to Common Bacterial Blight and Charcoal Rot Disease for the Production of Mesoamerican Common Beans in the Southern United States

The gene pool of Mesoamerican common beans (Phaseolus vulgaris L.) includes genotypes in the small-to-medium-size seeded dry beans, as well as some snap beans from hotter environments adapted to the Southeastern United States. However, the warm and humid climate of the Southeastern United States is conducive to diseases such as Common Bacterial Blight (CBB) and Charcoal Rot (CR). The pathogens for these two diseases can survive long periods in infested soil or on seeds and are difficult to control through pesticides. Hence, field-level resistance would be the best management strategy for these diseases. The goals of this study were (1) to evaluate field-level resistance from the various commercial classes and subgroups represented in the Mesoamerican gene pool as sources for breeding beans for the region and (2) to evaluate genome-wide marker × trait associations (GWAS) using genetic markers for the genotypes. A total of 300 genotypes from the Mesoamerican Diversity Panel (MDP) were evaluated for CBB and CR in field experiments for three years. CBB resistance was also tested with a field isolate in controlled greenhouse conditions. The analysis of variance revealed the presence of variability in the MDP for the evaluated traits. We also identified adapted common bean genotypes that could be used directly in Southeastern production or that could be good parents in breeding programs for CBB and CR resistance. The GWAS detected 14 significant Single-Nucleotide Polymorphism (SNP) markers associated with CBB resistance distributed on five chromosomes, namely Pv02, Pv04, Pv08, Pv10, and Pv11, but no loci for resistance to CR. A total of 89 candidate genes were identified in close vicinity (±100 kb) to the significant CBB markers, some of which could be directly or indirectly involved in plant defense to diseases. These results provide a basis to further understand the complex inheritance of CBB resistance in Mesoamerican common beans and show that this biotic stress is unrelated to CR resistance, which was evident during a drought period. Genotypes with good yield potential for the Southeastern U.S. growing conditions were found with resistant to infection by the two diseases, as well as adaptation to the hot and humid conditions punctuated by droughts found in this region.

Towards the Development, Maintenance, and Standardized Phenotypic Characterization of Single-Seed-Descent Genetic Resources for Common Bean

The optimal use of legume genetic resources represents a key prerequisite for coping with current agriculture-related societal challenges, including conservation of agrobiodiversity, agricultural sustainability, food security, and human health. Among legumes, the common bean (Phaseolus vulgaris) is the most economically important for human consumption, and its evolutionary trajectories as a species have been crucial to determining the structure and level of its present and available genetic diversity. Genomic advances are considerably enhancing the characterization and assessment of important genetic variants. For this purpose, the development and availability of, and access to, well-described and efficiently managed genetic resource collections that comprise pure lines derived by single-seed-descent cycles will be paramount for the use of the reservoir of common bean variability and for the advanced breeding of legume crops. This is one of the main aims of the new and challenging European project INCREASE, which is the implementation of Intelligent Collections with appropriate standardized protocols that must be characterized, maintained, and made available, along with the related data, to users such as breeders and researchers. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Characterizing common bean seeds for seed trait descriptors Basic Protocol 2: Bean seed imaging Basic Protocol 3: Characterizing bean lines for plant trait descriptors specific for common bean Primary Seed Increase.

Physiological and genetic characterization of heat stress effects in a common bean RIL population

Heat stress is a major abiotic stress factor reducing crop productivity and climate change models predict increasing temperatures in many production regions. Common bean (Phaseolus vulgaris L.) is an important crop for food security in the tropics and heat stress is expected to cause increasing yield losses. To study physiological responses and to characterize the genetics of heat stress tolerance, we evaluated the recombinant inbred line (RIL) population IJR (Indeterminate Jamaica Red) x AFR298 of the Andean gene pool. Heat stress (HS) conditions in the field affected many traits across the reproductive phase. High nighttime temperatures appeared to have larger effects than maximum daytime temperatures. Yield was reduced compared to non-stress conditions by 37% and 26% in 2016 and 2017 seasons, respectively. The image analysis tool HYRBEAN was developed to evaluate pollen viability (PolVia). A significant reduction of PolVia was observed in HS and higher viability was correlated with yield only under stress conditions. In susceptible lines the reproductive phase was extended and defects in the initiation of seed, seed fill and seed formation were identified reducing grain quality. Higher yields under HS were correlated with early flowering, high pollen viability and effective seed filling. Quantitative trait loci (QTL) analysis revealed a QTL for both pod harvest index and PolVia on chromosome Pv05, for which the more heat tolerant parent IJR contributed the positive allele. Also, on chromosome Pv08 a QTL from IJR improved PolVia and the yield component pods per plant. HS affected several traits during the whole reproductive development, from floral induction to grain quality traits, indicating a general heat perception affecting many reproductive processes. Identification of tolerant germplasm, indicator traits for heat tolerance and molecular tools will help to breed heat tolerant varieties to face future climate change effects.

Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean

The bean fly (Ophiomyia spp) is a key insect pest causing significant crop damage and yield loss in common bean (Phaseolus vulgaris L., 2n = 2x = 22). Development and deployment of agronomic superior and bean fly resistant common bean varieties aredependent on genetic variation and the identification of genes and genomic regions controlling economic traits. This study's objective was to determine the population structure of a diverse panel of common bean genotypes and deduce associations between bean fly resistance and agronomic traits based on single nucleotide polymorphism (SNP) markers. Ninety-nine common bean genotypes were phenotyped in two seasons at two locations and genotyped with 16 565 SNP markers. The genotypes exhibited significant variation for bean fly damage severity (BDS), plant mortality rate (PMR), and pupa count (PC). Likewise, the genotypes showed significant variation for agro-morphological traits such as days to flowering (DTF), days to maturity (DTM), number of pods per plant (NPP), number of seeds per pod (NSP), and grain yield (GYD). The genotypes were delineated into two populations, which were based on the Andean and Mesoamerican gene pools. The genotypes exhibited a minimum membership coefficient of 0.60 to their respective populations. Eighty-three significant (P<0.01) markers were identified with an average linkage disequilibrium of 0.20 at 12Mb across the 11 chromosomes. Three markers were identified, each having pleiotropic effects on two traits: M100049197 (BDS and NPP), M3379537 (DTF and PC), and M13122571 (NPP and GYD). The identified markers are useful for marker-assisted selection in the breeding program to develop common bean genotypes with resistance to bean fly damage.

Ancient relaxation of an obligate short-day requirement in common bean through loss of CONSTANS-like gene function

Common bean (Phaseolus vulgaris L.) is a major global food staple and source of dietary protein that was domesticated independently in Mexico and Andean South America. Its subsequent development as a crop of importance worldwide has been enabled by genetic relaxation of the strict short-day requirement typical of wild forms, but the genetic basis for this change is not well understood. Recently, a loss of photoperiod sensitivity was shown to result from mutations in the phytochrome photoreceptor gene Ppd/PHYA3 that arose independently within the two major domesticated lineages. Here, we define a second major photoperiod sensitivity locus, at which recessive alleles associate with deleterious mutations affecting the CONSTANS-like gene COL2. A wider survey of sequence variation in over 800 diverse lines, including wild, landrace, and domesticated accessions, show that distinct col2 haplotypes are associated with early flowering in Andean and Mesoamerican germplasm. The relative frequencies and distributions of COL2 and PHYA3 haplotypes imply that photoperiod adaptation developed in two phases within each gene pool: an initial reduction in sensitivity through impairment of COL2 function and subsequent complete loss through PHYA3. Gene expression analyses indicate that COL2 functions downstream of PHYA3 to repress expression of FT genes and may function in parallel with PvE1, the bean ortholog of a key legume-specific flowering repressor. Collectively, these results define the molecular basis for a key phenological adaptation, reveal a striking convergence in the naturally replicated evolution of this major crop, and further emphasize the wider evolutionary lability of CONSTANS effects on flowering time control.

Pod shattering in grain legumes: emerging genetic and environment-related patterns

A reduction in pod shattering is one of the main components of grain legume domestication. Despite this, many domesticated legumes suffer serious yield losses due to shattering, particularly under arid conditions. Mutations related to pod shattering modify the twisting force of pod walls or the structural strength of the dehiscence zone in pod sutures. At a molecular level, a growing body of evidence indicates that these changes are controlled by a relatively small number of key genes that have been selected in parallel across grain legume species, supporting partial molecular convergence. Legume homologs of Arabidopsis thaliana silique shattering genes play only minor roles in legume pod shattering. Most domesticated grain legume species contain multiple shattering-resistance genes, with mutants of each gene typically showing only partial shattering resistance. Hence, crosses between varieties with different genes lead to transgressive segregation of shattering alleles, producing plants with either enhanced shattering resistance or atavistic susceptibility to the trait. The frequency of these resistance pod-shattering alleles is often positively correlated with environmental aridity. The continued development of pod-shattering-related functional information will be vital for breeding crops that are suited to the increasingly arid conditions expected in the coming decades.

NAC Candidate Gene Marker for bgm-1 and Interaction With QTL for Resistance to Bean Golden Yellow Mosaic Virus in Common Bean

Genetic resistance is the primary means for control of Bean golden yellow mosaic virus (BGYMV) in common bean (Phaseolus vulgaris L.). Breeding for resistance is difficult because of sporadic and uneven infection across field nurseries. We sought to facilitate breeding for BGYMV resistance by improving marker-assisted selection (MAS) for the recessive bgm-1 gene and identifying and developing MAS for quantitative trait loci (QTL) conditioning resistance. Genetic linkage mapping in two recombinant inbred line populations and genome-wide association study (GWAS) in a large breeding population and two diversity panels revealed a candidate gene for bgm-1 and three QTL BGY4.1, BGY7.1, and BGY8.1 on independent chromosomes. A mutation (5 bp deletion) in a NAC (No Apical Meristem) domain transcriptional regulator superfamily protein gene Phvul.003G027100 on chromosome Pv03 corresponded with the recessive bgm-1 resistance allele. The five bp deletion in exon 2 starting at 20 bp (Pv03: 2,601,582) is expected to cause a stop codon at codon 23 (Pv03: 2,601,625), disrupting further translation of the gene. A T m -shift assay marker named PvNAC1 was developed to track bgm-1. PvNAC1 corresponded with bgm-1 across ∼1,000 lines which trace bgm-1 back to a single landrace "Garrapato" from Mexico. BGY8.1 has no effect on its own but exhibited a major effect when combined with bgm-1. BGY4.1 and BGY7.1 acted additively, and they enhanced the level of resistance when combined with bgm-1. T m -shift assay markers were generated for MAS of the QTL, but their effectiveness requires further validation.

Using Breeding Populations With a Dual Purpose: Cultivar Development and Gene Mapping-A Case Study Using Resistance to Common Bacterial Blight in Dry Bean (Phaseolus vulgaris L.)

Dry bean (Phaseolus vulgaris L.) is an important worldwide legume crop with low to moderate levels of resistance to common bacterial blight (CBB) caused by Xanthomonas axonopodis pv. phaseoli. A total of 852 genotypes (cultivars, preliminary and advanced breeding lines) from the North Dakota State University dry bean breeding program were tested for their effectiveness as populations for genome-wide association studies (GWAS) to identify genomic regions associated with resistance to CBB, to exploit the associated markers for marker-assisted breeding (MAB), and to identify candidate genes. The genotypes were evaluated in a growth chamber for disease resistance at both the unifoliate and trifoliate stages. At the unifoliate stage, 35% of genotypes were resistant, while 25% of genotypes were resistant at the trifoliate stage. Libraries generated from each genotype were sequenced using the Illumina platform. After filtering for sequence quality, read depth, and minor allele frequency, 41,998 single-nucleotide polymorphisms (SNPs) and 30,285 SNPs were used in GWAS for the Middle American and Andean gene pools, respectively. One region near the distal end of Pv10 near the SAP6 molecular marker from the Andean gene pool explained 26.7-36.4% of the resistance variation. Three to seven regions from the Middle American gene pool contributed to 25.8-27.7% of the resistance, with the most significant peak also near the SAP6 marker. Six of the eight total regions associated with CBB resistance are likely the physical locations of quantitative trait loci identified from previous genetic studies. The two new locations associated with CBB resistance are located at Pv10:22.91-23.36 and Pv11:52.4. A lipoxgenase-1 ortholog on Pv10 emerged as a candidate gene for CBB resistance. The state of one SNP on Pv07 was associated with susceptibility. Its subsequent use in MAB would reduce the current number of lines in preliminary and advanced field yield trial by up to 14% and eliminate only susceptible genotypes. These results provide a foundational SNP data set, improve our understanding of CBB resistance in dry bean, and impact resource allocation within breeding programs as breeding populations may be used for dual purposes: cultivar development as well as genetic studies.

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.

Field Screen and Genotyping of Phaseolus vulgaris against Two Begomoviruses in Georgia, USA

Simple Summary

Snap bean (Phaseolus vulgaris) production and quality have been negatively impacted by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. However, there is no information available in terms of resistance to these two viruses in commercial cultivars/genotypes. Hence, commercially available snap bean varieties/genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years (with a few exceptions) were screened in two field seasons of 2018 and 2019. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. As a result of this screening, we identified twenty Phaseolus genotypes with high-to-moderate levels of resistance and twenty-one genotypes with high levels of susceptibility. While there were differences among the Phaseolus spp. in severity of viral symptoms, suggesting differential susceptibility to viruses (CuLCrV and SiGMFV) and potential field resistance, the resistance mechanism is yet to be characterized. However, based on the greenhouse evaluation with two genotypes-each (susceptible vs. resistant) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistant might potentially be against the virus complex than against the whiteflies. Adult whitefly counts differed among the Phaseolus spp. in both the years, indicating variability in host preference. We further sequenced 82 genotypes (80 snap bean and two Lima bean) to unravel the variations within the genomes. Genome sequencing followed by bioinformatic analyses revealed a considerable number of sequence variants, single nucleotide polymorphisms (SNPs), and insertions and deletions (InDels) in the genomes. Considering the variations in disease response and the underlying variations in the sequenced genomes, it can be speculated that some of the phenotypic variations (against CuLCrV and SiGMFV) could be due to a high level of genomic variation in the host. Future genome-wide association studies with the identified genomic variants may shed some light on this.

Abstract

The production and quality of Phaseolusvulgaris (snap bean) have been negatively impacted by leaf crumple disease caused by two whitefly-transmitted begomoviruses: cucurbit leaf crumple virus (CuLCrV) and sida golden mosaic Florida virus (SiGMFV), which often appear as a mixed infection in Georgia. Host resistance is the most economical management strategy against whitefly-transmitted viruses. Currently, information is not available with respect to resistance to these two viruses in commercial cultivars. In two field seasons (2018 and 2019), we screened Phaseolus spp. genotypes (n = 84 in 2018; n = 80 in 2019; most of the genotypes were common in both years with a few exceptions) for resistance against CuLCrV and/or SiGMFV. We also included two commonly grown Lima bean (Phaseolus lunatus) varieties in our field screening. Twenty Phaseolus spp. genotypes with high to moderate-levels of resistance (disease severity ranging from 5%–50%) to CuLCrV and/or SiGMFV were identified. Twenty-one Phaseolus spp. genotypes were found to be highly susceptible with a disease severity of ≥66%. Furthermore, based on the greenhouse evaluation with two genotypes-each (two susceptible and two resistant; identified in field screen) exposed to viruliferous whiteflies infected with CuLCrV and SiGMFV, we observed that the susceptible genotypes accumulated higher copy numbers of both viruses and displayed severe crumple severity compared to the resistant genotypes, indicating that resistance might potentially be against the virus complex rather than against the whiteflies. Adult whitefly counts differed significantly among Phaseolus genotypes in both years. The whole genome of these Phaseolus spp. [snap bean (n = 82); Lima bean (n = 2)] genotypes was sequenced and genetic variability among them was identified. Over 900 giga-base (Gb) of filtered data were generated and >88% of the resulting data were mapped to the reference genome, and SNP and Indel variants in Phaseolus spp. genotypes were obtained. A total of 645,729 SNPs and 68,713 Indels, including 30,169 insertions and 38,543 deletions, were identified, which were distributed in 11 chromosomes with chromosome 02 harboring the maximum number of variants. This phenotypic and genotypic information will be helpful in genome-wide association studies that will aid in identifying the genetic basis of resistance to these begomoviruses in Phaseolus spp.

Genome-Wide Association Mapping of bc-1 and bc-u Reveals Candidate Genes and New Adjustments to the Host-Pathogen Interaction for Resistance to Bean Common Mosaic Necrosis Virus in Common Bean

Bean common mosaic necrosis virus (BCMNV) is a major disease in common bean (Phaseolus vulgaris L.). Host plant resistance is the primary disease control. We sought to identify candidate genes to better understand the host-pathogen interaction and develop tools for marker-assisted selection (MAS). A genome-wide association study (GWAS) approach using 182 lines from a race Durango Diversity Panel (DDP) challenged by BCMNV isolates NL-8 [Pathogroup (PG)-III] and NL-3 (PG-VI), and genotyped with 1.26 million single-nucleotide polymorphisms (SNPs), revealed significant peak regions on chromosomes Pv03 and Pv05, which correspond to bc-1 and bc-u resistance gene loci, respectively. Three candidate genes were identified for NL-3 and NL-8 resistance. Side-by-side receptor-like protein kinases (RLKs), Phvul.003G038700 and Phvul.003G038800 were candidate genes for bc-1. These RLKs were orthologous to linked RLKs associated with virus resistance in soybean (Glycine max). A basic Leucine Zipper (bZIP) transcription factor protein is the candidate gene for bc-u. bZIP protein gene Phvul.005G124100 carries a unique non-synonymous mutation at codon 14 in the first exon (Pv05: 36,114,516 bases), resulting in a premature termination codon that causes a nonfunctional protein. SNP markers for bc-1 and bc-u and new markers for I and bc-3 genes were used to genotype the resistance genes underpinning BCMNV phenotypes in the DDP, host group (HG) differentials, and segregating F₃ families. Results revealed major adjustments to the current host-pathogen interaction model: (i) there is only one resistance allele bc-1 for the Bc-1 locus, and differential expression of the allele is based on presence vs. absence of bc-u; (ii) bc-1 exhibits dominance and incomplete dominance; (iii) bc-1 alone confers resistance to NL-8; (iv) bc-u was absent from HGs 2, 4, 5, and 7 necessitating a new gene symbol bc-u ^d to reflect this change; (v) bc-u ^d alone delays susceptible symptoms, and when combined with bc-1 enhanced resistance to NL-3; and (vi) bc-u ^d is on Pv05, not Pv03 as previously thought. These candidate genes, markers, and adjustments to the host-pathogen interaction will facilitate breeding for resistance to BCMNV and related Bean common mosaic virus (BCMV) in common bean.

Genome-Wide Association Studies of Mineral Content in Common Bean

Micronutrient malnutrition is one of the main public health problems in many parts of the world. This problem raises the attention of all valuable sources of micronutrients for the human diet, such as common bean (Phaseolus vulgaris L.). In this research, a panel of 174 accessions representing Croatian common bean landraces was phenotyped for seed content of eight nutrients (N, P, K, Ca, Mg, Fe, Zn, and Mn), and genotyped using 6,311 high-quality DArTseq-derived SNP markers. A genome-wide association study (GWAS) was then performed to identify new genetic sources for improving seed mineral content. Twenty-two quantitative trait nucleotides (QTN) associated with seed nitrogen content were discovered on chromosomes Pv01, Pv02, Pv03, Pv05, Pv07, Pv08, and Pv10. Five QTNs were associated with seed phosphorus content, four on chromosome Pv07, and one on Pv08. A single significant QTN was found for seed calcium content on chromosome Pv09 and for seed magnesium content on Pv08. Finally, two QTNs associated with seed zinc content were identified on Pv06 while no QTNs were found to be associated with seed potassium, iron, or manganese content. Our results demonstrate the utility of GWAS for understanding the genetic architecture of seed nutritional traits in common bean and have utility for future enrichment of seed with macro- and micronutrients through genomics-assisted breeding.

Toward the introgression of PvPdh1 for increased resistance to pod shattering in common bean

A common bean shattering-resistance allele of PvPdh1 reduces pod twists during dehiscence, shows dominance that varies by phenotyping method, is part of a selective sweep, and can be introgressed using CAPS markers. Some varieties of common bean (Phaseolus vulgaris L.) suffer from pod shattering, which can severely reduce yields, especially in arid conditions. The PvPdh1 locus on chromosome Pv03 has recently been described as a major locus controlling pod shattering in common bean and could be used to mitigate pod shattering in the future. Despite this, the role of a possible second locus on chromosome Pv08 remains unclear and patterns of dominance and epistasis between alleles of these genes have not been resolved. This information will be vital for efficient selection to decrease pod shattering. Further, the genetic diversity around the PvPdh1 gene has not yet been thoroughly explored, and there are not yet genetic screens that can be used to evaluate pod shattering in segregating populations. Here, we have developed a recombinant inbred population to determine the roles of genes implicated in pod shattering and evaluate the patterns of dominance among the relevant alleles. Our results suggest that a PvPdh1 allele reduces pod valve twisting, and its dominance varies by phenotyping method. This allele is the only genetic variant that provides environmentally stable and widespread resistance to pod shattering in Middle American common beans grown for grain. Further analyses identified a selective sweep around PvPdh1 with greater nucleotide diversity in individuals with the ancestral, shattering-susceptible allele. Finally, we developed simple, effective CAPS markers to facilitate the introgression of PvPdh1 into new varieties of common bean. These genetic resources will be critical for improving the aridity resilience of a major global staple.

Genetic Factors Associated With Nodulation and Nitrogen Derived From Atmosphere in a Middle American Common Bean Panel

Among grain legume crops, common beans (Phaseolus vulgaris L.) are considered to have poor biological nitrogen (N₂) fixation (BNF) capabilities although variation in N₂ fixing capabilities exists within the species. The availability of genetic panel varying in BNF capacity and a large-scale single nucleotide polymorphism (SNP) data set for common bean provided an opportunity to discover genetic factors associated with N₂ fixation among genotypes in the Middle American gene pool. Using nodulation and percentage of N₂-derived from atmosphere (%NDFA) data collected from field trials, at least 11 genotypes with higher levels of BNF capacity were identified. Genome-wide association studies (GWASs) detected both major and minor effects that control these traits. A major nodulation interval at Pv06:28.0-28.27 Mbp was discovered. In this interval, the peak SNP was located within a small GTPase that positively regulates cellular polarity and growth of root hair tips. Located 20 kb upstream of this peak SNP is an auxin-responsive factor AUX/indole acetic auxin (IAA)-related gene involved in auxin transportation during root nodulation. For %NDFA, nitrate (NO₃ ^-) transporters, NRT1:2 and NRT1.7 (Pv02:8.64), squamosa promoter binding transcriptome factor (Pv08:28.42), and multi-antimicrobial extrusion protein (MATE) efflux family protein (Pv06:10.91) were identified as candidate genes. Three additional QTLs were identified on chromosomes Pv03:5.24, Pv09:25.89, and Pv11: 32.89 Mbp. These key candidate genes from both traits were integrated with previous results on N₂ fixation to describe a BNF pathway.

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.

GCViT: a method for interactive, genome-wide visualization of resequencing and SNP array data

Background

Large genotyping datasets have become commonplace due to efficient, cheap methods for SNP identification. Typical genotyping datasets may have thousands to millions of data points per accession, across tens to thousands of accessions. There is a need for tools to help rapidly explore such datasets, to assess characteristics such as overall differences between accessions and regional anomalies across the genome.

Results

We present GCViT (Genotype Comparison Visualization Tool), for visualizing and exploring large genotyping datasets. GCViT can be used to identify introgressions, conserved or divergent genomic regions, pedigrees, and other features for more detailed exploration. The program can be used online or as a local instance for whole genome visualization of resequencing or SNP array data. The program performs comparisons of variants among user-selected accessions to identify allele differences and similarities between accessions and a user-selected reference, providing visualizations through histogram, heatmap, or haplotype views. The resulting analyses and images can be exported in various formats.

Conclusions

GCViT provides methods for interactively visualizing SNP data on a whole genome scale, and can produce publication-ready figures. It can be used in online or local installations. GCViT enables users to confirm or identify genomics regions of interest associated with particular traits.

GCViT is freely available at https://github.com/LegumeFederation/gcvit. The 1.0 version described here is available at 10.5281/zenodo.4008713.

Genetic mapping for agronomic traits in a MAGIC population of common bean (Phaseolus vulgaris L.) under drought conditions

BACKGROUND

Common bean is an important staple crop in the tropics of Africa, Asia and the Americas. Particularly smallholder farmers rely on bean as a source for calories, protein and micronutrients. Drought is a major production constraint for common bean, a situation that will be aggravated with current climate change scenarios. In this context, new tools designed to understand the genetic basis governing the phenotypic responses to abiotic stress are required to improve transfer of desirable traits into cultivated beans.

RESULTS

A multiparent advanced generation intercross (MAGIC) population of common bean was generated from eight Mesoamerican breeding lines representing the phenotypic and genotypic diversity of the CIAT Mesoamerican breeding program. This population was assessed under drought conditions in two field trials for yield, 100 seed weight, iron and zinc accumulation, phenology and pod harvest index. Transgressive segregation was observed for most of these traits. Yield was positively correlated with yield components and pod harvest index (PHI), and negative correlations were found with phenology traits and micromineral contents. Founder haplotypes in the population were identified using Genotyping by Sequencing (GBS). No major population structure was observed in the population. Whole Genome Sequencing (WGS) data from the founder lines was used to impute genotyping data for GWAS. Genetic mapping was carried out with two methods, using association mapping with GWAS, and linkage mapping with haplotype-based interval screening. Thirteen high confidence QTL were identified using both methods and several QTL hotspots were found controlling multiple traits. A major QTL hotspot located on chromosome Pv01 for phenology traits and yield was identified. Further hotspots affecting several traits were observed on chromosomes Pv03 and Pv08. A major QTL for seed Fe content was contributed by MIB778, the founder line with highest micromineral accumulation. Based on imputed WGS data, candidate genes are reported for the identified major QTL, and sequence changes were identified that could cause the phenotypic variation.

CONCLUSIONS

This work demonstrates the importance of this common bean MAGIC population for genetic mapping of agronomic traits, to identify trait associations for molecular breeding tool design and as a new genetic resource for the bean research community.

European landrace diversity for common bean biofortification: a genome-wide association study

Mineral deficiencies represent a global challenge that needs to be urgently addressed. An adequate intake of iron and zinc results in a balanced diet that reduces chances of impairment of many metabolic processes that can lead to clinical consequences. In plants, bioavailability of such nutrients is reduced by presence of compounds such as phytic acid, that can chelate minerals and reduce their absorption. Biofortification of common bean (Phaseolus vulgaris L.) represents an important strategy to reduce mineral deficiencies, especially in areas of the world where this crop plays a key role in the diet. In this study, a panel of diversity encompassing 192 homozygous genotypes, was screened for iron, zinc and phytate seed content. Results indicate a broad variation of these traits and allowed the identification of accessions reasonably carrying favourable trait combinations. A significant association between zinc seed content and some molecular SNP markers co-located on the common bean Pv01 chromosome was detected by means of genome-wide association analysis. The gene Phvul001G233500, encoding for an E3 ubiquitin-protein ligase, is proposed to explain detected associations. This result represents a preliminary evidence that can foster future research aiming at understanding the genetic mechanisms behind zinc accumulation in beans.

Genome-Wide Association Studies Detect Multiple QTLs for Productivity in Mesoamerican Diversity Panel of Common Bean Under Drought Stress

Drought stress is an important abiotic factor limiting common bean yield, with great impact on the production worldwide. Understanding the genetic basis regulating beans' yield and seed weight (SW) is a fundamental prerequisite for the development of superior cultivars. The main objectives of this work were to conduct genome-wide marker discovery by genotyping a Mesoamerican panel of common bean germplasm, containing cultivated and landrace accessions of broad origin, followed by the identification of genomic regions associated with productivity under two water regimes using different genome-wide association study (GWAS) approaches. A total of 11,870 markers were genotyped for the 339 genotypes, of which 3,213 were SilicoDArT and 8,657 SNPs derived from DArT and CaptureSeq. The estimated linkage disequilibrium extension, corrected for structure and relatedness (r ² _sv ), was 98.63 and 124.18 kb for landraces and breeding lines, respectively. Germplasm was structured into landraces and lines/cultivars. We carried out GWASs for 100-SW and yield in field environments with and without water stress for 3 consecutive years, using single-, segment-, and gene-based models. Higher number of associations at high stringency was identified for the SW trait under irrigation, totaling ∼185 QTLs for both single- and segment-based, whereas gene-based GWASs showed ∼220 genomic regions containing ∼650 genes. For SW under drought, 18 QTLs were identified for single- and segment-based and 35 genes by gene-based GWASs. For yield, under irrigation, 25 associations were identified, whereas under drought the total was 10 using both approaches. In addition to the consistent associations detected across experiments, these GWAS approaches provided important complementary QTL information (∼221 QTLs; 650 genes; r ² from 0.01% to 32%). Several QTLs were mined within or near candidate genes playing significant role in productivity, providing better understanding of the genetic mechanisms underlying these traits and making available molecular tools to be used in marker-assisted breeding. The findings also allowed the identification of genetic material (germplasm) with better yield performance under drought, promising to a common bean breeding program. Finally, the availability of this highly diverse Mesoamerican panel is of great scientific value for the analysis of any relevant traits in common bean.

North-Western Himalayan Common Beans: Population Structure and Mapping of Quantitative Anthracnose Resistance Through Genome Wide Association Study

Common bean (Phaseolus vulgaris L.) is an important legume crop of north-western (NW) Himalayan region and the major disease that causes catastrophic loss to the crop is anthracnose, which is caused by Colletotrichum lindemuthianum. The pathogen is highly diverse and most of the commercial cultivars are susceptible to different races prevalent in the region. The lack of information on the genomic regions associated with anthracnose resistance in NW Himalayan common bean population prompted us to dissect Quantitative Resistance Loci (QRLs) against major anthracnose races. In this study, 188 common bean landraces collected from NW region were screened against five important anthracnose races and 113 bean genotypes showed resistance to one or multiple races. Genotyping by sequencing (GBS) was performed on a panel of 192 bean lines (4 controls plus 188 Indian beans) and 22,589 SNPs were obtained that are evenly distributed. Population structure analysis of 192 bean genotypes categorized 188 Indian beans into two major clusters representing Andean and Mesoamerican gene pools with obvious admixtures. Many QRLs associated with anthracnose resistance to Indian C. lindemuthianum virulences (race 3, 87, and 503) are located at Pv04 within the gene models that encode typical resistance gene signatures. The QRLs associated with race 73 are located on Pv08 and overlaps with Co-4 anthracnose resistance gene. A SNP located at distal end of Pv11 in a gene model Phvul.011G202300 which encodes a LRR with a typical NB-ARC domain showed association with race 73 resistance. Common bean genomic regions located at Pv03, Pv09, and Pv11 showed association with resistance to anthracnose race 2047. The present study showed presence of many novel bean genomic regions associated with anthracnose resistance. The presence of Co-4 and Co-2 genes in our material is encouraging for breeding durable anthracnose resistant cultivars for the region.

Genetic Diversity, Nitrogen Fixation, and Water Use Efficiency in a Panel of Honduran Common Bean (Phaseolus vulgaris L.) Landraces and Modern Genotypes

Common bean (Phaseolus vulgaris L.) provides critical nutrition and a livelihood for millions of smallholder farmers worldwide. Beans engage in symbiotic nitrogen fixation (SNF) with Rhizobia. Honduran hillside farmers farm marginal land and utilize few production inputs; therefore, bean varieties with high SNF capacity and environmental resiliency would be of benefit to them. We explored the diversity for SNF, agronomic traits, and water use efficiency (WUE) among 70 Honduran landrace, participatory bred (PPB), and conventionally bred bean varieties (HON panel) and 6 North American check varieties in 3 low-N field trials in Ontario, Canada and Honduras. Genetic diversity was measured with a 6K single nucleotide polymorphism (SNP) array, and phenotyping for agronomic, SNF, and WUE traits was carried out. STRUCTURE analysis revealed two subpopulations with admixture between the subpopulations. Nucleotide diversity was greater in the landraces than the PPB varieties across the genome, and multiple genomic regions were identified where population genetic differentiation between the landraces and PPB varieties was evident. Significant differences were found between varieties and breeding categories for agronomic traits, SNF, and WUE. Landraces had above average SNF capacity, conventional varieties showed higher yields, and PPB varieties performed well for WUE. Varieties with the best SNF capacity could be used in further participatory breeding efforts.

Marker association study of yield attributing traits in common bean (Phaseolus vulgaris L.)

Common bean is gaining acceptance as one of the most valuable major food consumed worldwide owing to innumerable nutritional and therapeutic benefits. Comparatively less productivity in underdeveloped countries encouraged us to proceed for QTL mining of yield traits in common bean. Heretofore, multiple yield associated markers have been detected all over the world; even so, the present work is looked on as the first report on identification of novel/new potent markers by exploiting the germplasm of Northern India. A panel of one hundred and thirty five genotypes was used for morphological studies and based on preliminary molecular evaluation; a set of ninety six diverse common bean genotypes (core set) was selected for association analysis. Molecular data generated by a total of ninety eight microsatellite markers (53 genomic and 45 genic SSRs) revealed high estimation of polymorphism among the genotypes that were observed to be divided into two major sub-populations and varying levels of admixtures based on population structure analyses. By employing both MLM and GLM analysis approaches, we identified 46 and 16 significant marker-trait associations (p ≤ 0.005) respectively, few of which have already been reported and hence validate our results. PVBR213 marker was found to be strongly associated with days to bud initiation trait when analyzed with both the approaches. Phenotypic variation of identified significant markers ranged from 3.1% to 32.7% where PVBR87, PVBR213, X96999 and X57022 explain more than 30% of phenotypic variation for 100 seed weight, days to bud initiation, pods per plant and pod length traits respectively. These findings introduce highly informative markers to aid marker-assisted selection program in common bean for high yield performance along with good agronomic merit.

QTL analysis of cooking time and quality traits in dry bean (Phaseolus vulgaris L.)

Three robust QTL for dry bean cooking time shortened cooking time 11-26 min and co-localized with QTL for increased cooked seed protein concentration. Cooking time is a major factor associated with consumer preference of dry beans (Phaseolus vulgaris L.). The genetic control of cooking time was investigated with a quantitative trait loci (QTL) study on a recombinant inbred line (RIL) population developed from TZ-27 (slow cooking) and TZ-37 (fast cooking). The RIL population of 146 lines was grown on research farms over 2 years in Arusha and Morogoro, Tanzania. Arusha is an important mid-altitude bean-growing region, with moderate temperatures and reliable rainfall, whereas the low altitude and high temperatures in Morogoro make it unfavorable for bean production. The population exhibited large variation for cooking time with a range of 22-98 min. On average, beans grown in Arusha cooked 15 min faster than those grown in Morogoro. A linkage map developed with 1951 SNP markers was used for QTL analysis. Ten QTL were identified for cooking time, three of which were found in multiple environments. RILs with all three QTL (CT3.1, CT6.1, and CT11.2) cooked on average 11 min faster in Arusha and 26 min faster in Morogoro than RILs with none. Seed attributes were related to cooking time such that seeds with greater seed mass and less seed coat percentage cooked faster. Cooked seed protein concentration ranged from 17.8 to 30.8% across the years and locations. All three of the most robust cooking time QTL co-localized with QTL for protein concentration, and TZ-37 always contributed faster cooking time and increased protein concentration.

Sources of Resistance to Fusarium solani and Associated Genomic Regions in Common Bean Diversity Panels

Common bean (Phaseolus vulgaris L.) production worldwide is hampered by Fusarium root rot (FRR), which is caused by Fusarium solani. Screening for FRR resistance on a large scale is notoriously difficult and often yields inconsistent results due to variability within the environment and pathogen biology. A greenhouse screening assay was developed incorporating multiple isolates of F. solani to improve assay reproducibility. The Andean (ADP; n = 270) and Middle American (MDP; n = 280) Diversity Panels were screened in the greenhouse to identify genetic factors associated with FRR resistance. Forty-seven MDP and 34 ADP lines from multiple market classes were identified as resistant to FRR. Greenhouse phenotyping repeatability was confirmed via five control lines. Genome-wide association mapping using ∼200k SNPs was performed on standard phenotyping score 1–9, as well as binary and polynomial transformation of score data. Sixteen and seven significant genomic regions were identified for ADP and MDP, respectively, using all three classes of phenotypic data. Most candidate genes were associated with plant immune/defense mechanisms. For the ADP population, ortholog of glucan synthase-like enzyme, senescence-associated genes, and NAC domain protein, associated with peak genomic region Pv08:0.04–0.18 Mbp, were the most significant candidate genes. For the MDP population, the peak SNPs Pv07:15.29 Mbp and Pv01:51 Mbp mapped within gene models associated with ethylene response factor 1 and MAC/Perforin domain-containing gene respectively. The research provides a basis for bean improvement through the use of resistant genotypes and genomic regions for more durable root rot resistance.

Determining the Genetic Control of Common Bean Early-Growth Rate Using Unmanned Aerial Vehicles

Vigorous early-season growth rate allows crops to compete more effectively against weeds and to conserve soil moisture in arid areas. These traits are of increasing economic importance due to changing consumer demand, reduced labor availability, and climate-change-related increasing global aridity. Many crop species, including common bean, show genetic variation in growth rate, between varieties. Despite this, the genetic basis of early-season growth has not been well-resolved in the species, in part due to historic phenotyping challenges. Using a range of UAV- and ground-based methods, we evaluated the early-season growth vigor of two populations. These growth data were used to find genetic regions associated with several growth parameters. Our results suggest that early-season growth rate is the result of complex interactions between several genetic and environmental factors. They also highlight the need for high-precision phenotyping provided by UAVs. The quantitative trait loci (QTLs) identified in this study are the first in common bean to be identified remotely using UAV technology. These will be useful for developing crop varieties that compete with weeds and use water more effectively. Ultimately, this will improve crop productivity in the face of changing climatic conditions and will mitigate the need for water and resource-intensive forms of weed control.

Genetic Associations in Four Decades of Multienvironment Trials Reveal Agronomic Trait Evolution in Common Bean

Multienvironment trials (METs) are widely used to assess the performance of promising crop germplasm. Though seldom designed to elucidate genetic mechanisms, MET data sets are often much larger than could be duplicated for genetic research and, given proper interpretation, may offer valuable insights into the genetics of adaptation across time and space. The Cooperative Dry Bean Nursery (CDBN) is a MET for common bean (Phaseolus vulgaris) grown for > 70 years in the United States and Canada, consisting of 20-50 entries each year at 10-20 locations. The CDBN provides a rich source of phenotypic data across entries, years, and locations that is amenable to genetic analysis. To study stable genetic effects segregating in this MET, we conducted genome-wide association studies (GWAS) using best linear unbiased predictions derived across years and locations for 21 CDBN phenotypes and genotypic data (1.2 million SNPs) for 327 CDBN genotypes. The value of this approach was confirmed by the discovery of three candidate genes and genomic regions previously identified in balanced GWAS. Multivariate adaptive shrinkage (mash) analysis, which increased our power to detect significant correlated effects, found significant effects for all phenotypes. Mash found two large genomic regions with effects on multiple phenotypes, supporting a hypothesis of pleiotropic or linked effects that were likely selected on in pursuit of a crop ideotype. Overall, our results demonstrate that statistical genomics approaches can be used on MET phenotypic data to discover significant genetic effects and to define genomic regions associated with crop improvement.

Exploration of the Yield Potential of Mesoamerican Wild Common Beans From Contrasting Eco-Geographic Regions by Nested Recombinant Inbred Populations

Genetic analyses and utilization of wild genetic variation for crop improvement in common bean (Phaseolus vulgaris L.) have been hampered by yield evaluation difficulties, identification of advantageous variation, and linkage drag. The lack of adaptation to cultivation conditions and the existence of highly structured populations make association mapping of diversity panels not optimal. Joint linkage mapping of nested populations avoids the later constraint, while populations crossed with a common domesticated parent allow the evaluation of wild variation within a more adapted background. Three domesticated by wild backcrossed-inbred-line populations (BC₁S₄) were developed using three wild accessions representing the full range of rainfall of the Mesoamerican wild bean distribution crossed to the elite drought tolerant domesticated parent SEA 5. These populations were evaluated under field conditions in three environments, two fully irrigated trials in two seasons and a simulated terminal drought in the second season. The goal was to test if these populations responded differently to drought stress and contained progenies with higher yield than SEA 5, not only under drought but also under water-watered conditions. Results revealed that the two populations derived from wild parents of the lower rainfall regions produced lines with higher yield compared to the domesticated parent in the three environments, i.e., both in the drought-stressed environment and in the well-watered treatments. Several progeny lines produced yields, which on average over the three environments were 20% higher than the SEA 5 yield. Twenty QTLs for yield were identified in 13 unique regions on eight of the 11 chromosomes of common bean. Five of these regions showed at least one wild allele that increased yield over the domesticated parent. The variation explained by these QTLs ranged from 0.6 to 5.4% of the total variation and the additive effects ranged from -164 to 277 kg ha^-1, with evidence suggesting allelic series for some QTLs. Our results underscore the potential of wild variation, especially from drought-stressed regions, for bean crop improvement as well the identification of regions for efficient marker-assisted introgression.

Pod indehiscence is a domestication and aridity resilience trait in common bean

Plant domestication has strongly modified crop morphology and development. Nevertheless, many crops continue to display atavistic characteristics that were advantageous to their wild ancestors but are deleterious under cultivation, such as pod dehiscence (PD). Here, we provide the first comprehensive assessment of the inheritance of PD in the common bean (Phaseolus vulgaris), a major domesticated grain legume. Using three methods to evaluate the PD phenotype, we identified multiple, unlinked genetic regions controlling PD in a biparental population and two diversity panels. Subsequently, we assessed patterns of orthology among these loci and those controlling the trait in other species. Our results show that different genes were selected in each domestication and ecogeographic race. A chromosome Pv03 dirigent-like gene, involved in lignin biosynthesis, showed a base-pair substitution that is associated with decreased PD. This haplotype may underlie the expansion of Mesoamerican domesticates into northern Mexico, where arid conditions promote PD. The rise in frequency of the decreased-PD haplotype may be a consequence of the markedly different fitness landscape imposed by domestication. Environmental dependency and genetic redundancy can explain the maintenance of atavistic traits under domestication.