Skip to main content

Chickpea shows genotype-specific nodulation responses across soil nitrogen environment and root disease resistance categories



The ability of chickpea to obtain sufficient nitrogen via its symbiotic relationship with Mesorhizobium ciceri is of critical importance in supporting growth and grain production. A number of factors can affect this symbiotic relationship including abiotic conditions, plant genotype, and disruptions to host signalling/perception networks. In order to support improved nodule formation in chickpea, we investigated how plant genotype and soil nutrient availability affect chickpea nodule formation and nitrogen fixation. Further, using transcriptomic profiling, we sought to identify gene expression patterns that characterize highly nodulated genotypes.


A study involving six chickpea varieties demonstrated large genotype by soil nitrogen interaction effects on nodulation and further identified agronomic traits of genotypes (such as shoot weight) associated with high nodulation. We broadened our scope to consider 29 varieties and breeding lines to examine the relationship between soilborne disease resistance and the number of nodules developed and real-time nitrogen fixation. Results of this larger study supported the earlier genotype specific findings, however, disease resistance did not explain differences in nodulation across genotypes. Transcriptional profiling of six chickpea genotypes indicates that genes associated with signalling, N transport and cellular localization, as opposed to genes associated with the classical nodulation pathway, are more likely to predict whether a given genotype will exhibit high levels of nodule formation.


This research identified a number of key abiotic and genetic factors affecting chickpea nodule development and nitrogen fixation. These findings indicate that an improved understanding of genotype-specific factors affecting chickpea nodule induction and function are key research areas necessary to improving the benefits of rhizobial symbiosis in chickpea.

Peer Review reports


The majority of legume and pulse crops are able to supplement their nitrogen (N) nutrition by forming a symbiotic association with N-fixing rhizobial bacteria. Farmers utilize this symbiotic relationship with rhizobia to counter nutrient limitations in soils as these bacteria can provide up to 97% of a plants total N requirements [1]. Further, because these symbioses are responsible for returning significant N into agricultural soils globally [2], they reduce the need for costly fertiliser applications. Unlike other diazotrophic bacteria that fix N while living freely in the soil, rhizobia must first colonise plant roots, forming specialised root organs called nodules. The formation of these nodules is controlled both by a series of abiotic cues and a number of important signalling events between the bacteria and the plant host. Chickpea (Cicer arietinum) is a legume of major economic value for food production, with chickpeas consistently being the second or third largest pulse crop grown around the world with nearly 11 million hectares in production [3]. Despite the formation of a functional symbiotic relationship between chickpea and the rhizobial species Mesorhizobium ciceri, N-fixation in chickpeas is limited as compared to other pulses. Field studies show that the chickpea-M. ciceri symbiosis often only fixes sufficient N to support plant growth with little or no enhancement of soil N concentrations after harvest [4, 5]. In addition, competition between inoculated and endemic but less efficient Mesorhizobium strains, as well as soil environment effects both affect successful symbiosis with chickpea [6]. Therefore, a better understanding of the factors that could improve the benefit of this symbiotic association would be advantageous in agriculture.

Abiotic factors such as soil moisture, pH and plant nutrient availability in soils have been identified as affecting nodule formation and N-fixation. With regards to soil nutrients, chickpea nodulation and N-fixation responses are highly dependent on plant available N concentrations in soils. For example, the chickpea var. Tyson, which receives between 60–80% of its N requirements from symbiotic fixation at soil nitrate concentrations between 20–90 kg nitrate/ha, shows a logarithmic decline to less than 25% dependency on its bacterial partner at > 125 kg nitrate/ha [7]. This research suggests that if chickpea crops are planted in moderate to high nitrogen containing agricultural soils, this negative feedback of nitrogen on N-fixation will lead to greater crop reliance on soil N. This will limit the ability of chickpea to enhance soil N for a following crop and reduces the sustainability of this particular legume in agriculture rotations. Therefore, it will become increasingly desirable to identify or breed genotypes that are less sensitive to high levels of plant-available N in soils to maintain the ability to form and acquire N from rhizobial symbiosis.

Understanding how genotype supports nodulation and N-fixation in chickpea is also important to maximise the benefit of N-fixation and reduce reliance on nitrogenous fertilisers. Surprisingly there have been few large-scale studies comparing chickpea nodulation and nitrogen production across multiple genotypes, with a recent review on this topic only citing two studies [8]. One study compared three American and three Middle Eastern chickpea varieties while the second compared 39 chickpea lines from a USDA collection [9, 10]. This latter study demonstrated that genotype vigour, as indicated by total plant biomass, was one of the most important factors affecting N-fixation, accounting for one third of the variation in N-fixation across the accessions [10]. This vigour was considered to be related to a larger plant biomass providing a stronger N sink [10]. Genotypes bred to exhibit increased disease resistance may also affect rhizobial symbiosis in chickpea. The productivity and profitability of chickpea crops are undermined each year by a range of root diseases, including Phytophthora root rot (PRR; caused by Phytophthora medicaginis) and Fusarium wilt (FW, caused by Fusarium oxysporum f.sp. ciceris). While the breeding of elite disease-resistant varieties is one of the most effective means used to combat diseases, especially root diseases [11, 12], genotype-based differences in nodulation and the proportion of plant-N derived from fixation were shown between PRR susceptible and moderately resistant chickpea varieties [13]. These differences were attributed to increases in the general basal defense processes of the moderately resistant varieties that, inadvertently, also compromised mutualist symbiosis [13]. Therefore, plant genotypes that support rhizobial symbiosis in chickpea need to be further investigated.

Until recently, sufficient genetic tools have not been available in chickpea to refine our understanding of the molecular bases of nodulation, but the release of 1 full and 29 re-sequenced genomes of chickpea [14] can facilitate the identification of genes related to the establishment of nodulation, and help understand how these interact with disease resistance response pathways, nitrogen perception and transport. The initial steps of nodule formation have been well studied in model legumes such as Medicago truncatula and Lotus japonicus [15]. From these systems, the general model is that upon perception of rhizobia, plant secondary metabolites, dominated by flavonoid compounds, are secreted by the host root and taken up by rhizobia leading to the production of secreted bacterial Nod-factors [16, 17]. These bacterial signals are perceived by host plant LysM receptor-like kinases, which induce a signalling cascade in the plant root called the Common Symbiosis Signalling Pathway (CSSP) and the induction of plant hormone pathways including cytokinin and auxin [18, 19]. Interruptions in these signalling steps can reduce the efficiency of nodulation or halt it altogether. Legumes also have a feedback signalling loop to control nodulation as, while an increased number of nodules may improve host plant N, excessive numbers of nodules negatively impact plant health due to the energy costs to the plant of both nodule formation and maintenance. Called the Autoregulation Of Nodulation (AON) feedback mechanism (reviewed by Ferguson et al. [20]), this pathway is regulated both by internal factors to the plant as well as by abiotic factors, particularly soil plant-available N [21,22,23]. It is likely that as we expand the number of model legume systems studied, that further genetic factors supporting nodulation will be found and, as such, the factors that affect these signalling events are an area of active research globally.

In this current study, we sought to further understand how Mesorhizobium ciceri (CC1192 [24]) inoculated chickpea genotypes differed in nodulation and nodule function. First, we consider how external factors (i.e. plant-available N) interact with genotype, as well as how intrinsic factors such as quantitative and qualitative disease resistance may alter nodulation potential. Using transcriptomic profiling of a core set of high- and low-nodulating chickpea genotypes, we then sought to understand what differential gene transcription supports high nodulation. Altogether, the results of this study inform our greater understanding of how host genotype and gene regulation affect symbiotic nodule formation and function and is an important step forward in identifying what genetic factors enable chickpea to benefit from rhizobia.


Nodule formation in chickpea varies with soil conditions in a genotype-specific manner

Six genotypes of chickpea (Sonali, PBA HatTrick, Jimbour, Moti, Kyabra and Yorker) were grown in soils with varying plant-available nitrogen (N): paddock soil with nitrate-based fertilizer (highest plant-available N; HN), paddock soil without N fertilizer (moderate plant-available N; MN) and paddock soil cut with sand (2 parts sand to 1 part soil; lowest plant-available N; LN). We found that chickpea genotype significantly interacted with soil N treatment for nodule formation (P < 0.001; Table 1a). The differential impact of genotype on how soil N treatment affected nodulation was seen in varieties such as PBA HatTrick, which had no significant change in nodule numbers across the N treatments, whereas others like Moti and Yorker had significantly more nodules in MN compared to HN soil. Genotype root dry weight values also significantly interacted with soil N treatment (P < 0.001). PBA HatTrick had a higher root dry weight than Yorker in both LN and HN soil, but in the MN soil Yorker had a higher root dry weight than PBA HatTrick (Table 1b). Similarly, a genotype by soil N interaction was observed for nodule numbers normalised for differences in root dry weight (P < 0.05; Table 1c) whereby all genotypes had significant increases in the number of nodules/mg root from the HN to the MN soil. Conversely, nodule numbers/mg root significantly decreased from the MN to the LN soil treatments for all genotypes with the exception of Jimbour and Sonali, indicating genotype specific thresholds of soil plant available N required to support optimal nodulation in chickpea.

Table 1 a) Number of nodules per plant b) root dry weight (g) and c) number of nodules per g of root tissue from the six genotype experiment with three soil nitrogen treatments: high nitrogen (HN, paddock soil + N), medium nitrogen (MN, paddock soil no N) and low nitrogen (LN, sand:soil 2:1). Measurements taken 6 weeks post-inoculation with rhizobia. Standard error of the differences among means (SED) and least significant difference (LSD) values included

The relationship between measured plant metrics (root and shoot dry weights, total plant weight and root:shoot ratio) and the number of nodules showed the strongest correlation coefficient (0.640) for shoot dry weight (Supplemental Table SM1). Further regression based evaluation, however, showed that shoot dry weight alone was not a strong predictor of nodule number (r2 39.5), but that shoot dry weight in conjunction with soil N as factors accounted for a substantial amount of the variation in nodule numbers (r2 71.9; Fig. 1). The N treatment with the smallest slope incline was the HN treatment, suggesting that nodule number is less crucial to plant dry weight in well fertilised conditions, while the MN treatment slope had the steepest incline.

Fig. 1

Chickpea shoot dry weight and soil N environment account for the number of nodules per plant. Linear regression plot of shoot dry weight (DW) for the six genotypes of chickpea grown with three levels of available N (HN (paddock soil + N), MN (paddock soil—N) and LN (sand:soil 2:1)) as predictors of the numbers of nodules per plant, r2 = 71.9, standard error of observations 9.97

Genotype dependent nodule formation in chickpea is not significantly related to root disease resistance rating

Given the significant effect of chickpea genotype on nodule formation and previous research suggesting that breeding for root disease resistance may impact nodulation [13], we expanded our study to a set of 29 chickpea genotypes sourced from either Australia or India with differing root disease resistance levels and type (i.e. quantitative for PRR vs. qualitative for FW) and grown at a single N level, results for all genotypes and the corresponding root disease resistance categories are presented in Supplemental Table SM2. Similar to earlier experiments with fewer varieties, the average number of nodules per plant differed significantly among chickpea genotypes (Fig. 2a). Comparison of the genotypes according to their respective FW or PRR disease resistance categories showed that for the FW resistance category comparisons (i.e. based on qualitative resistance) there were no significant differences in nodulation parameters and nodule performance (P > 0.05; Fig. 2b, Table 2a). Similarly, for the PRR resistance category comparisons (i.e. based on quantitative resistance) there was no evidence for significant differences in nodulation parameters and nodule performance (P > 0.05), although it was notable that the PRR moderately resistant group had lower nodule/g root values than the other two categories at the P < 0.10 level (Table 2b). Comparisons of the two geographic sources of genotypes showed contrasting trends (Table 2c). For the plant nodulation parameter, the Indian genotypes as a group had significantly more nodules per plant than the Australian bred genotypes (P < 0.01), but for the nodule performance results the Australian bred genotypes had significantly higher real-time N-fixation values per nodule than the Indian sourced genotypes (P < 0.01). To allow comparison of the range of genotype values in each geographic source group to previous studies, the ratio between the largest and smallest value (L:S ratio) was calculated for each group (Supplemental Table SM2). The Australian and Indian groups had similar a L:S ratio for the number of nodules per plant (values of 3.3 and 4.7, respectively) and nodules/g root (values of 4.1 and 4.2, respectively), but for shoot and root weights the Indian group had double the L:S ratio of the Australian group. Relationships between nodulation and plant growth related variates were evaluated across all 29 genotypes, the strongest correlation (0.568, P < 0.05) was between shoot DW and the number of nodules per plant (Supplemental Table SM3).

Fig. 2

Number of nodules formed per plant varies by genotype but is not linked to root disease resistance categories. A Average number of nodules formed per g root for 29 genotypes of chickpea, with SED, LSD and P values presented, N = 5 per genotype. Genotypes in bold were selected for further transcriptomic study. B Genotype root disease resistance categories by the average number of nodules per plant for 22 lines with differing Fusarium wilt (FW) resistance categories (MR moderately resistant, R resistant, AS asymptomatic) and seven lines with differing Phytophthora root rot (PRR) resistance (S susceptible, MS moderately susceptible, MR moderately resistant). Two FW MR genotypes also have resistance to Dry Root Rot (DRR) and two FW AS genotypes also have resistance to Botrytis Grey Mould (BGM). The LSD value of 21.7 is presented as an error bar, N = 5. All specific genotypes values in each figure are provided in Supplementary Table SM2

Table 2 Results for disease resistance category and genotype source comparisons for 29 chickpea genotypes for the variates of number of nodules/plant, number of nodules/g root, and ethylene production per nodule (ppm/g) for comparisons of a) Fusarium wilt (FW) disease resistance categories: moderately resistant (MR); resistant (R) and asymptomatic (AS); b) Phytophthora root rot (PRR) disease resistance categories: susceptible (S); moderately susceptible (MS), and moderately resistant (MR); and c) genotype source group of Indian (IND) and Australian genotypes (AUS). P values, standard error of differences (SED) and least significant difference (LSD) values

Chickpea genotypes exhibit similar transcriptomic profiles prior to rhizobial exposure

Six chickpea genotypes evenly split between high (ICC9002, ICC1510, ICC2072) and low (ICC11664, Jimbour, PBA Boundary) levels of nodulation were chosen based on the data above to study the difference in pre-symbiosis relative gene expression. Principle components analysis (PCA) of count data for all expressed genes (normalized and rlog transformed) from roots growing in the absence of rhizobia demonstrated that, with the exception of ICC1510 that partially separated from the group, all tested chickpea genotypes had a similar gene expression profile, with variability between genotypes being similar to that found within replicates (Fig. 3a). Similarly, relative gene expression in three different pathways associated with host rhizobial signalling and the initiation of nodulation—the CSSP, N-sensing/transport pathways, and flavonoid production pathways—was relatively consistent across all genotypes, regardless of each genotype’s nodulation potential (Fig. 3b-d). There are a few exceptions to this general observation, however. The genotypes forming more nodules also have lower levels of expression of chalcone synthase genes (Ca_13666, Ca_22255, Ca_22256, and Ca_22257; Fig. 3d), which are near the top of the flavonoid biosynthetic pathways, and a flavonoid methyltransferase (Ca_24494; Fig. 3d). Flavonoids, while being chemi-attractants for rhizobial bacteria, are also involved in plant defence pathways [25]. The homologue of DMI1 (Ca_02780), part of the CSSP, was most highly expressed in the two varieties forming the fewest nodules (Fig. 3b). Expression of the high affinity nitrate importers Ca_18399 and Ca_19740 (NRT2.1 and NRT1.2 homologues respectively) were also variable across the genotypes (Fig. 3c).

Fig. 3

Six chickpea genotypes with differing nodulation ability have similar transcriptomic profiles. A PCA plot of normalized and rlog transformed RNAseq count data for the roots of 6 chickpea genotypes grown in sterile soil. 95% confidence ellipses for each genotype are shown. B-D Normalized gene expression for genes known to be involved in B the common symbiotic signaling pathway (CSSP), C nitrogen transport and D flavonoid biosynthesis pathways in the roots of six chickpea genotypes grown in sterile soil. Genes indicated with arrows are (i) DMI1 (ii) high-affinity nitrate transporters, (iii) chalcone-synthase homologues or (iv) a flavonoid methyltransferase

Alterations to transcription in nodulation pathways following M. ciceri exposure do not mirror observed nodulation potentials in chickpea

Exposure to M. ciceri caused a number of genes to be differentially regulated in the chickpea roots 3 days post inoculation. Generally, ICC1510 was the most responsive genotype to the presence of the rhizobia with 2,383 genes being significantly differentially regulated (log2FC > 1, p > 0.05) compared to the sterile controls (Supplemental Figure 1). Early responses of legume roots to rhizobia have been well studied in model systems (Fig. 4a). Using specific genes in the chickpea CSSP and AON pathways, cytokinin and auxin responses, and flavonoid pathways, we considered the change in gene expression in genotypes with high or low nodulation after inoculation with rhizobia, compared to mock inoculated controls (Fig. 4b-d). High nodule forming ICC1510 and ICC2072 showed the most dramatic alterations in gene expression, particularly in the CSSP and AON pathways, however, the gene expression patterns showed high variability across genotypes with little to no link to nodulation potential. For example, Ca_22860, a homologue to Nod factor perception5 (NFR5), is up-regulated in ICC2072, but down-regulated in the two other high nodule forming genotypes ICC1510 and ICC9002. Nodule Inception (NIN; Ca_03225), which is a crucial transcriptional regulator controlling both nodule development and the AON pathway (Fig. 4a), also has inconsistent differential regulation among the genotypes, being upregulated in ICC2072, ICC16644 and Jimbour, but unchanged or downregulated in the other genotypes. Thus, there is no clear transcriptional pattern in the classical nodulation pathways distinguishing the response of high and low nodulation chickpea genotypes to rhizobia at this timepoint.

Fig. 4

Six chickpea genotypes with differing nodulation ability have inconsistent changes in gene expression in known nodulation pathways after inoculation with rhizobia. A Simplified schematic of early signalling events and interactions upon perception of rhizobia by legume roots. B-D Log2 fold change in gene expression for genes involved in B the common symbiotic signalling pathway (CSSP) and autoregulation of nodulation (AON), C cytokinin and auxin hormone responses and D flavonoid biosynthesis pathways in the roots of six chickpea genotypes 3 days after inoculation with M. ciceri

Altered expression of genes in signalling, cellular localization and N compound transport characterize high nodulating chickpea varieties

Given the lack of genes in classical nodulation pathways clearly differentiating the response of high or low nodule forming varieties to rhizobia, we used a supervised method (PLS-DA) to determine those genes whose expression most delineated the two nodulation phenotypes (Fig. 5). This resulted in 533 genes associated with the high nodulation varieties and showing, on average, increased expression in high nodulation genotypes as compared to low nodulation genotypes. An additional 861 genes were associated with low nodulation varieties (Supplemental Table SM4). GO enrichment analysis of these two groups of genes resulted in no significantly enriched GO terms for those genes associated with low nodulation, however, 42 terms were significantly enriched in genes associated with high nodulation (p < 0.01). These enriched GO terms generally fall under the categories of signalling, cellular localization, N compound transport and macromolecule or cellular component assembly (Table 3).

Fig. 5

A subset of genes expressed in chickpea roots after inoculation with Mesorhizobium ciceri differentiate varieties with high nodule formation from those with low nodule formation. A PLS-DA plot of normalized and rlog transformed RNAseq count data for the roots of 6 chickpea genotypes 3 days after inoculation with M. ciceri, with groupings assigned based on nodulation ability (high or low). 95% confidence ellipses for each nodulation level are shown. B Plot of loadings values from A for all expressed genes. Genes having coordinates with r > 0.013 were considered to contribute to the separation of high or low nodulation varieties

Table 3 Significantly enriched GO terms (p < 0.01) found in the set of 533 genes most differentiating gene expression in high nodulating varieties from low nodulating varieties following inoculation with rhizobia based on PLS-DA analysis


Biological N-fixation by rhizobia in legume crops, like chickpea, is important for plant nutrition and a reduced reliance on nitrogenous fertilizers. However, as our study reaffirms, the extent of root nodule formation is dependent both on the genotype of the plant and its abiotic environment [8, 26]. We demonstrate here that nodulation varies in a genotype-dependent manner across 29 genotypes of chickpea, with genotype specific responses to N-availability being demonstrated in an experiment using a smaller set of 6 varieties. We show that nodule number, while linked to plant growth traits such as above ground biomass, is not consistently related to the disease resistance rating of the genotype. This latter observation was most marked in plants known to encode qualitative single pathogen resistance. Finally, using transcriptomic analysis during the early stages of rhizobial colonization, we show that genotypes with opposing levels of nodule formation, while having generally overlapping gene expression profiles, show a sub-set of genes involved in signalling and transport pathways that vary by nodulation level, suggesting a potential role for these genes in chickpea rhizobial symbiosis.

It is well established that endemic Mesorhizobium species can affect the extent of nitrogen fixation by inoculated species and the soil environment effects rhizobium populations [6, 27]. We therefore carried out genotype based studies in a sterile non-field environment to provide control over potentially confounding factors. These experiments demonstrated that, despite a strong effect of plant-available N in influencing nodulation in legumes, not all plant genotypes have the same response to changing N levels in the soil. Among the six genotypes tested, the number of nodules differed significantly, both between genotypes and across each of the three soil N levels tested. Similar genotype specific differences in nodulation and N-fixation have previously been reported when considering restricted sets of chickpea genotypes (n < 7) [9, 28]. Given the role of this symbiosis in supplementing plant N, it could be assumed that the lowest concentration of plant-available N in the soil would be most favourable to the establishment of rhizobial interaction with the plant host. However, in our study and in others, nodule numbers were found to be highest at moderate levels of N and lowest at either reduced or high N [29]. As the formation of nodules is an energetically expensive investment on the part of the plant, it is possible that at very low soil N concentration, the plant lacks the resources to form large numbers of nodules. The perception of nitrogen status and availability by the plant host is dependent on a number of sensory systems and signalling pathways, which in turn affect pathways necessary for nodulation [30]. For example, the presence of high levels of nitrate in soils can induce the expression of CLE peptides in legume roots, triggering the AON pathway, thereby suppressing nodulation [22, 23]. While we were unable to test these mechanisms, sensitivity differences among chickpea genotypes to soil plant-available N explained our observed genotype by soil N results. Further investigation of these pathways may result in the identification of genotypic markers that predict low plant sensitivity to nitrate and, therefore, may promote high levels of nodulation in agricultural fields with high residual levels of plant-available N.

In our previous work, in which we analysed a small set of moderately PRR resistant varieties, we demonstrated a possible link between quantitative resistance to PRR, lower nodule numbers and N-fixation [13]. This was consistent with other reports of how certain plant genotypes can simultaneously impact pathogenic and mutualistic symbioses [31,32,33,34,35,36,37]. However, while the moderately resistant PBA HatTrick again had lower nodulation values than PRR-susceptible Sonali in this study, other PRR susceptible/moderately susceptible varieties such as PBA Boundary formed fewer nodules again. When we considered different genotypes with a range of FW disease resistance categories, we found both exceptional and poor performance of specific genotypes within each resistance category in terms of the number of nodules formed on the root system. This lack of consistent trend between disease resistance and nodulation within the FW genotype set is likely due to the fact that resistance to Fusarium oxysporum f.sp. ciceris is qualitatively based on a small number of disease-specific genes [12]. Therefore, the likelihood that one of these genes may also negatively affect beneficial rhizobia was low and the finding that there was no consistent pattern between FW resistance and repression of nodulation potential was not unexpected. Together, these findings demonstrate that, while disease resistance may affect nodulation, it is not the only indicator of nodulation potential.

Improved selection methods for high nodulation legumes are required. The association between nodulation, or nodule N-fixation, and agronomic traits such as root weight and shoot weight has led to suggestions that breeders can use these in-direct indicators of N associated traits, to achieve ‘simultaneous improvement’ within a particular legume crop [8, 38]. Simultaneous improvement, defined as the selection of one plant trait linked pleiotropically to multiple valued agronomic phenotypes, as a method for varietal development is undoubtedly highly important, and would improve efficiency in breeding programs. It is also clear, however, that an improved understanding of the genetic basis of nodule formation and enhancement of N-fixation across a broader range of legume crops is required before more targeted selection techniques can easily be adopted during the breeding process. Chickpea, known to have relatively limited genomic diversity [39, 40] (Fig. 4), may seem unlikely to show enough differences in the establishment of rhizobial symbiosis to be a candidate crop where this may be achieved. However, our study using putatively genetically similar FW resistant genotypes, and another by Biabani et al. [10], using genetically diverse chickpea genotypes from 12 countries, both found large genotype based differences in nodulation.

In a consideration of gene expression in known nodulation pathways for genotypes with high or low nodulation potential, most genes were consistently expressed across the genotypes in the absence of rhizobia, however, there were notable differences between the high and low nodulating genotypes. For example, chalcone synthase gene expression was lower in genotypes forming the greatest number of nodules without rhizobia present. Chalcones sit at the top of the flavonoid pathway and play a role both in plant defence and as chemi-attractants for rhizobial bacteria [25] and proper balance of chalcone synthase activity is necessary for nodulation [41,42,43]. Upon addition of rhizobia, gene expression in known nodulation pathways was not altered in a manner consistent with nodulation potential. While this may, in part, stem from the early time point used, earlier research suggests that chickpea is already responding to rhizobia with altered gene transcription 3 days after inoculation in this system [13, 44]. The expression of Nodule inception (NIN) is down-regulated in two of the high nodulation genotypes, despite commonly being up-regulated upon perception of rhizobia in other model legume systems [45]. This transcription factor, while essential for nodule organogenesis, also sits at the head of the AON pathway and its activity can inhibit nodule initiation [46, 47]. Its activity is also tissue specific with distinct functions in root epidermis versus cortex [45], therefore it is possible that while the expression of the gene is generally down regulated in some genotypes, it could still be locally accumulated in specific tissues. Overall, however, the gene expression patterns seen in model legumes are not consistently reflected in the chickpea roots tested here. Thus, rather than one key pathway across which differential gene regulation describes differences in nodulation between genotypes, cross talk between multiple signalling pathways based on differences in a small set of genes may co-ordinately explain the observed nodulation patterns. Further examination of gene expression following inoculation with rhizobia resulted in a set of 533 genes associated with the high nodulation phenotype. Although we do not know at this time if and how each of these integrate into the known pathways supporting rhizobial symbiosis, they provide new avenues for future research and may provide the necessary markers for simultaneous improvement during the breeding of chickpea.


Overall, the results of this paper give evidence to the possible scope for using genotypic screens in non-model crops such as chickpea to identify novel factors affecting nodule formation and function. Despite a narrow genetic pool, we have highlighted specific, discrete aspects by which chickpea genotype and transcriptional profiles impact rhizobial symbiosis. In future extending the results here may include consideration of the differences in the flavonoid profile produced by chickpea varieties with high and low nodule formation to identify specific compounds linked to nodule numbers. Further, the pathways controlled by the identified genes associated with the high nodulation should be refined with a larger set of genotypes and considered for their tractability as robust markers of nodulation potential. Together, results proceeding from the work described here as well as these proposed studies will generate the genetic framework necessary for the simultaneous improvement of non-model pulse crops.


Chickpea growth experiments

Six varieties of chickpea, Sonali, PBA HatTrick, Jimbour, Moti, Kyabra, and Yorker, were grown in three different soil types according to the methods of Plett et al. [13]. The source of all seed for the experiments was the New South Wales Department of Primary Industries led Pulse Breeding Australia chickpea breeding program (Tamworth, NSW, Australia). Briefly, seeds of each variety were surface sterilized with 10% bleach prior to planting in sterilized vermiculite. The seeds were allowed to germinate and grow for 2 weeks prior to transplanting into the test soils. Each plant was inoculated upon transplant with 10,000 CFU of M. ciceri (CC1192) and grown for another 6 weeks in growth chambers (18/10 °C day/night temperature; 12 h light cycle; BioChambers, Saskatchewan, Canada). The soils used for the experiments were all sterilized by gamma irradiation prior to use and included chickpea paddock soil with N fertilization (HN; 71.8 mg/kg NO3 + 7.5 mg/kg NH4), chickpea paddock soil without N fertilization (MN; 65 mg/kg NO3 + 7.5 mg/kg NH4) and paddock soil cut two times with sand (LN; 21.7 mg/kg NO3 + 2.5 mg/kg NH4). Plant available N was determined at harvest using KCl extraction as per Plett et al. [13]. Plants were fertilized three times after transplanting with an N free fertilizer and HN treatments were given an additional N fertilization treatment (3 × 1 mL of 0.1 M KNO3).

In the second series of experiments, 22 Desi lines were selected from International Crops Research Institute for the Semi-Arid Tropics (ICRISAT) mini-core collection of 211 accessions. The selections were based on the disease resistance information, seed colour and size, and origin information from Pande et al. [48]. To minimise potential origin based variation, all ICRISAT lines were chosen to originate from India. In addition only Desi lines and lines that were non-segregating for seed colour were selected. The final selection criteria was for Fusarium wilt resistance, with 6 lines rated as asymptomatic, 6 as moderately resistant and 6 as resistant to Fusarium wilt, and four lines with FW resistance that included multiple disease resistance as shown in Table SM-1. For comparison to the ICRISAT Fusarium wilt lines, from Australian sources, seven Australian bred Desi varieties were selected with differing PRR disease resistance ratings, two susceptible (Sonali, PBA Boundary), three moderately susceptible (Kyabra, Moti, Jimbour) and two moderately resistant (PBA HatTrick and Yorker). The chickpea plants were grown as for the first experiment but in one soil type (Sand:paddock soil 1:1, sterilized by gamma irradiation, 16 mg/kg NO3 + 1.2 mg/kg NH4) and given N free fertilizer only.

Chickpea experiment harvest and analysis

At harvest for both experiments, chickpea root systems were rinsed of soil and the number of nodules per root system were counted. About 10–20 nodules were carefully removed from the root system, weighed and placed in a sealed 20-mL GC vial for an acetylene reduction assay to assess nodule activity through the production of ethylene as per Plett et al. [13]. The plant above and below ground portions were separated and dried at 400C prior to weighing for biomass.

RNA sequencing and analysis

Chickpea seeds from three genotypes forming low numbers of nodules (Jimbour, Boundary, ICC11664) and three genotypes with high nodulation (ICC9002, ICC2072 and ICC1510) were sterilized in 10% bleach, rinsed in sterile water and planted into vermiculite for 2 weeks as above. Seedlings were transplanted to sterilized soil (same as 29 genotype experiment) and acclimatised to the new soil for a day prior to inoculation. Actively growing cultures of M. ciceri (~ 10,000 CFU per plant) was applied to half of the chickpeas while a mock solution without rhizobia was applied to the other half. Chickpea roots were harvested and frozen in liquid nitrogen 3 days after addition of live M. ciceri cultures or mock controls.

RNA was extracted from the chickpea roots (3 biological replicates per condition) using the Isolate II RNA Plant kit (Bioline, London, U.K.) according to manufacturer’s instructions. RNA libraries were constructed using the SENSE total RNA-Seq library prep kit (Lexogen, Vienna, Austria) and NovaSeq sequencing was performed at the Ramaciotti Centre for Genomics (University of New South Wales, Sydney, Australia). Raw reads from the RNA-sequencing were trimmed and mapped to the C. arietinum genome (desi type—ICC4958; [49]) using CLC Workbench v.10. Genes that were not expressed or lowly expressed (< 5 counts on average across each set of 3 replicates) were eliminated from the data set. Raw gene count data was normalized using the DESeq2 package in R (DeSEQ2 v.1.24.0; Love et al. [50]; R v.1.2.1335; [51]). For differential gene expression analysis, log2 fold change was calculated for the rhizobia inoculated chickpea roots as compared to the uninoculated control roots. Adjusted p-values were calculated using a Benjamini–Hochberg FDR correction. Genes up- or down- regulated more than 2 times and with p < 0.05 were considered significantly differentially regulated. Gene annotations were assigned based on the desi-type chickpea annotation [52].

Identification of genes in specific flavonoid production, nodulation or nitrogen transportation pathways was accomplished either by directly searching the C. arietinum genome annotation, or by BLAST analysis of corresponding protein sequences from the model legume M. truncatula or from A. thaliana against the C. arietinum (ICC4958) proteome [same ref as above] (Supplemental Table SM5). Determination of flavonoid biosynthesis enzymes was based off of Gifford et al. [16]. M. truncatula and A. thaliana protein sequences were obtained from an open commons database [53].

Statistical analyses

Data from the two chickpea controlled environment experiments were analysed with ANOVA, using an appropriate transformation prior to the analysis, to stabilise the variances. Means comparison tests used the Fisher’s LSD method. Correlations between all variates in both studies were completed, in addition selected variates were compared with linear regression. For the 29-genotype experiment, additional comparisons of values between the Indian and Australian genotype groups were made using the ratio between the highest and lowest value in each group for each variate was calculated and called the largest to smallest (L:S) ratio. All analyses were carried out in Genstat.

Principal components analysis (PCA) and partial least squares discriminate analysis (PLS-DA) was conducted in R using the mixOmics package [54]. Input transcriptomic data was first normalized and rlog transformed using the R DeSEQ2 package. Loadings data for all genes in the PLS-DA analysis were plotted and those genes falling outside of a radius of 0.013 from the origin were considered to be contributing towards the separation of high and low nodulation phenotypes. GO enrichment analyses were conducted in R using the topGO package [55] and GO annotations based on the C. arietinum (ICC4958) genome. GO terms were assigned as significantly enriched based on a classic Fisher test (p < 0.01). Heat maps of normalized count data and log2 fold change data were generated using the Morpheus online tool [56].

Availability of data and materials

The transcriptomic datasets generated and analysed during the current study along with count data tables are available at GEO submission number GSE162321:



Autoregulation of nodulation


Colony Forming Units


Common symbiotic signalling pathway


Dry weight


Fusarium wilt


Gene ontology


Least significant difference




Principal Components analysis


Partial Least Squares Discriminate analysis


Phytophthora root rot


Standard error of difference


  1. 1.

    Peoples MB, Craswell ET. Biological nitrogen fixation: Investments, expectations and actual contributions to agriculture. Plant Soil. 1992;141(1):13–39.

    CAS  Article  Google Scholar 

  2. 2.

    Herridge DF, Peoples MB, Boddey RM. Global inputs of biological nitrogen fixation in agricultural systems. Plant Soil. 2008;311(1–2):1–18.

    CAS  Article  Google Scholar 

  3. 3.

    Knights EJ, Açıkgöz N, Warkentin T, Bejiga G, Yadav SS, Sandhu JS. Chickpea breeding and management. In: Yadav SS, Redden RJ, Chen W, Sharma B, editors. Area, production and distribution. Wallingford: CABI; 2007. p. 167–78.

    Google Scholar 

  4. 4.

    Schwenke GD, Peoples MB, Turner GL, Herridge DF. Does nitrogen fixation of commercial, dryland chickpea and faba bean crops in north-west New South Wales maintain or enhance soil nitrogen? Aust J Exp Agric. 1998;38(1):61–70.

    Article  Google Scholar 

  5. 5.

    Gan Y, Johnston AM, Knight JD, McDonald C, Stevenson C. Nitrogen dynamics of chickpea: effects of cultivar choice, N fertilization, rhizobium inoculation, and cropping systems. Can J Plant Sci. 2010;90(5):655–66.

    CAS  Article  Google Scholar 

  6. 6.

    Greenlon A, Chang PL, Damtew ZM, Muleta A, Carrasquilla-Garcia N, Kim D, Nguyen HP, Suryawanshi V, Krieg CP, Yadav SK, et al. Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria. Proc Natl Acad Sci U S A. 2019;116(30):15200–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Doughton JA, Vallis I, Saffigna PG. Nitrogen fixation in chickpea. I. Influence of prior cropping or fallow, nitrogen fertilizer and tillage. Aust J Agric Res. 1993;44:1403–13.

    Article  Google Scholar 

  8. 8.

    Dwivedi SL, Sahrawat KL, Upadhyaya HD, Mengoni A, Galardini M, Bazzicalupo M, Biondi EG, Hungria M, Kaschuk G, Blair MW, et al. Advances in host plant and rhizobium genomics to enhance symbiotic nitrogen fixation in grain legumes. In: Sparks DL, editor. Advances in agronomy, vol. 129. San Diego: Elsevier Academic Press Inc; 2015. p. 1–116.

    Google Scholar 

  9. 9.

    Abi-Ghanem R, Carpenter-Boggs L, Smith J, Vandemark G. Nitrogen fixation by US and middle eastern chickpeas with commercial and wild middle eastern inocula. IRSN Soil Sci. 2012.

  10. 10.

    Biabani A, Carpenter-Boggs L, Coyne CJ, Taylor L, Smith JL, Higgins S. Nitrogen fixation potential in global chickpea mini-core collection. Biol Fertil Soils. 2011;47(6):679–85.

    Article  Google Scholar 

  11. 11.

    Singh KB, Malhotra RS, Halila MH, Knights EJ, Verma MM. Current status and future strategy in breeding chickpea for resistance to biotic and abiotic stresses. Euphytica. 1994;73(1–2):137–49.

    Article  Google Scholar 

  12. 12.

    Infantino A, Kharrat M, Riccioni L, Coyne CJ, McPhee KE, Grunwald NJ. Screening techniques and sources of resistance to root diseases in cool season food legumes. Euphytica. 2006;147(1–2):201–21.

    Article  Google Scholar 

  13. 13.

    Plett JM, Plett KL, Bithell SL, Mitchell CM, Moore K, Powell JR, Anderson IC. Improved Phytophthora resistance in commercial chickpea (Cicer arietinum) varieties negatively impacts symbiotic gene signalling and symbiotic potential in some varieties. Plant Cell Environ. 2016;39:1858–69.

  14. 14.

    Varshney RK, Song C, Saxena RK, Azam S, Yu S, Sharpe AG, Cannon S, Baek J, Rosen BD, Tar’an B, et al. Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat Biotechnol. 2013;31(3):240–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  15. 15.

    Young ND, Debellé F, Oldroyd GED, Geurts R, Cannon SB, Udvardi MK, Benedito VA, Mayer KFX, Gouzy J, Schoof H, et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature. 2011;480(7378):520–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Gifford I, Battenberg K, Vaniya A, Wilson A, Tian L, Fiehn O, Berry AM. Distinctive patterns of flavonoid biosynthesis in roots and nodules of Datisca glomerata and Medicago spp. revealed by metabolomic and gene expression profiles. Front Plant Sci. 2018;9:1463.

    PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Cooper JE. Early interactions between legumes and rhizobia: disclosing complexity in a molecular dialogue. J Appl Microbiol. 2007;103(5):1355–65.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  18. 18.

    Oldroyd GED, Murray JD, Poole PS, Downie JA. The rules of engagement in the legume-rhizobial symbiosis. In: Bassler BL, Lichten M, Schupbach G, editors. Annual review of genetics, vol. 45. Palo Alto: Annual Reviews; 2011. p. 119–44.

    Google Scholar 

  19. 19.

    Buhian WP, Bensmihen S. Mini-review: nod factor regulation of phytohormone signaling and homeostasis during rhizobia-legume symbiosis. Front Plant Sci. 2018;9:1247.

    PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Ferguson BJ, Mens C, Hastwell AH, Zhang M, Su H, Jones CH, Chu X, Gresshoff PM. Legume nodulation: the host controls the party. Plant Cell Environ. 2019;42(1):41–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Soyano T, Shimoda Y, Hayashi M. NODULE INCEPTION antagonistically regulates gene expression with nitrate in Lotus japonicus. Plant Cell Physiol. 2015;56(2):368–76.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Reid DE, Ferguson BJ, Hayashi S, Lin Y-H, Gresshoff PM. Molecular mechanisms controlling legume autoregulation of nodulation. Ann Bot. 2011;108(5):789–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Murray JD, Liu CW, Chen Y, Miller AJ. Nitrogen sensing in legumes. J Exp Bot. 2017;68(8):1919–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Haskett T, Wang P, Ramsay J, O’Hara G, Reeve W, Howieson J, Terpolilli J. Complete genome sequence of Mesorhizobium ciceri strain CC1192, an efficient nitrogen-fixing microsymbiont of Cicer arietinum. Genome Announc. 2016;4(3):e00516.

    PubMed  PubMed Central  Google Scholar 

  25. 25.

    Hassan S, Mathesius U. The role of flavonoids in root-rhizosphere signalling: opportunities and challenges for improving plant-microbe interactions. J Exp Bot. 2012;63(9):3429–44.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  26. 26.

    Peoples MB, Bowman AM, Gault RR, Herridge DF, McCallum MH, McCormick KM, Norton RM, Rochester IJ, Scammell GJ, Schwenke GD. Factors regulating the contributions of fixed nitrogen by pasture and crop legumes to different farming systems of eastern Australia. Plant Soil. 2001;228(1):29–41.

    CAS  Article  Google Scholar 

  27. 27.

    Elias NV, Herridge DF. Naturalised populations of mesorhizobia in chickpea (Cicer arietinum L.) cropping soils: effects on nodule occupancy and productivity of commercial chickpea. Plant Soil. 2015;387(1–2):233–49.

    CAS  Article  Google Scholar 

  28. 28.

    Rao DLN, Giller KE, Yeo AR, Flowers TJ. The effects of salinity and sodicity upon nodulation and nitrogen fixation in chickpea (Cicer arietinum). Ann Bot. 2002;89(5):563–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Xia X, Ma C, Dong S, Xu Y, Gong Z. Effects of nitrogen concentrations on nodulation and nitrogenase activity in dual root systems of soybean plants. Soil Sci Plant Nutr. 2017;63(5):470–82.

    CAS  Article  Google Scholar 

  30. 30.

    Gent L, Forde BG. How do plants sense their nitrogen status? J Exp Bot. 2017;68(10):2531–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Ané J-M, Lévy J, Thoquet P, Kulikova O, de Billy F, Penmetsa V, Kim D-J, Debellé F, Rosenberg C, Cook DR, et al. Genetic and cytogenetic mapping of DMI1, DMI2, and DMI3 genes of Medicago truncatula involved in nod factor transduction, nodulation, and mycorrhization. Mol Plant-Microbe Interact. 2002;15(11):1108–18.

    PubMed  Article  PubMed Central  Google Scholar 

  32. 32.

    Arrighi J-F, Barre A, Ben Amor B, Bersoult A, Soriano LC, Mirabella R, de Carvalho-Niebel F, Journet E-P, Ghérardi M, Huguet T, et al. The Medicago truncatula lysin motif-receptor-like kinase gene family includes NFP and new nodule-expressed genes. Plant Physiol. 2006;142(1):265.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Genre A, Ortu G, Bertoldo C, Martino E, Bonfante P. Biotic and abiotic stimulation of root epidermal cells reveals common and specific responses to abuscular mycorrhizal fungi. Plant Physiol. 2009;149(3):1424.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Wang E, Schornack S, Marsh John F, Gobbato E, Schwessinger B, Eastmond P, Schultze M, Kamoun S, Oldroyd Giles ED. A common signaling process that promotes mycorrhizal and oomycete colonization of plants. Curr Biol. 2012;22(23):2242–6.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Rey T, Chatterjee A, Buttay M, Toulotte J, Schornack S. Medicago truncatula symbiosis mutants affected in the interaction with a biotrophic root pathogen. New Phytol. 2015;206(2):497–500.

    PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Rey T, Nars A, Bonhomme M, Bottin A, Huguet S, Balzergue S, Jardinaud M-F, Bono J-J, Cullimore J, Dumas B, et al. NFP, a LysM protein controlling Nod factor perception, also intervenes in Medicago truncatula resistance to pathogens. New Phytol. 2013;198(3):875–86.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Zhang X, Dong W, Sun J, Feng F, Deng Y, He Z, Oldroyd GED, Wang E. The receptor kinase CERK1 has dual functions in symbiosis and immunity signalling. Plant J. 2015;81(2):258–67.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Wong-Bajracharya J, Plett J. Root renovation: how an improved understanding of basic root biology could inform the development of elite crops that foster sustainable soil health. Funct Plant Biol. 2019;46(7):597–612.

    Article  Google Scholar 

  39. 39.

    Singh R, Sharma P, Varshney RK, Sharma SK, Singh NK. Chickpea improvement: role of wild species and genetic markers. Biotechnol Genet Eng Rev. 2008;25(1):267–314.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  40. 40.

    Abbo S, Berger J, Turner N. Evolution of cultivated chickpea: four bottlenecks limit diversity and constrain adaptation. Funct Plant Biol. 2003;30(10):1081–7.

    PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Cooper J. Multiple responses of rhizobia to flavonoids during legume root infection. Adv Bot Res. 2004;41:1–62.

    CAS  Article  Google Scholar 

  42. 42.

    Zuanazzi JAS, Clergeot PH, Quirion J-C, Husson H-P, Kondorosi A, Ratet P. Production of Sinorhizobium meliloti nod gene activator and repressor flavonoids from Medicago sativa roots. Mol Plant Microbe Interact. 1998;11(8):784–94.

    CAS  Article  Google Scholar 

  43. 43.

    Abdel-Lateif K, Vaissayre V, Gherbi H, Verries C, Meudec E, Perrine-Walker F, Cheynier V, Svistoonoff S, Franche C, Bogusz D, et al. Silencing of the chalcone synthase gene in Casuarina glauca highlights the important role of flavonoids during nodulation. New Phytol. 2013;199(4):1012–21.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Plett JM, Solomon J, Snijders F, Marlow-Conway J, Plett KL, Bithell SL. Order of microbial succession affects rhizobia-mediated biocontrol efforts against Phytophthora root rot. Microbiol Res. 2020;242:126628.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  45. 45.

    Vernié T, Kim J, Frances L, Ding Y, Sun J, Guan D, Niebel A, Gifford ML, de Carvalho-Niebel F, Oldroyd GED. The NIN transcription factor coordinates diverse nodulation programs in different tissues of the Medicago truncatula root. Plant Cell. 2015;27(12):3410–24.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  46. 46.

    Yoro E, Suzaki T, Toyokura K, Miyazawa H, Fukaki H, Kawaguchi M. A positive regulator of nodule organogenesis, NODULE INCEPTION, acts as a negative regulator of rhizobial infection in Lotus japonicus. Plant Physiol. 2014;165(2):747–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Soyano T, Hirakawa H, Sato S, Hayashi M, Kawaguchi M. NODULE INCEPTION creates a long-distance negative feedback loop involved in homeostatic regulation of nodule organ production. Proc Natl Acad Sci. 2014;111(40):14607–12.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Pande S, Kishore GK, Upadhyaya HD, Rao JN. Identification of sources of multiple disease resistance in mini-core collection of chickpea. Plant Dis. 2006;90(9):1214–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  49. 49.

    Parween S, Nawaz K, Roy R, Pole AK, Suresh BV, Misra G, Jain M, Yadav G, Parida SK, Tyagi AK, et al. An advanced draft genome assembly of a desi type chickpea (Cicer arietinum L.). Sci Rep. 2015;5:14.

    Article  CAS  Google Scholar 

  50. 50.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):38.

    Article  CAS  Google Scholar 

  51. 51.

    R-Core-Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2016.

    Google Scholar 

  52. 52.

    Anon. Legume Information System. Accessed 9 Sept 2019.

  53. 53.

    Anon. Phytozome. Accessed 9 Sept 2019.

  54. 54.

    Le Cao K-A, Rohart F, Gonzalez I, Dejean S, Gautier B, Bartolo F, Monget P, Coquery J, Yao F, Liquet B. mixOmics: Omics data integration project. R package version 6.1.1. 2016.

    Google Scholar 

  55. 55.

    Alexa A, Rahnenfuhrer J. topGO: enrichment analysis for gene ontology. R package version 2.42.0. 2020.

    Google Scholar 

  56. 56.

    Anon. Morpheus. Accessed 26 Oct 2019.

Download references


The authors acknowledge the Ramaciotti Centre for Genomics at University of New South Wales for RNA sequencing and N. Plett and B. Laugesen for technical assistance.


We would like to acknowledge co-funded research by NSW DPI and the Grains Research and Development Corporation (GRDC) that supported SLB under project DAN00213. We also acknowledge funding support from the former Department of Agriculture, Fisheries and Forestry to JMP from a Science and Innovation grant.

Author information




KLP, SLB and JMP designed the research, KLP, AD and JMP performed experiments and all authors were involved in data analysis. KLP and SLB wrote the manuscript and all authors read and approved the final version.

Corresponding author

Correspondence to Jonathan M. Plett.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Figure SM1.

Number of genes significantly differentially expressed in chickpea roots with differing nodulation after inoculation with Mesorhizobium ciceri. Number of genes significantly up- (blue) or down- (orange) regulated (log2FC >1, p < 0.05) in the roots of six chickpea genotypes three days after inoculation with M. ciceri as compared to sterile controls. Average number of nodules formed by plants at harvest are indicated (grey line). Supplemental Table SM1. Correlation coefficient results for the six genotype experiment with soil nitrogen treatments for measurements after four weeks of plant height, leaf and branch number, number of nodules, number of nodules per g of root, root dry weight (DW), root to shoot DW ratio, shoot DW and total plant DW. *indicates P < 0.05 for a two-sided test of correlations different from zero. Supplemental Table SM2. Varietal comparison results for the controlled environment study of 29 chickpea genotypes with a range of root disease resistance ratings, for the variates: number of nodules (no. nod.); root (g.), shoot (g.), total plant (tot. pl, g) dry weights (DW); root:shoot (R:S) DW ratio per plant; number of nodules/g of root (no. nod/g); and nodule ethylene production (eth ppm per g nod). Root disease resistance abbreviations Fusarium wilt (FW), Phytophthora root rot (PRR), resistant (R), asymptomatic (AS), moderately resistant (MR), Dry root rot (DRR), Botrytis grey mould (BGM), moderately susceptible (MS), susceptible (S). P values, standard error of differences (SED) and least significant difference (LSD) values from ANOVA are included. For the Indian (ICC lines) and Australian (all non-ICC lines; PRR) lines the ratio between the highest and lowest value (L:S ratio) for each variate is also presented. *Genotypes studied for transcriptomic responses. Supplemental Table SM3. Correlation coefficient results for the 29 genotype experiment for measurements of number of nodules, number of nodules per g of root, average weight of nodules, ethylene production per gram of nodule, root dry weight (DW) shoot DW, root:shoot DW ratio and total plant DW. *indicates P < 0.05 for a two-sided test of correlations different from zero. Supplemental Table SM4. List of chickpea genes contributing to the separation of high and low nodulating chickpea genotypes in PLS-DA analysis (based on loadings coordinates with r > 0.013) of gene expression 3 days after inoculation with Mesorhizobium ciceri. Loadings coordinates on components 1 and 2 and radius are indicated, along with the gene and GO annotation. Supplemental Table SM5. List of chickpea genes and annotations in known pathways for early signalling and responses to rhizobia. Annotations were assigned based on either protein BLAST to known Medicago truncatula or Arabidopsis homologues if indicated, or gene assignments based on the Cicer arietinum (ICC4958) genome. * taken from Plett et al. 2016 [13] **genes identified based on Gifford et al. 2018 [16].

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Plett, K.L., Bithell, S.L., Dando, A. et al. Chickpea shows genotype-specific nodulation responses across soil nitrogen environment and root disease resistance categories. BMC Plant Biol 21, 310 (2021).

Download citation


  • Genotype specificity
  • Soil N environment
  • Rhizobia
  • Nodulation regulation
  • Nodule nitrogen production
  • Transcriptomic analysis
  • Root disease resistance