Linking photosynthesis and yield reveals a strategy to improve light use efficiency in a climbing bean breeding population.

Photosynthesis drives plant physiology, biomass accumulation, and yield. Photosynthetic efficiency, specifically the operating efficiency of PSII (Fq'/Fm'), is highly responsive to actual growth conditions, especially to fluctuating photosynthetic photon fluence rate (PPFR). Under field conditions, plants constantly balance energy uptake to optimize growth. The dynamic regulation complicates the quantification of cumulative photochemical energy uptake based on the intercepted solar energy, its transduction into biomass, and the identification of efficient breeding lines. Here, we show significant effects on biomass related to genetic variation in photosynthetic efficiency of 178 climbing bean (Phaseolus vulgaris L.) lines. Under fluctuating conditions, the Fq'/Fm' was monitored throughout the growing period using hand-held and automated chlorophyll fluorescence phenotyping. The seasonal response of Fq'/Fm' to PPFR (ResponseG:PPFR) achieved significant correlations with biomass and yield, ranging from 0.33 to 0.35 and from 0.22 to 0.31 in two glasshouse and three field trials, respectively. Phenomic yield prediction outperformed genomic predictions for new environments in four trials under different growing conditions. Investigating genetic control over photosynthesis, one single nucleotide polymorphism (Chr09_37766289_13052) on chromosome 9 was significantly associated with ResponseG:PPFR in proximity to a candidate gene controlling chloroplast thylakoid formation. In conclusion, photosynthetic screening facilitates and accelerates selection for high yield potential.

Interspecific common bean population derived from Phaseolus acutifolius using a bridging genotype demonstrate useful adaptation to heat tolerance.

Common bean (Phaseolus vulgaris L.) is an important legume crop worldwide and is a major nutrient source in the tropics. Common bean reproductive development is strongly affected by heat stress, particularly overnight temperatures above 20°C. The desert Tepary bean (Phaseolus acutifolius A. Gray) offers a promising source of adaptative genes due to its natural acclimation to arid conditions. Hybridization between both species is challenging, requiring in vitro embryo rescue and multiple backcrossing cycles to restore fertility. This labor-intensive process constrains developing mapping populations necessary for studying heat tolerance. Here we show the development of an interspecific mapping population using a novel technique based on a bridging genotype derived from P. vulgaris, P. Acutifolius and P. parvifolius named VAP1 and is compatible with both common and tepary bean. The population was based on two wild P. acutifolius accessions, repeatedly crossed with Mesoamerican elite common bush bean breeding lines. The population was genotyped through genotyping-by-sequencing and evaluated for heat tolerance by genome-wide association studies. We found that the population harbored 59.8% introgressions from wild tepary, but also genetic regions from Phaseolus parvifolius, a relative represented in some early bridging crosses. We found 27 significative quantitative trait loci, nine located inside tepary introgressed segments exhibiting allelic effects that reduced seed weight, and increased the number of empty pods, seeds per pod, stem production and yield under high temperature conditions. Our results demonstrate that the bridging genotype VAP1 can intercross common bean with tepary bean and positively influence the physiology of derived interspecific lines, which displayed useful variance for heat tolerance.

Genetic analysis of resistance to bean leaf crumple virus identifies a candidate LRR-RLK gene.

Bean leaf crumple virus (BLCrV) is a novel begomovirus (family Geminiviridae, genus Begomovirus) infecting common bean (Phaseolus vulgaris L.), threatening bean production in Latin America. Genetic resistance is required to ensure yield stability and reduce the use of insecticides, yet the available resistance sources are limited. In this study, three common bean populations containing a total of 558 genotypes were evaluated in different yield and BLCrV resistance trials under natural infection in the field. A genome-wide association study identified the locus BLC7.1 on chromosome Pv07 at 3.31 Mbp, explaining 8 to 16% of the phenotypic variation for BLCrV resistance. In comparison, whole-genome regression models explained 51 to 78% of the variation and identified the same region on Pv07 to confer resistance. The most significantly associated markers were located within the gene model Phvul.007G040400, which encodes a leucine-rich repeat receptor-like kinase subfamily III member and is likely to be involved in the innate immune response against the virus. The allelic diversity within this gene revealed five different haplotype groups, one of which was significantly associated with BLCrV resistance. As the same genome region was previously reported to be associated with resistance against other geminiviruses affecting common bean, our study highlights the role of previous breeding efforts for virus resistance in the accumulation of positive alleles against newly emerging viruses. In addition, we provide novel diagnostic single-nucleotide polymorphism markers for marker-assisted selection to exploit BLC7.1 for breeding against geminivirus diseases in one of the most important food crops worldwide.

Mr.Bean: a comprehensive statistical and visualization application for modeling agricultural field trials data.

Crop improvement efforts have exploited new methods for modeling spatial trends using the arrangement of the experimental units in the field. These methods have shown improvement in predicting the genetic potential of evaluated genotypes. However, the use of these tools may be limited by the exposure and accessibility to these products. In addition, these new methodologies often require plant scientists to be familiar with the programming environment used to implement them; constraints that limit data analysis efficiency for decision-making. These challenges have led to the development of Mr.Bean, an accessible and user-friendly tool with a comprehensive graphical visualization interface. The application integrates descriptive analysis, measures of dispersion and centralization, linear mixed model fitting, multi-environment trial analysis, factor analytic models, and genomic analysis. All these capabilities are designed to help plant breeders and scientist working with agricultural field trials make informed decisions more quickly. Mr.Bean is available for download at https://github.com/AparicioJohan/MrBeanApp.

Improving Association Studies and Genomic Predictions for Climbing Beans With Data From Bush Bean Populations.

Common bean (Phaseolus vulgaris L.) has two major origins of domestication, Andean and Mesoamerican, which contribute to the high diversity of growth type, pod and seed characteristics. The climbing growth habit is associated with increased days to flowering (DF), seed iron concentration (SdFe), nitrogen fixation, and yield. However, breeding efforts in climbing beans have been limited and independent from bush type beans. To advance climbing bean breeding, we carried out genome-wide association studies and genomic predictions using 1,869 common bean lines belonging to five breeding panels representing both gene pools and all growth types. The phenotypic data were collected from 17 field trials and were complemented with 16 previously published trials. Overall, 38 significant marker-trait associations were identified for growth habit, 14 for DF, 13 for 100 seed weight, three for SdFe, and one for yield. Except for DF, the results suggest a common genetic basis for traits across all panels and growth types. Seven QTL associated with growth habits were confirmed from earlier studies and four plausible candidate genes for SdFe and 100 seed weight were newly identified. Furthermore, the genomic prediction accuracy for SdFe and yield in climbing beans improved up to 8.8% when bush-type bean lines were included in the training population. In conclusion, a large population from different gene pools and growth types across multiple breeding panels increased the power of genomic analyses and provides a solid and diverse germplasm base for genetic improvement of common bean.

Multivariate genomic analysis and optimal contributions selection predicts high genetic gains in cooking time, iron, zinc, and grain yield in common beans in East Africa.

Common bean (Phaseolus vulgaris L.) is important in African diets for protein, iron (Fe), and zinc (Zn), but traditional cultivars have long cooking time (CKT), which increases the time, energy, and health costs of cooking. Genomic selection was used to predict genomic estimated breeding values (GEBV) for grain yield (GY), CKT, Fe, and Zn in an African bean panel of 358 genotypes in a two-stage analysis. In Stage 1, best linear unbiased estimates (BLUE) for each trait were obtained from 898 genotypes across 33 field trials in East Africa. In Stage 2, BLUE in a training population of 141 genotypes were used in a multivariate genomic analysis with genome-wide single nucleotide polymorphism data from the African bean panel. Moderate to high genomic heritability was found for GY (0.45 ± 0.10), CKT (0.50 ± 0.15), Fe (0.57 ± 0.12), and Zn (0.61 ± 0.13). There were significant favorable genetic correlations between Fe and Zn (0.91 ± 0.06), GY and Fe (0.66 ± 0.17), GY and Zn (0.44 ± 0.19), CKT and Fe (-0.57 ± 0.21), and CKT and Zn (-0.67 ± 0.20). Optimal contributions selection (OCS), based on economic index of weighted GEBV for each trait, was used to design crossing within four market groups relevant to East Africa. Progeny were predicted by OCS to increase in mean GY by 12.4%, decrease in mean CKT by 9.3%, and increase in mean Fe and Zn content by 6.9 and 4.6%, respectively, with low achieved coancestry of 0.032. Genomic selection with OCS will accelerate breeding of high-yielding, biofortified, and rapid cooking African common bean cultivars.

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

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

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

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

Genetic Analyses and Genomic Predictions of Root Rot Resistance in Common Bean Across Trials and Populations.

Root rot in common bean is a disease that causes serious damage to grain production, particularly in the upland areas of Eastern and Central Africa where significant losses occur in susceptible bean varieties. Pythium spp. and Fusarium spp. are among the soil pathogens causing the disease. In this study, a panel of 228 lines, named RR for root rot disease, was developed and evaluated in the greenhouse for Pythium myriotylum and in a root rot naturally infected field trial for plant vigor, number of plants germinated, and seed weight. The results showed positive and significant correlations between greenhouse and field evaluations, as well as high heritability (0.71-0.94) of evaluated traits. In GWAS analysis no consistent significant marker trait associations for root rot disease traits were observed, indicating the absence of major resistance genes. However, genomic prediction accuracy was found to be high for Pythium, plant vigor and related traits. In addition, good predictions of field phenotypes were obtained using the greenhouse derived data as a training population and vice versa. Genomic predictions were evaluated across and within further published data sets on root rots in other panels. Pythium and Fusarium evaluations carried out in Uganda on the Andean Diversity Panel showed good predictive ability for the root rot response in the RR panel. Genomic prediction is shown to be a promising method to estimate tolerance to Pythium, Fusarium and root rot related traits, indicating a quantitative resistance mechanism. Quantitative analyses could be applied to other disease-related traits to capture more genetic diversity with genetic models.

Introgression of bruchid (Zabrotes subfasciatus) resistance into small red common bean (Phaseolus vulgaris) background and validation of the BRU_00261 (snpPV0007) resistance marker.

Bruchids are a major storage pest of common bean. Genetic resistance is a suitable method to avoid grain losses during storage. The objective of the study was to introgress the arcelin-based resistance locus into selected advanced breeding line and to validate the molecular marker BRU_00261. A total of 208 progeny F4 families were phenotyped using a randomized complete block design, with three replications. Highly significant differences (P < .001) among the entries, parents and offspring were recorded for almost all traits. There was no significant difference between the two parents in the number of eggs laid. The progenies were grouped as highly resistant (34.3%), resistant (11.9%), moderately resistant (21.4%) and susceptible (32.4%). The levels of broad sense heritability ranged from 68.5%-93.9% for all the traits. Eighty-three most resistant lines and the parental lines were genotyped with the marker BRU_00261 (snpPV0007). The marker segregation deviated significantly from the expected independent segregation towards a strong enrichment for the resistant marker in the selected families. This marker will be useful for selecting promising materials in early generations and phenotypic confirmation of positive lines in later generations.

Genetic Correlation Between Fe and Zn Biofortification and Yield Components in a Common Bean (Phaseolus vulgaris L.).

Common bean (Phaseolus vulgaris L.) is the most important legume for direct human consumption worldwide. It is a rich and relatively inexpensive source of proteins and micronutrients, especially iron and zinc. Bean is a target for biofortification to develop new cultivars with high Fe/Zn levels that help to ameliorate malnutrition mainly in developing countries. A strong negative phenotypic correlation between Fe/Zn concentration and yield is usually reported, posing a significant challenge for breeders. The objective of this study was to investigate the genetic relationship between Fe/Zn. We used Quantitative Trait Loci (QTLs) mapping and Genome-Wide Association Studies (GWAS) analysis in three bi-parental populations that included biofortified parents, identifying genomic regions associated with yield and micromineral accumulation. Significant negative correlations were observed between agronomic traits (pod harvest index, PHI; pod number, PdN; seed number, SdN; 100 seed weight, 100SdW; and seed per pod, Sd/Pd) and micronutrient concentration traits (SdFe and SdZn), especially between pod harvest index (PHI) and SdFe and SdZn. PHI presented a higher correlation with SdN than PdN. Seventy-nine QTLs were identified for the three populations: 14 for SdFe, 12 for SdZn, 13 for PHI, 11 for SdN, 14 for PdN, 6 for 100SdW, and 9 for Sd/Pd. Twenty-three hotspot regions were identified in which several QTLs were co-located, of which 13 hotpots displayed QTL of opposite effect for yield components and Fe/Zn accumulation. In contrast, eight QTLs for SdFe and six QTLs for SdZn were observed that segregated independently of QTL of yield components. The selection of these QTLs will enable enhanced levels of Fe/Zn and will not affect the yield performance of new cultivars focused on biofortification.

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

BACKGROUND: Common bean is an important staple crop in the tropics of Africa, Asia and the Americas. Particularly smallholder farmers rely on bean as a source for calories, protein and micronutrients. Drought is a major production constraint for common bean, a situation that will be aggravated with current climate change scenarios. In this context, new tools designed to understand the genetic basis governing the phenotypic responses to abiotic stress are required to improve transfer of desirable traits into cultivated beans. RESULTS: A multiparent advanced generation intercross (MAGIC) population of common bean was generated from eight Mesoamerican breeding lines representing the phenotypic and genotypic diversity of the CIAT Mesoamerican breeding program. This population was assessed under drought conditions in two field trials for yield, 100 seed weight, iron and zinc accumulation, phenology and pod harvest index. Transgressive segregation was observed for most of these traits. Yield was positively correlated with yield components and pod harvest index (PHI), and negative correlations were found with phenology traits and micromineral contents. Founder haplotypes in the population were identified using Genotyping by Sequencing (GBS). No major population structure was observed in the population. Whole Genome Sequencing (WGS) data from the founder lines was used to impute genotyping data for GWAS. Genetic mapping was carried out with two methods, using association mapping with GWAS, and linkage mapping with haplotype-based interval screening. Thirteen high confidence QTL were identified using both methods and several QTL hotspots were found controlling multiple traits. A major QTL hotspot located on chromosome Pv01 for phenology traits and yield was identified. Further hotspots affecting several traits were observed on chromosomes Pv03 and Pv08. A major QTL for seed Fe content was contributed by MIB778, the founder line with highest micromineral accumulation. Based on imputed WGS data, candidate genes are reported for the identified major QTL, and sequence changes were identified that could cause the phenotypic variation. CONCLUSIONS: This work demonstrates the importance of this common bean MAGIC population for genetic mapping of agronomic traits, to identify trait associations for molecular breeding tool design and as a new genetic resource for the bean research community.

Biosynthesis of cyanogenic glucosides in Phaseolus lunatus and the evolution of oxime-based defenses.

Lima bean, Phaseolus lunatus, is a crop legume that produces the cyanogenic glucosides linamarin and lotaustralin. In the legumes Lotus japonicus and Trifolium repens, the biosynthesis of these two α-hydroxynitrile glucosides involves cytochrome P450 enzymes of the CYP79 and CYP736 families and a UDP-glucosyltransferase. Here, we identify CYP79D71 as the first enzyme of the pathway in P. lunatus, producing oximes from valine and isoleucine. A second CYP79 family member, CYP79D72, was shown to catalyze the formation of leucine-derived oximes, which act as volatile defense compounds in Phaseolus spp. The organization of the biosynthetic genes for cyanogenic glucosides in a gene cluster aided their identification in L. japonicus. In the available genome sequence of P. vulgaris, the gene orthologous to CYP79D71 is adjacent to a member of the CYP83 family. Although P. vulgaris is not cyanogenic, it does produce oximes as volatile defense compounds. We cloned the genes encoding two CYP83s (CYP83E46 and CYP83E47) and a UDP-glucosyltransferase (UGT85K31) from P. lunatus, and these genes combined form a complete biosynthetic pathway for linamarin and lotaustralin in Lima bean. Within the genus Phaseolus, the occurrence of linamarin and lotaustralin as functional chemical defense compounds appears restricted to species belonging to the closely related Polystachios and Lunatus groups. A preexisting ability to produce volatile oximes and nitriles likely facilitated evolution of cyanogenesis within the Phaseolus genus.

Genomic Prediction of Agronomic Traits in Common Bean (Phaseolus vulgaris L.) Under Environmental Stress.

In plant and animal breeding, genomic prediction models are established to select new lines based on genomic data, without the need for laborious phenotyping. Prediction models can be trained on recent or historic phenotypic data and increasingly available genotypic data. This enables the adoption of genomic selection also in under-used legume crops such as common bean. Beans are an important staple food in the tropics and mainly grown by smallholders under limiting environmental conditions such as drought or low soil fertility. Therefore, genotype-by-environment interactions (G × E) are an important consideration when developing new bean varieties. However, G × E are often not considered in genomic prediction models nor are these models implemented in current bean breeding programs. Here we show the prediction abilities of four agronomic traits in common bean under various environmental stresses based on twelve field trials. The dataset includes 481 elite breeding lines characterized by 5,820 SNP markers. Prediction abilities over all twelve trials ranged between 0.6 and 0.8 for yield and days to maturity, respectively, predicting new lines into new seasons. In all four evaluated traits, the prediction abilities reached about 50-80% of the maximum accuracies given by phenotypic correlations and heritability. Predictions under drought and low phosphorus stress were up to 10 and 20% improved when G × E were included in the model, respectively. Our results demonstrate the potential of genomic selection to increase the genetic gain in common bean breeding. Prediction abilities improved when more phenotypic data was available and G × E could be accounted for. Furthermore, the developed models allowed us to predict genotypic performance under different environmental stresses. This will be a key factor in the development of common bean varieties adapted to future challenging conditions.

Iron, Zinc and Phytic Acid Retention of Biofortified, Low Phytic Acid, and Conventional Bean Varieties When Preparing Common Household Recipes.

Biofortification is an effective method to improve the nutritional content of crops and nutritional intake. Breeding for higher micronutrient mineral content in beans is correlated with an increase in phytic acid, a main inhibitor of mineral absorption in humans. Low phytic acid (lpa) beans have a 90% lower phytic acid content compared to conventional beans. This is the first study to investigate mineral and total phytic acid retention after preparing common household recipes from conventional, biofortified and lpa beans. Mineral retention was determined for two conventional, three biofortified and two lpa bean genotypes. Treatments included soaking, boiling (boiled beans) and refrying (bean paste). The average true retention of iron after boiling was 77.2-91.3%; for zinc 41.2-84.0%; and for phytic acid 49.9-85.9%. Soaking led to a significant decrease in zinc and total phytic acid after boiling and refrying, whereas for iron no significant differences were found. lpa beans did not exhibit a consistent pattern of difference in iron and phytic acid retention compared to the other groups of beans. However, lpa beans had a significantly lower retention of zinc compared to conventional and biofortified varieties (p < 0.05). More research is needed to understand the underlying factors responsible for the differences in retention between the groups of beans, especially the low retention of zinc. Combining the lpa and biofortification traits could further improve the nutritional benefits of biofortified beans, by decreasing the phytic acid:iron and zinc ratio in beans.

Genetic Architecture and Genomic Prediction of Cooking Time in Common Bean (Phaseolus vulgaris L.).

Cooking time of the common bean is an important trait for consumer preference, with implications for nutrition, health, and environment. For efficient germplasm improvement, breeders need more information on the genetics to identify fast cooking sources with good agronomic properties and molecular breeding tools. In this study, we investigated a broad genetic variation among tropical germplasm from both Andean and Mesoamerican genepools. Four populations were evaluated for cooking time (CKT), water absorption capacity (WAC), and seed weight (SdW): a bi-parental RIL population (DxG), an eight-parental Mesoamerican MAGIC population, an Andean (VEF), and a Mesoamerican (MIP) breeding line panel. A total of 922 lines were evaluated in this study. Significant genetic variation was found in all populations with high heritabilities, ranging from 0.64 to 0.89 for CKT. CKT was related to the color of the seed coat, with the white colored seeds being the ones that cooked the fastest. Marker trait associations were investigated by QTL analysis and GWAS, resulting in the identification of 10 QTL. In populations with Andean germplasm, an inverse correlation of CKT and WAC, and also a QTL on Pv03 that inversely controls CKT and WAC (CKT3.2/WAC3.1) were observed. WAC7.1 was found in both Mesoamerican populations. QTL only explained a small part of the variance, and phenotypic distributions support a more quantitative mode of inheritance. For this reason, we evaluated how genomic prediction (GP) models can capture the genetic variation. GP accuracies for CKT varied, ranging from good results for the MAGIC population (0.55) to lower accuracies in the MIP panel (0.22). The phenotypic characterization of parental material will allow for the cooking time trait to be implemented in the active germplasm improvement programs. Molecular breeding tools can be developed to employ marker-assisted selection or genomic selection, which looks to be a promising tool in some populations to increase the efficiency of breeding activities.

Genomics, genetics and breeding of common bean in Africa: A review of tropical legume project.

Common bean (Phaseolus vulgaris L.) is an important legume crop worldwide. The International Centre for Tropical Agriculture (CIAT) and its national partners in Africa aim to overcome production constraints of common bean and address the food, nutrition needs and market demands through development of multitrait bean varieties. Breeding is guided by principles of market-driven approaches to develop client-demanded varieties. Germplasm accessions from especially two sister species, P. coccineus and P. acutifolius, have been utilized as sources of resistance to major production constraints and interspecific lines deployed. Elucidation of plant mechanisms governing pest and disease resistance, abiotic stress tolerance and grain nutritional quality guides the selection methods used by the breeders. Molecular markers are used to select for resistance to key diseases and insect pests. Efforts have been made to utilize modern genomic tools to increase scale, efficiency, accuracy and speed of breeding. Through gender-responsive participatory variety selection, market-demanded varieties have been released in several African countries. These new bean varieties are a key component of sustainable food systems in the tropics.

Haplotypes at the Phg-2 Locus Are Determining Pathotype-Specificity of Angular Leaf Spot Resistance in Common Bean.

Angular leaf spot (ALS) is one of the most devastating diseases of common bean (Phaseolus vulgaris L.) and causes serious yield losses worldwide. ALS resistance is reportedly pathotype-specific, but little is known about the efficacy of resistance loci against different pathotypes. Here, we report on ALS resistance evaluations of 316 bean lines under greenhouse and field conditions at multiple sites in Colombia and Uganda. Surprisingly, genome-wide association studies revealed only two of the five previously described resistance loci to be significantly associated with ALS resistance. Phg-2 on chromosome eight was crucial for ALS resistance in all trials, while the resistance locus Phg-4 on chromosome 4 was effective against one particular pathotype. Further dissection of Phg-2 uncovered an unprecedented diversity of functional haplotypes for a resistance locus in common bean. DNA sequence-based clustering identified eleven haplotype groups at Phg-2. One haplotype group conferred broad-spectrum ALS resistance, six showed pathotype-specific effects, and the remaining seven did not exhibit clear resistance patterns. Our research highlights the importance of ALS pathotype-specificity for durable resistance management strategies in common bean. Molecular markers co-segregating with resistance loci and haplotypes will increase breeding efficiency for ALS resistance and allow to react faster to future changes in pathogen pressure and composition.

Single and Multi-trait GWAS Identify Genetic Factors Associated with Production Traits in Common Bean Under Abiotic Stress Environments.

The genetic improvement of economically important production traits of dry bean (Phaseolus vulgaris L.), for geographic regions where production is threatened by drought and high temperature stress, is challenging because of the complex genetic nature of these traits. Large scale SNP data sets for the two major gene pools of bean, Andean and Middle American, were developed by mapping multiple pools of genotype-by-sequencing reads and identifying over 200k SNPs for each gene pool against the most recent assembly of the P. vulgaris genome sequence. Moderately sized B ean A biotic S tress E valuation (BASE) panels, consisting of genotypes appropriate for production in Central America and Africa, were assembled. Phylogenetic analyses demonstrated the BASE populations represented broad genetic diversity for the appropriate races within the two gene pools. Joint mixed linear model genome-wide association studies with data from multiple locations discovered genetic factors associated with four production traits in both heat and drought stress environments using the BASE panels. Pleiotropic genetic factors were discovered using a multi-trait mixed model analysis. SNPs within or near candidate genes associated with hormone signaling, epigenetic regulation, and ROS detoxification under stress conditions were identified and can be used as genetic markers in dry bean breeding programs.

Fine-mapping of angular leaf spot resistance gene Phg-2 in common bean and development of molecular breeding tools.

KEY MESSAGE: The Common Bean Angular Leaf Spot Resistance Gene Phg-2 was fine-mapped to a 409-Kbp region, and molecular markers for breeders were developed and validated in field experiments. Common bean (Phaseolus vulgaris L.) is an important food legume in Latin America, Asia and Africa. It is an important source of protein, carbohydrates and micro-minerals, particularly for smallholder farmers. Common bean productivity is affected by angular leaf spot (ALS) disease caused by the pathogenic fungus Pseudocercospora griseola, resulting in significant yield losses, particularly in low-input smallholder farming systems in the tropics. The ALS resistance gene Phg-2, which was found in several highly resistant common bean genotypes, was investigated in crosses between Mesoamerican pre-breeding lines and elite Andean breeding lines. Next-generation sequencing (NGS) data sets were used to design new SNP-based molecular markers. The Phg-2 locus was confirmed to be the major locus providing ALS resistance in these crosses. The locus was fine-mapped to a 409-Kbp region on chromosome 8. Two clusters of highly related LRR genes were identified in this region, which are the best candidate genes for Phg-2. Molecular markers were identified that are closely linked to the Phg-2 resistance gene and also highly specific to the donor germplasm. Marker-assisted selection (MAS) was used to introgress the Phg-2 resistance locus into Andean breeding germplasm using MAB lines. The usefulness of molecular markers in MAS was confirmed in several field evaluations in complex breeding crosses, under inoculation with different ALS pathotypes. This project demonstrates that NGS data are a powerful tool for the characterization of genetic loci and can be applied in the development of breeding tools.