- Open Access
A remarkable synergistic effect at the transcriptomic level in peach fruits doubly infected by prunus necrotic ringspot virus and peach latent mosaic viroid
Virology Journalvolume 10, Article number: 164 (2013)
Microarray profiling is a powerful technique to investigate expression changes of large amounts of genes in response to specific environmental conditions. The majority of the studies investigating gene expression changes in virus-infected plants are limited to interactions between a virus and a model host plant, which usually is Arabidopsis thaliana or Nicotiana benthamiana. In the present work, we performed microarray profiling to explore changes in the expression profile of field-grown Prunus persica (peach) originating from Chile upon single and double infection with Prunus necrotic ringspot virus (PNRSV) and Peach latent mosaic viroid (PLMVd), worldwide natural pathogens of peach trees.
Upon single PLMVd or PNRSV infection, the number of statistically significant gene expression changes was relatively low. By contrast, doubly-infected fruits presented a high number of differentially regulated genes. Among these, down-regulated genes were prevalent. Functional categorization of the gene expression changes upon double PLMVd and PNRSV infection revealed protein modification and degradation as the functional category with the highest percentage of repressed genes whereas induced genes encoded mainly proteins related to phosphate, C-compound and carbohydrate metabolism and also protein modification. Overrepresentation analysis upon double infection with PLMVd and PNRSV revealed specific functional categories over- and underrepresented among the repressed genes indicating active counter-defense mechanisms of the pathogens during infection.
Our results identify a novel synergistic effect of PLMVd and PNRSV on the transcriptome of peach fruits. We demonstrate that mixed infections, which occur frequently in field conditions, result in a more complex transcriptional response than that observed in single infections. Thus, our data demonstrate for the first time that the simultaneous infection of a viroid and a plant virus synergistically affect the host transcriptome in infected peach fruits. These field studies can help to fully understand plant-pathogen interactions and to develop appropriate crop protection strategies.
Plant viruses and viroids are obligate parasites that, consequently, depend on host factors to complete their life cycle. Subversion of host factors for viral replication and spread often disrupts cellular homeostasis and leads to pathogenesis and disease symptoms . Among the most important consequences of viral pathogenesis are changes in the expression of host genes [2–4]. However, the functional significance of transcript changes in specific plant-pathogen interactions is still not well understood [2, 5, 6]. Thus, an improved understanding of these interactions should lead to new and creative methods to control plant virus/viroid diseases.
Stone fruit viruses and viroids can cause significant crop damage and crop losses . Viroids are pathogens of food, industrial and ornamental plants. PLMVd, a member of the family Avsunviroidae, consists of a circular monocatenary RNA of 335-351 nucleotides and is distributed worldwide on peach . Although PLMVd was initially considered as a latent pathogen it has been described that several isolates of PLMVd can cause deformation, discolored spots and cracked sutures in infected peach fruit . PNRSV is a member of the genus Ilarvirus in the family Bromoviridae. PNRSV has been reported to cause a reduction of 12% to 70% in tree growth and 5% to 70% yield losses depending on the cultivar [7, 9].
Microarray technology is a powerful method to simultaneously evaluate changes in the expression of large amounts of genes in response to specific conditions. Currently, an increasing number of studies employ microarray profiling and de novo RNA sequencing to describe gene expression changes in virus/viroid-infected plants . These studies reveal a significant impact of the viral/viroidal infection on host gene expression. Among the cellular processes affected upon infection are defense/stress responses, cell wall structure, chloroplast function, hormone signaling, protein metabolism, silencing and diverse other functions [2, 5, 10–17]. Moreover, infection leads to the modulated expression of various key regulators such as numerous transcription factors, antioxidants, metabolic enzymes and transporters . However, the majority of studies analyzing gene expression changes upon infection are limited to specific host-pathogen interactions involving a given virus and a model host plant, usually Arabidopsis thaliana or Nicotiana benthamiana. Transcriptome analysis upon PNRSV infection, for example has only been performed in N. benthamiana plants . PLMVd infection has not been investigated at the transcript level. Here, we explore and compare the gene expression profiles of peach, a natural host of PNRSV and PLMVd, upon single and double infection with the virus and the viroid in field conditions. Mixed infections of two pathogens occur in natural situations and have biological and epidemiological implications: frequently, simultaneous infections cause symptoms with a higher severity than those induced in single infections by either of the two pathogens. This phenomenon is known as synergism in pathology  and has been mainly studied between two different viruses as, for example, mixed infection of Potato virus X (PVX) and Potato virus Y (PVY) [21–24].
Double infection of peach with PNRSV and PLMVd frequently occurs in the field. Although the absence of symptoms was a common characteristic for all the samples analyzed in this study, we wondered whether the double infection with PNRSV and PLMVd could affect host gene expression differently than single infection. To approach this question, we conducted microarray profiling of peach fruits originating from Chile, singly and doubly infected with PNRSV and PLMVd using microarray slides containing 4,261 unigenes obtained from mesocarp and exocarp tissues of two full-sib progeny contrasting for chilling injury (CI) . Our results clearly demonstrate, for the first time, a synergistic effect at the transcriptomic level between a viroid and a plant virus. The synergistic effect is represented by considerably increased amplitude of differentially expressed genes. The functional identity of those genes with significant alterations upon viral and/or viroidal infection and their putative role in the disease process are discussed.
Sampling and experimental design
Between 2006 and 2008 a number of orchards of different size were visited to assess their sanitary status (Fiore et al., unpublished data) in austral hemisphere spring (September to November) in Chile. Fruits of a total of 2,456 peach trees were sampled during the whole survey period and tested for the presence of viruses and viroids known to affect stone fruit plants in Chile. Among others, samples were tested positive for PNRSV and Prune dwarf virus (PDV), which are the viruses most frequently encountered in southeastern peach, Tomato ringspot virus (TomRSV), Plum pox virus (PPV) and PLMVd. Interestingly, the analysis also revealed that the most recurrent mixed infection was PLMVd-PNRSV (Fiore et al., unpublished data).
For microarray analysis, peach fruits infected either with PNRSV, PLMVd or simultaneously with both pathogens were collected at developmental stage (S4), at approximately 100 days after boom, period where the fruit reaches the final full size and enters the fruit ripening or climacteric stage according with the description provided by Zanchin et al. . The tissue used for RNA extraction was exocarp and mesocarp. The analysis of the PLMVd sequences recovered from the infected trees revealed a mix of quasispecies (e.g. PLMVd-4tun, ZZ46, lib P7, PC-C29 etc.) none of them associated to symptomatic strains. For the PNRSV all the sequence variants recovered belong to the PV32 group for which no specific symptomatology has been observed . To obtain robust and statistically accurate data, four biological replicates and two technical, dye swapped replicates of healthy, virus infected, viroid infected, and doubly infected samples were analyzed. Each biological replicate consisted of a pool of four peach fruits from four different trees previously tested positive (for the infected samples) or negative (for the healthy controls) (Figure 1). RNA of the different sets of samples was hybridized to cDNA microarrays representing all the unigenes in the ChillPeach (http://bioinfo.ibmcp.upv.es/genomics/ChillPeachDB) database.
Comparison of the PNRSV and PLMVd concentrations in infected peach fruits
To determine the concentration of PNRSV and/or PLMVd in peach fruits, serial dilutions of each total RNA preparation were blotted onto a nylon membrane. The RNA samples were analyzed using specific digoxigenin-labelled riboprobes [27, 28]. The concentration of either PLMVd or PNRSV (Figure 2: samples A, B, C, D and I, J, K and L respectively) was very similar among the singly infected samples reaching in both cases a detection limit of 0.3 ng/μL. In the mixed infections (samples E, F, G and H) a slight increase in the accumulation of both pathogens was observed (Figure 2 and Additional file 1: Figure S1).
Cumulative changes in gene expression
Microarray analysis of the samples was performed using the ChillPeach array containing 4,261 peach unigenes. The number of genes with statistically significant expression changes of at least 1.5 fold relative to healthy plants was 16 and 82 for fruits infected with PLMVd and PNRSV, respectively (Figure 3a). Out of these genes, 13 and 66 have orthologs in A. thaliana. Interestingly, fruits simultaneously infected with both pathogens presented a significantly higher number of differentially regulated genes (685, out of which 602 have orthologs in A. thaliana). The majority of differentially regulated genes exhibited fold changes between 1.5 and 2.0 (Figure 3b) upon infection, however, in contrast to singly, doubly infected fruits displayed a considerable number of genes with fold changes above 2. Comparison of induced and repressed genes revealed similar results as the analysis of the total number of genes (Figure 3c and d). Moreover, this analysis revealed that virus, viroid, and double infection with both pathogens mainly lead to the repression of host gene expression. In detail, approximately 70% of genes differentially regulated upon PNRSV, PLMVd, or double infection, respectively were downregulated, while the remaining 30% were induced (Figure 3c and d). Out of the genes with orthologs in Arabidopsis only a few were induced in PNRSV (4 genes) and PLMVd (9 genes) infected fruits whereas 202 genes exhibited increased expression levels compared to the healthy control in peaches infected with both pathogens (Figure 4a and Additional file 2: Table S1). 400, 62 and 4 genes with orthologs in Arabidopis displayed significantly reduced expression upon infection, in doubly, PNRSV and PLMVd infected samples, respectively (Figure 4b and Additional file 2: Table S1). Interestingly 70% (44 genes) of the PNRSV infected genes were found in PNRSV-PLMVd double infection (meaning 10% for double infected peaches), thus indicating that PNRSV and PNRSV-PLMVd infection deregulate a largely similar set of genes.
Taken together, transcript analysis revealed a higher number of repressed than induced genes upon infection. In addition, double infection with PNRSV and PLMVd resulted in many more expression changes compared to single infection with either pathogen. Thus, PNRSV and PLMVd appear to exhibit a synergistic effect with respect to the regulation of host gene expression.
Validation of expression changes by quantitative RT-PCR
As reference genes for quantitative RT-PCR we chose peach orthologs of Arabidopsis F-box family protein (F-BOX) and elongation factor 1-α (EF1-α), which were previously reported to exhibit stable expression upon virus infection in Arabidopsis  and the clathrin adaptor complex (CAC) as internal control of the microarray. All reference genes exhibited similar expression levels in virus, viroid infected or healthy peach fruits, demonstrating their suitability as reference genes in peach. To validate the results obtained by microarray analysis we performed quantitative RT-PCR with nine randomly selected genes displaying induced or reduced expression upon infection plus the CAC. The nine selected genes were Phytoene synthase, Auxin response protein (IAA9), Expansin (EXP8), Glutamate descarboxilase, Glutamate dehydrogenase 2 (GDH2), Cysteine proteinase (RD21A), Invertase/pectin methylesterase inhibitor family protein, Universal stress protein (USP) family protein, CBL-interacting protein kinase 6 (CIPK6) and Phosphate-responsive protein, putative (E). Out of these nine genes analyzed, seven displayed consistent gene expression changes with both methods (Figure 5). Phytoene synthase, Auxin response protein, Glutamate descarboxilase, Glutamate dehydrogenase 2 (GDH2), Cysteine proteinase (RD21A), Invertase/pectin methylesterase inhibitor family protein, CBL-interacting protein kinase 6 (CIPK6) exhibited reduced expression in PNRSV, PLMVd and PNRSV-PLMVd infected samples compared to the healthy control plants, and Glutamate decarboxylase was induced upon infection. As already seen by microarray analysis, all expression changes were more pronounced upon double infection with PNRSV and PLMVd. Expression levels for Expansin, however, differed in Microarray and quantitative RT-PCR experiments. While the gene was upregulated in microarray analyses, it was downregulated by qRT-PCR in PNRSV infected samples. Expansins exist in a multigene family in Arabidopsis . It is likely, that also peach has several expansin isoforms. Thus, the differences in expression of Expansin upon infection may result from several isoforms recognized by the cDNA array and the specific detection of only one isoform in our quantitative RT-PCR experiments. In a similar way, expression levels for Universal stress protein (USP) family protein are distinct in the two different analysis. A similar situation to that observed for Expansins occurs with this gene and thus, the same explanation could be valid for this protein family. Taken together, the results of qRT-PCR analyses confirmed the gene expression changes seen by microarrays analysis.
Functional categorization of genes with induced and repressed expression upon double infection with PNRSV and PLMVd
As stated above, a clear synergistic effect on the peach transcriptome was observed in samples simultaneously infected with both pathogens. To gain insight into functions of genes exhibiting modulated expression upon infection we grouped the genes with significantly modulated expression upon double infection with PNRSV and PLMVd into 23 functional categories using the Munich Information Center for Protein Sequences Functional Catalogue Database (MIPS FunCatDB, http://mips.helmholtz-muenchen.de/proj/funcatDB/) (Figure 6). Most repressed or induced genes can be assigned to three functional categories: “metabolism”, “protein with binding function or cofactor requirement” and “subcellular localization”. “Protein modification and degradation” was the functional category to which most repressed genes could be assigned, whereas induced genes were mainly predicted to encode proteins in the category “related to phosphate, C-compound and carbohydrate metabolism” and also “protein modification” (Table 1).
Statistical analysis of functional categories
To identify overrepresented functional categories among induced or repressed genes, we subjected genes to analysis by Babelomics 4.2 FatiGO and single enrichment analysis (SEA). FatiGO uses Fisher’s exact test for 2×2 contingency tables to check for significant over-representation of GO terms. Singular enrichment analysis identifies enriched GO terms in a list of microarray probe sets or gene identifiers. We used this test to identify over- or underrepresented GO terms in the set of genes with significant changes upon infection with both pathogens compared to the whole set of genes represented on the microarray. This analysis should enable us to predict a role of a certain biological processes in response to infection rather than merely calculate a number of up-regulated and down-regulated genes.
For over-representation analysis, each gene was assigned to one of the three main gene ontologies (GO): Biological Process, Cellular Component and Molecular Function. For the 198 up-regulated genes unique to double infection, no functional category was over- or under- represented with respect to the functional category distribution in the microarray. Table 2 shows the non-redundant functional categories significantly over- and under-represented among those genes that were significantly down-regulated upon double infection. The enrichment analysis resulted in a total of 6, 7 and 5 functional categories in GO Biological Process, Molecular Function and Cellular Component respectively. To estimate the number of times that a certain functional category was over- or under- represented, we calculated the ratio between the percentage of genes in each functional category and the percentage of genes corresponding to the same functional category in the whole array (Table 2). We found higher ratios for “alcohol metabolic process” and “response to cadmium ion” in GO Biological Process, “oxidoreductase activity, acting on NADH or NADPH” in GO Molecular Function and “peroxisome” in GO Cellular Component.
To gain a more precise idea about functional classes of proteins with altered expression upon double infection with PNRSV and PLMVd, we increased the threshold of genes subjected to over-representation analysis to at least 2-fold changes in expression level. With these settings, the majority of the functional groups over- represented within GO Biological Process was “related with cellular amino acid metabolic process” and “response to abiotic stresses”. Interestingly, “response to water deprivation” was overrepresented five times compared to the normal distribution of this functional class on the array (Table 2). No functional category was over- or under- represented in GO Molecular Function and only one category with 19 genes, “vacuole”, was overrepresented in GO Cellular Component.
Viroids and viruses depend on host components to complete their life cycle and thus interfere with host processes, such as host gene expression during infection. The identification of genes with altered expression in response to a certain virus/viroid infection can provide clues about the viral life cycle and requirements which need to be met by the host cell. Moreover, a better understanding of host responses to infection facilitates the development of methods to control pathogen infection. To date, the effects of viroid infection on the host transcriptome have been reported for four viroids. These studies used PCR-based cDNA library subtraction ; differential display RT-PCR and quantitative real-time RT-PCR [31, 32] or cDNA microarray [17, 33] to detect host gene expression changes. Currently, a number of studies explore gene expression changes in virus-infected plants using DNA microarrays. All these studies involve a given virus and a model host plant, usually A. thaliana or N. benthamiana. In the present work, we collected peach fruits, a natural host of both PNRSV and PLMVd, at development stage (S4) and studied changes in their transcriptome upon infection with the two pathogens. Today, no array including the whole peach fruit genome is available. Thus, we used a cDNA microarray containing probes for 4,261 genes expressed during chilling injury (CI) development .
During infection, the concentration of PNRSV and PLMVd in fruits is significantly higher (100:1) than in leaves [35, 36]. In addition, the occurrence of peach trees doubly infected with both, PNRSV and PLMVd in the field is considerable (e.g. Fiore et al., unpublished data); . Here, we compared transcript changes in asymptomatic singly and doubly PNRSV or PLMVd infected peach fruits. Overall, our gene expression analysis revealed a relatively low number of differentially expressed genes in single PNRSV infections compared to the healthy control plants (82 genes). This is consistent with the results obtained by Dardick  in which only 89 genes were significantly differentially expressed in PNRSV-infected N. benthamiana plants compared to mock controls. Peach fruits infected with PLMVd also mounted only a moderate response at transcript level (16 genes). The PLMVd variants detected in our study are all considered as latent. In future studies it would be interesting to compare the host transcripts changes caused by variants causing the “calico”  syndrome with those caused by milder variants. PNRSV infected samples were also asymptomatic and all the sequence variants recovered belonged to the PV32 group . Remarkably, in the doubly infected peach fruits a significant synergistic effect on the host transcriptome was observed. Considering the three different infections, the total number of genes with significantly altered expression (at least 1.5 fold change in expression level) was 783, which represents 18% of the whole array. From these 783 genes and eliminating those with significant expression changes occurring in common upon the three scenarios, 627 have orthologs in Arabidopsis (Additional file 3: Table S2). Among those genes 211 exhibited induced and 416 reduced expression. Thus, the total number of genes with reduced expression exceeded that of genes with induced expression by two-fold. This is surprising as comparative analysis of the Arabidopsis transcriptome during compatible interaction with plant viruses  revealed that there was a greater number of up-regulated than repressed genes in the course of viral pathogenesis. However, the authors demonstrated that each virus-host interaction is unique in terms of the genes with altered expression levels and to find a common pattern among different viruses is difficult . In addition, in one of the few studies in which a temporal analysis after virus/viroid infection was carried out, Rizza et al.,  showed that the pattern of up-regulated vs down regulated genes can change in pre-symptomatic when compared to post-symptomatic stages of Etrog citron infected with Citrus exocortis viroid (CEVd).
The overlap in significantly altered gene expression among the two single and the double infections was low. Common genes that were differentially expressed during each infection scenario are shown in Additional file 4: Table S3. This is not surprising considering the small number of genes with significant changes in both single infections. Consistently, Postnikova and Nemchinov  have recently shown that the number of host genes commonly affected during infection of Arabidopsis with either of twelve different viruses studied is very limited. Due to the nature of this specialized database it is difficult to establish a relationship between genes involved in peach ripening and those affected by the pathogens. Nevertheless in our study we observed some genes previously described as related with ripening with altered expression. Most of them did it in the, a priori, unexpected wayside slowing down the ripening (see Additional file 5: Table S4 and references [38–41]). Although we did not follow with great detail whether or not mixed infections had modified the ripening it is noteworthy to emphasize that alterations in ripening date have been described in several virus-host interactions regarding fruit trees (see  for review).
The higher number of genes significantly altered exclusively upon double infection with PLMVd and PNRSV led us to investigate whether some specific functional categories were over- or under- represented among those genes. Interestingly, among the 198 significantly induced genes unique to double infection (Figure 4) none of the functional categories appeared to be enriched in comparison with the normal distribution of genes in the whole microarray. As the array used in this study is enriched with sequences of genes that are implicated in chilling injury development, the lack of overrepresentation of functional categories among the induced genes upon infection compared to the genes represented on the array may indicate a similarity between the functional classes induced upon the two types of stresses, chilling injury and virus/viroid infection. By contrast, we identified over-represented functional categories and subcategories among the set of repressed genes belonging to each of the main GO domains (Table 2, threshold > 1.5 fold). We found genes involved in “response to external stimulus”, “defense response”, “catabolic processes” and “binding function or cofactor requirement”. Downregulation of these functions is presumably due to the virus counterattack against host defense-related pathways. Moreover, these host functions are frequently found in gene expression studies in response to virus infection .
By increasing the threshold to at least 2 fold changes we identified enriched functional categories only belonging to two of the main GO domains: Biological Process and Cellular Component (Table 2). Using this more stringent analysis we found groups of genes related with “response to stress”, to “external stimulus”, “amino acid metabolism” and “vacuole”. Importantly, the most overrepresented functional category was “response to water deprivation”. Among the genes assigned to this functional group genes related with “plant hormones” were found. Plant hormones play important roles in regulating developmental processes and signalling networks involved in plant responses to a wide range of biotic and abiotic stresses. For instance we identified a tonoplast resident H+-pump (At1g15690 H+-PPase AVP1) that contributes to vacuolar acidification, regulation of apoplastic pH and auxin transport . Auxin acts as an important component of hormone signaling network involved in the regulation of defense responses against various biotrophic and necrotrophic pathogens. Additionally, this hormone regulates the expression of genes associated with the biosynthesis, catabolism and signaling pathways of other hormones . Viral pathogens manipulate auxin signaling components to promote virulence and cause disease. One example is the interaction of Tobacco mosaic virus (TMV) replicase with Aux/IAA proteins. Development of symptoms promoted by this interaction has been described in Arabidopsis and tomato [45–47] and in addition disrupts the nuclear localization of several Arabidopsis Aux/IAA proteins. A second example is the protein phosphatase AtPP2CA (At3g11410) which acts as strong negative regulator of ABA signal transduction during seed germination and the regulation of stomatal closure . Tobacco plants infected with TMV showed increased ABA levels and treatment with ABA enhanced TMV resistance in tobacco . Some pathogens might have evolved the ability to produce ABA or ABA-mimicking substances to interfere with host defence. In any case, the role of ABA during plant-pathogen interactions depends on the individual plant-pathogen combination .
As members of the functional group “response to water deprivation” we also found proteins directly related with protection against oxidative stress like the Aldehyde dehydrogenase (At1g54100) and the Late embryogenesis abundant (LEA)-like protein (AtLEA5) (At4g02380). Aldehyde dehydrogenases (ALDHs) regulate the level of aldehydes by metabolizing excessive amounts of aldehyde molecules, which accumulate as a result of perturbed environmental conditions. Aldehydes are involved in different cellular metabolic processes but in excess they can have toxic effects on the cells. It is important then to maintain an appropriate balance of these molecules to avoid cellular damages. It has been described that over-expression of different ALDHs in A. thaliana confers tolerance to abiotic stress and protects plants against lipid peroxidation and oxidative stress . (LEA)-like proteins protect other proteins from aggregation from desiccation. AtLEA5 is the unique which is induced specifically by reactive oxygen species (ROS) as well as by ABA but it is unlikely to act as an antioxidant enzyme. It has been suggested that LEA5 may cooperate with other factors to protect cellular components against ROS-induced damage or indeed to enhance the turnover of specific proteins during stress to enable rapid acclimation to the prevailing conditions .
In the few cases in which changes in the host transcriptome have been investigated upon double virus infections (e.g. the PVX-PVY interaction in N. benthamiana) a severe oxidative stress was inferred to be induced in infected plant leaves, as increased transcript levels of genes encoding proteins important for lipid peroxidation for the generation of ROS were observed. We did not observe any significant over- or underrepresentation of genes involved in the antioxidative system in the doubly infected peach fruits. This difference between PLMVd-PNRSV-infected peach fruits and PVX-PVY infected N. benthamiana plants may be explained by the nature of the interaction between the two pathogens and the pathogens and the host. While co-infection of N. benthamiana with PVX-PVY can be considered as a true biological synergistic interaction leading to enhanced disease phenotypes compared to single infection with either pathogen, the doubly-infected PNRSV-PLMVd peach fruits did not exhibit an enhanced disease phenotype compared to single infection with either PNRSV or PLMVd. Interestingly, PNRSV was able to induce significant oxidative stress and an imbalance in the antioxidant systems in infected apricot-seeds  resulting in a decrease in seed germination. This observation indicates that host responses are not only specific for the host plant and the pathogen, but also tissue and/or organ-specific.
In the functional category “response to dehydration” which was strongly overrepresented upon double infection with PNRSV and PLMVd, we also found other important players with known functions in plant immunity. Among those was Cysteine protease RD21 (At1g47128), a Papain-like cysteine protease (PLCP). These proteins participate in immune responses and are targeted by pathogen- derived inhibitors. PLCPs are also required to trigger the hypersensitive response (HR) [54, 55]. Consistently, PLCP RD19 is required in Arabidopsis for RRS1-R- mediated resistance against the bacterial pathogen Ralstonia solanacearum producing effector PopP2 . A crucial role of PLCPs in disease immunity is also indicated by the observation that many pathogens produce effectors that manipulate these proteases [56–61].
Our data reveal for the first time a clear synergistic effect between a viroid and plant virus at the transcriptome level. We demonstrate that mixed infections, which occur frequently in the field, can result in a more complex transcriptional response than that observed in single infections under the same experimental conditions. Thus, our analysis takes first steps to illuminate the mechanistic basis of synergistic mixed infections in peach trees and reveals candidate genes which can be tried as targets for crop protection. Future additional studies will further elucidate facilitative or antagonistic interactions between plant viruses/viroids in mixed infections.
A total of 2,456 stone fruit trees were sampled during the whole survey period. Three trees were randomly sampled within each orchard. Twenty leaves per plant were collected from several points of the canopy and pooled leaves of each plant were tested for the presence of viruses and viroids previously reported from Chile, i.e. (PPV), (PDV), (PNRSV), (ACLSV), (ToRSV), (PLMVd) [62–65] and for viruses and viroids until now unrecorded in the country, i.e. Apple mosaic virus (ApMV), American plum line pattern virus (APLPV) and Hop stunt viroid (HSVd). Detection was carried out by nonradioactive molecular hybridization (MH) using the corresponding probes for the different viral/viroidal sequences [27, 66]. Reverse transcription-polymerase chain reaction (RT-PCR) was used on selected positive and negative samples, and in all the samples which showed questionable and weak reactions using MH. Infected and healthy fruits were collected from the selected positive and negative samples and analysed using the same detection methods described above.
Due to the absence of a unique peach tree variety infected with all the suggested combinations of virus/viroid for this study, we used three different peach tree varieties. Kawea: for healthy and PLMVd infected samples, White Lady: for PNRSV infected samples and Rosario Red for doubly infected samples. Before proceeding with the microarray analysis, we checked whether the distinct varieties could influence in the transcriptome analysis. For this purpose, we used Real Time qRT-PCR and amplified several genes in different infection scenarios against the same healthy variety. No influence of the different varieties used on the transcriptome profile was observed.
Total RNA from peaches for subsequent functional analyses was extracted as described by Meisel et al. . RNA quality and quantity was analyzed spectrophotometrically by evaluating the absorption ratios: A260/230 and A260/280 and confirmed with an EtBr stained 1.5% agarose gel containing 3% formaldehyde. Presence or absence of the different virus and viroids detected in above leaves was confirmed by MH. Positive tested samples infected either with PNRSV, PLMVd, or both pathogens, as well as healthy controls were selected for microarray profiling.
Microarray hybridization, scanning and data analyses
For each biological replicate equal amounts of samples from four peach fruits belonging to four different trees were pooled and RNA was extracted. To obtain robust and statistically accurate data, four biological replicates and two technical replicates (dye swaps: Cy3-healthy/Cy5-infected or Cy3-infected/Cy5-healthy) were used for all the different samples.
Amplification of RNA samples for microarray hybridization was carried out using the method of Van Gelder et al., 1990 . 1.25 μg of total RNA of each sample were amplified and aminoallyl-labelled using the MessageAmp® II aRNA kit (Ambion, http://www.ambion.com) and 5-(3-aminoallyl)-2´-deoxyuridine-5´-thiphosphate (aa-Dutp, Ambion), following manufacturer’s instructions. A rough amount of 20 μg of amplified RNA (aRNA) was obtained and 7.5 μg from each sample of aminoallyl-labelled aRNA were re-suspended in coupling buffer labelled with either Cy3 or Cy5 Mono NHS Ester (CyTM Dye Postlabelling Reactive Dye Pack, Amersham). Sample purification was carried out with MegaclearTM (Ambion) following manufacturer instructions. Incorporation of Cy3 and Cy5 was measured using 1 μl of the probe in a Nanodrop spectrophotometer (Nanodrop Technologies Inc.; http://www.nanodrop.com/).
Hybridization of samples and reference pool (RP) to the microarray slides was performed manually using Telechem Hybridization Chambers (Corning), according with the protocol described by Ogundiwin et al.  and manufacturer’s instructions. To detect differentially expressed genes in PNRSV, PLMVd or doubly infected samples in comparison to healthy samples, the data were analyzed with the SAM package (Significance Analysis of Microarray, . Statistical significance was assessed using SAM analyses with a false discovery rate of 5% with no fold change cut-off. To study the functional category distribution of genes with expression changes we used the Munich Information Center for Protein Sequences (MIPS FunCatDB, http://mips.helmholtz-muenchen.de/proj/funcatDB/). Gene lists were further analyzed with FatiGO to find differential distributions of gene ontology (GO) terms between statistically differential expressed genes and the rest of genes present on the microarray (Fisher’s exact test in 2 x 2 contingency tables), with P values adjusted after correcting for multiple testing.
Real-time qRT-PCR analysis
Nine genes found to be differentially expressed upon infection in the microarray experiments were selected for qRT-PCR analysis. qRT-PCR was performed with 100 ng total RNA using one step SYBR PrimeScript RT-PCR Kit II (Takara) following the manufacturer’s recommendations in an Applied Biosystems 7500 Real-Time PCR System. Each biological replicate was assayed in triplicate. Gene-specific oligonucleotide primers were designed using Primer Express® version 3.0 software (Applied Biosystems). Primer information is shown in Table 3. Expression levels for target genes were normalized to Elongation Factor 1-alpha (EF1 α) and fold expression changes compared to the healthy controls calculated using the DDthreshold cycle (Ct) method.
Total nucleic acids were first denatured with formaldehyde, serially diluted  and then applied to nylon membranes (positive charged, ROCHE), air dried and covalently UV cross-linked to the membrane (700 × 100 μJ/cm2). Prehybridizations and hybridizations with specific probes for the specific viruses and viroids were conducted as described previously [28, 70]. Chemiluminiscent detection with CSPD reagent as the substrate was performed as recommended by the manufacturer (ROCHE).
Pallas V, Garcia JA: How do plant viruses induce disease? Interactions and interference with host components. J Gen Virol 2011, 92: 2691-2705.
Whitham SA, Yang C, Goodin MM: Global impact: elucidating plant responses to viral infection. Mol Plant Microbe Interact 2006, 19: 1207-1215.
Havelda Z, Varallyay E, Valoczi A, Burgyan J: Plant virus infection-induced persistent host gene downregulation in systemically infected leaves. Plant J 2008, 55: 278-288.
Aranda M, Maule A: Virus-induced host gene shutoff in animals and plants. Virology 1998, 243: 261-267.
Whitham SA, Quan S, Chang HS, Cooper B, Estes B, Zhu T, Wang X, Hou YM: Diverse RNA viruses elicit the expression of common sets of genes in susceptible Arabidopsis thaliana plants. Plant J 2003, 33: 271-283.
Liu Y, Ren D, Pike S, Pallardy S, Gassmann W, Zhang S: Chloroplast-generated reactive oxygen species are involved in hypersensitive response-like cell death mediated by a mitogen-activated protein kinase cascade. Plant J 2007, 51: 941-954.
Hadidi A, Barba M: Economic impact of pome and stone fruit viruses and viroids. In Virus and Virus Like Diseases of Pome and Stone Fruits. Edited by: Hadidi A, Barba M, Candresse T, Jelkmann W. St Paul, MN: American Phytopathological Society; 2011:1-8.
Flores R, Delgado S, Rodio ME, Ambros S, Hernandez C, Serio FD: Peach latent mosaic viroid: not so latent. Mol Plant Pathol 2006, 7: 209-221.
Pallas V, Aparicio F, Herranz MC, Amari K, Sanchez-Pina MA, Myrta A, Sanchez-Navarro JA: Ilarviruses of Prunus spp.: A continued concern for fruit trees. Phytopathology 2012,102(12):1108-1120.
Rowland O, Jones JD: Unraveling regulatory networks in plant defense using microarrays. Genome Biol 2001,2(1):1001.1-1001.3.
Trinks D, Rajeswaran R, Shivaprasad PV, Akbergenov R, Oakeley EJ, Veluthambi K, Hohn T, Pooggin MM: Suppression of RNA silencing by a geminivirus nuclear protein, AC2, correlates with transactivation of host genes. J Virol 2005, 79: 2517-2527.
Senthil G, Liu H, Puram VG, Clark A, Stromberg A, Goodin MM: Specific and common changes in Nicotiana benthamiana gene expression in response to infection by enveloped viruses. J Gen Virol 2005, 86: 2615-2625.
Marathe R, Guan Z, Anandalakshmi R, Zhao H, Dinesh-Kumar SP: Study of Arabidopsis thaliana resistome in response to cucumber mosaic virus infection using whole genome microarray. Plant Mol Biol 2004, 55: 501-520.
Agudelo-Romero P, Carbonell P, de la Iglesia F, Carrera J, Rodrigo G, Jaramillo A, Perez-Amador MA, Elena SF: Changes in the gene expression profile of Arabidopsis thaliana after infection with Tobacco etch virus. Virol J 2008, 5: 92.
Itaya A, Matsuda Y, Gonzales RA, Nelson RS, Ding B: Potato spindle tuber viroid strains of different pathogenicity induces and suppresses expression of common and unique genes in infected tomato. Mol Plant Microbe Interact 2002, 15: 990-999.
Huang Z, Yeakley JM, Garcia EW, Holdridge JD, Fan JB, Whitham SA: Salicylic acid-dependent expression of host genes in compatible Arabidopsis-virus interactions. Plant Physiol 2005, 137: 1147-1159.
Rizza S, Conesa A, Juarez J, Catara A, Navarro L, Duran-Vila N, Ancillo G: Microarray analysis of Etrog citron (Citrus medica L.) reveals changes in chloroplast, cell wall, peroxidase and symporter activities in response to viroid infection. Mol Plant Pathol 2012,13(8):852-864.
Golem S, Culver JN: Tobacco mosaic virus induced alterations in the gene expression profile of Arabidopsis thaliana. Mol Plant Microbe Interact 2003, 16: 681-688.
Dardick C: Comparative expression profiling of Nicotiana benthamiana leaves systemically infected with three fruit tree viruses. Mol Plant Microbe Interact 2007, 20: 1004-1017.
Hull R: In Matthews’ Plant Virology. London: Edited by Academic Press; 2002.
Gonzalez-Jara P, Tenllado F, Martinez-Garcia B, Atencio FA, Barajas D, Vargas M, Diaz-Ruiz J, Diaz-Ruiz JR: Host-dependent differences during synergistic infection by Potyviruses with potato virus X. Mol Plant Pathol 2004, 5: 29-35.
Gonzalez-Jara P, Atencio FA, Martinez-Garcia B, Barajas D, Tenllado F, Diaz-Ruiz JR: A Single Amino Acid Mutation in the Plum pox virus Helper Component-Proteinase Gene Abolishes Both Synergistic and RNA Silencing Suppression Activities. Phytopathology 2005, 95: 894-901.
Vance VB: Replication of potato virus X RNA is altered in coinfections with potato virus Y. Virology 1991, 182: 486-494.
Garcia-Marcos A, Pacheco R, Martianez J, Gonzalez-Jara P, Diaz-Ruiz JR, Tenllado F: Transcriptional changes and oxidative stress associated with the synergistic interaction between Potato virus X and Potato virus Y and their relationship with symptom expression. Mol Plant Microbe Interact 2009, 22: 1431-1444.
Postnikova OA, Nemchinov LG: Comparative analysis of microarray data in Arabidopsis transcriptome during compatible interactions with plant viruses. Virol J 2012, 9: 101.
Zanchin A, Bonghi C, Casadoro G, Ramina A, Rascio N: Cell enlargement and cell separation during peach fruit development. International Journal of Plant Science 1994, 155: 49-56.
Herranz MC, Sanchez-Navarro JA, Aparicio F, Pallas V: Simultaneous detection of six stone fruit viruses by non-isotopic molecular hybridization using a unique riboprobe or ‘polyprobe’. J Virol Methods 2005, 124: 49-55.
Pallas V, Mas P, Sanchez-Navarro JA: Detection of plant RNA viruses by nonisotopic dot-blot hybridization. Methods Mol Biol 1998, 81: 461-468.
Lilly ST, Drummond RS, Pearson MN, MacDiarmid RM: Identification and validation of reference genes for normalization of transcripts from virus-infected Arabidopsis thaliana. Mol Plant Microbe Interact 2011, 24: 294-304.
Cosgrove JD: Expansive growth of plant cell walls. Plant Physiol Biochem 2000,38(1–2):109-124.
Tessitori M, Maria G, Capasso C, Catara G, Rizza S, De Luca V, Catara A, Capasso A, Carginale V: Differential display analysis of gene expression in Etrog citron leaves infected by Citrus viroid III. Biochim Biophys Acta 2007, 1769: 228-235.
Rizza S, Capasso C, Catara A, Capasso A, Conte E, Catara A Proceedings of the 17th Conference of the International Organization of Citrus Virologists-IOCV, pp. XVII. In Transcriptional response of Troyer citrange, sour orange and alemow rootstocks to Citrus viroid IIIb (CVd-IIIb) infection. Adana, Turkey: Conference of the International Organization of Citrus Virologists; 2010:142-149. http://www.ivia.es/iocv/
Owens RA, Tech KB, Shao JY, Sano T, Baker CJ: Global analysis of tomato gene expression during Potato spindle tuber viroid infection reveals a complex array of changes affecting hormone signaling. Mol Plant Microbe Interact 2012, 25: 582-598.
Ogundiwin EA, Marti C, Forment J, Pons C, Granell A, Gradziel TM, Peace CP, Crisosto CH: Development of ChillPeach genomic tools and identification of cold-responsive genes in peach fruit. Plant Mol Biol 2008, 68: 379-397.
Sánchez-Navarro JA FA, Rowhani A, Pallás V: Comparative analysis of ELISA, nonradioactive molecular hybridization and PCR for the detection of Prunus necrotic ringspot virus in herbaceous and prunus host. Plant Pathol 1998, 47: 780-786.
Astruc N, Marcos JF, Macquaire G, Candresse T, Pallas V: Studies on the diagnosis of hop stunt viroid in fruit trees: Identification of new hosts and application of a nucleic acid extraction procedure based on non-organic solvents. Eur J Plant Pathol 1996, 102: 837-846.
Myrta A, Di Terlizzi B, Pallas V, Savino V: Viruses and viroids of apricot in the Mediterranean: incidence and biodiversity. Acta Horticulturae 2006, 701: 409-417.
Bouzayen M, Latché A, Nath P, Pech JC: Mechanism of fruit ripening. In Plant Developmental Biology- Biotechnological Perspectives: Volume I Edited by: Pua EC, Darvey MR. 2010, 319-339. Chapter 16
Trainotti L, Bonghi C, Ziliotto F, Zanin D, Rasori A, Casadoro G, Ramina A, T P: The use of microarray mPEACH 1.0 to investigate transcriptome changes during transition from pre-climateric to climacteric phase in peach fruit. Plant Sci 2006, 170: 606-613.
Lombardo VA, Osorio S, Borsani J, Lauxmann MA, Bustamante CA, Budde CO, Andreo CS, Lara MV, Fernie AR MFD: Metabolic profiling during peach fruit development and ripening reveals the metabolic networks that underpin each developmental stage. Plant Physiol 2011,157(4):1696-1710.
Manganaris GA RA, Bassi D, Geuna F, Ramina A, Tonutti P, Bonghi C: Comparative transcript profiling of apricot (Prunus armeniaca L.) fruit development and on-tree ripening. Tree Genet Genomes 2011, 7: 609-616.
Uyemoto JK, Scott SW: Important diseases of Prunus caused by viruses and other graft-transmissible pathogens in California and South Carolina. Plant Dis 1992, 76: 5-11.
Li J, Yang H, Peer WA, Richter G, Blakeslee J, Bandyopadhyay A, Titapiwantakun B, Undurraga S, Khodakovskaya M, Richards EL, et al.: Arabidopsis H+-PPase AVP1 regulates auxin-mediated organ development. Science 2005, 310: 121-125.
Paponov IA, Paponov M, Teale W, Menges M, Chakrabortee S, Murray JA, Palme K: Comprehensive transcriptome analysis of auxin responses in Arabidopsis. Mol Plant 2008, 1: 321-337.
Padmanabhan MS, Goregaoker SP, Golem S, Shiferaw H, Culver JN: Interaction of the tobacco mosaic virus replicase protein with the Aux/IAA protein PAP1/IAA26 is associated with disease development. J Virol 2005, 79: 2549-2558.
Padmanabhan MS, Shiferaw H, Culver JN: The Tobacco mosaic virus replicase protein disrupts the localization and function of interacting Aux/IAA proteins. Mol Plant Microbe Interact 2006, 19: 864-873.
Padmanabhan MS, Kramer SR, Wang X, Culver JN: Tobacco mosaic virus replicase-auxin/indole acetic acid protein interactions: reprogramming the auxin response pathway to enhance virus infection. J Virol 2008, 82: 2477-2485.
Kuhn JM, Boisson-Dernier A, Dizon MB, Maktabi MH, Schroeder JI: The protein phosphatase AtPP2CA negatively regulates abscisic acid signal transduction in Arabidopsis, and effects of abh1 on AtPP2CA mRNA. Plant Physiol 2006, 140: 127-139.
Whenham RJ, Fraser RSS, Brown LP, Payne JA: Tobacco-mosaic-virus-induced increase in abscisic-acid concentration in tobacco leaves: Intracellular location in light and dark-green areas, and relationship to symptom development. Planta 1986, 168: 592-598.
Bari R, Jones JD: Role of plant hormones in plant defence responses. Plant Mol Biol 2009, 69: 473-488.
Kotchoni SO, Kuhns C, Ditzer A, Kirch HH, Bartels D: Over-expression of different aldehyde dehydrogenase genes in Arabidopsis thaliana confers tolerance to abiotic stress and protects plants against lipid peroxidation and oxidative stress. Plant Cell Environ 2006, 29: 1033-1048.
Mowla SB, Cuypers A, Driscoll SP, Kiddle G, Thomson J, Foyer CH, Theodoulou FL: Yeast complementation reveals a role for an Arabidopsis thaliana late embryogenesis abundant (LEA)-like protein in oxidative stress tolerance. Plant J 2006, 48: 743-756.
Amari K, Diaz-Vivancos P, Pallas V, Sanchez-Pina MA, Hernandez JA: Oxidative stress induction by Prunus necrotic ringspot virus infection in apricot seeds. Physiol Plant 2007, 131: 302-310.
Gilroy EM, Hein I, van der Hoorn R, Boevink PC, Venter E, McLellan H, Kaffarnik F, Hrubikova K, Shaw J, Holeva M, et al.: Involvement of cathepsin B in the plant disease resistance hypersensitive response. Plant J 2007, 52: 1-13.
Kruger J, Thomas CM, Golstein C, Dixon MS, Smoker M, Tang S, Mulder L, Jones JD: A tomato cysteine protease required for Cf-2-dependent disease resistance and suppression of autonecrosis. Science 2002, 296: 744-747.
Bernoux M, Timmers T, Jauneau A, Briere C, De Wit PJ, Marco Y, Deslandes L: RD19, an Arabidopsis cysteine protease required for RRS1-R-mediated resistance, is relocalized to the nucleus by the Ralstonia solanacearum PopP2 effector. Plant Cell 2008, 20: 2252-2264.
Shabab M, Shindo T, Gu C, Kaschani F, Pansuriya T, Chintha R, Harzen A, Colby T, Kamoun S, van der Hoorn RA: Fungal effector protein AVR2 targets diversifying defense-related cys proteases of tomato. Plant Cell 2008, 20: 1169-1183.
Van Esse HP, Van’t Klooster JW, Bolton MD, Yadeta KA, Van Baarlen P, Boeren S, Vervoort J, De Wit PJ, Thomma BP: The Cladosporium fulvum virulence protein Avr2 inhibits host proteases required for basal defense. Plant Cell 2008, 20: 1948-1963.
Song J, Win J, Tian M, Schornack S, Kaschani F, Ilyas M, van der Hoorn RA, Kamoun S: Apoplastic effectors secreted by two unrelated eukaryotic plant pathogens target the tomato defense protease Rcr3. Proc Natl Acad Sci U S A 2009, 106: 1654-1659.
Tian M, Win J, Song J, van der Hoorn R, van der Knaap E, Kamoun S: A Phytophthora infestans cystatin-like protein targets a novel tomato papain-like apoplastic protease. Plant Physiol 2007, 143: 364-377.
Rooney H, Van’t Klooster J, Van der Hoorn R, Joosten M, Jones J: Cladosporium Avr2 inhibits tomato Rcr3 protease required for Cf-2-dependent disease resistance. Science 2005, 308: 1783-1786.
Auger AJ: Tomato ringspot virus associated with brownline disease on prune trees in Chile. Acta Horticulturae 1989, 235: 197-204.
Herrera G: Enfermedades causadas por virus en frutales en Chile. Santiago, Chile: Instituto de Investigación Agropecuaria; 2001. Boletín INIA N°52. 65p
Fiore N, Abou Ghanem-Sabanadzovic N, Infante R, Myrta A, Pallás V: Detection of Peach latent mosaic viroid in stone fruits from Chile. In Option Méditerranéennes, Sér. B/n°45 –Virus ad virus-like disease of stone fruits, with particular reference to the Mediterranean region Edited by: Myrta A, Di Terlizzi B, Savino V. 2003, 143-145.
Torres H, Gómez G, Pallás V, Stamo B, Shalaby A, Aouane B, Gavriel I, Kominek P, Caglayan K, Sipahioglu M, et al.: Detection by tissue printing of stone fruit viroids, from europe, the mediterranean and north and south America. Acta Horticulturae 2004, 657: 379-383.
Peiró A, Pallás V, Sánchez-Navarro JA: Simultaneous detection of eight viruses and two viroids affecting stone fruit trees by using a unique polyprobe. Eur J Plant Pathol 2012,132(4):469-475.
Meisel L, Fonseca B, Gonzalez S, Baeza-Yates R, Cambiazo V, Campos R, Gonzalez M, Orellana A, Retamales J, Silva H: A rapid and efficient method for purifying high quality total RNA from peaches (Prunus persica) for functional genomics analyses. Biol Res 2005, 38: 83-88.
Van Gelder RN, Von Zastrow ME, Yool A, Dement WC, Barchas JD JHE: Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Natl Acad Sci U S A 1990,87(5):1663-1667.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A 2001, 98: 5116-5121.
Sanchez-Navarro JA, Canizares MC, Cano EA, Pallas V: Simultaneous detection of five carnation viruses by non-isotopic molecular hybridization. J Virol Methods 1999, 82: 167-175.
We thank Drs M.A. Pérez-Amador y J. Gadea for helping in the result analysis. This work was supported by grant BIO2011-25018 from the Spanish granting agency Direccion General de Investigacion Cientifica for the transcriptomic analyses and from the grant 2009CL0020 from the bilateral project INIA-Chile/CSIC-Spain for the phytosanitary evaluation. MC Herranz was the recipient of a contract from the Juan de la Cierva program of the Ministerio de Educación y Ciencia of Spain.
The authors declare that they have no competing interests.
Conceived and designed the experiments: MCH and VP. Performed the transcriptomic experiments: MCH. Performed the viral diagnosis: AG, NF, AZ, MCH and MR. Analyzed the data: MCH, AN and VP. Wrote the paper: MCH, AN and VP. All authors discussed the results and commented on the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Comparison between amounts of PLMVd and PNRSV in infected samples used for profiling. Average ± standard deviation (SD) of pixel intensity corresponding to PLMVd or PNRSV titer at a concentration of 0.3ng RNA (see Figure 2). Virus/viroid titer in samples infected with PLMVd (left graph) or PNRSV (right graph) was slightly lower compared to virus/viroid titer in samples infected with both pathogens simultaneously. (PPT 111 KB)
Additional file 4: Table S3: Common genes differentially expressed upon each infection scenario. According to the Venn diagrams (Figure 4) the table shows the description of the genes with altered expression which are common between the different combinations of infection. Between brackets the number of genes for each scenario is shown. (XLSX 13 KB)
Additional file 5: Table S4: Genes related with ripening process and with significantly altered expression upon the three different infections. Expansin putative (EXP1), At1g69530; Glutamate dehydrogenase 2 (GDH2), At5g07440; Expansin, putative (EXP8), At2g40610; Auxin-responsive protein / indoleacetic acid-induced protein 9 (IAA9), At5g65670; bZIP transcription factor family protein, similar to common plant regulatory factor 7, At4g34590; Glutamate dehydrogenase 2 (GDH2), At5g07440; Auxin-responsive protein / indoleacetic acid-induced protein 17 (IAA17), At1g04250; Invertase/pectin methylesterase inhibitor family protein, At5g62350; Agamous-like MADS box protein AGL8, At5g60910; Zinc finger (C3HC4-type RING finger) family protein, At5g25560. (PPTX 63 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.