Transcript dynamics at early stages of molecular interactions of MYMIV with resistant and susceptible genotypes of the leguminous host, Vigna mungo

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.

Alternative splicing shapes the transcriptome complexity in blackgram [Vigna mungo (L.) Hepper]

Vigna mungo, a highly consumed crop in the pan-Asian countries, is vulnerable to several biotic and abiotic stresses. Understanding the post-transcriptional gene regulatory cascades, especially alternative splicing (AS), may underpin large-scale genetic improvements to develop stress-resilient varieties. Herein, a transcriptome based approach was undertaken to decipher the genome-wide AS landscape and splicing dynamics in order to establish the intricacies of their functional interactions in various tissues and stresses. RNA sequencing followed by high-throughput computational analyses identified 54,526 AS events involving 15,506 AS genes that generated 57,405 transcripts isoforms. Enrichment analysis revealed their involvement in diverse regulatory functions and demonstrated that transcription factors are splicing-intensive, splice variants of which are expressed differentially across tissues and environmental cues. Increased expression of a splicing regulator NHP2L1/SNU13 was found to co-occur with lower intron retention events. The host transcriptome is significantly impacted by differential isoform expression of 1172 and 765 AS genes that resulted in 1227 (46.8% up and 53.2% downregulated) and 831 (47.5% up and 52.5% downregulated) transcript isoforms under viral pathogenesis and Fe²⁺ stressed condition, respectively. However, genes experiencing AS operate differently from the differentially expressed genes, suggesting AS is a unique and independent mode of regulatory mechanism. Therefore, it can be inferred that AS mediates a crucial regulatory role across tissues and stressful situations and the results would provide an invaluable resource for future endeavours in V. mungo genomics.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

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.

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

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

Differential responses of Phaseolus vulgaris cultivars following mungbean yellow mosaic India virus infection

Phaseolus vulgaris, commonly known as French bean is a vital leguminous crop worldwide and India stood 1st rank in dry bean and 4th rank in green bean production worldwide (FAOSTAT 2017). However, this production is severely affected by Mungbean yellow mosaic India virus (MYMIV) infection. Hence it is very important to identify MYMIV tolerant P. vulgaris cultivars. MYMIV infection results in the production of reactive oxygen species and plant cells have evolved complex defense mechanisms at different levels to overcome the damage. Our study for the first time focused on the changes at the morphological and biochemical level, as well as on the relative quantification of MYMIV genes in nine cultivars of P. vulgaris after MYMIV infection. Highest growth and the highest accumulation of four antioxidants of cv. 'Anupam' after MYMIV infection, established that cv. 'Anupam' was less affected by MYMIV infection amongst all nine cultivars. Relative quantification studies also correlated well with these results. Additionally, there is a consistent level of photosynthetic pigments content in mock- and MYMIV-treated seedlings of cv. 'Anupam' over early infection period. Combining all the results we conclude that cv. 'Anupam' is a MYMIV tolerant cultivar.

Plant virus interaction mechanism and associated pathways in mosaic disease of small cardamom (Elettaria cardamomum Maton) by RNA-Seq approach

Small cardamom (Elettaria cardamomum), grown in limited coastal tropical countries is one of the costliest and widely exported agri-produce having global turnover of >10 billion USD. Mosaic/marble disease is one of the major impediments that requires understanding of disease at molecular level. Neither whole genome sequence nor any genomic resources are available, thus RNA seq approach can be a rapid and economical alternative. De novo transcriptome assembly was done with Illumina Hiseq data. A total of 5317 DEGs, 2267 TFs, 114 pathways and 175,952 genic region putative markers were obtained. Gene regulatory network analysis deciphered molecular events involved in marble disease. This is the first transcriptomic report revealing disease mechanism mediated by perturbation in auxin homeostasis and ethylene signalling leading to senescence. The web-genomic resource (SCMVTDb) catalogues putative molecular markers, candidate genes and transcript information. SCMVTDb can be used in germplasm improvement against mosaic disease in endeavour of small cardamom productivity. Availability of genomic resource, SCMVTDb: http://webtom.cabgrid.res.in/scmvtdb/.

Complex molecular mechanisms underlying MYMIV-resistance in Vigna mungo revealed by comparative transcriptome profiling

Mungbean Yellow Mosaic India Virus (MYMIV)-infection creates major hindrance in V. mungo cultivation and poses significant threat to other grain legume production. Symptoms associated include severe patho-physiological alterations characterized by chlorotic foliar lesion accompanied by reduced growth. However, dissection of the host's defense machinery remains a tough challenge due to limited of host's genomic resources. A comparative RNA-Seq transcriptomes of resistant (VM84) and susceptible (T9) plants was carried out to identify genes potentially involved in V. mungo resistance against MYMIV. Distinct gene expression landscapes were observed in VM84 and T9 with 2158 and 1679 differentially expressed genes (DEGs), respectively. Transcriptomic responses in VM84 reflect a prompt and intense immune reaction demonstrating an efficient pathogen surveillance leading to activation of basal and induced immune responses. Functional analysis of the altered DEGs identified multiple regulatory pathways to be activated or repressed over time. Up-regulation of DEGs including NB-LRR, WRKY33, ankyrin, argonaute and NAC transcription factor revealed an insight on their potential roles in MYMIV-resistance; and qPCR validation shows a propensity of their accumulation in VM84. Analyses of the current RNA-Seq dataset contribute immensely to decipher molecular responses that underlie MYMIV-resistance and will aid in the improvement strategy of V. mungo and other legumes through comparative functional genomics.

Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus

BACKGROUND

Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus associated with betasatellite (Tobacco curly shoot betasatellite, TbCSB), which causes serious leaf curl disease on tomato and tobacco in China. It is interesting that TbCSV induced severe upward leaf curling in Nicotiana benthamiana, but in the presence of TbCSB, symptoms changed to be downward leaf curling. However, the mechanism of interactions between viral pathogenicity, host defense, viral-betasatellite interactions and virus-host interactions remains unclear.

METHODS

In this study, RNA-seq was used to analyze differentially expressed genes (DEGs) in N. benthamiana plants infected by TbCSV (Y35A) and TbCSV together with TbCSB (Y35AB) respectively.

RESULTS

Through mapping to N. benthamiana reference genome, 59,814 unigenes were identified. Transcriptome analysis revealed that a total of 4081 and 3196 DEGs were identified in Y35AB vs CK (control check) and Y35A vs CK, respectively. Both GO and KEGG analyses were conducted to classify the DEGs. Ten of the top 15 GO terms were enriched in both DEGs of Y35AB vs CK and Y35A vs CK, and these enriched GO terms mainly classified into three categories including biological process, cellular component and molecular function. KEGG pathway analysis indicated that 118 and 111 pathways were identified in Y35AB vs CK and Y35A vs CK, respectively, of which nine and six pathways were significantly enriched. Three major pathways in Y35AB vs CK involved in metabolic pathways, carbon metabolism and photosynthesis, while those in Y35A vs CK were related to Ribosome, Glyoxylate and dicarboxylate metabolism and DNA replication. We observed that 8 PR genes were significantly up-regulated and 44 LRR-RLK genes were significantly differentially expressed in Y35A treatment or in Y35AB treatment. In addition, 7 and 13 genes were identified to be significantly changed in biosynthesis and signal transduction pathway of brassinosteroid (BR) and jasmonic acid (JA) respectively.

CONCLUSIONS

These results presented here would be particularly useful to further elucidate the response of the host plant against virus infection.

Screening and Identification of putative long non coding RNAs from transcriptome data of a high yielding blackgram (Vigna mungo), Cv. T9

Blackgram (Vigna mungo) is one of primary legumes cultivated throughout India, Cv.T9 being one of its common high yielding cultivar. This article reports RNA sequencing data and a pipeline for prediction of novel long non-coding RNAs from the sequenced data. The raw data generated during sequencing are available at Sequence Read Archive (SRA) of NCBI with accession number- SRX1558530

Exogenous application of methyl jasmonate induces defense response and develops tolerance against mungbean yellow mosaic India virus in Vigna mungo

Vigna mungo (L.)Hepper is an economically important leguminous crop in south-east Asia. However, its production is severely affected by Mungbean yellow mosaic India virus (MYMIV). It is well established that methyl jasmonate (MeJA) is effective in inducing resistance against pathogens in several plants. To assess the role of MeJA in developing MYMIV tolerance in V. mungo, we analysed time-dependent biochemical and molecular responses of MYMIV susceptible V. mungo after exogenous application of different MeJA concentrations, followed by MYMIV infection. Our analysis revealed that exogenous application of different concentrations of MeJA resulted in decreased levels of malondialdehyde with higher membrane stability index values in MYMIV susceptible V. mungo, suggesting the protective role of MeJA through restoring the membrane stability. Moreover, the level of expression of different antioxidative enzymes revealed that exogenous MeJA is also very effective in ROS homeostasis maintenance. Enhanced expressions of the defence marker genes lipoxygenase and phenylalanine ammonia-lyase and the reduced expression of the MYMIV coat-protein encoding gene in all MeJA treated plants post MYMIV infection revealed that exogenous application of MeJA is effective for MYMIV tolerance in V. mungo. Our findings provide new insights into the physiological and molecular mechanisms of MYMIV tolerance in Vigna induced by MeJA.

High throughput transcriptome analysis of coffee reveals prehaustorial resistance in response to Hemileia vastatrix infection

We provide a transcriptional profile of coffee rust interaction and identified putative up regulated resistant genes Coffee rust disease, caused by the fungus Hemileia vastatrix, is one of the major diseases in coffee throughout the world. The use of resistant cultivars is considered to be the most effective control strategy for this disease. To identify candidate genes related to different mechanism defense in coffee, we present a time-course comparative gene expression profile of Caturra (susceptible) and Híbrido de Timor (HdT, resistant) in response to H. vastatrix race XXXIII infection. The main objectives were to obtain a global overview of transcriptome in both interaction, compatible and incompatible, and, specially, analyze up-regulated HdT specific genes with inducible resistant and defense signaling pathways. Using both Coffea canephora as a reference genome and de novo assembly, we obtained 43,159 transcripts. At early infection events (12 and 24 h after infection), HdT responded to the attack of H. vastatrix with a larger number of up-regulated genes than Caturra, which was related to prehaustorial resistance. The genes found in HdT at early hours were involved in receptor-like kinases, response ion fluxes, production of reactive oxygen species, protein phosphorylation, ethylene biosynthesis and callose deposition. We selected 13 up-regulated HdT-exclusive genes to validate by real-time qPCR, which most of them confirmed their higher expression in HdT than in Caturra at early stage of infection. These genes have the potential to assist the development of new coffee rust control strategies. Collectively, our results provide understanding of expression profiles in coffee-H. vastatrix interaction over a time course in susceptible and resistant coffee plants.

Molecular modeling and simulation of three important components of Plant Pathogen Interaction cascade in Vigna mungo

Plant pathogen interaction plays a great role in plant immunity. The regulation of various components of plant pathogen interactions is quite complicated and is very important in establishing relationship among components of this system. Yellow Mosaic Disease is common among legumes such as Vigna mungo. Mungbean Yellow Mosaic India Virus (MYMIV) and whitefly (Bemisia tabaci) is a vector causing the disease. Therefore, it is of interest to document the molecule models of three different components of Plant Pathogen interaction cascade- MAP kinase1, MAP kinase 2 and WRKY33 from V. mungo resistant to MYMIV. Both the MAP kinases were sequenced for this study while WRKY 33 was extracted and modeled from transcripts generated from two different transcriptome libraries, one set MYMIV- challenged, the other fed with aviruliferous whitefly. Post simulation studies revealed that MAPKs contained less percentage of disordered residues and were structurally more stable and than WRKY33.

Development of transcriptome based web genomic resources of yellow mosaic disease in Vigna mungo

Vigna mungo (Urdbean) is cultivated in the tropical and sub-tropical continental region of Asia. It is not only important source of dietary protein and nutritional elements, but also of immense value to human health due to medicinal properties. Yellow mosaic disease caused by Mungbean Yellow Mosaic India Virus is known to incur huge loss to crop, adversely affecting crop yield. Contrasting genotypes are ideal source for knowledge discovery of plant defence mechanism and associated candidate genes for varietal improvement. Whole genome sequence of this crop is yet to be completed. Moreover, genomic resources are also not freely accessible, thus available transcriptome data can be of immense use. V. mungo Transcriptome database, accessible at http://webtom.cabgrid.res.in/vmtdb/ has been developed using available data of two contrasting varieties viz., cv. VM84 (resistant) and cv. T9 (susceptible). De novo assembly was carried out using Trinity and CAP3. Out of total 240,945 unigenes, 165,894 (68.8%) showed similarity with known genes against NR database, and remaining 31.2% were found to be novel. We found 22,101 differentially expressed genes in all datasets, 44,335 putative genic SSR markers, 4105 SNPs and Indels, 64,964 transcriptional factor, 546 mature miRNA target prediction in 703 differentially expressed unigenes and 137 pathways. MAPK, salicylic acid-binding protein 2-like, pathogenesis-related protein and NBS-LRR domain were found which may play an important role in defence against pathogens. This is the first web genomic resource of V. mungo for future genome annotation as well as ready to use markers for future variety improvement program.

Molecular and biochemical characterization of a Vigna mungo MAP kinase associated with Mungbean Yellow Mosaic India Virus infection and deciphering its role in restricting the virus multiplication

Yellow Mosaic Disease caused by the begomovirus Mungbean Yellow Mosaic India Virus (MYMIV) severely affects many economically important legumes. Recent investigations in Vigna mungo - MYMIV incompatible interaction identified a MAPK homolog in the defense signaling pathway. An important branch of immunity involves phosphorylation by evolutionary conserved Mitogen-activated protein kinases (MAPK) that transduce signals of pathogen invasion to downstream molecules leading to diverse immune responses. However, most of the knowledge of MAPKs is derived from model crops, and functions of these versatile kinases are little explored in legumes. Here we report characterization of a MAP kinase (VmMAPK1), which was induced upon MYMIV-inoculation in resistant V. mungo. Phylogenetic analysis revealed that VmMAPK1 is closely related to other plant-stress-responsive MAPKs. Both mRNA and protein of VmMAPK1 were accumulated upon MYMIV infection. The VmMAPK1 protein localized in the nucleus as well as cytoplasm and possessed phosphorylation activity in vitro. A detailed biochemical characterization of purified recombinant VmMAPK1 demonstrated an intramolecular mechanism of autophosphorylation and self-catalyzed phosphate incorporation on both threonine and tyrosine residues. The V_max and K_m values of recombinant VmMAPK1 for ATP were 6.292nmol/mg/min and 0.7978μM, respectively. Furthermore, the ability of VmMAPK1 to restrict MYMIV multiplication was validated by its ectopic expression in transgenic tobacco. Importantly, overexpression of VmMAPK1 resulted in the considerable upregulation of defense-responsive marker PR genes. Thus, the present data suggests the critical role of VmMAPK1 in suppressing MYMIV multiplication presumably through SA-mediated signaling pathway and inducing PR genes establishing the significant implications in understanding MAP kinase gene function during Vigna-MYMIV interaction; and hence paves the way for introgression of resistance in leguminous crops susceptible to MYMIV.

Bruchid egg induced transcript dynamics in developing seeds of black gram (Vigna mungo)

Black gram (Vigna mungo) seeds are a rich source of digestible proteins, however, during storage these seeds are severely damaged by bruchids (Callosobruchus spp.), reducing seed quality and yield losses. Most of the cultivated genotypes of black gram are susceptible to bruchids, however, few tolerant genotypes have also been identified but the mechanism of tolerance is poorly understood. We employed Suppression Subtractive Hybridization (SSH) to identify specifically, but rarely expressed bruchid egg induced genes in black gram. In this study, Suppression Subtractive Hybridization (SSH) library was constructed to study the genes involved in defense response in black gram against bruchid infestation. An EST library of 277 clones was obtained for further analyses. Based on CAP3 assembly, 134 unigenes were computationally annotated using Blast2GOPRO software. In all, 20 defense related genes were subject to quantitative PCR analysis (qPCR) out of which 12 genes showed up-regulation in developing seeds of the pods oviposited by bruchids. Few major defense genes like defensin, pathogenesis related protein (PR), lipoxygenase (LOX) showed high expression levels in the oviposited population when compared with the non-oviposited plants. This is the first report on defense related gene transcript dynamics during the bruchid-black gram interaction using SSH library. This library would be useful to clone defense related gene(s) such as defensin as represented in our library for crop improvement.

High throughput sequencing reveals modulation of microRNAs in Vigna mungo upon Mungbean Yellow Mosaic India Virus inoculation highlighting stress regulation

MicroRNAs (miRNAs) are 20-24 nucleotides long non-coding RNAs known to play important regulatory roles during biotic and abiotic stresses by controlling gene expression. Blackgram (Vigna mungo), an economically important grain legume is highly susceptible to pathogenic begomovirus Mungbean Yellow Mosaic India Virus (MYMIV) and resulting in high yield loss. In this study two different leaf-small-RNA libraries were prepared from the pooled RNA at three different time points of resistant V. mungo inbred line VM84 inoculated either with viruliferous or non-viruliferous whiteflies carrying MYMIV and performed high-throughput Illumina sequencing. Sequencing followed by bioinformatics analysis of the small RNA reads indicated that the expression patterns of most of the known and novel miRNAs were altered in resistant line over mock-inoculated sample during the plant virus incompatible interaction. Highly altered miRNAs belong to the families of miR156, miR159, miR160, miR166, miR398, miR1511, miR1514, miR2118 and novel vmu-miRn7, vmu-miRn8, vmu-miRn13 and vmu-miRn14. These results were validated using qPCR, and most of the miRNAs showed similar pattern of expression like that of Illumina reads. The expression patterns of some selected known and novel miRNAs were also compared between the infected MYMIV-resistant and -susceptible genotypes and most of these were modulated after MYMIV-inoculation. Target transcripts like NB-LRR, NAC, MYB, Zinc finger, CCAAT-box transcription factor, fructose 2-6 bisphosphate, HDZIP protein that confers immune response were predicted as targets amongst identified miRNAs using psRNATarget server. Some selected target transcripts including NB-LRR, ARF, SOD, SPB, Basic blue copper protein were validated and their differential expression were demonstrated between MYMIV-resistant and -susceptible V. mungo by qPCR data analyses. In the present study we have identified miRNAs that implicate in the regulation of MYMIV-induced stress response in V. mungo; and generated genomic resources for a non-model legume with the aid of bioinformatics tools supplemented by experimental validation.

Analyses of MYMIV-induced transcriptome in Vigna mungo as revealed by next generation sequencing

Mungbean Yellow Mosaic Virus (MYMIV) is the viral pathogen that causes yellow mosaic disease to a number of legumes including Vigna mungo. VM84 is a recombinant inbred line resistant to MYMIV, developed in our laboratory through introgression of resistance trait from V. mungo line VM-1. Here we present the quality control passed transcriptome data of mock inoculated (control) and MYMIV-infected VM84, those have already been submitted in Sequence Read Archive (SRX1032950, SRX1082731) of NCBI. QC reports of FASTQ files generated by ‘SeqQC V2.2’ bioinformatics tool.