Transcriptomic profiling of Melon necrotic spot virus-infected melon plants revealed virus strain and plant cultivar-specific alterations

Cristina Gómez-Aix, Laura Pascual, Joaquín Cañizares, María Amelia Sánchez-Pina, Miguel A. Aranda

Viruses are among the most destructive and difficult to control plant pathogens. Melon (Cucumis melo L.) has become the model species for the agriculturally important Cucurbitaceae family. Approaches that take advantage of recently developed genomic tools in melon have been extremely useful for understanding viral pathogenesis and can contribute to the identification of target genes for breeding new resistant cultivars. In this work, we have used a recently described melon microarray for transcriptome profiling of two melon cultivars infected with two strains of Melon necrotic spot virus (MNSV) that only differ on their 3'-untranslated regions. Melon plant tissues from the cultivars Tendral or Planters Jumbo were locally infected with either MNSV-Mα5 or MNSV-Mα5/3'264 and analysed in a time-course experiment. Principal component and hierarchical clustering analyses identified treatment (healthy vs. infected) and sampling date (3 vs. 5 dpi) as the primary and secondary variables, respectively. Out of 7566 and 7074 genes deregulated by MNSV-Mα5 and MNSV-Mα5/3'264, 1851 and 1356, respectively, were strain-specific. Likewise, MNSV-Mα5/3'264 specifically deregulated 2925 and 1618 genes in Tendral and Planters Jumbo, respectively. The GO categories that were significantly affected were clearly different for the different virus/host combinations. Grouping genes according to their patterns of expression allowed for the identification of two groups that were specifically deregulated by MNSV-Mα5/3'264 with respect to MNSV-Mα5 in Tendral, and one group that was antagonistically regulated in Planters Jumbo vs. Tendral after MNSV-Mα5/3'264 infection. Genes in these three groups belonged to diverse functional classes, and no obvious regulatory commonalities were identified. When data on MNSV-Mα5/Tendral infections were compared to equivalent data on cucumber mosaic virus or watermelon mosaic virus infections, cytokinin-O-glucosyltransferase2 was identified as the only gene that was deregulated by all three viruses, with infection dynamics correlating with the amplitude of transcriptome remodeling. Strain-specific changes, as well as cultivar-specific changes, were identified by profiling the transcriptomes of plants from two melon cultivars infected with two MNSV strains. No obvious regulatory features shared among deregulated genes have been identified, pointing toward regulation through differential functional pathways.