A genetic linkage map for azuki bean [Vigna angularis (Willd.) Ohwi & Ohashi]

Progress in Adzuki Bean Seed Coat Colour Studies

Seed coat colour is an important quality trait, domestication trait, and morphological marker, and it is closely associated with flavonoid and anthocyanin metabolism pathways. The seed coat colour of the adzuki bean, an important legume crop, influences the processing quality, the commodity itself, and its nutritional quality. In this review, a genetic analysis of different seed coat colours, gene mapping, metabolite content determination, and varietal improvement in adzuki bean are summarized. It provides further insight into gene mapping and cloning of seed coat colour genes and varietal improvements in adzuki beans.

Exploiting genetic and genomic resources to enhance productivity and abiotic stress adaptation of underutilized pulses

Underutilized pulses and their wild relatives are typically stress tolerant and their seeds are packed with protein, fibers, minerals, vitamins, and phytochemicals. The consumption of such nutritionally dense legumes together with cereal-based food may promote global food and nutritional security. However, such species are deficient in a few or several desirable domestication traits thereby reducing their agronomic value, requiring further genetic enhancement for developing productive, nutritionally dense, and climate resilient cultivars. This review article considers 13 underutilized pulses and focuses on their germplasm holdings, diversity, crop-wild-crop gene flow, genome sequencing, syntenic relationships, the potential for breeding and transgenic manipulation, and the genetics of agronomic and stress tolerance traits. Recent progress has shown the potential for crop improvement and food security, for example, the genetic basis of stem determinacy and fragrance in moth bean and rice bean, multiple abiotic stress tolerant traits in horse gram and tepary bean, bruchid resistance in lima bean, low neurotoxin in grass pea, and photoperiod induced flowering and anthocyanin accumulation in adzuki bean have been investigated. Advances in introgression breeding to develop elite genetic stocks of grass pea with low β-ODAP (neurotoxin compound), resistance to Mungbean yellow mosaic India virus in black gram using rice bean, and abiotic stress adaptation in common bean, using genes from tepary bean have been carried out. This highlights their potential in wider breeding programs to introduce such traits in locally adapted cultivars. The potential of de-domestication or feralization in the evolution of new variants in these crops are also highlighted.

Unique Salt-Tolerance-Related QTLs, Evolved in Vigna riukiuensis (Na+ Includer) and V. nakashimae (Na+ Excluder), Shed Light on the Development of Super-Salt-Tolerant Azuki Bean (V. angularis) Cultivars

Wild relatives of crops have the potential to improve food crops, especially in terms of improving abiotic stress tolerance. Two closely related wild species of the traditional East Asian legume crops, Azuki bean (Vigna angularis), V. riukiuensis "Tojinbaka" and V. nakashimae "Ukushima" were shown to have much higher levels of salt tolerance than azuki beans. To identify the genomic regions responsible for salt tolerance in "Tojinbaka" and "Ukushima", three interspecific hybrids were developed: (A) azuki bean cultivar "Kyoto Dainagon" × "Tojinbaka", (B) "Kyoto Dainagon" × "Ukushima" and (C) "Ukushima" × "Tojinbaka". Linkage maps were developed using SSR or restriction-site-associated DNA markers. There were three QTLs for "percentage of wilt leaves" in populations A, B and C, while populations A and B had three QTLs and population C had two QTLs for "days to wilt". In population C, four QTLs were detected for Na⁺ concentration in the primary leaf. Among the F₂ individuals in population C, 24% showed higher salt tolerance than both wild parents, suggesting that the salt tolerance of azuki beans can be further improved by combining the QTL alleles of the two wild relatives. The marker information would facilitate the transfer of salt tolerance alleles from "Tojinbaka" and "Ukushima" to azuki beans.

Vignette of Vigna domestication: From archives to genomics

The genus Vigna comprises fast-growing, diploid legumes, cultivated in tropical and subtropical parts of the world. It comprises more than 200 species among which Vigna angularis, Vigna radiata, Vigna mungo, Vigna aconitifolia, Vigna umbellata, Vigna unguiculata, and Vigna vexillata are of enormous agronomic importance. Human selection along with natural variability within these species encompasses a vital source for developing new varieties. The present review convokes the early domestication history of Vigna species based on archeological pieces of evidence and domestication-related traits (DRTs) together with genetics of domestication. Traces of early domestication of Vigna have been evidenced to spread across several temperate and tropical regions of Africa, Eastern Asia, and few parts of Europe. Several DRTs of Vigna species, such as pod shattering, pod and seed size, dormancy, seed coat, seed color, maturity, and pod dehiscence, can clearly differentiate wild species from their domesticates. With the advancement in next-generation high-throughput sequencing techniques, exploration of genetic variability using recently released reference genomes along with de novo sequencing of Vigna species have provided a framework to perform genome-wide association and functional studies to figure out different genes related to DRTs. In this review, genes and quantitative trait loci (QTLs) related to DRTs of different Vigna species have also been summarized. Information provided in this review will enhance the in-depth understanding of the selective pressures that causes crop domestication along with nature of evolutionary selection made in unexplored Vigna species. Furthermore, correlated archeological and domestication-related genetic evidence will facilitate Vigna species to be considered as suitable model plants.

The bZIP Transcription Factor Family in Adzuki Bean (Vigna Angularis): Genome-Wide Identification, Evolution, and Expression Under Abiotic Stress During the Bud Stage

Adzuki bean (Vigna angularis) is an important dietary legume crop that was first cultivated and domesticated in Asia. Currently, little is known concerning the evolution and expression patterns of the basic leucine zipper (bZIP) family transcription factors in the adzuki bean. Through the PFAM search, 72 bZIP members of adzuki bean (VabZIP) were identified from the reference genome. Most of them were located on 11 chromosomes and seven on an unknown chromosome. A comprehensive analysis, including evolutionary, motifs, gene structure, cis-elements, and collinearity was performed to identify VabZIP members. The subcellular localization results showed VabZIPs might locate on the nuclear. Quantitative real-time PCR (qRT-PCR) analysis of the relative expression of VabZIPs in different tissues at the bud stage revealed that VabZIPs had a tissue-specific expression pattern, and its expression was influenced by abiotic stress. These characteristics of VabZIPs provide insights for future research aimed at developing interventions to improve abiotic stress resistance.

BAC- and oligo-FISH mapping reveals chromosome evolution among Vigna angularis, V. unguiculata, and Phaseolus vulgaris

Cytogenomic resources have accelerated synteny and chromosome evolution studies in plant species, including legumes. Here, we established the first cytogenetic map of V. angularis (Va, subgenus Ceratotropis) and compared this new map with those of V. unguiculata (Vu, subgenus Vigna) and P. vulgaris (Pv) by BAC-FISH and oligopainting approaches. We mapped 19 Vu BACs and 35S rDNA probes to the 11 chromosome pairs of Va, Vu, and Pv. Vigna angularis shared a high degree of macrosynteny with Vu and Pv, with five conserved syntenic chromosomes. Additionally, we developed two oligo probes (Pv2 and Pv3) used to paint Vigna orthologous chromosomes. We confirmed two reciprocal translocations (chromosomes 2 and 3 and 1 and 8) that have occurred after the Vigna and Phaseolus divergence (~9.7 Mya). Besides, two inversions (2 and 4) and one translocation (1 and 5) have occurred after Vigna and Ceratotropis subgenera separation (~3.6 Mya). We also observed distinct oligopainting patterns for chromosomes 2 and 3 of Vigna species. Both Vigna species shared similar major rearrangements compared to Pv: one translocation (2 and 3) and one inversion (chromosome 3). The sequence synteny identified additional inversions and/or intrachromosomal translocations involving pericentromeric regions of both orthologous chromosomes. We propose chromosomes 2 and 3 as hotspots for chromosomal rearrangements and de novo centromere formation within and between Vigna and Phaseolus. Our BAC- and oligo-FISH mapping contributed to physically trace the chromosome evolution of Vigna and Phaseolus and its application in further studies of both genera.

Construction of a High-Density Genetic Map and Its Application for QTL Mapping of Leaflet Shapes in Mung Bean (Vigna radiata L.)

Mung bean (Vigna radiata L.) is an important but understudied food legume in Asia and now worldwide. Genetic studies may help to accelerate the exploitation of new genes for breeding in this crop. Here, we used a recombination inbred line population to construct an SNP genetic linkage map by genome sequencing technology. We obtained 21,508 high-quality SNP markers integrated into 1,946 bin markers that were mapped onto 11 linkage groups (LGs) with 99-258 bin markers per LG. The total genetic length of the map was 1060.2 cM (38.76-168.03 cM per LG), with an average distance between markers of 0.54 cM. However, there were 18 gaps >5 cM, distribution on LG1, 3, 5, 7, and 9. Gene mapping for lobed and indented leaflets was conducted using the map. A major quantitative trait locus (QTL) associated with indented leaflets was detected on chromosome 10, with phenotypic variation explained (PVE) values of 39.7% and 45.4% under two different environments. Several QTLs for lobed leaflets were detected and most of them were tightly linked together on Chromosome 3. However, only one major QTL, which explained the largest phenotypic variation (27.7-69.5%), was stably detected under two different environments using both R and Q methods. In the two main stable QTLs regions on chromosomes 3 and 10, candidate genes for regulating the molecular mechanism of different leaflet shapes were detected by functional annotation. The overlap of major QTLs under different environments indicated that the present map would be good enough for precisely mapping genes, and both the QTL analysis and gene prediction were useful for investigating the mechanism of leaf development in mung bean or legumes.

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 F₂ 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.

QTL analysis of domestication syndrome in zombi pea (Vigna vexillata), an underutilized legume crop

Zombi pea (Vigna vexillata (L.) A. Rich) is an underutilized crop belonging to the genus Vigna. Two domesticated forms of zombi pea are cultivated as crop plants; seed and tuber forms. The cultivated seed form is present in Africa, while the cultivated tuber form is present in a very limited part of Asia. Genetics of domestication have been investigated in most of cultivated Vigna crops by means of quantitative trait locus (QTL) mapping. In this study, we investigated genetics of domestication in zombi pea by QTL analysis using an F2 population of 139 plants derived from a cross between cultivated tuber form of V. vexillata (JP235863) and wild V. vexillata (AusTRCF66514). A linkage map with 11 linkage groups (LGs) was constructed from this F2 population using 145 SSR, 117 RAD-seq and 2 morphological markers. Many highly segregation distorted markers were found on LGs 5, 6, 7, 8, 10 and 11. Most of the distorted markers were clustered together and all the markers on LG8 were highly distorted markers. Comparing this V. vexillata linkage map with linkage maps of other four Vigna species demonstrated several genome rearrangements in V. vexillata. QTL analysis for 22 domestication-related traits was investigated by inclusive composite interval mapping in which 37 QTLs were identified for 18 traits; no QTL was detected for 4 traits. Number of QTLs detected in each trait ranged from 1 to 5 with an average of only 2.3. Five QTLs for tuber width and three QTLs for tuber weight. Interestingly, 2 QTLs each for tuber width and tuber weight detected on LG2 and LG4 were located at similar position and wild allele increased tuber width and weight. This indicated wild germplasm having small tuber have potential to increase yield of large tuber cultivated type. Large-effect QTLs (PVE > 20%) were on LG4 (pod length), LG5 (leaf size and seed thickness), and LG7 (for seed-related traits). Comparison of domestication-related QTLs of the zombi pea with those of cowpea (Vigna unguiculata), azuki bean (Vigna angularis), mungbean (Vigna radiata) and rice bean (Vigna umbellata) revealed that there was conservation of some QTLs for seed size, pod size and leaf size between zombi pea and cowpea and that QTLs associated with seed size (weight, length, width and thickness) in each species were clustered on same linkage.

Identification of QTL and Qualitative Trait Loci for Agronomic Traits Using SNP Markers in the Adzuki Bean

The adzuki bean (Vigna angularis) is an important grain legume. Fine mapping of quantitative trait loci (QTL) and qualitative trait genes plays an important role in gene cloning, molecular-marker-assisted selection (MAS), and trait improvement. However, the genetic control of agronomic traits in the adzuki bean remains poorly understood. Single-nucleotide polymorphisms (SNPs) are invaluable in the construction of high-density genetic maps. We mapped 26 agronomic QTLs and five qualitative trait genes related to pigmentation using 1,571 polymorphic SNP markers from the adzuki bean genome via restriction-site-associated DNA sequencing of 150 members of an F₂ population derived from a cross between cultivated and wild adzuki beans. We mapped 11 QTLs for flowering time and pod maturity on chromosomes 4, 7, and 10. Six 100-seed weight (SD100WT) QTLs were detected. Two major flowering time QTLs were located on chromosome 4, firstly VaFld4.1 (PEVs 71.3%), co-segregating with SNP marker s690-144110, and VaFld4.2 (PEVs 67.6%) at a 0.974 cM genetic distance from the SNP marker s165-116310. Three QTLs for seed number per pod (Snp3.1, Snp3.2, and Snp4.1) were mapped on chromosomes 3 and 4. One QTL VaSdt4.1 of seed thickness (SDT) and three QTLs for branch number on the main stem were detected on chromosome 4. QTLs for maximum leaf width (LFMW) and stem internode length were mapped to chromosomes 2 and 9, respectively. Trait genes controlling the color of the seed coat, pod, stem and flower were mapped to chromosomes 3 and 1. Three candidate genes, VaAGL, VaPhyE, and VaAP2, were identified for flowering time and pod maturity. VaAGL encodes an agamous-like MADS-box protein of 379 amino acids. VaPhyE encodes a phytochrome E protein of 1,121 amino acids. Four phytochrome genes (VaPhyA1, VaPhyA2, VaPhyB, and VaPhyE) were identified in the adzuki bean genome. We found candidate genes VaAP2/ERF.81 and VaAP2/ERF.82 of SD100WT, VaAP2-s4 of SDT, and VaAP2/ERF.86 of LFMW. A candidate gene VaUGT related to black seed coat color was identified. These mapped QTL and qualitative trait genes provide information helpful for future adzuki bean candidate gene cloning and MAS breeding to improve cultivars with desirable growth periods, yields, and seed coat color types.

Genetic mapping and QTL analysis of agronomic traits in Indian Mucuna pruriens using an intraspecific F₂population

Mucuna pruriens is a well-recognized agricultural and horticultural crop with important medicinal use. However, antinutritional factors in seed and adverse morphological characters have negatively affected its cultivation. To elucidate the genetic control of agronomic traits, an intraspecific genetic linkage map of Indian M. pruriens has been developed based on amplified fragment length polymorphism (AFLP) markers using 200 F₂ progenies derived from a cross between wild and cultivated genotypes. The resulting linkage map comprised 129 AFLP markers dispersed over 13 linkage groups spanning a total distance of 618.88 cM with an average marker interval of 4.79 cM. For the first time, three QTLs explaining about 6.05-14.77% of the corresponding total phenotypic variation for three quantitative (seed) traits and, eight QTLs explaining about 25.96% of the corresponding total phenotypic variation for three qualitative traits have been detected on four linkage groups. The map presented here will pave a way for mapping of genes/QTLs for the important agronomic and horticultural traits contrasting between the parents used in this study.

Development of a high-density genetic linkage map and identification of flowering time QTLs in adzuki bean (Vigna angularis)

A high-density linkage map is crucial for the identification of quantitative trait loci (QTLs), positional cloning, and physical map assembly. Here, we report the development of a high-density linkage map based on specific length amplified fragment sequencing (SLAF-seq) for adzuki bean and the identification of flowering time-related QTLs. Through SLAF library construction and Illumina sequencing of a recombinant inbred line (RIL) population, a total of 4425 SLAF markers were developed and assigned to 11 linkage groups (LGs). After binning the SLAF markers that represented the same genotype, the final linkage map of 1628.15 cM contained 2032 markers, with an average marker density of 0.80 cM. Comparative analysis showed high collinearity with two adzuki bean physical maps and a high degree of synteny with the reference genome of common bean (Phaseolus vulgaris). Using this map, one major QTL on LG03 and two minor QTLs on LG05 associated with first flowering time (FLD) were consistently identified in tests over a two-year period. These results provide a foundation that will be useful for future genomic research, such as identifying QTLs for other important traits, positional cloning, and comparative mapping in legumes.

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

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.

The Vigna Genome Server, 'VigGS': A Genomic Knowledge Base of the Genus Vigna Based on High-Quality, Annotated Genome Sequence of the Azuki Bean, Vigna angularis (Willd.) Ohwi & Ohashi

The genus Vigna includes legume crops such as cowpea, mungbean and azuki bean, as well as >100 wild species. A number of the wild species are highly tolerant to severe environmental conditions including high-salinity, acid or alkaline soil; drought; flooding; and pests and diseases. These features of the genus Vigna make it a good target for investigation of genetic diversity in adaptation to stressful environments; however, a lack of genomic information has hindered such research in this genus. Here, we present a genome database of the genus Vigna, Vigna Genome Server ('VigGS', http://viggs.dna.affrc.go.jp), based on the recently sequenced azuki bean genome, which incorporates annotated exon-intron structures, along with evidence for transcripts and proteins, visualized in GBrowse. VigGS also facilitates user construction of multiple alignments between azuki bean genes and those of six related dicot species. In addition, the database displays sequence polymorphisms between azuki bean and its wild relatives and enables users to design primer sequences targeting any variant site. VigGS offers a simple keyword search in addition to sequence similarity searches using BLAST and BLAT. To incorporate up to date genomic information, VigGS automatically receives newly deposited mRNA sequences of pre-set species from the public database once a week. Users can refer to not only gene structures mapped on the azuki bean genome on GBrowse but also relevant literature of the genes. VigGS will contribute to genomic research into plant biotic and abiotic stresses and to the future development of new stress-tolerant crops.

The power of single molecule real-time sequencing technology in the de novo assembly of a eukaryotic genome

Second-generation sequencers (SGS) have been game-changing, achieving cost-effective whole genome sequencing in many non-model organisms. However, a large portion of the genomes still remains unassembled. We reconstructed azuki bean (Vigna angularis) genome using single molecule real-time (SMRT) sequencing technology and achieved the best contiguity and coverage among currently assembled legume crops. The SMRT-based assembly produced 100 times longer contigs with 100 times smaller amount of gaps compared to the SGS-based assemblies. A detailed comparison between the assemblies revealed that the SMRT-based assembly enabled a more comprehensive gene annotation than the SGS-based assemblies where thousands of genes were missing or fragmented. A chromosome-scale assembly was generated based on the high-density genetic map, covering 86% of the azuki bean genome. We demonstrated that SMRT technology, though still needed support of SGS data, achieved a near-complete assembly of a eukaryotic genome.

Construction of an SSR and RAD-Marker Based Molecular Linkage Map of Vigna vexillata (L.) A. Rich

Vigna vexillata (L.) A. Rich. (tuber cowpea) is an underutilized crop for consuming its tuber and mature seeds. Wild form of V. vexillata is a pan-tropical perennial herbaceous plant which has been used by local people as a food. Wild V. vexillata has also been considered as useful gene(s) source for V. unguiculata (cowpea), since it was reported to have various resistance gene(s) for insects and diseases of cowpea. To exploit the potential of V. vexillata, an SSR-based linkage map of V. vexillata was developed. A total of 874 SSR markers successfully amplified single DNA fragment in V. vexillata among 1,336 SSR markers developed from Vigna angularis (azuki bean), V. unguiculata and Phaseolus vulgaris (common bean). An F2 population of 300 plants derived from a cross between salt resistant (V1) and susceptible (V5) accessions was used for mapping. A genetic linkage map was constructed using 82 polymorphic SSR markers loci, which could be assigned to 11 linkage groups spanning 511.5 cM in length with a mean distance of 7.2 cM between adjacent markers. To develop higher density molecular linkage map and to confirm SSR markers position in a linkage map, RAD markers were developed and a combined SSR and RAD markers linkage map of V. vexillata was constructed. A total of 559 (84 SSR and 475 RAD) markers loci could be assigned to 11 linkage groups spanning 973.9 cM in length with a mean distance of 1.8 cM between adjacent markers. Linkage and genetic position of all SSR markers in an SSR linkage map were confirmed. When an SSR genetic linkage map of V. vexillata was compared with those of V. radiata and V. unguiculata, it was suggested that the structure of V. vexillata chromosome was considerably differentiated. This map is the first SSR and RAD marker-based V. vexillata linkage map which can be used for the mapping of useful traits.

Reciprocal translocation identified in Vigna angularis dominates the wild population in East Japan

Using an F2 population derived from cultivated and wild azuki bean, we previously detected a reciprocal translocation and a seed size QTL near the translocation site. To test the hypothesis that the translocation in the cultivated variety contributed to the larger seed size, we performed further linkage analyses with several cross combinations between cultivated and wild azuki beans. In addition, we visually confirmed the translocation by cytogenetic approach using 25 wild and cultivated accessions. As a result, we found the translocation-type chromosomes in none of the cultivated accessions, but in a number of the wild accessions. Interestingly, all the wild accessions with the translocation were originally collected from East Japan, while all the accessions with normal chromosomes were from West Japan or the Sea of Japan-side region. Such biased geographical distribution could be explained by the glacial refugium hypothesis, and supported narrowing down the domestication origin of cultivated azuki bean.

Genomics and molecular breeding in lesser explored pulse crops: current trends and future opportunities

Pulses are multipurpose crops for providing income, employment and food security in the underprivileged regions, notably the FAO-defined low-income food-deficit countries. Owing to their intrinsic ability to endure environmental adversities and the least input/management requirements, these crops remain central to subsistence farming. Given their pivotal role in rain-fed agriculture, substantial research has been invested to boost the productivity of these pulse crops. To this end, genomic tools and technologies have appeared as the compelling supplement to the conventional breeding. However, the progress in minor pulse crops including dry beans (Vigna spp.), lupins, lablab, lathyrus and vetches has remained unsatisfactory, hence these crops are often labeled as low profile or lesser researched. Nevertheless, recent scientific and technological breakthroughs particularly the next generation sequencing (NGS) are radically transforming the scenario of genomics and molecular breeding in these minor crops. NGS techniques have allowed de novo assembly of whole genomes in these orphan crops. Moreover, the availability of a reference genome sequence would promote re-sequencing of diverse genotypes to unlock allelic diversity at a genome-wide scale. In parallel, NGS has offered high-resolution genetic maps or more precisely, a robust genetic framework to implement whole-genome strategies for crop improvement. As has already been demonstrated in lupin, sequencing-based genotyping of the representative sample provided access to a number of functionally-relevant markers that could be deployed straight away in crop breeding programs. This article attempts to outline the recent progress made in genomics of these lesser explored pulse crops, and examines the prospects of genomics assisted integrated breeding to enhance and stabilize crop yields.

Detection of genome donor species of neglected tetraploid crop Vigna reflexo-pilosa (créole bean), and genetic structure of diploid species based on newly developed EST-SSR markers from azuki bean (Vigna angularis)

Vigna reflexo-pilosa, which includes a neglected crop, is the only one tetraploid species in genus Vigna. The ancestral species that make up this allotetraploid species have not conclusively been identified, although previous studies suggested that a donor genome of V. reflexo-pilosa is V. trinervia. In this study, 1,429 azuki bean EST-SSR markers were developed of which 38 EST-SSR primer pairs that amplified one product in diploid species and two discrete products in tetraploid species were selected to analyze 268 accessions from eight taxa of seven Asian Vigna species including V. reflexo-pilosa var. glabra, V. reflexo-pilosa var. reflexo-pilosa, V. exilis, V. hirtella, V. minima, V. radiata var. sublobata, V. tenuicaulis and V. trinervia to identify genome donor of V. reflexo-pilosa. Since both diploid and tetraploid species were analyzed and each SSR primer pair detected two loci in the tetraploid species, we separated genomes of the tetraploid species into two different diploid types, viz. A and B. In total, 445 alleles were detected by 38 EST-SSR markers. The highest gene diversity was observed in V. hirtella. By assigning the discrete PCR products of V. reflexo-pilosa into two distinguished genomes, we were able to identify the two genome donor parents of créole bean. Phylogenetic and principal coordinate analyses suggested that V. hirtella is a species complex and may be composed of at least three distinct taxa. Both analyses also clearly demonstrated that V. trinervia and one taxon of V. hirtella are the genome donors of V. reflexo-pilosa. Gene diversity indicates that the evolution rate of EST-SSRs on genome B of créole bean might be faster than that on genome A. Species relationship among the Vigna species in relation to genetic data, morphology and geographical distribution are presented.

Detection of quantitative trait loci for mungbean yellow mosaic India virus (MYMIV) resistance in mungbean (Vigna radiata (L.) Wilczek) in India and Pakistan

Yellow mosaic disease (YMD) is one of the major diseases affecting mungbean (Vigna radiata (L.) Wilczek). In this study, we report the mapping of the quantitative trait locus (QTL) for mungbean yellow mosaic India virus (MYMIV) resistance in mungbean. An F8 recombinant inbred line (RIL) mapping population was generated in Thailand from a cross between NM10-12-1 (MYMIV resistance) and KPS2 (MYMIV susceptible). One hundred and twenty-two RILs and their parents were evaluated for MYMIV resistance in infested fields in India and Pakistan. A genetic linkage map was developed for the RIL population using simple sequence repeat (SSR) markers. Composite interval mapping identified five QTLs for MYMIV resistance: three QTLs for India (qYMIV1, qYMIV2 and qYMIV3) and two QTLs for Pakistan (qYMIV4 and qYMIV5). qYMIV1, qYMIV2, qYMIV3, qYMIV4 and qYMIV5 explained 9.33%, 10.61%, 12.55%, 21.93% and 6.24% of variation in disease responses, respectively. qYMIV1 and qYMIV4 appeared to be the same locus and were common to a major QTL for MYMIV resistance in India identified previously using a different resistant mungbean.

QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields

The objective of this study was to identify quantitative trait loci (QTL) affecting fitness of hybrids between wild soybean (Glycine soja) and cultivated soybean (Glycine max). Seed dormancy and seed number, both of which are important for fitness, were evaluated by testing artificial hybrids of G. soja × G. max in a multiple-site field trial. Generally, the fitness of the F₁ hybrids and hybrid derivatives from self-pollination was lower than that of G. soja due to loss of seed dormancy, whereas the fitness of hybrid derivatives with higher proportions of G. soja genetic background was comparable with that of G. soja. These differences were genetically dissected into QTL for each population. Three QTLs for seed dormancy and one QTL for total seed number were detected in the F₂ progenies of two diverse cross combinations. At those four QTLs, the G. max alleles reduced seed number and severely reduced seed survival during the winter, suggesting that major genes acquired during soybean adaptation to cultivation have a selective disadvantage in natural habitats. In progenies with a higher proportion of G. soja genetic background, the genetic effects of the G. max alleles were not expressed as phenotypes because the G. soja alleles were dominant over the G. max alleles. Considering the highly inbreeding nature of these species, most hybrid derivatives would disappear quickly in early self-pollinating generations in natural habitats because of the low fitness of plants carrying G. max alleles.

Construction of a genetic linkage map and genetic analysis of domestication related traits in mungbean (Vigna radiata)

The genetic differences between mungbean and its presumed wild ancestor were analyzed for domestication related traits by QTL mapping. A genetic linkage map of mungbean was constructed using 430 SSR and EST-SSR markers from mungbean and its related species, and all these markers were mapped onto 11 linkage groups spanning a total of 727.6 cM. The present mungbean map is the first map where the number of linkage groups coincided with the haploid chromosome number of mungbean. In total 105 QTLs and genes for 38 domestication related traits were identified. Compared with the situation in other Vigna crops, many linkage groups have played an important role in the domestication of mungbean. In particular the QTLs with high contribution were distributed on seven out of 11 linkage groups. In addition, a large number of QTLs with small contribution were found. The accumulation of many mutations with large and/or small contribution has contributed to the differentiation between wild and cultivated mungbean. The useful QTLs for seed size, pod dehiscence and pod maturity that have not been found in other Asian Vigna species were identified in mungbean, and these QTLs may play the important role as new gene resources for other Asian Vigna species. The results provide the foundation that will be useful for improvement of mungbean and related legumes.

New Sources of Resistance to Cadophora gregata f. sp. adzukicola and Fusarium oxysporum f. sp. adzukicola in Vigna spp

Adzuki bean brown stem rot (BSR), caused by Cadophora gregata f. sp. adzukicola, and adzuki bean Fusarium wilt (AFW), caused by Fusarium oxysporum f. sp. adzukicola, are serious problems in Hokkaido, Japan, and have been managed using cultivars with multiple resistance. However, a limited number of adzuki bean varieties are resistant to these pathogens because of the frequent appearance of new races; thus, new sources of resistance have been sought in related Vigna spp., particularly in the section Angulares to which adzuki bean belongs. An analysis of selected Vigna accessions (JP81231 to JP235420) conserved in the Genebank of the National Institute of Agrobiological Sciences, Tsukuba, Japan, revealed wide variation in resistance spectra and resistance combinations, and eight disease response groups (A to H) were identified. Four of eight were newly detected response groups, suggesting the existence of novel resistance genes. Of 252 accessions from 26 species, 28 accessions in Vigna angularis var. nipponensis, V. hirtella, V. minima, and V. tenuicaulis (section Angulares) from group D, which were cross-compatible with adzuki bean, are expected to be potential sources of multiple resistance genes. They were resistant to all races of BSR and AFW pathogens examined.

Mapping of quantitative trait loci for phytic acid and phosphorus contents in seed and seedling of mungbean (Vigna radiata (L.) Wilczek)

Phytic acid (PA) is the storage form of phosphorus (P) in seeds and plays an important role in the nutritional quality of food crops. There is little information on the genetics of seed and seedling PA in mungbean [Vigna radiata (L.) Wilczek]. Quantitative trait loci (QTL) were identified for phytic acid P (PAP), total P (TP), and inorganic P (IP) in mungbean seeds and seedlings, and for flowering, maturity and seed weight, in an F(2) population developed from a cross between low PAP cultivated mungbean (V1725BG) and high PAP wild mungbean (AusTRCF321925). Seven QTLs were detected for P compounds in seed; two for PAP, four for IP and one for TP. Six QTLs were identified for P compounds in seedling; three for PAP, two for TP and one for IP. Only one QTL co-localized between P compounds in seed and seedling suggesting that low PAP seed and low PAP seedling must be selected for at different QTLs. Seed PAP and TP were positively correlated with days to flowering and maturity, indicating the importance of plant phenology to seed P content.

A model to predict the frequency of integration of fitness-related QTLs from cultivated to wild soybean

With the proliferation of genetically modified (GM) products and the almost exponential growth of land use for GM crops, there is a growing need to develop quantitative approaches to estimating the risk of escape of transgenes into wild populations of crop relatives by natural hybridization. We assessed the risk of transgene escape by constructing a population genetic model based on information on fitness-related QTLs obtained from an F (2) population of wild soybean G. soja × cultivated soybean Glycine max. Simulation started with ten F (1) and 990 wild soybeans reproducing by selfing or outcrossing. Seed production was determined from the genetic effects of two QTLs for number of seeds (SN). Each seed survived winter according to the maternal genotype at three QTLs for winter survival (WS). We assumed that one neutral transgene was inserted at various sites and calculated its extinction rate. The presence of G. max alleles at SN and WS QTLs significantly decreased the probability of introgression of the neutral transgene at all insertion sites equally. The presence of G. max alleles at WS QTLs lowered the risk more than their presence at SN QTLs. Although most model studies have concentrated only on genotypic effects of transgenes, we show that the presence of fitness-related domestication genes has a large effect on the risk of transgene escape. Our model offers the advantage of considering the effects of both domestication genes and a transgene, and they can be widely applied to other wild × crop relative complexes.

An SSR-based linkage map of yardlong bean (Vigna unguiculata (L.) Walp. subsp. unguiculata Sesquipedalis Group) and QTL analysis of pod length

Yardlong bean (Vigna unguiculata (L.) Walp. subsp. unguiculata Sesquipedalis Group) (2n = 2x = 22) is one of the most important vegetable legumes of Asia. The objectives of this study were to develop a genetic linkage map of yardlong bean using SSR makers from related Vigna species and to identify QTLs for pod length. The map was constructed from 226 simple sequence repeat (SSR) markers from cowpea (Vigna unguiculata (L.) Walp. subsp. unguiculata Unguiculata Group), azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi), and mungbean (Vigna radiata (L.) Wilczek) in a BC1F1 ((JP81610 × TVnu457) × JP81610) population derived from the cross between yardlong bean accession JP81610 and wild cowpea (Vigna unguiculata subsp. unguiculata var. spontanea) accession TVnu457. The markers were clustered into 11 linkage groups (LGs) spanning 852.4 cM in total length with a mean distance between adjacent markers of 3.96 cM. All markers on LG11 showed segregation distortion towards the homozygous yardlong bean JP81610 genotype. The markers on LG11 were also distorted in the rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) map, suggesting the presence of common segregation distortion factors in Vigna species on this LG. One major and six minor QTLs were identified for pod length variation between yardlong bean and wild cowpea. Using flanking markers, six of the seven QTLs were confirmed in an F2 population of JP81610 × TVnu457. The molecular linkage map developed and markers linked to pod length QTLs would be potentially useful for yardlong bean and cowpea breeding.

Genetic diversity of the Bambara groundnut (Vigna subterranea (L.) Verdc.) as assessed by SSR markers

Bambara groundnut ( Vigna subterranea (L.) Verdc.) is an important African legume crop. In this study, a collection consisting of 240 accessions was analyzed using 22 simple sequence repeat (SSR) markers. In total, 166 alleles were detected, with a mean of 7.59 alleles per locus. Allelic and gene diversities were higher in the west African and Cameroon/Nigeria regions with 6.68 and 6.18 alleles per locus, and 0.601 and 0.571, respectively. The genetic distance showed high similarity between west African and Cameroon/Nigeria accessions. Principal coordinate analyses and neighbor-joining analysis consistently revealed that the majority of west African accessions were grouped with Cameroon/Nigeria accessions, but they were differentiated from east African, central African, and southeast Asian accessions. Population structure analysis showed that two subpopulations existed, and most of the east African accessions were restricted to one subpopulation with some Cameroon/Nigeria accessions, whereas most of the west African accessions were associated with most of the Cameroon/Nigeria accessions in the other subpopulation. Comparison with SSR analysis of other Vigna cultigens, i.e., cultivated azuki bean ( Vigna angularis ) and mungbean ( Vigna radiata ), reveals that the mean gene diversity of Bambara groundnut was lower than azuki bean but higher than mungbean.

The genetically unstable dwarf locus in azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi)

We characterized a spontaneous dwarf mutant showing extremely short internodes and dark green leaves originating from azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi) cultivar "Erimo-shouzu." F(1) plants of 3 cross combinations between the dwarf mutant and several representative wild-type plants, Erimo-shouzu, V. angularis accession Acc2265 and wild relative V. riukiuensis accession Acc2482, supported the dwarf genotype being recessive. In a total of 3328 F(2) progeny of these 3 crosses, 65 dwarfs (2.0%) and 5 chimeric dwarfs (0.2%) segregated and the remainder were wild-type plants (97.8%). In F(3) progeny derived from self-pollinated dwarf F(2) plants, we observed wild type (54.3%), dwarf (39.1%), and chimeric dwarf (6.5%) plants. Two types of chimeric plants were observed: dwarf branches on the axils of wild-type plant stems and wild-type branches on the axils of dwarf stems. In 21 dwarf F(2) plants, the dwarf trait cosegregated with simple sequence repeat marker CEDG154 on chromosome 4. Conversely, homozygote F(2) plants at this chromosomal segment from the dwarf mutant frequently (>90%) expressed the wild-type phenotype. We concluded that the dwarf phenotype is mitotically and meiotically inheritable and controlled by a single genetically unstable locus, designated Azuki Dwarf1 (AD1), which converts between 2 phenotypic states bidirectionally.

A SNP and SSR based genetic map of asparagus bean (Vigna. unguiculata ssp. sesquipedialis) and comparison with the broader species

Asparagus bean (Vigna. unguiculata ssp. sesquipedialis) is a distinctive subspecies of cowpea [Vigna. unguiculata (L.) Walp.] that apparently originated in East Asia and is characterized by extremely long and thin pods and an aggressive climbing growth habit. The crop is widely cultivated throughout Asia for the production of immature pods known as 'long beans' or 'asparagus beans'. While the genome of cowpea ssp. unguiculata has been characterized recently by high-density genetic mapping and partial sequencing, little is known about the genome of asparagus bean. We report here the first genetic map of asparagus bean based on SNP and SSR markers. The current map consists of 375 loci mapped onto 11 linkage groups (LGs), with 191 loci detected by SNP markers and 184 loci by SSR markers. The overall map length is 745 cM, with an average marker distance of 1.98 cM. There are four high marker-density blocks distributed on three LGs and three regions of segregation distortion (SDRs) identified on two other LGs, two of which co-locate in chromosomal regions syntenic to SDRs in soybean. Synteny between asparagus bean and the model legume Lotus. japonica was also established. This work provides the basis for mapping and functional analysis of genes/QTLs of particular interest in asparagus bean, as well as for comparative genomics study of cowpea at the subspecies level.

The genetics of domestication of rice bean, Vigna umbellata

BACKGROUND AND AIMS

The Asian genus Vigna, to which four cultivated species (rice bean, azuki bean, mung bean and black gram) belong, is suitable for comparative genomics. The aims were to construct a genetic linkage map of rice bean, to identify the genomic regions associated with domestication in rice bean, and to compare these regions with those in azuki bean.

METHODS

A genetic linkage map was constructed by using simple sequence repeat and amplified fragment length polymorphism markers in the BC(1)F(1) population derived from a cross between cultivated and wild rice bean. Using this map, 31 domestication-related traits were dissected into quantitative trait loci (QTLs). The genetic linkage map and QTLs of rice bean were compared with those of azuki bean.

KEY RESULTS

A total of 326 markers converged into 11 linkage groups (LGs), corresponding to the haploid number of rice bean chromosomes. The domestication-related traits in rice bean associated with a few major QTLs distributed as clusters on LGs 2, 4 and 7. A high level of co-linearity in marker order between the rice bean and azuki bean linkage maps was observed. Major QTLs in rice bean were found on LG4, whereas major QTLs in azuki bean were found on LG9.

CONCLUSIONS

This is the first report of a genetic linkage map and QTLs for domestication-related traits in rice bean. The inheritance of domestication-related traits was so simple that a few major QTLs explained the phenotypic variation between cultivated and wild rice bean. The high level of genomic synteny between rice bean and azuki bean facilitates QTL comparison between species. These results provide a genetic foundation for improvement of rice bean; interchange of major QTLs between rice bean and azuki bean might be useful for broadening the genetic variation of both species.

NIASGBdb: NIAS Genebank databases for genetic resources and plant disease information

The National Institute of Agrobiological Sciences (NIAS) is implementing the NIAS Genebank Project for conservation and promotion of agrobiological genetic resources to contribute to the development and utilization of agriculture and agricultural products. The project’s databases (NIASGBdb; http://www.gene.affrc.go.jp/databases_en.php) consist of a genetic resource database and a plant diseases database, linked by a web retrieval database. The genetic resources database has plant and microorganism search systems to provide information on research materials, including passport and evaluation data for genetic resources with the desired properties. To facilitate genetic diversity research, several NIAS Core Collections have been developed. The NIAS Rice (Oryza sativa) Core Collection of Japanese Landraces contains information on simple sequence repeat (SSR) polymorphisms. SSR marker information for azuki bean (Vigna angularis) and black gram (V. mungo) and DNA sequence data from some selected Japanese strains of the genus Fusarium are also available. A database of plant diseases in Japan has been developed based on the listing of common names of plant diseases compiled by the Phytopathological Society of Japan. Relevant plant and microorganism genetic resources are associated with the plant disease names by the web retrieval database and can be obtained from the NIAS Genebank for research or educational purposes.

Extension of the core map of common bean with EST-SSR, RGA, AFLP, and putative functional markers

Microsatellites and gene-derived markers are still underrepresented in the core molecular linkage map of common bean compared to other types of markers. In order to increase the density of the core map, a set of new markers were developed and mapped onto the RIL population derived from the 'BAT93' x 'Jalo EEP558' cross. The EST-SSR markers were first characterized using a set of 24 bean inbred lines. On average, the polymorphism information content was 0.40 and the mean number of alleles per locus was 2.7. In addition, AFLP and RGA markers based on the NBS-profiling method were developed and a subset of the mapped RGA was sequenced. With the integration of 282 new markers into the common bean core map, we were able to place markers with putative known function in some existing gaps including regions with QTL for resistance to anthracnose and rust. The distribution of the markers over 11 linkage groups is discussed and a newer version of the common bean core linkage map is proposed.

Genetic diversity analysis in blackgram (Vigna mungo (L.) Hepper) using AFLP and transferable microsatellite markers from azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi)

Genetic diversity in 20 elite blackgram (Vigna mungo (L.) Hepper) genotypes was studied using microsatellite and AFLP markers. Thirty-six microsatellite markers from azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi) were successfully amplified across the 20 blackgram genotypes and 33 microsatellite markers showed polymorphism. A total of 137 microsatellite alleles were generated with an average of 4.1 alleles per locus. The number of alleles ranged from two to nine and the polymorphic information content value for the microsatellite markers varied from 0.10 to 0.87 with an average of 0.49. Microsatellite markers were highly informative and a combination of only three microsatellite markers (CEDG264, CEDG173, and CEDG044) was sufficient to discriminate all 20 blackgram genotypes. In the case of AFLP, 11 primer pairs generated 324 polymorphic marker fragments. The polymorphic information content values for AFLP primer combinations ranged from 0.21 to 0.34 with an average of 0.29. Similarity measures and clustering analyses were made using microsatellite and AFLP data separately. The resulting dendrograms distributed the 20 blackgram genotypes into five main clusters. The dendrograms were comparable with each other with the Mantel test between the cophenetic matrices of microsatellite data and AFLP data showing moderate correlation (r = 0.64). The results of the principal components analysis were well congruent with the dendrograms. In the dendrograms as well as in the principal components analyses, genotype Trombay wild (Vigna mungo var. silvestris) was placed separately from rest of the genotypes. This study demonstrated that the azuki bean microsatellite markers are highly polymorphic and informative and can be successfully used for genome analysis in blackgram. Results indicate that sufficient variability is present in the blackgram genotypes and would be helpful in the selection of suitable parents for breeding purposes and gene mapping studies.

Genetic diversity of the azuki bean (Vigna angularis (Willd.) Ohwi & Ohashi) gene pool as assessed by SSR markers

To facilitate the wider use of genetic resources including newly collected cultivated and wild azuki bean germplasm, the genetic diversity of the azuki bean complex, based on 13 simple sequence repeat (SSR) primers, was evaluated and a core collection was developed using 616 accessions originating from 8 Asian countries. Wild germplasm from Japan was highly diverse and represented much of the allelic variation found in cultivated germplasm. The SSR results together with recent archaeobotanical evidence support the view that Japan is one center of domestication of azuki bean, at least for the northeast Asian azuki bean. Cultivated azuki beans from China, Korea, and Japan were the most diverse and were genetically distinct from each other, suggesting a long and relatively isolated history of cultivation in each country. Cultivated azuki beans from eastern Nepal and Bhutan were similar to each other and quite distinct from others. For two primers, most eastern Nepalese and Bhutanese cultivated accessions had null alleles. In addition, wild accessions from the Yangtze River region of China and the Himalayan region had a null allele for one or the other of these primers. Whether the distinct diversity of azuki bean in the Himalayan region is due to introgression or separate domestication events requires further study. In contrast, western Nepalese azuki beans showed an SSR profile similar to that of Chinese azuki beans. The genetic distinctness of cultivated azuki beans from Vietnam has been revealed for the first time. The specific alleles indicate that Vietnamese azuki beans have been cultivated in isolation from Chinese azuki beans for a long time. Wild germplasm from the Himalayan region showed the highest level of variation. Based on the results, Himalayan germplasm could be considered a novel gene source for azuki bean breeding. A comparison with mungbean SSR analysis revealed that the mean gene diversity of cultivated azuki bean (0.74) was much higher than that of cultivated mungbean (0.41). The reduction in gene diversity due to domestication, the domestication bottleneck, in azuki bean is not strong compared with that in mungbean.

Construction of a genetic linkage map of black gram, Vigna mungo (L.) Hepper, based on molecular markers and comparative studies

A genetic linkage map of black gram, Vigna mungo (L.) Hepper, was constructed with 428 molecular markers using an F9 recombinant inbred population of 104 individuals. The population was derived from an inter-subspecific cross between a black gram cultivar, TU94-2, and a wild genotype, V. mungo var. silvestris. The linkage analysis at a LOD score of 5.0 distributed all 428 markers (254 AFLP, 47 SSR, 86 RAPD, and 41 ISSR) into 11 linkage groups. The map spanned a total distance of 865.1 cM with an average marker density of 2 cM. The largest linkage group spanned 115 cM and the smallest linkage group was of 44.9 cM. The number of markers per linkage group ranged from 11 to 86 and the average distance between markers varied from 1.1 to 5.6 cM. Comparison of the map with other published azuki bean and black gram maps showed high colinearity of markers, with some inversions. The current map is the most saturated map for black gram to date and will provide a useful tool for identification of QTLs and for marker-assisted selection of agronomically important characters in black gram.

Mapping of quantitative trait loci for a new source of resistance to bruchids in the wild species Vigna nepalensis Tateishi & Maxted (Vigna subgenus Ceratotropis)

Azuki bean breeders have long been interested in producing azuki bean [Vigna angularis (Willd.) Ohwi & Ohashi] varieties with bruchid resistance. A new bruchid (Callosobruchus spp.) resistance source was found in V. nepalensis Tateishi & Maxted, a species that is cross compatible with azuki bean. Quantitative trait loci (QTLs) analysis for resistance to C. chinensis (L.) and C. maculatus (F.) was conducted using F2 (V. nepalensisxV. angularis) and BC1F1 [(V. nepalensisxV. angularis)xV. angularis] populations derived from crosses between the bruchid resistant species V. nepalensis and bruchid susceptible species V. angularis. Resistance was measured using two traits, percentage of seeds damaged by bruchids and the time taken for adult bruchids to emerge from seeds. Based on the results from both populations seven QTLs were detected for bruchid resistance; five QTLs for resistance to C. chinensis and two QTLs for resistance to C. maculatus. The different locations found for some resistance QTL to the two bruchid species suggests different resistance mechanisms. QTLs on linkage group (LG) 1 and LG2 for bruchid resistance to C. chinensis co-localized with seed size QTLs suggesting that incremental increase in seed size accompanied susceptibility to C. chinensis. Based on linked markers the QTL on these two linkage groups appear to be the same as previously reported in other Asian Vigna. However, several other QTLs were newly detected including one on LG4 that appears unrelated to seed size. Transfer of these new sources of bruchid resistance from V. nepalensis to azuki bean will be aided by the progress being made in azuki genome mapping.

The genetics of domestication of the azuki bean (Vigna angularis)

Genetic differences between azuki bean (Vigna angularis var. angularis) and its presumed wild ancestor (V. angularis var. nipponensis) were resolved into QTL for traits associated with adaptation to their respective distinct habits. A genetic linkage map constructed using progenies from a cross between Japanese cultivated and wild azuki beans covers 92.8% of the standard azuki bean linkage map. A reciprocal translocation between cultivated and wild azuki bean parents was identified on the basis of the linkage map having a pseudolinkage group and clustering of seed productivity-related QTL with large effect near the presumed breakpoints. In total, 162 QTL were identified for 46 domestication-related traits. Domestication of azuki bean has involved a trade-off between seed number and seed size: fewer but longer pods and fewer but larger seeds on plants with shorter stature in cultivated azuki bean being at the expense of overall seed yield. Genes found related to germination and flowering time in cultivated azuki bean may confer a selective advantage to the hybrid derivatives under some ecological conditions and may explain why azuki bean has evolved as a crop complex in Japan.

Genome dissection of traits related to domestication in azuki bean (Vigna angularis) and comparison with other warm-season legumes

BACKGROUND

The objective of this study was to dissect into quantitative trait loci (QTLs) the large morphological and physiological differences between cultivated azuki bean (Vigna angularis) and a wild relative and to infer the commonalities of the QTLs for domestication-related traits across the Asian Vigna and with other warm-season legumes.

METHODS

Two linkage maps, for the BC(1)F(1) and F(2) populations, respectively, from the same cross between azuki bean and V. nepalensis were developed. Using these linkage maps QTLs for 33 domestication-related traits were analysed and mapped. The location of mapped QTLs was compared with locations of similar QTLs in other warm-season legumes.

KEY RESULTS

QTLs were detected for seed-, pod-, stem- and leaf-related traits. Most traits were controlled by between two and nine QTLs but several traits, such as pod dehiscence, were controlled by single genes. QTLs for domestication-related traits were restricted to particular regions of the azuki bean genome, especially linkage groups 1, 2, 4, 7 and 9. Linkage groups 1 and 2 had QTLs for a suite of traits including pod size, germination, seed size and lower stem length. QTLs on linkage groups 7 and 9 were associated with upper stem length, maximum leaf size and pod and seed size. Pleiotropy or close linkage of genes for domestication-related traits is suggested in these regions. While some QTLs are common to azuki bean and other warm-season legumes, many are recorded for the first time in azuki bean.

CONCLUSIONS

QTLs for a large number of domestication-related traits have been mapped for the first time in azuki bean. QTLs with unexpected effect and new QTLs for traits such as seed size have been found. The results provide a foundation that will be useful for improvement of azuki bean and related legumes.