De novo Analysis of the Epiphytic Transcriptome of the Cucurbit Powdery Mildew Fungus Podosphaera xanthii and Identification of Candidate Secreted Effector Proteins

The cucurbit powdery mildew fungus Podosphaera xanthii is a major limiting factor for cucurbit production worldwide. Despite the fungus’s agronomic and economic importance, very little is known about fundamental aspects of P. xanthii biology, such as obligate biotrophy or pathogenesis. To design more durable control strategies, genomic information about P. xanthii is needed. Powdery mildews are fungal pathogens with large genomes compared with those of other fungi, which contain vast amounts of repetitive DNA sequences, much of which is composed of retrotransposons. To reduce genome complexity, in this work we aimed to obtain and analyse the epiphytic transcriptome of P. xanthii as a starting point for genomic research. Total RNA was isolated from epiphytic fungal material, and the corresponding cDNA library was sequenced using a 454 GS FLX platform. Over 676,562 reads were obtained and assembled into 37,241 contigs. Annotation data identified 8,798 putative genes with different orthologues. As described for other powdery mildew fungi, a similar set of missing core ascomycete genes was found, which may explain obligate biotrophy. To gain insight into the plant-pathogen relationships, special attention was focused on the analysis of the secretome. After this analysis, 137 putative secreted proteins were identified, including 53 candidate secreted effector proteins (CSEPs). Consistent with a putative role in pathogenesis, the expression profile observed for some of these CSEPs showed expression maxima at the beginning of the infection process at 24 h after inoculation, when the primary appressoria are mostly formed. Our data mark the onset of genomics research into this very important pathogen of cucurbits and shed some light on the intimate relationship between this pathogen and its host plant.

[1]  P. Spanu The genomics of obligate (and nonobligate) biotrophs. , 2012, Annual review of phytopathology.

[2]  P. Auvinen,et al.  SNP Design from 454 Sequencing of Podosphaera plantaginis Transcriptome Reveals a Genetically Diverse Pathogen Metapopulation with High Levels of Mixed-Genotype Infection , 2012, PloS one.

[3]  Patrick Schweizer,et al.  HIGS: Host-Induced Gene Silencing in the Obligate Biotrophic Fungal Pathogen Blumeria graminis[W][OA] , 2010, Plant Cell.

[4]  S. Reissmann,et al.  Fungal effectors and plant susceptibility. , 2015, Annual review of plant biology.

[5]  James K. M. Brown,et al.  Coevolution between a Family of Parasite Virulence Effectors and a Class of LINE-1 Retrotransposons , 2009, PloS one.

[6]  G. Dimopoulos,et al.  Transcript profiles of Blumeria graminis development during infection reveal a cluster of genes that are potential virulence determinants. , 2005, Molecular plant-microbe interactions : MPMI.

[7]  H. Giese,et al.  A Blumeria graminis gene family encoding proteins with a C-terminal variable region with homologues in pathogenic fungi. , 2003, Gene.

[8]  N. Saitou,et al.  The neighbor-joining method: a new method for reconstructing phylogenetic trees. , 1987, Molecular biology and evolution.

[9]  P. Dodds,et al.  Recent progress in discovery and functional analysis of effector proteins of fungal and oomycete plant pathogens. , 2009, Current opinion in plant biology.

[10]  R. Schneiter,et al.  The CAP protein superfamily: function in sterol export and fungal virulence , 2013, Biomolecular concepts.

[11]  R. Dean,et al.  The adenylate cyclase gene MAC1 of Magnaporthe grisea controls appressorium formation and other aspects of growth and development. , 1997, The Plant cell.

[12]  Cathryn J. Rehmeyer,et al.  The genome sequence of the rice blast fungus Magnaporthe grisea , 2005, Nature.

[13]  A. Di Pietro,et al.  Fusarium oxysporum gas1 encodes a putative beta-1,3-glucanosyltransferase required for virulence on tomato plants. , 2005, Molecular plant-microbe interactions : MPMI.

[14]  Mukesh Jain,et al.  A next-generation approach to the characterization of a non-model plant transcriptome * , 2011 .

[15]  E. Selker,et al.  DNA methylation associated with repeat-induced point mutation in Neurospora crassa , 1995, Molecular and cellular biology.

[16]  X. Huang,et al.  CAP3: A DNA sequence assembly program. , 1999, Genome research.

[17]  A. Pérez-García,et al.  Occurrence of races and pathotypes of cucurbit powdery mildew in southeastern Spain , 2002, Phytoparasitica.

[18]  G. Camici,et al.  Purification, Characterization, and Amino Acid Sequence of Cerato-platanin, a New Phytotoxic Protein from Ceratocystis fimbriata f. sp. platani * , 1999, The Journal of Biological Chemistry.

[19]  A. Pérez-García,et al.  Long‐term Preservation of Podosphaera fusca Using Silica Gel , 2006 .

[20]  Katherine C. H. Amrine,et al.  Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator , 2014, BMC Genomics.

[21]  M. Pfaffl,et al.  A new mathematical model for relative quantification in real-time RT-PCR. , 2001, Nucleic acids research.

[22]  R. Dean,et al.  An eight-cysteine-containing CFEM domain unique to a group of fungal membrane proteins. , 2003, Trends in biochemical sciences.

[23]  S. Stock,et al.  Breast Cancer Attributable Costs in Germany: A Top-Down Approach Based on Sickness Funds Data , 2012, PloS one.

[24]  M. Stumpf,et al.  Gene Expression Profiles of Blumeria graminis Indicate Dynamic Changes to Primary Metabolism during Development of an Obligate Biotrophic Pathogenw⃞ , 2005, The Plant Cell Online.

[25]  Marta Matvienko,et al.  De novo assembly and characterization of the carrot transcriptome reveals novel genes, new markers, and genetic diversity , 2011, BMC Genomics.

[26]  Tom Hsiang,et al.  Specific adaptation of Ustilaginoidea virens in occupying host florets revealed by comparative and functional genomics , 2014, Nature Communications.

[27]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[28]  N. Talbot,et al.  Comparative genomic analysis of phytopathogenic fungi using expressed sequence tag (EST) collections. , 2006, Molecular plant pathology.

[29]  D. Glawe,et al.  The powdery mildews: a review of the world's most familiar (yet poorly known) plant pathogens. , 2008, Annual review of phytopathology.

[30]  F. López-Ruiz,et al.  The Podosphaera fusca TUB2 gene, a molecular "Swiss Army knife" with multiple applications in powdery mildew research. , 2014, Fungal biology.

[31]  R. Panstruga,et al.  Powdery mildew genomes in the crosshairs. 2nd International Powdery Mildew Workshop and 3rd New Phytologist Workshop, in Zürich, Switzerland, February 2012. , 2012, The New phytologist.

[32]  K. Mendgen,et al.  High level activation of vitamin B1 biosynthesis genes in haustoria of the rust fungus Uromyces fabae. , 2000, Molecular plant-microbe interactions : MPMI.

[33]  J. Sweigard,et al.  Magnaporthe grisea Pth11p Is a Novel Plasma Membrane Protein That Mediates Appressorium Differentiation in Response to Inductive Substrate Cues , 1999, Plant Cell.

[34]  Pari Skamnioti,et al.  Genome Expansion and Gene Loss in Powdery Mildew Fungi Reveal Tradeoffs in Extreme Parasitism , 2010, Science.

[35]  B. Barrell,et al.  Insight into the genome of Aspergillus fumigatus: analysis of a 922 kb region encompassing the nitrate assimilation gene cluster. , 2004, Fungal genetics and biology : FG & B.

[36]  B. Álvarez,et al.  Cultivo in vitro de Sphaerotheca fuliginea (Schlecht. ex Fr.), efecto de diferentes fuentes de carbono sobre su desarrollo , 1997 .

[37]  Liam J McGuffin,et al.  Structure and evolution of barley powdery mildew effector candidates , 2012, BMC Genomics.

[38]  J. Felsenstein CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP , 1985, Evolution; international journal of organic evolution.

[39]  Gertraud Burger,et al.  AutoFACT: An Automatic Functional Annotation and Classification Tool , 2005, BMC Bioinformatics.

[40]  Jeppe Emmersen,et al.  Powdery mildew fungal effector candidates share N-terminal Y/F/WxC-motif , 2010, BMC Genomics.

[41]  Alejandro Pérez-García,et al.  The powdery mildew fungus Podosphaera fusca (synonym Podosphaera xanthii), a constant threat to cucurbits. , 2009, Molecular plant pathology.

[42]  F. López-Ruiz,et al.  Occurrence and distribution of resistance to QoI fungicides in populations of Podosphaera fusca in south central Spain , 2006, European Journal of Plant Pathology.

[43]  F. Sbrana,et al.  Atomic force microscopy images suggest aggregation mechanism in cerato-platanin , 2007, European Biophysics Journal.

[44]  Leighton Pritchard,et al.  A translocation signal for delivery of oomycete effector proteins into host plant cells , 2007, Nature.

[45]  J. Rodriguez,et al.  De Novo Assembly and Functional Annotation of the Olive (Olea europaea) Transcriptome , 2013, DNA research : an international journal for rapid publication of reports on genes and genomes.

[46]  A. Ramos,et al.  The structure of the C-terminal KH domains of KSRP reveals a noncanonical motif important for mRNA degradation. , 2007, Structure.

[47]  S. Christiansen,et al.  Gene identification in the obligate fungal pathogen Blumeria graminis by expressed sequence tag analysis. , 2001, Fungal genetics and biology : FG & B.

[48]  H. Quesneville,et al.  The wheat powdery mildew genome shows the unique evolution of an obligate biotroph , 2013, Nature Genetics.

[49]  T. Cooper,et al.  Allantoin transport in Saccharomyces cerevisiae is regulated by two induction systems , 1987, Journal of bacteriology.

[50]  P. Pevzner,et al.  An Eulerian path approach to DNA fragment assembly , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[51]  Juan Miguel García-Gómez,et al.  BIOINFORMATICS APPLICATIONS NOTE Sequence analysis Manipulation of FASTQ data with Galaxy , 2005 .

[52]  R. Reinhardt,et al.  Transcriptome analysis of enriched Golovinomyces orontii haustoria by deep 454 pyrosequencing. , 2012, Fungal genetics and biology : FG & B.

[53]  M. Bölker,et al.  A H2O2-producing glyoxal oxidase is required for filamentous growth and pathogenicity in Ustilago maydis , 2005, Molecular Genetics and Genomics.

[54]  D. Roby,et al.  Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum , 2014, BMC Genomics.

[55]  James K. Hane,et al.  Repeat-Induced Point Mutation: A Fungal-Specific, Endogenous Mutagenesis Process , 2015 .

[56]  K. Hammond-Kosack,et al.  The Predicted Secretome of the Plant Pathogenic Fungus Fusarium graminearum: A Refined Comparative Analysis , 2012, PloS one.

[57]  Thomas Wetter,et al.  Genome Sequence Assembly Using Trace Signals and Additional Sequence Information , 1999, German Conference on Bioinformatics.

[58]  M. O’Bryan,et al.  The CAP superfamily: cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins--roles in reproduction, cancer, and immune defense. , 2008, Endocrine reviews.

[59]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[60]  Jaideep P. Sundaram,et al.  Genomic Islands in the Pathogenic Filamentous Fungus Aspergillus fumigatus , 2008, PLoS genetics.