Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Gene expression profiling by cDNA-AFLP reveals potential candidate genes for partial resistance of ‘Président Roulin’ against Venturia inaequalis

BMC Genomics volume 15, Article number: 1043 (2014) Cite this article

Abstract

Background

Scab, caused by the fungus Venturia inaequalis, is one of the most important diseases of cultivated apple. While a few scab resistance genes (R genes) governing qualitative resistance have been isolated and characterized, the biological roles of genes governing quantitative resistance, supposed to be more durable, are still unknown. This study aims to investigate the molecular mechanisms involved in the partial resistance of the old Belgian apple cultivar ‘Président Roulin’ against V. inaequalis.

Results

A global gene expression analysis was conducted in ‘Président Roulin’ (partially resistant) and in ‘Gala’ (susceptible) challenged by V. inaequalis by using the cDNA-AFLP method (cDNA-Amplified Fragment Length Polymorphism). Transcriptome analysis revealed significant modulation (up- or down-regulation) of 281 out of approximately 20,500 transcript derived fragments (TDFs) in ‘Président Roulin’ 48 hours after inoculation. Sequence annotation revealed similarities to several genes encoding for proteins belonging to the NBS-LRR and LRR-RLK classes of plant R genes and to other defense-related proteins. Differentially expressed genes were sorted into functional categories according to their gene ontology annotation and this expression signature was compared to published apple cDNA libraries by Gene Enrichment Analysis. The first comparison was made with two cDNA libraries from Malus x domestica uninfected leaves, and revealed in both libraries a signature of enhanced expression in ‘Président Roulin’ of genes involved in response to stress and photosynthesis. In the second comparison, the pathogen-responsive TDFs from the partially resistant cultivar were compared to the cDNA library from inoculated leaves of Rvi6 (HcrVf2)-transformed ‘Gala’ lines (complete disease resistance) and revealed both common physiological events, and notably differences in the regulation of defense response, the regulation of hydrolase activity, and response to DNA damage. TDFs were in silico mapped on the ‘Golden Delicious’ apple reference genome and significant co-localizations with major scab R genes, but not with quantitative trait loci (QTLs) for scab resistance nor resistance gene analogues (RGAs) were found.

Conclusions

This study highlights possible candidate genes that may play a role in the partial scab resistance mechanisms of ‘Président Roulin’ and increase our understanding of the molecular mechanisms involved in the partial resistance against apple scab.

Background

Apple scab caused by the hemi-biotrophic ascomycete Venturia inaequalis (Cke.) Wint. is one of the most serious diseases of apple (Malus x domestica, Borkh.) worldwide, causing huge economic losses. Scab infection leads to deformation in the shape and size of fruits, premature leaf and fruit drop, and enhances susceptibility of the tree to chilling and freezing injuries [1]. Currently, multiple applications of fungicides are required for effective control in commercial orchards planted with susceptible cultivars. Depending on the year and region, as many as 18 to 29 fungicide treatments may be necessary in one season to control the disease. For apple orchards in France, pesticide treatments costs account for about 10% of the fixed production expenses, representing a substantial cost per kg of apple (0.031 €) [2]. This intensive use of fungicides in orchards raises ecological problems and human health concerns in addition to the economic cost.

An effective alternative to chemical control is the use of scab-resistant apple cultivars. Phenotypically, the effect of resistance genes against V. inaequalis has been showed to cover a continuum from complete immunity to near-susceptibility depending on genetic background, pathogen and environment [3]. Despite the great deal of gray area between the extremes [4] hampering the classification of some genotypes, apple scab resistance is often qualified either as complete, when the pathogen growth is fully inhibited (complete or qualitative resistance), or partial, when the resistance allows limited but significantly reduced pathogen growth as compared to a susceptible genotype (partial or quantitative resistance). Based upon the extent of pathogen growth and nature of symptoms, Chevalier et al. [5] classified the macroscopic foliar symptoms into four classes; classes from 0 to 3 are categorized as resistance responses while class 4 is considered to be a susceptible response. It is usually referred that complete resistance is determined by major resistance gene (R gene) while incomplete resistance involved multiple genes or loci of partial effect (QTLs, Quantitative Trait Loci).

R genes, typically providing high levels of resistance, are relatively easy to manipulate. For these reasons, they were extensively used in both basic research and applied breeding programs. To date, at least 17 major scab resistance genes have been identified and mapped across nine linkage groups of the apple genome [3]. For more than 50 years, one of these R genes, the Rvi6 (Vf) gene from Malus floribunda 821, has provided effective resistance against apple scab by allowing a reduction of 75% in the number of fungicide treatments [6]. Nevertheless, the use of single R gene-mediated resistance for crop protection is hampered by a lack of durability, particularly with pathogens having high evolutionary potential, as with V. inaequalis[7]. This ephemeral nature of R gene-mediated resistance is highlighted by the recent emergence of some races of V. inaequalis that are virulent on cultivars carrying the widely-deployed Rvi6 (Vf) gene [8]. In contrast to major genes, the performance of partially resistant cultivars in the orchard is a function of the gene effects and spectra, which is thought to be more durable than single R gene-mediated resistance [9]. This durability could be explained by the smaller effects of partial resistance that result in a lower selection pressure on the pathogen and/or its presumed broader spectrum. Also, because partial resistance is controlled by multiple genes, pathogen isolates that overcome one of the genes would gain only a marginal advantage [4]. In apple, partial resistance has been mapped as QTLs to 10 out of the 17 linkage groups of the genome [10–14].

Extensive efforts have been made to clone and characterize major scab R genes, and the downstream responses that they trigger have become better understood. R genes typically encode proteins that recognize pathogen effectors or modifications of plant proteins that are targets of those effectors [4]. In this respect, the Rvi6 (Vf) resistance locus revealed the presence of a cluster of four resistance gene paralogs (called HcrVf genes), similar to the tomato Cf resistance gene, encoding leucine-rich repeats receptor-like proteins (LRR-RLP) [15] and the Rvi15 (Vr2) was reported to contain three TIR-NBS-LRR genes [16]. The function of all these genes was analyzed and only two of them, HcrVf2[17, 18] and Vr2-C[19] for Rvi6 and Rvi15, respectively, were proven to be functional against V. inaequalis. Literature is not in agreement on Vf1a (syn. HcrVf1) function [18, 20]. Recognition of V. inaequalis by these proteins triggers downstream defense reactions involving putative LRR receptor-like protein kinases [21, 22] and several defense-related proteins, such as b-1,3-glucanase, ribonuclease-like PR-10, cysteine protease inhibitor, endochitinase, ferrochelatase, and ADP-ribosylation factor [23, 24]. Methallothionein may also play a role in plant defense against V. inaequalis as it is present in large amounts in the apoplast of the Rvi6 cultivar ‘Remo’. Finally, recent publications explored the network of defense response triggered in Rvi6 (HcrVf2)-transformed apple plants using wide genome expression techniques, such as the PCR-based suppression subtractive hybridization and the dHPLC for cDNA- Amplified Fragment Length Polymorphism (cDNA-AFLP) transcript profiling [25, 26]. These studies gave new insight into the understanding of the plant pathogen-interaction that results in complete scab resistance. Nevertheless, the function of genes underlying the QTLs for partial apple scab resistance remains unknown. They are believed not to be based on pathogen recognition systems, as it is the case with most major R genes, but the possibility cannot be excluded [27].

The partially resistant apple cultivar ‘Président Roulin’ is an old Belgian cultivar that is used in apple breeding programs of the Walloon Agricultural Research Center (CRA-W) to broaden genetic apple scab resistance, and therefore reduce the risk of resistance breakdown. Under heavy scab infection, this cultivar shows typical resistance symptoms of chlorosis and necrosis with only slight sporulation (Chevalier class 3a). Its resistance against V. inaequalis has been durable for over 25 years in scab evaluation in different orchards of the CRA-W in Belgium without any fungicide treatment [28]. The partial resistance has been shown to be polygenic, but highly heritable [9].

In this study we investigated the defense response of ‘Président Roulin’ by identifying genes differentially expressed between this cultivar and the susceptible ‘Gala’ cultivar after pathogen challenge. For this purpose, cDNA-AFLP technology was chosen as it allowed a survey of transcriptional changes with no prior assumptions about which genes might be involved in the plant response [29]. This technique constitutes a robust solution for differential display, detecting changes in gene expression between samples, and it has been successfully applied in several quantitative expression studies in apple [26, 30–33]. The genes identified in this study were annotated and sorted into Gene Ontology (GO) categories. Comparisons were then made with various Malus x domestica libraries: two EST libraries from uninfected young expanding leaves [34, 35] and a cDNA library from an apple susceptible line transformed with the major scab R gene Rvi6 (HcrVf2) [25]. Common as well as different defense pathways were revealed and are discussed. Finally, we checked for co-localizations of our differentially expressed genes with resistance gene analogues (RGAs), QTLs and mapped R genes for apple scab resistance published in the literature.

Results

Fungal development across post inoculation time points

Microscopic observations revealed no significant difference between ‘Président Roulin’ (partially resistant) and ‘Gala’ (susceptible) for pathogen development at the early stages of infection. This comprises conidia germination, formation of appressoria and cuticle penetration (Figure 1). On both cultivars, conidia germination began within 4 hours post inoculation (hpi), and at 16 hpi most of the conidia had produced germ tubes. At this time, the rudimentary germ tubes had bulged at the tip to form characteristic appressoria adhering to the cuticle. At the end of the 24 hour period, the process was further advanced. Shortly after the formation of the appressoria, infection hyphae appeared and penetrated the host. At 48 hpi, 80% of the appressoria were formed and invasion of the host plant started with the expansion of the primary stroma. With the staining method used in this study, the subcuticular growth of the fungus was difficult to observe because it was poorly stained. Nevertheless, we can assume that, when the stroma was visible, at 120 hpi, subcuticular development on the susceptible host, ‘Gala’, had significantly exceeded that of the partially resistant ‘Président Roulin’ (data not shown). This difference in the extent of colonization of the fungus between the two cultivars remained throughout the whole infection cycle of the fungus. Between 7 and 12 days after inoculation, apple scab symptoms became macroscopically visible on ‘Président Roulin’ and ‘Gala’. After 21 days, 90% of the leaf surface of ‘Gala’ (susceptible) was covered by sporulating apple scab lesions (class 4) [5]. However, typical chlorotic and necrotic lesions with slight sporulation were observed on leaves of ‘Président Roulin’ (partially resistant), covering less than 15% of the leaf surface. These symptoms were considered as resistance responses and were classified in the class 3b, as described by Chevalier et al. [5].

Figure 1
figure 1

Kinetics of V. inaequalis conidial development on ‘Président Roulin’ and ‘Gala’ leaves. Germination of conidia and formation of appressoria were observed under light microscopy over time post inoculation. Fungal tissues were stained on whole leaves with periodic acid-basic fuchsin according to the method of Preece [97].

Full size image

cDNA-AFLP fingerprints: optimization and identification of differentially expressed transcripts

We used the AFLPinSilico application to choose the optimal restriction enzyme combination for the cDNA-AFLP experiments. The results of the different parameters for each enzyme pair combination are shown in Table 1. Enzymes with 4- or 5-base recognition sites yielded the highest number of Transcript Derived Fragments (TDFs), although these were generally relatively short and highly redundant, with up to 5 cleavage sites per cDNA. In opposite, the EcoRI/MseI recognized cleavage sites on less than half (34%) of the apple full-length cDNA tested, but this combination of restriction enzymes generated more informative TDFs than all the other enzyme pairs tested. In fact, EcoRI provided the longest TDF (mean size of 234 bp) with the lowest redundancy (0.5 restriction sites per cDNA) and derived at least partially from coding regions (683 bp from the poly(A) + tail). This enzyme was therefore the most appropriate and it was chosen in combination with MseI for cDNA-AFLP analysis.

Table 1 Suitability of restriction enzymes for use in cDNA-AFLP analysis in apple
Full size table

The cDNA-AFLP analysis using 141 primer combinations in ‘Président Roulin’ and ‘Gala’ resulted in 30 to 100 TDFs per primer combination, depending on the number of additional bases used for the selective amplification step, and a total of about 10,250 TDFs in each cultivar (representing a total of 20,500 for both cultivars). TDFs ranged in size from 30 to 800 bp. Figure 2 shows a typical cDNA-AFLP profile of the two apple genotypes challenged by V. inaequalis or water. Considering that 123 primer pairs with two additional selective nucleotides (EcoRI + 2/MseI + 2) were tested out of the 256 possible primer pair combinations, and taking into account that about 40% of the apple cDNA could potentially be visualized with the restriction enzyme EcoRI and MseI, we estimated that we analyzed a representative sample of approximatively 19% of the apple genes expressed in the tissues.

Figure 2
figure 2

Expression patterns of apple genes displayed by cDNA-AFLP fingerprints. The cDNA-AFLP compares transcriptional profiles from ‘Président Roulin’ (partially resistant) and ‘Gala’ (susceptible) mock-inoculated or challenged by V. inaequalis at 48 hpi. The 32 samples are arranged in 8 groups according to the different specific primers tested during the selective amplification step of the AFLP procedure. E and M refer to the EcoRI and MseI primers, followed by the selective nucleotides used. Within each of the 8 groups samples are ordered as follows: ‘Président Roulin’ infected (Ri) and mock-inoculated (Rm), and ‘Gala’ infected (Si) and mock-inoculated (Sm). Differentially expressed TDFs were classified into 2 categories: genotype-specific TDFs (group I) and pathogen-responsive TDFs (group II), further divided into two sub-groups; pathogen-responsive TDFs expressed in common by both genotypes (sub-group IIa) and pathogen-responsive and genotype specific TDFs (sub-group IIb). Illustrations are given.

Full size image

At the stringent threshold used, only 4.6% of all the 20,500 generated fragments exhibited significant differences in intensity among the different genotypes or treatments. According to their banding patterns, the differentially expressed TDFs were classified into two groups: (I) genotype-specific TDF, for all banding patterns differing between the genotypes but not affected by fungal infection and (II) pathogen-responsive TDFs representing all the TDFs showing differential expression during infection. Group II was further divided into two subgroups: (IIa) pathogen-responsive TDFs, representing the TDFs with differential expression upon pathogen challenge in both genotypes and (IIb) pathogen-responsive and genotype-specific TDFs, representing the TDFs showing differential expression induced by the pathogen in one of the two cultivars.

The first group (group I, genotype-specific) accounted for about 1.7% of all generated fragments (232 TDFs present only in ‘Président Roulin’ and 115 in ‘Gala’). In contrast, 281 (230 up-and 51 down-regulated) and 311 TDFs (241 up- and 70-down regulated), respectively, were identified in ‘Président Roulin’ and ‘Gala’ as significantly differentially expressed after fungal infection (Figure 3). These pathogen-responsive TDFs (group II) accounted for about 2.9% of all the 20,500 TDFs analyzed for both cultivars. Among them, 125 (111 up-regulated and 14 down-regulated) overlapped between the two genotypes (subgroup IIa). The remaining 156 and 186 TDFs, for ‘Président Roulin’ and ‘Gala’, respectively, displayed differential expression after fungal attack that depended on the genotype. They were classified into subgroup IIb (pathogen-responsive and genotype-specific) and represented about 1.7% of the fragments. From our point of view, these bands differentially expressed in ‘Président Roulin’ only, are the most interesting as they could be involved in specific plant defense reaction against V. inaequalis. In fact, genes that showed only a genotype effect may reflect the effect of the genetic background, whereas genes exhibiting only a treatment effect may represent a general plant response to pathogen challenge that does not determine the final different phenotype.

Figure 3
figure 3

Venn diagram showing number of pathogen-responsive TDFs in ‘Président Roulin’ (partially resistant) and/or ‘Gala’ (susceptible). Group II of TDFs was classified into two sub-groups: pathogen-responsive TDFs expressed in common by both genotypes (sub-group IIa) and pathogen-responsive and genotype specific TDFs (sub-group IIb). ‘+’ and ‘-’ represent up- and down-regulation, respectively.

Full size image

We then considered the amplitude and direction of expression changes for all pathogen-responsive TDFs (group II). We plotted the log10 transformed expression ratios of ‘Président Roulin’ against ‘Gala’ and distinguished the TDFs that were differentially expressed in only one of the two lines (blue squares for ‘Président Roulin’ and red squares for ‘Gala’, Figure 4) from those being differentially expressed by both cultivars (green triangle, intersect of the two circles in the Venn diagram of Figure 3). This graph showed the overall similarity and specificity of gene expression in the TDFs differentially expressed in common by ‘Président Roulin’ and ‘Gala’: most of these TDFs were differentially regulated in the same direction (up or down-regulation) by both cultivars. Fold-changes from 2 to 70 were observed, with the majority of the TDFs showing a difference in expression less than 7-fold.

Figure 4
figure 4

Scatter plot of log10-gene expression fold changes of pathogen-responsive TDFs from ‘Président Roulin’ and ‘Gala’. Fold changes are relative to mock inoculation. Colour-coded plots represent TDFs differentially regulated in one of the cultivars (cultivars-specific) or in both cultivars (common). Log ratios >0 or <0 indicate up- or down-regulation, respectively, dashed lines set at 0.3 and −0.3 (+ − log102) correspond to the threshold of a two-fold change in gene expression that was used as the cut-off value for biological significance.

Full size image

Sequences annotation

Two-hundred and fifty nine TDFs out of the 281 pathogen-responsive TDFs from the resistant cultivar and part of the genotype-specific TDFs (131 TDFs) were excised from the gels, re-amplified and cloned (390 TDFs in total). No bands belonging to the susceptible cultivar ‘Gala’ were excised and a priority was given to TDFs of sufficient length (upper part of the gel), to facilitate their functional characterization. Two clones were sequenced for each re-amplified TDF and for 38% of them (146 TDFs) a different sequence between the two clones was found, leading to a total of 536 sequences. To limit redundancy, sequences were clustered using the software Egassembler [36, 37] and then compared to the whole apple genome sequence assembly in order to identify the unigenes [38]. This returned 497 unique sequences (29 contigs and 468 singletons) from 53 bp to 803 bp that were annotated and submitted to the NCBI database with accession numbers assigned (Additional file 1).

Among the 497 sequences, 69% (344 TDFs) were similar to known expressed sequences in public databases (319 could be annotated, 25 were similar to encoded proteins with unclear function), while 153 sequences had no matches in the NCBI database. With the exception of 15 sequences similar to genes belonging to the apple stem pitting virus (ASPV), all TDFs were similar to sequences from model plants, with Vitis vinifera, Ricinus communis and Populus trichocarpa being the three most represented species. As expected, most TDFs (80%) were found to be similar (BLAST E-value <1e-4) to Malus x domestica sequences derived from the whole genome sequence of apple [38]. The ASPV is a widespread pathogen of pome fruit usually transmitted by the rootstock and frequently latent in the host [39]. As the cDNA encoding for ASPV coat proteins was found in both inoculated and controlled plants and as no symptoms of ASPV were observed in apple trees used for the experiments we believe that the virus was latent and did not interfere with gene expression data.

Putative functions of the apple TDFs were then classified according to the GO vocabulary. Figure 5 shows the percentages of apple genes assigned to the biological processes. About 27% of the annotated sequences have metabolic roles (particularly in protein and carbohydrate metabolism) and 24% are involved in cellular processes. Other relevant groups, accounting for 12, 10 and 7% of the TDFs, respectively, include single-organism process, response to stimulus (particularly response to biotic stress, oxidative stress, response to wounding and defense response, listed in Table 2) and biological regulation. In addition to TDFs related to plant defense response according to the GO classification, we also found TDFs showing similarities with genes encoding for proteins reported to have a potential role in the general defense response pathway through careful analysis of the scientific literature. These TDFs are also potential candidate genes for partial resistance of ‘Président Roulin’ against V. inaequalis. These TDFs encoded proteins involved in the recognition of the pathogen, signal transduction, oxidation reduction process, photosynthesis, transport, response to environmental stress, metabolism, transcription, and cell wall organization (Table 3). Expression of TDFs of interest was further investigated by qRT-PCR and putative functions of some of these proteins in the partial resistance against V. inaequalis are discussed.

Figure 5
figure 5

Distribution of differentially expressed TDFs within the GO categories of biological processes. GO annotations were made according to the International Gene Ontology Consortium using the automatic bioinformatics software Blast2GO.

Full size image
Table 2 TDFs associated with plant defense response, response to oxidative stress and response to wounding
Full size table
Table 3 Gene expression analysis of selected TDFs by qRT-PCR in ‘Président Roulin’ (resistant) and ‘Gala’ (susceptible)
Full size table

Gene enrichment analysis

The gene enrichment analysis revealed that genes accounting for the different functional GO categories were similarly represented (P > 0.05) in two subsets of our cDNA-AFLP: the TDFs being up- or down-regulated in ‘Président Roulin’ after pathogen attack (group II + and II - ). Nevertheless, when compared to two different collections of ESTs [34, 35] from uninfected expanding apple leaves (AELA and Mdstw), some biological pathways appeared to be significantly over-represented (P < 0.05) in the ‘Président Roulin’ cDNA-AFLP library (Figure 6). In both comparisons, it was notably the case for genes involved in photosynthesis and in response to stress. Additionally, other GO categories were significantly over-represented in ‘Président Roulin’ but depended on the library being compared: e.g. oxidation reduction process, regulation of gene expression in the AELA library, and homeostatic process in the Mdstw library. All these genes could be involved in the general defense response pathway of ‘Président Roulin’ against V. inaequalis. No biological pathways appeared to be significantly under-represented in our collection.

Figure 6
figure 6

Over-representation of GO categories in ‘Président Roulin’ cDNA-AFLP library compared to non-infected EST apple libraries. The comparison has been made by gene enrichment analysis for Biological Process GO categories between our cDNA-AFLP library from scab-infected leaves of ‘Président Roulin’ (partially resistant) and two EST libraries from uninfected actively growing shoot of: (A) cultivar ‘Royal Gala’ in the library AELA [34] and (B) cultivar ‘Wijcik’ in the library Mdstw [35]. Gene enrichment analysis was conducted with the software Blast2Go using Fisher’s Exact Test at a p-value <0.05. No GO categories were shown to be under-represented.

Full size image

In order to assess if mechanisms leading to partial resistance against V. inaequalis might involve biological pathways different to a complete resistance response, we compared by gene enrichment analysis our ‘Président Roulin’ TDFs dataset to the SSH library from Rvi6 (HcrVf2)-transformed ‘Gala’ line challenged by V. inaequalis[25]. Some biological processes were over-represented in our partial resistance compared to the complete resistance mediated by Rvi6 (HcrVf2) gene (Figure 7, Additional file 2). This was the case for genes involved in the plant response to DNA damage stimulus, particularly DNA repair, and those involved in the regulation of hydrolase activity. A significant over-representation of differentially regulated genes classified as ‘regulation of defense response’ was also found in ‘Président Roulin’. Again, no under-representation of any biological pathways was found.

Figure 7
figure 7

Over-representation of GO categories in ‘Président Roulin’cDNA-AFLP library compared to Rvi6 ( HcrVf2 )-’Gala’ transformed library. The comparison has been made by gene enrichment analysis for the Biological Process GO category between our cDNA-AFLP library (scab-infected leaves of ‘Président Roulin’, partially resistant) compared with a cDNA library from completely resistant Rvi6 (HcrVf2)-transformed ‘Gala’ lines challenged with V. inaequalis[25]. Gene enrichment analysis was conducted with the software Blast2Go using Fisher’s Exact Test at a p-value <0.05.

Full size image

Quantitative real-time reverse transcriptase PCR

To validate the reliability of our cDNA-AFLP analysis, qRT-PCR was performed on 24 pathogen-responsive TDFs representative of almost all functional categories identified, with a preference for defense-related genes and genes possibly involved in pathogenesis. Table 3 presents separately the qRT-PCR results carried out on the same RNA samples used for the cDNA-AFLP analysis and on a biological replication of the experiment. Results were expressed as fold-change of ‘Président Roulin’ (partially resistant) and ‘Gala’ (susceptible) after pathogen attack, in respect to mock-inoculated leaves. In general, expression data provided by qRT-PCR were in good agreement with profiles detected by cDNA-AFLP. When performed on RNA samples used for cDNA-AFLP, qRT-PCR confirmed the pattern of gene expression of 20 TDFs (83%). Expression of 22 of these RNA samples was then further verified by qRT-PCR on a biological replicate of the experiments, and the pattern of expression of 13 of these samples was confirmed (59%). TDFs that were not in accordance with the cDNA-AFLP showed no significant changes in expression in ‘Président Roulin’ (i.e. 42BUHcrVf, 43CU118, 39 AU13). For three TDFs, strong changes in gene expression (more than five-fold) were detected in infected leaves of ‘Président Roulin’ in both experiments (51HU129’, 44GU182 and 44EU122). Two different TDFs (43DU149 and 43DU149’) cloned from the same band and showing a significant increase in intensity after pathogen challenge were both confirmed to be up-regulated after pathogen attack. In most cases, no significant changes of expression were detected in infected leaves of the susceptible ‘Gala’ cultivar.

Co-localization of the TDFs with RGA, QTL or apple scab major R genes

Approximatively, 40% of the TDFs anchored in silico on one of the 17 chromosomes of apple were localized in the proximity (within 250 kb) of RGAs clusters, QTLs or major R genes for apple scab resistance (Table 4). Nevertheless, this frequency was not significantly greater than those calculated for random ESTs derived from an uninfected apple leaves library [34]. So co-localization could have occurred purely by chance. However, considering separately the three classes of TDFs (group I, IIa and IIb), we found that group I and IIb of TDFs mapped at a greater frequency in the vicinity of major R genes than EST from AELA library (5% instead of 1%).

Table 4 Frequency of co-localization of TDFs from ‘Président Roulin’ with RGAs, QTLs and major apple scab R genes
Full size table

Discussion

Plant disease resistance and susceptibility are governed by the combined host and pathogen genotypes, and depend on a complex exchange of signals and responses occurring under given environmental conditions. During the long process of host-pathogen co-evolution, plants have developed various elaborate mechanisms to ward off pathogen attack [40]. In addition to constitutive defense, it has been postulated that a key difference between resistant and susceptible plants is the timely recognition of the invading pathogen, and the rapid and effective activation of host defense mechanisms. Such induced resistance mechanisms have been demonstrated at transcriptional level in numerous studies on plant-pathogen interactions involving either complete [26, 41–44] or partial disease resistance [45, 46]. In both types of resistance, pathogen attack was accompanied by activation of host plant response genes and accumulation of corresponding gene products. Based on these findings, our cDNA-AFLP study attempted to elucidate the molecular mechanisms underlying partial resistance against apple scab by the identification of differentially expressed transcripts between ‘Président Roulin’ (partially resistant) and ‘Gala’ (susceptible) after pathogen attack.

Effective identification of differentially expressed transcripts after scab infection requires the determination of the stage of pathogen development at which it is first affected by the host’s defense mechanisms. This is important particularly in investigations built upon the hypothesis that resistance involves the induction of defense pathways. To identify the appropriate time-point for the extraction of RNA, we compared the kinetic of pathogen development on leaves of ‘Président Roulin’ (partially resistant) and ‘Gala’ (susceptible). Light microscopic assessment identified 48 hpi as the critical time-point when the pathogen is first affected by the host’s responses during formation of appressoria and stroma. In fact, whereas the early stages of pathogen growth occurred at the same rate in both cultivars, significant differences were observed in later stages, when development of subcuticular stroma subsequent to cuticle penetration was reduced by greater extent in the partially resistant cultivar compared to the susceptible one (Figure 1). These differences are concordant with the literature with various susceptible and resistant apple cultivars [47, 48] and led to a reduced disease severity and sporulation on ‘Président Roulin’ 21 days after infection. This also suggests that no particular pre-formed chemical or structural defense barriers are present at the leaf surface of ‘Président Roulin’. The choice of 48 hpi as the key time for our gene expression study is further motivated by the fact that in various plant-pathogen interactions the accumulation of defense-related transcripts in the host is often concomitant to the formation of appressoria and the development of the intracellular hyphae [49–52]. When this cytological event occurs, the molecular dialog between the pathogen and the host presumably begins. The host detects the pathogen and initiates defense-related transcriptional responses [53].

As for any transcriptome study, a genome-wide screen for differentially expressed genes based on cDNA-AFLP requires that as many transcripts as possible are analyzed in a unique way, and that the data obtained are informative enough to allow characterization of the transcripts. Considering criteria such as the genomic region of the restriction sites (coding or non-coding region), the redundancy and the length of the TDFs being produced, we showed in silico that the restriction enzyme pair EcoRI/MseI was the most appropriate for our cDNA-AFLP analysis. This restriction enzyme combination was successfully used in similar experiments on apple [25, 26, 32, 54, 55]. In our study, this enzyme combination yielded a very high informative content (69% of the TDFs could be effectively annotated) with limited redundancy (sequence assembly yielded only 29 contigs out of the 536 unclustered sequences).

The cDNA-AFLP method has the advantage that all treatments, time points and genotypes under investigation can be concurrently compared [56, 57]. In our study, we estimated that we analyzed a representative sample of approximatively 19% of the genes expressed in the leaves of ‘Président Roulin’ and ‘Gala’ in response to V. inaequalis infection at 48 hpi. Our estimation corroborates with the data from the recent sequencing of the apple genome project since our 10,250 bands would represent about 18% of the 57,386 putative genes in apple [38]. We highlighted that the differences between partially resistant and susceptible responses to V. inaequalis, were more quantitative (different regulation of expression of the same genes in both cultivars) than qualitative (different genes being expressed in both cultivars), and that quantitative differences mainly involved up-regulation of genes (71% of up-regulations). As similarly demonstrated in the literature [58, 59], we can argue that, despite the very clear phenotypic differences observed between the partially resistant and susceptible responses, the timing and amplitude of gene induction could have a greater influence on the disease resistance status of the plant than the identity of genes activated. Under-representation of down-regulated TDFs has also been reported in apple - V. inaequalis interaction [25] and in other pathosystems [60–62] but it has been suggested that such under-representation could be mainly due to the two-step PCR of the cDNA-AFLP method [61]. The overall similarity of transcriptional signature between ‘Président Roulin’ and ‘Gala’ is further demonstrated by the observation that genes were regulated in the same direction in both cultivars: when a TDF was up-regulated in one cultivar, up-regulation was also observed in the other one (and inversely) (Figure 4). Finally, all ‘Président Roulin’ pathogen-responsive TDFs were also found to be expressed in the mock-control from high, to intermediate or low levels, supporting the idea that many genes involved in the defense reaction were already expressed in resistant genotype before pathogen challenge [26, 63].

The potential presence of multiple sequences in one band cut from the gel is the major disadvantage of cDNA-AFLP fingerprinting. In our study, this mis-cloning has been observed for 38% of the excised bands and could be an underestimate of the multiplicity of sequences per band since only two clones per band were analyzed. Each restriction enzyme pair may in fact result in different fragments of the same size which are visualized as a single band on a cDNA-AFLP gel autoradiography. Therefore, confirmation of gene expression by an independent technique is required, as for hybridization-based assay like microarray. In this study, expression profiles observed by cDNA-AFLP were confirmed by qRT-PCR in 83% of the 24 TDFs under examination (on the same cDNA samples used for the AFLP, Table 3). Results were similar to those reported by Baldo et al. (79% of 28 ESTs) [55] and Paris et al. (81% of 32 TDFs) [26]. The percentage of confirmation decreased when using an independent biological replication (59%). According to these results it can be assumed that the modulation of about 83% of the observed TDFs is confirmed; however extending data to at least a biological repetition is highly recommended to further confirm that gene modulation strictly depends on the biological system under study (i.e. the plant-pathogen interaction). This consideration is valid not only for cDNA-AFLP but for all transcriptional analysis as modulation of gene expression can be influenced by other factors (environment, biotic stresses…). However, in our opinion, two or more sequences from one band can contribute to the change of band intensity observed on the autoradiography, as suggested by the fact that two different sequences cloned from the same band were both confirmed to be up-regulated after pathogen attack (43DU149 and 43DU149’).

The GO analysis of differentially expressed TDFs revealed that they are represented by a high diversity of functional categories (Figure 5). This is not surprising since, with the emerging of genome-wide gene expression profiling technologies, it is now clear that plant response to pathogens is associated with massive changes in gene expression. For example, in an Arabidopsis microarray, more than 2000 genes (out of 8000 genes) involved in a broad range of biological responses were regulated in response to the bacterial pathogen Pseudomonas syringae[64]. In our study, the annotation of ‘Président Roulin’ pathogen-responsive TDFs indicates that they may act in early events of plant defense response such as pathogen recognition (e.g. some TDFs encoded for putative nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins [65]), or in signal transduction (e.g., mitogen-activated protein (MAP) kinase [66]). Other TDFs identified may play a role in the later stages of the ‘Président Roulin’ defense response with the induction of genes aiming to stop or reduce the invasion of the host by pathogens. In plants, this often leads to the hypersensitive response, a form of localized programmed cell death orchestrated by the oxidative burst. TDFs involved in such reactions encoded for gluthatione-S-transferase [67], peroxidase [68], E3 ubiquitin protein ligase [69] and cysteine protease/cysteine protease inhibitor [70].

In order to understand the global molecular pathways involved in partial resistance, gene enrichment analysis were conducted to statistically determine whether specific GO terms were enriched in different sets of cDNA. First we demonstrated that the same functional categories were involved in our up- and down-regulated set of TDFs (group II + and II -). These results differ from the findings of Paris et al. [26] who found that TDFs similar to genes putatively involved in defense responses were generally up-regulated in R gene-mediated resistance (Rvi6/HcrVf2), while those putatively involved in general metabolism were down-regulated. In contrast, we demonstrated that genes involved in stress response (including disease/defense response) and in photosynthesis were preferentially regulated in the partially resistant cultivar ‘Président Roulin’ after pathogen attack, as compared to two published EST libraries from Malus x domestica uninfected young expanding leaves [34, 35] (Figure 6). When stress response genes could participate in the partial resistance of ‘Président Roulin’ against V. inaequalis, regulation of genes of the photosynthetic pathway might be the first step towards the appearance of chlorotic spots on infected leaves, which in the apple - V. inaequalis incompatible interaction appear about 8–12 days after inoculation [25]. It is also well known that defense responses is energy intensive [73] and requires transcriptional activation of genes [71]. These could be the reason why we found in ‘Président Roulin’ an over-representation of genes involved in different catabolic/metabolic processes, that might ‘fuel’ the implementation of downstream defense response, and in the regulation of gene expression. Nevertheless, as these later GO categories appear to be significantly over-represented in only one of the two libraries being compared, their involvement in the resistance of ‘Président Roulin’ against V. inaequalis still has to be confirmed. Finally, our comparison between the ‘Président Roulin’ cDNA library with the completely resistant Rvi6 (HcrVf2)-transformed ‘Gala’ lines [25] revealed both similar events (e.g. transport, photosynthesis functions), and differences in some biological process categories (Figure 7). Among the differences (Additional file 2) we noticed an over-representation in ‘Président Roulin’ of genes involved in (1) the regulation of defense response (e.g. E3 sumo ligase SIZ1 proteins [72] and tocopherol cyclase [74, 75]; (2) proteins involved in the regulation of hydrolase activity, also recognized to be key enzymes in the regulation of programmed cell death in incompatible plant-pathogen interactions (e.g. cysteine protease inhibitor and a 13-fold repressed cysteine protease [70, 76, 77]); (3) proteins involved in the response to DNA damage stimuli. This latter category is thought to be involved in the plant response to abiotic stress such as UV-B [78] and osmotic stress [79], but to our knowledge, these proteins were not yet known to be involved in quantitative disease resistance.

Beside these slight differences between complete and partial resistance, large parts of the transcriptional signatures did not demonstrate enrichment for genes in particular functional categories. This finding is consistent with the hypothesis that partial resistance could be due in part to the same genes governing complete resistance. To illustrate this hypothesis, some classical R genes encoding for the NBS-LRR family protein were found to be up-regulated in the partially resistant cultivar ‘Président Roulin’ after pathogen attack. In that context, partial resistance could be due to defective R genes that recognize with low efficiency pathogens and trigger weak defense response [65]. This may result either from mutation in the R genes themselves or at the corresponding avirulence locus of the pathogen. In fact, there are compelling lines of evidence that allelic variants of R genes account for quantitative disease resistance in plants (e.g. Xa21 in rice-blast interaction [80]). In the same way, when a pathogen strain overcomes an R gene due to a mutation at the avirulence locus, it has been proved that the “defeated” R gene still has a residual effect and can act as a QTL against virulent strains of the pathogen. This phenomenon has been observed for the major resistance genes Rvi4 (Vh4) [3] and Rvi6 (Vf) [81] in the V. inaequalis - apple interaction, and also in rice bacterial blight [82], wheat stem rust [83] and powdery mildew pathosystems [84]. Likewise, in the same way as major resistance genes, several QTLs have been shown to be isolate-specific [12, 85, 86] and co-localization of QTLs and R genes has been noted in several species including apple [62, 87–89]. Another evidence suggesting that genes controlling partial resistance in ‘Président Roulin’ could share structural and functional similarities with R genes resides in the fact that some subsets of our cDNA library mapped at greater frequency in apple genomic regions known to carry major scab R genes [3, 14] (as compared to random EST from uninfected apple library). The co-localization of ESTs with genomic regions carrying disease resistance factors (R genes, QTLs or RGAs) has been already reported in various genome-wide analyses studies [90–93]. No significant co-localization of our cDNA library with apple scab QTLs [3, 13, 14] nor apple RGAs [94] has been found. Obviously, we could only check for co-localizations with R genes and QTLs that have been detected so far. Moreover, information on the genomic loci that can regulate the expression level of the TDF of interest is still lacking from this analysis. In fact, the measured mRNA levels can either be the product of regulation of the parent gene or of another gene, mapping somewhere else in the genome (cis- or trans- regulatory elements) [95].

Conclusions

In conclusion, this study provides a wide transcriptional profile analysis for the comprehension of key events in partial resistance of ‘Président Roulin’ and highlights possible candidate resistance genes. We found altered gene expression in resistant and susceptible plants in response to V. inaequalis that involved many functional categories. Genes acting in pathogen recognition (NBS-LRR) as well functioning downstream of the initiated defense signaling pathways were identified. Biological processes related to stress response and photosynthesis were found to be over-represented in infected leaves of the partially resistant cultivar compared with two published libraries of uninfected apple leaves. In addition, through comparison between partial and complete resistance, the pathogen-responsive cDNA library revealed common physiological events, but differences in regulation of defense response, in the regulation of hydrolase activity, and in response to DNA damage stimuli. Finally, TDFs from ‘Président Roulin’ mapped more frequently in the vicinity of major R genes for apple scab resistance, suggesting that quantitative and complete resistance could be governed by the same types of genes. A functional analysis of the differentially expressed genes will allow more insights into their possible role in the quantitative resistance reaction of ‘Président Roulin’ against V. inaequalis. For example, an assessment of the differential expression of candidate resistance genes over different time points after infection could be investigated to find out how resistance is regulated by quantitative and/or kinetic enhancements. Also, analysis of candidate gene expression data in a segregating population could infer causal relationship between the differential expression of the genes and the resistance phenotype of the progeny. Finally, these candidate resistance genes might be at the basis of the development of molecular marker tools to be used in a genome-informed breeding program to speed-up the selection process of resistant plants.

Methods

Plant material and inoculation with Venturia inaequalis

Plants of the partially resistant Belgian cultivar ‘Président Roulin’ and the susceptible cultivar ‘Gala’ were grafted on M9 rootstocks and grown in a greenhouse at 20°C under 16 hours of illumination by daylight-incandescent lights. In the frame of the DARE European project, ‘Président Roulin’ has been shown to be resistant to a large range of inocula, including local mix inocula and monoconidal V. inaequalis strains belonging to the race 1, 6, 7 [96] and 2, 8, 9 (data not showed).

In this study we used six monoconidial strains of race 1 V. inaequalis originating from the INRA collection at Angers, France (EU-B04, EU-B16, EU-D49, EU-F05, EU-F11 and EU-I09) to prepare the inoculum. Each strain, first grown in Petri dishes for 10 days on malt extract agar and covered by a cellophane membrane, was multiplied separately on young seedlings raised from open-pollinated ‘Golden Delicious’. Infected leaves were dried and stored at −20°C for not more than 1 year before use. Conidia were harvested by shaking the leaves in sterile water. A conidial suspension was prepared by mixing conidia of the 6 strains at a final concentration of 2.5 × 105 conidia per milliliter. The conidial suspension was sprayed onto young leaves of actively growing ‘Gala’ and ‘Président Roulin’ plants in quantities sufficient to form small droplets on the leaf surface before run off. The inoculated plants were incubated at 20°C under maximum relative humidity (RH) for two days and were then transferred to the greenhouse (20°C, 60-80% RH). Control plants were inoculated with sterile water. Conidia vitality was verified after 24 h by determining the fraction of germinated conidia (~70% for all the isolates). Levels of scab infection were scored for each plant 21 days after inoculation. To reduce the level of biological variation among samples, two plants per treatment were used. For each plant, the three youngest leaves were collected at 48 hpi, immediately frozen in liquid nitrogen and stored at −80°C until RNA extraction.

Microscopic investigation of fungal development

To identify the optimal timing of sampling for the subsequent cDNA-AFLP experiment, progress of pathogen development was followed on ‘Président Roulin’ and ‘Gala’ using light microscopy. Plants were spray-inoculated and incubated as described above. At 1, 4, 6, 16, 20, 24, 32, 48, 54 and 120 hpi, the youngest expanding leaf of each cultivar was sampled (one leaf/sampling time/cultivar), cleared overnight in 99% methanol and stained with periodic acid-basic fuchsine according to the method by Preece [97]. Samples were thoroughly rinsed with water and mounted onto glass slides. Pathogen development stages were examined at the different time points by observing at least 200 conidia per sampling time for each of the cultivars.

In silicocDNA-AFLP simulations

A total of 450 full-length apple cDNAs sequences from the study by Newcomb et al. [34] were analyzed with the AFLPinSilico program [98]. The combinations of restriction enzymes used here were Sau3A, TaqI, DdeI, EcoRI, EcoRII, and ApoI as first sites in combination with the restriction enzyme MseI as second site, and vice-versa. The following parameters were estimated: percentage of cleaved cDNA, percentage of cDNA visualized on a gel, number of cleavage sites per cDNA, average distance between the last recognition site and the polyadenylation site, and the mean size of the restriction fragments.

RNA extraction and cDNA-AFLP

Total RNA was isolated from 100 mg of leaf material collected at 48 hpi (two plants per treatment), using the extraction method described by Gasic et al. [99]. After DNase I treatment (Fermentas Inc), purification of mRNA was performed starting from 250 μg total RNA using the Qiagen Oligotex mRNA kit (Qiagen Inc.). RNA purity and concentration was measured with the Nanodrop technology (Thermo Scientific Inc.). Double stranded cDNA was finally obtained starting from 500 ng mRNA following the instructions of the Superscript Double Stranded cDNA Synthesis kit (Invitrogen Inc.).

cDNA-AFLP analysis was performed with the AFLP Core Reagent kit (Invitrogen Inc.) as recommended by the manufacturer. The double-stranded cDNA was digested with EcoRI and MseI and ligated to the corresponding EcoRI and MseI adapters. The pre-amplification step was carried out using 20 cycles of amplification (94°C for 30 s; 56°C for 1 min; 72°C for 1 min) starting from a 5 μl aliquot of a 1:2 dilution of the ligation reaction and 75 ng of primer corresponding to the MseI and EcoRI adapter sequence without any extension, in 50 μl total volume. After 10-fold dilution of the PCR fragments, specific amplifications were carried out with a total of 141 EcoRI and MseI primer combinations containing two (123 primer combinations with EcoRI + 2/MseI + 2), or three additional selective bases at the 3’ end (18 primer combinations with EcoRI + 2/MseI + 3, EcoRI + 3/MseI + 3 or EcoRI + 3/MseI + 1). The Eco primers were labeled with [γ33P] dATP. Amplification products were separated by electrophoresis at 60 W on a vertical denaturing polyacrylamide gel (5%) containing 7 M urea for 3 hours 30 minutes. Gels were transferred onto Whatman 3MM paper before drying. Repeatability of the technique was checked trough an experimental replication of the selective amplification step (for a few selective primer pairs) starting from the same pre-amplified cDNA samples.

Bands corresponding to the TDFs were visualized on the polyacrylamide gel by autoradiography using X-ray films. Band intensities were digitized using the PhosphorImager tool (Biorad) and were quantified using QuantityOne software (Biorad). For each primer pair combination tested, only cDNA-AFLP profiles with the same global band intensity among genotypes and treatments were compared. This is presumed to reflect the fact that equivalent amounts of amplified cDNA were compared. We considered all the TDFs whose expression ratio between inoculated and non-inoculated control treatments was above the threshold of two (fold-change >2) as significantly up-regulated, and all those whose ratio was below the inverse threshold (fold-change <1/2) as significantly down-regulated.

Bands of interest belonging to the cultivar ‘Président Roulin’ were cut from the polyacrylamide gels, eluted in 50 μl of water, re-amplified with the selective primers and cloned into Pjet 1.2 vector (Fermentas). Two transformed colonies per TDF with an insert of the expected size were subsequently sequenced (Macrogen Inc.). cDNA sequences were manually checked for overall quality, and vector and primer contaminations were trimmed.

Bioinformatics analysis

In order to limit redundancy and to produce a longer consensus sequence, the sequences were clustered using Egassembler [36, 37]. We then compared the contigs and singletons identified with the whole apple genome sequence assembly v1.0 [38] by BLAST sequence similarity search. This second step allowed us to identify the TDFs that originated from the same original contig (and thus from the same gene) of the apple genome assembly. The functional annotation of this non-redundant database was performed using the automatic Blast2GO bioinformatics tool [100]. Basically, input sequences were queried by BLAST-X against the GenBank non-redundant sequences and ESTs database at the National Center for Biotechnology Information (NCBI) database, taking similarities with an E-value < 10−3 as significant matches. Then, the program extracts the GO terms associated to each of the obtained hits and returns an evaluated GO annotation for the query sequences (E-value < 1e−6). These E-values were the default values proposed by Blast2GO.

In order to detect GO annotations whose abundance was significantly different between different sets of annotated genes, a gene function enrichment analysis was performed using Blast2GO [100]. This software employs Fisher’s exact test to estimate the significance of associations between two categorical variables using a single test p-value. A set of GO terms that are under- or over-represented at a specified significance value (P < 0.05) were obtained as a result of performing the enrichment analysis. For this analysis, two subsets of our cDNA-AFLP library were compared: the TDFs being up- and down-regulated in ‘Président Roulin’ after pathogen infection (group II + and II -). We also compared our cDNA-AFLP library (group I and group II) from infected leaves of ‘Président Roulin’ to: (1) two EST libraries from uninfected young expanding leaves of Malus × domestica ‘Royal Gala’ (AELA) [34] and ‘Wijcik’ (Mdstw) [35], and (2) a cDNA library [25] of differentially expressed transcripts from Rvi6 (HcrVf2)-resistant transgenic ‘Gala’ lines challenged by V. inaequalis (obtained by suppression subtractive hybridization between Rvi6 (HcrVf2)-transformed ‘Gala’ lines and susceptible ‘Gala’ plants). The last three libraries were assembled in contigs and singletons using the EGassembler tool [37] before performing the analysis. Sequences were annotated with the bioinformatics tool Blast2GO using the same parameters as those applied for the annotation of the cDNA-AFLP fragments.

Quantitative real-time reverse transcription PCR (qRT-PCR)

qRT-PCR was carried out on the same RNA samples used for the cDNA-AFLP analysis and on RNA derived from one independent biological replication of the experiment. Specific TDF primers were designed with the software Primer3 [101] and qRT-PCR was performed using Biorad CFX96 and Maxima SYBR Green qPCR master mix (Fermentas Inc.), following the instruction of the manufacturer. A list of the specific primer pairs used for each TDF and product lengths is given in Additional file 3. PCR conditions were the same for all primer pairs: initial denaturation at 95°C for 10’ followed by 40 cycles of denaturation at 95°C for 15”, annealing at 60°C for 30” and extension at 72°C for 30”. Fluorescence data were collected at the end of the annealing step. Following cycling, samples were denatured at 95°C for 10”. The melting curve was analyzed to differentiate between the desired amplicons and any primer dimers or DNA contaminants (in the range 65°–95°C, with a temperature increment of 0.5°C for 5″). Each reaction was run in duplicate. The LinRegPCR software was used to confirm that individual PCR efficiencies were about 2 for each primer pair [102]. The relative expression ratio of the target genes between scab-inoculated and water-treated plants was evaluated using the ΔΔCt method described by Applied Biosystems (Relative expression ratio = 2ΔΔCt), with the glyceraldehyde 3-phosphate dehydrogenase gene (GAPDH) as the internal reference (primers sequence F-5′CAAGGTCATCCATGACAACTTTG3′, R-5′ GTCCACCACCCTGTTGCTGTAG3′). In fact, as it was the case in other qRT-PCR studies conducted on apple [33, 103], the GAPDH gene appeared to be the best housekeeping gene in our experimental conditions. Contrary to the elongation factor gene (EF, primers sequence F-5′TACTGGAACATCACAGGCTGAC3′, R-5′TGGACCTCTCAATCATGTTGTC3′), expression of the GAPDH was stable in scab-inoculated and water-treated leaf samples (Additional file 4). Individual relative expression values were then subjected to the ANOVA procedure, using Minitab 16 software, at a statistical significance level of P < 0.01.

In Silicomapping and co-localization of the TDFs with RGAs, QTL or apple scab major R genes

TDFs sequences were searched by BLAST-N within the whole genome sequence assembly v1.0 of apple [38] on the Genome Database of Rosaceae (GDR) [104] taking into account the best BLAST result (E-value < 0.001). Clusters of apple resistance genes analogues (RGAs) [90], QTLs and major scab resistance genes already identified in apple [3, 13, 14] were also anchored in silico to the apple genome sequence. Molecular markers flanking the major scab R genes and QTLs were obtained from the HIDRAS (High-quality Disease Resistant Apples for Sustainable Agriculture) website [105] and searched in the apple genome sequence by BLAST-N. Only SSR markers having an E-value ≤3e-3, a ratio of matched bases to marker sequence equal to 100% and a position on the expected chromosome on the ‘Golden Delicious genome sequence assembly (according to their genetic position) were anchored to the physical map. Clusters of RGAs and their physical position on the ‘Golden delicious’ genome assembly were retrieved from the publication of Perazzolli et al., [90]. Only TDFs mapping inside a QTL confidence interval or mapping within 250 kb from any cluster of RGAs or major scab R gene were considered to co-localize in the genomic regions involved in resistance. This distance has been used in previous publication as the average distance separating genes inside a cluster [90]. In order to determine if TDFs co-localized only by chance, frequencies of co-localization were compared to frequencies observed for ESTs from an uninfected apple library [34] using a χ2 test employing the statistical software Minitab 16.

Availability of supporting data

The whole cDNA-AFLP library isolated from ‘Président Roulin’ has been deposited at DDBJ/EMBL/GenBank under the accession numbers JZ719314 to JZ719813. Other supporting data are included as Additional files 1, 2, 3 and 4.

References

  1. Jha G, Thakur K, Thakur P: The Venturia apple pathosystem: pathogenicity mechanisms and plant defense responses. J Biomed Biotechnol. 2009, 2009: 1-10.

    Google Scholar 

  2. Sauphanor B, Dirwimmer C: Ecophyto R&D vers des systèmes de culture économes en produits phytosanitaires. Analyse comparative de différents systèmes en arboriculture fruitière. Tome IV. Edited by: MEEDDAT-MAP-INRA. 2009, France

    Google Scholar 

  3. Bus VGM, Rikkerink EHA, Caffier V, Durel CE, Plummer KM: Revision of the nomenclature of the differential host-pathogen interactions of Venturia inaequalis and Malus. Annu Rev Phytopathol. 2011, 49: 391-413. 10.1146/annurev-phyto-072910-095339.

    CAS  PubMed  Google Scholar 

  4. Poland JA, Balint-Kurti PJ, Wisser RJ, Pratt RC, Nelson RJ: Shades of gray: the world of quantitative disease resistance. Trends Plant Sci. 2009, 14: 21-29. 10.1016/j.tplants.2008.10.006.

    CAS  PubMed  Google Scholar 

  5. Chevalier M, Lespinasse Y, Renaudin S: A microscopic study of different classes of symptoms by the Vf gene in apple for resistance to scab (Venturia inaequalis). Plant Pathol. 1991, 40: 249-256. 10.1111/j.1365-3059.1991.tb02374.x.

    Google Scholar 

  6. Parisi L, Orts R, Rivenez-Damboise MO, Lefeuvre M, Lagarde MP: Protection intégrée du verger de pommiers de l’an 2000. Tavelure et oïdium: variétés résistantes et lutte raisonnée. Arboriculture Fruitière. 1995, 486: 25-29.

    Google Scholar 

  7. Lê Van A, Gladieux P, Lemaire C, Cornille A, Giraud T, Durel CE, Caffier V, Le Cam B: Evolution of pathogenicity traits in the apple scab fungal pathogen in response to the domestication of its host. Evol Appl. 2012, 5: 694-704. 10.1111/j.1752-4571.2012.00246.x.

    PubMed Central  PubMed  Google Scholar 

  8. Parisi L, Lespinasse Y, Guillaumes J, Krüger J: A new race of Venturia inaequalis virulent to apples with resistance due to the Vf gene. Phytopathology. 1993, 83: 533-537. 10.1094/Phyto-83-533.

    Google Scholar 

  9. Lateur M, Wagemans C, Populer C: Evaluation of fruit tree genetic resources as sources of polygenic scab resistance in an apple breeding programme. Acta Hortic. 1998, 484: 35-42.

    Google Scholar 

  10. Durel CE, Parisi L, Laurens F, Van de Weg WE, Liebhard R, Jourjon MF: Genetic dissection of partial resistance to race 6 of Venturia inaequalis in apple. Genome. 2003, 46: 224-234. 10.1139/g02-127.

    CAS  PubMed  Google Scholar 

  11. Liebhard R, Koller B, Patocchi A, Kellerhals M, Pfammatter W, Jermini M, Gessler C: Mapping quantitative field resistance against apple scab in a ‘Fiesta’ x ‘Discovery’ progeny. Phytopathology. 2003, 93: 493-501. 10.1094/PHYTO.2003.93.4.493.

    CAS  PubMed  Google Scholar 

  12. Calenge F, Faure A, Goerre M, Gebhardt C, Van de Weg WE, Parisi L, Durel CE: Quantitative Trait Loci (QTL) analysis reveals both broad-spectrum and isolate-specific QTL for scab resistance in an apple progeny challenged with eight isolates of Venturia inaequalis. Phytopathology. 2004, 94: 370-379. 10.1094/PHYTO.2004.94.4.370.

    CAS  PubMed  Google Scholar 

  13. Soufflet-Freslon V, Gianfranceschi L, Patocchi A, Durel CE: Inheritance studies of apple scab resistance and identification of Rvi14, a new major gene that acts together with other broad-spectrum QTL. Genome. 2008, 51: 657-667. 10.1139/G08-046.

    CAS  PubMed  Google Scholar 

  14. Gessler C, Patocchi A, Sansavini S, Tartarini S, Gianfranceschi L: Venturia inaequalis resistance in apple. Crit Rev Plant Sci. 2006, 25: 473-503. 10.1080/07352680601015975.

    CAS  Google Scholar 

  15. Vinatzer BA, Patocchi A, Gianfranceschi L, Tartarini S, Zhang HB, Gessler C, Sansavini S: Apple contains receptor-like genes homologous to the Cladiosporium fulvum resistance gene family of tomato with a cluster of genes cosegregating with Vf apple scab resistance. Mol Plant Microbe In. 2001, 14: 508-515. 10.1094/MPMI.2001.14.4.508.

    CAS  Google Scholar 

  16. Galli P, Patocchi A, Broggini GA, Gessler C: The Rvi15 (Vr2) apple scab resistance locus contains three TIR-NBS-LRR genes. Mol Plant Microbe In. 2010, 23: 608-617. 10.1094/MPMI-23-5-0608.

    CAS  Google Scholar 

  17. Belfanti E, Silfverberg-Dilworth E, Tartarini S, Patocchi A, Barbieri M, Zhu J, Vinatzer BA, Gianfranceschi L, Gessler C, Sansavini S: The HcrVf2 gene from a wild apple confers scab resistance to a transgenic cultivated variety. Proc Natl Acad Sci USA. 2004, 101: 886-890. 10.1073/pnas.0304808101.

    CAS  PubMed Central  PubMed  Google Scholar 

  18. Joshi SG, Schaart JG, Groenwold R, Jacobsen E, Schouten HJ, Krens FA: Functional analysis and expression profiling of HcrVf1 and HcrVf2 for development of scab resistant cisgenic and intragenic apples. Plant Mol Biol. 2011, 75: 579-591. 10.1007/s11103-011-9749-1.

    CAS  PubMed Central  PubMed  Google Scholar 

  19. Schouten HJ, Brinkhuis J, van der Burgh A, Schaart JG, Groenwold R, Broggini GA, Gessler C: Cloning and functional characterization of the Rvi15 (Vr2) gene for apple scab resistance. Tree Genet Genomes. 2014, 10: 251-260. 10.1007/s11295-013-0678-9.

    Google Scholar 

  20. Malnoy M, Xu M, Borejsza-Wysocka E, Korban SS, Aldwinckle HS: Two receptor-like genes, Vfa1 and Vfa2, confer resistance to the fungal pathogen Venturia inaequalis inciting apple scab disease. Mol Plant Microbe In. 2008, 21: 448-458. 10.1094/MPMI-21-4-0448.

    CAS  Google Scholar 

  21. Komjanc M, Festi S, Rizzotti L, Cervone F, De Lorenzo G: A leucine-rich repeat receptor-like protein kinase (LRPKm1) gene is induced in Malus x domestica by Venturia inaequalis infection and salicylic acid treatment. Plant Mol Biol. 1999, 40: 945-957. 10.1023/A:1006275924882.

    CAS  PubMed  Google Scholar 

  22. Cova V, Paris R, Passerotti S, Zini E, Gessler C, Pertot I, Loi N, Musetti R, Komjanc M: Mapping and functional analysis of four apple receptor-like protein kinases related to LRPKm1 in HcrVf2-transgenic and wild-type apple plants. Tree Genet Genomes. 2010, 6: 389-403. 10.1007/s11295-009-0257-2.

    Google Scholar 

  23. Gau AE, Koutb M, Piotrowski M, Kloppstech K: Accumulation of pathogenesis-related proteins in the apoplast of a susceptible cultivar of apple (Malus domestica cv. Elstar) after infection by Venturia inaequalis and constitutive expression of PR genes in the resistant cultivar Remo. Eur J Plant Pathol. 2004, 110: 703-711.

    CAS  Google Scholar 

  24. Degenhardt J, Al-Masri NA, Kurkcuoglu S, Szankowski I, Gau AE: Characterization by suppression subtractive hybridization of transcripts that are differentially expressed in leaves of apple scab resistant and susceptible cultivars of Malus domestica. Mol Genet Genomics. 2005, 273: 326-335. 10.1007/s00438-005-1136-7.

    CAS  PubMed  Google Scholar 

  25. Paris R, Cova V, Pagliarani G, Tartarini S, Komjanc M, Sansavini S: Expression profiling in HcrVf2-transformed apple plants in response to Venturia inaequalis. Tree Genet Genomes. 2009, 5: 81-91. 10.1007/s11295-008-0177-6.

    Google Scholar 

  26. Paris R, Dondini L, Zannini G, Bastia D, Marasco E, Gualdi V, Rizzi V, Piffanelli P, Mantovani V, Tartarini S: dHPLC efficiency for semi-automated cDNA-AFLP analyses and fragment collection in the apple scab-resistance gene model. Planta. 2012, 235: 1065-1080. 10.1007/s00425-012-1589-y.

    CAS  PubMed  Google Scholar 

  27. Vanderplank JE: Horizontal and vertical resistance. Disease Resistance in Plants. Edited by: Vanderplank JE. 1984, Orlando, Florida: Academic Press, 57-81.

    Google Scholar 

  28. Lateur M: Fruit tree genetic resources and Integrated Fruit Production. Acta Hortic. 2000, 525: 317-323.

    Google Scholar 

  29. Bachem CWB, Oomen RJFJ, Visser RGF: Transcript imaging with cDNA-AFLP: a step-by-step protocol. Plant Mol Biol Rep. 1998, 16: 157-173. 10.1023/A:1007468801806.

    CAS  Google Scholar 

  30. Qubbaj T, Reineke A, Zebitz CPW: Molecular interactions between rosy apple aphids, Dysaphis plantaginea, and resistant and susceptible cultivars of its primary host Malus domestica. Entomol Exp Appl. 2004, 115: 145-152.

    Google Scholar 

  31. Massart S, Jijakli MH: Identification of differentially expressed genes by cDNA-amplified fragment length polymorphism in the biocontrol agent Pichia anomala (Strain Kh5). Biol control. 2006, 96: 80-86.

    CAS  Google Scholar 

  32. Yao YX, Li M, Liu Z, Hao YJ, Zhai H: A novel gene, screened by cDNA-AFLP approach, contributes to lowering the acidity of fruit in apple. Plant Physiol Bioch. 2007, 45: 139-145. 10.1016/j.plaphy.2007.01.010.

    CAS  Google Scholar 

  33. Aldaghi M, Bertaccini A, Lepoivre P: cDNA-AFLP analysis of gene expression changes in apple trees induced by phytoplasma infection during compatible interaction. Eur J Plant Pathol. 2012, 134: 117-130. 10.1007/s10658-012-9970-z.

    CAS  Google Scholar 

  34. Newcomb RD, Crowhurst RN, Gleave AP, Rikkerink EHA, Allan AC, Beuning LL, Bowen JH, Gera E, Jamieson KR, Janssen BJ, Laing WA, McArtney S, Nain B, Ross GS, Snowden KC, Souleyre EJF, Walton EF, Yauk YK: Analyses of expressed sequence tags from apple. Plant Physiol. 2006, 141: 147-166. 10.1104/pp.105.076208.

    PubMed Central  PubMed  Google Scholar 

  35. Gasic K, Gonzalez DO, Thimmapuram J, Liu L, Malnoy M, Gong G, Han Y, Vodkin LO, Aldwinckle HS, Carroll NJ, Orvis KS, Goldsbrough P, Clifton S, Pape D, Fulton L, Martin J, Theising B, Wisniewski ME, Fazio G, Feltus FA, Korban SS: Comparative analysis and functional annotation of a large expressed sequence tag collection of apple. Plant Gen. 2009, 2: 23-38. 10.3835/plantgenome2008.11.0014.

    CAS  Google Scholar 

  36. The EG assembler website. http://egassembler.hgc.jp/,

  37. Masoudi-Nejad A, Tonomura K, Kawashima S, Moriya Y, Suzuki M, Itoh M, Kanehisa M, Endo T, Goto S: EGassembler: online bioinformatics service for large-scale processing, clustering and assembling ESTs and genomic DNA fragments. Nucleic Acids Res. 2006, 34: 459-462.

    Google Scholar 

  38. Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, Kalyanaraman A, Fontana P, Bhatnagar SK, Troggio M, Pruss D, Salvi S, Pindo M, Baldi P, Castelletti S, Cavaiuolo M, Coppola G, Costa F, Cova V, Dal Ri A, Goremykin V, Komjanc M, Longhi S, Magnago P, Malacarne G, Malnoy M, Micheletti D, Moretto M, Perazzolli M, Si-Ammour A, Vezzulli S, et al: The genome of the domesticated apple (Malus x domestica Borkh). Nat Genet. 2010, 42: 833-839. 10.1038/ng.654.

    CAS  PubMed  Google Scholar 

  39. Mathioudakis MM, Maliogka VI, Katsiani AT, Katis NI: Incidence and molecular variability of apple stem pitting and apple chlorotic leaf spot viruses in apple and pear orchards in Greece. J Plant Pathol. 2010, 92: 139-147.

    CAS  Google Scholar 

  40. Rausher MD: Co-evolution and plant resistance to natural enemies. Nature. 2001, 411: 857-864. 10.1038/35081193.

    CAS  PubMed  Google Scholar 

  41. Durrant WE, Rowland O, Piedras P, Hammond-Kosack KE, Jones JDG: cDNA-AFLP Reveals a striking overlap in race-specific resistance and wound response gene expression profiles. The Plant Cell. 2000, 12: 963-977. 10.1105/tpc.12.6.963.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. Gabriëls SHEJ, Takken FL, Vossen JH, de Jong CF, Liu Q, Turk SCHJ, Wachowski LK, Peters J, Witsenboer HMA, de Wit PJGM, Joosten MHJ: cDNA-AFLP combined with functional analysis reveals novel genes involved in the hypersensitive response. Mol Plant Microbe In. 2006, 19: 567-576. 10.1094/MPMI-19-0567.

    Google Scholar 

  43. Wang X, Liu W, Chen X, Tang C, Dong Y, Ma J, Huang X, Wei G, Han Q, Huang L, Kang Z: Differential gene expression in incompatible interaction between wheat and stripe rust fungus revealed by cDNA-AFLP and comparison to compatible interaction. BMC Plant Biol. 2010, 10: 9-10.1186/1471-2229-10-9.

    PubMed Central  PubMed  Google Scholar 

  44. Shi C, Chaudhary S, Yu K, Park SJ, Navabi A, McClean PE: Identification of candidate genes associated with CBB resistance in common bean HR45 (Phaseolus vulgaris L.) using cDNA-AFLP. Mol Biol Rep. 2011, 38: 75-81. 10.1007/s11033-010-0079-1.

    CAS  PubMed  Google Scholar 

  45. Li C, Faino L, Dong L, Fan J, Kiss L, De Giovanni C, Lebeda A, Scott J, Matsuda Y, Toyoda H, Lindhout P, Visser RGF, Bonnema G, Bai Y: Characterization of polygenic resistance to powdery mildew in tomato at cytological, biochemical and gene expression level. Mol Plant Pathol. 2012, 13: 148-159. 10.1111/j.1364-3703.2011.00737.x.

    CAS  PubMed  Google Scholar 

  46. Alignan M, Hewezi T, Petitprez M, Dechamp-Guillaume G, Gentzbittel L: A cDNA microarray approach to decipher sunflower (Helianthus annuus) responses to the necrotrophic fungus Phoma macdonaldii. New Phytol. 2006, 170: 523-536. 10.1111/j.1469-8137.2006.01696.x.

    CAS  PubMed  Google Scholar 

  47. Valsangiacomo C, Gessler C: Role of the cuticular membrane in ontogenic and Vf resistance of apple leaves against Venturia inaequalis. Phytopathology. 1988, 78: 1066-1069. 10.1094/Phyto-78-1066.

    Google Scholar 

  48. Silfverberg-Dilworth E, Patocchi A, Gessler C: Evaluation of in vitro grown apple shoot sensitivity to Venturia inaequalis using a detached leaf assay. IOBC-WPRS Bulletin. 2006, 29: 67-74.

    Google Scholar 

  49. Clark TA, Zeyen RJ, Smith AG, Bushnell WR, Szabo LJ, Vance CP: Host response gene transcript accumulation in relation to visible cytological events during Erysiphe graminis attack in isogenic barley lines differing at the Ml-a locus. Physiol Mol Plant P. 1993, 43: 283-298. 10.1006/pmpp.1993.1058.

    CAS  Google Scholar 

  50. Ruiz-Lozano JM, Roussel H, Gianinazzi S, Gianinazzi-Pearson V: Defense genes are differentially induced by a mycorrhizal fungus and Rhizobium sp. in wild-type and symbiosis-defective pea genotypes. Mol Plant Microbe In. 1999, 12: 976-984. 10.1094/MPMI.1999.12.11.976.

    CAS  Google Scholar 

  51. Lo SCC, Hipskind JD, Nicholson RL: cDNA cloning of a sorghum pathogenesis-related protein (PR-10) and differential expression of defense-related genes following inoculation with Cochliobolus heterostrophus or Colletotrichum sublineolum. Mol Plant Microbe In. 1999, 12: 479-489. 10.1094/MPMI.1999.12.6.479.

    CAS  Google Scholar 

  52. Kleemann J, Rincon-Rivera LJ, Takahara H, Neumann U, Ver Loren Van Themaat E, van der Does HC, Hacquard S, Stüber K, Will I, Schmalenbach W, Schmelzer E, O’Connell RJ: Sequential delivery of host-induced virulence effectors by appressoria and intracellular hyphae of the phytopathogen Colletotrichum higginsianum. PLOS Pathogens. 2012, 8: 1-15.

    Google Scholar 

  53. Chisholm S, Coaker G, Day B, Staskawicz BJ: Host-microbe interactions: shaping the evolution of the plant immune response. Cell. 2006, 124: 803-814. 10.1016/j.cell.2006.02.008.

    CAS  PubMed  Google Scholar 

  54. Dal Cin V, Barbaro E, Danesin M, Murayama H, Velasco R, Ramina A: Fruitlet abscission: a cDNA-AFLP approach to study genes differentially expressed during shedding of immature fruits reveals the involvement of a putative auxin hydrogen symporter in apple (Malus domestica L. Borkh). Gene. 2009, 442: 26-36. 10.1016/j.gene.2009.04.009.

    CAS  PubMed  Google Scholar 

  55. Baldo A, Norelli JL, Farrell RE, Bassett CL, Aldwinckle HS, Malnoy M: Identification of genes differentially expressed during interaction of resistant and susceptible apple cultivars (Malus X domestica) with Erwinia amylovora. BMC Plant Biol. 2010, 10: 1-10. 10.1186/1471-2229-10-1.

    PubMed Central  PubMed  Google Scholar 

  56. Breyne P, Zabeau M: Genome-wide expression analysis of plant cell cycle modulated genes. Curr Opin Plant Biol. 2001, 4: 136-142. 10.1016/S1369-5266(00)00149-7.

    CAS  PubMed  Google Scholar 

  57. Breyne P, Dreesen R, Cannoot B, Rombaut D, Vandepoele K, Rombauts S, Vanderhaeghen R, Inzé D, Zabeau M: Quantitative cDNA-AFLP analysis for genome-wide expression studies. Mol Genet Genomics. 2003, 269: 173-179.

    CAS  PubMed  Google Scholar 

  58. Maleck K, Levine A, Eulgem T, Morgan A, Schmid J, Lawton KA, Dangl JL, Dietrich RA: The transcriptome of Arabidopsis thaliana during systemic acquired resistance. Nature. 2000, 26: 403-410.

    CAS  Google Scholar 

  59. Katagiri F: A global view of defense gene expression regulation – a highly interconnected signaling network. Curr Opin Plant Biol. 2004, 7: 506-511. 10.1016/j.pbi.2004.07.013.

    CAS  PubMed  Google Scholar 

  60. Steiner B, Kurz H, Lemmens M, Buerstmayr H: Differential gene expression of related wheat lines with contrasting levels of head blight resistance after Fusarium graminearum inoculation. Theor Appl Genet. 2008, 118: 753-764.

    PubMed Central  PubMed  Google Scholar 

  61. Eckey C, Korell M, Leib K, Biedenkopf D, Janses C, Langen G, Kogel K-H: Identification of powdery mildew-induced barley genes by cDNA-AFLP: functional assessment of an early expressed MAP kinase. Plant Mol Biol. 2004, 55: 1-15.

    CAS  PubMed  Google Scholar 

  62. Wang GL, Mackill DJ, Bonman JM, McCouch SR, Champoux MC, Nelson RJ: RFLP mapping of genes conferring complete and partial resistance to blast in a durably resistant rice cultivar. Genetics. 1994, 136: 1421-1434.

    CAS  PubMed Central  PubMed  Google Scholar 

  63. Gilroy EM, Hein I, van der Hoorn R, Boevink PC, Venter E, McLellan H, Kaffarnik F, Hrubikova K, Shaw J, Holeva M, Lòpez EC, Borras-Hidalgo O, Pritchard L, Loake GJ, Lacomme C, Birch PRJ: Involvement of cathepsin B in the plant disease resistance hypersensitive response. Plant J. 2007, 52: 1-13. 10.1111/j.1365-313X.2007.03226.x.

    CAS  PubMed  Google Scholar 

  64. Tao Y, Xie Z, Chen W, Glazebrook J, Chang HS, Han B, Zhu T, Zou G, Katagiri F: Quantitative nature of Arabidopsis responses during compatible and incompatible interactions with the bacterial pathogen Pseudomonas syringae. Plant Cell. 2003, 15: 317-330. 10.1105/tpc.007591.

    CAS  PubMed Central  PubMed  Google Scholar 

  65. McHale L, Tan X, Koehl P, Michelmore RW: Plant NBS-LRR proteins: adaptable guards. Genome Biology. 2006, 7: 212-

    PubMed Central  PubMed  Google Scholar 

  66. Asai T, Tena G, Plotnikova J, Willmann MR, Chiu WL, Gomez-Gomez L, Boller T, Ausubel FM, Sheen J: MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002, 415: 977-983. 10.1038/415977a.

    CAS  PubMed  Google Scholar 

  67. Marrs KA: The functions and regulation of glutathione-S-transferase in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 127-158. 10.1146/annurev.arplant.47.1.127.

    CAS  PubMed  Google Scholar 

  68. Bindschedler LV, Dewdney J, Blee KA, Stone JM, Asai T, Plotnikov J, Denoux C, Hayes T, Gerrish C, Davies DR, Ausubel FM, Bolwell GP: Peroxidase-dependent apoplastic oxidative burst in Arabidopsis required for pathogen resistance. Plant J. 2006, 47: 851-863. 10.1111/j.1365-313X.2006.02837.x.

    CAS  PubMed Central  PubMed  Google Scholar 

  69. Yang CW, Gonzalez-Lamothe R, Ewan RA, Rowland O, Yoshioka H, Shenton M, Ye H, O’Donnell E, Jones JDG, Sadanandom A: The E3 ubiquitin ligase activity of Arabidopsis PLANT U-BOX17 and its functional tobacco homolog ACRE276 are required for cell death and defense. Plant Cell. 2006, 18: 1084-1098. 10.1105/tpc.105.039198.

    CAS  PubMed Central  PubMed  Google Scholar 

  70. Solomon M, Belenghia B, Delledonne M, Menachema E, Levine A: The involvement of cysteine proteases and protease inhibitor genes in the regulation of programmed cell death in plants. Plant Cell. 1999, 11: 431-443. 10.1105/tpc.11.3.431.

    CAS  PubMed Central  PubMed  Google Scholar 

  71. Bolton MD: Primary metabolism and plant defense: fuel for the fire. Mol Plant-Microbe In. 2009, 22: 487-497. 10.1094/MPMI-22-5-0487.

    CAS  Google Scholar 

  72. Rushton PJ, Somssich IE: Transcriptional control of plant genes responsive to pathogens. Curr Opin Plant Biol. 1998, 1: 311-315. 10.1016/1369-5266(88)80052-9.

    CAS  PubMed  Google Scholar 

  73. Lee J, Miura K, Bressan RA, Hasegawa PM, Yun DJ: Regulation of plant innate immunity by SUMO E3 ligase. Plant Signal Behav. 2007, 2: 253-254. 10.4161/psb.2.4.3867.

    PubMed Central  PubMed  Google Scholar 

  74. Sattler SE, Mene-Saffrane L, Farmer EE, Krischke M, Mueller MJ, DellaPenna D: Nonenzymatic lipid peroxidation reprograms gene expression and activates defense markers in Arabidopsis tocopherol-deficient mutants. Plant Cell. 2006, 18: 3706-3720. 10.1105/tpc.106.044065.

    CAS  PubMed Central  PubMed  Google Scholar 

  75. Maeda H, DellaPenna D: Tocopherol functions in photosynthetic organisms. Curr Opin Plant Biol. 2007, 10: 260-265. 10.1016/j.pbi.2007.04.006.

    CAS  PubMed  Google Scholar 

  76. Avrova AO, Stewart HE, De Jong W, Heilbronn J, Lyon GD, Birch PRJ: A cysteine protease gene is expressed early in resistant potato interactions with Phytophthora infestans. Mol Plant-Microbe In. 1999, 12: 1114-1119. 10.1094/MPMI.1999.12.12.1114.

    CAS  Google Scholar 

  77. Hao L, Hsiang T, Goodwin PH: Role of two cysteine proteinases in the susceptible response of Nicotiana benthamiana to Colletotrichum destructivum and the hypersensitive response to Pseudomonas syringae pv. tomato. Plant Sci. 2006, 170: 1001-1009. 10.1016/j.plantsci.2006.01.011.

    CAS  Google Scholar 

  78. Jansen MAK, Gaba V, Greenberg BM: Higher plants and UV-B radiation: balancing damage, repair and acclimatation. Trends Plant Sci. 1998, 3: 131-135. 10.1016/S1360-1385(98)01215-1.

    Google Scholar 

  79. Xiong L, Zhu JK: Molecular and genetic aspects of plant responses to osmotic stress. Plant Cell Environ. 2002, 25: 131-139. 10.1046/j.1365-3040.2002.00782.x.

    CAS  PubMed  Google Scholar 

  80. Andaya CB, Ronald PC: A catalytically impaired mutant of the rice Xa21 receptor kinase confers partial resistance to Xanthomonas oryzae pv oryzae. Physiol Mol Plant P. 2003, 62: 203-208. 10.1016/S0885-5765(03)00060-2.

    CAS  Google Scholar 

  81. Durel CE, Parisi L, Laurens F, Venisse JS, Jourjon MF: Does the Vf gene maintain a residual resistance to apple scab despite its breakdown by Venturia inaequalis race 6 strains. Acta Hortic. 2000, 538: 575-580.

    CAS  Google Scholar 

  82. Li ZK, Luo LJ, Mei HW, Paterson AH, Zhao XH, Zhong DB, Wang YP, Yu XQ, Zhu L, Tabien R, Stansel JW, Ying CS: A “defeated” rice resistance gene acts as a QTL against a virulent strain of Xanthomonas oryzae pv. oryzae. Mol Genet Genomics. 1999, 261: 58-63. 10.1007/s004380050941.

    CAS  Google Scholar 

  83. Brodny U, Nelson RR, Gregory LV: Residual and interactive expression of “defeated” wheat stem rust resistance genes. Phytopathology. 1986, 76: 546-549. 10.1094/Phyto-76-546.

    Google Scholar 

  84. Nass HA, Pedersen WL, MacKenzie DR, Nelson RR: The residual effect of some defeated powdery mildew Erysiphe graminis f.sp. tritici resistance genes in isolines of winter wheat. Phytopathology. 1981, 71: 1315-1348.

    Google Scholar 

  85. Talukder ZI, Tharreau D, Price AH: Quantitative trait loci analysis suggests that partial resistance to rice blast is mostly determined by race-specific interactions. New Phytol. 2004, 162: 197-209. 10.1111/j.1469-8137.2004.01010.x.

    CAS  Google Scholar 

  86. Perchepied L, Dogimont C, Pitrat M: Strain-specific and recessive QTLs involved in the control of partial resistance to Fusarium oxysporum f. sp. melonis race 1.2 in a recombinant inbred line population of melon. Theor and Appl Genet. 2005, 111: 65-74. 10.1007/s00122-005-1991-y.

    CAS  Google Scholar 

  87. Gebhardt C, Valkonen JPT: Organization of genes controlling disease resistance in the potato genome. Annu Rev Phytopathol. 2001, 39: 79-102. 10.1146/annurev.phyto.39.1.79.

    CAS  PubMed  Google Scholar 

  88. Calenge F, Van der Linden CG, Van de Weg E, Schouten HJ, Van Arkel G, Denancé C, Durel CE: Resistance gene analogues identified through the NBS-profiling method map close to major genes and QTL for disease resistance in apple. Theor Appl Genet. 2005, 110: 660-668. 10.1007/s00122-004-1891-6.

    CAS  PubMed  Google Scholar 

  89. Xiao W, Zhao J, Fan S, Li L, Dai J, Xu M: Mapping of genome-wide resistance gene analogs (RGAs) in maize (Zea mays l.). Theor App Genet. 2007, 115: 501-508. 10.1007/s00122-007-0583-4.

    CAS  Google Scholar 

  90. Wang Z, Taramino G, Yang D, Liu G, Tingey SV, Miao GH, Wang GL: Rice ESTs with disease-resistance gene-or defense-response gene-like sequences mapped to regions containing major resistance genes or QTLs. Mol Genet Genomics. 2001, 265: 302-310. 10.1007/s004380000415.

    CAS  PubMed  Google Scholar 

  91. Chu Z, Ouyang Y, Zhang J, Yang H, Wang S: Genome-wide analysis of defense-responsive genes in bacterial blight resistance of rice mediated by the recessive R gene xa13. Mol Genet Genomics. 2004, 271: 111-120. 10.1007/s00438-003-0964-6.

    CAS  PubMed  Google Scholar 

  92. Shi C, Ingvardsen C, Thümmler F, Melchinger AE, Wenzel G, Lübberstedt T: Identification by suppression subtractive hybridization of genes that are differentially expressed between near-isogenic maize lines in association with sugarcane mosaic virus resistance. Mol Genet Genomics. 2005, 273: 450-461. 10.1007/s00438-004-1103-8.

    CAS  PubMed  Google Scholar 

  93. Jensen PJ, Fazio G, Altman N, Praul C, McNellis TW: Mapping in an apple (Malus x domestica) F1 segregating population based on physical clustering of differentially expressed genes. BMC Genomics. 2014, 15: 261-10.1186/1471-2164-15-261.

    PubMed Central  PubMed  Google Scholar 

  94. Perazzolli M, Malacarne G, Baldo A, Righetti L, Bailey A, Fontana P, Velasco R, Malnoy M: Characterization of resistance gene analogues (RGAs) in apple (Malus x domestica Borkh.) and their evolutionary history of the Rosaceae family. PLoS One. 2014, 9: e83844-10.1371/journal.pone.0083844.

    PubMed Central  PubMed  Google Scholar 

  95. Gilad Y, Rifkin SA, Pritchard JK: Revealing the architecture of gene regulation: the promise of eQTL studies. Trends Genet. 2008, 24: 408-415. 10.1016/j.tig.2008.06.001.

    CAS  PubMed Central  PubMed  Google Scholar 

  96. Laurens F, Chevalier M, Dolega E, Gennari F, Goerre M, Fischer C, Kellerhals M, Lateur M, Lefrancq B, Parisi L, Schouten HJ, Tartarini S: Local European cultivars as sources of durable scab resistance in apple. Acta Hort. 2004, 663: 115-122.

    Google Scholar 

  97. Preece TF: A staining method for the study of apple scab infections. Plant Pathol. 1959, 8: 127-129. 10.1111/j.1365-3059.1959.tb00898.x.

    Google Scholar 

  98. Rombauts S, van de Peer Y, Rouzé P: AFLPinSilico, simulating AFLP fingerprints. Bioinformatics. 2003, 19: 776-777. 10.1093/bioinformatics/btg090.

    CAS  PubMed  Google Scholar 

  99. Gasic K, Hernandez A, Korban SS: RNA extraction from different apple tissues rich in polyphenols and polysaccharides for cDNA library construction. Plant Mol Biol Rep. 2004, 22: 437a-437g. 10.1007/BF02772687.

    Google Scholar 

  100. Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.

    CAS  PubMed  Google Scholar 

  101. Rozen S, Skaletsky HJ: Primer3 on the www for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, New Jersey: Humana Press, 365-386.

    Google Scholar 

  102. Ramakers C, Ruijter JM, Lekanne-Deprez RH, Moorman AFM: Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett. 2003, 339: 62-66. 10.1016/S0304-3940(02)01423-4.

    CAS  PubMed  Google Scholar 

  103. Gadiou S, Kundu JK: Evaluation of reference genes for the relative quantification of apple stem grooving virus and apple mosaic virus in apple trees. Indian J Virol. 2012, 23: 39-41. 10.1007/s13337-012-0065-4.

    CAS  PubMed Central  PubMed  Google Scholar 

  104. Genome Database for Rosaceae. http://www.rosaceae.org/,

  105. HIDRAS (High-quality Disease Resistant Apples for Sustainable Agriculture). website: http://www.hidras.unimi.it/HiDRAS-SSRdb/pages/index.php

Download references

Acknowledgements

We thank Valérie Caffier from INRA, France, for supplying the V. inaequalis isolates. We are also grateful to Amy Watson and Vincent Bus from the New Zealand Institute for Plant and Food Research for helping us in the revision of the draft manuscript. Many thanks also to the ‘RGF’ team from CRA-W (particularly Laurent Delpierre, Thibaut Donis, Pascal Dupont and Alain Rondia) for their valuable help in the grafting and raising of the apple trees used in this experiment. This work was supported by the Moerman funds (CRA-W) within the framework of the POMINNO project entitled “ Recherche de méthodes rapides de sélection de variétés innovantes de pommes de qualités différenciées et adaptées à une agriculture durable”.

Author information

Authors and Affiliations

  1. Life Sciences Department, Breeding and Biodiversity Unit, Walloon Agricultural Research Center, Rue de Liroux, 4, 5030, Gembloux, Belgium

    Héloïse Bastiaanse, Yordan Muhovski & Marc Lateur

  2. Life Sciences Department, Bioengineering Unit, Walloon Agricultural Research Center, Chaussée de Charleroi, 234, 5030, Gembloux, Belgium

    Dominique Mingeot

  3. Consiglio per la Ricerca e la sperimentazione in Agricoltura, Centro di Ricerca per le Colture Industriali, CRA-CIN, via di Corticella 133, 40128, Bologna, Italy

    Roberta Paris

  4. Plant Pathology Unit, Gembloux Agro-Bio Tech, University of Liège, Passage des déportés 2, 5030, Gembloux, Belgium

    Héloïse Bastiaanse & Olivier Parisi

Authors
  1. Héloïse Bastiaanse
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yordan Muhovski
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Olivier Parisi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Roberta Paris
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Dominique Mingeot
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Marc Lateur
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Héloïse Bastiaanse.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HB, DM, YM, OP and ML conceptualized and designed the study. ML and DM coordinated the study. HB did the inoculation and the phenotyping of the plants with the advice of ML. OP helped in the setting up of the cDNA-AFLP protocol. HB carried out the cDNA-AFLP and qRT-PCR experiments under the supervision of YM and DM. HB annotated the sequences and analyzed the results. HB, DM, YM, RP, OP and ML participated in the results discussion. HB drafted the manuscript. DM, YM and RP reviewed the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2014_6851_MOESM1_ESM.xlsx

Additional file 1: Sequence annotation using the bioinformatics software Blast2GO. This Table provides a full list of differentially expressed TDFs isolated from ‘Président Roulin’, their expression pattern according to the cDNA-AFLP and their annotation using the bioinformatics software Blast2GO. (XLSX 101 KB)

12864_2014_6851_MOESM2_ESM.xlsx

Additional file 2: Over-representation of GO categories in our cDNA-AFLP library compared to Rvi6 (HcrVf2 )-’Gala’ transformed library [25]. This table provides a list of TDFs and their annotation belonging to Biological Process GO categories that were over-represented in the ‘Président Roulin’ cDNA-AFLP library (partially resistant) compared with a cDNA library from completely resistant Rvi6 (HcrVf2)-transformed ‘Gala’ lines challenged with V. inaequalis[25]. Gene enrichment analysis was conducted with the software Blast2Go using Fisher’s Exact Test at a p-value <0.05. (XLSX 13 KB)

12864_2014_6851_MOESM3_ESM.xlsx

Additional file 3: List of primer pairs used for qRT-PCR validation of differential gene expression revealed by cDNA-AFLP. This table provides the primer sequences used for the qRT-PCR analyses. (XLSX 13 KB)

12864_2014_6851_MOESM4_ESM.docx

Additional file 4: Expression variation of candidate housekeeping genes in apple leaves challenged by V. inaequalis and mock-inoculated. This figure shows that the GAPDH gene appeared to be the best housekeeping gene in our experimental conditions. Contrary to the EF gene, expression of the GAPDH was stable in scab-inoculated and water-treated leaf samples. (DOCX 29 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bastiaanse, H., Muhovski, Y., Parisi, O. et al. Gene expression profiling by cDNA-AFLP reveals potential candidate genes for partial resistance of ‘Président Roulin’ against Venturia inaequalis. BMC Genomics 15, 1043 (2014). https://doi.org/10.1186/1471-2164-15-1043

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-1043

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us