Introduction

Among the several soft rot enterobacterial necrotrophic pathogens of the two genera Pectobacterium and Dickeya, which cause blackleg in potato plants and soft rot of potato tubers, Dickeya solani has received the most attention due to highly aggressive infection caused by this pathogen in potato plants and the ability of this species to initiate disease occurrence at low inoculum levels (Czajkowski et al. 2009; Toth et al. 2011). Potato tubers can be infected by bacteria through wounds, natural openings (stomata, lenticels), or stolons. Seed tubers infected with bacteria are the main sources of inocula. Bacteria can move throughout the plant via the water-conducting tissues and can infect growing plants and progeny tubers through stolons. Infected seed tubers may also rot in the soil and reinfect growing potato plants (through roots) and progeny tubers (directly through wounds and lenticels or through stolons). Potato tubers can be contaminated with pectinolytic bacteria and infected through wounds during harvest, transport and grading. Elphinstone and Perombelon (1986) stated that a single rotting tuber contaminated c. 100 kg of potatoes during mechanical grading, c. 50% being contaminated with 104–105 bacteria per tuber. There is no chemical control that can be applied, and extreme resistance has not been observed in the genus Solanum. Mapping studies in a diploid population originating from hybrids of Solanum chacoense and Solanum yungasense have shown the complexity of the inheritance of resistance to Pectobacterium atrosepticum (Zimnoch-Guzowska et al. 2000). Broad studies on soft rot resistance have shown wide-ranging variations of partial resistance in potato cultivars and candidate breeding lines in France (Pasco et al. 2006). Tuber infection occurs in two phases: an asymptomatic phase, in which bacteria do not express genes encoding plant cell wall-degrading enzymes (Toth and Birch 2005), and a late phase, in which cell wall-degrading enzymes genes are activated by various means, such as quorum sensing, temperature, lack of nutrients and low oxygen conditions (Pérombelon 2002; Sepulchre et al. 2007; Davidsson et al. 2013). The products of pectinolytic enzymes activate plant defense (Sepulchre et al. 2007). In oxygen-rich conditions, a barrier-like mechanism of bacterial resistance is observed in tubers as a result of a rapid form of wound healing (Chung et al. 2017).

Proteomics is a powerful tool for the study of defense mechanisms associated with plant-pathogen interactions. Proteomics could be useful for determination of quantitative and qualitative changes in proteins involved in the defense response in plants and to understand the relationships between these proteins and biological resistance (Fang et al. 2015). Potato tuber extract increased the expression more than 40 proteins in the secretome of the soft rot enterobacterium P. atrosepticum, including known virulence factors, such as pectic enzymes and proteases (Mattinen et al. 2007). The first proteomic study of potato tissue compared the protein expression profiles of P. atrosepticum- and mock-inoculated tubers of two potato cultivars (Barzic and Com 2012). Thirteen proteins were upregulated 48 h post inoculation in the partially resistant potato cultivar Kerpondy, but not in the susceptible cultivar Bintje, in comparison to the mock-inoculated tubers.

In our study we report, for the first time, the detection of proteins associated with increased resistance of tubers to D. solani in tetraploid and diploid potato plants. Nano-LC-MS-MS/MS (nano-liquid chromatography coupled with tandem mass spectrometry) was used to identify differentially expressed proteins between the resistant and susceptible potato plants following D. solani invasion at an early symptomatic phase of infection. In addition, differentially expressed proteins between the resistant and susceptible potato plants in intact tubers and after mock inoculation are presented.

Materials and methods

Plant material

Four potato cultivars—Katahdin and three Katahdin-derived cultivars, namely, Bea, Humalda, and Ulster Supreme—were obtained from the potato collection at the Plant Breeding and Acclimatization Institute, National Research Institute, Bonin, Poland. The cultivar Irys, used as standard of high susceptibility of tubers to pectinolytic bacteria, scored 2 on a scale of 1 to 9, where 9 is the most resistant, according to the descriptions of potato cultivars registered in Poland. Two diploid clones, namely, the highly resistant DG 00–270 and the susceptible DG 08–305, were obtained from the Młochów collection of diploid potato plants. DG 00–270 was derived by intercrossing Solanum tuberosum L with Solanum phureja Juz. and Bukasov and other wild species: S. chacoense Bitter, S. yungasense Hawk, Solanum verrucosum Schltdl, and Solanum microdontum Bitter. In addition, the moderately resistant cultivar Glada, which scored 5 on a scale of 1 to 9, where 9 is the most resistant, was used to test the kinetics of symptom development under the applied experimental conditions.

Bacterial inoculum

The bacterial strain of D. solani, IFB0099 syn. IPO2276 (IPO: the collection of Plant Research International, Wageningen, The Netherlands), was used in this study. This strain was kindly provided by Prof. E. Lojkowska, from the collection of the Intercollegiate Faculty of Biotechnology University of Gdańsk and Medical University of Gdańsk, Poland. The bacterium (stored at −70 °C) was first grown on LB agar medium for 24 h at 27 °C and then resuspended in sterile deionized water, and the suspension was adjusted to an optical density (OD600) of 1.0, equivalent to 109 CFU ml−1, using a spectrophotometer (Hitachi U-1900, Tokyo, Japan).

Testing of potato cultivars/clones for tuber soft rot resistance

Tuber resistance to pectinolytic bacteria was tested by the method described by Lebecka et al. (2018). Tubers were stored in the dark at 8 °C for 3 months before the experiments were conducted. Each tuber surface-disinfected in 1% sodium hypochlorite solution for 15 min, rinsed in distilled water and air dried at room temperature 1 day before inoculation, was wounded with a steel rod (10 mm, Ø 2 mm) and inoculated with 10 μL of bacterial suspension, which was pipetted into a wound hole. The holes were sealed with Vaseline and a piece of a parafilm. The tubers were externally sprayed with water, enclosed in boxes and kept for 3 days at 26 °C. Six experiments were conducted in January of 3 years (2015–2017). In total, during the 3 years of testing, 58 tubers were tested for each potato cultivar, and 40 tubers were tested for diploid clones. The resistance of potato tubers can be expressed as the number of successfull infected tubers and the level of infection in infected tubers. Therefore two components of resistance were assessed: disease severity (DS), expressed as mass of macerated tissue from tubers with symptoms of infection, and disease incidence (DI), expressed as the average percentage of tubers with symptoms of rot. Statistical analyses were performed using Statistica 10 software (Statsoft Inc. 2011). To test the effect of the potato cultivar/clone, the year and interactions with DS values, two-way factorial ANOVA was used. To test the significance of differences among cultivars, Tukey’s test was used.

Determining the time of symptoms occurrence in wounded potato tubers inoculated with bacteria

To identify the early phase of symptom development following the inoculation of potato tubers with the bacterium D. solani IFB0099 for the applied experimental conditions, the kinetics of symptom development were analyzed. At one-hour intervals starting 5 h after inoculation, the rotten tissue was weighed for 6 tubers from each of the two potato cultivars, namely, Irys and Glada.

Inoculation of potato tubers for protein extraction

Tubers were stored in the dark at 8 °C for 6 months before the experiments were conducted. For proteomic analyses, tubers were inoculated as described above. Control tubers were kept intact or wounded and treated with sterilized water instead of the bacterial suspension. Two samples from each of two different tubers on each of two different experimental dates were collected 8 h post inoculation. In total, 8 samples for each of five potato cultivars for every three combinations (of bacteria-inoculated, mock-inoculated and intact tubers) were examined. Accordingly, 8 samples for each of two diploid clones for every two combinations, bacteria-inoculated and mock-inoculated tubers and 4 samples for intact-tubers, were examined.

Protein extraction

The total proteomes of the tuber tissues were extracted as described by Murawska et al. (2017). Tuber tissue fragments of up to 250 mg were cut out from the edges of inoculated or symptom-exhibiting tubers using a corkborer (Ø 4 mm). The samples were immediately frozen in liquid nitrogen and then kept in a freezer at −80 °C before freeze-drying for 2 days. Lyophilized samples were homogenized using micropestles, and proteins were solubilized in 135 μl of 25 mM ammonium bicarbonate and 15 μl of 4% SDS. After 15 min of sonication, 100 μl of 25 mM ammonium bicarbonate was added, and the samples were vortexed and then centrifuged for 10 min at 12000×g. The supernatant was collected, and the total protein content in the sample was measured using the bicinchoninic acid assay (Smith et al. 1985). The protein concentration was estimated twice: once after protein extraction from homogenized tissue, to minimize variations resulting from material isolation and tissue homogenization, and once after acetone precipitation, to minimize the effect of the resolubilization of different proteins among samples. For digestion of peptide bonds carboxyterminal to lysine under high urea concentrations, the total protein content (and hence the amount of enzyme) in the samples was adjusted to 50 μg with 20 μl of 100 mM ammonium bicarbonate containing 0.1 μg of endoproteinase Lys-C (Pierce™ Lys-C Protease, MS-grade, 90,307). After 90 min of incubation at room temperature, to disrupt disulfide bonds, 1 μl of 0.5 M Tris (2-carboxyethyl)phosphine hydrochloride was added, and the samples were incubated for 20 min at 60 °C. To prevent the regeneration of disulfide bridges, 3 μl of 0.2 MS-methyl methanethiosulfonate was added. Overnight tryptic hydrolysis of 50 μg of total protein in the samples was performed at 37 °C using 20 μl of sequencing-grade modified trypsin (Promega, V5117). Two microliters of 1% TFA (trifluoroacetic acid) was added to the samples, which were then stored in refrigerator until nano-LC-MS-MS/MS analysis.

Mass spectrometry

Peptide mixtures were analyzed by nano-LC-MS-MS/MS using a nanoACQUITY (Waters) LC system and a Q-Exactive mass spectrometer (Thermo Electron Corp., San Jose, CA). The peptide mixture was applied to an RP-18 precolumn (nanoACQUITY Symmetry® C18, Waters 186,003,514) using water containing 0.1% TFA as a mobile phase and then transferred to a nano-HPLC RP-18 column (nanoACQUITY BEH C18, Waters 186,003,545) using an acetonitrile gradient (5–35% AcN in 180 min) in the presence of 0.05% formic acid with a flow rate of 250 nl/min. The column outlet was directly coupled to the electrospray ion source of the spectrometer functioning under a data-dependent MS-to-MS/MS switch. A blank run ensuring lack of cross contamination from previous samples preceded each analysis.

The raw data acquired were processed by Mascot Distiller followed by a Mascot Search (Matrix Science, London, UK, on-site license) against the UniProt Solanum tuberosum database (June 2015 release) restricted to human sequences. The search parameters were as follows: precursor and product ion mass tolerance of 30 ppm and 0.1 Da, respectively; enzyme specificity, trypsin; missed cleavages allowed, 1; fixed modification of cysteine by a methylthio group; and variable modification by methionine oxidation. Peptides with Mascot scores exceeding the threshold value corresponding to <5% of the expected value, as calculated by Mascot, were considered to be positively identified.

The mass spectrometry data have been deposited to the ProteomeXchange Consortium via the PRIDE [1] partner repository with the dataset identifier PXD013009 and https://doi.org/10.6019/PXD013009.

Quantitative analysis was performed as described by Bakun et al. (2012). Briefly, the mass calibration and data filtering described above were carried out with MScan software, which was developed in-house (http://proteom.ibb.waw.pl/). Only proteins represented by at least two peptides were considered, and proteins identified by a subset of peptides from another protein were excluded. Peptides lists from the LC-MS/MS runs that matched the acceptance criteria were merged and overlaid onto 2-D heat maps generated from the LC-MS/MS datasets using m/z and retention time coordinates and validated by predicted and observed isotopic pattern correlation (acceptance criteria: RMS > 90%) using MSparky (http://proteom.ibb.waw.pl/). Peptide intensities were determined as the surfaces of the tagged isotopic envelopes. Finally, lists of identified peptides with corresponding quantitative values were exported for statistical analysis, which was carried out with the Diffprot software program (version 1.5.19; 3.01.2013) developed in-house (Malinowska et al. 2012). Prior to analysis, quantitative values were normalized with LOWESS. Proteins with more than 90% common peptides were clustered, and only peptides unique to the cluster were used for statistical analysis. Only proteins with q-values ≤0.1 were taken into consideration during further analysis.

Principal component analysis was performed for proteins analysed from tubers of the five potato cultivars 8 h post mock- and bacteria-inoculation. The features used for PCA were normalized and log-transformed peptide intensitie, with median value imputation of missing data. Peptides missing in more than 30% of samples were removed from the data set. The analysis itself was performed using Scikit-learn 0.19.2 (Pedregosa et al. 2011) and the results visualized with ggplot 2 (Wickham 2017).

Results

Soft rot resistance of tubers of potato cultivars and diploid hybrids of Solanum spp.

Resistance to tuber soft rot was evaluated under laboratory conditions using point inoculation of whole tubers. The DS data used for analysis included only data for tubers with symptoms of infection. Two-way ANOVA indicated the significant effects of potato cultivar/clone (P < 0.001) and year of testing (P = 0.0034), but not of the interaction of cultivar and year (P = 0.072), on DS values. The ratings for the responses of potato cultivars and diploid clones to soft rot bacteria are shown in Figs. 1 and 2. Post hoc analysis of the DS data has shown that two potato cultivars, namely, Humalda and Bea, and one diploid clone, namely, DG 00–270, were significantly more resistant than the three potato cultivars Katahdin, Ulster Supreme, and Irys and the diploid clone DG 08–305. For DI, the percentage of tubers with symptoms of rot ranged from 63 to 84% in a group of resistant cultivars/clones and from 95 to 100% in a group of susceptible cultivars/clones. The correlation coefficient between DS and DI for tested cultivars/clones was r = 0.729 with P = 0.06.

Fig. 1
figure 1

Disease severity (DS) and disease incidence (DI) of potato cultivars/diploid clones after inoculation of tubers with D. solani. DS is expressed as weight of rotten tissue [g], and DI is expressed as percentage (%) of tubers with symptoms of rot. Each value represents the mean number of successfully infected tubers out of 58 tubers for the cultivars and 40 tubers for the clones after incubation at 26 °C for 3 days. The tubers were assayed in six independent experiments over 3 years. Means followed by the same letter do not differ significantly according to Tukey’s range test (P < 0.05). Error bars indicate the standard error

Fig. 2
figure 2

Reaction of potato cultivars and diploid clones to inoculation of tubers with bacteria with D. solani and incubation at temperature 26 °C for 3 days. High resistance is expressed in the form of wound healing (the second tuber of cultivar Humalda and tubers of the diploid clone DG 00–270) and was used for calculation the disease incidence. The level of infection (disease severity) in infected tubers was expressed as weight of rotten tissue [g], which was removed from cut tubers

Determining the time of symptoms occurrence in wounded potato tubers inoculated with bacteria

Symptoms of rot in tubers inoculated with D. solani under the conditions applied for evaluation of tuber resistance to soft rot at one-hour intervals were not observed 7 h post inoculation. However, the symptoms were visible 8 h after inoculation, regardless of the level of resistance to bacteria in the two tested potato cultivars (Fig. S1).

Proteins associated with D. solani resistance

Eight total samples from four different tubers from each of the cultivars were tested. Increase in protein abundance in the resistant cultivars (16 samples) compared to the susceptible group (24 samples) were analyzed for each of three treatments. For each of the highly resistant (DG 00–270) and the susceptible (DG 08–305) diploid clones, from four to eight samples for each of three treatments were compared. On average, in a single experiment, after filtering out proteins with low Mascot scores or those identified by one peptide only, approximately 1450 (ranging from 1233 to 1733) proteins were selected for quantification. From this data, differentially expressed proteins with ratios ≥1.9 and q-values ≤0.1 between the resistant and susceptible samples were selected. There were 8 constitutively expressed proteins that exhibited significant fold changes between the cultivar groups following bacterial infection: six exhibited increased abundance - probable inactive patatin-03-Kuras 1 (PT3K1), proteinase inhibitor PTI (PI), proteinase inhibitor type 2 (PI2), chloroplastic thiamine thiazole synthase (THI1), 4 subunits of chymotrypsin inhibitor (CI1), 4 isoforms of patatin (PT) and two exhibited decreased abundance – 4 other isoforms of patatin and defensing-like protein. For PI2, PT3K1, CI1 and PT, significant changes were observed after mock inoculation. PT also exhibited differential abundance in intact tubers. The fold change values are shown in Table 1. When the proteins of each of the resistant cultivars were compared with those of the susceptible group individually, PT3K1 and PI exhibited significantly increased protein abundance after D. solani inoculation in both resistant cultivars (Fig. 3). PI2, THI1, CI1 and PT exhibited increased abundance in D. solani-infected Humalda tubers and not in the Bea cultivar (Table 2).

Table 1 List of proteins in the resistant group that significantly differed from those in the susceptible group
Fig. 3
figure 3

Box plots of relative abundance of proteins in the tested cultivars: a the proteinase inhibitor PTI, b the probable inactive patatin-3-Kuras 1. The box represents quartiles, and the band inside the box is the median. The whiskers extend to data extremes. No outlier values more than 1.5 times the interquartile range from the box were observed

Table 2 List of proteins that exhibited significant differences in expression between the resistant cultivars Bea or Humalda and the susceptible group

Principal component analysis of data obtained from potato tubers of the five cultivars, both bacteria- and mock-inoculated, showed that potato cultivars were not well separated from each other in the dimension of the first component (Fig. 4). However, the cultivars were well separated in the dimension of the second component. PC1 explained 11% and PC2–8% of total variance.

Fig. 4
figure 4

Principal component analysis investigates whether the quantified protein patterns cluster according to cultivar and inoculation treatment. The variation explained by PC1 is 12% and by PC2–8%. Each point corresponds to the protein pattern of each replicate sample, projected onto a two-dimensional principal component space. The analysis shows separation of resistant and susceptible potato cultivars along the second component axis (PC2) indicating differences in protein abundance patterns

In the diploid clones, 21 proteins with significant fold changes were identified between DG 00–270 and DG 08–305 after D. solani infection. However, only metallocarboxypeptidase inhibitor (MCPI) and MCPI-like proteins exhibited much higher fold changes following pathogenic invasion (274.4- and 368.6-fold, respectively) than after mock inoculation (165.5- and 130.7-fold, respectively) and in intact tubers (nonsignificant and 49.2-fold, respectively) (Table 3).

Table 3 List of proteins that exhibited significant differences in expression between the resistant diploid clone DG 00–270 and the susceptible DG 08–305

Discussion

For the proteomics studies, five potato cultivars and two diploid clones of Solanum were selected based on the results of six biological tests, conducted in three consecutive years, in which DI and DS were evaluated in wounded tubers inoculated with the bacterium D. solani. An aggressive strain of D. solani IFB0099 in studies of Golanowska et al. (2015) has been selected for inoculation of potato in tuber soft rot assays and in protein extraction experiments. The wound inoculation method for testing tuber resistance to pectinolytic bacteria gives similar results as immersion of mechanically harvested tubers in a bacterial inoculum (Tzeng et al. 1990). Two potato cultivars, namely, Humalda and Bea, and one diploid clone, namely, DG 00–270, exhibited significantly more resistance to tuber maceration than the three potato cultivars, namely, Katahdin, Ulster Supreme, and Irys, and the diploid clone DG 08–305 used in this study. The clone DG 00–270 exhibited the lowest DS (0.7 g) and DI (63%); however, the DS of this clone did not differ significantly from those of the two resistant potato cultivars (1.6 and 2.0 g) but differed considerably from that of the most susceptible cultivar, Irys (12.9 g).

Pectinolytic bacteria are characterized by high levels of parasitism. These bacteria can coexist with their host plants without causing damage. The duration of latent infection is dependent on many factors. We determined the time of symptom occurrence under our experimental conditions, and accordingly, the tuber tissue samples used for proteomic studies were collected 8 h post inoculation, which reflects an early symptomatic phase of infection. In this phase, bacteria produce an arsenal of plant cell wall-degrading enzymes and trigger the resistance response in potato plants (Sepulchre et al. 2007). Simultaneously, synthesis of lignin and the prosuberin lamella in potato tubers occur 4 h post wounding (Lapierre et al. 1996).

Comparison of the protein expression profiles of pathogen-inoculated plants and controls provides information on only the response of plants to pathogenic attack, not on the mechanisms by which plants resist pathogenic infection (Zhang et al. 2013). In our study, to identify proteins associated with resistance processes, the protein profiles of tubers from potato cultivars or diploid clones that differed in D. solani resistance were analyzed. The cultivated potato is a highly heterozygous species with complex genetics. To minimize the possibility of false positive selection of proteins underlying plant resistance to the pathogen, we used samples from two groups of potato cultivars that differed in D. solani resistance. Six proteins, namely, PI2, PT3K1, THI1, PTI, CI1 and PT, were upregulated in tubers from the resistant group compared to the susceptible group after D. solani inoculation (Table 1). Four of these proteins, namely, PI2, THI1, CI1 and PT, exhibited significantly higher fold change values in only the resistant cultivar Humalda and not in the resistant cultivar Bea (Table 2). This result confirms the complexity of the defensive response in potato tubers against bacterial invasion.

Two proteins, namely, PT3K1 and PI, were revealed to be upregulated in both resistant cultivars (Fig. 3). PT3K1 levels were increased to 2.3- and 20.9-fold in the tubers of the Bea and Humalda cultivars, respectively, upon inoculation with D. solani. The corresponding fold change values for PI were 3.2 and 4.4. In potato tubers, patatin isoforms play dual roles as storage proteins and as acyl hydrolases and esterases, indicating the contribution of this protein to the resistance to pests, pathogens and abiotic stress (Shewry 2003; Barta and Bartova 2008). Proteases, also called peptidases or proteinases, are enzymes that are responsible for peptide bond hydrolysis. These enzymes play roles in the regulation of many different processes in plant physiology and development, including an important role in pathogenesis (Van der Hoorn 2008). In potato tubers, bacterial pathogens degrade host tissue by producing cell wall-degrading enzymes, such as pectin and pectate lyases, polygalacturonases, cellulases, and proteases (Pérombelon 2002; Toth et al. 2003). The strain D. solani IFB0099 (IPO 2276) used in our study is most closely related to the strain IPO 2222, which is most closely related to Dickeya dadantii 3947 (Van der Wolf et al. 2014; Pritchard et al. 2013), which produces a range of plant cell wall-degrading enzymes (Py et al. 1991; Kazemi-Pour et al. 2004). Protease activity is regulated by protease inhibitors (PIs). PIs are among the important defense related proteins and are divided into four classes: serine protease, metalloprotease, cysteine protease, and aspartic protease inhibitors (Van der Hoorn 2008; Rehman et al. 2017). PIs are associated with resistance and can inhibit proteinases of various plant pathogens (Ryan 1990). Mature potato tubers are rich sources of several PI families (Fischer et al. 2015). Valueva et al. (1998) have shown that infection of potato tubers by Phytophthora infestans induces the accumulation of chymotrypsin inhibitors, which significantly inhibit mycelial growth of this pathogen. Senser et al. (1974) demonstrated the inhibitory effects of PIs from potato tubers on the diversity of proteolytic microorganisms isolated from rotten potato tubers. Significant changes in protein abundance were also observed when comparing potato cultivars that were resistant or susceptible to bacteria after mock inoculation. Regardless of inoculation treatment (inoculation with bacteria and mock-inoculation), the five potato cultivars had different protein pattern when analyzed by principal component analysis. The effect of treatments was smaller than the effect of potato cultivars, which could be separated along with level of resistance based on PC2. These results indicated that in the resistance tests, protein expression is partially associated with the response of tuber tissue to mechanical wounding. Ten proteins exhibited significantly increased levels in the resistant clone DG 00–270 compared to the susceptible clone DG 08–305, both after D. solani inoculation and mock inoculation (Table 3). Only two of these proteins, namely, PI2 and PT3K1, were selected in the tetraploid cultivars. The most prominent fold change between the diploid clones was observed for MCPI and MCPI-like proteins. These protein exhibited much higher fold changes following pathogenic invasion than after mock inoculation and in intact tubers. MCPIs are known to be wound-inducible inhibitors in plants (Molesini et al. 2018). In potato tubers, these proteins are among the most abundant PIs (Fischer et al. 2015). MCPI accumulation in wounded leaves of potato plants has been reported (Graham and Ryan 1981). In transgenic rice, a potato carboxypeptidase inhibitor gene was shown to provide resistance against the economically important pathogens Magnaporthe oryzae and Fusarium verticillioides (Quilis et al. 2007). The higher abundance of MCPIs in DG 00–270 than in DG 08–305 indicates the defensive role of these proteins in the response to bacterial inoculation. Notably, in our study, changes in MCPI abundance between the potato cultivars were not statistically significant. This result suggests that different mechanisms of resistance to D. solani are present in the cultivated and diploid potato plants.

Protein abundances reflect a dynamic balance among post-transcriptional, translational and protein modification/destruction processes (Vogel and Marcotte 2012). Complementary approaches, such as quantitative genomics, allow the identification of genetic factors involved in the defensive response of potato to D. solani, which is the aim of our future studies.