Search

Flavonoids are secondary metabolites associated with plant seed coat and flower color. These compounds provide health benefits to humans as anti-inflammatory and antioxidant compounds. The expression of the late biosynthetic genes in the flavonoid pathway is controlled by a ternary MBW protein complex consisting of interfacing MYB, beta-helix-loop-helix (bHLH), and WD40 Repeat (WDR) proteins. P, the master regulator gene of the flavonoid expression in common bean (Phaseolus vulgaris L.), was recently determined to encode a bHLH protein. The T and Z genes control the distribution of color in bean seeds and flowers and have historically been considered regulators of the flavonoid gene expression. T and Z candidates were identified using reverse genetics based on genetic mapping, phylogenetic analysis, and mutant analysis. Domain and AlphaFold2 structure analyses determined that T encodes a seven-bladed β-propeller WDR protein, while Z encodes a R2R3 MYB protein. Deletions and SNPs in T and Z mutants, respectively, altered the 3D structure of these proteins. Modeling of the Z MYB/P bHLH/T WDR MBW complex identified interfacing sequence domains and motifs in all three genes that are conserved in dicots. One Z MYB motif is a possible beta-molecular recognition feature (β-MoRF) that only appears in a structured state when Z MYB is modeled as a component of a MBW complex. Complexes containing mutant T and Z proteins changed the interaction of members of the complex in ways that would alter their role in regulating the expression of genes in the flavonoid pathway., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Soler-Garzón, Mulube, Kamfwa, Lungu, Hamabwe, Roy, Salegua, Fourie, Porch, McClean and Miklas.)

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

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-u d gene as revealed by long-read sequencing, develop an efficient DNA marker to assist selection of bc-u d in snap and dry beans, and examine the interactions between the bc-u d 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-u d . An 84-bp repetitive insertion variant within the gene, exhibited 100% co-segregation with the bc-u d resistance allele across 264 common bean accessions. The 84-bp repetitive insertion was labeled with an indel marker IND&#95;05&#95;36225873 which was useful for tracking the bc-u d allele across diverse germplasm. A different single nucleotide polymorphism variant within the same candidate gene was associated with the bc-4 gene. Segregation in F 2 populations confirmed bc-u d and bc-4 were alleles, so bc-4 was renamed bc-u r to fit gene nomenclature guidelines. The interactions of bc-u d and bc-u r 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., (© 2023 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

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., (© 2023 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Oladzad, Roy, Mamidi, Miklas, Lee, Clevenger, Myers, Korani and McClean.)

Anthracnose, caused by the fungal pathogen Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.-Scrib., is one of the most devastating diseases in dry bean ( Phaseolus vulgaris L.) with seed yield losses up to 100%. Most anthracnose resistance genes thus far identified behave in a dominant manner and were identified by seedling screening. The Middle American Diversity Panel (MDP; n=266) was screened with a modified greenhouse screening method to evaluate the response to anthracnose race 73. Thirty MDP genotypes exhibited resistance to the race of which 16 genotypes were not known to contain anthracnose resistance genes to race 73. GWAS with ~93,000 SNP markers identified four genomic regions, two each on Pv01 and Pv10, associated race 73 resistance. A likelihood-ratio-based R2 analysis indicated the peak four SNP markers are responsible for 26% of the observed phenotypic variation, where one SNP, S10&#95;072250, explains 23% of the total variation. SNP S10&#95;072250 is associated with a new region of anthracnose resistance and is in an intron of a ZPR1 -like gene. Further greenhouse testing of the 16 resistant lines without previously known resistance to race 73 revealed various levels of resistance under various levels of disease pressure. Disease resistance was further characterized in the field using four representative genotypes. GTS-900 and Remington exhibited field resistance while Merlot and Maverick were susceptible. Field testing with two different fungicide regimes revealed the resistant genotypes had no significant disease differences. The results suggest resistance to anthracnose may differ at various growth stages and that breeders have been selecting for major genes at early seedling stages while ignoring the effect of alternative genes that may be active at later stages. The newly identified resistant lines may be related to Age Related Resistance (ARR) and could be exploited as parental sources of anthracnose resistance in addition to already known major genes. The physical localization of the multiple regions of resistance confirms the presence of two clusters of disease resistance genes on Pv01 and identifies two new regions of anthracnose resistance on Pv10 possibly associated with ARR. Future research should look at the mode of inheritance of this resistance and its effect when combined with other anthracnose resistance loci., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Simons, Schröder, Oladzad, McClean, Conner, Penner, Stoesz and Osorno.)

Key Message: GWAS detected ninety-eight significant SNPs associated with Sclerotinia sclerotiorum resistance. Six statistical models resulted in medium to high predictive ability, depending on trait, indicating potential of genomic prediction for disease resistance breeding. The lack of complete host resistance and a complex resistance inheritance nature between rapeseed/canola and Sclerotinia sclerotiorum often limits the development of functional molecular markers that enable breeding for sclerotinia stem rot (SSR) resistance. However, genomics-assisted selection has the potential to accelerate the breeding for SSR resistance. Therefore, genome-wide association (GWA) mapping and genomic prediction (GP) were performed using a diverse panel of 337 rapeseed/canola genotypes. Three-week-old seedlings were screened using the petiole inoculation technique (PIT). Days to wilt (DW) up to 2 weeks and lesion phenotypes (LP) at 3, 4, and 7 days post-inoculation (dpi) were recorded. A strong correlation (r = - 0.90) between DW and LP&#95;4dpi implied that a single time point scoring at four days could be used as a proxy trait. GWA analyses using single-locus (SL) and multi-locus (ML) models identified a total of 41, and 208 significantly associated SNPs, respectively. Out of these, ninety-eight SNPs were identified by a combination of the SL model and any of the ML models, at least two ML models, or two traits. These SNPs explained 1.25-12.22% of the phenotypic variance and considered as significant, could be associated with SSR resistance. Eighty-three candidate genes with a function in disease resistance were associated with the significant SNPs. Six GP models resulted in moderate to high (0.42-0.67) predictive ability depending on SSR resistance traits. The resistant genotypes and significant SNPs will serve as valuable resources for future SSR resistance breeding. Our results also highlight the potential of genomic selection to improve rapeseed/canola breeding for SSR resistance., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 McClean, Lee, Howe, Osborne, Grimwood, Levy, Haugrud, Plott, Robinson, Skiba, Tanha, Zamani, Thannhauser, Glahn, Schmutz, Osorno and Miklas.)

Dry bean (Phaseolus vulgaris L.) production in many regions is threatened by white mold (WM) [Sclerotinia sclerotiorum (Lib.) de Bary]. Seed yield losses can be up to 100% under conditions favorable for the pathogen. The low heritability, polygenic inheritance, and cumbersome screening protocols make it difficult to breed for improved genetic resistance. Some progress in understanding genetic resistance and germplasm improvement has been accomplished, but cultivars with high levels of resistance are yet to be released. A WM multiparent advanced generation inter-cross (MAGIC) population (n = 1060) was developed to facilitate mapping and breeding efforts. A seedling straw test screening method provided a quick assay to phenotype the population for response to WM isolate 1980. Nineteen MAGIC lines were identified with improved resistance. For genome-wide association studies (GWAS), the data was transformed into three phenotypic distributions-quantitative, polynomial, and binomial-and coupled with ∼52,000 single-nucleotide polymorphisms (SNPs). The three phenotypic distributions identified 30 significant genomic intervals [-log 10 (P value) ≥ 3.0]. However, across distributions, four new genomic regions as well as two regions previously reported were found to be associated with resistance. Cumulative R 2 values were 57% for binomial distribution using 13 genomic intervals, 41% for polynomial using eight intervals, and 40% for quantitative using 11 intervals. New resistant germplasm as well as new genomic regions associated with resistance are now available for further investigation., (© 2022 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Soler-Garzón, McClean and Miklas.)

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 F 3 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Soler-Garzón, McClean and Miklas.)

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

Dry bean ( Phaseolus vulgaris L.) is an important worldwide legume crop with low to moderate levels of resistance to common bacterial blight (CBB) caused by Xanthomonas axonopodis pv. phaseoli . A total of 852 genotypes (cultivars, preliminary and advanced breeding lines) from the North Dakota State University dry bean breeding program were tested for their effectiveness as populations for genome-wide association studies (GWAS) to identify genomic regions associated with resistance to CBB, to exploit the associated markers for marker-assisted breeding (MAB), and to identify candidate genes. The genotypes were evaluated in a growth chamber for disease resistance at both the unifoliate and trifoliate stages. At the unifoliate stage, 35% of genotypes were resistant, while 25% of genotypes were resistant at the trifoliate stage. Libraries generated from each genotype were sequenced using the Illumina platform. After filtering for sequence quality, read depth, and minor allele frequency, 41,998 single-nucleotide polymorphisms (SNPs) and 30,285 SNPs were used in GWAS for the Middle American and Andean gene pools, respectively. One region near the distal end of Pv10 near the SAP6 molecular marker from the Andean gene pool explained 26.7-36.4% of the resistance variation. Three to seven regions from the Middle American gene pool contributed to 25.8-27.7% of the resistance, with the most significant peak also near the SAP6 marker. Six of the eight total regions associated with CBB resistance are likely the physical locations of quantitative trait loci identified from previous genetic studies. The two new locations associated with CBB resistance are located at Pv10:22.91-23.36 and Pv11:52.4. A lipoxgenase-1 ortholog on Pv10 emerged as a candidate gene for CBB resistance. The state of one SNP on Pv07 was associated with susceptibility. Its subsequent use in MAB would reduce the current number of lines in preliminary and advanced field yield trial by up to 14% and eliminate only susceptible genotypes. These results provide a foundational SNP data set, improve our understanding of CBB resistance in dry bean, and impact resource allocation within breeding programs as breeding populations may be used for dual purposes: cultivar development as well as genetic studies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Simons, Oladzad, Lamppa, Maniruzzaman, McClean, Osorno and Pasche.)

Background: Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses., Results: Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max., Conclusions: The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Background: Physical seed dormancy is an important trait in legume domestication. Although seed dormancy is beneficial in wild ecosystems, it is generally considered to be an undesirable trait in crops due to reduction in yield and / or quality. The physiological mechanism and underlying genetic factor(s) of seed dormancy is largely unknown in several legume species. Here we employed an integrative approach to understand the mechanisms controlling physical seed dormancy in common bean (Phaseolus vulgaris L.)., Results: Using an innovative CT scan imaging system, we were able to track water movements inside the seed coat. We found that water uptake initiates from the bean seed lens. Using a scanning electron microscopy (SEM) we further identified several micro-cracks on the lens surface of non-dormant bean genotypes. Bulked segregant analysis (BSA) was conducted on a bi-parental RIL (recombinant inbred line) population, segregating for seed dormancy. This analysis revealed that the seed water uptake is associated with a single major QTL on Pv03. The QTL region was fine-mapped to a 118 Kb interval possessing 11 genes. Coding sequence analysis of candidate genes revealed a 5-bp insertion in an ortholog of pectin acetylesterase 8 that causes a frame shift, loss-of-function mutation in non-dormant genotype. Gene expression analysis of the candidate genes in the seed coat of contrasting genotypes indicated 21-fold lower expression of pectin acetylesterase 8 in non-dormant genotype. An analysis of mutational polymorphism was conducted among wild and domesticated beans. Although all the wild beans possessed the functional allele of pectin acetylesterase 8, the majority (77%) of domesticated beans had the non-functional allele suggesting that this variant was under strong selection pressure through domestication., Conclusions: In this study, we identified the physiological mechanism of physical seed dormancy and have identified a candidate allele causing variation in this trait. Our findings suggest that a 5-bp insertion in an ortholog of pectin acetylesterase 8 is likely a major causative mutation underlying the loss of seed dormancy during domestication. Although the results of current study provide strong evidences for the role of pectin acetylesterase 8 in seed dormancy, further confirmations seem necessary by employing transgenic approaches.

Background: In plants, the plasma membrane is enclosed by the cell wall and anchors RLK and RLP proteins, which play a fundamental role in perception of developmental and environmental cues and are crucial in plant development and immunity. These plasma membrane receptors belong to large gene/protein families that are not easily classified computationally. This detailed analysis of these plasma membrane proteins brings a new source of information to the legume genetic, physiology and breeding research communities., Results: A computational approach to identify and classify RLK and RLP proteins is presented. The strategy was evaluated using experimentally-validated RLK and RLP proteins and was determined to have a sensitivity of over 0.85, a specificity of 1.00, and a Matthews correlation coefficient of 0.91. The computational approach can be used to develop a detailed catalog of plasma membrane receptors (by type and domains) in several legume/crop species. The exclusive domains identified in legumes for RLKs are WaaY, APH Pkinase&#95;C, LRR&#95;2, and EGF, and for RLP are L-lectin LPRY and PAN&#95;4. The RLK-nonRD and RLCK subclasses are also discovered by the methodology. In both classes, less than 20% of the total RLK predicted for each species belong to this class. Among the 10-species evaluated ~ 40% of the proteins in the kinome are RLKs. The exclusive legume domain combinations identified are B-Lectin/PR5K domains in G. max, M. truncatula, V. angularis, and V. unguiculata and a three-domain combination B-lectin/S-locus/WAK in C. cajan, M. truncatula, P. vulgaris, V. angularis. and V. unguiculata., Conclusions: The analysis suggests that about 2% of the proteins of each genome belong to the RLK family and less than 1% belong to RLP family. Domain diversity combinations are greater for RLKs compared with the RLP proteins and LRR domains, and the dual domain combination LRR/Malectin were the most frequent domain for both groups of plasma membrane receptors among legume and non-legume species. Legumes exclusively show Pkinase extracellular domains, and atypical domain combinations in RLK and RLP compared with the non-legumes evaluated. The computational logic approach is statistically well supported and can be used with the proteomes of other plant species.

Common bean ( Phaseolus vulgaris L.) production worldwide is hampered by Fusarium root rot (FRR), which is caused by Fusarium solani . Screening for FRR resistance on a large scale is notoriously difficult and often yields inconsistent results due to variability within the environment and pathogen biology. A greenhouse screening assay was developed incorporating multiple isolates of F. solani to improve assay reproducibility. The Andean (ADP; n = 270) and Middle American (MDP; n = 280) Diversity Panels were screened in the greenhouse to identify genetic factors associated with FRR resistance. Forty-seven MDP and 34 ADP lines from multiple market classes were identified as resistant to FRR. Greenhouse phenotyping repeatability was confirmed via five control lines. Genome-wide association mapping using ∼200k SNPs was performed on standard phenotyping score 1-9, as well as binary and polynomial transformation of score data. Sixteen and seven significant genomic regions were identified for ADP and MDP, respectively, using all three classes of phenotypic data. Most candidate genes were associated with plant immune/defense mechanisms. For the ADP population, ortholog of glucan synthase-like enzyme, senescence-associated genes, and NAC domain protein, associated with peak genomic region Pv08:0.04-0.18 Mbp, were the most significant candidate genes. For the MDP population, the peak SNPs Pv07:15.29 Mbp and Pv01:51 Mbp mapped within gene models associated with ethylene response factor 1 and MAC/Perforin domain-containing gene respectively. The research provides a basis for bean improvement through the use of resistant genotypes and genomic regions for more durable root rot resistance., (Copyright © 2020 Zitnick-Anderson, Oladzadabbasabadi, Jain, Modderman, Osorno, McClean and Pasche.)

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., (Copyright © 2020 MacQueen et al.)

Key Message: Two stem rust resistance genes identified on chromosome arms 2BL and 6AL of the cultivated emmer wheat accession PI 193883 can be used for protecting modern varieties against Ug99 strains. The wheat research community consistently strives to identify new genes that confer resistance to stem rust caused by the fungal pathogen Puccinia graminis f. sp. tritici Eriks & E. Henn (Pgt). In the current study, our objective was to identify and genetically characterize the stem rust resistance derived from the cultivated emmer accession PI 193883. A recombinant inbred line population developed from a cross between the susceptible durum wheat line Rusty and PI 193883 was genotyped and evaluated for reaction to Pgt races TTKSK, TRTTF, and TMLKC. Two QTLs conferring resistance were identified on chromosome arms 2BL (QSr.fcu-2B) and 6AL (QSr.fcu-6A). The stem rust resistance gene (Sr883-2B) underlying QSr.fcu-2B was recessive, and based on its physical location it is located proximal to the Sr9 region. QSr.fcu-6A was located in the Sr13 region, but PI 193883 is known to carry the susceptible haplotype S4 for Sr13, indicating that the gene underlying QSr.fcu-6A (Sr883-6A) is likely a new allele of Sr13 or a gene residing close to Sr13. Three IWGSC scaffold-based simple sequence repeat (SSR) and two SNP-based semi-thermal asymmetric reverse PCR (STARP) markers were developed for the Sr883-2B region, and one STARP marker was developed for Sr883-6A. Sr883-2B was epistatic to Sr883-6A for reaction to TTKSK and TRTTF, and the two genes had additive effects for TMLKC. These two genes and the markers developed in this research provide additional resources and tools for the improvement in stem rust resistance in durum and common wheat breeding programs.

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&#95;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.

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&#95;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.

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., (Copyright © 2019 Oladzad et al.)

Genetic mutations in genes governing wheat threshability were critical for domestication. Knowing when these genes mutated during wheat evolution will provide more insight into the domestication process and lead to further exploitation of primitive alleles for wheat improvement. We evaluated a population of recombinant inbred lines derived from a cross between the durum variety Rusty and the cultivated emmer accession PI 193883 for threshability, rachis fragility, and other spike-related traits. Quantitative trait loci (QTL) associated with spike length, spikelets per spike, and spike compactness were primarily associated with known genes such as the pleiotropic domestication gene Q. Interestingly, rachis fragility was not associated with the Q locus, suggesting that this trait, usually a pleiotropic effect of the q allele, can be influenced by the genetic background. Threshability QTL were identified on chromosome arms 2AS, 2BS, and 5AL corresponding to the tenacious glume genes Tg 2A and Tg 2B as well as the Q gene, respectively, further demonstrating that cultivated emmer harbors the primitive non-free-threshing alleles at all three loci. Genetic analysis indicated that the effects of the three genes are mostly additive, with Q having the most profound effects on threshability, and that free-threshing alleles are necessary at all three loci to attain a completely free-threshing phenotype. These findings provide further insight into the timeline and possible pathways of wheat domestication and evolution that led to the formation of modern day domesticated wheats.

Common bean (Phaseolus vulgaris L.) is an important high protein crop grown worldwide. North Dakota and Minnesota are the largest producers of common beans in the USA, but crop production is threatened by soybean cyst nematode (SCN; Heterodera glycines Ichinohe) because most current cultivars are susceptible. Greenhouse screening data using SCN HG type 0 from 317 plant introductions (PI's) from the USDA core collection was used to conduct a genome wide association study (GWAS). These lines were divided into two subpopulations based on principal component analysis (Middle American vs. Andean). Phenotypic results based on the female index showed that accessions could be classified as highly resistant (21% and 27%), moderately resistant (51% and 48%), moderately susceptible (27% and 22%) and highly susceptible (1% and 3%) for Middle American and Andean gene pools, respectively. Mixed models with two principal components (PCs) and kinship matrix for Middle American genotypes and Andean genotypes were used in the GWAS analysis using 3,985 and 4,811 single nucleotide polymorphic (SNP) markers, respectively which were evenly distributed across all 11 chromosomes. Significant peaks on Pv07, and Pv11 in Middle American and on Pv07, Pv08, Pv09 and Pv11 in Andean group were found to be associated with SCN resistance. Homologs of soybean rhg1, a locus which confers resistance to SCN in soybean, were identified on chromosomes Pv01 and Pv08 in the Middle American and Andean gene pools, respectively. These genomic regions may be the key to develop SCN-resistant common bean cultivars., Competing Interests: The authors have declared that no competing interests exist.

Cochliobolus victoria, the causal agent of Victoria blight, is pathogenic due to its production of a toxin called victorin. Victorin sensitivity in oats, barley, Brachypodium spp., and Arabidopsis has been associated with nucleotide-binding site leucine-rich repeat (NLR) genes, a class of genes known for conferring disease resistance. In this work, we investigated the sensitivity of Phaseolus vulgaris to victorin. We found that victorin sensivity in Phaseolus vulgaris is a developmentally regulated, quantitative trait. A single quantitative trait locus (QTL) accounted for 34% of the phenotypic variability in victorin sensitivity among Stampede × Red Hawk (S×R) recombinant inbred lines. We cloned two NLR-encoding genes within this QTL and showed one, Phvul05G031200 (PvLOV), confers victorin-dependent cell death when overexpressed in Nicotiana benthamiana. Protein sequences of PvLOV from victorin-sensitive and the victorin-resistant bean parents differ by two amino acids in the leucine-rich repeat region, but both proteins confer victorin-dependent cell death when overexpressed in N. benthamiana.

The presence of seed color in common bean (Phaseolus vulgaris) requires the dominant-acting P (pigment) gene, and white seed is a recessive phenotype in all domesticated races of the species. P was classically associated with seed size, thus describing it as the first genetic marker for a quantitative trait. The molecular structure of P was characterized to understand the selection of white seeds during bean diversification and the relationship of P to seed weight. P was identified by homology searches, a genome-wide association study (GWAS) and gene remodeling, and confirmed by gene silencing. Allelic variation was assessed by a combination of resequencing and marker development, and the relationship between P and seed weight was assessed by a GWAS study. P is a member of clade B of subclass IIIf of plant basic helix-loop-helix (bHLH) proteins. Ten race-specific P alleles conditioned the white seed phenotype, and each causative mutation affected at least one bHLH domain required for color expression. GWAS analysis confirmed the classic association of P with seed weight. In common bean, white seeds are the result of convergent evolution and, among plant species, orthologous convergence on a single transcription factor gene was observed., (© 2018 The Authors. New Phytologist © 2018 New Phytologist Trust.)

Common bean ( L.) is the most consumed edible grain legume worldwide and contains a wide range of nutrients for human health including dietary fiber. Diets high in beans are associated with lower rates of chronic diseases such as obesity and type 2 diabetes, and the content of dietary fibers varies among different market classes of dry bean. In this study, we evaluated the dietary fiber content in a Middle American diversity panel (MDP) of common bean and evaluated the genetic architecture of the various dietary fiber components. The dietary fiber components included insoluble and soluble dietary fibers as well as the antinutritional raffinose family of oligosaccharides (RFOs; raffinose, stachyose, and verbascose). All variables measured differed among market classes and entries. Colored bean seeds had higher levels of insoluble dietary fibers with the black market class showing also the highest raffinose and stachyose content. Cultivars and lines released since 1997 had higher insoluble dietary fibers and RFO content in race Durango. Higher levels of RFOs were also observed in cultivars with type II growth habit that was a recent breeding target in Durango race germplasm. Candidate genes for dietary fiber traits, especially homologs to two main genes in the RFO biosynthesis pathway, were identified. The knowledge of diversity of dietary fibers in the MDP accompanied with the identification of candidate genes could effectively improve dietary fiber components in common bean., (Copyright © 2018 Crop Science Society of America.)

Flooding is a devastating abiotic stress that endangers crop production in the twenty-first century. Because of the severe susceptibility of common bean ( Phaseolus vulgaris L.) to flooding, an understanding of the genetic architecture and physiological responses of this crop will set the stage for further improvement. However, challenging phenotyping methods hinder a large-scale genetic study of flooding tolerance in common bean and other economically important crops. A greenhouse phenotyping protocol was developed to evaluate the flooding conditions at early stages. The Middle-American diversity panel ( n = 272) of common bean was developed to capture most of the diversity exits in North American germplasm. This panel was evaluated for seven traits under both flooded and non-flooded conditions at two early developmental stages. A subset of contrasting genotypes was further evaluated in the field to assess the relationship between greenhouse and field data under flooding condition. A genome-wide association study using ~150 K SNPs was performed to discover genomic regions associated with multiple physiological responses. The results indicate a significant strong correlation ( r > 0.77) between greenhouse and field data, highlighting the reliability of greenhouse phenotyping method. Black and small red beans were the least affected by excess water at germination stage. At the seedling stage, pinto and great northern genotypes were the most tolerant. Root weight reduction due to flooding was greatest in pink and small red cultivars. Flooding reduced the chlorophyll content to the greatest extent in the navy bean cultivars compared with other market classes. Races of Durango/Jalisco and Mesoamerica were separated by both genotypic and phenotypic data indicating the potential effect of eco-geographical variations. Furthermore, several loci were identified that potentially represent the antagonistic pleiotropy. The GWAS analysis revealed peaks at Pv08/1.6 Mb and Pv02/41 Mb that are associated with root weight and germination rate, respectively. These regions are syntenic with two QTL reported in soybean ( Glycine max L.) that contribute to flooding tolerance, suggesting a conserved evolutionary pathway involved in flooding tolerance for these related legumes.

Common bean ( Phaseolus vulgaris L.) is an important legume, useful for its high protein and dietary fiber. The fungal pathogen Uromyces appendiculatus (Pers.) Unger can cause major loss in susceptible varieties of the common bean. The Ur-3 locus provides race specific resistance to virulent strains or races of the bean rust pathogen along with Crg , (Complements resistance gene), which is required for Ur-3 -mediated rust resistance. In this study, we inoculated two common bean genotypes (resistant "Sierra" and susceptible crg) with rust race 53 of U. appendiculatus , isolated leaf RNA at specific time points, and sequenced their transcriptomes. First, molecular markers were used to locate and identify a 250 kb deletion on chromosome 10 in mutant crg (which carries a deletion at the Crg locus). Next, we identified differential expression of several disease resistance genes between Mock Inoculated (MI) and Inoculated (I) samples of "Sierra" leaf RNA within the 250 kb delineated region. Both marker assisted molecular profiling and RNA-seq were used to identify possible transcriptomic locations of interest regarding the resistance in the common bean to race 53. Identification of differential expression among samples in disease resistance clusters in the bean genome may elucidate significant genes underlying rust resistance. Along with preserving favorable traits in the crop, the current research may also aid in global sustainability of food stocks necessary for many populations., Competing Interests: The authors declare no conflict of interest.

Common bean ( L.) breeding programs aim to improve both agronomic and seed characteristics traits. However, the genetic architecture of the many traits that affect common bean production are not completely understood. Genome-wide association studies (GWAS) provide an experimental approach to identify genomic regions where important candidate genes are located. A panel of 280 modern bean genotypes from race Mesoamerica, referred to as the Middle American Diversity Panel (MDP), were grown in four US locations, and a GWAS using >150,000 single-nucleotide polymorphisms (SNPs) (minor allele frequency [MAF] ≥ 5%) was conducted for six agronomic traits. The degree of inter- and intrachromosomal linkage disequilibrium (LD) was estimated after accounting for population structure and relatedness. The LD varied between chromosomes for the entire MDP and among race Mesoamerica and Durango-Jalisco genotypes within the panel. The LD patterns reflected the breeding history of common bean. Genome-wide association studies led to the discovery of new and known genomic regions affecting the agronomic traits at the entire population, race, and location levels. We observed strong colocalized signals in a narrow genomic interval for three interrelated traits: growth habit, lodging, and canopy height. Overall, this study detected ∼30 candidate genes based on a priori and candidate gene search strategies centered on the 100-kb region surrounding a significant SNP. These results provide a framework from which further research can begin to understand the actual genes controlling important agronomic production traits in common bean., (Copyright © 2016 Crop Science Society of America.)

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., (Copyright © 2016 Crop Science Society of America.)

A total of 992,682 single-nucleotide polymorphisms (SNPs) was identified as ideal for Illumina Infinium II BeadChip design after sequencing a diverse set of 17 common bean (Phaseolus vulgaris L) varieties with the aid of next-generation sequencing technology. From these, two BeadChips each with >5000 SNPs were designed. The BARCBean6K&#95;1 BeadChip was selected for the purpose of optimizing polymorphism among market classes and, when possible, SNPs were targeted to sequence scaffolds in the Phaseolus vulgaris 14× genome assembly with sequence lengths >10 kb. The BARCBean6K&#95;2 BeadChip was designed with the objective of anchoring additional scaffolds and to facilitate orientation of large scaffolds. Analysis of 267 F2 plants from a cross of varieties Stampede × Red Hawk with the two BeadChips resulted in linkage maps with a total of 7040 markers including 7015 SNPs. With the linkage map, a total of 432.3 Mb of sequence from 2766 scaffolds was anchored to create the Phaseolus vulgaris v1.0 assembly, which accounted for approximately 89% of the 487 Mb of available sequence scaffolds of the Phaseolus vulgaris v0.9 assembly. A core set of 6000 SNPs (BARCBean6K&#95;3 BeadChip) with high genotyping quality and polymorphism was selected based on the genotyping of 365 dry bean and 134 snap bean accessions with the BARCBean6K&#95;1 and BARCBean6K&#95;2 BeadChips. The BARCBean6K&#95;3 BeadChip is a useful tool for genetics and genomics research and it is widely used by breeders and geneticists in the United States and abroad., (Copyright © 2015 Song et al.)

Gene expression related to drought response in the leaf tissues of two Brazilian upland cultivars, the drought-tolerant Douradão and the drought-sensitive Primavera, was analyzed. RNA-seq identified 27,618 transcripts in the Douradão cultivar, with 24,090 (87.2%) homologous to the rice database, and 27,221 transcripts in the Primavera cultivar, with 23,663 (86.9%) homologous to the rice database. Gene-expression analysis between control and water-deficient treatments revealed 493 and 1154 differentially expressed genes in Douradão and Primavera cultivars, respectively. Genes exclusively expressed under drought were identified for Douradão, including two genes of particular interest coding for the protein peroxidase precursor, which is involved in three distinct metabolic pathways. Comparisons between the two drought-exposed cultivars revealed 2314 genes were differentially expressed (978 upregulated, 1336 downregulated in Douradão). Six genes distributed across 4 different transcription factor families (bHLH, MYB, NAC, and WRKY) were identified, all of which were upregulated in Douradão compared to Primavera during drought. Most of the genes identified in Douradão activate metabolic pathways responsible for production of secondary metabolites and genes coding for enzymatically active signaling receptors. Quantitative PCR validation showed that most gene expression was in agreement with computational prediction of these transcripts. The transcripts identified here will define molecular markers for identification of Cis-acting elements to search for allelic variants of these genes through analysis of polymorphic SNPs in GenBank accessions of upland rice, aiming to develop cultivars with the best combination of these alleles, resulting in materials with high yield potential in the event of drought during the reproductive phase.

In wheat, exotic genotypes harbor a broad range of spike-related traits, and can be used as a source of new genes for germplasm enhancement in wheat breeding programs. In the present study, a population of 163 recombinant inbred lines was derived from a cross between an elite line (WCB414) and an exotic line (WCB617) with branched spike (supernumerary spikelet; SS) head morphology. The population was evaluated over four to six environments to identify quantitative trait loci (QTL) associated with nine spike-related traits and 10 agronomic traits. A genetic map consisting of 939 diversity arrays technology (DArT) markers was constructed. Composite interval mapping identified a total of 143 QTL located on 17 different wheat chromosomes and included 33 consistent and definitive QTL. The amount of phenotype variation explained (PVE) by individual QTL ranged from 0.61 to 91.8%. One major QTL for glume pubescence was located in a QTL-rich region on the short arm of chromosome 1A, where loci for other traits such as for kernels per spike (KS) and spike length (SL) were also identified. Similarly, a cluster of QTL associated with yield-related, agronomic and spike-related traits contributing up to 40.3% of PVE was found on the short arm of chromosome 2D, in the vicinity of a major QTL for SS-related traits. Consistent and major QTL identified in the present study may be useful in marker-assisted breeding programs to facilitate transfer of desirable alleles into other germplasm. Desirable QTL alleles were also contributed by the exotic line, suggesting the possibility of enriching the breeding germplasm with alleles from SS genotypes., (© 2015 The Authors.)

Leaf rust, caused by Puccinia triticina (Pt), and stripe rust, caused by P. striiformis f. sp. tritici (Pst), are destructive foliar diseases of wheat worldwide. Breeding for disease resistance is the preferred strategy of managing both diseases. The continued emergence of new races of Pt and Pst requires a constant search for new sources of resistance. Here we report a genome-wide association analysis of 567 winter wheat (Triticum aestivum) landrace accessions using the Infinium iSelect 9K wheat SNP array to identify loci associated with seedling resistance to five races of Pt (MDCL, MFPS, THBL, TDBG, and TBDJ) and one race of Pst (PSTv-37) frequently found in the Northern Great Plains of the United States. Mixed linear models identified 65 and eight significant markers associated with leaf rust and stripe rust, respectively. Further, we identified 31 and three QTL associated with resistance to Pt and Pst, respectively. Eleven QTL, identified on chromosomes 3A, 4A, 5A, and 6D, are previously unknown for leaf rust resistance in T. aestivum.

Key Message: A population developed from an exotic line with supernumerary spikelets was genetically dissected for eight quality traits, discovering new genes/alleles with potential use in wheat breeding programs. Identifying new QTLs and alleles in exotic germplasm is paramount for further improvement of quality traits in wheat. In the present study, an RIL population developed from a cross of an elite wheat line (WCB414) and an exotic genotype with supernumerary spikelets (SS) was used to identify QTLs and new alleles for eight quality traits. Composite interval mapping for 1,000 kernels weight (TKW), kernel volume weight (KVW), grain protein content (GPC), percent of flour extraction (FE) and four mixograph-related traits identified a total of 69 QTLs including 19 stable QTLs. These QTLs were located on 18 different chromosomes (except 4D, 5D, and 6D). Thirteen of these QTLs explained more than 15% of phenotypic variation (PV) and were considered as major QTLs. In this study, we identified 11 QTLs for TKW (R (2) = 7.2-17.1 %), 10 for KVW (R (2) = 6.7-22.5%), 11 for GPC (R (2) = 4.7-16.9%), 6 for FE (R (2) = 4.8-19%) and 31 for mixograph-related traits (R (2) = 3.2-41.2%). In this population, several previously identified QTLs for SS, nine spike-related and ten agronomic traits were co-located with the quality QTLs, suggesting pleiotropic effects or close linkage among loci. The traits GPC and mixogram-related traits were positively correlated with SS. Indeed, several loci for quality traits were co-located with QTL for SS. The exotic parent contributed positive alleles that increased PV of the traits at 56% of loci demonstrating the suitability of germplasm with SS to improve quality traits in wheat.

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&#95;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.

Background: Common bean (Phaseolus vulgaris) is grown throughout the world and comprises roughly 50% of the grain legumes consumed worldwide. Despite this, genetic resources for common beans have been lacking. Next generation sequencing, has facilitated our investigation of the gene expression profiles associated with biologically important traits in common bean. An increased understanding of gene expression in common bean will improve our understanding of gene expression patterns in other legume species., Results: Combining recently developed genomic resources for Phaseolus vulgaris, including predicted gene calls, with RNA-Seq technology, we measured the gene expression patterns from 24 samples collected from seven tissues at developmentally important stages and from three nitrogen treatments. Gene expression patterns throughout the plant were analyzed to better understand changes due to nodulation, seed development, and nitrogen utilization. We have identified 11,010 genes differentially expressed with a fold change ≥ 2 and a P-value < 0.05 between different tissues at the same time point, 15,752 genes differentially expressed within a tissue due to changes in development, and 2,315 genes expressed only in a single tissue. These analyses identified 2,970 genes with expression patterns that appear to be directly dependent on the source of available nitrogen. Finally, we have assembled this data in a publicly available database, The Phaseolus vulgaris Gene Expression Atlas (Pv GEA), http://plantgrn.noble.org/PvGEA/ . Using the website, researchers can query gene expression profiles of their gene of interest, search for genes expressed in different tissues, or download the dataset in a tabular form., Conclusions: These data provide the basis for a gene expression atlas, which will facilitate functional genomic studies in common bean. Analysis of this dataset has identified genes important in regulating seed composition and has increased our understanding of nodulation and impact of the nitrogen source on assimilation and distribution throughout the plant.

Iron deficiency chlorosis (IDC) is a yield limiting problem in soybean (Glycine max (L.) Merr) production regions with calcareous soils. Genome-wide association study (GWAS) was performed using a high density SNP map to discover significant markers, QTL and candidate genes associated with IDC trait variation. A stepwise regression model included eight markers after considering LD between markers, and identified seven major effect QTL on seven chromosomes. Twelve candidate genes known to be associated with iron metabolism mapped near these QTL supporting the polygenic nature of IDC. A non-synonymous substitution with the highest significance in a major QTL region suggests soybean orthologs of FRE1 on Gm03 is a major gene responsible for trait variation. NAS3, a gene that encodes the enzyme nicotianamine synthase which synthesizes the iron chelator nicotianamine also maps to the same QTL region. Disease resistant genes also map to the major QTL, supporting the hypothesis that pathogens compete with the plant for Fe and increase iron deficiency. The markers and the allelic combinations identified here can be further used for marker assisted selection.

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.

Next generation sequence data provides valuable information and tools for genetic and genomic research and offers new insights useful for marker development. This data is useful for the design of accurate and user-friendly molecular tools. Common bean (Phaseolus vulgaris L.) is a diverse crop in which separate domestication events happened in each gene pool followed by race and market class diversification that has resulted in different morphological characteristics in each commercial market class. This has led to essentially independent breeding programs within each market class which in turn has resulted in limited within market class sequence variation. Sequence data from selected genotypes of five bean market classes (pinto, black, navy, and light and dark red kidney) were used to develop InDel-based markers specific to each market class. Design of the InDel markers was conducted through a combination of assembly, alignment and primer design software using 1.6× to 5.1× coverage of Illumina GAII sequence data for each of the selected genotypes. The procedure we developed for primer design is fast, accurate, less error prone, and higher throughput than when they are designed manually. All InDel markers are easy to run and score with no need for PCR optimization. A total of 2687 InDel markers distributed across the genome were developed. To highlight their usefulness, they were employed to construct a phylogenetic tree and a genetic map, showing that InDel markers are reliable, simple, and accurate.

Legumes contain a variety of phytochemicals derived from the phenylpropanoid pathway that have important effects on human health as well as seed coat color, plant disease resistance and nodulation. However, the information about the genes involved in this important pathway is fragmentary in common bean (Phaseolus vulgaris L.). The objectives of this research were to isolate genes that function in and control the phenylpropanoid pathway in common bean, determine their genomic locations in silico in common bean and soybean, and analyze sequences of the 4CL gene family in two common bean genotypes. Sequences of phenylpropanoid pathway genes available for common bean or other plant species were aligned, and the conserved regions were used to design sequence-specific primers. The PCR products were cloned and sequenced and the gene sequences along with common bean gene-based (g) markers were BLASTed against the Glycine max v.1.0 genome and the P. vulgaris v.1.0 (Andean) early release genome. In addition, gene sequences were BLASTed against the OAC Rex (Mesoamerican) genome sequence assembly. In total, fragments of 46 structural and regulatory phenylpropanoid pathway genes were characterized in this way and placed in silico on common bean and soybean sequence maps. The maps contain over 250 common bean g and SSR (simple sequence repeat) markers and identify the positions of more than 60 additional phenylpropanoid pathway gene sequences, plus the putative locations of seed coat color genes. The majority of cloned phenylpropanoid pathway gene sequences were mapped to one location in the common bean genome but had two positions in soybean. The comparison of the genomic maps confirmed previous studies, which show that common bean and soybean share genomic regions, including those containing phenylpropanoid pathway gene sequences, with conserved synteny. Indels identified in the comparison of Andean and Mesoamerican common bean 4CL gene sequences might be used to develop inter-pool phenylpropanoid pathway gene-based markers. We anticipate that the information obtained by this study will simplify and accelerate selections of common bean with specific phenylpropanoid pathway alleles to increase the contents of beneficial phenylpropanoids in common bean and other legumes.

A small fast neutron (FN) mutant population has been established from Phaseolus vulgaris cv. Red Hawk. We leveraged the available P. vulgaris genome sequence and high throughput next generation DNA sequencing to examine the genomic structure of five P. vulgaris cv. Red Hawk FN mutants with striking visual phenotypes. Analysis of these genomes identified three classes of structural variation (SV); between cultivar variation, natural variation within the FN mutant population, and FN induced mutagenesis. Our analyses focused on the latter two classes. We identified 23 large deletions (>40 bp) common to multiple individuals, illustrating residual heterogeneity and regions of SV within the common bean cv. Red Hawk. An additional 18 large deletions were identified in individual mutant plants. These deletions, ranging in size from 40 bp to 43,000 bp, are potentially the result of FN mutagenesis. Six of the 18 deletions lie near or within gene coding regions, identifying potential candidate genes causing the mutant phenotype.

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π(sil)=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F(st)=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

Necrotrophic effectors (also known as host-selective toxins) are important determinants of disease in the wheat-Stagonospora nodorum pathosystem. To date, five necrotrophic effector-host gene interactions have been identified in this system. Most of these interactions have additive effects while some are epistatic. The Snn4-SnTox4 interaction was originally identified in a recombinant-inbred population derived from a cross between the Swiss winter wheat cultivars 'Arina' and 'Forno' using the S. nodorum isolate Sn99CH 1A7a. Here, we used a recombinant-inbred population consisting of 121 lines developed from a cross between the hexaploid land race Salamouni and the hexaploid wheat 'Katepwa' (SK population). The SK population was used for the construction of linkage maps and quantitative trait loci (QTL) detection using the Swiss S. nodorum isolate Sn99CH 1A7a. The linkage maps developed in the SK population spanned 3,228 centimorgans (cM) and consisted of 441 simple-sequence repeats, 9 restriction fragment length polymorphisms, 29 expressed sequence tag sequence-tagged site markers, and 5 phenotypic markers. The average marker density was 6.7 cM/marker. Two QTL, designated QSnb.fcu-1A and QSnb.fcu-7A on chromosome arms 1AS and 7AS, respectively, were associated with disease caused by the S. nodorum isolate Sn99CH 1A7a. The effects of QSnb.fcu-1A were determined by the Snn4-SnTox4 interaction and accounted for 23.5% of the phenotypic variation in this population, whereas QSnb.fcu-7A accounted for 16.4% of the phenotypic variation for disease but was not associated with any known effector sensitivity locus. The effects of both QTL were largely additive and collectively accounted for 35.7% of the total phenotypic variation. The results of this research validate the effects of a compatible Snn4-SnTox4 interaction in a different genetic background, and it provides knowledge regarding genomic regions and molecular markers that can be used to improve Stagonospora nodorum blotch resistance in wheat germplasm.

Climate change and global population increase are two converging forces that will jointly challenge researchers to design programs that ensure crop production systems meet the world's food demand. Climate change will potentially reduce productivity while a global population increase will require more food. If productivity is not improved for future climatic conditions, food insecurity may foster major economic and political uncertainty. Given the importance of grain legumes in general - common bean (Phaseolus vulgaris L.) in particular - a workshop entitled 'Improving Tolerance of Common Bean to Abiotic Stresses' was held with the goal of developing an interdisciplinary research agenda designed to take advantage of modern genotyping and breeding approaches that are coupled with large scale phenotyping efforts to improve common bean. Features of the program included a multinational phenotyping effort to evaluate the major common bean core germplasm collections and appropriate genetic populations. The phenotyping effort will emphasise the response of root and shoot traits to individual and combined stress conditions. These populations would also be genotyped using newly emerging high density single nucleotide polymorphism (SNP) marker arrays or next generation sequencing technology. Association analysis of the core collections aims to identify key loci associated with the response to the stress conditions. Companion bi-parental quantitative trait loci (QTL) experiments will act as confirmation experiments for the association analysis. The upcoming release of the genome sequence of common bean will be leveraged by utilising population genomic approaches to discover genomic regions that differentiate stress-responsive and non-responsive genotypes. The genome sequence will also enable global gene expression studies that will highlight specific molecular-based stress responses. This collective knowledge will inform the selection of parental lines to improve the efficiency of common bean improvement programs.

Common bacterial blight (CBB), incited by Xanthomonas axonopodis pv. phaseoli (Xap), is a serious seed-borne disease of common bean (Phaseolus vulgaris L.) in both temperate and tropical production zones. The line HR45 is highly resistant to Xap infection on leaves and pods in both field and greenhouse. To understand the molecular mechanisms underlying CBB resistance in HR45, cDNA-amplified fragment length polymorphism (AFLP) technique was used to identify the genes that are differentially expressed in the leaves of HR45 at different time-periods after inoculation. Selective amplifications with 34 primer combinations allowed the visualization of 2,448 transcript-derived fragments (TDFs) in infected leaves, and 259 (10.6%) of them were differentially expressed TDFs (DE-TDFs). Seventy-seven of the DE-TDFs were cloned and sequenced. Thirty-nine of the 77 (50.6%) DE-TDFs representing bean transcripts were not previously reported in any EST database. The expression patterns of 10 representative DE-TDFs were further confirmed by real-time RT-PCR. BLAST analysis suggested that 40% (31 of 77) of the DE-TDFs were homologous to the genes related to metabolism, photosynthesis, and cellular transport, whereas 28% (22 of 77) of the DE-TDFs showed homology to the genes involved in defence response, response to stimulus, enzyme regulation, and transcription regulation. Thus, the 22 pathogenesis-related DE-TDFs were selected as potential functional candidate genes (FCGs) in association with CBB resistance. Meanwhile, six of the DE-TDFs (1FCG and five other DE-TDFs) were in silico mapped to the distal region of the bean linkage group B6 (the genomic location containing the major CBB resistance QTL in HR45) and, therefore, were considered as positional candidate genes (PCGs). This study represents a first step towards the discovery of bean genes expressed upon Xap infection. This information will be useful for elucidating the molecular basis of the resistance response process and identifying the genes that underlie the CBB-resistance.

Map-based cloning and fine mapping to find genes of interest and marker assisted selection (MAS) requires good genetic maps with reproducible markers. In this study, we saturated the linkage map of the intra-gene pool population of common bean DOR364 × BAT477 (DB) by evaluating 2,706 molecular markers including SSR, SNP, and gene-based markers. On average the polymorphism rate was 7.7% due to the narrow genetic base between the parents. The DB linkage map consisted of 291 markers with a total map length of 1,788 cM. A consensus map was built using the core mapping populations derived from inter-gene pool crosses: DOR364 × G19833 (DG) and BAT93 × JALO EEP558 (BJ). The consensus map consisted of a total of 1,010 markers mapped, with a total map length of 2,041 cM across 11 linkage groups. On average, each linkage group on the consensus map contained 91 markers of which 83% were single copy markers. Finally, a synteny analysis was carried out using our highly saturated consensus maps compared with the soybean pseudo-chromosome assembly. A total of 772 marker sequences were compared with the soybean genome. A total of 44 syntenic blocks were identified. The linkage group Pv6 presented the most diverse pattern of synteny with seven syntenic blocks, and Pv9 showed the most consistent relations with soybean with just two syntenic blocks. Additionally, a co-linear analysis using common bean transcript map information against soybean coding sequences (CDS) revealed the relationship with 787 soybean genes. The common bean consensus map has allowed us to map a larger number of markers, to obtain a more complete coverage of the common bean genome. Our results, combined with synteny relationships provide tools to increase marker density in selected genomic regions to identify closely linked polymorphic markers for indirect selection, fine mapping or for positional cloning.

Determinacy is an agronomically important trait associated with the domestication in soybean (Glycine max). Most soybean cultivars are classifiable into indeterminate and determinate growth habit, whereas Glycine soja, the wild progenitor of soybean, is indeterminate. Indeterminate (Dt1/Dt1) and determinate (dt1/dt1) genotypes, when mated, produce progeny that segregate in a monogenic pattern. Here, we show evidence that Dt1 is a homolog (designated as GmTfl1) of Arabidopsis terminal flower 1 (TFL1), a regulatory gene encoding a signaling protein of shoot meristems. The transition from indeterminate to determinate phenotypes in soybean is associated with independent human selections of four distinct single-nucleotide substitutions in the GmTfl1 gene, each of which led to a single amino acid change. Genetic diversity of a minicore collection of Chinese soybean landraces assessed by simple sequence repeat (SSR) markers and allelic variation at the GmTfl1 locus suggest that human selection for determinacy took place at early stages of landrace radiation. The GmTfl1 allele introduced into a determinate-type (tfl1/tfl1) Arabidopsis mutants fully restored the wild-type (TFL1/TFL1) phenotype, but the Gmtfl1 allele in tfl1/tfl1 mutants did not result in apparent phenotypic change. These observations indicate that GmTfl1 complements the functions of TFL1 in Arabidopsis. However, the GmTfl1 homeolog, despite its more recent divergence from GmTfl1 than from Arabidopsis TFL1, appears to be sub- or neo-functionalized, as revealed by the differential expression of the two genes at multiple plant developmental stages and by allelic analysis at both loci.