Skip to main content

Comparative analysis of microarray data in Arabidopsis transcriptome during compatible interactions with plant viruses

Abstract

Background

At the moment, there are a number of publications describing gene expression profiling in virus-infected plants. Most of the data are limited to specific host-pathogen interactions involving a given virus and a model host plant – usually Arabidopsis thaliana. Even though several summarizing attempts have been made, a general picture of gene expression changes in susceptible virus-host interactions is lacking.

Methods

To analyze transcriptome response to virus infection, we have assembled currently available microarray data on changes in gene expression levels in compatible Arabidopsis-virus interactions. We used the mean r (Pearson’s correlation coefficient) for neighboring pairs to estimate pairwise local similarity in expression in the Arabidopsis genome.

Results

Here we provide a functional classification of genes with altered expression levels. We also demonstrate that responsive genes may be grouped or clustered based on their co-expression pattern and chromosomal location.

Conclusions

In summary, we found that there is a greater variety of upregulated genes in the course of viral pathogenesis as compared to repressed genes. Distribution of the responsive genes in combined viral databases differed from that of the whole Arabidopsis genome, thus underlining a role of the specific biological processes in common mechanisms of general resistance against viruses and in physiological/cellular changes caused by infection. Using integrative platforms for the analysis of gene expression data and functional profiling, we identified overrepresented functional groups among activated and repressed genes. Each virus-host interaction is unique in terms of the genes with altered expression levels and the number of shared genes affected by all viruses is very limited. At the same time, common genes can participate in virus-, fungi- and bacteria-host interaction. According to our data, non-homologous genes that are located in close proximity to each other on the chromosomes, and whose expression profiles are modified as a result of the viral infection, occupy 12% of the genome. Among them 5% form co-expressed and co-regulated clusters.

Background

Viruses are among the most agriculturally important groups of plant pathogens, causing serious economic losses in many major crops by reducing yield and quality [1]. Although viruses have relatively simple genetic structure, the detailed mechanisms of their interaction with host plants and means by which they manipulate a plant’s physiology toward their needs and trigger antiviral responses in hosts are still not well-defined [25]. Among the most important consequences of viral pathogenesis are changes in expression of host genes that define both the fate of the virus and the host’s survival chances. If plants are capable of efficiently fighting infections by inherited genetic tools, such as resistance (R) genes that are abundant in every plant species [6], they immediately initiate general resistance pathways leading to a hypersensitive response (HR). In susceptible plants lacking R genes to a specific viral pathogen, viruses induce a variety of responses to prime and elevate their infections. These include expression changes associated with cellular processes redirected by viruses for their demands and host defensive reactions to the pathogenesis [3]. Understanding the balance and interplay between these two types of responses would bring light to poorly characterized molecular mechanisms of viral comprehensive control of host immune system and to the counteracting host signaling pathways. It will also help to explain continuous and interconnected genetic variability in viral and host populations, that is, co-evolution of plants and viruses.

At the moment, there are a number of publications describing gene expression profiling in virus-infected plants that are derived mostly from DNA microarrays. They indicate a significant impact of viral infection on a wide array of cellular processes [7]. Usually, altered functional categories include responses to biotic and abiotic stresses, changes in basal plant metabolism, protein synthesis, developmental and photosynthetic processes [710].

Most of the data are limited to specific host-pathogen interactions involving a given virus and a model host plant, which usually is Arabidopsis thaliana. In spite of several efforts to summarize general changes in plant gene expression (due to viral, bacterial and fungal infections, insect attack, other biotic and abiotic stresses) having been made [35, 11], a general picture of gene expression changes in susceptible virus-host interactions is missing. Detailed knowledge about the groups of host genes participating in and/or responsive to viral pathogenesis may lead to new assertions on how host cells are controlled by infection, which defense and stress mechanisms are deployed, and why disease symptoms or deviation from normal in the growth of a plant are developed [1, 5].

In this work, in order to analyze transcriptome response to virus infection, we have assembled currently available microarray data on changes in gene expression levels in compatible Arabidopsis-virus interactions and attempted to create a functional classification of the genes with altered expression levels. We conclude that each virus-host interaction is unique in terms of the genes with altered expression levels and the number of shared genes affected by all viruses is very limited. Importantly, we also demonstrate that responsive genes may be grouped or clustered based on their co-expression pattern and chromosomal location.

Methods

Data source, microarray data

Arabidopsis expression data were obtained from Nottingham Arabidopsis Stock Centre microarray (NASC), ArrayExpress from the European Bioinformatics Institute database and the Gene Expression Omnibus database. Additional data were retrieved from supplementary material of published papers [see Additional file 1]. The data sets were log transformed (when needed) and significant genes were selected according to P < 0.05. We selected only those genes that significantly changed their expression level in response to pathogen attack by at least two-fold. Tandem duplicates were removed from the resulting profile. The total number of collected genes across all experiments was 52488. These data represent 44 experiments with 3 different types of pathogens: virus, bacteria and fungus. Among them there were 11 viruses and the total number of genes with significantly altered expression elicited by these viruses was 16816. This number included many identical genes (with the same ID) recorded in different experiments. After subtraction of the repeating genes, a list of 7639 unique genes was obtained. The same data set was used to obtain data for bacteria 17734 (11409 unique genes) and 15426 for fungi (among them 11047 unique genes).

Data analysis

We performed a meta-analysis of all the collected data on compatible virus-host interactions and also on the whole database representing viral, bacterial and fungal interactions with the host plant. We used tools from TAIR to search GO annotations and functionally classify Arabidopsis genes. To find over-represented functional groups among activated or repressed genes during virus-host interactions we used Babelomics 4 FatiGO [12] and SAE from agriGO [13]. To visualize this data we used REViGO software [14].

Clustering pathogen related genes

The level of co-expression between two genes was defined as the Pearson’s correlation coefficient (r) of the expression level for these genes. To test for pairwise local similarity in expression in the Arabidopsis genome, the mean r of the expression profiles for neighboring pairs of genes was calculated [15].

The mean r calculated from the real data set was then compared with the mean r calculated from 1000 data sets in which the order of genes in the Arabidopsis genome was randomized. We generated the stochastic distribution using a function that generates an even distribution of stochastic numbers. The proportion of genes found in clusters and the size distribution of clusters were calculated, and the values were averaged for 1000 iterations.

Results and discussion

Broad changes in gene expression during susceptible virus-host interactions

To analyze plant response to virus infection, we have assembled currently available microarray data on changes of gene expression levels in Arabidopsis thaliana in response to infection with various plant viruses: Cabbage leaf curl virus (CalCuV), Cauliflower mosaic virus (CAMV), Cucumber mosaic virus (CMV), Lettuce mosaic virus (LMV), Plum pox virus (PPV), Turnip crinkle virus (TCV), Tobacco etch virus (TEV), Tobacco mosaic virus (TMV and TMV-Cg), Tobacco rattle virus (TRV), Turnip mosaic virus (TuMV) and Oilseed rape mosaic tobamovirus (ORMV) [7, 8, 10, 1622].

The total number of genes in the assembled experiments with significantly altered expression elicited by these viruses was 16816. Among them 8684 were upregulated and 8132 were downregulated (a threshold of at least 2-fold change in expression level). However, this number included many identical genes (with the same ID) recorded in different experiments. After subtraction of the repeated genes, a list of 7639 unique genes was obtained [see Additional file 2], which represents 23% of the whole Arabidopsis genome. These are the genes either needed for the host to defend itself against the virus or for the virus to re-arrange host cellular machinery for its own needs. More than two thirds of these genes (69%) were always upregulated and only 13% were always downregulated. A sizeable portion of the genes (17%) had differential expression in response to infection with different viruses (Figure 1). Thus, the total number of induced genes (5282) exceeds that of repressed genes (1056) more than five-fold in our reduced (unique IDs) database. Approximately 15.5% of responsive genes had previously been described as involved in plant defense.

Figure 1
figure 1

Venn diagram depicting the distribution of 7639 unique genes in response to infection with different viruses. The yellow circle represents induced genes and the green circle represents repressed genes. The radius of each circle corresponds to the number of genes in the group of induced or repressed genes.

Does the larger number of activated genes as compared to repressed genes correspond to a general trend of plant response to virus, reflecting a greater diversity of upregulated genes? In other words, based on this information, can we conclude as other authors have done [23], that there is a widespread induction of the host’s biological processes due to the virus infection? De facto, it depends on several conditions. First, individual databases available for different viruses differ extensively in the number of repressed or induced genes and combined analysis is greatly influenced by this ratio in the most comprehensive databases. Second, as mentioned above, the pool of upregulated genes is larger because of the greater variety of affected genes whereas the quantity of downregulated genes is limited. Otherwise stated, more diverse genes are upregulated during different virus infections whereas downregulated genes tend to be common regardless of the particular virus. Lastly, a more accurate illustration of the general status of gene expression changes can be derived from the analysis of their proportional representation in the sets of induced or repressed genes within each functional category, which is the subject of the following section.

Using TAIR’s functional categorization, we first assigned each gene to one of the three main gene ontologies (GO) - Biological Process, Cellular Components and Molecular Function, and next to a specific functional category (FC). It is important to emphasize three key points when relying on the GO terms in analyzing expression profiles: their generality, their obvious redundancy and their incompleteness. Redundant annotations and multiple descriptions of the same biological mechanisms represent special concern undermining an effort to address consistency in characterization of gene products. Still, the Gene Ontology project [24] currently provides the most constructive way to find functionally equivalent terms for the purpose of classifying gene product properties.

Figure 2 and Table 1 show distributions of the responsive genes in different FC with respect to the total number of genes in each given GO domain of the assembled viral database. Noticeably, some of the key functions with a large number of affected genes (upregulated- or downregulated) are in the following categories: chloroplast (21% of total genes in FC), nucleus (15%) and cytosol (15%) in GO Cellular Components; hydrolase and transferase activity (13% both), protein and DNA or RNA binding (16% and 10%, respectively) in GO Molecular Function; protein metabolism (18%) and response to stress (16%) in GO Biological Process.

Figure 2
figure 2

Distributions of the responsive genes by FC in each given GO domain of the assembled viral database.

Table 1 Distributions of genes in the three main GO domains for the whole genome and assembled viral database

It was essential to determine the extent of involvement of specific FC that represents groups of genes implicated in a particular biological mechanism in the host reaction to infection. Therefore we compared the distribution of genes assigned to different FC on the whole genome of Arabidopsis with the corresponding distribution within our database of genes that are involved in response to virus infection (Table 1), [see Additional file 3]. Presumably, the greater the share each category occupies in the virus database versus in the whole genome, the greater this FC participates in host response. We found that the percentage of genes covered by several categories, such as cell wall, cytosol, extracellular, ribosome, electron transport or energy pathways, was twice as much as the normal distribution in the whole genome, thus emphasizing the important role of these functions in host-viral interactions. A share of genes in the FC “response to stress” was also 1.7 times higher in the viral database as compared to the whole genome (Table 1) and[see Additional file 3].

Common responses to different viruses

To find common responses to different viruses, we compared patterns of gene response in individual susceptible interactions. In order to do this, we used the number of shared genes among every pair of viruses to compute a similarity matrix between them according to the formula S ij  = 2n ij /(n i  + n j ), where n i and n j are the number of genes with altered expression level belonging to the database for virus i and virus j, respectively, and n ij represents the number of genes shared between both viruses [9]. Next, we arranged the computed data in accordance with the value of Sij; the higher this value is, the more similarity that exists between the two compared virus-host interactions.

As presented in Figure 3, changes of the expression pattern in response to infection with the majority of viruses were similar to the ones associated with TMV and TMV-Cg. Responses to TEV and LMV potyviruses also showed significant similarity to each other. On the other hand, even responses to RNA versus DNA viruses can be quite similar as exemplified by CMV and CalCuV.

Figure 3
figure 3

Similarity matrix ( S ij ) reflecting changes in the Arabidopsis transcriptome in response to different Arabidopsis -virus interactions. The higher the value of the S ij , the more similarity exists between two compared virus-host interactions. Each cell represents an individual virus-host interaction S ij , color-delineated according to the level of similarity ranging from considerable (red) to weak (dark blue).

Although the number of common genes affected by all viruses is very limited, each virus-host interaction is unique in terms of which genes have altered expression levels. Among them are several pathogenesis-related (PR) genes, albeit they seem to be specifically upregulated in response to particular viruses: CalCuV, TMV and TEV (PR1), CMV, CalCuV, TMV and TRV (PR2 and PR4), LMV and TRV (PR5), CMV and TMV (PR3).

Overall, we have found only 198 genes that frequently change their expression in response to the majority of the viruses [see Additional file 4. Those include genes participating in the defense or immune pathways as well as genes of catabolic and regulation processes. One of the most frequently induced genes in all interactions is AT5G38530 of the tryptophan biosynthetic pathway. The tryptophan pathway provides precursors for the synthesis of key secondary metabolites such as auxin, indole-3-acetic acid (IAA), and other molecules that help protect plants against pathogens and herbivores [25].

Proportional representation of different functional categories in the sets of induced or repressed genes

To identify overrepresented functional groups among activated or repressed genes, we subjected genes to analysis by Babelomics 4 FatiGO [12] and SEA from agriGO [13]. FatiGO uses Fisher’s exact test for 2 × 2 contingency tables to scan for significant over-representation of GO terms in one set with respect to the other. Singular enrichment analysis (SEA) identifies enriched GO terms in a list of microarray probe sets or gene identifiers. Using both FatiGO and SEA ensures finding accurate, condensed biological data by comparing a query list to a background population from which it is derived [12]. That is, such analyses predict a role of a certain biological processes in total response to infection rather than merely calculate a number of upregulated and downregulated genes.

When implemented with the assembled data set, FatiGO and SEA identified over-represented FC and subcategories among the sets of induced or repressed genes belonging to each of the main GO domains. Overrepresented in the set of repressed genes were those involved in defense response, hormone signaling (JA, ABA), response to external stimulus, photosynthesis and bioenergetics processes (encoding photosystem I and II proteins and electron transport chain) [see Additional file 5. Downregulation of these functions is presumably due to the virus overtaking host defense-related pathways and causing physiological changes associated with the disease symptoms [22].

Overrepresented in the set of induced genes were those participating in response to abiotic stimulus, responses to organic and inorganic substances, nitrogen component metabolic processes and protein transport (Golgi vesicle transport, protein targeting, and cytoskeletal protein binding). The latter host pathways are essential for facilitating virus intracellular movement. Another example of a biological process that was found only in the upregulated gene set is chromatin organization (histone modification). Chromatin structural features and posttranslational modifications play a crucial role in the regulation of gene expression [26]. Epigenetic ‘marks’ generated by modifications of histones and DNA are spread over vast regions of chromosomes and can be altered in response to stress.

Assembling information from multiple sources, such as different microarray platforms, experimental conditions, stages of infection when samples were collected, etc. raises a question of the integrity of the combined data, since it is hardly possible to eliminate the “batch effect” from influencing final results. Even so, these disparities are not likely to change the biological truth. For instance, as presented in Figure 3, the values for TMV and TMV-Cg (a crucifer-infecting strain of TMV) are very similar to each other even though they were obtained by different authors using different platforms in totally different environments. To take into account some of the determining factors (such as infection stages at the time of analyses), when different genes may become activated and/or repressed, we also looked into combined microarray data on early and late responses to virus infection and analyzed it as much as statistically possible.

Applying agriGO tools, we found that early, non-symptomatic, phases of infection are characterized by massive induction of genes belonging to both common and stress-responsive pathways. Overrepresented in the set of activated genes were amine biosynthetic processes, aromatic amino acid family metabolic processes, photosynthetic activity and responses to biotic and abiotic stresses (Figure 4). Late stages of pathogenesis, when plants are systemically infected, are characterized by repression of the majority of the stress-responsive pathways activated at the early phases, such as response to abscisic acid stimulus, response to wounding, innate immune response, response to oxidative stress, response to auxin stimulus, callose deposition in cell wall and glycosinolate metabolic process.

Figure 4
figure 4

Overrepresented GO terms in induced and repressed gene sets in the early and late stages of infection. The ontology table displays significant GO categories in response to infection with a compatible virus as determined by Babelomics 4 FatiGO and agriGO SEA. Box color reflects the flash discovery rate (FDR) to which the given GO category belongs, as shown on the scale.

One of the main features of the late host response is repression of photosynthetic and energy pathways: photosystem I and II assembly, pentose-phosphate cycle, etc. (Figure 4). Chlorosis, or yellowing of normally green plant tissue, because of the disruption of chloroplast structure and function and a decreased amount of chlorophyll is often a direct result of these deficiencies. On the contrary, cellular respiration, catabolic processes, proteolysis and senescence are overrepresented in the induced gene set at the late stages of virus infection. Among the common characteristics of both the early and late responses is negative regulation of developmental processes. In essence, sets of host genes affected at the late stage of infection closely resemble the general picture of gene expression changes caused by viral pathogenesis.

Pathogenesis-related two- to four-gene clusters in the genome of Arabidopsis thaliana

While looking at the data derived from the analysis of publicly available microarray repositories, we noticed that genes with expression profiles modified as a result of viral infection were often (12% of genome) located in close proximity to each other on the same chromosomes. Analyzing only close proximity and differential response to the infection (repressed or activated genes) we discovered 1594 such groups of genes (Figure 5A). Among them were 5 groups consisting of 8 genes, 7 groups of 7 genes, 20 groups of 6 genes, and 35 groups of 5 genes. Assuming that the order of genes with altered expression patterns along the chromosomes is not accidental [27] but reflects their functional role, we hypothesized that these groups of neighboring genes distributed across the Arabidopsis genome may further be divided into co-regulated and co-expressed blocks of genes or clusters. Since microarray data sets on susceptible host-virus interactions were not large enough to statistically predict clusters of genes with similar expression changes, we combined them with data from microarray experiments representing bacterial-host and fungal-host interactions and then used “viral sets” as a base for filtering out only analogous genes. This way, we were able to compose groups of co-expressed genes.

Figure 5
figure 5

Clusters of non-homologous, defense-related genes in the genome of A . thaliana. (A) Chromosomal distribution of 1594 groups of neighboring genes, which responded differentially (repressed or activated) to various viral infections. (B) Distribution of pathogen-related, co-expressed clusters on the chromosomes of A. thaliana. (C) Comparison between the number of clustered genes obtained from the random distribution data set or found experimentally. The height of the bars represents the number of defense-related gene clusters of the corresponding size either revealed in the experiment (black columns) or estimated by stochastic distribution (unfilled columns).

We found 207 neighboring co-expressed genes which fall into 98 clusters under conditions of pathogenesis [see Additional file 6]. These clusters consist of groups of physically linked and functionally related genes (response to pathogen) that are co-expressed (correlation coefficient r ≥ 0.7) and possibly co-regulated but share no sequence homology. Although most of them were differentially expressed, 22 clusters were always upregulated and only 2 clusters were always downregulated. Among all identified clusters only two contained four genes, nine were composed of three genes, and eighty-six contained two genes (Figure 5B).

Differences between the stochastic distribution and the actual distribution revealed in this experiment were observed both in the number of genes included in clusters and in cluster size (Figure 5C). The number of two-gene and three-gene clusters in the experiment was almost 2.5 times and 4 times higher, respectively, than expected by chance. Four-gene clusters were obtained in the experiment only; clusters of this size were not predicted to form by chance. We found 16 overlapping clusters between our pathogen-response clusters and those predicted by Zhan et al. using microarray data representing 128 experimental conditions [28]. Apparently, genes forming these clusters are broadly co-expressed in a wide range of conditions.

To find out if there are any functional relationships between locally co-expressed genes, we used TAIR’s GO for Arabidopsis. We found 7 molecular functions, which are shared for each of 8 gene pairs as well as 13 cellular components that are common for 23 gene pairs. As revealed by the AraCyc database files from the Plant Metabolic Network, none of our clusters belongs to the same pathway [see Additional file 6].

Therefore, co-expressed neighbors do not seem to be associated with a particular GO [29]. However, it does not mean that clustered genes are not functionally related to each other. That is, in spite of belonging to different GO categories, co-expressed groups of genes are affiliated with the same function – stress response. Plants re-arrange their metabolism upon recognition of pathogen-associated molecular patterns (PAMP) so that genes of different GO categories that are involved in defense mechanisms are engaged [30].

Interestingly, one of the clusters includes three genes encoding leucine-rich repeat (LRR) family proteins: AT1G33590, AT1G33600 and AT1G33610. Two more genes that are absent in available microarray data sets, AT1G33612 and AT1G33670, also encode LRR proteins and are located in the same chromosomal region. In addition, we found three other clusters containing genes with common domain structure and functional characteristics: i. cluster with a Toll-Interleukin-Resistance (TIR) domain (AT1G72900 and AT1G72910); ii. cluster with genes encoding SAM superfamily proteins (S-adenosyl-L-methionine-dependent methyltransferases superfamily, AT4G00740, and AT4G00750); and iii. cluster with genes encoding histone superfamily proteins (AT4G40030 and AT4G40040).

Previously, we reported on the clustering of pathogen-response genes in the genome of Arabidopsis thaliana[31]. That study was based on the profiling of EST databases derived from different plant species infected with fungi, bacteria, and viruses [32]. While comparing gene clusters revealed by broad EST mining with analysis of microarray data sets specific for compatible virus-host interactions (this investigation), we found that most groups of neighboring genes determined in the former study could be included in the clusters identified in this work, providing that both up and downregulated genes derived from different experiments are counted (Figure 5A). However, if only co-expressed genes are considered (Figure 5B), overlap between these two data sets is quite low, which could possibly be explained by a unique pattern of chromosomal gene clustering characteristic for different types and/or individual pathogens.

Common genes participate in virus-, fungi- and bacteria-host interaction

As mentioned above, we tried to assemble the largest possible database of Arabidopsis genes responsive to viral infections using currently available microarray data. For comparison, we also put together genes derived from microarray experiments with bacteria and fungi. Since there is significantly more information on changes in plant gene expression due to infection with bacteria and fungi, we limited the number of genes to approximately the same number as was compiled for viruses: 17734 for bacteria (11409 unique genes) and 15426 for fungi (among them 11047 unique genes) [11, 3336].

Gene expression changes in response to all pathogens were very similar. In spite of specific interactions between host plants and each of the pathogens, nearly half of the genes associated with viral infections in susceptible hosts were also involved in response to bacterial or fungal infections (Figure 6). Next, we selected genes which were induced or repressed in all three types of interactions [see Additional file 7]. Most of these genes belong to co-expressed chromosomal regions, or clusters: there were 79 two-gene clusters and 4 three-gene clusters. This suggests that common genes participating in response to biotic stress may be co-regulated and organized in clusters. A small cluster containing non-homologous genes AT1G20100 and AT1G20110 is especially interesting since it is engaged in response to the majority of plant viruses. One of these genes encoding a RING/FYVE/PHD zinc finger superfamily protein participates in signal transduction pathways and another one is a protein of unknown function.

Figure 6
figure 6

Venn diagram depicting similarity between responses to viral, bacterial and fungal pathogens. The numbers in parentheses represent a number of unique genes responsive to a particular type of pathogen.

While determining significantly over-represented functional groups that include genes activated during virus-host interactions, we have also found groups of genes whose biological functions are common for all three pathogens (virus, bacteria, and fungus) regardless of whether there is a susceptible or resistant type of interaction. To obtain significantly over-represented GO terms from gene sets upregulated during infection with these pathogens, we used the SAE tool from agriGO [13]. In order to visualize data we used REViGO [14]. We found that some of the GO categories, such as sulfur metabolism, cellular aromatic compound metabolism, and cell wall modification were enriched in the downregulated genes under compatible virus infections. However, when we considered resistant-type interactions of bacteria and fungi with host plants, we found that the same GO categories were enriched with activated genes (Figure 7). In addition, when genes involved in general immune responses were analyzed using SAE and REViGO in both susceptible and resistant types of interactions, they were also found to be induced. Unfortunately, the limited amount of microarray data did not allow a full-scale comparison between susceptible and resistant types of interactions, which would be useful in terms of understanding the mechanisms of R gene-mediated resistance.

Figure 7
figure 7

Scatter plot of over-represented GO terms from upregulated gene sets of viral, bacterial and fungal pathogens. GO terms were obtained using the SAE tool from agriGO and data were visualized by REViGO. The axes have no intrinsic meaning; the guiding principle is that semantically similar GO terms remain close together in the plot. Color scale indicates log10 p-value (red is higher and blue is lower). Disc size is proportional to the log number of genes in the category.

Conclusions

We have assembled currently available microarray data on changes in gene expression levels in compatible Arabidopsis-virus interactions. In summary, we found that there is a greater variety of upregulated genes in the course of viral pathogenesis as compared to repressed genes. Distribution of the responsive genes in combined viral databases differed from that of the whole Arabidopsis genome, thus underlining a role of the specific FC in common mechanisms of general resistance against viruses and in physiological/cellular changes caused by infection. Using integrative platforms for the analysis of gene expression data and functional profiling, we identified overrepresented functional groups among activated and repressed genes, which provided an in-depth view of the role of certain biological processes in response to infection. Each virus-host interaction was found to be unique in terms of the genes with altered expression levels, and the number of common genes affected by all viruses was very limited. We discovered that genes with expression profiles modified as a result of viral infection were often located in close proximity to each other on the same chromosomes forming a multiple clusters, consisting of physically linked and functionally related genes. Finally, combining genes derived from microarray experiments with bacteria and fungi with a viral data set, we observed that gene expression changes in response to all pathogens were very similar and that nearly half of the genes associated with viral infections in susceptible hosts were also involved in response to bacterial or fungal infections.

Abbreviations

HR:

Hypersensitive response

CalCuV:

Cabbage leaf curl virus

CAMV:

Cauliflower mosaic virus

CMV:

Cucumber mosaic virus

LMV:

Lettuce mosaic virus

PPV:

Plum pox virus

TCV:

Turnip crinkle virus

TEV:

Tobacco etch virus

TMV and TMV-Cg:

Tobacco mosaic virus

TRV:

Tobacco rattle virus

TuMV:

Turnip mosaic virus

ORMV:

Oilseed rape mosaic tobamovirus

GO:

Gene ontology annotations

ET:

Ethylene

JA:

Jasmonic acid

ABA:

Abscisic acid.

References

  1. Hull R: Matthews’ plant virology. 2002, Academic, San Diego

    Google Scholar 

  2. Kang BC, Yeam I, Jahn MM: Genetics of plant virus resistance. Annu Rev Phytopathol. 2005, 43: 581-621. 10.1146/annurev.phyto.43.011205.141140.

    Article  PubMed  CAS  Google Scholar 

  3. Whitham SA, Quan S, Chang HS, Cooper B, Estes B, Zhu T, Wang X, Hou YM: Diverse RNA viruses elicit the expression of common sets of genes in susceptible Arabidopsis thaliana plants. Plant J. 2003, 33: 271-283. 10.1046/j.1365-313X.2003.01625.x.

    Article  PubMed  CAS  Google Scholar 

  4. Ma S, Bohnert HJ: Integration of Arabidopsis thaliana stress-related transcript profiles, promoter structures, and cell-specific expression. Genome Biol. 2007, 8: R49-10.1186/gb-2007-8-4-r49.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Whitham SA, Yang C, Goodin MM: Global impact: elucidating plant responses to viral infection. Mol Plant Microbe Interact. 2006, 19: 1207-1215. 10.1094/MPMI-19-1207.

    Article  PubMed  CAS  Google Scholar 

  6. Bakker EG, Toomajian C, Kreitman M, Bergelson J: A genome-wide survey of R gene polymorphisms in Arabidopsis. Plant Cell. 2006, 18: 1803-1818. 10.1105/tpc.106.042614.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Golem S, Culver JN: Tobacco mosaic virus induced alterations in the gene expression profile of Arabidopsis thaliana. Mol Plant Microbe Interact. 2003, 16: 681-688. 10.1094/MPMI.2003.16.8.681.

    Article  PubMed  CAS  Google Scholar 

  8. Espinoza C, Medina C, Somerville S, Arce-Johnson P: Senescence-associated genes induced during compatible viral interactions with grapevine and Arabidopsis. J Exp Bot. 2007, 58: 3197-3212. 10.1093/jxb/erm165.

    Article  PubMed  CAS  Google Scholar 

  9. Agudelo-Romero P, Carbonell P, de la Iglesia F, Carrera J, Rodrigo G, Jaramillo A, Perez-Amador MA, Elena SF: Changes in the gene expression profile of Arabidopsis thaliana after infection with Tobacco etch virus. Virol J. 2008, 5: 92-10.1186/1743-422X-5-92.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Babu M, Griffiths JS, Huang TS, Wang A: Altered gene expression changes in Arabidopsis leaf tissues and protoplasts in response to Plum pox virus infection. BMC Genomics. 2008, 9: 325-10.1186/1471-2164-9-325.

    Article  PubMed  PubMed Central  Google Scholar 

  11. De Vos M, Van Oosten VR, Van Poecke RM, Van Pelt JA, Pozo MJ, Mueller MJ, Buchala AJ, Metraux JP, Van Loon LC, Dicke M, Pieterse CM: Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant Microbe Interact. 2005, 18: 923-937. 10.1094/MPMI-18-0923.

    Article  PubMed  CAS  Google Scholar 

  12. Medina I, Carbonell J, Pulido L, Madeira SC, Goetz S, Conesa A, Tarraga J, Pascual-Montano A, Nogales-Cadenas R, Santoyo J, et al: Babelomics: an integrative platform for the analysis of transcriptomics, proteomics and genomic data with advanced functional profiling. Nucleic Acids Res. 2010, 38: W210-213. 10.1093/nar/gkq388.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: W64-70. 10.1093/nar/gkq310.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Supek F, Skunca N, Repar J, Vlahovicek K, Smuc T: Translational selection is ubiquitous in prokaryotes. PLoS Genet. 2010, 6: e1001004-10.1371/journal.pgen.1001004.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Williams EJ, Bowles DJ: Coexpression of neighboring genes in the genome of Arabidopsis thaliana. Genome Res. 2004, 14: 1060-1067. 10.1101/gr.2131104.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Ishihara T, Sakurai N, Sekine KT, Hase S, Ikegami M, Shibata D, Takahashi H: Comparative analysis of expressed sequence tags in resistant and susceptible ecotypes of Arabidopsis thaliana infected with cucumber mosaic virus. Plant Cell Physiol. 2004, 45: 470-480. 10.1093/pcp/pch057.

    Article  PubMed  CAS  Google Scholar 

  17. Ascencio-Ibanez JT, Sozzani R, Lee TJ, Chu TM, Wolfinger RD, Cella R, Hanley-Bowdoin L: Global analysis of Arabidopsis gene expression uncovers a complex array of changes impacting pathogen response and cell cycle during geminivirus infection. Plant Physiol. 2008, 148: 436-454. 10.1104/pp.108.121038.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Marathe R, Guan Z, Anandalakshmi R, Zhao H, Dinesh-Kumar SP: Study of Arabidopsis thaliana resistome in response to cucumber mosaic virus infection using whole genome microarray. Plant Mol Biol. 2004, 55: 501-520. 10.1007/s11103-004-0439-0.

    Article  PubMed  CAS  Google Scholar 

  19. Yu W, Murfett J, Schoelz JE: Differential induction of symptoms in Arabidopsis by P6 of Cauliflower mosaic virus. Mol Plant Microbe Interact. 2003, 16: 35-42. 10.1094/MPMI.2003.16.1.35.

    Article  PubMed  CAS  Google Scholar 

  20. Agudelo-Romero P, Carbonell P, Perez-Amador MA, Elena SF: Virus adaptation by manipulation of host’s gene expression. PLoS One. 2008, 3: e2397-10.1371/journal.pone.0002397.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Dempsey DA, Pathirana MS, Wobbe KK, Klessig DF: Identification of an Arabidopsis locus required for resistance to turnip crinkle virus. Plant J. 1997, 11: 301-311. 10.1046/j.1365-313X.1997.11020301.x.

    Article  PubMed  CAS  Google Scholar 

  22. Yang C, Guo R, Jie F, Nettleton D, Peng J, Carr T, Yeakley JM, Fan JB, Whitham SA: Spatial analysis of arabidopsis thaliana gene expression in response to Turnip mosaic virus infection. Mol Plant Microbe Interact. 2007, 20: 358-370. 10.1094/MPMI-20-4-0358.

    Article  PubMed  CAS  Google Scholar 

  23. Dardick C: Comparative expression profiling of Nicotiana benthamiana leaves systemically infected with three fruit tree viruses. Mol Plant Microbe Interact. 2007, 20: 1004-1017. 10.1094/MPMI-20-8-1004.

    Article  PubMed  CAS  Google Scholar 

  24. [http://www.geneontology.org/GO.doc.shtml] GOp

  25. Niyogi KK, Fink GR: Two anthranilate synthase genes in Arabidopsis: defense-related regulation of the tryptophan pathway. Plant Cell. 1992, 4: 721-733.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128: 693-705. 10.1016/j.cell.2007.02.005.

    Article  PubMed  CAS  Google Scholar 

  27. Lee JM, Sonnhammer EL: Genomic gene clustering analysis of pathways in eukaryotes. Genome Res. 2003, 13: 875-882. 10.1101/gr.737703.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Zhan S, Horrocks J, Lukens LN: Islands of co-expressed neighbouring genes in Arabidopsis thaliana suggest higher-order chromosome domains. Plant J. 2006, 45: 347-357. 10.1111/j.1365-313X.2005.02619.x.

    Article  PubMed  CAS  Google Scholar 

  29. Mentzen WI, Wurtele ES: Regulon organization of Arabidopsis. BMC Plant Biol. 2008, 8: 99-10.1186/1471-2229-8-99.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  CAS  Google Scholar 

  31. Postnikova OA MN, Boutanaev AM, Nemchinov LG: Clustering Of Pathogen-Response Genes In The Genome Of Arabidopsis Thaliana. JIPB. 2011, in press

    Google Scholar 

  32. Boutanaev AM, Postnikova OA, Nemchinov LG: Mapping of heterologous expressed sequence tags as an alternative to microarrays for study of defense responses in plants. BMC Genomics. 2009, 10: 273-10.1186/1471-2164-10-273.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Bethke G, Unthan T, Uhrig JF, Poschl Y, Gust AA, Scheel D, Lee J: Flg22 regulates the release of an ethylene response factor substrate from MAP kinase 6 in Arabidopsis thaliana via ethylene signaling. Proc Natl Acad Sci U S A. 2009, 106: 8067-8072. 10.1073/pnas.0810206106.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Ferrari S, Galletti R, Denoux C, De Lorenzo G, Ausubel FM, Dewdney J: Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol. 2007, 144: 367-379. 10.1104/pp.107.095596.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Wang L, Mitra RM, Hasselmann KD, Sato M, Lenarz-Wyatt L, Cohen JD, Katagiri F, Glazebrook J: The genetic network controlling the Arabidopsis transcriptional response to Pseudomonas syringae pv. maculicola: roles of major regulators and the phytotoxin coronatine. Mol Plant Microbe Interact. 2008, 21: 1408-1420. 10.1094/MPMI-21-11-1408.

    Article  PubMed  CAS  Google Scholar 

  36. Wang L, Tsuda K, Sato M, Cohen JD, Katagiri F, Glazebrook J: Arabidopsis CaM binding protein CBP60g contributes to MAMP-induced SA accumulation and is involved in disease resistance against Pseudomonas syringae. PLoS Pathog. 2009, 5: e1000301-10.1371/journal.ppat.1000301.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was supported by the United States Department of Agriculture, Agricultural Research Service. We thank Wesley Schonborn for critical reading of the manuscript and language editing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lev G Nemchinov.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

OP collected and analyzed the data. LN wrote the manuscript. All authors discussed the results and commented on the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Postnikova, O.A., Nemchinov, L.G. Comparative analysis of microarray data in Arabidopsis transcriptome during compatible interactions with plant viruses. Virol J 9, 101 (2012). https://doi.org/10.1186/1743-422X-9-101

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1743-422X-9-101

Keywords