Narrowing down a major QTL region reveals Phytochrome E (PHYE) as the candidate gene controlling flowering time in mungbean (Vigna radiata).

Flowering time is an important agronomic trait that is highly correlated with plant height, maturity time and yield in mungbean. Up to present, however, molecular basis of flowering time in mungbean is poorly understood. Previous studies demonstrated that flowering time in mungbean is largely controlled by a major QTL on linkage group 2 (LG2). In this study, the QTL on the LG2 in mungbean was investigated using F2 and F2:3 populations derived from a cross between mungbean cultivar Kamphaeng Saen 2 (KPS2) and wild mungbean accession ACC41. The QTL was narrowed down to a genome region of 164.87 Kb containing a phytochrome gene, designated VrPHYE, encoding phytochrome E (phyE), a known photoreceptor modulating flowering time. Compared to VrPHYE of the wild ACC41, VrPHYE of KPS2 contained several single nucleotide polymorphisms (SNPs) causing amino acid changes. Those SNPs were also found in other mungbean cultivars. Some amino acid changes were predicted to occur in the regulatory region of phytochromes. Gene expression analysis revealed that VrPHYE in KPS2 was expressed significantly higher than that in ACC41. These results showed that VrPHYE is the candidate gene controlling flowering time in the mungbean.

Two genes encoding caffeoyl coenzyme A O-methyltransferase 1 (CCoAOMT1) are candidate genes for physical seed dormancy in cowpea (Vigna unguiculata (L.) Walp.).

KEY MESSAGE: The major QTL Sdp1.1+ controlling seed dormancy in cowpea was finely mapped, and two CCoAOMT1 genes were identified as candidate genes for the dormancy. Seed dormancy in wild cowpea may be useful in breeding cultivated cowpea with pre-harvest sprouting resistance. A previous study identified a major quantitative trait locus (QTL) for seed dormancy, Sdp1.1+ , using the population of the cross between cultivated cowpea 'JP81610' and wild cowpea 'JP89083.' However, the molecular basis of seed dormancy in cowpea is not yet known. In this study, we aimed to finely map the locus Sdp1.1+ and identify candidate gene(s) for it. Germination tests demonstrated that the seed coat is the major factor controlling seed dormancy in the wild cowpea JP89083. Microscopic observations revealed that wild cowpea seeds, unlike cultivated cowpea seeds, possessed a palisade cuticle layer. Fine mapping using a large F2 population of the cross JP81610 × JP89083 grown in Thailand revealed a single QTL, Sdp1.1+ , controlling seed dormancy. The Sdp1.1+ was confirmed using a small F2 population of the same cross grown in Japan. The Sdp1.1+ was mapped to a 37.34-Kb region containing three genes. Two closely linked genes, Vigun03g278900 (VuCCoAOMT1a) and Vigun03g290000 (VuCCoAOMT1b), located 4.844 Kb apart were considered as candidate genes for seed dormancy. The two genes encoded caffeoyl coenzyme A O-methyltransferase 1 (CCoAOMT1). DNA sequencing and alignment of VuCCoAOMT1a and VuCCoAOMT1b between JP89083 and JP81610 revealed a single nucleotide polymorphism (SNP) causing an amino acid change in VuCCoAOMT1a and several SNPs leading to six amino acid changes in VuCCoAOMT1b. Altogether, these results indicate that VuCCoAOMT1a and VuCCoAOMT1b are candidate genes controlling physical seed dormancy in the wild cowpea JP89083.

A Gene Encoding Xylanase Inhibitor Is a Candidate Gene for Bruchid (Callosobruchus spp.) Resistance in Zombi Pea (Vigna vexillata (L.) A. Rich).

Two bruchid species, Callosobruchus maculatus and Callosobruchus chinensis, are the most significant stored insect pests of tropical legume crops. Previously, we identified a major QTL, qBr6.1, controlling seed resistance to these bruchids in the cultivated zombi pea (Vigna vexillata) accession 'TVNu 240'. In this study, we have narrowed down the qBr6.1 region and identified a candidate gene conferring this resistance. Fine mapping using F2 and F2:3 populations derived from a cross between TVNu 240 and TVNu 1623 (susceptible) revealed the existence of two tightly linked QTLs, designated qBr6.1-A and qBr6.1-B, within the qBr6.1. The QTLs qBr6.1-A and qBr6.1-B explained 37.46% and 10.63% of bruchid resistance variation, respectively. qBr6.1-A was mapped to a 28.24 kb region containing four genes, from which the gene VvTaXI encoding a xylanase inhibitor was selected as a candidate gene responsible for the resistance associated with the qBr6.1-A. Sequencing and sequence alignment of VvTaXI from TVNu 240 and TVNu 1623 revealed a 1-base-pair insertion/deletion and five single-nucleotide polymorphisms (SNPs) in the 5' UTR and 11 SNPs in the exon. Alignment of the VvTAXI protein sequences showed five amino acid changes between the TVNu 240 and TVNu 1623 sequences. Altogether, these results demonstrated that the VvTaXI encoding xylanase inhibitor is the candidate gene conferring bruchid resistance in the zombi pea accession TVNu 240. The gene VvTaXI will be useful for the molecular breeding of bruchid resistance in the zombi pea.

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.

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.

The First Genetic Linkage Map of Winged Bean [Psophocarpus tetragonolobus (L.) DC.] and QTL Mapping for Flower-, Pod-, and Seed-Related Traits.

Winged bean [Psophocarpus tetragonolobus (L.) DC.] (2n = 2× = 18) is a tropical legume crop with multipurpose usages. Recently, the winged bean has regained attention from scientists as a food protein source. Currently, there is no breeding program for winged bean cultivars. All winged bean cultivars are landraces or selections from landraces. Molecular markers and genetic linkage maps are pre-requisites for molecular plant breeding. The aim of this study was to develop a high-density linkage map and identify quantitative trait loci (QTLs) for pod and seed-related traits of the winged bean. An F2 population of 86 plants was developed from a cross between winged bean accessions W054 and TPT9 showing contrasting pod length, and pod, flower and seed colors. A genetic linkage map of 1384 single nucleotide polymorphism (SNP) markers generated from restriction site-associated DNA sequencing was constructed. The map resolved nine haploid chromosomes of the winged bean and spanned the cumulative length of 4552.8 cM with the number of SNPs per linkage ranging from 36 to 218 with an average of 153.78. QTL analysis in the F2 population revealed 31 QTLs controlling pod length, pod color, pod anthocyanin content, flower color, and seed color. The number of QTLs per trait varied between 1 (seed length) to 7 (banner color). Interestingly, the major QTLs for pod color, anthocyanin content, and calyx color, and for seed color and flower wing color were located at the same position. The high-density linkage map QTLs reported in this study will be useful for molecular breeding of winged beans.

Genetic diversity and structure of landrace of lablab (Lablab purpureus (L.) Sweet) cultivars in Thailand revealed by SSR markers.

Lablab (Lablab purpureus (L.) Sweet) is a legume crop widely cultivated in tropical and subtropical regions of Africa and Asia. In this study, we assessed genetic diversity and population structure of 299 individuals of subspecies purpureus and bengalensis of lablab from Thailand using 13 simple sequence repeat (SSR) markers. The SSR markers detected only 34 alleles in total with a mean of 2.6 alleles per locus. Overall gene diversity was 0.360. Gene diversity (H E) and allelic richness (A R) in different geographic regions was comparable. Similarly, both H E and A R between subspecies purpureus and bengalensis were similar. STRUCTURE and neighbor-joining (NJ) analyses revealed that the 299 individuals were clustered into two major groups. In contrast, principal coordinate analysis (PCoA) revealed admixture of the lablab germplasm. STRUCTURE, NJ and PCoA analyses also revealed that the subspecies purpureus and bengalensis are not genetically differentiated. Although the number of individuals from the west of Thailand was small and all of them were collected from the same province, they possessed comparable gene diversity with those from the other geographic regions. These results demonstrated that there is moderately low genetic diversity of lablab in Thailand and the west of the country possesses high diversity of lablab.

Mapping and Functional Characterization of Stigma Exposed 1, a DUF1005 Gene Controlling Petal and Stigma Cells in Mungbean (Vigna radiata).

Flowers with exposed stigma increase the outcrossing rate and are useful in developing improved hybrid crop cultivars. This exposure results mainly from the cellular morphology of the petal and pistil, but what affects the formation of the petal and pistil in the late developmental stages is less understood. Here, we characterized a novel floral mutant in mungbean (Vigna radiata), stigma exposed 1 (se1), which displays irregular petals and pistils. Floral organ initiation in the se1 mutant was normal, but petal and pistil growth malfunctioned during late development. A histological analysis revealed that the se1 mutant had wrinkled petals with knotted structures and elongated styles. The cellular morphology of the epidermal layers of the se1 petals was deformed, while the cell lengths in the styles increased. A genetic analysis indicated that the se1 phenotype is controlled by a single recessive gene, and it was mapped to chromosome 11. A sequence analysis suggested that a DUF1005-encoding gene, LOC106777793, is the gene controlling the se1 phenotype. The se1 mutant possessed a single-nucleotide polymorphism that resulted in an amino acid change in VrDUF1005. Overexpression of VrDUF1005 in Arabidopsis resulted in rolling leaves and reduced floral size. Consequently, we proposed that VrSE1 functions to modulate cell division in petals and cell expansion in styles during the late developmental stages in mungbean. The se1 mutant is a new genetic resource for mung bean hybrid breeding.

Identification of QTLs for Domestication-Related Traits in Zombi Pea [Vigna vexillata (L.) A. Rich], a Lost Crop of Africa.

Zombi pea [Vigna vexillata (L.) A. Rich] is a legume crop found in Africa. Wild zombi pea is widely distributed throughout the tropical and subtropical regions, whereas domesticated zombi pea is rarely cultivated. Plant domestication is an evolutionary process in which the phenotypes of wild species, including seed dormancy, pod shattering, organ size, and architectural and phenological characteristics, undergo changes. The molecular mechanism underlying the domestication of zombi pea is relatively unknown. In this study, the genetic basis of the following 13 domestication-related traits was investigated in an F2 population comprising 198 individuals derived from a cross between cultivated (var. macrosperma) and wild (var. vexillata) zombi pea accessions: seed dormancy, pod shattering, days-to-flowering, days-to-maturity, stem thickness, stem length, number of branches, leaf area, pod length, 100-seed weight, seed width, seed length, and seeds per pod. A genetic map containing 6,529 single nucleotide polymorphisms constructed for the F2 population was used to identify quantitative trait loci (QTLs) for these traits. A total of 62 QTLs were identified for the 13 traits, with 1-11 QTLs per trait. The major QTLs for days-to-flowering, stem length, number of branches, pod length, 100-seed weight, seed length, and seeds per pod were clustered in linkage group 5. In contrast, the major QTLs for seed dormancy and pod shattering belonged to linkage groups 3 and 11, respectively. A comparative genomic analysis with the cowpea [Vigna unguiculata (L.) Walp.] genome used as the reference sequence (i.e., the genome of the legume species most closely related to zombi pea) enabled the identification of candidate genes for the major QTLs. Thus, we revealed the genomic regions associated with domestication-related traits and the candidate genes controlling these traits in zombi pea. The data presented herein may be useful for breeding new varieties of zombi pea and other Vigna species.

Construction of a high density linkage map and genome dissection of bruchid resistance in zombi pea (Vigna vexillata (L.) A. Rich).

Zombi pea (Vigna vexillata) is a legume crop that is resistant to several biotic and abiotic stresses. Callosobruchus maculatus and Callosobruchus chinensis are serious stored-insect pests of legume crops. We constructed a high-density linkage map and performed quantitative trait loci (QTLs) mapping for resistance to these insect species in zombi pea. An F2 population of 198 individuals from a cross between 'TVNu 240' (resistant) and 'TVNu 1623' (susceptible) varieties was used to construct a linkage map of 6,529 single nucleotide polymorphism markers generated from sequencing amplified fragments of specific loci. The map comprised 11 linkage groups, spanning 1,740.9 cM, with an average of 593.5 markers per linkage group and an average distance of 0.27 cM between markers. High levels of micro-synteny between V. vexillata and cowpea (Vigna unguiculata), mungbean (Vigna radiata), azuki bean (Vigna angularis) and common bean (Phaseolus vulgaris) were found. One major and three minor QTLs for C. chinensis resistance and one major and one minor QTLs for C. maculatus resistance were identified. The major QTLs for resistance to C. chinensis and C. maculatus appeared to be the same locus. The linkage map developed in this study will facilitate the identification of useful genes/QTLs in zombi pea.

Construction of genetic linkage map and genome dissection of domestication-related traits of moth bean (Vigna aconitifolia), a legume crop of arid areas.

The moth bean (Vigna aconitifolia), possibly the most primitive crop of the genus Vigna, is a highly drought- and heat-resistant legume grown in arid areas. Moth bean domestication involved phenotypic changes, including reduction of seed dormancy and pod shattering, increased organ size, and earlier flowering and maturity. However, the genetics of the domestication process in moth bean is not known. In this study, we constructed a genetic linkage map for moth bean and used the map to identify quantitative trait loci (QTL) for domestication-related traits of an F2 population of 188 individuals produced from a cross of wild moth bean (TN67) and cultivated moth bean (ICPMO056). The genetic linkage map comprised 11 linkage groups (LG) of 172 simple sequence repeat markers and spanned a total length of 1016.8 centiMorgan (cM), with an average marker distance of 7.34 cM. A comparative genome analysis showed high genome synteny between moth bean and mungbean (Vigna radiata), adzuki bean (Vigna angularis), rice bean (Vigna umbellata), and yardlong bean (Vigna unguiculata). In total, 50 QTLs and 3 genes associated with 20 domestication-related traits were identified. Most of the QTLs belonged to five LGs (1, 2, 4, 7, and 10). Key traits related to domestication such as seed dormancy and pod shattering were controlled by large-effect QTLs (PVE > 20%) with one or two minor QTLs, whereas all other traits were controlled by one-seven minor QTLs, apart from seed weight, which was controlled by one major and seven minor QTLs. These results suggest that a small number of mutations with large phenotypic effects have contributed to the domestication of the moth bean. Comparative analysis of QTLs with related Vigna crops revealed that there are several domestication-related large-effect QTLs that had not been used in moth bean domestication. This study provides a basic genetic map and identified genome regions associated with domestication-related traits, which will be useful for the genetic improvement of the moth bean and related Vigna species.

Mapping of QTLs for Seed Phorbol Esters, a Toxic Chemical in Jatropha curcas (L.).

Jatropha (Jatropha curcas L.) is an oil-bearing plant that has potential to be cultivated as a biodiesel crop. The seed cake after oil extraction has 40-50% protein that can be used in animal feeds. A major limitation in utilizing the cake is the presence of phorbol esters (PE), a heat-tolerant toxic chemical. To identify the quantitative trait loci (QTLs) for PE, we constructed a genetic linkage map from an F2 population of 95 individuals from a cross "Chai Nat" × "M10" using 143 simple sequence repeat (SSR) markers. M10 is low in seed PE while Chai Nat is high. Seeds from each F2 individual were quantified for PE content by high performance liquid chromatography. A single marker analysis revealed five markers from linkage group 3 (LG3) and nine markers from LG8 associated with seed PE. Inclusive composite interval mapping identified two QTLs, each on LG3 (qPE3.1) and LG8 (qPE8.1) responsible for the PE. qPE3.1 and qPE8.1 accounted for 14.10%, and 15.49% of total variation in seed PE, respectively. Alelle(s) from M10 at qPE3.1 increased seed PE, while at qPE8.1 decreased seed PE. qPE3.1 is a new loci for PE, while qPE8.1 is the same locus with that reported recently for PE.