Genome-Wide Linkage and Association Mapping of Halo Blight Resistance in Common Bean to Race 6 of the Globally Important Bacterial Pathogen

GWAS of resistance to three bacterial diseases in the Andean common bean diversity panel

Bacterial brown spot (BBS) caused by Pseudomonas syringae pv. syringae (Pss), common bacterial blight (CBB) caused by Xanthomonas axonopodis pv. phaseoli (Xap) and Xanthomonas fuscans subsp. fuscans (Xff), and halo bacterial blight (HBB), caused by Pseudomonas syringae pv. phaseolicola (Psph), are major bacterial diseases that severely affect common bean yields and global food security. Andean-origin dry beans, representing large-seeded market classes, are particularly susceptible. Using 140,325 SNPs, a multi-locus GWAS was conducted on subsets of the Andean diversity panel (ADP) phenotyped for BBS in South Africa, CBB in Puerto Rico, South Africa, and Zambia, and HBB in South Africa, through natural infection, artificial inoculation, or both. Twenty-four QTL associated with resistance were identified: nine for BBS, eight for CBB, and seven for HBB. Four QTL intervals on Pv01, Pv03, Pv05, and Pv08 overlapped with BBS and HBB resistance. A genomic interval on Pv01, near the fin gene, which determines growth habit, was linked to resistance to all three pathogens. Different QTLs were detected for BBS and CBB resistance when phenotyped under natural infection versus artificial inoculation. These results underscore the importance of combining phenotyping methods in multi-GWAS to capture the full genetic spectrum. Previously recognized CBB resistance QTL SAP6 and SU91 and HBB resistance QTL HB4.2, and HB5.1, were observed. Other common (MAF >0.25) and rare (MAF <0.05) resistance QTL were also detected. Overall, these findings enhance the understanding and utilization of bacterial resistance present in ADP for the development of common beans with improved resistance.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Genome-wide association analysis of anthracnose resistance in the Yellow Bean Collection of Common Bean

Anthracnose caused by Colletotrichum lindemuthianum is a major disease of common bean (Phaseolus vulgaris) worldwide. Yellow beans are a major market class of common bean especially in eastern and southern Africa. The Yellow Bean Collection (YBC), which is comprised of 255 genotypes, and has not been used previously in genetic studies on anthracnose, is an excellent genetic resource for understanding the extent of anthracnose resistance and its genetic architecture in the yellow bean market class. The objectives of this study were i) evaluate the YBC for resistance to races 5, 19, 39, 51, 81, 183, 1050 and 1105 of C. lindemuthianum. and ii) conduct genome-wide association analysis to identify genomic regions and candidate genes associated with resistance to C. lindemuthianum. The YBC was genotyped with 72,866 SNPs, and genome-wide association analysis was conducted using Mixed Linear Model in TASSEL. Andean and Middle American genotypes with superior levels of resistance to the eight races were identified. YBC278 was the only one among 255 genotypes that was highly resistant to all eight races. Resistance to anthracnose in the YBC was controlled by major-effect loci on chromosomes Pv01, Pv03, Pv04, Pv05 and Pv07. The genomic region on Pv01, which overlapped with the Andean locus Co-1 provided resistance to races 81, 1050 and 1105. Significant SNPs for resistance to race 39 were identified on Pv02. The genomic region on Pv04, which overlapped with known major-effect loci Co-3, Co-15, Co-16, Co-y and Co-z, provided resistance to races 5, 19, 51 and 183. Novel genomic regions for resistance to race 39 were identified on Pv05 and Pv07. Plant resistance genes (R genes) with NB-ARC and LRR domains, which occurred in clusters, were identified as positional candidate genes for genomic regions on Pv02 and Pv04.

Genome-wide association mapping for component traits of drought tolerance in dry beans (Phaseolus vulgaris L.)

Understanding the genetic basis of traits of economic importance under drought stressed and well-watered conditions is important in enhancing genetic gains in dry beans (Phaseolus vulgaris L.). This research aims to: (i) identify markers associated with agronomic and physiological traits for drought tolerance and (ii) identify drought-related putative candidate genes within the mapped genomic regions. An andean and middle-american diversity panel (AMDP) comprising of 185 genotypes was screened in the field under drought stressed and well-watered conditions for two successive seasons. Agronomic and physiological traits, viz., days to 50% flowering (DFW), plant height (PH), days to physiological maturity (DPM), grain yield (GYD), 100-seed weight (SW), leaf temperature (LT), leaf chlorophyll content (LCC) and stomatal conductance (SC) were phenotyped. Principal component and association analysis were conducted using the filtered 9370 Diversity Arrays Technology sequencing (DArTseq) markers. The mean PH, GYD, SW, DPM, LCC and SC of the panel was reduced by 12.1, 29.6, 10.3, 12.6, 28.5 and 62.0%, respectively under drought stressed conditions. Population structure analysis revealed two sub-populations, which corresponded to the andean and middle-american gene pools. Markers explained 0.08-0.10, 0.22-0.23, 0.29-0.32, 0.43-0.44, 0.65-0.66 and 0.69-0.70 of the total phenotypic variability (R2) for SC, LT, PH, GYD, SW and DFW, respectively under drought stressed conditions. For well-watered conditions, R2 varied from 0.08 (LT) to 0.70 (DPM). Overall, 68 significant (p < 10-03) marker-trait associations (MTAs) and 22 putative candidate genes were identified across drought stressed and well-watered conditions. Most of the identified genes had known biological functions related to regulating the response to drought stress. The findings provide new insights into the genetic architecture of drought stress tolerance in common bean. The findings also provide potential candidate SNPs and putative genes that can be utilized in gene discovery and marker-assisted breeding for drought tolerance after validation.

Genetic mapping of the Andean anthracnose resistance gene present in the common bean cultivar BRSMG Realce

The rajado seeded Andean bean (Phaseolus vulgaris L.) cultivar BRSMG Realce (striped seed coat) developed by Embrapa expressed a high level of anthracnose resistance, caused by Colletotrichum lindemuthianum, in field and greenhouse screenings. The main goal of this study was to evaluate the inheritance of anthracnose resistance in BRSMG Realce, map the resistance locus or major gene cluster previously named as Co-Realce, identify resistance-related positional genes, and analyze potential markers linked to the resistance allele. F₂ plants derived from the cross BRSMG Realce × BRS FC104 (Mesoamerican) and from the cross BRSMG Realce × BRS Notável (Mesoamerican) were inoculated with the C. lindemuthianum races 475 and 81, respectively. The BRSMG Realce × BRS FC104 F₂ population was also genotyped using the DArTseq technology. Crosses between BRSMG Realce and BAT 93 (Mesoamerican) were also conducted and resulting F₂ plants were inoculated with the C. lindemuthianum races 65 and 1609, individually. The results shown that anthracnose resistance in BRSMG Realce is controlled by a single locus with complete dominance. A genetic map including 1,118 SNP markers was built and shown 78% of the markers mapped at a distances less than 5.0 cM, with a total genetic length of 4,473.4 cM. A major locus (Co-Realce) explaining 54.6% of the phenotypic variation of symptoms caused by the race 475 was identified in Pv04, flanked by the markers snp1327 and snp12782 and 4.48 cM apart each other. These SNPs are useful for marker-assisted selection, due to an estimated selection efficiency of 99.2%. The identified resistance allele segregates independently of the resistance allele Co-3³ (Pv04) present in BAT 93. The mapped genomic region with 704,867 bp comprising 63 putative genes, 44 of which were related to the pathogen-host interaction. Based on all these results and evidence, anthracnose resistance in BRSMG Realce should be considered as monogenic, useful for breeding purpose. It is proposed that locus Co-Realce is unique and be provisionally designated as CoPv04^R until be officially nominated in accordance with the rules established by the Bean Improvement Cooperative Genetics Committee.

Characterization of fungal pathogens and germplasm screening for disease resistance in the main production area of the common bean in Argentina

The common bean (Phaseolus vulgaris L.) is the most important grain legume in the human diet, mainly in Africa and Latin America. Argentina is one of the five major producers of the common bean in the world, and the main cultivation areas are concentrated in the northwestern provinces of this country. Crop production of the common bean is often affected by biotic factors like some endemic fungal diseases, which exert a major economic impact on the region. The most important fungal diseases affecting the common bean in Argentina are white mold caused by Sclerotinia sclerotiorum, angular leaf spot caused by Pseudocercospora griseola, web blight and root rot caused by Rhizoctonia solani, which can cause production losses of up to 100% in the region. At the present, the most effective strategy for controlling these diseases is the use of genetic resistance. In this sense, population study and characterization of fungal pathogens are essential for developing cultivars with durable resistance. In this review we report diversity studies carried out on these three fungal pathogens affecting the common bean in northwestern Argentina, analyzing more than 200 isolates by means of molecular, morphological and pathogenic approaches. Also, the screening of physiological resistance in several common bean commercial lines and wild native germplasm is reviewed. This review contributes to the development of sustainable management strategies and cultural practices in bean production aimed to minimize yield losses due to fungal diseases in the common bean.

Pseudomonas phaseolicola preferentially modulates genes encoding leucine-rich repeat and malectin domains in the bean landrace G2333

Candidate resistance genes encoding malectin-like and LRR domains mapped to halo blight resistance loci throughout the common bean genome are co-expressed to fight a range of Pph races. Common bean (Phaseolus vulgaris L.) is an important crop both as a source of protein and other nutrients for human nutrition and as a nitrogen fixer that benefits sustainable agriculture. This crop is affected by halo blight disease, caused by the bacterium Pseudomonas syringae pv. phaseolicola (Pph), which can lead to 45% yield losses. Common bean resistance to Pph is conferred by six loci (Pse-1 to Pse-6) and minor-effect quantitative trait loci (QTLs); however, information is lacking on the molecular mechanisms implicated in this resistance. Here, we describe an in-depth RNA-sequencing (RNA-seq) analysis of the tolerant G2333 bean line in response to the Pph strain NPS3121. We identified 275 upregulated and 357 downregulated common bean genes in response to Pph infection. These differentially expressed genes were mapped to all 11 chromosomes of P. vulgaris. The upregulated genes were primarily components of plant immune responses and negative regulation of photosynthesis, with enrichment for leucine-rich repeat (LRRs) and/or malectin-like carbohydrate-binding domains. Interestingly, LRRs and malectin genes mapped to the same location as previously identified Pph resistance loci or QTLs. For instance, the major loci Pse-6/HB4.2 involved in broad-resistance to many Pph races co-located with induced LRR-encoding genes on Pv04. These findings indicate a coordinated modulation of genes involved in pathogen perception and signal transduction. In addition, the results further support these LRR/malectin loci as resistance genes in response to halo blight. Thus, these genes are potential targets for future genetic manipulation, enabling the introduction of resistance to Pph into elite cultivars of common bean.

Genome-Wide Association Study Identifies Genomic Regions for Important Morpho-Agronomic Traits in Mesoamerican Common Bean

The population growth trend in recent decades has resulted in continuing efforts to guarantee food security in which leguminous plants, such as the common bean (Phaseolus vulgaris L.), play a particularly important role as they are relatively cheap and have high nutritional value. To meet this demand for food, the main target for genetic improvement programs is to increase productivity, which is a complex quantitative trait influenced by many component traits. This research aims to identify Quantitative Trait Nucleotides (QTNs) associated with productivity and its components using multi-locus genome-wide association studies. Ten morpho-agronomic traits [plant height (PH), first pod insertion height (FPIH), number of nodules (NN), pod length (PL), total number of pods per plant (NPP), number of locules per pod (LP), number of seeds per pod (SP), total seed weight per plant (TSW), 100-seed weight (W100), and grain yield (YLD)] were evaluated in four environments for 178 Mesoamerican common bean domesticated accessions belonging to the Brazilian Diversity Panel. In order to identify stable QTNs, only those identified by multiple methods (mrMLM, FASTmrMLM, pLARmEB, and ISIS EM-BLASSO) or in multiple environments were selected. Among the identified QTNs, 64 were detected at least thrice by different methods or in different environments, and 39 showed significant phenotypic differences between their corresponding alleles. The alleles that positively increased the corresponding traits, except PH (for which lower values are desired), were considered favorable alleles. The most influenced trait by the accumulation of favorable alleles was PH, showing a 51.7% reduction, while NN, TSW, YLD, FPIH, and NPP increased between 18 and 34%. Identifying QTNs in several environments (four environments and overall adjusted mean) and by multiple methods reinforces the reliability of the associations obtained and the importance of conducting these studies in multiple environments. Using these QTNs through molecular techniques for genetic improvement, such as marker-assisted selection or genomic selection, can be a strategy to increase common bean production.

Towards the Development, Maintenance, and Standardized Phenotypic Characterization of Single-Seed-Descent Genetic Resources for Common Bean

The optimal use of legume genetic resources represents a key prerequisite for coping with current agriculture-related societal challenges, including conservation of agrobiodiversity, agricultural sustainability, food security, and human health. Among legumes, the common bean (Phaseolus vulgaris) is the most economically important for human consumption, and its evolutionary trajectories as a species have been crucial to determining the structure and level of its present and available genetic diversity. Genomic advances are considerably enhancing the characterization and assessment of important genetic variants. For this purpose, the development and availability of, and access to, well-described and efficiently managed genetic resource collections that comprise pure lines derived by single-seed-descent cycles will be paramount for the use of the reservoir of common bean variability and for the advanced breeding of legume crops. This is one of the main aims of the new and challenging European project INCREASE, which is the implementation of Intelligent Collections with appropriate standardized protocols that must be characterized, maintained, and made available, along with the related data, to users such as breeders and researchers. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Characterizing common bean seeds for seed trait descriptors Basic Protocol 2: Bean seed imaging Basic Protocol 3: Characterizing bean lines for plant trait descriptors specific for common bean Primary Seed Increase.

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.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Genome-Wide Association Studies of Mineral Content in Common Bean

Micronutrient malnutrition is one of the main public health problems in many parts of the world. This problem raises the attention of all valuable sources of micronutrients for the human diet, such as common bean (Phaseolus vulgaris L.). In this research, a panel of 174 accessions representing Croatian common bean landraces was phenotyped for seed content of eight nutrients (N, P, K, Ca, Mg, Fe, Zn, and Mn), and genotyped using 6,311 high-quality DArTseq-derived SNP markers. A genome-wide association study (GWAS) was then performed to identify new genetic sources for improving seed mineral content. Twenty-two quantitative trait nucleotides (QTN) associated with seed nitrogen content were discovered on chromosomes Pv01, Pv02, Pv03, Pv05, Pv07, Pv08, and Pv10. Five QTNs were associated with seed phosphorus content, four on chromosome Pv07, and one on Pv08. A single significant QTN was found for seed calcium content on chromosome Pv09 and for seed magnesium content on Pv08. Finally, two QTNs associated with seed zinc content were identified on Pv06 while no QTNs were found to be associated with seed potassium, iron, or manganese content. Our results demonstrate the utility of GWAS for understanding the genetic architecture of seed nutritional traits in common bean and have utility for future enrichment of seed with macro- and micronutrients through genomics-assisted breeding.

Genome-Wide Association Studies Detect Multiple QTLs for Productivity in Mesoamerican Diversity Panel of Common Bean Under Drought Stress

Drought stress is an important abiotic factor limiting common bean yield, with great impact on the production worldwide. Understanding the genetic basis regulating beans' yield and seed weight (SW) is a fundamental prerequisite for the development of superior cultivars. The main objectives of this work were to conduct genome-wide marker discovery by genotyping a Mesoamerican panel of common bean germplasm, containing cultivated and landrace accessions of broad origin, followed by the identification of genomic regions associated with productivity under two water regimes using different genome-wide association study (GWAS) approaches. A total of 11,870 markers were genotyped for the 339 genotypes, of which 3,213 were SilicoDArT and 8,657 SNPs derived from DArT and CaptureSeq. The estimated linkage disequilibrium extension, corrected for structure and relatedness (r ² _sv ), was 98.63 and 124.18 kb for landraces and breeding lines, respectively. Germplasm was structured into landraces and lines/cultivars. We carried out GWASs for 100-SW and yield in field environments with and without water stress for 3 consecutive years, using single-, segment-, and gene-based models. Higher number of associations at high stringency was identified for the SW trait under irrigation, totaling ∼185 QTLs for both single- and segment-based, whereas gene-based GWASs showed ∼220 genomic regions containing ∼650 genes. For SW under drought, 18 QTLs were identified for single- and segment-based and 35 genes by gene-based GWASs. For yield, under irrigation, 25 associations were identified, whereas under drought the total was 10 using both approaches. In addition to the consistent associations detected across experiments, these GWAS approaches provided important complementary QTL information (∼221 QTLs; 650 genes; r ² from 0.01% to 32%). Several QTLs were mined within or near candidate genes playing significant role in productivity, providing better understanding of the genetic mechanisms underlying these traits and making available molecular tools to be used in marker-assisted breeding. The findings also allowed the identification of genetic material (germplasm) with better yield performance under drought, promising to a common bean breeding program. Finally, the availability of this highly diverse Mesoamerican panel is of great scientific value for the analysis of any relevant traits in common bean.

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.

The Proteomics of Resistance to Halo Blight in Common Bean

Halo blight disease of beans is caused by a gram-negative bacterium, Pseudomonas syringae pv. phaseolicola. The disease is prevalent in South America and Africa and causes crop loss for indigent people who rely on beans as a primary source of daily nutrition. In susceptible beans, P. syringae pv. phaseolicola causes water-soaking at the site of infection and produces phaseolotoxin, an inhibitor of bean arginine biosynthesis. In resistant beans, P. syringae pv. phaseolicola triggers a hypersensitive response that limits the spread of infection. Here, we used high-throughput mass spectrometry to interrogate the responses to two different P. syringae pv. phaseolicola isolates on a single line of common bean, Phaseolus vulgaris PI G19833, with a reference genome sequence. We obtained quantitative information for 4,135 bean proteins. A subset of 160 proteins with similar accumulation changes during both susceptible and resistant reactions included salicylic acid responders EDS1 and NDR1, ethylene and jasmonic acid biosynthesis enzymes, and proteins enabling vesicle secretion. These proteins revealed the activation of a basal defense involving hormonal responses and the mobilization of extracellular proteins. A subset of 29 proteins specific to hypersensitive immunity included SOBIR1, a G-type lectin receptor-like kinase, and enzymes needed for glucoside and phytoalexin production. Virus-induced gene silencing revealed that the G-type lectin receptor-like kinase suppresses bacterial infection. Together, the results define the proteomics of disease resistance to P. syringae pv. phaseolicola in beans and support a model whereby the induction of hypersensitive immunity reinstates defenses targeted by P. syringae pv. phaseolicola.

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.

Natural Variation in Portuguese Common Bean Germplasm Reveals New Sources of Resistance Against Fusarium oxysporum f. sp. phaseoli and Resistance-Associated Candidate Genes

Common bean (Phaseolus vulgaris) is one of the most consumed legume crops in the world, and Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. phaseoli, is one of the major diseases affecting its production. Portugal holds a very promising common bean germplasm with an admixed genetic background that may reveal novel genetic resistance combinations between the original Andean and Mesoamerican gene pools. To identify new sources of Fusarium wilt resistance and detect resistance-associated single-nucleotide polymorphisms (SNPs), we explored, for the first time, a diverse collection of the underused Portuguese common bean germplasm by using genome-wide association analyses. The collection was evaluated for Fusarium wilt resistance under growth chamber conditions, with the highly virulent F. oxysporum f. sp. phaseoli strain FOP-SP1 race 6. Fourteen of the 162 Portuguese accessions evaluated were highly resistant and 71 intermediate. The same collection was genotyped with DNA sequencing arrays, and SNP-resistance associations were tested via a mixed linear model accounting for the genetic relatedness between accessions. The results from the association mapping revealed nine SNPs associated with resistance on chromosomes Pv04, Pv05, Pv07, and Pv08, indicating that Fusarium wilt resistance is under oligogenic control. Putative candidate genes related to phytoalexin biosynthesis, hypersensitive response, and plant primary metabolism were identified. The results reported here highlight the importance of exploring underused germplasm for new sources of resistance and provide new genomic targets for the development of functional markers to support selection in future disease resistance breeding programs.

Improving the Health Benefits of Snap Bean: Genome-Wide Association Studies of Total Phenolic Content

Snap beans are a significant source of micronutrients in the human diet. Among the micronutrients present in snap beans are phenolic compounds with known beneficial effects on human health, potentially via their metabolism by the gut-associated microbiome. The genetic pathways leading to the production of phenolics in snap bean pods remain uncertain. In this study, we quantified the level of total phenolic content (TPC) in the Bean Coordinated Agriculture Program (CAP) snap bean diversity panel of 149 accessions. The panel was characterized spectrophotometrically for phenolic content with a Folin-Ciocalteu colorimetric assay. Flower, seed and pod color were also quantified, as red, purple, yellow and brown colors are associated with anthocyanins and flavonols in common bean. Genotyping was performed through an Illumina Infinium Genechip BARCBEAN6K_3 single nucleotide polymorphism (SNP) array. Genome-Wide Association Studies (GWAS) analysis identified 11 quantitative trait nucleotides (QTN) associated with TPC. An SNP was identified for TPC on Pv07 located near the P gene, which is a major switch in the flavonoid biosynthetic pathway. Candidate genes were identified for seven of the 11 TPC QTN. Five regulatory genes were identified and represent novel sources of variation for exploitation in developing snap beans with higher phenolic levels for greater health benefits to the consumer.

Winged bean (Psophocarpus tetragonolobus (L.) DC.) for food and nutritional security: synthesis of past research and future direction

Winged bean is popularly known as "One Species Supermarket" for its nutrient-dense green pods, immature seeds, tubers, leaves, and mature seeds. This underutilised crop has potential beneficial traits related to its biological nitrogen-fixation to support low-input farming. Drawing from past knowledge, and based on current technologies, we propose a roadmap for research and development of winged bean for sustainable food systems. Reliance on a handful of "major" crops has led to decreased diversity in crop species, agricultural systems and human diets. To reverse this trend, we need to encourage the greater use of minor, "orphan", underutilised species. These could contribute to an increase in crop diversity within agricultural systems, to improve human diets, and to support more sustainable and resilient food production systems. Among these underutilised species, winged bean (Psophocarpus tetragonolobus) has long been proposed as a crop for expanded use particularly in the humid tropics. It is an herbaceous perennial legume of equatorial environments and has been identified as a rich source of protein, with most parts of the plant being edible when appropriately prepared. However, to date, limited progress in structured improvement programmes has restricted the expansion of winged bean beyond its traditional confines. In this paper, we discuss the reasons for this and recommend approaches for better use of its genetic resources and related Psophocarpus species in developing improved varieties. We review studies on the growth, phenology, nodulation and nitrogen-fixation activity, breeding programmes, and molecular analyses. We then discuss prospects for the crop based on the greater understanding that these studies have provided and considering modern plant-breeding technologies and approaches. We propose a more targeted and structured research approach to fulfil the potential of winged bean to contribute to food security.

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.

A Review of Angular Leaf Spot Resistance in Common Bean

Angular leaf spot (ALS), caused by Pseudocercospora griseola, is one of the most devastating diseases of common bean (Phaseolus vulgaris L.) in tropical and subtropical production areas. Breeding for ALS resistance is difficult due to the extensive virulence diversity of P. griseola and the recurrent appearance of new virulent races. Five major loci, Phg-1 to Phg-5, conferring ALS resistance have been named, and markers tightly linked to these loci have been reported. Quantitative trait loci (QTLs) have also been described, but the validation of some QTLs is still pending. The Phg-1, Phg-4, and Phg-5 loci are from common bean cultivars of the Andean gene pool, whereas Phg-2 and Phg-3 are from beans of the Mesoamerican gene pool. The reference genome of common bean and high-throughput sequencing technologies are enabling the development of molecular markers closely linked to the Phg loci, more accurate mapping of the resistance loci, and the comparison of their genomic positions. The objective of this report is to provide a comprehensive review of ALS resistance in common bean. Furthermore, we are reporting three case studies of ALS resistance breeding in Latin America and Africa. This review will serve as a reference for future resistance mapping studies and as a guide for the selection of resistance loci in breeding programs aiming to develop common bean cultivars with durable ALS resistance.

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.

Evidence for Strong Kinship Influence on the Extent of Linkage Disequilibrium in Cultivated Common Beans

Phaseolus vulgaris is an important grain legume for human consumption. Recently, association mapping studies have been performed for the species aiming to identify loci underlying quantitative variation of traits. It is now imperative to know whether the linkage disequilibrium (LD) reflects the true association between a marker and causative loci. The aim of this study was to estimate and analyze LD on a diversity panel of common beans using ordinary r² and r2 extensions which correct bias due to population structure (rS²), kinship (rV²), and both (rVS²). A total of 10,362 single nucleotide polymorphisms (SNPs) were identified by genotyping by sequencing (GBS), and polymorphisms were found to be widely distributed along the 11 chromosomes. In terms of r2, high values of LD (over 0.8) were identified between SNPs located at opposite chromosomal ends. Estimates for rV² were lower than those for rS². Results for rV² and rVS² were similar, suggesting that kinship may also include information on population structure. Over genetic distance, LD decayed to 0.1 at a distance of 1 Mb for rVS². Inter-chromosomal LD was also evidenced. This study showed that LD estimates decay dramatically according to the population structure, and especially the degree of kinship. Importantly, the LD estimates reported herein may influence our ability to perform association mapping studies on P. vulgaris.

Dissection of Resistance Genes to Pseudomonas syringae pv. phaseolicola in UI3 Common Bean Cultivar

Few quantitative trait loci have been mapped for resistance to Pseudomonas syringae pv. phaseolicola in common bean. Two F₂ populations were developed from the host differential UI3 cultivar. The objective of this study was to further characterize the resistance to races 1, 5, 7 and 9 of Psp included in UI3. Using a QTL mapping approach, 16 and 11 main-effect QTLs for pod and primary leaf resistance were located on LG10, explaining up to 90% and 26% of the phenotypic variation, respectively. The homologous genomic region corresponding to primary leaf resistance QTLs detected tested positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL), Natural Resistance Associated Macrophage (NRAMP) and Pentatricopeptide Repeat family (PPR) proteins. It is worth noting that the main effect QTLs for resistance in pod were located inside a 3.5 Mb genomic region that included the Phvul.010G021200 gene, which encodes a protein that has the highest sequence similarity to the RIN4 gene of Arabidopsis, and can be considered an important candidate gene for the organ-specific QTLs identified here. These results support that resistance to Psp from UI3 might result from the immune response activated by combinations of R proteins, and suggest the guard model as an important mechanism in pod resistance to halo blight. The candidate genes identified here warrant functional studies that will help in characterizing the actual defense gene(s) in UI3 genotype.