BACKGROUND
currently, a prevailing view is that land plants coevolved with pathogens since their first emergence on earth, 500– <dig> million years ago  <cit> .
this coevolution created complex relationships involving various mechanisms of plant defenses and diverse strategies used by pathogens to overcome these barriers  <cit> .
a comprehensive analysis and comparison of gene repertoires associated with host-pathogen interactions are of significant interest.
a difficulty of this task, however, lies in the fact that nucleotide sequences and functional annotation of plant genomes are at different stages of completion.
therefore, a straight sequence comparison of homologous genes belonging to distinct genomes and displaying differential activity is challenging.
heterologous est mapping, allowing indirect evaluation of similarity between gene sets of various genomes, represents an alternative solution  <cit> .
using this approach, we show that gene expression changes in response to diverse biological agents may be evolutionarily conserved in higher plant species and such conservation involves a wide range of inducible genes composing defense-related regulatory networks and up-regulation of genes controlling general plant metabolic pathways.
RESULTS
profiling of pathogenesis-related genes in arabidopsis
for heterologous est mapping we used the most representative est sets of wheat , tomato , potato , and soybean  infected with various plant pathogens  as well as est from the same species  that were not infected with any microorganisms .
among pathogenic species were fungi, bacteria, and viruses .
using the blast search engine  <cit>  together with custom-developed est-mapping software, that attributes each est from the database to a gene in the annotated arabidopsis genome  <cit> , we mapped "infectious" and "healthy" est of each origin  to a. thaliana genome, which resulted in several expression profiles of arabidopsis genes corresponding to the four selected genomes.
"healthy" expression profiles were sequentially deduced from the "infectious" profiles of the same genomes and therefore served as filters to eliminate genes that do not participate in defense response.
resultant individual profiles of defense-related a. thaliana genes corresponding to the four other plant species were added together to yield a generalized profile, reflecting a global pattern of a. thaliana gene expression in response to infection.
furthermore, at the next level of selection and using the same technique of profile building, this previously obtained comprehensive profile was additionally filtered through  <dig>  est derived from a. thaliana "healthy" libraries.
eventually, a final expression profile of  <dig> genes  presumably participating in defense response of a. thaliana and exhibiting positive profile values  has been created.
we used individual arabidopsis gene expression profiles obtained in the first round of selection to evaluate similarity between defense-related gene sets of corresponding species.
this was done by locating common genes in pairs of the resultant expression profiles of arabidopsis.
generally, a number of upregulated genes in a profile depends on a quantitative representation of est libraries, therefore a number of common genes had to be normalized by the sum of "infection-related" est derived from a pair of analyzed species.
we obtained the following normalized values: s. tuberosum/s.
lycopersicum , s. tuberosum/g.
max , s. lycopersicum/g.
max , t. aestivum/s.
tuberosum , t. aestivum/s.
lycopersicum , t. aestivum/g.
max .
these numbers have been mostly identical among dicotyledonous species and at least four times lower between the monocotyledon wheat and any of the four dicots.
since arabidopsis has the most comprehensive functional annotation of any plant genome, we were able to join pathogenesis-related genes identified by expression profiling into functional groups .
these groups could further be characterized by a degree of their engagement in defense reactions as well as by the average expression level of constituent genes.
the contribution of each individual functional group to defense reaction was evaluated as follows.
first, we found what fraction of genes is attributed to the group on the scale of the profile.
next, we found what fraction of genes is attributed to the group on the scale of the whole genome.
finally, a ratio between the first value and the second one showed possible involvement of each of the groups in defense response .
in general, the ratio indicates overrepresentation in the expression profile if greater than  <dig> , or underrepresentation, if less than  <dig> .
in other words, this is a degree to which different functional groups of arabidopsis genes related to stress responses and to general plant metabolism are involved in defense reactions.
although both abiotic and biotic stress response genes contributed considerably to the expression profile, degree of such contribution varied from  <dig>   to  <dig>   in case of abiotic stress, whereas in case of biotic stress it was  <dig>   and  <dig>  , respectively.
interestingly, genes of general plant metabolism were also substantially involved in defense against pathogen challenge.
implication of photosynthesis-related genes was especially high  whereas involvement of mitochondrial genes and genes related to protein synthesis was  <dig>  and  <dig>  respectively.
we next compared average values of expression levels in groups of functionally related genes with that of the whole profile .
the average values for the genes related to photosynthesis, oxidative stress and cold were significantly higher than the average value for the whole profile while the corresponding values for the rest of the functional groups were not.
among other notable groups of genes related to defense response were protein kinases and transcription factors .
we found that among  <dig> kinases,  <dig>  were related to endomembrane system with  <dig> annotated to encode membrane receptors and  <dig> described as defense response genes .
among  <dig> tf  <dig>  were annotated as genes responsive to both abiotic and biotic stimuli.
experimental assessment of computer-generated expression profile
to support est-driven predictions, we took advantage of a previously described experimental model, which is based on the interaction between a. thaliana and its compatible pathogen, yellow strain of cucumber mosaic virus )  <cit> .
cmv-infected a. thaliana plants ecotype columbia  <dig>  developed chlorotic symptoms on inoculated leaves  <dig> days post inoculation .
symptoms became systemic  <dig> to  <dig> dpi .
virus presence in the inoculated arabidopsis plants was confirmed by reverse transcription polymerase chain reaction  with cmv-specific primers .
expression of the genes encoding tf was of special interest because of their basic role in regulating plant responses to environmental stresses.
since the entire expression profile varied from  <dig>  to  <dig>  with an average of  <dig> , we have arbitrarily selected  <dig> tf with profile values in the range of  <dig> – <dig> , which approximately represented the middle of the profile.
out of twenty studied tf, thirteen  had elevated expression levels in response to the viral infection .
figure  <dig> shows some of the tf, upregulated in response to cmv.
a complete overlap between computer data and its experimental validation could hardly be expected since expression profiles were deduced from plants infected with a variety of different pathogens whereas our model system included a single plant virus.
needless to say, these data demonstrate a high stringency of computer-generated gene profiling even when applied to a random model of host-pathogen interaction.
comparison between computer-generated profile and microarray analysis of gene expression changes
to quantitatively compare the depth of the pathogen response repertoire found by a computer-generated suit of upregulated genes with arabidopsis pathogen array data, we used the shannon index, which is frequently employed to measure diversity in categorical data.
global gene expression changes assessed by the affymetrix ath <dig> array in a. thaliana plants infected with plum pox virus   <cit>  were used as a representative example of pathogen array experiment.
the whole genome ath <dig> chip contains  <dig> genes, which allowed us to compare these data with our global profile.
the index values were  <dig>  in the case of est mapping and  <dig>  in the microarray experiment.
however, for the genes encoding transcription factors, index values calculated from est mapping were higher when compared to microarray data .
this indicates a similarity of gene coverage between whole sets of induced genes identified by est mapping and microarrays.
it also points out the advantage of est-generated profiling for uncovering low-expression level genes, such as tf.
in both cases , the total amount of common induced genes was  <dig>  which corresponded to 51% and 21%, respectively .
in the microarray experiment, a number of defense-related genes specifically responsive to ppv pathogenesis in arabidopsis plants were found to be induced.
since corresponding est, obtained from ppv-infected arabidopsis, were not used to build a computer-generated profile, it is not surprising that the same genes were not revealed by est mapping.
on the other hand, the low-expression level genes were likely underrepresented in the microarray.
discussion
microarray technology helped to accumulate an immense pool of data on gene expression changes in response to different environmental factors.
yet, a computer-generated gene profiling, although limited by availability of est collections in databases, represents a valuable alternative to microarrays, which allows efficient discovery of homologous sequences in evolutionarily different species and comparison of gene sets on the whole genome scale.
importantly, computer profiling is especially sensitive to low-abundance transcripts that are normally underrepresented in cellular mrna, such as transcription factors.
furthermore, the method does not require a significant statistical support because of a large number of available transcripts  and reliable selection of upregulated genes in expression profile.
previously, we successfully used computer-assisted gene profiling for identification of testis-specific genes in drosophila  <cit> .
thus, an expansion of data obtained by microarrays with the help of est-based gene profiling as well as comparative analysis of information collected by these two approaches would undoubtedly be beneficial.
indeed, our quantitative comparison of gene profiling data using the shannon index demonstrated that est mapping and microarrays can efficiently complement each other.
application of computer-assisted gene profiling toward study of gene expression induced by plant pathogens was the main goal of this study.
in general, our results are in a good agreement with those generated by genome chip technology.
as a whole, two independent approaches  reveal the same classes of functionally related genes with expression levels elevated as a result of defense response.
these are many groups of genes involved in general plant metabolism  along with the genes of host stress and defense responses  <cit> .
in the est-derived profile, genes annotated as defense-related and responsive to abiotic stress represented 20% and 23%, respectively, of all the stress-responsive genes in the arabidopsis genome.
other classes found included putative membrane-associated receptors and many transcription factors with unknown regulatory functions.
according to our estimate, contribution of different groups of functionally-related genes to plant defense responses varied significantly.
a share of genes associated with photosynthesis, which under the same conditions is a rather constitutive and stable process, was especially high .
this may indicate their role in basal defense responses rather than a rapid switch toward viral needs  <cit> .
other researchers also found considerable variation in the degree of involvement of different groups of genes in host responses  <cit> .
in addition, computer-generated profile of gene expression changes during pathogen attack uncovered activation of many genes normally induced in response to abiotic factors.
engagement and expression levels of cold stress-responsive genes were the largest.
even though it is hard to suggest their specific role in defense reactions, it is known that plant defense mechanisms may be temperature-dependent.
for instance, at low temperatures rna silencing-mediated antiviral defense is inhibited  <cit> .
the experimental system of compatible virus-host interaction between a. thaliana and cmv, which was used in this work to evaluate reliability of computer-generated profiling, is based on the fact that the plant is not equipped with powerful means of genetic resistance such as the presence of r genes and therefore will not develop hypersensitive response or systemic acquired resistance to the pathogen.
observed expression changes in the genes, pre-selected according to the created profiles, demonstrated high stringency of est profiling even when applied to a random model of host-pathogen interaction.
activation of transcription factors, genes of general metabolism and stress genes may correspond to the initial, basal response to the pathogen, which is eventually overtaken by infection.
profiles of defense-related genes, obtained by mapping of heterologous est, represent putative arabidopsis homologs of the corresponding species.
comparison of these profiles in pairs and locating common genes allowed indirect estimation of the similarity between defense-related gene sets of different plant species based on arabidopsis homologs.
this similarity was high among three dicotyledonous species and rather low in monocotyledonous t. aestivum vs three dicots.
this suggests that the repertoire of genes participating in defense reactions in dicots and monocots, although similar to some degree, has nevertheless considerably diverged.
it is, apparently, generally conserved in dicots.
since arabidopsis  and glycine max  are in the rosids and both lycopersicum esculentum and solanum tuberosum  are in the asterids, it appears that evolutionary conservation of defense responses in these two groups may be traced to their common ancestor as far as  <dig> million years ago  <cit> .
differences in defense mechanisms between dicots and monocots could in theory date back to their split ~ <dig> million years ago  <cit> .
in summary, we have demonstrated that computer-assisted gene profiling based on heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens.
unlike microarrays, it permitted indirect comparison of the complete sets of functionally related genes in different species on the whole genome scale.
the method allows effective identification of tissue-specific and organ-specific genes, genes expressed at a particular developmental stage, and genes responsive to internal and external stimuli.
moreover, est profiling permits to further narrow down these groups to individual biological processes, cellular compartments, or molecular functions.
for instance, in this work, est mining identified a large group of defense-related genes encoding calcium-binding proteins.
since calcium is an important second messenger, playing an essential role in plant defence responses  <cit> , these data may eventually be applied to identify a role of ca2+-binding proteins in specific plant-pathogen interactions.
for the quantitative representation of expression level in profile, only est derived from non-normalized cdna libraries can be used, even though low-level expression genes in such libraries are underrepresented.
on the contrary, in normalized libraries where specific genes are better represented, their est cannot be used for quantitative analysis of expression level.
presently, est mining can easily be adopted to reveal specific genes in arabidopsis, drosophila, and mouse genomes and, to a lesser degree, in the human genome.
databases such as genbank are constantly supplied with new est and annotated genome sequences.
soon, vast datasets will be available from deep rna sequencing on next generation platforms.
because of this, we believe that the approach conveyed in this study may be successfully applied to further increase our knowledge and understanding of transcriptom dynamics in general and more specifically, its role in host-pathogen interactions and plant defense mechanisms.
CONCLUSIONS
we demonstrated that heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens.
upregulated genes identified in this study substantially overlap with those previously obtained by microarrays.
as a whole, two independent approaches  reveal the same classes of functionally related genes with expression levels elevated as a result of defense response.
these include groups of genes involved in general plant metabolism  alongside genes of host stress and defense responses.

 currently, a prevailing view is that land plants coevolved with pathogens since their first emergence on earth, 500– million years ago. this coevolution created complex relationships involving various mechanisms of plant defenses and diverse strategies used by pathogens to overcome these barriers. a comprehensive analysis and comparison of gene repertoires associated with host-pathogen interactions are of significant interest. a difficulty of this task, however, lies in the fact that nucleotide sequences and functional annotation of plant genomes are at different stages of completion. therefore, a straight sequence comparison of homologous genes belonging to distinct genomes and displaying differential activity is challenging. heterologous est mapping, allowing indirect evaluation of similarity between gene sets of various genomes, represents an alternative solution. using this approach, we show that gene expression changes in response to diverse biological agents may be evolutionarily conserved in higher plant species and such conservation involves a wide range of inducible genes composing defense-related regulatory networks and up-regulation of genes controlling general plant metabolic pathways.  profiling of pathogenesis-related genes in arabidopsis for heterologous est mapping we used the most representative est sets of wheat, tomato, potato, and soybean infected with various plant pathogens as well as est from the same species that were not infected with any microorganisms. among pathogenic species were fungi, bacteria, and viruses. using the blast search engine together with custom-developed est-mapping software, that attributes each est from the database to a gene in the annotated arabidopsis genome, we mapped "infectious" and "healthy" est of each origin to a. thaliana genome, which resulted in several expression profiles of arabidopsis genes corresponding to the four selected genomes. "healthy" expression profiles were sequentially deduced from the "infectious" profiles of the same genomes and therefore served as filters to eliminate genes that do not participate in defense response. resultant individual profiles of defense-related a. thaliana genes corresponding to the four other plant species were added together to yield a generalized profile, reflecting a global pattern of a. thaliana gene expression in response to infection. furthermore, at the next level of selection and using the same technique of profile building, this previously obtained comprehensive profile was additionally filtered through est derived from a. thaliana "healthy" libraries. eventually, a final expression profile of genes presumably participating in defense response of a. thaliana and exhibiting positive profile values has been created. we used individual arabidopsis gene expression profiles obtained in the first round of selection to evaluate similarity between defense-related gene sets of corresponding species. this was done by locating common genes in pairs of the resultant expression profiles of arabidopsis. generally, a number of upregulated genes in a profile depends on a quantitative representation of est libraries, therefore a number of common genes had to be normalized by the sum of "infection-related" est derived from a pair of analyzed species. we obtained the following normalized values: s. tuberosum/s. lycopersicum, s. tuberosum/g. max, s. lycopersicum/g. max, t. aestivum/s. tuberosum, t. aestivum/s. lycopersicum, t. aestivum/g. max. these numbers have been mostly identical among dicotyledonous species and at least four times lower between the monocotyledon wheat and any of the four dicots. since arabidopsis has the most comprehensive functional annotation of any plant genome, we were able to join pathogenesis-related genes identified by expression profiling into functional groups. these groups could further be characterized by a degree of their engagement in defense reactions as well as by the average expression level of constituent genes. the contribution of each individual functional group to defense reaction was evaluated as follows. first, we found what fraction of genes is attributed to the group on the scale of the profile. next, we found what fraction of genes is attributed to the group on the scale of the whole genome. finally, a ratio between the first value and the second one showed possible involvement of each of the groups in defense response. in general, the ratio indicates overrepresentation in the expression profile if greater than, or underrepresentation, if less than. in other words, this is a degree to which different functional groups of arabidopsis genes related to stress responses and to general plant metabolism are involved in defense reactions. although both abiotic and biotic stress response genes contributed considerably to the expression profile, degree of such contribution varied from to in case of abiotic stress, whereas in case of biotic stress it was and, respectively. interestingly, genes of general plant metabolism were also substantially involved in defense against pathogen challenge. implication of photosynthesis-related genes was especially high whereas involvement of mitochondrial genes and genes related to protein synthesis was and respectively. we next compared average values of expression levels in groups of functionally related genes with that of the whole profile. the average values for the genes related to photosynthesis, oxidative stress and cold were significantly higher than the average value for the whole profile while the corresponding values for the rest of the functional groups were not. among other notable groups of genes related to defense response were protein kinases and transcription factors. we found that among kinases, were related to endomembrane system with annotated to encode membrane receptors and described as defense response genes. among tf were annotated as genes responsive to both abiotic and biotic stimuli. experimental assessment of computer-generated expression profile to support est-driven predictions, we took advantage of a previously described experimental model, which is based on the interaction between a. thaliana and its compatible pathogen, yellow strain of cucumber mosaic virus ). cmv-infected a. thaliana plants ecotype columbia developed chlorotic symptoms on inoculated leaves days post inoculation. symptoms became systemic to dpi. virus presence in the inoculated arabidopsis plants was confirmed by reverse transcription polymerase chain reaction with cmv-specific primers. expression of the genes encoding tf was of special interest because of their basic role in regulating plant responses to environmental stresses. since the entire expression profile varied from to with an average of, we have arbitrarily selected tf with profile values in the range of –, which approximately represented the middle of the profile. out of twenty studied tf, thirteen had elevated expression levels in response to the viral infection. figure shows some of the tf, upregulated in response to cmv. a complete overlap between computer data and its experimental validation could hardly be expected since expression profiles were deduced from plants infected with a variety of different pathogens whereas our model system included a single plant virus. needless to say, these data demonstrate a high stringency of computer-generated gene profiling even when applied to a random model of host-pathogen interaction. comparison between computer-generated profile and microarray analysis of gene expression changes to quantitatively compare the depth of the pathogen response repertoire found by a computer-generated suit of upregulated genes with arabidopsis pathogen array data, we used the shannon index, which is frequently employed to measure diversity in categorical data. global gene expression changes assessed by the affymetrix ath array in a. thaliana plants infected with plum pox virus were used as a representative example of pathogen array experiment. the whole genome ath chip contains genes, which allowed us to compare these data with our global profile. the index values were in the case of est mapping and in the microarray experiment. however, for the genes encoding transcription factors, index values calculated from est mapping were higher when compared to microarray data. this indicates a similarity of gene coverage between whole sets of induced genes identified by est mapping and microarrays. it also points out the advantage of est-generated profiling for uncovering low-expression level genes, such as tf. in both cases, the total amount of common induced genes was which corresponded to 51% and 21%, respectively. in the microarray experiment, a number of defense-related genes specifically responsive to ppv pathogenesis in arabidopsis plants were found to be induced. since corresponding est, obtained from ppv-infected arabidopsis, were not used to build a computer-generated profile, it is not surprising that the same genes were not revealed by est mapping. on the other hand, the low-expression level genes were likely underrepresented in the microarray. discussion microarray technology helped to accumulate an immense pool of data on gene expression changes in response to different environmental factors. yet, a computer-generated gene profiling, although limited by availability of est collections in databases, represents a valuable alternative to microarrays, which allows efficient discovery of homologous sequences in evolutionarily different species and comparison of gene sets on the whole genome scale. importantly, computer profiling is especially sensitive to low-abundance transcripts that are normally underrepresented in cellular mrna, such as transcription factors. furthermore, the method does not require a significant statistical support because of a large number of available transcripts and reliable selection of upregulated genes in expression profile. previously, we successfully used computer-assisted gene profiling for identification of testis-specific genes in drosophila. thus, an expansion of data obtained by microarrays with the help of est-based gene profiling as well as comparative analysis of information collected by these two approaches would undoubtedly be beneficial. indeed, our quantitative comparison of gene profiling data using the shannon index demonstrated that est mapping and microarrays can efficiently complement each other. application of computer-assisted gene profiling toward study of gene expression induced by plant pathogens was the main goal of this study. in general, our results are in a good agreement with those generated by genome chip technology. as a whole, two independent approaches reveal the same classes of functionally related genes with expression levels elevated as a result of defense response. these are many groups of genes involved in general plant metabolism along with the genes of host stress and defense responses. in the est-derived profile, genes annotated as defense-related and responsive to abiotic stress represented 20% and 23%, respectively, of all the stress-responsive genes in the arabidopsis genome. other classes found included putative membrane-associated receptors and many transcription factors with unknown regulatory functions. according to our estimate, contribution of different groups of functionally-related genes to plant defense responses varied significantly. a share of genes associated with photosynthesis, which under the same conditions is a rather constitutive and stable process, was especially high. this may indicate their role in basal defense responses rather than a rapid switch toward viral needs. other researchers also found considerable variation in the degree of involvement of different groups of genes in host responses. in addition, computer-generated profile of gene expression changes during pathogen attack uncovered activation of many genes normally induced in response to abiotic factors. engagement and expression levels of cold stress-responsive genes were the largest. even though it is hard to suggest their specific role in defense reactions, it is known that plant defense mechanisms may be temperature-dependent. for instance, at low temperatures rna silencing-mediated antiviral defense is inhibited. the experimental system of compatible virus-host interaction between a. thaliana and cmv, which was used in this work to evaluate reliability of computer-generated profiling, is based on the fact that the plant is not equipped with powerful means of genetic resistance such as the presence of r genes and therefore will not develop hypersensitive response or systemic acquired resistance to the pathogen. observed expression changes in the genes, pre-selected according to the created profiles, demonstrated high stringency of est profiling even when applied to a random model of host-pathogen interaction. activation of transcription factors, genes of general metabolism and stress genes may correspond to the initial, basal response to the pathogen, which is eventually overtaken by infection. profiles of defense-related genes, obtained by mapping of heterologous est, represent putative arabidopsis homologs of the corresponding species. comparison of these profiles in pairs and locating common genes allowed indirect estimation of the similarity between defense-related gene sets of different plant species based on arabidopsis homologs. this similarity was high among three dicotyledonous species and rather low in monocotyledonous t. aestivum vs three dicots. this suggests that the repertoire of genes participating in defense reactions in dicots and monocots, although similar to some degree, has nevertheless considerably diverged. it is, apparently, generally conserved in dicots. since arabidopsis and glycine max are in the rosids and both lycopersicum esculentum and solanum tuberosum are in the asterids, it appears that evolutionary conservation of defense responses in these two groups may be traced to their common ancestor as far as million years ago. differences in defense mechanisms between dicots and monocots could in theory date back to their split ~ million years ago. in summary, we have demonstrated that computer-assisted gene profiling based on heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens. unlike microarrays, it permitted indirect comparison of the complete sets of functionally related genes in different species on the whole genome scale. the method allows effective identification of tissue-specific and organ-specific genes, genes expressed at a particular developmental stage, and genes responsive to internal and external stimuli. moreover, est profiling permits to further narrow down these groups to individual biological processes, cellular compartments, or molecular functions. for instance, in this work, est mining identified a large group of defense-related genes encoding calcium-binding proteins. since calcium is an important second messenger, playing an essential role in plant defence responses, these data may eventually be applied to identify a role of ca2+-binding proteins in specific plant-pathogen interactions. for the quantitative representation of expression level in profile, only est derived from non-normalized cdna libraries can be used, even though low-level expression genes in such libraries are underrepresented. on the contrary, in normalized libraries where specific genes are better represented, their est cannot be used for quantitative analysis of expression level. presently, est mining can easily be adopted to reveal specific genes in arabidopsis, drosophila, and mouse genomes and, to a lesser degree, in the human genome. databases such as genbank are constantly supplied with new est and annotated genome sequences. soon, vast datasets will be available from deep rna sequencing on next generation platforms. because of this, we believe that the approach conveyed in this study may be successfully applied to further increase our knowledge and understanding of transcriptom dynamics in general and more specifically, its role in host-pathogen interactions and plant defense mechanisms.  we demonstrated that heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens. upregulated genes identified in this study substantially overlap with those previously obtained by microarrays. as a whole, two independent approaches reveal the same classes of functionally related genes with expression levels elevated as a result of defense response. these include groups of genes involved in general plant metabolism alongside genes of host stress and defense responses. 
est-driven prediction identified genes which, according to the origin of est sets, were associated with defense responses in the reference genome. we demonstrated that heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens. profiles of defense-related genes, obtained by mapping of heterologous est, represent putative arabidopsis homologs of the corresponding species. yet, computer- generated gene profiling using expressed sequence tags represents a valuable alternative to microarrays, which allows efficient discovery of homologous sequences in evolutionarily different species and comparison of gene sets on the whole genome scale. microarray technology helped to accumulate an immense pool of data on gene expression changes in response to different environmental factors. in this study, we used publicly available est database derived from different plant species infected with a variety of pathogens, to generate an expression profile of homologous genes involved in defense response of a model organism, arabidopsis thaliana. comparison of these profiles in pairs and locating common genes allowed estimating similarity between defense-related gene sets of different plant species. their expression levels were examined by real-time polymerase chain reaction during infection with yellow strain of cucumber mosaic virus, a compatible virus systemically infecting arabidopsis.
 microarray technology helped to accumulate an immense pool of data on gene expression changes in response to different environmental factors. yet, computer- generated gene profiling using expressed sequence tags represents a valuable alternative to microarrays, which allows efficient discovery of homologous sequences in evolutionarily different species and comparison of gene sets on the whole genome scale. in this study, we used publicly available est database derived from different plant species infected with a variety of pathogens, to generate an expression profile of homologous genes involved in defense response of a model organism, arabidopsis thaliana.  est-driven prediction identified genes which, according to the origin of est sets, were associated with defense responses in the reference genome. profiles of defense-related genes, obtained by mapping of heterologous est, represent putative arabidopsis homologs of the corresponding species. comparison of these profiles in pairs and locating common genes allowed estimating similarity between defense-related gene sets of different plant species. to experimentally support computer data, we arbitrarily selected a number of transcription factor genes detected by est mapping. their expression levels were examined by real-time polymerase chain reaction during infection with yellow strain of cucumber mosaic virus, a compatible virus systemically infecting arabidopsis. we observed that 65% of the designated tf were upregulated in accordance with the est-generated profile.  we demonstrated that heterologous est mapping may be efficiently used to reveal genes involved in host defense responses to pathogens. upregulated genes identified in this study substantially overlap with those previously obtained by microarrays. 