Insights into organ-specific pathogen defense responses in plants: RNA-seq analysis of potato tuber-Phytophthora infestans interactions
© Gao et al.; licensee BioMed Central Ltd. 2013
Received: 17 December 2012
Accepted: 17 May 2013
Published: 23 May 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.
KeywordsSolanum tuberosum Next-generation sequencing (NGS) RNA-seq Phytophthora infestans FPKM (fragments per kilobase of exon per million mapped reads) JA (jasmonic acid) SA (salicylic acid) ET (Ethylene) qRT-PCR
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
RNA-seq reads aligned well with the potato reference genome sequence
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
Differentially expressed (DE) genes
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
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.
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.
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.
- Stevenson WR, Loria R, Frank GD, Weingartner DP: Compedium of potato diseases. 2001, Saint Paul, MN, USA: APS PressGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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-Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- Metzker ML: Applications of next generation sequencing sequencing technologies - the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.19220.127.116.116.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Saldanha AJ: Java treeview-extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- 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-0Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.