RNA-Seq profile of flavescence dorée phytoplasma in grapevine

Abstract

Background: The phytoplasma-borne disease flavescence dorée is still a threat to European viticulture, despite mandatory control measures and prophylaxis against the leafhopper vector. Given the economic importance of grapevine, it is essential to find alternative strategies to contain the spread, in order to possibly reduce the current use of harmful insecticides. Further studies of the pathogen, the vector and the mechanisms of phytoplasma-host interactions could improve our understanding of the disease. In this work, RNA-Seq technology followed by three de novo assembly strategies was used to provide the first comprehensive transcriptomics landscape of flavescence dorée phytoplasma (FD) infecting field-grown Vitis vinifera leaves.

Results: With an average of 8300 FD-mapped reads per library, we assembled 347 sequences, corresponding to 215 annotated genes, and identified 10 previously unannotated genes, 15 polycistronic transcripts and three genes supposedly localized in the gaps of the FD92 draft genome. Furthermore, we improved the annotation of 44 genes with the addition of 5'/3' untranslated regions. Functional classification revealed that the most expressed genes were either related to translation and protein biosynthesis or hypothetical proteins with unknown function. Some of these hypothetical proteins were predicted to be secreted, so they could be bacterial effectors with a potential role in modulating the interaction with the host plant. Interestingly, qRT-PCR validation of the RNA-Seq expression values confirmed that a group II intron represented the FD genomic region with the highest expression during grapevine infection. This mobile element may contribute to the genomic plasticity that is necessary for the phytoplasma to increase its fitness and endorse host-adaptive strategies.

Conclusions: The RNA-Seq technology was successfully applied for the first time to analyse the FD global transcriptome profile during grapevine infection. Our results provided new insights into the transcriptional organization and gene structure of FD. This may represent the starting point for the application of high-throughput sequencing technologies to study differential expression in FD and in other phytoplasmas with an unprecedented resolution.

Figure 1

Distribution of assembled transcripts into functional categories according to KEGG classification. Number of expressed genes (y-axis) associated to KEGG pathways (x-axis).

Figure 2

Enrichment analysis of GO terms calculated by Fisher’s exact test. Reference set (red) is represented by all the FD expressed genes, whereas the test set (blue) is represented by the top 20 highly expressed FD genes.

Figure 3

Correlation between RNA-Seq and qRT-PCR expression data. Plot of the transcriptional levels of the 15 selected genes expressed as RPK values (y-axis) and Expression Index values (x-axis). R² = 0.98 (P < 0.05).