Transcriptomic profiling of Melon necrotic spot virus-infected melon plants revealed virus strain and plant cultivar-specific alterations
© The Author(s). 2016
Received: 10 November 2015
Accepted: 25 May 2016
Published: 7 June 2016
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.
KeywordsMNSV CMV Cucurbits Cytokinin-O-glucosyltransferase Resistance WMV
Viruses change the physiology and metabolism of infected plants, and can directly or indirectly influence the host’s gene expression patterns [1–4]. High-throughput technologies such as expressed sequence tags (ESTs), microarrays and next-generation sequencing, have made possible the simultaneous analysis of functional data for many genes and the study of the plant’s transcriptomic remodelling in response to virus infections . In the past few years, DNA microarrays have become popular tools for comparative high-throughput gene expression analysis, and microarray platforms have become available for both model and non-model crop species. Melon (Cucumis melo L.), in addition to its agronomic importance, has biological features that make it an interesting experimental model, favouring the development of a growing number of genetic and molecular tools for this species, including large ESTs collections [6, 7], TILLING platforms [8, 9] and the sequencing of its genome . More specifically, EST sequencing has allowed the development of a melon-specific microarray , which has been used for transcriptomic profiling of Cucumber mosaic virus (CMV), Watermelon mosaic virus (WMV) and Monosporacus cannonballus-infected plants [11–13]. In this work, we have used the melon microarray to profile the melon transcriptome after infection with Melon necrotic spot virus (MNSV).
MNSV (genus Carmovirus, family Tombusviridae) is endemic in cucurbit crops worldwide, often causing significant economic losses due to epidemic outbreaks. The MNSV genome is composed of a 4.3Kb, single-stranded positive-sense RNA containing at least five open reading frames (ORFs)  which are flanked by two untranslated regions (UTRs) at their 5′ and 3′ termini. The 3′ ORF encodes the capsid protein (CP) which has a structural role, is necessary for vascular transport of the virus, plays a role in suppression of RNA silencing  and is involved in virus transmission . It also contains a double gene block (DGB), typical of carmoviruses, consisting of two small, centrally located ORFs, which encode two consecutive 7 kDa proteins (p7A and p7B) involved in the cell-to-cell movement of the virus [15, 17]. The 5′ ORF can either encode a 29 kDa protein (p29) ending in an amber codon, or a larger 89 kDa gene product (p89) if it is read-through, which contains the RNA-dependent RNA polymerase (RdRp) domain. The p29 and p89 proteins are involved in viral replication [14, 15], which takes place in virus-altered mitochondria . The 3′-UTR of genomic MNSV RNAs, which are identical to those of sub-genomic RNAs, has been shown to contain sequences that act as cap-independent translational enhancers (3′-CITEs) [19, 20]. Depending on the specific nature of these 3′-CITEs, MNSV can infect N. benthamiana and the otherwise resistant melon plants that carry the recessive eIF4E 228L allele at the nsv locus [20–22]. Interestingly, 3′-CITEs exhibit a modular nature, as they can be exchanged among viral strains or even viral species through recombination [19, 23].
In this work, we have used two MNSV strains that only differed in their 3-UTRs, namely, MNSV-Mα5 and a chimera with its 3′-UTR from MNSV-264 (MNSV-Mα5/3’264) for infection profiling. MNSV-264 is a strain that is able to break the resistance controlled by nsv [21, 23]. The characterization of melon cultivar-specific responses was also investigated, and two melon cultivars were used for this purpose. These were: cv. Tendral, which is fully susceptible to MNSV, and cv. Planters Jumbo, which is homozygous for the recessive eIF4E 228L resistance allele and is therefore resistant to most MNSV isolates but not to those carrying the MNSV-264 3′-CITE , as is the case for MNSV-Mα5/3’264 . Locally-infected tissues were analysed in a time-course experiment and the melon microarray  was used for describing differential alterations of the melon transcriptome associated with: (i) the presence of one or another 3′-UTR in the MNSV RNA, (ii) the melon genotype and (iii) the type of plant tissue infected.
Identification of the main sources of variability
According to these results, MNSV-Mα5 induced faster and more marked changes in Tendral as compared to MNSV-Mα5/3’264, an effect that could also be seen in inoculated leaves. Among cultivars, the course of the infection resulted in greater differentiation among the samples of the different cultivars inoculated with the same virus.
Transcriptomic remodeling in inoculated cotyledons
Progression of MNSV accumulation
Differentially expressed genes during infection progression
Virus-specific transcriptomic alterations
Cultivar-specific transcriptomic alterations
Comparison of changes induced by MNSV with those induced by Cucumber mosaic virus (CMV) or Watermelon mosaic virus (WMV)
Functional analysis of the genes shared by MNSV-Mα5 and CMV did not identify categories that were statistically significant. Nevertheless, manual exploration of the 318 shared genes identified transcripts that were annotated as WRKY transcription factors, peroxidases, β-1, 3 glucanases, lipoxigenases, histones or heat-shock proteins, among others (Additional file 5). A comparison of functional categories deregulated by each virus showed important differences among them. MNSV-Mα5 deregulated a great number of GO categories, with emphasis on processes linked to peroxidation, carbohydrate metabolism and responses to various types of stress. However, on the list of CMV-deregulated genes, only two GO categories were found, among them, “sequence-specific DNA binding transcription factor activity” (GO:0003700), which includes a great number of transcription factors related to hormonal response regulation (Additional file 5). Among these transcription factors, a few of the Basic Leucine Zipper Domain (bZIP) type, involved in the mediation of the response to salicylic acid, were activated, while many WRKY and Ethylene response factors (ERF) were repressed (Additional file 5), evidencing the early regulation on the hormonal responses that CMV exerts. Among the genes regulated by WMV, we identified processes related to oxidative stress as being over-represented (Additional file 5).
Comparison of transcriptomic changes induced by MNSV in directly-inoculated leaves vs cotyledons
In this work, we have compared transcriptomic profiles of melon plants from two different cultivars infected with two different strains of MNSV. Melon cultivars differed on their alleles at the nsv locus, which encodes the eIF4E gene that controls susceptibility to most MNSV strains. The viral strains differed on their 3′-UTRs, which have been shown to control translational efficiency of MNSV RNAs and, hence, resistance-breaking capabilities.
Comparison of profiles associated with two strains of MNSV that differed in their 3-UTRs
Our results showed that accumulation of MNSV-Mα5/3’264 was lower to that of MNSV-Mα5 in Tendral tissues, in agreement with previous observations [19, 20, 23]. As the accumulation levels of viral RNAs have been related to the intensity of the transcriptomic changes induced by certain viruses [29, 30], the reduced transcriptomic impact induced by MNSV-Mα5/3’264 at 3 dpi (Fig. 3b) could be thus explained. This is an important aspect, because other observations with both strains at discrete time points could be, at least partially, due to differences in their infection dynamics. For this reason we decided to do our study at two different time points and compare the expression patterns as a function of time for both strains.
Both viruses activated genes involved in defense responses, oxidative stress and secondary metabolism or ubiquitin-dependent catabolic processes. Primary metabolism processes, mainly photosynthesis and genes related to the chloroplasts and the organization of the cell wall, were over-represented among the genes inhibited by both viruses, leading to changes that have consistently been described for other viruses [5, 31–33]. However, there were striking differences that were mainly found in two groups of genes that had a level of activation or inhibition in response to MNSV-Mα5/3’264 that were well above the levels of those induced by MNSV-Mα5 infection. Among these genes, various transcripts of unknown unigenes (sSSH1G12_c) were found, as well as transcripts that coded for resistance-related proteins (cCL2022Contig2), proteins involved in lipid transfer (cCL5847Contig1), cytochrome p450 CYP87A3 (cCL2810Contig1), or membrane steroid binding proteins (c46d_14-E05-M13R_c), all belonging to a diversity of functional classes. These differences suggested differential interactions of both viral strains with the corresponding factors of the host, or at least with the routes and/or processes where these factors were involved. It is important to note that the region that was exchanged between the two viral strains does not code for any protein, but contains RNA structural elements that have critical regulatory functions in a variety of viral processes, including translation, replication and transcription of sub-genomic RNAs [19–21, 23, 34]. The potential functions of the altered genes for either the virus or the host are unknown, as are the regulatory mechanisms that control the expression of these genes. In this regard, our attempts to identify common sequences shared by transcripts and viral 3′-UTRs that may point toward regulation through small RNAs have not been successful.
Comparison of profiles associated with different varieties of melon
In agreement to previous data [20, 23], MNSV-Mα5/3’264 reached higher levels of accumulation in Planters Jumbo than in Tendral. However, the amplitude of the transcriptomic response of Tendral was greater than in Planters Jumbo. The genes deregulated only in Tendral were mostly related with defense processes and functions involved in photosynthesis in general and photosystem II specifically. The decrease in photosynthetic activity has been related to defense mechanisms through the production of reactive oxygen species (ROS) derived from the chloroplasts [35, 36], and this has also been related to inhibition of photosystem II proteins with the increase in concentration of specific viruses, suggesting their involvement in basal defense responses . In this sense, the inhibition of these processes in Tendral and in Planters Jumbo could be related to the launching of different basal defense responses by the different cultivars. The list of genes that were specifically deregulated in Planters Jumbo, however, was enriched in functions related to potassium ion transport, which could be involved in early signaling of events that take place during viral infection in this cultivar .
Planters Jumbo is a cultivar that is susceptible to MNSV-264 and MNSV-Mα5/3’264 but resistant to other strains [23, 39], with the resistance being due to a mutation in eIF4E . The two cultivars tested here differ in this resistance gene, but presumably also in several other genes, as they are not isogenic lines. Thus, perhaps not surprisingly, both cultivars responded differently to infection by MNSV-Mα5/3’264. In this regard, the most remarkable finding was the unequivocal identification of a set of genes that showed antagonistic expression tendencies between both cultivars. Notably, clues of cultivar-viral strain double interactions were found: Among the deregulated genes, many of them coincided with those inhibited in Tendral that did not suffer changes in plants infected with MNSV-Mα5, which suggested the importance of those genes in the infection processes by MNSV-Mα5/3’264 and its differential behavior depending on the melon genotype. Among the different genes included on the list, we found Myb factors (cCL4076Contig1), which are involved, together with WRKY transcription factors, in the modulation of the plant’s hormonal expression. These factors are frequently modified as a consequence of viral infections in relation to the alteration of hormonal expression of the plant in favor of the virus or as an integrated part of the plant’s defense system . Other annotated genes included diverse protein kinases, which are an essential part of the signaling events required during defense responses, as well as related to cellular death associated to resistance [41–43]. The unigene cCL555Contig1, annotated as coding for the L31 protein of the 60S ribosome subunit, showed the greatest activation in Planters Jumbo, and was found among those that were most inhibited in Tendral, becoming a very interesting candidate for the analysis of its implication in the MNSV cycle of infection . Altogether, these genes represent potential targets for functional studies during infection with MNSV-Mα5/3’264 and point to the differential involvement and regulation of metabolic processes between both cultivars. Further research could include transcriptomic profiling of melon isogenic lines that only differ on eIF4E after MNSV-Mα5/3’264 infection.
Comparison of profiles associated with different melon tissues
A possible criticism of the generic analysis of transcriptomic profiles could come from the assumption that the different cells or tissues respond in similar ways to viral infections, without taking into account specific alterations in the tissue as well as space-time variations of lesser intensity [30, 45]. In this sense, the transcriptomic deregulation induced by MNSV in leaves was lesser than that in cotyledons. This was probably related to the progression of infection, which was slower in leaves than in cotyledons as shown by the levels of viral accumulation. However, there was a big overlap in the nature of the genes deregulated in leaves with those from the cotyledons. Likewise, the biological functions and metabolic processes among the deregulated genes in infected leaves mostly identified the same over-represented functions in both of these tissues. These results provide validity to the general view of MNSV infections obtained in melon cotyledons.
Comparison of profiles associated to infection by three different viruses
This comparison was possible due to the existence of previous data from research on infections of Tendral melon plants by CMV and WMV [11, 12]. As each data set come from independent experiments and with the objective of minimizing variations that could be attributed to each process of analysis, the data sets were normalized and analyzed separately following the same criteria for all three cases. After the analysis, a single gene was identified, annotated as cytokinin-O-glucosyltransferase2, which was deregulated by all three viruses, but with an expression pattern that differed according to the virus, that is, activated by MNSV and WMV but inhibited by CMV. Although the exact roles that cytokynins may play in plant-pathogen interactions are unknown, in Arabidopsis the involvement of these hormonal routes in some responses mediated by R proteins have been identified . The common deregulation of this gene by the different viruses could indicate its involvement in the response to viruses in melon plants .
Among MNSV and CMV, however, numerous shared genes were identified, among them, genes related to response to stress and general defense, heat shock proteins (cCL5861Contig1), glutathione S-transferases (cA_23-D09-M13R_c), transcripts that code for resistance proteins (cCL1320Contig1), or many WRKY transcription factors, which have been identified as a response to various viruses in other hosts [5, 48]. However, each virus-host interaction was unique in terms of modified biological functions, as well as in the levels of genetic deregulation. For example, CMV caused a rapid deregulation of genes related to hormonal routes, while MNSV induced a rapid defense response and the activation of oxidative stress routes in infected plants. Unlike the other two viruses, the amplitude of the transcriptomic response induced by WMV was small, only activating oxidative-stress genes. A hypothesis that could explain the scarce genetic deregulation by WMV could be that the degree of response was related to the levels of viral accumulation. The quantification by RT-qPCR of a group of stress-response genes as a function of time showed the existence of correlation between the viral accumulation and the values of gene deregulation. These results suggest the existence of potential methodological errors when performing comparative analysis that use a single time point of observation or narrow temporal windows to compare different infections. Besides the dynamics of viral accumulation, the common deregulation of specific genes by different viruses has also been related to the phylogenetic distance of the viruses compared , so that new analysis and comparisons through the use of wider temporal windows and genetically-related viruses could provide more information on the processes that are commonly deregulated as a response to viral infection in melon plants.
By comparing transcriptomic profiles of plants from the same cultivar infected with each of the two viral strains, we have shown that there are common but also strain-specific changes, the latter referring to a variety of genes with very different functions that were affected. No obvious regulatory features shared among deregulated genes were identified. Similarly, by comparing transcriptomic profiles of plants from each of the two cultivars infected with the same viral strain, common but also cultivar-specific changes were identified. Again, no obvious features among deregulated genes arose, but our analysis suggested the launching of different basal defense responses resulting in differential involvement of hormonal and stress response processes. An important methodological aspect emerging from this work is the influence of infection dynamics in transcriptome profiling. When comparing different viruses or viral strains, observations done at single time points could be significantly influenced by different infection dynamics. Biologically relevant data can be obtained by performing observations at different time points and comparing expression patterns as a function of time.
Plant material, viral isolates and virus inoculation
Melon (Cucumis melo L.) seeds from the cultivars Tentral (Fitó Seeds, Barcelona, Spain) and Planters Jumbo (accession C46 from the Experimental Station of “La Mayora”-CSIC, Malaga, Spain) were used. The seeds were germinated in Petri dishes for 48 h at 25 °C. After germination, the seedlings were transplanted onto 35-cell trays with soil, and grown at 27/19 °C day/night conditions. Mechanical inoculations were done on fully-expanded cotyledons after 7 days in the greenhouse or on the first fully-expanded leaves after 15 days. For these inoculations, a mix of fresh inoculum in a 30 mM potassium phosphate buffer (pH 8.0) with active charcoal and 0.037 mm Carborundum particles was used. The inocula used were from infected melon plants which had been infected with lyophilized material from MNSV-Mα5 and MNSV-Mα5/3’264 .
Experimental design and sampling
The biological assay consisted of cv. Tendral plants inoculated either with MNSV-Mα5 or MNSV-Mα5/3’264, and Planters Jumbo plants inoculated with MNSV-Mα5/3’264. For the healthy controls, plants from each cultivar were treated with virus-free buffer as used for the inoculations. The sampling was done at 1, 3 and 5 days post-inoculation (dpi) in cotyledons and at 5dpi in inoculated leaves. We used three biological replicates that were composed of a pool of three different plants. Therefore, a single treatment (virus/cultivar/sampling time point) included at least 9 different plants. Viral accumulation was quantified in cotyledon samples inoculated with the different viruses (MNSV-Mα5, CMV, WMV) at 1, 3, 5 and 9 dpi. At 15 dpi, the second leaf showing systemic infection from plants inoculated with CMV-fny  and WMV-M116  was also sampled.
RNA extraction and microarray hybridization
All the samples were independently harvested and frozen in liquid N2 and stored at -80 °C. The RNA extractions were performed with Tri-Reagent (Sigma-Aldrich, St. Louis), according to the manufacturer’s instructions. After the extraction, the RNA was analyzed by dot-blot to check for the presence of the virus in the infected samples. To eliminate traces of genomic DNA, total RNA was incubated with DNAse I (New England Biolabs, London) for 10 min at 37 °C. The reaction volume was adjusted to 100 μl, and the aqueous phase was extracted with phenol/chloroform/isoamyl alcohol (25:24:1). Lastly, the RNA was precipitated with 10 % (v/v) NaOAc (3 M) and 2.5 volumes of absolute alcohol by centrifugation (12,000 x g, 20 min at 4 °C). The quality and quantity of RNA was verified with a ND-1000 spectrometer (Nano Drop Technologies, Wilmington, DE, USA) and a Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). For this work we used the melon microarray  adding 244 new unigenes , already used in other research studies [12, 13, 52]. In total, the melon microarray contain 17,443 melon unigenes that represent 10,649 genes in the melon genome, with additional 2,021 unigenes with no assigned hit in the annotated melon genome (Additional file 7: Table S7) . Hybridizations were performed by NimbleGen’s microarray hybridization service (IRB Functional Genomics Core, BaldiriReizac, 10-12, 08028, Barcelona, Spain).
The hybridization data were obtained from our experiment and from public repositories. Our experiment data provided by NimbleGen, were grouped for normalization into two groups: cotyledon and leaf. As the CMV data came from an older platform, each hybridization data group (CMV, WMV and MNSV-α5 at 3dpi) was independently normalized for later analysis and comparison. Each group of data were normalized and transformed to a log2 scale using the RMA (Robust Multi-array Average) algorithm found in the oligo package  of Bioconductor (http://www.bioconductor.org). For the cotyledon time-course experiment, the maSigPro package  was used to identify differentially-expresed genes. This program uses a two-regression step strategy. In the first step, a general regression model is defined. Then, the defined model is adjusted to the data through least squares, and the genes that significantly differ from this regression model are identified by correcting with a specific false discovery rate (FDR) of 1 % (Q = 0.01). In the second step, a stepwise regression is employed, and a probability (p) is calculated for each variable, showing the probability that causes the deviation. After the analysis, a list of the differentially expressed genes is obtained according to each variable (“TIME”, “TIME x Virus” and “Virus vs. Control”). We discard the differentially expressed genes associated with only the variable “TIME” in order to select those genes that were deregulated with time and associated with the virus in each cultivar. A maSigPro analysis was conducted for each cultivar. For the single stage experiment, the identification of differentially expressed genes was done through the SAM (Significant analysis of microarrays) module  found in the Multi Experimental Viewer (MeV, v. 4.9.0) program , using a FDR = 0. Genes with a fold change smaller than 2 (FC ≤2, cut-off of log2 ≤ 1) were filtered out.
Samples were grouped with the PCA module from MeV . Clustering of the samples was done with Euclidean distance by hierachical clustering  and the bootstrap was done by Support trees  (bootstrap 100 replicates). Genes were clustered by their expression pattern by using the k-means clustering method  and Pearson’s correlation for the calculation of distances. Lastly, the functional analysis was done with the Blast2GO program , extracting the over- or under-represented GO terms among the differentially-expressed genes from each condition by the application of Fisher’s test (p-value <0.05).
Microarray validation and real time quantitative reverse transcription PCR (RT-qPCR)
The melon microarray was validated in previous works [11–13] and we undertook further verification by comparing microarray and RT-qPCR expression patterns of a pathogen response protein, a calmodulin-binding protein, a lipoxigenase and a glucosyl transferase transcript. The same RNA samples from cotyledon hybridized to the microarray were used for this purpose. Data from RT-qPCR were transformed to a log2 scale to make the data comparable with microarray results. A strong positive correlation was found between the two sets of values (R 2 = 0.86; correlation coefficient of 0.93) (Additional file 8), confirming previous results [12, 13].
For real time quantitative PCR, the first strand cDNA was synthesized using 1.5 μg of total RNA, following the directions of the reverse transcriptase manufacturer (Roche) with an oligo-dT(16) as reverse primer. As MNSV and CMV do not have a poly(A) tail, reverse primers for the respective viruses (CE-948, 5′-CCCACTATCATCACGATCTTTAC-3′, and CE-169, 5′-CCGCTTACGATTCCCAACTGT-3′) were added for transcription of the viral RNAs. The qPCR for the quantification of messenger RNA, as well as viral accumulation was performed on an AB7500 System (Applied Biosystems), using SYBR Green PCR Master Mix (Applied Biosystems) as the detector and ROX as the passive reference. All the reactions (final volume of 20 μl) contained 10 μl Master Mix, 0.15 μl of each primer (100 mM) and 60 ng of cDNA. Each reaction was done in triplicate, along with controls without DNA (NTC), using a two-step amplification protocol and adding a melting curve. The analysis of the melting curves and the NTC were done in order to ensure the specific amplification of the product and the absence of dimerization of the primers. The primers used for amplification of the target and reference genes are listed in Additional file 9.
For calculating the relative quantification of each transcript, we used the 2ΔΔct method. The relative expression levels were determined through the normalization of the samples with mRNA from cyclophilin (cCL3169Contig1) as an internal control and relating it to the expression values of the healthy controls. The analysis was carried out with the SDS-7500 software and exported to a spread sheet for further calculations. The specific primer pairs were designed with Primer Express software v3.0 (Applied Biosystems). The efficiency of each primer pair was calculated through the equation: Efficiency (%) = (10[-1/slope] - 1) x 100 (Guide to performing relative quantitation of gene expression using real-time quantitative PCR, Applied Biosystems).
We thank Mari Carmen Montesinos and Blanca Gosalvez for their excellent technical assistance. M. Fon (email@example.com) edited the manuscript.
This work was supported by grants AGL2012-37390 and PCIN-2013-043 (Ministerio de Economía y Competitividad, Spain)
Availability of data and materials
The data reported in this publication have been deposited in NCBI’s Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE74840 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE74840).
CGA performed the experiments. CGA, LP and JC analysed the data. MASP and MAA conceived the study. CGA, MASP and MAA wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Maule A, Leh V, Lederer C. The dialogue between viruses and hosts in compatible interactions. Curr Opin Plant Biol. 2002;5(4):279–84.View ArticlePubMedGoogle Scholar
- Pallás V, Garcia JA. How do plant viruses induce disease? Interactions and interference with host components. J Gen Virol. 2011;92(Pt 12):2691–705.View ArticlePubMedGoogle Scholar
- Wang D, Maule AJ. Inhibition of host gene expression associated with plant virus replication. Science. 1995;267(5195):229–31.View ArticlePubMedGoogle Scholar
- Aranda MA, Escaler M, Wang D, Maule AJ. Induction of HSP70 and polyubiquitin expression associated with plant virus replication. Proc Natl Acad Sci U S A. 1996;93(26):15289–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Whitham SA, Yang CL, Goodin MM. Global impact: Elucidating plant responses to viral infection. Mol Plant Microbe Interact. 2006;19(11):1207–15.View ArticlePubMedGoogle Scholar
- Clepet C, Joobeur T, Zheng Y, Jublot D, Huang M, Truniger V, et al. Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. BMC Genomics. 2011;12:252.View ArticlePubMedPubMed CentralGoogle Scholar
- González-Ibeas D, Blanca J, Roig C, Gonzalez-To M, Pico B, Truniger V, et al. MELOGEN: an EST database for melon functional genomics. BMC Genomics. 2007;8:306.View ArticlePubMedPubMed CentralGoogle Scholar
- Dahmani-Mardas F, Troadec C, Boualem A, Leveque S, Alsadon AA, Aldoss AA, et al. Engineering melon plants with improved fruit shelf life using the TILLING approach. Plos One. 2010;5(12):e15776.View ArticlePubMedPubMed CentralGoogle Scholar
- González M, Xu M, Esteras C, Roig C, Monforte AJ, Troadec C, et al. Towards a TILLING platform for functional genomics in Piel de Sapo melons. BMC Res Notes. 2011;4:289.View ArticlePubMedPubMed CentralGoogle Scholar
- García MJ. The genome of melon (Cucumis melo L.). Proc Natl Acad Sci U S A. 2012;109(29):11872–7.View ArticleGoogle Scholar
- Mascarell-Creus A, Cañizares J, Vilarrasa-Blasi J, Mora-García S, Blanca J, González-Ibeas D, et al. An oligo-based microarray offers novel transcriptomic approaches for the analysis of pathogen resistance and fruit quality traits in melon (Cucumis melo L.). BMC Genomics. 2009;10:467.View ArticlePubMedPubMed CentralGoogle Scholar
- González-Ibeas D, Canizares J, Aranda MA. Microarray analysis shows that recessive resistance to Watermelon mosaic virus in melon is associated with the induction of defense response genes. Mol Plant Microbe Interact. 2012;25(1):107–18.View ArticlePubMedGoogle Scholar
- Roig C, Fita A, Rios G, Hammond JP, Nuez F, Pico B. Root transcriptional responses of two melon genotypes with contrasting resistance to Monosporascus cannonballus (Pollack et Uecker) infection. BMC Genomics. 2012;13:601.View ArticlePubMedPubMed CentralGoogle Scholar
- Riviere CJ, Rochon DM. Nucleotide-Sequence and genomic organization of melon necrotic spot virus. J Gen Virol. 1990;71:1887–96.View ArticlePubMedGoogle Scholar
- Genovés A, Navarro JA, Pallás V. Functional analysis of the five melon necrotic spot virus genome-encoded proteins. J Gen Virol. 2006;87(Pt 8):2371–80.View ArticlePubMedGoogle Scholar
- Ohki T, Akita F, Mochizuki T, Kanda A, Sasaya T, Tsuda S. The protruding domain of the coat protein of Melon necrotic spot virus is involved in compatibility with and transmission by the fungal vector Olpidium bornovanus. Virology. 2010;402(1):129–34.View ArticlePubMedGoogle Scholar
- Navarro JA, Genovés A, Climent J, Sauri A, Martinez-Gil L, Mingarro I, et al. RNA-binding properties and membrane insertion of Melon necrotic spot virus (MNSV) double gene block movement proteins. Virology. 2006;356(1-2):57–67.View ArticlePubMedGoogle Scholar
- Gómez-Aix C, García-García M, Aranda MA, Sánchez-Pina MA. Melon necrotic spot virus Replication Occurs in Association with Altered Mitochondria. Mol Plant Microbe Interact. 2015;28(4):387–97.View ArticlePubMedGoogle Scholar
- Miras M, Sempere RN, Kraft JJ, Miller WA, Aranda MA, Truniger V. Interfamilial recombination between viruses led to acquisition of a novel translation-enhancing RNA element that allows resistance breaking. New Phytologist. 2014;202(1):233–46.View ArticlePubMedGoogle Scholar
- Truniger V, Nieto C, Gonzalez-Ibeas D, Aranda M. Mechanism of plant eIF4E-mediated resistance against a Carmovirus (Tombusviridae): cap-independent translation of a viral RNA controlled in cis by an (a)virulence determinant. Plant J. 2008;56(5):716–27.View ArticlePubMedGoogle Scholar
- Nieto C, Morales M, Orjeda G, Clepet C, Monfort A, Sturbois B, et al. An eIF4E allele confers resistance to an uncapped and non-polyadenylated RNA virus in melon. Plant J. 2006;48(3):452–62.View ArticlePubMedGoogle Scholar
- Nieto C, Rodriguez-Moreno L, Rodriguez-Hernandez AM, Aranda MA, Truniger V. Nicotiana benthamiana resistance to non-adapted Melon necrotic spot virus results from an incompatible interaction between virus RNA and translation initiation factor 4E. Plant J. 2011;66(3):492–501.View ArticlePubMedGoogle Scholar
- Díaz JA, Nieto C, Moriones E, Truniger V, Aranda MA. Molecular characterization of a Melon necrotic spot virus strain that overcomes the resistance in melon and nonhost plants. Mol Plant Microbe Interact. 2004;17(6):668–75.View ArticlePubMedGoogle Scholar
- Conesa A, Nueda MJ, Ferrer A, Talon M. maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments. Bioinformatics. 2006;22(9):1096–102.View ArticlePubMedGoogle Scholar
- Conesa A, Gotz S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.View ArticlePubMedGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95(25):14863–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response (vol 98, pg 5116, 2001). Proc Natl Acad Sci U S A. 2001;98(18):10515.View ArticleGoogle Scholar
- Whitham SA, Quan S, Chang HS, Cooper B, Estes B, Zhu T, et al. Diverse RNA viruses elicit the expression of common sets of genes in susceptible Arabidopsis thaliana plants. Plant J. 2003;33(2):271–83.View ArticlePubMedGoogle Scholar
- Hillung J, Cuevas JM, Elena SF. Transcript Profiling of Different Arabidopsis thaliana Ecotypes in Response to Tobacco etch potyvirus Infection. Front Microbiol. 2012;3:229.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang CL, Guo R, Jie F, Nettleton D, Peng JQ, Carr T, et al. Spatial analysis of Arabidopsis thaliana gene expression in response to Turnip mosaic virus infection. Mol Plant Microbe Interact. 2007;20(4):358–70.View ArticlePubMedGoogle Scholar
- Agudelo-Romero P, Carbonell P, de la Iglesia F, Carrera J, Rodrigo G, Jaramillo A, et al. Changes in the gene expression profile of Arabidopsis thaliana after infection with Tobacco etch virus. Virol J. 2008;5:92.View ArticlePubMedPubMed CentralGoogle Scholar
- Hanssen IM, van Esse HP, Ballester AR, Hogewoning SW, Parra NO, Paeleman A, et al. Differential tomato transcriptomic responses induced by pepino mosaic virus isolates with differential aggressiveness. Plant Physiol. 2011;156(1):301–18.View ArticlePubMedPubMed CentralGoogle Scholar
- Wise RP, Moscou MJ, Bogdanove AJ, Whitham SA. Transcript profiling in host-pathogen interactions. Annu Rev Phytopathol. 2007;45:329–69.View ArticlePubMedGoogle Scholar
- Miller WA, White KA. Long-distance RNA-RNA interactions in plant virus gene expression and replication. Annu Rev Phytopathol. 2006;44:447–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Díaz-Vivancos P, Clemente-Moreno MJ, Rubio M, Olmos E, Garcia JA, Martinez-Gomez P, et al. Alteration in the chloroplastic metabolism leads to ROS accumulation in pea plants in response to plum pox virus. J Exp Bot. 2008;59(8):2147–60.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(6):941–54.View ArticlePubMedGoogle Scholar
- Abbink TE, Peart JR, Mos TN, Baulcombe DC, Bol JF, Linthorst HJ. Silencing of a gene encoding a protein component of the oxygen-evolving complex of photosystem II enhances virus replication in plants. Virology. 2002;295(2):307–19.View ArticlePubMedGoogle Scholar
- Shabala S, Babourina O, Rengel Z, Nemchinov LG. Non-invasive microelectrode potassium flux measurements as a potential tool for early recognition of virus-host compatibility in plants. Planta. 2010;232(4):807–15.View ArticlePubMedGoogle Scholar
- Coudriet DL, Kishaba AN, Bohn GW. Inheritance of resistance to muskmelon. J Am Soc Hortic Sci. 1981;106:789–91.Google Scholar
- Alazem M, Lin NS. Roles of plant hormones in the regulation of host-virus interactions. Mol Plant Pathol. 2015;16:529–540.Google Scholar
- Takabatake R, Ando Y, Seo S, Katou S, Tsuda S, Ohashi Y, et al. MAP kinases function downstream of HSP90 and upstream of mitochondria in TMV resistance gene N-mediated hypersensitive cell death. Plant Cell Physiol. 2007;48(3):498–510.View ArticlePubMedGoogle Scholar
- Tena G, Boudsocq M, Sheen J. Protein kinase signaling networks in plant innate immunity. Curr Opin Plant Biol. 2011;14(5):519–29.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang KY, Liu Y, Zhang S. Activation of a mitogen-activated protein kinase pathway is involved in disease resistance in tobacco. Proc Natl Acad Sci U S A. 2001;98(2):741–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang C, Zhang C, Dittman JD, Whitham SA. Differential requirement of ribosomal protein S6 by plant RNA viruses with different translation initiation strategies. Virology. 2009;390(2):163–73.View ArticlePubMedGoogle Scholar
- Aranda MA, Escaler M, Thomas CL, Maule AJ. A heat shock transcription factor in pea is differentially controlled by heat and virus replication. Plant J. 1999;20(2):153–61.View ArticlePubMedGoogle Scholar
- Igari K, Endo S, Hibara K, Aida M, Sakakibara H, Kawasaki T, et al. Constitutive activation of a CC-NB-LRR protein alters morphogenesis through the cytokinin pathway in Arabidopsis. Plant J. 2008;55(1):14–27.View ArticlePubMedGoogle Scholar
- Bari R, Jones JD. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88.View ArticlePubMedGoogle Scholar
- Postnikova OA, Nemchinov LG. Comparative analysis of microarray data in Arabidopsis transcriptome during compatible interactions with plant viruses. Virol J. 2012;9:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Rodrigo G, Carrera J, Ruiz-Ferrer V, del Toro FJ, Llave C, Voinnet O, et al. A meta-analysis reveals the commonalities and differences in Arabidopsis thaliana response to different viral pathogens. Plos One. 2012;7(7):e40526.View ArticlePubMedPubMed CentralGoogle Scholar
- Rizzo TM, Palukaitis P. Nucleotide sequence and evolutionary relationships of cucumber mosaic virus (CMV) strains: CMV RNA 1. J Gen Virol. 1989;70(Pt 1):1–11.View ArticlePubMedGoogle Scholar
- Díaz-Pendón JA, Fernández-Muñoz R, Gómez-Guillamon ML, Moriones E. Inheritance of Resistance to Watermelon mosaic virus in Cucumis melo that Impairs Virus Accumulation, Symptom Expression, and Aphid Transmission. Phytopathology. 2005;95(7):840–6.View ArticlePubMedGoogle Scholar
- Saladié M, Cañizares J, Phillips MA, Rodríguez-Concepción M, Larrigaudiere C, Gibon Y, et al. Comparative transcriptional profiling analysis of developing melon (Cucumis melo L.) fruit from climacteric and non-climacteric varieties. BMC Genomics. 2015;16:440.View ArticlePubMedPubMed CentralGoogle Scholar
- Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26(19):2363–7. doi:https://doi.org/10.1093/bioinformatics/btq431.View ArticlePubMedPubMed CentralGoogle Scholar
- Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, et al. TM4 microarray software suite. Methods Enzymol. 2006;411:134–93.View ArticlePubMedGoogle Scholar
- Graur D, Li WH. Fundamentals of Molecular Evolution. 2nd ed. Sunderland: Sinauer Associates; 2000. p. 209–10.Google Scholar
- Soukas A, Cohen P, Socci ND, Friedman JM. Leptin-specific patterns of gene expression in white adipose tissue. Genes Dev. 2000;14(8):963–80.PubMedPubMed CentralGoogle Scholar