- Research article
- Open Access
Root transcriptional responses of two melon genotypes with contrasting resistance to Monosporascus cannonballus (Pollack et Uecker) infection
BMC Genomics volume 13, Article number: 601 (2012)
Monosporascus cannonballus is the main causal agent of melon vine decline disease. Several studies have been carried out mainly focused on the study of the penetration of this pathogen into melon roots, the evaluation of symptoms severity on infected roots, and screening assays for breeding programs. However, a detailed molecular view on the early interaction between M. cannonballus and melon roots in either susceptible or resistant genotypes is lacking. In the present study, we used a melon oligo-based microarray to investigate the gene expression responses of two melon genotypes, Cucumis melo ‘Piel de sapo’ (‘PS’) and C. melo ‘Pat 81’, with contrasting resistance to the disease. This study was carried out at 1 and 3 days after infection (DPI) by M. cannonballus.
Our results indicate a dissimilar behavior of the susceptible vs. the resistant genotypes from 1 to 3 DPI. ‘PS’ responded with a more rapid infection response than ‘Pat 81’ at 1 DPI. At 3 DPI the total number of differentially expressed genes identified in ‘PS’ declined from 451 to 359, while the total number of differentially expressed transcripts in ‘Pat 81’ increased from 187 to 849. Several deregulated transcripts coded for components of Ca2+ and jasmonic acid (JA) signalling pathways, as well as for other proteins related to defence mechanisms. Transcriptional differences in the activation of the JA-mediated response in ‘Pat 81’ compared to ‘PS’ suggested that JA response might be partially responsible for their observed differences in resistance.
As a result of this study we have identified for the first time a set of candidate genes involved in the root response to the infection of the pathogen causing melon vine decline. This information is useful for understanding the disease progression and resistance mechanisms few days after inoculation.
Monosporascus cannonballus (Pollack et Uecker) is an ascomycete soil-borne pathogen adapted to hot, arid and semi-arid areas that causes root/rot vine decline. This disease has severe economic impact on global melon (Cucumis melo L) and watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai) production . Infection of roots may occur via germinating ascospores and/or active mycelium in infested soils. When roots are infected, they become necrotic with numerous discrete lesions causing the loss of most of secondary and tertiary roots. During the growing season, symptoms are characterized by reduced plant growth, progressive defoliation and partial or complete canopy collapse, resulting in fruit sunburn and total crop loss at harvest.
Ascospores probably function as the primary survival structure, as well as the primary inoculum for root infection [2, 3]. Studies on the specificity of M. cannonballus showed that germination of ascospores is extremely host specific and occurs only in the rhizosphere of certain genera and species of plants belonging exclusively to the Cucurbitaceae family . In a histological study using artificial inoculation with active mycelium, Alfaro-Fernandez and García-Luis  examined the early colonization of M. cannonballus in two cucurbit species that differed in their sensitivity to this disease: a highly sensitive muskmelon (C. melo) and a tolerant squash (Cucurbita maxima). Results showed that M. cannonballus was capable of infecting the tissue of both host plants, colonizing the epidermis and cortex with decreased density of mycelium at the endodermis level. Differences in sensitivity of muskmelon and squash seemed to be due to the differential resistance to the initial penetration of the fungus.
Studies show that M. cannonballus is a facultative saprophyte fungus, does not form special penetration structures such as appressoria and pathogen reproduction in infected roots occurs primarily after plant death . It behaves in a manner most similar to that of vascular wilt pathogens, but differs in that it is not systemic and cannot be isolated from aerial portions of infected plants .
Strategies for controlling vine decline based on the development of resistant cultivars have been initiated in most affected countries. Since large genetic variability is displayed within C. melo varieties, the selection of resistant material has become a major objective for plant breeding approaches . C. melo germplasm is commonly classified into two subspecies, ssp. melo including the main commercial types, and ssp. agrestis with important sources of resistance and quality traits . The genotype ‘Pat 81’ of C. melo ssp. agrestis showed resistance to M. cannonballus under field conditions [10, 11], which was employed to initiate a breeding program aimed at introducing resistance into Spanish melon cultivars such as C. melo ssp. melo ‘Piel de sapo’ (‘PS’) [12–14]. The resistance to vine decline of ‘Pat 81’ was expressed as a delay in the appearance of root lesions and as a slow rate of disease development with a low percentage of wilted plants at the end of the growing cycle. Root lesions caused by M. cannonballus infection in ‘Pat 81’ were less severe than in susceptible genotypes, being limited to small damaged areas without loss of root biomass. A vigorous and branched root structure with high regeneration potential could also help to maintain a good hydraulic conductivity and improve resistance to the infection towards the end of the season .
During the last decade, relevant genetic and genomic tools have been developed in melon. Available resources include new mapping populations [16, 17], genetic maps [18–21], BAC-based physical maps , melon transcriptome characterization , TILLING and ECOTILLING platforms [24–26] and large EST collections [27, 28]. Recently, an oligonucleotide-based microarray has been developed from a dataset of 17,510 unigenes obtained from normalized cDNA libraries, representing different melon genotypes, tissues and physiological conditions, including M. cannonballus infected root tissue . A preliminary comparison between ‘Pat 81’ infected vs. non-infected roots was used for the validation of this microarray platform . Inoculation was carried out on artificially infected soil and samples were collected 14 days after inoculation. At this late stage of the infection, only a low number of genes were found to be differentially regulated between M. cannonballus infected and non-infected roots of ‘Pat 81’.
To date, a detailed molecular view on the early interaction between M. cannonballus and melon roots in either susceptible or resistant genotypes is lacking. To improve our knowledge on the genetic regulation of root/rot vine decline resistance in melon, the transcriptional changes associated with M. cannonballus inoculation of the susceptible ‘PS’ vs. the resistant ‘Pat 81’ are studied here at two time points of early infection.
Results and discussion
Root infection with M. cannonballus occurs rapidly in ‘PS’
The study of the early response of melon roots to fungal infection by M. cannonballus required a method for rapid infection of roots different from the traditional soil inoculation procedures employed previously . A method based on the direct contact of melon roots with M. cannonballus mycelium, and subsequent growth in a hydroponic system optimally met these requirements.
Infected roots of ‘PS’ and their respective mock-inoculated controls were collected at 1, 3 and 5 days post inoculation (DPI). There were no visible symptoms of disease at these time points; therefore, quantitative PCR was used to detect the pathogen in these roots . The amount of M. cannonballus DNA detected in root samples increased significantly from 1 to 3 DPI (2 to 500 pg of M. cannonballus DNA per ng of total DNA), and subsequently remained constant until 5 DPI (Figure 1). Consequently, samples from 1 and 3 DPI were chosen for transcriptional profiling of the early stages of M. cannonballus infection in the roots of the susceptible genotype ‘PS’ and the resistant genotype ‘Pat 81’.
Global gene expression trends reflect susceptibility to M. cannonballus root infection
Changes in the expression of genes during M. cannonballus root infection were identified using a melon oligonucleotide microarray (Roche NimbleGen, Madison, WI, USA). This microarray was designed using 17,443 unigenes assembled from different melon cDNA libraries, including two libraries obtained from roots infected with M. cannonballus[27, 29]. Inoculated and mock-inoculated roots of the susceptible ‘PS’ and the resistant ‘Pat 81’ were collected at 1 DPI and 3 DPI. Three biological replicates of cDNA obtained for each condition were labelled and used for microarray hybridization.
Principal component analysis (PCA) was used for preliminary data analysis and quality control. The first and second components explained 29.7% and 20.7% of total variance, respectively. Biological replicates from each accession and stage (inoculated and mock-inoculated) clustered together in all cases (Figure 2).
Both root developmental stage and fungal infection seem to be associated with the first component (Figure 2). The differences observed between roots collected 1 day and 3 days in control plants, similar in both genotypes, suggest that there exist a large number of genes deregulated during root development. We also observed striking differences between inoculated and mock-inoculated roots (both at 1 and 3 DPI), with differential patterns in ‘PS’ and ‘Pat 81’. In fact, at 3 DPI, the infection produced similar plot shifts in both genotypes. However, the global gene expression differences between 1 DPI-infected and mock sample in ‘Pat 81’ background were lower when compared to the equivalent ’PS’ samples (Figure 2). This observation is consistent with slower M. cannonballus infection of ‘Pat 81’. The second component explained genotypic differences between ‘PS’and ‘Pat 81’ samples, suggesting the influence of genotype on hybridization signals.
At 1 DPI, the abundance of 165 transcripts significantly increased and 286 decreased in ‘PS’, whereas the abundance of 56 transcripts increased and 131 decreased in ‘Pat 81’ (Additional file 1 and Additional file 2). Thus, most of modulated transcripts at 1 DPI reduced their abundance following M. cannonballus infection (Figure 3). A similar enrichment in down-regulated genes was observed during the susceptible plant-pathogen interaction of cotton roots with soil-borne, vascular wilt fungus Fusarium oxysporum f. sp. vasinfectum, particularly at early stages of infection . However, at 3 DPI, the ratio of induced/repressed transcripts was close to one and similar in both genotypes (Figure 3; Additional file 3 and Additional file 4). The total number of differentially expressed genes identified in ‘PS’ declined from 451 (1 DPI) to 359 (3 DPI), while the total number of differentially expressed transcripts in ‘Pat 81’ increased from 187 (1 DPI) to 849 (3 DPI), which emphasizes the dissimilar behaviour of the susceptible and the resistant genotypes from the initial stages of infection.
The overlap of differentially expressed transcripts shared by two or more pairs of genotype-time groups, was small for most combinations (Figure 4). The percentage of total (induced and repressed) unique transcripts in ‘Pat 81’-1DPI, ‘Pat 81’-3DPI, ‘PS’-1DPI and ‘PS’-3DPI were 54%, 63%, 56% and 27% respectively. Transcriptional profiles for ‘PS’ at 3 DPI and ‘Pat 81’ at 3 DPI had a high number of common induced (113) and repressed transcripts (65), suggesting that at 3 DPI both susceptible and resistant genotypes, had initiated a set of common transcriptional responses to M. cannonballus infection. In contrast, the number of common transcripts induced (13) or repressed (31) at 1DPI between both genotypes was smaller (Figure 4), consistently with the differential response of these two genotypes at the very early stages of infection.
When comparing these data with previous results studying the transcriptomic response of melon to M. cannonballus we found only 11 coincident genes, 4 of them repressed and 7 induced (Additional file 5). None of these genes was directly related to known plant defence mechanisms. Such divergence between experiments could be due to some differences in infection time (14 days in  versus 1–3 days in this work) and the inoculation procedure. Whereas in the previous work the plants were grown on soil infected with 50 colony-forming units of M. cannonballus per gram, in this work roots were inoculated by direct contact with M. cannonballus mycelium and plants were grown on liquid medium, allowing a faster and more uniform infection.
Quantitative reverse transcription-PCR validation of microarray data
In order to confirm microarray results, nine differentially expressed transcripts with different expression profiles were selected and tested by quantitative reverse transcription-PCR (qRT-PCR) analysis in RNA samples from the four assayed genotype/ time combinations (primers in Additional file 6). Selected genes were mainly involved in defence (induced and repressed in different conditions), signal transduction, metabolism, and stress response. All transcript profiles corresponding to differentially expressed genes in each condition showed similar magnitudes and directions of change, with a correlation coefficient of 0.83 (Figure 5). These results confirm the reliability of our microarray data.
Differentially expressed transcripts were assigned to functional categories using MELOGEN database information , published data, and functional annotation performed with the BLAST2go package . The percentage of transcripts assigned to each functional category in each treatment is summarized in Figure 6. “Signal transduction” category represented 9% of total genes in ‘PS’ and 10% in ‘Pat 81’ at 1 DPI, with only one common gene. At 3 DPI, the ratio of differentially expressed genes in this category slightly increased to about 11%, but the number of common induced genes to both genotypes was higher (17 transcripts). In the “stress response” category 8% and 14% of differentially expressed genes were observed in ‘PS’ and ‘Pat 81’ 1 DPI, respectively, with the majority down-regulated. At 3 DPI the percentage of genes categorized under ‘stress response’ was 8% in both genotypes, however the number of induced genes was twice the number of repressed genes in ‘Pat 81’. The number of genes classified in the “defence” group ranged from 2.7% to 3.8% and most of them were down-regulated. Transcripts categorized under “hormone metabolism” were mainly down-regulated in both genotypes and sampling dates, although some of them were up-regulated at 3 DPI in both genotypes. Genes belonging to these four categories are discussed in more detail below.
Transcriptional evidence for Ca2+ and jasmonic acid in signalling pathogen responses
Calcium (Ca2+) plays a crucial role as a messenger in transduction during various developmental processes and in response to stress, including biotic stresses caused by fungi, bacteria and viruses [34, 35]. In this study we identified seven transcripts up-regulated in response to M. cannonballus infection with similarity to Ca2+ sensor genes (Table 1). Three of them, coding for calmodulin (CaM)-like proteins (cCL1433Contig1, cCL417Contig1 and c46d_38-D02-M13R), had great transcript abundance exclusively in ‘Pat 81’ 1 DPI. At 3 DPI, two transcripts with homology to calcineurin B-like (CBL) protein (cCL3052Contig1) and CBL-interacting serine/threonine protein kinases (cCL3317 contig1) also had greater transcript abundance specifically in ‘Pat 81’ (Table 1), whereas the CBL-interacting serine/threonine protein kinase coded by cHS_34-D08-M13R was found to be up-regulated in both ‘Pat 81’ and ‘PS’. An Additional transcript coding for a calcium-dependent protein kinase (CDPK, cCL_57-B01-M13R) was found up-regulated in ‘PS’3 DPI. Recent evidence suggests that Ca2+ signalling via CDPKs, CBL/CBL-interacting protein kinases and CaM are involved in different aspects of biotic defence response [36, 37]. Up-regulation of these genes mainly in the ‘Pat 81’ genotype suggests that Ca2+ signalling may play a part in the melon-M. cannonballus interaction.
Members of the complex family of WRKY transcriptional factors play a broad and pivotal role in regulating defence . Most of the WRKY-like transcripts identified in this work were found to have lower transcript abundance in infected samples. The transcript cPSI_30-E06-M13R, which was down-regulated in ‘Pat 81’ 3 DPI, showed high homology to ClWRKY70 from Citrullus lanatus (Table 1), a salicylic acid (SA)-inducible gene that increases resistance to pathogens by overexpression in Arabidopsis. WRKY70 from Arabidopsis represses jasmonic acid (JA)-responsive genes and is itself down-regulated by JA, acting at a convergence point determining the balance between SA- and JA-dependent defence pathways. WRKY70 is also required for R gene-mediated resistance [40–42]. Down-regulation of this transcript in ‘Pat 81’ points to the prevalence of JA signalling cascade over SA pathway for activating the early defence mechanisms triggered by M. cannonballus infection. In general, SA-dependent gene expression responses are effective against biotrophic pathogens, whereas defence genes activated by JA are effective against necrotrophic pathogens. Some exceptions to this statement have been documented, and frequently pathogens cannot be merely classified as biotrophic or necrotrophic as in the case of M. cannonballus[43, 44].
Additional transcriptional changes indicate a role of JA pathways in the response to M. cannonballus. For example, among genes specifically down-regulated in ‘Pat 81’ 1 and 3 DPI, transcripts cCL935Contig1 and cCL935Contig2 had homology to JASMONATE-ZIM DOMAIN 10 (JAZ10) of Arabidopsis (Table 1). JAZ proteins are negative regulators of the JA signal cascade through the interaction with certain transcription factors such as MYC2. Fast down-regulation of these genes upon M. cannonballus infection suggests the possibility of proteasome-independent pathways for the activation of JA-mediated plant defence responses in ‘Pat 81’ genotype .
The transcripts cPS_05-A06-M13R and cA_30-B04-M13R, encoding two related cytochrome P450 proteins, were induced 2 to 4-fold in ‘PS’ and 3 to 9-fold in ‘Pat 81’ 3 DPI. These transcripts show a high similarity to CYP82C2 from Arabidopsis, involved in JA-induced expression of defence genes and accumulation of indole glucosinolates. A mutation in CYP82C2 gene reduces plant resistance to the necrotrophic fungus Botrytis cinerea and alters root growth sensitivity to exogenous JA, whereas CYP82C2 overexpression improves resistance to B. cinerea. These data suggest the preferential activation of the JA signalling pathway 3 days after M. cannonballus infection of melon roots, instead of SA-dependent cascades. Interestingly, the differential quantitative and qualitative expression of CYP82C2-like and JAZ10-like genes in ‘Pat 81’ and ‘PS’ genotypes suggest that JA response might be partially responsible for their observed differences in resistance, although sequence differences between the two genotypes affecting array hybridization may also account for part of this variation.
Different studies indicate that JA- and ethylene-signalling frequently operate synergistically to induce the effector genes of defence responses , however we found no transcriptional evidence under our experimental conditions: several transcripts (cA_02-A09-M13R, cFR15J17, cCL451Contig1) coding for ACC synthase (ACS) and ACC oxidase (ACO), the enzymes catalysing the last two steps in the ethylene biosynthetic pathway, were found repressed in ‘Pat 81’ and Piel de sapo’ at different infection times (Table 1). Other authors have found also different transcripts with similarity to ACC oxidase genes differentially expressed after infection of melon with Fusarium oxysporum f. sp. melonis. However, this does not preclude post-transcriptional and translational processes altering the activity of these enzymes and the production of ethylene. Further work is required to elucidate the roles of these hormones during the melon-M. cannonballus interaction.
Defence responses show transcriptional similarity to previously identified pathogen responses
In ‘Pat 81’ a transcript (cCL3498Contig1) homologous to AtSYP121/PEN1, encoding the protein syntaxin 121 from Arabidopsis, was induced at early stages of infection (1 DPI), but did not change its expression significantly in the susceptible ‘PS’ (Table 1). PEN1 was identified in a screening for penetration (pen) mutants, required for the resistance to fungal penetration in the non host interaction between Arabidopsis and Blumeria graminis f. sp. hordei (Bgh) . PEN1 protein is a constituent of a SNARE complex that contributes to the formation of cell wall appositions [49, 50]. The ortholog of PEN1 in barley (ROR2) was described as required for basal penetration resistance against Bgh in mlo (mildew resistance locus o) mutants [49, 51]. In addition to barley, loss-of-function mutations of MLO genes conferred broad-spectrum resistance to powdery mildew in Arabidopsis, tomato and pea (Pisum sativum) [52–56]. In our study, 2 transcripts (cPS_18-B02-M13R and cFR17N13) with similarity to MLO genes were found to be down-regulated in ‘Pat 81’ at 3 DPI. A Blastx comparison of these ESTs against the Arabidopsis, tomato and pea protein databases using an E-value cut-off of 10-5 found only components of the MLO gene family from these three species. The first Arabidopsis blastx hits of cPS_18-B02-M13R and cFR17N13 were respectively MLO4 and MLO6. MLO4 was recently described to affect growth responses of the Arabidopsis root in response to mechanical stimuli . While the disruption of MLO6, together with mutants in the related MLO2 and MLO12 genes, was required for the resistance to powdery mildew .
To date, there are no evidences supporting the occurrence of cell wall appositions in melon roots infected with M. cannonballus; however it is tempting to speculate that an early expression of a PEN1-like could contribute to delay and prevent to some extent the penetration of the fungus in ‘Pat 81’, resulting in an altered development of the infection. The subsequent repression of MLO-like genes, detected two days later, could tune up this specific response in ‘Pat 81’. Recently, a MLO-like gene has been described in melon . The expression of this gene, designated CmMlo1, was up-regulated under cadmium exposure, which suggested its participation in abiotic stress responses, but this transcript is different from the transcripts identified in this work.
A gene specifically expressed in ‘Pat 81’ (cPS_23-C06-M13R) had similarity to GhZFP1, encoding a CCCH-type zinc-finger transcription factor of Gossypium hirsutum. Recently, this protein was characterized as a relevant positive regulator conferring salt tolerance and fungal pathogen resistance to plants . Overexpression of GhZFP1 in transgenic tobacco enhanced tolerance to salt stress and resistance to Rhizoctonia solani. Two possible interactors of GhZFP1 protein: GZIRD21A, similar to responsive to dehydration protein 21A, and GZIPR5, a pathogenesis-related protein 5 (PR5)-like were also identified . Interestingly, a transcript with high similarity to PR5 of Arabidopsis was found up-regulated in ‘Pat 81’.
Several significantly differentially expressed transcripts with high similarity to pathogenesis-related (PR) genes (glucanases and chitinases among others) were down-regulated in both genotypes except for a transcript with homology to PR5 of Arabidopsis (cCL4557Contig1), coding for a thaumatin-like protein, which was up-regulated in ‘Pat 81’ at 3 DPI. Such reduction of PR-related transcripts suggests that PR-specific defence mechanisms are not activated within the first 3 days after M. cannonballus inoculation. The down-regulation of PR genes was also reported by Schlink  in infected roots of Fagus sylvatica by Phytophthora citricola (hemibiotrophic oomycete), where they hypothesized that down-regulation would alter the pathogens’ chance to escape recognition. Nevertheless, in our system, additional studies on the late pathogenic response are required to elucidate the role of specific PR proteins in melon-M. cannonballus interaction.
Non-pathogen related stress responses are consistent between genotypes
Several transcripts related to non-pathogen stress responses were differentially expressed after fungal infection. These included members of the small and large heat-shock protein families (HSP), which were up-regulated in both genotypes (cSSH1P11, cPSI_23-F06-M13R, cPSI_32-H04-M13R, cCL172Contig1, cCL3362Contig1, cCL5902Contig1). We also found transcripts coding for proteins related to the oxidative stress and the regulation of reactive oxygen species (ROS) as peroxidase-like and glutathione S-transferase (cCL3733Contig1, cPSI_21-C11-M13R, cPS_30-C04-M13R, cAI_08-H07-M13R) to be significantly down-regulated in both genotypes.
The results show common and divergent responses of the susceptible and resistant melon genotypes to infection with M. cannonballus. Transcriptomic differences are more apparent at an early stage of infection. Transcriptional differences in the activation of the JA-mediated response in ‘Pat 81’ compared to ‘PS’ suggest that JA response may be partially responsible for their observed differences in resistance. Several transcripts, previously implicated in plant fungal resistance, were also significantly differentially expressed in ‘Pat 81’, also potentially resulting in an altered infection development. Further studies are needed to quantify differences in tissue hormone concentrations between the two genotypes, as implicated in the differential expression of JA regulated genes, and identify the functional roles of many of the transcripts observed to be expressed more abundantly in ‘Pat 81’ melon compared to the susceptible ‘PS’ genotype. Recently the genome sequence of melon has been reported . The authors predicted 27,427 protein-coding genes. Thus, this work offers a partial view on the whole picture of the transcriptomic changes occurring in our experimental model. Nevertheless these data along with future functional studies could lead to the identification/characterization of defence genes involved in resistance of melon to M. cannonballus vine decline disease.
Two melon accessions, Cucumis melo spp. melo ‘Piel de sapo cv piñonet’ (‘PS’), fully susceptible to the infection by M. cannonballus and C. melo spp. agrestis ‘Pat 81’, resistant to the infection by M. cannonballus were used in this study. These accessions are maintained by Cucurbits Breeding group at COMAV-UPV.
In vitro inoculation of C. melo roots with M. cannonballus
In a preliminary study, seeds from ‘PS’ were surface-sterilized with 20% bleach and a drop of Tween-20 for a minute and after rinsing, the seeds were placed in Petri dishes with wet filter paper under sterile conditions. After 6 days, seedlings were inoculated by direct contact with M. cannonballus mycelium grown on PDA (potato dextrose agar). To ensure the correct inoculation each root was rolled in germination paper with 2 discs (1 cm2 aprox.) of PDA with active growing mycelium. Mock treatments were prepared in the same way using PDA without mycelium. Plants were placed in a container containing a half diluted nutrient solution composed of: 3 mM KNO3, 2 mM Ca(NO3)2.4H2O, 0.5 mM MgSO4.7H2O, 0.5 mM (NH4)H2PO4, 25 μM KCl, 12.5 μM H3BO3, 1 μM MnSO4.H2O, 1μM ZnSO4.7H2O, 0.25 μM CuSO4.5H2O, 1.3 μM (NH4)6 Mo7O24.4H2O and 25 μM Fe-NaEDTA. Axenic conditions were maintained. Plants were grown in a climatic chamber (28°C, 16/8 h light/dark). Infected roots of ‘PS’ melon (four plants per time-point) and their respective mock-inoculated controls (one plant per time-point) were collected at 1, 3 and 5 days post inoculation (DPI). Total DNA was extracted using the DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) and the presence of the fungus and the infection levels were assessed by quantitative PCR with specific M. cannonballus primers as described previously .
Sample collection and RNA isolation
Root samples for the extraction of total RNA used to hybridize to the melon microarray were taken from plants grown using the inoculation method described above. For each biological sample we collected roots from 4 to 6 plants per genotype (‘PS’ and ‘Pat 81’), treatment (inoculated and mock-inoculated as control) and time post inoculation (1DPI and 3DPI). Three biological replicates were used for each genotype/treatment/time combination and independently hybridized to the melon microarray. Total RNA was isolated using TRI Reagent (Sigma-Aldrich Corporation, St. Louis, MO, USA) according to manufacturer’s protocols and further purified using RNeasy Mini Kit (Qiagen). RNA integrity and quality was checked on agarose electrophoresis. Quantity and purity of total RNA were determined by Nano-Drop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Total RNA samples were sent to Roche NimbleGen Systems where cDNA synthesis, Cy3 labelling and hybridizations were performed following the manufacturer’s procedures.
Microarray data analysis
The melon microarray is an oligo-based (60-mer) microarray representing a total of 17,443 unigenes derived from 33,418 high-quality melon ESTs . Sequences of these unigenes are listed in Additional file 7. Hybridization signal intensity was calculated using a GenePix 4000B (Molecular Devices, Sunnyvale, CA, USA) and the data were extracted using NimbleScan software (Roche NimbleGen). The intensity values obtained from the array scanning were normalized using the robust multiarray average (RMA) . Normalized probe set data, provided by Roche NimbleGen Systems in RMA calls files, were imported into ArrayStar software 3.0 version (Dnastar, Madison, WI, USA), where statistical analysis was performed. Data from infected samples were normalized to their respective controls. Data were log2 transformed, thus normalized values are the log2 of the ratio between infected and control samples at a given time-point. Significantly differentially expressed genes were identified using an unpaired t-test with a Benjamini-Hochberg multiple testing corrected p-value cut-off of 0.05  and a fold change cut-off of 2. The microarray data were deposited at ArrayExpress (http://www.ebi.ac.uk/microarray-as/ae/) under the accession number E-MEXP-3732. Transcripts differentially expressed were annotated based on the MELOGEN database , and genes discussed in detail were re-annotated using Cucurbit Genomic Database Melon Unigene v. 4.0 . Additionally, we performed a functional classification of transcripts following the Gene Ontology (GO) scheme with Blast2GO package . This information and previously published data allowed us to classify manually the genes in functional groups.
Principal component analysis (PCA) of all samples was generated using TMeV 4.0 software from TIGR . The Venn diagrams were made manually using the output lists of the statistical analysis.
To validate the microarray experiments, the transcript levels of nine selected genes were quantified using qRT-PCR. First strand cDNA was synthesized from 1 μg of total RNA with the Oligo (dT)20 (50 μM) primer using the Expand Reverse Transcriptase (Roche Applied Science, Penzberg, Germany), according to the manufacturer’s instructions. Quantitative PCR was performed with an ABI PRISM 7000 Sequence Detector System (Applied Biosystems, Foster City, CA, U.S.A), using FastStart Universal SYBR Green Master (ROX) (Roche Applied Science) and 2 μl of diluted 1:10 cDNA for each PCR reaction. The relative expression level was determined using the cyclophilin (cCL1375) housekeeping gene from melon as reference . The gene specific primers for PCR amplification were designed using Primer3 v.0.4.0  (Additional file 6). The fold changes in each infected sample compared to the expression level detected in the corresponding sample under control conditions were calculated using the 2-ΔΔCT method . Intra-assay variation was evaluated by performing all amplification reactions in duplicate.
Martyn RD, Miller ME: Monosporascus root rot and vine decline, an emerging disease of melon worldwide. Plant Dis. 1996, 80: 716-725. 10.1094/PD-80-0716.
Stanghellini ME, Kim DH, Rasmussen SL: Ascospores of Monosporascus cannonballus. Germination and distribution in cultivated and desert soils in Arizona. Phytopathology. 1996, 86: 509-514. 10.1094/Phyto-86-509.
Stanghellini ME, Kim DH, Waugh M: Microbe-mediated germination of ascospores of Monosporascus cannonballus. Phytopathology. 2000, 90: 243-247. 10.1094/PHYTO.2000.90.3.243.
Stanghellini ME, Alcantara TP, Ferrin DM: Germination of Monosporascus cannonballus ascospores in the rhizosphere: a hostspecific response. Can J Plant Pathol. 2010, 32: 402-405. 10.1080/07060661.2010.499270.
Alfaro-Fernandez A, Garcia-Luis A: Colonisation and histological changes in muskmelon and autumn squash tissues infected by Acremonium cucurbitacearum or Monosporascus cannonballus. Eur J Plant Pathol. 2009, 125: 73-85. 10.1007/s10658-009-9460-0.
Waugh MM, Kim DH, Ferrin DM, Stanghellini ME: Reproductive potential of Monosporascus cannonballus. Plant Dis. 2003, 87: 45-50. 10.1094/PDIS.2003.87.1.45.
Waugh MM, Ferrin DM, Stanghellini ME: Colonization of cantaloupe roots by Monosporascus cannonballus. Mycol Res. 2005, 109: 1297-1301. 10.1017/S0953756205003722.
Fita A, Picó B, Nuez F: Melon roots under stress: melon vine decline. Plant Stress. 2007, 1: 93-104.
Pitrat M: Melon (Cucumis melo L.). Handbook of Crop Breeding Vol I. Vegetables. Edited by: Prohens J, Nuez F. 2008, New York: Springer, 283-315.
Esteva J, Nuez F: Field resistance to melon dieback in Cucumis melo L. Cucurbit Genet Cooperative. 1994, 17: 76-77.
Iglesias A, Nuez F: Caracterización de diversas entradas de melón frente al colapso o muerte súbita. Actas Horticultura. 1998, 22: 139-147.
Iglesias A, Picó B, Nuez F: A temporal genetic analysis of disease resistance genes: resistance to melon vine decline derived from Cucumis melo var agrestis. Plant Breed. 2000, 118: 1-6.
Dias R, De CS, Picó B, Espinós A, Nuez F: Resistance to melon vine decline from Cucumis melo subs agrestis: genetic analysis for root structure and rot response. Plant Breed. 2004, 123: 66-72. 10.1046/j.1439-0523.2003.00944.x.
Fita A, Picó B, Dias RCS, Nuez F: ‘Piel de Sapo’ breeding lines tolerant to melon vine decline. HortSci. 2009, 44: 1458-1460.
Fita A, Picó B, Dias R, Nuez F: Effects of root architecture on response to melon vine decline. J Hortic Sci Biotech. 2008, 83: 616-623.
Perin C, Hagen S, De Conto V, Katzir N, Danin-Poleg Y, Portnoy V, Baudracco-Arnas S, Chadoeuf J, Dogimont C, Pitrat M: A reference map of Cucumis melo based on two recombinant inbred line populations. Theor Appl Genet. 2002, 104: 1017-1034. 10.1007/s00122-002-0864-x.
Eduardo I, Arús P, Monforte AJ: Development of a genomic library of near isogenic lines (NILs) in melon (Cucumis melo L.) from the exotic accession PI161375. Theor Appl Genet. 2005, 112: 139-148. 10.1007/s00122-005-0116-y.
Gonzalo MJ, Oliver M, Garcia-Mas J, Monfort A, Dolcet-Sanjuan R, Katzir N, Arus P, Monforte AJ: Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.). Theor Appl Genet. 2005, 110: 802-811. 10.1007/s00122-004-1814-6.
Fernandez-Silva I, Eduardo I, Blanca J, Esteras C, Pico B, Nuez F, Arus P, Garcia-Mas J, Monforte AJ: Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.). Theor Appl Genet. 2008, 118: 139-150. 10.1007/s00122-008-0883-3.
Deleu W, Esteras C, Roig C, González-To M, Fernández-Silva I, González-Ibeas D, Blanca J, Aranda MA, Arús P, Nuez F, Monforte AJ, Picó MB, Garcia-Mas J: A set of EST-SNPs for map saturation and cultivar identification in melon. BMC Plant Biol. 2009, 9: 90-10.1186/1471-2229-9-90.
Harel-Beja R, Tzuri G, Portnoy V, Lotan-Pompan M, Lev S, Cohen S, Dai N, Yeselson L, Meir A, Libhaber SE, Avisar E, Melame T, van Koert P, Verbakel H, Hofstede R, Volpin H, Oliver M, Fougedoire A, Stalh C, Fauve J, Copes B, Fei Z, Giovannoni J, Ori N, Lewinsohn E, Sherman A, Burger J, Tadmor Y, Schaffer AA, Katzir N: A genetic map of melon highly enriched with fruit quality QTLs and EST markers, including sugar and carotenoid metabolism genes. Theor Appl Genet. 2010, 121: 511-533. 10.1007/s00122-010-1327-4.
González V, Garcia-Mas J, Arús P, Puigdomènech P: Generation of a BAC based physical map of the melon genome. BMC Genomics. 2010, 11: 339-10.1186/1471-2164-11-339.
Blanca J, Cañizares J, Ziarsolo P, Esteras C, Mir G, Nuez F, Garcia-Mas J, Picó B: Melon transcriptome characterization. SSRs and SNPs discovery for high throughput genotyping across the species. The Plant Genome. 2011, 4: 118-131. 10.3835/plantgenome2011.01.0003.
Nieto C, Piron F, Dalmais M, Marco CF, Moriones E, Gómez-Guillamón ML, Truniger V, Gómez P, Garcia-Mas J, Aranda MA, Bendahmane A: EcoTILLING for the identification of allelic variants of melon eIF4E, a factor that controls virus susceptibility. BMC Plant Biol. 2007, 7: 34-10.1186/1471-2229-7-34.
Dahmani-Mardas F, Troadec C, Boualem A, Leveque S, Alsadon AA, Aldoss AA, Dogimont C, Bendahman A: Engineering melon plants with improved fruit shelf life using the TILLING approach. PLoS One. 2010, 5: e15776-10.1371/journal.pone.0015776.
González M, Xu M, Esteras C, Roig C, Monforte AJ, Troadec C, Pujol M, Nuez F, Bendahmane A, Garcia-Mas J, Picó B: Towards a TILLING platform for functional genomics in Piel de Sapo melons. BMC Res Notes. 2011, 4: 289-10.1186/1756-0500-4-289.
González-Ibeas D, Blanca J, Roig C, González-To M, Picó B, Truniger V, Gómez P, Deleu W, Caño-Delgado A, Arús P, Nuez F, García-Mas J, Puigdomènech P, Aranda MA: MELOGEN: an EST database for melon functional genomics. BMC Genomics. 2007, 8: 306-10.1186/1471-2164-8-306.
Clepet C, Joobeur T, Zheng Y, Jublot D, Huang M, Truniger V, Boualem A, Hernandez-Gonzalez ME, Dolcet-Sanjuan R, Portnoy V, Mascarell-Creus A, Cano-Delgado AI, Katzir N, Bendahmane A, Giovannoni JJ, Aranda MA, Garcia-Mas J, Fei Z: Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. BMC Genomics. 2011, 12: 252-10.1186/1471-2164-12-252.
Mascarell-Creus A, Cañizares J, Vilarrasa J, Mora-García S, Blanca J, González-Ibeas D, Saladié M, Roig C, Deleu W, Picó B, López-Bigas N, Aranda MA, Garcia-Mas J, Nuez F, Puigdomènech P, Caño-Delgado A: 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-10.1186/1471-2164-10-467.
Roig C, Picó B, Fita A, Nuez F: Phenotyping root response to complex diseases: the Monosporascus root rot/vine decline of melons. Proceedings of IXth EUCARPIA Meeting on Genetics and Breeding of Cucurbitaceae: 21-24 May 2008. Edited by: Pitrat M. 2008, Avignon: INRA, 423-428.
Picó B, Roig C, Fita A, Nuez F: Quantitative detection of Monosporascus cannonballus in infected melon roots. Eur J Plant Pathol. 2008, 120: 147-156. 10.1007/s10658-007-9203-z.
Dowd C, Wilson LW, McFadden H: Gene expression profile changes in cotton root and hypocotyl tissues in response to infection with Fusarium oxysporum f. sp vasinfectum. Mol Plant Microbe In. 2004, 17: 654-667. 10.1094/MPMI.2004.17.6.654.
Conesa A, Götz S: Blast2GO: a suite for functional analysis in plant genomics. Int J Plant Genomics. 2008, 200: 619832-
Reddy VS, Ali GS, Reddy AS: Characterization of a pathogen-induced calmodulin-binding protein: mapping of four Ca2+−dependent calmodulin-binding domains. Plant Mol Biol. 2003, 52: 143-159. 10.1023/A:1023993713849.
Reddy ASN, Ali GS, Celesnik H, Day IS: Coping with stresses: roles of calcium- and calcium/calmodulin-regulated gene expression. Plant Cell. 2011, 23: 2010-2032. 10.1105/tpc.111.084988.
DeFalco TA, Bendre KW, Snedden WA: Breaking the code: Ca2+ sensors in plant signaling. Biochem J. 2010, 425: 27-40. 10.1042/BJ20091147.
Tena G, Boudsocq M, Sheen J: Protein kinase signaling networks in plant innate immunity. Curr Opin Plant Biol. 2011, 14: 519-529. 10.1016/j.pbi.2011.05.006.
Eugelm T, Somssich E: Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007, 10: 366-371. 10.1016/j.pbi.2007.04.020.
Cho SM, Kang EY, Min KH, Lee YK, Kim YC, Yang KY, Kim KS, Choi YS, Cho BH: A positive role of watermelon ClWRKY70 gene for disease resistance in transgenic Arabidopsis thaliana. Biol Plantarum. 2012, 56: 560-565. 10.1007/s10535-012-0070-x.
Li J, Brader G, Kariola T, Tapio Palva E: WRKY70 modulates the selection of signaling pathways in plant defense. Plant J. 2006, 46: 477-491. 10.1111/j.1365-313X.2006.02712.x.
Li J, Brader G, Palva ET: The WRKY70 transcription factor: a node of convergence for jasmonate-mediated and salycilate-mediated signals in plant defense. Plant Cell. 2004, 16: 319-331. 10.1105/tpc.016980.
Knoth C, Ringler J, Dangl JL, Eulgem T: Arabidopsis WRKY70 is required for full RPP4-mediated disease resistance and basal defense against Hyaloperonospora parasitica. Mol Plant Microbe In. 2007, 20: 120-128. 10.1094/MPMI-20-2-0120.
Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Ann Rev Phytopathology. 2005, 43: 205-227. 10.1146/annurev.phyto.43.040204.135923.
Thaler JS, Owen B, Higgins VJ: The role of the jasmonate response in plant susceptibility to diverse pathogens with a range of lifestyles. Plant Physiol. 2004, 135: 1-9. 10.1104/pp.900109.
Pauwels L, Goossens A: The JAZ proteins: a crucial interface in the jasmonate signaling cascade. Plant Cell. 2011, 23: 3089-3100. 10.1105/tpc.111.089300.
Liu F, Jiang H, Ye S, Chen WP, Liang W, Xu Y, Sun B, Sun J, Wang Q, Cohen JD, Li C: The Arabidopsis P450 protein CYP82C2 modulates jasmonate-induced root growth inhibition, defense gene expression and indole glucosinolate biosynthesis. Cell Res. 2010, 20: 539-552. 10.1038/cr.2010.36.
Broekaert WF, Delauré SL, De Bolle MFC, Cammue BPA: The role of ethylene in host-pathogen interactions. Ann Rev Phytopathol. 2006, 44: 393-416. 10.1146/annurev.phyto.44.070505.143440.
Sestili S, Polverari A, Luongo L, Ferrarini A, Scotton M, Hussain J, Delledonne M, Ficcadenti N, Belisario A: Distinct colonization patterns and cDNA-AFLP transcriptome profiles in compatible and incompatible interactions between melon and different races of Fusarium oxysporum f. sp. Melonis. BMC Genomics. 2011, 12: 122-10.1186/1471-2164-12-122.
Collins NC, Thordal-Christensen H, Lipka V, Bau S, Kombrink E, Qiu J, Hückelhoven R, Stein M, Freialdenhoven A, Somerville SC, Schulze-Lefert P: SNARE-protein-mediated disease resistance at the plant cell wall. Nature. 2003, 425: 973-977. 10.1038/nature02076.
Schulze-Lefert P: Knocking on the heaven's wall: pathogenesis of and resistance to biotrophic fungi at the cell wall. Curr Opin Plant Biol. 2004, 7: 377-383. 10.1016/j.pbi.2004.05.004.
Frialdenhoven A, Peterhänsel C, Kurth J, Kreuzaler F, Schulze-Lefert P: Identification of genes required for the function of non-race-specific mlo resistance to powdery mildew in barley. Plant Cell. 1996, 8: 5-14.
Büschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, van der Lee T, Diergaarde P, Groenendijk J, Töpsch S, Vos P, Salamini F, Schulze-Lefert P: The barley Mlo gene: a novel control element of plant pathogen resistance. Cell. 1997, 88: 695-705. 10.1016/S0092-8674(00)81912-1.
Consonni C, Humphry M, Hartmann HA, Livaja M, Durner J, Westphal L, Vogel J, Lipka V, Kemmerling B, Schulze-Lefert P: Conserved requirement for a plant host cell protein in powdery mildew pathogenesis. Nat Genet. 2006, 38: 716-720. 10.1038/ng1806.
Bai YL, Pavan S, Zheng Z, Zappel NF, Reinstädler A, Lotti C, De Giovanni C, Ricciardi L, Lindhout P, Visser R, Theres K, Panstruga R: Naturally occurring broad-spectrum powdery mildew resistance in a central American tomato accession is caused by loss of Mlo function. Mol Plant Microbe In. 2008, 21: 30-39. 10.1094/MPMI-21-1-0030.
Humphry M, Reinstädler A, Ivanov S, Bisseling T, Panstruga R: Durable broad-spectrum powdery mildew resistance in pea er1 plants is conferred by natural loss-of-function mutations in PsMLO1. Mol Plant Pathol. 2011, 12: 866-878. 10.1111/j.1364-3703.2011.00718.x.
Humphry M, Bednarek P, Kemmerling B, Koh S, Stein M, Gobel U, Stüber K, Pislewska-Bednarek M, Loraine A, Schulze-Lefert P, Somerville S, Panstruga R: A regulon conserved in monocot and dicot plants defines a functional module in antifungal plant immunity. PNAS. 2010, 107: 21896-21901. 10.1073/pnas.1003619107.
Chen Z, Noir S, Kwaaitaal M, Hartmann A, Wu J, Mudgil Y, Sukuman P, Muday G, Panstruga R, Jones AM: Two seven-transmembrane domain MILDEW RESISTANCE LOCUS O proteins cofunction in Arabidopsis root thigmomorphogenesis. Plant Cell. 2009, 21: 1972-1991. 10.1105/tpc.108.062653.
Cheng H, Kun W, Liu D, Su Y, He Q: Molecular cloning and expression analysis of CmMlo in melon. Mol Biol Rep. 2012, 39: 1903-1907. 10.1007/s11033-011-0936-6.
Guo Y, Yu Y, Wang D, Wu C, Yang G, Huang J, Zheng C: GhZFP1, a novel CCCH-type zinc finger protein from cotton, enhances salt stress tolerance and fungal disease resistance in transgenic tobacco by interacting with GZIRD21A and GZIPR5. New Phytol. 2009, 183: 62-75. 10.1111/j.1469-8137.2009.02838.x.
Schlink K: Down-regulation of defense genes and resource allocation into infected roots as factor for compatibility between Fagus sylvatica and Phytophthora citricola. Funct Integr Genomics. 2010, 10: 253-264. 10.1007/s10142-009-0143-x.
Garcia-Mas J, Benjak A, Sanseverino W, Bourgeois M, Mir G, González VM, Hénaff E, Câmara F, Cozzuto L, Lowy E, Alioto T, Capella-Gutiérrez S, Blanca J, Cañizares J, Ziarsolo P, Gonzalez-Ibeas D, Rodríguez-Moreno L, Droege M, Du L, Alvarez-Tejado M, Lorente-Galdos B, Melé M, Yang L, Weng Y, Navarro A, Marques-Bonet T, Aranda MA, Nuez F, Picó B, Gabaldón T, Roma G, Guigó R, Casacuberta JM, Arús P, Puigdomènech P: The genome of melon (Cucumis melo L.). PNAS. 2012, 109: 11872-11877. 10.1073/pnas.1205415109.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.
Benjamini Y, Hochberg Y: Controlling the false rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995, 57: 289-300.
The International Cucurbit Genomics Initiative (ICuGI):http://www.icugi.org]website,
Saeed AI, Bhagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li J, Thiagarajan M, White JA, Quackenbush J: TM4 microarray software suite. Method Enzymol. 2006, 411: 134-193.
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, NJ: Humana Press, Totowa, 365-386.
Livak KJ, Schmitgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
This work was supported by a project of Universitat Politècnica de València (PAID-06-09). We also acknowledge the grant from the Spanish Ministry of Education and Science GEN2003-20237-C06 in the frame of which the microarray was developed and Plant KBBE project (SAFQIM: PIM2010PKB-00691).
The authors have declared that no competing interest exists.
CR, BP, AF and FN were involved in the conception and design of the study. CR and AF were responsible for growing melon plants and performed melon infection. CR prepared RNA extraction and qRT-PCR. CR and JH performed microarray statistical analysis. GR performed functional classification of transcripts and revised the manuscript. CR was responsible for drafting and revising the manuscript with contributions from co-authors. All read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Roig, C., Fita, A., Ríos, G. et al. Root transcriptional responses of two melon genotypes with contrasting resistance to Monosporascus cannonballus (Pollack et Uecker) infection. BMC Genomics 13, 601 (2012) doi:10.1186/1471-2164-13-601
- Jasmonic Acid
- Infected Root
- Roche NimbleGen
- Robust Multiarray Average