Genetic variability and genome-wide association analysis of flavor and texture in cooked beans (Phaseolus vulgaris L.).

KEY MESSAGE: Cooked bean flavor and texture vary within and across 20 Andean seed types; SNPs are significantly associated with total flavor, beany, earthy, starchy, bitter, seed-coat perception, and cotyledon texture. Common dry beans are a nutritious food recognized as a staple globally, but their consumption is low in the USA. Improving bean flavor and texture through breeding has the potential to improve consumer acceptance and suitability for new end-use products. Little is known about genetic variability and inheritance of bean sensory characteristics. A total of 430 genotypes of the Andean Diversity Panel representing twenty seed types were grown in three locations, and cooked seeds were evaluated by a trained sensory panel for flavor and texture attribute intensities, including total flavor, beany, vegetative, earthy, starchy, sweet, bitter, seed-coat perception, and cotyledon texture. Extensive variation in sensory attributes was found across and within seed types. A set of genotypes was identified that exhibit extreme attribute intensities generally stable across all three environments. seed-coat perception and total flavor intensity had the highest broad-sense heritability (0.39 and 0.38, respectively), while earthy and vegetative intensities exhibited the lowest (0.14 and 0.15, respectively). Starchy and sweet flavors were positively correlated and highest in white bean genotypes according to principal component analysis. SNPs associated with total flavor intensity (six SNPs across three chromosomes), beany (five SNPs across four chromosomes), earthy (three SNPs across two chromosomes), starchy (one SNP), bitter (one SNP), seed-coat perception (three SNPs across two chromosomes), and cotyledon texture (two SNPs across two chromosomes) were detected. These findings lay a foundation for incorporating flavor and texture in breeding programs for the development of new varieties that entice growers, consumers, and product developers alike.

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.

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

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

Uneven recombination rate and linkage disequilibrium across a reference SNP map for common bean (Phaseolus vulgaris L.).

Recombination (R) rate and linkage disequilibrium (LD) analyses are the basis for plant breeding. These vary by breeding system, by generation of inbreeding or outcrossing and by region in the chromosome. Common bean (Phaseolus vulgaris L.) is a favored food legume with a small sequenced genome (514 Mb) and n = 11 chromosomes. The goal of this study was to describe R and LD in the common bean genome using a 768-marker array of single nucleotide polymorphisms (SNP) based on Trans-legume Orthologous Group (TOG) genes along with an advanced-generation Recombinant Inbred Line reference mapping population (BAT93 x Jalo EEP558) and an internationally available diversity panel. A whole genome genetic map was created that covered all eleven linkage groups (LG). The LGs were linked to the physical map by sequence data of the TOGs compared to each chromosome sequence of common bean. The genetic map length in total was smaller than for previous maps reflecting the precision of allele calling and mapping with SNP technology as well as the use of gene-based markers. A total of 91.4% of TOG markers had singleton hits with annotated Pv genes and all mapped outside of regions of resistance gene clusters. LD levels were found to be stronger within the Mesoamerican genepool and decay more rapidly within the Andean genepool. The recombination rate across the genome was 2.13 cM / Mb but R was found to be highly repressed around centromeres and frequent outside peri-centromeric regions. These results have important implications for association and genetic mapping or crop improvement in common bean.

Genome sequence data from 17 accessions of Ensete ventricosum, a staple food crop for millions in Ethiopia.

We present raw sequence reads and genome assemblies derived from 17 accessions of the Ethiopian orphan crop plant enset (Ensete ventricosum (Welw.) Cheesman) using the Illumina HiSeq and MiSeq platforms. Also presented is a catalogue of single-nucleotide polymorphisms inferred from the sequence data at an average density of approximately one per kilobase of genomic DNA.

Trait Associations in Diversity Panels of the Two Common Bean (Phaseolus vulgaris L.) Gene Pools Grown under Well-watered and Water-Stress Conditions.

Common beans are a warm-season, food legume cultivated in areas prone to water limitation throughout their growing season. This study assessed the magnitude and pattern of trait associations for a total of 202 common bean genotypes divided into panels of 81 Andean and 121 Mesoamerican gene pool accessions grown under contrasting treatments of well-watered, non-stress, and water-limited, terminal drought-stress conditions. Linear correlation, complex path coefficient, and genetic divergence analyses were used to dissect the relationship dynamics between traits and the relative contribution of adaptive traits to differentiation among gene pools and genotypes based on drought stress. Drought severity level for the trial was high and created the ideal condition to reveal genotypic differences, as seen by the differential response of the genotypes for the various traits measured. The value for phenotypic coefficients of variation for all traits was higher than the corresponding genotypic values. Seed yield had positive and strong genotypic and phenotypic correlation with pods per plant across gene pools and stress levels. The overall amount of genetic correlation was greater than the corresponding phenotypic correlation matrix for all the traits within the gene pool and across stress levels. Moreover, the results depicted the phenotypic correlation as equal or better than its genotypic counterpart in estimating drought tolerance in common bean plants. Clustering analysis with Mahanalobis's coefficient of generalized distance grouped genotypes with a differential level of drought adaptation into different classes within each panel. This indicates drought tolerance involves different mechanisms of plant response and is present separately in each gene pool panel. Pods per plant, seed weight, pod partitioning index, and harvest index are useful selection objectives to improve drought adaptation in common bean, but must be differentially weighted in each gene pool. The analysis of genetic variation and association between adaptive traits on the two panels provided useful insights on which traits could be used to improve common bean adaptation to low water availability during the growth season.

Morphological Variation and Inter-Relationships of Quantitative Traits in Enset (Ensete ventricosum (welw.) Cheesman) Germplasm from South and South-Western Ethiopia.

Enset (Ensete ventricosum (Welw.) Cheesman) is Ethiopia's most important root crop. A total of 387 accessions collected from nine different regions of Ethiopia were evaluated for 15 quantitative traits at Areka Agricultural Research Centre to determine the extent and pattern of distribution of morphological variation. The variations among the accessions and regions were significant (p ≤ 0.01) for all the 15 traits studied. Mean for plant height, central shoot weight before grating, and fermented squeezed kocho yield per hectare per year showed regional variation along an altitude gradient and across cultural differences related to the origin of the collection. Furthermore, there were significant correlations among most of the characters. This included the correlation among agronomic characteristics of primary interest in enset breeding such as plant height, pseudostem height, and fermented squeezed kocho yield per hectare per year. Altitude of the collection sites also significantly impacted the various characteristics studied. These results reveal the existence of significant phenotypic variations among the 387 accessions as a whole. Regional differentiations were also evident among the accessions. The implication of the current results for plant breeding, germplasm collection, and in situ and ex situ genetic resource conservation are discussed.