Identification of genes differentially expressed in a resistant reaction to Mycosphaerella pinodes in pea using microarray technology
© Fondevilla et al; licensee BioMed Central Ltd. 2011
Received: 20 August 2010
Accepted: 13 January 2011
Published: 13 January 2011
Ascochyta blight, caused by Mycosphaerella pinodes is one of the most important pea pathogens. However, little is known about the genes and mechanisms of resistance acting against M. pinodes in pea. Resistance identified so far to this pathogen is incomplete, polygenic and scarce in pea, being most common in Pisum relatives. The identification of the genes underlying resistance would increase our knowledge about M. pinodes- pea interaction and would facilitate the introgression of resistance into pea varieties. In the present study differentially expressed genes in the resistant P. sativum ssp. syriacum accession P665 comparing to the susceptible pea cv. Messire after inoculation with M. pinodes have been identified using a M. truncatula microarray.
Of the 16,470 sequences analysed, 346 were differentially regulated. Differentially regulated genes belonged to almost all functional categories and included genes involved in defense such as genes involved in cell wall reinforcement, phenylpropanoid and phytoalexins metabolism, pathogenesis- related (PR) proteins and detoxification processes. Genes associated with jasmonic acid (JA) and ethylene signal transduction pathways were induced suggesting that the response to M. pinodes in pea is regulated via JA and ET pathways. Expression levels of ten differentially regulated genes were validated in inoculated and control plants using qRT-PCR showing that the P665 accession shows constitutively an increased expression of the defense related genes as peroxidases, disease resistance response protein 39 (DRR230-b), glutathione S-transferase (GST) and 6a-hydroxymaackiain methyltransferase.
Through this study a global view of genes expressed during resistance to M. pinodes has been obtained, giving relevant information about the mechanisms and pathways conferring resistance to this important disease. In addition, the M. truncatula microarray represents an efficient tool to identify candidate genes controlling resistance to M. pinodes in pea.
Legumes are a versatile and inexpensive source of protein for human food and animal feed. In addition, legumes provide numerous environmental benefits that could contribute to the sustainability of agriculture. Legumes are able to symbiotically fix atmospheric nitrogen, improving soil fertility and decreasing N fertilizers needs .
Dry pea is the most produced legume in Europe and the fourth most in the world  and one of the most productive. However, the instability of pea yields, caused mainly by the occurrence of diseases, hampers the expansion of this legume.
Ascochyta blight, caused by Mycosphaerella pinodes (Berk & Blox) Vesterg, the teleomorph of Ascochyta pinodes (Berk & Blox) Jones, is one of the most important pea pathogens . It is widespread throughout the major pea growing areas [4, 5], especially in temperate regions of Europe, North America, Australia and New Zealand  and constitutes the major constraint for the crop after broomrape in the Mediterranean basin . The disease causes 10% yield losses as an average and can reach 50% under some conditions . The use of resistant varieties would be the most efficient, economical and ecologically strategy to control the disease. However, pea varieties resistant to M. pinodes are not available.
Complete resistance to M. pinodes has not been identified so far. Although extensive searches have been carried out, only moderate resistance has been reported in the cultivated pea [8–10] and this has been inadequate to control the disease. Higher levels of resistance have been identified in wild species of Pisum[8, 11, 12], but their use in breeding programs is hampered by the polygenic nature of resistance.
The identification of the genes controlling resistance to M. pinodes in these wild resistant accessions would facilitate their introgression into pea varieties but these genes are difficult to identify by traditional approaches. Quantitative Trait Loci (QTL) analysis have identified numerous genomic regions involved in resistance to this disease in pea [13–17]. In addition, candidate genes approaches and comparative mapping have revealed the co-localization of QTLs for resistance to M. pinodes and resistance gene analogs, the putative transcription factor PsDof1 and the pea defensin DRR230-b [18, 16], but still very little is known about the mechanisms of response to M. pinodes in pea at the histological, molecular and biochemical level.
Large scale expression studies would allow the establishment of a global and detailed picture of all genes and metabolic pathways expressed or differentially regulated during M. pinodes-pea interaction and would contribute to the identification of candidate genes implicated in ascochyta blight resistance. However, this approach has never been performed in this pathosystem.
The goal of this study was to identify genes and mechanisms of resistance underlying phenotypic variation in resistance to M. pinodes in pea using microarray technology. The advent of microarray technology has enabled large-scale surveys leading to a more integrated view of gene expression responses . In plant-pathogen interactions microarray studies allow a more comprehensive understanding of molecular responses in the infection process making the elucidation of mechanisms involved in resistance possible . The microarray technology requires prior knowledge of the sequence of the genome, but sequence information of pea is at the moment limited. Therefore, this study has taken advantage of the knowledge and tools developed in the model legume Medicago truncatula. A microarray (Mt16KOLI1Plus)  containing 16,470 different 70 mer oligonucleotides from M. truncatula, that represent all tentative consensus sequences (TCs) of the TIGR M. truncatula Gene Index 5 (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=medicago) is available. In the present study cDNA obtained from resistant and susceptible pea plants inoculated with M. pinodes has been hybridised to this microarray and genes differentially expressed in the resistant genotype during M. pinodes infection have been identified.
In the case of genes less expressed in P665 than in Messire (called down-regulated), the most abundant category was 'Primary metabolism' (25.5%). Categories ' Signal transduction and post-translational regulation', 'Gene expression and RNA metabolism', 'Defense and cell rescue' and 'Cell Wall' formed a second group representing around 6% each. The percentage of genes included in the categories 'Protein synthesis and processing', 'Secondary metabolism and hormone metabolism', 'Membrane transport' or not classified in any category, ranged from 2.9 to 4.4 while the remaining categories were weakly represented.
Most relevant genes differentially expressed in P665 compared to Messire after inoculation with M. pinodes
Repetitive proline-rich cell wall protein 2 precursor
Repetitive proline-rich cell wall protein 1 precursor
Caffeic acid O-methyltransferase
Probable xyloglucan endotransglycosylase
Putative proline-rich protein APG isolog
Peroxisomal copper-containing amine oxidase
Cytochrome P450 78A3
Putative amine oxidase
UDP-glycose flavonoid glycosyltransferase
4-coumarate--CoA ligase-like protein
Cytochrome P-450LXXIA1 (cyp71A1)
Flavanone 3 beta-hydroxylase
Hyoscyamine 6 beta-hydroxylase
Flavonol synthase-like protein
Cytochrome P450 71A21
Tropinone reductase homolog
Amygdalin hydrolase isoform AH I precursor
NAC domain protein NAC1
Ethylene responsive element binding factor-like
EBNA-1 nuclear protein
Promoter-binding factor-like protein
Putative steroid membrane binding protein
GRAS family transcription factor
Transcription factor JERF1
Pathogenesis related transcriptional factor and ERF
Probable cysteinyl-tRNA synthetase
Probable CCCH-type zinc finger protein
Probable C2H2 type zinc finger protein ID1 like
Chromosome chr7 scaffold_31 whole genome shotgun sequence
Putative transcription factor APFI
GATA-binding transcription factor-like protein
Putative CTP synthase
bHLH transcription factor GBOF-1
Glycine-rich RNA binding protein
Probable homeobox protein T9L24.43
Transcriptional regulator AraC family
12-oxophytodienoate reductase (OPR2)
12-oxophytodienoate reductase (OPR2)
Putative wall associated serine/theorine kinase
Calcium/calmodulin-dependent protein kinase CaMK2
Putative protein kinase
Signal recognition particle 54 kDa subunit precursor
Mitogen-activated protein kinase
Bacterial-induced peroxidase precursor
Cationic peroxidase 2 precursor
Lipid transfer protein SDi-9 drought-induced
Bacterial-induced peroxidase precursor
Syringolide-induced protein B13-1-9
Disease resistance response protein 39 precursor
Bacterial-induced peroxidase precursor
Probable glutathione S-transferase
Cyanogenic Beta-Glucosidase Molid 1
Elicitor inducible gene product Nt-SubE80
Disease resistance protein-like
TIR-similar-domain-containing protein TSDC
Putative NBS-LRR type disease resistance protein
Putative resistance protein
Syringolide-induced protein 19-1-5
Epoxide hydrolase homolog
CIC protein cold-inducible
Auxin influx carrier protein
Probable wound-induced protein T9A4.6
Putative 16.9 kDa heat shock protein
Sixteen sequences involved in cell wall were differentially expressed in inoculated P665 plants comparing to Messire. Genes implicated in cell wall reinforcement were in general more expressed in P665 than in Messire while genes involved in cell elongation, wall expansion and wall degradation were less expressed.
Up-regulated genes belonging to this category included those associated with protein and amino acid transport and a putative Na+/H+ antiporter. In contrast, several putative membrane transporter proteins and a gene involved in potassium transport were down regulated.
A high percentage of genes differentially expressed in P665 comparing to Messire corresponded to genes involved in primary metabolism. Most of them (73%) showed down regulation. Of them, numerous sequences corresponded to genes participating in photosynthesis. Other down-regulated genes of this category were involved in mobilization and degradation of carbohydrates, degradation of storage oil and nitrogen metabolism. Interestingly, two sequences encoding NADH-plastoquinone oxidoreductase chain 1 chloroplast were also down regulated.
Genes up-regulated included glucosyltransferases, a probable anthocyanin 5-aromatic acyltransferase, a putative acid phosphatase and genes involved in amino acid and phosphor metabolism.
Secondary metabolism and hormone metabolism
The secondary metabolism plays an important role in the response against pathogens. Several genes associated with the synthesis of antimicrobial compounds and involved in defense were more expressed in P665 compared to Messire after inoculation with M. pinodes. Among them, genes involved in phenylpropanoid, alkaloid and flavonoid metabolism, a sequence encoding phenylalanine ammonia-lyase, 2 sequences encoding lipoxygenases and genes involved in H2O2 production were included. However, a sequence showing similarities with a lipoxygenase, another encoding tropinone reductase homolog and a flavonol synthase-like protein were down regulated.
Gene expression and RNA metabolism
Several transcription factors and binding proteins were differentially regulated in P665 comparing to Messire. Genes up regulated encoded a NAC domain protein, an ethylene responsive element binding factor and transcriptions factors belonging to ERF and GRAS families. Transcription factors belonging to bHLH and GATA family were down regulated.
Signal transduction and post-translational regulation
Interestingly, two sequences encoding a 12-oxophytodienoate reductase (OPR2) were more expressed in P665 than in Messire. Several protein kinases involved in different processes were less expressed.
Eleven sequences included in the 'Miscellaneous' category, according to Journet et al.  showed a different regulation in P665 compared to Messire. Among them there were three up-regulated sequences encoding lectins, proteins that can be involved in defense. Other up-regulated genes encoded a beta-glucosidase and an early light inducible protein. A CLC-b chloride channel protein was down regulated. Genes similar to dermal glycoproteins and a legumin were represented by different sequences that were up-regulated in some cases and down regulated in others.
Defense and cell rescue
As expected, numerous sequences corresponding to genes involved in defense were more expressed in P665 comparing to Messire. Those included, among others, sequences encoding peroxidases, a disease resistance response protein 39 precursor, nine-cis-epoxycarotenoid dioxygenase and glutathione S-transferases. However other proteins with possible roles in defense were less expressed in P665. Those included GA protein, PR-10, ascorbate peroxidase and a putative NBS-LRR type disease resistance protein.
Abiotic stimuli and development
Different proteins related to the response to 'Abiotic stimuli and development' category were also found to be differentially regulated in P665 comparing to Messire. Up-regulated proteins included ripening-related proteins, cold- and wound-inducible proteins. An auxin-induced protein and a putative 16.9 kDa heat shock protein were down-regulated.
qRT-PCR validation of 10 differentially expressed genes according to microarray experiment
Disease resistance response protein 39 precursor
12-oxophytodienoate reductase (OPR2)
chlorophyll a/b-binding protein
ribulose 1 5-bisphosphate carboxylase small subunit
Quantitative RT-PCR techniques was used not only to validate the microarray data, but also to investigate the expression level of 10 selected genes in non-inoculated P665 and Messire plants. That allowed a calculation of the fold changes of different genotype and treatment combinations and hence provided interesting extra-information about the regulation of these genes.
Log2 of normalized expression ratios according to qRT-PCR
Disease resistance response protein 39 (DRR230-b)
12-oxophytodienoic acid 10,10-reductase (OPR1)
Nine-cis-epoxycarotenoid dioxygenase 4 (nced4)
6a-hydroxymaackiain methyltransferase (hmm6)
Ferrodoxin NADP oxidoreductase
Chlorophyll a/b-binding protein
Ribulose 1 5-bisphosphate carboxylase small subunit
12-oxophytodienoic acid 10,10-reductase was down-regulated in P665 as compared to Messire after inoculation according to qRT-PCR. That was caused by a high induction of this gene after inoculation in Messire, while the expression of this sequence was almost not induced after inoculation with M. pinodes in P665. In addition, this sequence was constitutively less expressed in P665 than in Messire.
A glutathione S-transferase like gene was constitutively around 4 times more expressed in P665 than in Messire and was slightly induced after inoculation in both genotypes. As a result, P665 inoculated plants showed a higher level of expression of this gene than those of Messire.
The nced4 gene, encoding a nine-cis-epoxycarotenoid dioxygenase 4, had a similar regulation pattern in both genotypes. Thus, both genotypes possessed a constitutively similar level of expression of the gene and in both genotypes the gene was repressed after inoculation with M. pinodes. Consequently, the level of expression of this gene in P665 was similar to that of Messire after inoculation.
hmm6 gene, encoding a 6a-hydroxymaackiain methyltransferase, was constitutively around 2 times more expressed in P665 than in Messire. After inoculation with M. pinodes this gene was over expressed in both genotypes but more strongly in Messire. As a result this gene showed a similar level of expression in both genotypes after inoculation.
Ferrodoxin NADP oxidoreductase was less expressed in P665 than in Messire after inoculation. This gene was repressed in P665 after inoculation with M. pinodes but induced in Messire. In addition, in P665 control plants this gene was less expressed than in Messire ones.
Messire plants inoculated with M. pinodes showed a higher level of expression of chlorophyll a/b-binding protein than P665 plants. Constitutively, both genotypes had a similar level of expression of the gene but this protein was repressed after inoculation in P665 while it was not differentially regulated in Messire.
Messire plants inoculated with M. pinodes showed also a higher expression level of a GA protein encoding gene than P665. In this case, the gene was also repressed after inoculation in P665 and only slightly induced in Messire. In addition, the level of expression of the gene in control plants was lower in P665 than in Messire.
The gene encoding ribulose 1 5-bisphosphate carboxylase small subunit showed almost no regulation after inoculation with M. pinodes in both genotypes and was not differentially expressed in control plants of both genotypes.
Resistance to M. pinodes in pea is a complex trait. Only incomplete resistance to this disease has been identified and genetic analyses have shown that numerous genomic regions are involved in resistance [13–17]. In addition, the necrotrophic nature of M. pinodes complicates the performance of histological studies to elucidate the mechanisms of resistance acting to this pathogen. As a consequence, very little is known about the genes and mechanisms of resistance conferring resistance to this important disease. The present study offers a global view of genes and metabolic pathways expressed in a resistant interaction with M. pinodes and hence provides an excellent tool to increase our knowledge about pea-M. pinodes interaction and to identify candidate genes useful for marker assisted selection.
Previous studies have given some insight into defence responses induced after infection with M. pinodes or treatment with a M. pinodes elicitor. These studies were performed in susceptible pea lines and showed that the M. pinodes elicitor induced the production of the phytoalexin pisatin, the enzymes chalcone sintase and phenylalanine ammonia-lyase, PR proteins as chitinase and endo-b-1,3-glucanase and the generation of superoxide anion. ATPase activity and polyphosphoinositide metabolisms were also activated. On the other hand, M. pinodes produces two suppressors that inhibit these defence responses [23–28]. The present study is the first report on genes differentially expressed after infection with M. pinodes in a resistant line. We compared gene expression profiling in this resistant line with that of a susceptible line using the microarray technology. This approach can contribute to the identification of the specific genes and mechanisms conferring resistance to M. pinodes in pea.
The microarray technology allows the simultaneous assessment of the expression of thousands of genes, being an excellent tool to characterize, at the transcription level, several processes such as defence response to pathogens. As sequence information in pea is limited, we used a microarray containing 70-mer oligos representing all tentative consensus sequences (TCs) of the TIGR M. truncatula Gene Index 5. We obtained a successful cross-hybridization between pea targets and M. truncatula probes. That result was expected due to the high level of homology and syntheny between these two species [29–31]. A high level of data quality and reproducibility was achieved through the use of tree independent biological replicates and two technical replicates, the use of negative controls and a strict statistical analysis to select the genes differentially expressed.
Comparison between microarray and qRT-PCR results showed common expression kinetics for many of the genes indicating that this microarray experiment is a useful tool to select candidates genes potentially involved in resistance to M. pinodes in pea. However, our results also indicate that the involvement of these candidate genes in resistance to M. pinodes must be verified by qRT-PCR using Pisum sequences. Differences observed between microarray and qRT-PCR may be due to the presence of different gene isoforms or to the cross-hybridization between M. truncatula probes and different Pisum genes having similar sequences. For example, microarray experiment showed that a M. truncatula sequence showing similarities to the gene 12-oxophytodienoate reductase (OPR2) was more expressed in P665 than in Mesire. OPR genes are highly similar to each other. Therefore, to investigate by qRT-PCR which OPR gene was up regulated we used a primer pair based on the sequence of the Pisum gene OPR1. But these primers were also able to amplify the genes OPR2, OPR3, OPR4 and OPR6. So, is possible that we have amplified by qRT-PCR sequences corresponding to different genes belonging to this family showing different expression levels. Therefore, further experiments using primers specific for each OPR gene are needed to discern which of them is more expressed in P665 than in Messire after inoculation with M. pinodes. In the case of the gene 'hmm6, 6a-hydroxymaackiain methyltransferase,'in the microarray there were a M. truncatula sequence (MT000707) showing similarities to this gene that was strongly down regulated in P665 comparing to Messire. However, when primers were designed according to the pea sequence of this gene (NCBI accession U69554.1) and the gene was amplified by qRT-PCR, results obtained showed that this gene was not differentially regulated in P665 comparing to Messire. Therefore it is possible that other pea genes, different from '6a-hydroxymaackiain methyltransferase' but having a sequence similar to the probe MT000707, or other unknown isoforms of the gene are also able to hybridise to this probe interfering in the results. In addition, the probe MT000707 has homology with the fragment of the gene located from 1041 to 1094 pb (accession U69554.1), while for qRT-PCR we used a pair of primers that amplified the region of the gene located between 730 and 809 pb. Therefore, as we have checked different fragments by microarray and qRT-PCR, another possibility is that P665 differs from Messire in the fragment of the gene corresponding to the probe MT000707. That hypothesis could be further clarified by amplifying by qRT-PCR the fragment of the gene located between 1041 and 1094 pb.
Plants express a wide range of defence responses that can contribute to resistance to pathogens. These include preformed structural and chemical components, activation of the phytoalexin biosynthetic pathway, production of PR proteins, cell wall reinforcement mediated by hydrogen peroxide and detoxification of fungal toxins. Our results suggest that several of these mechanisms may contribute to resistance to M. pinodes in pea accession P665.
In interactions with susceptible genotypes, M. pinodes spores germinate producing a germ tube and penetrate the pea cuticle directly through the wall of the epidermal cells. Beneath the cuticle, hyphae grow within the outer wall of the epidermis being predominantly aligned with the longitudinal axis of the epidermal cells. Subsequently hyphae grow further within the periplasmic space between plasmalemma and wall, displacing the cell contents, but not causing collapse of the protoplast. This probably biotrophic phase is followed by a necrotrophic one in which M. pinodes causes necrotic lesions in the pea mesophyll. These necrotic lesions rapidly spread in susceptible genotypes and hyphae can grow beyond the necrotic zone [32, 33]. Previous histological studies performed by the authors  indicated that in P665 a proportion of M. pinodes infection units were stopped at the epidermal cells not being able to reach the mesophyl. This lower success in establishing colonies was associated with a rapid death of the epidermal cell that was being infected, resembling a hypersensitive response. In addition, those infection units that succeeded penetrating the epidermis and reached the mesophyll in P665 resulted in lesions significantly smaller than those formed in the susceptible line Messire. These results suggest that a battery of resistance mechanisms are acting in P665, starting from a barrier stopping the infection of M. pinodes at the epidermis and further barriers acting after the pathogen has penetrated epidermis and reached the mesophyll restricting the growth of M. pinodes in the mesophyll. Several genes involved in cell wall fortification were found to be more expressed in P665 than in Messire in the microarray experiment. The involvement of wall reinforcement in the resistance to M. pinodes in pea has been suggested by Clulow and Lewis  and Wroth . This cell wall reinforcement could contribute to the development of physical barriers hampering the expansion of M. pinodes within the P665 tissues or reducing the diffusion of pathogenic toxins. In our microarray experiment 'repetitive proline-rich cell wall proteins', that are structural proteins of the primary cell wall involved in cell wall strengthening, and a 'caffeic acid O-methyltransferase' involved in the lignin synthesis were up regulated in P665 comparing to Messire. In addition, enzymes involved in the production of H2O2, such as 'peroxidases' and 'amine oxidase' were found also to be up-regulated. H2O2 is thought to be required for lignification of the cell wall and for the oxidative crosslinking of hydroxyproline-rich glycoproteins in the cell wall [35, 36]. The accumulation of reactive oxygen species (ROS) is also associated with the occurrence of hypersensitive response  that may play a role in the resistance of line P665 to M. pinodes. ROS can also be toxic and inhibit fungal growth  and act as signaling agents in plant defense .
In addition to ROS, other compounds have antimicrobial properties and can contribute to the inhibition of pathogens development. The recognition of a pathogen by the plant activates several defensive responses including the activation of the phenylpropanoid metabolism and the production of phytoalexins . Phenylpropanoids are a group of plant secondary metabolites derived from phenylalanine that have a wide variety of functions, including defense against microbial attacks or other sources of injury. Phytoalexins are plant antibiotics that are synthesized after the plant tissue is exposed to microbial infection. The production of these antifungal compounds has a relevant role in plant defence. Our results indicate that many enzymes involved in their synthesis were higher expressed in P665 than in Messire after inoculation with M. pinodes. Those included 'UDP-glycose flavonoid glycosiltransferase', 'coumarate CoA ligase like protein', 'phenylalanine ammonia-lyase', 'flavanone 3 beta-hydroxylase' and 'hyoscyamine 6 beta-hydroxylase'. However, other enzymes of these pathways as 'flavonol synthase', 'tropinone reductase" and "amygdalin hydrolase isoform AH I precursor" were down regulated.
In addition, qRT-PCR data showed that the enzyme '6a-hydroxymaackiain methyltransferase', which catalyses the last step of the synthesis of pisatin, the main pea phytoalexin, is constitutively at a higher concentration in P665 than in the susceptible cultivar Messire. This enzyme was activated after inoculation with M. pinodes in both genotypes and both genotypes showed similar amount of this enzyme 16hai. The constitutively higher expression of this enzyme in P665, suggests that in P665 pisatin can start acting earlier against the pathogen and can reach the same final level as Messire with a lower effort by the plant.
Pathogenesis- related (PR) proteins are also induced during infection by pathogens and several of them possess antimicrobial properties . The PR14 are 'lipid transfer protein' and a sequence encoding such a protein was up regulated in P665. Another up-regulated PR protein encoded a precursor of the defensin 'disease resistance response 39 (DRR230-b)'. This gene was induced after infection with M. pinodes and showed a constitutively higher expression in P665 comparing to Messire. DRR230-b defensin was first identified by Chian and Hadwiger  from pea pods in response to infection by the fungal pathogen Fusarium solani. The gene encoding 'disease resistance response 39 precursor' was also present in a cDNA library obtained from a resistant Lathyrus sativus accession inoculated with M. pinodes. More recently, this defensin was found to co-localize with the QTL mpIII-4 involved in field resistance to M. pinodes in pea . In addition, the related defensins DRR230-a and DRR230-c were also found to be induced after infection with several pathogens including A. pinodes (the teleomorph of M. pinodes) . Our results reinforce these recent studies suggesting the important role of this protein in resistance to diseases in pea, and specially in resistance to M. pinodes.
Necrotrophic fungi, as M. pinodes, kill host tissues during infection, usually through the secretion of toxic substances. Therefore, the ability of a plant to detoxify these fungal toxins may contribute to resistance to necrotic pathogens. Thus, chickpea cultivars with higher sensitivity to the phytoxins produced by Ascochyta rabiei are more susceptible to this pathogen . Our results show that P665 may posses a higher ability to detoxify M. pinodes toxins as two genes involved in detoxification processes, the 'glutatione S-transferase' and 'ABC transporter', were found to be more expressed in P665 than on Messire. Glutatione S-transferases are involved in several metabolic processes and in the detoxification of a wide variety of compounds including microbial toxins . ATP-binding cassette transporters (ABC-transporter) are transmembrane proteins that function in the transport of a wide variety of substrates across extra- and intracellular membranes including toxins, drugs, glutatione conjugates, peptides and secondary metabolites [47, 48].
In addition to the genes reported above, other genes involved in defence were also up-regulated in P665. These included a "syringolide-induced protein" that have been found to be induced after treatment with the syringolide elicitors produced by the bacteria Pseudomonas syringae. Other up-regulated genes with a possible involvement in defence were genes encoding lectins, as several plant lectins have been shown to induce the production of pisatin .
Perception of both general and specific pathogen associated molecules triggers defence responses via signal transduction cascades and transcriptional activation of numerous genes . The expression of transcription factors and proteins kinases, as well as elevation of cytosolic calcium, is integral to the signalling of these defences . We identified several genes involved in signal recognition and transduction pathways, such as kinases, CCCH-type zinc finger protein and transcription factors, that were differently expressed in P665 comparing to Messire after inoculation with M. pinodes. Among the differentially regulated transcriptional factors there were some associated with Jasmonic Acid (JA) and Ethylene (ET): "ethylene responsive element binding factor-like", "transcription factor JERF1" and "pathogenesis related transcriptional factor ERF". This suggests that the response to M. pinodes in pea is regulated via JA and ET pathways. This is in agreement with the predominant necrotrophic nature of M. pinodes, as gene-for gene resistance and SA signalling are generally effective against biotrophic pathogens whereas JA/ET signalling is generally effective against necrotrophs .
In addition to genes involved in defence against pathogens also genes involved in response to abiotic stresses and development such as 'ripening-related protein-like', 'CIC protein cold-inducible' and 'wound-induced protein T9A4.6' were more expressed in P665 than in Messire showing that response to abiotic and biotic stresses and proteins involved in development are interlinked, as many other studies suggest. For example a "ripening related protein" was found also to be expressed in a pea line resistant to Erysiphe pisi and in the model legume legume Medicago truncatula in response to the parasitic plant Orobanche crenata.
Our results suggest that resistance to M. pinodes in P665 is in part due to a constitutively higher expression of genes involved in defense such as peroxidases, DRR230-b, GST and 6a-hydroxymaackiain methyltransferase. The first step in the response to a pathogen is the recognition of the pathogen by the plant. This recognition leads to the induction of the defence responses. In gene-for-gene resistance, early recognition of specific pathogen strains, a key step in a successfully defense, depends on complementary pairs of dominant genes, one in the host and one in the pathogen. Gene-for-gene resistance is common in interactions with many biotrophic pathogens . In contrast, resistance mediated by a single host resistance gene is uncommon in the case of necrotrophic fungal pathogens. In the case of necrotrophic pathogens plants usually recognize non-specific elicitors that activate a battery of basal defense responses that act against a wide range of pathogens. In this case, as is the case of resistance to the necrotrophic fungi M. pinodes, the recognition of a pathogen is not so fast and a preformed higher expression of genes with antimicrobial properties can be an advantage to get a fast and effective defence response.
In this study, we have obtained a global view of genes expressed during resistance to M. pinodes. This gave us information about the possible mechanisms and pathways involved in the resistance to this important disease such as cell wall reinforcement, production of phytoalexins, phenylpropanoids and PR proteins and detoxification of fungal toxins. This study is also an useful tool to identify candidates genes involved in the control of resistance to M. pinodes in pea useful for marker assisted selection. Further studies will include the mapping of the most relevant genes identified in this study in a RIL population derived from the cross P665 x Messire where QTLs associated with resistance to M. pinodes have been identified and functional analysis to discern the role of these genes in resistance.
Plant material and inoculation
Two pea genotypes, P665 and Messire, the parental lines of a RIL population previously used to identify QTLs associated with resistance to M. pinodes were used in the experiment. Messire is a commercial Pisum sativum ssp. sativum cultivar highly susceptible to M. pinodes. P665 is a P. sativum ssp. syriacum accession displaying incomplete resistance to M. pinodes. Previous histological studies revealed that resistance to M. pinodes in accession P665 was characterized by a lower succeed in colony establishment, associated with the rapid death of the epidermal cell being attacked by M. pinodes and by a smaller colony size .
For inoculation plants were grown until the fifth leave stage in a growth chamber (20 ± 2°C with a 12 h dark/12 h light photoperiod, at 250 μmol m-2 sec-1). Plants were inoculated with the monoconidial M. pinodes isolate C0-99, obtained from infected pea material collected in commercial fields at Córdoba, Spain. The isolate was multiplied in Petri dishes containing V8 juice medium located in a growth chamber at 21 ± 2°C with a 12 h dark/12 h light photoperiod, at 106 lmol/m2 s. A spore suspension was prepared by flooding the surface of 12 days old cultures with sterile water, scraping the colony with a needle and filtering the suspension through two layers of sterile cheesecloth. The concentration of spores in the solution obtained was further determined with a haemocytometer and adjusted to 350.000 spores per ml. Finally, Tween-20 (120 μl per 100 ml of suspension) was added as a wetting agent and the spore suspension was applied with a sprayer at a rate of 1 ml per plant. After inoculation high humidity was ensured during the first 24 h by ultrasonic humidifiers operating for 15 minutes every two hours. After that period the humidifiers were turned off.
The experiment was performed in three independent replicates, each having 3 to 5 plants per genotype (Messire/P665), treatment (inoculated/control) and time of harvesting (16, 24 and 48 hours after inoculation). In each replicate, plants grown under the same conditions but not inoculated were used as control.
Sample collection and RNA extraction
At 16, 24 and 48 hours after inoculation (hai) leaflets of control and inoculated plants were harvested, immediately frozen in liquid nitrogen and stored at -80°C. RNA was isolated using Trizol reagent (Invitrogen, Karlsruhe, Germany) according to manufacture's protocols. Integrity of total RNA was checked on agarose gels and its quantity, as well as purity, was determined using NanoDrop ND1000 (NanoDrop Technologies, Inc., Wilmington, USA). RNA from infected plants was further purified and concentrated to 0.8 μg/μl using Microcon-30 YM columns (Millipore, Schwalbach, Germany).
Microarray experiment was performed at the Institute for Genome Research of Bielefeld University, Germany. For each time of harvesting and replicate, Cy-labelled cDNA samples from resistant and susceptible inoculated plants labelled with different Cy dies were co-hybridized to Mt16kOLI1Plus microarray as described by Küster et al. . The experiment included three biological and two technical replicates incorporating one dye swap. The resulting images were analysed using the ImaGene 5.5 software (Bio-Discovery, Los Angeles) as described by Hohnjec et al., . Data files were imported into the EMMA1.1 array analysis software  and normalized using Lowess normalization. To identify the genes differentially expressed in the inoculated resistant genotype compared to the susceptible one a t-test followed by FDR correction was performed. Genes were considered differentially regulated when p ≤ 0.05 and M ≤ -0.8 or M ≥ 0.8, being M = Log2 (red/green). The microarray data have been deposited into the public data base ArrayExpress (E-TABM-1084).
Data validation by quantitative real time Reverse Transcription PCR (qRT-PCR)
The expression profiles of 10 genes differentially expressed according to the microarray experiment were validated in inoculated and control plants using two steps qRT-PCR. Total RNA was extracted from different samples obtained from the same three replicates used for microarray study using TRISure (Bioline, London, UK). After checking its quality, any possible residual genomic DNA was removed using RQ1 RNase-Free Dnase (Promega, Madison, USA). RNA was further purified using RNeasy Plant Mini Kit (Quiagen, Hilden, Germany). The absence of genomic DNA was checked by PCR using specific primers that amplify and intron-exon-intron sequence of the P. sativum gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Fw: 5'-3': GTGGTCTCCACTGACTTTATTGGT/Rv 5'-3': TTCCTGCCTTGGCATCAAA, Die et al., 2010). Total RNA (5 μg) was reverse-transcribed using SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen, Karlsruhe, Germany).
In order to ensure equal starting cDNA amounts, real-time PCR amplification of α-Tubuline (TUB) was run for all the different templates and, depending on the CT (threshold cycle) number, cDNA samples were diluted to obtain similar CT values. In addition, to check the quality of the reverse transcription, specific primers were used to amplify in each template two fragments of the gene GAPDH located 915 bp apart at the 5'or 3'end of the transcript (GAPDH1 Fw 5'-3': ctccactgactttattggtgaca/Rv 5'-3': caaacttgtcatttaaggcaattc; GAPDH2 Fw 5'-3': tcaagatcggaatcaacggatt/Rv 5'-3': cgagttcaacatcatctctcttcaa).
Polymerase chain reactions were performed in a 96-well plate with a 7500 Real Time PCR System (Applied Biosystems, Foster City, CA, USA), using SYBR Green to monitor dsDNA synthesis. Reactions contained 0.5 μl of Fast Start Universal SYBR Green Master (ROX), 1 μl of cDNA, and 0.3 μM of each gene-specific primer in a final volume of 10 μl. The following standard thermal profile was used for all PCR reactions: polymerase activation (95°C for 10 min), amplification and quantification cycles repeated 40 times (95°C for 15 seconds, 60°C for 1 min) and dissociation curve (95°C 15 seconds, 60°C 1 min, 95°C 30 seconds).
Primers used to amplify 10 candidate genes by qRT-PCR
Forward primer (5'-3')
Reverse primer (5'-3')
Disease resistance response protein 39 (DRR230-b)
12-oxophytodienoic acid 10,10-reductase (OPR1)
Nine-cis-epoxycarotenoid dioxygenase 4 (nced4)
6a-hydroxymaackiain methyltransferase (hmm6)
Ferrodoxin NADP oxidoreductase
Chlorophyll a/b-binding protein
Ribulose 1 5-bisphosphate carboxylase small subunit
The genes TUB, histone H3 and GAPDH  were used as reference genes for normalization.
The PCR efficiency of each primer pair in each individual reaction was calculated using LingRegPCR 7.5 software and used to calculate an average efficiency (E) per primer pair. This average efficiency was used to calculate the expression in each reaction using the formula Expresion = ECT. A normalization index was calculated for each plate as the geometric mean of the expression of the reference genes TUB, GAPDH and histone H3 and a relative expression was calculated for each reaction as the ratio of the gene expression of the gene of interest in each reaction against the normalization index.
The authors thank the Spanish Ministry of Science and Innovation for providing a 'Juan de la Cierva' contract to the senior author and AGL2008-01239 and FP6-CT-2004-FOOD-1-506223 projects for financial support.
- Köpje U, Nemecek T: Ecological services of faba bean. Field Crops Res. 2010, 115: 217-233.View ArticleGoogle Scholar
- FAOSTAT: [http://faostat.fao.org/site/567/default.aspx#ancor]
- Moussart A, Tivoli B, Lemarchand E, Deneufbourg F, Roi S, Sicard G: Role of seed infection by the Ascochyta blight pathogen of dried pea (Mycosphaerella pinodes) in seedling emergence, early disease development and transmission of the disease to aerial plant parts. Eur J Plant Pathol. 1998, 104: 93-102. 10.1023/A:1008673914537.View ArticleGoogle Scholar
- Lawyer SA: Diseases caused by Ascochyta spp. Compendium of pea diseases. Edited by: Hargedon DJ. 1984, St Paul, Minnesota: APS Press, 11-15.Google Scholar
- Wallen VR: Field evaluation of the importance of the Ascochyta complex of peas. Can J Plant Sci. 1965, 45: 27-33. 10.4141/cjps65-004.View ArticleGoogle Scholar
- Rubiales D, Pérez-de-Luque A, Cubero JI, Sillero JC: Crenate broomrape (Orobanche crenata) infection in field pea cultivars. Cro Prot. 2003, 22: 865-872. 10.1016/S0261-2194(03)00070-X.View ArticleGoogle Scholar
- Xue AG, Warkentin TD, Kenaschuk EO: Effect of timings of inoculation with Mycosphaerella pinodes on yield and seed infection on field pea. Can J Plant Sci. 1997, 77: 685-689.View ArticleGoogle Scholar
- Clulow SA, Lewis BG, Matthews P: A pathotype classification for Mycosphaerella pinodes. Phytopathology. 1991, 131: 322-332. 10.1111/j.1439-0434.1991.tb01203.x.View ArticleGoogle Scholar
- Kraft JM: A search for resistance in peas to Mycosphaerella pinodes. Plant Dis. 1998, 82: 251-253. 10.1094/PDIS.19188.8.131.52.View ArticleGoogle Scholar
- Wroth JM: Host- pathogen relationship of the ascochyta bligt (Mycosphaerella pinodes (Berk & Blox) Vesterg ) disease of field pea (Pisum sativum L ). PhD. thesis. 1996, University of Western Australia, PerthGoogle Scholar
- Fondevilla S, Ávila CM, Cubero JI, Rubiales D: Response to Mycosphaerella pinodes in a germplasm collection of Pisum spp. Plant Breeding. 2005, 124: 313-315. 10.1111/j.1439-0523.2005.01104.x.View ArticleGoogle Scholar
- Wroth JM: Possible role for wild genotypes of Pisum spp. to enhance ascochyta bligt resistance in pea. Aust J Exp Agr. 1998, 38: 469-479. 10.1071/EA98024.View ArticleGoogle Scholar
- Fondevilla S, Satovic Z, Rubiales D, Moreno MT, Torres AM: Mapping of quantitative trait loci for resistance to Mycosphaerella pinodes in Pisum sativum subsp. syriacum. Mol Breeding. 2008, 21: 439-454. 10.1007/s11032-007-9144-4.View ArticleGoogle Scholar
- Prioul S, Frankewitz A, Deniot G, Morin G, Baranger A: Mapping of quantitative trait loci for partial resistance to Mycosphaerella pinodes in pea (Pisum sativum L.) at the seedling and adult plant stages. Theor Appl Genet. 2004, 108: 1322-1334. 10.1007/s00122-003-1543-2.PubMedView ArticleGoogle Scholar
- Tar'an B, Warkentin T, Somers DJ, Miranda D, Vandenberg A, Balde S, Woods S, Bing D, Xue A, DeKoeyer D, Penner G: Quantitative trait loci for lodging resistance, plant height and partial resistance to mycosphaerella blight in field pea (Pisum sativum L.). Theor Appl Genet. 2003, 107: 1482-1491.PubMedView ArticleGoogle Scholar
- Timmerman-Vaughan GM, Frew TJ, Russell AC, Khan T, Butler R, Gilpin M, Murray S, Falloon K: QTL mapping of partial resistance to field epidemics of ascochyta blight of pea. Crop Sci. 2002, 42: 2100-2111. 10.2135/cropsci2002.2100.View ArticleGoogle Scholar
- Timmerman-Vaughan GM, Frew TJ, Butler R, Murray S, Gilpin M, Falloon K, Johnston P, Lakeman MB, Russell AC, Khan T: Validation of quantitative trait loci for Ascochyta blight resistance in pea (Pisum sativum L.), using populations from two crosses. Theor Appl Genet. 2004, 109: 1620-1631. 10.1007/s00122-004-1779-5.PubMedView ArticleGoogle Scholar
- Prioul-Gervais S, Deniot G, Receveur EM, Frankewitz A, Fourmann M, Rameau C, Pilet-Nayel ML, Baranger A: Candidate genes for quantitative resistance to Mycosphaerella pinodes in pea (Pisum sativum L.). Theor Appl Genet. 2007, 114: 971-984. 10.1007/s00122-006-0492-y.PubMedView ArticleGoogle Scholar
- Dita MA, Die JV, Román B, Krajinski F, Küster H, Moreno MT, Cubero JI, Rubiales D: Gene expression profiling of Medicago truncatula roots in response to the parasitic plant Orobanche crenata. Weed Res. 2009, 49 (Suppl. 1): 66-80. 10.1111/j.1365-3180.2009.00746.x.View ArticleGoogle Scholar
- Bar-Or C, Kapulnik Y, Koltai H: A broad characterization of the transcriptional profile of the compatible tomato response to the plant parasitic root net nematode Meloidogyne javanica. Eur J Plant Pathol. 2005, 111: 181-192. 10.1007/s10658-004-2134-z.View ArticleGoogle Scholar
- Küster H, Becker A, Firnhaber C, Hohnjec N, Manthey K, Perlick AM, Bekel T, Dondrup M, Henckel K, Goesmann A, Meyer F, Wipf D, Requena N, Hildebrandt I, Hampp R, Nehls U, Krajinski F, Franken P, Pühler A: Development of bioinformatic tools to support EST-sequencing, in silico- and microarray-based transcriptome profiling in mycorrhizal symbioses. Phytochemistry. 2007, 68: 19-32.PubMedView ArticleGoogle Scholar
- Journet EP, van Tuinen D, Gouzy J, Crespeau H, Carreau V, Farmer M, Niebel A, Schiex T, Jaillon O, Chatagnier O, Godiard L, Micheli F, Kahn D, Gianinazzi-Pearson V, Gamas P: Exploring root symbiotic programs in the model legume Medicago truncatula using EST analysis. Nucleic Acids Res. 2002, 30: 5579-5592. 10.1093/nar/gkf685.PubMed CentralPubMedView ArticleGoogle Scholar
- Kiba A, Miyake C, Toyoda K, Ichinose Y, Yamada T, Shiraishi T: Superoxide generation in extracts from isolated plant cell walls is regulated by fungal signal molecules. Phytopathol. 1997, 87: 846-852. 10.1094/PHYTO.19184.108.40.2066.View ArticleGoogle Scholar
- Shiraishi T, Oku H, Tsuji Y, Ouchi S: Inhibitory effect of pisatin on infection process of Mycosphaerella pinodes on pea. Ann Phytopathol Soc Jpn. 1978, 44: 641-645.View ArticleGoogle Scholar
- Shiraishi T, Saitoh K, Mo Kim H, Kato T, Tahara M, Oku H, Yamada T, Ichinose Y: Two suppressors, supprescins A and B, secreted by a pea pathogen, Mycosphaerella pinodes. Plant Cell Physiol. 1992, 33: 663-667.Google Scholar
- Toyoda K, Koyama M, Mizukoshi R, Ichinose Y, Yamada T, Shiraishi T: Phosphorylation of phosphatidylinositols and production of lysophospholipid in pea plasma membrane are coordinately regulated by elicitor and suppressor from Mycosphaerella pinodes. Sci Rep Fac Agr Okayama Univ. 1998, 87: 109-116.Google Scholar
- Yamada T, Shiraishi T, Ichinose Y, Kato H, Seki H, Murakami Y: Regulation of genes for phenylpropanoid synthesis in pea elicitor and suppressor. Molecular aspects of phatogenicity and resistance: requirement for signal transduction. Edited by: Mills D, Kunoh H, Keen NT, Mayama S. 1996, St Paul: American Phytopathological Society, 151-162.Google Scholar
- Yoshioka H, Shiraishi T, Nasu K, Yamada T, Ichinose Y, Oku H: Suppression of activation of chitinase and ß-1,3-glucanase in pea epicotyls by orthovanadate and suppressor from Mycosphaerella pinodes. Ann Phytopathol Soc Jpn. 1992, 58: 405-410.View ArticleGoogle Scholar
- Choi HK, Mun JH, Kim DJ, Zhu H, Baek JM, Mudge J, Roe B, Ellis N, Doyle J, Kiss GB, Young ND, Cook DR: Estimating genome conservation between crop and model legume species. Proc Natl Acad Sci USA. 2004, 101: 15289-15294. 10.1073/pnas.0402251101.PubMed CentralPubMedView ArticleGoogle Scholar
- Gualtieri G, Kulikova O, Limpens E, Kim DJ, Cook DR, Bisseling T, Geurts R: Microsynteny between pea and Medicago truncatula in the SYM2 region. Plant Mol Biol. 2002, 50: 225-235. 10.1023/A:1016085523752.PubMedView ArticleGoogle Scholar
- Kalo P, Seres A, Taylor SA, Jakab J, Kevei Z, Kereszt A, Endre G, Ellis TH, Kiss GB: Comparative mapping between Medicago sativa and Pisum sativum. Mol Genet Genomics. 2004, 272: 235-246. 10.1007/s00438-004-1055-z.PubMedView ArticleGoogle Scholar
- Clulow SA, Lewis BG, Matthews P: Infection of pea epicotyls by Mycosphaerella pinodes. Mycol Res. 1991, 95 (7): 817-820. 10.1016/S0953-7562(09)80044-0.View ArticleGoogle Scholar
- Clulow SA, Lewis BG, Matthews P: Expression of resistance to Mycosphaerella pinodes in Pisum sativum. Plant Pathol. 1992, 41: 362-369. 10.1111/j.1365-3059.1992.tb02359.x.View ArticleGoogle Scholar
- Fondevilla S: Identification and characterization of sources of resistance to Mycosphaerella pinodes in Pisum spp. MSC thesis. 2000, University of Córdoba, Genetic DepartmentGoogle Scholar
- Brady JD, Fry SC: Formation of di-isodityrosine and loss of isodityrosine in the cell walls of tomato cell-suspension cultures treated with fungal elicitors or H2O2. Plant Physiol. 1997, 115: 87-92.PubMed CentralPubMedGoogle Scholar
- Brisson LF, Tenhaken R, Lamb C: Function of oxidative cross-linking of cell wall structural proteins in plant disease resistance. The Plant Cell. 1994, 6: 1703-1712. 10.2307/3869902.PubMed CentralPubMedView ArticleGoogle Scholar
- Torres MA, Jones JDG, Dangl JL: Pathogen-induced, NADPH oxidase-derived reactive oxygen intermediates suppress spread of cell death in Arabidopsis thaliana. Nat Genet. 2005, 37: 1130-1134. 10.1038/ng1639.PubMedView ArticleGoogle Scholar
- Lu H, Higgins VJ: The effect of hydrogen peroxide on the viability of tomato cells and of the fungal pathogen Cladosporium fulvum. Physiol Mol Plant Pathol. 1999, 54: 131-143. 10.1006/pmpp.1998.0195.View ArticleGoogle Scholar
- Lamb C, Dixon RA: The oxidative burst in plant disease resistance. Annu Rev Plant Physiol Plant Mol Biol. 1997, 48: 251-275. 10.1146/annurev.arplant.48.1.251.PubMedView ArticleGoogle Scholar
- Dixon R, Achnine L, Kota P, Lui C, Reddy M, Wang L: The phenylpropanoid pathway and plant defense- a genomics perspective. Mol Plant Pathol. 2002, 3: 371-390. 10.1046/j.1364-3703.2002.00131.x.PubMedView ArticleGoogle Scholar
- van Loon L, van Strien EA: The families of pathogenesis related proteins, their activities and comparative analysis of PR-1 type proteins. Physiol Mol Plant Pathol. 1999, 55: 85-97. 10.1006/pmpp.1999.0213.View ArticleGoogle Scholar
- Chiang CC, Hadwiger LA: The Fusarium solani -induced expression of a pea gene family encoding high cysteine content proteins. Mol Plant Microbe Interact. 1991, 4: 324-331. 10.1094/MPMI-4-324.PubMedView ArticleGoogle Scholar
- Skiba B, Ford R, Pang ECK: Construction of a cDNA library of Lathyrus sativus inoculated with Mycosphaerella pinodes and the expression of potential defence-related expressed sequence tags (ESTs). Physiol Mol Plant Pathol. 2005, 66: 55-67. 10.1016/j.pmpp.2005.04.005.View ArticleGoogle Scholar
- Lai FM, DeLong C, Mei K, Wignes T, Fobert PR: Analysis of the DRR230 family of pea defensins: gene expression pattern and evidence of broad host-range antifungal activity. Plant Sci. 2002, 163: 855-864. 10.1016/S0168-9452(02)00230-3.View ArticleGoogle Scholar
- Jayakumar P, Gossen BD, Gan YT, Wankentin TD, Banniza S: Ascochyta blight of chickpea: infection and host resistance mechanisms. Can J Plant Pathol. 2005, 27: 499-509. 10.1080/07060660509507251.View ArticleGoogle Scholar
- Marrs KA: The functions and regulation of glutathione S-transferases in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 127-158. 10.1146/annurev.arplant.47.1.127.PubMedView ArticleGoogle Scholar
- Theodoulou FL: Plant ABC transporters. Biochim Biophys Acta. 2000, 1465: 79-103. 10.1016/S0005-2736(00)00132-2.PubMedView ArticleGoogle Scholar
- Yazaki K: ABC transporters involved in the transport of plant secondary metabolites. FEBS Letters. 2006, 580: 1183-1191. 10.1016/j.febslet.2005.12.009.PubMedView ArticleGoogle Scholar
- Hagihara T, Hashi M, Takeuchi Y, Yamaoka N: Cloning of soybean genes induced during hypersensitive cell death caused by syringolide elicitor. Planta. 2004, 218: 1432-2048. 10.1007/s00425-003-1136-y.View ArticleGoogle Scholar
- Toyoda K, Miki K, Ichinose Y, Yamada T, Shiraishi T: Plants lectins induce the production of phytoalexin in Pisum sativum. Plant Cell Physiol. 1995, 36 (5): 799-807.Google Scholar
- Coram TE, Pang ECK: Transcriptional profiling of chickpea genes differentially regulated by salicylic acid, methyl jasmonate and aminocyclopropane carboxylic acid to reveal pathways of defense-related gene regulation. Func Plant Biol. 2007, 34: 52-64. 10.1071/FP06250.View ArticleGoogle Scholar
- Grant M, Mansfield J: Early events in host-pathogen interactions. Curr Opin Plant Biol. 2007, 2: 312-319. 10.1016/S1369-5266(99)80055-7.View ArticleGoogle Scholar
- Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005, 43: 205-227. 10.1146/annurev.phyto.43.040204.135923.PubMedView ArticleGoogle Scholar
- Curto M, Camafeita E, Lopez JA, Maldonado AM, Rubiales D, Jorrín V: A proteomic approach to study pea (Pisum sativum) responses to powdery mildew (Erysiphe pisi). Proteomics. 2006, 6: 163-174. 10.1002/pmic.200500396.View ArticleGoogle Scholar
- Die JV, Dita MA, Krajinski F, González-Verdejo CI, Rubiales D, Moreno MT, Román B: Identification by suppression subtractive hybridization and expression analysis of Medicago truncatula putative defense genes in response to Orobanche crenata parasitization. Physiol Mol Plant Pathol. 2007, 70: 49-59. 10.1016/j.pmpp.2007.06.001.View ArticleGoogle Scholar
- Hammond-Kosack KE, Jones JDG: Plant disease resistance genes. Annu Rev Plant Physiol Plant Mol Biol. 1997, 48: 575-607. 10.1146/annurev.arplant.48.1.575.PubMedView ArticleGoogle Scholar
- Küster H, Hohnjec N, Krajinski F, El Yahyaoui F, Manthey K, Gouzy J, Dondrup M, Meyer F, Kalinowski J, Brechenmacher L, van Tuinen D, Gianinazzi-Pearson V, Pühler A, Gamas P, Becker A: Construction and validation of cDNA-based Mt6k-RIT macro- and microarrays to explore root endosymbioses in the model legume Medicago truncatula. J Biotech. 2004, 108: 95-113.View ArticleGoogle Scholar
- Hohnjec N, Vieweg MF, Phuler A, Becker A, Kuster H: Overlaps in the transcriptional profiles of Medicago truncatula roots inoculated with two different glomus fungi provide insights into the genetic program activated during arbuscular mycorriza. Plant Physiol. 2005, 137: 1283-1301. 10.1104/pp.104.056572.PubMed CentralPubMedView ArticleGoogle Scholar
- Dondrup M, Goesmann A, Bartels D, Kalinowski J, Krause L, Linke B, Rupp O, Sczyrba A, Puhler A, Meyer F: EMMA: a platform for consistent storage and efficient analysis of microarray data. J Biotech. 2003, 106: 135-146. 10.1016/j.jbiotec.2003.08.010.View ArticleGoogle Scholar
- Die JV, Román B, Nadal S, González-Verdejo CI: Evaluation of candidate reference genes for expression studies in Pisum sativum under different experimental conditions. Planta. 2010, 232: 145-153. 10.1007/s00425-010-1158-1.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.