The alleles bc-ud and bc-ur (previously bc-4 gene), representing coding mutations within Vps4 AAA+ ATPase ESCRT protein, interact with other genes to condition resistance to BCMV and BCMNV in common bean.

Bean common mosaic virus (BCMV) and bean common mosaic necrosis virus (BCMNV) have a damaging impact on global common bean (Phaseolus vulgaris L.) cultivation, causing potential yield losses of over 80%. The primary strategy for controlling these viruses is through host plant resistance. This research aimed to identify and validate structural variations for the bc-ud gene as revealed by long-read sequencing, develop an efficient DNA marker to assist selection of bc-ud in snap and dry beans, and examine the interactions between the bc-ud allele and other BCMV resistance genes. A gene (Phvul.005G125100) model on chromosome Pv05, encoding a vacuolar protein-sorting 4 (Vps4) AAA+ ATPase endosomal sorting complexes required for transport (ESCRT) protein, was identified as the best candidate gene for bc-ud. An 84-bp repetitive insertion variant within the gene, exhibited 100% co-segregation with the bc-ud resistance allele across 264 common bean accessions. The 84-bp repetitive insertion was labeled with an indel marker IND_05_36225873 which was useful for tracking the bc-ud allele across diverse germplasm. A different single nucleotide polymorphism variant within the same candidate gene was associated with the bc-4 gene. Segregation in F2 populations confirmed bc-ud and bc-4 were alleles, so bc-4 was renamed bc-ur to fit gene nomenclature guidelines. The interactions of bc-ud and bc-ur with other resistance genes, such as bc-1 (receptor-like kinase on Pv03) and bc-2 (Vps4 AAA+ ATPase ESCRT protein on Pv11), validated gene combinations in the differential "host groups" effective against specific BCMV/BCMNV "pathogroups." These findings increase our understanding of the Bc-u locus, and enhance our ability to develop more resilient bean varieties through marker-assisted selection, reducing the impact of BCMV and BCMNV.

Exploring the Mycovirus Sclerotinia sclerotiorum Hypovirulence-Associated DNA Virus 1 as a Biocontrol Agent of White Mold Caused by Sclerotinia sclerotiorum.

Sclerotinia sclerotiorum is a necrotrophic fungal pathogen causing white mold on many important economic crops. Recently, some mycoviruses such as S. sclerotiorum hypovirulence-associated DNA virus 1 (SsHADV-1) converted S. sclerotiorum into a beneficial symbiont that helps plants manage pathogens and other stresses. To explore the potential use of SsHADV-1 as a biocontrol agent in the United States and to test the efficacy of SsHADV-1-infected United States isolates in managing white mold and other crop diseases, SsHADV-1 was transferred from the Chinese strain DT-8 to United States isolates of S. sclerotiorum. SsHADV-1 is readily transmitted horizontally among United States isolates of S. sclerotiorum and consistently conferred hypovirulence to its host strains. Biopriming of dry bean seeds with hypovirulent S. sclerotiorum strains enhanced resistance to white mold, gray mold, and Rhizoctonia root rot. To investigate the underlying mechanisms, endophytic growth of hypovirulent S. sclerotiorum in dry beans was confirmed using PCR, and the expression of 12 plant defense-related genes were monitored before and after infection. The results indicated that the endophytic growth of SsHADV-1-infected strains in plants stimulated the expression of plant immunity pathway genes that assisted a rapid response from the plant to fungal infection. Finally, application of the seed biopriming technology with SsHADV-1-infected hypervirulent strain has promise for the biological control of several diseases of wheat, pea, and sunflower.

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.

Automating high-throughput screening for anthracnose resistance in common bean using allele specific PCR.

BACKGROUND: Common beans (Phaseolus vulgaris L.) provide important protein and calories globally. Anthracnose (Colletotrichum lindemuthianum (Sacc. & Magnus) Briosi & Cavara, 1889) is a major disease in common bean and causes significant yield losses in bean production areas. Screening for markers linked to known disease resistance genes provides useful information for plant breeders to develop improved common bean varieties. The Kompetitive Allele Specific PCR (KASP) assay is an affordable genetic screening technique that can be used to accelerate breeding programs, but manual DNA extraction and KASP assay preparation are time-consuming. Several KASP markers have been developed for genes involved in resistance to bean anthracnose, which can reduce yield by up to 100%, but their usefulness is hindered by the labor required to screen a significant number of bean lines. Our research objective was to develop publicly available protocols for DNA extraction and KASP assaying using a liquid handling robot (LHR) which would facilitate high-throughput genetic screening with less active human time required. Anthracnose resistance markers were used to compare manual and automated results. RESULTS: The 12 bean anthracnose differential cultivars were screened for four anthracnose KASP markers linked to the resistance genes Co-1, Co-3 and Co-42 both by hand and with the use of an LHR. A protocol was written for DNA extraction and KASP assay thermocycling to implement the LHR. The LHR protocol reduced the active human screening time of 24 samples from 3h44 to 1h23. KASP calls were consistent across replicates but not always accurate for their known linked resistance genes, suggesting more specific markers still need to be developed. Using an LHR, information from KASP assays can be accumulated with little active human time. CONCLUSION: Results suggest that LHRs can be used to expedite time-consuming and tedious lab work such as DNA extraction or PCR plate filling. Notably, LHRs can be used to prepare KASP assays for large sample sizes, facilitating higher throughput use of genetic marker screening tools.

Identification of quantitative trait loci for drought tolerance in Bukoba/Kijivu Andean mapping population of common bean.

KEY MESSAGE: Quantitative Trait Loci "hotspots" for drought tolerance were identified on chromosomes Pv06, Pv07 and Pv10 of common bean. Drought is a major production constraint of common bean (Phaseolus vulgaris L.) worldwide. The objective of this study was to identify the Quantitative Trait Loci (QTL) for drought tolerance in an Andean population of Recombinant Inbred Lines (RILs). A total of 155 F5:7 RILs derived from a cross between Kijivu (drought tolerant) and Bukoba (drought susceptible) were evaluated for drought tolerance in field and pot experiments. Four field experiments were conducted at three locations in Zambia in 2020 and 2021. All field trials were conducted in the dry season under irrigation. The 155 RILs were genotyped with 11,292 SNPs, and composite interval mapping was conducted to identify QTL for drought tolerance. Seed yield for Kijivu under drought stress was consistently higher than for Bukoba across all four field trials. A total of 60 QTL were identified for morphological, agronomic, and physiological traits under drought stress and non-stress conditions. However, the majority of these QTL were specific to drought stress. QTL "hotspots" for drought tolerance were identified on chromosomes Pv06, Pv07, and Pv10. Extensive co-localizations for agronomic and morpho-physiological traits under drought stress were observed at the three drought-tolerance QTL hotspots. Additionally, these three QTL hotspots overlapped with previously identified QTL for drought tolerance, while several others identified QTL are novel. The three identified QTL hotspots could be used in future marker-assisted selection for drought tolerance in common bean.

Linked candidate genes of different functions for white mold resistance in common bean (Phaseolus vulgaris L) are identified by multiple QTL mapping approaches.

White mold (WM) is a major disease in common bean (Phaseolus vulgaris L.), and its complex quantitative genetic control limits the development of WM resistant cultivars. WM2.2, one of the nine meta-QTL with a major effect on WM tolerance, explains up to 35% of the phenotypic variation and was previously mapped to a large genomic interval on Pv02. Our objective was to narrow the interval of this QTL using combined approach of classic QTL mapping and QTL-based bulk segregant analysis (BSA), and confirming those results with Khufu de novo QTL-seq. The phenotypic and genotypic data from two RIL populations, 'Raven'/I9365-31 (R31) and 'AN-37'/PS02-029C-20 (Z0726-9), were used to select resistant and susceptible lines to generate subpopulations for bulk DNA sequencing. The QTL physical interval was determined by considering overlapping interval of the identified QTL or peak region in both populations by three independent QTL mapping analyses. Our findings revealed that meta-QTL WM2.2 consists of three regions, WM2.2a (4.27-5.76 Mb; euchromatic), WM 2.2b (12.19 to 17.61 Mb; heterochromatic), and WM2.2c (23.01-25.74 Mb; heterochromatic) found in both populations. Gene models encoding for gibberellin 2-oxidase 8, pentatricopeptide repeat, and heat-shock proteins are the likely candidate genes associated with WM2.2a resistance. A TIR-NBS-LRR class of disease resistance protein (Phvul.002G09200) and LRR domain containing family proteins are potential candidate genes associated with WM2.2b resistance. Nine gene models encoding disease resistance protein [pathogenesis-related thaumatin superfamily protein and disease resistance-responsive (dirigent-like protein) family protein etc] found within the WM2.2c QTL interval are putative candidate genes. WM2.2a region is most likely associated with avoidance mechanisms while WM2.2b and WM2.2c regions trigger physiological resistance based on putative candidate genes.

Integrating de novo QTL-seq and linkage mapping to identify quantitative trait loci conditioning physiological resistance and avoidance to white mold disease in dry bean.

White mold (WM), caused by the ubiquitous fungus Sclerotinia sclerotiorum, is a devastating disease that limits production and quality of dry bean globally. In the present study, classic linkage mapping combined with QTL-seq were employed in two recombinant inbred line (RIL) populations, "Montrose"/I9365-25 (M25) and "Raven"/I9365-31 (R31), with the initial goal of fine-mapping QTL WM5.4 and WM7.5 that condition WM resistance. The RILs were phenotyped for WM reactions under greenhouse (straw test) and field environments. The general region of WM5.4 and WM7.5 were reconfirmed with both mapping strategies within each population. Combining the results from both mapping strategies, WM5.4 was delimited to a 22.60-36.25 Mb interval in the heterochromatic regions on Pv05, while WM7.5 was narrowed to a 0.83 Mb (3.99-4.82 Mb) region on the Pv07 chromosome. Furthermore, additional QTL WM2.2a (3.81-7.24 Mb), WM2.2b (11.18-17.37 Mb, heterochromatic region), and WM2.2c (23.33-25.94 Mb) were mapped to a narrowed genomic interval on Pv02 and WM4.2 in a 0.89 Mb physical interval at the distal end of Pv04 chromosome. Gene models encoding gibberellin 2-oxidase proteins regulating plant architecture are likely candidate genes associated with WM2.2a resistance. Nine gene models encoding a disease resistance protein (quinone reductase family protein and ATWRKY69) found within the WM5.4 QTL interval are putative candidate genes. Clusters of 13 and 5 copies of gene models encoding cysteine-rich receptor-like kinase and receptor-like protein kinase-related family proteins, respectively, are potential candidate genes associated with WM7.5 resistance and most likely trigger physiological resistance to WM. Acquired knowledge of the narrowed major QTL intervals, flanking markers, and candidate genes provides promising opportunities to develop functional molecular markers to implement marker-assisted selection for WM resistant dry bean cultivars.

A robust SNP-haplotype assay for Bct gene region conferring resistance to beet curly top virus in common bean (Phaseolus vulgaris L.).

Beet curly top virus (BCTV), which is synonymous with curly top virus (CTV), causes significant yield loss in common bean (snap and dry beans) cultivars and several other important crops. Common bean cultivars have been found to be resistant to CTV, but screening for resistance is challenging due to the cyclical nature of epidemics and spotty feeding by the leafhopper that vectors the virus. We used an SNP dataset for the Snap Bean Association Panel (SnAP) agro-inoculated with CTV-Logan (CA/Logan) strain to locate the Bct gene region to a 1.7-Mb interval on chromosome Pv07 using genome-wide association study (GWAS) analysis. Recombinant lines from the SnAP were used to further narrow the Bct region to a 58.0-kb interval. A missense SNP (S07_2970381) in candidate gene Phvul.007G036300 Exonuclease V (EXO5) was identified as the most likely causal mutation, and it was the most significant SNP detected by GWAS in a dry bean population (DBP) naturally infected by the CTV-Worland (Wor) strain. Tm-shift assay markers developed for SNP S07_2970381 and two linked SNPs, S07_2970276 and S07_2966197, were useful for tracking different origins of the Bct EXO5 candidate gene resistance to CTV in common bean. The three SNPs identified four haplotypes, with haplotype 3-1 (Haplo3-1) of Middle American origin associated with the highest levels of CTV resistance. This SNP-haplotype assay will enable breeders to track resistance sources and to develop cultivars with better CTV resistance.

The Common Bean V Gene Encodes Flavonoid 3'5' Hydroxylase: A Major Mutational Target for Flavonoid Diversity in Angiosperms.

The classic V (violet, purple) gene of common bean (Phaseolus vulgaris) functions in a complex genetic network that controls seed coat and flower color and flavonoid content. V was cloned to understand its role in the network and the evolution of its orthologs in the Viridiplantae. V mapped genetically to a narrow interval on chromosome Pv06. A candidate gene was selected based on flavonoid analysis and confirmed by recombinational mapping. Protein and domain modeling determined V encodes flavonoid 3'5' hydroxylase (F3'5'H), a P450 enzyme required for the expression of dihydromyricetin-derived flavonoids in the flavonoid pathway. Eight recessive haplotypes, defined by mutations of key functional domains required for P450 activities, evolved independently in the two bean gene pools from a common ancestral gene. V homologs were identified in Viridiplantae orders by functional domain searches. A phylogenetic analysis determined F3'5'H first appeared in the Streptophyta and is present in only 41% of Angiosperm reference genomes. The evolutionarily related flavonoid pathway gene flavonoid 3' hydroxylase (F3'H) is found nearly universally in all Angiosperms. F3'H may be conserved because of its role in abiotic stress, while F3'5'H evolved as a major target gene for the evolution of flower and seed coat color in plants.

Coding Mutations in Vacuolar Protein-Sorting 4 AAA+ ATPase Endosomal Sorting Complexes Required for Transport Protein Homologs Underlie bc-2 and New bc-4 Gene Conferring Resistance to Bean Common Mosaic Virus in Common Bean.

Bean common mosaic virus (BCMV) is a major disease in common bean (Phaseolus vulgaris L.). Host plant resistance is the most effective strategy to minimize crop damage against BCMV and the related Bean common mosaic necrosis virus (BCMNV). To facilitate breeding for resistance, we sought to identify candidate genes and develop markers for the bc-2 gene and the unknown gene with which it interacts. Genome-wide association study (GWAS) of the Durango Diversity Panel (DDP) identified a peak region for bc-2 on chromosome Pv11. Haplotype mapping narrowed the bc-2 genomic interval and identified Phvul.011G092700, a vacuolar protein-sorting 4 (Vps4) AAA+ ATPase endosomal sorting complexes required for transport (ESCRT) protein, as the bc-2 candidate gene. The race Durango Phvul.011G092700 gene model, bc-2 [UI 111], contains a 10-kb deletion, while the race Mesoamerican bc-2 [Robust] consists of a single nucleotide polymorphism (SNP) deletion. Each mutation introduces a premature stop codon, and they exhibit the same interaction with the pathogroups (PGs) tested. Phvul.005G125100, another Vps4 AAA+ ATPase ESCRT protein, was identified as the candidate gene for the new recessive bc-4 gene, and the recessive allele is likely an amino acid substitution in the microtubule interacting and transport (MIT) domain. The two Vps4 AAA+ ATPase ESCRT proteins exhibit high similarity to the Zym Cucsa.385040 candidate gene associated with recessive resistance to Zucchini yellow mosaic virus in cucumber. bc-2 alone has no resistance effect but, when combined with bc-4, provides resistance to BCMV (except PG-V) but not BCMNV, and, when combined with bc-u d, provides resistance to BCMV (except BCMV PG-VII) and BCMNV. So instead of different resistance alleles (i.e., bc-2 and bc-2 2), there is only bc-2 with a differential reaction based on whether it is combined with bc-4 or bc-u d , which are tightly linked in repulsion. The new tools and enhanced understanding of this host-virus pathogen interaction will facilitate breeding common beans for resistance to BCMV and BCMNV.

Genome-Wide Association Mapping of bc-1 and bc-u Reveals Candidate Genes and New Adjustments to the Host-Pathogen Interaction for Resistance to Bean Common Mosaic Necrosis Virus in Common Bean.

Bean common mosaic necrosis virus (BCMNV) is a major disease in common bean (Phaseolus vulgaris L.). Host plant resistance is the primary disease control. We sought to identify candidate genes to better understand the host-pathogen interaction and develop tools for marker-assisted selection (MAS). A genome-wide association study (GWAS) approach using 182 lines from a race Durango Diversity Panel (DDP) challenged by BCMNV isolates NL-8 [Pathogroup (PG)-III] and NL-3 (PG-VI), and genotyped with 1.26 million single-nucleotide polymorphisms (SNPs), revealed significant peak regions on chromosomes Pv03 and Pv05, which correspond to bc-1 and bc-u resistance gene loci, respectively. Three candidate genes were identified for NL-3 and NL-8 resistance. Side-by-side receptor-like protein kinases (RLKs), Phvul.003G038700 and Phvul.003G038800 were candidate genes for bc-1. These RLKs were orthologous to linked RLKs associated with virus resistance in soybean (Glycine max). A basic Leucine Zipper (bZIP) transcription factor protein is the candidate gene for bc-u. bZIP protein gene Phvul.005G124100 carries a unique non-synonymous mutation at codon 14 in the first exon (Pv05: 36,114,516 bases), resulting in a premature termination codon that causes a nonfunctional protein. SNP markers for bc-1 and bc-u and new markers for I and bc-3 genes were used to genotype the resistance genes underpinning BCMNV phenotypes in the DDP, host group (HG) differentials, and segregating F3 families. Results revealed major adjustments to the current host-pathogen interaction model: (i) there is only one resistance allele bc-1 for the Bc-1 locus, and differential expression of the allele is based on presence vs. absence of bc-u; (ii) bc-1 exhibits dominance and incomplete dominance; (iii) bc-1 alone confers resistance to NL-8; (iv) bc-u was absent from HGs 2, 4, 5, and 7 necessitating a new gene symbol bc-u d to reflect this change; (v) bc-u d alone delays susceptible symptoms, and when combined with bc-1 enhanced resistance to NL-3; and (vi) bc-u d is on Pv05, not Pv03 as previously thought. These candidate genes, markers, and adjustments to the host-pathogen interaction will facilitate breeding for resistance to BCMNV and related Bean common mosaic virus (BCMV) in common bean.

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.

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum.

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

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.

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.

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

Pseudomonas syringae pv. phaseolicola (Psph) Race 6 is a globally prevalent and broadly virulent bacterial pathogen with devastating impact causing halo blight of common bean (Phaseolus vulgaris L.). Common bean lines PI 150414 and CAL 143 are known sources of resistance against this pathogen. We constructed high-resolution linkage maps for three recombinant inbred populations to map resistance to Psph Race 6 derived from the two common bean lines. This was complemented with a genome-wide association study (GWAS) of Race 6 resistance in an Andean Diversity Panel of common bean. Race 6 resistance from PI 150414 maps to a single major-effect quantitative trait locus (QTL; HB4.2) on chromosome Pv04 and confers broad-spectrum resistance to eight other races of the pathogen. Resistance segregating in a Rojo × CAL 143 population maps to five chromosome arms and includes HB4.2. GWAS detected one QTL (HB5.1) on chromosome Pv05 for resistance to Race 6 with significant influence on seed yield. The same HB5.1 QTL, found in both Canadian Wonder × PI 150414 and Rojo × CAL 143 populations, was effective against Race 6 but lacks broad resistance. This study provides evidence for marker-assisted breeding for more durable halo blight control in common bean by combining alleles of race-nonspecific resistance (HB4.2 from PI 150414) and race-specific resistance (HB5.1 from cv. Rojo).

Meta-QTL for resistance to white mold in common bean.

White mold, caused by the fungus Sclerotinia sclerotiorum (Lib.) de Bary, is a major disease that limits common bean production and quality worldwide. The host-pathogen interaction is complex, with partial resistance in the host inherited as a quantitative trait with low to moderate heritability. Our objective was to identify meta-QTL conditioning partial resistance to white mold from individual QTL identified across multiple populations and environments. The physical positions for 37 individual QTL were identified across 14 recombinant inbred bi-parental populations (six new, three re-genotyped, and five from the literature). A meta-QTL analysis of the 37 QTL was conducted using the genetic linkage map of Stampede x Red Hawk population as the reference. The 37 QTL condensed into 17 named loci (12 previously named and five new) of which nine were defined as meta-QTL WM1.1, WM2.2, WM3.1, WM5.4, WM6.2, WM7.1, WM7.4, WM7.5, and WM8.3. The nine meta-QTL had confidence intervals ranging from 0.65 to 9.41 Mb. Candidate genes shown to express under S. sclerotiorum infection in other studies, including cell wall receptor kinase, COI1, ethylene responsive transcription factor, peroxidase, and MYB transcription factor, were found within the confidence interval for five of the meta-QTL. The nine meta-QTL are recommended as potential targets for MAS for partial resistance to white mold in common bean.

Sequence-Based Introgression Mapping Identifies Candidate White Mold Tolerance Genes in Common Bean.

White mold, caused by the necrotrophic fungus (Lib.) de Bary, is a major disease of common bean ( L.). WM7.1 and WM8.3 are two quantitative trait loci (QTL) with major effects on tolerance to the pathogen. Advanced backcross populations segregating individually for either of the two QTL, and a recombinant inbred (RI) population segregating for both QTL were used to fine map and confirm the genetic location of the QTL. The QTL intervals were physically mapped using the reference common bean genome sequence, and the physical intervals for each QTL were further confirmed by sequence-based introgression mapping. Using whole-genome sequence data from susceptible and tolerant DNA pools, introgressed regions were identified as those with significantly higher numbers of single-nucleotide polymorphisms (SNPs) relative to the whole genome. By combining the QTL and SNP data, WM7.1 was located to a 660-kb region that contained 41 gene models on the proximal end of chromosome Pv07, while the WM8.3 introgression was narrowed to a 1.36-Mb region containing 70 gene models. The most polymorphic candidate gene in the WM7.1 region encodes a BEACH-domain protein associated with apoptosis. Within the WM8.3 interval, a receptor-like protein with the potential to recognize pathogen effectors was the most polymorphic gene. The use of gene and sequence-based mapping identified two candidate genes whose putative functions are consistent with the current model of pathogenicity.

Application of in silico bulked segregant analysis for rapid development of markers linked to Bean common mosaic virus resistance in common bean.

BACKGROUND: Common bean was one of the first crops that benefited from the development and utilization of molecular marker-assisted selection (MAS) for major disease resistance genes. Efficiency of MAS for breeding common bean is still hampered, however, due to the dominance, linkage phase, and loose linkage of previously developed markers. Here we applied in silico bulked segregant analysis (BSA) to the BeanCAP diversity panel, composed of over 500 lines and genotyped with the BARCBEAN_3 6K SNP BeadChip, to develop codominant and tightly linked markers to the I gene controlling resistance to Bean common mosaic virus (BCMV). RESULTS: We physically mapped the genomic region underlying the I gene. This locus, in the distal arm of chromosome Pv02, contains seven putative NBS-LRR-type disease resistance genes. Two contrasting bulks, containing BCMV host differentials and ten BeanCAP lines with known disease reaction to BCMV, were subjected to in silico BSA for targeting the I gene and flanking sequences. Two distinct haplotypes, containing a cluster of six single nucleotide polymorphisms (SNP), were associated with resistance or susceptibility to BCMV. One-hundred and twenty-two lines, including 115 of the BeanCAP panel, were screened for BCMV resistance in the greenhouse, and all of the resistant or susceptible plants displayed distinct SNP haplotypes as those found in the two bulks. The resistant/susceptible haplotypes were validated in 98 recombinant inbred lines segregating for BCMV resistance. The closest SNP (~25-32 kb) to the distal NBS-LRR gene model for the I gene locus was targeted for conversion to codominant KASP (Kompetitive Allele Specific PCR) and CAPS (Cleaved Amplified Polymorphic Sequence) markers. Both marker systems accurately predicted the disease reaction to BCMV conferred by the I gene in all screened lines of this study. CONCLUSIONS: We demonstrated the utility of the in silico BSA approach using genetically diverse germplasm, genotyped with a high-density SNP chip array, to discover SNP variation at a specific targeted genomic region. In common bean, many disease resistance genes are mapped and their physical genomic position can now be determined, thus the application of this approach will facilitate further development of codominant and tightly linked markers for use in MAS.

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

Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement.