Integrating genetic and physical positions of the anthracnose resistance genes described in bean chromosomes Pv01 and Pv04

The Kirkhouse Trust: Successes and Challenges in Twenty Years of Supporting Independent, Contemporary Grain Legume Breeding Projects in India and African Countries

This manuscript reviews two decades of projects funded by the Kirkhouse Trust (KT), a charity registered in the UK. KT was established to improve the productivity of legume crops important in African countries and in India. KT's requirements for support are: (1) the research must be conducted by national scientists in their home institution, either a publicly funded agricultural research institute or a university; (2) the projects need to include a molecular biology component, which to date has mostly comprised the use of molecular markers for the selection of one or more target traits in a crop improvement programme; (3) the projects funded are included in consortia, to foster the creation of scientific communities and the sharing of knowledge and breeding resources. This account relates to the key achievements and challenges, reflects on the lessons learned and outlines future research priorities.

Screening for resistance to four fungal diseases and associated genomic regions in a snap bean diversity panel

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.

Mapping of adult plant recessive resistance to anthracnose in Indian common bean landrace Baspa/KRC 8

BACKGROUND

The common bean (Phaseolus vulgaris) has become the food of choice owing to its wealthy nutritional profile, leading to a considerable increase in its cultivation worldwide. However, anthracnose has been a major impediment to production and productivity, as elite bean cultivars are vulnerable to this disease. To overcome barriers in crop production, scientists worldwide are working towards enhancing the genetic diversity of crops. One way to achieve this is by introducing novel genes from related crops, including landraces like KRC 8. This particular landrace, found in the North Western Himalayan region, has shown adult plant resistance against anthracnose and also possesses a recessive resistance gene.

METHODS AND RESULTS

In this study, a population of 179 F2:9 RIL individuals (Jawala × KRC 8) was evaluated at both phenotypic and genotypic levels using over 830 diverse molecular markers to map the resistance gene present in KRC 8. We have successfully mapped a resistance gene to chromosome Pv01 using four SSR markers, namely IAC 238, IAC 235, IAC 259, and BM 146. The marker IAC 238 is closely linked to the gene with a distance of 0.29 cM, while the other markers flank the recessive resistance gene at 10.87 cM (IAC 259), 17.80 cM (BM 146), and 25.22 cM (IAC 235). Previously, a single recessive anthracnose resistance gene (co-8) has been reported in the common bean accession AB 136. However, when we performed PCR amplification with our tightly linked marker IAC 238, we got different amplicons in AB 136 and KRC 8. Interestingly, the susceptible cultivar Jawala produced the same amplicon as AB 136. This observation indicated that the recessive gene present in KRC 8 is different from co-8. As the gene is located far away from the Co-1 locus, we suggest naming the recessive gene co-Indb/co-19. Fine mapping of co-Indb in KRC 8 may provide new insights into the cloning and characterization of this recessive gene so that it can be incorporated into future bean improvement programs. Further, the tightly linked marker IAC 238 can be utilized in marker assisted introgression in future bean breeding programs.

CONCLUSION

The novel co-Indb gene present in Himalayan landrace KRC 8, showing adult plant resistance against common bean anthracnose, is independent from all the resistance genes previously located on chromosome Pv01.

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5

KEY MESSAGE

Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

Modified screening method of middle american dry bean genotypes reveals new genomic regions on Pv10 associated with anthracnose resistance

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

Genetic mapping of the Andean anthracnose resistance gene present in the common bean cultivar BRSMG Realce

The rajado seeded Andean bean (Phaseolus vulgaris L.) cultivar BRSMG Realce (striped seed coat) developed by Embrapa expressed a high level of anthracnose resistance, caused by Colletotrichum lindemuthianum, in field and greenhouse screenings. The main goal of this study was to evaluate the inheritance of anthracnose resistance in BRSMG Realce, map the resistance locus or major gene cluster previously named as Co-Realce, identify resistance-related positional genes, and analyze potential markers linked to the resistance allele. F₂ plants derived from the cross BRSMG Realce × BRS FC104 (Mesoamerican) and from the cross BRSMG Realce × BRS Notável (Mesoamerican) were inoculated with the C. lindemuthianum races 475 and 81, respectively. The BRSMG Realce × BRS FC104 F₂ population was also genotyped using the DArTseq technology. Crosses between BRSMG Realce and BAT 93 (Mesoamerican) were also conducted and resulting F₂ plants were inoculated with the C. lindemuthianum races 65 and 1609, individually. The results shown that anthracnose resistance in BRSMG Realce is controlled by a single locus with complete dominance. A genetic map including 1,118 SNP markers was built and shown 78% of the markers mapped at a distances less than 5.0 cM, with a total genetic length of 4,473.4 cM. A major locus (Co-Realce) explaining 54.6% of the phenotypic variation of symptoms caused by the race 475 was identified in Pv04, flanked by the markers snp1327 and snp12782 and 4.48 cM apart each other. These SNPs are useful for marker-assisted selection, due to an estimated selection efficiency of 99.2%. The identified resistance allele segregates independently of the resistance allele Co-3³ (Pv04) present in BAT 93. The mapped genomic region with 704,867 bp comprising 63 putative genes, 44 of which were related to the pathogen-host interaction. Based on all these results and evidence, anthracnose resistance in BRSMG Realce should be considered as monogenic, useful for breeding purpose. It is proposed that locus Co-Realce is unique and be provisionally designated as CoPv04^R until be officially nominated in accordance with the rules established by the Bean Improvement Cooperative Genetics Committee.

Phaseolus vulgaris-Colletotrichum lindemuthianum Pathosystem in the Post-Genomic Era: An Update

Phaseolus vulgaris-Colletotrichum lindemuthianum is one among the oldest host and pathogen interface. Researchers have taken painstaking efforts across the world for understanding the dialogue during early and late phases of interaction. Collectively, these efforts resulted in the deluge of information that helped the researchers to underpin the interface. The latest molecular biology techniques furnished novel detection methods for the anthracnose pathogen, refined the understanding of pathogen population dynamics, and provided the insights on co-evolutionary common bean resistance and C. lindemuthianum virulence dynamics. One of the important breakthroughs came when the Phaseolus vulgaris and its corresponding anthracnose pathogen (C. lindemuthianum) genomes were decoded in 2014 and 2017, respectively. Availability of both the genomes yielded a significant genomic information that helped bean communities to fine map the economically important traits and to identify the pathogenicity determinants and effector molecules. The interface is in a continuous development as knowledge of the anthracnose resistance genes, their precise physical locations, and the identification of effector proteins; the fungus arsenals are being routinely updated. Hence, we revisited the interface and tried to provide an overview of host pathogen dialogue in the genomic era. Additionally, we compiled the sporadic information on this pathosystem from India and provided its futuristic road map to shape its research in the world and northern India, the major dry bean area in the country.

Genome-wide association mapping reveals race-specific SNP markers associated with anthracnose resistance in carioca common beans

Brazil is the largest consumer of dry edible beans (Phaseolus vulgaris L.) in the world, 70% of consumption is of the carioca variety. Although the variety has high yield, it is susceptible to several diseases, among them, anthracnose (ANT) can lead to losses of up to 100% of production. The most effective strategy to overcome ANT, a disease caused by the fungus Colletotrichum lindemuthianum, is the development of resistant cultivars. For that reason, the selection of carioca genotypes resistant to multiple ANT races and the identification of loci/markers associated with genetic resistance are extremely important for the genetic breeding process. Using a carioca diversity panel (CDP) with 125 genotypes and genotyped by BeadChip BARCBean6K_3 and a carioca segregating population AM (AND-277 × IAC-Milênio) genotyped by sequencing (GBS). Multiple interval mapping (MIM) and genome-wide association studies (GWAS) were used as mapping tools for the resistance genes to the major ANT physiological races present in the country. In general, 14 single nucleotide polymorphisms (SNPs) showed high significance for resistance by GWAS, and loci associated with multiple races were also identified, as the Co-3 locus. The SNPs ss715642306 and ss715649427 in linkage disequilibrium (LD) at the beginning of chromosome Pv04 were associated with all the races used, and 16 genes known to be related to plant immunity were identified in this region. Using the resistant cultivars and the markers associated with significant quantitative resistance loci (QRL), discriminant analysis of principal components (DAPC) was performed considering the allelic contribution to resistance. Through the DAPC clustering, cultivar sources with high potential for durable anthracnose resistance were recommended. The MIM confirmed the presence of the Co-1⁴locus in the AND-277 cultivar which revealed that it was the only one associated with resistance to ANT race 81. Three other loci were associated with race 81 on chromosomes Pv03, Pv10, and Pv11. This is the first study to identify new resistance loci in the AND-277 cultivar. Finally, the same Co-1⁴locus was also significant for the CDP at the end of Pv01. The new SNPs identified, especially those associated with more than one race, present great potential for use in marker-assisted and early selection of inbred lines.

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

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

Different loci control resistance to different isolates of the same race of Colletotrichum lindemuthianum in common bean

Linkage and genome-wide association analyses using high-throughput SNP genotyping revealed different loci controlling resistance to different isolates of race 65 of Colletotrichum lindemuthianum in common bean. Development of varieties with durable resistance to anthracnose is a major challenge in common bean breeding programs because of the extensive virulence diversity of Colletotrichum lindemuthianum fungus. We used linkage and genome-wide association analyses to tap the genomic regions associated with resistance to different isolates of race 65. Linkage mapping was done using an F₂ population derived from the cross between the Mesoamerican common beans BRS Estilo x Ouro Vermelho, inoculated with two different isolates of race 65. Association genetics relied on a diversity common bean panel containing 189 common bean accessions inoculated with five different isolates of race 65 as an attempt to validate the linkage analysis findings and, eventually, identify other genomic regions associated with resistance to race 65. The F₂ population and diversity panel were genotyped with the BARCBean6K_3 Illumina BeadChip containing 5398 SNP markers. Both linkage and genome-wide association analyses identified different loci controlling resistance to different isolates of race 65 on linkage group Pv04. Genome-wide association analysis also detected loci on Pv05, Pv10 and Pv11 associated with resistance to race 65. These findings indicate that resistance to race 65 can be overcome by the virulence diversity among different isolates of the same race and could lead to the loss of resistance after cultivar release. We identified 25 resistant common bean cultivars to all five isolates of race 65 in the diversity panel. The accessions should be useful to develop cultivars combining different resistance genes that favor durable resistance to anthracnose in common bean.

Comprehensive genomic resources related to domestication and crop improvement traits in Lima bean

Lima bean (Phaseolus lunatus L.), one of the five domesticated Phaseolus bean crops, shows a wide range of ecological adaptations along its distribution range from Mexico to Argentina. These adaptations make it a promising crop for improving food security under predicted scenarios of climate change in Latin America and elsewhere. In this work, we combine long and short read sequencing technologies with a dense genetic map from a biparental population to obtain the chromosome-level genome assembly for Lima bean. Annotation of 28,326 gene models show high diversity among 1917 genes with conserved domains related to disease resistance. Structural comparison across 22,180 orthologs with common bean reveals high genome synteny and five large intrachromosomal rearrangements. Population genomic analyses show that wild Lima bean is organized into six clusters with mostly non-overlapping distributions and that Mesomerican landraces can be further subdivided into three subclusters. RNA-seq data reveal 4275 differentially expressed genes, which can be related to pod dehiscence and seed development. We expect the resources presented here to serve as a solid basis to achieve a comprehensive view of the degree of convergent evolution of Phaseolus species under domestication and provide tools and information for breeding for climate change resiliency.

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

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.

North-Western Himalayan Common Beans: Population Structure and Mapping of Quantitative Anthracnose Resistance Through Genome Wide Association Study

Common bean (Phaseolus vulgaris L.) is an important legume crop of north-western (NW) Himalayan region and the major disease that causes catastrophic loss to the crop is anthracnose, which is caused by Colletotrichum lindemuthianum. The pathogen is highly diverse and most of the commercial cultivars are susceptible to different races prevalent in the region. The lack of information on the genomic regions associated with anthracnose resistance in NW Himalayan common bean population prompted us to dissect Quantitative Resistance Loci (QRLs) against major anthracnose races. In this study, 188 common bean landraces collected from NW region were screened against five important anthracnose races and 113 bean genotypes showed resistance to one or multiple races. Genotyping by sequencing (GBS) was performed on a panel of 192 bean lines (4 controls plus 188 Indian beans) and 22,589 SNPs were obtained that are evenly distributed. Population structure analysis of 192 bean genotypes categorized 188 Indian beans into two major clusters representing Andean and Mesoamerican gene pools with obvious admixtures. Many QRLs associated with anthracnose resistance to Indian C. lindemuthianum virulences (race 3, 87, and 503) are located at Pv04 within the gene models that encode typical resistance gene signatures. The QRLs associated with race 73 are located on Pv08 and overlaps with Co-4 anthracnose resistance gene. A SNP located at distal end of Pv11 in a gene model Phvul.011G202300 which encodes a LRR with a typical NB-ARC domain showed association with race 73 resistance. Common bean genomic regions located at Pv03, Pv09, and Pv11 showed association with resistance to anthracnose race 2047. The present study showed presence of many novel bean genomic regions associated with anthracnose resistance. The presence of Co-4 and Co-2 genes in our material is encouraging for breeding durable anthracnose resistant cultivars for the region.

New Andean source of resistance to anthracnose and angular leaf spot: Fine-mapping of disease-resistance genes in California Dark Red Kidney common bean cultivar

Anthracnose (ANT) and angular leaf spot (ALS) caused by Colletotrichum lindemuthianum and Pseudocercospora griseola, respectively, are devastating diseases of common bean around the world. Therefore, breeders are constantly searching for new genes with broad-spectrum resistance against ANT and ALS. This study aimed to characterize the genetic resistance of California Dark Red Kidney (CDRK) to C. lindemuthianum races 73, 2047, and 3481 and P. griseola race 63-39 through inheritance, allelism testing, and molecular analyses. Genetic analysis of response to ANT and ALS in recombinant inbred lines (RILs) from a CDRK × Yolano cross (CY) showed that the resistance of CDRK cultivar is conferred by a single dominant loci, which we named CoPv01CDRK/PhgPv01CDRK. Allelism tests performed with race 3481showed that the resistance gene in CDRK is independent of the Co-1 and Co-AC. We conducted co-segregation analysis in genotypes of 110 CY RILs and phenotypes of the RILs in response to different races of the ANT and ALS pathogens. The results revealed that CoPv01CDRK and PhgPv01CDRK are coinherited, conferring resistance to all races. Genetic mapping of the CY population placed the CoPv01CDRK/PhgPv01CDRK loci in a 245 Kb genomic region at the end of Pv01. By genotyping 19 RILs from the CY population using three additional markers, we fine-mapped the CoPv01CDRK/PhgPv01CDRK loci to a smaller genomic region of 33 Kb. This 33 Kb region harbors five predicted genes based on the common bean reference genome. These results can be applied in breeding programs to develop bean cultivars with ANT and ALS resistance using marker-assisted selection.