Pseudoperonospora cubensis, an oomycete, is the causal agent of cucurbit downy mildew, and is responsible for significant losses on cucurbit crops worldwide. While other oomycete plant pathogens have been extensively studied at the molecular level, Ps. cubensis and the molecular basis of its interaction with cucurbit hosts has not been well examined. Here, we present the first large-scale global gene expression analysis of Ps. cubensis infection of a susceptible Cucumis sativus cultivar, ‘Vlaspik’, and identification of genes with putative roles in infection, growth, and pathogenicity. Using high throughput whole transcriptome sequencing, we captured differential expression of 2383 Ps. cubensis genes in sporangia and at 1, 2, 3, 4, 6, and 8 days post-inoculation (dpi). Additionally, comparison of Ps. cubensis expression profiles with expression profiles from an infection time course of the oomycete pathogen Phytophthora infestans on Solanum tuberosum revealed similarities in expression patterns of 1,576–6,806 orthologous genes suggesting a substantial degree of overlap in molecular events in virulence between the biotrophic Ps. cubensis and the hemi-biotrophic P. infestans. Co-expression analyses identified distinct modules of Ps. cubensis genes that were representative of early, intermediate, and late infection stages. Collectively, these expression data have advanced our understanding of key molecular and genetic events in the virulence of Ps. cubensis and thus, provides a foundation for identifying mechanism(s) by which to engineer or effect resistance in the host.
References
[1]
Savory EA, Granke LL, Quesada-Ocampo LM, Varbanova M, Hausbeck MK, et al. (2011) The cucurbit downy mildew pathogen Pseudoperonospora cubensis. Mol Plant Pathol 12: 217–226.
[2]
Lebeda A, Cohen Y (2011) Cucurbit downy mildew (Pseudoperonospora cubensis)—biology, ecology, epidemiology, host-pathogen interaction and control. Eur J Plant Pathol 129: 157–192.
[3]
Thomas C, Inaba T, Cohen Y (1987) Physiological specialization in Pseudoperonospora cubensis. Phytopathol 77: 1621–1624.
[4]
Sarris P, Abdelhalim M, Kitner M, Skandalis N, Panopoulos N, et al. (2009) Molecular polymorphisms between populations of Pseudoperonospora cubensis from Greece and the Czech Republic and the phytopathological and phylogenetic implications. Plant Pathol 58: 933–943.
[5]
Lebeda A, Widrlechner MP (2003) A set of Cucurbitaceae taxa for differentiation of Pseudoperonespora cubensis pathotypes. J Plant Dis Prot 110: 337–349.
[6]
Runge F, Choi Y-J, Thines M (2011) Phylogenetic investigations in the genus Pseudoperonospora reveal overlooked species and cryptic diversity in the P. cubensis species cluster. Eur J Plant Pathol 129: 135–146.
[7]
Choi Y (2005) A re-consideration of Pseudoperonospora cubensis and P. humuli based on molecular and morphological data. Mycological Res 109: 841–848.
[8]
Savory EA, Zou C, Adhikari BN, Hamilton JP, Buell CR, et al. (2012) Alternative splicing of a multi-drug transporter from Pseudoperonospora cubensis generates an RXLR effector protein that elicits a rapid cell death. PLoS ONE. 10.1371/journal.pone.0034701.
[9]
Tian M, Win J, Savory E, Burkhardt A, Held M, et al. (2011) 454 Genome sequencing of Pseudoperonospora cubensis reveals effector proteins with a QXLR translocation motif. Mol Plant Microbe Interact 24: 543–553.
[10]
Haas BJ, Kamoun S, Zody MC, Jiang RH, Handsaker RE, et al. (2009) Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature 461: 393–398.
[11]
Kamoun S (2006) A catalogue of the effector secretome of plant pathogenic oomycetes. Ann Rev Phytopathol 44: 41–60.
[12]
Cabral A, Stassen JH, Seidl MF, Bautor J, Parker JE, et al. (2011) Identification of Hyaloperonospora arabidopsidis transcript sequences expressed during Infection reveals isolate-specific effectors. PLoS ONE 6: e19328.
[13]
Sierra R, Rodriguez RL, Chaves D, Pinzon A, Grajales A, et al. (2010) Discovery of Phytophthora infestans genes expressed in planta through mining of cDNA libraries. PLoS ONE 5: e9847.
[14]
Torto-Alalibo TA, Tripathy S, Smith BM, Arredondo FD, Zhou L, et al. (2007) Expressed sequence tags from Phytophthora sojae reveal genes specific to development and infection. Mol Plant Microbe Interact 20: 781–793.
[15]
Randall TA, Dwyer RA, Huitema E, Beyer K, Cvitanich C, et al. (2005) Large-scale gene discovery in the oomycete Phytophthora infestans reveals likely components of phytopathogenicity shared with true fungi. Mol Plant Microbe Interact 18: 229–243.
[16]
Bos JI, Kanneganti TD, Young C, Cakir C, Huitema E, et al. (2006) The C-terminal half of Phytophthora infestans RXLR effector AVR3a is sufficient to trigger R3a-mediated hypersensitivity and suppress INF1-induced cell death in Nicotiana benthamiana. Plant J 48: 165–176.
[17]
Fabro G, Steinbrenner J, Coates M, Ishaque N, Baxter L, et al. (2011) Multiple Candidate Effectors from the Oomycete Pathogen Hyaloperonospora arabidopsidis Suppress Host Plant Immunity. PLoS Pathog 7: e1002348.
[18]
Armstrong MR, Whisson SC, Pritchard L, Bos JI, Venter E, et al. (2005) An ancestral oomycete locus contains late blight avirulence gene Avr3a, encoding a protein that is recognized in the host cytoplasm. Proc Natl Acad Sci U S A 102: 7766–7771.
[19]
Allen RL (2004) Host-Parasite Coevolutionary Conflict Between Arabidopsis and Downy Mildew. Science 306: 1957–1960.
[20]
Money NP, Davis CM, Ravishankar JP (2004) Biomechanical evidence for convergent evolution of the invasive growth process among fungi and oomycete water molds. Fung Gen Biol 41: 872–876.
[21]
Dong S, Qutob D, Tedman-Jones J, Kuflu K, Wang Y, et al. (2009) The Phytophthora sojae avirulence locus Avr3c encodes a multi-copy RXLR effector with sequence polymorphisms among pathogen strains. PLoS ONE 4: e5556.
[22]
Torto TA (2003) EST mining and functional expression fssays identify extracellular effector proteins from the plant pathogen Phytophthora. Genome Res 13: 1675–1685.
[23]
Schornack S, van Damme M, Bozkurt TO, Cano LM, Smoker M, et al. (2010) Ancient class of translocated oomycete effectors targets the host nucleus. Proc Natl Acad Sci U S A 107: 17421–17426.
[24]
Damasceno CM, Bishop JG, Ripoll DR, Win J, Kamoun S, et al. (2008) Structure of the glucanase inhibitor protein (GIP) family from phytophthora species suggests coevolution with plant endo-beta-1,3-glucanases. Mol Plant Microbe Interact 21: 820–830.
[25]
Rose JK, Ham KS, Darvill AG, Albersheim P (2002) Molecular cloning and characterization of glucanase inhibitor proteins: coevolution of a counterdefense mechanism by plant pathogens. Plant Cell 14: 1329–1345.
[26]
Tian M, Benedetti B, Kamoun S (2005) A Second Kazal-like protease inhibitor from Phytophthora infestans inhibits and interacts with the apoplastic pathogenesis-related protease P69B of tomato. Plant Physiol 138: 1785–1793.
[27]
Tian M, Huitema E, Da Cunha L, Torto-Alalibo T, Kamoun S (2004) A Kazal-like extracellular serine protease inhibitor from Phytophthora infestans targets the tomato pathogenesis-related protease P69B. J Biol Chem 279: 26370–26377.
[28]
Tian M, Win J, Song J, van der Hoorn R, van der Knaap E, et al. (2007) A Phytophthora infestans cystatin-like protein targets a novel tomato papain-like apoplastic protease. Plant Physiol 143: 364–377.
[29]
Liu Z, Bos JI, Armstrong M, Whisson SC, da Cunha L, et al. (2005) Patterns of diversifying selection in the phytotoxin-like scr74 gene family of Phytophthora infestans. Mol Biol Evol 22: 659–672.
[30]
Levesque CA, Brouwer H, Cano L, Hamilton JP, Holt C, et al. (2010) Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire. Genome Biol 11: R73.
[31]
Fellbrich G, Romanski A, Varet A, Blume B, Brunner F, et al. (2002) NPP1, a Phytophthora-associated trigger of plant defense in parsley and Arabidopsis. Plant J 32: 375–390.
[32]
Cabral A, Oome S, Sander N, Kuefner I, Nürnberger T, et al. (2012) Non-toxic Nep1-like proteins of the downy mildew pathogen Hyaloperonospora arabidopsidis; repression of necrosis-inducing activity by a surface-exposed region. Mol Plant Microbe Interact.
[33]
Gaulin E, Drame N, Lafitte C, Torto-Alalibo T, Martinez Y, et al. (2006) Cellulose binding domains of a Phytophthora cell wall protein are novel pathogen-associated molecular patterns. Plant Cell 18: 1766–1777.
[34]
Gaulin E, Jauneau A, Villalba F, Rickauer M, Esquerré-Tugayé MT, et al. (2002) The CBEL glycoprotein of Phytophthora parasitica var-nicotianae is involved in cell wall deposition and adhesion to cellulosic substrates. J Cell Sci 115: 4565.
[35]
Tyler BM (2006) Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis. Science 313: 1261–1266.
[36]
Baxter L, Tripathy S, Ishaque N, Boot N, Cabral A, et al. (2010) Signatures of adaptation to obligate biotrophy in the Hyaloperonospora arabidopsidis genome. Science 330: 1549–1551.
[37]
Moy P, Qutob D, Chapman BP, Atkinson I, Gijzen M (2004) Patterns of gene expression upon infection of soybean plants by Phytophthora sojae. Mol Plant Microbe Interact 17: 1051–1062.
[38]
Costanzo S, Ospina-Giraldo MD, Deahl KL, Baker CJ, Jones RW (2006) Gene duplication event in family 12 glycosyl hydrolase from Phytophthora spp. Fung Gen Biol 43: 707–714.
[39]
Yamaguchi Y, Huffaker A (2011) Endogenous peptide elicitors in higher plants. Current opinion in plant biology 14: 351–357.
[40]
Fernandez D, Tisserant E, Talhinhas P, Azinheira H, Vieira ANA, et al. (2011) 454-pyrosequencing of Coffea arabica leaves infected by the rust fungus Hemileia vastatrix reveals in planta-expressed pathogen-secreted proteins and plant functions in a late compatible plant–rust interaction. Mol Plant Pathol 13: 17–37.
[41]
Joly DL, Feau N, Tanguay P, Hamelin RC (2010) Comparative analysis of secreted protein evolution using expressed sequence tags from four poplar leaf rusts (Melampsora spp.). BMC Gen 11: 422.
[42]
Miranda M, Ralph SG, Mellway R, White R, Heath MC, et al. (2007) The transcriptional response of hybrid poplar (Populus trichocarpa×P. deltoides) to infection by Melampsora medusae leaf rust involves induction of flavonoid pathway genes leading to the accumulation of proanthocyanidins. Mol Plant Microbe Interact 20: 816–931.
[43]
Duplessis S, Hacquard S, Delaruelle C, Tisserant E, Frey P, et al. (2011) Melampsora larici-populina tanscript profiling during germination and timecourse infection of poplar leaves reveals dynamic expression patterns associated with virulence and biotrophy. Mol Plant Microbe Interact 24: 808–818.
[44]
Mosquera G, Giraldo M, Khang CH, Coughlan S, Valent B (2009) Interaction transcriptome analysis identifies Magnaporthe oryzae BAS1-4 as biotrophy-associated secreted proteins in rice blast disease. Plant Cell 21: 1273–1290.
[45]
Polesani M, Desario F, Ferrarini A, Zamboni A, Pezzotti M, et al. (2008) cDNA-AFLP analysis of plant and pathogen genes expressed in grapevine infected with Plasmopara viticola. BMC Gen 9: 142–156.
[46]
Adhikari B, Savory E, Vaillancourt B, Childs KL, Hamilton JP, et al. (2012) Expression profiling of Cucumis sativus in response to infection by Pseudoperonospora cubensis. PLOS ONE. in press.
[47]
Ospina-Giraldo M, Griffith J, Laird E, Mingora C (2010) The CAZyome of Phytophthora spp.: A comprehensive analysis of the gene complement coding for carbohydrate-active enzymes in species of the genus Phytophthora. BMC Gen 11: 1–16.
[48]
Ihmels J, Bergmann S, Berman J, Barkai N (2005) Comparative gene expression analysis by a differential clustering approach: Application to the Candida albicans transcription program. PLoS Genet 1: e39.
[49]
Langfelder P, Zhang B, Horvath S (2008) Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics 24: 719–720.
[50]
Langfelder P, Horvath S (2007) Eigengene networks for studying the relationships between co-expression modules. BMC Sys Biol 1: 54.
[51]
Iyer LM, Anantharaman V, Wolf MY, Aravind L (2008) Comparative genomics of transcription factors and chromatin proteins in parasitic protists and other eukaryotes. In J Parasitol 38: 1–31.
[52]
Wang Y, Dou D, Wang X, Li A, Sheng Y, et al. (2009) The PsCZF1 gene encoding a C2H2 zinc finger protein is required for growth, development and pathogenesis in Phytophthora sojae. Micro Path 47: 78–86.
[53]
Blanco FA, Judelson HS (2005) A bZIP transcription factor from Phytophthora interacts with a protein kinase and is required for zoospore motility and plant infection. Mol Micro 56: 638–648.
[54]
Judelson H, Ah-Fong A (2010) The kinome of Phytophthora infestans reveals oomycete-specific innovations and links to other taxonomic groups. BMC Gen 11: 700.
[55]
Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.
[56]
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28: 511–515.
Bullard J, Purdom E, Hansen K, Dudoit S (2010) Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics 11: 94.
[59]
Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30: 207–210.
[60]
Irizarry R, Hobbs B, Collin F, Beazer-Barclay Y, Antonellis K, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4: 249–264.
[61]
Chen F, Mackey AJ, Vermunt JK, Roos DS (2007) Assessing performance of orthology detection strategies applied to eukaryotic genomes. PLoS ONE 2: e383.
Bateman A, Birney E, Durbin R, Eddy SR, Howe KL, et al. (2000) The Pfam protein families database. Nucleic Acids Res 28: 263–266.
[64]
Eddy SR (2009) A new generation of homology search tools based on probabilistic inference. Genome Inform 23: 205–211.
[65]
Childs KL, Davidson RM, Buell CR (2011) Gene coexpression network analysis as a source of functional annotation for rice genes. PLoS ONE 6: e22196.
[66]
Park BH, Karpinets TV, Syed MH, Leuze MR, Uberbacher EC (2010) CAZymes Analysis Toolkit (CAT): Web-service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. Glycobiology.
[67]
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, et al. (2009) The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res 37: D233–D238.
