Skip to main content
  • Research article
  • Open access
  • Published:

Analysis of root proteome unravels differential molecular responses during compatible and incompatible interaction between chickpea (Cicer arietinum L.) and Fusarium oxysporum f. sp. ciceri Race1 (Foc1)



Vascular wilt caused by Fusarium oxysporum f. sp. ciceri Race 1 (Foc1) is a serious disease of chickpea (Cicer arietinum L.) accounting for approximately 10-15% annual crop loss. The fungus invades the plant via roots, colonizes the xylem vessels and prevents the upward translocation of water and nutrients, finally resulting in wilting of the entire plant. Although comparative transcriptomic profiling have highlighted some important signaling molecules, but proteomic studies involving chickpea-Foc1 are limited. The present study focuses on comparative root proteomics of susceptible (JG62) and resistant (WR315) chickpea genotypes infected with Foc1, to understand the mechanistic basis of susceptibility and/or resistance.


The differential and unique proteins of both genotypes were identified at 48 h, 72 h, and 96 h post Foc1 inoculation. 2D PAGE analyses followed by MALDI-TOF MS and MS/MS identified 100 differentially (>1.5 fold<, p < 0.05) or uniquely expressed proteins. These proteins were further categorized into 10 functional classes and grouped into GO (gene ontology) categories. Network analyses of identified proteins revealed intra and inter relationship of these proteins with their neighbors as well as their association with different defense signaling pathways. qRT-PCR analyses were performed to correlate the mRNA and protein levels of some proteins of representative classes.


The differential and unique proteins identified indicate their involvement in early defense signaling of the host. Comparative analyses of expression profiles of obtained proteins suggest that albeit some common components participate in early defense signaling in both susceptible and resistant genotypes, but their roles and regulation differ in case of compatible and/or incompatible interactions. Thus, functional characterization of identified PR proteins (PR1, BGL2, TLP), Trypsin protease inhibitor, ABA responsive protein, cysteine protease, protein disulphide isomerase, ripening related protein and albumins are expected to serve as important molecular components for biotechnological application and development of sustainable resistance against Foc1.


Plants are often challenged by different types of biotic and abiotic stress factors. Their immobile nature precludes escape from these stress causing agents. Therefore, they possess preformed and inducible defensive strategies to overcome these stresses. In most cases, the host arrests the invading rival at the site of penetration [1]. Such immune response adapted by the host is termed as pattern triggered immunity (PTI) which include reprogramming of host cellular metabolism, reinforcement of cell wall by callose occlusions and production of antimicrobial compounds that act directly to prevent pathogen invasion [2, 3]. However, in some selected cases the invading pathogens secrete effector molecules that try to overcome host immunity, which in the absence of cognate host resistant protein/proteins (R-proteins) lead to effector triggered susceptibility (ETS) [4]. On the other hand, in the presence of cognate R protein/proteins the host mounts a defense response of much greater amplitude known as the effector triggered immunity (ETI), which largely overlaps with that of PTI [5]. However, the defense mechanisms of both PTI and/or ETI are regulated by altered protein synthesis and their time dependent degradation. Hence, qualitative and quantitative changes in protein levels are believed to be probable indicators of the ultimate outcome of any plant-pathogen interaction.

Amongst agronomical important crop plants, legume crops are known for their nutritive value that play very important roles in human nutrition as well as serve as supplement to improve growth of livestock [6]. Besides, they also fix atmospheric nitrogen enhancing soil fertility and boosting the yield of subsequently grown crops [7]. These crops are equally vulnerable towards pathogen. But studies on the molecular interaction involving legume-pathogen case study are significantly limited. Chickpea is the third most important legume crop in the world and the most important one in India (FAO). It is a rich source of digestible protein, and hence is considered globally as a valuable crop. However it is found that it accounts for 10 to 15% of yield loss worldwide by wilt causing fungus Fusarium oxysporum f. sp. ciceri (Foc). This seed or soil borne fungus has two different pathotypes, a yellowing pathotype and a wilt causing pathotype [8]. Amongst eight pathogenic races of Foc (Races 0, 1, 1B/C, 2, 3, 4, 5, and 6) Race 1, known to show a wide geographic distribution throughout India has received major scientific concern. The fungus invades the plant via roots, colonizes the xylem vessels and prevents the upward translocation of water and nutrients, finally resulting in wilting [9]. Yellowing of rootlets, chlorosis of basal leaflets and drooping of lower branches are the initial symptoms of pathogenic infection in chickpea plants [10]. Until recently, Fusarium wilt was being managed by resistance breeding programs. But the main hurdle faced by plant breeders was pathogenic variability and mutability that resulted in breakdown of natural resistance over prolonged period of time and generations [11]. Although, chemical fungicides are used as alternatives under such circumstances [12], but high cost and environmental safety issues are known to raise social concerns regarding its long term utilization. Therefore a proper understanding of the molecular mechanism involved in chickpea-Foc1 interaction could suggest effective measures for developing sustainable resistance.

Previous studies conducted on understanding the molecular interaction of chickpea- Fusarium oxysporum was based mainly on transcriptomic studies taking lead from the model plant Arabidopsis thaliana and tomato (Solanum lycopersicum) [13, 14]. Moreover, previously published histopathological reports suggested that Foc1 enters the roots through the breaches of root hairs and colonizes the xylem vessel of compatible host at about 4dpi (days post inoculation). Rapid establishment of Foc1 coupled with massive tissue disintegration led to the total collapse of root architecture ultimate causing wilting of susceptible plants at about 12dpi, whereas resistant plants showed minimal signs of stress even at later stages of infection [1517]. Besides, reports based on transcriptomic studies suggested early recognition of wound inducing Foc1 by the host. Such early recognition triggered reprogramming of the primary metabolism of the host where ROS (reactive oxygen species), cellular transporters, transcription factors and sugar molecules acted as signal modulators [17, 18]. Apart from these, biochemical analyses and analytical studies on molecular markers and molecular linkages relating to wilt disease were also performed [1921]. But the inferences drawn from transcriptomic studies are rather inadequate without proteomic support, as there are reports of huge numbers of genes with no assigned functions at their protein level. Additionally, the correlation between mRNAs and protein levels are remarkably low and fail to provide indications about post translational modifications or protein-protein interaction that are believed to have significant regulatory effect on defense responses [22]. Thus, in order to predict the actual scenario of pathogen driven molecular signaling within the host, the knowledge of defense responsive proteins are strongly desirable.

The present study involves understanding of chickpea-Fusarium interaction using proteomic techniques like two dimensional electrophoresis (2DE) and mass spectrometry (MALDI-TOF MS and MS/MS) followed by high throughput data base search. These techniques are effectively used nowadays to identify and analyze differentially expressed proteins involved in plant pathogen interaction and also their post translational modifications [23]. Barring a few, most of these proteomic studies are performed on model plants like Arabidopsis or Medicago[24, 25]. Information gathered from these model plants definitely boost up knowledge of plant immunity but biological interpretation of this knowledge in crop models require experimental substantiation. Moreover, some features and processes are likely to be unique for crop plants and hence cannot be approached via model plant in totality [26]. The present study focuses on the legume crop chickpea and its early response to infection by Foc1. This study aims to understand the mechanistic basis of susceptibility and/or resistance offered by two different genotypes (JG62 wilt susceptible, WR315 wilt resistant) respectively. Approximately 100 proteins were significantly identified by MALDI-TOF MS and MS/MS which included differentially regulated as well as unique proteins identified from both resistant and susceptible genotype of chickpea at different time points of 48 h, 72 h and 96 h after infection with Foc1. These identified proteins are categorized and their probable roles in plant defense are illustrated through interaction network based studies.


Plant growth and fungal treatment

Chickpea (Cicer arietinum L.) genotypes JG62 (wilt susceptible) and WR315 (wilt resistant), obtained from ICRISAT (International Crops Research Institute for Semi Arid Tropics), Hyderabad, India were used for experimental analysis. Seeds of both genotypes were grown in a mixture of soil and sand (1:1) under natural green house conditions of 22 to 28°C, 35 to 40% relative humidity and 16 h:8 h photoperiod of day and night respectively [15].

F. oxysporum f. sp. ciceri Race1 (Foc1) was obtained from ICRISAT and further purified according to the protocol of Summerell et al[27]. Spores obtained were harvested and stored at -80°C until further use. Two week old seedlings of both genotypes were inoculated with Foc1 using sick soil method as described by Gupta et al[15]. Plants of both genotypes grown on inoculum free soil served as control samples. Both control and infected plants were kept under same growth conditions. Root samples from control and infected plants at 48, 72 and 96 h post inoculation (hpi) were harvested, instantly frozen in liquid nitrogen and stored at -80°C for further analysis. Proteins were extracted from pooled tissue to run triplicate gels of each time points [28]. The entire experiment of plant growth and fungal treatment was repeated three times to generate three biological replicate.

Protein extraction and quantification

Chickpea root proteins were obtained from one gram of root tissue by following Phenol-SDS buffer extraction method with sonication [29]. One gram of root tissue was pulverized in mortar and pestle with liquid nitrogen and homogenized with 3ml of SDS buffer (30% sucrose, 2% SDS, 0.1M Tris-Cl, 5% β-mercaptoethanol and 1 mM phenyl methyl sulfonyl fluoride (PMSF), pH 8.0). The extract was sonicated (60 amps, 15 secs, 6 times) and further treated with Tris buffered phenol. The phenolic phase obtained by centrifugation at 8000 g for 10 min at 4°C was rinsed with SDS buffer. This final phenolic phase was collected and precipitated overnight with four volumes of 0.1M ammonium acetate in methanol at -20°C. Precipitate was obtained at 10,000 g for 30 min.Washing of protein pellet was performed thrice at 8,000 g for 10 min with cold 0.1 M ammonium acetate and finally washed with cold 80% acetone. The pellet was then dried and resuspended in 100 μl sample buffer (Biorad) for further analysis. Extracted proteins were quantified using Bradford protein assay method using BSA as standard [30].

Two dimensional polyacrylamide gel electrophoresis (2D-PAGE)

Isoelectric focusing (IEF) was carried out on PROTEAN IEF cell (Bio-Rad, USA) using immobilized pH gradient (IPG) strips. Two hundred fifty micrograms of each sample protein dissolved in 185 μl of rehydration sample buffer (8M urea, 2% CHAPS, 50 mM DTT, 0.2% Biolyte ampholytes) was loaded onto 11 cm immobilized pH 3-10 nonlinear (NL) gradient strips (Bio-Rad, USA) and was passively rehydrated overnight at room temperature. IEF was conducted at field strength of 600 V/cm and 50 mA/IPG strip. The strips were focused at 250 V for 20 min, 8000 V for 2 h 30 min with linear voltage amplification and finally to 20,000 Volt hour with rapid amplification. After focusing the strips were reduced and alkylated using 135 mM DTT and 135 mM iodoacetamide respectively, in 4 ml of equilibration buffer (20% v/v glycerol, 0.375M tris- Cl, 6M urea, 2% w/v SDS, pH8.8) for 15 min. Second dimensional electrophoresis was run with strips transferred to 12% SDS polyacrylamide gels (13.8 cm × 13.0 cm × 1 mm) in an AE-6200 slab electrophoresis chamber (Atto Biosciences and Technology, China) at a constant volt (200 V) for 3 h 30 mins in tris-glycine SDS running buffer. The gels were stained with 0.1% (w/v) coomassie brilliant blue R-250 (Sigma) overnight, destained and stored in 5% acetic acid at 4°C. 2D-PAGE gel separation was performed with both technical and biological replications of three.

Image acquisition and analysis

Coomassie stained 2-D gel images were captured with Versa Doc Imaging system (Model 4000, Bio-Rad, USA) and analyzed with PD Quest Advanced 2-D gel analysis software (version 8.0.1, Bio-Rad, USA). For this study in total 72 reproducible gels were generated (three replicates, four time points, two genotypes and three biological replicates). Three technical replicates from three biological replicates at different time points (control, 48 h, 72 h, 96 h) for both genotypes (JG62,WR315) were assembled to create the master gel image (match set). Replicate gels used for making the match set had correlation coefficient value of at least 0.8. Background subtraction between the gels was done using floating ball method. Spots were detected automatically by the spot detection parameter wizard using Gaussian model with advance settings, by choosing faint spot, small spot and large spot cluster. Detected spots were visually checked and manually added when required [31]. Each spot included for analysis were present at least in two of the three replicate gels for a particular time point and also was of high quality. Detected spot volumes were normalized by the spot volume of the entire gel and used as a parameter for quantifying protein abundance. The differential spots which showed statistical significance level of p < 0.05 (Student’s t-test) were selected for analyses. However, the spots selected for downstream MALDI-TOF MS and MS/MS analyses fell under three main categories. Firstly it included the spots showing 1.5 fold changes (above or below) in protein abundance level in infected samples at least in any of the time points as compared to the comparable protein level of both the controls. Second category included spots which were accumulated after infection and present in more than one time point in infected samples but absent in controls. Third category included qualitative spots which are reproducibly present only in one infected variety for a particular time point. Spots which were present only in one replicate were not considered for analysis to minimize the interference of missing value. Experimental molecular mass and pI were calculated using 2D-PAGE gel images of standard molecular mass and pI markers. Data were further analyzed using Statistica v10.0 software (Statsoft Inc) through coefficient of variance calculation (CV), followed by comparison of control and treated values to find out statistical differences by multivariate analysis of variance (MANOVA) and Duncan’s multiple range test (DMRT), at p value 0.05. Protein spots that showed significant difference between treatments through DMRT were further processed for downstream MALDI-TOF MS and MS/MS analyses.

Protein identification using MALDI-TOF MS and MS/MS

Protein spots were manually excised from 2D-PAGE gels, destained and in gel digested according to the protocol mentioned by Shevchenko et al.[32] with minor modifications. In gel digestion of proteins were carried out with porcine trypsin (Promega, USA) and peptides were extracted with 25% acetonitrile and 1% trifluroacetic acid. One microlitre of sample was loaded along with matrix (1 μl, α-cyano-4-hydroxy cinnamic acid, HCCA) (Bruker Daltonics, Germany) in an Anchor Chip MALDI Plate (Bruker Daltonics, Germany).

Mass spectra were generated in an Autoflex II MALDI TOF/TOF (Bruker Daltonics, Germany) mass spectrometer equipped with a pulsed nitrogen laser (λ-337 nm, 50 Hz) in the m/z range from 500 to 3500 Da. The enzyme used was trypsin with one missed cleavage. The spectra obtained were analyzed with Flex Analysis Software (version 2.4, Bruker Daltonics, Germany) for deletion of matrix peaks and tryptic autolysis peaks. Processed spectra were then searched using MS Biotools (version 3.2) program against the taxonomy Viridiplantae (Green plants) in the MSDB 20060831 (3239079 sequences; 1079594700 residues), NCBInr 20140323 (38032689 sequences; 13525028931 residues), SwissProt 2013_12 (541954 sequences; 192668437 residues) databases using MASCOT search engine (version 2.2). The standard parameters used in the search included peptide mass tolerance (±0.5 Da); fragment mass tolerance (±0.8 Da); proteolytic enzyme (trypsin); global modification (caramidomethyl, Cys); variable modification (oxidation, Met); peptide charge state (1+) and maximum missed cleavage of 1, for MALDI-TOF MS minimum S/N = 10 and for MS/MS minimum S/N =3. The significance threshold was set to a maximum of 95% (p</= 0.05). The criteria used to accept protein identification were based on molecular weight search (MOWSE) score, and the percentage of sequence coverage. From each samples most intense m/z values were chosen for further fragmentation (MS/MS). Automatic decoy database search was performed by choosing the decoy checkbox on MASCOT search engine. Decoy search was performed to avoid false identification of peptide by matching it to a random sequence from a decoy database. Only the results with 0% false discovery rate were accepted. Final protein identification was done by a combined search of PMF (Peptide Mass Fingerprint) and MS/MS data in MASCOT search engine.

Protein interaction network generation and analysis

Pathway Studio software (version 7.1) (Ariadne Genomics, USA) and Res Net database (version 3.0) was used to study the biological interactions [33] between the identified proteins of the present study. These differential and unique protein sequences identified by PMF and MS/MS studies were subjected to BLAST analyses at TAIR database (The Arabidopsis Information Resource) and their homologous genes (bearing TAIR gene IDs, Additional file 1) used as inputs for network generation. Ambiguities and components without any interactive neighbors were eliminated from the import list. Interaction network was generated using the neighbor joining method with a degree of correlation as 1 (only the immediate upstream and downstream neighbors having direct relationship to the protein/protein products were considered for analyses). In addition, standard filter parameters and relation types were selected for interaction map generation. Presence of the identified proteins in known biological pathways was analysed using AraCyc and Ariadne Pathway data list. Functional classification of the identified proteins based on gene ontologies (GO) were also studied using Pathway Studio software. In both cases statistical significance (p < 0.05) of the pathway locations and GO classification of the identified protein were calculated.

Quantitative real time pcr (qRT-PCR)

Total RNA was extracted from one gram root tissues of infected and uninfected plants of both genotypes at different time points of 48 h, 72 h and 96 h post infection. RNA was extracted using TRI reagent kit (Himedia, India) as per manufacturers’ instruction. For avoiding any DNA contamination RNA samples were treated with RNase free DNase (Fermentas, USA). cDNA synthesized using Revert Aid first strand cDNA synthesis kit (Fermentas, USA), was further used for qRT-PCR. Specific primers were designed based on the corresponding nucleotide sequence of identified proteins from CTDB (Chickpea Transcriptomic Database), DFCI (Medicago trancatula database), PDB (Protein data Bank) and NCBI database using Gene Runner software (version 3.1) and listed in Additional file 2. qRT-PCR was performed on Biorad i cycler (Bio-Rad I-Q5, USA) using SyBr green super mix. A reaction mix of 20 μl was prepared containing 25 ng cDNA, 0.3 μM of forward and reverse primers. The PCR conditions used were 95°C for 5 mins, followed by 40 cycles at 95°C for 30 sec, 50°C-55°C for 30 sec and 72°C for 30 sec [16]. A melt curve was also generated at the end of each PCR cycle to verify primer specificity. Sample variation was minimized by normalization using actin as internal standard [34]. Mean fold change was calculated using 2 -ΔΔct method [35]. All experiments were repeated three times and standard error was calculated.

Results and discussion

Analysis of chickpea root proteome

Chickpea root proteome was studied with a view to understand the molecular mechanism governing the susceptibility and/or resistance of chickpea plant upon pathogen infection. Previous results based on histopathological and transcriptomic analyses performed by our research group as well as others, suggested the time points of 48 h, 72 h and 96 h to be crucial for delineating the early defense responses of chickpea during Foc1 attack [15, 16, 18]. These previous reports stated 96 h as the onset for xylem vessel colonization in compatible roots, while significant differential transcriptomic alterations were detected at as early as 48 h in both the susceptible and resistant genotypes [15, 16, 18]. An estimated protein yield for all the samples are provided in Additional file 3. Total root proteins were resolved onto 11 cm IPG strip (pH 3-10 NL). Figure 1(A and B) shows representative 2D experimental gel profiles corresponding to control and infected samples at different time points for both the genotypes, JG62 and WR315 respectively. The experimental design is shown in Additional file 4. Three independent experiments were performed to ensure that the changes in protein abundance at each time point were reproducible and significant. Two dimensional gel analyses indicated differential protein profiles for JG62 and WR315 plants upon Foc1 infection. Further PD Quest software analysis detected a total of 274 spots in the master gel (Figure 1C).The number of total spots detected and the differential spots (quantitative and qualitative) obtained post inoculation with Foc1 for each sample is provided in Additional file 3. To assess the reproducibility of the corresponding protein quantification, the CV was calculated for all protein spots, at all time points examined. The CV of protein spots for each sample type and time points was within 21% which is in accordance with other plant stress related studies [36] indicating stability and reproducibility of the present data. Among the total 206 differential spots obtained 163 spots which fell under the previously described three categories were processed for downstream MALDI-TOF MS and MS/MS analysis. MS/MS analyses was performed with 137 spots of which 100 spots that showed significant scores were taken into consideration for further functional clustering. Differential spots obtained due to differences in genotypes, depicting the natural variation between the susceptible and resistant genotypes (i.e differentially abundant between control samples of JG62 and WR315) were excluded from further downstream analyses in the present study (data not shown). Relevance of such differences between both genotypes that could also add significantly to the understanding of chickpea-Foc1 interaction shall be dealt separately in future studies. Selected protein spots were found to be interspersed at and around the median region of IPG strip suggesting the critical pH range for resolving the differential proteins to be around pH 4-7 (Figure 1C). Finally MS/MS analyses using mascot search engine in the available databases (NCBI, MSDB, Swissprot) led to the successful identification of 100 spots (Figure 1C). The details of these proteins and their peptides identified by MS/MS is provided as a table in Additional file 5. Among these 100 spots, 65 spots showed significant (1.5 fold change) quantitative changes in infected genotypes (JG62 and WR315) as compared to comparable protein level in control and 35 spots showed qualitative changes. Out of these 35 spots, 28 were accumulated after infection in more than one time point of either/or both infected genotypes, absent in controls and 7 spots were unique for any one time point and genotypes. MANOVA followed by DMRT indicated the statistical significance of the data provided in Additional file 6. Means that do not share any common alphabet differ significantly by DMRT at 5% level.

Figure 1
figure 1

Root proteome expression profile of control and Foc1 ( F. oxysporum f.sp. ciceri Race1) infected chickpea genotypes JG62 (A), WR315 (B). Root proteins (250 μg) of control (Jc) and infected JG62 extracted at different time intervals of 48 h, 72 h, 96 h (J48, J72, J96) post infection (A). Root proteins (250 μg) of control (Wc) and infected WR315 extracted at different time intervals 48 h, 72 h, 96 h (W48, W72, W96) post infection (B). Total proteins separated in first dimension (IEF,11 cm IPG strip, 3-10 NL) followed by second dimension in 12% SDS-PAGE, stained with coomassie brilliant blue R-250. Master gel generated by PD Quest (C). Identified spots are encircled and the number corresponds to the spot number mentioned in Additional file 5.

Identification and classification of differential and unique proteins in chickpea during Foc1 infection

The identified proteins were classified into nine functional categories based on their putative biological functions and proteins with unassigned functions were categorized as unclassified group. Metabolism related protein (36%) constitute the most abundant group followed by proteins related to scavenging of reactive oxygen species (ROS) (16%), protein synthesis and degradation related proteins (11%), defense related proteins (7%), signaling proteins (7%), storage proteins (6%), transport proteins (4%), developmental proteins (3%), structural proteins (1%) (Figure 2A). The unclassified group accounts for 9% of total identified proteins. Metabolism related proteins were further classified into glycolysis related proteins (31%), proteins of TCA cycle (17%), ATP synthesis and degradation regulating proteins (14%), proteins related to amino acid metabolism (19%), secondary metabolism (8%) and sugar metabolism (5%). Moreover, 3% proteins were found to be related to electron transport and another 3% were related to cell wall metabolism and transport (Figure 2B). Many defense related proteins identified had scores below 70. Complete draft genome sequence of chickpea has been recently reported, but functional annotations of genes and gene products are still at initial stages. Chickpea is a legume and its closest completely sequenced neighbor legumes are model plants Medicago and Lotus. However, the functional annotation of these neighbor model legumes are also underway and constantly being updated. Besides, chickpea being a crop legume is expected to have some distinct differences with these model legumes. Such differences are likely to be reflected in the protein identification scores of chickpea when subjected to homology matches with these model legumes. Hence, all the important defense related proteins obtained from the present study were discussed even though their scores were in the range of 40-60. Previous studies conducted with chickpea also reported similar identification scores for protein identification [36, 37]. In most of the cases each protein spots were identified as a single, unique protein but in some cases the identified protein spots contained more than a single protein; in such cases, the first hit with maximum score was considered for their protein IDs [38]. In addition to this, multiple spots were also found which were identified as the same protein. The appearance of such proteins probably suggests them being chemically and/or molecularly different products of a single gene and referred to as protein species [39] (Additional file 5). They basically fall under three main categories (i) with same molecular mass and different pI; for example, Kunitz proteinase inhibitor (sp 2, sp 13), Annexin (sp107, sp 499), Glyceraldehyde 3 phosphate dehydrogenase EC (GAPC) (sp505, 509), (ii) with different molecular mass but same pI; for example Cysteine proteinase (sp 99, sp 45),( iii) or with different molecular mass and pI; for example, Superoxide dismutase EC (sp 67, sp103, sp401), Triose phosphate isomerase EC (sp 109, sp 40, sp 85). The differences in Mr and pI values, suggest that these changes in the proteome are probably due to the post-transcriptional modification. They may belong to different members of the same functional family, indicated by small shift in the pI or are degraded protein products as suggested by significant differences between theoretical and observed Mr values. The slight differences in pI and Mr values probably reflect post translational modifications (like phosphorylation, acetylation, glcosylation, methylation) occurring in vivo or may be the result of modifications such as deamidation of the proteins during sample preparation and processing [25]. It is known that the same protein may have different functions in different subcellular compartments. In the present study superoxide dismutase (sp 103,sp 401), triose phosphate isomerase (sp 109,sp 85) and GAPC (sp 505, sp 509) were identified as protein variants present in different cellular compartments like mitochondria, chloroplast or cytosol. Hence their multiple forms may be attributed to their multiple cellular locations [39]. In most of the stress related studies GAPC showed post translational modification like phosphorylation and was found to be present as multiple protein species. But whether the same observation in the present study indicates same modifications needs validation [40].

Figure 2
figure 2

Distribution of functional classification of Identified proteins. Functional classification and relative distribution of proteins altogether identified in JG62 and WR315 chickpea genotypes after infection (A). Classification and categorization of metabolism related proteins (B). The proteins belonging to different categories and their expression level at different time intervals are mentioned in Additional file 5 in detail

Proteins related to direct defense responses against Foc1

Defense related proteins contribute to about 7% of total identified proteins. They include PR1 (pathogenesis related protein 1), BGL, EC (glucan endo 1-3 beta glucosidase), TLP (thaumatin like protein) and TPI (trypsin protease inhibitor) (Figure 3, Additional files 5 and 7). Pathway analysis showed the association of these proteins with defense and hypersensitive response related pathways (Additional file 8). Gene ontology (GO) based classification showed their relation with biological processes, molecular function and their cellular location (Additional files 9 and 10). Schematic network showed the interaction of these components with other Foc1 inducible proteins (Figure 4). PR1(sp 145) protein known to be directly involved in plant defense against pathogen attack was found to be accumulated at 48 h and 72 h post infection in resistant plants while in case of susceptible plants protein level was not detectable after infection (Figure 3, Additional files 5 and 7). PR1 expression known to be regulated by salicylic acid (SA) is positively regulated by NPR1 (Non expressor of PR genes1) during defense [41]. Besides, ACD (accelerated cell death), known to accelerate cell death in Arabidopsis is also a positive regulator of PR1 [42]. MAP kinase (Mitogen activated protein kinase), EDS4 (Enhanced disease susceptibility 4), PAD2 (Phytoalexin deficient 2) linked to fungal defense response also regulate PR1 expression [43, 44]. On the other hand, studies conducted on Arabidopsis thaliana reported EDR2 (Enhanced disease resistance 2), NPR3 and NPR4 to be negative regulators of PR1 [45, 46]. PR1 expression is also reported to be altered by phospholipase C and fatty acids [47, 48]. In the present study the increase of PR1 protein in resistant plants suggests its direct role in Foc1 induced defense, although the role of SA in modulating resistance in the present case study is still speculative. BGL also known as PR2, are enzymes which mainly act by hydrolyzing 1-3 β D glucosidic linkage of fungal cell wall and hence known to provide resistance in plants. BGL (sp 239) was found to be up accumulated in response to fungal attack in both genotypes. However, the susceptible plants showed highest accumulation at 72h (Figure 3, Additional files 5 and 7) that decreased later. Both BGL and PR1 are known to have SA dependent expressional regulation [49]. Both PR1 and BGL are reported to be upregulated in over expression lines containing EIL (ethylene-insensitive3-like) transcription factor in Vigna mungo indicating a positive role of ethylene in regulating defense response [50]. TLPs are pathogenesis related proteins having antifungal activity. TLP, also known as PR5 (sp 83,129) was found to be significantly increased in response to Foc1 in both genotypes (Figure 3, Additional files 5 and 7). However, in resistant plants it showed uniform accumulation while in susceptible plants (sp 129) it was found to be absent at later time points (72 h and 96 h). TLP was found to be up accumulated in Medicago trancatula during Orobanche crenata infection indicating that it may eventually take part in defense mechanism against parasitic infection [25]. TPI are known to participate in the wound induced defense response of plants against herbivores and pathogens. TPI (sp 2, 13, 81) were found to be uniformly enhanced in response to Foc1 induction in resistant plants while susceptible ones showed protein level undulations (except for sp 81, which showed uniform protein accumulation). (Figure 3, Additional files 5 and 7). TPI is positively regulated by JA signaling [51]. WRKY transcription factors coordinating herbivory are also known to regulate TPI expression [52]. The induction of TPI probably indicates the involvement of SA/JA mediated hormonal crosstalk which needs further experimentation. Role of PR proteins (PR1, PR2 and PR5) in modulating defense network were also elaborated by transcriptomic as well as proteomic studies involving wheat (Triticum aestivum L.) and stripe rust fungus Puccinia striiformis f.sp. tritici Eriks. (Pst) [53, 54].

Figure 3
figure 3

Heat map representation of differentially expressed proteins of JG62 and WR315 chickpea genotypes on infection with Foc1. Heat map was generated with the fold change values considering infected/control ratios. Each column represents a particular time point of infection and each row represents corresponding proteins with their identities. Up regulation or down regulation is indicated by the above scale which shows pale to saturated colors of green and red respectively. Yellow color represents mid-value and white represents no expression.

Figure 4
figure 4

Schematic representation showing the location and interaction between the different Foc1 induced proteins in chickpea roots. Representation shows the intra and inter relationship between the Foc1 induced proteins and their regulatory biological processes. (Complete names of abbreviated proteins are provided in Additional file 1).

Role of ROS scavengers/regulators

Sixteen percent of total proteins were classified as ROS scavengers/regulators. Superoxide dismutases, EC (SOD), Peroxiredoxin proteins, Ascorbate peroxidase, EC (APX), Ferric reductase EC, Glutathione S transferase, EC (GST), Peroxidase, Thioredoxin (NTRA, NTRB), Monodehydroascorbate reductase, EC (MDHAR, MDAR), Quinone oxidoreductase, EC etc (Figure 3, Additional files 5 and 7) are the proteins included in this class. Pathway analysis showed association of some proteins (SOD, APX, NTRA and NTRB, MDHAR and MDAR etc) with ROS regulatory pathways (Figure 4, Additional file 8). GO classification illustrated the roles of these proteins according to their biological processes, molecular functions and cellular components (Additional files 9 and 10). Figure 4 showed their cellular location and their interaction with other Foc1 induced proteins. SODs (sp 67, 103, 401) showed oscillations in protein accumulation post infection in both genotypes (Figure 3, Additional files 5 and 7). SODs are known to provide the first line of defense to infected hosts by scavenging the pathogen triggered ROS [55]. SODs are also reported to induce ROS mediated PR1 expression in Nicotiana[56]. APX (sp 134, 87, 323) (Figure 3, Additional files 5 and 7) is an important enzyme participating in anti oxidation metabolism in plants [57]. Besides, they are also reported to be upregulated during heat stress in Arabidopsis[58]. Differential induction of SOD and APX in the present study indicated the role of Foc1 induced ROS in triggering defense responses in chickpea. These observations support previous reports based on transcriptomic studies [16, 18]. Ferric reductase plays an important role in maintaining iron homeostasis, disruption of which may lead to generation of toxic free radicals. Ferric reductase (sp 232) also known to be an antioxidant for peroxides, showed enhanced protein level in susceptible plants compared to resistant ones (Figure 3, Additional file 5). GST (sp 324,211,317) showed marginal changes in protein accumulation in resistant plants while susceptible plants showed relatively sharp increments and decrements in protein level post Foc1 induction (Additional files 5 and 7). GSTs are reported to reduce oxidative stress inductive organic hydroperoxides in Nicotiana benthamiana following Colletotrichum destructivum infection [59]. In the present study, steady state protein level of GST in resistant plants may indicate lesser accumulation of oxidative stress components as compared to susceptible plants. NTR (sp 78) showed increment in resistant plants following Foc1 infection while susceptible plants showed sharp decline after 72 h of infection (Figure 3, Additional file 5). Such up accumulation of NTR only in resistant plants indicated their efficient role in regulating oxidative stress tolerance [60]. Previous proteomic studies conducted on wheat showed enhanced accumulation of GST and NTR during incompatible interaction with Puccinia striiformis f.sp. tritici Eriks. (Pst). Besides, level of peroxiredoxin was also found to be induced [54]. In addition transcriptomic studies showed the enhancement of peroxidase transcripts in wheat following Puccinia striiformis f.sp. tritici Eriks. (Pst) infection [53]. MDHAR (sp 199) showed similar protein accumulation levels in both plants post infection (Figure 3, Additional files 5 and 7). Such increment indicated role of MDHAR in JA mediated antioxidation metabolism in the present case study that was found to be similar to previous results reported on Arabidopsis thaliana[61]. Besides, increment of MDHAR also linked to increased lipid peroxidation which is marked as a feature during pathogen mediated membrane injury [62]. Quinone oxidoreductase, known to act as detoxifier of ROS induced oxidative stress along with GST was found to be up accumulated at later time points of infection in susceptible plants compared to resistant ones (Figure 3, Additional file 5).

Role of signaling proteins

Signaling proteins constitute about 7% of total identified proteins. Guanine nucleotide binding protein (AGB), Annexins (ANNATs), ABA responsive protein (RAB), Ran binding protein (RANBP), Auxin induced protein and Zinc binding dehydrogenase are classified under this category (Additional file 5). Pathway analysis based on Arabidopsis homologues showed only the association of AGB with signaling pathway (Figure 4). While GO classified all the proteins in this category (AGB, ANNATs, RAB18 and RANBP) according to their relation with biological processes, molecular function and cellular components (Additional files 9 and 10). Network map showed their interaction with other Foc1 induced proteins (Additional file 8). AGB (sp 233) coupled with other G proteins and GPCRs are known to modulate defense responses in Arabidopsis[63]. Besides, AGB are also known to modulate ABA driven K+ and anion channels thus regulating stomatal movement [64]. In the present study, similar protein accumulation pattern of AGB (Figure 3, Additional file 5) in both plants indicate a common regulation of AGB that is probably directed towards stomatal movement, a significant phenomenon observed during vascular wilt. Annexins (sp 107, 499) (Figure 3, Additional files 5 and 7) are reported to regulate pH mediated cellular responses that are directly influenced by ABA and calcium conductance during stress in Arabidopsis and Zea mays[65, 66]. The up accumulation of annexins in both plants probably directs the role of Foc1 in triggering pH alterations as well as ABA driven calcium oscillations during infection that needs to be investigated. RAB (sp 71) was found to be up accumulated only in resistant plants post infection. RAB was reported to be induced during ABA perception that activated calcium influx in Arabidopsis thaliana suspension culture cells [67]. Such induction was further known to be mediated by phospholipase D activation [68]. In the present study induction of RAB only in resistant plants directs towards role of ABA and calcium signaling in modulating defense in chickpea during Foc1 infection. RANBP (Sp16) known to regulate nucleocytoplasmic transport under the control of hormones and light, was found to be uniquely expressed at 72 h post infection in susceptible plants [69]. The relevance of such selective induction in the present study requires further investigation.

Role of metabolism related proteins

Majority of the proteins identified (36%) fell under metabolism related proteins (Figure 2A). This category was further re-categorized into several sub classes (Figure 2B, Figure 3, Additional files 5 and 7). Such large assemblage of metabolism related proteins indicates that pathogens usually target the host metabolism for self survival and reproduction, while on the other hand host puts forth complete effort in shielding their primary metabolism from the devastations of pathogen attack [16]. Pathway analysis showed the association of some of these proteins with metabolic pathways (Figure 4). GO classification grouped them according to their biological processes, molecular functions and cellular components (Additional files 9 and 10). Interaction map further showed the location and interaction of some of these proteins with their neighbors as well as within themselves (Additional files 8). Glycolytic enzymes triose phosphate isomerase, EC (TIM) (sp 109, 85, 40) and glyceraldehyde dehydrogenase phosphate, EC (GAPC) (sp 212,505,509) were found to show similar pattern of protein level undulations in both compatible and incompatible interaction suggesting the common role of glycolytic ATP on pathogen triggered immune response of host [70]. However, enolase EC (LOS) (sp 198) showed sharp decline at later time points of infection in susceptible plants while resistant plants showed steady state protein level (sp 198) or sharp induction (sp 182, 351) at different time points of infection. Enzymes of TCA cycle such as isocitrate dehydrogenase, EC (ICDH) (sp 165), malate dehydrogenase, EC (sp 156) and fumarase, EC (FUM) (sp 271,272) showed elevated or stable protein accumulation in resistant plants as compared to susceptible plants suggesting a constant energy supply, which is required for different processes like photosynthesis, respiration and photorespiration during stress [7173]. ATP synthase (sp 36, sp 196) and ATPases (VHA) (sp 480,124,168) however showed similar protein level patterns in both plants after infection. This may indicate the need for maintaining energy and solute homeostasis necessary for protein sorting and cell wall repair that probably aid to cell protection during pathogen progression [74]. Similar interaction studies involving wheat and stripe rust fungus reported the increment of ATP synthase both at transcriptomic as well as proteomic levels [53, 54]. Cytochrome c oxidase, EC (COX) (sp 11) which is known for translocation of protons to drive aerobic respiration as well as to regulate stress mediated signals [75] showed elevated protein levels in resistant plant at later time period as compared to susceptible plants. These findings suggest that even though, energy requirement is necessary for both the genotypes during stress, but proper channelization of energy needed for running basic metabolic activities controls resistance, which perhaps is efficiently maintained by the resistant plants. Cysteine protease (RD) (sp 45) showed up accumulation of protein only after Foc1 infection in both plants. RD is known to be important players in plant immunity, especially in regulating resistance response against necrotrophic pathogen [76]. In the present study the selective up accumulation of RD after Foc1 infection predicts the role of RD in regulating biotrophic interaction also. However, such assumption requires further experimental support. Phosphoserine amino transferase, EC (PSAT) (sp 264) was found to be absent at 96 h in susceptible plants while resistant plants maintained a moderate protein accumulation level even after infection suggesting the need of serine biosynthesis which is known to be associated with photorespiration [77]. S-adenosyl methionine synthetase, EC (SAM) (sp 122) a direct product of methionine catabolism acts as substrate for several transmethylation reactions including those that occur during lignin biosynthesis [78]. SAM was found to be absent at 96 h post inoculation in susceptible plants while resistant plants regained the protein accumulation at 96 h suggesting the role of transmethylation and lignin biosynthesis in somehow regulating repair mechanisms caused by pathogen invasion. Proteins related to secondary metabolism (sp 207, 173, 164) showed differential abundance level post infection in both plants. They are known to regulate defense response during biotic stress [79]. Methylesterase (sp 290) were found to be selectively enhanced at 72 h post infection in both the plants. Methylesterases are known to be directly or indirectly associated with defense reactions by regulating the degree of methyl esterification of pectin that is known as essential cell wall components [80]. Selective accumulation of methyl esterase after infection in both plants suggests a possible cell wall repair mechanism to be operational, which however may be more efficient in resistant plants as indicated by its elevated level.

Role of proteins involved in its folding, synthesis and degradation

This group of proteins accounts for about 11% of total identified proteins. Pathway analysis showed the association of adenylate kinase EC (ADK, AMK2) with protein synthesis and purine biosynthetic pathways (Figure 4). While all other proteins related to protein synthesis, folding and degradations showed enlistment under categories of GO (biological function, molecular function and cellular component) (Additional files 9 and 10). Network analyses also showed the intra and inter relationship of these proteins with other Foc1 induced proteins (Additional file 8). 26S proteasome subunits are known to contribute to both basal defense as well as R gene mediated defense in Arabidopsis. Activation of these proteins are known to regulate innate immunity both positively and negatively as appropriate protein degradation are necessary for mounting defense [81]. Besides, studies on Nicotiana reported the induction of 20S proteasome subunits that was found to be linked to HR and SAR [82]. In the present study the differential accumulation of 26S proteasome subunits EC (sp 57, 10, 240) in both the plants post infection suggests the role of protein degradation in regulating defense (Figure 3, Additional files 5 and 7). However, whether such regulation is directed towards positive and/or negative influences needs to be investigated in detail separately for both compatible and incompatible interaction. Adenylate kinase EC (ADK, AMK2) (sp 234) was found to be up accumulated in 48 h in resistant plants while susceptible plants maintained an overall low protein level. ADK, known to be involved in salvage pathways of adenine and adenosine also convert cytokinin and ribosides to corresponding nucleotides. Such cytokinin conversion regulates the hormonal level of plants [83]. Absence of ADKs is known to cause chloroplastic deformity in Arabidopsis[84]. In the present case study overall down accumulation of ADKs (Figure 3, Additional files 5 and 7) probably indicates pathogen mediated chloroplastic damage and hormonal alteration. Eukaryotic translation initiation factor (elF5lpha) (sp 106) protein was uniquely accumulated at 48 h in resistant plants while susceptible plants showed no accumulation. Studies conducted on Arabidopsis showed the involvement of elF5alpha in controlling resistance by preventing pathogen growth and development of Pseudomonas syringae[85]. Besides, elF5alpha was also up accumulated during infection with stripe rust fungus in resistant wheat plants [54]. However, whether the accumulation of elF5alpha protein at a specific time point post infection in resistant plants has a similar role in restricting the pathogen progression needs to be experimented. Protein disulphide isomerase (PDIL) (sp 102) was found to be up accumulated at 48 h in resistant plants which gradually declined at later time points however maintaining a moderate level compared to control plants even after 96 h of infection. In susceptible plants the accumulation level of PDIL were greater compared to control samples only at 72 h post infection and absent in other time points (Figure 3, Additional files 5 and 7). PDIL acts as chaperones of cysteine proteases, thus regulating their trafficking from endoplasmic reticulum to vacuole prior to PCD [86]. Besides, PDIL are also known to be reduced by thioredoxin reductases and actin and removing aberrant disulphides formed by oxidative stress [87]. Level of PDIL was found to be elevated in wheat following inoculation with stripe rust fungus [54]. The abundance of PDIL in resistant plants in the present case study suggests the operation of antioxidant defense machinery during incompatible interaction in chickpea against Foc1 attack.

Role of developmental, structural, channel and storage proteins in Foc1 induced defense

Developmental, structural, channel and storage proteins contributes to about 14% in total.

Functional classification identified developmental proteins such as ripening related protein (RLP) (sp 137) and germin (sp 556, 478) to be differentially expressed in both the plants after infection (Figure 3, Additional files 5 and 7). RLP contains the conserved Bet v fold domain also present in major latex proteins (MLPs) and PR10 group of allergen proteins. These proteins are associated with fruit and flower development as well as defense. However their role in defense is not well characterized [88]. In the present study the protein abundance of RLP at 48 h in resistant plants suggests this protein somehow modulate initial defense response which requires further characterization. Germins, known to have roles in plant development and defense, are associated with extra cellular manganese-SOD activity [89]. The up accumulation of sp 556 protein at later time points (96 h) in resistant plants and protein level undulations of sp 478 in both plants post infection suggests a differential operation of antioxidant defense mechanism in controlling pathogen invasion in both plants. However, transcriptomic based studies reported the increment of germin like transcripts in response to stripe rust fungus specifically in resistant genotypes of wheat [53]. Structural protein profilins (PRFs) are actin monomer binding proteins that regulate the assembly-disassembly of uncapped-capped actin molecules in forming cytoskeletal filaments [90]. Profilin (sp 3) protein was found to be uniquely accumulated at 48h in resistant plants post infection. Such selective accumulation probably indicated the need of cytoskeletal assembly to strengthen the cell and prevent further fungal ingress. However, such assumption needs further experimental support. Channel proteins porin (sp 267) and plasma membrane intrinsic protein (PMIP/PIP) (sp 19) belong to the aquaporin family of proteins that are known to regulate hydraulic conductance during cold and oxidative stress [91]. The present study showed differential protein accumulation profiles in both genotypes after Foc1 induction, which suggested that probably the channel proteins regulated water transport differently during incompatible and compatible interaction. The enhanced level of these proteins at 96 h post infection in resistant plants compared to susceptible ones suggested proper water conductance in resistant plants when susceptible plants succumbed to wilting symptoms. Plant albumins are known to serve as storage proteins as well as defense responsive proteins possessing insecticidal and antimicrobial properties that are induced in response to stress [92]. The present study showed up accumulation of albumin (sp 244, 485) in resistant plants compared to susceptible ones indicating the role of storage proteins in controlling defense against fungal attack.

Unclassified proteins

This group mainly includes proteins with unknown functions (Figure 3, Additional file 5). They contribute to about 9% of total proteins identified (Figure 2). Recent availability of chickpea whole genome sequences and updating of functional annotations is believed to provide proper naming and functional designations to these unclassified proteins [93].

Probable roles of identified proteins in imparting resistance against Foc1

To understand the mechanism of resistance in plants it is important to know what are the different proteins involved and how they come into play during the pathogen attack. Pathogenesis related proteins are defense related proteins which are induced on pathogen attack and have a direct role in plant defense, but how these proteins operate or accumulate in compatible and incompatible interaction actually decides the sustainability of resistance. In the present study three important PR proteins were identified, PR1 (pathogenesis related protein 1), PR5b (Thaumatin like protein), PR2 (β-1, 3-glucanases) and their accumulation at different time points post infection were studied. All these PR proteins were found to show a stable level of accumulation in resistant plants after infection where as in susceptible interaction although the proteins appear in early time points but at later time point they either decrease or disappear (Additional file 5). More specifically PR1 which has antifungal activity showed high level accumulation in resistant chickpea plants post infection (Figure 3, Additional file 7). Both PR1 and PR2 are also known to be associated with salicylic acid and ethylene signaling indicating their probable roles in modulating defense [49, 50]. In addition, uniform accumulation of TPI (trypsin protease inhibitor) in resistant plants pointed towards the role of JA in regulating defense [51]. PR5b (Thaumatin like proteins) are known to be induced exclusively in response to wounding or pathogen infection. This protein exhibit a balanced accumulation in resistant plant which indicated its role in disease resistance (Additional file 5). This protein is known to inhibit hyphal growth and reduce spore germination probably by altering membrane permealization or by interacting with pathogen receptors [94]. PR2, involved in cleavage of the β-1, 3-glucosidic bonds of β-1, 3-glucan of fungal call wall, was found at elevated level post infection in susceptible plants which decreased at later time point (Additional file 5). This probably indicated an initial struggle between susceptible plants and pathogen, which was followed by pathogen overpowering the host. However in resistant plants a high and stable PR2 protein accumulation (Additional file 5) indicated that oligosaccharides from fungal cell wall probably acted as an activator, for other PR proteins or antifungal compounds, such as phytoalexins [95]. Besides, this study also mentions high level accumulation of a developmental protein namely ripening related protein (RLP) in incompatible interaction after infection (Additional file 5). Ripening related proteins are found to share homology with some defense responsive proteins in plants. Defense responsive proteins are often found to express during fruit ripening, suggesting that both these processes possibly share a common regulator [96]. Ethylene acts as a key regulator in fruit ripening and in response to stresses caused by pathogens and wounding [97]. Moreover accumulation of S-adenosyl methionine synthetase at later time period in resistant plant also indicates a role of ethylene in plant defense (Additional file 5). Hence these findings indicate that ripening related protein may have a dual role during plant defense and fruit ripening. In addition, this study also provides evidence of hormonal cross talks in host chickpea involving SA, JA, ABA and ethylene during Foc1 invasion.

Salicylic acid is also associated with ROS generation. In the present study the proteins related to ROS scavenging showed abrupt increase and downfall in their accumulation level in susceptible plants after infection whereas in resistant plants the accumulation of protein was found to be in a synchronized way and maintained stability (Additional file 5). It may be assumed that the ROS machinery gets activated in both resistant and susceptible plants after infection but acts differently. The sudden generation of ROS and lack of proper ROS scavenging machinery leads to oxidative stress in case of compatible interaction while in case of incompatible interaction they are efficiently detoxified by scavenging machinery. Hence balanced ROS generation in resistant plants act as signaling molecules and communicate downstream defense signals. Previous studies based on transcriptomic profiling indicate the significance of several ROS regulators to act as the initial trigger communicating downstream defense signals [16]. Interestingly the present case study identified similar set of ROS regulating proteins that not only provided correlation between the transcriptomic and proteomic studies, but also highlighted the conservation of ROS components in regulating host defense during Foc1 infection.

The high accumulation of signaling protein like, ABA responsive protein (RAB) only in resistant genotype in present study predicts the involvement of ABA-mediated signaling in plant defense (Additional file 5). ABA responsive protein was reported to be involved in PR-protein induction and disease resistance in other related studies [98]. Besides, the similar accumulation pattern of AGB (Guanine nucleotide binding protein) and annexinD1 in both genotypes after infection further highlighted the role of ABA and calcium in regulating defense signals.

The accumulation of proteins related to energy metabolism, ATP synthesis and degradation, amino acid metabolism, secondary metabolism etc, in both genotypes; although of different levels suggest that in both cases the pathogen targets the primary metabolism of the host (Additional file 5). Resistant plants probably safeguard their essential metabolic elements from the fungal catastrophe while the susceptible plants fail to do so and submit to pathogenic endeavors.

Initially the events of PTI and ETI were thought to be distinct, but recent studies revealed that the components of PTI and ETI overlap [5]. In present study the initial accumulation of ROS scavengers and regulators direct towards possible responses related to PTI. However categorization of other proteins under the categories of PTI and/or ETI could prove to be erroneous without further experimentation. Even then, all these findings as a whole indicate that plant defenses are controlled by complex signaling pathways which are interconnected to each other.

Correlation between protein and mRNA levels

To correlate the protein levels with mRNA levels eleven representative genes corresponding to MS/MS identified proteins were selected and their transcript accumulation versus protein abundance analyzed (Figure 5). The genes of corresponding proteins selected for transcript accumulation were pathogenesis related protein 1(PR1) (sp 145), thaumatin like protein (TLP) (sp 83), glucan-endo-1,3-beta-glucosidase EC (BGL) (sp239), elongation factor 1 (EF1) (sp 275), protein disulfide isomerase (PDI) (sp 102), guanine nucleotide binding protein (GNBP) (sp 233), triose phosphate isomerase EC (TIM) (sp 109), fructokinase like protein (FLP) (sp54), isocitrate dehydrogenase EC (ICDH) (sp 165), plasma membrane intrinsic protein (PMIP) (sp19) and superoxide dismutase EC (SOD) (sp 103). In general the abundance of mRNA differed from that of protein levels suggesting that the fold increment and/or decrement in mRNA accumulation do not correlate with the protein fold changes. Except for PR1, TLP and ICDH, the other eight proteins and their corresponding transcripts showed a similar qualitative trend in their accumulation patterns. However, the profiles did not provide quantitative similarity. PR, TLP and ICDH showed dissimilar patterns suggesting that mRNA and protein levels often exhibit different profiles. Transcript to protein production involves several regulatory factors which are spatially and temporally regulated due to which there are seldom profile matches between mRNA and protein levels. Moreover this disparity between mRNA and protein level might be due to posttranscriptional or posttranslational modifications, complexities of protein expression or presence of multigene families [39]. Similar results were reported in proteomic analysis of strawberry during Colletotrichum fragariae infection [99].

Figure 5
figure 5

Comparison of mRNA and protein expression levels of eleven representative genes. Quantitative real time PCR was performed using gene specific primers (Additional file 2). The log10 transformed fold change values (infected/control) of protein spot intensities and mRNA expression level were plotted at different time intervals (48 h, 72 h and 96 h) after infection with Foc1 for both chickpea genotypes. JG62 represented by grey color bars and WR315 represented by black color bars. The proteins selected are SOD (superoxide dismutase), PMIP (Pasma membrane intrinsic protein), ICDH (Isocitrate dehydrogenase), FLP (Fructokinase like protein), TIM (Triose phosphate isomerase), GNBP (Guanine nucleotide binding protein), PDI (Protein disulfide isomerase), EF1 (Elongation factor1), BGL2 (Glucan-endo-1, 3-beta-glucosidase), TLP (Thaumatin like protein), PR1 (Pathogenesis related protein1).


The present study was an attempt to investigate the differential root proteome and identify defense related proteins in chickpea during Foc1 infection. Previous report based on proteome studies involving chickpea Foc5 and root knot nematode Meloidogyne artiellia highlighted the presence of several defense responsive proteins [100]. But the difference in pathogenic race is expected to yield some case specific results and hence needs to be studied as an individual case study. The findings of this study suggests that albeit some common proteins are accumulated in response to Foc1 infection in both compatible and/or incompatible chickpea genotypes, but their differential temporal accumulation and regulation probably governs the net outcome of the interaction. The present study highlights the role of several important proteins like PR proteins (PR1, BGL2, TLP), Trypsin protease inhibitors (TPI), ABA responsive protein (RAB18), cysteine proteases (RD19, RD21), methylesterases, 26S proteasome subunits, protein disulphide isomerase (PDIL), ripening related protein (RLP), profilins (PFRs) and albumins and their varied accumulation in susceptible and resistant plants. The functional characterization of these proteins could not only yield important new findings in reevaluating the resistance mechanism of chickpea during Foc1 infection but also help in directing crop improvement programs by using breeding and genetic engineering techniques. Therefore further experiments are necessitated to strengthen the knowledge and understanding through detailed investigations.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files. The protein and peptide data sets supporting the results are presented in Additional file 5.



Fusarium oxysporum f. sp ciceri Race 1


Two-dimensional polyacrylamide gel electrophoresis


Matrix-assisted laser desorption/ionization time of flight tandem mass spectrometry


Gene ontology


Real-time quantitative reverse transcriptional PCR


Pathogenesis related


Pattern triggered immunity


Effector triggered susceptibility


Effector triggered immunity


Food and Agriculture Organization


Isoelectric focusing


Multivariate analysis of variance


Duncan’s multiple range test


Salicylic acid


Jasmonic acid


abscisic acid


Non expressor of PR genes1


Accelerated cell death

MAP kinase:

Mitogen activated protein kinase


Enhanced disease susceptibility 4


Phytoalexin deficient 2


Enhanced disease resistance 2




Reactive oxygen species


G-protein coupled receptors


Hypersensitive response


Systemic acquired resistance.


  1. Richard J, Panstruga R: Tête à tête inside a plant cell:establishing compatibility between plants and biotrophic fungi and oomycetes. New Phytol. 2006, 171: 699-718. 10.1111/j.1469-8137.2006.01829.x.

    Article  Google Scholar 

  2. Tsuda K, Katagiri F: Comparing signaling mechanisms engaged in pattern triggered & effector triggered immunity. Curr Opin Plant Biol. 2010, 13: 459-465. 10.1016/j.pbi.2010.04.006.

    Article  CAS  PubMed  Google Scholar 

  3. Buchanan BB, Gruissem W, Jones RL: Biochemistry and Molecular Biology of Plants. 2007, New Delhi: IK International

    Google Scholar 

  4. Jones JDG, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.

    Article  CAS  PubMed  Google Scholar 

  5. Thomma BPHJ, Nu¨rnberger T, Joosten MHAJ: Of PAMPs and effectors: The blurred PTI-ETI dichotomy. Plant Cell. 2011, 23: 4-15. 10.1105/tpc.110.082602.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Lavin M, Herendeen PS, Wojciechowski MF: Evolutionary rates analysis of Leguminosae implicates a rapid diversification of lineages during the tertiary. Syst Biol. 2005, 54: 575-594. 10.1080/10635150590947131.

    Article  PubMed  Google Scholar 

  7. Ferguson BJ, Indrasumunar A, Hayashi S, Lin MH, Reid DE, Gresshoff PM: Molecular analysis of legume nodule development and autoregulation. J. Integrative Plant Biol. 2010, 52: 61-76. 10.1111/j.1744-7909.2010.00899.x.

    Article  CAS  Google Scholar 

  8. Haware MP, Nene YL: Races of Fusarium oxysporum. Plant Dis. 1982, 66: 809-810. 10.1094/PD-66-809.

    Article  Google Scholar 

  9. Jime’nez-Gasco MM, Navas-Cortes JA, Jime’nez-Diaz RM: The Fusarium oxysporum f. sp. ciceris/Cicer arietinum pathosystem: a case study of the evolution of plant-pathogenic fungi into races and pathotypes. Internat Microbiol. 2004, 7: 95-104.

    Google Scholar 

  10. Gupta S, Bhar A, Das S: Understanding the molecular defence responses of host during chickpea–Fusarium interplay: where do we stand?. Funct Plant Biol. 2013,

    Google Scholar 

  11. Nimalkar SB, Harsulkar AM, Giri AP, Sainani MN, Franceshi V, Gupta VS: Differentially expressed gene transcripts in roots of resistant and susceptible chickpea plant (Cicer arietinum L.) upon Fusarium oxysporum infection. Physiol Mol Plant Pathol. 2006, 68: 176-188. 10.1016/j.pmpp.2006.10.003.

    Article  Google Scholar 

  12. Hossain MM, Hossain N, Sultana F, Islam SMN, Islam MS, Bhuiyan MKA: Integrated management of Fusarium wilt of chickpea (Cicer arietinum L.) caused by Fusarium oxysporum f. sp. ciceris with microbial antagonist, botanical extract and fungicide. Afr J Biotechnol. 2013, 12: 4699-4706.

    Article  Google Scholar 

  13. Berrocal-Lobo M, Molina A: Arabidopsis defense response against Fusarium oxysporum. Cell. 2007, 13: 145-150.

    Google Scholar 

  14. Van der Does HC, Duyvesteijn RG, Goltstein PM, van Schie CC, Manders EM, Cornelissen BJ, Rep M: Expression of effector gene SIX1 of Fusarium oxysporum requires living plant cells. Fungal Genet Biol. 2008, 45: 1257-1264. 10.1016/j.fgb.2008.06.002.

    Article  CAS  PubMed  Google Scholar 

  15. Gupta S, Chakraborti D, Rangi RK, Basu D, Das S: A molecular insight into the early events of chickpea (Cicer arietinum) and Fusarium oxysporum f. sp. ciceri (race 1) interaction through cDNA AFLP analysis. Phytopathology. 2009, 99: 1245-1257. 10.1094/PHYTO-99-11-1245.

    Article  CAS  PubMed  Google Scholar 

  16. Gupta S, Chakraborti D, Sengupta A, Basu D, Das S: Primary metabolism of chickpea is the initial target of wound inducing early sensed Fusarium oxysporum f. sp. ciceri race 1. PLoS One. 2010, 5: e9030-10.1371/journal.pone.0009030. doi:10.1371/journal.pone.0009030

    Article  PubMed Central  PubMed  Google Scholar 

  17. Gupta S, Bhar A, Chatterjee M, Das S: Fusarium oxysporum f.sp. ciceri Race 1 induced redox state alterations are coupled to downstream defense signaling in root tissues of chickpea (Cicer arietinum L.). PLoS One. 2013, 8: e73163-10.1371/journal.pone.0073163. doi:10.1371/journal.pone.0073163

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Ashraf N, Ghai D, Barman P, Basu S, Gangisetty N, Mandal MK, Chakraborty N, Datta A, Chakraborty S: Comparative analyses of genotype dependent expressed sequence tags and stress-responsive transcriptome of chickpea wilt illustrate predicted and unexpected genes and novel regulators of plant immunity. BMC Genomics. 2009, 10: e415-10.1186/1471-2164-10-415.

    Article  Google Scholar 

  19. Giri AP, Harsulkar AM, Patankar AG, Gupta VS, Sainani MN, Deshpande VV, Ranjeka PK: Association of induction of protease and chitinase in chickpea roots with resistance to Fusarium oxysporum. f.sp. ciceris. Plant Pathol. 1998, 47: 693-699. 10.1046/j.1365-3059.1998.00299.x.

    Article  CAS  Google Scholar 

  20. Cho S, Muehlbauer FJ: Genetic effect of differentially regulated fungal response genes on resistance to necrotrophic fungal pathogens in chickpea (Cicer arietinum L.). Physiol Mol Plant Pathol. 2004, 64: 57-66. 10.1016/j.pmpp.2004.07.003.

    Article  CAS  Google Scholar 

  21. Flandez-Galvez H, Ford R, Pang ECK, Taylor PWJ: An intraspecific linkage map of chickpea (Cicer arietinum L.) genome based on sequence tagged microsatellite site and resistant gene analogue markers. Theo Appl Genet. 2003, 106: 1447-1456.

    CAS  Google Scholar 

  22. Jorrín JV, Rubiales D, Dumas-Gaudot E, Recorbet G, Maldonado A, Castillejo MA, Curto M: Proteomics: a promising approach to study biotic interaction in legumes.A review. Euphytica. 2006, 147: 37-47. 10.1007/s10681-006-3061-1.

    Article  Google Scholar 

  23. Padliya ND, Cooper B: Mass spectrometry-based proteomics for the detection of plant pathogens. Proteomics. 2006, 6: 4069-4075. 10.1002/pmic.200600146.

    Article  CAS  PubMed  Google Scholar 

  24. Asano T, Makota K, Takumi N: The defense response in Arabidopsis thaliana against Fusarium sporotrichoides. Proc Natl Acad Sci U S A. 2012, 10: 61-

    CAS  Google Scholar 

  25. Castillejo MA, Maldonado AM, Dumas-Gaudot E, Fernandez-Aparico M, Susin R, Diego R, Jorrin JV: Differential expression proteomics to investigate responses and resistance to Orobanche crenata in Medicago trancatula. BMC Genomics. 2009, 10: 294-10.1186/1471-2164-10-294.

    Article  PubMed Central  PubMed  Google Scholar 

  26. Carpentier SC, Panis B, Vertommen A, Swennen R, Sergeant K, Renaut J, Laukens K, Witters E, Samyn B, Devreese B: Proteome analysis of non-model plants: a challenging but powerful approach. Mass Spectrom Rev. 2008, 27: 354-377. 10.1002/mas.20170.

    Article  CAS  PubMed  Google Scholar 

  27. Summerell BA, Sallen B, Leslie JF: A utilitarian approach to Fusarium identification. Plant Dis. 2003, 87: 117-128. 10.1094/PDIS.2003.87.2.117.

    Article  Google Scholar 

  28. Subba P, Barua P, Kumar R, Datta A, Soni KK, Chakraborty S, Chakraborty N: Phosphoproteomic dynamics of chickpea (Cicer arietinum L.) reveals shared and distinct components of dehydration response. J Proteome Res. 2013, 12: 5025-5047. 10.1021/pr400628j.

    Article  CAS  PubMed  Google Scholar 

  29. Chatterjee M, Gupta S, Bhar A, Das S: Optimization of an efficient protein extraction protocol compatible with two-dimensional electrophoresis and mass spectrometry from recalcitrant phenolic rich roots of chickpea (Cicer arietinum L.). Int J Proteomics. 2012, 2012: 536963-doi:10.1155/2012/536963

    Article  PubMed Central  PubMed  Google Scholar 

  30. Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254. 10.1016/0003-2697(76)90527-3.

    Article  CAS  PubMed  Google Scholar 

  31. Valledor L, Jorrín J: Back to the basics: Maximizing the information obtained by quantitative two dimensional gel electrophoresis analyses by an appropriate experimental design and statistical analyses. J Proteomics. 2011, 74: 1-18. 10.1016/j.jprot.2010.07.007.

    Article  CAS  PubMed  Google Scholar 

  32. Shevchenko A, Tomas H, Havliˇs J, Olsen JV, Mann M: In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. 2007, 1: 2856-2860. 10.1038/nprot.2006.468.

    Article  Google Scholar 

  33. Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio–the analysis and navigation of molecular networks. Bioinformatics. 2003, 19: 2155-2157. 10.1093/bioinformatics/btg290.

    Article  CAS  PubMed  Google Scholar 

  34. Garg R, Sahoo A, Tyagi AK, Jain M: Validation of internal control genes for quantitative gene expression studies in chickpea (Cicer arietinum L.). Biochem Biophys Res Comm. 2010, 396: 283-288. 10.1016/j.bbrc.2010.04.079.

    Article  CAS  PubMed  Google Scholar 

  35. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C (T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.

    Article  CAS  PubMed  Google Scholar 

  36. Chattopadhyay A, Subba P, Pandey A, Bhushan D, Kumar R, Datta A, Chakraborty S, Chakraborty N: Analysis of the grasspea proteome and identification of stress-responsive proteins upon exposure to high salinity, low temperature, and abscisic acid treatment. Phytochemistry. 2011, 72: 1293-1307. 10.1016/j.phytochem.2011.01.024.

    Article  CAS  PubMed  Google Scholar 

  37. Subba P, Kumar R, Gayali S, Shekhar S, Parveen S, Pandey A, Datta A, Chakraborty S, Chakraborty N: Characterisation of the nuclear proteome of a dehydration-sensitive cultivar of chickpea and comparative proteomic analysis with a tolerant cultivar. Proteomics. 2013, 13: 1973-1992. 10.1002/pmic.201200380.

    Article  CAS  PubMed  Google Scholar 

  38. Schlüter H, Apweiler R, Holzhütter HG, Jungblut PR: Finding one’s way in proteomics: a protein species nomenclature. Chem Cent J. 2009, 3: 11-10.1186/1752-153X-3-11.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Zelko IN, Mariani TJ, Folz RJ: Superoxide dismutase multigene family: a comparison of the CuZn-SOD (SOD1), Mn-SOD (SOD2), and EC-SOD (SOD3) gene structures, evolution, and expression. Free Radic Biol Med. 2002, 33: 337-349. 10.1016/S0891-5849(02)00905-X.

    Article  CAS  PubMed  Google Scholar 

  40. Pandey A, Choudhary MK, Bhushan D, Chattopadhyay A, Chakraborty S, Datta A, Chakraborty N: The nuclear proteome of chickpea (Cicer arietinum L.) reveals predicted and unexpected proteins. J Proteome Res. 2006, 5: 3301-3311. 10.1021/pr060147a.

    Article  CAS  PubMed  Google Scholar 

  41. Aboul-Soud MAM, Chen X, Kang JG, Yun BW, Raja MU, Malik SI, Loake GJ: Activation tagging of ADR2 conveys a spreading lesion phenotype and resistance to biotrophic pathogens. New Phytol. 2009, 183: 1163-1175. 10.1111/j.1469-8137.2009.02902.x.

    Article  CAS  PubMed  Google Scholar 

  42. Lu H, Rate DN, Song JT, Greenberg JT: ACD6, a novel ankyrin protein, is a regulator and an effector of salicylic acid signaling in the Arabidopsis defense response. Plant Cell. 2003, 15: 2408-2420. 10.1105/tpc.015412.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Qiu JL, Zhou L, Yun BW, Nielsen HB, Fiil BK, Petersen K, MacKinlay J, Loake GJ, Mundy J, Morris PC: Arabidopsis mitogen-activated protein kinase kinases MKK1 and MKK2 have overlapping functions in defense signaling mediated by MEKK1, MPK4, and MKS1. Plant Physiol. 2008, 148: 212-222. 10.1104/pp.108.120006.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Ferrari S, Plotnikova JM, Lorenzo GD, Ausubel FM: Arabidopsis local resistance to Botrytis cinerea involves salicylic acid and camalexin and requires EDS4 and PAD2, but not SID2, EDS5 or PAD4. Plant J. 2003, 35: 193-205. 10.1046/j.1365-313X.2003.01794.x.

    Article  CAS  PubMed  Google Scholar 

  45. Vorwerk S, Schiff C, Santamaria M, Koh S, Nishimura M, Vogel J, Somerville C, Somerville S: EDR2 negatively regulates salicylic acid-based defenses and cell death during powdery mildew infections of Arabidopsis thaliana. BMC Plant Biol. 2007, 7: 35-10.1186/1471-2229-7-35. doi:10.1186/1471-2229-7-35

    Article  PubMed Central  PubMed  Google Scholar 

  46. Zhang Y, Cheng YT, Qu N, Zhao Q, Bi D, Li X: Negative regulation of defense responses in Arabidopsis by two NPR1 paralogs. Plant J. 2006, 48: 647-656. 10.1111/j.1365-313X.2006.02903.x.

    Article  CAS  PubMed  Google Scholar 

  47. Laxalt AM, Raho N, Have AT, Lamattina L: Nitric oxide is critical for inducing phosphatidic acid accumulation in xylanase-elicited tomato cells. J Biol Chem. 2007, 282: 21160-21168. 10.1074/jbc.M701212200.

    Article  CAS  PubMed  Google Scholar 

  48. Nandi A, Krothapalli K, Buseman CM, Li M, Welti R, Enyedi A, Shah J: Arabidopsis sfd mutants affect plastidic lipid composition and suppress dwarfing, cell death, and the enhanced disease resistance phenotypes resulting from the deficiency of a fatty acid desaturase. Plant Cell. 2003, 15: 2383-2398. 10.1105/tpc.015529.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Sandhu D, Tasma IM, Frasch R, Bhattacharyya MK: Systemic acquired resistance in soybean is regulated by two proteins, orthologous to Arabidopsis NPR1. BMC Plant Biol. 2009, 9: 105-10.1186/1471-2229-9-105. doi:10.1186/1471-2229-9-105

    Article  PubMed Central  PubMed  Google Scholar 

  50. Lee JH, Kim WT: Molecular and biochemical characterization of VR-EILs encoding mung bean ETHYLENE INSENSITIVE3-LIKE proteins. Plant Physiol. 2003, 132: 1475-1488. 10.1104/pp.103.022574.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Demkura PV, Abdala G, Baldwin IT, Ballare CL: Jasmonate-dependent and independent pathways mediate specific effects of solar ultraviolet-B radiation on leaf phenolics and antiherbivore defense. Plant Physiol. 2010, 152: 1084-1095. 10.1104/pp.109.148999.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  52. Skibbe M, Qu N, Galis I, Baldwin IT: Induced plant defenses in the natural environment: Nicotiana attenuata WRKY3 and WRKY6 coordinate responses to herbivory. Plant Cell. 2008, 20: 1984-2000. 10.1105/tpc.108.058594.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  53. Coram TE, Settles ML, Chen X: Transcriptome analysis of high temperature adult plant resistance conditioned by Yr39 during wheat-Pucinia striiformis f.sp. tritici interaction. Molecular Plant Pathol. 2008, 9: 479-493. 10.1111/j.1364-3703.2008.00476.x.

    Article  CAS  Google Scholar 

  54. Matyalman D, Mert Z, Baykal AT, Inan C, Gunel A, Hasancebi S: Proteomic analysis of early responsive resistance proteins of wheat (Triticum aestivum) to yellow rust (Puccinia striiformis F.Sp. tritici) using ProteomeLab PF2D. Plant Omics J. 2013, 6: 24-35.

    Google Scholar 

  55. Alscher RG, Erturk N, Heath LS: Role of superoxide dismutases (SODs) in controlling oxidative stress in plants. J Exp Bot. 2002, 53: 1331-1341. 10.1093/jexbot/53.372.1331.

    Article  CAS  PubMed  Google Scholar 

  56. Takemoto D, Hardham AR, Jones DA: Differences in cell death induction by phytophthora elicitins are determined by signal components downstream of MAP Kinase Kinase in different species of Nicotiana and Cultivars of Brassica rapa and Raphanus sativus. Plant Physiol. 2005, 138: 1491-1504. 10.1104/pp.104.058388.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Narendra S, Venkataramani S, Shen G, Wang J, Pasapula V, Lin Y, Kornyeyev D, Holaday AS, Zhang H: The Arabidopsis ascorbate peroxidase 3 is a peroxisomal membrane-bound antioxidant enzyme and is dispensable for Arabidopsis growth and development. J Exp Bot. 2006, 57: 3033-3042. 10.1093/jxb/erl060.

    Article  CAS  PubMed  Google Scholar 

  58. Panchuk II, Volkov RA, Schoffl F: Heat stress-and heat shock transcription factor-dependent expression and activity of ascorbate peroxidase in Arabidopsis. Plant Physiol. 2002, 129: 838-853. 10.1104/pp.001362.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Dean JD, Goodwin PH, Hsiang T: Induction of glutathione S-transferase genes of Nicotiana benthamiana following infection by Colletotrichum destructivum and C. orbiculare and involvement of one in resistance. J Exp Bot. 2005, 56: 1525-1533. 10.1093/jxb/eri145.

    Article  CAS  PubMed  Google Scholar 

  60. Vieira Dos Santos C, Rey P: Plant thioredoxins are key actors in the oxidative stress response. Trends Plant Sci. 2006, 11: 329-334. 10.1016/j.tplants.2006.05.005.

    Article  CAS  PubMed  Google Scholar 

  61. Sasaki-Sekimoto Y, Taki N, Obayashi T, Aono M, Matsumoto F, Sakurai N, Suzuki H, Hirai MY, Noji M, Saito K, Masuda T, Takamiya K, Shibata D, Ohta H: Coordinated activation of metabolic pathways for antioxidants and defence compounds by jasmonates and their roles in stress tolerance in Arabidopsis. Plant J. 2005, 44: 653-668. 10.1111/j.1365-313X.2005.02560.x.

    Article  CAS  PubMed  Google Scholar 

  62. Hodges DM, Forney CF: The effects of ethylene, depressed oxygen and elevated carbon dioxide on antioxidant profiles of senescing spinach leaves. J Exp Bot. 2000, 51: 645-655. 10.1093/jexbot/51.344.645.

    Article  CAS  PubMed  Google Scholar 

  63. Zhu H, Li GJ, Ding L, Cui X, Berg H, Assmann SM, Xia Y: Arabidopsis extra large G-protein 2 (XLG2) interacts with the G beta subunit of heterotrimeric G protein and functions in disease resistance. Mol Plant. 2009, 2: 513-525. 10.1093/mp/ssp001.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Fan LM, Zhang W, Chen JG, Taylor JP, Jones AM, Assmann SM: Abscisic acid regulation of guard-cell K+ and anion channels in G beta- and RGS-deficient Arabidopsis lines. Proc Natl Acad Sci. 2008, 105: 8476-8481. 10.1073/pnas.0800980105.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Gorecka KM, Thouverey C, Buchet R, Pikula S: Potential role of annexin AnnAt1 from Arabidopsis thaliana in pH-mediated cellular response to environmental stimuli. Plant Cell Physiol. 2007, 48: 792-803. 10.1093/pcp/pcm046.

    Article  CAS  PubMed  Google Scholar 

  66. Laohavisit A, Mortimer JC, Demidchik V, Coxon KM, Stancombe MA, Macpherson N, Brownlee C, Hofmann A, Webb AAR, Miedema H, Battey NH, Davies JM: Zea mays annexins modulate cytosolic free Ca and generate a Ca2+ -permeable conductance. Plant Cell. 2009, 21: 479-493. 10.1105/tpc.108.059550.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Ghelis T, Dellis O, Jeannette E, Bardat F, Miginiac E, Sotta B: Abscisic acid plasmalemma perception triggers a calcium influx essential for RAB18 gene expression in Arabidopsis thaliana suspension cells. FEBS Lett. 2000, 483: 67-70. 10.1016/S0014-5793(00)02088-3.

    Article  CAS  PubMed  Google Scholar 

  68. Hallouin M, Ghelis T, Brault M, Bardat F, Cornel D, Miginiac E, Rona JP, Sotta B, Jeannette E: Plasmalemma abscisic acid perception leads to RAB18 expression via phospholipase D activation in Arabidopsis suspension cells. Plant Physiol. 2002, 130: 265-272. 10.1104/pp.004168.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  69. Lee Y, Lee HS, Lee JS, Kim SK, Kim SH: Hormone and light-regulated nucleocytoplasmic transport in plants: current status. J Exp Bot. 2008, 59: 3229-3245. 10.1093/jxb/ern200.

    Article  CAS  PubMed  Google Scholar 

  70. Dhar-Chowdhury P, Harrell MD, Han SY, Jankowska D, Parachuru L, Morrissey A, Srivastava S, Liu W, Malester B, Yoshida H, Coetzee WA: The glycolytic enzymes, glyceraldehyde-3-phosphate dehydrogenase, triose-phosphate isomerase, and pyruvate kinase are components of the K (ATP) channel macromolecular complex and regulate its function. J Biol Chem. 2005, 280: 38464-38470. 10.1074/jbc.M508744200.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Tovar-Mendez A, Miernyk JA, Randall DD: Regulation of pyruvate dehydrogenase complex activity in plant cells. Eur J Biochem. 2003, 270: 1043-1049. 10.1046/j.1432-1033.2003.03469.x.

    Article  CAS  PubMed  Google Scholar 

  72. Ga’lvez L, Gonza’lez EM, Arrese-Igor C: Evidence for carbon flux shortage and strong carbon/nitrogen interactions in pea nodules at early stages of water stress. J Exp Bot. 2005, 56: 2551-2561. 10.1093/jxb/eri249.

    Article  Google Scholar 

  73. Nunes-Nesi A, Fernie AR, Stitt M: Metabolic and signaling aspects underpinning the regulation of plant carbon nitrogen interactions. Mol Plant. 2010, 3: 973-996. 10.1093/mp/ssq049.

    Article  CAS  PubMed  Google Scholar 

  74. Padmanaban S, Lin X, Perera I, Kawamura Y, Sze H: Differential expression of vacuolar H+ -ATPase subunit c genes in tissues active in membrane trafficking and their roles in plant growth as revealed by RNAi. Plant Physiol. 2004, 134: 1514-1526. 10.1104/pp.103.034025.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  75. Lee A, Kirichenko A, Vygodina T, Siletsky SA, Das TK, Rousseau DL, Gennis R, Konstantinov AA: Ca(2+)-binding site in Rhodobacter sphaeroides cytochrome C oxidase. Biochemistry. 2002, 41: 8886-8898. 10.1021/bi020183x.

    Article  CAS  PubMed  Google Scholar 

  76. Shindo T, Misas-Villamil JC, Horger AC, Song J, van der Hoorn RA: A role in immunity for Arabidopsis cysteine protease RD21, the ortholog of the tomato immune protease C14. PLoS One. 2012, 7: e29317-10.1371/journal.pone.0029317. doi:10.1371/journal.pone.0029317

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  77. Ho C-L, Noji M, Saito M, Yamazaki M, Saito K: Molecular characterization of plastidic phosphoserine aminotransferase in serine biosynthesis from Arabidopsis. Plant J. 1998, 16: 443-452. 10.1046/j.1365-313x.1998.00313.x.

    Article  CAS  PubMed  Google Scholar 

  78. Shen B, Li C, Tarczynski MC: High free-methionine and decreased lignin content result from a mutation in the Arabidopsis S-adenosyl-Lmethionine synthetase 3 gene. Plant J. 2002, 29: 371-380. 10.1046/j.1365-313X.2002.01221.x.

    Article  CAS  PubMed  Google Scholar 

  79. Aoki T, Akashi T, Ayabe SI: Flavonoids of leguminous plants: structure, biological activity, and biosynthesis. J Plant Res. 2000, 113: 475-488. 10.1007/PL00013958.

    Article  Google Scholar 

  80. Pelloux J, Ruste’ rucci C, Mellerowicz EJ: New insights into pectin methylesterase structure and function. Trends Plant Sci. 2007, 12: 267-277. 10.1016/j.tplants.2007.04.001.

    Article  CAS  PubMed  Google Scholar 

  81. Yao C, Wu Y, Nie H, Tang D: RPN1a, a 26S proteasome subunit, is required for innate immunity in Arabidopsis. Plant J. 2012, 71: 1015-1028. 10.1111/j.1365-313X.2012.05048.x.

    Article  CAS  PubMed  Google Scholar 

  82. Suty L, Lequeu J, Lançon A, Etienne P, Petitot AS, Blein JP: Preferential induction of 20S proteasome subunits during elicitation of plant defense reactions: towards the characterization of “plant defense proteasomes”. Int J Biochem Cell Biol. 2003, 35: 637-650. 10.1016/S1357-2725(02)00386-2.

    Article  CAS  PubMed  Google Scholar 

  83. Moffatt BA, Wang L2, Allen MS, Stevens YY, Qin W, Snider J, von Schwartzenberg K: Adenosine kinase of Arabidopsis. Kinetic properties and gene expression. Plant Physiol. 2000, 124: 1775-1785. 10.1104/pp.124.4.1775.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  84. Lange PR, Geserick C, Tischendorf G, Zrenner R: Functions of chloroplastic adenylate kinases in Arabidopsis. Plant Physiol. 2008, 146: 492-504.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  85. Hopkins MT, Lampi Y, Wang TW, Liu Z, Thompson JE: Eukaryotic translation initiation factor 5A is involved in pathogen-induced cell death and development of disease symptoms in Arabidopsis. Plant Physiol. 2008, 148: 479-489. 10.1104/pp.108.118869.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  86. Ondzighi CA, Christopher DA, Cho EJ, Chang SC, Staehelin LA: Arabidopsis protein disulfide isomerase-5 inhibits cysteine proteases during trafficking to vacuoles before programmed cell death of the endothelium in developing seeds. Plant Cell. 2008, 20: 2205-2220. 10.1105/tpc.108.058339.

    Article  CAS  Google Scholar 

  87. Rigobello MP, Donella-Deana A, Cesaro L, Bindoli A: Distribution of protein disulphide isomerase in rat liver mitochondria. Biochem J. 2001, 356: 567-570. 10.1042/0264-6021:3560567.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  88. Lytle BL, Song J, de la Cruz NB, Peterson FC, Johnson KA, Bingman CA, Phillips GN, Volkman BF: Structures of two Arabidopsis thaliana major latex proteins represent novel helix-grip folds. Proteins. 2009, 76: 237-243. 10.1002/prot.22396.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  89. Carter C, Thornburg RW: Tobacco nectarin I. Purification and characterization as a germin-like, manganese superoxide dismutase implicated in the defense of floral reproductive tissues. J Biol Chem. 2000, 275: 36726-36733. 10.1074/jbc.M006461200.

    Article  CAS  PubMed  Google Scholar 

  90. Day B, Henty JL, Porter KJ, Staiger CJ: The pathogen-actin connection: a platform for defense signaling in plants. Annu Rev Phytopathol. 2011, 49: 483-506. 10.1146/annurev-phyto-072910-095426.

    Article  CAS  PubMed  Google Scholar 

  91. Aroca R, Amodeo G, Ferna’ndez-Illescas S, Herman EM, Chaumont F, Chrispeels MJ: The role of aquaporins and membrane damage in chilling and hydrogen peroxide induced changes in the hydraulic conductance of maize roots. Plant Physiol. 2005, 137: 341-353. 10.1104/pp.104.051045.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  92. Cândido Ede S, Pinto MF, Pelegrini PB, Lima TB, Silva ON, Pogue R, Grossi-de-Sá MF, Franco OL: Plant storage proteins with antimicrobial activity: novel insights into plant defense mechanisms. FASEB J. 2011, 25: 3290-3305. 10.1096/fj.11-184291.

    Article  PubMed  Google Scholar 

  93. Varshney RK, Song C, Saxena RK, Azam S, Yu S, Sharpe AG, Cannon S, Baek J, Rosen BD, Tar’an B, Millan T, Zhang X, Ramsay LD, Iwata A, Wang Y, Nelson W, Farmer AD, Gaur PM, Soderlund C, Penmetsa RV, Xu C, Bharti AK, He W, Winter P, Zhao S, Hane JK, Carrasquilla-Garcia N, Condie JA, Upadhyaya HD, Luo MC, et al: Draft genome sequence of chickpea (Cicer arietinum) provides a resource for trait improvement. Nat Biotechnol. 2013, 3: 240-246.

    Article  Google Scholar 

  94. Thompson CE, Fernandes CL, de Souza ON, Salzano FM, Bonatto SL, Freitas LB: Molecular modeling of pathogenesis-related proteins of family 5. Cell Biochem Biophys. 2006, 44: 385-394. 10.1385/CBB:44:3:385.

    Article  CAS  PubMed  Google Scholar 

  95. Shibuya N, Minami E: Oligosaccharide signalling for defence responses in plant. Physiol Mol Plant P. 2001, 59: 223-233. 10.1006/pmpp.2001.0364.

    Article  CAS  Google Scholar 

  96. Dopico B1, Lowe AL, Wilson ID, Merodio C, Grierson D: Cloning and characterization of avocado fruit mRNAs and their expression during ripening and low-temperature storage. Plant Mol Biol. 1993, 21: 437-449. 10.1007/BF00028802.

    Article  CAS  PubMed  Google Scholar 

  97. Hadfield KA1, Dang T, Guis M, Pech JC, Bouzayen M, Bennett AB: Characterization of ripening-regulated cDNAs and their expression in ethylene-suppressed charentais melon fruit. Plant Physiol. 2000, 122: 977-983. 10.1104/pp.122.3.977.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  98. Choi DS, Hwang BK: Proteomics and functional analyses of pepper abscisic acid-responsive 1 (ABR1), which is involved in cell death and defense signaling. Plant Cell. 2011, 23: 823-842. 10.1105/tpc.110.082081.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  99. Fang X, Chen W, Xin Y, Zhang H, Yan C, Yu H, Liu H, Xiao W, Wang S, Zheng G, Liu H, Jin L, Ma H, Ruan S: Proteomic analysis of strawberry leaves infected with Colletotrichum fragariae. J Proteomics. 2012, 75: 4074-4090. 10.1016/j.jprot.2012.05.022.

    Article  CAS  PubMed  Google Scholar 

  100. Palomares-Rius JE, Castillo P, Navas-Cortés JA, Jiménez-Díaz RM, Tena M: A proteomic study of in-root interactions between chickpea pathogens: The root-knot nematode Meloidogyne artiellia and the soil-borne fungus Fusarium oxysporum f. sp. ciceris race 5. J Proteomics. 2011, 74: 2034-2051. 10.1016/j.jprot.2011.05.026.

    Article  CAS  PubMed  Google Scholar 

Download references


Authors are thankful to Dr. Suresh C Pande, Dr. Kiran Kumar Sharma and Dr. Pooja Bhatnagar (International Crops Research Institute for the Semi-Arid Tropics, Patancheru, Hyderabad, India) for providing fungal culture and chickpea seeds. All the authors are thankful to Bose Institute for infrastructure. Special thanks are offered to Rajesh Vashisth (Service manager, Bruker Daltonics) for providing technical help for conducting mass spectrometry. The help provided by the Central Instrumentation Facility, Bose Institute on proteomic services is duly acknowledged. Sincere thanks are also extended to Mr. Arup Kumar Dey for his help during green house experiments. Authors also thank Swarnava Das and Sudipta Basu for technical support. Special thanks are reserved for Mr. Binoy Krishna Modak for his help in manuscript formatting. M.C is thankful to Department of Biotechnology, Government of India for financial support. S.G is thankful to Bose Institute for financial assistance. A.B is thankful to Council of Scientific and Industrial Research for financial support.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Sampa Das.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MC and SG contributed in designing the experiments. MC carried out protein extraction, 2-DE gel analysis, and MALDI-TOF MS and MS/MS experiments. AB conducted real time experiments. DC performed the statistical analysis. MC, SG, AB, DC, DB and SD analyzed the data. MC, SG and AB drafted the manuscript. SD and DB edited the manuscript and supervised the work. All authors read and approved the final manuscript.

Moniya Chatterjee, Sumanti Gupta contributed equally to this work.

Electronic supplementary material


Additional file 1: Protein names, abbreviations and TAIR gene IDs. Table containing list of proteins, their abbreviations used for pathway construction and qRT-PCR and TAIR homologous IDs of the identified proteins used as input for network generation. (XLSX 12 KB)


Additional file 2: List of primers designed for the qRT-PCR. List includes the primer pair sequences used for qRT-PCR for identified proteins wth their respective spot IDs. (XLSX 10 KB)


Additional file 3: Comparative analysis of differentially accumulated protein spots in infected chickpea genotypes. Includes details of protein yield, average number of spots, variable spots (Quantitative and qualitative) obtained in control and infected chickpea genotypes (JG62 and WR315) at different time points post Foc1 infection. (XLSX 11 KB)


Additional file 4: Schematic representation of experimental design. A flow chart depicting the experimental design of Foc1 infected root proteome in chickpea plants. Two weeks old seedlings were infected with Foc1. Root tissues were harvested and 250 μg of proteins were extracted from pooled root tissue to run gels for each time points. The experiments were repeated three times to generate three biological replicates. The gels were stained with coomassie blue and further processed for downstream analyses (In total 72 reproducible 2DE gels were generated). Three technical replicates from three biological replicates were used for PD quest analysis. Differential spots were picked, trypsinized and processed for MALDI- TOF MS and MS/MS. (TIFF 7 MB)


Additional file 5: Protein spots identified by MALDI-TOF MS AND MS/MS. Includes details of differential and unique proteins and their peptides identified by MS and MS/MS. The expression pattern of these proteins in control and infected chickpea genotypes (JG62 and WR315) at different time points post Foc1 infection are also illustrated. (XLSX 654 KB)


Additional file 6: MANOVA Table. Table includes mean protein spot intensities for identified protein spots for control and infected chickpea cultivars (JG62 and WR315) at different time points upon Foc1 infection. Each value represents mean of three repeated experiments each with three replications. The means followed by the same letters within a row do not differ statistically according to Duncan’s multiple range tests at a 5% probability level. (XLSX 20 KB)


Additional file 7: Representative cropped gel images of protein spots belonging to different functional categories. Images show quantitative changes among control and infected plants of both (JG62) and (WR315) chickpea genotypes at different time intervals of 48 h, 72 h, and 96 h after Foc1infection.The number and name indicates the spot identity and name of the proteins mentioned in Additional file 5. Arrow represents the presence of spots. The proteins represented are Pathogenesis related Protein (PR 1), Thaumatin like protein PR- 5b (TLP), Glucan-endo-1,3-beta-glucosidase (BGL2), Trypsin protein inhibitor 3(TrpI-3), Superoxide dismutase (SOD; Mitochondrila manganese SOD), Ascorbate peroxidase (APX), Glutathione S transferase parA(GST), Monodehydro ascorbate reductase (MDAR), Annexin (ANX), ABA-responsive protein (RAB/ABARE), Auxin - induced protein PCNT 115 (Aux ind pro), Triose phosphate isomerase(TIM), Enolase, Isocitrate dehydogenase [NADP] chloroplastic (ICDH), ATP synthase, sub unit D chain (ATPase sub D), S-adenosyl methionine synthetase (SAM), Cysteine proteinase (Cys Pro), Chalcone isomerase(CI), Fructokinase-like protein (FLP), Cytochrome C oxidase subunit 6b-1(COX), Methylesterase1(MER), Adenylate kinase (ADK), 20S proteasome alpha subunit D (20S Prot alpha-D), Protein disulfide-isomerase A6 (PDI), Ripening related protein (RLP/RRP), Germin-like protein (GLP), Profilin-1, Outer plastidial membrane protein porin (Porin), Chain A, Crystal Structure Of A Plant Albumin (Albumin ). (TIFF 5 MB)


Additional file 8: Network showing the total interaction of different upregulated proteins in chickpea obtained after 48 h, 72 h and 96 h post infection with Foc1. Green and pink highlighted components represent the upregulated proteins of WR315 plants and JG62 plants respectively obtained after 48h of infection with Foc1. Blue and orange highlighted components represent the upregulated proteins of WR315 plants and JG62 plants respectively obtained after 72 h of infection with Foc1. Yellow highlighted components represent the upregulated proteins of WR315 plants obtained after 96 h of infection with Foc1. Complete names of protein abbreviations are provided in Additional file 1. (TIFF 12 MB)


Additional file 9: GO classification (biological process and cellular components). Graphical representation of differentially expressed protein spots in chickpea roots (JG62 and WR315) based on network derived (pathway studio version 7.1 software) gene ontology classification. Graphs represent up regulated and down regulated proteins at different time points under different biological processes and cellular components. (TIFF 1 MB)


Additional file 10: GO classification (molecular function). Graphical representation of differentially expressed protein spots in chickpea roots (JG62 and WR315) based on network derived (pathway studio version 7.1 software) gene ontology classification. Graphs represent upregulated and downregulated proteins at different time points under different molecular functions. (TIFF 2 MB)

Authors’ original submitted files for images

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chatterjee, M., Gupta, S., Bhar, A. et al. Analysis of root proteome unravels differential molecular responses during compatible and incompatible interaction between chickpea (Cicer arietinum L.) and Fusarium oxysporum f. sp. ciceri Race1 (Foc1). BMC Genomics 15, 949 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: