<title>Background</title> <p>Systemic virus infection in plants relies on complex molecular interactions between the invading virus and the host plant [##REF##12179959##1##, ####REF##17073303##2##, ##UREF##0##3##, ##REF##17417941##4####17417941##4##]. Through such interactions, the virus recruits diverse host gene products (host factors) and metabolites for the translation and replication of its genome, which disrupts the normal biological processes of the host [##REF##17417941##4##, ####REF##15231257##5##, ##UREF##1##6####1##6##]. On the other hand, the expression of resistance genes in the infected host is also triggered by specific molecular interactions to combat the infection. Thus, virus infection induces the activation and suppression of global gene expressions in the host [##REF##12906112##7##, ####REF##12535341##8##, ##REF##15111722##9##, ##REF##15604696##10##, ##REF##16099921##11##, ##UREF##2##12##, ##REF##17722703##13##, ##REF##16775684##14##, ##REF##17378427##15##, ##REF##17427806##16####17427806##16##]. These gene expression changes are the molecular basis of general stress and defence-like responses, viral pathogenesis and host symptom development.</p> <p><italic>Plum pox virus </italic>(PPV) is a member of potyviruses, and is considered the most devastating viral pathogen of many stone-fruit species such as <italic>Prunus persica</italic>, <italic>P. persica var. nucipersica</italic>, and <italic>P. domestica </italic>[##UREF##1##6##]. Despite numerous extensive screenings for genetic resistance against the virus, no natural resistant germplams are currently available for conventional breeding programs. Very little is known about the effects of PPV on host cells at the molecular level. A better understanding of molecular PPV-plant interactions would help identify key host factors required for PPV infection and assist in the development of novel recessive resistance to PPV [##UREF##1##6##].</p> <p>To date, a major problem in dissecting molecular virus-host interactions is the progressive and asynchronous nature of the infection in leaves or in whole plants [##REF##12179959##1##,##REF##17427806##16##, ####UREF##3##17##, ##REF##17480183##18####17480183##18##]. In systemically infected plants, some cells may be infected earlier, while others may be infected later. Some parts may have higher virus titres while others may have fewer viruses or even virus-free. In the whole plant, it may take from days to weeks for the invading virus to spread from the infection site and develop systemic symptoms in the intact plant. Thus, the infected tissues contain cells at various infection stages or with viral concentrations in a wide range. As a result, temporal and spatial expression changes of host genes are very difficult to investigate. In the absence of a synchronous virus infection system and a global gene expression profiling technology, Maule and colleagues elegantly assessed the expression of selected host genes in individual cells across an advancing infection front using a cell biological approach (immunocytochemistry and <italic>in situ </italic>hybridization) (see review by Maule et al.) [##UREF##3##17##]. To date, microarray technology has been adapted to profile genome-wide gene expression changes in plants in response to infections of a number of positive sense RNA viruses [##REF##12906112##7##, ####REF##12535341##8##, ##REF##15111722##9##, ##REF##15604696##10####15604696##10##,##UREF##2##12##, ####REF##17722703##13##, ##REF##16775684##14##, ##REF##17378427##15##, ##REF##17427806##16####17427806##16##]. More recently, Yang <italic>et al</italic>. [##REF##17427806##16##] creatively dissected infection foci away from the non-infected tissue into four distinct zones using an infectious clone of <italic>Turnip mosaic virus </italic>(TuMV) tagged with green fluorescence protein (GFP) (TuMV-GFP) and profiled global gene expression changes in each of these zones with microarray, making it possible to spatially and temporally correlate virus infection and accumulation with host gene expression. Collectively these works largely advance our understanding of virus-host interactions.</p> <p>Considering the availability of whole genome microarrays, and increasing use of genetic and biochemical approaches in archetypal plants like <italic>Arabidopsis thaliana </italic>to explore the molecular interactions between viruses and their host plants [##UREF##4##19##], we used an <italic>Arabidopsis </italic>ATH1 Affymetrix array and analyzed global gene expression changes in <italic>Arabidopsis </italic>leaves infected with PPV at 17 days post inoculation (dpi). Further, we generated <italic>Arabidopsis </italic>synchronized single-cell protoplasts, transfected with a PPV infectious clone and assayed global gene expression changes in the protoplasts at 3, 6 and 12 hr post transfection (hpt). The aim of this study was to investigate general gene expression changes in infected leaf tissues and their association with PPV infection and symptom development, and to identify putative host gene factors differentially regulated upon PPV infection, particularly those at early time points involved in defence and viral pathogenesis.</p>
<title>Methods</title> <title>Plant growth and virus inoculation</title> <p><italic>Arabidopsis thaliana </italic>accession Col-0 was grown in a growth chamber at 80% RH with a day/night regime of 16 h light at 22°C followed by 8 h dark at 18°C. Prior to the preparation of viral inoculum, leaf tissues collected from <italic>Prunus glandulosa </italic>infected with a Canadian PPV-D isolate were tested positive by ELISA and by RT-PCR. The inoculum was then prepared by homogenizing approximately 1 g of <italic>P. glandulosa </italic>leaf tissues infected with PPV-D in 10 ml inoculation buffer (PO<sub>4 </sub>buffer, pH 7.4 containing 2% PVP-40 (w/v) and 0.2% diethyl dithiocarbonic acid). Viral inoculum was mechanically rub-inoculated onto <italic>Arabidopsis </italic>leaves at 5–6 leaf stage. Control plants were mock inoculated with extracts from healthy <italic>P. glandulosa </italic>leaf tissues ground in inoculation buffer. All the mock- and PPV-infected <italic>Arabidopsis </italic>plants used in the microarray experiment were randomized and grown in different locations within a greenhouse under the same growth conditions. Fifty fully expanded rosette leaves above the inoculated leaves were harvested from 25 infected or mock inoculated control plants at 17 dpi. The leaves were immediately frozen with liquid nitrogen and stored at -80°C prior to RNA extraction. Six independent hybridizations were performed using total RNA isolated from three independent biological replicates of the virus-infected or mock-inoculated control samples.</p> <title>PPV infectious cDNA clone and its deletion mutant</title> <p>The PPV infectious cDNA clone pPPV-SK68 derived from an M strain was a kind gift from Drs. Laszlo Palkovics and Ervin Balazs, Agricultural Biotechnology Center, Hungary. A control plasmid, pPPV-SK68Δ, was generated by deleting 1489 nt Sma I cDNA fragment from nt 4272 to 5760 (inclusive). This deletion resulted in the frame-shift of the portion of PPV genome after nt 5760 [##REF##16777241##58##]. The deletion region abolishes the translation of the viral protein domains CI-6k2-NIa/VPg-NIa/Pro-NIB-CP which is required for viral replication. Thus, clone pPPV-SK68Δ is non-infectious.</p> <title>Protoplast isolation and transfection</title> <p>Protoplasts were isolated using fully expanded leaves of 3-week-old <italic>Arabidopsis </italic>accession Col-0 and transfected with the PPV infectious cDNA clone or the mutated PPV cDNA clone following the protocol we have developed recently [##REF##16774842##54##]. Transfected protoplasts were incubated at room temperature for 5 min, and then transferred to 25°C and incubated in the dark for 3, 6 and 12 h, respectively. At selected time points, transfected protoplast cells were harvested, immediately frozen in liquid nitrogen and stored at -80°C prior to RNA extractions. Gene expression differences were averaged out by pooling six RNA samples isolated independently from protoplasts transfected with pPPV-SK68 and pPPV-SK68Δ at each time point after transfection. For time course comparison, three independent biological replicates were maintained. Thus, for three different time points, 18 gene chips in total were used to hybridize with RNA isolated from protoplasts transfected with pPPV-SK68 and pPPV-SK68Δ.</p> <title>RNA extraction, purification and virus detection</title> <p>Total RNA from virus-infected and mock-inoculated leaves of <italic>Arabidopsis </italic>plants and from protoplasts transfected with pPPV-SK68 and pPPV-SK68Δ were isolated using Trizol reagent (Invitrogen) following the manufacturer's suggestion. Total RNA was purified using RNeasy Mini columns (Qiagen) after DNase I treatment (Invitrogen). Integrity of total RNA was assessed using Bioanalyzer 2100 (Agilent Technologies, Boblingen, Germany). Primers OL-PPVDF3 (5'TCAATGGAATGTGGGTGATGA 3') and OL-PPVDR6 (5' GATTAGACTCTCACCCAGGTA 3'); and PPVMF18 (5' GGAAATTGAGAGATACCT 3') and PPVMR16 (5' TTTGTCTAAAAGTGGGTTT 3') were used to detect PPV in the virus-infected or mock-inoculated control samples and in pPPV-SK68- and pPPV-SK68Δ-transfected protoplasts, respectively.</p> <title>Target preparation and hybridization</title> <p>Standard RNA processing and hybridization protocols were followed as recommended by the Affymetrix manual (Santa Clara, CA, USA). A double-amplification approach (Affymetrix) was used to generate cRNA for hybridization. Total RNA isolated from virus-infected or mock-inoculated leaves were subjected to double-amplification. Labelling, hybridization and detection of hybridization signals were performed at the McGill University and Genome Quebec Innovation Centre (please see Availability &amp; requirements for more information).</p> <title>Data analysis</title> <p>Hybridization signals obtained from the scanned array images were quantified using MAS5.0 expression analysis algorithm, implemented in the Affymetrix GCOS (GeneChip<sup>® </sup>Operating) software version 1.2. Raw MAS5.0 data files obtained from scanned array images were imported into GeneSpring microarray analysis software version 7.3.1 (Silicon Genetics, Redwood City, CA, USA). Systemic variations across the chips were normalized using per gene and per chip normalization methods available in the GeneSpring software. Mean normalized value of the signal intensity for each gene in three biological replicate experiments was adopted as the expression value of the gene. The significant changes in gene expression during PPV infection in infected leaves and in protoplasts were determined using <italic>P </italic>values (<italic>P </italic>≤ 0.05) derived from one-way ANOVA and the significant changes in the values of the identified genes were further adjusted using multiple testing correction of the Benjamini and Hochberg false discovery rate procedure [##UREF##5##21##] which applied a stringent cut off of 0.05 that corresponds to <italic>Q </italic>≤ 0.05 [##REF##12883005##22##]. 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 order to identify genes that show significant changes of transcript levels during virus infection, we further searched for genes that were significantly up- (≥ 2.5-fold) or downregulated (≤ -2.5-fold) relative to the control. For clustering analysis, gene expression profiles from the protoplasts transfected with the PPV infectious clone at each time point were clustered based on their similarity in expression using hierarchical average linkage algorithm and Pearson correlation distance metric implemented in the Cluster 3.0. The clustering results were visualized by a gene tree heat map using the TreeView program [##REF##9843981##82##]. K-means clustering was employed to cluster the gene expression data from the PPV transfected protoplasts into several distinct expression profiles using Pearson correlation distance metric. The microarray data from this work was submitted to the NCBI database with Gene Expression Omnibus (GEO) accession number [<ext-link ext-link-type="gen" xlink:href="GSE11217">GSE11217</ext-link>].</p> <title>Functional categorization of genes</title> <p>Each gene was assigned with a putative function by blasting the AGI probe identifier sequences against the <italic>Arabidopsis </italic>database deposited at TAIR (the <italic>Arabidopsis </italic>Information Resources) using the BLASTX program. Functional classification was performed by assigning the function of each gene based on protein sequence similarity using MIPS (Munich Information Centre for Protein Sequences; please see Availability &amp; requirements for more information) <italic>Arabidopsis </italic>functional classification schema. Genes without a category name or with unclear classification were labeled as "unknowns". Over-or under-representation of genes involved in each functional group was evaluated using Fisher's exact test function [##REF##14668247##23##,##REF##12620386##24##].</p> <title>Determination of <italic>P. persica </italic>orthologs to <italic>Arabidopsis </italic>genes induced by PPV</title> <p>The PPV-infected <italic>P. persica </italic>leaf (PPL) EST library [##UREF##15##64##], consisting of 2156 uniEST sequences extracted from NCBI dbESTs [GenBank: <ext-link ext-link-type="gen" xlink:href="DN552797">DN552797</ext-link>, <ext-link ext-link-type="gen" xlink:href="DN556565">DN556565</ext-link>] were formatted so that they were searchable using the tBLASTx algorithm. Sequences of 2013 PPV induced and 1457 PPV repressed <italic>Arabidopsis </italic>transcripts in infected leaf tissues respectively were aligned to this collection to search for the presence of <italic>P. persica </italic>orthologs using the tBLASTx algorithm with default settings [##REF##9254694##83##]. The BLAST output for each transcript was then parsed for high-scoring pair (HSP) and the associated E-value. Any hits with HSP 100 and E-value 10<sup>-20 </sup>were taken as an indicative of significant similarity [##REF##10731134##84##].</p> <title>Cross comparison of genes differentially regulated by PPV infection against the <italic>Arabidopsis </italic>genes differentially regulated by infections of other positive sense RNA viruses</title> <p>The <italic>Arabidopsis </italic>genes differentially regulated in PPV-infected leaf tissues were cross compared by searching against the <italic>Arabidopsis </italic>gene set that were differentially expressed in response to the infection by other positive sense RNA viruses [##REF##12906112##7##, ####REF##12535341##8##, ##REF##15111722##9##, ##REF##15604696##10####15604696##10##,##REF##17427806##16##].</p> <title>Northern hybridizations and semi quantitative RT-PCR (sqRT-PCR)</title> <p>Northern hybridizations and sqRT-PCR were carried out to verify the microarray data. Probes for sqRT-PCR and Northern hybridizations were generated by PCR amplification of <italic>Arabidopsis </italic>cDNA with gene specific primers (Table ##TAB##1##2##). <italic>Arabidopsis </italic>Actin 2 (ACT2; At3g18780) was used as an internal control with the following primer sets; 5' GCCATCCAAGC TGTTCTCTC 3' and 5' GAACCACCGATCCAGACACT 3' to normalize small differences in template amounts. For Northern hybridization, 10 μg of total RNA was denatured with glyoxal, electrophoretically separated on 1.5% agarose gels, blotted onto a nylon membrane and UV cross-linked (1200 μJ/cm<sup>2</sup>). For sqRT-PCR, amplicons of RT-PCR were electrophoretically separated on a 1% agarose gel, transferred to a nylon membrane by capillary blotting and cross-linked as above. Hybridizations for Northern blots and sqRT-PCR were performed with at least three independent biological replicates. The probes were radiolabeled with [α-<sup>32</sup>P] dCTP by Ready-To-GO DNA labeling kit following the manufacturer's instructions (Amersham Biosciences). A Phosphor Imager (Imaging Screen K, Personal Molecular Imager<sup>® </sup>FX, Bio-Rad) and Quality One Quantitation software v. 4.2 (Bio-Rad) were used for visualization and quantification of radioactive signals.</p>
<title>Results and discussion</title> <title>Symptom development in PPV-infected <italic>Arabidopsis </italic>plants</title> <p>Susceptibility of <italic>Arabidopsis </italic>accession Col-0 to the infection of a PPV-D strain isolated in Canada was monitored by recording the development of symptoms and the virus titre. No typical viral symptoms such as mosaic, chlorosis, ringspot, vein yellowing, or leaf curling were visible in plants inoculated with PPV at 17 dpi and up to 42 dpi. This is consistent with the recent finding that <italic>Arabidopsis </italic>plants infected by two PPV-D isolates displayed mild or no symptoms [##REF##16673941##20##]. At 17 dpi, the most obvious phenotypes observed in this study were mild growth retardation, slightly delayed flowering and bolting (Figure ##FIG##0##1A##). These phenotypes were not documented for <italic>Arabidopsis </italic>accession Col-0 infected with two PPV-D isolates but were very similar to accession Bay-O inoculated with one of the two PPV-D isolates [##REF##16673941##20##]. Considering the diverse phenotypes of different accession-isolate combinations, the observed phenotypic variation between the two studies is likely due to the PPV isolates used. Using reverse transcription-polymerase chain reaction (RT-PCR) and enzyme-linked immunosorbent assay (ELISA), PPV was detected in newly developed rosette leaves above the inoculated leaves (Figure ##FIG##0##1B##) at 17 dpi. Since PPV-induced symptoms appeared and the virus was detected by RT-PCR and ELISA at 17 dpi in the infected leaf tissues, it was expected that changes in host gene expression would be altered at this time point. In addition, we also inoculated <italic>Nicotiana benthamiana </italic>plants with PPV. Infected leaves showed chlorosis and severe mosaic symptoms as early as 7 dpi (data not shown). At 17 dpi, the virus titre in <italic>N. benthamiana </italic>was at least 100 times or higher than that in <italic>Arabidopsis </italic>(data not shown), suggesting that <italic>Arabidopsis </italic>may have stronger basal or inducible resistance to PPV than <italic>N. benthamiana</italic>.</p> <title>Over- or under-represented functional categories of differentially expressed <italic>Arabidopsis </italic>genes in PPV-infected leaf tissues</title> <p>At 17 dpi, newly expanded PPV-infected rosette leaves from the inoculated plants and the corresponding leaves from the mock-inoculated plants were sampled for microarray hybridizations. Consistency of the microarray data was assessed by comparing the mean normalized signal intensities obtained from the mock-inoculated control samples across three independent biological replicates. A high correlation (0.96 – 0.98) among the mock inoculated control samples was observed, indicating low biological variability between the replicates.</p> <p>Gene expression data were analyzed by one-way ANOVA to identify differentially expressed genes in PPV-infected rosette leaves from the inoculated plants relative to leaves from the mock-inoculated plants. Using stringent selection criteria with the Benjamini and Hochberg false discovery rate (FDR) of 5% [##UREF##5##21##], corresponding to <italic>Q </italic>≤ 0.05 [##REF##12883005##22##], ~31.4% of the non-redundant (nr) genes (7,151 out of the 22,810 <italic>Arabidopsis </italic>genes printed on the chip) were identified as significantly differentially expressed (either induced or repressed) in response to PPV infection. These 7,151 genes were further filtered using a 2.5-fold increase or decrease in signal intensity. Using this fold change cut-off, we identified 2,013 and 1,457 <italic>Arabidopsis </italic>genes that were significantly induced and repressed by PPV in systemically infected tissues, respectively.</p> <p>Each identified gene was then assigned to a functional class according to the <italic>Arabidopsis </italic>MIPS (Munich Information Centre for Protein Sequences) functional classification scheme (Figure ##FIG##1##2##). With the information gathered from TAIR (The <italic>Arabidopsis </italic>Information Resource) database (please see Availability &amp; requirements for more information), we were able to assign a putative function for 54.3% of the total 22,810 genes on the array, while the remaining 45.7% were undefined as they were annotated as "unknowns or orphans" (Table ##TAB##0##1##). The expression levels of <italic>Arabidopsis </italic>genes significantly induced or repressed by PPV were assigned to 12 different major functional groups. In some major functional groups, genes were further separated into subcategories to simplify the data for biological interpretation [see Additional file ##SUPPL##0##1##]. The number of genes in each functional category was subjected to Fisher's exact test [##REF##14668247##23##,##REF##12620386##24##] to determine if genes involved in certain biological processes were over- or under-represented (<italic>P </italic>≤ 0.001). The functional groups significantly over-represented in the set of upregulated genes were carbohydrate/soluble sugar/starch/amino acid metabolism (<italic>P </italic>&lt; 7E-12), primary and secondary metabolism (<italic>P </italic>&lt; 2.2E-07), intracellular membrane/membrane-bound organelle (<italic>P </italic>&lt; 4.5E-06), clathrin binding (<italic>P </italic>&lt; 1.2E-05), glycolysis (<italic>P </italic>&lt; 2.0E-05), protein folding/heat shock/chaperone activity (<italic>P </italic>&lt; 5.2E-06), cell wall-associated transcripts (<italic>P </italic>&lt; 1.6E-05), RNA metabolism (<italic>P </italic>&lt; 1.9E-05), chloroplast (<italic>P </italic>&lt; 1.8E-06), protein fate (<italic>P </italic>&lt; 3.5E-07), ubiquitin-like conjugating enzyme activity (<italic>P </italic>&lt; 4.1E-05) and development/storage proteins (<italic>P </italic>&lt; 4.8E-08). Among the downregulated genes, the over-represented functional groups were metal ion binding/metal binding (<italic>P </italic>&lt; 1.6E-07), development/storage proteins (<italic>P </italic>&lt; 3.3E-06), cell wall-associated transcripts (<italic>P </italic>&lt; 0.00013) and protein synthesis and translation (<italic>P </italic>&lt; 1.4E-43). Genes involved in primary/secondary metabolism (<italic>P </italic>&lt; 1.4E-06) and intracellular membrane/membrane-bound organelle (<italic>P </italic>&lt; 1.3E-05) functional groups were significantly under-represented.</p> <p>The analysis showed that that the number of upregulated genes was more than that of downregulated genes in PPV-infected <italic>Arabidopsis </italic>leaves, probably reflecting more diverse genes induced by PPV infection. This notion is supported by a disproportionate allocation in the number of altered genes in the infected leaf tissues belonging to major functional groups of metabolism, transcription/splicing/RNA processing proteins, defence and development/storage proteins.</p> <title>Gene expression changes in metabolism, chloroplast and protein fate in PPV-infected <italic>Arabidopsis </italic>leaf tissues likely involved in symptom development</title> <p>Virus infections often disturb the general metabolism in the host. This includes the accumulation of soluble sugars (e.g., sucrose synthase) and starch (e.g., starch synthase), an increase in the rate of respiration, a decrease in the rate of photosynthesis and elevated levels of amino acid synthesis (e.g., amino acid permease and aminotransferase) [##UREF##6##25##, ####UREF##7##26##, ##REF##12226342##27##, ##UREF##8##28##, ##REF##17170466##29##, ##UREF##9##30####9##30##]. In this study, we found that in PPV-infected leaves, 23% (124/548) of genes associated with metabolism of starch, soluble sugars and amino acids and 2.7% (15/548) of genes related to glycolysis showed elevated expression levels [see Additional file ##SUPPL##0##1##]. Appearance of disease symptoms such as chlorosis on the <italic>Cucumber mosaic virus </italic>(CMV)-infected leaf tissues in cucumber was correlated with increased glycolysis, phosphoenol pyruvate phosphate, respiration and starch accumulation, decreased photosynthesis and reduction in total protein synthesis [##REF##12226342##27##]. The increased expression of these genes may result from the physical blockage of the transport path, the inactivation of sugar transport proteins or the induction of cell wall invertase [##UREF##10##31##, ####REF##10508760##32##, ##REF##12663220##33####12663220##33##]. Indeed, several cell wall invertase genes were also induced in PPV-infected leaves in this study. These data indicate that a series of intercrossing pathways involved in sugar metabolism was differentially regulated by PPV infection. It is possible that altered gene expression pathways involved in sugar metabolism might contribute to the mild symptom development observed in PPV-infected <italic>Arabidopsis </italic>leaves.</p> <p>In addition to soluble sugars, starch and amino acids, PPV also induced 16% (87/548) of membrane-associated proteins, including proteins associated with integral membrane, plasma membrane, endoplasmic reticulum (ER), vacuole, mitochondria and peroxisomes [see Additional File ##SUPPL##0##1##]. Recently, membrane proteins have also been found to be upregulated in <italic>N. benthamiana </italic>plants infected with <italic>Tomato ringspot virus </italic>(ToRSV) and PPV [##REF##17722703##13##]. For example, genes encoding proteins such as plasma membrane polypeptide, putative integral membrane, COP-coated vesicle membrane, clathrin protein and mitochondrial malate dehydrogenase were induced in PPV-infected <italic>Arabidopsis </italic>leaves that have similar functional identity to those upregulated in <italic>N. benthamiana </italic>plants infected with ToRSV and PPV. It is believed that many positive-strand RNA viruses replicate their genome in association with cellular membranes. Viral proteins that are required for the formation of viral replication complexes (VRCs) are often shown to be located in ER [##REF##17292435##34##, ####REF##10579716##35##, ##UREF##11##36####11##36##]. It would be interesting to investigate whether these upregulated membrane proteins are co-localized to VRCs. Characterization of such host membrane factors will provide useful information regarding molecular mechanisms governing membrane specificity for viral RNA replication.</p> <p>Of the upregulated genes coding for nuclear or cytosolic proteins in PPV-infected <italic>Arabidopsis </italic>leaves, two OGTs (O-linked N-acetyl glucosamine transferase), named SEC (SECRET AGENT) and SPY (SPINDLY, a gibberellin response protein), which can potentially posttranslationally modify the viral capsid protein (CP), are very interesting. Study by Chen et al. [##REF##16014901##37##] has shown that the N-terminus of the CP of PPV is modified with O-linked N-acetyl glucosamine in <italic>N. clevelandii</italic>, <italic>N. benthamiana </italic>and <italic>Arabidopsis</italic>, indicating that the modifications are not host specific. In <italic>Arabidopsis</italic>, this modification is catalyzed by SEC but not by SPY. Such modification is not essential but plays a role in the infection process [##REF##16014901##37##]. SEC and SPY have partial functional redundancy in plants [##REF##16014901##37##,##REF##12136030##38##]. It is not clear whether the upregulation of SPY in PPV-infected <italic>Arabidopsis </italic>is a compensation for the loss of SEC that is used for the modification of the PPV CP.</p> <p>We also observed as many as 52 chloroplast genes that were upregulated in PPV-infected leaf tissues [see Additional file ##SUPPL##0##1##]. These represent approximately 60% of the chloroplast genes in the <italic>Arabidopsis </italic>genome [##REF##10574454##39##]. Several studies suggest that differential regulation of chloroplast gene expression is correlated to viral symptoms such as mosaic and chlorosis [##REF##17722703##13##,##REF##17427806##16##]. Recently, Dardick [##REF##17722703##13##] has reported that ToRSV or PPV infection in <italic>N. benthamiana </italic>leaves induces severe mosaic and chlorosis symptoms. However, the infection of this host with <italic>Prunus necrotic ringspot virus </italic>(PNRSV) hardly induces any symptoms. Using potato cDNA chips, it has been shown that infections of ToRSV and PPV result in the suppression of 34% (83/243) of chloroplast-related genes, while no common genes are among the induced genes in ToRSV- and PPV-infected <italic>N. benthamiana</italic>. In contrast to other two viruses, PNRSV infection downregulates 67% (10/15) and upregulates 33% (5/15) of the chloroplast-related genes [##REF##17722703##13##]. Analogous to PNRSV-infected <italic>N. benthamiana</italic>, PPV-infected <italic>Arabidopsis </italic>leaves did not show typical viral symptoms such as mosaic and chlorosis. The induction of such a large number of chloroplast-related genes in <italic>Arabidopsis </italic>may help maintain the normal function and morphology of the infected leaf and counteract possible detrimental effects resulting from PPV infection.</p> <p>In virus infected plants, protein activities are often altered through posttranslational modifications such as ubiquitination [##REF##17417941##4##]. In the current study, about 9% (60 genes) of the upregulated protein fate genes encoding proteins associated with the 26S proteasome were very significant (<italic>Q </italic>≤ 0.05) and over-represented (<italic>P </italic>&lt; 3.5E-07). Of the seven upregulated 26S proteasome genes, two genes code for AAA-ATPase subunits including RPT3 that reside in the mature 26S proteasome [##REF##17427806##16##] and five genes including <italic>Rpn1</italic>, <italic>3 </italic>and <italic>6 </italic>are involved in ubiquitination. Fundamentally, the 26S proteasome is known to control many biochemical processes by programmed degradation of regulatory protein targets [##REF##16905670##40##]. Upregulation of these genes may interfere with ubiquitin-dependent reactions in plants [##REF##11473704##41##]. It has been reported that, the 26S proteasome degrades the virus-encoded movement protein that facilitates the viral cell-to-cell spread [##REF##10708450##42##]. Therefore, the upregulation of the 26S proteasome genes may be a resistance response to PPV infection.</p> <p>On the other hand, in <italic>N. benthamiana </italic>plants infected with ToRSV, about 72% (33/46) of the proteasome genes have been shown to be upregulated [##REF##17722703##13##]. A few archetype proteasome proteins such as subtilisin-like protease, aspartyl protease, ubiquitin conjugating enzyme, F-box protein and 26S proteasome regulatory subunits induced by PPV in <italic>Arabidopsis </italic>leaves are upregulated by ToRSV in <italic>N. benthamiana </italic>plants [##REF##17722703##13##]. Upregulation of such protein fate and ubiquitin genes has also been observed in other studies during viral infection [##REF##15604696##10##,##REF##17427806##16##,##REF##8986804##43##,##REF##11041891##44##]. Study by Marathe <italic>et al</italic>. [##REF##15604696##10##] has demonstrated that the induction of F-box proteins may mediate the degradation components of the CMV-Y resistance pathway. Though the above mentioned protease genes induced by PPV infection are known to function in the protein degradation, plant defence responses and hypersensitive reactions [##REF##11976458##45##, ####REF##11896264##46##, ##REF##8843949##47####8843949##47##], their exact roles in PPV infection are to be elucidated.</p> <title>Association of the repression of development/storage proteins, protein synthesis and cell wall-related genes in <italic>Arabidopsis </italic>leaf tissues with PPV infection and symptom development</title> <p>Many virus-induced symptoms have been linked to alterations in plant hormone synthesis [##REF##17417941##4##,##REF##15681455##48##,##REF##16903352##49##]. Of the 96 development/storage genes identified in this study, 10 auxin, 5 ethylene and 4 gibberellin-related genes were repressed in PPV-infected leaves [see Additional file ##SUPPL##1##2##]. Recent studies [##REF##15681455##48##,##REF##16903352##49##] on TMV infection in <italic>Arabidopsis </italic>have shown that during the process of infection the TMV replicase protein interacts with a subset of Aux/IAA proteins, resulting in the activation of auxin responsive transcription factors (ARF). Such an alteration in the transcription levels of auxin-responsive genes and the disruption of host normal physiological process have shown to be associated with disease symptoms. In this study, PPV infection suppressed the expression of two such genes encoding auxin IAA proteins (At5g57420 and At1g52830). With respect to the downregulation of ethylene responsive genes, though there is no direct evidence suggesting interactions between virus and ethylene signalling components, specific viral proteins have been suggested to alter the ethylene signalling pathway [##REF##17417941##4##]. For instance, the P6 protein of <italic>Cauliflower mosaic virus </italic>(CaMV) is associated with the disruption of the ethylene response pathway and this may play a key role in the induction of symptoms such as mild vein chlorosis to severe chlorosis and stunting [##REF##9390424##50##]. In the case of downregulation of gibberellic acid (GA), a direct interaction between the <italic>Rice dwarf virus </italic>(RDV) P2 protein and the rice ent-kaurene oxidase, a key enzyme in the synthesis of GA in plants was found to be responsible for the significant reduction of GA1 in RDV-infected plants, and in the appearance of dwarf phenotype [##REF##16299167##51##]. Collectively these examples demonstrate that the suppression of hormone-related genes by viruses induces symptoms through altering the host hormone and developmental signalling pathway. It would be interesting to know whether any of the observed downregulated hormone-related genes is associated with the reduced growth and delayed development of PPV-infected <italic>Arabidopsis</italic>.</p> <p>Another significantly over-represented subset of genes suppressed in PPV-infected tissues includes genes coding for the protein synthesis and translation machinery [see Additional file ##SUPPL##1##2##]. This functional category contains 33 ribosomal subunits: 21 large, 11 small and 1 plastid specific ribosomal proteins. In the <italic>Arabidopsis </italic>genome, there are ~227 ribosomal genes [##REF##17722703##13##,##REF##14681386##52##]. About 15% of these ribosomal genes were printed on the array. Based on the analysis, we observed that approximately 81% (27/33) of the small and large ribosomal subunit genes were significantly repressed in PPV-infected leaves.</p> <p>However, an opposite trend was observed in PPV-infected <italic>N. benthamiana </italic>plants [##REF##17722703##13##], where 80% (105/131) of potato ribosomal genes were upregulated. Comparison of these potato upregulated ribosomal genes in <italic>Arabidopsis </italic>through BLAST (E-value &lt; -10) sequence similarity searches identified 55 nuclear ribosomal subunits [##REF##17722703##13##]. Comparative analysis of 27 <italic>Arabidopsis </italic>nuclear ribosomal subunits that were downregulated by PPV in our study with those in Dardick's study (using Table S5 from Dardick's study for comparison) showed 8 (7 large and 1 small) <italic>Arabidopsis </italic>ribosomal subunits with an opposing trend. Since plant viruses lack protein synthesis ability, they rely on the host cell translation machinery to produce viral proteins. Thus, PPV might benefit from the upregulation of ribosomal genes in <italic>N. benthamiana </italic>which showed more susceptibility to PPV and displayed much more severe symptoms than <italic>Arabidopsis</italic>. However, at present there is no evidence to discern whether the observed downregulation of ribosomal genes is due to the stronger basal resistance of <italic>Arabidopsis </italic>to PPV infection or due to the natural intricacy of the different host responses to PPV that are diverse depending on the tissue type, time post infection, temperature, host and virus strain, or simply due to an artifact of hybridizations, given that two different platforms were used in the studies.</p> <p>Cell wall-related genes are major determinants of cell morphogenesis in plants, encoding several hundreds of different structural proteins and cell wall-related enzymes [##REF##11554474##53##,##REF##16774842##54##]. In this study, 75 cell wall-related genes encoding protein for structure (23 genes), reassembly (34 genes), expansins (9 genes) and matrix polymers (9 genes) were significantly downregulated in PPV-infected leaves [see Additional file ##SUPPL##1##2##]. This is consistent with results observed in <italic>Arabidopsis </italic>plants infected by TuMV [##REF##17427806##16##] and in rice plants infected with RDV [##REF##17378427##15##]. In these studies, repression of cell wall related genes is correlated with symptom development. Of these downregulated genes, extensin is a cell wall-localized hydroxyproline-rich glycoprotein that usually forms a cross-linked network with pectin to create a highly impassable barrier against pathogens [##REF##7659756##55##]. Downregulation of extensin during PPV infection suggests a decreased extensin deposition [##UREF##12##56##]. Conceivably, PPV exploits the plasmodesmata to move from initially infected cells to the neighbouring healthy cells. It is not likely that PPV cell-to-cell movement is related to the deposition of extensin. Nonetheless, the molecular mechanisms and the biological connotation underlying the downregulation of cell wall-related genes in the infected leaf remain to be understood.</p> <p>In addition to the above described genes, many genes in the category of metal ion binding or metal binding were also downregulated. It is not clear if the downregulation of these genes can lead to physiological changes that are associated with PPV infection.</p> <title>Identification of genes differentially regulated by PPV in infected <italic>Arabidopsis </italic>protoplasts</title> <p>Resistance induced through signal transduction pathways in <italic>Arabidopsis </italic>may take place immediately after PPV infection. The <italic>Arabidopsis </italic>leaves infected with PPV at 17 dpi contained a mixture of uninfected, infected, and post-infected cells. Thus, the spatial and temporal information on gene expression was lost. Sampling at earlier time points such as 1 or 2 dpi can reduce the percentage of post-infected cells but will largely increase the percentage of non-infected cells. So, protoplasts have been suggested to overcome this barrier [##UREF##3##17##]. There are several advantages in using <italic>Arabidopsis </italic>protoplasts to study plant-virus interactions. First, in contrary to unpredictable infections on the leaf, synchronized protoplast cells allows for a more precise determination of gene expression. Second, a much higher proportion of the cells are infected than that in infected leaf tissues [##UREF##13##57##]. For example, we have shown that in <italic>Arabidopsis </italic>protoplasts transfected with plasmid pPPV-SK68, a PPV infectious clone, the transfection efficiency is approximately 35% [##REF##16777241##58##]. Hence, the protoplast system makes it possible to study molecular events as early as just a few hours post infection. Like other experimental systems, the protoplast system also has some limitations. Protoplasts lack the cell wall so the physiological integrity of the plant organ is disrupted. During isolation, transfection and culture of protoplasts, gene expression changes induced by harsh treatments such as osmotic stress, enzymatic digestion and centrifugation may mask those induced by virus infection [##UREF##3##17##]. Physiological milieu that a virus encounters in protoplast cells might be quite different from that of an intact leaf cell. Sample bias may result from the minute amounts of RNA isolated from the protoplasts and from the subsequent amplification step. However, improved technologies and experimental designs may be used to minimize the negative impact of the system [##REF##7731989##59##, ####REF##3409786##60##, ##REF##16027030##61##, ##UREF##14##62####14##62##]. Recently, we reported a protocol describing the isolation and transfection of <italic>Arabidopsis </italic>protoplasts for studying PPV-host interactions [##REF##16777241##58##]. In order to identify early time points for profiling global gene expression changes in infected protoplast cells, a time course experiment was performed. <italic>Arabidopsis </italic>protoplasts were transfected with a PPV infectious cDNA clone, pPPV-SK68 derived from an M strain and with a non-infectious deletion mutant, pPPV-SK68Δ as control, respectively [##REF##16777241##58##]. In the mutant clone, the deletion abolishes the translation of viral proteins CI, 6K2 and NIa/VPg required for PPV replication. The quantitative results showed that the viral RNA in the protoplasts was detectable at 3 hpt, accumulated at 6 hpt and peaked at 12 hpt (Figure ##FIG##0##1C##), followed by a decline at 24 hpt [##REF##16777241##58##]. The time course progression of PPV infection in <italic>Arabidopsis </italic>protoplasts is similar to that of <italic>Cymbidium mosaic virus </italic>(CyMV) in orchid protoplasts [##REF##8249447##63##]. Thus, three time points, 3, 6 and 12 hpt were selected for this study.</p> <p>Gene expression in PPV-infected <italic>Arabidopsis </italic>protoplasts was analyzed by one-way ANOVA to identify genes differentially expressed at each time point. The degree of biological variability and reproducibility of the microarray data was assessed by comparing the mean values of the normalized signal intensity from the protoplasts transfected with a non-infectious mutant, pPPV-SK68Δ with all three independent biological replicates at each time point. The analysis resulted with a high correlation coefficient of 0.90 – 0.95, indicating low technical and biological variability between the samples used in the hybridizations. Using FDR at 5%, we identified ~8.4% (1,907 out of 22,810) of non-redundant genes that were significantly (<italic>Q </italic>≤ 0.05) differentially regulated in protoplasts transfected with the PPV infectious clone, pPPV-SK68 at three different time points relative to the control transfected with the non-infectious deletion mutant, pPPV-SK68Δ. These genes were further filtered using a 2.5-fold increase or decrease in signal intensity, resulting in 263 nr genes that were induced (≥ 2.5-fold), and 304 genes that were repressed (≤ -2.5-fold), as illustrated in the Venn diagrams (Figure ##FIG##2##3A## and ##FIG##2##3B##). A comparison analysis between 263 induced and 304 suppressed genes identified 107 genes induced at one time point and were not repressed at other time points, 148 genes repressed at one time point and were not induced at other time points, and 156 genes induced and repressed at different time points (Figure ##FIG##2##3C##) [see Additional file ##SUPPL##2##3##]. The analysis also revealed that the number of repressed genes was greater than that of the induced genes at all three time points analyzed. The number of genes that were commonly repressed among the time points was 14 between 3 and 12 hpt, 10 between 3 and 6 hpt, and 12 between 6 and 12 hpt. We further examined each of these 263 induced and 304 repressed genes to see if they were also consistently induced or repressed at all three different time points of infection. Fifteen genes were found to be constantly repressed and 4 were consistently induced in the PPV transfected protoplasts, suggesting that the proportion of genes repressed at all three time points were greater than that of the genes induced.</p> <p>The genes induced or repressed in PPV-infected protoplasts were annotated, assigned to a predicted function and classified into 12 different major functional categories using the <italic>Arabidopsis </italic>MIPS functional classification scheme (Figure ##FIG##3##4##). Fisher's exact test [##REF##14668247##23##,##REF##12620386##24##] was used to determine whether the representation of the number of genes in each functional category was significant (<italic>P </italic>≤ 0.001) or by random chance. Of the upregulated genes, transcripts involved in primary and secondary metabolism (<italic>P </italic>&lt; 6.2E-06), cellular signalling (<italic>P </italic>&lt; 4.4E-11), defence response (<italic>P </italic>&lt; 2.7E-06), transcription factors (<italic>P </italic>&lt; 1.1E-11) and transporters (<italic>P </italic>&lt; 1E-12) were over-represented, while in the case of downregulation, genes involved in transcription factors (<italic>P </italic>&lt; 5E-09), primary and secondary metabolism (<italic>P </italic>&lt; 6.4E-05), kinases or cellular signalling (<italic>P </italic>&lt; 1.3E-05), and defence responses (<italic>P </italic>&lt; 5.6E-05) were over-represented. A complete list of induced or repressed genes in PPV-infected protoplasts in different functional groups is provided [see Additional file ##SUPPL##2##3##].</p> <p>As evident [see Additional file ##SUPPL##3##4##; Additional file ##SUPPL##4##5##; Additional file ##SUPPL##5##6##], comparison of the differentially expressed genes between PPV-infected <italic>Arabidopsis </italic>protoplasts (263 up- and 304 downregulated genes) and PPV-infected leaf tissues (2,013 induced genes and 1,457 suppressed genes) resulted in only 32 and 33 common genes that were significantly induced (≥ 2.5-fold) and repressed (≤ -2.5-fold), respectively. A few vital genes such as genes encoding heavy metal protein, disease resistance protein, FtsH protease, cytochrome P450 and β-1, 3-glucanase were significantly differentially regulated by PPV in the infected leaves and transfected protoplasts. The smaller percentage of significantly differentially regulated genes observed in the PPV-infected protoplasts than in the PPV-infected leaf suggests a remarkable gene expression differences between the two systems. However, as discussed above, the differences in the alteration of gene expression changes in the protoplasts induced by PPV infection could have been masked by harsh treatments during protoplast isolation, transfection and culture or might result from the two different PPV strains used. Though the PPV-D and -M strains used in this study share more than 95% sequence identity [##UREF##15##64##], these two strains may induce different symptoms on the same host [##REF##16673941##20##,##UREF##16##65##].</p> <p>To visualize gene expression profiling at all three time points, a hierarchical clustering analysis was carried out on the 263 induced and 304 repressed genes (Figure ##FIG##4##5A##). A gene tree heat map was built using Pearson correlation distance metric with an average linkage algorithm. To understand further the temporal relationship of the host gene expression in PPV-infected protoplasts at each time point, twelve distinct groups of expression pattern clusters were generated on the significant differential gene expression datasets using K-means clustering with Pearson correlation distance metric (Figure ##FIG##4##5B##) [see Additional file ##SUPPL##6##7##]. A complete list of genes induced or repressed at each time point in different cluster groups along with their expression levels and putative functions is provided [see Additional file ##SUPPL##7##8##]. The clustering analysis also revealed greater changes in gene expression (236 genes) at 6 hpt than at 3 hpt (163 genes) and at 12 hpt (201 genes) in PPV-infected <italic>Arabidopsis </italic>protoplasts. These observations suggest that early gene expressions in response to PPV infection are active and dynamic.</p> <title>Early upregulated genes related to disease resistance or defence responses in PPV-infected <italic>Arabidopsis </italic>protoplasts</title> <p>PPV infection in protoplasts induced the expression of 263 genes. In contrast to about 11% of the suppressed genes that were common at different time points, there were no induced genes shared at different time points, suggesting induction of host gene expression was much more transient than suppression. Genes related to defence, cellular signalling, primary and secondary metabolism, transcription and transporters were upregulated and over-represented. Within these five categories, ~52% (48/93) of genes were predominantly upregulated (≥ 2.5-fold) at 6 hpt [see Additional file ##SUPPL##2##3##]. Examples of the upregulated genes in the above functional categories include MADS box protein, zinc finger family protein, WRKY family transcription factor, avirulence (Avr)-induced protein, hydrophobic protein, disease resistance TIR-NBS-LRR (toll-interleukin-1-receptor/nucleotide-binding site/leucine-rich repeat), protein kinases and phosphatases.</p> <p>Another interesting gene upregulated by PPV encodes a disease resistance protein, TIR-NBS-LRR. This gene was upregulated with 6.39-fold in PPV-infected protoplasts at 6 hpt and 2.12-fold at 12 hpt, and 10.84-fold in PPV-infected leaves [see Additional file ##SUPPL##3##4##]. In the <italic>Arabidopsis </italic>(Col-0 accession) genome, there are 149 NBS-LRR-like genes [##UREF##17##66##,##REF##12671079##67##]. The resistance gene <italic>NBS-LRR </italic>mostly confers resistance to viral diseases but not necessarily associates with cell death or tissue necrosis [##UREF##17##66##]. Genes <italic>Sw5</italic>, <italic>N </italic>and <italic>Rx </italic>are homologous to the <italic>NBS-LRR </italic>gene. <italic>Rx </italic>resistance to <italic>Potato virus X </italic>(PVX; Potexvirus) is the most classical example, which shows an extreme resistance that inhibits virus replication without hypersensitive cell death [##REF##10330465##68##]. The <italic>Sw5 </italic>gene confers resistance to several Tospoviruses [##REF##11043474##69##], and <italic>N </italic>and <italic>Rx </italic>are effective against many natural variants of TMV and PVX [##UREF##18##70##,##UREF##19##71##]. Since the specificity determinants subsist in LRR domains [##UREF##17##66##], it would be interesting to find out if the NBS-LRR gene detected in the PPV-infected leaf plays a role in resistance responses to PPV infection.</p> <p>The comparison of induced genes at different time points in PPV-infected protoplasts revealed that genes upregulated at 3 hpt might regulate their target genes at 12 hpt. For instance, two genes coding for cnd41 (At5g24820), a plastid specific DNA binding factor, and calmodulin (At3g15050, At3g16490), a signalling binding protein were induced at 3 hpt and their target genes encoding metabolic enzymes such as serine carboxypeptidase (At1g11080) were upregulated at 12 hpt [##UREF##13##57##,##REF##16777241##58##]. Similarly, transcription factors such as zinc finger proteins (At5g54020; At3g45480; At4g12140) was upregulated at 3 hpt and the genes encoding enzymes such as GDSL-motif lipase (At3g43550) and glycosyl hydrolases (At3g26140, At4g39010) were upregulated at 12 hpt [##REF##9426602##72##]. These results suggest that virus infection may induce expression of specific upstream regulator genes first which then activate the expression of downstream target genes.</p> <title>Inhibition of host mRNA accumulation in PPV-infected protoplasts</title> <p>Based on the expression analysis (Figure ##FIG##2##3##), the number of repressed genes in the transfected protoplasts was greater than that of the genes induced at all time points, suggesting that at early infection stages, suppression is more active than induction of host gene expression. Moreover, about 11% (33 out of 304) of the suppressed genes overlapped at different time points, indicating a relatively stable suppression of genes. Consistent with this, in pea cotyledons infected with <italic>Pea seed borne mosaic virus </italic>(PSbMV), suppression of diverse host genes including seed storage proteins (e.g., vicilin and convicilin), and starch synthase occurs only in a zone of six to eight cells immediately behind the infection front, which coincided with the onset of viral replication [##REF##17791345##73##]. In the present study, a gene coding for the seed storage protein cupin family (At2g18540) that shares sequence homology to the pea vicilin and convicilin was significantly repressed by PPV infection at 12 hpt in <italic>Arabidopsis </italic>protoplasts. Suppression of such host gene expression or "host shutoff" has been suggested to be due to stress responses, induced defence states or a mechanism employed by viruses to outcompete against host mRNAs for host resources for their translation and replication [##REF##17427806##16##,##UREF##3##17##,##REF##8986804##43##,##REF##11104816##74##]. Study by Rottier [##UREF##20##75##] has demonstrated that protein synthesis in protoplasts is blocked in cowpea protoplasts transfected with <italic>Cowpea mosaic virus</italic>. Despite this inhibition, viral mRNAs are preferentially translated. Since both viruses and hosts share the same translational apparatus, viruses must have a mechanism to avoid this inhibition. It is not clear whether the disappearance of host mRNAs is the consequence of mRNA destabilization triggered by stress response, or host transcription inhibition by the virus, or targeted degradation through a virus-induced ribonuclease activity.</p> <title>Identification of a common set of genes in general stress and defence-related groups differentially regulated by infections of PPV and other positive sense RNA viruses</title> <p>We performed a cross sequence comparison analysis by blasting individual sequence of all the genes significantly differentially regulated by PPV in infected <italic>Arabidopsis </italic>leaves against 2,156 uniEST sequences derived from PPV-infected <italic>P. persica </italic>leaves [##UREF##15##64##]. High-scoring pair (HSP; ≥ 100 score) and E-value (E-value ≤ 10–20) were used to identify orthologs. A total of 170 induced [see Additional File ##SUPPL##8##9##] and 18 repressed <italic>Arabidopsis </italic>genes showed high sequence similarity to <italic>P. persica </italic>uniESTs. Of these 170 induced genes, approximately ~39% (67/170) and ~16% (28/170) of the genes were involved in metabolism and defence, respectively. The defence-related genes include those coding for pathogenesis related (PR)-proteins, heat shock and germin-like proteins (Ger2) [see Additional file ##SUPPL##8##9##]. An increased expression of PR proteins requires salicylic acid (SA) signifying that compatible host-virus interactions induce defence-like responses through a SA dependent pathway [##REF##15728340##76##]. More recently, it has been shown that PPV infection in the apoplastic space of <italic>P. persica </italic>leaves induces a thaumatin-like protein (PR-5) [##REF##17043083##77##]. Induction of such protein is thought to be mediated by an increase in the level of peroxidase accumulation [##REF##17043083##77##]. A similar trend was also observed in PPV infected <italic>Arabidopsis </italic>leaves where genes encoding PR-5 and peroxidase were significantly induced and their counterparts were also found in the uniEST sequences derived from PPV-infected <italic>P. persica</italic>.</p> <p>To determine whether the diverse RNA viruses have the ability to elicit common gene expression changes, we cross-compared the genes that were significantly differentially regulated by PPV infection in this study [Additional file ##SUPPL##0##1##] [Additional file ##SUPPL##1##2##] to those <italic>Arabidopsis </italic>genes that were differentially regulated by infections of other positive sense RNA viruses including CMV [##REF##15111722##9##,##REF##15604696##10##]; TMV [##REF##12906112##7##]; TuMV, and <italic>Oil seed rape mosaic virus </italic>(ORMV) [##REF##17427806##16##], and to five diverse RNA viruses, i.e., TuMV, ORMV, <italic>Turnip vein clearing virus </italic>(TVCV), <italic>Potato virus X </italic>(PVX), and CMV [##REF##12535341##8##]. Collectively, 416 common genes [see Additional file ##SUPPL##9##10##] were identified from this analysis. These genes were further classified into 6 different categories. Group A contains 202 genes induced by infections of PPV and other positive sense RNA viruses. About ~40% (80/202) are metabolic genes and ~8% (17/202) are defence-related genes. Group B includes 41 <italic>Arabidopsis </italic>genes that were induced during PPV infection but repressed by infections of other positive sense RNA viruses. About ~27% (11/41) are involved in metabolism. Group C comprises of 119 <italic>Arabidopsis </italic>genes repressed during PPV infection but induced by infections of other positive sense RNA viruses. For example, a defence-related gene encoding glycosyl transferase (At1g70090), which was repressed in PPV-infected leaves at 17 dpi was observed to be induced more than 1.5-fold in TMV-infected leaves [##REF##12906112##7##] at 14 dpi. Group D consists of 7 <italic>Arabidopsis </italic>genes repressed by infections of PPV and other positive sense RNA viruses. Group E comprises 25 genes that were induced during PPV infection but either induced or suppressed by infections of other positive sense RNA viruses. Group F contains 22 <italic>Arabidopsis </italic>genes repressed during PPV infection but induced during infections of more than one positive sense RNA virus. Previous studies have revealed that in virus-infected <italic>Arabidopsis </italic>and tomato leaves as many as one-third of the genes commonly induced by viruses are known or predicted to be involved in plant defence- and stress-related genes [##REF##15231257##5##,##REF##17378427##15##]. In the current study, about 52% (217/416) of the functional proteins differentially regulated in response to infections of PPV and other positive sense RNA viruses are involved in defence, cellular signalling and metabolism. Twenty-nine cellular signalling genes involved in defence mechanisms including protein kinase (At3g08760, At4g32250), protein phosphatase 2C (At5g02760, At3g05640), pyruvate kinase (At2g36580), serine/threonine protein kinase (At3g17410) and calmodulin-binding family protein (At5g62390) [##UREF##0##3##,##UREF##21##78##,##REF##16838784##79##] were induced by infections of PPV and other positive sense RNA viruses. Substantial induction of these defence, cellular signalling and metabolism genes by PPV and other positive sense RNA viruses in infected leaf tissues suggests that there is a shared pathway modulating their expression.</p> <title>Validation of microarray data by semi-quantitative RT-PCR (sqRT-PCR) and Northern hybridizations</title> <p>Since most changes in gene expression occurred at 6 hpt, 9 transcripts significantly upregulated and 1 transcript significantly downregulated at 6 hpt were selected to validate the microarray data using sqRT-PCR. We also selected 3 other transcripts that were significantly upregulated at 3 or 12 hpt for validating the microarray data through sqRT-PCR. We further confirmed the microarray data of 5 <italic>Arabidopsis </italic>genes that were differentially regulated in both PPV-infected protoplasts and leaves using sqRT-PCR or Northern hybridizations. Overall, the expression patterns of all the 13 genes analyzed by sqRT-PCR were consistent with those by microarray hybridizations (Figures ##FIG##5##6## and ##FIG##6##7##) [see Additional file ##SUPPL##10##11##]. However, the relative levels of gene expression obtained from sqRT-PCR or Northern hybridizations were lower than those from microarray hybridizations. This inconsistency could have been attributed to several factors: the differences in the efficiencies of reverse transcriptase, the low copy of an mRNA transcript and the different priming methods or the fundamental differences in data normalization procedures between microarray and Northern hybridizations or sqRT-PCR. For example, in this study, the former was achieved through the global normalization, while the later used the expression of the housekeeping gene <italic>Actin </italic>to normalize the sqRT-PCR or Northern blots. So, application of normalization criteria might have an effect in the correlation of fold changes between these two methods [##REF##11161795##80##,##REF##14706121##81##]. Though the housekeeping gene <italic>Actin </italic>was considered to be constitutively expressed in tissues, due to the nature of the processing of hybridizing samples to Northern blots, one could reasonably assume a possible inconsistency resulting from these two methods. Certainly, other sensitive techniques such as quantitative RT-PCR (qRT-PCR) and nuclear run-on assays can be used further to validate the results generated from these expression studies.</p>
<title>Results and discussion</title> <title>Symptom development in PPV-infected <italic>Arabidopsis </italic>plants</title> <p>Susceptibility of <italic>Arabidopsis </italic>accession Col-0 to the infection of a PPV-D strain isolated in Canada was monitored by recording the development of symptoms and the virus titre. No typical viral symptoms such as mosaic, chlorosis, ringspot, vein yellowing, or leaf curling were visible in plants inoculated with PPV at 17 dpi and up to 42 dpi. This is consistent with the recent finding that <italic>Arabidopsis </italic>plants infected by two PPV-D isolates displayed mild or no symptoms [##REF##16673941##20##]. At 17 dpi, the most obvious phenotypes observed in this study were mild growth retardation, slightly delayed flowering and bolting (Figure ##FIG##0##1A##). These phenotypes were not documented for <italic>Arabidopsis </italic>accession Col-0 infected with two PPV-D isolates but were very similar to accession Bay-O inoculated with one of the two PPV-D isolates [##REF##16673941##20##]. Considering the diverse phenotypes of different accession-isolate combinations, the observed phenotypic variation between the two studies is likely due to the PPV isolates used. Using reverse transcription-polymerase chain reaction (RT-PCR) and enzyme-linked immunosorbent assay (ELISA), PPV was detected in newly developed rosette leaves above the inoculated leaves (Figure ##FIG##0##1B##) at 17 dpi. Since PPV-induced symptoms appeared and the virus was detected by RT-PCR and ELISA at 17 dpi in the infected leaf tissues, it was expected that changes in host gene expression would be altered at this time point. In addition, we also inoculated <italic>Nicotiana benthamiana </italic>plants with PPV. Infected leaves showed chlorosis and severe mosaic symptoms as early as 7 dpi (data not shown). At 17 dpi, the virus titre in <italic>N. benthamiana </italic>was at least 100 times or higher than that in <italic>Arabidopsis </italic>(data not shown), suggesting that <italic>Arabidopsis </italic>may have stronger basal or inducible resistance to PPV than <italic>N. benthamiana</italic>.</p> <title>Over- or under-represented functional categories of differentially expressed <italic>Arabidopsis </italic>genes in PPV-infected leaf tissues</title> <p>At 17 dpi, newly expanded PPV-infected rosette leaves from the inoculated plants and the corresponding leaves from the mock-inoculated plants were sampled for microarray hybridizations. Consistency of the microarray data was assessed by comparing the mean normalized signal intensities obtained from the mock-inoculated control samples across three independent biological replicates. A high correlation (0.96 – 0.98) among the mock inoculated control samples was observed, indicating low biological variability between the replicates.</p> <p>Gene expression data were analyzed by one-way ANOVA to identify differentially expressed genes in PPV-infected rosette leaves from the inoculated plants relative to leaves from the mock-inoculated plants. Using stringent selection criteria with the Benjamini and Hochberg false discovery rate (FDR) of 5% [##UREF##5##21##], corresponding to <italic>Q </italic>≤ 0.05 [##REF##12883005##22##], ~31.4% of the non-redundant (nr) genes (7,151 out of the 22,810 <italic>Arabidopsis </italic>genes printed on the chip) were identified as significantly differentially expressed (either induced or repressed) in response to PPV infection. These 7,151 genes were further filtered using a 2.5-fold increase or decrease in signal intensity. Using this fold change cut-off, we identified 2,013 and 1,457 <italic>Arabidopsis </italic>genes that were significantly induced and repressed by PPV in systemically infected tissues, respectively.</p> <p>Each identified gene was then assigned to a functional class according to the <italic>Arabidopsis </italic>MIPS (Munich Information Centre for Protein Sequences) functional classification scheme (Figure ##FIG##1##2##). With the information gathered from TAIR (The <italic>Arabidopsis </italic>Information Resource) database (please see Availability &amp; requirements for more information), we were able to assign a putative function for 54.3% of the total 22,810 genes on the array, while the remaining 45.7% were undefined as they were annotated as "unknowns or orphans" (Table ##TAB##0##1##). The expression levels of <italic>Arabidopsis </italic>genes significantly induced or repressed by PPV were assigned to 12 different major functional groups. In some major functional groups, genes were further separated into subcategories to simplify the data for biological interpretation [see Additional file ##SUPPL##0##1##]. The number of genes in each functional category was subjected to Fisher's exact test [##REF##14668247##23##,##REF##12620386##24##] to determine if genes involved in certain biological processes were over- or under-represented (<italic>P </italic>≤ 0.001). The functional groups significantly over-represented in the set of upregulated genes were carbohydrate/soluble sugar/starch/amino acid metabolism (<italic>P </italic>&lt; 7E-12), primary and secondary metabolism (<italic>P </italic>&lt; 2.2E-07), intracellular membrane/membrane-bound organelle (<italic>P </italic>&lt; 4.5E-06), clathrin binding (<italic>P </italic>&lt; 1.2E-05), glycolysis (<italic>P </italic>&lt; 2.0E-05), protein folding/heat shock/chaperone activity (<italic>P </italic>&lt; 5.2E-06), cell wall-associated transcripts (<italic>P </italic>&lt; 1.6E-05), RNA metabolism (<italic>P </italic>&lt; 1.9E-05), chloroplast (<italic>P </italic>&lt; 1.8E-06), protein fate (<italic>P </italic>&lt; 3.5E-07), ubiquitin-like conjugating enzyme activity (<italic>P </italic>&lt; 4.1E-05) and development/storage proteins (<italic>P </italic>&lt; 4.8E-08). Among the downregulated genes, the over-represented functional groups were metal ion binding/metal binding (<italic>P </italic>&lt; 1.6E-07), development/storage proteins (<italic>P </italic>&lt; 3.3E-06), cell wall-associated transcripts (<italic>P </italic>&lt; 0.00013) and protein synthesis and translation (<italic>P </italic>&lt; 1.4E-43). Genes involved in primary/secondary metabolism (<italic>P </italic>&lt; 1.4E-06) and intracellular membrane/membrane-bound organelle (<italic>P </italic>&lt; 1.3E-05) functional groups were significantly under-represented.</p> <p>The analysis showed that that the number of upregulated genes was more than that of downregulated genes in PPV-infected <italic>Arabidopsis </italic>leaves, probably reflecting more diverse genes induced by PPV infection. This notion is supported by a disproportionate allocation in the number of altered genes in the infected leaf tissues belonging to major functional groups of metabolism, transcription/splicing/RNA processing proteins, defence and development/storage proteins.</p> <title>Gene expression changes in metabolism, chloroplast and protein fate in PPV-infected <italic>Arabidopsis </italic>leaf tissues likely involved in symptom development</title> <p>Virus infections often disturb the general metabolism in the host. This includes the accumulation of soluble sugars (e.g., sucrose synthase) and starch (e.g., starch synthase), an increase in the rate of respiration, a decrease in the rate of photosynthesis and elevated levels of amino acid synthesis (e.g., amino acid permease and aminotransferase) [##UREF##6##25##, ####UREF##7##26##, ##REF##12226342##27##, ##UREF##8##28##, ##REF##17170466##29##, ##UREF##9##30####9##30##]. In this study, we found that in PPV-infected leaves, 23% (124/548) of genes associated with metabolism of starch, soluble sugars and amino acids and 2.7% (15/548) of genes related to glycolysis showed elevated expression levels [see Additional file ##SUPPL##0##1##]. Appearance of disease symptoms such as chlorosis on the <italic>Cucumber mosaic virus </italic>(CMV)-infected leaf tissues in cucumber was correlated with increased glycolysis, phosphoenol pyruvate phosphate, respiration and starch accumulation, decreased photosynthesis and reduction in total protein synthesis [##REF##12226342##27##]. The increased expression of these genes may result from the physical blockage of the transport path, the inactivation of sugar transport proteins or the induction of cell wall invertase [##UREF##10##31##, ####REF##10508760##32##, ##REF##12663220##33####12663220##33##]. Indeed, several cell wall invertase genes were also induced in PPV-infected leaves in this study. These data indicate that a series of intercrossing pathways involved in sugar metabolism was differentially regulated by PPV infection. It is possible that altered gene expression pathways involved in sugar metabolism might contribute to the mild symptom development observed in PPV-infected <italic>Arabidopsis </italic>leaves.</p> <p>In addition to soluble sugars, starch and amino acids, PPV also induced 16% (87/548) of membrane-associated proteins, including proteins associated with integral membrane, plasma membrane, endoplasmic reticulum (ER), vacuole, mitochondria and peroxisomes [see Additional File ##SUPPL##0##1##]. Recently, membrane proteins have also been found to be upregulated in <italic>N. benthamiana </italic>plants infected with <italic>Tomato ringspot virus </italic>(ToRSV) and PPV [##REF##17722703##13##]. For example, genes encoding proteins such as plasma membrane polypeptide, putative integral membrane, COP-coated vesicle membrane, clathrin protein and mitochondrial malate dehydrogenase were induced in PPV-infected <italic>Arabidopsis </italic>leaves that have similar functional identity to those upregulated in <italic>N. benthamiana </italic>plants infected with ToRSV and PPV. It is believed that many positive-strand RNA viruses replicate their genome in association with cellular membranes. Viral proteins that are required for the formation of viral replication complexes (VRCs) are often shown to be located in ER [##REF##17292435##34##, ####REF##10579716##35##, ##UREF##11##36####11##36##]. It would be interesting to investigate whether these upregulated membrane proteins are co-localized to VRCs. Characterization of such host membrane factors will provide useful information regarding molecular mechanisms governing membrane specificity for viral RNA replication.</p> <p>Of the upregulated genes coding for nuclear or cytosolic proteins in PPV-infected <italic>Arabidopsis </italic>leaves, two OGTs (O-linked N-acetyl glucosamine transferase), named SEC (SECRET AGENT) and SPY (SPINDLY, a gibberellin response protein), which can potentially posttranslationally modify the viral capsid protein (CP), are very interesting. Study by Chen et al. [##REF##16014901##37##] has shown that the N-terminus of the CP of PPV is modified with O-linked N-acetyl glucosamine in <italic>N. clevelandii</italic>, <italic>N. benthamiana </italic>and <italic>Arabidopsis</italic>, indicating that the modifications are not host specific. In <italic>Arabidopsis</italic>, this modification is catalyzed by SEC but not by SPY. Such modification is not essential but plays a role in the infection process [##REF##16014901##37##]. SEC and SPY have partial functional redundancy in plants [##REF##16014901##37##,##REF##12136030##38##]. It is not clear whether the upregulation of SPY in PPV-infected <italic>Arabidopsis </italic>is a compensation for the loss of SEC that is used for the modification of the PPV CP.</p> <p>We also observed as many as 52 chloroplast genes that were upregulated in PPV-infected leaf tissues [see Additional file ##SUPPL##0##1##]. These represent approximately 60% of the chloroplast genes in the <italic>Arabidopsis </italic>genome [##REF##10574454##39##]. Several studies suggest that differential regulation of chloroplast gene expression is correlated to viral symptoms such as mosaic and chlorosis [##REF##17722703##13##,##REF##17427806##16##]. Recently, Dardick [##REF##17722703##13##] has reported that ToRSV or PPV infection in <italic>N. benthamiana </italic>leaves induces severe mosaic and chlorosis symptoms. However, the infection of this host with <italic>Prunus necrotic ringspot virus </italic>(PNRSV) hardly induces any symptoms. Using potato cDNA chips, it has been shown that infections of ToRSV and PPV result in the suppression of 34% (83/243) of chloroplast-related genes, while no common genes are among the induced genes in ToRSV- and PPV-infected <italic>N. benthamiana</italic>. In contrast to other two viruses, PNRSV infection downregulates 67% (10/15) and upregulates 33% (5/15) of the chloroplast-related genes [##REF##17722703##13##]. Analogous to PNRSV-infected <italic>N. benthamiana</italic>, PPV-infected <italic>Arabidopsis </italic>leaves did not show typical viral symptoms such as mosaic and chlorosis. The induction of such a large number of chloroplast-related genes in <italic>Arabidopsis </italic>may help maintain the normal function and morphology of the infected leaf and counteract possible detrimental effects resulting from PPV infection.</p> <p>In virus infected plants, protein activities are often altered through posttranslational modifications such as ubiquitination [##REF##17417941##4##]. In the current study, about 9% (60 genes) of the upregulated protein fate genes encoding proteins associated with the 26S proteasome were very significant (<italic>Q </italic>≤ 0.05) and over-represented (<italic>P </italic>&lt; 3.5E-07). Of the seven upregulated 26S proteasome genes, two genes code for AAA-ATPase subunits including RPT3 that reside in the mature 26S proteasome [##REF##17427806##16##] and five genes including <italic>Rpn1</italic>, <italic>3 </italic>and <italic>6 </italic>are involved in ubiquitination. Fundamentally, the 26S proteasome is known to control many biochemical processes by programmed degradation of regulatory protein targets [##REF##16905670##40##]. Upregulation of these genes may interfere with ubiquitin-dependent reactions in plants [##REF##11473704##41##]. It has been reported that, the 26S proteasome degrades the virus-encoded movement protein that facilitates the viral cell-to-cell spread [##REF##10708450##42##]. Therefore, the upregulation of the 26S proteasome genes may be a resistance response to PPV infection.</p> <p>On the other hand, in <italic>N. benthamiana </italic>plants infected with ToRSV, about 72% (33/46) of the proteasome genes have been shown to be upregulated [##REF##17722703##13##]. A few archetype proteasome proteins such as subtilisin-like protease, aspartyl protease, ubiquitin conjugating enzyme, F-box protein and 26S proteasome regulatory subunits induced by PPV in <italic>Arabidopsis </italic>leaves are upregulated by ToRSV in <italic>N. benthamiana </italic>plants [##REF##17722703##13##]. Upregulation of such protein fate and ubiquitin genes has also been observed in other studies during viral infection [##REF##15604696##10##,##REF##17427806##16##,##REF##8986804##43##,##REF##11041891##44##]. Study by Marathe <italic>et al</italic>. [##REF##15604696##10##] has demonstrated that the induction of F-box proteins may mediate the degradation components of the CMV-Y resistance pathway. Though the above mentioned protease genes induced by PPV infection are known to function in the protein degradation, plant defence responses and hypersensitive reactions [##REF##11976458##45##, ####REF##11896264##46##, ##REF##8843949##47####8843949##47##], their exact roles in PPV infection are to be elucidated.</p> <title>Association of the repression of development/storage proteins, protein synthesis and cell wall-related genes in <italic>Arabidopsis </italic>leaf tissues with PPV infection and symptom development</title> <p>Many virus-induced symptoms have been linked to alterations in plant hormone synthesis [##REF##17417941##4##,##REF##15681455##48##,##REF##16903352##49##]. Of the 96 development/storage genes identified in this study, 10 auxin, 5 ethylene and 4 gibberellin-related genes were repressed in PPV-infected leaves [see Additional file ##SUPPL##1##2##]. Recent studies [##REF##15681455##48##,##REF##16903352##49##] on TMV infection in <italic>Arabidopsis </italic>have shown that during the process of infection the TMV replicase protein interacts with a subset of Aux/IAA proteins, resulting in the activation of auxin responsive transcription factors (ARF). Such an alteration in the transcription levels of auxin-responsive genes and the disruption of host normal physiological process have shown to be associated with disease symptoms. In this study, PPV infection suppressed the expression of two such genes encoding auxin IAA proteins (At5g57420 and At1g52830). With respect to the downregulation of ethylene responsive genes, though there is no direct evidence suggesting interactions between virus and ethylene signalling components, specific viral proteins have been suggested to alter the ethylene signalling pathway [##REF##17417941##4##]. For instance, the P6 protein of <italic>Cauliflower mosaic virus </italic>(CaMV) is associated with the disruption of the ethylene response pathway and this may play a key role in the induction of symptoms such as mild vein chlorosis to severe chlorosis and stunting [##REF##9390424##50##]. In the case of downregulation of gibberellic acid (GA), a direct interaction between the <italic>Rice dwarf virus </italic>(RDV) P2 protein and the rice ent-kaurene oxidase, a key enzyme in the synthesis of GA in plants was found to be responsible for the significant reduction of GA1 in RDV-infected plants, and in the appearance of dwarf phenotype [##REF##16299167##51##]. Collectively these examples demonstrate that the suppression of hormone-related genes by viruses induces symptoms through altering the host hormone and developmental signalling pathway. It would be interesting to know whether any of the observed downregulated hormone-related genes is associated with the reduced growth and delayed development of PPV-infected <italic>Arabidopsis</italic>.</p> <p>Another significantly over-represented subset of genes suppressed in PPV-infected tissues includes genes coding for the protein synthesis and translation machinery [see Additional file ##SUPPL##1##2##]. This functional category contains 33 ribosomal subunits: 21 large, 11 small and 1 plastid specific ribosomal proteins. In the <italic>Arabidopsis </italic>genome, there are ~227 ribosomal genes [##REF##17722703##13##,##REF##14681386##52##]. About 15% of these ribosomal genes were printed on the array. Based on the analysis, we observed that approximately 81% (27/33) of the small and large ribosomal subunit genes were significantly repressed in PPV-infected leaves.</p> <p>However, an opposite trend was observed in PPV-infected <italic>N. benthamiana </italic>plants [##REF##17722703##13##], where 80% (105/131) of potato ribosomal genes were upregulated. Comparison of these potato upregulated ribosomal genes in <italic>Arabidopsis </italic>through BLAST (E-value &lt; -10) sequence similarity searches identified 55 nuclear ribosomal subunits [##REF##17722703##13##]. Comparative analysis of 27 <italic>Arabidopsis </italic>nuclear ribosomal subunits that were downregulated by PPV in our study with those in Dardick's study (using Table S5 from Dardick's study for comparison) showed 8 (7 large and 1 small) <italic>Arabidopsis </italic>ribosomal subunits with an opposing trend. Since plant viruses lack protein synthesis ability, they rely on the host cell translation machinery to produce viral proteins. Thus, PPV might benefit from the upregulation of ribosomal genes in <italic>N. benthamiana </italic>which showed more susceptibility to PPV and displayed much more severe symptoms than <italic>Arabidopsis</italic>. However, at present there is no evidence to discern whether the observed downregulation of ribosomal genes is due to the stronger basal resistance of <italic>Arabidopsis </italic>to PPV infection or due to the natural intricacy of the different host responses to PPV that are diverse depending on the tissue type, time post infection, temperature, host and virus strain, or simply due to an artifact of hybridizations, given that two different platforms were used in the studies.</p> <p>Cell wall-related genes are major determinants of cell morphogenesis in plants, encoding several hundreds of different structural proteins and cell wall-related enzymes [##REF##11554474##53##,##REF##16774842##54##]. In this study, 75 cell wall-related genes encoding protein for structure (23 genes), reassembly (34 genes), expansins (9 genes) and matrix polymers (9 genes) were significantly downregulated in PPV-infected leaves [see Additional file ##SUPPL##1##2##]. This is consistent with results observed in <italic>Arabidopsis </italic>plants infected by TuMV [##REF##17427806##16##] and in rice plants infected with RDV [##REF##17378427##15##]. In these studies, repression of cell wall related genes is correlated with symptom development. Of these downregulated genes, extensin is a cell wall-localized hydroxyproline-rich glycoprotein that usually forms a cross-linked network with pectin to create a highly impassable barrier against pathogens [##REF##7659756##55##]. Downregulation of extensin during PPV infection suggests a decreased extensin deposition [##UREF##12##56##]. Conceivably, PPV exploits the plasmodesmata to move from initially infected cells to the neighbouring healthy cells. It is not likely that PPV cell-to-cell movement is related to the deposition of extensin. Nonetheless, the molecular mechanisms and the biological connotation underlying the downregulation of cell wall-related genes in the infected leaf remain to be understood.</p> <p>In addition to the above described genes, many genes in the category of metal ion binding or metal binding were also downregulated. It is not clear if the downregulation of these genes can lead to physiological changes that are associated with PPV infection.</p> <title>Identification of genes differentially regulated by PPV in infected <italic>Arabidopsis </italic>protoplasts</title> <p>Resistance induced through signal transduction pathways in <italic>Arabidopsis </italic>may take place immediately after PPV infection. The <italic>Arabidopsis </italic>leaves infected with PPV at 17 dpi contained a mixture of uninfected, infected, and post-infected cells. Thus, the spatial and temporal information on gene expression was lost. Sampling at earlier time points such as 1 or 2 dpi can reduce the percentage of post-infected cells but will largely increase the percentage of non-infected cells. So, protoplasts have been suggested to overcome this barrier [##UREF##3##17##]. There are several advantages in using <italic>Arabidopsis </italic>protoplasts to study plant-virus interactions. First, in contrary to unpredictable infections on the leaf, synchronized protoplast cells allows for a more precise determination of gene expression. Second, a much higher proportion of the cells are infected than that in infected leaf tissues [##UREF##13##57##]. For example, we have shown that in <italic>Arabidopsis </italic>protoplasts transfected with plasmid pPPV-SK68, a PPV infectious clone, the transfection efficiency is approximately 35% [##REF##16777241##58##]. Hence, the protoplast system makes it possible to study molecular events as early as just a few hours post infection. Like other experimental systems, the protoplast system also has some limitations. Protoplasts lack the cell wall so the physiological integrity of the plant organ is disrupted. During isolation, transfection and culture of protoplasts, gene expression changes induced by harsh treatments such as osmotic stress, enzymatic digestion and centrifugation may mask those induced by virus infection [##UREF##3##17##]. Physiological milieu that a virus encounters in protoplast cells might be quite different from that of an intact leaf cell. Sample bias may result from the minute amounts of RNA isolated from the protoplasts and from the subsequent amplification step. However, improved technologies and experimental designs may be used to minimize the negative impact of the system [##REF##7731989##59##, ####REF##3409786##60##, ##REF##16027030##61##, ##UREF##14##62####14##62##]. Recently, we reported a protocol describing the isolation and transfection of <italic>Arabidopsis </italic>protoplasts for studying PPV-host interactions [##REF##16777241##58##]. In order to identify early time points for profiling global gene expression changes in infected protoplast cells, a time course experiment was performed. <italic>Arabidopsis </italic>protoplasts were transfected with a PPV infectious cDNA clone, pPPV-SK68 derived from an M strain and with a non-infectious deletion mutant, pPPV-SK68Δ as control, respectively [##REF##16777241##58##]. In the mutant clone, the deletion abolishes the translation of viral proteins CI, 6K2 and NIa/VPg required for PPV replication. The quantitative results showed that the viral RNA in the protoplasts was detectable at 3 hpt, accumulated at 6 hpt and peaked at 12 hpt (Figure ##FIG##0##1C##), followed by a decline at 24 hpt [##REF##16777241##58##]. The time course progression of PPV infection in <italic>Arabidopsis </italic>protoplasts is similar to that of <italic>Cymbidium mosaic virus </italic>(CyMV) in orchid protoplasts [##REF##8249447##63##]. Thus, three time points, 3, 6 and 12 hpt were selected for this study.</p> <p>Gene expression in PPV-infected <italic>Arabidopsis </italic>protoplasts was analyzed by one-way ANOVA to identify genes differentially expressed at each time point. The degree of biological variability and reproducibility of the microarray data was assessed by comparing the mean values of the normalized signal intensity from the protoplasts transfected with a non-infectious mutant, pPPV-SK68Δ with all three independent biological replicates at each time point. The analysis resulted with a high correlation coefficient of 0.90 – 0.95, indicating low technical and biological variability between the samples used in the hybridizations. Using FDR at 5%, we identified ~8.4% (1,907 out of 22,810) of non-redundant genes that were significantly (<italic>Q </italic>≤ 0.05) differentially regulated in protoplasts transfected with the PPV infectious clone, pPPV-SK68 at three different time points relative to the control transfected with the non-infectious deletion mutant, pPPV-SK68Δ. These genes were further filtered using a 2.5-fold increase or decrease in signal intensity, resulting in 263 nr genes that were induced (≥ 2.5-fold), and 304 genes that were repressed (≤ -2.5-fold), as illustrated in the Venn diagrams (Figure ##FIG##2##3A## and ##FIG##2##3B##). A comparison analysis between 263 induced and 304 suppressed genes identified 107 genes induced at one time point and were not repressed at other time points, 148 genes repressed at one time point and were not induced at other time points, and 156 genes induced and repressed at different time points (Figure ##FIG##2##3C##) [see Additional file ##SUPPL##2##3##]. The analysis also revealed that the number of repressed genes was greater than that of the induced genes at all three time points analyzed. The number of genes that were commonly repressed among the time points was 14 between 3 and 12 hpt, 10 between 3 and 6 hpt, and 12 between 6 and 12 hpt. We further examined each of these 263 induced and 304 repressed genes to see if they were also consistently induced or repressed at all three different time points of infection. Fifteen genes were found to be constantly repressed and 4 were consistently induced in the PPV transfected protoplasts, suggesting that the proportion of genes repressed at all three time points were greater than that of the genes induced.</p> <p>The genes induced or repressed in PPV-infected protoplasts were annotated, assigned to a predicted function and classified into 12 different major functional categories using the <italic>Arabidopsis </italic>MIPS functional classification scheme (Figure ##FIG##3##4##). Fisher's exact test [##REF##14668247##23##,##REF##12620386##24##] was used to determine whether the representation of the number of genes in each functional category was significant (<italic>P </italic>≤ 0.001) or by random chance. Of the upregulated genes, transcripts involved in primary and secondary metabolism (<italic>P </italic>&lt; 6.2E-06), cellular signalling (<italic>P </italic>&lt; 4.4E-11), defence response (<italic>P </italic>&lt; 2.7E-06), transcription factors (<italic>P </italic>&lt; 1.1E-11) and transporters (<italic>P </italic>&lt; 1E-12) were over-represented, while in the case of downregulation, genes involved in transcription factors (<italic>P </italic>&lt; 5E-09), primary and secondary metabolism (<italic>P </italic>&lt; 6.4E-05), kinases or cellular signalling (<italic>P </italic>&lt; 1.3E-05), and defence responses (<italic>P </italic>&lt; 5.6E-05) were over-represented. A complete list of induced or repressed genes in PPV-infected protoplasts in different functional groups is provided [see Additional file ##SUPPL##2##3##].</p> <p>As evident [see Additional file ##SUPPL##3##4##; Additional file ##SUPPL##4##5##; Additional file ##SUPPL##5##6##], comparison of the differentially expressed genes between PPV-infected <italic>Arabidopsis </italic>protoplasts (263 up- and 304 downregulated genes) and PPV-infected leaf tissues (2,013 induced genes and 1,457 suppressed genes) resulted in only 32 and 33 common genes that were significantly induced (≥ 2.5-fold) and repressed (≤ -2.5-fold), respectively. A few vital genes such as genes encoding heavy metal protein, disease resistance protein, FtsH protease, cytochrome P450 and β-1, 3-glucanase were significantly differentially regulated by PPV in the infected leaves and transfected protoplasts. The smaller percentage of significantly differentially regulated genes observed in the PPV-infected protoplasts than in the PPV-infected leaf suggests a remarkable gene expression differences between the two systems. However, as discussed above, the differences in the alteration of gene expression changes in the protoplasts induced by PPV infection could have been masked by harsh treatments during protoplast isolation, transfection and culture or might result from the two different PPV strains used. Though the PPV-D and -M strains used in this study share more than 95% sequence identity [##UREF##15##64##], these two strains may induce different symptoms on the same host [##REF##16673941##20##,##UREF##16##65##].</p> <p>To visualize gene expression profiling at all three time points, a hierarchical clustering analysis was carried out on the 263 induced and 304 repressed genes (Figure ##FIG##4##5A##). A gene tree heat map was built using Pearson correlation distance metric with an average linkage algorithm. To understand further the temporal relationship of the host gene expression in PPV-infected protoplasts at each time point, twelve distinct groups of expression pattern clusters were generated on the significant differential gene expression datasets using K-means clustering with Pearson correlation distance metric (Figure ##FIG##4##5B##) [see Additional file ##SUPPL##6##7##]. A complete list of genes induced or repressed at each time point in different cluster groups along with their expression levels and putative functions is provided [see Additional file ##SUPPL##7##8##]. The clustering analysis also revealed greater changes in gene expression (236 genes) at 6 hpt than at 3 hpt (163 genes) and at 12 hpt (201 genes) in PPV-infected <italic>Arabidopsis </italic>protoplasts. These observations suggest that early gene expressions in response to PPV infection are active and dynamic.</p> <title>Early upregulated genes related to disease resistance or defence responses in PPV-infected <italic>Arabidopsis </italic>protoplasts</title> <p>PPV infection in protoplasts induced the expression of 263 genes. In contrast to about 11% of the suppressed genes that were common at different time points, there were no induced genes shared at different time points, suggesting induction of host gene expression was much more transient than suppression. Genes related to defence, cellular signalling, primary and secondary metabolism, transcription and transporters were upregulated and over-represented. Within these five categories, ~52% (48/93) of genes were predominantly upregulated (≥ 2.5-fold) at 6 hpt [see Additional file ##SUPPL##2##3##]. Examples of the upregulated genes in the above functional categories include MADS box protein, zinc finger family protein, WRKY family transcription factor, avirulence (Avr)-induced protein, hydrophobic protein, disease resistance TIR-NBS-LRR (toll-interleukin-1-receptor/nucleotide-binding site/leucine-rich repeat), protein kinases and phosphatases.</p> <p>Another interesting gene upregulated by PPV encodes a disease resistance protein, TIR-NBS-LRR. This gene was upregulated with 6.39-fold in PPV-infected protoplasts at 6 hpt and 2.12-fold at 12 hpt, and 10.84-fold in PPV-infected leaves [see Additional file ##SUPPL##3##4##]. In the <italic>Arabidopsis </italic>(Col-0 accession) genome, there are 149 NBS-LRR-like genes [##UREF##17##66##,##REF##12671079##67##]. The resistance gene <italic>NBS-LRR </italic>mostly confers resistance to viral diseases but not necessarily associates with cell death or tissue necrosis [##UREF##17##66##]. Genes <italic>Sw5</italic>, <italic>N </italic>and <italic>Rx </italic>are homologous to the <italic>NBS-LRR </italic>gene. <italic>Rx </italic>resistance to <italic>Potato virus X </italic>(PVX; Potexvirus) is the most classical example, which shows an extreme resistance that inhibits virus replication without hypersensitive cell death [##REF##10330465##68##]. The <italic>Sw5 </italic>gene confers resistance to several Tospoviruses [##REF##11043474##69##], and <italic>N </italic>and <italic>Rx </italic>are effective against many natural variants of TMV and PVX [##UREF##18##70##,##UREF##19##71##]. Since the specificity determinants subsist in LRR domains [##UREF##17##66##], it would be interesting to find out if the NBS-LRR gene detected in the PPV-infected leaf plays a role in resistance responses to PPV infection.</p> <p>The comparison of induced genes at different time points in PPV-infected protoplasts revealed that genes upregulated at 3 hpt might regulate their target genes at 12 hpt. For instance, two genes coding for cnd41 (At5g24820), a plastid specific DNA binding factor, and calmodulin (At3g15050, At3g16490), a signalling binding protein were induced at 3 hpt and their target genes encoding metabolic enzymes such as serine carboxypeptidase (At1g11080) were upregulated at 12 hpt [##UREF##13##57##,##REF##16777241##58##]. Similarly, transcription factors such as zinc finger proteins (At5g54020; At3g45480; At4g12140) was upregulated at 3 hpt and the genes encoding enzymes such as GDSL-motif lipase (At3g43550) and glycosyl hydrolases (At3g26140, At4g39010) were upregulated at 12 hpt [##REF##9426602##72##]. These results suggest that virus infection may induce expression of specific upstream regulator genes first which then activate the expression of downstream target genes.</p> <title>Inhibition of host mRNA accumulation in PPV-infected protoplasts</title> <p>Based on the expression analysis (Figure ##FIG##2##3##), the number of repressed genes in the transfected protoplasts was greater than that of the genes induced at all time points, suggesting that at early infection stages, suppression is more active than induction of host gene expression. Moreover, about 11% (33 out of 304) of the suppressed genes overlapped at different time points, indicating a relatively stable suppression of genes. Consistent with this, in pea cotyledons infected with <italic>Pea seed borne mosaic virus </italic>(PSbMV), suppression of diverse host genes including seed storage proteins (e.g., vicilin and convicilin), and starch synthase occurs only in a zone of six to eight cells immediately behind the infection front, which coincided with the onset of viral replication [##REF##17791345##73##]. In the present study, a gene coding for the seed storage protein cupin family (At2g18540) that shares sequence homology to the pea vicilin and convicilin was significantly repressed by PPV infection at 12 hpt in <italic>Arabidopsis </italic>protoplasts. Suppression of such host gene expression or "host shutoff" has been suggested to be due to stress responses, induced defence states or a mechanism employed by viruses to outcompete against host mRNAs for host resources for their translation and replication [##REF##17427806##16##,##UREF##3##17##,##REF##8986804##43##,##REF##11104816##74##]. Study by Rottier [##UREF##20##75##] has demonstrated that protein synthesis in protoplasts is blocked in cowpea protoplasts transfected with <italic>Cowpea mosaic virus</italic>. Despite this inhibition, viral mRNAs are preferentially translated. Since both viruses and hosts share the same translational apparatus, viruses must have a mechanism to avoid this inhibition. It is not clear whether the disappearance of host mRNAs is the consequence of mRNA destabilization triggered by stress response, or host transcription inhibition by the virus, or targeted degradation through a virus-induced ribonuclease activity.</p> <title>Identification of a common set of genes in general stress and defence-related groups differentially regulated by infections of PPV and other positive sense RNA viruses</title> <p>We performed a cross sequence comparison analysis by blasting individual sequence of all the genes significantly differentially regulated by PPV in infected <italic>Arabidopsis </italic>leaves against 2,156 uniEST sequences derived from PPV-infected <italic>P. persica </italic>leaves [##UREF##15##64##]. High-scoring pair (HSP; ≥ 100 score) and E-value (E-value ≤ 10–20) were used to identify orthologs. A total of 170 induced [see Additional File ##SUPPL##8##9##] and 18 repressed <italic>Arabidopsis </italic>genes showed high sequence similarity to <italic>P. persica </italic>uniESTs. Of these 170 induced genes, approximately ~39% (67/170) and ~16% (28/170) of the genes were involved in metabolism and defence, respectively. The defence-related genes include those coding for pathogenesis related (PR)-proteins, heat shock and germin-like proteins (Ger2) [see Additional file ##SUPPL##8##9##]. An increased expression of PR proteins requires salicylic acid (SA) signifying that compatible host-virus interactions induce defence-like responses through a SA dependent pathway [##REF##15728340##76##]. More recently, it has been shown that PPV infection in the apoplastic space of <italic>P. persica </italic>leaves induces a thaumatin-like protein (PR-5) [##REF##17043083##77##]. Induction of such protein is thought to be mediated by an increase in the level of peroxidase accumulation [##REF##17043083##77##]. A similar trend was also observed in PPV infected <italic>Arabidopsis </italic>leaves where genes encoding PR-5 and peroxidase were significantly induced and their counterparts were also found in the uniEST sequences derived from PPV-infected <italic>P. persica</italic>.</p> <p>To determine whether the diverse RNA viruses have the ability to elicit common gene expression changes, we cross-compared the genes that were significantly differentially regulated by PPV infection in this study [Additional file ##SUPPL##0##1##] [Additional file ##SUPPL##1##2##] to those <italic>Arabidopsis </italic>genes that were differentially regulated by infections of other positive sense RNA viruses including CMV [##REF##15111722##9##,##REF##15604696##10##]; TMV [##REF##12906112##7##]; TuMV, and <italic>Oil seed rape mosaic virus </italic>(ORMV) [##REF##17427806##16##], and to five diverse RNA viruses, i.e., TuMV, ORMV, <italic>Turnip vein clearing virus </italic>(TVCV), <italic>Potato virus X </italic>(PVX), and CMV [##REF##12535341##8##]. Collectively, 416 common genes [see Additional file ##SUPPL##9##10##] were identified from this analysis. These genes were further classified into 6 different categories. Group A contains 202 genes induced by infections of PPV and other positive sense RNA viruses. About ~40% (80/202) are metabolic genes and ~8% (17/202) are defence-related genes. Group B includes 41 <italic>Arabidopsis </italic>genes that were induced during PPV infection but repressed by infections of other positive sense RNA viruses. About ~27% (11/41) are involved in metabolism. Group C comprises of 119 <italic>Arabidopsis </italic>genes repressed during PPV infection but induced by infections of other positive sense RNA viruses. For example, a defence-related gene encoding glycosyl transferase (At1g70090), which was repressed in PPV-infected leaves at 17 dpi was observed to be induced more than 1.5-fold in TMV-infected leaves [##REF##12906112##7##] at 14 dpi. Group D consists of 7 <italic>Arabidopsis </italic>genes repressed by infections of PPV and other positive sense RNA viruses. Group E comprises 25 genes that were induced during PPV infection but either induced or suppressed by infections of other positive sense RNA viruses. Group F contains 22 <italic>Arabidopsis </italic>genes repressed during PPV infection but induced during infections of more than one positive sense RNA virus. Previous studies have revealed that in virus-infected <italic>Arabidopsis </italic>and tomato leaves as many as one-third of the genes commonly induced by viruses are known or predicted to be involved in plant defence- and stress-related genes [##REF##15231257##5##,##REF##17378427##15##]. In the current study, about 52% (217/416) of the functional proteins differentially regulated in response to infections of PPV and other positive sense RNA viruses are involved in defence, cellular signalling and metabolism. Twenty-nine cellular signalling genes involved in defence mechanisms including protein kinase (At3g08760, At4g32250), protein phosphatase 2C (At5g02760, At3g05640), pyruvate kinase (At2g36580), serine/threonine protein kinase (At3g17410) and calmodulin-binding family protein (At5g62390) [##UREF##0##3##,##UREF##21##78##,##REF##16838784##79##] were induced by infections of PPV and other positive sense RNA viruses. Substantial induction of these defence, cellular signalling and metabolism genes by PPV and other positive sense RNA viruses in infected leaf tissues suggests that there is a shared pathway modulating their expression.</p> <title>Validation of microarray data by semi-quantitative RT-PCR (sqRT-PCR) and Northern hybridizations</title> <p>Since most changes in gene expression occurred at 6 hpt, 9 transcripts significantly upregulated and 1 transcript significantly downregulated at 6 hpt were selected to validate the microarray data using sqRT-PCR. We also selected 3 other transcripts that were significantly upregulated at 3 or 12 hpt for validating the microarray data through sqRT-PCR. We further confirmed the microarray data of 5 <italic>Arabidopsis </italic>genes that were differentially regulated in both PPV-infected protoplasts and leaves using sqRT-PCR or Northern hybridizations. Overall, the expression patterns of all the 13 genes analyzed by sqRT-PCR were consistent with those by microarray hybridizations (Figures ##FIG##5##6## and ##FIG##6##7##) [see Additional file ##SUPPL##10##11##]. However, the relative levels of gene expression obtained from sqRT-PCR or Northern hybridizations were lower than those from microarray hybridizations. This inconsistency could have been attributed to several factors: the differences in the efficiencies of reverse transcriptase, the low copy of an mRNA transcript and the different priming methods or the fundamental differences in data normalization procedures between microarray and Northern hybridizations or sqRT-PCR. For example, in this study, the former was achieved through the global normalization, while the later used the expression of the housekeeping gene <italic>Actin </italic>to normalize the sqRT-PCR or Northern blots. So, application of normalization criteria might have an effect in the correlation of fold changes between these two methods [##REF##11161795##80##,##REF##14706121##81##]. Though the housekeeping gene <italic>Actin </italic>was considered to be constitutively expressed in tissues, due to the nature of the processing of hybridizing samples to Northern blots, one could reasonably assume a possible inconsistency resulting from these two methods. Certainly, other sensitive techniques such as quantitative RT-PCR (qRT-PCR) and nuclear run-on assays can be used further to validate the results generated from these expression studies.</p>
<title>Conclusion</title> <p>Microarray analyses conducted in this study allowed us to gain insights into global gene expression changes in <italic>Arabidopsis </italic>associated with PPV infection, symptom development and defence responses. At 17 dpi, PPV displayed mild growth retardation, slightly delayed flowering and bolting. PPV infection was associated with the downregulation of genes related to development/storage proteins, protein synthesis and translation and cell wall-related genes and the upregulation of genes involved in soluble sugars, starch, amino acid, intracellular membrane/membrane-bound organelles, chloroplast and protein fate. Transcriptional profiling of PPV-infected protoplasts revealed that at early infection stages, downregulation of host gene expression was more prominent than upregulation. Concomitant with this downregulation, several defence and cellular signalling genes were induced during PPV infection. These genes may be involved in signal transduction pathways to induce resistance to combat the virus. Further targeted functional studies of these early defence-related genes will help unravel their biological function during the course of PPV infection.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Virus infection induces the activation and suppression of global gene expression in the host. Profiling gene expression changes in the host may provide insights into the molecular mechanisms that underlie host physiological and phenotypic responses to virus infection. In this study, the <italic>Arabidopsis </italic>Affymetrix ATH1 array was used to assess global gene expression changes in <italic>Arabidopsis thaliana </italic>plants infected with <italic>Plum pox virus </italic>(PPV). To identify early genes in response to PPV infection, an <italic>Arabidopsis </italic>synchronized single-cell transformation system was developed. <italic>Arabidopsis </italic>protoplasts were transfected with a PPV infectious clone and global gene expression changes in the transfected protoplasts were profiled.</p> <title>Results</title> <p>Microarray analysis of PPV-infected <italic>Arabidopsis </italic>leaf tissues identified 2013 and 1457 genes that were significantly (<italic>Q </italic>≤ 0.05) up- (≥ 2.5 fold) and downregulated (≤ -2.5 fold), respectively. Genes associated with soluble sugar, starch and amino acid, intracellular membrane/membrane-bound organelles, chloroplast, and protein fate were upregulated, while genes related to development/storage proteins, protein synthesis and translation, and cell wall-associated components were downregulated. These gene expression changes were associated with PPV infection and symptom development. Further transcriptional profiling of protoplasts transfected with a PPV infectious clone revealed the upregulation of defence and cellular signalling genes as early as 6 hours post transfection. A cross sequence comparison analysis of genes differentially regulated by PPV-infected <italic>Arabidopsis </italic>leaves against uniEST sequences derived from PPV-infected leaves of <italic>Prunus persica</italic>, a natural host of PPV, identified orthologs related to defence, metabolism and protein synthesis. The cross comparison of genes differentially regulated by PPV infection and by the infections of other positive sense RNA viruses revealed a common set of 416 genes. These identified genes, particularly the early responsive genes, may be critical in virus infection.</p> <title>Conclusion</title> <p>Gene expression changes in PPV-infected <italic>Arabidopsis </italic>are the molecular basis of stress and defence-like responses, PPV pathogenesis and symptom development. The differentially regulated genes, particularly the early responsive genes, and a common set of genes regulated by infections of PPV and other positive sense RNA viruses identified in this study are candidates suitable for further functional characterization to shed lights on molecular virus-host interactions.</p>
<title>Availability &amp; requirements</title> <p>The <italic>Arabidopsis </italic>Information Resource: <ext-link ext-link-type="uri" xlink:href="http://www.arabidopsis.org/"/></p> <p>McGill University and Genome Quebec Innovation Centre: <ext-link ext-link-type="uri" xlink:href="http://genomequebec.mcgill.ca"/></p> <p>Munich Information Centre for Protein Sequences: <ext-link ext-link-type="uri" xlink:href="http://mips.gsf.de/proj/thal/db/"/></p> <title>Authors' contributions</title> <p>MB carried out most of the experimental work including design of experiments, preparation and transfection of protoplasts, RNA isolation, analysis and interpretation of microarray data, and verification experiments. JSG participated in preparation and transfection of protoplasts, and isolation of RNA. T–SH prepared and inoculated <italic>Arabidopsis </italic>plants. AW conceived the work and interpreted the data. MB and AW wrote the paper. All authors have read and approved the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank Drs. Laszlo Palkovics and Ervin Balazs (Agricultural Biotechnology Center, Hungary) for providing the PPV infectious clone; Drs. Hélène Sanfaçon, Brian Miki and Dan Brown (Agriculture and Agri-Food Canada) for very helpful discussion; Drs. Antonet Svircev and Lorne Stobbs (AAFC) for assistance in live virus work; and Dr. Daowen Wang (Chinese Academy of Sciences, China) and Dr. Hélène Sanfaçon for critical reading of the manuscript. This work was supported by grants to AW from AAFC and from the Natural Science and Engineering Research Council of Canada.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Systemic infection and detection of PPV in <italic>Arabidopsis thaliana </italic>plants</bold>. (A) Symptoms on mock-inoculated (left) and PPV-infected (right) plants 17 days post inoculation (dpi). (B) RT-PCR and ELISA analysis of PPV-infected and mock-inoculated leaf tissues at 17 dpi. BR1, BR2, BR3, represents three independent biological replicates of PPV-infected and mock-inoculated leaf tissues, respectively. RT-PCR amplification of a 467 bp fragment corresponding to a segment of the PPV genome (nt 9140 to 9607), as tested using 25 cycles of amplification. +, ELISA positive; -; ELISA negative. (C) sqRT-PCR amplification of a cDNA fragment of the PPV genome isolated from <italic>Arabidopsis </italic>protoplasts transfected with pPPV-SK68 and pPPV-SK68Δ at 3, 6 and 12 hours post transfection (hpt), respectively. Histone 3 gene was used as a loading control. RT-PCR amplification for both panel B and C were carried out using the same RNA samples that were used in microarray hybridizations.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Functional distribution of <italic>Arabidopsis </italic>genes significantly induced and repressed in PPV-infected leaves</bold>. The genes were grouped following the methods of the MIPS <ext-link ext-link-type="uri" xlink:href="http://mips.gsf.de/proj/thal/db/"/><italic>Arabidopsis </italic>classification scheme. Genes whose function has not been determined were grouped under "unknown function". Number of genes identified in each functional group is indicated on the x-axis. Genes that belong to major functional categories are highlighted in bold and the subcategories within a major functional category are highlighted in italics.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Venn diagrams depicting the distribution of induced (≥ 2.5 fold) and repressed (≤ -2.5 fold) genes in PPV-infected protoplasts at three different time points</bold>. A statistical cut-off with a FDR at 5% that corresponds to <italic>Q </italic>≤ 0.05, after Benjamini and Hochberg's correction was used to determine genes significantly differentially regulated in protoplasts transfected with a PPV infectious cDNA clone, pPPV-SK68 at 3, 6 and 12 hours post transfection (hpt). The number of genes in the non-overlapping sector represents unique significant genes at each time point, while the overlapping sectors represent genes that are in common at time points indicated. The number in parentheses in the Venn diagram circles corresponds to the genes induced or repressed in the PPV infected leaves.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Functional distribution of <italic>Arabidopsis </italic>genes significantly induced and repressed in PPV-infected protoplasts</bold>. The genes were grouped following the method of the MIPS <ext-link ext-link-type="uri" xlink:href="http://mips.gsf.de/proj/thal/db/"/><italic>Arabidopsis </italic>classification scheme. Genes whose function has not been determined were grouped under "unknown function". Number of genes identified in each functional group is indicated on the x-axis. Genes that belong to major functional categories are highlighted in bold and the subcategories within a major functional category are highlighted in italics.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Clustering analysis of differentially regulatedgenes in PPV-infected <italic>Arabidopsis </italic>protoplasts</bold>. (A) Hierarchical clustering and changes in gene expression of 411 significantly (<italic>Q </italic>≤ 0.05) differentially regulated genes in PPV-infected <italic>Arabidopsis </italic>protoplasts at 3, 6 and 12 hours post transfection (hpt). Gene tree map generated from 411 significantly differentially regulated <italic>Arabidopsis </italic>genes are derived from 263 upregulated and 304 downregulated genes. (B) Using k-mean clustering, the gene expression profiles were grouped into twelve major cluster groups. The expression pattern of a gene in a cluster group is indicated in parentheses. Expression levels are color coded with red indicating upregulation by PPV infection; green indicating downregulation by PPV infection; and black indicating no change in expression. The intensity of color represents the degree of gene expression levels. List of genes induced or repressed at each time point in different cluster groups along with their expression levels and putative functions is provided in [Additional file ##SUPPL##7##8##].</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Confirmation of relative expression levels of the transcripts selected from the microarray analysis with semi quantitative RT-PCR (sqRT-PCR)</bold>. Expression changes of 8 selected genes were determined by sqRT-PCR and microarray. The signal intensity of each transcript was normalized using At3g18780 (Actin 2). The <italic>Arabidopsis </italic>Genome Initiative (AGI) locus identifier of each gene is provided. The y-axis indicates the normalized expression level of the transcript. The x-axis represents hours post transfection. Expression ratios are the average of three independent hybridizations ± standard deviation (SD). hpt, hours post transfection; dpi, days post inoculation. The error bars represent the standard deviation for the signals from each of the three independent hybridizations.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Confirmation of relative expression levels of the transcripts selected from the microarray analysis with semi quantitative RT-PCR (sqRT-PCR) or Northern blot hybridizations</bold>. Expression changes of 5 selected genes were determined by sqRT-PCR, microarray and Northern blot. The signal intensity of each transcript was normalized using At3g18780 (Actin 2). The <italic>Arabidopsis </italic>Genome Initiative (AGI) locus identifier of each gene is provided. The y-axis indicates the normalized expression level of the transcript. The x-axis represents hours post transfection or days post inoculation. Expression ratios are the average of three independent hybridizations ± standard deviation (SD). hpt, hours post transfection; dpi, days post inoculation. The error bars represent the standard deviation for the signals from each of the three independent hybridizations.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Assignment of upregulated and downregulated genes in each functional group identified in PPV-infected <italic>Arabidopsis </italic>leaves and protoplasts.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">S.No</td><td align="left">Functional category</td><td align="center">Number of genes on the array <sup>a</sup></td><td align="center">PPV-infected leaves Upregulated/Downregulated [Number (%)/Number (%)] <sup>b</sup></td><td align="center">pPPV-SK68-transfected protoplasts Upregulated/Downregulated [Number (%)/Number (%)] <sup>b</sup></td></tr></thead><tbody><tr><td align="center">A</td><td align="left">Defence and virulence</td><td align="center">1488</td><td align="center">204 (13.7)/165 (11.1)</td><td align="center">32 (2.15)/33 (2.22)</td></tr><tr><td align="center">B</td><td align="left">Cellular communication/Signal transduction mechanism/transmembrane signal transduction</td><td align="center">1657</td><td align="center">94 (5.67)/74 (4.47)</td><td align="center">16 (0.97)/18 (1.09)</td></tr><tr><td align="center">C</td><td align="left">Cell Cycle/DNA processing/chromatin regulation and cytoskeleton reorganization</td><td align="center">771</td><td align="center">87 (11.3)/30 (3.89)</td><td align="center">10 (1.30)/9 (1.17)</td></tr><tr><td align="center">D</td><td align="left">Development/storage proteins</td><td align="center">657</td><td align="center">138 (21)/100 (15.2)</td><td align="center">20 (3.04)/20 (3.04)</td></tr><tr><td align="center">E</td><td align="left">Metabolism</td><td align="center">4238</td><td align="center">548 (12.9)/226 (5.3)</td><td align="center">35 (0.83)/41 (0.97)</td></tr><tr><td align="center">F</td><td align="left">Proteins with binding function or cofactor requirement</td><td align="center">921</td><td align="center">210 (22.8)/139 (15.1)</td><td align="center">20 (2.17)/22 (2.39)</td></tr><tr><td align="center">G</td><td align="left">Protein synthesis and translation</td><td align="center">314</td><td align="center">37 (11.8)/39 (12.4)</td><td align="center">1 (0.32)/1 (0.32)</td></tr><tr><td align="center">H</td><td align="left">Protein fate</td><td align="center">647</td><td align="center">77 (11.9)/69 (10.6)</td><td align="center">9 (1.39)/12 (1.85)</td></tr><tr><td align="center">I</td><td align="left">Transcription/splicing/RNA processing/modification</td><td align="center">1013</td><td align="center">175 (17.3)/102 (10.1)</td><td align="center">32 (3.16)/29 (2.86)</td></tr><tr><td align="center">J</td><td align="left">Transporters</td><td align="center">563</td><td align="center">78 (13.9)/19 (3.4)</td><td align="center">7 (1.24)/7 (1.24)</td></tr><tr><td align="center">K</td><td align="left">Transposable elements, viral and plasmid proteins</td><td align="center">124</td><td align="center">22 (17.7)/44 (35.5)</td><td align="center">8 (6.45)/21 (16.94)</td></tr><tr><td align="center">L</td><td align="left">Unknown proteins</td><td align="center">10417</td><td align="center">343 (3.3)/450 (4.3)</td><td align="center">73 (0.70)/91 (0.87)</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Primers used in semi quantitative RT-PCR<sup>a </sup>and in Northern hybridzations<sup>b</sup></p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Probe set IDs </bold><sup>c</sup></td><td align="left"><bold>AGI</bold><sup>d</sup><bold>locus</bold></td><td align="left"><bold>Description</bold></td><td align="left"><bold>Forward Primer (5'---3')</bold></td><td align="left"><bold>Reverse Primer (5'---3')</bold></td></tr></thead><tbody><tr><td align="left" colspan="5"><bold>A. Genes induced only in PPV-infected <italic>Arabidopsis </italic>protoplasts</bold></td></tr><tr><td align="left">256487_at</td><td align="left">At1g31540</td><td align="left">Disease resistance protein, TIR-NBS-LRR class</td><td align="left">GCCAGTAGAGCTCGAGGTAAAATG</td><td align="left">AGGCTCTTCACTGTGTTCAATCTC</td></tr><tr><td align="left">246338_s_at</td><td align="left">At3g44930</td><td align="left">Cation/hydrogen exchanger CHX10, putative</td><td align="left">TAACCCGAGGATCCGTGTAAC</td><td align="left">ATCTCTAGCCGCCAGCAAATC</td></tr><tr><td align="left">265682_at</td><td align="left">At2g24390</td><td align="left">Avirulence induced gene (AIG) protein, related</td><td align="left">CAATGCACAGTGTCTTCGTCTA</td><td align="left">TTCGGTAAATTCAACTCTTCAGC</td></tr><tr><td align="left">254582_at</td><td align="left">At4g19470</td><td align="left">Disease resistance protein, related</td><td align="left">CAGAAACGAAGACGGTGATTAC</td><td align="left">TCTTCCTTTATTCTTTTCCTGTTC</td></tr><tr><td align="left">267056_at</td><td align="left">At2g32470</td><td align="left">F-box family protein, related</td><td align="left">TGGGGTTTCACAGTCAATAATAGT</td><td align="left">TTCGCATACCAATCCCACAAC</td></tr><tr><td align="left">248188_at</td><td align="left">At5g54070</td><td align="left">Heat shock transcription factor</td><td align="left">TGAGTCATTGAAGGAGGAACAGAG</td><td align="left">GAGCCACCGTCAACAAGTAGG</td></tr><tr><td align="left">253366_at</td><td align="left">At4g33180</td><td align="left">Hydrolase family protein</td><td align="left">TCGCAGACGCTTTCCATAGACT</td><td align="left">GAAGAAGACAAGATCGGGAGAATA</td></tr><tr><td align="left">260699_at</td><td align="left">At1g32320</td><td align="left">MAP kinase, putative</td><td align="left">GACATAAAGCCTTCAAATCTCC</td><td align="left">ATCTCCCAACAAAATCTCTAAACT</td></tr><tr><td align="left">249440_at</td><td align="left">At5g40030</td><td align="left">Protein kinase PK1- like</td><td align="left">CTTCAAGCCTCGGTTCCTCAA</td><td align="left">GATCGGTTTCGGTATTTCTGGTG</td></tr><tr><td align="left" colspan="5"><bold>B. Genes induced in PPV-infected <italic>Arabidopsis </italic>protoplasts and leaf tissues</bold></td></tr><tr><td align="left">260921_at</td><td align="left">At1g21540</td><td align="left">AMP-binding protein, putative</td><td align="left">TTAGGGACCCGAGAACC</td><td align="left">GCAGAGACTTAGCCATTTGT</td></tr><tr><td align="left">258962_at</td><td align="left">At3g10570</td><td align="left">Cytochrome P450, putative</td><td align="left">TCGTGAGGTTGAGGAAAAAGATG</td><td align="left">GATAAGCGCTCCACTCAAACTC</td></tr><tr><td align="left">253502_at</td><td align="left">At4g31940</td><td align="left">Cytochrome P450-like protein</td><td align="left">TGCACGGTCGCTGGTTACT</td><td align="left">AGGCGTGGCTTTAGGAATG</td></tr><tr><td align="left">248989_at</td><td align="left">At5g45200</td><td align="left">Disease resistance protein, TIR-NBS-LRR</td><td align="left">ATTCTTGGTGGTTGGACTG</td><td align="left">CACTTCTCTTTGGACTTTC</td></tr><tr><td align="left">257229_at</td><td align="left">At3g16490</td><td align="left">Calmodulin-binding family protein</td><td align="left">AATGAATGTGGCTGTATCC</td><td align="left">CGTTTGCGGTTGTGGTTGT</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Supplemental Table 1. Expression levels of genes induced in PPV-infected <italic>Arabidopsis </italic>leaf tissues 17 days post inoculation.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Supplemental Table 2. Expression levels of genes repressed in PPV-infected <italic>Arabidopsis </italic>leaf tissues 17 days post inoculation.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Supplemental Table 3. Expression levels of differentially regulated genes in PPV-infected <italic>Arabidopsis </italic>protoplasts.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Supplemental Table 4. Expression levels of common genes induced in PPV-infected <italic>Arabidopsis </italic>protoplasts and in infected leaf tissues.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>Supplemental Table 5. Expression levels of common genes repressed in PPV-infected <italic>Arabidopsis </italic>protoplasts and in infected leaf tissues.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional file 6</title><p>Supplemental Figure 1. Hierarchical clustering of common genes that were significantly differentially regulated by PPV in infected leaf tissues and in transfected protoplast cells. Expression levels are color-coded with red indicating upregulation by pPPV-SK68; green indicating downregulation by pPPV-SK68 and black indicating no change in expression. The intensity of color represents the degree of gene expression levels. The putative function of each gene is shown on the right side of the cluster.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional file 7</title><p>Supplemental Figure 2. K-means profiling of gene expression using 411 <italic>Arabidopsis </italic>genes differentially regulated by PPV in transfected protoplasts at three different time points. The expression profiles were grouped into twelve distinct cluster groups. The AGI locus identifier of each gene differentially regulated by PPV in the transfected protoplasts at different time points in each cluster group is shown on the right side of the cluster. Values on the y-axis indicate the relative expression level of the gene, while the x-axis represents hours post transfection. Number of genes belonging to each cluster is shown in the cluster inset.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional file 8</title><p>Supplemental Table 6. Gene expression profiles of <italic>Arabidopsis </italic>genes differentially regulated by PPV infectious clone in transfected protoplasts belonging to twelve distinct cluster groups.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional file 9</title><p>Supplemental Table 7. Identification of <italic>Prunus persica </italic>orthologs to <italic>Arabidopsis </italic>genes induced by PPV infection in the leaf tissues at 17 days post inoculation.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional file 10</title><p>Supplemental Table 8. Cross comparison of genes significantly differentially regulated by PPV (≥ 2.5- or ≤ -2.5-fold) in this study with the <italic>Arabidopsis </italic>genes regulated by other positive sense RNA viruses.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S11"><caption><title>Additional file 11</title><p>Supplemental Figure 3. Confirmation of microarray data by sqRT-PCR and Northern hybridizations. Panel A shows confirmation of microarray data using sqRT-PCR for genes induced in PPV-infected <italic>Arabidopsis </italic>protoplasts. Panel B shows confirmation of microarray data using sqRT-PCR (left panel) and Northern hybridizations (right panel) for <italic>Arabidopsis </italic>genes differentially regulated in PPV-infected protoplasts and in PPV-infected leaves. Probes for sqRT-PCR and Northern hybridizations were generated by PCR amplification of <italic>Arabidopsis </italic>cDNA using gene specific primers shown in Table 3. sqRT-PCR of the constitutively expressed Actin 2 gene (At3g18780) was used as a loading control. pPPV-SK68, a PPV infectious cDNA clone used to transfect protoplasts; pPPV-SK68Δ, a mutant non-infectious clone of pPPV-SK68 was used as a control; hpt, hours post transfection.</p></caption></supplementary-material>
<table-wrap-foot><p><sup>a </sup>Number of genes in a given category represented on the chip.</p><p><sup>b </sup>Percentage was calculated based on the number of differentially regulated genes divided by the total number of genes in the functional category on the chip.</p></table-wrap-foot> <table-wrap-foot><p><sup>a</sup>Genes shown in panels "A" and "B" were validated using semi quantitative RT-PCR.</p><p><sup>b</sup>Genes shown in panel "B" were validated using Northern hybridizations.</p><p><sup>c</sup>Probe set ID represents Affymeterix probe set number.</p><p><sup>d</sup>AGI represents <italic>Arabidopsis </italic>Genome Initiative (AGI) locus identifier that corresponds to each gene represented on the array.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-325-1"/> <graphic xlink:href="1471-2164-9-325-2"/> <graphic xlink:href="1471-2164-9-325-3"/> <graphic xlink:href="1471-2164-9-325-4"/> <graphic xlink:href="1471-2164-9-325-5"/> <graphic xlink:href="1471-2164-9-325-6"/> <graphic xlink:href="1471-2164-9-325-7"/>
<media xlink:href="1471-2164-9-325-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S5.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S6.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S7.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S8.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S9.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S10.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-325-S11.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Stange"], "given-names": ["C"], "article-title": ["Plant-virus interactions during the infective process"], "source": ["Cien Inv Agr"], "year": ["2006"], "volume": ["33"], "fpage": ["1"], "lpage": ["18"]}, {"surname": ["Wang", "Sanfa\u00e7on", "Stobbs", "James", "Thompson", "Svircev", "Brown"], "given-names": ["A", "H", "LW", "D", "D", "AM", "DCW"], "italic": ["Plum pox virus "], "source": ["Can J Plant Pathol"], "year": ["2006"], "volume": ["28"], "fpage": ["182"], "lpage": ["196"]}, {"surname": ["Pompe-Novak", "Gruden", "Baebler", "Kre\u010di\u010dStres", "Kova\u010d", "Jongsma", "Ravnikar"], "given-names": ["M", "K", "\u0160", "H", "M", "M", "M"], "italic": ["Potato virus Y ", "Solanum tuberosum "], "source": ["Physiol Mol Plant Pathol"], "year": ["2006"], "volume": ["67"], "fpage": ["237"], "lpage": ["247"]}, {"surname": ["Maule", "Escaler", "Aranda"], "given-names": ["AJ", "M", "MA"], "article-title": ["Programmed responses to virus replication in plants"], "source": ["Mol Plant Pathol"], "year": ["2000"], "volume": ["1"], "fpage": ["9"], "lpage": ["15"]}, {"surname": ["Carr", "Whitham"], "given-names": ["T", "SA"], "article-title": ["An emerging model system: "], "italic": ["Arabidopsis "], "source": ["Plant Cell Monogr"], "year": ["2007"], "volume": ["2"], "fpage": ["159"], "lpage": ["183"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate-a practical and powerful approach to multiple testing"], "source": ["JR Stat Soc Ser"], "year": ["1995"], "volume": ["B57"], "fpage": ["289"], "lpage": ["300"]}, {"surname": ["Balachandran", "Hull", "Vaadia", "Wolf", "Lucas"], "given-names": ["S", "RJ", "Y", "S", "WJ"], "article-title": ["Alteration of carbon partitioning induced by the movement protein of "], "italic": ["Tobacco mosaic virus "], "source": ["Plant Cell Environ"], "year": ["1995"], "volume": ["18"], "fpage": ["1301"], "lpage": ["1310"]}, {"surname": ["T\u00e9csi", "Maule", "Smith", "Leegood"], "given-names": ["LI", "AJ", "AM", "RC"], "article-title": ["Complex, localized changes in CO2 assimilation and starch content associated with the susceptible interaction between "], "italic": ["Cucumber mosaic virus "], "source": ["Plant J"], "year": ["1994"], "volume": ["6"], "fpage": ["837"], "lpage": ["847"]}, {"surname": ["Handford", "Carr", "Loebenstein G, Carr JP"], "given-names": ["MG", "JP"], "article-title": ["Plant metabolism associated with resistance and susceptibility"], "source": ["Natural Resistance Mechanisms of Plants to Viruses"], "year": ["2006"], "publisher-name": ["Dordrecht: Springer"], "fpage": ["315"], "lpage": ["340"]}, {"surname": ["Loebenstein", "Akad", "Loebenstein G, Carr JP"], "given-names": ["G", "F"], "article-title": ["The local lesion response"], "source": ["Natural Resistance Mechanisms of Plants to Viruses"], "year": ["2006"], "publisher-name": ["Dordrecht: Springer"], "fpage": ["99"], "lpage": ["124"]}, {"surname": ["Herbers", "Takahata", "Melzer", "Mock", "Hajirezaei", "Sonnewald"], "given-names": ["K", "Y", "M", "HP", "M", "U"], "article-title": ["Regulation of carbohydrate partitioning during the interaction of potato virus Y with tobacco"], "source": ["Mol Plant Pathol"], "year": ["2000"], "volume": ["1"], "fpage": ["51"], "lpage": ["59"]}, {"surname": ["Sanfa\u00e7on"], "given-names": ["H"], "article-title": ["Replication of positive-strand RNA viruses in plants: contact points between plant and virus components"], "source": ["Can J Bot"], "year": ["2005"], "volume": ["83"], "fpage": ["1529"], "lpage": ["1549"]}, {"surname": ["Smith", "Rodriguez-Buey", "Karlsson", "Campbell"], "given-names": ["CM", "M", "J", "MM"], "article-title": ["The response of the poplar transcriptome to wounding and subsequent infection by a viral pathogen"], "source": ["New Phytologist"], "year": ["2004"], "volume": ["164"], "fpage": ["123"], "lpage": ["136"]}, {"surname": ["Takebe"], "given-names": ["I"], "article-title": ["The use of protoplasts in plant virology"], "source": ["Ann Rev Phytopathol"], "year": ["1975"], "volume": ["13"], "fpage": ["105"], "lpage": ["125"]}, {"surname": ["Ramsay", "Jones", "Wang"], "given-names": ["K", "MGK", "Z"], "article-title": ["Laser capture microdissection: a novel approach to microanalysis of plant-microbe interactions"], "source": ["Mol Plant Pathol"], "year": ["2006"], "volume": ["7"], "fpage": ["429"], "lpage": ["435"]}, {"surname": ["Wang", "Chapman", "Chen", "Stobbs", "Brown", "Brandle"], "given-names": ["A", "P", "L", "LW", "DCW", "JE"], "article-title": ["A comparative survey, by expressed sequence tag analysis, of genes expressed in peach leaves infected with "], "italic": ["Plum pox virus "], "source": ["Can J Plant Pathol"], "year": ["2005"], "volume": ["27"], "fpage": ["410"], "lpage": ["419"]}, {"surname": ["Ll\u00e1cer"], "given-names": ["G"], "article-title": ["Hosts and symptoms of "], "italic": ["plum pox virus"], "source": ["OEPP/EPPO Bulletin"], "year": ["2006"], "volume": ["36"], "fpage": ["227"], "lpage": ["228"]}, {"surname": ["Maule", "Caranta", "Boulton"], "given-names": ["AJ", "C", "MI"], "article-title": ["Sources of natural resistance to plant viruses: status and prospects"], "source": ["Mol Plant Pathol"], "year": ["2007"], "volume": ["8"], "fpage": ["223"], "lpage": ["231"]}, {"surname": ["Querci", "Baulcombe", "Goldbach", "Salazar"], "given-names": ["M", "DC", "RW", "LF"], "article-title": ["Analysis of resistance breaking determinants of potato virus \u00d7 (PVX) strain HB on different potato genotypes expressing extreme resistance to PVX"], "source": ["Phytopathology"], "year": ["1995"], "volume": ["85"], "fpage": ["1003"], "lpage": ["1011"]}, {"surname": ["Dunigan", "Golemboski", "Zaitlin", "Evered D, Harnett S"], "given-names": ["DD", "DB", "M"], "article-title": ["Analysis of the N gene of Nicotiana"], "source": ["Plant Resistance to Viruses"], "year": ["1987"], "publisher-name": ["Chichester, UK: Wiley"], "fpage": ["120"], "lpage": ["135"]}, {"surname": ["Rottier"], "given-names": ["P"], "article-title": ["Viral protein synthesis in cowpea mosaic virus infected protoplasts"], "source": ["PhD thesis"], "year": ["1980"], "publisher-name": ["University of Wageningen, The Netherlands"]}, {"surname": ["Hammond-Kosack", "Jones", "Buchanan BB, Gruissem W, Jones RL"], "given-names": ["K", "JDG"], "article-title": ["Responses to plant pathogens"], "source": ["Biochemistry and Molecular Biology of plants"], "year": ["2000"], "publisher-name": ["Maryland: American Society of Plant Physiology"], "fpage": ["1102"], "lpage": ["1156"]}]
BMC Genomics. 2008 Jul 9; 9:325