Investigation of the domestication of common bean (Phaseolus vulgaris) using multilocus sequence data

Morphological and molecular characterization of variation in common bean (Phaseolus vulgaris L.) germplasm from Azad Jammu and Kashmir, Pakistan

Phaseolus vulgaris, an essential food and source of protein, is cultivated across the world. This study was carried out to investigate the diversity and population structure of 34 P. vulgaris landrace accessions collected from the Azad Jammu and Kashmir (AJ&K) regions of Pakistan. The samples were analyzed both morphologically and using genetic variation identified through RNA sequencing. Our results indicated that most genetic variation occurs among local accessions, with little genetic variation occurring between geographical regions. In addition, the accessions fell into two major genetic groups. Morphological analysis revealed that these two genetic groups differ in a number of quantitative traits, including seed length, seed width, and seed weight. One accession, DUD-11, appears to be a mixture of the two major groups genetically as well as morphologically. Among the other accessions, DUD-8, RWK-2, and NGD-1 depicted particularly high seed weight along with higher seed length, seed width, and seed yield per plant. We suggest focusing on these accessions in future breeding programs. More generally, our results provide baseline data that will be useful for crop improvement and effective cultivation practices in Pakistan.

DISTRIBUTION AND BIOLOGY OF INVASIVE SPECIES OF BEAN BRUCHID <i>ACANTHOSCELIDES OBTECTUS</i> (INSECTA, COLEOPTERA, BRUCHIDAE)

The review of literary sources on ecology, biology, distribution of bean bruchid ( Acanthoscelides obtectus ) and its main food plant - Phaseolus vulgaris in North and South America; Europe, Africa, Asia, Australia and more details in Russia; the influence of abiotic, biotic and anthropogenic factors on the invasive process, phytosanitary condition of common bean crops in Russia is presented. Some aspects of the invader management are shown. The main stages and areas of cultivation of common bean and invasion of bean bruchid from their primary habitat in South America and in the south of North America are traced; the vectors and reasons causing them are considered. In Russia, the economic importance of bean bruchid has increased since the mid-1980s, which coincided with the climate warming; there was an expansion of its distribution in the eastern and north-western directions. At the last decades of the 20th century, it had penetrated in Smolensk and in the south part of the Tver and the Tomsk regions. With the increase in production of beans in Russia, the lack of systemic protection from bean bruchid and further increase of climate warming will contribute to the extension of its range to the north in the European part of Russia and the Urals to 57-58° N. Lat., where the conditions of the summer period are favorable for development of common bean and bean bruchid. To the east, it may spread to Tyva, Buryatia, the Trans-Baikal territory, the Amur region, the Jewish Autonomous region, and the southern part of the Khabarovsk territory. With the introduction of strict internal quarantine and a system of protection of common bean from this pest, which prevents the spread of infected dry bean, on the contrary, it is possible to reduce the distribution range of the bean bruchid, with its disappearance in the Siberian, Ural districts, Bashkortostan and Tatarstan.

Genome-Environment Association Analysis for Bio-Climatic Variables in Common Bean (Phaseolus vulgaris L.) from Brazil

Abiotic stress is a limiting factor for common bean (Phaseolus vulgaris L.) production globally. The study of the genotypic, phenotypic, and bio-climatic variables in a broad set of accessions may assist the identification of genomic regions involved in the climatic adaptation of the common bean. We conducted a genotyping-by-sequencing analysis using 28,823 SNPs on 110 georeferenced common bean accessions from Brazil to discover associations between SNPs and bio-climatic indexes. The population structure analysis clustered the accessions into two groups corresponding to the Andean and Mesoamerican gene pools. Of the 19 bioclimatic variables, 17 exhibited a significant association with SNPs on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, Pv09, Pv10, and Pv11 of common bean. Ten candidate genes were associated with specific bio-climatic variables related to temperature and precipitation. The candidate genes associated with this significant Pv09 region encode a Platz transcription factor family protein previously reported to be an essential regulator of drought stress. The SNP markers and candidate genes associated with the bio-climatic variables should be validated in segregating populations for water stress, which could further be used for marker-assisted selection. As a result, bean breeding programs may be able to provide advances in obtaining drought-tolerant cultivars.

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

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

Genome-Wide Association Study Reveals Genomic Regions Associated with Fusarium Wilt Resistance in Common Bean

Fusarium wilt (Fusarium oxysporum f. sp. phaseoli, Fop) is one of the main fungal soil diseases in common bean. The aim of the present study was to identify genomic regions associated with Fop resistance through genome-wide association studies (GWAS) in a Mesoamerican Diversity Panel (MDP) and to identify potential common bean sources of Fop's resistance. The MDP was genotyped with BARCBean6K_3BeadChip and evaluated for Fop resistance with two different monosporic strains using the root-dip method. Disease severity rating (DSR) and the area under the disease progress curve (AUDPC), at 21 days after inoculation (DAI), were used for GWAS performed with FarmCPU model. The p-value of each SNP was determined by resampling method and Bonferroni test. For UFV01 strain, two significant single nucleotide polymorphisms (SNPs) were mapped on the Pv05 and Pv11 for AUDPC, and the same SNP (ss715648096) on Pv11 was associated with AUDPC and DSR. Another SNP, mapped on Pv03, showed significance for DSR. Regarding IAC18001 strain, significant SNPs on Pv03, Pv04, Pv05, Pv07 and on Pv01, Pv05, and Pv10 were observed. Putative candidate genes related to nucleotide-binding sites and carboxy-terminal leucine-rich repeats were identified. The markers may be important future tools for genomic selection to Fop disease resistance in beans.

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.

Genetic diversity and structure of landrace accessions, elite lineages and cultivars of common bean estimated with SSR and SNP markers

Common bean (Phaseolus vulgaris L.) is an important source of proteins, fibers and minerals for humans, being grown mainly in developing countries and representing a source of income for small farmers. In this work, a set of 206 Brazilian landraces and 59 elite lineages and cultivars were genotyped with 23 SSR (Simple Sequence Repeats) and 251 SNPs (Single-Nucleotide Polymorphism) markers. The ideal number of groups, according to STRUCTURE, was K = 2 for both SNPs and SSRs. This could be expected considering the two original gene pools-Andean (AND) and Mesoamerican (MES). The matrices of genetic simple matching dissimilarity for SSRs and SNPs were highly correlated; therefore, the allelic data of the markers was combined and analyzed to understand the genetic relationships of the studied collection. The neighbor-joining analysis considering the genetic distance of simple matching grouped the 265 genotypes into 17 subgroups. The markers SSR and SNP presented high power to discriminate among the genotypes. The ample genetic diversity observed in the work collection makes it a valuable source for the conservation, sustainable management and exploration in breeding programs of the crop.

Marker-Trait Association Analysis of Seed Traits in Accessions of Common Bean (Phaseolus vulgaris L.) in China

Seed weight and seed size are the key agronomic traits that determine yield in common bean. To investigate the genetic architecture of four seed traits (100-seed weight, seed length, seed width, and seed height) of common bean in China, marker-trait association analysis of these seed traits was performed in a nationwide population of 395 common bean accessions using 116 polymorphic SSR markers. The four seed traits were evaluated in six trials across three environments. Seed size varied among the environments. Population structure was evaluated based on SSR markers and phaseolin, which divided the accessions into two main subpopulations representing the two known gene pools. Seed weight and seed size had a strong relationship with population clustering. In addition, in a Genome-wide association studies (GWAS), 21 significantly associated markers were identified for the four seed traits with two models, namely, general linear model (GLM) and mixed linear model (MLM). Some markers had pleiotropic effects, i.e., controlled more than one trait. The significant quantitative trait loci identified in this study could be used in marker-assisted breeding to accelerate the genetic improvement of yield in common bean.

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.

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.

Current State and Perspectives in Population Genomics of the Common Bean

* Correspondence: r [...].

Root and shoot variation in relation to potential intermittent drought adaptation of Mesoamerican wild common bean (Phaseolus vulgaris L.)

BACKGROUND

Wild crop relatives have been potentially subjected to stresses on an evolutionary time scale prior to domestication. Among these stresses, drought is one of the main factors limiting crop productivity and its impact is likely to increase under current scenarios of global climate change. We sought to determine to what extent wild common bean (Phaseolus vulgaris) exhibited adaptation to drought stress, whether this potential adaptation is dependent on the climatic conditions of the location of origin of individual populations, and to what extent domesticated common bean reflects potential drought adaptation.

METHODS

An extensive and diverse set of wild beans from across Mesoamerica, along with a set of reference Mesoamerican domesticated cultivars, were evaluated for root and shoot traits related to drought adaptation. A water deficit experiment was conducted by growing each genotype in a long transparent tube in greenhouse conditions so that root growth, in addition to shoot growth, could be monitored.

RESULTS

Phenotypic and landscape genomic analyses, based on single-nucleotide polymorphisms, suggested that beans originating from central and north-west Mexico and Oaxaca, in the driest parts of their distribution, produced more biomass and were deeper-rooted. Nevertheless, deeper rooting was correlated with less root biomass production relative to total biomass. Compared with wild types, domesticated types showed a stronger reduction and delay in growth and development in response to drought stress. Specific genomic regions were associated with root depth, biomass productivity and drought response, some of which showed signals of selection and were previously related to productivity and drought tolerance.

CONCLUSIONS

The drought tolerance of wild beans consists in its stronger ability, compared with domesticated types, to continue growth in spite of water-limited conditions. This study is the first to relate bean response to drought to environment of origin for a diverse selection of wild beans. It provides information that needs to be corroborated in crosses between wild and domesticated beans to make it applicable to breeding programmes.

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

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

Agronomic Performance and Nitrogen Fixation of Heirloom and Conventional Dry Bean Varieties Under Low-Nitrogen Field Conditions

Common beans (Phaseolus vulgaris) form a relationship with nitrogen-fixing rhizobia and through a process termed symbiotic nitrogen fixation (SNF) which provides them with a source of nitrogen. However, beans are considered poor nitrogen fixers, and modern production practices involve routine use of N fertilizer, which leads to the down-regulation of SNF. High-yielding, conventionally bred bean varieties are developed using conventional production practices and selection criteria, typically not including SNF efficiency, and may have lost this trait over decades of modern breeding. In contrast, heirloom bean genotypes were developed before the advent of modern production practices and may represent an underutilized pool of genetics which could be used to improve SNF. This study compared the SNF capacity under low-N field conditions, of collections of heirloom varieties with and conventionally bred dry bean varieties. The heirloom-conventional panel (HCP) consisted of 42 genotypes from various online seed retailers or from the University of Guelph Bean Breeding program seedbank. The HCP was genotyped using a single nucleotide polymorphism (SNP) array to investigate genetic relatedness within the panel. Field trials were conducted at three locations in ON, Canada from 2014 to 2015 and various agronomic and seed composition traits were measured, including capacity for nitrogen fixation (using the natural abundance method to measure seed N isotope ratios). Significant variation for SNF was found in the panel. However, on average, heirloom genotypes did not fix significantly more nitrogen than conventionally bred varieties. However, five heirloom genotypes fixed >60% of their nitrogen from the atmosphere. Yield (kg ha-1) was not significantly different between heirloom and conventional genotypes, suggesting that incorporating heirloom genotypes into a modern breeding program would not negatively impact yield. Nitrogen fixation was significantly higher among Middle American genotypes than among Andean genotypes, confirming previous findings. The best nitrogen fixing line was Coco Sophie, a European heirloom white bean whose genetic makeup is admixed between the Andean and Middle American genepools. Heirloom genotypes represent a useful source of genetics to improve SNF in modern bean breeding.

Genotypes and Genomic Regions Associated With Rhizoctonia solani Resistance in Common Bean

Rhizoctonia solani Kühn (teleomorph Thanatephorus cucumeris) is an important root rot pathogen of common bean (Phaseolus vulgaris L.). To uncover genetic factors associated with resistance to the pathogen, the Andean (ADP; n = 273) and Middle American (MDP; n = 279) diversity panels, which represent much of the genetic diversity known in cultivated common bean, were screened in the greenhouse using R. solani anastomosis group 2-2. Repeatability of the assay was confirmed by the response of five control genotypes. The phenotypic data for both panels were normally distributed. The resistance responses of ∼10% of the ADP (n = 28) and ∼6% of the MDP (n = 18) genotypes were similar or higher than that of the resistant control line VAX 3. A genome-wide association study (GWAS) was performed using ∼200k single nucleotide polymorphisms to discover genomic regions associated with resistance in each panel, For GWAS, the raw phenotypic score, and polynomial and binary transformation of the scores, were individually used as the input data. A major QTL peak was observed on Pv02 in the ADP, while a major QTL was observed on Pv01 with the MDP. These regions were associated with clusters of TIR-NB_ARC-LRR (TNL) gene models encoding proteins similar to known disease resistance genes. Other QTL, unique to each panel, were mapped within or adjacent to a gene model or cluster of related genes associated with disease resistance. This is a first case study that provides evidence for major as well as minor genes involved in resistance to R. solani in common bean. This information will be useful to integrate more durable root rot resistance in common bean breeding programs and to study the genetic mechanisms associated with root diseases in this important societal legume.

Multi-tissue integration of transcriptomic and specialized metabolite profiling provides tools for assessing the common bean (Phaseolus vulgaris) metabolome

Common bean (Phaseolus vulgaris L.) is an important legume species with a rich natural diversity of landraces that originated from the wild forms following multiple independent domestication events. After the publication of its genome, several resources for this relevant crop have been made available. A comprehensive characterization of specialized metabolism in P. vulgaris, however, is still lacking. In this study, we used a metabolomics approach based on liquid chromatography-mass spectrometry to dissect the chemical composition at a tissue-specific level in several accessions of common bean belonging to different gene pools. Using a combination of literature search, mass spectral interpretation, 13 C-labeling, and correlation analyses, we were able to assign chemical classes and/or putative structures for approximately 39% of all measured metabolites. Additionally, we integrated this information with transcriptomics data and phylogenetic inference from multiple legume species to reconstruct the possible metabolic pathways and identify sets of candidate genes involved in the biosynthesis of specialized metabolites. A particular focus was given to flavonoids, triterpenoid saponins and hydroxycinnamates, as they represent metabolites involved in important ecological interactions and they are also associated with several health-promoting benefits when integrated into the human diet. The data are presented here in the form of an accessible resource that we hope will set grounds for further studies on specialized metabolism in legumes.

Genome wide association study discovers genomic regions involved in resistance to soybean cyst nematode (Heterodera glycines) in common bean

Common bean (Phaseolus vulgaris L.) is an important high protein crop grown worldwide. North Dakota and Minnesota are the largest producers of common beans in the USA, but crop production is threatened by soybean cyst nematode (SCN; Heterodera glycines Ichinohe) because most current cultivars are susceptible. Greenhouse screening data using SCN HG type 0 from 317 plant introductions (PI's) from the USDA core collection was used to conduct a genome wide association study (GWAS). These lines were divided into two subpopulations based on principal component analysis (Middle American vs. Andean). Phenotypic results based on the female index showed that accessions could be classified as highly resistant (21% and 27%), moderately resistant (51% and 48%), moderately susceptible (27% and 22%) and highly susceptible (1% and 3%) for Middle American and Andean gene pools, respectively. Mixed models with two principal components (PCs) and kinship matrix for Middle American genotypes and Andean genotypes were used in the GWAS analysis using 3,985 and 4,811 single nucleotide polymorphic (SNP) markers, respectively which were evenly distributed across all 11 chromosomes. Significant peaks on Pv07, and Pv11 in Middle American and on Pv07, Pv08, Pv09 and Pv11 in Andean group were found to be associated with SCN resistance. Homologs of soybean rhg1, a locus which confers resistance to SCN in soybean, were identified on chromosomes Pv01 and Pv08 in the Middle American and Andean gene pools, respectively. These genomic regions may be the key to develop SCN-resistant common bean cultivars.

Characterization of genetic diversity in Turkish common bean gene pool using phenotypic and whole-genome DArTseq-generated silicoDArT marker information

Turkey presents a great diversity of common bean landraces in farmers' fields. We collected 183 common bean accessions from 19 different Turkish geographic regions and 5 scarlet runner bean accessions to investigate their genetic diversity and population structure using phenotypic information (growth habit, and seed weight, flower color, bracteole shape and size, pod shape and leaf shape and color), geographic provenance and 12,557 silicoDArT markers. A total of 24.14% markers were found novel. For the entire population (188 accessions), the expected heterozygosity was 0.078 and overall gene diversity, Fst and Fis were 0.14, 0.55 and 1, respectively. Using marker information, model-based structure, principal coordinate analysis (PCoA) and unweighted pair-group method with arithmetic means (UPGMA) algorithms clustered the 188 accessions into two main populations A (predominant) and B, and 5 unclassified genotypes, representing 3 meaningful heterotic groups for breeding purposes. Phenotypic information clearly distinguished these populations; population A and B, respectively, were bigger (>40g/100 seeds) and smaller (<40g/100 seeds) seed-sized. The unclassified population was pure and only contained climbing genotypes with 100 seed weight 2-3 times greater than populations A and B. Clustering was mainly based on A: seed weight, B: growth habit, C: geographical provinces and D: flower color. Mean kinship was generally low, but population B was more diverse than population A. Overall, a useful level of gene and genotypic diversity was observed in this work and can be used by the scientific community in breeding efforts to develop superior common bean strains.

White seed color in common bean (Phaseolus vulgaris) results from convergent evolution in the P (pigment) gene

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.

Direct archaeological evidence for Southwestern Amazonia as an early plant domestication and food production centre

Southwestern Amazonia is considered an early centre of plant domestication in the New World, but most of the evidence for this hypothesis comes from genetic data since systematic archaeological fieldwork in the area is recent. This paper provides first-hand archaeobotanical evidence of food production from early and middle Holocene (ca. 9,000-5000 cal. BP) deposits at Teotonio, an open-air site located on a 40 m-high bluff on the south bank of the Madeira river. Such evidence includes the presence of local and exotic domesticates such as manioc (Manihot esculenta), squash (Cucurbita sp.) and beans (Phaseolus sp.), alongside edible fruits such as pequiá (Caryocar sp.) and guava (Psidium sp.) that point to the beginnings of landscape domestication. The results contribute to an ever-growing number of studies that posit southwest Amazonia as an important centre for early crop domestication and experimentation, and which highlight the longue-durée of human impacts on tropical forest biodiversity around the world.

Resequencing of Common Bean Identifies Regions of Inter-Gene Pool Introgression and Provides Comprehensive Resources for Molecular Breeding

Common bean ( L.) is the most important grain legume for human consumption and is a major nutrition source in the tropics. Because bean production is reduced by both abiotic and biotic constraints, current breeding efforts are focused on the development of improved varieties with tolerance to these stresses. We characterized materials from different breeding programs spanning three continents to understand their sequence diversity and advance the development of molecular breeding tools. For this, 37 varieties belonging to , (A. Gray), and L. were sequenced by whole-genome sequencing, identifying more than 40 million genomic variants. Evaluation of nuclear DNA content and analysis of copy number variation revealed important differences in genomic content not only between and the two other domesticated species, but also within , affecting hundreds of protein-coding genomic regions. A large number of inter-gene pool introgressions were identified. Furthermore, interspecific introgressions for disease resistance in breeding lines were mapped. Evaluation of newly developed single nucleotide polymorphism markers within previously discovered quantitative trait loci for common bacterial blight and angular leaf spot provides improved specificity to tag sources of resistance to these diseases. We expect that this dataset will provide a deeper molecular understanding of breeding germplasm and deliver molecular tools for germplasm development, aiming to increase the efficiency of bean breeding programs.

Testing Domestication Scenarios of Lima Bean (Phaseolus lunatus L.) in Mesoamerica: Insights from Genome-Wide Genetic Markers

Plant domestication can be seen as a long-term process that involves a complex interplay among demographic processes and evolutionary forces. Previous studies have suggested two domestication scenarios for Lima bean in Mesoamerica: two separate domestication events, one from gene pool MI in central-western Mexico and another one from gene pool MII in the area Guatemala-Costa Rica, or a single domestication from gene pool MI in central-western Mexico followed by post-domestication gene flow with wild populations. In this study we evaluated the genetic structure of the wild gene pool and tested these two competing domestication scenarios of Lima bean in Mesoamerica by applying an ABC approach to a set of genome-wide SNP markers. The results confirm the existence of three gene pools in wild Lima bean, two Mesoamerican gene pools (MI and MII) and the Andean gene pool (AI), and suggest the existence of another gene pool in central Colombia. The results indicate that although both domestication scenarios may be supported by genetic data, higher statistical support was given to the single domestication scenario in central-western Mexico followed by admixture with wild populations. Domestication would have involved strong founder effects reflected in loss of genetic diversity and increased LD levels in landraces. Genomic regions affected by selection were detected and these may harbor candidate genes related to domestication.

In-depth genome characterization of a Brazilian common bean core collection using DArTseq high-density SNP genotyping

BACKGROUND

Common bean is a legume of social and nutritional importance as a food crop, cultivated worldwide especially in developing countries, accounting for an important source of income for small farmers. The availability of the complete sequences of the two common bean genomes has dramatically accelerated and has enabled new experimental strategies to be applied for genetic research. DArTseq has been widely used as a method of SNP genotyping allowing comprehensive genome coverage with genetic applications in common bean breeding programs.

RESULTS

Using this technology, 6286 SNPs (1 SNP/86.5 Kbp) were genotyped in genic (43.3%) and non-genic regions (56.7%). Genetic subdivision associated to the common bean gene pools (K = 2) and related to grain types (K = 3 and K = 5) were reported. A total of 83% and 91% of all SNPs were polymorphic within the Andean and Mesoamerican gene pools, respectively, and 26% were able to differentiate the gene pools. Genetic diversity analysis revealed an average H E of 0.442 for the whole collection, 0.102 for Andean and 0.168 for Mesoamerican gene pools (F ST  = 0.747 between gene pools), 0.440 for the group of cultivars and lines, and 0.448 for the group of landrace accessions (F ST  = 0.002 between cultivar/line and landrace groups). The SNP effects were predicted with predominance of impact on non-coding regions (77.8%). SNPs under selection were identified within gene pools comparing landrace and cultivar/line germplasm groups (Andean: 18; Mesoamerican: 69) and between the gene pools (59 SNPs), predominantly on chromosomes 1 and 9. The LD extension estimate corrected for population structure and relatedness (r2SV) was ~ 88 kbp, while for the Andean gene pool was ~ 395 kbp, and for the Mesoamerican was ~ 130 kbp.

CONCLUSIONS

For common bean, DArTseq provides an efficient and cost-effective strategy of generating SNPs for large-scale genome-wide studies. The DArTseq resulted in an operational panel of 560 polymorphic SNPs in linkage equilibrium, providing high genome coverage. This SNP set could be used in genotyping platforms with many applications, such as population genetics, phylogeny relation between common bean varieties and support to molecular breeding approaches.

Fast diffusion of domesticated maize to temperate zones

Adaptation to a temperate climate was a prerequisite for the spread of maize across a broad geographical range. To explicitly explore the demographic process underlying maize adaptation, we used a diffusion-based method to model the differentiation between temperate and tropical populations using the Non-Stiff Stalk group as a proxy for temperate maize. Based on multiple sequential Markovian coalescent approaches, we estimate that tropical and temperate maize diverged approximately 3'000 to 5'000 years ago and the population size shrank after the split. Using composite likelihood approaches, we identified a distinct tropical-temperate divergence event initiated 4'958 years ago (95% confidence interval (CI): 4'877-5'039) from an ancestral population whose effective size was 24,162 (95% CI: 23,914-24,409). We found that continuous gene flow between tropical and temperate maize accompanied the differentiation of temperate maize. Long identical-by-descent tracts shared by tropical and temperate inbred lines have been identified, which might be the result of gene flow between tropical and temperate maize or artificial selection during domestication and crop improvement. Understanding the demographic history of maize diffusion not only provides evidence for population dynamics of maize, but will also assist the identification of regions under selection and the genetic basis of complex traits of agronomic importance.

Beans (Phaseolus ssp.) as a Model for Understanding Crop Evolution

Here, we aim to provide a comprehensive and up-to-date overview of the most significant outcomes in the literature regarding the origin of Phaseolus genus, the geographical distribution of the wild species, the domestication process, and the wide spread out of the centers of origin. Phaseolus can be considered as a unique model for the study of crop evolution, and in particular, for an understanding of the convergent phenotypic evolution that occurred under domestication. The almost unique situation that characterizes the Phaseolus genus is that five of its ∼70 species have been domesticated (i.e., Phaseolus vulgaris, P. coccineus, P. dumosus, P. acutifolius, and P. lunatus), and in addition, for P. vulgaris and P. lunatus, the wild forms are distributed in both Mesoamerica and South America, where at least two independent and isolated episodes of domestication occurred. Thus, at least seven independent domestication events occurred, which provides the possibility to unravel the genetic basis of the domestication process not only among species of the same genus, but also between gene pools within the same species. Along with this, other interesting features makes Phaseolus crops very useful in the study of evolution, including: (i) their recent divergence, and the high level of collinearity and synteny among their genomes; (ii) their different breeding systems and life history traits, from annual and autogamous, to perennial and allogamous; and (iii) their adaptation to different environments, not only in their centers of origin, but also out of the Americas, following their introduction and wide spread through different countries. In particular for P. vulgaris this resulted in the breaking of the spatial isolation of the Mesoamerican and Andean gene pools, which allowed spontaneous hybridization, thus increasing of the possibility of novel genotypes and phenotypes. This knowledge that is associated to the genetic resources that have been conserved ex situ and in situ represents a crucial tool in the hands of researchers, to preserve and evaluate this diversity, and at the same time, to identify the genetic basis of adaptation and to develop new improved varieties to tackle the challenges of climate change, and food security and sustainability.

High Level of Nonsynonymous Changes in Common Bean Suggests That Selection under Domestication Increased Functional Diversity at Target Traits

Crop species have been deeply affected by the domestication process, and there have been many efforts to identify selection signatures at the genome level. This knowledge will help geneticists to better understand the evolution of organisms, and at the same time, help breeders to implement successful breeding strategies. Here, we focused on domestication in the Mesoamerican gene pool of Phaseolus vulgaris by sequencing 49 gene fragments from a sample of 45 P. vulgaris wild and domesticated accessions, and as controls, two accessions each of the closely related species Phaseolus coccineus and Phaseolus dumosus. An excess of nonsynonymous mutations within the domesticated germplasm was found. Our data suggest that the cost of domestication alone cannot explain fully this finding. Indeed, the significantly higher frequency of polymorphisms in the coding regions observed only in the domesticated plants (compared to noncoding regions), the fact that these mutations were mostly nonsynonymous and appear to be recently derived mutations, and the investigations into the functions of their relative genes (responses to biotic and abiotic stresses), support a scenario that involves new functional mutations selected for adaptation during domestication. Moreover, consistent with this hypothesis, selection analysis and the possibility to compare data obtained for the same genes in different studies of varying sizes, data types, and methodologies allowed us to identify four genes that were strongly selected during domestication. Each selection candidate is involved in plant resistance/tolerance to abiotic stresses, such as heat, drought, and salinity. Overall, our study suggests that domestication acted to increase functional diversity at target loci, which probably controlled traits related to expansion and adaptation to new agro-ecological growing conditions.

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

Common bean ( 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 Durango-Jalisco 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.

To Have and to Hold: Selection for Seed and Fruit Retention During Crop Domestication

Crop domestication provides a useful model system to characterize the molecular and developmental bases of morphological variation in plants. Among the most universal changes resulting from selection during crop domestication is the loss of seed and fruit dispersal mechanisms, which greatly facilitates harvesting efficiency. In this review, we consider the molecular genetic and developmental bases of the loss of seed shattering and fruit dispersal in six major crop plant families, three of which are primarily associated with seed crops (Poaceae, Brassicaceae, Fabaceae) and three of which are associated with fleshy-fruited crops (Solanaceae, Rosaceae, Rutaceae). We find that the developmental basis of the loss of seed/fruit dispersal is conserved in a number of independently domesticated crops, indicating the widespread occurrence of developmentally convergent evolution in response to human selection. With regard to the molecular genetic approaches used to characterize the basis of this trait, traditional biparental quantitative trait loci mapping remains the most commonly used strategy; however, recent advances in next-generation sequencing technologies are now providing new avenues to map and characterize loss of shattering/dispersal alleles. We anticipate that continued application of these approaches, together with candidate gene analyses informed by known shattering candidate genes from other crops, will lead to a rapid expansion of our understanding of this critical domestication trait.

The Domestication Syndrome in Phoenix dactylifera Seeds: Toward the Identification of Wild Date Palm Populations

Investigating crop origins is a priority to understand the evolution of plants under domestication, develop strategies for conservation and valorization of agrobiodiversity and acquire fundamental knowledge for cultivar improvement. The date palm (Phoenix dactylifera L.) belongs to the genus Phoenix, which comprises 14 species morphologically very close, sometimes hardly distinguishable. It has been cultivated for millennia in the Middle East and in North Africa and constitutes the keystone of oasis agriculture. Yet, its origins remain poorly understood as no wild populations are identified. Uncultivated populations have been described but they might represent feral, i.e. formerly cultivated, abandoned forms rather than truly wild populations. In this context, this study based on morphometrics applied to 1625 Phoenix seeds aims to (1) differentiate Phoenix species and (2) depict the domestication syndrome observed in cultivated date palm seeds using other Phoenix species as a "wild" reference. This will help discriminate truly wild from feral forms, thus providing new insights into the evolutionary history of this species. Seed size was evaluated using four parameters: length, width, thickness and dorsal view surface. Seed shape was quantified using outline analyses based on the Elliptic Fourier Transform method. The size and shape of seeds allowed an accurate differentiation of Phoenix species. The cultivated date palm shows distinctive size and shape features, compared to other Phoenix species: seeds are longer and elongated. This morphological shift may be interpreted as a domestication syndrome, resulting from the long-term history of cultivation, selection and human-mediated dispersion. Based on seed attributes, some uncultivated date palms from Oman may be identified as wild. This opens new prospects regarding the possible existence and characterization of relict wild populations and consequently for the understanding of the date palm origins. Finally, we here describe a pipeline for the identification of the domestication syndrome in seeds that could be used in other crops.

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Landrace Germplasm for Improving Yield and Abiotic Stress Adaptation

Plant landraces represent heterogeneous, local adaptations of domesticated species, and thereby provide genetic resources that meet current and new challenges for farming in stressful environments. These local ecotypes can show variable phenology and low-to-moderate edible yield, but are often highly nutritious. The main contributions of landraces to plant breeding have been traits for more efficient nutrient uptake and utilization, as well as useful genes for adaptation to stressful environments such as water stress, salinity, and high temperatures. We propose that a systematic landrace evaluation may define patterns of diversity, which will facilitate identifying alleles for enhancing yield and abiotic stress adaptation, thus raising the productivity and stability of staple crops in vulnerable environments.

Identification of novel drought-tolerant-associated SNPs in common bean (Phaseolus vulgaris)

Common bean (Phaseolus vulgaris L.) is a leguminous in high demand for human nutrition and a very important agricultural product. Production of common bean is constrained by environmental stresses such as drought. Although conventional plant selection has been used to increase production yield and stress tolerance, drought tolerance selection based on phenotype is complicated by associated physiological, anatomical, cellular, biochemical, and molecular changes. These changes are modulated by differential gene expression. A common method to identify genes associated with phenotypes of interest is the characterization of Single Nucleotide Polymorphims (SNPs) to link them to specific functions. In this work, we selected two drought-tolerant parental lines from Mesoamerica, Pinto Villa, and Pinto Saltillo. The parental lines were used to generate a population of 282 families (F3:5) and characterized by 169 SNPs. We associated the segregation of the molecular markers in our population with phenotypes including flowering time, physiological maturity, reproductive period, plant, seed and total biomass, reuse index, seed yield, weight of 100 seeds, and harvest index in three cultivation cycles. We observed 83 SNPs with significant association (p < 0.0003 after Bonferroni correction) with our quantified phenotypes. Phenotypes most associated were days to flowering and seed biomass with 58 and 44 associated SNPs, respectively. Thirty-seven out of the 83 SNPs were annotated to a gene with a potential function related to drought tolerance or relevant molecular/biochemical functions. Some SNPs such as SNP28 and SNP128 are related to starch biosynthesis, a common osmotic protector; and SNP18 is related to proline biosynthesis, another well-known osmotic protector.

Inference of domestication history and differentiation between early- and late-flowering varieties in pearl millet

Pearl millet (Pennisetum glaucum) is a staple crop in Sahelian Africa. Farmers usually grow varieties with different cycle lengths and complementary functions in Sahelian agrosystems. Both the level of genetic differentiation of these varieties and the domestication history of pearl millet have been poorly studied. We investigated the neutral genetic diversity and population genetic structure of early- and late-flowering domesticated and wild pearl millet populations using 18 microsatellite loci and 8 nucleotide sequences. Strikingly, early- and late-flowering domesticated varieties were not differentiated over their whole distribution area, despite a clear difference in their isolation-by-distance pattern. Conversely, our data brought evidence for two well-differentiated genetic pools in wild pearl millet, allowing us to test scenarios with different numbers and origins of domestication using approximate Bayesian computation (ABC). The ABC analysis showed the likely existence of asymmetric migration between wild and domesticated populations. The model choice procedure indicated that a single domestication from the eastern wild populations was the more likely scenario to explain the polymorphism patterns observed in cultivated pearl millet.

Genotyping-by-Sequencing Enabled Mapping and Marker Development for the By-2 Potyvirus Resistance Allele in Common Bean

Since its emergence in 2001, an aphid-transmitted virus disease complex has caused substantial economic losses to snap bean (Phaseolus vulgaris L.) production and processing in the Great Lakes Region of the United States. The general ineffectiveness of chemical control measures for nonpersistently transmitted viruses established an urgent need for the development and deployment of cultivars with resistance to the component viruses. Our objectives were to further characterize the inheritance of resistance to Bean yellow mosaic virus (BYMV), which is conditioned by the By-2 allele, to adapt genotyping-by-sequencing (GBS) to common bean to discover and genotype genome-wide single nucleotide polymorphisms (SNPs) in a set of recombinant inbred lines (RILs) derived from an introgression program, and to enable and validate marker-assisted selection for By-2. We optimized ApeKI for GBS in common bean and retained 7530 high-quality SNPs that segregated in our introgression RILs. A case-control genome-wide association study (GWAS) was used to discover 44 GBS SNPs that were strongly associated with the resistance phenotype and which delimited a 974 kb physical interval on the distal portion of chromosome 2. Seven of these SNPs were converted to single-marker Kompetitive Allele-Specific Polymerase chain reaction (KASP) assays and were demonstrated to be tightly linked to BYMV resistance in an F₂ population of 185 individuals. This research enables marker-assisted selection of By-2, provides enhanced resolution for fine mapping, and demonstrates the potential of GBS as a highly efficient, high-throughput genotyping platform for common bean breeding and genetics.

A reference genome for common bean and genome-wide analysis of dual domestications

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.

Developing market class specific InDel markers from next generation sequence data in Phaseolus vulgaris L

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.

Storytelling and story testing in domestication

The domestication of plants and animals marks one of the most significant transitions in human, and indeed global, history. Traditionally, study of the domestication process was the exclusive domain of archaeologists and agricultural scientists; today it is an increasingly multidisciplinary enterprise that has come to involve the skills of evolutionary biologists and geneticists. Although the application of new information sources and methodologies has dramatically transformed our ability to study and understand domestication, it has also generated increasingly large and complex datasets, the interpretation of which is not straightforward. In particular, challenges of equifinality, evolutionary variance, and emergence of unexpected or counter-intuitive patterns all face researchers attempting to infer past processes directly from patterns in data. We argue that explicit modeling approaches, drawing upon emerging methodologies in statistics and population genetics, provide a powerful means of addressing these limitations. Modeling also offers an approach to analyzing datasets that avoids conclusions steered by implicit biases, and makes possible the formal integration of different data types. Here we outline some of the modeling approaches most relevant to current problems in domestication research, and demonstrate the ways in which simulation modeling is beginning to reshape our understanding of the domestication process.

The contribution of genetic and genomic approaches to plant domestication studies

The application of genomic approaches to the phenomenon of plant domestication promises a better understanding of the origins of agriculture, but also of the way plant genomes in general are organized and expressed. Building on earlier genetic research, more detailed information has become available on the organization of genetic diversity at the genome level and the effects of gene flow on diversity in different regions of the genome. In addition, putative domestication genes have been identified through population genomics approaches (selective sweeps or divergence scanning). Further information has been obtained on the origin of domestication syndrome mutations and the dispersal and adaptation of crops after domestication. For the future, increasingly multidisciplinary approaches using combinations of genomics and other approaches will prevail.

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L

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 (π(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 ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% 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.

Molecular analysis of the parallel domestication of the common bean (Phaseolus vulgaris) in Mesoamerica and the Andes

We have studied the nucleotide diversity of common bean, Phaseolus vulgaris, which is characterized by two independent domestications in two geographically distinct areas: Mesoamerica and the Andes. This provides an important model, as domestication can be studied as a replicate experiment. We used nucleotide data from five gene fragments characterized by large introns to analyse 214 accessions (102 wild and 112 domesticated). The wild accessions represent a cross-section of the entire geographical distribution of P. vulgaris. A reduction in genetic diversity in both of these gene pools was found, which was three-fold greater in Mesoamerica compared with the Andes. This appears to be a result of a bottleneck that occurred before domestication in the Andes, which strongly impoverished this wild germplasm, leading to the minor effect of the subsequent domestication bottleneck (i.e. sequential bottleneck). These findings show the importance of considering the evolutionary history of crop species as a major factor that influences their current level and structure of genetic diversity. Furthermore, these data highlight a single domestication event within each gene pool. Although the findings should be interpreted with caution, this evidence indicates the Oaxaca valley in Mesoamerica, and southern Bolivia and northern Argentina in South America, as the origins of common bean domestication.

Why not beans?

Changes in climate and urbanisation rapidly affecting human livelihood are particularly threatening to developing nations in tropical regions. Food production crises have focused the global development agenda on agricultural research, a proven approach for increasing crop yield. A few crops benefit from private investment, but improvement of most crops will rely on limited public funding that must be deployed strategically, pushing forward both proven approaches and new ideas. Why not invest in beans? More than 300 million people rely on this crop, considered to be the most important grain legume for human consumption. Yet the yield of beans, especially in poor regions or marginal soils, is reduced by abiotic stresses such as phosphorus deficiency, aluminum toxicity and especially drought. Is it possible to assemble resources, including genetic diversity in beans, breeding expertise, genomic information and tools, and physiological insight to generate rapid progress in developing new lines of beans more tolerant to abiotic stress? A workshop to address this question was held in November 2010 at the International Center for Tropical Agriculture (CIAT) in Colombia. The resulting 'call to action' is presented in this issue which also includes research papers focused on tolerance of beans to stress.