Abstract
Effectors are secreted by plant-associated microorganisms to modify the host cell physiology. As effectors, the Necrosis- and Ethylene-inducing peptide 1-like proteins (NLPs) are involded in the early phases of plant infection and may trigger host immune responses. Corynespora cassiicola is a polyphagous plant pathogen that causes target spot on many agriculturally important crops. Using genome assembly, gene prediction, and proteome annotation tools, we retrieved 135 NLP-encoding genes from proteomes of 44 isolates. We explored the evolutionary history of NLPs using Bayesian phylogeny, gene genealogies, and selection analyses. We accessed the expression profiles of the NLP genes during the early phase of C. cassiicola–soybean interaction. Three NLP putative-effector genes (Cc_NLP1.1, Cc_NLP1.2A, and Cc_NLP1.2B) were maintained in the genomes of all isolates tested. An NLP putative-non-effector gene (Cc_NLP1.3) was found in three isolates that had been originally obtained from soybean. Putative-effector NLPs were under different selective constraints: Cc_NLP1.1 was under stronger selective pressure, while Cc_NLP1.2A was under a more relaxed constraint. Meanwhile, Cc_NLP1.2B likely evolved under either positive or balancing selection. Despite highly divergent, the putative-effector NLPs maintain conserved the residues necessary to trigger plant immune responses, suggesting they are potentially functional. Only the Cc_NLP1.1 putative-effector gene was significantly expressed at the early hours of soybean colonization, while Cc_NLP1.2A and Cc_NLP1.2B showed much lower levels of gene expression.
Similar content being viewed by others
Data availability statement
The authors confirm that the data supporting the findings of this study are available in the supplementary material. Additional data are available from the corresponding author upon request.
References
Bae H, Kim MS, Sicher RC et al (2006) Necrosis- and Ethylene-Inducing Peptide from Fusarium oxysporum induces a complex cascade of transcripts associated with signal transduction and cell death in Arabidopsis. Plant Physiol 141:1056–1067. https://doi.org/10.1104/pp.106.076869.1056
Bailey TL, Boden M, Buske FA et al (2009) MEME Suite: tools for motif discovery and searching. Nucleic Acids Res 37:202–208. https://doi.org/10.1093/nar/gkp335
Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16:37–48
Bankevich A, Nurk S, Antipov D et al (2012) SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455–477. https://doi.org/10.1089/cmb.2012.0021
Bashi ZD, Hegedus DD, Buchwaldt L et al (2010) Expression and regulation of Sclerotinia sclerotiorum Necrosis and Ethylene-inducing peptides (NEPs). Mol Plant Pathol 11:43–53. https://doi.org/10.1111/j.1364-3703.2009.00571.x
Bhatt S, Katzourakis A, Pybus OG (2010) Detecting natural selection in RNA virus populations using sequence summary statistics. Infect Genet Evol 10:421–430. https://doi.org/10.1016/j.meegid.2009.06.001
Böhm H, Albert I, Oome S et al (2014) A conserved peptide pattern from a widespread microbial virulence factor triggers pattern-induced immunity in Arabidopsis. PLoS Pathog. https://doi.org/10.1371/journal.ppat.1004491
Bolger AM, Lohse M, Usadel B (2014) Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120. https://doi.org/10.1093/bioinformatics/btu170
Breton F, Sanier C, D’Auzac J (2000) Role of cassiicolin, a host-selective toxin in pathogenicity of C. cassiicola, causal agent of leaf disease of Hevea. J Nat Rubber Res 3:115–128
Buchfink B, Xie C, Huson DH (2014) Fast and sensitive protein alignment using DIAMOND. Nat Methods 12:59–60. https://doi.org/10.1038/nmeth.3176
Cabral A, Oome S, Sander N et al (2012) Nontoxic Nep1-like proteins of the downy mildew pathogen Hyaloperonospora arabidopsidis: repression of necrosis-inducing activity by a surface-exposed region. Mol Plant-Microbe Interact 25:697–708. https://doi.org/10.1094/MPMI-10-11-0269
Dal’Sasso TCS, Rody HVS, Grijalba PE, Oliveira LO, (2021) Genome sequences and in silico effector mining of Corynespora cassiicola CC_29 and Corynespora olivacea CBS 114450. Arch Microbiol 203:5257–5265. https://doi.org/10.1007/s00203-021-02456-7
Dal’Sasso TCS, Rody HVS, Oliveira LO (2022) Genome-wide analysis and evolutionary history of the Necrosis and Ethylene-inducing peptide 1-like protein (NLP) superfamily across the Dothideomycetes class of fungi. bioRxiv. https://doi.org/10.1101/2022.02.07.479250
Darriba D, Taboada GL, Doallo R, Posada D (2012) jModelTest 2: more models, new heuristics and high-performance computing Europe PMC Funders Group. Nat Methods 9:772. https://doi.org/10.1038/nmeth.2109.jModelTest
de Lamotte F, Duviau MP, Sanier C et al (2007) Purification and characterization of cassiicolin, the toxin produced by Corynespora cassiicola, causal agent of the leaf fall disease of rubber tree. J Chromatogr B Anal Technol Biomed Life Sci 849:357–362. https://doi.org/10.1016/j.jchromb.2006.10.051
Déon M, Scomparin A, Tixier A et al (2012) First characterization of endophytic Corynespora cassiicola isolates with variant cassiicolin genes recovered from rubber trees in Brazil. Fungal Divers 54:87–99. https://doi.org/10.1007/s13225-012-0169-6
Déon M, Fumanal B, Gimenez S et al (2014) Diversity of the cassiicolin gene in Corynespora cassiicola and relation with the pathogenicity in Hevea brasiliensis. Fungal Biol 118:32–47. https://doi.org/10.1016/j.funbio.2013.10.011
Dixon LJ, Schlub RL, Pernezny K, Datnoff LE (2009) Host specialization and phylogenetic diversity of Corynespora cassiicola. Phytopathology 99:1015–1027. https://doi.org/10.1094/PHYTO-99-9-1015
Dong S, Kong G, Qutob D et al (2012) The NLP toxin family in Phytophthora sojae includes rapidly evolving groups that lack necrosis-inducing activity. Mol Plant-Microbe Interact 25:896–909. https://doi.org/10.1094/MPMI-01-12-0023-R
Duhan D, Gajbhiye S, Jaswal R et al (2021) Functional characterization of the Nep1-Like protein effectors of the necrotrophic pathogen – Alternaria brassicae. Front Microbiol 12:1–11. https://doi.org/10.3389/fmicb.2021.738617
El-Gebali S, Mistry J, Bateman A et al (2019) The Pfam protein families database in 2019. Nucleic Acids Res 47:D427–D432. https://doi.org/10.1093/nar/gky995
Emanuelsson O, Brunak S, von Heijne G, Nielsen H (2007) Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc 2:953–971. https://doi.org/10.1038/nprot.2007.131
Emms DM, Kelly S (2015) OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol 16:1–14. https://doi.org/10.1186/s13059-015-0721-2
Fang YL, Peng YL, Fan J (2017) The Nep1-like protein family of Magnaporthe oryzae is dispensable for the infection of rice plants. Sci Rep 7:1–10. https://doi.org/10.1038/s41598-017-04430-0
Fellbrich G, Romanski A, Varet A et al (2002) NPP1, a Phytophthora-associated trigger of plant defense in parsley and Arabidopsis. Plant J 32:375–390. https://doi.org/10.1046/j.1365-313X.2002.01454.x
Fouché S, Plissonneau C, Croll D (2018) The birth and death of effectors in rapidly evolving filamentous pathogen genomes. Curr Opin Microbiol 46:34–42. https://doi.org/10.1016/j.mib.2018.01.020
Franceschetti M, Maqbool A, Jiménez-Dalmaroni MJ et al (2017) Effectors of filamentous plant pathogens: commonalities amid diversity. Microbiol Mol Biol Rev 81:1–17. https://doi.org/10.1128/mmbr.00066-16
Fu YX, Li WH (1993) Statistical tests of neutrality of mutations. Genetics 133:693 LP – 709
Gao S, Zeng R, Xu L et al (2020) Genome sequence and spore germination-associated transcriptome analysis of Corynespora cassiicola from cucumber. BMC Microbiol 20:1–20. https://doi.org/10.1186/s12866-020-01873-w
Gijzen M, Nürnberger T (2006) Nep1-like proteins from plant pathogens: recruitment and diversification of the NPP1 domain across taxa. Phytochemistry 67:1800–1807. https://doi.org/10.1016/j.phytochem.2005.12.008
Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by Maximum Likelihood. Syst Biol 52:696–704. https://doi.org/10.1080/10635150390235520
Irieda H, Inoue Y, Mori M et al (2019) Conserved fungal effector suppresses PAMP-triggered immunity by targeting plant immune kinases. Proc Natl Acad Sci U S A 116:496–505. https://doi.org/10.1073/pnas.1807297116
Jones P, Binns D, Chang HY et al (2014) InterProScan 5: Genome-scale protein function classification. Bioinformatics 30:1236–1240. https://doi.org/10.1093/bioinformatics/btu031
Jones DA, Bertazzoni S, Turo CJ et al (2018) Bioinformatic prediction of plant–pathogenicity effector proteins of fungi. Curr Opin Microbiol 46:43–49. https://doi.org/10.1016/j.mib.2018.01.017
Katoh K, Standley DM (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol 30:772–780. https://doi.org/10.1093/molbev/mst010
Krogh A, Larsson B, Von Heijne G, Sonnhammer ELL (2001) Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol 305:567–580. https://doi.org/10.1006/jmbi.2000.4315
Lenarčič T, Albert I, Böhm H et al (2017) Eudicot plant-specific sphingolipids determine host selectivity of microbial NLP cytolysins. Science (80-) 358:1431–1434
Levin E, Raphael G, Ma J et al (2019) Identification and functional analysis of NLP-encoding genes from the postharvest pathogen Penicillium expansum. Microorganisms. https://doi.org/10.3390/microorganisms7060175
Liu L, Xu L, Jia Q et al (2019) Arms race: diverse effector proteins with conserved motifs. Plant Signal Behav. https://doi.org/10.1080/15592324.2018.1557008
Lo Presti L, Lanver D, Schweizer G et al (2015) Fungal effectors and plant susceptibility. Annu Rev Plant Biol 66:513–545. https://doi.org/10.1146/annurev-arplant-043014-114623
Looi HK, Toh YF, Yew SM et al (2017) Genomic insight into pathogenicity of dematiaceous fungus Corynespora cassiicola. PeerJ 5:1–28. https://doi.org/10.7717/peerj.2841
Lopez D, Ribeiro S, Label P et al (2018) Genome-wide analysis of Corynespora cassiicola leaf fall disease putative effectors. Front Microbiol 9:1–21. https://doi.org/10.3389/fmicb.2018.00276
McDonald J, Kreitman H (1991) Adaptive protein evolution at the Adh locus in Drosophila. Nature 351:652–654
Menardo F, Praz CR, Wicker T, Keller B (2017) Rapid turnover of effectors in grass powdery mildew (Blumeria graminis). BMC Evol Biol 17:1–14. https://doi.org/10.1186/s12862-017-1064-2
Mikheenko A, Prjibelski A, Saveliev V et al (2018) Versatile genome assembly evaluation with QUAST-LG. Bioinformatics 34:i142–i150. https://doi.org/10.1093/bioinformatics/bty266
Ngo KX, Nguyen PDN, Furusho H et al (2022) Unraveling the host-selective toxic interaction of cassiicolin with lipid membranes and its cytotoxicity. Phytopathology. https://doi.org/10.1094/PHYTO-09-21-0397-R
Ono E, Mise K, Takano Y (2020) RLP23 is required for Arabidopsis immunity against the grey mould pathogen Botrytis cinerea. Sci Rep 10:1–12. https://doi.org/10.1038/s41598-020-70485-1
Oome S, Van Den Ackerveken G (2014) Comparative and functional analysis of the widely occurring family of Nep1-like proteins. Mol Plant-Microbe Interact 27:1081–1094. https://doi.org/10.1094/MPMI-04-14-0118-R
Oome S, Raaymakers TM, Cabral A et al (2014) Nep1-like proteins from three kingdoms of life act as a microbe-associated molecular pattern in Arabidopsis. Proc Natl Acad Sci U S A 111:16955–16960. https://doi.org/10.1073/pnas.1410031111
Ottmann C, Luberacki B, Küfner I et al (2009) A common toxin fold mediates microbial attack and plant defense. Proc Natl Acad Sci U S A 106:10359–10364. https://doi.org/10.1073/pnas.0902362106
Petersen TN, Brunak S, Von Heijne G, Nielsen H (2011) SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods 8:785–786. https://doi.org/10.1038/nmeth.1701
Pour FN, Cobos R, Coque JJR et al (2020) Toxicity of recombinant necrosis and ethylene-inducing proteins (NLPs) from Neofusicoccum parvum. Toxins (basel) 12:1–17. https://doi.org/10.3390/toxins12040235
Qutob D, Kemmerling B, Brunner F et al (2006) Phytotoxicity and innate immune responses induced by Nep1-like proteins. Plant Cell 18:3721–3744. https://doi.org/10.1105/tpc.106.044180
Rambaut A, Drummond AJ, Xie D et al (2018) Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol 67:901–904. https://doi.org/10.1093/sysbio/syy032
Ramírez-Soriano A, Ramos-Onsins SE, Rozas J et al (2008) Statistical power analysis of neutrality tests under demographic expansions, contractions and bottlenecks with recombination. Genetics 179:555–567. https://doi.org/10.1534/genetics.107.083006
Rody HVS, Oliveira LO (2018) Evolutionary history of the cobalamin-independent methionine synthase gene family across the land plants. Mol Phylogenet Evol 120:33–42. https://doi.org/10.1016/j.ympev.2017.12.003
Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572–1574
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC et al (2017) DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol 34:3299–3302. https://doi.org/10.1093/molbev/msx248
Sánchez-Vallet A, Fouché S, Fudal I et al (2018) The genome biology of effector gene evolution in filamentous plant pathogens. Annu Rev Phytopathol. https://doi.org/10.1146/annurev-phyto-080516-035303
Seidl MF, Faino L, Shi-Kunne X et al (2015) The genome of the saprophytic fungus Verticillium tricorpus reveals a complex effector repertoire resembling that of its pathogenic relatives. Mol Plant-Microbe Interact 28:362–373. https://doi.org/10.1094/MPMI-06-14-0173-R
Seidl MF, Van den Ackerveken G (2019) Activity and phylogenetics of the broadly occurring family of microbial Nep1-like proteins. Annu Rev Phytopathol 57:367–386. https://doi.org/10.1146/annurev-phyto-082718-100054
Shrestha SK, Lamour K, Young-Kelly H (2018) Genome sequences and SNP analyses of Corynespora cassiicola from cotton and soybean in the southeastern United States reveal limited diversity. PLoS ONE 12:6–14. https://doi.org/10.1371/journal.pone.0184908
Simão FA, Waterhouse RM, Ioannidis P et al (2015) BUSCO: Assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31:3210–3212. https://doi.org/10.1093/bioinformatics/btv351
Sonah H, Deshmukh RK, Bélanger RR (2016) Computational prediction of effector proteins in fungi: opportunities and challenges. Front Plant Sci 7:1–14. https://doi.org/10.3389/fpls.2016.00126
Sperschneider J, Dodds PN (2022) EffectorP 3.0: prediction of apoplastic and cytoplasmic effectors in fungi and oomycetes. Mol Plant Microbe Interact 35:146–156. https://doi.org/10.1094/MPMI-08-21-0201-R
Staats M, van Baarlen P, Schouten A et al (2007) Positive selection in phytotoxic protein-encoding genes of Botrytis species. Fungal Genet Biol 44:52–63. https://doi.org/10.1016/j.fgb.2006.07.003
Stanke M, Morgenstern B (2005) AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res 33:465–467. https://doi.org/10.1093/nar/gki458
Sumabat LG, Kemerait RC, Brewer MT (2018) Phylogenetic diversity and host specialization of Corynespora cassiicola responsible for emerging target spot disease of cotton and other crops in the southeastern United States. Phytopathology 108:892–901. https://doi.org/10.1094/PHYTO-12-17-0407-R
Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595
Terauchi R, Yoshida K (2010) Towards population genomics of effector-effector target interactions. New Phytol 187:929–939. https://doi.org/10.1111/j.1469-8137.2010.03408.x
Wu J, Xie X, Shi Y et al (2018) Variation of cassiicolin genes among Chinese isolates of Corynespora cassiicola. J Microbiol 56:634–647. https://doi.org/10.1007/s12275-018-7497-5
Yang Z (2007) PAML 4: phylogenetic analysis by Maximum likelihood. Mol Biol Evol 24:1586–1591
Acknowledgements
We thank Dr. Eveline Caixeta and the Laboratório de Biotecnologia do Cafeeiro for the use of the Real-Time PCR System. This work was supported by The Minas Gerais State Foundation of Research Aid—FAPEMIG (grant number APQ-00150-17) and by The National Council of Scientific and Technological Development – CNPq (fellowship number PQ 302336/2019-2) to LOO. TCSD received student fellowships from the CAPES Foundation (PROEX—0487 No. 1684083) and CNPq (GM/GD 142400/2018-1).
Author information
Authors and Affiliations
Contributions
TCSD performed data assembly and analyses. TCSD, HVSR, and LOO contributed to the study conception. VDR and MDLC contributed to the gene expression analyses. TCSD and LOO wrote the manuscript. LOO and HVSR supervised the research. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Conflict of interest
The authors have declared that no competing interests exist.
Additional information
Communicated by M. Polymenis.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
294_2022_1252_MOESM1_ESM.pdf
Fig. S1. Sliding window neutrality tests of three predicted Necrosis- and Ethylene-inducing peptide 1-like protein (NLP)-encoding genes (CDS) of Corynespora cassiicola. The x axis shows the midpoints of contiguous windows (25 bases, with steps of 5 bases). The gray vertical bar indicates the relative position of the GHRHDWE motif. The arrow indicates sites with significant departures from zero (p-value < 0.05)
294_2022_1252_MOESM2_ESM.pdf
Fig. S2. Placement of nucleotide substitutions along the haplotypes of each of the putative-effector three NLP-encoding genes (Cc_NLP1.1, Cc_NLP1.2A, and Cc_NLP1.2B) (CDS) in Corynespora cassiicola. The putative-effector NLP-encoding genes from isolate CC_29 were used as reference for haplotype sequence alignment, according to the haplogroups in Fig. 2. Upper protein representations indicate portions of the CDS that encode for the signal peptide (SP), the NPP1 domain, and the conserved GHRHDWE motif in the putative-effector NLPs of isolate CC_29. A. Cc_NLP1.1 compared to CC_29_g12228 (haplotype 3). B. Cc_NLP1.2A compared to CC_29_g13117 (haplotype 15). C. Cc_NLP1.2B compared to CC_29_g13473 (haplotype 8). Plus (+) signs indicate substitutions that resulted in amino acids replacements predicted to be positively selected, according to codeml/PAML analysis. Types of nucleotide substitutions as indicated
294_2022_1252_MOESM3_ESM.pdf
Fig. S3. Expression patterns of the Necrosis- and Ethylene-inducing peptide 1-like (NLP) genes of Corynespora cassiicola (isolate CC_29) after inoculation in soybean plants. Gene expression was inferred by RT-qPCR analyses. The time (x axis) is given in hours post inoculation (hpi). The fold variation of gene expression of each NLP gene (y axis) was calculated using the 2-ΔCtmethod, the values were normalized to an endogenous control gene (β-tubulin). Letters show significant differences (p-value < 0.05) after a one-way ANOVA followed by Tukey’s test for each NLP gene individually. Error bars represent the standard deviation of three biological replicates
Rights and permissions
Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Dal’Sasso, T.C.d., Rocha, V.D., Rody, H.V.S. et al. The necrosis- and ethylene-inducing peptide 1-like protein (NLP) gene family of the plant pathogen Corynespora cassiicola. Curr Genet 68, 645–659 (2022). https://doi.org/10.1007/s00294-022-01252-0
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00294-022-01252-0