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.

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.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding.

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Thirty Years of Mungbean Genome Research: Where Do We Stand and What Have We Learned?

Mungbean is a socioeconomically important legume crop in Asia that is currently in high demand by consumers and industries both as dried beans and in plant-based protein foods. Marker-assisted and genomics-assisted breeding are promising approaches to efficiently and rapidly develop new cultivars with improved yield, quality, and resistance to biotic and abiotic stresses. Although mungbean was at the forefront of research at the dawn of the plant genomics era 30 years ago, the crop is a "slow runner" in genome research due to limited genomic resources, especially DNA markers. Significant progress in mungbean genome research was achieved only within the last 10 years, notably after the release of the VC1973A draft reference genome constructed using next-generation sequencing technology, which enabled fast and efficient DNA marker development, gene mapping, and identification of candidate genes for complex traits. Resistance to biotic stresses has dominated mungbean genome research to date; however, research is on the rise. In this study, we provide an overview of the past progress and current status of mungbean genomics research. We also discuss and evaluate some research results to provide a better understanding of mungbean genomics.

The mungbean VrP locus encoding MYB90, an R2R3-type MYB protein, regulates anthocyanin biosynthesis.

Anthocyanins are water-soluble pigments present in several tissues/parts of plants. The pigments provide color and are wildly known for health benefits for human, insect attraction for plant pollination, and stress resistance in plants. Anthocyanin content variations in mungbean [Vigna radiata (L.) Wilczek] were first noticed a long time ago, but the genetic mechanism controlling the anthocyanins in mungbean remains unknown. An F2 population derived from the cross between purple-hypocotyl (V2709) and green-hypocotyl (Sulv1) mungbeans was used to map the VrP locus controlling purple hypocotyl. The VrP locus was mapped to a 78.9-kb region on chromosome 4. Sequence comparison and gene expression analysis identified an R2R3-MYB gene VrMYB90 as the candidate gene for the VrP locus. Haplotype analysis using 124 mungbean accessions suggested that 10 single nucleotide polymorphisms (SNPs) in exon 3 may lead to an abolished expression of VrMYB90 and an absence of anthocyanin accumulation in the hypocotyl of Sulv1 and KPS2. The overexpression of VrMYB90 in mungbean hairy root, tobacco leaf, and Arabidopsis resulted in anthocyanin accumulation (purple color). Gene expression analysis demonstrated that VrMYB90 regulated anthocyanin accumulation in the hypocotyl, stem, petiole, and flowers, and the expression was sensitive to light. VrMYB90 protein may upregulate VrDFR encoding dihydroflavonol 4-reductase at the late biosynthesis step of anthocyanins in mungbeans. These results suggest that VrMYB90 is the dominator in the spatiotemporal regulation of anthocyanin biosynthesis. Our results provide insight into the biosynthesis mechanism of anthocyanin and a theoretical basis for breeding mungbeans.

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.

Identification of Pathogenicity Loci in Magnaporthe oryzae Using GWAS with Neck Blast Phenotypic Data.

Magnaporthae oryzae (M. oryzae) is the most destructive disease of rice worldwide. In this study, one hundred and two isolates of M. oryzae were collected from rice (Oryzae sativa L.) from 2001 to 2017, and six rice varieties with resistance genes Pizt, Pish, Pik, Pib, and Pi2 were used in a genome-wide association study to identify pathogenicity loci in M. oryzae. Genome-wide association analysis was performed using 5338 single nucleotide polymorphism (SNPs) and phenotypic data of neck blast screening by TASSEL software together with haplotype block and SNP effect analysis. Twenty-seven significant SNPs were identified on chromosomes 1, 2, 3, 4, 5, 6, and 7. Many predicted genes (820 genes) were found in the target regions of six rice varieties. Most of these genes are described as putative uncharacterized proteins, however, some genes were reported related to virulence in M. oryzae. Moreover, this study revealed that R genes, Pik, Pish, and Pi2, were broad-spectrum resistant against neck blast disease caused by Thai blast isolate. Haplotype analysis revealed that the combination of the favorable alleles causing reduced virulence of isolates against IRBLz5-CA carrying Pi2 gene contributes 69% of the phenotypic variation in pathogenicity. The target regions and information are useful to develop marker-specific genes to classify blast fungal isolates and select appropriate resistance genes for rice cultivation and improvement.

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.

A chromosome-scale assembly of the black gram (Vigna mungo) genome.

Black gram (Vigna mungo) is an important short duration grain legume crop. Black gram seeds provide an inexpensive source of dietary protein. Here, we applied the 10X Genomics linked-read technology to obtain a de novo whole genome assembly of V. mungo cultivated variety Chai Nat 80 (CN80). The preliminary assembly contained 12,228 contigs and had an N50 length of 5.2 Mb. Subsequent scaffolding using the long-range Chicago and HiC techniques yielded the first high-quality, chromosome-level assembly of 499 Mb comprising 11 pseudomolecules. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that black gram and mungbean (Vigna radiata) diverged about 2.7 million years ago . The transversion rate (4DTv) analysis in V. mungo revealed no evidence supporting a recent genome-wide duplication event observed in the tetraploid créole bean (Vigna reflexo-pilosa). The proportion of repetitive elements in the black gram genome is slightly lower than the numbers reported for related Vigna species. The majority of long terminal repeat retrotransposons appeared to integrate into the genome within the last five million years. We also examined alternative splicing events in V. mungo using full-length transcript sequences. While intron retention was the most prevalent mode of alternative splicing in several plant species, alternative 3' acceptor site selection represented the majority of events in black gram. Our high-quality genome assembly along with the genomic variation information from the germplasm provides valuable resources for accelerating the development of elite varieties through marker-assisted breeding and for future comparative genomics and phylogenetic studies in legume species.

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.

The Genome and Transcriptome Analysis of the Vigna mungo Chloroplast.

Vigna mungo is cultivated in approximately 5 million hectares worldwide. The chloroplast genome of this species has not been previously reported. In this study, we sequenced the genome and transcriptome of the V. mungo chloroplast. We identified many positively selected genes in the photosynthetic pathway (e.g., rbcL, ndhF, and atpF) and RNA polymerase genes (e.g., rpoC2) from the comparison of the chloroplast genome of V. mungo, temperate legume species, and tropical legume species. Our transcriptome data from PacBio isoform sequencing showed that the 51-kb DNA inversion could affect the transcriptional regulation of accD polycistronic. Using Illumina deep RNA sequencing, we found RNA editing of clpP in the leaf, shoot, flower, fruit, and root tissues of V. mungo. We also found three G-to-A RNA editing events that change guanine to adenine in the transcripts transcribed from the adenine-rich regions of the ycf4 gene. The edited guanine bases were found particularly in the chloroplast genome of the Vigna species. These G-to-A RNA editing events were likely to provide a mechanism for correcting DNA base mutations. The V. mungo chloroplast genome sequence and the analysis results obtained in this study can apply to phylogenetic studies and chloroplast genome engineering.

Genome sequence of Jatropha curcas L., a non-edible biodiesel plant, provides a resource to improve seed-related traits.

Jatropha curcas (physic nut), a non-edible oilseed crop, represents one of the most promising alternative energy sources due to its high seed oil content, rapid growth and adaptability to various environments. We report ~339 Mbp draft whole genome sequence of J. curcas var. Chai Nat using both the PacBio and Illumina sequencing platforms. We identified and categorized differentially expressed genes related to biosynthesis of lipid and toxic compound among four stages of seed development. Triacylglycerol (TAG), the major component of seed storage oil, is mainly synthesized by phospholipid:diacylglycerol acyltransferase in Jatropha, and continuous high expression of homologs of oleosin over seed development contributes to accumulation of high level of oil in kernels by preventing the breakdown of TAG. A physical cluster of genes for diterpenoid biosynthetic enzymes, including casbene synthases highly responsible for a toxic compound, phorbol ester, in seed cake, was syntenically highly conserved between Jatropha and castor bean. Transcriptomic analysis of female and male flowers revealed the up-regulation of a dozen family of TFs in female flower. Additionally, we constructed a robust species tree enabling estimation of divergence times among nine Jatropha species and five commercial crops in Malpighiales order. Our results will help researchers and breeders increase energy efficiency of this important oil seed crop by improving yield and oil content, and eliminating toxic compound in seed cake for animal feed.

Genetics of resistance to Cercospora leaf spot disease caused by Cercospora canescens and Psuedocercospora cruenta in yardlong bean (Vigna unguiculata ssp. sesquipedalis) × grain cowpea (V. unguiculata ssp. unguiculata) populations.

Yardlong bean (Vigna unguiculata ssp. sesquipedalis), a type of cowpea, is an important vegetable legume of Asia. Cercospora leaf spot (CLS) caused by Cercospora canescens and Psuedocercospora cruenta is an important phytopathological problem of the yardlong bean grown in tropical regions. The objectives of this study were to (i) determine mode of inheritance of resistance to CLS caused by C. canescens and P. cruenta, (ii) estimate the heritability of the resistance, (iii) estimate genetic effects on the resistance using six basic populations generated from the cross between the susceptible yardlong bean 'CSR12906' and the resistant grain cowpea (V.unguiculata spp. unguiculata) 'IT90K-59-120'. Segregation for the resistance to both fungi in the F2 population fitted both 3 : 1 ratio and 13 : 3 ratio of susceptible:resistant, while that in the BC2 ((CSR12906×IT90K-59-120)×IT90K- 59-120) population fitted a 1 : 1 ratio, suggesting one recessive gene or two genes with inhibitory gene action control the resistance. Generation mean analysis showed that a simple additive-dominance model was adequate to explain the genetic control of CLS disease resistance, indicating that a single gene controls the resistance. The average number of major genes (effective factors) controlling the resistance was estimated to be 1.05 and 0.92 for C. canescens and P. cruenta, respectively. The broad-sense heritability calculated for resistance to both diseases was higher than 0.90. Altogether, these results indicated that the resistance to CLS disease caused by C. canescens and P. cruenta in grain cowpea IT90K-59-120 is a highly heritable trait governed by a single major recessive gene.

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.

Draft genome sequence of adzuki bean, Vigna angularis.

Adzuki bean (Vigna angularis var. angularis) is a dietary legume crop in East Asia. The presumed progenitor (Vigna angularis var. nipponensis) is widely found in East Asia, suggesting speciation and domestication in these temperate climate regions. Here, we report a draft genome sequence of adzuki bean. The genome assembly covers 75% of the estimated genome and was mapped to 11 pseudo-chromosomes. Gene prediction revealed 26,857 high confidence protein-coding genes evidenced by RNAseq of different tissues. Comparative gene expression analysis with V. radiata showed that the tissue specificity of orthologous genes was highly conserved. Additional re-sequencing of wild adzuki bean, V. angularis var. nipponensis, and V. nepalensis, was performed to analyze the variations between cultivated and wild adzuki bean. The determined divergence time of adzuki bean and the wild species predated archaeology-based domestication time. The present genome assembly will accelerate the genomics-assisted breeding of adzuki bean.