Ana Campa
Valérie Geffroy2,3
Elena Bitocchi
Juan José Ferreira- 1Plant Genetic Group, Regional Service for Agrofood Research and Development (SERIDA), Villaviciosa, Asturias, Spain
- 2Université Paris-Saclay, CNRS, INRAE, Université Evry, Institute of Plant Sciences Paris-Saclay (IPS2), Gif sur Yvette, France
- 3Université Paris-Cité, CNRS, INRAE, Institute of Plant Sciences Paris-Saclay (IPS2), Gif sur Yvette, France
- 4Department of Agricultural, Food and Environmental Sciences, Marche Polytechnic University, Via Brecce Bianche, Ancona, Italy
Anthracnose, white mold, powdery mildew, and root rot caused by Colletotrichum lindemuthianum, Scletorinia sclerotiorum, Erysiphe spp., and Pythium ultimum, respectively, are among the most frequent diseases that cause significant production losses worldwide in common bean (Phaseolus vulgaris L.). Reactions against these four fungal diseases were investigated under controlled conditions using a diversity panel of 311 bean lines for snap consumption (Snap bean Panel). The genomic regions involved in these resistance responses were identified based on a genome-wide association study conducted with 16,242 SNP markers. The highest number of resistant lines was observed against the three C. lindemuthianum isolates evaluated: 156 lines were resistant to CL124 isolate, 146 lines resistant to CL18, and 109 lines were resistant to C531 isolate. Two well-known anthracnose resistance clusters were identified, the Co-2 on chromosome Pv11 for isolates CL124 and CL18, and the Co-3 on chromosome Pv04 for isolates CL124 and C531. In addition, other lesser-known regions of anthracnose resistance were identified on chromosomes Pv02, Pv06, Pv08, and Pv10. For the white mold isolate tested, 24 resistant lines were identified and the resistance was localized to three different positions on chromosome Pv08. For the powdery mildew local isolate, only 12 resistant lines were identified, and along with the two previous resistance genes on chromosomes Pv04 and Pv11, a new region on chromosome Pv06 was also identified. For root rot caused by Pythium, 31 resistant lines were identified and two main regions were located on chromosomes Pv04 and Pv05. Relevant information for snap bean breeding programs was provided in this work. A total of 20 lines showed resistant or intermediate responses against four or five isolates, which can be suitable for sustainable farm production and could be used as resistance donors. Potential genes and genomic regions to be considered for targeted improvement were provided, including new or less characterized regions that should be validated in future works. Powdery mildew disease was identified as a potential risk for snap bean production and should be considered a main goal in breeding programs.
1 Introduction
Snap bean (green bean or French bean) is a group of common beans (Phaseolus vulgaris L.) whose immature pods can be consumed as a vegetable. The origin of the snap bean is still uncertain. The very complex evolutionary history that characterizes P. vulgaris, including the spread and adaptation processes outside the centers of domestication (Bellucci et al., 2023), makes it difficult to obtain a clear picture. Snap beans may have arisen in Europe as a result of selective pressure on the pod characteristics of dry cultivars (Wade, 1937; Myers and Bagget, 1999). However, an independent derivation from dry beans may have also occurred in the Americas, with snap beans subsequently brought to Europe (Wallace et al., 2018). Most snap beans cultivated today are elite cultivars obtained through intensive breeding programs.
Plant diseases cause significant losses in food crop production worldwide (Ristaino et al., 2021; Singh et al., 2023). Anthracnose (ANT), white mold (WM), powdery mildew (PM), or root rot diseases are the most frequent fungal diseases present throughout bean production areas worldwide (Schwartz et al., 2005).
ANT is caused by the fungus Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.- Scrib., a highly variable pathogen with several reported pathotypes (Padder et al., 2017). Resistance generally follows a qualitative mode of inheritance, with a very specific plant-pathogen interaction. Several ANT resistance genes (Co-) have been described and mapped on all the P. vulgaris chromosomes, except Pv05, Pv06, and Pv10 (Kelly and Vallejo, 2004; Ferreira et al., 2013; Campa et al., 2014);. Three chromosome regions including ANT resistance genes have been identified as complex clusters of resistant genes including closely linked race-specific genes, different allelic series of a Co- gene, or even genes conferring resistance to a diverse selection of pathogens (Kelly and Vallejo, 2004; Miklas et al., 2006; Ferreira et al., 2013; Meziadi et al., 2016). Some of these clusters have been physically positioned in the bean reference genome (Schmutz et al., 2014): the commonly known Co-1 cluster was located at the end position of chromosome Pv01, approximately at 50 Mbp (Richard et al., 2014; 2021; Vaz Bisneta and Gonçalves-Vidigal, 2020). On Pv04 a well-known anthracnose resistance cluster (David et al., 2008, David et al, 2009), usually identified as the Co-3 cluster, was located in a physical position of 1.28–2.04 Mbp (Murube et al., 2019; Vaz Bisneta et al., 2021). Finally, the first anthracnose resistance cluster identified by Creusot et al. (1999) was physically located at the end position of chromosome Pv11, around 45,15–48,66 Mbp (Chen et al., 2018). Using a genome-wide association study (GWAS) a broader range of genetic variation was explored, and new chromosome regions have been recently identified on chromosomes Pv02, Pv05, Pv06, and Pv10 (Wu et al., 2017; Banoo et al., 2020; Vidigal Filho et al., 2020; Almeida et al., 2021; Alvarez-Diaz et al., 2021; Costa et al., 2021; Vaz Bisneta et al., 2021; Simons et al., 2022).
The causal agent of WM in common bean is the fungus Sclerotinia sclerotiorum (Lib.) de Bary. Defense against this pathogen is conditioned by avoidance and physiological resistance mechanisms, both being quantitatively inherited (Miklas et al., 2001). Substantial progress has been made in understanding pathogenic variation, developing screening methods, understanding the genetics of resistance, and developing resistance-related molecular markers (Schwartz and Singh, 2013). Several quantitative trait loci (QTL) for resistance and avoidance have been identified based on bi-parental populations, with most having small to moderate effects, and located in all chromosomes except Pv10 (Schwartz and Singh, 2013; Vasconcellos et al., 2017; Escobar et al., 2022). Recently, three GWAS against WM were conducted and new chromosome regions associated with the resistance were identified, including chromosome Pv10 (Campa et al., 2020; Arkwazee et al., 2022; Escobar et al., 2022).
A broad range of genera and species of the order Erysiphales (Glawe, 2008) can cause PM. Two different species were identified as the common bean PM´s causal agent, Erysiphe poligony DC (Rezende et al., 1999; Schwartz et al., 2005), and Erysiphe diffusa (Cooke & Peck) U. Braun & S. Taka (Almeida et al., 2008). Genetic analyses have shown a qualitative mode of inheritance, with two independent epistatic genes (Pm genes) having been identified (Trabanco et al., 2012; Perez-Vega et al., 2013). One gene confers total resistance and masks the effect of the other gene which confers partial resistance. These genes have been located on chromosomes Pv04 (0,08–0,23 Mbp, candidate gene Phvul.004G00150; Campa and Ferreira, 2017) and Pv11, in a physical position corresponding to that of the ANT resistance cluster Co-2 (Perez-Vega et al., 2013). To date, one GWAS against PM has been conducted in common bean (Binagwa et al., 2021) and confirmed the locations of resistance genes on chromosomes Pv04, and Pv11, as well as identifying a new region on chromosome Pv10 (40,50–40,9 Mbp).
Root rot diseases can be caused by a complex group of fungi species, among which Fusarium solani, Rhizoctonia solani and Pythium ultimum are very frequent. Pythium ultimum (Py) is a soil pathogen that can cause seed rot before germination, pre-emergence and post-emergence damping-off, or root rot (Schwartz et al., 2005). Root rot resistance, including to Py, has been a main target trait for breeding programs in the East African region, where this disease is a primary cause of crop failures among subsistence farmers (Abawi and Pastor Corrales, 1990). Genetic resistance to this disease has been difficult to increase due to (i) its negative association with agronomically accepted white-seeded cultivars; (ii) the polygenic nature of this resistance, being identified as having qualitative and quantitative inheritance; (iii) and its linkage with the P gene responsible for the white seed color (McClean et al., 2018) or with undesirable traits, such as late maturity, low yield potential, and poor pod quality (Deakin, 1974; Schwartz et al., 2005; Campa et al., 2010). White-seeded cultivars with high levels of resistance have not been identified, but within colored seeds, two QTL associated with the emergence rate and with seedling vigor were identified in linkage groups Pv03 and Pv06 (Campa et al., 2010). A GWAS against Py (Dramadri et al., 2020) led to the identification of new genomic regions associated with this resistance on chromosomes Pv01 (6,27 Mbp), Pv02 (25,48 and 44,90 Mbp), Pv04 (41,96 Mbp), Pv05 (29,48 Mbp), and Pv09 (6,25 Mbp).
The control of plant diseases using classical pesticides raises serious concerns about food safety, environmental quality, and pesticide resistance. In the context of sustainable agriculture, which aims to reduce the environmental and climatic impact of primary production, or in the context of the organic farming system, which aims to use only inputs of natural origin, the control of diseases through plant breeding is the most effective strategy. The objective of this work was to investigate the reactions against four fungal diseases, ANT, WM, PM, and root rot caused by Py in a set of 311 snap bean lines (Snap Bean Panel-SBP; García-Fernández et al., 2022) to identify the most resistant materials that will be the more resilient under sustainable farming conditions. A GWAS was also conducted against each of these diseases to identify the genomic regions involved in the resistance. The results provide relevant information regarding the need for snap bean breeding programs in the snap bean market.
2 Material and methods
2.1 Plant material
The SBP contains 311 lines derived from accessions for snap consumption collected in European gene banks, working collections, and seed companies (Ferreira et al., 2021a, Ferreira et al., 2021b; García-Fernández et al., 2022; Supplementary Table S1). A set of 14 common bean lines derived from elite bean varieties including well-known resistance sources against the different diseases, was also used (Supplementary Table S2).
The 311 SBP lines were morphologically characterized in the field by García-Fernández et al. (2022) recording two qualitative plant traits (growth habit, wing flower color, standard flower color), four qualitative pod traits (pod length, pod cross-section, pod color, curvature) and six qualitative seed traits (primary seed color, secondary seed color, shine, longitudinal shape, transversal shape, seed size).
2.2 Resistance test
2.2.1 ANT
Monosporic isolates CL18 and CL124 from local bean crops in northern Spain are maintained at the SERIDA pathogen collection (Ferreira et al., 2008). Monosporic isolate C531 is maintained at the Orsay fungal library and was collected in Costa Rica (Geffroy et al., 1999). Resistance tests were conducted according to standardized resistance tests (Schwartz et al., 1981 Each isolate was evaluated independently in separate experiments, CL18 and CL124 were evaluated at the SERIDA and C531 at the INRAE. Each experiment consisted of all SBP and check lines, with 8–10 seedlings per line sown in trays randomly distributed in the greenhouse. Plant responses were recorded per line 7 days after inoculation. Plants without visible symptoms or showing limited necrotic lesions were scored as resistant, R. Plants with large, sporulating lesions and dead plants were scored as susceptible, S (Van Schoonhoven and Pastor-Corrales, 1987). Resistant lines and those showing unclear reactions were evaluated in a second resistance test.
2.2.2 PM
A local isolate of PM is preserved at the SERIDA, on plants of the S bean susceptible cultivar Xana, grown in spore-proof chambers. Resistance tests were conducted following Trabanco et al. (2012). Two independent resistance tests were carried out in the greenhouse, each one including half of the SBP and the check lines. In each resistant test, two pots per line were sown, with 4–5 plants per pot, randomly distributed in the greenhouse. Lines showing resistant and unclear responses were confirmed in a third resistance test. Plant responses were recorded per line 8–10 days after inoculation as infection type (IT) following a 0–4 scale: plants showing the IT0 response were considered R, plants showing IT1, IT2, or IT3 were considered intermediate (I), and plants showing IT4 were considered S (Trabanco et al., 2012).
2.2.3 WM
The local isolate, WM2, maintained at the SERIDA pathogen collection was used (Pascual et al., 2010). Resistance tests were conducted using the straw method (Petzoldt and Dickson, 1996) and the same experimental design as that used with PM. Disease progression was evaluated 8–10 days after inoculation based on the level of invasion of the main stem, using a 1–9 severity scale, where 1 represents no symptoms and 9 represents an invasion of the third node and plant death. Because WM resistance is quantitatively inherited, each plant in the pot was rated separately and the disease value per pot was calculated as the arithmetic mean of the 4–5 plants. The disease value per line was calculated as the arithmetic mean of the two pots. Outliers were identified using the coefficient of variation (coefficient of variation = (standard deviation/mean) *100). Values equal to or less than 4.5 were considered highly R, values between 4.5 and 7 were considered I, and values above 7 were considered S.
2.2.4 Py
A local Py isolate is maintained at the SERIDA pathogen collection on infected soil maintained at 4°C. Resistance tests were conducted in accordance with Campa et al. (2010). Two independent resistance tests were carried out, each one including half of the SBP and the check lines randomly distributed in the greenhouse. For each line, 10 seeds were inoculated with the fungus and 10 seeds (used as a control) were inoculated with agar plugs only. Disease assessment was evaluated 15 to 18 days after planting (first trifoliate leaf opened) as follows:
Values of 100% of germination were considered highly R, between 90%-30% of germination were considered I, and values equal to or less than 30% of germination were considered S.
2.3 Genotyping
A total of 311 SBP lines were genotyped using the Genotyping by Sequencing (GBS) method (Elshire et al., 2011) optimized following Schröder et al. (2016), based on the use of two enzymes (Taqα1 and MseI). Genomic DNAs from the lines were digested individually with the restriction enzyme combination and post-PCR purified with AMPure XP Beads. Paired-end library sequencing was carried out using the Illumina platform. For the bioinformatics analyses carried out for mapping and variant calling. Cutadapt v3.2 (Martin, 2011) was used to filter and trim low-quality reads and FastQC (Andrews, 2010) was used to check the final quality of the reads. Data from line SBP365 were not available. Each cleaned library was then mapped to the G19833 P. vulgaris v2.0 reference genome (Schmutz et al., 2014) using BWA-mem v0.7.15 (Li, 2013). The Genome Analysis ToolKit v4.1.9.0 (McKenna et al., 2010) was used for variant calling. Duplicated reads were sorted and filtered with Picard (2022) v2.4.1 (http://broadinstitute.github.io/picard) for sequencing duplicates. Variants were then discovered using “HaplotypeCaller” mode for each accession separately. Then the joint genotyping across all accessions was performed using “GenotypeGVCFs” mode. SNPs were extracted using the “SelectVariants” mode and pre-filtered using the “VariantFiltration” mode with the recommended parameters for hard filtering (–filterExpression “ QUAL< 60 || QD< 2.0 || MQ< 40.0 || FS > 60.0 || SOR > 3.0 || MQRankSum< -20.0 || ReadPosRankSum< -8.0”). SNPs were named based on their physical positions in the bean chromosomes (Pv) and their genomic positions in bp.
2.4 GWAS
SNPs were filtered using Tassel v5.2 software (Bradbury et al., 2007) for missing values (< 20%) and minor allele frequency (MAF> 0.05). A total of 16,242 SNPs were considered. The distribution of the SNPs along the 11 chromosomes was visualized with the “CMplot” package of the R project (R core Team 2022). GWAS was conducted based on the multi-locus model FASTmrEMMA (Wen et al., 2018) using the mrMLM package of the R project (Zhang et al., 2020; R Core Team, 2022). A Principal Component Analysis for two components and a Kinship matrix, obtained by the centered-IBS method, were considered to account for multiple levels of relatedness within the lines included in the panel. FASTmrEMMA includes a multi-locus model for true SNP detection at LOD≥ 3 from the associations showing p-values ≤ 0.005 (-log10(p) ≥ 2.3).
3 Results
3.1 Responses to four pathogens
The responses of the 311 SBP lines against the 6 isolates are summarized in Figure 1 (see Supplementary Table S1). The highest number of R lines was observed against the three ANT isolates: 156, 146, and 109 lines were R to CL124, CL18, and C531, respectively. For WM, 24 lines were considered R, 125 lines were I, and 148 were S. For PM, only 12 lines were R, 78 were I, and most lines, 212, were S. For Py, 31 lines were considered R (all with colored seed), 83 were I, and 190 were S.
Figure 1 Stacked bar graph comparing the number of lines identified as R (resistant), I (intermediate) or S (susceptible) for each disease using the SBP. CL124, CL18, and C531: Colletotrichum lindemuthianum isolates. WM, white mold. PM, powdery mildew. Py, Pythium ultimum isolate.
The evaluation of the six isolates was jointly considered in 282 lines (Supplementary Table S1). Considering the response to the three ANT isolates, 150 lines showed resistance to at least two ANT isolates, 53.3% of which were R to the three isolates, 35.3% to CL124 and CL18, 6.0% to CL124 and C531 and 5.3% to CL18 and C531. For resistance to ANT isolates and PM, 78% of the PM-R lines were also R to C531, and 89% were also R to CL18 and CL124. Considering responses to the ANT isolates and WM, 16%, 21%, and 33% of the lines showing R to WM were also R to CL18, C531, and CL124 isolates, respectively. Considering responses to Py and the ANT isolates, 17% of the lines were R to Py and CL124, and 10% were R to Py and CL18 and C531 isolates. No lines were R to WM and to PM, nor to PM and Py, and only 12% of the lines showed resistance against WM and Py.
3.2 Morphological characterizations of the most promising lines
In total, 20 lines were considered the most promising materials of the SBP based on their R or I responses against more than three isolates (Table 1). Following the results of García-Fernández et al. (2022), the morpho-agronomic traits of these lines are shown in Table 1 (Supplementary Table S1). Most lines showed determinate growth habits, and only five showed indeterminate climbing habits. For pod color, the most frequent phenotype was green or green-mottled color, with only one line having yellow pods (SBP372) and one line having purple pods (SBP096). Most lines (11) showed colored seeds, with only single lines having either red or black seeds, and nine lines were white-seeded. Different pod and seed sizes were homogeneously distributed (Table 1).
Table 1 Morpho-agronomic characterization of the SBP lines showing resistance or intermediate response to more than three isolates.
3.3 GWAS
The 16,242 SNPs considered after filtering were homogeneously distributed across the 11 bean chromosomes (Supplementary Figure S1). The GWAS conducted in the SBP led to the identification of 20 significantly associated SNPs on chromosomes Pv02, Pv04, Pv05, Pv06, Pv08, Pv10, and Pv11 (Figures 2, 3; Table 2). Four SNPs were identified for the ANT isolate CL124 on chromosomes Pv02 (S02_6165835), Pv04 (S04_1022360), Pv10 (S10_3752313), and Pv11 (S11_49692310). For the ANT isolate CL18, two SNPs were identified, one at chromosome Pv08 (S08_58508678) and the other at Pv11 (S11_49692310), at the same position as for CL124. For the ANT isolate C531, five SNPs were identified on chromosomes Pv04 (S04_1022420), Pv06 (S06_23139662), and Pv10 (S10_2713358, S10_21453937, and S10_41922734). For the WM isolate, three SNPs were identified on chromosome Pv08 (S08_1528630, S08_25049002, S08_59153321). For PM, four SNPs were identified on chromosomes Pv04 (S04_235050), Pv06 (S06_16484868), Pv08 (S08_2467368), and Pv11 (S11_52547132). The SNP S08_2467368 was considered a false positive because it was associated with the allele N, which corresponded to a missing value (Table 2). Finally, two SNPs were identified for Py on chromosomes Pv04 (S04_25514460) and Pv05 (S05_29460608).
Figure 2 Manhattan and QQ plots obtained in the GWAS conducted in the SBP for anthracnose isolates (CL124, CL18, C531), white mold (WM), powdery mildew (PM), and Pythium ultimum (Py). The two threshold lines considered in the FASTmrEMMA method are indicated, in orange color for -log(p)≥ 2.3 and in grey color for LOD score ≥ 3. Pink dots indicate the SNPs significantly associated with each isolate. Correspondence with the previously described Co-2 and Co-3 clusters and with the Pm genes is indicated.
Figure 3 SNPs per chromosome associated with the different isolates evaluated in the Snap Bean Panel: anthracnose isolates (CL18, CL124, C531; blue color) white mold (WM; green color), powdery mildew (PM; red color), and Pythium ultimum. (Py; orange color). The specific isolate associated with each SNP is indicated in superscript. Scale bar represents chromosome size in Mbp.
Table 2 SNPs that were significantly associated with each isolate in the GWAS performed with the FASTmrEMMA method and using the SBP.
4 Discussion
The SBP includes a wide diversity of genotypes collected from European gene banks, working collections, and seed companies, and it can be considered to be representative of the level of diversity in European snap beans (García-Fernández et al., 2022). Disease resistance is a relevant goal for the adaptation of snap bean crops to more sustainable agriculture, which aims to reduce the environmental and climatic impact of primary production. The control of diseases through plant breeding is the most effective strategy, but further knowledge is still required to meet this objective. Improving our understanding of resistance mechanisms and the identification of new R genes are important to mitigate the risk of future disease outbreaks related to climate scenarios that threaten food security in many areas of the world (Ristaino et al., 2021; Singh et al., 2023).
The largest number of R lines was observed against the three ANT isolates evaluated, CL18, CL124, and C531. Traditionally, major achievements in snap bean breeding programs have been the modification of pod traits, seed sizes, and yield traits, but also resistance to diseases. ANT is regarded as a serious disease of common bean due to its seed-borne nature and pathogenic variability (Schwartz et al., 2005). Breeding to protect common bean against ANT has been intense since the beginning of the 20th century, therefore, it is not surprising to identify a high prevalence of resistance in the SBP.
In the GWAS for ANT isolates CL124 and CL18, the same SNP at the end of chromosome Pv11 was identified (S11_49692310). In common bean, this chromosomal region has been described as a complex cluster of R genes encoding nucleotide-binding leucine-rich repeat proteins (NB-LRR) that protect against different diseases or different strains of the same pathogen (Creusot et al., 1999; Campa et al., 2014; Meziadi et al., 2016). One of the first ANT resistance genes identified in common bean, the Co-2 gene, has been mapped to this region (Mastenbroek, 1960). This was the first ANT resistance locus to be ‘tagged’ with molecular markers (Adam-Blondon et al., 1994), and has been extensively used in breeding programs since 1960 (Mastenbroek, 1960; Ferreira et al., 2012). For this reason, it was highly represented within the SBP lines, with several of them being derived from commercial varieties. Isolates CL124 and CL18 correspond to race 3 and race 38, respectively (Ferreira et al., 2008). ANT resistance genes against these two races have been located previously in the Co-2 cluster using different R genotypes (Rodriguez-Suarez et al., 2007; Campa et al., 2011, Campa et al, 2014). This is in agreement with the results obtained in the present work. For ANT isolate C531, the SNP located at Pv04 is close (60 bp apart) to the SNP identified for CL124 isolate, therefore they may correspond to closely linked genes of the same cluster. In this region of the Pv04 chromosome, the ANT resistance cluster Co-3 has also been located (Meziadi et al., 2016; Murube et al., 2019). The Co-3 gene, originally known as Mexique1, was first described by Bannerot (1965) and was one of the first loci having described alleles that conditioned resistance to ANT. The breeding value of Co-3 has been described as very specific for certain locations and was proposed to be included in pyramiding breeding programs (Kelly and Vallejo, 2004). One response to the short-term protection conferred by single genes is to pyramid specific R genes into one plant (Mundt, 1990). Here, the genes Co-2 and Co-3 were identified against the isolate CL124 and have probably been introgressed and pyramided in different breeding programs against ANT. Especially remarkable is the implication of chromosome Pv10 in the resistance to isolate C531 at physical positions 2.71, 21.5, and 41.92 Mbp. The Pv10 chromosome has not been previously identified in the resistance to ANT when genetic analyses were conducted using biparental populations. However, since 2016, the use of GWAS in the study of ANT resistance has led to the exploration of a broader range of genetic variation and the identification of new chromosomal regions (Zuiderveen et al., 2016). Recently, using GWAS, genomic regions on Pv10 were associated with ANT resistance (Banoo et al., 2020; Vidigal Filho et al., 2020; Simons et al., 2022), which is in agreement with the results obtained here. Furthermore, Wu et al. (2017) identified NB-LRR genes on chromosome Pv10 that are related to resistance against different diseases, including ANT. Other lesser-known regions for ANT resistance were identified in this work on chromosomes Pv02 (6.16 Mbp), Pv06 (23.13 Mbp), and Pv08 (58.50 Mbp) for CL124, C531, and CL18 isolates, respectively. At the start position of chromosome Pv02 (4.77 Mbp), an association with ANT resistance was recently identified and could correspond with the association identified here (Banoo et al., 2020). The region identified on chromosome Pv08 for the CL18 isolate may correspond to the one identified by Almeida et al. (2021) at 58.76 Mbp. In summary, some of the genomic regions identified in this work as being associated with resistance against ANT were previously identified in other GWAS. They represent minor R regions that could have important breeding values and should be validated in future works.
WM is endemic and widespread in North and South American countries, many of which have developed intensive breeding programs against the disease. The development of R bean cultivars against WM has been slow owing to the low number of R sources identified and the complex nature of the resistance, having both physiological and avoidance mechanisms. In this work, 24 lines were R to WM, and they can be used as R sources in snap breeding programs. Among them, only five lines were derived from commercial varieties (Supplementary Table S1), highlighting the importance of characterizing different genetic resources, such as landraces and/or wild materials, to exploit their diversity in breeding.
Three regions involved in the resistance to WM were identified on chromosome Pv08 (1.52, 25.04, 59.15 Mbp). Previously, several QTL and one meta-QTL (WM8.3; 37.32–46.73 Mbp) were identified on chromosome Pv08 based on biparental populations and GWAS (Vasconcellos et al., 2017; Campa et al., 2020). However, the physical positions of the regions identified in the present work do not correspond to those identified previously. The results may be difficult to compare owing to the different methods (fields or controlled conditions) used for the WM evaluations. In general, it has been assumed that WM evaluations under controlled conditions reveal physiological resistance, whereas field evaluations measure physiological resistance and disease avoidance mechanisms (Miklas et al., 2001). These new R regions could have important breeding values and should be validated in future works.
A low frequency of resistance against PM was observed, with only 12 SBP R lines identified. PM is a widespread plant disease with a rapidly increasing geographical distribution (Glawe, 2008). This expansion has been related to global climate change, because PM pathogen populations evolve quickly, showing high rates of variation due to the coexistence of sexual and asexual stages and the high dispersal capability. Here, the low frequency of observed resistance indicates the risk that can arise from this disease in the snap bean market and in sustainable farming systems. Increasing the resistance to PM should be a major goal of future breeding programs. The GWAS revealed three regions for PM resistance on chromosomes Pv04 (0.23 Mbp), Pv06 (16.49 Mbp), and Pv11 (52.54 Mbp). The regions on Pv04 and Pv11 corresponded with previously identified and mapped major resistance genes for PM (Pm genes; Trabanco et al., 2012; Perez-Vega et al., 2013; Campa et al., 2014). These regions at the beginning of chromosome Pv04 and the end of chromosome Pv11 are close to the physical positions of ANT clusters Co-3 and Co-2, respectively. This could explain the strong association identified in this work between the responses to PM and ANT isolates. However, the low number of R lines identified against PM, in contrast to the high number of ANT-R lines identified, also suggests that the genes conferring resistance against the pathogens are different. In fact, these two regions correspond to large clusters of NB-LRR genes, the major class of disease-resistance genes in plants (Meziadi et al., 2016). The new region identified on chromosome Pv06 (16.49 Mbp) against PM should be validated because it could have an important breeding value.
For Py disease, only 31 SBP lines showed high levels of resistance, all of them having colored seeds. Among the R lines, five were derived from commercial lines, whereas the remaining were derived from landraces or from unknown origins. Thus, they could be new interesting R sources to be exploited in snap bean breeding programs. The GWAS identified two regions on chromosomes Pv04 (25.51 Mbp) and Pv05 (29.46 Mbp). The physical position of the region identified on Pv05 corresponds to the one identified by Dramadri et al. (2020), confirming that this region plays a role in the genetic control of Py resistance. The other region at Pv04 was identified for the first time and should be confirmed in future works. Two other chromosomal regions on Pv02 and Pv07 have not been identified as significantly associated by the FASTmrEMMA method but are noteworthy regions showing log10(p) ≥ 3. On Pv02, a region between 43.92 and 44.21 Mbp was tagged with four SNPs having high -log(p) values in the GWAS: SNP S02_43921938 - log10(p)=4; S02_44213587 - log10(p)=4; S02_44217242 - log10(p)=3; S02_44217294 - log10(p)=4. This region was also identified as being associated with Py resistance in the GWAS conducted by Dramadri et al. (2020). Py resistance is often associated with colored seeds, and at this position of chromosome Pv02, between 33.7–48.5 Mbp, was located the seed coat color gene B, responsible for gray-brown colors and proposed to regulate simultaneously the flavonoid (color) and isoflavonoid (resistance) pathways (Beninger et al., 2000; Hagerty et al., 2016; Campa et al., 2023). Another seed coat color gene, the Rk, was located at this position, linked to B. Rk was described to be involved in the synthesis of red or garnet brown seed colors, and was proposed as a modifying color gene that affects the amount and type of tannins produced in the seed coat (Prakken, 1972). The results obtained in this work suggest that the seed coat color genes B or Rk could be involved in the resistance against Py. On Pv07, two SNPs showed high log(p) values (S07_28657174, log(p)=3; S07_28657186, log(p)=3) in a position corresponding to that of P gene (candidate gene Phvul.007G171333; location 28,752,131.28,766,155 bp; McClean et al., 2018) responsible for white seed color and linked to Py response (Deakin, 1974; Campa et al., 2010). White-seeded lines showing high levels of resistance have not been identified. Lines combining total or I resistance to WM and PM, or to PM and Py, have not been identified among all the isolates evaluated. This indicates that such resistance through breeding requires the pyramiding of genetic resistance to different pathogens and the introgression of resistance in specific pod phenotypes or snap bean market classes. For instance, in the groups of lines having flat cross-sections and very long green-colored pods (Type Romano), round cross-sections, or long yellow-colored pods (Yellow wax), a lack of genotypes R to PM was detected (Supplementary Table S1). In total, 20 lines were identified as the most promising, showing R or I responses to more than three isolates. Among them, the line SBP240 is very promising, with resistance to ANT and PM and moderate resistance to WM. This set of 20 lines provides different genotypes with high levels of resistance and wide phenotypic diversity that can be used in sustainable farm production systems with low input needs. These lines are also potential resistance sources for snap bean breeding programs. The data revealed by the GWAS can help breeders to identify useful gene/s or genomic regions for targeted breeding programs.
5 Conclusions
This work led to the identification of the most promising snap bean lines for sustainable farm production and potential resistance sources for breeding programs focused on developing varieties with broad disease resistance. The smallest number of R lines was found against PM which poses a high risk owing to the rapid and easy spread of the causal pathogen. Likewise, this study provides relevant information about the needs of future breeding programs for this group of common beans and the potential genes and genomic regions to be considered for targeted improvement. A GWAS of the SBP showed the involvement of well-known chromosomal regions in the R responses against the four pathogens and identified new, or less characterized genomic regions that should be validated in future works. This highlighted the importance of characterizing different genetic resources, such as landraces and/or wild materials, to exploit their diversity in breeding. Breeding resistance to ANT in snap beans is mainly conducted using the Co-2 and Co-3 clusters. For future breeding programs, the use of other ANT resistance genes or clusters should be considered.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: (https://www.ncbi.nlm.nih.gov/, PRJNA1076096).
Author contributions
AC: Writing – original draft, Methodology, Formal analysis, Data curation, Conceptualization. VG: Writing – review & editing, Data curation. EB: Writing – review & editing, Data curation. AN: Writing – review & editing, Data curation. RP: Writing – review & editing, Supervision, Project administration, Funding acquisition. JF: Writing – review & editing, Supervision, Funding acquisition, Data curation.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported in part by the two European Union’s Horizon 2020 Research and Innovation Projects BRESOV (Grant Agreement no. 774244) and INCREASE (Grant Agreement no. 862862).
Acknowledgments
The authors thank JA Poladura, F Díez, and Seyed Adib Sheikhahma for their assistance.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1386877/full#supplementary-material
References
Abawi, G. S., Pastor Corrales, M. A. (1990). Root rots of beans in Latin America and Africa: Diagnosis, research methodologies, and management strategies (Cali, Colombia: Centro Internacional de Agricultura Tropical (CIAT), 114.
Adam-Blondon, A. F., Sevignac, M., Bannerot, H., Dron, M. (1994). SCAR, RAPD and RFLP markers linked to a dominant gene (ARE) conferring resistance to anthracnose in common bean. Theor. Appl. Genet. 88, 865–870. doi: 10.1007/BF01253998
Almeida, A. M. R., Binneck, E., Piuga, F. F., Marin, S. R. R., Ribeiro do Valle, P. R. Z., Silveir, C. A., et al. (2008). Characterization of powdery mildews strains from soybean, bean, sunflower, and weeds in Brazil using rDNA-ITS sequences. Trop. Plant Pathol. 33, 20–26. doi: 10.1590/S1982-56762008000100004
Almeida, C. P., Paulino, J. F. C., Barbosa, C. C. F., Gonçalves, G. M. C., Fritsche-Neto, R., Carbonell, S. A. M., et al. (2021). Genome-wide association mapping reveals race-specific SNP markers associated with anthracnose resistance in carioca common beans. PloS One 16, e0251745. doi: 10.1371/journal.pone.0251745
Alvarez-Diaz, J. C., Richard, M. M. S., Thareau, V., Teano, G., Paysant-Le-Roux, C., Rigaill, G., et al. (2021). Genome-wide identification of key components of RNA silencing in two phaseolus vulgaris genotypes of contrasting origin and their expression analyses in response to fungal infection. Genes 2713, 64. doi: 10.3390/genes13010064
Andrews, S. (2010). FastQC: A Quality Control tool for High Throughput Sequence Data. Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
Arkwazee, H. A., Wallace, L. T., Hart, J. P., Griffiths, P. D., Myers, J. R. (2022). Genome-wide association study (GWAS) of white mold resistance in snap bean. Genes 613, 2297. doi: 10.3390/genes13122297
Bannerot, H. (1965). Résults de l’infection d’une collection de haricots par six races physiologiques d’anthracnose. Ann. Amélior. Des. Plantes (Paris) 15, 201–222.
Banoo, A., Nabi, A., Rasool, R. S., Mahiya-Farooq, S. M. D., Ahmad, M., Sofi, P. A., et al. (2020). North-western himalayan common beans: population structure and mapping of quantitative anthracnose resistance through genome-wide association study. Front. Plant Sci. 611. doi: 10.3389/fpls.2020.571618
Bellucci, E., Benazzo, A., Xu, C., Bitocchi, E., Rodriguez, M., Alseekh, A., et al. (2023). Selection and adaptive introgression guided the complex evolutionary history of the European common bean. Nat. Commun. 14, 1908. doi: 10.1038/s41467-023-37332-z
Beninger, C. W., Hosfield, G. L., Bassett, M. J., Owens, S. (2000). Chemical and morphological expression of the B and Asp seedcoat genes in Phaseolus vulgaris. J. Amer Soc. Hort Sci. 125, 52–58. doi: 10.21273/JASHS.125.1.52
Binagwa, P. H., Traore, S. M., Egnin, M., Bernard, G. C., Ritte, I., Mortley, D., et al. (2021). Genome-wide identification of powdery mildew resistance in common bean (Phaseolus vulgaris L.). Front. Genet. 2212. doi: 10.3389/fgene.2021.673069
Bradbury, P. J., Zhang, Z., Kroon, D. E., Casstevens, T. M., Ramdoss, Y., Buckler, E. S. (2007). TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics 23, 2633–2635. doi: 10.1093/bioinformatics/btm308
Campa, A., Ferreira, J. J. (2017). Gene coding for an elongation factor is involved in resistance against powdery mildew in common bean. Theor. Appl. Genet. 130, 849–860. doi: 10.1007/s00122-017-2864-x
Campa, A., García-Fernández, C., Ferreira, J. J. (2020). Genome-wide association study (GWAS) for resistance to sclerotinia sclerotiorum in common bean. Genes 11, 1496. doi: 10.3390/genes11121496
Campa, A., Giraldez, R., Ferreira, J. J. (2011). Genetic analysis of the resistance to eight anthracnose races in the common bean differential cultivar Kaboon. Phytopathology 101, 757–764. doi: 10.1094/PHYTO-11-10-0296
Campa, A., Perez-Vega, E., Pascual, A., Ferreira, J. J. (2010). Genetic analysis and molecular mapping of quantitative trait loci in common bean against Pythium ultimum. Phytopathology 100, 1315–1320. doi: 10.1094/PHYTO-06-10-0161
Campa, A., Rodríguez Madrera, R., Jurado, M., García-Fernández, C., Suárez Vallés, B. (2023). Genome-wide association study for the extractable phenolic profile and coat color of common bean seeds (Phaseolus vulgaris L.). BMC Plant Biol. 23, 158. doi: 10.1186/s12870-023-04177-z
Campa, A., Rodriguez-Suarez, C., Giraldez, R., Ferreira, J. J. (2014). Genetic analysis of the response to eleven Colletotrichum lindemuthianum races in a RIL population of common bean (Phaseolus vulgaris L.). BMC Plant Biol. 14, 115. doi: 10.1186/1471-2229-14-115
Chen, N. W. G., Thareau, V., Ribeiro, T., Magdelenat, G., Ashfield, T., Innes, R. W., et al. (2018). Common bean subtelomeres are hot spots of recombination and favor resistance gene evolution. Front. Plant Sci. 149. doi: 10.3389/fpls.2018.01185
Costa, L. C., Nalin, R. S., Dias, M. A., Ferreira, M. E., Song, Q., Pastor-Corrales, M. A., et al. (2021). Different loci control resistance to different isolates of the same race of Colletotrichum lindemuthianum in common bean. Theor. Appl. Genet. 134, 543–556. doi: 10.1007/s00122-020-03713-x
Creusot, F., Macadre, C., Cana, E. F., Riou, C., Geffroy, V., Sevignac, M., et al. (1999). Cloning and molecular characterization of three members of the NBS-LRR subfamily located in the vicinity of the Co-2 locus for anthracnose resistance in Phaseolus vulgaris. Genome 42, 254–264. doi: 10.1139/g98-134
David, P., Chen, N. W. G., Pedrosa-Harand, A., Thareau, V., Sevignac, M., Cannon, S. B., et al. (2009). A nomadic subtelomeric disease resistance gene cluster in common bean. Plant Physiol. 151, 1048–1065. doi: 10.1104/pp.109.142109
David, P., Sévignac, M., Thareau, V., Catillon, Y., Kami, J., Gepts, P., et al. (2008). BAC end sequences corresponding to the B4 resistance gene cluster in common bean: a resource for markers and synteny analyses. Mol. Genet. Genomics 280, 521–533. doi: 10.1007/s00438-008-0384-8
Deakin, J. R. (1974). Association of seed color with emergence and seed yield of snap beans. J. Am. Soc Hortic. Sci. 99, 110–114. doi: 10.21273/JASHS.99.2.110
Dramadri, I. O., Amongi, W., Kelly, J. D., Mugisha Mukankusi, C. (2020). Genome-wide association analysis of resistance to Pythium ultimum in common bean (Phaseolus vulgaris). Plant Breed. 139, 1168–1180. doi: 10.1111/pbr.12855
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PloS One 6, e19379. doi: 10.1371/journal.pone.0019379
Escobar, E., Oladzad, A., Simons, K., Miklas, P., Lee, R. K., Schroder, S., et al. (2022). New genomic regions associated with white mold resistance in dry bean using a MAGIC population. Plant Genome 15, e20190. doi: 10.1002/tpg2.20190
Ferreira, J. J., Campa., A., García-Fernández, C. (2021a). Pod phenotypic variation in the Snap Bean Panel (SBP). Zenodo. doi: 10.5281/zenodo.5557139
Ferreira, J. J., Campa, A., García-Fernández, C. (2021b). Seed diversity catalog of the Snap Bean Panel. Zenodo. doi: 10.5281/zenodo.5557174
Ferreira, J. J., Campa, A., Kelly, J. D. (2013). “Organization of genes conferring resistance to anthracnose in common bean,” in Translational genomic for crop breeding. Eds. Varshney, R. K., Tuberosa, R. (Wiley Blackwell, USA: John Wiley & Sons), 151–182. doi: 10.1002/9781118728475
Ferreira, J. J., Campa, A., Perez-Vega, E., Giraldez, R. (2008). Reaction of a bean germplasm collection against five races of colletotrichum lindemuthianum identified in northern Spain and implications for breeding. Plant Dis. 92, 705–708. doi: 10.1094/PDIS-92-5-0705
Ferreira, J. J., Campa, A., Perez-Vega, E., Rodriguez-Suarez, C., Giraldez, R. (2012). Introgression and pyramiding into common bean market class fabada of genes conferring resistance to anthracnose and potyvirus. Theor. Appl. Genet. 124, 777–788. doi: 10.1007/s00122-011-1746-x
García-Fernández, C., Jurado, M., Campa, A., Brezeanu, C., Geffroy, V., Bitocchi, E., et al. (2022). A core set of snap bean genotypes was established by phenotyping a large panel collected in Europe. Plants 2211, 577. doi: 10.3390/plants11050577
Geffroy, V., Sicard, D., de Oliveira, J. C., Sévignac, M., Cohen, S., Gepts, P., et al. (1999). Identification of an ancestral resistance gene cluster involved in the coevolution process between Phaseolus vulgaris and its fungal pathogen Colletotrichum lindemuthianum. Mol. Plant Microbe Interact. 12, 774–784. doi: 10.1094/MPMI.1999.12.9.774
Glawe, D. A. (2008). The powdery mildews. A review of the world’s most familiar (yet poorly known) plant pathogens. Ann. Rev. Phytopathol. 46, 27–51doi: 10.1371/journal.pone.0235215
Hagerty, C. H., Cuesta-Marcos, A., Cregan, P., Song, Q., McClean, P., Myers, J. R. (2016). Mapping snap bean pod and color traits, in a dry bean × Snap bean recombinant inbred population. J. Amer. Soc Hortic. Sci. 141, 131–138. doi: 10.21273/JASHS.141.2.131
Kelly, J. D., Vallejo, V. A. (2004). A comprehensive review of the major genes conditioning resistance to anthracnose in common bean. HortScience 39, 1196–1207. doi: 10.21273/HORTSCI.39.6.1196
Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv, 1303.3997v2. doi: 10.6084/M9.FIGSHARE.963153.V1
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17, 10. doi: 10.14806/ej.17.1.200
Mastenbroek, C. A. (1960). A breeding program for resistance to anthracnose in dry shell haricot beans, based on a new gene. Euphytica 9, 177–184. doi: 10.1007/BF00022219
McClean, P. E., Bett, K. E., Stonehouse, R., Lee, R., Pflieger, S., Moghaddam, S. M., et al. (2018). White seed color in common bean (Phaseolus vulgaris) results from convergent evolution in the (pigment) gene. New Phytol. 219, 1112–1123. doi: 10.1111/nph.15259
McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi: 10.1101/gr.107524.110
Meziadi, C., Richard, M. M. S., Derquennes, A., Thareau, V., Blanchet, S., Gratias, A., et al. (2016). Development of molecular markers linked to disease resistance genes in common bean on wholes genome sequence. Plant Sci. 242, 351–357. doi: 10.1016/j.plantsci.2015.09.006
Miklas, P. N., Delorme, R., Johnson, W. C., Gepts, P. (2001). QTL conditioning physiological resistance and avoidance to white mold in dry bean. Crop Sci. 41, 309–315. doi: 10.2135/cropsci2001.412309x
Miklas, P. N., Kelly, J. D., Beebe, S. E., Blair, M. W. (2006). Common bean breeding for resistance against biotic and abiotic stresses: From classical to MAS breeding. Euphytica 147, 105–131. doi: 10.1007/s10681-006-4600-5
Mundt, C. C. (1990). Probability of mutation to multiple virulence and durability of resistance gene pyramids: further comments. Phytopathology 80, 221–223. doi: 10.1094/Phyto-80-221
Murube, E., Campa, A., Ferreira, J. J. (2019). Integrating genetic and physical positions of the anthracnose resistance genes described in bean chromosomes Pv01 and Pv04. PloS One 14, e0212298. doi: 10.1371/journal.pone.0212298
Myers, J. R., Bagget, J. R. (1999). “Improvement of snap beans,” in Common bean improvement in the twenty-first century. Ed. Singh, S. P. (Germany: Springer), 289–329.
Padder, B. A., Sharma, P. N., Awale, H. E., Kelly, J. D. (2017). Colletotrichum lindemuthianum, the causal agent of bean anthracnose. J. Plant Pathol. 99, 317–330. doi: 10.4454/jpp.v99i2.3867
Pascual, A., Campa, A., Perez-Vega, E., Giraldez, R., Miklas, P. N., Ferreira, J. J. (2010). Screening common bean for resistance to four sclerotinia sclerotiorum isolates collected in northern Spain. Plant Dis. 947, 885–890. doi: 10.1094/PDIS-94-7-0885
Perez-Vega, E., Trabanco, N., Campa, A., Ferreira, J. J. (2013). Mapping of two genes conferring resistance to powdery mildew in common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 126, 1503–1512. doi: 10.1007/s00122-013-2068-y
Petzoldt, R., Dickson, M. H. (1996). Straw test for resistance to white mold in beans. Annu. Rep. Bean Improv. Coop. 39, 142–143.
Picard. (2022). Broad institute. GitHub. Available at: http://broadinstitute.github.io/picard.
Prakken, R. (1972). Inheritance of colours in Phaseolus vulgaris L. III. On genes for red seedcoat colour and a general synthesis. Meded. Landbouwhogeschool Wageningen 72–29, 1–82.
R Core Team. (2022). R: A language and environment for statistical computing (Vienna, Austria: Foundation for Statistical Computing). Available at: https://www.R-project.org/ (Accessed September 2022).
Rezende, V. F., Patto, R. M. A., Corte, H. (1999). Genetic control of common bean (Phaseolus vulgaris L.) resistance to powdery mildew (Erysiphe polygoni). Genet. Mol. Biol. 22, 233–236. doi: 10.1590/S1415-47571999000200017
Richard, M. M. S., Gratias, A., Alvarez Díaz, J. C., Thareau, V., Pflieger, S., Meziadi., C., et al. (2021). A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum. J. Exp. Bot. 472, 3569–3581. doi: 10.1093/jxb/erab082
Richard, M. M. S., Pflieger, S., Sévignac, M., Thareau, V., Blanchet, S., Li, Y., et al. (2014). Fine mapping of Co-x, an anthracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean. Theor. Appl. Genet. 127, 1653–1666. doi: 10.1007/s00122-014-2328-5
Ristaino, J. B., Anderson, P. K., Bebber, D. P., Brauman, K. A., Cunniffe, N. J., Fedoroff, N. V., et al. (2021). The persistent threat of emerging plant disease pandemics to global food security. Proc. Natl. Acad. Sci. U.S.A. 118, e2022239118. doi: 10.1073/pnas.2022239118
Rodriguez-Suarez, C., Méndez-Vigo, B., Pañeda, A., Ferreira, J. J., Giraldez, R. (2007). A genetic linkage map of Phaseolus vulgaris L. and localization of genes for specific resistance to six races of anthracnose (Colletotrichum lindemuthianum). Theor. Appl. Genet. 114, 713–722. doi: 10.1007/s00122-006-0471-3
Schmutz, J., McClean, P., Mamidi, S., Wu, G. A., Cannon, S. B., Grimwood, J., et al. (2014). A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 46, 707–713. doi: 10.1038/ng.3008
Schröder, S., Mamidi, S., Lee, R., McKain, M. R., McClean, P. E., Osorno, J. M. (2016). Optimization of genotyping by sequencing (GBS) data in common bean (Phaseolus vulgaris L.). Mol. Breed. 36, 6. doi: 10.1007/s11032-015-0431-1
Schwartz, H. F., Pastor Corrales, M. A., Correa Victoria, F. J., Ospina, O., Héctor, F., Flor, M., et al. (1981). La antracnosis del fríjol y su control (Cali, Colombia: CIAT (Centro Internacional de Agricultura Tropical).
Schwartz, H. F., Singh, S. P. (2013). Breeding common bean for resistance to white mold: A review. Crop Sci. 53, 1832–1844. doi: 10.2135/cropsci2013.02.0081
Schwartz, H. F., Steadman, J. R., Hall, R., Forster, R. L. (2005). Compendium of bean diseases. 2nd edition (USA: The American Phytopathological Society).
Simons, K. J., Schröder, S., Oladzad, A., McClean, P. E., Conner, R. L., Penner, W. C., et al. (2022). Modified screening method of middle American dry bean genotypes reveals new genomic regions on Pv10 associated with anthracnose resistance. Front. Plant Sci. 1513. doi: 10.3389/fpls.2022.1015583
Singh, B. K., Delgado-Baquerizo, M., Egidi, E., Guirado, E., Leach, J. E., Liu, H., et al. (2023). Climate change impacts on plant pathogens, food security and paths forward. Nat. Rev. Microbiol. 21, 640–656. doi: 10.1038/s41579-023-00900-7
Trabanco, N., Perez-Vega, E., Campa, A., Rubiales, D., Ferreira, J. J. (2012). Genetic resistance to powdery mildew in common bean. Euphytica 186, 875–882. doi: 10.1007/s10681-012-0663-7
Van Schoonhoven, A., Pastor-Corrales, M. A. (1987). Standard system for the evaluation of bean germplasm (Cali, Colombia: CIAT (Centro Internacional de Agricultura Tropical).
Vasconcellos, R. C. C., Oraguzie, O. B., Soler, A., Arkwazee, H., Myers, J. R., Ferreira, J. J., et al. (2017). Meta-QTL for resistance to white mold in common bean. PloS One 12, e0171685. doi: 10.1371/journal.pone.0171685
Vaz Bisneta, M., Gonçalves-Vidigal, M. C. (2020). Integration of anthracnose resistance loci and RLK and NBS-LRR-encoding genes in the Phaseolus vulgaris L. Genome. Crop Sci. 60, 2901–2918. doi: 10.1002/csc2.20288
Vaz Bisneta, M., Gonçalves-Vidigal, M. C., Filho, V., Ferreira Elias, J. C., Valentini, G., Lima, L. R. L., et al. (2021). New genomic regions for resistance to anthracnose (Colletotrichum lindemuthianum) through GBS-based genome-wide association study in common bean (Phaseolus vulgaris). WJARR 12, 020–040. doi: 10.30574/wjarr.2021.12.1.0493
Vidigal Filho, P. S., Gonçalves-Vidigal, M. C., Vaz Bisneta, M., Souza, V. B., Gilio, T. A. S., Calvi, A. A., et al. (2020). Genome-wide association study of resistance to the anthracnose and angular leaf spot diseases in Brazilian Mesoamerican and Andean common bean cultivars. Crop Sci. 12, 440–447. doi: 10.1002/csc2.20308
Wade, B. L. (1937). “Breeding and improvement of peas and beans,” in Yearbook of agriculture (United States Government Printing Office, Washington, DC, USA), 251–282.
Wallace, L., Arkwazee, H., Vining, K., Myers, J. R. (2018). Genetic diversity within snap beans and their relation to dry beans. Genes 9, 587. doi: 10.3390/genes9120587
Wen, Y. J., Zhang, H., Ni, Y. L., Huang, B., Zhang, J., Feng, J. Y., et al. (2018). Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief Bioinform. 19, 700–712. doi: 10.1093/bib/bbw145
Wu, J., Zhu, J., Wang, L., Wang, S. (2017). Genome-wide association study identifies NBS-LRR-encoding genes related with anthracnose and common bean bacterial blight in the common bean. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01398
Zhang, Y. W., Tamba, C. L., Wen, Y.-J., Li, P., Ren, W.-L., Ni, Y. L., et al. (2020). mrMLM v4.0.2: an R platform for multi-locus genome-wide association studies. Genom. Proteom. Bioinf. 18, 481–487. doi: 10.1016/j.gpb.2020.06.006
Keywords: Phaseolus vulgaris, resistance genes, anthracnose, powdery mildew, Pythium ultimum, white mold, genome-wide association mapping
Citation: Campa A, Geffroy V, Bitocchi E, Noly A, Papa R and Ferreira JJ (2024) Screening for resistance to four fungal diseases and associated genomic regions in a snap bean diversity panel. Front. Plant Sci. 15:1386877. doi: 10.3389/fpls.2024.1386877
Received: 16 February 2024; Accepted: 22 May 2024;
Published: 11 June 2024.
Edited by:
Frédéric Marsolais, Agriculture and Agri-Food Canada (AAFC), CanadaReviewed by:
James Myers, Oregon State University, United StatesMaria Celeste Gonçalves-Vidigal, Universidade Estadual de Maringá, Brazil
Copyright © 2024 Campa, Geffroy, Bitocchi, Noly, Papa and Ferreira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Campa, acampa@serida.org