Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement.

Common bean (Phaseolus vulgaris L.) breeding programs aim to improve both agronomic and seed characteristics traits. However, the genetic architecture of the many traits that affect common bean production are not completely understood. Genome-wide association studies (GWAS) provide an experimental approach to identify genomic regions where important candidate genes are located. A panel of 280 modern bean genotypes from race Mesoamerica, referred to as the Middle American Diversity Panel (MDP), were grown in four US locations, and a GWAS using >150,000 single-nucleotide polymorphisms (SNPs) (minor allele frequency [MAF] ³ 5%) was conducted for six agronomic traits. The degree of inter-and intrachromosomal linkage disequilibrium (LD) was estimated after accounting for population structure and relatedness. The LD varied between chromosomes for the entire MDP and among race Mesoamerica and DurangoJalisco genotypes within the panel. The LD patterns reflected the breeding history of common bean. Genome-wide association studies led to the discovery of new and known genomic regions affecting the agronomic traits at the entire population, race, and location levels. We observed strong colocalized signals in a narrow genomic interval for three interrelated traits: growth habit, lodging, and canopy height. Overall, this study detected ~30 candidate genes based on a priori and candidate gene search strategies centered on the 100-kb region surrounding a significant SNP. These results provide a framework from which further research can begin to understand the actual genes controlling important agronomic production traits in common bean.

Iron defi ciency chlorosis (IDC) is a signifi cant yield-limiting problem in several major soybean [Glycine max (L.) Merr.] production regions in the United States. Soybean plants display a variety of symptoms that range from a slight yellowing of the leaf to interveinal chlorosis, to stunted growth that reduces yield. The objective of this analysis was to employ single nucleotide polymorphism (SNP)-based genome-wide association mapping to uncover genomic regions associated with IDC tolerance. A number of the loci discovered here mapped at or near previously discovered IDC quantitative trait loci (QTL). A total of 15 genes known to be involved in iron metabolism mapped in the vicinity (<500 kb) of signifi cant markers in one or both populations.

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (p sil ¼ 0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F st ¼ 0.29). It is estimated that the two current wild gene pools diverged from a common ancestor B111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted B40 000 years. The Middle American bottleneck population size was B46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

The presence of seed color in common bean (Phaseolus vulgaris) requires the dominant-acting P (pigment) gene, and white seed is a recessive phenotype in all domesticated races of the species. P was classically associated with seed size, thus describing it as the first genetic marker for a quantitative trait. The molecular structure of P was characterized to understand the selection of white seeds during bean diversification and the relationship of P to seed weight. P was identified by homology searches, a genome-wide association study (GWAS) and gene remodeling, and confirmed by gene silencing. Allelic variation was assessed by a combination of resequencing and marker development, and the relationship between P and seed weight was assessed by a GWAS study. P is a member of clade B of subclass IIIf of plant basic helix-loop-helix (bHLH) proteins. Ten race-specific P alleles conditioned the white seed phenotype, and each causative mutation affected at least one bHLH domain required for color expression. GWAS analysis confirmed the classic association of P with seed weight. In common bean, white seeds are the result of convergent evolution and, among plant species, orthologous convergence on a single transcription factor gene was observed.

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

Abstract. Multilocus sequence data collected from domesticated and related wild relatives provides a rich source of information on the effect of human selection on the diversity and adaptability of a species to complex environments. To evaluate the domestication history of common bean (Phaseolus vulgaris L.), multilocus sequence data from landraces representing the various races within the Middle American (MA) and Andean gene pools was evaluated. Across 13 loci, nucleotide diversity was similar between landraces and wild germplasm in both gene pools. The diversity data were evaluated using the approximate Bayesian computation approach to test multiple domestication models and estimate population demographic parameters. A model with a single domestication event coupled with bidirectional migration between wild and domesticated genotypes fitted the data better than models consisting of two or three domestication events in each genepool. The effective bottleneck population size was~50% of the base population in each genepool. The bottleneck began~8200 and 8500 years before present and ended at~6300 and~7000 years before present in MA and Andean gene pools respectively. Linkage disequilibrium decayed to a greater extent in the MA genepool. Given the (1) geographical adaptation bottleneck in each wild gene pool, (2) a subsequent domestication bottleneck within each gene pool, (3) differentiation into gene-pool specific races and (4) variable extents of linkage disequilibrium, association mapping experiments for common bean would more appropriately be performed within each genepool.

Next generation sequence data provides valuable information and tools for genetic and genomic research and offers new insights useful for marker development. This data is useful for the design of accurate and user-friendly molecular tools. Common bean (Phaseolus vulgaris L.) is a diverse crop in which separate domestication events happened in each gene pool followed by race and market class diversification that has resulted in different morphological characteristics in each commercial market class. This has led to essentially independent breeding programs within each market class which in turn has resulted in limited within market class sequence variation. Sequence data from selected genotypes of five bean market classes (pinto, black, navy, and light and dark red kidney) were used to develop InDel-based markers specific to each market class. Design of the InDel markers was conducted through a combination of assembly, alignment and primer design software using 1.6× to 5.1× coverage of Illumina GAII sequence data for each of the selected genotypes. The procedure we developed for primer design is fast, accurate, less error prone, and higher throughput than when they are designed manually. All InDel markers are easy to run and score with no need for PCR optimization. A total of 2687 InDel markers distributed across the genome were developed. To highlight their usefulness, they were employed to construct a phylogenetic tree and a genetic map, showing that InDel markers are reliable, simple, and accurate.