Genome sequencing and transcript analysis of Hemileia vastatrix reveal expression dynamics of candidate effectors dependent on host compatibility

Coffee leaf rust caused by the fungus Hemileia vastatrix is one of the most important leaf diseases of coffee plantations worldwide. Current knowledge of the H. vastatrix genome is limited and only a small fraction of the total fungal secretome has been identified. In order to obtain a more comprehensive understanding of its secretome, we aimed to sequence and assemble the entire H. vastatrix genome using two next-generation sequencing platforms and a hybrid assembly strategy. This resulted in a 547 Mb genome of H. vastatrix race XXXIII (Hv33), with 13,364 predicted genes that encode 13,034 putative proteins with transcriptomic support. Based on this proteome, 615 proteins contain putative secretion peptides, and lack transmembrane domains. From this putative secretome, 111 proteins were identified as candidate effectors (EHv33) unique to H. vastatrix, and a subset consisting of 17 EHv33 genes was selected for a temporal gene expression analysis during infection. Five genes were significantly induced early during an incompatible interaction, indicating their potential role as pre-haustorial effectors possibly recognized by the resistant coffee genotype. Another nine genes were significantly induced after haustorium formation in the compatible interaction. Overall, we suggest that this fungus is able to selectively mount its survival strategy with effectors that depend on the host genotype involved in the infection process.

[1]  Keith Bradnam,et al.  CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes , 2007, Bioinform..

[2]  Julie Poulain,et al.  Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection , 2014, Front. Plant Sci..

[3]  Burkhard Morgenstern,et al.  AUGUSTUS: a web server for gene finding in eukaryotes , 2004, Nucleic Acids Res..

[4]  Ana Vieira,et al.  The coffee leaf rust pathogen Hemileia vastatrix: one and a half centuries around the tropics. , 2017, Molecular plant pathology.

[5]  A. Krogh,et al.  Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. , 2001, Journal of molecular biology.

[6]  J. Cornelis,et al.  Impact of rice cultivar and organ on elemental composition of phytoliths and the release of bio-available silicon , 2014, Front. Plant Sci..

[7]  Paul Horton,et al.  Nucleic Acids Research Advance Access published May 21, 2007 WoLF PSORT: protein localization predictor , 2007 .

[8]  R. Gibbs,et al.  Mind the Gap: Upgrading Genomes with Pacific Biosciences RS Long-Read Sequencing Technology , 2012, PloS one.

[9]  David A. Jones,et al.  The genome sequence and effector complement of the flax rust pathogen Melampsora lini , 2014, Front. Plant Sci..

[10]  S. Brunak,et al.  Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. , 2000, Journal of molecular biology.

[11]  E. Gichuru,et al.  Additional physiological races of coffee leaf rust (Hemileia vastatrix) identified in Kenya , 2012 .

[12]  Marco Cristancho,et al.  The coffee rust crises in Colombia and Central America (2008–2013): impacts, plausible causes and proposed solutions , 2015, Food Security.

[13]  O. Paulo,et al.  Validation of RT-qPCR reference genes for in planta expression studies in Hemileia vastatrix, the causal agent of coffee leaf rust. , 2011, Fungal biology.

[14]  M. C. Heath.  A comparative study of non-host interactions with rust fungi , 1977 .

[15]  A. Petitot,et al.  Cellular and molecular analyses of coffee resistance to Hemileia vastatrix and nonhost resistance to Uromyces vignae in the resistance-donor genotype HDT832/2 , 2011, European Journal of Plant Pathology.

[16]  B. Bertrand,et al.  Coffee resistance to the main diseases: leaf rust and coffee berry disease , 2006 .

[17]  S. Robatzek,et al.  Breaking the barriers: microbial effector molecules subvert plant immunity. , 2008, Annual review of phytopathology.

[18]  Ana Vieira,et al.  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. , 2012, Molecular plant pathology.

[19]  A. Lecouls,et al.  Monitoring of the early molecular resistance responses of coffee (Coffea arabica L.) to the rust fungus (Hemileia vastatrix) using real-time quantitative RT-PCR , 2006 .

[20]  M. Vainstein,et al.  Fungal zinc metabolism and its connections to virulence , 2013, Front. Cell. Infect. Microbiol..

[21]  R. Oliver,et al.  Effectors as tools in disease resistance breeding against biotrophic, hemibiotrophic, and necrotrophic plant pathogens. , 2014, Molecular plant-microbe interactions : MPMI.

[22]  R. Terauchi,et al.  Emerging concepts in effector biology of plant-associated organisms. , 2009, Molecular plant-microbe interactions : MPMI.

[23]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[24]  M. C. Heath.  Signalling between Pathogenic Rust Fungi and Resistant or Susceptible Host Plants , 1997 .

[25]  J. Jurka,et al.  Next Generation Sequencing Provides Rapid Access to the Genome of Puccinia striiformis f. sp. tritici, the Causal Agent of Wheat Stripe Rust , 2011, PloS one.

[26]  B. Bertrand,et al.  Next generation variety development for sustainable production of arabica coffee (Coffea arabica L.): a review , 2015, Euphytica.

[27]  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.

[28]  S. McCook,et al.  The Big Rust and the Red Queen: Long-Term Perspectives on Coffee Rust Research. , 2015, Phytopathology.

[29]  N. Friedman,et al.  Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data , 2011, Nature Biotechnology.

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

[31]  Ian Korf,et al.  MaskerAid : a performance enhancement to RepeatMasker , 2000, Bioinform..

[32]  C. J. Rodrigues,et al.  Hypersensitive cell death and post-haustorial defence responses arrest the orange rust (Hemileia vastatrix) growth in resistant coffee leaves , 2002 .

[33]  Fyodor A. Kondrashov,et al.  Sequence space and the ongoing expansion of the protein universe , 2010, Nature.

[34]  Yi Wang,et al.  OrthoVenn: a web server for genome wide comparison and annotation of orthologous clusters across multiple species , 2015, Nucleic Acids Res..

[35]  Jonathan D. G. Jones,et al.  The plant immune system , 2006, Nature.

[36]  Aaron A. Klammer,et al.  Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data , 2013, Nature Methods.

[37]  P. D. de Wit,et al.  Fungal effector proteins: past, present and future. , 2009, Molecular plant pathology.

[38]  Jian Wang,et al.  SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler , 2012, GigaScience.

[39]  M. Nicole,et al.  Involvement of peroxidases in the coffee resistance to orange rust (Hemileia vastatrix) , 2008 .

[40]  R. Voegele,et al.  Genome size analyses of Pucciniales reveal the largest fungal genomes , 2014, Front. Plant Sci..

[41]  R. Barreto,et al.  Coffee rust genome measured using flow cytometry: does size matter? , 2014 .

[42]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[43]  Christina A. Cuomo,et al.  Comparative genomics of Australian isolates of the wheat stem rust pathogen Puccinia graminis f. sp. tritici reveals extensive polymorphism in candidate effector genes , 2015, Front. Plant Sci..

[44]  S. Brommonschenkel,et al.  The Hemileia vastatrix effector HvEC‐016 suppresses bacterial blight symptoms in coffee genotypes with the SH1 rust resistance gene , 2016, The New phytologist.

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

[46]  T. D. de Kievit,et al.  Elucidating the Role of Effectors in Plant-Fungal Interactions: Progress and Challenges , 2016, Front. Microbiol..

[47]  Glenn Tesler,et al.  Mapping single molecule sequencing reads using basic local alignment with successive refinement (BLASR): application and theory , 2012, BMC Bioinformatics.

[48]  Christina A. Cuomo,et al.  Obligate Biotrophy Features Unraveled by the Genomic Analysis of the Rust Fungi, Melampsora larici-populina and Puccinia graminis f. sp. tritici , 2011 .

[49]  M. Nicole,et al.  Cytochemical Aspects of the Plant–Rust Fungus Interface during the Compatible Interaction Coffea arabica (cv. Caturra)–Hemileiavastatrix (race III) , 1999, International Journal of Plant Sciences.

[50]  M. Borodovsky,et al.  Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. , 2008, Genome research.

[51]  S. Brunak,et al.  SignalP 4.0: discriminating signal peptides from transmembrane regions , 2011, Nature Methods.

[52]  L. Zambolim,et al.  Current status and management of coffee leaf rust in Brazil , 2016, Tropical Plant Pathology.

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

[54]  O. Kohany,et al.  Repbase Update, a database of repetitive elements in eukaryotic genomes , 2015, Mobile DNA.

[55]  Sofia M. C. Robb,et al.  MAKER: an easy-to-use annotation pipeline designed for emerging model organism genomes. , 2007, Genome research.

[56]  P. D. de Wit,et al.  Fungal effector proteins. , 2009, Annual review of phytopathology.

[57]  E. Stukenbrock,et al.  The evolving fungal genome , 2014 .

[58]  Luis F. Rivera,et al.  Annotation of a hybrid partial genome of the coffee rust (Hemileia vastatrix) contributes to the gene repertoire catalog of the Pucciniales , 2014, Front. Plant Sci..

[59]  A. Petitot,et al.  Biphasic haustorial differentiation of coffee rust (Hemileia vastatrix race II) associated with defence responses in resistant and susceptible coffee cultivars , 2009 .

[60]  D. MacLean,et al.  Genome analyses of the wheat yellow (stripe) rust pathogen Puccinia striiformis f. sp. tritici reveal polymorphic and haustorial expressed secreted proteins as candidate effectors , 2013, BMC Genomics.