Genome-Wide Analysis of Corynespora cassiicola Leaf Fall Disease Putative Effectors

Corynespora cassiicola is an Ascomycetes fungus with a broad host range and diverse life styles. Mostly known as a necrotrophic plant pathogen, it has also been associated with rare cases of human infection. In the rubber tree, this fungus causes the Corynespora leaf fall (CLF) disease, which increasingly affects natural rubber production in Asia and Africa. It has also been found as an endophyte in South American rubber plantations where no CLF outbreak has yet occurred. The C. cassiicola species is genetically highly diverse, but no clear relationship has been evidenced between phylogenetic lineage and pathogenicity. Cassiicolin, a small glycosylated secreted protein effector, is thought to be involved in the necrotrophic interaction with the rubber tree but some virulent C. cassiicola isolates do not have a cassiicolin gene. This study set out to identify other putative effectors involved in CLF. The genome of a highly virulent C. cassiicola isolate from the rubber tree (CCP) was sequenced and assembled. In silico prediction revealed 2870 putative effectors, comprising CAZymes, lipases, peptidases, secreted proteins and enzymes associated with secondary metabolism. Comparison with the genomes of 44 other fungal species, focusing on effector content, revealed a striking proximity with phylogenetically unrelated species (Colletotrichum acutatum, Colletotrichum gloesporioides, Fusarium oxysporum, nectria hematococca, and Botrosphaeria dothidea) sharing life style plasticity and broad host range. Candidate effectors involved in the compatible interaction with the rubber tree were identified by transcriptomic analysis. Differentially expressed genes included 92 putative effectors, among which cassiicolin and two other secreted singleton proteins. Finally, the genomes of 35 C. cassiicola isolates representing the genetic diversity of the species were sequenced and assembled, and putative effectors identified. At the intraspecific level, effector-based classification was found to be highly consistent with the phylogenomic trees. Identification of lineage-specific effectors is a key step toward understanding C. cassiicola virulence and host specialization mechanisms.

[1]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[2]  Erik L. L. Sonnhammer,et al.  A Hidden Markov Model for Predicting Transmembrane Helices in Protein Sequences , 1998, ISMB.

[3]  P. Johnston,et al.  Colletotrichum – current status and future directions , 2012, Studies in mycology.

[4]  E. Birney,et al.  Velvet: algorithms for de novo short read assembly using de Bruijn graphs. , 2008, Genome research.

[5]  A. Clément-Vidal,et al.  Genetic Determinism of Sensitivity to Corynespora cassiicola Exudates in Rubber Tree (Hevea brasiliensis) , 2016, PloS one.

[6]  J. C. del Pozo,et al.  A bacterial cysteine protease effector protein interferes with photosynthesis to suppress plant innate immune responses , 2012, Cellular microbiology.

[7]  Xin Gao,et al.  Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups. , 2011, Current protocols in bioinformatics.

[8]  K. Nishimura,et al.  Rare case of fungal keratitis caused by Corynespora cassiicola , 2013, Journal of infection and chemotherapy : official journal of the Japan Society of Chemotherapy.

[9]  Keith Johnstone,et al.  Evidence for a role of cutinase in pathogenicity of Pyrenopeziza brassicae on brassicas , 2000 .

[10]  B. Lyon,et al.  Molecular, physiological and pathological characterization of Corynespora leaf spot fungi from rubber plantations in Sri Lanka , 1998 .

[11]  Sarah Calvo,et al.  Insights from the genome of the biotrophic fungal plant pathogen Ustilago maydis , 2006, Nature.

[12]  K. Hyde,et al.  Variation between freshwater and terrestrial fungal communities on decaying bamboo culms , 2006, Antonie van Leeuwenhoek.

[13]  Pari Skamnioti,et al.  Against the grain: safeguarding rice from rice blast disease. , 2009, Trends in biotechnology.

[14]  Pedro M. Coutinho,et al.  The carbohydrate-active enzymes database (CAZy) in 2013 , 2013, Nucleic Acids Res..

[15]  Ramón Doallo,et al.  ProtTest 3: fast selection of best-fit models of protein evolution , 2011, Bioinform..

[16]  Md. Sharifur Rahman,et al.  Tools to kill: Genome of one of the most destructive plant pathogenic fungi Macrophomina phaseolina , 2012, BMC Genomics.

[17]  Y. Benjamini,et al.  THE CONTROL OF THE FALSE DISCOVERY RATE IN MULTIPLE TESTING UNDER DEPENDENCY , 2001 .

[18]  N. Lemke,et al.  Cupin: A candidate molecular structure for the Nep1-like protein family , 2008, BMC Plant Biology.

[19]  P. Schulze-Lefert,et al.  Root Endophyte Colletotrichum tofieldiae Confers Plant Fitness Benefits that Are Phosphate Status Dependent , 2016, Cell.

[20]  W. Knogge,et al.  Amino acid alterations in isoforms of the effector protein NIP1 from Rhynchosporium secalis have similar effects on its avirulence- and virulence-associated activities on barley , 2002 .

[21]  V. Bouriotis,et al.  Purification and Characterization of Chitin Deacetylase from Colletotrichum lindemuthianum(*) , 1995, The Journal of Biological Chemistry.

[22]  A. Gnirke,et al.  High-quality draft assemblies of mammalian genomes from massively parallel sequence data , 2010, Proceedings of the National Academy of Sciences.

[23]  C. Sanier,et al.  Role of cassiicolin, a host-selective toxin, in pathogenicity of Corynespora cassiicola, causal agent of a leaf fall disease of Hevea. , 2000 .

[24]  R. Gazis,et al.  Diversity of fungal endophytes in leaves and stems of wild rubber trees (Hevea brasiliensis) in Peru , 2010 .

[25]  B. Henrissat,et al.  Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes , 2013, Biotechnology for Biofuels.

[26]  B. Dentinger,et al.  Assembling the Fungal Tree of Life: constructing the Structural and Biochemical Database , 2006, Mycologia.

[27]  Y. Zhang,et al.  A class-wide phylogenetic assessment of Dothideomycetes , 2009, Studies in mycology.

[28]  M. Kimura,et al.  A novel lactonohydrolase responsible for the detoxification of zearalenone: enzyme purification and gene cloning. , 2002, The Biochemical journal.

[29]  Kai Blin,et al.  antiSMASH 2.0—a versatile platform for genome mining of secondary metabolite producers , 2013, Nucleic Acids Res..

[30]  J. Woolliams,et al.  What is Genetic Diversity , 2007 .

[31]  Su Mei Yew,et al.  Genomic insight into pathogenicity of dematiaceous fungus Corynespora cassiicola , 2017, PeerJ.

[32]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[33]  G. Sherlock,et al.  Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads , 2010, BMC Genomics.

[34]  R. Kemerait,et al.  Clonality and geographic structure of host-specialized populations of Corynespora cassiicola causing emerging target spot epidemics in the southeastern United States , 2018, PloS one.

[35]  Y. Qi,et al.  Morphological and molecular analysis of genetic variability within isolates of Corynespora cassiicola from different hosts , 2011, European Journal of Plant Pathology.

[36]  P. Bowyer,et al.  Characterization of SNP1, a cell wall-degrading trypsin, produced during infection by Stagonospora nodorum. , 2000, Molecular plant-microbe interactions : MPMI.

[37]  E. Filho,et al.  An overview of mannan structure and mannan-degrading enzyme systems , 2008, Applied Microbiology and Biotechnology.

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

[39]  Averil E. Brown,et al.  Pathogenicity of Colletotrichum acutatum and C. gloeosporioides on leaves of Hevea spp. , 1994 .

[40]  J. Cairney,et al.  A simple and efficient method for isolating RNA from pine trees , 1993, Plant Molecular Biology Reporter.

[41]  Vijay Chandra Verma,et al.  Study of endophytic fungal community from different parts of Aegle marmelos Correae (Rutaceae) from Varanasi (India) , 2007 .

[42]  Priscila Chaverri,et al.  Species delimitation in fungal endophyte diversity studies and its implications in ecological and biogeographic inferences , 2011, Molecular ecology.

[43]  M. Hahn,et al.  Molecular and functional characterization of a secreted lipase from Botrytis cinerea. , 2005, Molecular plant pathology.

[44]  A. Salamov,et al.  Diverse Lifestyles and Strategies of Plant Pathogenesis Encoded in the Genomes of Eighteen Dothideomycetes Fungi , 2012, PLoS pathogens.

[45]  Seonju Lee,et al.  Diversity of saprobic hyphomycetes on Proteaceae and Restionaceae from South Africa , 2004 .

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

[47]  Neil D. Rawlings,et al.  MEROPS: the database of proteolytic enzymes, their substrates and inhibitors , 2013, Nucleic Acids Res..

[48]  S. Lumyong,et al.  A Phylogenetic Evaluation of Whether Endophytes Become Saprotrophs at Host Senescence , 2007, Microbial Ecology.

[49]  Jin-Rong Xu,et al.  Correction: Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi , 2014, BMC Genomics.

[50]  M. Delledonne,et al.  De novo genome assembly of the soil-borne fungus and tomato pathogen Pyrenochaeta lycopersici , 2014, BMC Genomics.

[51]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[52]  M. Morgante,et al.  The powdery mildew resistance gene REN1 co-segregates with an NBS-LRR gene cluster in two Central Asian grapevines , 2009, BMC Genetics.

[53]  F. de Lamotte,et al.  Characterization of a cassiicolin-encoding gene from Corynespora cassiicola, pathogen of rubber tree (Hevea brasiliensis). , 2012, Plant science : an international journal of experimental plant biology.

[54]  L. Lo Leggio,et al.  Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components , 2011, Proceedings of the National Academy of Sciences.

[55]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[56]  Brandi L. Cantarel,et al.  The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics , 2008, Nucleic Acids Res..

[57]  S. Mathivanan,et al.  Extracellular peptidases of the cereal pathogen Fusarium graminearum , 2015, Front. Plant Sci..

[58]  Neil Moore,et al.  Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses , 2012, Nature Genetics.

[59]  F. de Lamotte,et al.  Purification and characterization of cassiicolin, the toxin produced by Corynespora cassiicola, causal agent of the leaf fall disease of rubber tree. , 2007, Journal of chromatography. B, Analytical technologies in the biomedical and life sciences.

[60]  G. Felix,et al.  NPP1, a Phytophthora-associated trigger of plant defense in parsley and Arabidopsis. , 2002, The Plant journal : for cell and molecular biology.

[61]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[62]  Mei Wang,et al.  Chromone derivatives from a sponge-derived strain of the fungus Corynespora cassiicola. , 2015, Journal of natural products.

[63]  P. Carbonero,et al.  Softening-up mannan-rich cell walls. , 2012, Journal of experimental botany.

[64]  V. Pujade-Renaud,et al.  First characterization of endophytic Corynespora cassiicola isolates with variant cassiicolin genes recovered from rubber trees in Brazil , 2012, Fungal Diversity.

[65]  T. Lumley,et al.  gplots: Various R Programming Tools for Plotting Data , 2015 .

[66]  L. Dixon,et al.  Host specialization and phylogenetic diversity of Corynespora cassiicola. , 2009, Phytopathology.

[67]  J. van den Brink,et al.  Closely related fungi employ diverse enzymatic strategies to degrade plant biomass , 2015, Biotechnology for Biofuels.

[68]  C. Hsiao,et al.  Subcutaneous infection caused by Corynespora cassiicola, a plant pathogen. , 2010, The Journal of infection.

[69]  C. Eyles,et al.  A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by Zymoseptoria tritici , 2017, Fungal genetics and biology : FG & B.

[70]  Sébastien Lê,et al.  FactoMineR: An R Package for Multivariate Analysis , 2008 .

[71]  L. E. Gray,et al.  Isolation of the soybean pathogens Corynespora cassiicola and Phialophora gregata from cysts of Heterodera glycines in Illinois. , 1986 .

[72]  V. Chi,et al.  Genetic Diversity and Pathogenicity of Corynespora cassiicola Isolates from Rubber Trees and Other Hosts in Vietnam , 2014 .

[73]  C. Tsai,et al.  Purification and Characterization , 2006 .

[74]  E. Uberbacher,et al.  CAZymes Analysis Toolkit (CAT): web service for searching and analyzing carbohydrate-active enzymes in a newly sequenced organism using CAZy database. , 2010, Glycobiology.

[75]  Su Mei Yew,et al.  Dissecting the fungal biology of Bipolaris papendorfii: from phylogenetic to comparative genomic analysis , 2015, DNA research : an international journal for rapid publication of reports on genes and genomes.

[76]  P. Barthe,et al.  Structural analysis of cassiicolin, a host-selective protein toxin from Corynespora cassiicola. , 2007, Journal of molecular biology.

[77]  M. Wingfield,et al.  Botryosphaeriaceae as endophytes and latent pathogens of woody plants: diversity, ecology and impact , 2007 .

[78]  F. Bao,et al.  CARD9 mutation linked to Corynespora cassiicola infection in a Chinese patient , 2016, The British journal of dermatology.

[79]  Hanhong Bae,et al.  The diversity of fungal genome , 2015, Biological Procedures Online.

[80]  Wei Qian,et al.  Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. , 2000, Molecular biology and evolution.

[81]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[82]  Simon R. Harris,et al.  SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments , 2016, bioRxiv.

[83]  C. Stoeckert,et al.  OrthoMCL: identification of ortholog groups for eukaryotic genomes. , 2003, Genome research.

[84]  G. Platas,et al.  Geographical and seasonal influences on the distribution of fungal endophytes in Quercus ilex. , 1999, The New phytologist.

[85]  Inna Dubchak,et al.  MycoCosm portal: gearing up for 1000 fungal genomes , 2013, Nucleic Acids Res..

[86]  E. Mahgoub.  Corynespora cassiicola, a new agent of maduromycetoma. , 1969, The Journal of tropical medicine and hygiene.

[87]  R. Wijesundera,et al.  Genetic variation in Corynespora cassiicola: a possible relationship between host origin and virulence. , 2003, Mycological research.

[88]  J. V. van Kan,et al.  Cutinase A of Botrytis cinerea is expressed, but not essential, during penetration of gerbera and tomato. , 1997, Molecular plant-microbe interactions : MPMI.

[89]  P. Comménil,et al.  Antilipase antibodies prevent infection of tomato leaves byBotrytis cinerea , 1998 .

[90]  W. Friedt,et al.  Jasmonate and ethylene dependent defence gene expression and suppression of fungal virulence factors: two essential mechanisms of Fusarium head blight resistance in wheat? , 2012, BMC Genomics.

[91]  J. Dunwell.  Cupins: a new superfamily of functionally diverse proteins that include germins and plant storage proteins. , 1998, Biotechnology & genetic engineering reviews.

[92]  O. Gascuel,et al.  A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. , 2003, Systematic biology.

[93]  Anton J. Enright,et al.  An efficient algorithm for large-scale detection of protein families. , 2002, Nucleic acids research.

[94]  S. Basit.  Linkage analysis coupled with exome sequencing identified defects in gene ‘X’ causing premature ovarian insufficiency , 2014, BMC Genomics.

[95]  R. Sukumar,et al.  Endophytic fungal communities in woody perennials of three tropical forest types of the Western Ghats, southern India , 2011, Biodiversity and Conservation.

[96]  John S. Hu,et al.  Genome ‐ wide identification , classification and expression analysis in fungal – plant interactions of cutinase gene family and functional analysis of a putative ClCUT 7 in Curvularia lunata , 2016 .

[97]  J. Kadir,et al.  Morphological and Inter Simple Sequence Repeat (ISSR) Markers Analyses of Corynespora cassiicola Isolates from Rubber Plantations in Malaysia , 2008, Mycopathologia.

[98]  A. Covarrubias,et al.  Characterization of an Extracellular Serine Protease of Fusarium eumartii and its Action on Pathogenesis Related Proteins , 2004, European Journal of Plant Pathology.

[99]  A. Able,et al.  Secretome analysis of virulent Pyrenophora teres f. teres isolates , 2016, Proteomics.

[100]  Robert D. Finn,et al.  Dfam: a database of repetitive DNA based on profile hidden Markov models , 2012, Nucleic Acids Res..

[101]  Ole Tange,et al.  GNU Parallel: The Command-Line Power Tool , 2011, login Usenix Mag..

[102]  C. Jayasinghe,et al.  Corynespora leaf spot disease of rubber (Hevea brasiliensis) - a new record , 1986 .

[103]  N. L. Glass,et al.  Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. , 2014, Annual review of phytopathology.

[104]  Mira V. Han,et al.  Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. , 2013, Molecular biology and evolution.

[105]  John S. Hu,et al.  Genome-wide identification, classification and expression analysis in fungal–plant interactions of cutinase gene family and functional analysis of a putative ClCUT7 in Curvularia lunata , 2016, Molecular Genetics and Genomics.

[106]  C. Alabouvette,et al.  Fusarium oxysporum and its biocontrol. , 2003, The New phytologist.

[107]  Siliang Huang,et al.  Molecular and Pathogenic Variation Identified Among Isolates of Corynespora cassiicola , 2009, Molecular biotechnology.

[108]  C. Maurel,et al.  Insights into Populus XIP aquaporins: evolutionary expansion, protein functionality, and environmental regulation. , 2012, Journal of experimental botany.

[109]  G. Kingsland.  Pathogenicity and epidemiology of Corynespora cassiicola in the Republic of Seychelles , 1986 .

[110]  D. Bieysse,et al.  Diversity of the cassiicolin gene in Corynespora cassiicola and relation with the pathogenicity in Hevea brasiliensis. , 2014, Fungal biology.

[111]  Y. Hathout,et al.  Advances in the proteomic investigation of the cell secretome , 2012, Expert review of proteomics.