In silico comparison of transcript abundances during Arabidopsis thaliana and Glycine max resistance to Fusarium virguliforme

Background Sudden death syndrome (SDS) of soybean (Glycine max L. Merr.) is an economically important disease, caused by the semi-biotrophic fungus Fusarium solani f. sp. glycines, recently renamed Fusarium virguliforme (Fv). Due to the complexity and length of the soybean-Fusarium interaction, the molecular mechanisms underlying plant resistance and susceptibility to the pathogen are not fully understood. F. virguliforme has a very wide host range for the ability to cause root rot and a very narrow host range for the ability to cause a leaf scorch. Arabidopsis thaliana is a host for many types of phytopathogens including bacteria, fungi, viruses and nematodes. Deciphering the variations among transcript abundances (TAs) of functional orthologous genes of soybean and A. thaliana involved in the interaction will provide insights into plant resistance to F. viguliforme. Results In this study, we reported the analyses of microarrays measuring TA in whole plants after A. thaliana cv 'Columbia' was challenged with fungal pathogen F. virguliforme. Infection caused significant variations in TAs. The total number of increased transcripts was nearly four times more than that of decreased transcripts in abundance. A putative resistance pathway involved in responding to the pathogen infection in A. thaliana was identified and compared to that reported in soybean. Conclusion Microarray experiments allow the interrogation of tens of thousands of transcripts simultaneously and thus, the identification of plant pathways is likely to be involved in plant resistance to Fusarial pathogens. Dissection of the set functional orthologous genes between soybean and A. thaliana enabled a broad view of the functional relationships and molecular interactions among plant genes involved in F. virguliforme resistance.


Background
Transcriptional changes play a major role in many plant defense processes [1]. Investigation of alterations in transcript abundance in functional genomics has provided unique opportunities to delve into gene functions by the comparison of species, tissue and time specific transcript accumulation for thousands of genes simultaneously [2][3][4]. The transcript abundances of the annotated genes of Arabidopsis, soybean and many other crops can be evaluated in parallel using high-density microarrays of sequenced cDNAs (AGI, 2000) or oligomers [5]. Microarray experiments have enabled the detection of significant variation in mRNA abundance and improved the understanding of the molecular mechanism of partial defense responses [6][7][8][9]. The host-pathogen interaction involved in incomplete, quantitative and partial resistance of soybean roots to F. virguliforme has been intensively investigated [9][10][11]. Transcription factors, chromatin remodeling proteins and transcript stabilizing factors are likely candidates to be involved. Regulated pathways are expected to include the synthesis of phytoalexins, signal molecules, cell wall deposition and carbon (C) and nitrogen partitioning.
Several studies suggested that disease resistance genes shared the same specificity across distantly related plant species [12][13][14][15]. The specificity of response was maintained, perhaps because of balancing selection, in lineages leading to multiple plant species [16]. However, it is difficult to conclude that a unified model of host-pathogen interactions has been determined because many of the genes underlying pathogen recognition were functional orthologs rather than the closest sequence homologous in different species.
Phytoalexins, phytoanticipins and signal molecules are three major natural products involved in plant defense with common precursors [17]. Phenylalanine ammonialyase (PAL; EC 4.3.1.5) expression has been associated with resistance to fungal pathogens in many plant species [18,19]. PAL catalyzes the deamination of phenylalanine to produce trans-cinnamic acid, the first step in the phenylpropanoid pathway leading to phytoalexins, lignins or coumarins. Multiple isoforms of the pal gene were identified in plants [20]. Manipulation of PAL, the first enzyme of the phenylpropanoid pathway together with the downstream enzymes such as cinnamate 4-hydroxylase (C4H; EC 1.14.13.11), diphenol oxidase (laccase; EC 1.10.3.2) and 4-hydroxycinnamoyl CoA ligase (4CL; EC 6.2.1.12), revealed an association with resistance to viral and fungal infection [21][22][23]. Reduction of phenylpropanoid biosynthesis in tobacco via down-regulation of PAL reduced local and systemic acquired resistance to fungal or viral infection [24,25]. Phenylpropanoid derived polymers like lignin also play an important role as a physical barrier against pathogen invasion [26]. Lignin, a complex racemic aromatic heteropolymer is the second most abundant cell wall polymer (after cellulose) and provides rigidity for the cell wall and a physical barrier against pathogens [27]. Lignin is synthesized from the phenylpropanoid metabolism reactions. These series of hydroxylation and O-methylation and conversion of side-chain carboxyl to an alcohol result in the building blocks of lignin, which is initiated by deamination of phenylalanine by the enzyme PAL where hydroxycinnamic acid esters play a central role [28,29]. 4CL is responsible for the CoA esterification of p-coumaric acid. Down-regulation of isozymes of 4CLs may alter metabolite concentrations other than those involved in lignin production, with a secondary effect on growth as a consequence [30]. Laccase was the first enzyme demonstrated to be able to perform lignin polymerization in vitro. Over-expressed laccase in the roots caused cell wall lignin deposition increases in the developing xylem [31,32,23]. Therefore lignin might be involved in the disease resistance mechanism of plant cells where wall fortification occurred in response to many pathogens and would be especially useful against root rots.
In plant, SnRKs possess a catalytic domain similar to that of sucrose non-fermenting-1 (SNF1) of yeast (Saccharomyces cerevisiae) and AMP-activated protein kinase (AMPK) of animals. SnRK1, SnRK2 and SnRK3 are the key members of SnRKs family but the SnRK1 subfamily appears to share direct structural and functional homology with the SNF1/AMPK family. The SnRK1 protein complex can be further divided into 3 subunits: the α subunit which is SNF1-like protein; the β subunit that is composed of SIP1-, SIP2-and GAL83-like proteins; and the γ subunit that is SNF4-like protein based on sequence structure and expression patterns [33,34]. Homologues of SnRKs occur in all kingdoms and they appear to be highly conserved among yeasts, animals and plants suggesting that they may play very similar roles across species [35]. SNF1 modulates the phosphorylation state of a number of metabolic enzymes whilst SnRK1 regulates several enzymes involved in sugar metabolism and cell energy metabolism [36,33]. In yeast, the function of SNF1 is to coordinate about 600 genes to respond to lower cellular glucose concentrations. Hong et al. [33] identified GAL83 as mediating carbon partitioning during the plant response to herbivore Manduca sexta attack. GAL83, a β-subunit of a heterotrimeric SnRK1, showed a decrease in source leaves whereas the abundance of the catalytic α-subunit of SNF remained unaltered. The herbivore-induced changes in sink-source relations in Nicotiana attenuata was regulated by the β-subunit of SnRK1 (SNF1-related kinase) protein kinase, GAL83. GAL83 silenced plants were unable to enhance root reserves, delay senescence or prolonged flowering following herbivore attack during early stages of develop-ment. In turn, SnRK1 can be used to alter resource allocation thereby plants may be equipped to better tolerate the pest attack [37].
Partial resistance can often be subcategorized as rate reducing resistance for fungal infections of roots [38,39]. Partial resistance may result from a reduced infection frequency, an extended latent period of infection, a reduced sporulation of the pathogen, or a combination of these [40][41][42]. However, the defense pathways that are induced during partial resistance do not share the same temporal and spatial patterns of gene expression observed in complete resistance [9,10,40,42,43]. Sudden death syndrome (SDS) of soybean caused by F. virguliforme [44] results from two distinct interactions. SDS has the root infection component, where the fungus exogenously penetrates root cell walls and infect specific cells and causes root rot whilst leaf scorch component, where some of the toxins produced in the roots are translocated to leaves and cause the leaf scorch. The leaf symptoms only occur in soybean [45,46] but the root rot occurs in all legumes, most dicots and some cereals. Both root rot and leaf scorch contributes to yield losses [38,47]. Arabidopsis thaliana is a host for many types of phytopathogens. We have observed that Arabidopsis is also an excellent model plant for F. virguliforme resistance (authors unpublished data). Arabidopsis was a host for F. virguliforme and the responses to the pathogen with nicely respect to the spore concentration (Figure 1). Our date showed that this ecotype demonstrated a rate reduction resistance to the fungal pathogen. Thus, A. thaliana should be useful for studying the interactions between plant and F. virguliforme.
Iqbal et al. [9] measured changes in TAs of 192 known plant defense and biotic/abiotic stress related genes in soybean roots at five time points over a period of 10 days after F. virguliforme inoculation. The cDNA arrays used were chosen from a soybean root cDNA library [48] and a subtractive hybridization experiment [9]. The temporal and spatial response differed among soybean genotypes with different numbers of SDS resistance genes. The responses were different in the partial resistance and susceptible genotypes among genes involved in the plant defense, signal recognition and transduction and metabolic processes. For most of the responding genes and all genotypes, there was an initial decrease in TAs in the inoculated roots. However, by later stages of post inoculation, the inoculated roots of the partially susceptible cultivars failed to increase abundance of any transcripts of known defense-related genes. In the most resistant cultivar, a set of 35 genes maintained at least a two-fold higher abundance at all time points. The increase in TA in RIL23 was in contrast to that observed in Essex, susceptible parent, where most of the ESTs showed either no change or a decreased TA [9].
Here, we analyzed the changes that occurred in the abundance of transcripts corresponding to 10,560 A. thaliana expressed sequence tags (ESTs) after A. thaliana cv 'Columbia' was treated with F. virguliforme. Reverse labeled slides were used. An ortholog analysis was exploited to understand the evolutionary roles of the regulated genes based on investigation of orthologous relationship between soybean and Arabidopsis. A set of resistance pathways involved in response to the pathogen infection in A. thaliana has been proposed. The comparison of the transcriptional activity in the resistance pathways between soybean and Arabidopsis after F. virguliforme pathogenesis allows for the examination of evolution of the disease response in both species.

Transcript profile of Arabidopsis genes in response to F. virguliforme infection
Analysis of the microarray data demonstrated a significant variation within and between the slides after local and global normalization [49] and 6,109 genes correlated in both reverse labeled slides from the 10,560 EST arrayed. The position and label variations between two replications did not significantly alter the topography between slides. The fungal infestation caused 168 transcripts to increase in abundance more than 1 fold (>1 on Log2 scale). About 24 of the transcripts were increased in abundance more than 2 fold (Log2 scale). In contrast, only forty-two transcripts were observed to decrease in abun-Arabidopsis thaliana responded to F. virguliforme with respect to the spore concentration Figure 1 Arabidopsis thaliana responded to F. virguliforme with respect to the spore concentration. The plant exhibited a remarkable response to the concentration of the spore on the plant fresh weight (mean weight). The mean weights labeled with the same letter were not significantly different (P < 0.05). dance by more than 1 fold (Log2 scale) and 14 of them went down more than 2 fold (Log2 scale) following F. virguliform infection. In addition, the number of the up-regulated genes was nearly four times more than that of down-regulated genes. Therefore, A. thaliana cv Columbia responses to F. virguliforme were more similar to resistant than susceptible soybean cultivars [9].
After the Arabidopsis EST-based microarray data were converted into Affymetrix 22 K array annotation, the 6,109 genes on the arrays were subjected to pathway construction using the MapMan platform [50]http:// gabi.rzpd.de/projects/MapMan. The majority (3,541) of genes altered by the treatment (P < 0.05) could not be assigned to any of the known function (Table 1). In contrast, 2,568 genes were assembled into 14 major bins, each bin representing a set of related pathways ( Table 1). Classification of the transcript abundances changed by F. virguliforme infestation showed several interesting features. There were 571 genes distributed among five protein metabolism related bins: synthesis, activation, posttranslational modification, degradation and folding. Bin16, secondary metabolism, bin29, protein metabolism and bin30, signal transduction were subjected further investigation. A large proportion of these protein metabolism related genes was assigned to bin29.2 (synthesis) and bin29.4 (posttranslational modification). Bin29.4 (protein post-translational modification) contained many The bin number of the microarray was also denoted from 1-100 and it was only shown for the bins with significant changes in transcript abundance at P < 0.05.

genes involved in published plant defense schema [51].
The 186 genes encompassed in this bin category would allow various stages of plant response to pathogen challenge to be investigated. Of the 186 genes involved in protein posttranslational modification there were only 9 genes with significantly altered TAs. The synthesis of new proteins and the alteration of the activities of existing proteins by modification have been frequently reported to be important to plant pathogen resistance.
Further in silico analysis, the platform of the Arabidopsis Interactions Viewer [52] was used to investigate the potential molecular protein-protein interaction based on the observed transcript changes. When the 186 genes were deployed in the Arabidopsis Interactions Viewer, there were more than 745 interlogs obtained. However, only two hubs (AMP kinase and phosphatase associated protein 46 and TAP46) were identified if these 9 genes were loaded into the Arabidopsis Interactions Viewer under higher stringency.

Secondary metabolism and lignin biosynthesis pathways
The 31 genes which had been classified as potentially being involved in secondary metabolism (bin 16; Table 2; Figure 2) were prominent. Several other genes that were reported to be involved in plant responses to pathogen attack were found within the schema. They were subdivided into pathways leading to isoprenoid, phenylpropanoid and lignin biosynthesis. The 18 genes were mapped into bin 16.1.4 (isoprenoids, carotenoids) included genes involved in signaling, signaling-calcium, cell wall and cellular synthesis and 13 genes were mapped into bin16.2.1 (phenylpropanoids) and bin16.7 (waxes).
The transcriptional analysis of genes from the multiple branches of the phenylpropanoid pathway showed that 6 of the 31 genes were significantly altered in TA (Table 2). Further, TA changes inferred that synthesis of cinnamic acid in Arabidopsis was one of the early responses to F. viguliforme infection. Phytoalexins and phytoanticipins serve both structural and metabolic functions in disease resistance. Our microarray data showed that the resistance response of the phenylpropanoid pathway was different in Arabidopsis and soybean [9]. The microarray data showed that the phenylpropanoid pathway was activated during the resistance response in the Arabidopsis. In contrast to soybean [9], despite the fact that the PAL was induced in this metabolic pathway, other major genes such as cinnamic acid 4-hydroxylase (C4L), chalcone synthase (CHS) and cinnamyl-alcohol dehydrogenase (CAD) in the phenylpropanoid pathway were all suppressed ( Table 2), indicating that soybean and Arabidopsis did not share similar strategies in the specific pathway to resistance F. virguliforme. Based on the up-regulation of ESTs, total six genes in the phenylpropanoid pathway had more than 1-fold increase in abundances (Log2 scale). However, the rest of genes involved in major branches in the pathway were either down-regulated or no significant increase in TAs ( Table 2). The suppression of these branches of the phenylpropanoid pathway seemed somehow different from the soybean response to the pathogen [9]. Since transcripts of those enzymes leading to the synthesis of flavonols, terpenes and proanthocyanidins, decrease in abundances of down-stream genes in the procedure of the resistance response in the pathway but increase in PAL and cinnamoyl-CoA reductase transcripts during this same time suggesting the existence of a potential bypass to synthesize secondary metabolites in the defense response. After the 192 soybean genes have been converted into 158 functional orthologs of Arabidopsis genes, comparison of the protein-protein interaction network in both Arabidopsis and soybean demonstrated a high specificity trend in gene regulation of the two species ( Figure 3). There were 12 hubs identified in the 158 functional orthologous soybean genes whilst more than 15 network hubs were observed to mediate the resistance in the up-regulated Arabidopsis genes. Unfortunately, no resistance pathway was generated from our Arabidopsis TA changes in responses of Arabidopsis plant to F. virguliforme microarray data by the Arabidopsis Interactions Viewer program in the bin16 and no significant interlog among genes in bin16 of the data obtained.

SNF1 (sucrose non-fermenting-1)-related protein kinase (SnRK1, AKIN10) gene
When all 50 genes altered in TA from in bin30 (signaling) were subjected to Arabidopsis Interactions Viewer platform, more than 300 interlogs were identified. The genes clustered into three major categories based on both interlog confidence value and Pearson correlation coefficient.
The first group of the interlogs was mediated around MAP (mitogen-activated protein) kinases, the second group was centered on SNF1-related protein kinase 1 (SnRK1; Figure 4), and the third group was formed around ATMPK ( Table 3). The second group contained 32 genes altered in TAs and this group was composed of three putative major nodes with a total 111 hits and 443 interlog confidence values (Table 3). A putative signal transduction pathway was derived from the interaction network ( Figure 5). The results suggested that the SnRK1 gene could be important for coordinating the signal assembly of a cellular apparatus associated with the "endogenous fuel gauge" [53], since of the 31 genes in the reaction center may be differentially regulated by the endogenous AMP and sugar content in order to maintain cell defense. As shown in Figure  3, SnRK1 was centered among 31 interacting genes. The group was composed of 3 clades and most genes corresponded to the signal transduction cascade and cellular responses (data not shown).
Upstream of the SnRK1 signal cascade, the ATMPK cascade may be connected to cellular calcium (through calmodulin binding protein) and protein kinases (Figure 4, Network of protein-protein interaction of A. thaliana and soybean infested with F. virguliforme Figure 3 Network of protein-protein interaction of A. thaliana and soybean infested with F. virguliforme. The networks were generated by Cytoscape [67] and visual displays were saved as Cytoscape graphs. Interaction hubs were marked by red.   . The SnRK1 shared 33% identify with one of the mouse LKB1 paralogs (BAA76749) and there was a significant interlog predicted between the SnRK1 and LKB1. Pathogen elicitor or biotic and abiotic stresses may activate MAPK pathway through an unknown protein. On the other hand, cellular C and AMP concentrations may also trigger the AMPK to induce SnRK1via SNF4. However, the role of the AMPK of Arabidopsis and its molecular significance and molecular fragments located immediately in the upstream of SnRK1 are still unknown. The TAs of the AMPK (At5g21170) and cyclin family protein (At5g48630) were significantly increased by F. virguliform but ACC1 (acetyl-CoA carbox-ylase 1) was suppressed ( Table 3). Analysis of functional orthologs of the soybean genes in Arabidopsis indicated that significant variations existed between two species facing the same pathogen. Interestingly, SnRK1 and SNF4 were also found in the soybean resistance network ( Figure  3).

Comparison of the fungal resistance genes between Arabidopsis and soybean
168 genes with altered TAs in the Arabidopsis microarray data shared homology to the 192 soybean genes involved in resistance, signal transduction, plant defense and transport of metabolites. Those TAs from both species were functionally clustered into three major groups (data not shown) using CLUSTAL X  mum Parsimony analysis with Kimura two parameter distances. The reliabilities of each branch point were assessed by the analysis of 500 bootstrap replicates. The Maximum Parsimony (MP) analysis showed that these genes possessed homology by evolutionary descent. Maximum Parsimony phylogenetic analyses partitioned these resistance genes in Arabidopsis and soybean [9] jointly into 7 clades (denoted I-VII, data not shown). Each clade contained various members. Some clades could be further subdivided into subclasses. Clade III contained members that were the most divergent compared with those encoded by the other clades. Interestingly, there were several ESTs that were only found in one species. Majority families or subfamilies in the higher nodes tended to give low bootstrap values (data not shown).

Discussion
It is well-known that over-expression PAL in tobacco plant can effectively be used against a virulent fungal pathogen Cercospora nicotiana infection and also altered the lesion phenotype inoculated with TMV [57]. Reduction of phenylpropanoid biosynthesis in tobacco via down-regulation of PAL reduces local and systemic acquired resistance to fungal or viral infection [27,28]. The major function of PAL is to catalyze phytoalexins and phytoanticipins production and thereby reduce disease severity [58]. Constitutive over-expression of the PAL gene from the tropical pasture legume Stylosanthes humilis in tobacco plants provided resistance to Cercospora nicotianae and to pathogen Phytophthora parasitica pv. Nicotianae [59]. The key enzymes of phenylpropanoid pathway, PAL, C4L and CHS were all induced in soybean resistance to F. virguliforme [9]. The global expression analysis of the Arabidopsis plant challenged by F. virguliforme appeared to be different. The fungal pathogen affected the secondary metabolism not only in the phenylpropanoid pathway but also via pathways leading to other cellular functions (terpenes, phenolics, special N compounds metabolism). Transcripts of the PAL gene encode the key enzyme feeding all these branches. Therefore, the decreases in most key transcripts at the times when PAL and cinnamoyl-CoA reductase transcripts increases in abundance is a novel finding showing that the synthesis of PAL and cinnamoyl-CoA reductase transcripts is with negative correlation to the activity of the other genes during the pathogen stress response. Perhaps in A. thaliana, the phenylpropanoid pathway did not directly participate in the cascade of reactions elicited during the defense and probably, the biosynthesis of phytoalexins and phytoanticipins was regulated in different fashions in soybean and Arabidopsis. Transcripts of genes involved in the synthesis of lignin, flavonols, anthocyanins and proanthocyanidins (data not shown) seemed to be down-regulated during the periods of time when the genes involved in the synthesis of cinnamic acid was up-regulated. Arabidopsis lacks isoflavonoid and phytoalexins and produces mainly camalexin rather than a cocktail of isoflavonoids [60]. Therefore, in Arabidopsis the flavonoids (anthocyanins, proanthocyanidins, flavones and flavonols) may indeed not play a significant role in defense. In contrast, there was evidence that cinnamoyl alcohol dehydrogenase (CAD) was induced rapidly in Arabidopsis infected with Xanthomonas [61] and CAD is in the pathway that leads to lignin biosynthesis. In the phenylpropanoid pathway, key genes like PAL are encoded by a group of genes showing a great degree of sequence diversity. Significantly up-regulated PAL transcripts support for the hypothesis that in Arabidopsis expression of the PAL is significantly associated with pathogen infection. However, the PAL over-expression plants caused a significant reduction in growth and delayed flowering. These phenotypes may also be related to energy status or amino acid pool availability being altered by PAL.
It has been shown that the addition of sugar activates genes related to disease resistance [62]. The metabolic change in the pathogenesis could negatively affect cellular C content based on the expression of sugar-repressive Simplified diagram of putative SnRK1 pathway Figure 5 Simplified diagram of putative SnRK1 pathway. SnRK1 hub was marked by a green ring.   nuclear genes and may be involved in the regulation of these genes [63]. Therefore, there should be other mechanisms for maintaining the cellular C concentrations. Moreover, metabolic movement should provide some explanation. Our results infer that the metabolite profiles may affect the status of disease resistance in the Arabidopsis plant. The response to C stresses is probably the key role of SnRK1 in eukaryotes [64]. SnRK1 may represent a primordial protein that in plants plays an important role in all resistances. It should be noted that SnRK1 interacts with 31 proteins, and each is located up-and downstream of the gene in the interactome. Perhaps, the SnRK1 sequence has shared motifs with the rest of other proteins in the network. Interestedly, the presence of proteins in Arabidopsis related to three subunits of the SnRK1 protein might indicate an early duplication of these subunits in the process of evolution. To date, it is still not clear how, during the course of evolution of multicellular organisms, the SnRK1 system acquired the ability to be a signal system. One important implication is that SnRK1 serves as one of the hubs for the signal cascades in the plant.
Although the relative specificity of SnRK1 in plants is yet unknown, it is possible that further investigation of such role may allow us to draw new resistance pathways based on energy status. Indeed, plant SnRK1 associated with herbivore tolerance was recently identified [53]. Our results suggest that interaction between SnRK1 with its network partners may be induced by the fungal pathogen. Further studies of the gene and the interacting partners should help to uncover important specific aspects of SNF1-mediated signaling during pathogen infection.
Orthologous genes are defined by direct evolutionary descent and should play similar developmental or physiological roles. Several gene groups were identified according to the outcome of functional orthologous analysis in this study, showing same functional orthologous relationships in these two species. In our primary results, the 1 to 1 ratio of the orthologous relationship was found among those several paired groups, indicating these gene groups were descended from a common ancestor and corresponded to well-conserved functions. These 1 to 1 ortholog classes are presumed to represent conserved functions in Arabidopsis and soybean, but they shared diverse bootstrapping value. Diversification following gene duplication may have occurred but the degree is still unknown. Twenty-eight pairs of the genes from both species were taxoned together and each pair possessed unique bootstrapping value and Kimura phylogenetic distance (data not shown). Based on phylogenetic analysis, the organization of soybean phenylalanine ammonialyase (PAL) proteins were very similar to that of the Arabidopsis PAL proteins, implying that the soybean and Arabidopsis PAL proteins analyzed here were all derived from a common ancestor. The fact that these clades were formed by both Arabidopsis and soybean PAL genes gave a suggestion that these PAL genes existed before the divergence of monocots and dicots. Interestingly, fair numbers of EST sequences were not found in the other organism, even in closely related gene families that were associated with the fungal resistance. However, detailed screening phylogenetic relationship among those resistance genes is necessary for detection of the gene evolution.
Plants possess an ability of resistance to most potentially pathogenic microbes. Gene transcriptional changes are critical for many plant defense processes [1].  [66]. The spore suspension, at 5 × 10 4 spore ml-1 with sterile distilled water, was made by adding F. virguliforme spores from several F. virguliforme culture plates and was continuously stirred on a stir flask to keep a uniform suspension. The spore suspension was poured on the growth medium for the infested plants and the same volume of sterile distilled water was added to non-infested plants.

Data analysis
A visualization software, MapMan was used to perform a gene ontology where the set of Arabidopsis genes of the microarray was assigned based on the non-redundant and hierarchically categorized assignment of BINs and sub-BINs at TIGR (The Institute for Genomic Research). The ontology was derived from the Affymetrix 22 K array corresponding to similar or sub-modal biological functions [50]http://gabi.rzpd.de/projects/MapMan/. The classification of Image Annotator in the software was also used to diagram the data display. The changes were expressed relative to those in pathogen challenged roots. Transcripts that increased in abundance were denoted in blue, and transcripts were decreased in abundance were denoted in red. In the scale used for the visualized data, a 1-fold change (Log2 scale) was required to produce a visible coloration, and the scale saturates at a 3-fold (log2 scale) change. Blue and Red squares were denoted increased (positive) and decreased (negative) transcript abundances (TAs), respectively. The bar was shown on log2 scale. The Arabidopsis Interaction Viewer queries a database of 19,979 predicted and 1,499 confirmed interacting proteins. The predicted interactions (interologs) were generated by Geisler-Lee et al. [52]. Output of two interlogs was transferred to the Cytoscape software environment [67] for network visualization and modeling against each other in order to catalog all of their conserved pathways and gene interaction networks. The program was equipped with a plug-in architecture for customizing applications. The visual displays were saved as a Cytoscape graphs.
Stringent quality control measures were applied to all stages of data analysis. The Microarray data were normalized by local (local background value was subtracted from the intensity value of each spot) and global metrics. The procedures described by Pevsner [49] were followed to adjust for differences in the intensity of the two labels.
Coefficients of means and variances on the signal intensities in each channel and ratio of signals from two replicates were calculated by our C++ program (available on request), which was also used to handle the missing and extra data values. The average ratio for a signal microarray from two replicates was computed by the equation of [Ratio1st+ (1/Ratio2nd)]/2. The Student's t-test was used to determine the statistical significance for genes considered between and within Microarrays slides and the plant mean weight differences (P < 0.05).

Phylogenetic analyses
Functional sequence analysis was performed on amino acid sequences using Clustal X [55] with the default settings. Nucleotide sequences were aligned with Clustal W [68]. The results of multiple alignments were subjected to phylogenetic analysis using the algorithm of the MEGA package version 4.0 using the Maximum Parsimony analysis with Kimura two parameter distances [56]. The Maximum-Parsimony was assessed by 500 bootstrap replicates. Only nucleotide sequences were employed in phylogenetic analysis based on higher stringy consideration for the phylogenetic tree construction.