Abstract
Phosphate (Pi) deficiency is known to be a major limitation for symbiotic nitrogen fixation (SNF), and hence legume crop productivity globally. However, very little information is available on the adaptive mechanisms, particularly in the important legume crop chickpea (Cicer arietinum L.), which enable nodules to respond to low-Pi availability. Thus, to elucidate these mechanisms in chickpea nodules at molecular level, we used an RNA sequencing approach to investigate transcriptomes of the nodules in Mesorhizobium mediterraneum SWRI9-(MmSWRI9)-chickpea and M. ciceri CP-31-(McCP-31)-chickpea associations under Pi-sufficient and Pi-deficient conditions, of which the McCP-31-chickpea association has a better SNF capacity than the MmSWRI9-chickpea association during Pi starvation. Our investigation revealed that more genes showed altered expression patterns in MmSWRI9-induced nodules than in McCP-31-induced nodules (540 vs. 225) under Pi deficiency, suggesting that the Pi-starvation-more-sensitive MmSWRI9-induced nodules required expression change in a larger number of genes to cope with low-Pi stress than the Pi-starvation-less-sensitive McCP-31-induced nodules. The functional classification of differentially expressed genes (DEGs) was examined to gain an understanding of how chickpea nodules respond to Pi starvation, caused by soil Pi deficiency. As a result, more DEGs involved in nodulation, detoxification, nutrient/ion transport, transcriptional factors, key metabolic pathways, Pi remobilization and signalling were found in Pi-starved MmSWRI9-induced nodules than in Pi-starved McCP-31-induced nodules. Our findings have enabled the identification of molecular processes that play important roles in the acclimation of nodules to Pi deficiency, ultimately leading to the development of Pi-efficient chickpea symbiotic associations suitable for Pi-deficient soils.