- Research article
- Open Access
Insights into organ-specific pathogen defense responses in plants: RNA-seq analysis of potato tuber-Phytophthora infestans interactions
BMC Genomics volume 14, Article number: 340 (2013)
The late blight pathogen Phytophthora infestans can attack both potato foliage and tubers. Although interaction transcriptome dynamics between potato foliage and various pathogens have been reported, no transcriptome study has focused specifically upon how potato tubers respond to pathogen infection. When inoculated with P. infestans, tubers of nontransformed ‘Russet Burbank’ (WT) potato develop late blight disease while those of transgenic ‘Russet Burbank’ line SP2211 (+RB), which expresses the potato late blight resistance gene RB (Rpi-blb1), do not. We compared transcriptome responses to P. infestans inoculation in tubers of these two lines.
We demonstrated the practicality of RNA-seq to study tetraploid potato and present the first RNA-seq study of potato tuber diseases. A total of 483 million paired end Illumina RNA-seq reads were generated, representing the transcription of around 30,000 potato genes. Differentially expressed genes, gene groups and ontology bins that exhibited differences between the WT and +RB lines were identified. P. infestans transcripts, including those of known effectors, were also identified.
Faster and stronger activation of defense related genes, gene groups and ontology bins correlate with successful tuber resistance against P. infestans. Our results suggest that the hypersensitive response is likely a general form of resistance against the hemibiotrophic P. infestans—even in potato tubers, organs that develop below ground.
Cultivated potato is the world’s third most important human food crop and the number one non-grain food commodity (FAOSTAT 2010). Unfortunately, potato is also host of a broad range of pathogens . Late blight disease is caused by Phytophthora infestans, resulted in the Irish Potato Famine in the 1840s, and still today results in multi-billion dollar losses worldwide annually. P. infestans is a notorious plant destroyer with the capacity to attack both potato foliage and tubers. Foliage resistance against late blight does not guarantee tuber resistance—contrasting disease resistance phenotypes can be evident in comparing foliage and tubers from a single genotype .
Gene RB (Rpi-blb1) [3, 4], a disease resistance gene cloned from a wild potato species, confers broad-spectrum foliar resistance against all major late blight pathogen isolates. The gene has been introduced into several potato cultivars using transgenic approaches . Transgenic lines have been or are currently being tested for eventual commercial release in Europe, India, Bangladesh, US and other places. Halterman et al.  reported that the four foliar late blight resistant +RB transgenic potato lines examined in their study lacked statistically significant tuber blight resistance. However, in a more extensive survey of +RB transgenic potato lines, we identified two transgenic lines with unusually high RB transcript levels that were late blight resistant in both foliage and tuber . Thus, RB has the potential to function in the tuber and the RB-potato tuber-P. infestans interaction provides a tractable system to study how potato tubers defend against plant pathogens.
Next-generation sequencing (NGS) technologies are fast evolving and are transforming biology research [7, 8]. Genome sequences of potato and P. infestans have been published [9, 10], making sequencing-based transcriptome studies (RNA-seq) more accessible to potato late blight researchers. RNA-seq is a relatively new approach towards study of the transcriptome [11, 12]. To our knowledge, there has been no RNA-seq study focused on potato-microbe interactions. Indeed, one NGS-based transcriptome study focused on the interaction between potato and P. infestans has been published . However, that study utilized the DeepSAGE method, not RNA-seq. DeepSAGE differs substantially from RNA-seq and researchers in that study relied heavily on assembled tags; the newly available genome sequence data were mostly not utilized. Furthermore, that and all other published studies focused on potato foliage transcriptome dynamics [13–15]; no previous study has reported transcriptome dynamics of potato tubers in response to pathogen attack.
In this study, we employed RNA-seq to study the transcriptome dynamics of potato tuber- P. infestans interactions in compatible and incompatible potato genotypes. We employed genome-wide sequence data from both potato and P. infestans in our analyses. Differentially expressed genes and ontology bins were identified that distinguish compatible and incompatible interactions. Transcripts of P. infestans, including those from candidate effectors, were also identified. Our study has important implications for potato R gene deployment and contributes to scientific understanding of organ-specific defense regulation in plants.
The +RBtransgenic line exhibits enhanced tuber late blight resistance
In previous field-based evaluations, the +RB line examined in this study, SP2211, was ranked as the third most foliar late blight resistant of 57 transgenic lines tested and was rated as “Resistant” to foliar late blight . In contrast, the WT line (nontransformed ‘Russet Burbank’) was rated as “Susceptible” to foliar late blight . Here, in replicated whole tuber assays performed six weeks after harvest, the +RB line showed no tuber late blight disease after normalization to the water-inoculated controls, whereas the WT line showed clear tuber late blight disease development (Figure 1). Results from two sets of whole tuber experiments were in agreement (32 of 33 P. infestans inoculations of WT tubers resulted in disease, none of the 60 inoculation sites in +RB tubers developed late blight disease). Thus, while WT is susceptible to P. infestans infection in both foliage and tuber, the RB transgene renders the +RB line SP2211 resistant to P. infestans infection in both foliage and tuber. RB gene transcription in the tubers of the +RB line was unaltered by inoculation with P. infestans at 0, 24, and 48 hpi (p>0.05).
RNA-seq reads aligned well with the potato reference genome sequence
A total of 36 RNA samples, collected from three biological replicates of P. infestans- or water-inoculated tuber tissues of two potato genotypes (WT and +RB) at three time points post inoculation, were subjected to RNA-seq. Approximately 483 million paired end reads were generated, yielding an average of 13.4 million paired end reads per sample. A total of 436.3 million Illumina reads (90.3%) passed quality filtering. The majority of reads (380.2 million or 78.7% of all reads) could be mapped uniquely to one location within the doubled monoploid (DM) potato reference genome sequence  (Figure 2). An additional 16.6 million reads (3.5%) were mapped to multiple locations within the reference genome sequence (Figure 2).
We detected transcription of 29,319 potato genes based on cufflinks FPKM (fragments per kilobase of exon per million mapped reads) information and gene models reported by the Potato Genome Sequencing Consortium (PGSC) . Wang et al.  demonstrated that a data set of 10 million RNA-seq reads represented about 80% of annotated chicken genes, concluding that RNA-seq is a viable alternative to the microarray for transcriptome study. In this study, we detected more than 75% of PGSC potato gene models—a percentage similar to that reported by Wang et al. . This scale of gene detection is higher than reported for any previous potato-P. infestans transcriptome study. Importantly, 19.8% of RNA-seq reads that passed quality filters were mapped to regions outside of the current potato reference genome gene models (or predicted exons), suggesting further refinement of gene annotations associated with the potato reference genome sequence may be warranted. Nonetheless, qRT-PCR results correlated well with RNA-seq data with an average Pearson correlation coefficient of r=0.86 (Additional file 1). Together, our results suggest that RNA-seq and qPCR approaches can be cross-validated, confirming that current potato gene models  are appropriate for functional genomics studies of potato tubers.
Overall transcriptome dynamics
A principal component analysis (PCA) of log2 transformed FPKM values for 29,319 genes from the 36 RNA-seq samples is shown in Figure 3. Together, PC1 and PC2 explained >40% of the total variance of this dataset. Samples collected at 0 hpi were distinct from all other samples, with PC1 providing clear separation of 0 hpi samples from 24 hpi and 48 hpi samples. PC1 also differentiates water- and P. infestans-inoculated samples at 48 hpi. Overall, the distribution of +RB samples is similar to that of corresponding WT samples. These results indicate that time and treatments (P. infestans- vs. water-inoculation) play a greater role than genotype in defining overall transcriptome dynamics. This results is expected given that the +RB line was created by transformation of the WT line with the RB transgene.
Differentially expressed (DE) genes
We also identified 2,531 genes that are differentially expressed (DE) during at least one of the between time transitions in P. infestans-inoculated samples (Additional file 2). Transcriptome dynamics of both genotypes in P. infestans- compared to water-inoculated samples during transitions from 0 hpi to 24 hpi and from 24 hpi to 48 hpi were visualized using MA plots (hybridization intensity plotted against the fold change in expression; Figure 4). Interestingly, at later stages of infection (24 hpi – 48 hpi), WT tubers exhibited mostly up-regulation of genes, whereas the +RB line displayed approximately equal numbers of up- and down-regulated genes (Figure 4), potentially reflecting contrasting rates of pathogen proliferation (as indicated by total pathogen reads counts, see results sections below) and corresponding host response patterns.
Transcriptome dynamics of between line comparisons were also visualized using MA plots (Additional file 3). Consistent with PCA which shows that genotype has a lesser role (compared to time or treatment) in defining overall transcriptome dynamics, fewer genes were determined to be DE during between line comparisons, despite a trend of increasing numbers of DE genes over time (Additional file 3). This suggests that both the WT and the +RB line use that same set of genes to respond to P. infestans attack but that there are significant differences in the temporal regulation of these genes.
Next we employed hierarchical clustering to visualize the DE genes identified through between time point comparisons. Despite a preponderance of shared DE gene regulation patterns, groups of genes that distinguish compatible and incompatible interactions were still identified (Additional files 3 and 4). These include WRKY (PGSC0003DMG400008188) and ethylene response (PGSC0003DMG400025989, ERF1) transcription factors. The ERF1 transcription factor gene is a likely key element in the integration of JA (jasmonic acid) and ET (ethylene) signals for the regulation of defense responses .
The +RBtransgenic line has faster and stronger induction of DE genes
A total of 1,102 DE genes were identified when water- and P. infestans-inoculated samples were compared (between treatment comparisons) at the same time points (24 hpi or 48 hpi). Of these, 959 were DE only in the +RB line, 34 were DE only in the WT line, and 109 were DE in both lines. We identified two representative hierarchical clusters, comprising a total of 31 genes, showing interesting regulation patterns at 48 hpi (Figure 5). In WT, these 31 genes were suppressed (compared to water-inoculated); in +RB, these 31 genes were induced (compared to water-inoculated). These genes include ethylene response factors (ERF2), cellulose synthase, chitinase, and elicitor inducible cytochrome P450. Two genes specifically associated with the hypersensitive response (HR), PGSC0003DMG400025335 [Hypersensitive-induced reaction protein (HIR); Figure 5] and PGSC0003DMG400027473 (Hypersensitive response assisting protein) were significantly induced in the +RB line but not in the WT line at 48 hpi. It is important to highlight that the +RB line displayed 903 DE genes at 48 hpi, accounting for about 82% of the total DE genes identified during between treatment comparisons (Figure 6). These results suggest that the +RB line responds to P. infestans attack by rapid differential regulation of large sets of genes, whereas the WT line is slower to respond (e.g., 48 hpi).
Ontology bins and gene groups distinguishing compatible and incompatible interactions
Mapman ontology enrichment analysis of DE genes (Figure 6) revealed that WRKY transcription factors are up-regulated in +RB but not in WT at 48 hpi. Strikingly, 13 of the 14 DE (during water- vs. P. infestans-inoculated comparisons) WRKY transcription factor genes show significant induction in the +RB line at 48 hpi (Additional file 5). In contrast, none of the 14 WRKY transcription factor genes showed statistically significant (FDR<0.001) induction in the WT line at 48 hpi. Homologs of several of these WRKY transcription factors have been previously reported to be involved in defense against Phytophthora spp. and other pathogens [18, 19]. We hypothesize that DE WRKY transcription factors identified in this study are downstream regulatory components of RB-mediated defenses in potato tubers.
Very importantly, we discovered that when WT and +RB were compared directly, the +RB line displays stronger transcription of genes within receptor kinase and defense related bins, even at 0 hpi (Figure 7 and Additional file 6). It is especially relevant to note that many of the bins commonly associated with plant defense, including receptor kinases and PR (pathogenesis-related) protein genes are more highly expressed in the +RB line. Interestingly, WT actually has much higher transcription of similar defense-related bins at 24 hpi (but not 48 hpi) if an indirect comparison (adjusting FPKM values of P. infestans-inoculated samples by FPKM values of water-inoculated samples of the same genotype) was made, indicating that WT and +RB probably share large sets of defense components. Thus, WT responds to P. infestans by up-regulation of sets of defense response genes that are constitutively transcribed at higher levels in the +RB line even in the absence of the pathogen (i.e., 0 hpi). But slower activation of defense related genes in WT is insufficient to prevent disease development, resulting in a tuber blight susceptible phenotype. These results were corroborated by an independent pilot RNA-seq study (Illumina GAIIx single end reads) using field grown tubers of the same genotypes from a different year (Additional file 7). Collectively, these findings suggest that faster and stronger expression of defense related genes plays a role in enhanced disease resistance in the +RB line.
RNA-seq reads mapping to Phytophthora infestansreference transcripts
A total of 17,353 paired end sequence reads were mapped to the P. infestans reference transcript set , representing the transcription of over 4,600 genes. Ninety-six P. infestans transcripts are predicted to encode RxLR effectors .
Throughout the infection time course, a trend of increasing mappable P. infestans reads is evident for WT samples. This increase of mappable reads was not seen in the +RB line from 24 hpi to 48 hpi (Additional file 8) and a majority (>82%) of all mapped reads originated from infected WT potato tuber samples collected at 48 hpi. This is in agreement with tuber late blight disease development observed in WT but not +RB tubers. In total, a majority of all P. infestans reads (>95%) were from infected WT potato tuber samples; less than 5% of the reads were from infected +RB potato tuber samples (Additional file 9).
Class I and II ipiO proteins, members of the RxLR family of putative effectors, induce avirulence responses in potato lines containing RB[20–22]. After inoculation (24, 48 hpi), class I ipiO transcripts were detected (though in low abundance) in five of six tuber samples collected from the P. infestans-inoculated WT line but from no samples derived from the +RB line (Additional file 10). This is probably due to suppression of P. infestans proliferation in the +RB line.
Faster and stronger activation of defense related genes, gene groups and ontology bins correlates with successful tuber defense
Direct comparison between WT and +RB lines at 0 hpi revealed a total of 11 DE genes (Additional file 3). Seven of the eleven DE genes were more highly transcribed in +RB. These include Carbonic anhydrase (PGSC0003DMG400006956), Aquaglyceroporin (PGSC0003DMG400009604), Cytochrome P450 (PGSC0003DMG400015185), Malic enzyme (PGSC0003DMG401026923), and 1-aminocyclopropane-1-carboxylate oxidase (ACO) (PGSC0003DMG401026923). Interestingly, these genes or homologs of these genes were previously reported to be involved in basal or pathogen-induced defense responses in plants [14, 23, 24]. The four remaining DE genes [protoporphyrinogen oxidase (PGSC0003DMG400001117), phenylcoumaran benzylic ether reductase (PGSC0003DMG400003691), DNA binding protein (PGSC0003DMG400003902), Major Latex (PGSC0003DMG400008811)] are transcribed at higher levels in WT than in the +RB line. Each of these genes has also been associated with plant defense [25–28]. Together, our results demonstrate that the +RB line differs from WT in the constitutive transcriptional regulation of defense related genes.
GO enrichment analysis of DE genes identified during between treatment comparisons (Figure 6) suggests that +RB has faster response to pathogen attack, as evidenced by higher and statistically significant induction of ontology bins commonly associated with plant defense (e.g., ET, WRKY, Signaling receptor kinases) at 48 hpi. One notable example is WRKY transcription factors (Additional file 5), with 13 out of 14 WRKY transcription factors significantly induced in +RB. These same transcription factors are also induced, but at nonsignificant levels, in the WT line, documenting that more potent and rapid activation of defense related genes, rather than novel resistance responses per se, might play a role in enhanced tuber resistance in the +RB line.
Furthermore, when a more holistic approach (using all available data points: 39,031 gene transcript values per sample) was employed to identify differentially regulated ontology bins, various defense related ontology bins were discovered (Figure 7). +RB potato tubers displayed faster and stronger transcription of genes within bins commonly associated with plant defense such as receptor kinases and pathogenesis-related (PR) genes (Figure 7). Together, these results suggest that the +RB line is tuber late blight resistant due to faster and stronger activation of defense related components.
Cao et al.  found that increased transcription of the rice Xa3 resistance gene correlated with enhanced expression of defense-responsive genes and an enlarged resistance spectrum to Xanthomonas oryzae pv. oryzae (Xoo). Previously we demonstrated that higher transgene RB copy numbers correlated with both higher RB transcript levels and enhanced late blight resistance in the foliage . The +RB line employed in this study, SP2211, displayed the highest tuber transgene transcription levels among 11 transgenic lines examined . Thus, the unusually high levels of RB gene transcription in the tubers of SP2211 correlates with faster and stronger activation of defense related components (Figures 6, 7). The faster and stronger up-regulation of these defense-related components very likely contributed to the observed successful defense against P. infestans (Figure 1) in this line. One interesting observation is that WT actually has much higher transcription of similar defense-related bins at 24 hpi (but not 48 hpi) in indirect (adjusting FPKM values of P. infestans-inoculated samples by FPKM values of water-inoculated samples of the same genotype) comparisons. Collectively, our data support a model of tuber blight resistance resulting as a function of faster and stronger expression of defense related genes due to high RB transcript levels.
HR is likely a general form of resistance to P. infestans, even in potato tubers
Prevailing scientific thought suggests that the phyllosphere frequently encounters biotrophic or hemibiotrophic pathogens, against which the HR is effective. In contrast, roots are more likely to encounter necrotrophic microbes that would theoretically benefit from the HR . Consistently, working in Arabidopsis, Hermanns et al.  showed that HR is present in above ground incompatible host-pathogen interactions, but absent during the same host-pathogen interaction in the roots—despite R gene transcription in both leaves and roots. These authors concluded that R gene function is modulated in an organ-specific but mechanistically poorly defined manner for the purpose of suppressing HR in the roots. Observations of differential late blight resistance levels in potato foliage and tubers within a single genotype  suggest that R gene function in potato may also be modulated in an organ-specific manner. But the potato tuber, although it develops below ground, is a modified stem, not a root. Previous studies [32, 33] have suggested that HR plays a role in tuber defense. The identification of SP2211, a +RB transgenic line with both foliar and tuber late blight resistance provided opportunity to examine the role of HR in defending non-root organs against pathogen attack.
P. infestans is a hemibiotrophic pathogen and Kamoun et al.  argued that HR is likely a general form of resistance against Phytophthora pathogens. In agreement, Chen and Halterman  demonstrated that RB triggers an HR in potato foliage challenged with P. infestans. In examining RB-mediated tuber responses to P. infestans, we identified DE genes associated with the HR (potato HIR and hypersensitive response assisting genes) that are highly up-regulated in pathogen-inoculated tuber samples of SP2211 compared to water-inoculated tuber samples at 48 hpi. Previous studies have shown that the pepper homolog of HIR is a positive regulator of hypersensitive cell death [36, 37]. Qi et al.  also demonstrated that the NBS-LRR protein RPS2 forms complexes with AtHIR proteins in Arabidopsis and tobacco plants, providing mechanistic insight into R protein function and the HR. We hypothesize that RB and HIR might form similar complexes in potato, a topic that warrants further experimentation. Rapid production of reactive oxygen species by plants is a hallmark of pathogen recognition and correlates with the HR. Working with disks cut from potato tubers that carry the R1 late blight resistance gene, Doke  demonstrated that incompatible but not compatible races of P. infestans triggered production of reactive oxygen species. In our study, a gene (PGSC0003DMG400024754) encoding respiratory burst oxidase homolog protein B (NADPH oxidase RBOHB) was induced (based on between treatment comparisons) 6.6 fold in the +RB line but not in the WT line 48 hours after P. infestans-inoculation. This gene is known to be involved in the massive phase II oxidative burst induced in potato by pathogen infection  and offers further support that RB-mediated resistance to tuber late blight likely entails an HR or HR-like phenomenon.
Interestingly, HR is commonly associated with SA (salicylic acid) mediated responses. But our ontology bin analysis didn’t reveal SA metabolism as a differentiating factor for the tuber blight resistant and susceptible lines. Instead, our study reveals that the ET (ethylene) bin is predominantly associated with the incompatible interaction. Leon-Reyes et al.  reported that the antagonistic relationship between SA and JA  only exists when there is no strong production of ET. Nunez-Pastrana et al.  reported that ethylene but neither SA nor Methyl JA (jasmonic acid) induces resistance response against Phytophthora capsici in pepper. Extensive hormonal crosstalk  during plant responses to pathogens is a topic of extensive ongoing research and future results will likely yield new insights that will allow fine-tuning of our understanding of potato tuber-pathogen interactions and organ-specific defense responses in plants.
We presented the first RNA-seq (or transcriptome dynamics) study focused on potato tuber responses to pathogen attack. We mapped approximately 400 million RNA-seq reads onto the recently published potato reference genome sequence, documenting the utility of RNA-seq for biological study of tetraploid cultivated potato. We identified sets of DE genes that distinguish resistant and susceptible lines. Our data suggest that potent regulation of defense genes and gene groups or ontology bins (e.g., HIR and WRKY transcription factors) plays a role in RB-mediated tuber defense. In particular, faster and stronger expression of defense related genes correlates with enhanced tuber blight resistance in the +RB transgenic line. In agreement with Kamoun et al.  and Doke , our data suggest that HR is likely a general form of resistance against late blight even in the potato tuber, an organ that develops below ground.
Plant material and RNA preparation
The tuber late blight susceptible nontransformed ‘Russet Burbank’ (WT, provided by Dr. Carl Rosen, UMN) and SP2211, a tuber late blight resistant transformed ‘Russet Burbank’ line carrying the RB transgene (+RB) were examined in this study. Tubers were produced under standard cultural practices at the University of Minnesota Sand Plain Research Farm (Becker, MN). Tubers of each line were harvested and held for three days at room temperature before storage for six weeks at 11–13 degrees Celsius. P. infestans US8 isolate US940480  was maintained on Rye A medium . Prior to inoculations, sporangia were harvested from plates by physical scraping into distilled water. Inoculum was adjusted to 6.75*104 sporangia/ml. Prepared inoculum was incubated for 1 hour at 4 degrees Celsius and then at room temperature for 30 minutes. Six week old WT and +RB tubers were inoculated using a modified whole tuber assay , as detailed below.
For RNA extractions and RNA-Seq, six tubers of each genotype were randomly selected, washed with deionized water, and allowed to air dry at room temperature for 24 hours. Tubers were wounded [0.2 cm*0.3 cm (depth*diameter)] at six sites spaced uniformly (~3 cm apart) across the tuber surface. Each wound was inoculated with 10 μl sporangial suspension (3 tubers per genotype) or water (mock treatment; 3 tubers per genotype). Inoculated tubers were stored in tightly sealed dark boxes under high (~95%) humidity and at room temperature for 72 hours (during which time all RNA-seq samples were collected as specified below) and then moved to 11–13 Celsius for disease development and phenotyping at 11 days post inoculation. From each of three replicate tubers per genotype, tuber tissue (0.7 cm*(0.5-0.8)cm) was collected using a cork borer from an inoculation site at time point 0 (pre-inoculation) and at 6, 12, 24, and 48 hours post inoculation (hpi). The sampled tissue includes cells from the periderm, cortex and medulla layers. The sixth inoculation site on each tuber was left intact and the tubers were phenotyped for late blight disease development 11 days later by multiplying measured length, width, and depth of disease lesions. Collected tissue samples (one tuber core per replicate tuber for each time point) were immediately frozen in liquid nitrogen and stored at −80 degrees Celsius. To reduce possible experimental variation, tissues collected for RNA-seq originated from a common experiment involving a single P. infestans inoculum preparation or water. In this study, each tuber is considered a biological replicate of the corresponding genotype by treatment combination. In total, 36 tissue samples from the two plant genotypes (WT and +RB) x three time points (0, 24, and 48 hpi) x two inocula (P. infestans or water) x three replicates were employed for RNA extraction. [Note: We collected tissue at 0, 6, 12, 24, and 48 hpi, but only 0, 24, and 48 hpi samples were subjected to RNA-seq]. Total RNA was extracted from frozen tissue using the SV Total RNA Isolation System (Promega Corporation, Madison, WI) according to manufacturer’s instructions. The quantity and quality of RNA samples were assessed using a Nanodrop 1000 machine (Thermo Fisher Scientific Inc., Wilmington, DE). High quality total RNA (5ug, 100 ng/μl) samples were sent to the University of Minnesota BioMedical Genomics Center (BMGC) for RNA-seq library preparation using the TruSeq SBS Kit (50 Cycles) (Illumina Inc., San Diego, CA) and paired end sequencing using an Illumina Hi-Seq 2000 machine (Illumina).
A separate, additional set of tubers (nine for WT and eighteen for +RB) were simultaneously subjected to the whole tuber assay (using the same inoculum and protocol described above), but without tissue collection, for the purposes of disease phenotyping.
RNA-seq reads mapping and DE genes clustering
RNA-seq reads were quality filtered using SolexaQA packages  with default parameters and a length filter of greater than 27 bp for both ends of each read pair. Sequence data have been submitted to the NCBI Sequence Read Archive [accession number SRP022916]. Quality filtered RNA-seq reads were analyzed using the “Tuxedo Suite” software packages : Bowtie  v0.12.7, Tophat  v1.3.2, Cufflinks  v1.1.0 and a reference genome of potato . We suppressed Tophat  from identifying novel junctions. A total of eleven pair-wise comparisons [between time (4), between line (3), between treatment (4)] were made using Cuffdiff of the Cufflinks software packages. Differentially expressed (DE) genes are those genes that showed significantly (FDR adjusted p-value <0.001) different transcript levels among comparisons. False discovery rate correction (FDR) was done using the Benjamini Hochberg method .
Quality filtered RNA-seq reads were also mapped to reference transcript sequences of P. infestans and different classes (see Champouret et al.  for details of classification) of ipiO (potential cognate avr effector of RB) gene sequences. The mapping was done using Bowtie  v0.12.7. Cufflinks software suites were not used in this analysis as underlying normalization assumptions were not met. Potato gene expression fold change values were used to cluster and partition genes into groups using hierarchical clustering and the complete linkage method in Cluster 3.0  and visualized in Treeview software  using adjusted pixel settings (threshold 1.0).
The UNIX command line options for Tophat and Cuffdiff were as follows: “tophat –G genes.gff –o sample.out --no-novel-juncs –r 100 genome.fasta read.left.fq read.right.fq; cuffdiff –N –u –O cuffdiff.out genes.gff 1.bam,2.bam,3.bam 4.bam,5.bam,6.bam”. The UNIX command line options for Bowtie were as follows: “Bowtie -a -v 2 --best –M 1 --fr ref_bowtie_build −1 read1.fq −2 read2.fq”. Custom perl scripts were used to parse the mapping/alignment results into tabular formats. R statistical software  was used to generate various plots.
Functional assignment and MapMan analysis of potato genes
To assign potato genes into functional categories (bins), we adopted MapMan ontology . The Mercator (Max Planck Institute, Germany, Dr. Mark Lohse, personal communication) annotation pipeline was used to assign potato genes into functional bins by searching a variety of reference databases. The mapping file was manually curated to remove suspicious assignments. Water- and P. infestans-inoculated samples were compared using Fisher’s exact test to identify ontology bins enriched with significantly up or down regulated DE (FDR <0.001) genes. WT and +RB lines were compared using the Wilcoxon rank-sum test to identify functional bins that behaved differently. Resulting p-values from Fisher’s exact tests or Wilcoxon rank-sum tests were subjected to multiple test correction using the Benjamini Hochberg method .
Real-time quantitative RT-PCR
For determining the transcript levels of the transgene RB, we used a previously established qRT-PCR protocol [5, 55]. For RNA-seq validation, Primer Express 3.0 (Applied Biosystems, Foster City, CA, USA) was used with default parameters to generate primer pairs for selected transcripts (Additional file 11). Total RNA (150 ng, 30 ng/ul) from diluted stocks of the same RNA that was subjected to RNA-seq was used in each reverse transcription reaction using the SuperScript III First Strand Synthesis Kit (Life Technologies Inc., Carlsband CA) according to manufacturer’s instructions. All qPCRs were performed using the Power SYBR Green Master Mix (Life Technologies) and an ABI 7500 Real Time System (Applied Biosystems). Target gene transcript levels were normalized to EF1α using 2-ΔCT values. Normalized gene expression levels were compared with RNA-seq FPKM values derived from Cufflinks.
Stevenson WR, Loria R, Frank GD, Weingartner DP: Compedium of potato diseases. 2001, Saint Paul, MN, USA: APS Press
Kirk WW, Felcher KJ, Douches DS, Niemira BA, Hammerschmidt R: Susceptibility of potato (Solanum tuberosum L.) foliage and tubers to the US8 genotype of Phytophthora infestans. Am J Potato Res. 2001, 78 (4): 319-322. 10.1007/BF02875697.
Song JQ, Bradeen JM, Naess SK, Raasch JA, Wielgus SM, Haberlach GT, Liu J, Kuang HH, Austin-Phillips S, Buell CR: Gene RB cloned from Solanum bulbocastanum confers broad spectrum resistance to potato late blight. Proc Natl Acad Sci USA. 2003, 100 (16): 9128-9133. 10.1073/pnas.1533501100.
van der Vossen E, Sikkema A, Hekkert BTL, Gros J, Stevens P, Muskens M, Wouters D, Pereira A, Stiekema W, Allefs S: 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. 2003, 36 (6): 867-882. 10.1046/j.1365-313X.2003.01934.x.
Bradeen JM, Lorizzo M, Mollov DS, Raasch J, Kramer LC, Millett BP, Austin-Phillips S, Jiang JM, Carputo D: Higher copy numbers of the potato RB transgene correspond to enhanced transcript and late blight resistance levels. Mol Plant Microbe Interact. 2009, 22 (4): 437-446. 10.1094/MPMI-22-4-0437.
Millett BP: Characterization of transcription and function of foliar blight resistance gene RB in foliage and tubers of transgenic cultivated potato at different physiological ages. 2008, University of Minnesota: PhD Thesis, 178-
Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92 (5): 255-264. 10.1016/j.ygeno.2008.07.001.
Metzker ML: Applications of next generation sequencing sequencing technologies - the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.
Xu X, Pan SK, Cheng SF, Zhang B, Mu DS, Ni PX, Zhang GY, Yang S, Li RQ, Wang J: Genome sequence and analysis of the tuber crop potato. Nature. 2011, 475 (7355): 189-U194. 10.1038/nature10158.
Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, Grabherr M, Kodira CD, Raffaele S, Torto-Alalibo T: Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature. 2009, 461 (7262): 393-398. 10.1038/nature08358.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.
Gyetvai G, Sonderkaer M, Gobel U, Basekow R, Ballvora A, Imhoff M, Kersten B, Nielsen KL, Gebhardt C: The transcriptome of compatible and incompatible interactions of potato (Solanum tuberosum) with Phytophthora infestans revealed by DeepSAGE analysis. PLoS One. 2012, 7 (2): e31526-10.1371/journal.pone.0031526.
Restrepo S, Myers KL, del Pozo O, Martin GB, Hart AL, Buell CR, Fry WE, Smart CD: Gene profiling of a compatible interaction between Phytophthora infestans and Solanum tuberosum suggests a role for carbonic anhydrase. Mol Plant Microbe Interact. 2005, 18 (9): 913-922. 10.1094/MPMI-18-0913.
Birch PRJ, Avrova AO, Duncan JM, Lyon GD, Toth RL: Isolation of potato genes that are induced during an early stage of the hypersensitive response to Phytophthora infestans. Mol Plant Microbe Interact. 1999, 12 (4): 356-361. 10.1094/MPMI.19184.108.40.2066.
Wang Y, Ghaffari N, Johnson CD, Braga-Neto UM, Wang H, Chen R, Zhou HJ: Evaluation of the coverage and depth of transcriptome by RNA-Seq in chickens. BMC Bioinforma. 2011, 12 (suppl 10): S5-10.1186/1471-2105-12-S10-S5.
Lorenzo O, Piqueras R, Sanchez-Serrano JJ, Solano R: Ethylene response factor1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell. 2003, 15 (1): 165-178. 10.1105/tpc.007468.
Nakazawa-Ueji YE, Nunez-Pastrana R, Souza-Perera RA, Santana-Buzzy N, Zuniga-Aguilar JJ: Mycelium homogenates from a virulent strain of Phytophthora capsici promote a defence-related response in cell suspensions from Capsicum chinense. Eur J Plant Pathol. 2010, 126 (3): 403-415. 10.1007/s10658-009-9544-x.
Katou S, Yoshioka H, Kawakita K, Rowland O, Jones JDG, Mori H, Doke N: Involvement of PPS3 phosphorylated by elicitor-responsive mitogen-activated protein kinases in the regulation of plant cell death. Plant Physiol. 2005, 139 (4): 1914-1926. 10.1104/pp.105.066795.
Pieterse CMJ, Vanwest P, Verbakel HM, Brasse P, Vandenbergvelthuis GCM, Govers F: Structure and genomic organization of the ipib and ipio gene clusters of Phytophthora-infestans. Gene. 1994, 138 (1–2): 67-77.
Vleeshouwers VGAA, Rietman H, Krenek P, Champouret N, Young C, Oh S-K, Wang M, Bouwmeester K, Vosman B, Visser RGF: Effector genomics accelerates discovery and functional profiling of potato disease resistance and Phytophthora infestans avirulence genes. PLoS One. 2008, 3 (8): e2875-10.1371/journal.pone.0002875.
Champouret N, Bouwmeester K, Rietman H, van der Lee T, Maliepaard C, Heupink A, van de Vondervoort PJI, Jacobsen E, Visser RGF, van der Vossen EAG: Phytophthora infestans isolates lacking class I ipiO variants are virulent on Rpi-blb1 potato. Mol Plant Microbe Interact. 2009, 22 (12): 1535-1545. 10.1094/MPMI-22-12-1535.
Yu MM, Shen L, Zhang AJ, Sheng JP: Methyl jasmonate-induced defense responses are associated with elevation of 1-aminocyclopropane-1-carboxylate oxidase in Lycopersicon esculentum fruit. J Plant Physiol. 2011, 168 (15): 1820-1827. 10.1016/j.jplph.2011.05.021.
Ishihama N, Yamada R, Yoshioka M, Katou S, Yoshioka H: Phosphorylation of the Nicotiana benthamiana WRKY8 transcription factor by MAPK functions in the defense response. Plant Cell. 2011, 23 (3): 1153-1170. 10.1105/tpc.110.081794.
Molina A, Volrath S, Guyer D, Maleck K, Ryals J, Ward E: Inhibition of protoporphyrinogen oxidase expression in arabidopsis causes a lesion-mimic phenotype that induces systemic acquired resistance. Plant J. 1999, 17 (6): 667-678. 10.1046/j.1365-313X.1999.00420.x.
Gang DR, Kasahara H, Xia ZQ, Vander Mijnsbrugge K, Bauw G, Boerjan W, Van Montagu M, Davin LB, Lewis NG: Evolution of plant defense mechanisms - relationships of phenylcoumaran benzylic ether reductases to pinoresinol-lariciresinol and isoflavone reductases. J Biol Chem. 1999, 274 (11): 7516-7527. 10.1074/jbc.274.11.7516.
Korfhage U, Trezzini GF, Meier I, Hahlbrock K, Somssich IE: Plant homeodomain protein involved in transcriptional regulation of a pathogen defense-related gene. Plant Cell. 1994, 6 (5): 695-708.
Konno K: Plant latex and other exudates as plant defense systems: roles of various defense chemicals and proteins contained therein. Phytochemistry. 2011, 72 (13): 1510-1530. 10.1016/j.phytochem.2011.02.016.
Cao YL, Ding XH, Cai M, Zhao J, Lin YJ, Li XH, Xu CG, Wang SP: Expression pattern of a rice disease resistance gene Xa3/Xa26 is differentially regulated by the genetic backgrounds and developmental stages that influence its function. Genetics. 2007, 177 (1): 523-533. 10.1534/genetics.107.075176.
Glazebrook J: Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005, 43: 205-227. 10.1146/annurev.phyto.43.040204.135923.
Hermanns M, Slusarenko AJ, Schlaich NL: Organ-specificity in a plant disease is determined independently of R gene signaling. Mol Plant Microbe Interact. 2003, 16 (9): 752-759. 10.1094/MPMI.2003.16.9.752.
Doke N: Involvement of superoxide anion generation in the hypersensitive response of potato - tuber tissues to infection with an incompatible race of Phytophthora infestans and to the hyphal wall components. Physiological Plant Pathology. 1983, 23 (3): 345-357. 10.1016/0048-4059(83)90019-X.
Tomiyama K: Further observation on the time requirement for hypersensitive cell death of potatoes infected by Phytophthora infestans and its relation to metabolic activity. Phytopathologische Zeitschrift. 1967, 58 (4): 367-378. 10.1111/j.1439-0434.1967.tb02321.x.
Kamoun S, Huitema E, Vleeshouwers V: Resistance to oomycetes: a general role for the hypersensitive response?. Trends Plant Sci. 1999, 4 (5): 196-200. 10.1016/S1360-1385(99)01404-1.
Chen Y, Halterman DA: Phenotypic characterization of potato late blight resistance mediated by the broad-spectrum resistance gene RB. Phytopathology. 2011, 101 (2): 263-270. 10.1094/PHYTO-04-10-0119.
Jung HW, Hwang BK: The leucine-rich repeat (LRR) protein, CaLRR1, interacts with the hypersensitive induced reaction (HIR) protein, CaHIR1, and suppresses cell death induced by the CaHIR1 protein. Mol Plant Pathol. 2007, 8 (4): 503-514. 10.1111/j.1364-3703.2007.00410.x.
Jung HW, Lim CW, Lee SC, Choi HW, Hwang CH, Hwang BK: Distinct roles of the pepper hypersensitive induced reaction protein gene CaHIR1 in disease and osmotic stress, as determined by comparative transcriptome and proteome analyses. Planta. 2008, 227 (2): 409-425.
Qi YP, Tsuda K, Nguyen LV, Wang X, Lin JS, Murphy AS, Glazebrook J, Thordal-Christensen H, Katagiri F: Physical association of arabidopsis hypersensitive induced reaction proteins (HIRs) with the immune receptor RPS2. J Biol Chem. 2011, 286 (36): 31297-31307. 10.1074/jbc.M110.211615.
Yoshioka H, Sugie K, Park HJ, Maeda H, Tsuda N, Kawakita K, Doke N: Induction of plant gp91 phox homolog by fungal cell wall, arachidonic acid, and salicylic acid in potato. Mol Plant Microbe Interact. 2001, 14 (6): 725-736. 10.1094/MPMI.2001.14.6.725.
Leon-Reyes A, Spoel SH, De Lange ES, Abe H, Kobayashi M, Tsuda S, Millenaar FF, Welschen RAM, Ritsema T, Pieterse CMJ: Ethylene modulates the role of nonexpressor of pathogenesis related Genes1 in cross talk between salicylate and jasmonate signaling. Plant Physiol. 2009, 149 (4): 1797-1809. 10.1104/pp.108.133926.
Niki T, Mitsuhara I, Seo S, Ohtsubo N, Ohashi Y: Antagonistic effect of salicylic acid and jasmonic acid on the expression of pathogenesis-related (PR) protein genes in wounded mature tobacco leaves. Plant Cell Physiol. 1998, 39 (5): 500-507. 10.1093/oxfordjournals.pcp.a029397.
Nunez-Pastrana R, Arcos-Ortega GF, Souza-Perera RA, Sanchez-Borges CA, Nakazawa-Ueji YE, Garcia-Villalobos FJ, Guzman-Antonio AA, Zuniga-Aguilar JJ: Ethylene, but not salicylic acid or methyl jasmonate, induces a resistance response against Phytophthora capsici in habanero pepper. Eur J Plant Pathol. 2011, 131 (4): 669-683. 10.1007/s10658-011-9841-z.
Robert-Seilaniantz A, Grant M, Jones JDG: Hormone crosstalk in plant disease and defense: more than just jasmonate-salicylate antagonism. Annual review of phytopathology, Vol 49. Edited by: VanAlfen NK, Bruening G, Leach JE. 2011, Palo Alto: Annual Reviews, 317-343.
Tumwine J, Frinking HD, Jeger MJ: Isolation techniques and cultural media for Phytophthora infestans from tomatoes. Mycologist. 2000, 14 (3): 137-139. 10.1016/S0269-915X(00)80096-8.
Cox M, Peterson D, Biggs P: SolexaQA: at-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinforma. 2010, 11 (1): 485-10.1186/1471-2105-11-485.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012, 7 (3): 562-578. 10.1038/nprot.2012.016.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-U174. 10.1038/nbt.1621.
Benjamini Y, Hochberg Y: Controlling false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Ser B-Methodol. 1995, 57 (1): 289-300.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Saldanha AJ: Java treeview-extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248. 10.1093/bioinformatics/bth349.
R Development Core Team: R: a language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, URL http://www.R-project.org/, 3-900051-07-0
Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.
Millett BP, Bradeen JM: Development of allele-specific PCR and RT-PCR assays for clustered resistance genes using a potato late blight resistance transgene as a model. Theor Appl Genet. 2007, 114 (3): 501-513. 10.1007/s00122-006-0449-1.
Nicot N, Hausman JF, Hoffmann L, Evers D: Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot. 2005, 56 (421): 2907-2914. 10.1093/jxb/eri285.
The work was funded by USDA NIFA, EPA, MN area II Potato Research and Promotion Council and a Thesis Research Grant to LLG from the Graduate School, University of Minnesota. The authors thank Dr. Mark Lohse of Max Planck Institute (Germany) for generating a draft PGSC mapping file for MapMan analysis in potato. The authors thank Dr. Fumiaki Katagiri for thoughtful comments on statistical analysis and draft versions of our manuscript. Computing resources from the Minnesota Supercomputing Institute at the University of Minnesota are greatly appreciated.
The authors declare that they have no competing interests.
LLG was responsible of all bench work and informatics analyses of RNA-seq reads filtering, mapping, DE genes detection, and ontology regulation analysis. ZJT provided some informatics support, especially in relation to a pilot RNA-seq project. BPM identified transgenic line SP2211 (+RB), pioneered tuber resistance phenotyping, and developed molecular assays for transgene transcription measurement. JMB provided the conceptual impetus for this work, suggestions on all aspects of experimental design and data analysis, and secured funding for this project. LLG wrote this manuscript and ZJT, BPM, and JMB provided comments and suggestions to improve the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: RNA-seq FPKM and qPCR correlation. First column is potato tuber sample/IDs. Remaining columns indicate qPCR results and RNA-seq FPKM values for each potato gene, as indicated. The per gene and average correlations between RNA-seq and qPCR are listed at the bottom of the page. (XLS 34 KB)
MA plots (hybridization intensity plotted against the fold change in expression) of between line comparisons.
Additional file 3: Each x-axis indicates mean gene expression levels [log2(FPKM1)+log2(FPKM2)] across the two time points selected for comparison. Each y-axis indicates fold change values [log2(FPKM1/FPKM2)]. (A) WT compared to +RB at 0 hpi; (B): WT compared to +RB at 24 hpi; (C): WT compared to +RB at 48 hpi. (PDF 485 KB)
Hierarchical clustering and Treeview visualization of 1,767 DE genes (determined based on between time point comparisons; Note: genes that are DE in between time point comparisons in water-inoculated samples were exclude from this analysis) show that DE gene regulation patterns in compatible and incompatible interactions are predominantly similar.
Additional file 4: Each column represents a comparison between two time points. Column one: 0 hpi to 24 hpi in WT; column two: 0 hpi to 24 hpi in +RB; column three: 24 hpi to 48 hpi in WT; column four: 24 hpi to 48 hpi in +RB. Red indicates up-regulation, green indicates down regulation. Left panel: the overall pattern of the 1,767 genes. Right panel: Magnified images of small gene clusters. (PDF 675 KB)
Additional file 5: +RB but not WT. Column one and two are gene ID and PGSC annotation descriptions. Columns three to six are log2 fold change values derived from between treatment comparisons (P. infestans- vs. water-inoculated). Values highlighted in red are statistically different between the +RB and WT lines. (PDF 78 KB)
Additional file 6: +RB ) were inoculated with P. infestans and water. We compared the RNA-seq FPKM counts for WT and +RB using all 39,031 gene models included in the Potato Genome Sequencing Consortium (PGSC) v3 dataset (i.e., all genes were included, regardless of whether or not a given gene was DE). Genes were grouped into ontology bins using a Mapman mapping file. Each column represents a comparison between the two genotypes at a defined time point post inoculation, as indicated. Bins in blue are transcribed at higher levels in WT than in +RB; bins in red are transcribed at higher levels in +RB than in WT; bins in white did not significantly differ in transcript levels between WT and +RB. Results indicate that stronger activation of defense bins, including stress responses and receptor kinases, occurred in +RB (the tuber blight resistant line). (PDF 2 MB)
Additional file 7: RB samples at 0, 24 and 48 hours post P. infestans inoculation were collected from three bio-reps and pooled into a single composite sample for RNA-seq. A total of 146.6 million single end Illumina reads (51 bp) were filtered and mapped to the reference potato genome using SolexaQA and Tuxedo software suite packages . Cuffdiff was used to generate log2 transformed fold change values for each between genotype comparisons at 0, 24 and 48 hpi. Mapman analyses and Wilcoxon rank sum tests were performed (see methods). Columns one, two, and three represent +RB vs. WT comparisons at 0, 24 and 48 hpi, respectively. Blue bins show higher transcription in WT; red bins show higher transcription in +RB. Note that +RB has faster and stronger activation of defense related bins (stress and receptor kinases) at 24 and 48 hpi. (PDF 1 MB)
Additional file 8: Phytophthora infestans transcripts. The X-axis indicates different time points (0, 24, 48 hpi) post P. infestans inoculation. The Y-axis indicates log2 transformed total mapped reads count. Results indicate an increase in P. infestans RNA-seq reads in the WT but not the +RB line over time. (PNG 33 KB)
Additional file 9: Phytophthora infestans genes. The first column contains gene ID, the second column indicates whether the gene is a predicted RxLR effector gene, each remaining column contains mapping summary results for each of the 36 potato tuber samples. (XLS 1 MB)
Additional file 10: ipiO genes. The first and second columns are ref id and short IDs for the ipiO genes adopted from Champouret et al. . Each remaining column contains a bowtie mapping summary for each sample following tuber inoculation with Phytophthora infestans.(XLS 34 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Gao, L., Tu, Z.J., Millett, B.P. et al. Insights into organ-specific pathogen defense responses in plants: RNA-seq analysis of potato tuber-Phytophthora infestans interactions. BMC Genomics 14, 340 (2013). https://doi.org/10.1186/1471-2164-14-340
- Solanum tuberosum
- Next-generation sequencing (NGS)
- Phytophthora infestans
- FPKM (fragments per kilobase of exon per million mapped reads)
- JA (jasmonic acid)
- SA (salicylic acid)
- ET (Ethylene)