Analysis of gene expression changes in peach leaves in response to Plum pox virus infection using RNA-Seq

Differences in gene expression were studied after Plum pox virus (PPV, sharka disease) infection in peach GF305 leaves with and without sharka symptoms using RNA-Seq. For each sample, more than 80% of 100-nucleotide paired-end (PE) Illumina reads were aligned on the peach reference genome. In the symptomatic sample, a significant proportion of reads were mapped to PPV reference genomes (1.04% compared with 0.00002% in non-symptomatic leaves), allowing for the ultra-deep assembly of the complete genome of the PPV isolate used (9775 nucleotides, missing only 11 nucleotides at the 5' genome end). In addition, significant alternative splicing events were detected in 359 genes and 12 990 single nucleotide polymorphisms (SNPs) were identified, 425 of which could be annotated. Gene ontology annotation revealed that the high-ranking mRNA target genes associated with the expression of sharka symptoms are mainly related to the response to biotic stimuli, to lipid and carbohydrate metabolism and to the negative regulation of catalytic activity. A greater number of differentially expressed genes were observed in the early asymptomatic phase of PPV infection in comparison with the symptomatic phase. These early infection events were associated with the induction of genes related to pathogen resistance, such as jasmonic acid, chitinases, cytokinin glucosyl transferases and Lys-M proteins. Once the virus had accumulated, the overexpression of Dicer protein 2a genes suggested a gene silencing plant response that was suppressed by the virus HCPro and P1 proteins. These results illustrate the dynamic nature of the peach-PPV interaction at the transcriptome level and confirm that sharka symptom expression is a complex process that can be understood on the basis of changes in plant gene expression.