A Cluster of Peronospora parasitica 13-like (NBS-LRR) Genes Is Associated with Powdery Mildew (Erysiphe polygoni) Resistance in Mungbean (Vigna radiata).

Powdery mildew (PM) caused by Erysiphe polygoni is an important foliar disease in mungbean (Vigna radiata). A previous study showed that QTL qPMRUM5-2 is a major locus for PM resistance in mungbean accession RUM5 (highly resistant). Bioinformatics analysis revealed that flanking markers of the qPMRUM5-2 covered a region of 1.93 Mb. In this study, we conducted fine mapping for the qPMRUM5-2 using the F2 population of 1156 plants of the cross between Chai Nat 60 (CN60; highly susceptible) and RUM5. PM resistance evaluation was performed under field conditions using F2:3 lines grown in three different environments. QTL analyses consistently located the qPMRUM5-2 to a 0.09 cm interval on linkage group 6 between InDel markers VrLG6-InDel05 and VrLG6-InDel10, which corresponded to a 135.0 kb region on chromosome 8 containing nine predicted genes of which five were NBS-LRR-type genes Recognition of Peronospora parasitica 13-like protein (RPP13L). Whole-genome re-sequencing of RUM5 and CN60 showed polymorphisms in four RPP13L genes predictively cause substantial amino acid changes, rendering them important candidate genes for PM resistance. The InDel markers VrLG6-InDel05 and VrLG6-InDel10 flanking to the qPMRUM5-2 would be useful for marker-assisted breeding of PM resistance in the mungbean.

Fine mapping of QTL conferring resistance to calcareous soil in mungbean reveals VrYSL3 as candidate gene for the resistance.

Iron is a crucial nutrient for biological functions in plants. High-pH and calcareous soil is a major stress causing iron deficiency chlorosis (IDC) symptoms and yield losses in crops. Use of calcareous soil-tolerance genetic resources is the most effective preventative method to combat the effects of high-pH and calcareous soils. A previous study using a mungbean recombinant inbred line (RIL) population of the cross Kamphaeg Saen 2 (KPS2; IDC susceptible) × NM-10-12 identified a major quantitative trait locus (QTL), qIDC3.1, which controls resistance and explains more than 40% of IDC variation. In this study, we fine-mapped qIDC3.1 and identified an underlying candidate gene. A genome wide association analysis (GWAS) using 162 mungbean accessions identified single nucleotide polymorphisms (SNPs) on chromosome 6; several SNPs were associated with soil plant analysis development (SPAD) values and IDC visual scores of mungbeans planted on calcareous soil, respectively. These SNPs corresponded to qIDC3.1. Using the same RIL population as in the previous study and an advanced backcross population developed from KPS2 and IDC-resistant inbred line RIL82, qIDC3.1 was further confirmed and fine-mapped to an interval of 217 kilobases harboring five predicted genes, including LOC106764181 (VrYSL3), which encodes a yellow stripe1-like-3 (YSL3) protein, YSL3 is involved in iron deficiency resistance. Gene expression analysis revealed that VrYSL3 was highly expressed in mungbean roots. In calcareous soil, expression of VrYSL3 was significantly up-regulated, and it was more obviously upregulated in the roots of RIL82, than in those of KPS2. Sequence comparison of VrYSL3 between the RIL82 and KPS2 revealed four SNPs that result in amino acid changes in the VrYSL3 protein and a 20-bp insertion/deletion in the promoter where a cis-regulatory element resides. Transgenic Arabidopsis thaliana plants overexpressing VrYSL3 showed enhanced iron and zinc contents in the leaves. Taken together, these results indicate that VrYSL3 is a strong candidate gene responsible for calcareous soil resistance in mungbean.

Tandemly duplicated genes encoding polygalacturonase inhibitors are associated with bruchid (Callosobruchus chinensis) resistance in moth bean (Vigna aconitifolia).

Bruchids are stored-grain insect pests responsible for serious seed loss in legume crops. A previous study using an F2 population (F2OA) derived from a cross between wild moth-bean (Vigna aconitifolia [Jacq.] Maréchal) accession TN67 (resistant) and cultivated moth-bean accession ICPMO056 (susceptible) revealed that resistance to the azuki bean weevil (Callosobruchus chinensis L.) in TN67 was regulated by a single gene located in the major quantitative trait locus-qVacBrc2.1. In this study, qVacBrc2.1 was finely mapped and candidate genes in this locus were identified using F2OA and another large F2 population (F2NB) derived from the cross mentioned previously. In contrast to the previous study, segregation analysis in the F2NB population revealed that resistance against this pest was controlled by two genes. Furthermore, the addition of novel markers to qVacBrc2.1 and reanalysis of the QTL in the F2OA population demonstrated that qVacBrc2.1 constituted two linked QTLs-qVacBrc2.1-A and qVacBrc2.1-B. The presence of qVacBrc2.1-B was verified using the population F2NB. Comparative genomics using three Vigna spp. strongly suggested the presence of two tandemly duplicated genes, VacPGIP1 and VacPGIP2, which encoded polygalacturonase inhibitors (polygalacturonase-inhibiting proteins) as the candidates for conferring resistance, but only VacPGIP1 could be successfully cloned and sequenced. The alignment of VacPGIP1 coding sequences of TN67 and ICPMO056 revealed eight single nucleotide polymorphisms, three of which altered the amino-acid sequence of the predicted domains of polygalacturonase inhibitors in ICPMO056. Overall, these findings indicate that VacPGIP1 and VacPGIP2 regulated C. chinensis resistance in TN67.

A Class II KNOX Gene, KNAT7-1, Regulates Physical Seed Dormancy in Mungbean [Vigna radiata (L.) Wilczek].

Seed dormancy in wild mungbean (Vigna radiata var. sublobata) may be useful for the breeding of cultivated mungbean (var. radiata) with pre-harvest sprouting resistance. Previous studies have identified two major quantitative trait loci (QTLs) for seed dormancy, HsA and Sdwa5.1.1+, in wild mungbean that are possibly having the same locus or linked. However, these QTLs have not been confirmed/verified and a molecular basis of seed dormancy in mungbean is not yet known. In this study, we aimed to finely map the Sdwa5.1.1+ and identify candidate gene(s) for this locus. Microscopic observations revealed that wild mungbean "ACC41" seeds had a palisade cuticle layer, while cultivated mungbean "Kamphaeng Saen 2" (KPS2) seeds lacked this layer. Fine mapping using an F2 population developed from a cross between ACC41 and KPS2 revealed two linked QTLs, Sdwa5.1.1+ and Sdwa5.1.2+, controlling seed dormancy. The Sdwa5.1.1+ was confirmed in an F2:3 population derived from the same cross and mapped to a 3.298-Kb region containing only one gene LOC106767068, designated as VrKNAT7-1, which encodes the transcription factor KNOTTED ARABIDOPSIS THALIANA7 (KNAT7), a class II KNOTTED1-LIKE HOMEOBOX (KNOX II) protein. VrKNAX7 sequence alignment between ACC41 and KPS2 revealed several polymorphisms in the coding, untranslated, and promoter regions. Quantitative real-time PCR (qRT-PCR) analysis revealed that the expression of VrKNAT7-1 and VrCYP86A, a putative downstream regulation of VrKNAT7-1, in the seed coat of ACC41 is statistically much higher than that of KPS2. Altogether, these results indicate that VrKNAT7-1 controls physical seed dormancy in the wild mungbean ACC41.

Mapping QTLs Controlling Soybean Rust Disease Resistance in Chiang Mai 5, an Induced Mutant Cultivar.

Soybean rust (SBR) caused by the fungus Phakopsora pachyrhizi is an important folia disease of soybean (Glycine max). In this study, we identified QTLs controlling SBR in Chiang Mai 5 (CM5), an SBR-resistant cultivar developed by induced mutation breeding. A recombinant inbred line (RIL) population of 108 lines developed from a cross between Sukhothai 2 (SKT2, a susceptible cultivar) and CM5 was evaluated for SBR resistance under field conditions in Thailand. QTL analysis for the resistance in the RIL population identified a single QTL, qSBR18.1, for resistance. qSBR18.1 was mapped to a 212-kb region on chromosome 18 between simple sequence repeat markers Satt288 and sc21_3420 and accounted for 21.31-35.09% depending on the traits evaluated for resistance. The qSBR18.1 interval overlapped with genomic regions containing resistance to P. pachyrhizi 4 (Rpp4), a locus for SBR resistance. Three tightly linked genes, Glyma.18G226250, Glyma.18G226300, and Glyma.18G226500, each encoding leucine-rich repeat-containing protein, were identified as candidate genes for SBR resistance at the qSRB18.1. The qSBR18.1 would be useful for breeding of SBR resistance.

QTL Mapping for Agronomic and Adaptive Traits Confirmed Pleiotropic Effect of mog Gene in Black Gram [Vigna mungo (L.) Hepper].

Organ size and architecture of plants are important traits affecting crop yield and agronomic practices. An induced mutant, multiple-organ gigantism (MOG), of black gram (Vigna mungo) has been obtained, which shows gigantic leaves, fruit, seed, and architecture (plant height) but lower number of pods per plant. These traits are a pleiotropic effect of a single recessive gene, mog. In this study, we investigated variation of 16 agronomic and adaptive traits in a recombinant inbred line (RIL) population derived from a cross between the MOG mutant (V. mungo var. mungo) and wild black gram (V. mungo var. silvestris) accession TC2210 and identified quantitative trait loci (QTLs) controlling those traits to gain a better understanding of the effect of the mog gene on breeding. The results showed that most of the traits (100-seed weight, leaf size, and plant height) showed moderate narrow-sense heritability (h 2) (45-65%), while pod size and seed length (SDL) showed high h 2 (>75%) and pod dehiscence (shattering), and seed width (SDW) and days to flowering showed low h 2 (<35%). The QTLs for the traits were mapped onto a high-density linkage map developed for the RIL population. Inclusive composite interval mapping identified 42 QTLs in total for the 16 traits with number of QTLs per trait ranging from one to six. Major QTLs for the MOG phenotypes were clustered on linkage group (LG) 6, confirming the pleiotropic effect of the mog gene. Effect of the mog gene/QTL for the MOG phenotypic variations was not high, ranging from about 15% in plant height to 40% in leaf size. For 100-seed weight, which is the most interesting trait, the mog gene/QTL contributed about 30% of the total trait variation and showed an additive effect of only 0.51 g, which is only about 1.5-fold higher than that of the other five QTLs detected for this trait. These results indicated that mog gene expression is highly affected by environment and the effect of the gene toward organ size and plant height is not extraordinarily high. Implications of the findings of this study and exploiting of the MOG mutant in breeding were also discussed.

Development of an SNP-based high-density linkage map and QTL analysis for bruchid (Callosobruchus maculatus F.) resistance in black gram (Vigna mungo (L.) Hepper).

Black gram (Vigna mungo var. mungo) is an important pulse crop in Asia. The cowpea weevil (Callosobruchus maculatus) is a stored-seed insect pest (seed weevil/bruchid) that causes serious postharvest losses in pulse crops, including black gram. In this study, we constructed a high-density linkage map for black gram and identified quantitative trait loci (QTLs) for C. maculatus resistance. A recombinant inbred line (RIL) population of 150 lines from a cross between BC48 [cultivated black gram (var. mungo); bruchid-susceptible] and TC2210 [wild black gram (var. silvestris); bruchid-resistant] were used to construct a linkage map of 3,675 SNP markers from specific-locus amplified fragment sequencing. The map comprised 11 linkage groups spanning 1,588.7 cM with an average distance between adjacent markers of 0.57 cM. Seeds of the RIL population grown in 2016 and 2017 were evaluated for C. maculatus resistance through two traits; the percentage of damaged seeds (PDS) and infestation severity progress (AUDPS). Inclusive composite interval mapping identified three QTLs each for PDS and AUDPS. Two QTLs, qVmunBr6.1 and qVmunBr6.2, mapped about 10 cM apart on linkage group 6 were common between PDS and AUDPS. Comparative genome analysis revealed that qVmunBr6.1 and qVmunBr6.2 are new loci for C. maculatus resistance in Vigna species and that genes encoding a lectin receptor kinase and chitinase are candidates for qVmunBr6.2. The high-density linkage map constructed and QTLs for bruchid resistance identified in this study will be useful for molecular breeding of black gram.

Novel Alleles of Two Tightly Linked Genes Encoding Polygalacturonase-Inhibiting Proteins (VrPGIP1 and VrPGIP2) Associated with the Br Locus That Confer Bruchid (Callosobruchus spp.) Resistance to Mungbean (Vigna radiata) Accession V2709.

Nearly all mungbean cultivars are completely susceptible to seed bruchids (Callosobruchus chinensis and Callosobruchus maculatus). Breeding bruchid-resistant mungbean is a major goal in mungbean breeding programs. Recently, we demonstrated in mungbean (Vigna radiata) accession V2802 that VrPGIP2, which encodes a polygalacturonase inhibiting protein (PGIP), is the Br locus responsible for resistance to C. chinensis and C. maculatus. In this study, mapping in mungbean accession V2709 using a BC11F2 population of 355 individuals revealed that a single major quantitative trait locus, which controlled resistance to both C. chinensis and C. maculatus, was located in a 237.35 Kb region of mungbean chromosome 5 that contained eight annotated genes, including VrPGIP1 (LOC106760236) and VrPGIP2 (LOC106760237). VrPGIP1 and VrPGIP2 are located next to each other and are only 27.56 Kb apart. Sequencing VrPGIP1 and VrPGIP2 in "V2709" revealed new alleles for both VrPGIP1 and VrPGIP2, named VrPGIP1-1 and VrPGIP2-2, respectively. VrPGIP2-2 has one single nucleotide polymorphism (SNP) at position 554 of wild type VrPGIP2. This SNP is a guanine to cystine substitution and causes a proline to arginine change at residue 185 in the VrPGIP2 of "V2709". VrPGIP1-1 has 43 SNPs compared with wild type and "V2802", and 20 cause amino acid changes in VrPGIP1. One change is threonine to proline at residue 185 in VrPGIP1, which is the same as in VrPGIP2. Sequence alignments of VrPGIP2 and VrPGIP1 from "V2709" with common bean (Phaseolus vulgaris) PGIP2 revealed that residue 185 in VrPGIP2 and VrPGIP1 contributes to the secondary structures of proteins that affect interactions between PGIP and polygalacturonase, and that some amino acid changes in VrPGIP1 also affect interactions between PGIP and polygalacturonase. Thus, tightly linked VrPGIP1 and VrPGIP2 are the likely genes at the Br locus that confer bruchid resistance in mungbean "V2709".

Genetic diversity and structure of the zombi pea (Vigna vexillata (L.) A. Rich) gene pool based on SSR marker analysis.

Zombi pea (Vigna vexillata (L.) A. Rich) is an underutilized legume species and a useful gene source for resistance to biotic and abiotic stresses, although there is little understanding on its genetic diversity and structure. In this study, 422 (408 wild and 14 cultivated) accessions of zombi pea from diverse origins (201 from Africa, 126 from America, 85 from Australia, 5 from Asia and 5 from unknown origin) were analyzed with 20 simple sequence repeat (SSR) markers to determine its genetic diversity and genetic structure. The SSR markers detected 273 alleles in total with a mean of 13.6 alleles per locus. Polymorphism information content values of the markers varied from 0.58 to 0.90 with an average of 0.76. Overall gene diversity was 0.715. Gene diversity and average allelic richness was highest in Africa (0.749 and 8.08, respectively) and lowest in America (0.435 and 4.10, respectively). Nei's genetic distance analysis revealed that the highest distance was between wild Australia and cultivated Africa (0.559), followed by wild West Africa and wild Australia (0.415). STRUCTURE, neighbor-joining (NJ), and principal coordinate analyses consistently showed that these zombi pea accessions were clustered into three major groups, viz. America, Africa and Asia, and Australia. NJ tree also suggested that American and Australian accessions are originated from East African zombi peas, and that the cultivated accessions from Africa and Asia were genetically distinct, while those from America were clustered with some cultivated accessions from Africa. These results suggest that Africa is the center of origin and diversity of zombi pea, and that domestication of this pea took place more than once in different regions.

A gene encoding a polygalacturonase-inhibiting protein (PGIP) is a candidate gene for bruchid (Coleoptera: bruchidae) resistance in mungbean (Vigna radiata).

KEY MESSAGE: The Br locus confers bruchid resistance in mungbean; VrPGIP2 (encoding a polygalacturonase inhibitor) is a strong candidate gene for this resistance. The VrPGIP2 sequence differs between resistant and susceptible lines. Azuki bean weevil (Callosobruchus chinensis) and cowpea weevil (Callosobruchus maculatus) are serious insect pests of mungbean during storage. Bruchid resistance in mungbean is controlled by a single dominant locus, Br. Although the Br locus has been located on a genetic map, molecular basis and function of the gene remain unknown. In this study, high-resolution mapping using a BC11F2 population of 418 plants derived from a cross between 'Kamphaeng Saen 1' (KPS1; susceptible) and 'V2802' (resistant) using simple sequence repeat (SSR) markers delimited the Br locus to a genomic region of 38 Kb of chromosome 5 containing two annotated genes. EST-SSR marker DMB-SSR158 co-segregated perfectly with the Br locus. Bioinformatics analyses revealed that DMB-SSR158 corresponds to a gene encoding a polygalacturonase inhibitor (polygalacturonase-inhibiting protein PGIP) and was designated as VrPGIP2. Comparison of VrPGIP2 coding sequences between four bruchid-resistant (V2802, V1128, V2817 and TC1966) and four bruchid-susceptible (KPS1, Sulu-1, CM and an unknown accession) mungbean lines revealed six single nucleotide polymorphisms (SNPs) between the resistant and susceptible groups. Three of the six SNPs resulted in amino acid changes; namely, alanine (A) to serine (S) at position 320, leucine (L) to proline (P) at position 332, and threonine (T) to P at position 335 of the VrPGIP2 sequence in resistant lines, compared with that in susceptible lines. The A to S change at position 320 may affect the interaction between PGIP and polygalacuronase. These results indicate that VrPGIP2 is very likely the gene at the Br locus responsible for bruchid resistance in mungbean.