RNA-Seq analysis of mung bean (Vigna radiata L.) roots shows differential gene expression and predicts regulatory pathways responding to taxonomically different rhizobia

Symbiotic interaction among legume and rhizobia is a complex phenomenon which results in the formation of nitrogen-fixing nodules. Mung bean is promiscuous host however expression profile of this important legume plant in response to rhizobial infection was particularly lacking and urgently needed. We have demonstrated the pattern of gene expression of mung bean roots inoculated with two symbionts Bradyrhizobium yuanmingense Vr50 and Sinorhizobium (Ensifer) aridi Vr33 and non-inoculated control (CK). The RNA-Seq data analyzed at two growth stages i.e., 1–3 h and 10–16 days post inoculation revealed significantly higher number of differentially expressed genes (DEGs) at nodulation stage. The DEGs encoding receptor kinases identified at early stage might be involved in perception of Nod factors produced by different rhizobia. At nodulation stage important genes involved in plant hormone signal transduction, nitrogen and sulfur metabolism were identified. KEGG pathway enrichment analysis showed that metabolic pathways were most prominent in both groups (Group 1: Vr33 vs CK; Group 2: Vr50 vs CK), followed by biosynthesis of secondary metabolites, plant hormone signal transduction and biosynthesis of amino acids. Furthermore, DEGs involved in cell communication and plant hormone signal transduction were found to be different among two symbiotic systems while DEGs involved in carbon, nitrogen and sulfur metabolism were similar but their expression varied in response to two rhizobial strains. This study provides the first insight into the mechanisms underlying interactions of mung bean host with two taxonomically different symbionts (Bradyrhizobium and Sinorhizobium) and the candidate genes for better understanding the mechanisms of symbiotic host-specificity.