Detection of QTLs associated with mungbean yellow mosaic virus (MYMV) resistance using the interspecific cross of Vigna radiata × Vigna umbellata

Multi-GWAS reveals significant genomic regions for Mungbean yellow mosaic India virus resistance in urdbean (Vigna mungo (L.) across multiple environments

KEY MESSAGE

Unraveling genetic markers for MYMIV resistance in urdbean, with 8 high-confidence marker-trait associations identified across diverse environments, provides crucial insights for combating MYMIV disease, informing future breeding strategies. Globally, yellow mosaic disease (YMD) causes significant yield losses, reaching up to 100% in favorable environments within major urdbean cultivating regions. The introgression of genomic regions conferring resistance into urdbean cultivars is crucial for combating YMD, including resistance against mungbean yellow mosaic India virus (MYMIV). To uncover the genetic basis of MYMIV resistance, we conducted a genome-wide association study (GWAS) using three multi-locus models in 100 diverse urdbean genotypes cultivated across six individual and two combined environments. Leveraging 4538 high-quality single nucleotide polymorphism (SNP) markers, we identified 28 unique significant marker-trait associations (MTAs) for MYMIV resistance, with 8 MTAs considered of high confidence due to detection across multiple GWAS models and/or environments. Notably, 4 out of 28 MTAs were found in proximity to previously reported genomic regions associated with MYMIV resistance in urdbean and mungbean, strengthening our findings and indicating consistent genomic regions for MYMIV resistance. Among the eight highly significant MTAs, one localized on chromosome 6 adjacent to previously identified quantitative trait loci for MYMIV resistance, while the remaining seven were novel. These MTAs contain several genes implicated in disease resistance, including four common ones consistently found across all eight MTAs: receptor-like serine-threonine kinases, E3 ubiquitin-protein ligase, pentatricopeptide repeat, and ankyrin repeats. Previous studies have linked these genes to defense against viral infections across different crops, suggesting their potential for further basic research involving cloning and utilization in breeding programs. This study represents the first GWAS investigation aimed at identifying resistance against MYMIV in urdbean germplasm.

Identification of quantitative trait loci (QTLs) regulating leaf SPAD value and trichome density in mungbean (Vigna radiata L.) using genotyping-by-sequencing (GBS) approach

Quantitative trait loci (QTL) mapping is used for the precise localization of genomic regions regulating various traits in plants. Two major QTLs regulating Soil Plant Analysis Development (SPAD) value (qSPAD-7-1) and trichome density (qTric-7-2) in mungbean were identified using recombinant inbred line (RIL) populations (PMR-1×Pusa Baisakhi) on chromosome 7. Functional analysis of QTL region identified 35 candidate genes for SPAD value (16 No) and trichome (19 No) traits. The candidate genes regulating trichome density on the dorsal leaf surface of the mungbean include VRADI07G24840, VRADI07G17780, and VRADI07G15650, which encodes for ZFP6, TFs bHLH DNA-binding superfamily protein, and MYB102, respectively. Also, candidate genes having vital roles in chlorophyll biosynthesis are VRADIO7G29860, VRADIO7G29450, and VRADIO7G28520, which encodes for s-adenosyl-L-methionine, FTSHI1 protein, and CRS2-associated factor, respectively. The findings unfolded the opportunity for the development of customized genotypes having high SPAD value and high trichome density having a possible role in yield and mungbean yellow vein mosaic India virus (MYMIV) resistance in mungbean.

Major viral diseases in grain legumes: designing disease resistant legumes from plant breeding and OMICS integration

Grain legumes play a crucial role in human nutrition and as a staple crop for low-income farmers in developing and underdeveloped nations, contributing to overall food security and agroecosystem services. Viral diseases are major biotic stresses that severely challenge global grain legume production. In this review, we discuss how exploring naturally resistant grain legume genotypes within germplasm, landraces, and crop wild relatives could be used as promising, economically viable, and eco-environmentally friendly solution to reduce yield losses. Studies based on Mendelian and classical genetics have enhanced our understanding of key genetic determinants that govern resistance to various viral diseases in grain legumes. Recent advances in molecular marker technology and genomic resources have enabled us to identify genomic regions controlling viral disease resistance in various grain legumes using techniques such as QTL mapping, genome-wide association studies, whole-genome resequencing, pangenome and 'omics' approaches. These comprehensive genomic resources have expedited the adoption of genomics-assisted breeding for developing virus-resistant grain legumes. Concurrently, progress in functional genomics, especially transcriptomics, has helped unravel underlying candidate gene(s) and their roles in viral disease resistance in legumes. This review also examines the progress in genetic engineering-based strategies, including RNA interference, and the potential of synthetic biology techniques, such as synthetic promoters and synthetic transcription factors, for creating viral-resistant grain legumes. It also elaborates on the prospects and limitations of cutting-edge breeding technologies and emerging biotechnological tools (e.g., genomic selection, rapid generation advances, and CRISPR/Cas9-based genome editing tool) in developing virus-disease-resistant grain legumes to ensure global food security.

Genetic basis of maize kernel protein content revealed by high-density bin mapping using recombinant inbred lines

Maize with a high kernel protein content (PC) is desirable for human food and livestock fodder. However, improvements in its PC have been hampered by a lack of desirable molecular markers. To identify quantitative trait loci (QTL) and candidate genes for kernel PC, we employed a genotyping-by-sequencing strategy to construct a high-resolution linkage map with 6,433 bin markers for 275 recombinant inbred lines (RILs) derived from a high-PC female Ji846 and low-PC male Ye3189. The total genetic distance covered by the linkage map was 2180.93 cM, and the average distance between adjacent markers was 0.32 cM, with a physical distance of approximately 0.37 Mb. Using this linkage map, 11 QTLs affecting kernel PC were identified, including qPC7 and qPC2-2, which were identified in at least two environments. For the qPC2-2 locus, a marker named IndelPC2-2 was developed with closely linked polymorphisms in both parents, and when tested in 30 high and 30 low PC inbred lines, it showed significant differences (P = 1.9E-03). To identify the candidate genes for this locus, transcriptome sequencing data and PC best linear unbiased estimates (BLUE) for 348 inbred lines were combined, and the expression levels of the four genes were correlated with PC. Among the four genes, Zm00001d002625, which encodes an S-adenosyl-L-methionine-dependent methyltransferase superfamily protein, showed significantly different expression levels between two RIL parents in the endosperm and is speculated to be a potential candidate gene for qPC2-2. This study will contribute to further research on the mechanisms underlying the regulation of maize PC, while also providing a genetic basis for marker-assisted selection in the future.

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.

Consensus genetic map and QTLs for morphological and agronomical traits in mungbean using DArT markers

Mungbean is an important but understudied food legume compared with other major grain crops. Genetic studies through development of high-through put markers, linkage map construction and QTL analysis can accelerate and improve the efficiency of mining for genes for breeding in this crop. This study used four mungbean F₅ recombinant inbred lines (RILs) from crosses of two wild types (ACC 1, ACC 87) and two cultivars (Berken, Kiloga) and DArT markers to construct individual and consensus linkage maps and to identify QTLs associated with 54 traits in mungbean. The number of polymorphic DArT markers identified among the four RIL populations varied from 1062 to 2013. The individual maps covered the lengths of 629.7-883.5 cM, comprising 672-981 DArT markers and 15-19 linkage groups (LG) with average distance between markers of 0.9-1.2 cM. The consensus map had the total length of 795.3 cM, comprising 1539 DArT markers and resolved 11 LGs with an average inter-marker distance of 0.65 cM. Sixty-two QTLs were identified for 39 traits across 10 LGs of the consensus map. Major QTLs were identified for two special traits, late flowering inherited from ACC 1 (6 QTLs, PVE of 11.2-29.9%) and perenniality inherited from ACC 87 (3 QTLs, PVE of 17.4-22.6%) in separate population analysis. Number of congruent QTLs across four mungbean populations and the consensus map was 18 for 13 traits. These results illustrated the high efficiency of DArT marker application in mungbean genetic dissection and suggested the future potential employment of identified QTLs for mungbean improvement.

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.

Dynamic Transcriptome Profiling of Mungbean Genotypes Unveil the Genes Respond to the Infection of Mungbean Yellow Mosaic Virus

Yellow mosaic disease (YMD), incited by mungbean yellow mosaic virus (MYMV), is a primary viral disease that reduces mungbean production in South Asia, especially in India. There is no detailed knowledge regarding the genes and molecular mechanisms conferring resistance of mungbean to MYMV. Therefore, disclosing the genetic and molecular bases related to MYMV resistance helps to develop the mungbean genotypes with MYMV resistance. In this study, transcriptomes of mungbean genotypes, VGGRU-1 (resistant) and VRM (Gg) 1 (susceptible) infected with MYMV were compared to those of uninfected controls. The number of differentially expressed genes (DEGs) in the resistant and susceptible genotypes was 896 and 506, respectively. Among them, 275 DEGs were common between the resistant and susceptible genotypes. Functional annotation of DEGs revealed that the DEGs belonged to the following categories defense and pathogenesis, receptor-like kinases; serine/threonine protein kinases, hormone signaling, transcription factors, and chaperons, and secondary metabolites. Further, we have confirmed the expression pattern of several DEGs by quantitative real-time PCR (qRT-PCR) analysis. Collectively, the information obtained in this study unveils the new insights into characterizing the MYMV resistance and paved the way for breeding MYMV resistant mungbean in the future.

Identification of a Locus Controlling Compound Raceme Inflorescence in Mungbean [Vigna radiata (L.) R. Wilczek]

Mungbean [Vigna radiata (L.) R. Wilczek] produces a compound raceme inflorescence that branches into secondary inflorescences, which produce flowers. This architecture results in the less-domesticated traits of asynchronous pod maturity and multiple harvest times. This study identified the genetic factors responsible for the compound raceme of mungbean, providing a unique biological opportunity to improve simultaneous flowering. Using a recombinant inbred line (RIL) population derived from VC1973A, an elite cultivar with a compound raceme type, and IT208075, a natural mutant with a simple raceme type, a single locus that determined the inflorescence type was identified based on 1:1 segregation ratio in the F8 generation, and designated Comraceme. Linkage map analysis showed Comraceme was located on chromosome 4 within a marker interval spanning 520 kb and containing 64 genes. RILs carrying heterozygous fragments around Comraceme produced compound racemes, indicating this form was dominant to the simple raceme type. Quantitative trait loci related to plant architecture and inflorescence have been identified in genomic regions of soybean syntenic to Comraceme. In IT208075, 15 genes were present as distinct variants not observed in other landrace varieties or wild mungbean. These genes included Vradi04g00002481, a development-related gene encoding a B3 transcriptional factor. The upstream region of Vradi04g00002481 differed between lines producing the simple and compound types of raceme. Expression of Vradi04g00002481 was significantly lower at the early vegetative stage and higher at the early reproductive stage, in IT208075 than in VC1973A. Vradi04g00002481 was therefore likely to determine inflorescence type in mungbean. Although further study is required to determine the functional mechanism, this finding provides valuable genetic information for understanding the architecture of the compound raceme in mungbean.

Comparative RNA-Seq analysis unfolds a complex regulatory network imparting yellow mosaic disease resistance in mungbean [Vigna radiata (L.) R. Wilczek]

Yellow Mosaic Disease (YMD) in mungbean [Vigna radiata (L.) R. Wilczek] is one of the most damaging diseases in Asia. In the northern part of India, the YMD is caused by Mungbean Yellow Mosaic India Virus (MYMIV), while in southern India this is caused by Mungbean Yellow Mosaic Virus (MYMV). The molecular mechanism of YMD resistance in mungbean remains largely unknown. In this study, RNA-seq analysis was conducted between a resistant (PMR-1) and a susceptible (Pusa Vishal) mungbean genotype under infected and control conditions to understand the regulatory network operating between mungbean-YMV. Overall, 76.8 million raw reads could be generated in different treatment combinations, while mapping rate per library to the reference genome varied from 86.78% to 93.35%. The resistance to MYMIV showed a very complicated gene network, which begins with the production of general PAMPs (pathogen-associated molecular patterns), then activation of various signaling cascades like kinases, jasmonic acid (JA) and brassinosteroid (BR), and finally the expression of specific genes (like PR-proteins, virus resistance and R-gene proteins) leading to resistance response. The function of WRKY, NAC and MYB transcription factors in imparting the resistance against MYMIV could be established. The string analysis also revealed the role of proteins involved in kinase, viral movement and phytoene synthase activity in imparting YMD resistance. A set of novel stress-related EST-SSRs are also identified from the RNA-Seq data which may be used to find the linked genes/QTLs with the YMD resistance. Also, 11 defence-related transcripts could be validated through quantitative real-time PCR analysis. The identified gene networks have led to an insight about the defence mechanism operating against MYMIV infection in mungbean which will be of immense use to manage the YMD resistance in mungbean.

High-Density GBS-Based Genetic Linkage Map Construction and QTL Identification Associated With Yellow Mosaic Disease Resistance in Bitter Gourd (Momordica charantia L.)

Yellow mosaic disease (YMD) in bitter gourd (Momordica charantia) is a devastating disease that seriously affects its yield. Although there is currently no effective method to control the disease, breeding of resistant varieties is the most effective and economic option. Moreover, quantitative trait locus (QTL) associated with resistance to YMD has not yet been reported. With the objective of mapping YMD resistance in bitter gourd, the susceptible parent "Punjab-14" and the resistant parent "PAUBG-6" were crossed to obtain F₄ mapping population comprising 101 individuals. In the present study, the genotyping by sequencing (GBS) approach was used to develop the genetic linkage map. The map contained 3,144 single nucleotide polymorphism (SNP) markers, consisted of 15 linkage groups, and it spanned 2415.2 cM with an average marker distance of 0.7 cM. By adopting the artificial and field inoculation techniques, F_4:5 individuals were phenotyped for disease resistance in Nethouse (2019), Rainy (2019), and Spring season (2020). The QTL analysis using the genetic map and phenotyping data identified three QTLs qYMD.pau_3.1, qYMD.pau_4.1, and qYMD.pau_5.1 on chromosome 3, 4, and 5 respectively. Among these, qYMD.pau_3.1, qYMD.pau_4.1 QTLs were identified during the rainy season, explaining the 13.5 and 21.6% phenotypic variance respectively, whereas, during the spring season, qYMD.pau_4.1 and qYMD.pau_5.1 QTLs were observed with 17.5 and 22.1% phenotypic variance respectively. Only one QTL qYMD.pau_5.1 was identified for disease resistance under nethouse conditions with 15.6% phenotypic variance. To our knowledge, this is the first report on the identification of QTLs associated with YMD resistance in bitter gourd using SNP markers. The information generated in this study is very useful in the future for fine-mapping and marker-assisted selection for disease resistance.

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.

V. R, Tripathi K, Kumar RR, Aski M, Singh A, Roy A, Priti, Kumari N, Dasgupta U, Kumar A, Praveen S, Nair RM. Yellow Mosaic Disease (YMD) of Mungbean (Vigna radiata (L.) Wilczek): Current Status and Management Opportunities

Globally, yellow mosaic disease (YMD) remains a major constraint of mungbean production, and management of this deadly disease is still the biggest challenge. Thus, finding ways to manage YMD including development of varieties possessing resistance against mungbean yellow mosaic virus (MYMV) and mungbean yellow mosaic India virus (MYMIV) is a research priority for mungbean crop. Characterization of YMD resistance using various advanced molecular and biochemical approaches during plant-virus interactions has unfolded a comprehensive network of pathogen survival, disease severity, and the response of plants to pathogen attack, including mechanisms of YMD resistance in mungbean. The biggest challenge in YMD management is the effective utilization of an array of information gained so far, in an integrated manner for the development of genotypes having durable resistance against yellow mosaic virus (YMV) infection. In this backdrop, this review summarizes the role of various begomoviruses, its genomic components, and vector whiteflies, including cryptic species in the YMD expression. Also, information about the genetics of YMD in both mungbean and blackgram crops is comprehensively presented, as both the species are crossable, and same viral strains are also found affecting these crops. Also, implications of various management strategies including the use of resistance sources, the primary source of inoculums and vector management, wide-hybridization, mutation breeding, marker-assisted selection (MAS), and pathogen-derived resistance (PDR) are thoroughly discussed. Finally, the prospects of employing various powerful emerging tools like translational genomics, and gene editing using CRISPR/Cas9 are also highlighted to complete the YMD management perspective in mungbean.