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

Fine mapping and identification of a Fusarium wilt resistance gene FwS1 in pea

KEY MESSAGE

A Fusarium wilt resistance gene FwS1 on pea chromosome 6 was identified and mapped to a 91.4 kb region by a comprehensive genomic-based approach, and the gene Psat6g003960 harboring NB-ARC domain was identified as the putative candidate gene. Pea Fusarium wilt, incited by Fusarium oxysporum f. sp. pisi (Fop), has always been a devastating disease that causes severe yield losses and economic damage in pea-growing regions worldwide. The utilization of pea cultivars carrying resistance gene is the most efficient approach for managing this disease. In order to finely map resistance gene, F2 populations were established through the cross between Shijiadacaiwan 1 (resistant) and Y4 (susceptible). The resistance genetic analysis indicated that the Fop resistance in Shijiadacaiwan 1 was governed by a single dominant gene, named FwS1. Based on the bulked segregant analysis sequencing analyses, the gene FwS1 was initially detected on chromosome 6 (i.e., linking group II, chr6LG2), and subsequent linkage mapping with 589 F2 individuals fine-mapped the gene FwS1 into a 91.4 kb region. The further functional annotation and haplotype analysis confirmed that the gene Psat6g003960, characterized by a NB-ARC (nucleotide-binding adaptor shared by APAF-1, R proteins, and CED-4) domain, was considered as the most promising candidate gene. The encoding amino acids were altered by a "T/C" single-nucleotide polymorphism (SNP) in the first exon of the Psat6g003960, and based on this SNP locus, the molecular marker A016180 was determined to be a diagnostic marker for FwS1 by validating its specificity in both pea accessions and genetic populations with different genetic backgrounds. The FwS1 with diagnostic KASP marker A016180 could facilitate marker-assisted selection in resistance pea breeding in pea. In addition, a comparison of the candidate gene Psat6g003960 in 74SN3B and SJ1 revealed the same sequences. This finding indicated that 74SN3B carried the candidate gene for FwS1, suggesting that FwS1 and Fwf may be closely linked or an identical resistant gene against Fusarium wilt.

Identification of novel sources of partial and incomplete hypersensitive resistance to rust and associated genomic regions in common bean

Common bean (Phaseolus vulgaris) is one of the legume crops most consumed worldwide and bean rust is one of the most severe foliar biotrophic fungal diseases impacting its production. In this work, we searched for new sources of rust resistance (Uromyces appendiculatus) in a representative collection of the Portuguese germplasm, known to have accessions with an admixed genetic background between Mesoamerican and Andean gene pools. We identified six accessions with incomplete hypersensitive resistance and 20 partially resistant accessions of Andean, Mesoamerican, and admixed origin. We detected 11 disease severity-associated single-nucleotide polymorphisms (SNPs) using a genome-wide association approach. Six of the associations were related to partial (incomplete non-hypersensitive) resistance and five to incomplete hypersensitive resistance, and the proportion of variance explained by each association varied from 4.7 to 25.2%. Bean rust severity values ranged from 0.2 to 49.1% and all the infection types were identified, reflecting the diversity of resistance mechanisms deployed by the Portuguese germplasm.The associations with U. appendiculatus partial resistance were located in chromosome Pv08, and with incomplete hypersensitive resistance in chromosomes Pv06, Pv07, and Pv08, suggesting an oligogenic inheritance of both types of resistance. A resolution to the gene level was achieved for eight of the associations. The candidate genes proposed included several resistance-associated enzymes, namely β-amylase 7, acyl-CoA thioesterase, protein kinase, and aspartyl protease. Both SNPs and candidate genes here identified constitute promising genomics targets to develop functional molecular tools to support bean rust resistance precision breeding.

Genome-Wide Association Study and Genomic Prediction of Fusarium Wilt Resistance in Common Bean Core Collection

The common bean (Phaseolus vulgaris L.) is a globally cultivated leguminous crop. Fusarium wilt (FW), caused by Fusarium oxysporum f. sp. phaseoli (Fop), is a significant disease leading to substantial yield loss in common beans. Disease-resistant cultivars are recommended to counteract this. The objective of this investigation was to identify single nucleotide polymorphism (SNP) markers associated with FW resistance and to pinpoint potential resistant common bean accessions within a core collection, utilizing a panel of 157 accessions through the Genome-wide association study (GWAS) approach with TASSEL 5 and GAPIT 3. Phenotypes for Fop race 1 and race 4 were matched with genotypic data from 4740 SNPs of BARCBean6K_3 Infinium Bea Chips. After ranking the 157-accession panel and revealing 21 Fusarium wilt-resistant accessions, the GWAS pinpointed 16 SNPs on chromosomes Pv04, Pv05, Pv07, Pv8, and Pv09 linked to Fop race 1 resistance, 23 SNPs on chromosomes Pv03, Pv04, Pv05, Pv07, Pv09, Pv10, and Pv11 associated with Fop race 4 resistance, and 7 SNPs on chromosomes Pv04 and Pv09 correlated with both Fop race 1 and race 4 resistances. Furthermore, within a 30 kb flanking region of these associated SNPs, a total of 17 candidate genes were identified. Some of these genes were annotated as classical disease resistance protein/enzymes, including NB-ARC domain proteins, Leucine-rich repeat protein kinase family proteins, zinc finger family proteins, P-loopcontaining nucleoside triphosphate hydrolase superfamily, etc. Genomic prediction (GP) accuracy for Fop race resistances ranged from 0.26 to 0.55. This study advanced common bean genetic enhancement through marker-assisted selection (MAS) and genomic selection (GS) strategies, paving the way for improved Fop resistance.

Genetic variation in a tepary bean (Phaseolus acutifolius A. Gray) diversity panel reveals loci associated with biotic stress resistance

Tepary bean (Phaseolus acutifolius A. Gray), indigenous to the arid climates of northern Mexico and the Southwest United States, diverged from common bean (Phaseolus vulgaris L.), approximately 2 million years ago and exhibits a wide range of resistance to biotic stressors. The tepary genome is highly syntenic to the common bean genome providing a foundation for discovery and breeding of agronomic traits between these two crop species. Although a limited number of adaptive traits from tepary bean have been introgressed into common bean, hybridization barriers between these two species required the development of bridging lines to alleviate this barrier. Thus, to fully utilize the extant tepary bean germplasm as both a crop and as a donor of adaptive traits, we developed a diversity panel of 422 cultivated, weedy, and wild tepary bean accessions which were then genotyped and phenotyped to enable population genetic analyses and genome-wide association studies for their response to a range of biotic stressors. Population structure analyses of the panel revealed eight subpopulations and the differentiation of botanical varieties within P. acutifolius. Genome-wide association studies revealed loci and candidate genes underlying biotic stress resistance including quantitative trait loci for resistance to weevils, common bacterial blight, Fusarium wilt, and bean common mosaic necrosis virus that can be harnessed not only for tepary bean but also common bean improvement.

Genetic architecture of inducible and constitutive metabolic profile related to drought resistance in qingke (Tibetan hulless barley)

Qingke (Tibetan hulless barley, Hordeum vulgare L. var. nudum) is the primary food crop on the Tibet Plateau, the long-term drought and other harsh environments makes qingke an important resource for the study of abiotic resistance. Here, we evaluated the drought sensitivity of 246 qingke varieties. Genome-wide association studies (GWAS) found that root-specific expressed gene CYP84 may be involved in the regulation of drought resistance. Based on widely targeted metabolic profiling, we identified 2,769 metabolites in qingke leaves, of which 302 were significantly changed in response to drought stress, including 4-aminobutyric acid (GABA), proline, sucrose and raffinose. Unexpectedly, these drought-induced metabolites changed more violently in drought-sensitive qingkes, while the constitutive metabolites that had little response to drought stress, such as C-glycosylflavonoids and some amino acids, accumulated excessively in drought-resistant qingkes. Combined with metabolite-based genome-wide association study (mGWAS), a total of 1,006 metabolites under optimal condition and 1,031 metabolites under mild drought stress had significant associated loci. As a marker metabolite induced by drought stress, raffinose was significantly associated with two conservatively adjacent α-galactosidase genes, qRT-PCR suggests that these two genes may jointly regulate the raffinose content in qingke. Besides, as constituent metabolites with stable differences between drought-sensitive and drought-resistant qingkes, a class of C-glycosylflavonoids are simultaneously regulated by a UDP-glucosyltransferase gene. Overall, we performed GWAS for sensitivity and widely targeted metabolites during drought stress in qingke for the first time, which provides new insights into the response mechanism of plant drought stress and drought resistance breeding.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Associated SNPs, Heritabilities, Trait Correlations, and Genomic Breeding Values for Resistance in Snap Beans (Phaseolus vulgaris L.) to Root Rot Caused by Fusarium solani (Mart.) f. sp. phaseoli (Burkholder)

Root rot is a major constraint to snap bean (Phaseolus vulgaris) production in the United States and around the world. Genetic resistance is needed to effectively control root rot disease because cultural control methods are ineffective, and the pathogen will be present at the end of one season of production on previously clean land. A diversity panel of 149 snap bean pure lines was evaluated for resistance to Fusarium root rot in Oregon. Morphological traits potentially associated with root rot resistance, such as aboveground biomass, adventitious roots, taproot diameter, basal root diameter, deepest root angle, shallowest root angle, root angle average, root angle difference, and root angle geometric mean were evaluated and correlated to disease severity. A genome wide association study (GWAS) using the Fixed and random model Circulating Probability Unification (FarmCPU) statistical method, identified five associated single nucleotide polymorphisms (SNPs) for disease severity and two SNPs for biomass. The SNPs were found on Pv03, Pv07, Pv08, Pv10, and Pv11. One candidate gene for disease reaction near a SNP on Pv03 codes for a peroxidase, and two candidates associated with biomass SNPs were a 2-alkenal reductase gene cluster on Pv10 and a Pentatricopeptide repeat domain on Pv11. Bean lines utilized in the study were ranked by genomic estimated breeding values (GEBV) for disease severity, biomass, and the root architecture traits, and the observed and predicted values had high to moderate correlations. Cross validation of genomic predictions showed slightly lower correlational accuracy. Bean lines with the highest GEBV were among the most resistant, but did not necessarily rank at the very top numerically. This study provides information on the relationship of root architecture traits to root rot disease reaction. Snap bean lines with genetic merit for genomic selection were identified and may be utilized in future breeding efforts.

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

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