Impact Factor 4.106 | CiteScore 4.47
More on impact ›

Original Research ARTICLE

Front. Plant Sci., 17 September 2015 |

A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones

  • Department of Plant Protection Biology, Swedish University of Agricultural Sciences, Alnarp, Sweden

Comparative transcriptomics between species can provide valuable understanding of plant-pathogen interactions. Here, we focus on wild Solanum species and potato clones with varying degree of resistance against Phytophthora infestans, which causes the devastating late blight disease in potato. The transcriptomes of three wild Solanum species native to Southern Sweden, Solanum dulcamara, Solanum nigrum, and Solanum physalifolium were compared to three potato clones, Desiree (cv.), SW93-1015 and Sarpo Mira. Desiree and S. physalifolium are susceptible to P. infestans whereas the other four have different degrees of resistance. By building transcript families based on de novo assembled RNA-seq across species and clones and correlating these to resistance phenotypes, we created a novel workflow to identify families with expanded or depleted number of transcripts in relation to the P. infestans resistance level. Analysis was facilitated by inferring functional annotations based on the family structure and semantic clustering. More transcript families were expanded in the resistant clones and species and the enriched functions of these were associated to expected gene ontology (GO) terms for resistance mechanisms such as hypersensitive response, host programmed cell death and endopeptidase activity. However, a number of unexpected functions and transcripts were also identified, for example transmembrane transport and protein acylation expanded in the susceptible group and a cluster of Zinc knuckle family proteins expanded in the resistant group. Over 400 expressed putative resistance (R-)genes were identified and resistant clones Sarpo Mira and SW93-1015 had ca 25% more expressed putative R-genes than susceptible cultivar Desiree. However, no differences in numbers of susceptibility (S-)gene homologs were seen between species and clones. In addition, we identified P. infestans transcripts including effectors in the early stages of P. infestans-Solanum interactions.


Wild Solanum species are alternative hosts for pathogens of cultivated potato, which is a crop attacked by many different microorganisms. Local populations of these often grow in close proximity of potato fields. Compared to cultivated potatoes, which are not well adapted to the local environment and also might have lost important genetic variation through domestication, wild species could be a source for desirable traits to breed disease-resistant varieties of the crop.

The oomycete Phytophthora infestans is the causing agent of late blight leading to annual losses in potato worth billions of US dollars (Haverkort et al., 2008). In the last few decades, many resistance (R-)genes against this devastating pathogen have been successfully transferred from various wild Solanum species to potato (Vleeshouwers et al., 2011; Rodewald and Trognitz, 2013). However, the pathogen contains fast-evolving effector genes localized in a highly dynamic region of its genome (Haas et al., 2009), enabling it to quickly overcome the R-gene based resistance in potato.

Functional stacking of R-genes has been employed to achieve resistance (Zhu et al., 2012) and new cultivars such as Sarpo Mira have a wider mix of R-genes introgressed (Rietman et al., 2012). The emergence of next-generation sequencing techniques enable inexpensive characterization of expressed R-genes also in non-model plants (Draffehn et al., 2013; Yang et al., 2014). Apart from R-gene conferred resistance, studying Solanum lineage-specific adaptive responses and comparing wild Solanum and potato interactions with pathogens can provide valuable insights into the underlying molecular events (Draffehn et al., 2013; Du et al., 2013; Ali et al., 2014).

Three wild Solanum species, which are potential hosts for P. infestans, grow in Sweden. Bittersweet nightshade (S. dulcamara) is a diploid semi-woody, outcrossing, perennial climber, while the di- tetra- or hexaploid European black nightshade (S. nigrum) and diploid hairy nightshade (S. physalifolium) are self-compatible, annual weeds. Populations of these species grow in the proximity of potato cultivations but differ in sensitivity to P. infestans. S. dulcamara and S. nigrum are only rarely infected by P. infestans but still can harbor the pathogen (Deahl et al., 2004, 2010), while S. physalifolium has been reported to be a sensitive host, which might influence pathogen aggressiveness in potato (Andersson et al., 2003; Grönberg et al., 2012). S. dulcamara has previously been shown to vary in resistance toward P. infestans, and R-genes, Rpi-dlc1 and Rpi-dlc2, have been mapped (Golas et al., 2010, 2013). A comprehensive characterization of the S. dulcamara transcriptome and genome-wide genetic map have also been presented (D'Agostino et al., 2013).

Coevolution between plants and their natural enemies is believed to be an important driver of biological diversity (Burdon and Thrall, 2009). Variation in transcriptomes among wild species can therefore be expected to be largely influenced by natural selection on defense to pests and pathogens, in addition to selection in relation to beneficial biotic interactions and abiotic factors (Oh et al., 2014). In crop species, on the other hand, artificial selection caused by the domestication process and subsequent breeding have been shown to have direct effect on the transcriptome (Swanson-Wagner et al., 2012; Bellucci et al., 2014).

The aim of this study was to characterize defense responses against P. infestans in wild Solanum species and to identify factors successfully conferring resistance. Using RNA-seq we therefore generated leaf transcriptome atlases of individuals of S. dulcamara (accession Sd 3:6), S. nigrum (accession Sn 4:3) and S. physalifolium (accession Sp 2:4) and compared these to the transcriptomes of three tetraploid potato clones: Sarpo Mira and SW93-1015, which exhibit strong P. infestans resistance (Ali et al., 2012), and Desiree (cv.), which is susceptible. We have previously explored the response to infection with P. infestans in these potato clones by gene expression microarrays and proteomics (Ali et al., 2014).

In order to capture lineage-specific events reflecting adaptive evolution and conferring resistance, we used a novel workflow incorporating a correlation-based method to identify clusters either expanded or depleted in transcript numbers dependent on the level of P. infestans resistance of the wild Solanum and potato clones studied. An overview of the workflow used is presented in Figure 1. We used class vectors describing the level of plant resistance either qualitatively, i.e., resistant vs. susceptible, or quantitatively, i.e., based on an experimentally determined resistance gradient. Clusters were analyzed functionally by gene ontology (GO) and clusters populated with R-genes studied in more detail. We also identified 100's of putative R-genes in the three wild Solanum species and three potato clones. Furthermore, P. infestans transcripts representing RXLR effectors, Crinklers (CRN) and elicitins were detected.


Figure 1. Overview the workflow used to identify resistance and susceptible factors in the host and effectors of the infecting pathogen based on RNA-seq. (A) Phenotypic analysis and RNA Sequencing (RNA-Seq). (B) Bioinformatic analysis of the RNA-seq data and transcriptomes.

Materials and Methods

Plant Material, Growth Conditions, and Resistance Assay

Seeds of wild Solanum species (S. dulcamara, accession Sd 3:6; S. nigrum, accession Sn 4:3; and S. physalifolium, accession Sp 2:4) were collected from natural habitats in Southern Sweden. Wild Solanum and potato clones Desiree (cv.), Sarpo Mira, and SW93-1015 were grown in a growth chamber set to 20°C, 16:8 h light:dark cycle and 70% relative humidity. Based on a detached leaf resistance assay with P. infestans (strain SE-03058) performed as described in Ali et al. (2012) one individual from each wild Solanum species and the three potato clones were selected for transcriptome profiling. Detached leaves without inoculation, hereafter referred to as at 0 h after inoculation (hai), at 24 and 48 hai were frozen in liquid nitrogen for RNA extraction.

For the resistance assay, detached leaf assays with five leaves per each clone/species plant were inoculated with P. infestans. Long and short axes of the lesion on each leaf were measured 7 days after inoculation, and then averaged to show the infection level quantitatively in each clone/species. Based on the lesion size, a gradient (from most resistant to most susceptible) was determined. The resistance gradient was statistically validated by a Kruskal-Wallis test with the PCMCR R-package (Pohlert, 2014) for pairwise multiple comparison between mean ranks of the lesion size.

RNA Isolation and Sequencing

RNA was isolated from fully expanded leaves detached from 6-week-old plants and incubated in conditions used for the resistance assay. For the potato clones leaves were collected before infection and at 24 hai, as described in Ali et al. (2014), whereas the wild Solanum leaves were harvested at 0, 24, and 48 hai.

The RNA extractions were performed using the RNeasy Plant Mini kit (Qiagen GmbH, Hilden, Germany). Samples were DNase treated and cleaned using the Qiagen RNA cleanup kit. RNA concentration and purity (260/280 nm > 1.8) was determined with an ND-1000 NanoDrop (Wilmington, USA) and integrity of the samples were analyzed with an Experion™ (Bio-Rad Laboratories, Hercules, USA). Equimolar amount of RNA from the leaves at 0, 24, and 48 hai was pooled for each wild Solanum before sequencing using Illumina HiSeq 2000 at Beijing Genomics Institute (BGI, Beijing, China). The RNA-seq data has been deposited as de novo assembled FASTA files in ArrayExpress: E-MTAB-2953.

Transcript Assembly, Identification of Orthologous Groups

The raw reads were processed with Trimmomatic v 0.3 (Bolger et al., 2014) in order to remove low-quality bases. The P. infestans reads were removed from the Solanum datasets by mapping the reads on the P. infestans genome using TopHat2 (version 0.8b) (Kim et al., 2013). Three independent de novo assemblies were then built from the remaining paired-end reads with Trinity (r2013_08_14) (Grabherr et al., 2011) using default parameters. To estimate the assemblies' quality bowtie2 (v. 2.1.0) was used to calculate the overall alignment rate of the quality-trimmed reads of the transcriptome assemblies. For each Trinity component a representative sequence was chosen, based on length, and Transdecoder was then used to obtain the longest protein sequence. The three protein sequence FASTA files obtained were then used as an input for OrthoMCL (Li et al., 2003) along with proteomes from Arabidopsis thaliana (TAIR10), S. lycopersicum (ITAG2.3), O. sativa (RGAP7), S. tuberosum group Phureja (PGSC3.4) and the International Tomato Annotation Group (ITAG1.0). Three additional proteomes of S. tuberosum clones Desirée (cv.), Sarpo Mira and SW93-1015 derived from either uninfected or P. infestans infected (isolate SE-03058) transcriptomes assembled as described in Ali et al. (2014) were finally added. OrthoMCL (v2.0.3) (Li et al., 2003) was run on these 11 proteomes with default parameters.

Identification of Expanded/Depleted OrthoMCL Clusters

In order to identify expanded or depleted OrthoMCL clusters, prototype vectors holding the differences in transcript numbers for the resistant and susceptible species and clones qualitatively (i.e., resistant vs. susceptible species and clones) or quantitatively (i.e., related to a gradient vector based on the results of the resistance assay, described above) were constructed. For both analyses a matrix of transcript numbers for each OrthoMCL clusters of all 6 clones/species was generated (Table S1). Additionally, two prototype vectors were created reflecting the qualitative and quantitative susceptibility-resistance relationship between the plants. Thereafter association between paired samples was tested by Spearman's rho statistics (test for rank-based measure of association) between each OrthoMCL cluster and each prototype vector in order to determine positive or negative associations. This gave four p-values for each OrthoMCL cluster that measure quantitative and qualitative positive and negative associations (Table S1), representing expanded and depleted clusters (p < 0.05), respectively.

Functional Annotation of the Expanded and Depleted OrthoMCL Clusters

GO enrichment analysis of the expanded and depleted OrthoMCL clusters was done by a hypergeometric test using GOEAST (Zheng and Wang, 2008) with default parameters. In order to create an annotation file, all OrthoMCL clusters were assigned GO terms based on the genomic sequences belonging to the cluster. In this way the potato clones and Solanum wild species sequences inherited the annotations of the genome sequences dependent on the OrthoMCL structure. The following functional annotations were used: TAIR10 (Arabidopsis), RGAP7 (rice), ITAG2.4. For the two potato gene models PGSC and ITAG GO terms were used as assigned in Amar et al. (2014). Related processes were displayed as significant bootstrap clusters (Suzuki and Shimodaira, 2006) of GO term semantic similarities calculated using the GOSimSem package (Yu et al., 2010).

Identification of Plant Resistance (R)-, Susceptibility (S)-genes, and P. infestans Effectors

Two complementing ways were used to identify R-genes. Firstly, a discriminatory psp-gen MEME (Bailey et al., 2009) analysis identified 20 NB-LRR protein sequence motifs (Table S3) using the same positive NB-LRR and negative non-NB-LRR sequence training sets as in Jupe et al. (2012) keeping their classification of the motifs. MAST (Bailey et al., 2009) was then used to search for occurrences of the 20 motifs in the wild Solanum species and potato clone transcript sequences. Transcripts were considered candidate NB-LRRs if the reported MAST E-values were lower than the smallest E-value observed for a member of the negative training set in Jupe et al. (2012) (i.e., E < 8.5e-24). The NB-LRRs were classified as CNLs or TNLs based on the presence of previously classified MEME motifs (Table S4). Alternatively, BLASTP searches (parameters: -evalue 1e-5 -outfmt 6 -max_target_seqs 1) with 112 reference resistance genes deposited in the Plant Resistance Genes database (Sanseverino et al., 2012) were performed (Table S5). For the BLASTP results a sequence cut-off of 30% identity and 30% coverage was applied and then classified according to the domain classification of their homologs in the Plant Resistance Genes database. The resulting R-genes were validated with a Pfam domain search (Finn et al., 2014) (Table S5). In order to identify specific R-gene members BLASTP and BLASTN searches (parameters: -evalue 1e-5 -outfmt 6 -max_target_seqs 1) were performed with R1, R2, R3, and Rpi-blb1 genes as baits (Table S6). S-genes were identified by BLASTP and BLASTN (parameters: -evalue 1e-5 -outfmt 6 -max_target_seqs 1) with A. thaliana mitogen-activated protein kinase 4 (MAPK 4), A. thaliana homoserine kinase 1 and 6 (dmr1, dmr6), and the Barley mlo together with its Capsicum and S. tuberosum homologs as baits (Table S7).

To identify effectors of P. infestans acting during the infection process in different species, sequence reads of each wild Solanum and potato clone were mapped and aligned using (TopHat v2.0.8b) on the P. infestans strain T30-4 reference genome ( Unique reads mapped to the P. infestans reference genome were counted using the htseq-count program (HTSeq-0.5.4p3). The P. infestans annotation file from The Broad institute (downloaded 08/07/2014; was used to identify the potential RxLR effectors, Crinklers (CRNs) and elicitins among the P. infestans reads (Table S8).

Handling and Visualization of Data

R (v3.0.3) and Perl (v5.18.2) were used to handle, parse, transform and statistically analyze the data. Venn diagrams was produced by VENNERABLE (Swinton, 2011). Qlucore Omics Explorer (v3.0) was used for construction of principle component analysis (PCA).

Results and Discussion

De novo Transcriptome Assembly

Between 54 and 90 Mb clean paired-end sequence reads of S. dulcamara, S. nigrum, S. physalifolium, Desiree, Sarpo Mira, and SW93-1015 were de novo assembled using Trinity (Table 1). The reads were of high and consistent quality for all species and clones (Figure S1). The amount of clean reads should be sufficient to generate in-depth, representative assemblies of leaf transcriptome atlases for the species and clones studied here, since 20 million reads for tissue samples has been shown to be sufficient for non-model mammal species (Francis et al., 2013). There was no big variation in the GC content of the assemblies, 40–41% (Table 1), which is similar to the content in the coding domain sequences predicted from the potato genome (43%).


Table 1. Summary of transcriptome assemblies for three wild Solanum species and three potato clones.

The N50 value for the contigs, generated using the longest representative, were larger than the values recently reported for closely related species S. melongena and S. torvum (Yang et al., 2014). This could be due the longer read length (90 bp) of the RNA-seq data we used and is indicating a high quality of the assemblies with a high level of represented leaf transcripts. In all the assemblies, the contig N50 value generated using the longest isoform was smaller (736–1012 bp) than the N50 derived from all the assembled transcripts (797–1402 bp) (Table 1). The lower N50 value in the assembly using only the longest isoform could be associated with the reduced total sequence length (Table 1) due to the removal of redundant transcripts. However, commonly used reference-free measures including median contig length, number of contigs and N50 can be misleading as measures of assembly quality. Therefore, we mapped all the transcripts back to the assemblies and found that overall alignment rates were constant and high for all six assemblies between 91.7 and 95.2% (Table 1).

The number of transcripts identified is within the range of what can be expected in plant leaf transcriptomes, granted that both tomato and potato have ca 35,000 predicted genes (Tomato Genome Consortium, 2012). There was a slightly higher number of transcripts based on the longest representative in the wild Solanum species (26,156–27,988) than in the potato clones (23,925–25,045), which is in line with previous reports on decreased diversity of gene expression in domesticated crops (Bellucci et al., 2014) (Table 1).

Resistance Levels of Species and Clones

Both the S. nigrum (Sn 4:3) and S. dulcamara (Sd 3:6) accessions used in this study are resistant, whereas the S. physalifolium (Sp 2:4) accession is susceptible to the P. infestans strain. Potato clones responded to P. infestans infection as expected. Desirée (cv.) is a well-studied food potato with a high degree of susceptibility to P. infestans infection. Sarpo Mira is a potato cultivar that recognizes five different effectors from P. infestans and shows a classical HR reaction in response to inoculation with the pathogen (Rietman et al., 2012). SW93-1015 is a breeding potato clone, which is resistant to Swedish P. infestans populations with reduced HR expansion (Ali et al., 2012).

We were only able to detect visible macroscopic symptom of the P. infestans inoculation after 48 hai as noted for potato in a review by Fry (2008). The 45 hai is referred as biotrophic phase (Vleeshouwers et al., 2000) whilst the 72 hai is considered the start of the necrotrophic phase of the pathogen (Birch et al., 2003). Consequently, RNA was isolated either up to 24 hpi (potato clones) or 48 hpi (wild Solanum species) before macroscopic lesions occurred in any of the species or clones and before the onset of the necrotrophic phase. Therefore, the samples in this study cover only the biotrophic phase of the P. infestans infection. Although our analysis groups the samples as either as susceptible or resistant to make general findings, the samples should not be regarded as proper biological replications. That sort of comparison is difficult, if not impossible, to accomplish, since progression of infection inevitably varies between species and clones with varying resistance levels. Based on the resistance assay using detached leaves a resistant to susceptible gradient was established in the following order: Sarpo Mira, SW93-1015, S. nigrum, S. dulcamara, Desiree, and S. physalifolium (Figure 2). Applying the Kruskal-Wallis test on mean rank of the lesion size [χ(5)2=26.1, P < 0.01] showed a resistance difference and the grouping was in line with the gradient.


Figure 2. Detached leaves from plants of three wild Solanum species, (S. nigrum, Sn 4:3; S. dulcamara, Sd 3:6; and S. physalifolium, Sp 2:4), and potato clones (Sarpo Mira, SW93-1015, and Desiree) inoculated with Phytophthora infestans. The pathogen grew on the susceptible Desiree and Sp 2:4 (gray bars) whilst its growth was restricted to the inoculation site on the resistant ones (white bars). Pictures were taken at 7 days after inoculation (dpi). Error bars indicate standard error of the mean. Bars that are not connected by the same letter are significantly different (P < 0.05).

In plant resistance biology detached leaf assays are widely used even if they do not reflect the physiology of the intact plant or environmental conditions in the field. In more in depth studies it is necessary to relate the results between the two systems to each other. Here, we have only studied the infection progress and transcriptomes under detached leaf conditions, to which the species and clones might react differently. This also overlooks reactions in other environments such as field conditions. However, the correlation-based method identifying expanded and depleted transcript families based on resistance level is dependent on relating the phenotypic response during the same conditions, justifying the capturing the transcriptomes during detached leaf conditions. Future more comprehensive studies should take the genotype x environment interactions into account by testing different environments and pathogen strains and their influence on resistance levels and transcriptome compositions. The overall workflow of the present study is available in Figure 1.

Comparative Analysis of OrthoMCL Cluster Composition

Using the predicted proteomes of the 11 species and clones in OrthoMCL, 38,890 clusters were generated, in which leaf transcripts of Desiree, Sarpo Mira, SW93-1015, S. nigrum, S. dulcamara, and S. physalifolium were represented in 27,979 (Table S2). In the analyses of the transcriptomes only the longest representative transcript was used to avoid allelic variants caused by different polyploidy, mis-assemblies and sequencing errors (D'Agostino et al., 2013). The OrthoMCL cluster structure was analyzed by PCA (Figure 3) and GO term enrichment analysis (Figure 4; Figure S2). Displaying the OrthoMCL clusters constructed by transcripts of Desiree, Sarpo Mira, SW93-1015, S. nigrum, S. dulcamara, and S. physalifolium and genes S. tuberosum (PGSC or ITAG gene models), S. lycopersicum and A. thaliana, there are clear differences between transcript and gene data (principle component 1, PC1) as can be expected since not all genes are present in the leaf transcriptomes studied here. A division of S. tuberosum clones and wild Solanum follows PC2. Furthermore, OrthoMCL clusters expanded or depleted as for numbers of transcripts were identified and either linked to a gradient depending on the degree of resistance to P. infestans (quantitative analysis) or grouped as resistant vs. susceptible (qualitative analysis). These analyses found 143 and 1291 OrthoMCL clusters significantly (p < 0.05) expanded quantitatively and qualitatively, respectively, and 249 and 525 significantly depleted quantitatively and qualitatively, respectively (Table S1). In the qualitative analysis as we contrast two susceptible with four resistant species and clones and therefore we have unequal sequence amounts for the resistant and susceptible entities, 158.8 and 77.7 Mb, respectively, which could influence results. However, this unequal relation is not seen in the quantitative analysis and we find a substantial overlap of 86% in the clusters identified by both analyses showing coherence between both approaches.


Figure 3. OrthoMCL principle component analysis (PCA). The PCA is based on number of transcripts for Solanum nigrum, S. dulcamara, and S. physalifolium, Desiree, Sarpo Mira, SW93-1015 as well as number of genes for S. tuberosum (PGSC or ITAG gene models), S. lycopersicum, and A. thaliana.


Figure 4. Clades of GO enriched terms for Biological Process (Siddappa et al., 2014) of the expanded and depleted OrthoMCL clusters determined quantitatively according to the resistance gradient (GO enriched terms based on the qualitative analysis is available in Figure S2). R and S denotes GO terms enriched in the resistant or susceptible groups, respectively. Broad, unspecific GO terms have been removed.

Accession- and Species-specific OrthoMCL Clusters

Under detached leaf conditions S. dulcamara had slightly fewer transcripts represented in the OrthoMCL clusters compared to the others, which had representatives in roughly 2/3 of all OrthoMCL clusters (Table 2). With 456 OrthoMCL clusters only containing SW93-1015 transcripts, this clone had the highest number of unique clusters (Table 2; Figure 5). However, since these observed differences in expressed transcripts are from detached leaf assays solely, to which the species/clones might react differently to the incubation condition, results on species and clone-specific features should be seen as preliminary until the transcriptomes from several environments and most importantly from field and natural habitat conditions have been obtained.


Table 2. OrthoMCL cluster summary.


Figure 5. OrthoMCL cluster overlap between Solanum nigrum (yellow), S. dulcamara (red), and S. physalifolium (lilac), Desiree (blue), Sarpo Mira (green), SW93-1015 (orange) based on detected transcripts.

A GO enrichment analysis of the unique OrthoMCL clusters (Table S9) found apoptotic process, response to oxidative stress, chitin binding and chromatin binding corresponding to clusters of transcripts with sequence similarity to Rpi-vnt1, other NBS-LRRs, 42 kDa chitin-binding protein and MYB factors in SW93-1015. The Rpi-vnt1 gene from S. venturii confers late-blight resistance in both potato and tomato (Foster et al., 2009). Interestingly, in Sarpo Mira the GO term terpenoid biosynthetic process was uniquely enriched. Among the wild Solanum enrichment of defense response to fungus and regulation of peptidase activity in S. dulcamara, peptidase activity acting on L-amino acid peptides and histone modification in S. nigrum and cell wall modification and epidermis development in S. physalifolium are noteworthy with regards to plant resistance. Enriched GO terms of unique OrthoMCL clusters specific for each accession are found in Table S9.

Transcript Clusters Expanded in the Resistant Group

The qualitative (resistant vs. susceptible species and clones) and quantitative (resistance gradient) analyses identified OrthoMCL clusters significantly expanded in the resistant group. GO term enrichment analysis of expanded groups identified terms such as protein phosphorylation, defense response, plant-type hypersensitive response, host programmed cell death induced by symbiont and aspartic-type endopeptidase activity (Figure 4; Figure S2). These are all associated to plant defense and were therefore studied in more detail.

A number of expanded transcript clusters were populated by transcripts encoded by resistance genes showing the usefulness of the method. Examples include several clusters with the NBS-LRR resistance family members and a large cluster of R1 Late blight resistance homologs. This specific OrthoMCL cluster (WILD0000054) had 13 and 11 expressed members in SW93-1015 and Sarpo Mira, respectively, and 10 members in S. nigrum, whereas susceptible Desiree and S. physalifolium only had seven and six members, respectively. Many of the clusters with NBS-LRR representatives were specific to the Solanum species lacking transcript members from Arabidopsis and rice. The protein of the NBS-LRR resistance gene (PGSC0003DMP400054560) shows high sequence identity (60%) to a potato Bacterial spot disease resistance protein 4 (Bs4) homolog. In tomato, BS4 (Solyc05g007850) gives resistance against Xanthomonas campestris causing bacterial leaf spots (Schornack et al., 2004). In S. demissum the BS4 homolog is flanking the functional R1 locus (Kuang et al., 2005).

An expanded defense-related cluster contained Glutamate-gated kainate-type ion channel receptor subunits (GluR5s), which are ligand-gated non-selective cation channels. Overexpression in Arabidopsis delays fungal infection of fungal pathogen Botrytis cinerea (Kang et al., 2006) and have been shown to mediate MAMP-induced calcium influx in A. thaliana (Mark et al., 2011).

Some expanded transcript clusters in the resistant group without clear connections to defense were identified. For example, a small cluster of Zinc knuckle family proteins, which has been shown to have RNA binding and splicing activities (Lopato et al., 1999), and Diminutivo 1 involved in Brassinosteroids biosynthesis and lignin formation (Hossain et al., 2012).

Transcript Clusters Expanded in the Susceptible Group

GO term enrichment analysis associated clusters expanded in the susceptible group to terms such as lipoprotein biosynthetic process, signal transduction, protein acylation, transmembrane transport and nitrogen biosynthesis (Figure 4; Figure S2). An example of a putative susceptibility factor present in one of the clusters is the cytokinin-regulated kinase 1 (CRK1). In tobacco (Schäfer and Schmülling, 2002), CRK1 of tobacco is a class I receptor kinase located in transmembrane, hypothesized to be involved in hormone signaling.

Transcripts related to ATPase activity, such as the vacuolar H+-pumping ATPase required for salinity stress tolerance in potato (Queirós et al., 2009), was overrepresented in the susceptible group. Histidine kinase is a vital component in signaling transduction and used by plants to sense and respond to biotic and abiotic stresses were also expanded. As were transcripts for sugar and carbohydrate transport, which could be integral part of the network distributing photosynthetically assimilated carbon and transcripts for development of embryo, reproductive structure, fruit and seed.

Even if fewer than in the resistant group, some GO terms associated to plant defense responses were enriched also in the susceptible group, and could be components of the basal disease resistance related to PAMP-triggered immunity (PTI) present also in susceptible hosts (Jones and Dangl, 2006). One cluster was populated with Avr9/Cf-9 rapidly elicited protein 189 (ACRE189). The ACRE189 encodes widely conserved F-box/leucine rich proteins, which are up-regulated upon elicitation by the race-specific avirulence peptide Avr9 in tomato and tobacco, and plays a crucial role in HR development (van den Burg et al., 2008). GO enrichment identified transcripts of serine/threonine protein kinases, and a member of this protein kinase class, Avr9/Cf-9 inducing kinase 1 (ACIK1), is important in the Cf9/Avr9- and Cf4/Avr4-mediated HR response against the tomato leaf mold fungus Cladiosporium fulvum, which contains the Avr9 (Rowland et al., 2005). Another expanded transcript contained zinc finger binding proteins, which could be integral components of some resistance genes in potato against potato virus X (PVX) and nematodes (Gupta et al., 2012).

Comparison of Expressed R-gene Homologs

The R-gene encoded R-proteins, most of which are nucleotide-binding-leucine-rich-repeat (NB-LRR) type, recognize avirulence effectors of the pathogen during the infection process (Hogenhout et al., 2009; Thomma et al., 2011; Giraldo and Valent, 2013) resulting in a hypersensitive response (HR), which restricts expansion of the pathogen. Most of the cloned NB-LRR proteins from Solanum species belong to the N-terminal coiled coil-nucleotide-binding site-leucine-rich repeat (CNL), Toll interleukin1 receptor-nucleotide-binding site-leucine-rich repeat (TNL) and receptor-like protein (RLP) classes (van Ooijen et al., 2007). A motif based search identified 438 NB-LRR type genes in the sequenced potato genome of S. tuberosum Group Phureja genome (Jupe et al., 2012). Still, Phureja is susceptible to P. infestans in spite of the large number of NB-LRR genes, and thus, these genes seem ineffective against P. infestans. Previously identified functional R-genes, belonging to the largest class CC-NB-LRR, conferring resistance against P. infestans includes the R1 from S. demissum (Ballvora et al., 2002), Rpi-sto2 from S. stoloniferum (Champouret, 2010), Rpi-blb1 and Rpi-blb2 of S. bulbocastanum (van der Vossen et al., 2003, 2005), and the Rpi-vnt1 of S. venturii (Foster et al., 2009; Pel et al., 2009).

R-gene homologs sharing NB-LRR domains in the wild Solanum were identified by the MEME/MAST analysis as in Jupe et al. (2012), and by blasting 112 reference R-genes from the Plant Resistance Genes database (Sanseverino et al., 2012) (Tables S4, S5). The identified motifs include CNLs, TNLs, RLPs, three truncated classes (kinase, which has a kinase domain and an extracellular leucine-rich repeat, NL and TN) and “Other” with no typical resistance-related domains, but which still have been described as conferring resistance through different molecular mechanisms, e.g., asc-1 (Table 3).


Table 3. Summary of expressed R-gene homologs in Solanum species.

The MEME/MAST analysis is powerful in identifying putative R-genes lacking high overall sequence identity with known R-genes, but which still contain conserved domains. However, with the training set used this analysis is restricted to NB-LRR protein identification (CNLs or TNLs) in contrast to identification by sequence identity of the whole genes. Furthermore, a MAST search with the 20 MEME motifs identified only 69 of the 112 reference R-genes in the Plant Resistance Genes database. The discrepancy between the two methods to find putative R-genes justifies the use of both approaches. The validation done with Pfam indeed showed that the vast majority contain the NB-ARC domain, indicative of putative R-genes (Tables S4, S5). Figure S3 shows the overlap of orthoMCL clusters containing R-genes between Solanum wild species and potato clones.

According to the MEME/MAST analysis, we found both TNL and CNL representatives in the transcript assemblies represented by 3–11 and 28–64 putative R-genes, respectively (Table 3), reflecting the fact that the largest part of NB-LRR genes belonged to the CNL class in the Phureja genome (Jupe et al., 2012). SW93-1015 had the largest number of expressed putative R-gene (107). Overall, SW93-1015 and Sarpo Mira had 45–65% more expressed putative R-gene than the wild Solanum and 25–43% more than the susceptible cultivar Desirée (Table 3).

Closer analysis of putative R1, R2, R3 and Rpi-blb1 representatives based on sequence identity revealed a small number of genes (Table 4; Table S6). Interestingly, S. dulcamara and Sarpo Mira, which are resistant against P. infestans, harbor unique OrthoMCL clusters for these specific classes of R-genes (WILD0014801 and WILD0031982, and WILD0036062, respectively).


Table 4. Summary of expressed putative R1, R2, R3, and Rpi-blb1 genes.

In addition to R-genes, we investigated four families of susceptibility (S)-gene homologs whose impairment results in disease resistance. For example, a loss-of-function mutation of the SlMlo1 in tomato was recently demonstrated to decrease susceptibility to fungal pathogen Leveillula taurica (Zheng et al., 2013). Expect for the classic example of mlo, we investigated the putative expansion of DMR6, DMR1, and MPK4 (Petersen et al., 2000; van Damme et al., 2008, 2009). However, all species and lines had the same number and composition of S-genes expressed (not shown). So at least for the four S-gene families tested, no link existed between (S)-gene transcript composition and P. infestans susceptibility.

Our results suggest that it is possible to scan the composition of R-genes expressed by RNA-seq and characterize a number of putative R-genes from wild Solanum species. By focusing on the transcriptome non-expressed genes are eliminated, which are likely to be unimportant in the resistance reaction and thus of little use as novel resources of resistance.

Analysis of Transcription of P. infestans Genes

Even if our methods where not optimized to identify P. infestans transcripts we detected transcripts including effectors expressed during the biotrophic phase of the P. infestans-plant interaction. A total of 7769, 2612, 1471, 892, and 73 transcripts from P. infestans were identified in inoculated samples from S. physalifolium, S. nigrum, SW93-1015, Desiree, and S. dulcamara, respectively (Figure S4). However, this number of P. infestans transcripts does not follow the resistance gradient. We hypothesize that this can either be due to differences in host resistance mechanisms that prevent the pathogen at different stages of the infection process reflecting different layers of resistance responses or due to low overall detection of P. infestans reads, which were substantially lower than plant reads forming only 0.35–0.0003% of total reads detected for each species or clone. Still, these results indicate that RNA-seq studies, if optimized, can be a powerful tool to monitor the progress of P. infestans during the infection, not least before macroscopic lesions occur in relation to susceptibility or resistance of the host.

The P. infestans transcripts represented putative RxLR effectors, Crinklers (CRNs) and elicitins (Table 5; Table S8). No P. infestans transcript was detected in the resistant potato clone Sarpo Mira (Table 5). The P. infestans genome consists of 563 RXLR and 196 CRN genes out of the total 17,797 predicted protein coding genes (Haas et al., 2009). We found no less than 134 RXLR effectors of the pathogen in the susceptible S. physalifolium (Table 5; Table S8), whereas no P. infestans effectors were detected in S. dulcamara and Sarpo Mira, indicating an efficient early defense response by both hosts upon perceiving the pathogen associated molecular patterns (PAMPs).


Table 5. Number of expressed P. infestans effector transcripts identified in each species and clone.

The secreted RxLR effector peptide avrblb1 (ipiO1; PITG_21388) was detected in all but S. dulcamara and Sarpo Mira. Identified in a screening of P. infestans effectors in resistant potato clones (Vleeshouwers et al., 2008), the avrblb1 leads to HR in plants containing the Rpi-blb1 R-gene from S. bulbocastanum. The avrblb2 (PITG_04090, PITG_04086, PITG_04085) and avr2 (PITG_22870, PITG_22870, PITG_08943, PITG_13940) family of secreted RxLR effector peptides were also found during the infection of both S. nigrum and S. physalifolium. Secreted RxLR effector peptide Avr4 (PITG_07387), which has been studied previously (van Poppel et al., 2008) was only identified in S. physalifolium.

Among the identified elicitins were the elicitin-like proteins INF4 (PITG_21410) and INF6 (PITG_12556), elicitin precursor INF1 (PITG_12551) and elicitin INF2A-like protein (PITG_12561).

Conclusions and Perspective

Compared to cultivated potatoes, which might have lost genetic variation through domestication, wild relatives could provide a source of traits that can be used to breed disease-resistant cultivated varieties. This study provides a first resource of comparative transcriptomics for three wild Solanum species. Furthermore, we have contrasted these transcriptomes to those of three potato clones. In our analysis we identified a number of expanded and depleted transcript families by applying a gradient based on the level of resistance of the species and clones to P. infestans. Additionally, we present a number of putative R-genes based on sequence identity and presence of specific domains, which could be future candidates in resistance breeding. This study also shows that transcripts in the P. infestans and Solanum pathosystem can be analyzed simultaneously.

It is an important but challenging task to link phenotypes to transcript profiles. In this study several candidate transcript with an effect on the resistance phenotype were identified, such as the members of the R1 transcript family (WILD0000054) expanded in resistant Sarpo Mira and S. dulcamara. Furthermore, two clusters (WILD0014801 and WILD0031982) were found to be populated with an R-gene only present in S. dulcamara. The importance of these R-genes should be further tested by generating transgenic plants. There are a number of elements that can be seen as novel in this study and the combination of these certainly is. Firstly, our analysis of the Solanum material is based on correlations between resistance phenotypes and transcript family compositions. It has been common to contrast expressional changes between phenotypes focusing on differentially expressed genes ever since the microarray technique was introduced, but using the correlations of phenotype-transcript family compositions remain little explored. Secondly, instead of using genomes as input for OrthoMCL to identify gene families, which is common, transcript sequences determined by RNA-seq was used to identify expanded transcript families associated to resistance or susceptibility. Lastly, in order to identify processes associated to resistance or susceptibility of infection with P. infestans, the network structure created by OrthoMCL was used to transfer annotations from annotated transcriptomes to the unannotated wild Solanum transcriptomes. This enables to do GO term enrichment analyses and subsequent clustering by semantic similarity to improve visualization and help interpretation.

The rapidly decreasing cost of sequencing is increasing the number of available transcriptomes from species and accessions with varying level of resistance to biotic as well as abiotic stresses, and will, thus, enable more comprehensive comparative analyses on the RNA-level. In addition, new emerging sequencing techniques with longer read length will lessen the obstacle of transcript assembly in the future. This study illustrates the power of comparative transcriptomics and presents a new correlation-based approach useful to analyze this type of data.

Author Contributions

IF, KA, EP, EAnd and EAlex planned the study. KA performed the lab experiment with assistance of EAlex. Bioinformatics analysis was done by IF and EP. All authors interpreted the data, wrote and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


We thank PlantLink, Swedish foundation for Strategic Research and the Crafoord foundation (20120533) for financial support. We thank Graham Etherington for providing support as well as the settings of MEME. We thank Markus Stenemo and Benjamin Almeida for their bioinformatics support.

Supplementary Material

The Supplementary Material for this article can be found online at:

Figure S1. QC plots for the species and clones. S. nigrum (A), S. dulcamara (B), S. physalifolium (C), Sarpo Mira (D), SW93-1015 (E), and Desiree (F). For each basepair position a Box-and-Whisker plot is drawn. The y-axis on the graph shows the quality scores. The higher the score the better the base call. The background of the graph divides the y-axis into calls of very good quality (green), calls of reasonable quality (orange), and calls of poor quality (red).

Figure S2. Clades of GO enriched terms for Biological Process (Siddappa et al., 2014) in of the expanded and depleted OrthoMCL clusters determined qualitatively (resistant vs. susceptible). Broad, unspecific GO terms have been removed.

Figure S3. Venn diagram showing the overlap of the OrthoMCL clusters containing identified R-genes in the studied Solanum species and potato clones.

Figure S4. Overlap of P. infestans transcripts detected in the wild Solanum species and potato clones.

Table S1. Transcript numbers for the Solanum species and potato clones for each OrthoMCL cluster including p-values (Spearman's rho test) for the degree of quantitative and qualitative positive or negative association.

Table S2. OrthoMCL clusters with transcript PGSC identifiers.

Table S3. NB-LRR-specific amino acid motifs identified with psp-gen MEME. Motif number, motif sequence, domain group (TNL or CNL) and whether the motif is similar to a known motif are indicated.

Table S4. Transcripts of Desiree, S. nigrum and S. physalifolium, S. dulcamara, Sarpo Mira, and SW93-1015 identified by MEME. The classification (TNL or CNL) is provided as well as NBS-LRR motif validation by pfam. For the potential R-gene transcripts the OrthoMCL cluster ID is given as well as number of transcripts in the clusters and p-values (Spearman's rho test) for the degree of quantitative and qualitative association.

Table S5. Potential R-gene homolog transcripts in Desiree, S. nigrum and S. physalifolium, S. dulcamara, Sarpo Mira, and SW93-1015 identified by BLASTP searches against the Plant Resistance Genes database. The classification is provided as well as NBS-LRR motif validation by pfam. For the potential R-gene transcripts the OrthoMCL cluster ID is given as well as number of transcripts in the clusters and p-values (Spearman's rho test) for the degree of quantitative and qualitative association.

Table S6. Putative R1, R2, R3, and blb1 transcripts in Desiree, S. nigrum and S. physalifolium, S. dulcamara, Sarpo Mira, and SW93-1015 identified by BLASTP and BLASTN searches with R1, R2, R3, and blb1 representatives. For the potential R-gene transcripts the OrthoMCL cluster ID is given as well as number of transcripts in the clusters and p-values (Spearman's rho test) for the degree of quantitative and qualitative association.

Table S7. S-gene homologs among the transcripts in the Solanum wild species and clones identified by BLASTP searches using A. thaliana mitogen-activated protein kinase 4 (MAPK 4), A. thaliana homoserine kinase 1 and 6 (dmr1, dmr6) and Barley mlo and its Capsicum and S. tuberosum homologs.

Table S8. The number of transcripts of CRN, RxLR effector and elicitin in Desiree, S. nigrum and S. physalifolium, S. dulcamara, Sarpo Mira and SW93-1015 together with the accession ID and annotation.

Table S9. (A) GO terms enriched in unique OrthoMCL clusters in Desiree, S. nigrum and S. physalifolium, S. dulcamara, Sarpo Mira and SW93-1015. (B) GO terms enriched in unique OrthoMCL clusters in potato clones vs. Solanum wild species. The hypergeometic tests p-value are provided.


Ali, A., Alexandersson, E., Sandin, M., Resjö, S., Lenman, M., Hedley, P., et al. (2014). Quantitative proteomics and transcriptomics of potato in response to Phytophthora infestans in compatible and incompatible interactions. BMC Genomics 15:497. doi: 10.1186/1471-2164-15-497

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, A., Moushib, L. I., Lenman, M., Levander, F., Olsson, K., Carlson-Nilson, U., et al. (2012). Paranoid potato phytophthora-resistant genotype shows constitutively activated defense. Plant Signal. Behav. 7, 400–408. doi: 10.4161/psb.19149

PubMed Abstract | CrossRef Full Text | Google Scholar

Amar, D., Frades, I., Danek, A., Goldberg, T., Sharma, S. K., Hedley, P. E., et al. (2014). Evaluation and integration of functional annotation pipelines for newly sequenced organisms: the potato genome as a test case. BMC Plant Biol. 14:329. doi: 10.1186/s12870-014-0329-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, B., Johansson, M., and Jonsson, B. (2003). First report of Solanum physalifolium as a host plant for Phytophthora infestans in Sweden. Plant Dis. 87, 1538–1538. doi: 10.1094/PDIS.2003.87.12.1538B

CrossRef Full Text | Google Scholar

Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, 202–208. doi: 10.1093/nar/gkp335

PubMed Abstract | CrossRef Full Text

Ballvora, A., Ercolano, M. R., Weiss, J., Meksem, K., Bormann, C. A., Oberhagemann, P., et al. (2002). The R1 gene for potato resistance to late blight (Phytophthora infestans) belongs to the leucine zipper/NBS/LRR class of plant resistance genes. Plant J. 30, 361–371. doi: 10.1046/j.1365-313X.2001.01292.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellucci, E., Bitocchi, E., Ferrarini, A., Benazzo, A., Biagetti, E., Klie, S., et al. (2014). Decreased nucleotide and expression diversity and modified coexpression patterns characterize domestication in the common bean. Plant Cell. 26, 1901–1912. doi: 10.1105/tpc.114.124040

PubMed Abstract | CrossRef Full Text | Google Scholar

Birch, P., Avrova, A., Armstrong, M., Venter, E., Taleb, N., Gilroy, E., et al. (2003). The potato - Phytophthora infestans interaction transcriptome. Can. J. Plant Pathol. 25, 226–231. doi: 10.1080/07060660309507074

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 170, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Burdon, J. J., and Thrall, P. H. (2009). Co-evolution of plants and their pathogens in natural habitats. Science 324, 755. doi: 10.1126/science.1171663

PubMed Abstract | CrossRef Full Text | Google Scholar

Champouret, N. (2010). Functional Genomics of Phytophthora infestans Effectors and Solanum Resistance Genes. Ph.D. thesis Wageningen University, Wageningen.

Consortium, T. G. (2012). The tomato genome sequence provides insights into fleshy fruit evolution. Nature 485, 635–641. doi: 10.1038/nature11119

PubMed Abstract | CrossRef Full Text | Google Scholar

D'Agostino, N., Golas, T., van de Geest, H., Bombarely, A., Dawood, T., Zethof, J., et al. (2013). Genomic analysis of the native European Solanum species, S. dulcamara. BMC Genomics 14:356. doi: 10.1186/1471-2164-14-356

PubMed Abstract | CrossRef Full Text | Google Scholar

Deahl, K. L., Perez, F., Baker, C. J., Jones, R. W., Cooke, L., and McGrath, M. (2010). Natural occurrence of Phytophthora infestans causing late blight on woody nightshade (Solanum dulcamara) in New York. Plant Dis. 94, 1063–1063. doi: 10.1094/PDIS-94-8-1063B

CrossRef Full Text | Google Scholar

Deahl, K. L., Shaw, D. S., and Cooke, L. R. (2004). Natural occurrence of Phytophthora infestans on black nightshade (Solanum nigrum) in Wales. Plant Dis. 88, 771–771. doi: 10.1094/PDIS.2004.88.7.771A

CrossRef Full Text | Google Scholar

Draffehn, A. M., Li, L., Krezdorn, N., Ding, J., Lübeck, J., Strahwald, J., et al. (2013). Comparative transcript profiling by SuperSAGE identifies novel candidate genes for controlling potato quantitative resistance to late blight not compromised by late maturity. Front. Plant Sci. 4:423. doi: 10.3389/fpls.2013.00423

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, J., Tian, Z., Liu, J., Vleeshouwers, V. G., Shi, X., and Xie, C. (2013). Functional analysis of potato genes involved in quantitative resistance to Phytophthora infestans. Mol. Biol. Rep. 40, 957–967. doi: 10.1007/s11033-012-2137-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Finn, R. D., Bateman, A., Clements, J., Coggill, P., Eberhardt, R. Y., Eddy, S. R., et al. (2014). Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230. doi: 10.1093/nar/gkt1223

PubMed Abstract | CrossRef Full Text | Google Scholar

Foster, S. J., Park, T.-H., Pel, M., Brigneti, G., Sliwka, J., Jagger, L., et al. (2009). Rpi-vnt1. 1, a Tm-22 homolog from Solanum venturii, confers resistance to potato late blight. Mol. Plant Microbe Interact. 22, 589–600. doi: 10.1094/MPMI-22-5-0589

PubMed Abstract | CrossRef Full Text | Google Scholar

Francis, W. R., Christianson, L. M., Kiko, R., Powers, M. L., Shaner, N. C., and Haddock, S. H. (2013). A comparison across non-model animals suggests an optimal sequencing depth for de novo transcriptome assembly. BMC Genomics 14:167. doi: 10.1186/1471-2164-14-167

PubMed Abstract | CrossRef Full Text | Google Scholar

Fry, W. (2008). Phytophthora infestans: the plant (and R gene) destroyer. Mol. Plant Pathol. 9, 385–402. doi: 10.1111/j.1364-3703.2007.00465.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Giraldo, M. C., and Valent, B. (2013). Filamentous plant pathogen effectors in action. Nat. Rev. Microbiol. 11, 800–814. doi: 10.1038/nrmicro3119

PubMed Abstract | CrossRef Full Text | Google Scholar

Golas, T., van de Geest, H., Gros, J., Sikkema, A., D'Agostino, N., Nap, J., et al. (2013). Comparative next-generation mapping of the Phytophthora infestans resistance gene Rpi-dlc2 in a European accession of Solanum dulcamara. Theor. Appl. Genet. 126, 59–68. doi: 10.1007/s00122-012-1959-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Golas, T. M., Sikkema, A., Gros, J., Feron, R. M., van den Berg, R. G., van der Weerden, G. M., et al. (2010). Identification of a resistance gene Rpi-dlc1 to Phytophthora infestans in European accessions of Solanum dulcamara. Theor. Appl. Genet. 120, 797–808. doi: 10.1007/s00122-009-1202-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883

PubMed Abstract | CrossRef Full Text | Google Scholar

Grönberg, L., Andersson, B., and Yuen, J. (2012). Can weed hosts increase aggressiveness of Phytophthora infestans on potato? Phytopathology 102, 429–433. doi: 10.1094/PHYTO-07-11-0192

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, S. K., Rai, A. K., Kanwar, S. S., and Sharma, T. R. (2012). Comparative analysis of zinc finger proteins involved in plant disease resistance. PLoS ONE 7:e42578. doi: 10.1371/journal.pone.0042578

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, B. J., Kamoun, S., Zody, M. C., Jiang, R. H., Handsaker, R. E., Cano, L. M., et al. (2009). Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature 461, 393–398. doi: 10.1038/nature08358

PubMed Abstract | CrossRef Full Text | Google Scholar

Haverkort, A., Boonekamp, P., Hutten, R., Jacobsen, E., Lotz, L., Kessel, G., et al. (2008). Societal costs of late blight in potato and prospects of durable resistance through cisgenic modification. Potato Res. 51, 47–57. doi: 10.1007/s11540-008-9089-y

CrossRef Full Text | Google Scholar

Hogenhout, S. A., van der Hoorn, R. A. L., Terauchi, R., and Kamoun, S. (2009). Emerging concepts in effector biology of plant-associated organisms. Mol. Plant Microbe Interact. 22, 115–122. doi: 10.1094/MPMI-22-2-0115

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, Z., McGarvey, B., Amyot, L., Gruber, M., Jung, J., and Hannoufa, A. (2012). DIMINUTO 1 affects the lignin profile and secondary cell wall formation in Arabidopsis. Planta 235, 485–498. doi: 10.1007/s00425-011-1519-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, J. D., and Dangl, J. L. (2006). The plant immune system. Nature 444, 323–329. doi: 10.1038/nature05286

PubMed Abstract | CrossRef Full Text | Google Scholar

Jupe, F., Pritchard, L., Etherington, G. J., Mackenzie, K., Cock, P. J., Wright, F., et al. (2012). Identification and localisation of the NB-LRR gene family within the potato genome. BMC Genomics 13:75. doi: 10.1186/1471-2164-13-75

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, S., Kim, H. B., Lee, H., Choi, J. Y., Heu, S., Oh, C. J., et al. (2006). Overexpression in Arabidopsis of a plasma membrane-targeting glutamate receptor from small radish increases glutamate-mediated Ca2+ influx and delays fungal infection. Mol. Cells 21, 418–427.

PubMed Abstract | Google Scholar

Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, R36. doi: 10.1186/gb-2013-14-4-r36

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, H., Wei, F., Marano, M. R., Wirtz, U., Wang, X., Liu, J., et al. (2005). The R1 resistance gene cluster contains three groups of independently evolving, type I R1 homologues and shows substantial structural variation among haplotypes of Solanum demissum. Plant J. 44, 37–51. doi: 10.1111/j.1365-313X.2005.02506.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Stoeckert, C. J., and Roos, D. S. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189. doi: 10.1101/gr.1224503

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopato, S., Gattoni, R., Fabini, G., Stevenin, J., and Barta, A. (1999). A novel family of plant splicing factors with a Zn knuckle motif: examination of RNA binding and splicing activities. Plant Mol. Biol. 39, 761–773. doi: 10.1023/A:1006129615846

PubMed Abstract | CrossRef Full Text | Google Scholar

Mark, K., Rik, H., Jens, M., Anja, R., and Ralph, P. (2011). Ionotropic glutamate receptor (iGluR)-like channels mediate MAMP-induced calcium influx in Arabidopsis thaliana. Biochem. J. 440, 355–365. doi: 10.1042/BJ20111112

PubMed Abstract | CrossRef Full Text

Oh, D.-H., Hong, H., Lee, S. Y., Yun, D.-J., Bohnert, H. J., and Dassanayake, M. (2014). Genome structures and transcriptomes signify niche adaptation for the multiple-ion-tolerant extremophyte Schrenkiella parvula. Plant Physiol. 164, 2123–2138. doi: 10.1104/pp.113.233551

PubMed Abstract | CrossRef Full Text | Google Scholar

Pel, M. A., Foster, S. J., Park, T. H., Rietman, H., van Arkel, G., Jones, J. D. G., et al. (2009). Mapping and cloning of late blight resistance genes from Solanum venturii using an interspecific candidate gene approach. Mol. Plant Microbe Interact. 22, 601–615. doi: 10.1094/MPMI-22-5-0601

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, M., Brodersen, P., Naested, H., Andreasson, E., Lindhart, U., Johansen, B., et al. (2000). Arabidopsis map kinase 4 negatively regulates systemic acquired resistance. Cell 103, 1111–1120. doi: 10.1016/S0092-8674(00)00213-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Pohlert, T. (2014). The Pairwise Multiple Comparison of Mean Ranks Package (PMCMR), R-package.

Queirós, F., Fontes, N., Silva, P., Almeida, D., Maeshima, M., Gerós, H., et al. (2009). Activity of tonoplast proton pumps and Na+/H+ exchange in potato cell cultures is modulated by salt. J. Exp. Bot. 60, 1363–1374. doi: 10.1093/jxb/erp011

PubMed Abstract | CrossRef Full Text | Google Scholar

Rietman, H., Bijsterbosch, G., Cano, L. M., Lee, H.-R., Vossen, J. H., Jacobsen, E., et al. (2012). Qualitative and quantitative late blight resistance in the potato cultivar Sarpo Mira is determined by the perception of five distinct RXLR effectors. Mol. Plant Microbe Interact. 25, 910–919. doi: 10.1094/MPMI-01-12-0010-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodewald, J., and Trognitz, B. (2013). Solanum resistance genes against Phytophthora infestans and their corresponding avirulence genes. Mol. Plant Pathol. 14, 740–757. doi: 10.1111/mpp.12036

PubMed Abstract | CrossRef Full Text | Google Scholar

Rowland, O., Ludwig, A. A., Merrick, C. J., Baillieul, F., Tracy, F. E., Durrant, W. E., et al. (2005). Functional analysis of Avr9/Cf-9 rapidly elicited genes identifies a protein kinase, ACIK1, that is essential for full Cf-9-dependent disease resistance in tomato. Plant Cell 17, 295–310. doi: 10.1105/tpc.104.026013

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanseverino, W., Hermoso, A., D'Alessandro, R., Vlasova, A., Andolfo, G., Frusciante, L., et al. (2012). PRGdb 2.0: towards a community-based database model for the analysis of R-genes in plants. Nucleic Acids Res. 41, 1167–1171. doi: 10.1093/nar/gks1183

PubMed Abstract | CrossRef Full Text | Google Scholar

Schäfer, S., and Schmülling, T. (2002). The CRK1 receptor-like kinase gene of tobacco is negatively regulated by cytokinin. Plant Mol. Biol. 50, 155–166. doi: 10.1023/A:1016087908746

PubMed Abstract | CrossRef Full Text | Google Scholar

Schornack, S., Ballvora, A., Gürlebeck, D., Peart, J., Ganal, M., Baker, B., et al. (2004). The tomato resistance protein Bs4 is a predicted non-nuclear TIR-NB-LRR protein that mediates defense responses to severely truncated derivatives of AvrBs4 and overexpressed AvrBs3. Plant J. 37, 46–60. doi: 10.1046/j.1365-313X.2003.01937.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Siddappa, S., Tiwari, J., Sindhu, R., Sharma, S., Bhardwaj, V., Chakrabarti, S., et al. (2014). Phytophthora infestans associated global gene expression profile in a late blight resistant Indian potato cv. Kufri Girdhari. AJCS 8, 215–222.

Google Scholar

Suzuki, R., and Shimodaira, H. (2006). Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics 22, 1540–1542. doi: 10.1093/bioinformatics/btl117

PubMed Abstract | CrossRef Full Text | Google Scholar

Swanson-Wagner, R., Briskine, R., Schaefer, R., Hufford, M. B., Ross-Ibarra, J., Myers, C. L., et al. (2012). Reshaping of the maize transcriptome by domestication. Proc. Natl. Acad. Sci. U.S.A. 109, 11878–11883. doi: 10.1073/pnas.1201961109

PubMed Abstract | CrossRef Full Text | Google Scholar

Swinton, J. (2011). Venn Diagrams in R with the Vennerable Package.

Thomma, B. P., Nürnberger, T., and Joosten, M. H. (2011). Of PAMPs and effectors: the blurred PTI-ETI dichotomy. Plant Cell 23, 4–15. doi: 10.1105/tpc.110.082602

PubMed Abstract | CrossRef Full Text | Google Scholar

van Damme, M., Huibers, R. P., Elberse, J., and van den Ackerveken, G. (2008). Arabidopsis DMR6 encodes a putative 2OG-Fe (II) oxygenase that is defense-associated but required for susceptibility to downy mildew. Plant J. 54, 785–793. doi: 10.1111/j.1365-313X.2008.03427.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van Damme, M., Zeilmaker, T., Elberse, J., Andel, A., de Sain-van der Velden, M., and van den Ackerveken, G. (2009). Downy mildew resistance in Arabidopsis by mutation of HOMOSERINE KINASE. Plant Cell 21, 2179–2189. doi: 10.1105/tpc.109.066811

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Burg, H. A., Tsitsigiannis, D. I., Rowland, O., Lo, J., Rallapalli, G., Maclean, D., et al. (2008). The F-box protein ACRE189/ACIF1 regulates cell death and defense responses activated during pathogen recognition in tobacco and tomato. Plant Cell 20, 697–719. doi: 10.1105/tpc.107.056978

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Vossen, E., Sikkema, A., Hekkert, B., Gros, J., Stevens, P., Muskens, M., et al. (2003). An ancient R gene from the wild potato species Solanum bulbocastanum confers broad-spectrum resistance to Phytophthora infestans in cultivated potato and tomato. Plant J. 36, 867–882. doi: 10.1046/j.1365-313X.2003.01934.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Vossen, E. A. G., Gros, J., Sikkema, A., Muskens, M., Wouters, D., Wolters, P., et al. (2005). The Rpi-blb2 gene from Solanum bulbocastanum is an Mi-1 gene homolog conferring broad-spectrum late blight resistance in potato. Plant J. 44, 208–222. doi: 10.1111/j.1365-313X.2005.02527.x

PubMed Abstract | CrossRef Full Text | Google Scholar

van Ooijen, G., van den Burg, H. A., Cornelissen, B. J. C., and Takken, F. L. W. (2007). Structure and function of resistance proteins in solanaceous plants. Annu. Rev. Phytopathol. 45, 43–72. doi: 10.1146/annurev.phyto.45.062806.094430

PubMed Abstract | CrossRef Full Text | Google Scholar

van Poppel, P. M., Guo, J., van de Vondervoort, P. J., Jung, M. W., Birch, P. R., Whisson, S. C., et al. (2008). The Phytophthora infestans avirulence gene Avr4 encodes an RXLR- dEER effector. Mol. Plant Microbe Interact. 21, 1460–1470. doi: 10.1094/MPMI-21-11-1460

PubMed Abstract | CrossRef Full Text | Google Scholar

Vleeshouwers, V. G. A. A., Raffaele, S., Vossen, J. H., Champouret, N., Oliva, R., Segretin, M. E., et al. (2011). Understanding and exploiting late blight resistance in the age of effectors. Annu. Rev. Phytopathol. 49, 507–531. doi: 10.1146/annurev-phyto-072910-095326

PubMed Abstract | CrossRef Full Text | Google Scholar

Vleeshouwers, V. G. A. A., Rietman, H., Krenek, P., Champouret, N., Young, C., Oh, S. K., et al. (2008). Effector genomics accelerates discovery and functional profiling of potato disease resistance and Phytophthora infestans Avirulence genes. PLoS ONE 3:e2875. doi: 10.1371/journal.pone.0002875

PubMed Abstract | CrossRef Full Text | Google Scholar

Vleeshouwers, V. G. A. A., van Dooijeweert, W., Govers, F., Kamoun, S., and Colon, L. T. (2000). The hypersensitive response is associated with host and nonhost resistance to Phytophthora infestans. Planta 210, 853–864. doi: 10.1007/s004250050690

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Cheng, Y. F., Deng, C., Ma, Y., Wang, Z. W., Chen, X. H., et al. (2014). Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis. BMC Genomics 15:412. doi: 10.1186/1471-2164-15-412

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Li, F., Qin, Y., Bo, X., Wu, Y., and Wang, S. (2010). GOSemSim: an R package for measuring semantic similarity among GO terms and gene products. Bioinformatics 26, 976–978. doi: 10.1093/bioinformatics/btq064

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Q., and Wang, X.-J. (2008). GOEAST: a web-based software toolkit for gene ontology enrichment analysis. Nucleic Acids Res. 36, W358–W363. doi: 10.1093/nar/gkn276

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Z., Nonomura, T., Appiano, M., Pavan, S., Matsuda, Y., Toyoda, H., et al. (2013). Loss of function in Mlo orthologs reduces susceptibility of pepper and tomato to powdery mildew disease caused by Leveillula taurica. PLoS ONE 8:e70723. doi: 10.1371/journal.pone.0070723

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, S., Li, Y., Vossen, J. H., Visser, R. G., and Jacobsen, E. (2012). Functional stacking of three resistance genes against Phytophthora infestans in potato. Transgenic Res. 21, 89–99. doi: 10.1007/s11248-011-9510-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Solanum, Phytophthora infestans, RNA-seq, R-genes, potato, resistance

Citation: Frades I, Abreha KB, Proux-Wéra E, Lankinen Å, Andreasson E and Alexandersson E (2015) A novel workflow correlating RNA-seq data to Phythophthora infestans resistance levels in wild Solanum species and potato clones. Front. Plant Sci. 6:718. doi: 10.3389/fpls.2015.00718

Received: 13 February 2015; Accepted: 27 August 2015;
Published: 17 September 2015.

Edited by:

Brigitte Mauch-Mani, Université de Neuchâtel, Switzerland

Reviewed by:

Weixing Shan, Northwest A&F University, China
Dennis Halterman, United States Department of Agriculture/Agricultural Research Service, USA

Copyright © 2015 Frades, Abreha, Proux-Wéra, Lankinen, Andreasson and Alexandersson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Erik Alexandersson, Department of Plant Protection Biology, Swedish University of Agricultural Sciences, Sundsvägen 14, SE-230 53 Alnarp, Sweden,

Present Address: Estelle Proux-Wéra, SciLifeLab, Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden

Joint first authors.