Next Article in Journal
Transcriptome Analyses Reveal Differential Transcriptional Profiles in Early- and Late-Dividing Porcine Somatic Cell Nuclear Transfer Embryos
Previous Article in Journal
Extracellular Vesicles-Mediated Transfer of miRNA Let-7b from PC3 Cells to Macrophages
Previous Article in Special Issue
Fine-Mapping of a Wild Genomic Region Involved in Pod and Seed Size Reduction on Chromosome A07 in Peanut (Arachis hypogaea L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 12
10.3390/genes11121496
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean

by
Ana Campa
*,
Carmen García-Fernández
and
Juan José Ferreira
Plant Genetics Group, Regional Service for Agrofood Research and Development (SERIDA), Ctra AS-267 PK 19, 33300 Asturias, Spain
*
Author to whom correspondence should be addressed.
Genes 2020, 11(12), 1496; https://doi.org/10.3390/genes11121496
Submission received: 28 October 2020 / Revised: 3 December 2020 / Accepted: 10 December 2020 / Published: 12 December 2020
(This article belongs to the Special Issue Food Legume Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
White mold (WM) is a devastating fungal disease affecting common bean (Phaseolus vulgaris L.). In this research, a genome-wide association study (GWAS) for WM resistance was conducted using 294 lines of the Spanish diversity panel. One single-locus method and six multi-locus methods were used in the GWAS. Response to this fungus showed a continuous distribution, and 28 lines were identified as potential resistance sources, including lines of Andean and Mesoamerican origin, as well as intermediate lines between the two gene pools. Twenty-two significant associations were identified, which were organized into 15 quantitative trait intervals (QTIs) located on chromosomes Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09. Seven of these QTIs were identified for the first time, whereas eight corresponded to chromosome regions previously identified in the WM resistance. In all, 468 genes were annotated in these regions, 61 of which were proposed potential candidate genes for WM resistance, based on their function related to the three main defense stages on the host: recognition (22), signal transduction (8), and defense response (31). Results obtained from this work will contribute to a better understanding of the complex quantitative resistance to WM in common bean and reveal information of significance for future breeding programs.

1. Introduction

Common bean, Phaseolus vulgaris L., is the most important grain legume for human consumption worldwide [1] and can be harvested fresh as snap beans (pods harvested before the seed development phase) or as dry beans (seeds harvested at complete maturity). It is a diploid species (2n = 2x = 22) native to the Americas, in which two major eco-geographically and genetically distinct gene pools, the Andean (AN) and the Mesoamerican (MA) have been described [2,3,4,5,6]. Common bean was introduced into Europe from the Americas in the 16th century, and since then, a considerable diversity of this crop has been traditionally grown in many Spanish regions. Part of this traditional Spanish diversity has been gathered in the Spanish Diversity Panel SDP [7].
Diseases are one of the main constraint factors in common bean production. Among them, white mold is responsible for important losses of bean yield production on a global scale. The causal agent of white mold in the common bean is the fungus Sclerotinia sclerotiorum (Lib.) de Bary, which is capable of infecting over 400 plant species worldwide [8,9]. Symptoms of the disease are easily identified with a white cottony appearance caused by the fungal mycelium on stems, leaves, flowers, or pods. Resting structures called sclerotia are produced during the infection process as hardened, dense mycelial bodies, that can survive in the soil for several years, increasing the incidence of the disease when susceptible varieties are grown. The wide host range and prolonged survival in the soil make conventional control measures against S. sclerotiorum more problematic, increasing the importance of developing resistant bean varieties. The fungus exhibits a sophisticated pathogenic process, showing an initial biotrophic lifestyle, followed by a subsequent necrotrophic lifestyle [10]. Recently, a model for the plant defense process against S. sclerotiorum infection has been proposed based on the current knowledge from different dicotyledonous hosts [11]. This plant defense model can be split into three main stages: (i) recognition of pathogen-associated molecules at the cell wall level, involving proteins for recognition, such as nucleotide-binding site-leucine-rich repeats (NBS-LRR), and different pathways, such as calcium signaling; (ii) signal-transduction at the cytoplasmic level, involving kinase cascades or signal molecules like salicylic acid, jasmonic acid, ethylene, auxin, abscisic acid, or nitric oxide, as well as reactive oxygen species; (iii) activation of specific defense responses, including transcription factors, protein kinases, pathogenesis-related proteins, reactive oxygen species production, oxidative protection, callose deposition, and production of camalexin and other specialized secondary plant metabolites [11].
Defense against this pathogen in the common bean is conditioned by both avoidance and physiological resistance mechanisms, both being quantitatively inherited [12]. Since 2001, several quantitative trait loci (QTLs) have been mapped for WM resistance and disease avoidance using different bi-parental mapping populations [12,13,14,15,16,17,18,19,20,21,22,23]. Due to the availability of an annotated reference genome sequence in common bean [24], it has been possible to define the physical position of several QTLs mapped in different populations and to establish a correspondence between them. Vasconcellos et al. [23] defined the physical position of 37 QTLs, which were condensed into 17 loci, of which nine were identified as consensus genomic regions, called meta-QTLs. However, bi-parental QTL mapping populations show some limitations, because only allelic diversity that segregates between the parents can be detected, or the limited number of recombination events decreases the mapping resolution [25].
A complementary strategy for the detection of QTLs is the use of diversity panels and genome-wide association studies (GWAS), which overcome the main limitations listed above for bi-parental mapping populations. Different methods can be used for conducting a GWAS, most being based on single-locus one-dimensional models, which do not provide true information about complex traits that are frequently controlled by multiple loci showing epistatic relationships. The single-locus model for GWAS used conventionally is the mixed linear model MLM [26,27]. Following the establishment of the MLM for single-locus analysis, several multi-locus-based methods have been proposed. The first was the multi-locus random-SNP-effect mixed linear model mrMLM [28], followed by pLARmEB [29], ISIS EM-BLASSO [30], FASTmrMLM [31], FASTmrEMMA [32], and pKWmEB [33]. To date, only one preliminary GWAS for WM resistance has been conducted in common bean, using the MLM and two diversity panels, namely the Snap Bean Diversity Panel of 150 cultivars and breeding lines, and the Snap Bean Association Panel, consisting of 382 cultivars and breeding lines [34]. This preliminary work described 39 chromosome regions involved in WM resistance, 25 of which were newly identified, with the remaining 14 corresponding to regions identified in other studies.
In the present study, a GWAS for WM resistance was conducted considering seven different GWAS methods, based on single-locus and multi-locus analyses, and using the SDP, a panel including wide phenotypic and genotypic diversity. Results will broaden the understanding of the genetic architecture of this complex quantitative resistance trait in common bean.

2. Materials and Methods

2.1. Plant Material

A total of 294 lines from the SDP were used in this study [7]. This panel was established from local Spanish germplasm, as well as including old and elite cultivars, mainly used for snap consumption. Most of the landraces included were derived from the Spanish common bean core collection [35]. According to Pascual et al. [36], lines Cornell49242 and AB136 were used as checks of susceptibility and high resistance responses, respectively.
Genotyping-by-sequencing GBS, [37] was carried out at BGI-Tech (Copenhagen, Denmark) using the ApeKI restriction enzyme. The sequencing reads obtained were aligned using the Phaseolus vulgaris L. genome v1 (www.ncbi.nlm.nih.gov; [24]) and a total of 9070 single-nucleotide polymorphism (SNPs) were mapped [7]. Data were filtered in Tassel v 5.2 software [38] for missing values (<10%) and minor allele frequencies (MAF ≥ 0.05), being considered for GWAS 5394 SNPs distributed across the eleven bean chromosomes (Figure S1). Each SNP was named to indicate the chromosome onto which it was mapped, followed by its physical location in base pairs.

2.2. White Mold Response Evaluation

Resistance tests were conducted using the straw method under controlled greenhouse conditions [39]. All lines were evaluated in two separate experiments (E1: sowing day 17 of April 2018; E2: sowing date 26 of April 2018). Each experiment consisted of two replicates (pots) per line, with 4–5 plants per pot, arranged in a randomized complete block design. Seeds were sown in 15 cm diameter plastic pots of 1.5 L volume containing 80% peats: 20% perlite and were watered and fertilized for normal plant growth. Each plant was inoculated when the third internode of the main stem was fully developed (21–28 days after planting).
The local S. sclerotiorum isolate, WM2, derived from a single hyphal tip and conserved in the form of sclerotia at −20 °C, was used for the resistance test [36]. One sclerotium was germinated on potato dextrose agar medium (Difco, Becton, Dickinson and Company, Spark, MD, USA) in a standard Petri plate and then subcultured to prepare plates of fresh and actively growing mycelium for the inoculations. A single mycelial plug was placed, mycelium-side down, on the cut main stem. The inoculated plants were maintained at moderate temperature (16 °C to 24 °C) in the shade and at a high relative humidity (70%). Disease progression was evaluated 8–10 days after inoculation (when the susceptible control was dead), based on the level of invasion of the main stem, using a 1–9 severity scale, where 1 = no symptoms and 9 = invasion of the third node and total plant collapse [40]. Values equal to or less than 4.5 were considered to be a resistant reaction (R), values between 4.5 and 7 were considered to be an intermediate reaction (I), and values equal to or greater than 7 were considered to be susceptible (S). Each plant in the pot was rated separately. The disease value per pot was calculated as the arithmetic mean of the 4–5 plants. Mean values were adjusted identifying outliers through the coefficient of variation (CV = (standard deviation/mean) × 100). Coefficient of variation over 50% was not accepted. In each evaluation, E1 and E2, the disease value per genotype was calculated as the average of the two replicates (pots). E1, E2, and the arithmetic mean of E1 and E2 (Em) were considered to be three sets of data.

2.3. Data Analysis

Statistical analyses of the genotypic and phenotypic data were conducted with the R Project for Statistical Computing [41]. The distribution of the SNPs along the eleven chromosomes was visualized with the “CMplot” package in R. Spearman’s correlation coefficient was calculated to estimate the association between results from the independent inoculations, E1 and E2.
Broad-sense heritability (H2) for the trait “WM response” within the 294 SDP lines was estimated using the function repeatability() of the “heritability” package in R. H2 was estimated at plot level in the E1 and E2 evaluations following the equation Vg/(Vg + Ve), where Vg = (MS(G) − MS(E)/r and Ve = MS(E), r is the number of replicates per genotype, and MS(G) an MS(E) are the mean sums of squares for genotype and residual error obtained from the analysis of variance. H2 was estimated at the genotype level in the Em data set following the equation Vg/(Vg + Ve/r).

2.4. Genome-Wide Association Study

To dissect the genetic architecture of a trait as complex as white mold resistance, two approaches were used in this study: (i) consistency between analysis in different data sets (E1, E2, and Em) and (ii) multiple GWAS methodologies and models.
Two different GWAS methodologies were considered: single-locus GWAS (SL-GWAS) based on MLM [26,27] and multi-locus GWAS (ML-GWAS) based on six models: mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, pKWmEB, and pLARmEB [28−33]. SL-GWAS was conducted in Tassel v 5.2 [38], while ML-GWAS was performed using the “mrMLM” package [42] in R. To account for the multiple levels of relatedness within the lines included in the panel, population structure and kinship matrix were considered. Population structure was measured by Principal Component Analysis (PCA) conducted on the genotypic data set, considering two principal components according to the results of Campa et al. [7]. The centered-IBS method was used to obtain the kinship matrix, which is an estimate of the additive genetic variance. Both, PCA and kinship matrices were obtained in Tassel v5.2 [7] and used in all the GWAS models. A threshold of −log10(p) equal to or greater than 3.0 was used to identify significant associations between a trait and an SNP.
Three data sets (E1, E2, and Em) were evaluated for each one of the seven GWAS methods (MLM, mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, pKWmEB, and pLARmEB). For each method, only significant associations identified in the mean data set (Em) and in, at least, one experimental environment (E1, E2) were considered. Each quantitative trait nucleotide (QTN) was named based on the chromosome and physical position in Mbp of the associated SNP (i.e., the SNP marker located at 42.46 Mbp of chromosome Pv01 would result in the WM1_42.46).
To verify the robustness of these QTNs, T-student tests were conducted to detect significant differences for the response to WM between the two alleles of the corresponding SNP. The observed variation was visualized in a boxplot.

2.5. Quantitative Trait Intervals

Quantitative trait intervals (QTIs) were defined as a group of QTNs located at less than 100-kbp upstream and downstream [43] from the significant SNP. In some cases, the physical positions of different QTIs overlapped or were situated very close to one another. Physical positions less than 400-kbp apart were considered to be overlapping QTIs. QTIs were named using the nomenclature chromosome and the lowest physical position in Mbp of the QTI (i.e., if the QTI expands from 42.46 to 43.02 Mbp of chromosome Pv01, it would be named as QTI1_42).
The “ShinyCircos” package [44] in R was used to visualize the position of each QTI and meta-QTLs/QTLs/QTNs previously reported in the bean genome from the underlying markers [23,34].

2.6. Candidate Gene Identification

Genes underlying each QTI were analyzed using the legume information system (https://legumeinfo.org/). Potential candidate genes were considered, based on their function and according to the three main stages proposed in the defense process against S. sclerotiorum: recognition, signal transduction, and defense response [11]. Genes coding for pathogenesis-related proteins or stress-antifungal proteins were also considered as potential candidate genes.

3. Results

3.1. White Mold Evaluation

A total of 294 common bean lines included in the SDP were evaluated against the local WM isolate in two independent evaluations, E1 and E2. E1 and E2 were significantly positively correlated (R = 0.68, p < 0.01; Figure S2). Figure 1 shows the distribution of the reaction to inoculation observed over the 1 to 9 disease severity scale. A wide and continuous distribution of disease reactions was observed. The broad-sense heritability of the trait WM response was estimated to be H2 = 0.69 in the E1, H2 = 0.78 in the E2, and H2 = 0.57 in the Em.
Most lines showed an intermediate or susceptible reaction on the disease severity scale. A total of 28 lines showed values less than 4.5 and were potential resistance sources against WM (Table 1). According to Campa et al. [7], seven of these lines came from the Andean gene pool of common bean, seven from the Mesoamerican, and 14 were recombinants between the two gene pools. Most of these potential resistance sources exhibited an indeterminate growth habit [7].

3.2. Genome-Wide Association Study

For SL-GWAS, 55 significant WM-SNP associations were identified, using the three WM data sets independently (Table S1A): 15 WM-SNP associations were identified for E1, 23 for E2, and 17 for Em. Of these associations, 13 were considered to be robust QTNs located on chromosomes Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09 (Table 2).
For ML-GWAS, 147 significant WM-SNP associations were identified in the three WM data sets (Table S1B). The mrMLM method reveals 34 significant associations, of which only three were consistent QTNs located on chromosomes Pv01, Pv02, and Pv09 (Table 2). For the FASTmrMLM method, 20 significant associations were identified, three of which were considered to be QTNs located on chromosomes Pv01, Pv03, and Pv08. For the FASTmrEMMA, 33 significant associations resulted in the identification of three QTNs on chromosomes Pv01 and Pv08. For the ISIS EM-BLASSO, 22 significant associations resulted in four QTNs on chromosomes Pv01, Pv04, and Pv08. For the pKWmEB, 27 significant associations resulted in three QTNs on chromosomes Pv04 and Pv08. Finally, for the pLARmEB 33 significant associations were identified, giving rise to eight QTNs located on Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09.
In sum, 22 robust QTNs, located on chromosomes Pv01, Pv02, Pv03, Pv04, Pv08, and Pv09, were obtained from SL-GWAS and ML-GWAS and the three data sets considered (Table 2). Among them, five QTNs were identified using more than one GWAS method. Special mention should be made for WM1_45.59 and WM8_52.39, which were identified by seven and six GWAS methods, respectively. To verify the robustness of these 22 QTNs, Student’s t-tests were carried out for the mean WM evaluation (Figure 2). Significant differences were observed in all cases.

3.3. Quantitative Trait Intervals

The 22 QTNs identified were grouped into 15 QTIs (Table 3). Three QTIs were located in a telomeric position of chromosome Pv01, two QTIs were on Pv02, two on Pv03, one on Pv04, four on Pv08, and three on Pv09. The physical positions of these QTIs were compared with those of WM resistance loci previously reported in common bean. Figure 3 shows the overlapping between chromosome regions identified in this work and the QTLs, meta-QTLs, or QTNs identified in previous studies [23,34]. In all, seven new QTIs were identified in the current study as regions involved in the response against WM: QTI1_42, QTI2_34, QTI3_0.9, QTI4_2, QTI8_52, QTI8_55, and QTI9_36. The remaining eight QTIs identified corresponded to regions previously associated with WM resistance: QTI1_45, QTI1_50, QTI2_4, QTI3_39, QTI8_13, QTI8_37, QTI9_11, and QTI9_16.

3.4. Candidate Genes

Genes underlying each QTI (spanning 5.5 Mbp) were investigated to identify potential candidate genes involved in the WM response. There were 468 genes annotated, 61 of which were identified as potential candidate genes based on their function (Table 3; Table S2). Twenty-two proteins were related to the recognition of the pathogen, eight were related to signal transduction, and 31 were related to defense response. No candidate genes were identified at QTI2_34, QTI8_13, and QTI8_37.

4. Discussion

The use of diversity panels has notably increased since the development of advanced genomic sequencing technologies and the availability of annotated reference genome sequences. In the current work, the resistance response against a local isolate of S. sclerotiorum was studied through GWAS, using the common bean SDP. The SDP was established from local Spanish germplasm, including the Spanish common bean core collection [35], and is considered to be representative of the Spanish genetic diversity for this species [7]. Using the straw test evaluation under greenhouse conditions, a moderate broad-sense heritability was observed for WM resistance. Heritability estimated for WM resistance in previous studies was also low to moderate [12,13,45,46], indicating the large environmental effect in the expression of this resistance.
A total of 28 SDP lines were identified as potential resistance sources against WM, most of them having indeterminate growth habit. Nine of these lines were obtained from the commercial snap varieties Sacha, Bilma, Helda, Planeta, Donna, Vitalis, Florencia, Marconi, and Tendergreen [7] so it could not be excluded that these resistances derived from different breeding programs in commercial varieties. The remaining 19 resistance sources were obtained from local Spanish landraces, which were assigned to different gene pools of common beans: six of them were genetically close to the Andean gene pool, six to the Mesoamerican gene pool and seven were recombinants between the two gene pools [7]. Most of the WM resistance sources identified to date are of Andean origin or from the secondary gene pool such as Phaseolus coccineus [23] so that the identification of six putative new resistance sources from the Mesoamerican gene pool is a significant finding.
Lines SDP035 and SDP071, identified as resistant in the current study, are derived from the Spanish landraces BGE003121 and BGE04000, respectively, and are both included in the Spanish common bean core collection [7,35]. This result agrees with that of Pascual et al. [36] who identified these Spanish landraces as potential sources of resistance against WM. The accession BGE003121 was identified as being of Andean origin based on 12 molecular markers, including the seed protein phaseolin [35]. Later, the use of 3099 SNP markers, obtained by mass genotyping, led to the identification of the line SDP035, derived from accession BGE003121, as a recombinant between the two gene pools, Andean and Mesoamerican [7]. The accession BGE04000 was identified as being of Andean origin, based on 12 molecular markers [35] and 3099 SNP markers [7].
In the present study, two approaches were used to improve the detection power and robustness of the GWAS: (i) analyses in multiple data sets and (ii) and multiple GWAS methodologies, based on both SL-GWAS and ML-GWAS. The control of false positives is very crucial in GWAS, but the effect of false negatives should not be ignored, which can occur if the cutoff value is too stringent. The low heritability observed for white mold resistance has led us to consider the less restrictive threshold value of −log(p) = 3. With these considerations, 22 WM QTNs were identified from the SDP. These QTNs were grouped into 15 QTIs: 3 QTIs on chromosome Pv01, 2 on Pv02, 2 on Pv03, 1 on Pv04, 4 on Pv08, and 3 on Pv09. No significant associations were identified on chromosomes Pv05, Pv06, or Pv07, on which other authors had described white mold QTLs [12,13,16,18,19,20,21,23]. It is important to note that different WM evaluation methods were used by different authors, so the results from different studies cannot always be compared. The seedling straw test under controlled conditions [39], as used in the present work, is the most frequently used method, but field evaluations [16,18,21,22] or a combination of both, controlled conditions and field, [12,13,17,19,23,40] were also used. Under field evaluations, there is a greater environmental influence and it can be difficult to distinguish between physiological resistance or avoidance mechanisms, such as growth habit, early development, the arrangement of leaves, or internode lengths.
Seven of the 15 QTIs described in the present research were identified for the first time at chromosomes Pv01 (QTI1_42), Pv02 (QTI2_34), Pv03 (QTI3_0.9), Pv04 (QTI4_2), Pv08 (QTI8_52, QTI8_55), and Pv09 (QTI9_36). The QTI4_2 was located at 2.02–2.47 Mbp on chromosome Pv04. It is not surprising that this chromosome region is involved in the resistant response against WM because it corresponds to a well-known cluster of resistance genes encoding protein kinases or proteins with NBS-LRR domains, which are typically involved in resistance response against pathogens [24]. In fact, resistance genes against multiple fungal diseases have been mapped at this position [47,48,49]. One of these genes is Co-3, which is a cluster of closely linked genes that confers resistance against Colletotrichum lindemuthianum, the causal agent of anthracnose of common bean. Co-3 has been identified mainly in the Mesoamerican gene pool. This result could explain why Pv04 is identified in the WM response only in GWAS, where the panels include materials from the Mesoamerican gene pool, whereas most of the resistant sources used in bi-parental populations are of Andean origin [23].
The remaining eight QTIs identified in this work co-localized to WM QTLs, meta-QTLs or QTNs previously reported [23,34]:
The QTI1_45 (45.36–45.69 Mbp) corresponds to several QTLs identified from different populations (WM1.1XC, WM1.1AO, WM1.1O83) which together extends from 45.15 to 50.22 Mbp of chromosome Pv01 [23]. A chromosome region comprising 43.71–49.94 Mbp of Pv01 was also significantly associated with WM response in the GWAS [34]. QTI1_45 could be related to morphological plant traits involved in pathogen avoidance response, because at this position (45.56 Mbp) the gene PHAVU_001G189200g has been proposed as a candidate of fin, the gene controlling growth habit with recessive genotypes showing determinacy [50,51]. The implication of this chromosome region in avoidance mechanisms have also been suggested by other authors [12,20].
The QTI1_50 (50.97–51.17 Mbp of Pv01) is located close to the meta-QTL WM1.1 (49.57–50.22 Mbp) define by [23] and could correspond to the same region.
Correspondence was observed between QTI2_4 (4.53–5.68 Mbp), the meta-QTL WM2.2 (3.57–5.24 Mbp), and QTLs WM2.2AN and WM2.2R31, which together extends from 3.91 to 21.76 Mbp of Pv02 [23]. A chromosome region between 3.14 and 4.61 Mbp of Pv02 was significantly associated with WM response in the GWAS [34].
The QTI3_39 (39.85–40.64 Mbp) co-localized with the meta-QTL WM3.1 (34.96–39.49 Mbp), with two QTLs WM3.1AP, WM3.1XC (34.33–48.32 Mbp) [23], and the region between 33.55–44.48 Mbp of Pv03 was also identified in the GWAS [34].
QTI8_13 (13.22–13.42 Mbp) and QTI8_37 (37.04–37.24 Mbp) overlap with QTLs WM8.3PS02−029C and WM8.3Z0725, which extend over most of chromosome Pv08, from 4.97 to 46.71 Mbp [23]. These regions (11.16–14.70 Mbp and 32.48–37.24 Mbp of Pv08) had also been identified in the GWAS [34]. Moreover, QTI8_37 is located less than 100 kbp from the meta-QTL WM8.3 (37.32–46.73 Mbp) so they could be considered the same region.
QTI9_11 (11.61–11.88 Mbp) and QTI9_16 (16.68–16.88 Mbp) extended from 8.19 to 17.73 Mbp, a region significantly associated with WM response in previous GWAS [34].
Among the genes underlying these QTIs, 61 were proposed as potential candidate genes, based on their function related to the three main stages proposed against the defense process: recognition, signal transduction, and defense response [11]. Among them, 15 Cytochrome P450 proteins were identified organized into two main clusters at QTI3_0.9 and QTI4_2, and 15 genes encoding typical R products like protein kinases or NBS-LRR proteins at QTI1_42, QTI1_45, QTI1_50, QTI2_4, QTI3_39, QTI4_2, QTI8_52, QTI9_11, and QTI9_16. Cytochromes P450s are a superfamily of enzymes involved in several plant functions, including response to biotic stresses [52]. For instance, transgenic plants of Nicotiana benthamiana overexpressing the gene GmCYP82A3 of the cytochrome P450 family CYP82 exhibited high resistance to gray mold (Botrytis cinerea) and black shank (Phytophthora nicotianae) [53].
Protein kinase and NBS-LRR genes constitute the largest plant disease resistance gene (R gene) family. Protein kinases are involved in defense mechanisms at different cellular levels, including elicitor recognition, as extracellular or intracellular receptors, in signal transduction, and the induction of transcriptional activation [54]. Furthermore, many NBS-LRR proteins recognize effectors secreted by pathogens directly or indirectly, which, in turn, activate downstream signaling pathways, leading to the activation of plant defense responses against various classes of pathogens, including bacterial, fungal, and viral, as well as nematode and insect pests [47,55]. Among the potential candidate genes identified, 28 encode for typical R genes products like protein kinases or NBS-LRR proteins and are located in all QTIs except for QTI8_55 and QTI9_36.
Genes Phvul.002G055700 and Phvul.002G055800, both ethylene-responsive transcription factors, have been proposed as candidate resistance genes against WM by Vasconcellos et al. [23]. Ethylene signaling has been shown to be involved in defense against S. sclerotiorum in Arabidopsis thaliana [56] and the closely related Brassica napus [57]. Both genes, Phvul.002G055700 and Phvul.002G055800, were identified in the present research at QTI2_4 and have been proposed as potential candidate genes against WM [23].
White mold produces many molecules and proteins, such as polygalacturonases (PGs), during infection and colonization to breach host cell walls [58]. Host plants produce polygalacturonase-inhibiting proteins (PGIPs) that recognize PGs and prevent enzymatic action during the invasion. In common bean, PGIPs consist of a complex locus of four genes organized in a cluster (PvPGIP1, PvPGIP2, PvPGIP3, and PvPGIP4) that spans a 50-kbp region at around 36 Mbp of chromosome Pv02 [59]. It has been demonstrated that this locus in the common bean is involved in defense against insects and fungi [59], including S. sclerotiorum infection [60,61].
Close to this PGIP cluster was located the QTI2_34 (34.54–34.84 Mbp), in which no candidate genes for WM resistance had previously been identified. It cannot be discarded that this cluster of PGIP genes may be the causal genes for WM response at this QTI. It is well known that SNPs may occasionally affect distant genes and the common practice of identifying candidate genes based on the nearest SNP association may lead to some difficulties in the identification of candidate genes [62,63].

5. Conclusions

The evaluation of the response to a local isolate of S. sclerotiorum in the SDP revealed 28 lines with high levels of resistance against this pathogen. Among these potential resistance sources, six were of Mesoamerican origin. This is a significant finding because most of the WM resistance sources identified to date are of Andean origin or from the secondary gene pool. The GWAS led to the validation of eight regions previously reported as being associated with WM resistance at chromosomes Pv01 (QTI1_45, QTI1_50), Pv02 (QTI2_4), Pv03 (QTI3_39), Pv08 (QTI8_13 and QTI8_37), and Pv09 (QTI9_11, QTI9_16). The GWAS also identified seven novel regions as being involved in the WM response located at chromosomes Pv01 (QTI1_42), Pv02 (QTI2_34), Pv03 (QTI3_0.9), Pv04 (QTI4_2), Pv08 (QTI8_52, QTI8_55), and Pv09 (QTI9_36). These results verify the complexity of this quantitative resistance, validate the implication of eight chromosome regions in WM response, and provide new targets that could be of major interest to bean breeders.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/12/1496/s1, Figure S1: Distribution along the eleven bean chromosomes of the 5384 single nucleotide polymorphisms (SNPs); Figure S2: Scatter plot obtained using the ggpubr package in R to display the relationship between the two variables WM evaluation 1 and WM evaluation 2. R, Spearman correlation coefficient; p, significance level; Table S1: (A) Significant SNP-WM response associations (−log10(p) ≥ 3) identified in the single-locus analysis conducted in Tassel. (B) Significant SNP-WM response associations (−log10(p) ≥ 3) identified in the multi-locus analysis conducted in mrMLM package of R software; Table S2: Genes underlying each QTI identified in this work. Potential candidate genes are indicated in red color. The potential implication of these candidate genes in recognition, signal transduction, or defense response is indicated.

Author Contributions

Data curation, A.C.; formal analysis, A.C. and C.G.-F.; funding acquisition, A.C. and J.J.F.; project administration, J.J.F.; supervision, J.J.F.; writing—original draft, A.C.; writing—review and editing, A.C., C.G.-F., and J.J.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Government, grant number AGL2017-87050-R.

Acknowledgments

The authors thank S. Panera, M. Bueno, J.A. Poladura, and F. Díaz for their technical assistance.

Conflicts of Interest

The authors declare that there are no competing interests.

References

  1. Food and Agriculture Organization of the United Nations. FAOSTAT Statistics Database; FAO: Rome, Italy, 1998. [Google Scholar]
  2. Gepts, P.; Osborne, T.C.; Rashka, K.; Bliss, F.A. Phaseolin protein variability in wild forms and landraces of the common bean (Phaseolus vulgaris L.): Evidence for multiple centers of domestication. Econ. Bot. 1986, 40, 451–468. [Google Scholar] [CrossRef]
  3. Koenig, R.L.; Gepts, P. Allozyme diversity in wild Phaseolus vulgaris: Further evidence for two major centers of genetic diversity. Theor. Appl. Genet. 1989, 78, 809–817. [Google Scholar] [CrossRef]
  4. Singh, S.P.; Gepts, P.; Debouck, D.G. Races of common bean (Phaseolus vulgaris, Fabaceae). Econ. Bot. 1991, 45, 379–396. [Google Scholar] [CrossRef]
  5. Kwak, M.; Gepts, P. Structure of genetic diversity in the two major gene pools of common bean (Phaseolus vulgaris L., Fabaceae). Theor. Appl. Genet. 2009, 118, 979–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rossi, M.; Bitocchi, E.; Bellucci, E.; Nanni, L.; Rau, D.; Attene, G.; Papa, R. Linkage disequilibrium and population structure in wild and domesticated populations of Phaseolus vulgaris L. Evol. Appl. 2009, 2, 504–522. [Google Scholar] [CrossRef] [PubMed]
  7. Campa, A.; Murube, E.; Ferreira, J.J. Genetic diversity, population structure, and linkage disequilibrium in a Spanish common bean diversity panel revealed through genotyping-by-sequencing. Genes 2018, 9, 518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bolland, G.J.; Hall, R. Index of plant host of Sclerotinia sclerotiorum. Can. J. Plant Pathol. 1994, 16, 93–108. [Google Scholar] [CrossRef]
  9. Schwartz, H.F.; Steadman, J.R.; Hall, R.; Forster, R.L. Compendium of Bean Diseases, 2nd ed.; APS Press: Fort Collins, CO, USA, 2005; p. 109. ISBN 0890543283. [Google Scholar]
  10. Kabbage, M.; Yarden, O.; Dickman, M.B. Pathogenic attributes of Sclerotinia sclerotiorum: Switching from a biotrophic to necrotrophic lifestyle. Plant Sci. 2015, 233, 53–60. [Google Scholar] [CrossRef]
  11. Wang, Z.; Ma, L.-Y.; Cao, J.; Li, Y.-L.; Ding, L.-N.; Zhu, K.-M.; Yang, Y.-H.; Tan, X.-L. Recent advances in mechanism of plant defense to Sclerotinia sclerotiorum. Front. Plant Sci. 2019, 10, 1314. [Google Scholar] [CrossRef] [PubMed]
  12. Miklas, P.N.; Delorme, R.; Johnson, W.C.; Gepts, P. QTL conditioning physiological resistance and avoidance to white mold in dry bean. Crop Sci. 2001, 41, 309–315. [Google Scholar] [CrossRef]
  13. Park, O.S.; Coyne, D.P.; Steadman, J.R.; Skroch, P.W. Mapping of QTL for resistance to white mold disease in common bean. Crop Sci. 2001, 41, 1253–1262. [Google Scholar] [CrossRef]
  14. Kolkman, J.M.; Kelly, J.D. QTL conferring resistance and avoidance to white mold in common bean. Crop Sci. 2003, 43, 539–548. [Google Scholar] [CrossRef]
  15. Miklas, P.N.; Delorme, R. Identification of QTL conditioning resistance to white mold in snap bean. J. Am. Soc. Hortic. Sci. 2003, 128, 564–570. [Google Scholar] [CrossRef] [Green Version]
  16. Ender, M.; Kelly, J.D. Identification of QTL associated with white mold resistance in common bean. Crop Sci. 2005, 45, 2482–2490. [Google Scholar] [CrossRef]
  17. Maxwell, J.J.; Brick, M.A.; Byrne, P.F.; Schwartz, H.F.; Shan, X.; Ogg, J.B.; Hensen, R.A. Quantitative trait loci linked to white mold resistance in common bean. Crop Sci. 2007, 47, 2285–2294. [Google Scholar] [CrossRef]
  18. Miklas, P.N.; Larsen, K.M.; Terpstra, K.A.; Hauf, D.C.; Grafton, K.F.; Kelly, J.D. QTL analysis of ICA Bunsi-derived resistance to white mold in a Pinto 3 navy bean cross. Crop Sci. 2007, 47, 174–179. [Google Scholar] [CrossRef] [Green Version]
  19. Soule, M.; Porter, L.; Medina, J.; Santana, G.P.; Blair, M.W.; Miklas, P.N. Comparative QTL map for white mold resistance in common bean, and characterization of partial resistance in dry bean lines VA19 and I9365-31. Crop Sci. 2011, 51, 123–139. [Google Scholar] [CrossRef] [Green Version]
  20. Pérez-Vega, E.; Pascual, A.; Campa, A.; Giraldez, R.; Miklas, P.N.; Ferreira, J.J. Mapping quantitative trait loci conferring partial physiological resistance to white mold in the common bean RIL population Xana x Cornell49242. Mol. Breed. 2012, 29, 31–41. [Google Scholar] [CrossRef]
  21. Miklas, P.N.; Porter, L.D.; Kelly, J.D.; Myers, J.R. Characterization of white mold disease avoidance in common bean. Eur. J. Plant Pathol. 2013, 135, 525–543. [Google Scholar] [CrossRef] [Green Version]
  22. Hoyos-Villegas, V.; Mkwaila, W.; Cregan, P.B.; Kelly, J.D. Quantitative trait loci analysis of white mold avoidance in pinto bean. Crop Sci. 2015, 55, 2116–2129. [Google Scholar] [CrossRef]
  23. Vasconcellos, R.C.C.; Oraguzie, O.B.; Soler, A.; Arkwazee, H.; Myers, J.R.; Ferreira, J.J.; Song, Q.; McClean, P.; Miklas, P.N. Meta-QTL for resistance to white mold in common bean. PLoS ONE 2017, 12, e0171685. [Google Scholar] [CrossRef] [PubMed]
  24. Schmutz, J.; McClean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Jenkins, J.; Shu, S.; Song, Q.; Chavarro, C.; et al. A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Korte, A.; Farlow, A. The advantages and limitations of trait analysis with GWAS: A review. Plant Methods 2013, 9, 29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zhang, Y.M.; Mao, Y.; Xie, C.; Smith, H.; Luo, L.; Xu, S. Mapping quantitative trait loci using naturally occurring genetic variance among commercial inbred lines of maize (Zea mays L.). Genetics 2005, 169, 2267–2275. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Yu, J.; Pressoir, G.; Briggs, W.H.; Bi, I.V.; Yamasaki, M.; Doebley, J.F.; McMullen, M.D.; Gaut, B.S.; Nielsen, D.M.; Holland, J.B.; et al. A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 2006, 38, 203–208. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, S.B.; Feng, J.Y.; Ren, W.L.; Huang, B.; Zhou, L.; Wen, Y.J.; Zhang, J.; Dunwell, J.M.; Xu, S.; Zhang, Y.M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 2016, 6, 19444. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, J.; Feng, J.Y.; Ni, Y.L.; Wen, Y.J.; Niu, Y.; Tamba, C.L.; Yue, C.; Song, Q.; Zhang, Y.M. pLARmEB: Integration of least angle regression with empirical bayes for multi-locus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef] [Green Version]
  30. Tamba, C.L.; Ni, Y.L.; Zhang, Y.M. Iterative sure independence screening EM-Bayesian LASSO algorithm for multilocus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef]
  31. Zhang, Y.M.; Tamba, C.L. A fast mrMLM algorithm for multi-locus genome-wide association studies. bioRxiv 2018. [Google Scholar] [CrossRef]
  32. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wang, S.B.; Dunwell, J.M.; Zhang, Y.M.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief. Bioinform. 2018, 19, 700–712. [Google Scholar] [CrossRef] [Green Version]
  33. Ren, W.L.; Wen, Y.J.; Dunwell, J.M.; Zhang, Y.M. pKWmEB: Integration of Kruskal-Wallis test with empirical bayes under polygenic background control for multi-locus genome-wide association study. Heredity 2018, 120, 208–218. [Google Scholar] [CrossRef] [PubMed]
  34. Arkwazee, H.; Hart, J.; Porch, T.; Griffiths, P.; Davis, J.; Myers, J.R. Genome wide association study (GWAS) for white mold resistance in snap bean. Annu. Rept. Bean Improv. Coop. 2018, 61, 85. [Google Scholar]
  35. Pérez-Vega, E.; Campa, A.; De la Rosa, L.; Giraldez, R.; Ferreira, J.J. Genetic diversity in a core collection established from the main genebank in Spain. Crop Sci. 2009, 49, 1377–1386. [Google Scholar] [CrossRef]
  36. Pascual, A.; Campa, A.; Pérez-Vega, E.; Giraldez, R.; Miklas, P.N.; Ferreira, J.J. Screening common bean for resistance to four Sclerotinia sclerotiorum isolates collected in northern Spain. Plant Dis. 2010, 94, 885–890. [Google Scholar] [CrossRef] [Green Version]
  37. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef] [Green Version]
  38. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  39. Petzoldt, R.; Dickson, M.H. Straw test for resistance to white mold in beans. Annu. Rep. Bean Improv. Coop. Meet 1996, 39, 142–143. [Google Scholar]
  40. Miklas, P.N. Marker-assisted backcrossing QTL for partial resistance to Sclerotinia white mold in dry bean. Crop Sci. 2007, 47, 935–942. [Google Scholar] [CrossRef] [Green Version]
  41. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; ISBN 3-900051-07-0. Available online: http://www.R-project.org (accessed on 1 April 2020).
  42. Zhang, Y.W.; Tamba, C.L.; Wen, Y.J.; Li, P.; Ren, W.L.; Ni, Y.L.; Gao, J.; Zhang, Y.M. mrMLM v4.0: An R platform for multi-locus genome-wide association studies. bioRxiv 2020. [Google Scholar] [CrossRef]
  43. Moghaddam, S.M.; Mamidi, S.; Osorno, J.M.; Lee, R.; Brick, M.; Kelly, J.; Miklas, P.; Urrea, C.; Song, Q.; Cregan, P.; et al. Genome-Wide association study identifies candidate loci underlying agronomic traits in a middle american diversity panel of common bean. Plant Genome 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  44. Yu, Y.; Ouyang, Y.; Yao, W. ShinyCircos: An R/shiny application for interactive creation of circos plot. Bioinformatics 2018, 34, 1229–1231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Fuller, P.A.; Coyne, D.P.; Steadman, J.R. Inheritance of resistance to white mold disease in a diallel cross of dry beans. Crop Sci. 1984, 24, 929–933. [Google Scholar] [CrossRef]
  46. Miklas, P.N.; Hauf, D.C.; Henson, R.A.; Grafton, K.F. Inheritance of ICA Bunsi derived resistance in a navy × pinto bean cross. Crop Sci. 2004, 44, 1584–1588. [Google Scholar] [CrossRef] [Green Version]
  47. Meziadi, C.; Richard, M.M.S.; Derquennes, A.; Thareau, V.; Blanchet, S.; Gratias, A.; Pflieger, S.; Geffroy, V. Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence. Plant Sci. 2016, 242, 351–357. [Google Scholar] [CrossRef] [PubMed]
  48. Campa, A.; Ferreira, J.J. Gene coding for an elongation factor is involved in resistant against powdery mildew in common bean. Theor. Appl. Genet. 2017, 130, 849–860. [Google Scholar] [CrossRef] [PubMed]
  49. Murube, E.; Campa, A.; Ferreira, J.J. Integrating genetic and physical positions of the anthracnose resistance genes described in bean chromosomes Pv01 and Pv04. PLoS ONE 2019, 14, e0212298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Kwak, M.; Velasco, D.; Gepts, P. Mapping homologous sequences for determinacy and photoperiod sensitivity in common bean (Phaseolus vulgaris). J. Hered. 2008, 99, 283–291. [Google Scholar] [CrossRef]
  51. Repinski, S.L.; Kwak, M.; Gepts, P. The common bean growth habit gene PvTFL1y is a functional homolog of Arabidopsis TFL1. Theor. Appl. Genet. 2012, 124, 1539–1547. [Google Scholar] [CrossRef]
  52. Pandian, B.A.; Sathishraj, R.; Djanaguiraman, M.; Prasad, P.V.V.; Jugulam, M. Role of cytochrome P450 enzymes in plant stress response. Antioxidants 2020, 9, 454. [Google Scholar] [CrossRef]
  53. Yan, Q.; Cui, X.; Lin, S.; Gan, S.; Xing, H.; Dou, D. GmCYP82A3, a soybean cytochrome P450 family gene involved in the jasmonic acid and ethylene signaling pathway, enhances plant resistance to biotic and abiotic stresses. PLoS ONE 2016, 11, e0162253. [Google Scholar] [CrossRef] [Green Version]
  54. Sessa, G.; Martin, G.B. Protein kinases in the plant defence response. Adv. Bot. Res. 2000, 32, 379–404. [Google Scholar] [CrossRef]
  55. Marone, D.; Russo, M.A.; Laido, G.; Leonardis, A.M.; Mastrangelo, A.M. Plant nucleotide binding site-leucine-rich repeat (NBS-LRR) genes: Active guardians in host defense response. Int J. Mol. Sci. 2013, 14, 7302–7326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Guo, X.; Stotz, H.U. Defense against Sclerotinia sclerotiorum in arabidopsis is dependent on jasmonic acid, salicylic acid, and ethylene signaling xiaomei. Mol. Plant Microbe. Interact. 2007, 20, 1384–1395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Zhao, J.; Wang, J.; An, L.; Doerge, R.W.; Chen, Z.J.; Grau, C.R.; Meng, J.; Osborn, T.C. Analysis of gene expression profiles in response to Sclerotinia sclerotiorum in Brassica napus. Planta 2007, 227, 13–24. [Google Scholar] [CrossRef] [PubMed]
  58. Kalunke, R.M.; Janni, M.; Sella, L.; David, P.; Geffroy, V.; Favaron, F.; D´Ovidio, R. Transcript analysis of the bean polygalacturonase-inhibiting protein gene family reveals that Pvpgip2 is expressed in the whole plant and is strongly induced by pathogen infection. J. Plant Pathol. 2011, 93, 141–148. [Google Scholar] [CrossRef]
  59. D´Ovidio, R.; Raiola, A.; Capodicasa, C.; Devoto, A.; Pontiggia, D.; Roberti, S.; Galletti, R.; Conti, E.; O’Sullivan, D.; De Lorenzo, G. Characterization of the complex locus of bean encoding Polygalacturonase-Inhibiting Proteins reveals subfunctionalization for defense against fungi and insects. Plant Physiol. 2004, 135, 2424–2435. [Google Scholar] [CrossRef] [Green Version]
  60. Vasconcellos, R.C.C.; Lima, T.F.C.; Fernandes-Brum, C.N.; Chalfun-Junior, A.; Santos, J.B. Expression and validation of PvPGIP genes for resistance to white mold (Sclerotinia sclerotiorum) in common beans (Phaseolus vulgaris L.). Genet. Mol. Res. 2016, 15. [Google Scholar] [CrossRef]
  61. Porto, A.C.M.; Cardon, C.H.; Vasconcellos, R.C.C.; Novaes, E.; Leite, M.E.; Chalfun-Junior, A.; Pereira, W.A.; Santos, J.B. Expression of candidate genes related to white mold resistance in common beans. Trop. Plant Pathol. 2019, 44, 483–493. [Google Scholar] [CrossRef]
  62. Visel, A.; Rubin, E.M.; Pennacchio, L.A. Genomic views of distant-acting enhancers. Nature 2009, 461, 199–205. [Google Scholar] [CrossRef] [Green Version]
  63. Brodie, A.; Azaria, J.R.; Ofran, Y. How far from the SNP may the causative genes be? Nucleic Acids Res. 2016, 44, 6046–6054. [Google Scholar] [CrossRef]
Genes 11 01496 g001 550
Figure 1. Histogram showing the distribution of the reaction against a local Sclerotinia sclerotiorum isolate in the 294 lines of the SDP (Spanish Diversity Panel) and the (A) E1, (B) E2, (C) Em. Arrows indicate the average value in the disease severity scale of resistant and susceptible checks.
Figure 1. Histogram showing the distribution of the reaction against a local Sclerotinia sclerotiorum isolate in the 294 lines of the SDP (Spanish Diversity Panel) and the (A) E1, (B) E2, (C) Em. Arrows indicate the average value in the disease severity scale of resistant and susceptible checks.
Genes 11 01496 g001
Genes 11 01496 g002 550
Figure 2. Box plot for the mean WM evaluation plotted according to the haplotype for the 22 SNPs associated with WM response. The X-axis represents the two alleles for each SNP, while the Y-axis corresponds to the phenotype (1 to 9 disease severity scale). T, T-student test value; ** p < 0.01; *** p < 0.001.
Figure 2. Box plot for the mean WM evaluation plotted according to the haplotype for the 22 SNPs associated with WM response. The X-axis represents the two alleles for each SNP, while the Y-axis corresponds to the phenotype (1 to 9 disease severity scale). T, T-student test value; ** p < 0.01; *** p < 0.001.
Genes 11 01496 g002
Genes 11 01496 g003 550
Figure 3. Circos plot showing the overlapping between WM quantitative trait intervals (QTIs) identified in this work (red track) and WM meta-quantitative trait loci (QTLs) (dark blue track), QTLs (orange track), and quantitative trait nucleotides (QTNs) (green track) previously described.
Figure 3. Circos plot showing the overlapping between WM quantitative trait intervals (QTIs) identified in this work (red track) and WM meta-quantitative trait loci (QTLs) (dark blue track), QTLs (orange track), and quantitative trait nucleotides (QTNs) (green track) previously described.
Genes 11 01496 g003
Table
Table 1. Lines included in the SDP (Spanish Diversity Panel) showing partial resistance against the Sclerotinia sclerotiorum local isolate. Resistant and susceptible checks are included. The severity scale value assigned in each evaluation (E1, E2, and Em) is indicated. According to Campa et al. [7], the accession/variety from which each SDP line was obtained, the gene pool (GP) assigned to each line (A, Andean; MA, Mesoamerican; M, recombinant between A and MA) and the growth habit (GH: det, determinate vs. ind, indeterminate) are indicated.
Table 1. Lines included in the SDP (Spanish Diversity Panel) showing partial resistance against the Sclerotinia sclerotiorum local isolate. Resistant and susceptible checks are included. The severity scale value assigned in each evaluation (E1, E2, and Em) is indicated. According to Campa et al. [7], the accession/variety from which each SDP line was obtained, the gene pool (GP) assigned to each line (A, Andean; MA, Mesoamerican; M, recombinant between A and MA) and the growth habit (GH: det, determinate vs. ind, indeterminate) are indicated.
SDPAccession/VarietyType of MaterialGPGHE1E2Em
SDP194BGE043359Spanish landraceAInd4.04.04.0
SDP217CN_227Spanish landraceMInd4.04.04.0
SDP289SachaCommercial varietyMAInd4.24.04.1
SDP203BILMACommercial varietyMAInd4.34.04.1
SDP216CN_226Spanish landraceMInd4.34.04.1
SDP246HeldaCommercial varietyMAInd4.14.14.1
SDP274PLANETACommercial varietyMAInd4.14.34.2
SDP226DONNACommercial varietyMAInd4.34.24.2
SDP108BGE011762Spanish landraceAInd4.34.34.3
SDP039BGE003139Spanish landraceMAInd4.24.44.3
SDP133BGE022831Spanish landraceMAInd4.74.04.3
SDP305VITALISCommercial varietyMAInd4.54.14.3
SDP234FLORENCIACommercial varietyMAInd4.44.34.4
SDP304V381Spanish landraceMAInd4.24.54.4
SDP106BGE011736Spanish landraceMDet4.84.04.4
SDP107BGE011758Spanish landraceMAInd4.63.74.4
SDP255MARCONICommercial varietyMAInd4.64.24.4
SDP189BGE039982Spanish landraceMInd4.64.14.4
SDP025BGE002196Spanish landraceMInd4.44.54.5
SDP035BGE003121Spanish landraceMInd4.54.54.5
SDP053BGE003482Spanish landraceAInd4.44.64.5
SDP071BGE004000Spanish landraceAInd5.04.04.5
SDP088BGE005475Spanish landraceMAInd5.23.84.5
SDP221CN_241Spanish landraceADet4.64.44.5
SDP295TendergreenCommercial varietyADet5.04.14.5
SDP110BGE013962Spanish landraceMInd4.94.04.5
SDP150BGE025745Spanish landraceAInd5.04.14.5
SDP024BGE002189Spanish landraceMAInd5.04.04.5
AB136Resistant checkMAInd4.04.04.0
Cornell49242Susceptible checkMAInd8.99.09.0
Table
Table 2. Significant white mold (WM)-SNP associations identified using one single locus-genome-wide association study (SL-GWAS) method (mixed linear model—MLM), 6 ML-GWAS methods (mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, pLARmEB, and pKWmEB), and three data sets (E1, E2, and Em). SNPs are named according to chromosome number and physical position in bp.
Table 2. Significant white mold (WM)-SNP associations identified using one single locus-genome-wide association study (SL-GWAS) method (mixed linear model—MLM), 6 ML-GWAS methods (mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, pLARmEB, and pKWmEB), and three data sets (E1, E2, and Em). SNPs are named according to chromosome number and physical position in bp.
E1E2Em
QTNSNPGWAS Method−log10(p)−log10(p)−log10(p)
WM1_42.55s1_42558186MLM3-3
WM1_42.85s1_42857105MLM4-3
pLARmEB7-6
WM1_42.91s1_42916300FASTmrEMMA7-5
WM1_45.45s1_45457672MLM-43
WM1_45.59s1_45594842MLM-64
mrMLM-75
FASTmrMLM-910
ISIS EM-BLASSO-1111
FASTmrEMMA-75
pLARmEB-68
pLARmEB5-7
WM1_51.06s1_51067156ISIS EM-BLASSO3-3
WM2_4.62s2_4625990MLM9-5
WM2_5.57s2_5576016pLARmEB3-3
WM2_34.64s2_34644216MLM-54
WM2_34.73s2_34738605mrMLM3-3
WM3_0.9s3_947347FASTmrMLM-45
WM3_39.95s3_39951398MLM7-7
pLARmEB3-3
WM3_40.53s3_40535470MLM-65
WM4_2.11s4_2115855ISIS EM-BLASSO-43
WM4_2.37s4_2371248MLM-76
pKWmEB-48
pLARmEB435
WM8_13.31s8_13317583MLM-33
WM8_37.13s8_37139214MLM4-4
WM8_52.39s8_52395252MLM7-7
FASTmrMLM667
ISIS EM-BLASSO10-10
FASTmrEMMA-55
pKWmEB5-4
pLARmEB-76
WM8_55.56s8_55569419pKWmEB3-3
WM9_11.70s9_11707260MLM-54
WM9_16.78s9_16783907mrMLM-88
WM9_36.67s9_36677215pLARmEB3-3
Table
Table 3. QTIs identified in the SDP for resistance to WM. The physical position of each QTI is indicated (chromosome, start, and end position in Mbp). The number of annotated genes and candidate genes underlying each QTI is indicated. Potential candidate genes were grouped in the main three stages described for WM response: recognition, signaling, and defense.
Table 3. QTIs identified in the SDP for resistance to WM. The physical position of each QTI is indicated (chromosome, start, and end position in Mbp). The number of annotated genes and candidate genes underlying each QTI is indicated. Potential candidate genes were grouped in the main three stages described for WM response: recognition, signaling, and defense.
QTIChrStart-EndQTNsNº Annotated GenesNº Candidate GenesStage of Resistance Response (Number of Genes)
QTI1_42Pv0142.46−43.02WM1_42.55, WM11_42.85, WM1_42.91468Recognition (6), Signaling (2)
QTI1_45Pv0145.36−45.69WM1_45.45, WM1_45.59364Recognition (1), Signaling (1), Defense (2)
QTI1_50Pv0150.97–51.17WM1_51.06282Recognition (1), Defense (1)
QTI2_4Pv024.53–5.68WM2_4.62, WM2_5.57747Recognition (3), Signaling (2), Defense (2)
QTI2_34Pv0234.54–34.84WM2_34.64, WM2_34.7323--
QTI3_0.9Pv030.80–1.05WM3_0.9296Recognition (1), Defense (5)
QTI3_39Pv0339.85–40.64WM3_39.95, WM3_40.53777Recognition (2), Signaling (1), Defense (4)
QTI4_2Pv042.02–2.47WM4_2.11, WM4_2.374212Recognition (1), Defense (11)
QTI8_13Pv0813.22–13.42WM8_13.316--
QTI8_37Pv0837.04–37.24WM8_37.132--
QTI8_52Pv0852.30–52.50WM8_52.39245Recognition (4), Defense (1)
QTI8_55Pv0855.47–55.67WM8_55.56221Signaling (1)
QTI9_11Pv0911.61–11.88WM9_11.70242Recognition (2)
QTI9_16Pv0916.68–16.88WM9_16.78153Recognition (1), Defense (2)
QTI9_36Pv0936.58–36.78WM9_36.67204Signaling (1), Defense (3)
46861
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Campa, A.; García-Fernández, C.; Ferreira, J.J. Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean. Genes 2020, 11, 1496. https://doi.org/10.3390/genes11121496

AMA Style

Campa A, García-Fernández C, Ferreira JJ. Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean. Genes. 2020; 11(12):1496. https://doi.org/10.3390/genes11121496

Chicago/Turabian Style

Campa, Ana, Carmen García-Fernández, and Juan José Ferreira. 2020. "Genome-Wide Association Study (GWAS) for Resistance to Sclerotinia sclerotiorum in Common Bean" Genes 11, no. 12: 1496. https://doi.org/10.3390/genes11121496

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop