| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1 Southern Crop Protection and Food Research Centre, Agriculture and Agri-Food Canada (AAFC), 1391 Sandford St, London, ON N5V 4T3, Canada
2 Department of Biology, The University of Western Ontario, Biological and Geological Building, 1151 Richmond St, London, ON N6A 5B7, Canada
Correspondence
Aiming Wang
wanga{at}agr.gc.ca
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTS AND DISCUSSIONREFERENCES |
|---|
2- or 
–2-fold change and a Q value of 0.05), 273 (1.5 %) and 173 (0.9 %) transcripts were identified to be up- and downregulated, respectively, from 18 613 soybean cDNAs on the array. The expression levels of many transcripts encoding proteins for hormone metabolism, cell-wall biogenesis, chloroplast functions and photosynthesis were repressed at 14 days p.i. and were associated with the highest levels of viral RNA in the host cells. A number of transcripts corresponding to genes involved in defence were either downregulated or not affected at the early stages of infection, but upregulated at the late stages, indicating that the plant immune response is not activated until the late time points of infection. Such a delayed defence response may be critical for SMV to establish its systemic infection.
Present address: Department of Molecular Genetics, University of Toronto, Toronto M5S 1A8, Canada. 
The NCBI Gene Expression Omnibus (GEO) accession numbers for the microarray data from this work are GSE9824, GPL6258, GSM247941, GSM247942 and GSM247943.
A supplementary figure showing confirmation of microarray data by Northern hybridizations, as well as five supplementary tables and references cited therein, are available with the online version of this paper.
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONREFERENCES |
|---|
; Golem & Culver, 2003; Whitham et al., 2003; Pompe-Novak et al., 2005). The magnitude of physiological and phenotypic changes in the host during viral infection suggests the involvement of a large number of host genes (Golem & Culver, 2003; Whitham et al., 2006). Thus, the intimate interaction between a plant virus and its host is complicated by the systemic nature of infection and global alterations in host gene expression (Maule et al., 2002; Whitham & Wang, 2004). Dissecting the plant gene-expression network that occurs in response to virus infection should assist in a better understanding of the infection process and ultimately in the control of plant viruses. To achieve this, microarray technology has been adopted to profile global gene expression of the host (mainly model plant species) infected with several positive-sense RNA viruses (Golem & Culver, 2003; Whitham et al., 2003; Ishihara et al., 2004; Marathe et al., 2004; Pompe-Novak et al., 2005; Senthil et al., 2005; Dardick, 2007; Shimizu et al., 2007; Yang et al., 2007).
Soybean mosaic virus (SMV) is a member of the genus Potyvirus in the family Potyviridae, which is the largest plant virus family. SMV is the most prevalent viral pathogen of soybean [Glycine max (L.) Merr.] in the world. Infection by SMV usually causes yield losses of between 35 and 50 % under natural field conditions and of up to 50–100 % in severe outbreaks (Arif & Hassan, 2002; Liao et al., 2002). SMV has a single-stranded, positive-sense [ss(+)] RNA genome approximately 9600 nt in length (Jayaram et al., 1992). Whilst the virus itself has been relatively well studied, the molecular mechanisms underlying soybean responses to SMV infection remain poorly understood.
In this work, we report analysis of gene expression in SMV-infected soybean plants during the course of infection. The most pronounced changes in gene expression occurred when viral RNA accumulation reached the highest level and the leaf underwent moderate symptom development. The SMV-induced and -repressed transcripts were classified based on their putative functions and clustered according to their expression patterns. A cross-sequence comparison of the transcriptional response to SMV and to other ss(+) RNA viruses identified transcripts with similar changes in gene expression. Cumulatively, the gene-expression analysis presented in this study reveals correlations between host gene-expression changes and disease-symptom development in soybean.
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONREFERENCES |
|---|
Microarray probe labelling, hybridization, scanning and data analysis.
Soybean cDNA microarrays (18K A series) containing 18 613 soybean cDNAs of low redundancy were obtained from Dr Lila Vodkin (University of Illinois, Urbana–Champaign, IL, USA). The GEO accession number of these expressed sequence tags (ESTs) is GPL3015
[NCBI GEO]
. Two separate cDNA-labelling reactions (Cy3 and Cy5), one for the SMV-infected leaf and the other for the mock-inoculated leaf, were carried out. In total, 18 microarray slides were used for this study: six array hybridizations, including three reciprocal labelling experiments from three independent biological replications of either SMV-infected or mock-inoculated leaf, at each of the three time points. Dye swaps were performed to ensure that the results were not biased by dye effects. Total RNA (10 µg), isolated from each of the virus-infected and mock-inoculated leaf tissues at different time points, was labelled with either Cy3 or Cy5 fluorescent dye by using a CyScribe post-labelling kit (Amersham Biosciences). Synthesized probes were purified by CyScribe GFX (Amersham Biosciences) and hybridized to soybean chips following the manufacturer's protocol. Subsequent processing of the slides was essentially as described by Moy et al. (2004). The quantified data extracted from the 16-bit TIFF images from ArrayVision v. 6.0 (Imaging Research) were background-subtracted and then analysed with GeneSpring microarray-analysis software version 7.3 (Silicon Genetics), where LOWESS normalization was used to correct for any spatial or intensity-dependent biases within each array. Data analysis was essentially as described by Senthil et al. (2005). In brief, the mean normalized signal-intensity values for each transcript were calculated from six replicate hybridizations (three biological replicates and a dye-swap hybridization for each biological replicate) for each time point. log2 ratios were then calculated as described in the GeneSpring microarray-analysis instruction manual, where the mean of normalized signal-intensity values from the virus-infected samples was divided by the mean of respective values from control samples (Cy5-infected/Cy3-control or Cy3-infected/Cy5-control in a dye-swap experiment). The fold changes of differentially regulated transcripts in virus-infected samples compared with the control samples were calculated based on these ratios. To select transcripts with significant changes during virus infection, P values (P0.05) derived from ANOVA were adjusted by using the multiple testing correction of Benjamini & Hochberg (1995) with a 5 % false-discovery rate (FDR), corresponding to a Q value of 0.05 (Storey & Tibshirani, 2003). This multiple testing correction procedure is used to correct for the occurrence of false positives and to maximize the likelihood of finding significant gene sets. In addition to these statistical criteria, we searched for transcripts that showed significant up- (2.0-fold) or down- (–2.0-fold) regulation at at least one time point. Expression profiles from each time point were clustered based on their similarity in expression pattern by using a hierarchical average linkage clustering algorithm and Pearson correlation distance metric, implemented in the GeneSpring v. 7.3 software. A gene tree heat map was also built based on this analysis.
Functional categorization of transcripts.
The 5' and 3' sequences of each differentially expressed soybean unique sequence of expressed sequence tag (uniEST) were queried against the non-redundant (nr) protein database by using the BLASTX algorithm (Altschul et al., 1997). In addition, we identified the protein in the Dana Farber Cancer Institute (DFCI) Soybean Gene Index database (http://compbio.dfci.harvard.edu/tgi/cgi-bin/tgi/gimain.pl?gudb=soybean), using the highest BLASTX score to identify the function of each nr uniEST. Significantly induced and repressed soybean transcripts were grouped into functional groups by performing searches using annotated soybean transcripts against the Arabidopsis MIPS (Munich Information Center for Protein Sequences) functional classification scheme to create a custom set of functional gene categories and subcategories. Transcripts with no matches in the database were labelled as ‘unknown’. Both 
2 and Fisher's exact tests were carried out to confirm the significance on the representation of the numbers of genes in each functional category at each time point (P0.05). Statistical measurement to calculate the significant category was performed as described by Draghici et al. (2003).
Cross-sequence comparison of transcripts regulated differentially in soybean infected with SMV and in plants infected with other ss(+) RNA viruses.
The EST sequences of soybean transcripts regulated differentially in SMV-infected leaf tissues were cross-compared by BLAST-searching against the sequences of the transcripts that were significantly differentially expressed in response to infection by other ss(+) RNA viruses (Golem & Culver, 2003; Whitham et al., 2003; Ishihara et al., 2004; Marathe et al., 2004; Dardick, 2007; Yang et al., 2007). Sequence searches were performed by using the TBLASTX algorithm with default settings. The resulting BLAST output for each transcript was then parsed for high-scoring pair (HSP) and E value. Any hits with the existence of HSP100 and E10–20 were indicative of significant similarity (Rubin et al., 2000).
Northern hybridization.
Three independent Northern blots were conducted to validate the microarray data. Total RNA (10 µg) was separated and transferred to a nylon membrane. RT-PCR products derived from primers (given in Supplementary Table S1, available in JGV Online) were used as a probe. The probes were labelled with [
-32P]dCTP by using a Ready-To-Go DNA-labelling kit (Amersham Biosciences). A phosphorimager and QualityOne quantification v. 4.2 software (Bio-Rad) were used for visualization and quantification of radioactive signals.
|
RESULTS AND DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONREFERENCES |
|---|
). Moderate mosaic and vein-clearing symptoms were evident in the trifoliate leaf at approximately 14 days p.i. Symptoms of the trifoliate leaf turned to severe mosaic and mild mottling at about 21 days p.i. (Fig. 1a). Accumulation of viral RNA was detected by Northern hybridizations in the first trifoliate leaf at 7 days p.i., reached its peak level at 14 days p.i. and then decreased by 1.7-fold (58 % of the maximum viral RNA level) at 21 days p.i. (Fig. 1b).
|
Identification of soybean transcripts regulated differentially over time
Gene-expression data were analysed by one-way ANOVA to identify differentially expressed transcripts. By using stringent criteria with a 5 % FDR that corresponds to a Q value of 0.05, we identified approximately 4.8 % nr transcripts (894 of the 18 613 soybean transcripts printed on the chip) that were significantly differentially expressed (induced or repressed) in response to SMV infection at three different time points. These significantly differentially regulated transcripts were filtered further by using a 2.0-fold increase or a –2.0-fold decrease in signal intensity. At 7 days p.i., 62 nr transcripts were induced and 46 were repressed by SMV; at 14 days p.i., 154 were induced and 90 repressed; at 21 days p.i., 100 were induced and 45 repressed. This final filtering resulted in 316 transcripts (representing 273 nr transcripts) that were induced by 2.0-fold, and 181 transcripts (representing 173 nr transcripts) that were repressed by –2.0-fold, as illustrated by Venn diagrams (Fig. 2a, b). Comparison analysis between the 273 induced and the 173 repressed nr transcripts identified 168 transcripts that were significantly induced only at one time point but not significantly repressed at other time points, 68 transcripts significantly repressed only at one time point but not significantly induced at other time points and 105 transcripts significantly induced and repressed at different time points (Fig. 2c).
|
). uniESTs were assigned with a putative function only if the best match had an E value of 10–20. Based on the putative functions, the differentially expressed transcripts were classified into 11 functional categories (Fig. 3; see Supplementary Tables S2, S3 and S4, available in JGV Online). In four instances, subcategories were created to interpret the data better. 2 and Fisher's exact tests suggested that the majority of transcripts are associated with defence, chromatin regulation and cytoskeleton reorganization, protein synthesis and translation, metabolism (including cell wall-related and sulphur-assimilation transcripts) and development/storage proteins (including hormone-, chloroplast- and photosynthesis-related transcripts), comprising approximately 60 % of all differentially regulated transcripts changed significantly at a P value of 0.05 (Table 2; Fig. 3).
|
|
|
). A gene tree heat map was built by using the Pearson correlation distance metric with an average linkage clustering algorithm. Horizontal time points, representing differences in expression of transcripts, were plotted against the vertical grouping, corresponding to expression similarity. Six distinct groups of expression pattern were identified from the clustering analysis. A complete list of genes induced or repressed in different cluster groups, along with their expression levels and putative functions, is provided in Supplementary Table S2 (available in JGV Online).
|
Confirmation of microarray data by Northern hybridization
As most gene changes occurred at 14 and 21 days p.i., nine transcripts either up- or downregulated significantly at 14 or 21 days p.i. were selected for Northern hybridizations to validate the microarray data. Overall expression patterns of all nine transcripts analysed by Northern blotting were consistent with those obtained by microarray hybridizations (Fig. 5; Supplementary Fig. S1, available in JGV Online). The relative expression ratios of the transcripts analysed by Northern hybridizations were lower than those from microarray hybridizations. This could be attributed to methodological variations, for example, normalization to actin expression in Northern blotting compared with a more global normalization method used in microarray hybridizations (Taniguchi et al., 2001; Yao et al., 2004).
|
). In virus–plant interactions, Aranda et al. (1996) showed that the level of LOX1 transcript was depleted significantly in pea embryo cells undergoing pea seed-borne mosaic virus replication, and recovered at a shorter distance behind the infection front when virus replication declined. Based on the above study, LOX is probably associated with virus infection. However, past microarray studies on soybean have shown that LOX transcripts change dramatically under diverse treatments (Thibaud-Nissen et al., 2003; Moy et al., 2004). Therefore, it is also possible that the expression changes of this gene family may simply result from biotic stress.
Another interesting metabolic transcript repressed at 7 days p.i., but induced at 14 days p.i., is a transcript encoding callose synthase (GenBank accession no. AW278243). This enzyme is deposited as a layer between the plasma membrane and the cell wall (Jacobs et al., 2003; Nishimura et al., 2003; Ueki & Citovsky, 2005). Beffa et al. (1996) showed that the increased deposition of callose in and around tobacco mosaic virus (TMV)-induced lesions in β-1,3-glucanase-deficient tobacco plants resulted in decreased susceptibility to TMV. In the current study, repression of callose synthase at 7 days p.i. and upregulation of this gene at 14 days p.i. may be associated with a delayed resistance response in soybean (see Supplementary Table S3, available in JGV Online).
An intriguing subset of transcripts that were induced substantially at 14 days p.i. includes enzymes critical in the Krebs cycle, e.g. NADP-dependent malate dehydrogenase, cytosolic malate dehydrogenase and acetyl-CoA carboxylase; the oxidative pentose phosphate pathway, e.g. phosphoenolpyruvate/phosphate translocator; amino acid synthesis, such as amine oxidase, adenosine 5-phosphosulphate reductase, adenosylhomocysteinase and S-adenosylmethionine decarboxylase; sugar metabolism, e.g. sucrose synthase; and carbohydrate synthesis, e.g. β-amylase 1, mannosyloligosaccharide 1,2--mannosidase, GDP-mannose 3,5-epimerase and glutamate carboxypeptidase. Previous studies have suggested that the appearance of disease symptoms such as chlorosis on cucumber mosaic virus (CMV)-infected leaf tissues in cucumber is correlated with increased glycolysis, phosphoenolpyuruvate phosphate, respiration and starch accumulation, decreased photosynthesis and reduction in total protein synthesis (Técsi et al., 1996). The accumulation of sugar and starch in virus-infected plants is correlated with symptom development, such as altered coloration in leaf, leaf distortion, chlorosis, mosaic and stunted phenotype (von Schaewen et al., 1990; Herbers et al., 1997). Based on the above findings, it is possible that host gene-expression changes related to sugar, starch and amino acid metabolism constitute part of the host response to infection and that these changes may contribute to symptom development in soybean.
Repression of transcripts related to hormone, anti-oxidative metabolism, cell wall, sulphur assimilation, chloroplast and photosynthesis with SMV infection and symptom development
A wide range of symptoms in virus-infected plants can be linked to plant hormones, such as auxin or gibberellin (Padmanabhan et al., 2005, 2006; Culver & Padmanabhan, 2007). Interaction between the TMV replicase and auxin/indole acetic acid proteins has been shown to affect auxin-mediated pathways and contribute to disease symptoms in Arabidopsis (Padmanabhan et al., 2005, 2006). Transcripts encoding auxin-repressed protein were downregulated significantly in potato plants infected with potato virus Y at 7 days p.i. (Pompe-Novak et al., 2005). In rice plants infected with rice dwarf virus, the expression of a gibberellin-biosynthetic enzyme, ent-kaurene oxidase, was repressed, leading to the appearance of dwarf phenotype (Itoh et al., 2004; Zhu et al., 2005). In this study, five auxin-associated transcripts and one gibberellin-associated transcript were repressed in SMV-infected soybean leaf tissues at 14 or 21 days p.i. (see Supplementary Table S4, available in JGV Online). It is not clear whether these hormone-related gene changes also contribute to SMV symptom development in soybean.
Reactive or activated oxygen species have been suggested to be key mediators of local and systemic resistance responses in incompatible plant–pathogen reactions and to be involved in symptom development and pathogenesis in compatible plant–virus interactions (Sandermann, 2000; Hernández et al., 2004). In this study, the expression of two anti-oxidative metabolism-related ESTs encoding catalase (CAT; GenBank accession no. AW349008) and glutathione S-transferase (GST; GenBank no. AW472161) was suppressed significantly in the SMV-infected leaf at 14 days p.i. Another key detoxification gene, superoxide dismutase (SOD), was downregulated significantly at 7 days p.i. In contrast, two ESTs (GenBank accession nos AW471843 and AI938378) encoding peroxidase (PO) were induced significantly at 7 days p.i. (by about 1.7-fold) and 14 days p.i. (by >4-fold). Consistent with these transcriptional profiles, Zhuang et al. (1993) reported that SOD activity decreased with an increase in PO activity in the SMV-infected soybean leaf. A decrease of CAT and SOD activity and an increase of PO activity were also found in many other compatible host–virus interactions, such as in Phaseolus vulgaris infected with white clover mosaic virus (Clarke et al., 2002). Upregulation of PO transcripts and downregulation of CAT, SOD and GST transcripts may induce an oxidative stress in the early infection process. Such an anti-oxidative metabolism imbalance may be associated with the progression of SMV infection and symptom development, as suggested for the plum pox virus (PPV)–peach interaction (Hernández et al., 2004).
The other interesting changes observed in this study were the significant repression of cell wall-related genes. This gene group, containing several hundred different structural proteins and cell wall-related enzymes, is known to be a major determinant of cell morphogenesis in plants (Milioni et al., 2001; Huckelhoven, 2007). In SMV-infected leaves, 10 cell wall-related genes encoding proteins for reassembly (six genes), matrix polymers (two genes) and expansins (two genes) (see Supplementary Table S4, available in JGV Online) were downregulated significantly at 14 days p.i. Downregulation of cell wall-related transcripts during SMV infection may be associated with a reduction in cell wall cross-linking that contributes to disease-symptom development. This is consistent with the recent finding that cell wall-related transcripts were downregulated significantly in Arabidopsis infected with turnip mosaic virus (TuMV) (Yang et al., 2007) and in rice infected with rice dwarf virus (Shimizu et al., 2007). In these two studies, repression of cell wall-related transcripts was correlated with symptom development.
Chlorosis is usually associated with susceptible interactions in which virus replicates and moves throughout the plant. Chlorotic symptoms appear as yellowed areas in expanded leaves. In fully developed leaves, virus is largely confined to light-green areas, where it seems to interfere with chloroplast structure, function and/or development (Culver et al., 1991). In this study, 18 chloroplast and four photosynthesis transcripts belonging to the ‘chloroplast and photosynthesis’ functional category were downregulated in the SMV-infected leaf at 14 days p.i. (see Supplementary Table S4, available in JGV Online). Recently, Yang et al. (2007) showed that a large fraction of the metabolic genes encoding chloroplast- and photosynthesis-related proteins were repressed significantly in Arabidopsis plants infected with TuMV. Downregulation of photosynthesis-related genes is also correlated with the development of infection symptoms, such as chlorosis, stunting or mosaic, in plants (Técsi et al., 1994, 1996; Maule et al., 2002; Pompe-Novak et al., 2005; Espinoza et al., 2007).
In this study, 10 transcripts involved in sulphur assimilation and utilization were upregulated significantly at 14 days p.i., followed by downregulation at 21 days p.i. These transcripts encode enzymes such as methionine synthase, adenosine 5-phosphosulphate reductase, adenosylhomocysteinase and adenosylmethionine decarboxylase, which play a major role in sulphur uptake. Recently, Yang et al. (2007) also reported that TuMV infection suppressed expression of genes involved in the sulphur-uptake pathway, affecting plant growth and development. Taken together, these results suggest that gene-expression changes in hormone, cell-wall biogenesis, chloroplast, photosynthesis and sulphur-assimilation pathways may contribute to symptom development in SMV-infected soybean plants.
Association of upregulated transcripts in protein synthesis and disease resistance with SMV infection
Translational regulation is a critical component of the cellular response to a variety of types of stress, such as viral infection, nutrient deprivation and heat shock. In the current study, 11 differentially regulated transcripts related to protein synthesis and translation were identified. These transcripts, including two translation-elongation factors and nine ribosomal transcripts, were affected slightly at 7 days p.i., but upregulated substantially at 14 days p.i. (see Supplementary Table S3, available in JGV Online). Similar induction of ribosomal genes was also observed in Nicotiana benthamiana infected with PPV (Dardick, 2007) and in Arabidopsis infected with TuMV (Yang et al., 2007). It is not known whether the increased expression of these ribosomal proteins is a simple stress response to compensate for the host cell that may lack sufficient translation components to maintain its viability, because many such components are hijacked by the virus for its genome translation and replication.
Transcripts encoding proteins related to defence and virulence were significantly over-represented at 14 and 21 days p.i., but not at 7 days p.i. (Table 2
), suggesting a general delayed defence in the SMV-infected leaf. Of the 24 upregulated defence-related transcripts, a subset of 17 defence-related transcripts, such as pathogenesis-related (PR) protein 3 (chitinase), GST10, HSP, SOD and PO, were either downregulated or slightly affected at 7 days p.i., but substantially upregulated at 14 or 21 days p.i. (see Supplementary Table S3, available in JGV Online). These transcripts have been found to be vital in disease signalling, plant defence and stress responses (Cardinale et al., 2002; Marathe et al., 2004; Garcia-Brugger et al., 2006; Whitham et al., 2006). It is possible that, at late infection stages, the soybean plant responds to SMV infection by expressing defence-related genes, as in the case of other virus-infected plants (Whitham et al., 2003; Senthil et al., 2005; Shimizu et al., 2007). Collectively, these data suggest that there is a delayed host defence response; the immune reponse in soybean plants is not activated until the relatively late stage of infection.
Association of a common set of induced transcripts in general stress- and defence-related categories with infection by SMV and other ss(+) RNA viruses
To examine whether the diverse ss(+) RNA viruses have the ability to elicit common gene-expression changes, a cross-sequence comparison analysis was performed. The uniEST sequences of 894 significantly differentially expressed (either up- or downregulated, regardless of fold change) soybean transcripts identified in this study were BLAST-searched against all of the differentially regulated transcripts identified in potato in response to infection by three distinct fruit tree viruses, namely PPV, tomato ringspot virus (ToRSV) and Prunus necrotic ringspot virus (Dardick, 2007), and in Arabidopsis in response to infection by CMV (Ishihara et al., 2004; Marathe et al., 2004), TuMV (Yang et al., 2007), TMV (Golem & Culver, 2003) and five positive-sense plant viruses including TuMV, oilseed rape mosaic virus, turnip vein clearing virus, potato virus X and CMV (Whitham et al., 2003). As a result, 107 unique transcripts were identified (see Supplementary Table S5, available in JGV Online).
The 107 unique transcripts were classified into four major groups. Group A contains 43 soybean transcripts induced by SMV infection whose sequences match those of transcripts induced by infection with other ss(+) RNA viruses. Six of them encode ribosomal proteins (three 40S and three 60S ribosomal proteins). Group B includes 34 soybean transcripts induced during SMV infection, but these transcripts are repressed in plants infected by other ss(+) RNA viruses. For example, a transcript encoding plasma membrane ATPase (GenBank accession no. BE020672) was induced 1.3-fold by SMV infection at 21 days p.i., but was downregulated 1.5-fold by CMV infection at 6 h p.i. Three histone proteins related to chromatin regulation were induced approximately 1.2-fold by SMV infection at 7 days p.i. and were downregulated approximately 1.2-fold by ToRSV infection at 14 days p.i. Group C contains 19 SMV-induced soybean transcripts, each of which shares sequence similarity with more than one transcript induced or repressed in plants infected with other ss(+) RNA viruses. One example is a transcript induced by SMV infection at 14 days p.i. encoding heat-shock protein 70 kDa (Hsp70) that shares sequence similarity with four Arabidopsis transcripts (AT5g02500, AT3g09440, AT3g12580 and AT5g02490) induced by TuMV infection at 10 days p.i. (Yang et al., 2007). Group D consists of 11 SMV-induced soybean transcripts. In this group, more than one SMV-induced soybean transcript matches the sequence of a transcript induced or repressed by infection with other ss(+) RNA viruses. For example, the sequences of two photosystem I transcripts (GenBank accession nos AI461105 and AI495711) induced by SMV infection at 21 days p.i. match the sequence of an Arabidopsis transcript (At4g12800) induced by CMV infection at 12 h p.i. (Ishihara et al., 2004).
Among the 107 unique soybean transcripts, 62 (approx. 58 %), 14 (13 %), 12 (approx. 11 %) and 35 (33 %) transcripts share sequence similarity with transcripts regulated differentially by infection with other viruses identified by Ishihara et al. (2004), Marathe et al. (2004), Yang et al. (2007) and Dardick (2007), respectively (see Supplementary Table S5, available in JGV Online). As diverse plant RNA viruses may elicit a general stress and defence response to plant virus infection (Whitham et al., 2003), we directed our searches specifically to a common subset of such transcripts. Indeed, soybean transcripts involved in general stress and defence, such as LRR protein kinase (GenBank accession no. AW186515), pyruvate kinase (AW830175), protein phosphatase 2C (AW830157), protein kinase (AW831515), calmodulin (AI441176), peroxidase (AI496108), 2-Cys peroxiredoxin (AI443769), hypersensitive-induced response protein (AT3G01290) and universal stress protein (USP; AI735896), were identified from this analysis. It seems that, despite the differences in symptoms, hosts and viruses, there is a commonly shared stress and defence response in plants to the infections of ss(+) RNA viruses. Interestingly, all of these transcripts were upregulated by SMV infection at 14 or 21 days p.i., but not at 7 days p.i., further suggesting a delayed defence response in SMV-infected soybean plants.
In conclusion, high-throughput microarray analysis used in this study permitted us to draw some general associations between SMV infection and gene-expression changes in the SMV-infected leaf. Even though large sets of informative data were generated, our study is limited in certain ways. First, the total RNA samples derived from SMV-infected soybean leaves contained a mixture of uninfected, infected and post-infected cells, and therefore spatial and temporal information on the patterns of gene expression was not optimized. Second, certain zones in the SMV-infected leaf tissue might contain little or no viral RNA, commonly referred to as dark-green islands (Atkinson & Matthews, 1970; Moore & MacDiarmid, 2006). These zones may dilute or mask expression changes in the infected cells. Despite these limitations, this analysis suggests clearly that a number of genes are probably associated with SMV-compatible infection and symptom development, and activation of defence-like genes seems not to occur until virus accumulation reaches its highest level. Such a delayed defence response may be critical for SMV to establish its systemic infection. Further targeted functional studies of these differentially expressed transcripts, particularly those involved in defence, cell wall, sulphur assimilation, hormone, chloroplast and photosynthesis pathways, will provide new insights into SMV–soybean interactions.
|
ACKNOWLEDGEMENTS |
|---|
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTS AND DISCUSSIONREFERENCES |
|---|
Aranda, M. A., Escaler, M., Wang, D. & Maule, A. J. (1996). Induction of HSP70 and polyubiquitin expression associated with plant virus replication. Proc Natl Acad Sci U S A 93, 15289–15293.
Arif, M. & Hassan, S. (2002). Evaluation of resistance in soybean germplasm to soybean mosaic potyvirus under field conditions. J Biol Sci 2, 601–604.
Atkinson, P. H. & Matthews, R. E. F. (1970). On the origin of dark green tissue in tobacco leaves infected with tobacco mosaic virus. Virology 40, 344–356.[CrossRef][Medline]
Beffa, R. S., Hofer, R.-M., Thomas, M. & Meins, F., Jr (1996). Decreased susceptibility to virus disease of β-1,3-glucanase-deficient plants generated by antisense transformation. Plant Cell 8, 1001–1011.[Abstract]
Benjamini, Y. & Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol 57, 289–300.
Cardinale, F., Meskiene, I., Ouaked, F. & Hirt, H. (2002). Convergence and divergence of stress-induced mitogen-activated protein kinase signaling pathways at the level of two distinct mitogen-activated protein kinases. Plant Cell 14, 703–711.
Clarke, S. F., Guy, P. L., Burritt, D. J. & Jameson, P. E. (2002). Changes in the activities of antioxidant enzymes in response to virus infection and hormone treatment. Physiol Plant 114, 157–164.[CrossRef][Medline]
Culver, J. N. & Padmanabhan, M. S. (2007). Virus-induced disease: altering host physiology one interaction at a time. Annu Rev Phytopathol 45, 221–243.[CrossRef][Medline]
Culver, J. N., Lindbeck, A. G. C. & Dawson, W. O. (1991). Virus-host interactions: induction of chlorotic and necrotic responses in plants by tobamoviruses. Annu Rev Phytopathol 29, 193–217.
Dardick, C. (2007). Comparative expression profiling of Nicotiana benthamiana leaves systemically infected with three fruit tree viruses. Mol Plant Microbe Interact 20, 1004–1017.[CrossRef][Medline]
Draghici, S., Khatri, P., Martins, R. P., Ostermeier, C. & Krawetz, S. A. (2003). Global functional profiling of gene expression. Genomics 81, 98–104.[CrossRef][Medline]
Espinoza, C., Vega, A., Medina, C., Schlauch, K., Cramer, G. & Arce-Johnson, P. (2007). Gene expression associated with compatible viral diseases in grapevine cultivars. Funct Integr Genomics 7, 95–110.[CrossRef][Medline]
Garcia-Brugger, A. G., Lamotte, O., Vandelle, E., Bourque, S., Lecourieux, D., Poinssot, B., Wendehenne, D. & Pugin, A. (2006). Early signaling events induced by elicitors of plant defenses. Mol Plant Microbe Interact 19, 711–724.[CrossRef][Medline]
Golem, S. & Culver, J. N. (2003). Tobacco mosaic virus induced alterations in the gene expression profile of Arabidopsis thaliana. Mol Plant Microbe Interact 16, 681–688.[Medline]
Herbers, K., Tacke, E., Hazirezaei, M., Krause, K. P., Melzer, M., Rohde, W. & Sonnewald, U. (1997). Expression of a luteoviral movement protein in transgenic plants leads to carbohydrate accumulation and reduced photosynthetic capacity in source leaves. Plant J 12, 1045–1056.[CrossRef][Medline]
Hernández, J. A., Rubio, M., Olmos, E., Ros-Barceló, A. & Martínez-Gómez, P. (2004). Oxidative stress induced by long-term plum pox virus infection in peach (Prunus persica). Physiol Plant 122, 486–495.[CrossRef]
Huckelhoven, R. (2007). Cell wall-associated mechanisms of disease resistance and susceptibility. Annu Rev Phytopathol 45, 101–127.[CrossRef][Medline]
Ishihara, T., Sakurai, N., Sekine, K. T., Hase, S., Ikegami, M., Shibata, D. & Takahashi, H. (2004). Comparative analysis of expressed sequence tags in resistance and susceptible ecotypes of Arabidopsis thaliana infected with cucumber mosaic virus. Plant Cell Physiol 45, 470–480.
Itoh, H., Tatsumi, T., Sakamoto, T., Otomo, K., Toyomasu, T., Kitano, H., Ashikari, M., Ichihara, S. & Matsuoka, M. (2004). A rice semi-dwarf gene, Tan-Ginbozu (D35), encodes the gibberellin biosynthesis enzyme, ent-kaurene oxidase. Plant Mol Biol 54, 533–547.[CrossRef][Medline]
Jacobs, A. K., Lipka, V., Burton, R. A., Panstruga, R., Strizhov, N., Lefert, P. S. & Fincher, G. B. (2003). An Arabidopsis callose synthase, GSL5, is required for wound and papillary callose formation. Plant Cell 15, 2503–2513.
Jayaram, C. H., Hill, J. H. & Miller, W. A. (1992). Complete nucleotide sequences of two soybean mosaic virus strains differentiated by response of soybean containing the Rsv resistance gene. J Gen Virol 73, 2067–2077.
Kolomiets, M. V., Chen, H., Gladon, R. J., Braun, E. J. & Hannapel, D. J. (2000). A leaf lipoxygenase of potato induced specifically by pathogen infection. Plant Physiol 124, 1121–1130.
Liao, L., Chen, P., Buss, G. R., Yang, Q. & Tolin, S. A. (2002). Inheritance and allelism of resistance to soybean mosaic virus in Zao18 soybean from China. J Hered 93, 447–452.
Marathe, R., Guan, Z., Anandalakshmi, R., Zhao, H. & Dinesh-Kumar, S. P. (2004). Study of Arabidopsis thaliana resistome in response to cucumber mosaic virus infection using whole genome microarray. Plant Mol Biol 55, 501–520.[CrossRef][Medline]
Maule, A. J., Escaler, M. & Aranda, M. A. (2000). Programmed responses to virus replication in plants. Mol Plant Pathol 1, 9–15.[CrossRef]
Maule, A., Leh, V. & Lederer, C. (2002). The dialogue between viruses and hosts in compatible interactions. Curr Opin Plant Biol 5, 279–284.[CrossRef][Medline]
Milioni, D., Sado, P. E., Stacey, N. J., Domingo, C., Roberts, K. & McCann, M. C. (2001). Differential expression of cell-wall-related genes during the formation of tracheary elements in the Zinnia mesophyll cell system. Plant Mol Biol 47, 221–238.[CrossRef][Medline]
Moore, C. J. & MacDiarmid, R. M. (2006). Dark green islands: the phenomenon. In Natural Resistance Mechanisms of Plants to Viruses, pp. 187–209. Edited by G. Loebenstein & J. P. Carr. Dordrecht: Springer.
Moy, P., Qutob, D., Chapman, B. P., Atkinson, I. & Gijzen, M. (2004). Patterns of gene expression upon infection of soybean plants by Phytophthora sojae. Mol Plant Microbe Interact 17, 1051–1062.[Medline]
Nishimura, M. T., Stein, M., Hou, B. H., Vogel, J. P., Edwards, H. & Somerville, S. C. (2003). Loss of callose synthase results in salicylic acid dependent disease resistance. Science 301, 969–972.
Padmanabhan, M. S., Goregaoker, S. P., Golem, S., Shiferaw, H. & Culver, J. N. (2005). Interaction of the tobacco mosaic virus replicase protein with the Aux/IAA protein PAP1/IAA26 is associated with disease development. J Virol 79, 2549–2558.
Padmanabhan, M. S., Shiferaw, H. & Culver, J. N. (2006). The Tobacco mosaic virus replicase protein disrupts the localization and function of interacting Aux/IAA proteins. Mol Plant Microbe Interact 19, 864–873.[CrossRef][Medline]
Pompe-Novak, M., Gruden, K., Baebler, S., Krecic-Stres, H., Kovac, M., Jongsma, M. & Ravnikar, M. (2005). Potato virus Y induced changes in the gene expression of potato (Solanum tuberosum L.). Physiol Mol Plant Pathol 67, 237–247.[CrossRef]
Rubin, G. M., Yandell, M. D., Wortman, J. R., Gabor Miklos, G. L., Nelson, C. R., Hariharan, I. K., Fortini, M. E., Li, P. W., Apweiler, R. & other authors (2000). Comparative genomics of the eukaryotes. Science 287, 2204–2215.
Sandermann, H. (2000). Active oxygen species as mediators of plant immunity: three case studies. Biol Chem 381, 649–653.[CrossRef][Medline]
Senthil, G., Liu, H., Puram, V. G., Clark, A., Stromberg, A. & Goodin, M. M. (2005). Specific and common changes in Nicotiana benthamiana gene expression in response to infection by enveloped viruses. J Gen Virol 86, 2615–2625.
Shimizu, T., Satoh, K., Kikuchi, S. & Omura, T. (2007). The repression of cell wall- and plastid-related genes and the induction of defense-related genes in rice plants infected with rice dwarf virus. Mol Plant Microbe Interact 20, 247–254.[CrossRef][Medline]
Storey, J. D. & Tibshirani, R. (2003). Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100, 9440–9445.
Taniguchi, M., Miura, K., Iwao, H. & Yamanaka, S. (2001). Quantitative assessment of DNA microarrays – comparison with northern blot analyses. Genomics 71, 34–39.[CrossRef][Medline]
Técsi, L. I., Maule, A. J., Smith, A. M. & Leegood, R. C. (1994). Complex, localized changes in CO2 assimilation and starch content associated with the susceptible interaction between cucumber mosaic virus and a cucurbit host. Plant J 5, 837–847.[CrossRef]
Técsi, L. I., Smith, A. M., Maule, A. J. & Leegood, R. C. (1996). A spatial analysis of physiological changes associated with infection of cotyledons of marrow plants with cucumber mosaic virus. Plant Physiol 111, 975–985.[Abstract]
Thibaud-Nissen, F., Shealy, R. T., Khanna, A. & Voldkin, L. O. (2003). Clustering of microarray data reveals transcript patterns associated with somatic embryogenesis in soybean. Plant Physiol 132, 118–136.
Ueki, S. & Citovsky, V. (2005). Identification of an interactor of cadmium ion-induced glycine-rich protein involved in regulation of callose levels in plant vasculature. Proc Natl Acad Sci U S A 102, 12089–12094.
von Schaewen, A., Stitt, M., Schmidt, R., Sonnewald, U. & Willmitzer, L. (1990). Expression of a yeast-derived invertase in the cell wall of tobacco and Arabidopsis plants leads to accumulation of carbohydrate and inhibition of photosynthesis and strongly influences growth and phenotype of transgenic tobacco plants. EMBO J 9, 3033–3044.[Medline]
Whitham, S. A. & Wang, Y. (2004). Roles for host factors in plant viral pathogenicity. Curr Opin Plant Biol 7, 365–371.[CrossRef][Medline]
Whitham, S. A., Quan, S., Chang, H. S., Cooper, B., Estes, B., Zhu, T., Wang, X. & Hou, Y. M. (2003). Diverse RNA viruses elicit the expression of common sets of genes in susceptible Arabidopsis thaliana plants. Plant J 33, 271–283.[CrossRef][Medline]
Whitham, S. A., Yang, C. & Goodin, M. M. (2006). Global impact: elucidating plant responses to viral infection. Mol Plant Microbe Interact 19, 1207–1215.[CrossRef][Medline]
Yang, C., Guo, R., Jie, F., Nettleton, D., Peng, J., Carr, T., Yeakley, J. M., Fan, J. B. & Whitham, S. A. (2007). Spatial analysis of Arabidopsis thaliana gene expression in response to turnip mosaic virus infection. Mol Plant Microbe Interact 20, 358–370.[CrossRef][Medline]
Yao, B., Rakhade, S. N., Li, Q., Ahmed, S., Krauss, R., Draghici, S. & Loeb, J. A. (2004). Accuracy of cDNA microarray methods to detect small gene expression changes induced by neuregulin on breast epithelial cells. BMC Bioinformatics 5, 99[CrossRef][Medline]
Zhu, S., Gao, F., Cao, X., Chen, M., Ye, G., Wei, C. & Li, Y. (2005). The rice dwarf virus P2 protein interacts with ent-kaurene oxidases in vivo, leading to reduced biosynthesis of gibberellins and rice dwarf symptoms. Plant Physiol 139, 1935–1945.
Zhuang, B. C., Xu, B. & Liao, L. (1993). Change of superoxide dismutase, peroxidase and storage protein in soybean leaves after inoculation with soybean mosaic virus. Acta Phytopathol Sin 23, 261–265.
Received 18 October 2007;
accepted 24 December 2007.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |
1.0) in expression. The intensity of colour represents the degree of gene-expression level. A complete list of genes in each cluster is given in Supplementary Table S2 (available in JGV Online).