Early developmental and stress responsive ESTs from mungbean, Vigna radiata (L.) Wilczek, seedlings

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.

Transcriptomic analysis of Vigna radiata in response to chilling stress and uniconazole application

Background

Chilling injury of mung bean (Vigna radiata (L.)) during the blooming and podding stages is a major agricultural threat in Northeast China. Uniconazole (UNZ) can alleviate water deficit stress in soybean and waterlogging stress in mung bean. However, there has been no report on the effect of UNZ application on the growth and transcriptomic profile of mung bean under chilling stress.

Results

UNZ application before chilling stress at the R1 stage alleviated the decline in mung bean yield. UNZ delayed the decrease in leaf chlorophyll content under chilling stress at the R1 stage and accelerated the increase in leaf chlorophyll content during the recovery period. Eighteen separate RNA-Seq libraries were generated from RNA samples collected from leaves exposed to six different treatment schemes. The numbers of DEGs specific for UNZ treatment between D1 + S vs. D1 and D4 + S vs. D4 were 708 and 810, respectively. GO annotations showed that photosynthesis genes were obviously enriched among the genes affected by chilling stress and UNZ application. KEGG pathway enrichment analysis indicated that 4 pathways (cutin, suberin and wax biosynthesis; photosynthesis; porphyrin and chlorophyll metabolism; and ribosome) were downregulated, while plant–pathogen interaction was upregulated, by chilling stress. UNZ application effectively prevented the further downregulation of the gene expression of members of these 4 KEGG pathways under chilling stress.

Conclusions

UNZ application effectively delayed the decrease in photosynthetic pigment content under chilling stress and accelerated the increase in photosynthetic pigment content during the recovery period, thus effectively limiting the decline in mung bean yield. UNZ application effectively prevented the further downregulation of the gene expression of members of 4 KEGG pathways under chilling stress and increased mung bean tolerance to chilling stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08443-6.

A near-complete genome sequence of mungbean (Vigna radiata L.) provides key insights into the modern breeding program

Mungbean (Vigna radiata L.), a fast-growing legume species, is an important source of carbohydrates and proteins in developing countries of Asia. Here, we constructed a near-complete genome sequence of mungbean with a scaffold N50 value of 5.2 Mb and only a 0.4% gap, with a total scaffold size of 475 Mb. We identified several misassembled pseudomolecules (Chr03, Chr04, Chr05, and Chr08) in the previous draft assembly; Chr03, Chr04, and Chr08 were assembled into one chromosome, and Chr05 was broken into two chromosomes in the improved reference genome assembly, thus providing more accurate linkage information to breeders. Additionally, using an ultra-high-resolution linkage map constructed based on resequencing data, we identified several quantitative trait loci (QTLs) and the underlying candidate genes affecting synchronous pod maturity (SPM). Mungbean homologs of two soybean ([Glycine max (L.) Merr.] flowering genes, E3 (phytochrome A) and J (early flowering 3), were identified as candidate genes for the QTLs, and the candidate genes for plant height, node number, and SPM showed critical nucleotide substitutions between the reference cultivar and other genotypes (landraces and wild accessions). Based on the analysis of genetic diversity among 276 accessions collected from 23 countries, we identified 36 selective sweep regions and observed that the overall genetic diversity of cultivars decreased to 30% of that in wild accessions postdomestication. The near-complete genome sequence of mungbean represents an important resource for genome-assisted improvement in the mungbean breeding program.

Proteomic analysis of the response of Funnelifor mismosseae/Medicago sativa to atrazine stress

BACKGROUND

Arbuscular mycorrhizal (AM) fungi form symbiotic associations with host plants can protect host plants against diverse biotic and abiotic stresses, and promote biodegradation of various contaminants. However, the molecular mechanisms of how the arbuscular mycorrhizal fungi and host plant association on atrazine stress were still poorly understood. To better characterize how arbuscular mycorrhizal fungi and host plant interactions increase atrazine stress, we performed physiological and proteomic analysis of Funneliformis mosseae (mycorrhizal fungi) and Medicago sativa (alfalfa) association under atrazine stress.

RESULTS

The results showed that in the Arbuscular mycorrhizal, protective enzymes were up regulated and the malondialdehyde content increased relative to those of non-mycorrhizal M.sativa. We also examined the atrazine degradation rates within the nutrient solution, and a 44.43% reduction was observed with the mycorrhizal M.sativa, with 30.83% of the reduction attributed to F. mosseae. The accumulation content in root and stem of mycorrhizal M.sativa were obviously increased 11.89% and 16.33% than those of non- mycorrhizal M.sativa. The activity of PPO, POD, CAT and SOD in mycorrhizal M.sativa were obviously higher than non mycorrhizal M.sativa under atrazine stess. We identified differential root proteins using isobaric tags for relative and absolute quantization coupled with liquid chromatography-mass spectrometry, with 533 proteins identified (276 unregulated and 257 downregulated). The differentially expressed proteins were further examined using GO, BLAST comparisons, and a literature inquiry and were classified into the categories of atrazine degradation (37.1%); atrazine stress response (28.6%); plant immune responses (14.3%); translation, synthesis, and processing (10%); and signal transduction and biological processes (10%). Furthermore, we identified glycosyl transferase, glutathione S-transferase, laccase, cytochrome P450 monooxygenase, peroxidase, and other proteins closely related to the degradation process.

CONCLUSIONS

Mycorrhizal Medicago showed improved atrazine degradation within the culturing medium and increased atrazine enrichment in the roots and stems. Additionally, AMF increased the plant root response to atrazine, with relevant enzymes up regulated and toxic effects alleviated. Overall, the findings of this study show that AMF played an important role in easing atrazine stress in plants and contributed to atrazine remediation and further contributed to the understanding of the molecular mechanism associated with atrazine stresses and potential mycorrhizal contributions in M.sativa.

Chilling susceptibility in mungbean varieties is associated with their differentially expressed genes

BACKGROUND

Mungbean (Vigna radiata L. Wilczek) is an economically important legume of high nutritional value, however, its cultivation is limited by susceptibility to chilling. Varieties NM94 and VC1973A, with differential susceptibility to stress, serve as good materials for uncovering how they differ in chilling tolerance. This study aimed to identify the ultrastructural, physiological and molecular changes to provide new insights on the differential susceptibility to chilling between varieties VC1973A and NM94.

RESULTS

Chilling stress caused a greater reduction in relative growth rate, a more significant decrease in maximum photochemical efficiency of PSII and DPPH scavenging activity and more-pronounced ultrastructural changes in VC1973A than in NM94 seedlings. Comparative analyses of transcriptional profiles in NM94 and VC1973A revealed that the higher expression of chilling regulated genes (CORs) in NM94. The transcript levels of lipid transfer protein (LTP), dehydrin (DHN) and plant defensin (PDF) in NM94 seedlings after 72 h at 4 °C was higher than that in its parental lines VC1973A, 6601 and VC2768A.

CONCLUSIONS

Our results suggested that LTP, DHN and PDF may mediate chilling tolerance in NM94 seedlings.

Gene expression associated with intersterility in Heterobasidion

Intersterility (IS) is thought to prevent mating compatibility between homokaryons that belong to different species. Although IS in Heterobasidion is regulated by the genes located at the IS loci, it is not yet known how the IS genes influence sexual compatibility and heterokaryon formation. To increase our understanding of the molecular events underlying IS, we studied mRNA abundance changes during IS compatible and incompatible interactions over time. The clustering of the transcripts into expression profiles, followed by the application of Gene Ontology (GO) enrichment pathway analysis of each of the clusters, allowed inference of biological processes participating in IS. These analyses identified events involved in mating and sexual development (i.e., linked with IS compatibility), which included processes associated with cell-cell adhesion and recognition, cell cycle control and signal transduction. We also identified events potentially involved in overriding mating between individuals belonging to different species (i.e., linked with IS incompatibility), which included reactive oxygen species (ROS) production, responses to stress (especially to oxidative stress), signal transduction and metabolic biosynthesis. Our findings thus enabled detection and characterization of gene expression changes associated with IS in Heterobasidion, as well as identification of important processes and pathways associated with this phenomenon. Overall, the results of this study increase current knowledge regarding the molecular mechanisms underpinning IS in Heterobasidion and allowed for the establishment of a vital baseline for further studies.

Cold response in Phalaenopsis aphrodite and characterization of PaCBF1 and PaICE1

Phalaenopsis is a winter-blooming orchid genus commonly cultivated in tropical Asian countries. Because orchids are one of the most economically important flower crops in Taiwan, it is crucial to understand their response to cold and other abiotic stresses. The present study focused on gene regulation of P. aphrodite in response to abiotic stress, mainly cold. Our results demonstrate that P. aphrodite is sensitive to low temperatures, especially in its reproductive stage. We found that after exposure to 4°C, plants in the vegetative stage maintained better membrane integrity and photosynthetic capacity than in the flowering stage. At the molecular level, C-repeat binding factor1 (PaCBF1) and its putative target gene dehydrin1 (PaDHN1) mRNAs were induced by cold, whereas inducer of CBF expression1 (PaICE1) mRNA was constitutively expressed. PaICE1 transactivated MYC motifs in the PaCBF1 promoter, indicating that up-regulation of PaCBF1 may be mediated by the binding of PaICE1 to MYC motifs. Overexpression of PaCBF1 in transgenic Arabidopsis induced AtCOR6.6 and RD29a without cold stimulus and maintained better membrane integrity after cold stress. Herein, we present evidence that cold induction of PaCBF1 transcripts in P. aphrodite may be transactivated by PaICE1 and consequently protect plants from cold damage through up-regulation of cold-regulated (COR) genes, such as DHN. To our knowledge, this study is the first report of the isolation and characterization of CBF, DHN and ICE genes in the Orchidaceae family.

Proteomic analysis of soybean defense response induced by cotton worm (prodenia litura, fabricius) feeding

BACKGROUND

Cotton worm is one of the main insects of soybean in southern China. Plants may acquire defense mechanisms that confer protection from predation by herbivores. Induced responses can lead to increased resistance against herbivores in many species. This study focuses on searching changed proteins in soybean defense response induced by cotton worm feeding.

RESULTS

Ten protein spots that are changed in abundance in response to cotton worm feeding were identified by Two-dimensional gel electrophoresis (2-DE). A total of 11 unique proteins from these spots were identified by MALDI-TOF MS. The mRNA and protein relative expression levels of most changed proteins were up-regulated. These proteins were mainly involved in physiological processes, including active oxygen removal, defense signal transduction, and metabolism regulation.

CONCLUSION

This is the first proteomic analysis of the soybean defense response induced by cotton worm. The differentially expressed proteins could work together to play a major role in the induced defense response. PAL and SAMS were up-regulated at both the protein and mRNA levels. These genes can be strongest candidates for further functional research.

Changes in the nature of phenolic deposits after re-warming as a result of melatonin pre-sowing treatment of Vigna radiata seeds

Changes in phenolics (PhC - phenolic compounds) measured as UV-absorbing compounds (UVAC) and their localization as well as growth, lipid peroxidation (TBARS level) and proline (Pro) level in three-day-old roots of seedlings (T(o) stage) obtained from hydroprimed (H) and hydroprimed with melatonin (H-MEL) seeds after 2 days of chilling (5°C) and 2 days of re-warming were examined. H and H-MEL resulted in inhibition of root growth under optimal conditions, but after re-warming, a positive effect of MEL was noted. The results also showed a positive MEL impact on TBARS level already after chilling and especially after re-warming. Exposure of Vigna radiata seedlings to chilling caused a significant increase in Pro level, especially in H-MEL roots, but after re-warming it drastically decreased. Under chilling stress, accumulation of UVAC also decreased. However, after re-warming it returned to the level observed in the roots grown constantly at 25°C. Even if after re-warming of V. radiata seedlings only slight accumulation of total PhC was observed, phenolic deposits accumulating in the vacuoles of H-MEL roots were completely different from those in the vacuoles of the control and H roots. H-MEL application to the seeds resulted in a significant increase in small granular composite materials, while in the control and H roots, large oval deposits prevailed. Taken together, it is probable that all of these differences were connected with positive effects of MEL on chilled V. radiata seedlings after re-warming.

A raison d'être for two distinct pathways in the early steps of plant isoprenoid biosynthesis?

When compared to other organisms, plants are atypical with respect to isoprenoid biosynthesis: they utilize two distinct and separately compartmentalized pathways to build up isoprene units. The co-existence of these pathways in the cytosol and in plastids might permit the synthesis of many vital compounds, being essential for a sessile organism. While substrate exchange across membranes has been shown for a variety of plant species, lack of complementation of strong phenotypes, resulting from inactivation of either the cytosolic pathway (growth and development defects) or the plastidial pathway (pigment bleaching), seems to be surprising at first sight. Hundreds of isoprenoids have been analyzed to determine their biosynthetic origins. It can be concluded that in angiosperms, under standard growth conditions, C₂₀-phytyl moieties, C₃₀-triterpenes and C₄₀-carotenoids are made nearly exclusively within compartmentalized pathways, while mixed origins are widespread for other types of isoprenoid-derived molecules. It seems likely that this coexistence is essential for the interaction of plants with their environment. A major purpose of this review is to summarize such observations, especially within an ecological and functional context and with some emphasis on regulation. This latter aspect still requires more work and present conclusions are preliminary, although some general features seem to exist.

Capturing cold-stress-related sequence diversity from a wild relative of common bean (Phaseolus angustissimus)

One restriction to the cultivation of common bean, Phaseolus vulgaris L., is its limited tolerance to low temperatures. In the present study, subtraction suppression hybridization was employed to enrich for stress responsive genes in both a chilling-susceptible common bean and a relatively more chilling-tolerant wild bean species, Phaseolus angustissimus. For each species, approximately 11 000 expressed sequence tags were generated. Comparative sequence analysis of the EST collection with the available annotated genome sequences of the model Fabaceae species Medicago truncatula and Glycine max identified protein homologues for approximately 65% and 80% of the Phaseolus sequences, respectively. This difference reflects the closer phylogenetic relationship between the genera Phaseolus and Glycine compared with Medicago. Annotation of the Phaseolus sequences was facilitated through this comparative analysis and indicated that several heat shock proteins, cytochrome P450s, and DNA binding factors were uniquely found among the sequences from the wild species P. angustissimus. The Phaseolus sequences have been made available on a GBrowse implementation using M. truncatula as the reference genome, providing rapid access to the sequence data and associated comparative genome data.

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.