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

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.

Genome-wide association analysis of anthracnose resistance in the Yellow Bean Collection of Common Bean

Anthracnose caused by Colletotrichum lindemuthianum is a major disease of common bean (Phaseolus vulgaris) worldwide. Yellow beans are a major market class of common bean especially in eastern and southern Africa. The Yellow Bean Collection (YBC), which is comprised of 255 genotypes, and has not been used previously in genetic studies on anthracnose, is an excellent genetic resource for understanding the extent of anthracnose resistance and its genetic architecture in the yellow bean market class. The objectives of this study were i) evaluate the YBC for resistance to races 5, 19, 39, 51, 81, 183, 1050 and 1105 of C. lindemuthianum. and ii) conduct genome-wide association analysis to identify genomic regions and candidate genes associated with resistance to C. lindemuthianum. The YBC was genotyped with 72,866 SNPs, and genome-wide association analysis was conducted using Mixed Linear Model in TASSEL. Andean and Middle American genotypes with superior levels of resistance to the eight races were identified. YBC278 was the only one among 255 genotypes that was highly resistant to all eight races. Resistance to anthracnose in the YBC was controlled by major-effect loci on chromosomes Pv01, Pv03, Pv04, Pv05 and Pv07. The genomic region on Pv01, which overlapped with the Andean locus Co-1 provided resistance to races 81, 1050 and 1105. Significant SNPs for resistance to race 39 were identified on Pv02. The genomic region on Pv04, which overlapped with known major-effect loci Co-3, Co-15, Co-16, Co-y and Co-z, provided resistance to races 5, 19, 51 and 183. Novel genomic regions for resistance to race 39 were identified on Pv05 and Pv07. Plant resistance genes (R genes) with NB-ARC and LRR domains, which occurred in clusters, were identified as positional candidate genes for genomic regions on Pv02 and Pv04.

Genome-wide association study as a powerful tool for dissecting competitive traits in legumes

Legumes are extremely valuable because of their high protein content and several other nutritional components. The major challenge lies in maintaining the quantity and quality of protein and other nutritional compounds in view of climate change conditions. The global need for plant-based proteins has increased the demand for seeds with a high protein content that includes essential amino acids. Genome-wide association studies (GWAS) have evolved as a standard approach in agricultural genetics for examining such intricate characters. Recent development in machine learning methods shows promising applications for dimensionality reduction, which is a major challenge in GWAS. With the advancement in biotechnology, sequencing, and bioinformatics tools, estimation of linkage disequilibrium (LD) based associations between a genome-wide collection of single-nucleotide polymorphisms (SNPs) and desired phenotypic traits has become accessible. The markers from GWAS could be utilized for genomic selection (GS) to predict superior lines by calculating genomic estimated breeding values (GEBVs). For prediction accuracy, an assortment of statistical models could be utilized, such as ridge regression best linear unbiased prediction (rrBLUP), genomic best linear unbiased predictor (gBLUP), Bayesian, and random forest (RF). Both naturally diverse germplasm panels and family-based breeding populations can be used for association mapping based on the nature of the breeding system (inbred or outbred) in the plant species. MAGIC, MCILs, RIAILs, NAM, and ROAM are being used for association mapping in several crops. Several modifications of NAM, such as doubled haploid NAM (DH-NAM), backcross NAM (BC-NAM), and advanced backcross NAM (AB-NAM), have also been used in crops like rice, wheat, maize, barley mustard, etc. for reliable marker-trait associations (MTAs), phenotyping accuracy is equally important as genotyping. Highthroughput genotyping, phenomics, and computational techniques have advanced during the past few years, making it possible to explore such enormous datasets. Each population has unique virtues and flaws at the genomics and phenomics levels, which will be covered in more detail in this review study. The current investigation includes utilizing elite breeding lines as association mapping population, optimizing the choice of GWAS selection, population size, and hurdles in phenotyping, and statistical methods which will analyze competitive traits in legume breeding.

Phylogeny and pathogenicity of Colletotrichum lindemuthianum causing anthracnose of Phaseolus vulgaris cv. Bhaderwah-Rajmash from northern Himalayas, India

With an annual loss of up to 100%, anthracnose caused by Colletotrichum is one of the most devastating diseases of common beans (Phaseolus vulgaris L.). Due to few distinctive morphological characters, Colletotrichum species are frequently misidentified. In India, several Colletotrichum species have been reported as pathogens of Phaseolus species, but none had previously been validated by means of molecular tools. In this study, we studied Colletotrichum strains from common beans cv. Bhaderwah-Rajmash from the northern Himalayas of India based on both morphological and DNA sequence data of six loci, namely ITS, gapdh, chs-1, his3, act, tub2. The strains were identified as C. lindemuthianum that belongs to the C. orbiculare species complex. Representative C. lindemuthianum strains tested on Phaseolus vulgaris cv. Bhaderwah-Rajmash were pathogenic and exhibited variation in symptomology and disease progression. By identifying the causal agent, we provided substantial information to develop the best control strategies for anthracnose of Phaseolus vulgaris from the northern Himalayas of India.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03216-0.

Morphological and molecular characterization of variation in common bean (Phaseolus vulgaris L.) germplasm from Azad Jammu and Kashmir, Pakistan

Phaseolus vulgaris, an essential food and source of protein, is cultivated across the world. This study was carried out to investigate the diversity and population structure of 34 P. vulgaris landrace accessions collected from the Azad Jammu and Kashmir (AJ&K) regions of Pakistan. The samples were analyzed both morphologically and using genetic variation identified through RNA sequencing. Our results indicated that most genetic variation occurs among local accessions, with little genetic variation occurring between geographical regions. In addition, the accessions fell into two major genetic groups. Morphological analysis revealed that these two genetic groups differ in a number of quantitative traits, including seed length, seed width, and seed weight. One accession, DUD-11, appears to be a mixture of the two major groups genetically as well as morphologically. Among the other accessions, DUD-8, RWK-2, and NGD-1 depicted particularly high seed weight along with higher seed length, seed width, and seed yield per plant. We suggest focusing on these accessions in future breeding programs. More generally, our results provide baseline data that will be useful for crop improvement and effective cultivation practices in Pakistan.

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.

Optimizing the Agrobacterium tumifaciens mediated transformation conditions in Colletotrichum lindemuthianum: A step forward to unravel the functions of pathogenicity arsenals

Colletotrichum lindemuthianum is a hemibiotrophic fungal pathogen that causes bean anthracnose and it is rated among the top 10 important diseases infecting beans. Currently our knowledge on molecular mechanisms underlying C. lindemuthianum pathogenesis is limited. About five pathogenicity genes have been identified in C. lindemuthianum using Restricted Enzyme Mediated Integration (REMI) and the transformation using Agroinfection has not been optimized. In this study, a series of experiments were conducted to optimize the key parameters affecting the Agrobacterium tumefaciens- mediated transformation (ATMT) for C. lindemuthianum. The transformation efficiency increased with increase in spore concentration and co-cultivation time. However, the optimum conditions that yielded significant number of transformants were 10⁶ ml^-1 spore concentration, co-cultivation time of 72 h, incubation at 25ºC and using a cellulose membrane filter for the co-cultivation. The optimized protocol resulted in establishment of large mutant library (2400). A few mutants were melanin deficient and a few were unable to produce conidia. To determine the altered pathogenicity, two new approaches such as detached leaf and twig techniques proved reliable and require fewer resources to screen the large mutant libraries in a short time. Among the 1200 transformants tested for virulence, 90% transformants were pathogenically similar to wild type (race 2047), 96 and 24 were reduced and impaired, respectively. The altered avirulent transformants can prove vital for understanding the missing link between growth and developmental stages of pathogen with virulence. This platform will help to develop strategies to determine the potential pathogenicity genes and to decipher molecular mechanisms of host-pathogen interactions in more detail.

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.

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

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