A methyl esterase 1 (PvMES1) promotes the salicylic acid pathway and enhances Fusarium wilt resistance in common beans.

KEY MESSAGE: Methyl esterase (MES), PvMES1, contributes to the defense response toward Fusarium wilt in common beans by regulating the salicylic acid (SA) mediated signaling pathway from phenylpropanoid synthesis and sugar metabolism as well as others. Common bean (Phaseolus vulgaris L.) is an important food legume. Fusarium wilt caused by Fusarium oxysporum f. sp. phaseoli is one of the most serious soil-borne diseases of common bean found throughout the world and affects the yield and quality of the crop. Few sources of Fusarium wilt resistance exist in legumes and most are of quantitative inheritance. In this study, we have identified a methyl esterase (MES), PvMES1, that contributes to plant defense response by regulating the salicylic acid (SA) mediated signaling pathway in response to Fusarium wilt in common beans. The result showed the role of PvMES1 in regulating SA levels in common bean and thus the SA signaling pathway and defense response mechanism in the plant. Overexpression of the PvMES1 gene enhanced Fusarium wilt resistance; while silencing of the gene caused susceptibility to the diseases. RNA-seq analysis with these transiently modified plants showed that genes related to SA level changes included the following gene ontologies: (a) phenylpropanoid synthesis; (b) sugar metabolism; and (c) interaction between host and pathogen as well as others. These key signal elements activated the defense response pathway in common bean to Fusarium wilt. Collectively, our findings indicate that PvMES1 plays a pivotal role in regulating SA biosynthesis and signaling, and increasing Fusarium wilt resistance in common bean, thus providing novel insight into the practical applications of both SA and MES genes and pathways they contribute to for developing elite crop varieties with enhanced broad-spectrum resistance to this critical disease.

Meta-QTL analysis of seed iron and zinc concentration and content in common bean (Phaseolus vulgaris L.).

KEY MESSAGE: Twelve meta-QTL for seed Fe and Zn concentration and/or content were identified from 87 QTL originating from seven population grown in sixteen field trials. These meta-QTL include 2 specific to iron, 2 specific to zinc and 8 that co-localize for iron and zinc concentrations and/or content. Common bean (Phaseolus vulgaris L.) is the most important legume for human consumption worldwide and it is an important source of microelements, especially iron and zinc. Bean biofortification breeding programs develop new varieties with high levels of Fe and Zn targeted for countries with human micronutrient deficiencies. Biofortification efforts thus far have relied on phenotypic selection of raw seed mineral concentrations in advanced generations. While numerous quantitative trait loci (QTL) studies have been conducted to identify genomic regions associated with increased Fe and Zn concentration in seeds, these results have yet to be employed for marker-assisted breeding. The objective of this study was to conduct a meta-analysis from seven QTL studies in Andean and Middle American intra- and inter-gene pool populations to identify the regions in the genome that control the Fe and Zn levels in seeds. Two meta-QTL specific to Fe and two meta-QTL specific to Zn were identified. Additionally, eight Meta QTL that co-localized for Fe and Zn concentration and/or content were identified across seven chromosomes. The Fe and Zn shared meta-QTL could be useful candidates for marker-assisted breeding to simultaneously increase seed Fe and Zn. The physical positions for 12 individual meta-QTL were identified and within five of the meta-QTL, candidate genes were identified from six gene families that have been associated with transport of iron and zinc in plants.

Dissecting the Genetic Diversity of USDA Cowpea Germplasm Collection Using Kompetitive Allele Specific PCR-Single Nucleotide Polymorphism Markers.

Cowpea (Vigna unguiculata L. Walp) is an important grain legume crop of the subtropics, particularly in West Africa, where it contributes to the livelihoods of small-scale farmers. Despite being a drought-resilient crop, cowpea production is hampered by insect pests, diseases, parasitic weeds, and various abiotic stresses. Genetic improvement can help overcome these limitations, and exploring diverse cowpea genetic resources is crucial for cowpea breeding. This study evaluated the genetic diversity of 361 cowpea accessions from the USDA core collection for the species using 102 Kompetitive Allele Specific PCR (KASP) single nucleotide polymorphism (SNP) markers. A total of 102 KASP-SNP was validated in the germplasm panel, and 72 showed polymorphism across the germplasm panel. The polymorphism information content (PIC) of all SNPs ranged from 0.1 to 0.37, with an average of 0.29, while the mean observed heterozygosity was 0.52. The population structure revealed three distinct populations that clustered into two major groups after phylogenetic analysis. Analysis of molecular variance (AMOVA) indicated greater genetic variation within populations than among populations. Although cowpea generally has a narrow genetic diversity, the accessions used in this study exhibited considerable variation across geographical regions, sub-species, and improvement status. These results indicated that the selected KASP genotyping assay can provide robust and accurate genotyping data for application in the selection and management of cowpea germplasm in breeding programs and genebanks.

A high-throughput SNP marker system for parental polymorphism screening, and diversity analysis in common bean (Phaseolus vulgaris L.).

Single nucleotide polymorphism (SNP) detection has become a marker system of choice, because of the high abundance of source polymorphisms and the ease with which allele calls are automated. Various technologies exist for the evaluation of SNP loci and previously we validated two medium throughput technologies. In this study, our goal was to utilize a 768 feature, Illumina GoldenGate assay for common bean (Phaseolus vulgaris L.) developed from conserved legume gene sequences and to use the new technology for (1) the evaluation of parental polymorphisms in a mini-core set of common bean accessions and (2) the analysis of genetic diversity in the crop. A total of 736 SNPs were scored on 236 diverse common bean genotypes with the GoldenGate array. Missing data and heterozygosity levels were low and 94 % of the SNPs were scorable. With the evaluation of the parental polymorphism genotypes, we estimated the utility of the SNP markers in mapping for inter-genepool and intra-genepool populations, the latter being of lower polymorphism than the former. When we performed the diversity analysis with the diverse genotypes, we found Illumina GoldenGate SNPs to provide equivalent evaluations as previous gene-based SNP markers, but less fine-distinctions than with previous microsatellite marker analysis. We did find, however, that the gene-based SNPs in the GoldenGate array had some utility in race structure analysis despite the low polymorphism. Furthermore the SNPs detected high heterozygosity in wild accessions which was probably a reflection of ascertainment bias. The Illumina SNPs were shown to be effective in distinguishing between the genepools, and therefore were most useful in saturation of inter-genepool genetic maps. The implications of these results for breeding in common bean are discussed as well as the advantages and disadvantages of the GoldenGate system for SNP detection.

Isolation and characterization of nucleotide-binding site resistance gene homologues in common bean (Phaseolus vulgaris).

Common bean production is constrained by many fungal, viral, and bacterial pathogens. Thus, the identification of resistance (R) genes is an important focal point of common bean research. The main goal of our study was to identify resistance gene homologues (RGH) in the crop, using degenerate primers designed from conserved sequences in the nucleotide-binding site (NBS) domains of R-genes from the model legume Medicago truncatula. Total DNA of the Andean common bean genotype G19833 was used for amplification of over 500 primer combinations. Sequencing of amplicons showed that 403 cloned fragments had uninterrupted open reading frames and were considered representative of functional RGH genes. The sequences were grouped at two levels of nucleotide identity (90 and 80%) and representative sequences of each group were used for phylogenetic analyses. The RGH sequence diversity of common bean was divided into TIR and non-TIR families, each with different clusters. The TIR sequences grouped into 14 clades while non-TIR sequences grouped into seven clades. Pairwise comparisons showed purifying selection, although some sequences may have been the result of diversifying selection. Knowledge about RGH genes in common bean can allow the design of molecular markers for pyramiding of resistance genes against various pathogens.

Differentiation of Andean and Mesoamerican accessions in a proposed core collection of grain amaranths.

Grain amaranths are made up of three New World species of pseudo-cereals with C4 photosynthesis from the dicotyledonous family Amaranthaceae and the genus Amaranthus. They originate in two ecoregions of the Americas, namely, the inter-Andean valleys of South America and the volcanic axis and lowlands of Mexico and Central America. These correspond to two centers of domestications for Andean and Mesoamerican crops, with one cultivated species found in the first region and two found in the latter region. To date, no core collection has been made for the grain amaranths in the United States Department of Agriculture (USDA) germplasm system. In this study, our objective was to create a core for the 2,899 gene bank accessions with collection site data by town or farm site of which 1,090 have current geo-referencing of latitude and longitude coordinates. We constituted the core with 260 genotypes of Amaranthus, which we evaluated with 90 single-nucleotide polymorphism markers. Our goal was to distinguish between Andean and Mesoamerican gene pools of amaranths, including the cultivated species and three possible progenitor or wild relative ancestors along with two more species in an outgroup. Population structure, clustering, and discriminant analysis for principal components showed that Andean species Amaranthus caudatus and Amaranthus quitensis shared fewer alleles with Mesoamerican species Amaranthus cruentus and Amaranthus hypochondriacus, compared to each group individually. Amaranthus hybridus was a bridge species that shared alleles with both regions. Molecular markers have the advantage over morphological traits at quickly distinguishing the Andean and Mesoamerican cultivars and have the added benefit of being useful for following inter-species crosses and introgression.

Use of the advanced backcross-QTL method to transfer seed mineral accumulation nutrition traits from wild to Andean cultivated common beans.

Iron deficiency anemia and zinc deficiency are major health concerns across the world and can be addressed by biofortification breeding of higher mineral concentration in staple crops, such as common bean. Wild common beans have for the most part had higher average seed mineral concentration than cultivars of this species but have small un-commercial seeds. A logical approach for the transfer of the seed mineral trait from wild beans to cultivated beans is through the advanced backcross breeding approach. The goal of this study was to analyze a population of 138 BC(2)F(3:5) introgression lines derived from the very high iron wild genotype G10022 backcrossed into the genetic background of the commercial-type variety 'Cerinza', a large-red seeded bush bean cultivar of the Andean genepool. In addition to measuring seed mineral accumulation traits and the quantitative trait loci (QTL) controlling these traits we were interested in simultaneously testing the adaptation of the introgression lines in two replicated yield trials. We found the cross to have high polymorphism and constructed an anchored microsatellite map for the population that was 1,554-cM long and covered all 11 linkage groups of the common bean genome. Through composite interval mapping (CIM) and single point analysis (SPA), we identified associations of markers and mineral traits on b01, b06, b07, b08, b10 and b11 for seed iron concentration, and markers on b01, b04 and b10 for seed zinc concentration. The b07 and b08 QTL aligned with previous QTL for iron concentration. A large number of QTL were found for seed weight (9 with CIM and 36 with SPA analysis) and correlations between seed size and mineral content affected the identification of iron and zinc contents' QTL on many linkage groups. Segregation distortion around domestication genes made some areas difficult to introgress. However, in conclusion, the advanced backcross program produced some introgression lines with high mineral accumulation traits using a wild donor parent.

First use of microsatellite markers in a large collection of cultivated and wild accessions of tepary bean (Phaseolus acutifolius A. Gray).

Tepary bean (Phaseolus acutifolius A. Gray) is a dry-land crop species that originated in the deserts of Mexico and the south-western United States and therefore is proposed as a source of drought and salt tolerance for related species and for production in marginal rainfall areas. Few genetic tools have been developed or tested for tepary bean but microsatellites from common bean are an obvious choice for diversity analysis in the crop. The first goal of this study was to validate a set of gene-derived and non-gene simple sequence repeat or microsatellite markers from common bean in tepary bean cultivars and wild relative accessions. The second and more extensive objective of this study was to evaluate the genetic diversity and population structure of the tepary bean accessions to determine if leaf-morphology variants are valid as separate sub-groups of wild tepary beans; if P. parvifolius exist as a separate variants or species; and if cultivated tepary beans originated from one domestication event or several events. Our analysis of 140 tepary bean genotypes showed that a single domestication was likely as the cultivars were most closely related to accessions from Sinaloa and northern Mexico and that diversity was much higher in the wild genotypes compared to the cultivated ones. Other results were that P. parvifolius was classified as a separate species by population structure analysis while the variants P. acutifolius var. acutifolius and var. tenuifolius were admixed and inter-crossed. P. latifolius is not a valid species or variant of P. acutifolius but represents a group of cultivars within tepary bean. This is the first analysis of microsatellite diversity in tepary beans and has implications for breeding and conservation of this crop and its wild relatives.

Nucleotide diversity patterns at the drought-related DREB2 encoding genes in wild and cultivated common bean (Phaseolus vulgaris L.).

Common beans are an important food legume faced with a series of abiotic stresses the most severe of which is drought. The crop is interesting as a model for the analysis of gene phylogenies due to its domestication process, race structure, and origins in a group of wild common beans found along the South American Andes and the region of Mesoamerica. Meanwhile, the DREB2 transcription factors have been implicated in controlling non-ABA dependent responses to drought stress. With this in mind our objective was to study in depth the genetic diversity for two DREB2 genes as possible candidates for association with drought tolerance through a gene phylogenetic analysis. In this genetic diversity assessment, we analyzed nucleotide diversity at the two candidate genes Dreb2A and Dreb2B, in partial core collections of 104 wild and 297 cultivated common beans with a total of 401 common bean genotypes from world-wide germplasm analyzed. Our wild population sample covered a range of semi-mesic to very dry habitats, while our cultivated samples presented a wide spectrum of low to high drought tolerance. Both genes showed very different patterns of nucleotide variation. Dreb2B exhibited very low nucleotide diversity relative to neutral reference loci previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, Dreb2A exhibited higher levels of nucleotide diversity, which is indicative of adaptive selection and population expansion. These patterns were more distinct in wild compared to cultivated common beans. These approximations suggested the importance of Dreb2 genes in the context of drought tolerance, and constitute the first steps towards an association study between genetic polymorphism of this gene family and variation in drought tolerance traits. We discuss the utility of allele mining in the DREB gene family for the discovery of new drought tolerance traits from wild common bean.

Molecular ecology and selection in the drought-related Asr gene polymorphisms in wild and cultivated common bean (Phaseolus vulgaris L.).

BACKGROUND: The abscisic acid (ABA) pathway plays an important role in the plants' reaction to drought stress and ABA-stress response (Asr) genes are important in controlling this process. In this sense, we accessed nucleotide diversity at two candidate genes for drought tolerance (Asr1 and Asr2), involved in an ABA signaling pathway, in the reference collection of cultivated common bean (Phaseolus vulgaris L.) and a core collection of wild common bean accessions. RESULTS: Our wild population samples covered a range of mesic (semi-arid) to very dry (desert) habitats, while our cultivated samples presented a wide spectrum of drought tolerance. Both genes showed very different patterns of nucleotide variation. Asr1 exhibited very low nucleotide diversity relative to the neutral reference loci that were previously surveyed in these populations. This suggests that strong purifying selection has been acting on this gene. In contrast, Asr2 exhibited higher levels of nucleotide diversity, which is indicative of adaptive selection. These patterns were more notable in wild beans than in cultivated common beans indicting that natural selection has played a role over long time periods compared to farmer selection since domestication. CONCLUSIONS: Together these results suggested the importance of Asr1 in the context of drought tolerance, and constitute the first steps towards an association study between genetic polymorphism of this gene family and variation in drought tolerance traits. Furthermore, one of our major successes was to find that wild common bean is a reservoir of genetic variation and selection signatures at Asr genes, which may be useful for breeding drought tolerance in cultivated common bean.

Gene-based single nucleotide polymorphism markers for genetic and association mapping in common bean.

BACKGROUND: In common bean, expressed sequence tags (ESTs) are an underestimated source of gene-based markers such as insertion-deletions (Indels) or single-nucleotide polymorphisms (SNPs). However, due to the nature of these conserved sequences, detection of markers is difficult and portrays low levels of polymorphism. Therefore, development of intron-spanning EST-SNP markers can be a valuable resource for genetic experiments such as genetic mapping and association studies. RESULTS: In this study, a total of 313 new gene-based markers were developed at target genes. Intronic variation was deeply explored in order to capture more polymorphism. Introns were putatively identified after comparing the common bean ESTs with the soybean genome, and the primers were designed over intron-flanking regions. The intronic regions were evaluated for parental polymorphisms using the single strand conformational polymorphism (SSCP) technique and Sequenom MassARRAY system. A total of 53 new marker loci were placed on an integrated molecular map in the DOR364 × G19833 recombinant inbred line (RIL) population. The new linkage map was used to build a consensus map, merging the linkage maps of the BAT93 × JALO EEP558 and DOR364 × BAT477 populations. A total of 1,060 markers were mapped, with a total map length of 2,041 cM across 11 linkage groups. As a second application of the generated resource, a diversity panel with 93 genotypes was evaluated with 173 SNP markers using the MassARRAY-platform and KASPar technology. These results were coupled with previous SSR evaluations and drought tolerance assays carried out on the same individuals. This agglomerative dataset was examined, in order to discover marker-trait associations, using general linear model (GLM) and mixed linear model (MLM). Some significant associations with yield components were identified, and were consistent with previous findings. CONCLUSIONS: In short, this study illustrates the power of intron-based markers for linkage and association mapping in common bean. The utility of these markers is discussed in relation with the usefulness of microsatellites, the molecular markers by excellence in this crop.

Genome-Environment Associations, an Innovative Tool for Studying Heritable Evolutionary Adaptation in Orphan Crops and Wild Relatives.

Leveraging innovative tools to speed up prebreeding and discovery of genotypic sources of adaptation from landraces, crop wild relatives, and orphan crops is a key prerequisite to accelerate genetic gain of abiotic stress tolerance in annual crops such as legumes and cereals, many of which are still orphan species despite advances in major row crops. Here, we review a novel, interdisciplinary approach to combine ecological climate data with evolutionary genomics under the paradigm of a new field of study: genome-environment associations (GEAs). We first exemplify how GEA utilizes in situ georeferencing from genotypically characterized, gene bank accessions to pinpoint genomic signatures of natural selection. We later discuss the necessity to update the current GEA models to predict both regional- and local- or micro-habitat-based adaptation with mechanistic ecophysiological climate indices and cutting-edge GWAS-type genetic association models. Furthermore, to account for polygenic evolutionary adaptation, we encourage the community to start gathering genomic estimated adaptive values (GEAVs) for genomic prediction (GP) and multi-dimensional machine learning (ML) models. The latter two should ideally be weighted by de novo GWAS-based GEA estimates and optimized for a scalable marker subset. We end the review by envisioning avenues to make adaptation inferences more robust through the merging of high-resolution data sources, such as environmental remote sensing and summary statistics of the genomic site frequency spectrum, with the epigenetic molecular functionality responsible for plastic inheritance in the wild. Ultimately, we believe that coupling evolutionary adaptive predictions with innovations in ecological genomics such as GEA will help capture hidden genetic adaptations to abiotic stresses based on crop germplasm resources to assist responses to climate change. "I shall endeavor to find out how nature's forces act upon one another, and in what manner the geographic environment exerts its influence on animals and plants. In short, I must find out about the harmony in nature" Alexander von Humboldt-Letter to Karl Freiesleben, June 1799.

Genetic differentiation of grain, fodder and pod vegetable type cowpeas (Vigna unguiculata L.) identified through single nucleotide polymorphisms from genotyping-by-sequencing.

The species Vigna unguiculata L. (Walp), commonly known as cowpea, is a multi-purpose legume that has been selected into three subspecies that are divided into grain, fodder and pod (yardlong bean) types. However, genetic bases for distinctions are not well understood. The purpose of this study was to apply genotyping-by-sequencing (GBS) and current reference genome for V. unguiculata to distinguish three subspecies and identify signatures of divergence. The collection of 130 accessions included 128 cultivated from: 1) ssp. cylindrica, fodder type; 2) ssp. sesquipedalis, pod vegetable type; and 3) ssp. unguiculata, grain type. Two wilds genotypes from spp. dekindtiana and spp. pubescens, were used to anchor phylogeny. A total of 11,083 highly informative single nucleotide polymorphisms (SNPs) were discovered. Wild accessions showed distinct genetic fingerprints and were separated from cultivated subspecies. Principal component analysis showed closer relationship between ssp. unguiculata and ssp. cylindrica compared to ssp. sesquipedalis. Relative differentiation of cultivated subspecies (with Fixation Index, FST) indicated the existence of discrete signatures of selection. This work clarifies the population structure, phylogeny, and domestication of cultivated cowpeas. Furthermore, significant genetic differences between grain and pod vegetable types can provide valuable information for future breeding in three cowpea groups.

Construction and EST sequencing of full-length, drought stress cDNA libraries for common beans (Phaseolus vulgaris L.).

BACKGROUND: Common bean is an important legume crop with only a moderate number of short expressed sequence tags (ESTs) made with traditional methods. The goal of this research was to use full-length cDNA technology to develop ESTs that would overlap with the beginning of open reading frames and therefore be useful for gene annotation of genomic sequences. The library was also constructed to represent genes expressed under drought, low soil phosphorus and high soil aluminum toxicity. We also undertook comparisons of the full-length cDNA library to two previous non-full clone EST sets for common bean. RESULTS: Two full-length cDNA libraries were constructed: one for the drought tolerant Mesoamerican genotype BAT477 and the other one for the acid-soil tolerant Andean genotype G19833 which has been selected for genome sequencing. Plants were grown in three soil types using deep rooting cylinders subjected to drought and non-drought stress and tissues were collected from both roots and above ground parts. A total of 20,000 clones were selected robotically, half from each library. Then, nearly 10,000 clones from the G19833 library were sequenced with an average read length of 850 nucleotides. A total of 4,219 unigenes were identified consisting of 2,981 contigs and 1,238 singletons. These were functionally annotated with gene ontology terms and placed into KEGG pathways. Compared to other EST sequencing efforts in common bean, about half of the sequences were novel or represented the 5' ends of known genes. CONCLUSIONS: The present full-length cDNA libraries add to the technological toolbox available for common bean and our sequencing of these clones substantially increases the number of unique EST sequences available for the common bean genome. All of this should be useful for both functional gene annotation, analysis of splice site variants and intron/exon boundary determination by comparison to soybean genes or with common bean whole-genome sequences. In addition the library has a large number of transcription factors and will be interesting for discovery and validation of drought or abiotic stress related genes in common bean.

Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series.

BACKGROUND: Sequencing of cDNA libraries for the development of expressed sequence tags (ESTs) as well as for the discovery of simple sequence repeats (SSRs) has been a common method of developing microsatellites or SSR-based markers. In this research, our objective was to further sequence and develop common bean microsatellites from leaf and root cDNA libraries derived from the Andean gene pool accession G19833 and the Mesoamerican gene pool accession DOR364, mapping parents of a commonly used reference map. The root libraries were made from high and low phosphorus treated plants. RESULTS: A total of 3,123 EST sequences from leaf and root cDNA libraries were screened and used for direct simple sequence repeat discovery. From these EST sequences we found 184 microsatellites; the majority containing tri-nucleotide motifs, many of which were GC rich (ACC, AGC and AGG in particular). Di-nucleotide motif microsatellites were about half as common as the tri-nucleotide motif microsatellites but most of these were AGn microsatellites with a moderate number of ATn microsatellites in root ESTs followed by few ACn and no GCn microsatellites. Out of the 184 new SSR loci, 120 new microsatellite markers were developed in the BMc (Bean Microsatellites from cDNAs) series and these were evaluated for their capacity to distinguish bean diversity in a germplasm panel of 18 genotypes. We developed a database with images of the microsatellites and their polymorphism information content (PIC), which averaged 0.310 for polymorphic markers. CONCLUSIONS: The present study produced information about microsatellite frequency in root and leaf tissues of two important genotypes for common bean genomics: namely G19833, the Andean genotype selected for whole genome shotgun sequencing from race Peru, and DOR364 a race Mesoamerica subgroup 2 genotype that is a small-red seeded, released variety in Central America. Both race Peru and Mesoamerica subgroup 2 (small red beans) have been understudied in comparison to race Nueva Granada and Mesoamerica subgroup 1 (black beans) both with regards to gene expression and as sources of markers. However, we found few differences between SSR type and frequency between the G19833 leaf and DOR364 root tissue-derived ESTs. Overall, our work adds to the analysis of microsatellite frequency evaluation for common bean and provides a new set of 120 BMc markers which combined with the 248 previously developed BMc markers brings the total in this series to 368 markers. Once we include BMd markers, which are derived from GenBank sequences, the current total of gene-based markers from our laboratory surpasses 500 markers. These markers are basic for studies of the transcriptome of common bean and can form anchor points for genetic mapping studies in the future.

SNP marker diversity in common bean (Phaseolus vulgaris L.).

Single nucleotide polymorphism (SNP) markers have become a genetic technology of choice because of their automation and high precision of allele calls. In this study, our goal was to develop 94 SNPs and test them across well-chosen common bean (Phaseolus vulgaris L.) germplasm. We validated and accessed SNP diversity at 84 gene-based and 10 non-genic loci using KASPar technology in a panel of 70 genotypes that have been used as parents of mapping populations and have been previously evaluated for SSRs. SNPs exhibited high levels of genetic diversity, an excess of middle frequency polymorphism, and a within-genepool mismatch distribution as expected for populations affected by sudden demographic expansions after domestication bottlenecks. This set of markers was useful for distinguishing Andean and Mesoamerican genotypes but less useful for distinguishing within each gene pool. In summary, slightly greater polymorphism and race structure was found within the Andean gene pool than within the Mesoamerican gene pool but polymorphism rate between genotypes was consistent with genepool and race identity. Our survey results represent a baseline for the choice of SNP markers for future applications because gene-associated SNPs could themselves be causative SNPs for traits. Finally, we discuss that the ideal genetic marker combination with which to carry out diversity, mapping and association studies in common bean should consider a mix of both SNP and SSR markers.

QTL analyses for seed iron and zinc concentrations in an intra-genepool population of Andean common beans (Phaseolus vulgaris L.).

Legumes provide essential micronutrients that are found only in low amounts in the cereals or root crops. An ongoing project at CIAT has shown that the legume common bean is variable in the amount of seed minerals (iron, zinc, and other elements), vitamins, and sulfur amino acids that they contain and that these traits are likely to be inherited quantitatively. In this study we analyzed iron and zinc concentrations in an Andean recombinant inbred line (RIL) population of 100 lines derived from a cross between G21242, a Colombian cream-mottled climbing bean with high seed iron/zinc and G21078, an Argentinean cream seeded climbing bean with low seed iron/zinc. The population was planted across three environments; seed from each genotype was analyzed with two analytical methods, and quantitative trait loci (QTL) were detected using composite interval mapping and single-point analyses. A complete genetic map was created for the cross using a total of 74 microsatellite markers to anchor the map to previously published reference maps and 42 RAPD markers. In total, nine seed mineral QTL were identified on five linkage groups (LGs) with the most important being new loci on b02 and other QTL on b06, b08, and b07 near phaseolin. Seed weight QTL were associated with these on b02 and b08. These Andean-derived QTL are candidates for marker-assisted selection either in combination with QTL from the Mesoamerican genepool or with other QTL found in inter and intra-genepool crosses, and the genetic map can be used to anchor other intra-genepool studies.

Biofortified red mottled beans (Phaseolus vulgaris L.) in a maize and bean diet provide more bioavailable iron than standard red mottled beans: studies in poultry (Gallus gallus) and an in vitro digestion/Caco-2 model.

BACKGROUND: Our objective was to compare the capacities of biofortified and standard colored beans to deliver iron (Fe) for hemoglobin synthesis. Two isolines of large-seeded, red mottled Andean beans (Phaseolus vulgaris L.), one standard ("Low Fe") and the other biofortified ("High Fe") in Fe (49 and 71 µg Fe/g, respectively) were used. This commercial class of red mottled beans is the preferred varietal type for most of the Caribbean and Eastern and Southern Africa where almost three quarters of a million hectares are grown. Therefore it is important to know the affect of biofortification of these beans on diets that simulate human feeding studies. METHODS: Maize-based diets containing the beans were formulated to meet the nutrient requirements for broiler except for Fe (Fe concentrations in the 2 diets were 42.9 ± 1.2 and 54.6 ± 0.9 mg/kg). One day old chicks (Gallus gallus) were allocated to the experimental diets (n = 12). For 4 wk, hemoglobin, feed-consumption and body-weights were measured. RESULTS: Hemoglobin maintenance efficiencies (HME) (means ± SEM) were different between groups on days 14 and 21 of the experiment (P < 0.05). Final total body hemoglobin Fe contents were different between the standard (12.58 ± 1.0 mg {0.228 ± 0.01 µmol}) and high Fe (15.04 ± 0.65 mg {0.273 ± 0.01 µmol}) bean groups (P < 0.05). At the end of the experiment, tissue samples were collected from the intestinal duodenum and liver for further analyses. Divalent-metal-transporter-1, duodenal-cytochrome-B, and ferroportin expressions were higher and liver ferritin was lower (P < 0.05) in the standard group vs. the biofortified group. In-vitro analysis showed lower iron bioavailability in cells exposed to standard ("Low Fe") bean based diet. CONCLUSIONS: We conclude that the in-vivo results support the in-vitro observations; biofortified colored beans contain more bioavailable-iron than standard colored beans. In addition, biofortified beans seems to be a promising vehicle for increasing intakes of bioavailable Fe in human populations that consume these beans as a dietary staple. This justifies further work on the large-seeded Andean beans which are the staple of a large-region of Africa where iron-deficiency anemia is a primary cause of infant death and poor health status.

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.

Relationship of Cultivated Grain Amaranth Species and Wild Relative Accessions.

Amaranthus is a genus of C4 dicotyledonous herbaceous plants, and three New World species have been domesticated to produce grain crops with light colored seed which are classified as pseudo-cereals rich in protein and minerals. A core collection of grain amaranths and immediate precursor species has been established, representing the closest related species. The goal of this study was to evaluate the genetic diversity in that collection of cultivated and wild species, using competitive allele single nucleotide polymorphism markers. A secondary objective was to determine the relationships among the three cultivated species and non-domesticated Amaranthus, while a third objective was to evaluate the utility of the markers in detecting diversity in the 276 genotypes. The markers were found to be highly variable with an average polymorphism information content of 0.365. All markers were bi-allelic; and the major allele frequency ranged from 0.388 to 0.871. Population structure analysis of the cultigens revealed the presence of two sub populations. Phylogeny confirmed that the two Mesoamerican species, Amaranthus cruentus and Amaranthus hypochondriacus, were related and distant from the South American species Amaranthus caudatus, which in turn was very closely clustered with Amaranthus quitensis, even though this is considered a weedy relative. The first pair of species were likely to have inter-crossed, while the latter two likely exist in a wild-cultivated hybrid state. In conclusion, the results of this SNP study provided insights on amaranth cultivars and their relationship to wild species, the probable domestication events leading to the cultivars, and possible crop breeding or germplasm conservation strategies.