Transcriptome-wide association mapping provides insights into the genetic basis and candidate genes governing flowering, maturity and seed weight in rice bean (Vigna umbellata)

BACKGROUND

Rice bean (Vigna umbellata), an underrated legume, adapts to diverse climatic conditions with the potential to support food and nutritional security worldwide. It is used as a vegetable, minor food crop and a fodder crop, being a rich source of proteins, minerals, and essential fatty acids. However, little effort has been made to decipher the genetic and molecular basis of various useful traits in this crop. Therefore, we considered three economically important traits i.e., flowering, maturity and seed weight of rice bean and identified the associated candidate genes employing an associative transcriptomics approach on 100 diverse genotypes out of 1800 evaluated rice bean accessions from the Indian National Genebank.

RESULTS

The transcriptomics-based genotyping of one-hundred diverse rice bean cultivars followed by pre-processing of genotypic data resulted in 49,271 filtered markers. The STRUCTURE, PCA and Neighbor-Joining clustering of 100 genotypes revealed three putative sub-populations. The marker-trait association analysis involving various genome-wide association study (GWAS) models revealed significant association of 82 markers on 48 transcripts for flowering, 26 markers on 22 transcripts for maturity and 22 markers on 21 transcripts for seed weight. The transcript annotation provided information on the putative candidate genes for the considered traits. The candidate genes identified for flowering include HSC80, P-II PsbX, phospholipid-transporting-ATPase-9, pectin-acetylesterase-8 and E3-ubiquitin-protein-ligase-RHG1A. Further, the WRKY1 and DEAD-box-RH27 were found to be associated with seed weight. Furthermore, the associations of PIF3 and pentatricopeptide-repeat-containing-gene with maturity and seed weight, and aldo-keto-reductase with flowering and maturity were revealed.

CONCLUSION

This study offers insights into the genetic basis of key agronomic traits in rice bean, including flowering, maturity, and seed weight. The identified markers and associated candidate genes provide valuable resources for future exploration and targeted breeding, aiming to enhance the agronomic performance of rice bean cultivars. Notably, this research represents the first transcriptome-wide association study in pulse crop, uncovering the candidate genes for agronomically useful traits.

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.

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.

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.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development–related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development–related traits in ricebeans. A total of 51 genes encoding SCF^TIR1/AFB, Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean–derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Rice bean: a lesser known pulse with well-recognized potential

Required genetic resources for the improvement of agronomic, nutritional and economic value of rice bean are available in the world collection. International cooperative effort is required to utilize and conserve them. Rice bean [Vigna umbellata (Thunb.) Ohwi and Ohashi], a lesser known pulse among the Asiatic Vigna, has long been considered as a food security crop of small and marginal farmers of Southeast Asia. Considered as a nutritionally rich food and fodder, it is also a source of genes for biotic and abiotic stress tolerance including drought, soil acidity and storage pest. Although it spread from its centre of domestication in the Indo-China region to other parts around the world, it never became an important crop anywhere probably because of agronomic disadvantages. Crop improvement for determinate nature, good yield, less variable seed colour, pleasant organoleptic properties and lower antinutrients is required. Scanning of scientific literature indicates that genetic resources with desirable agronomic and nutritional traits exist within the current collection but are spread across countries. Genomic studies in the species indicate that except for insect resistance and aluminium toxicity tolerance, not much attention has been paid to decipher and utilize other stress tolerance and nutritional quality traits. Collaborative efforts towards improving farming, food, trade value and off-farm conservation of rice bean would not only help marginal farmers but will also help to preserve the yet to be explored genomic resources available in this sturdy pulse.

Functional analysis of the promoter of the molt-inhibiting hormone (mih) gene in mud crab Scylla paramamosain

In this study, the 5'-flanking region of molt-inhibiting hormone (MIH) gene was cloned by Tail-PCR. It is 2024 bp starting from the translation initiation site, and 1818 bp starting from the predicted transcription start site. Forecast analysis results by the bioinformatics software showed that the transcription start site is located at 207 bp upstream of the start codon ATG, and TATA box is located at 240 bp upstream of the start codon ATG. Potential transcription factor binding sites include Sp1, NF-1, Oct-1, Sox-2, RAP1, and so on. There are two CpG islands, located at -25- +183 bp and -1451- -1316 bp respectively. The transfection results of luciferase reporter constructs showed that the core promoter region was located in the fragment -308 bp to -26 bp. NF-kappaB and RAP1 were essential for mih basal transcriptional activity. There are three kinds of polymorphism CA in the 5'-flanking sequence, and they can influence mih promoter activity. These findings provide a genetic foundation of the further research of mih transcription regulation.

Development and characterization of a new set of genomic microsatellite markers in rice bean (Vigna umbellata (Thunb.) Ohwi and Ohashi) and their utilization in genetic diversity analysis of collections from North East India

Rice bean [Vigna umbellate (Thumb.) Ohwi and Ohashi] is an underutilized crop believed to be domesticated in the Myanmar-Thailand region of Asia. In India, rice bean is mainly cultivated in the North-Eastern Hills, which is a hotspot for biological diversity. A 5' anchored PCR was used to develop microsatellite markers in rice bean. Twenty-eight specific primer pairs were designed and employed to characterize sixty five ricebean accessions collected from North East India. A total of 179 alleles were amplified with an average of 6.393 alleles per locus. The gene diversity was high (mean 0.534) in the accessions collected from Darjeeling, Nagaland and Manipur, which are bordering areas with East Nepal and Myanmar, respectively. Exceptionally high outcrossing rate was observed in the entire population. Population structure analysis identified three distinct clusters in which accessions collected from areas bordering Myanmar and East Nepal grouped separately. Using a combination of STRUCTURE and Principal Coordinate Analysis, relative affinity of the intermediate accessions could be established. However, differences in allelic counts among populations were non-significant. The results showed that there is a high level of genetic diversity within the accessions, with high outcrossing rate.

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

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

Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data

Rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) is a warm season annual legume mainly grown in East Asia. Only scarce genomic resources are currently available for this legume crop species and no simple sequence repeat (SSR) markers have been specifically developed for rice bean yet. In this study, approximately 26 million high quality cDNA sequence reads were obtained from rice bean using Illumina paired-end sequencing technology and assembled into 71,929 unigenes with an average length of 986 bp. Of these unigenes, 38,840 (33.2%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases. Furthermore, 30,170 (76.3%) could be classified into gene ontology categories, 25,451 (64.4%) into Swiss-Prot categories and 21,982 (55.6%) into KOG database categories (E-value < 1.0E-5). A total of 9,301 (23.5%) were mapped onto 118 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. A total of 3,011 genic SSRs were identified as potential molecular markers. AG/CT (30.3%), AAG/CTT (8.1%) and AGAA/TTCT (20.0%) are the three main repeat motifs. A total of 300 SSR loci were randomly selected for validation by using PCR amplification. Of these loci, 23 primer pairs were polymorphic among 32 rice bean accessions. A UPGMA dendrogram revealed three major clusters among 32 rice bean accessions. The large number of SSR-containing sequences and genic SSRs in this study will be valuable for the construction of high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

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.

A Homoploid Hybrid Between Wild Vigna Species Found in a Limestone Karst

Genus Vigna comprise several domesticated species including cowpea and mungbean, and diverse wild species. We found an introgressive hybrid population derived from two wild species, Vigna umbellata and Vigna exilis, in Ratchaburi district, Thailand. The hybrid was morphologically similar to V. umbellata but habituated in a limestone rock mountain, which is usually dominated by V. exilis. Analyzing simple sequence repeat loci indicated the hybrid has undergone at least one round of backcross by V. umbellata. We found the hybrid acquired vigorous growth from V. umbellata and drought tolerance plus early flowering from V. exilis, and thus has taken over some habitats of V. exilis in limestone karsts. Given the wide crossability of V. umbellata, the hybrid can be a valuable genetic resource to improve drought tolerance of some domesticated species.

Genetic diversity of the black gram [Vigna mungo (L.) Hepper] gene pool as revealed by SSR markers

In this study, 520 cultivated and 14 wild accessions of black gram (Vigna mungo (L.) Hepper) were assessed for diversity using 22 SSR markers. Totally, 199 alleles were detected with a mean of 9.05 alleles per locus. Wild black gram showed higher gene diversity than cultivated black gram. Gene diversity of cultivated accessions among regions was comparable, while allelic richness of South Asia was higher than that of other regions. 78.67% of the wild gene diversity presented in cultivated accessions, indicating that the domestication bottleneck effect in black gram is relatively low. Genetic distance analysis revealed that cultivated black gram was more closely related to wild black gram from South Asia than that from Southeast Asia. STRUCTURE, principal coordinate and neighbor-joining analyses consistently revealed that 534 black gram accessions were grouped into three major subpopulations. The analyses also revealed that cultivated black gram from South Asia was genetically distinct from that from West Asia. Comparison by SSR analysis with other closely related Vigna species, including mungbean, azuki bean, and rice bean, revealed that level of gene diversity of black gram is comparable to that of mungbean and rice bean but lower than that of 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.