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

Mapping yield and yield-related traits using diverse common bean germplasm

Common bean (bean) is one of the most important legume crops, and mapping genes for yield and yield-related traits is essential for its improvement. However, yield is a complex trait that is typically controlled by many loci in crop genomes. The objective of this research was to identify regions in the bean genome associated with yield and a number of yield-related traits using a collection of 121 diverse bean genotypes with different yields. The beans were evaluated in replicated trials at two locations, over two years. Significant variation among genotypes was identified for all traits analyzed in the four environments. The collection was genotyped with the BARCBean6K_3 chip (5,398 SNPs), two yield/antiyield gene-based markers, and seven markers previously associated with resistance to common bacterial blight (CBB), including a Niemann-Pick polymorphism (NPP) gene-based marker. Over 90% of the single-nucleotide polymorphisms (SNPs) were polymorphic and separated the panel into two main groups of small-seeded and large-seeded beans, reflecting their Mesoamerican and Andean origins. Thirty-nine significant marker-trait associations (MTAs) were identified between 31 SNPs and 15 analyzed traits on all 11 bean chromosomes. Some of these MTAs confirmed genome regions previously associated with the yield and yield-related traits in bean, but a number of associations were not reported previously, especially those with derived traits. Over 600 candidate genes with different functional annotations were identified for the analyzed traits in the 200-Kb region centered on significant SNPs. Fourteen SNPs were identified within the gene model sequences, and five additional SNPs significantly associated with five different traits were located at less than 0.6 Kb from the candidate genes. The work confirmed associations between two yield/antiyield gene-based markers (AYD1m and AYD2m) on chromosome Pv09 with yield and identified their association with a number of yield-related traits, including seed weight. The results also confirmed the usefulness of the NPP marker in screening for CBB resistance. Since disease resistance and yield measurements are environmentally dependent and labor-intensive, the three gene-based markers (CBB- and two yield-related) and quantitative trait loci (QTL) that were validated in this work may be useful tools for simplifying and accelerating the selection of high-yielding and CBB-resistant bean cultivars.

Novel candidate loci for morpho-agronomic and seed quality traits detected by targeted genotyping-by-sequencing in common bean

Phaseolus vulgaris L., known as common bean, is one of the most important grain legumes cultivated around the world for its immature pods and dry seeds, which are rich in protein and micronutrients. Common bean offers a cheap food and protein sources to ameliorate food shortage and malnutrition around the world. However, the genetic basis of most important traits in common bean remains unknown. This study aimed at identifying QTL and candidate gene models underlying twenty-six agronomically important traits in common bean. For this, we assembled and phenotyped a diversity panel of 200 P. vulgaris genotypes in the greenhouse, comprising determinate bushy, determinate climbing and indeterminate climbing beans. The panel included dry beans and snap beans from different breeding programmes, elite lines and landraces from around the world with a major focus on accessions of African, European and South American origin. The panel was genotyped using a cost-conscious targeted genotyping-by-sequencing (GBS) platform to take advantage of highly polymorphic SNPs detected in previous studies and in diverse germplasm. The detected single nucleotide polymorphisms (SNPs) were applied in marker-trait analysis and revealed sixty-two quantitative trait loci (QTL) significantly associated with sixteen traits. Gene model identification via a similarity-based approach implicated major candidate gene models underlying the QTL associated with ten traits including, flowering, yield, seed quality, pod and seed characteristics. Our study revealed six QTL for pod shattering including three new QTL potentially useful for breeding. However, the panel was evaluated in a single greenhouse environment and the findings should be corroborated by evaluations across different field environments. Some of the detected QTL and a number of candidate gene models only elucidate the understanding of the genetic nature of these traits and provide the basis for further studies. Finally, the study showed the possibility of using a limited number of SNPs in performing marker-trait association in common bean by applying a highly scalable targeted GBS approach. This targeted GBS approach is a cost-efficient strategy for assessment of the genetic basis of complex traits and can enable geneticists and breeders to identify novel loci and targets for marker-assisted breeding more efficiently.

Genome-wide association studies dissect the genetic architecture of seed shape and size in common bean

Seed weight and size are important yield components. Thus, selecting for large seeds has been a key objective in crop domestication and breeding. In common bean, seed shape is also important since it influences industrial processing and plays a vital role in determining the choices of consumers and farmers. In this study, we performed genome-wide association studies on a core collection of common bean accessions to dissect the genetic architecture and identify genomic regions associated with seed morphological traits related to weight, size, and shape. Phenotypic data were collected by high-throughput image-based approaches, and utilized to test associations with 10,362 single-nucleotide polymorphism markers using multilocus mixed models. We searched within genome-associated regions for candidate genes putatively involved in seed phenotypic variation. The collection exhibited high variability for the entire set of seed traits, and the Andean gene pool was found to produce larger, heavier seeds than the Mesoamerican gene pool. Strong pairwise correlations were verified for most seed traits. Genome-wide association studies identified marker-trait associations accounting for a considerable amount of phenotypic variation in length, width, projected area, perimeter, and circularity in 4 distinct genomic regions. Promising candidate genes were identified, e.g. those encoding an AT-hook motif nuclear-localized protein 8, type 2C protein phosphatases, and a protein Mei2-like 4 isoform, known to be associated with seed size and weight regulation. Moreover, the genes that were pinpointed are also good candidates for functional analysis to validate their influence on seed shape and size in common bean and other related crops.

Genome-wide association study reveals white lupin candidate gene involved in anthracnose resistance

GWAS identifies candidate gene controlling resistance to anthracnose disease in white lupin. White lupin (Lupinus albus L.) is a promising grain legume to meet the growing demand for plant-based protein. Its cultivation, however, is severely threatened by anthracnose disease caused by the fungal pathogen Colletotrichum lupini. To dissect the genetic architecture for anthracnose resistance, genotyping by sequencing was performed on white lupin accessions collected from the center of domestication and traditional cultivation regions. GBS resulted in 4611 high-quality single-nucleotide polymorphisms (SNPs) for 181 accessions, which were combined with resistance data observed under controlled conditions to perform a genome-wide association study (GWAS). Obtained disease phenotypes were shown to highly correlate with overall three-year disease assessments under Swiss field conditions (r > 0.8). GWAS results identified two significant SNPs associated with anthracnose resistance on gene Lalb_Chr05_g0216161 encoding a RING zinc-finger E3 ubiquitin ligase which is potentially involved in plant immunity. Population analysis showed a remarkably fast linkage disequilibrium decay, weak population structure and grouping of commercial varieties with landraces, corresponding to the slow domestication history and scarcity of modern breeding efforts in white lupin. Together with 15 highly resistant accessions identified in the resistance assay, our findings show promise for further crop improvement. This study provides the basis for marker-assisted selection, genomic prediction and studies aimed at understanding anthracnose resistance mechanisms in white lupin and contributes to improving breeding programs worldwide.

Genetic diversity and inter-gene pool introgression of Mesoamerican Diversity Panel in common beans

Brazil is among the largest producers and consumers of common bean (Phaseolus vulgaris L.) and can be considered a secondary center of diversity for the species. The aim of this study was to estimate the genetic diversity, population structure, and relationships among 288 common bean accessions in an American Diversity Panel (ADP) genotyped with 4,042 high-quality single nucleotide polymorphisms (SNPs). The results showed inter-gene pool hybridization (hybrids) between the two main gene pools (i.e., Mesoamerican and Andean), based on principal component analysis (PCA), discriminant analysis of principal components (DAPC), and STRUCTURE analysis. The genetic diversity parameters showed that the Mesoamerican group has higher values of diversity and allelic richness in comparison with the Andean group. Considering the optimal clusters (K), clustering was performed according to the type of grain (i.e., market group), the institution of origin, the period of release, and agronomic traits. A new subset was selected and named the Mesoamerican Diversity Panel (MDP), with 205 Mesoamerican accessions. Analysis of molecular variance (AMOVA) showed low genetic variance between the two panels (i.e., ADP and MDP) with the highest percentage of the limited variance among accessions in each group. The ADP showed occurrence of high genetic differentiation between populations (i.e., Mesoamerican and Andean) and introgression between gene pools in hybrids based on a set of diagnostic SNPs. The MDP showed better linkage disequilibrium (LD) decay. The availability of genetic variation from inter-gene pool hybridizations presents a potential opportunity for breeders towards the development of superior common bean cultivars.

Genome-wide association mapping reveals new loci associated with light-colored seed coat at harvest and slow darkening in carioca beans

BACKGROUND

Common bean (Phaseolus vulgaris L.) is a legume whose grain can be stored for months, a common practice among Brazilian growers. Over time, seed coats become darker and harder to cook, traits that are undesirable to consumers, who associate darker-colored beans with greater age. Like commercial pinto and cranberry bean varieties, carioca beans that have darker seeds at harvest time and after storage are subject to decreased market values.

RESULTS

The goal of our study was to identify the genetic control associated with lightness of seed coat color at harvest (HL) and with tolerance to post-harvest seed coat darkening (PHD) by a genome-wide association study. For that purpose, a carioca diversity panel previously validated for association mapping studies was used with 138 genotypes and 1,516 high-quality SNPs. The panel was evaluated in two environments using a colorimeter and the CIELAB scale. Shelf storage for 30 days had the most expressive results and the L* (luminosity) parameter led to the greatest discrimination of genotypes. Three QTL were identified for HL, two on chromosome Pv04 and one on Pv10. Regarding PHD, results showed that genetic control differs for L* after 30 days and for the ΔL* (final L*-initial L*); only ΔL* was able to properly express the PHD trait. Four phenotypic classes were proposed, and five QTL were identified through six significant SNPs.

CONCLUSIONS

Lightness of seed coat color at harvest showed an oligogenic inheritance corroborated by moderate broad-sense heritability and high genotypic correlation among the experiments. Only three QTL were significant for this trait - two were mapped on Pv04 and one on Pv10. Considering the ΔL, six QTL were mapped on four different chromosomes for PHD. The same HL QTL at the beginning of Pv10 was also associated with ΔL* and could be used as a tool in marker-assisted selection. Several candidate genes were identified and may be useful to accelerate the genetic breeding process.

Genome-Wide Association Study Reveals Genomic Regions Associated with Fusarium Wilt Resistance in Common Bean

Fusarium wilt (Fusarium oxysporum f. sp. phaseoli, Fop) is one of the main fungal soil diseases in common bean. The aim of the present study was to identify genomic regions associated with Fop resistance through genome-wide association studies (GWAS) in a Mesoamerican Diversity Panel (MDP) and to identify potential common bean sources of Fop's resistance. The MDP was genotyped with BARCBean6K_3BeadChip and evaluated for Fop resistance with two different monosporic strains using the root-dip method. Disease severity rating (DSR) and the area under the disease progress curve (AUDPC), at 21 days after inoculation (DAI), were used for GWAS performed with FarmCPU model. The p-value of each SNP was determined by resampling method and Bonferroni test. For UFV01 strain, two significant single nucleotide polymorphisms (SNPs) were mapped on the Pv05 and Pv11 for AUDPC, and the same SNP (ss715648096) on Pv11 was associated with AUDPC and DSR. Another SNP, mapped on Pv03, showed significance for DSR. Regarding IAC18001 strain, significant SNPs on Pv03, Pv04, Pv05, Pv07 and on Pv01, Pv05, and Pv10 were observed. Putative candidate genes related to nucleotide-binding sites and carboxy-terminal leucine-rich repeats were identified. The markers may be important future tools for genomic selection to Fop disease resistance in beans.

Genome-wide association analysis of bean fly resistance and agro-morphological traits in common bean

The bean fly (Ophiomyia spp) is a key insect pest causing significant crop damage and yield loss in common bean (Phaseolus vulgaris L., 2n = 2x = 22). Development and deployment of agronomic superior and bean fly resistant common bean varieties aredependent on genetic variation and the identification of genes and genomic regions controlling economic traits. This study's objective was to determine the population structure of a diverse panel of common bean genotypes and deduce associations between bean fly resistance and agronomic traits based on single nucleotide polymorphism (SNP) markers. Ninety-nine common bean genotypes were phenotyped in two seasons at two locations and genotyped with 16 565 SNP markers. The genotypes exhibited significant variation for bean fly damage severity (BDS), plant mortality rate (PMR), and pupa count (PC). Likewise, the genotypes showed significant variation for agro-morphological traits such as days to flowering (DTF), days to maturity (DTM), number of pods per plant (NPP), number of seeds per pod (NSP), and grain yield (GYD). The genotypes were delineated into two populations, which were based on the Andean and Mesoamerican gene pools. The genotypes exhibited a minimum membership coefficient of 0.60 to their respective populations. Eighty-three significant (P<0.01) markers were identified with an average linkage disequilibrium of 0.20 at 12Mb across the 11 chromosomes. Three markers were identified, each having pleiotropic effects on two traits: M100049197 (BDS and NPP), M3379537 (DTF and PC), and M13122571 (NPP and GYD). The identified markers are useful for marker-assisted selection in the breeding program to develop common bean genotypes with resistance to bean fly damage.

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

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

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

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

Genetic Diversity, Population Structure, and Andean Introgression in Brazilian Common Bean Cultivars after Half a Century of Genetic Breeding

Brazil is the largest consumer and third highest producer of common beans (Phaseolus vulgaris L.) worldwide. Since the 1980s, the commercial Carioca variety has been the most consumed in Brazil, followed by Black and Special beans. The present study evaluates genetic diversity and population structure of 185 Brazilian common bean cultivars using 2827 high-quality single-nucleotide polymorphisms (SNPs). The Andean allelic introgression in the Mesoamerican accessions was investigated, and a Carioca panel was tested using an association mapping approach. The results distinguish the Mesoamerican from the Andean accessions, with a prevalence of Mesoamerican accessions (94.6%). When considering the commercial classes, low levels of genetic differentiation were seen, and the Carioca group showed the lowest genetic diversity. However, gain in gene diversity and allelic richness was seen for the modern Carioca cultivars. A set of 1060 'diagnostic SNPs' that show alternative alleles between the pure Mesoamerican and Andean accessions were identified, which allowed the identification of Andean allelic introgression events and shows that there are putative introgression segments in regions enriched with resistance genes. Finally, genome-wide association studies revealed SNPs significantly associated with flowering time, pod maturation, and growth habit, showing that the Carioca Association Panel represents a powerful tool for crop improvements.

QTL analysis of cooking time and quality traits in dry bean (Phaseolus vulgaris L.)

Three robust QTL for dry bean cooking time shortened cooking time 11-26 min and co-localized with QTL for increased cooked seed protein concentration. Cooking time is a major factor associated with consumer preference of dry beans (Phaseolus vulgaris L.). The genetic control of cooking time was investigated with a quantitative trait loci (QTL) study on a recombinant inbred line (RIL) population developed from TZ-27 (slow cooking) and TZ-37 (fast cooking). The RIL population of 146 lines was grown on research farms over 2 years in Arusha and Morogoro, Tanzania. Arusha is an important mid-altitude bean-growing region, with moderate temperatures and reliable rainfall, whereas the low altitude and high temperatures in Morogoro make it unfavorable for bean production. The population exhibited large variation for cooking time with a range of 22-98 min. On average, beans grown in Arusha cooked 15 min faster than those grown in Morogoro. A linkage map developed with 1951 SNP markers was used for QTL analysis. Ten QTL were identified for cooking time, three of which were found in multiple environments. RILs with all three QTL (CT3.1, CT6.1, and CT11.2) cooked on average 11 min faster in Arusha and 26 min faster in Morogoro than RILs with none. Seed attributes were related to cooking time such that seeds with greater seed mass and less seed coat percentage cooked faster. Cooked seed protein concentration ranged from 17.8 to 30.8% across the years and locations. All three of the most robust cooking time QTL co-localized with QTL for protein concentration, and TZ-37 always contributed faster cooking time and increased protein concentration.

Recommendations for Choosing the Genotyping Method and Best Practices for Quality Control in Crop Genome-Wide Association Studies

High-throughput genotyping boosts genome-wide association studies (GWAS) in crop species, leading to the identification of single-nucleotide polymorphisms (SNPs) associated with economically important traits. Choosing a cost-effective genotyping method for crop GWAS requires careful examination of several aspects, namely, the purpose and the scale of the study, crop-specific genomic features, and technical and economic matters associated with each genotyping option. Once genotypic data have been obtained, quality control (QC) procedures must be applied to avoid bias and false signals in genotype–phenotype association tests. QC for human GWAS has been extensively reviewed; however, QC for crop GWAS may require different actions, depending on the GWAS population type. Here, we review most popular genotyping methods based on next-generation sequencing (NGS) and array hybridization, and report observations that should guide the investigator in the choice of the genotyping method for crop GWAS. We provide recommendations to perform QC in crop species, and deliver an overview of bioinformatics tools that can be used to accomplish all needed tasks. Overall, this work aims to provide guidelines to harmonize those procedures leading to SNP datasets ready for crop GWAS.

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

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

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

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