Chitosan induces differential transcript usage of chitosanase 3 encoding gene (csn3) in the biocontrol fungus Pochonia chlamydosporia 123
BMC Genomics volume 23, Article number: 101 (2022)
Pochonia chlamydosporia is an endophytic fungus used for nematode biocontrol that employs its cellular and molecular machinery to degrade the nematode egg-shell. Chitosanases, among other enzymes, are involved in this process. In this study, we improve the genome sequence assembly of P. chlamydosporia 123, by utilizing long Pacific Biosciences (PacBio) sequence reads. Combining this improved genome assembly with previous RNA-seq data revealed alternative isoforms of a chitosanase in the presence of chitosan. This study could open new insights into understanding fungal resistance to chitosan and root-knot nematode (RKN) egg infection processes.
The P. chlamydosporia 123 genome sequence assembly has been updated using long-read PacBio sequencing and now includes 12,810 predicted protein-coding genes. Compared with the previous assembly based on short reads, there are 701 newly annotated genes, and 69 previous genes are now split. Eight of the new genes were differentially expressed in fungus interactions with Meloidogyne javanica eggs or chitosan.
A survey of the RNA-seq data revealed alternative splicing in the csn3 gene that encodes a chitosanase, with four putative splicing variants: csn3_v1, csn3_v2, csn3_v3 and csn3_v4. When P. chlamydosporia is treated with 0.1 mg·mL− 1 chitosan for 4 days, csn3 is expressed 10-fold compared with untreated controls. Furthermore, the relative abundances of each of the four transcripts are different in chitosan treatment compared with controls. In controls, the abundances of each transcript are nil, 32, 55, and 12% for isoforms csn3_v1, csn3_v2, csn3_v3 and csn3_v4 respectively. Conversely, in chitosan-treated P. chlamydosporia, the abundances are respectively 80, 15%, 2—3%, 2—3%. Since isoform csn3_v1 is expressed with chitosan only, the putatively encoded enzyme is probably induced and likely important for chitosan degradation.
Alternative splicing events have been discovered and described in the chitosanase 3 encoding gene from P. chlamydosporia 123. Gene csn3 takes part in RKN parasitism process and chitosan enhances its expression. The isoform csn3_v1 would be related to the degradation of this polymer in bulk form, while other isoforms may be related to the degradation of chitosan in the nematode egg-shell.
Pochonia chlamydosporia (Goddard) Zare and Gams (Pc) is an endophytic fungus used for biocontrol of nematode eggs and females from plant parasitic nematodes, including root-knot nematodes (RKN) as Meloidogyne javanica  and M. incognita  and cyst nematodes such as Globodera spp. , among others. When infecting M. javanica eggs, Pc generates chitosan from chitin using chitin deacetylases . Chitosan functional importance is reflected in the observation that Pc has more genes encoding enzymes related to chitin and chitosan degradation than any other sequenced fungus [4, 17, 18]. Chitosanases (EC 220.127.116.11, glycoside hydrolase 75, www.cazy.org) hydrolyze chitosan to oligosaccharides . They have been detected in plants , where they have been considered as defense enzymes against pathogens . Chitosanolytic enzymes are also present in soil microorganisms [1, 11, 19, 34]; most of them are chitosan-resistant, including Pc . Chitosan enhances appressorium differentiation and RKN egg parasitism by Pc . Genes encoding chitosanases are expressed by Pc on its own, with RKN eggs or with chitosan, but mainly with the fungus, RKN eggs and chitosan together [4, 29].
Higher eukaryotic organisms have evolved mechanisms to increase the variability of proteins that are synthesized by a cell. Alternative splicing consists of the generation of different transcripts from the same single DNA strand, resulting in proteins of different conformation and length, which usually have different activities . Alternative splicing patterns can sometimes be detected in transcriptomic data, such as RNA-seq . In a system with multiple experimental conditions, a distinction can be made between differential transcript expression (DTE) and differential transcript usage (DTU ;). In DTE, it is possible to observe expression changes for at least one transcript between conditions. This implies gene overexpression or repression. In DTU, relative expression level of each isoform varies with experimental conditions. DTU implies DTE, but not vice versa . DTU analyses can reveal genes that express different isoforms under given conditions. This could be related to environmental adaptation. Although alternative splicing has been found in fungi , it has not been widely studied. In this work we report the alternative splicing and gene expression patterns of a chitosanase in Pc isolate 123 growing with chitosan.
Materials and methods
Pochonia chlamydosporia var. chlamydosporia (=Metacordyceps chlamydosporia var. chlamydosporium) isolate 123 (Pc123) (ATCC No. MYA-4875; CECT No. 20929) was selected for genome resequencing and alternative splicing analysis. Pc123 was isolated from Heterodera avenae infected eggs  in south-west Spain.
DNA isolation and Pacific biosciences sequencing
Pc123 conidia (final concentration 106 conidia·mL− 1) were inoculated into 250 mL flasks each containing 50 mL Potato Dextrose Broth medium (24 g·L− 1). Flasks were incubated at 25 °C with shaking at 120 rpm. After 5 days, mycelia were recovered by filtration through Miracloth (Calbiochem) and washed twice with sterile distilled water (SDW). DNA from Pc123 resulting fresh mycelia (ca. 0.5 g) was extracted using DNeasy Plant Mini Kit (Qiagen) following manufacturer’s instructions.
Pc123 DNA was sent to Macrogen Inc. to perform Pacific Biosciences (PacBio) sequencing. PacBio Sequel SMRT (20 Kb insert size) was used as library, with PacBio Sequel SMRT 1 cell Run as sequencing platform. Throughput was around 6–7 Gb/spl.
Genome sequence assembly and gene-calling
We used an assembly strategy that combined de novo assembly of the new, long reads with the previously assembled and annotated genome sequence for Pc 123 (GenBank: GCA_000411695.2). First, we refined the original Pc123 annotation using RNA-seq data from Suarez-Fernandez et al.  with Program to Assemble Spliced Alignments (PASA; Haas et al. ). PacBio data were then assembled de novo with the long-read sequence assembler Canu . Both assemblies were combined using RagTag, a tool for reference-guided genome assembly improvement that allows current annotation features to be preserved and updated for the new reference . The full details of command lines, with parameter values and options, can be found in the Supplementary Information. BLASTx (NCBI) with standard genetic code and non-redundant protein sequences as database was used to identify homologies of selected novel genes of interest with other organisms.
Genome annotation and alternative splicing analysis
PASA was then re-run on the new assembly to create a new annotation with updated gene models, whilst retaining all previous annotation information. Format conversion, data tidying and script preparation were done for alternative transcript usage testing.
Predicted polypeptide and transcript sequences were analysed using Pfam to detect conserved domains , SignalP 5.0 to detect secretion signals , CPC2 to assess coding potential  and NetSurfP-2 to detect intrinsically disordered regions (IDR) . The results of these searches were combined to analyse alternative transcript usage using IsoformSwitchAnalyzeR [32, 33], which enables identification and analysis of alternative splicing and isoform switches from RNA-seq data. Full details can be found in the Supplementary Information. PacBio and RNA-seq data have been submitted to NCBI Sequence Read Archive (BioProjects PRJNA68669 and PRJNA741387, respectively).
Identification of differentially spliced genes
After genome resequencing, RNA-seq data from Suarez-Fernandez et al.  were mapped against the updated genome using Salmon with a wrapper script (align_and_estimate_abundance.pl) from the Trinity software package . Full details are provided in the Supplementary Information.
Suarez-Fernandez et al.  determined the transcriptomic effect of chitosan on Pc123 root-knot nematode parasitism. We selected Pc (control, Pc growing in minimal medium for 4 days) and PcQ (Pc growing in minimal medium amended with 0,1 mg·mL− 1 chitosan for 4 days) treatments from that experiment.
PacBio sequencing reduces number of Pc123 scaffolds and predicts additional putative genes
PacBio reads (average length: 15.6 kb) were used to improve the previous assembly of the Pochonia chlamydosporia 123 (Pc123) genome  (GenBank: GCA_000411695.2). The new version of the assembly is deposited as GenBank accession GCA_000411695.4.
We transferred annotation, facilitating comparisons to previous experiments. Gene models were updated using transcriptomic data . These include 20 original gene models that were merged to create 10 new genes and 69 genes that were split. Furthermore, 701 putative genes that were not previously detected have also been identified, 499 of them non-overlapping with current models. After PacBio sequencing, it was possible to reduce the number of scaffolds from 956 to 121 and that of contigs from 9087 to 8409. Finally, after this new sequencing, 12,721 genes are now predicted, which is 770 more than in the previous Pc123 genome prediction (AOSW02000000). Pc123 resequencing and genome improvement is summarized in Table 1. PacBio data have been submitted to NCBI Sequence Read Archive database with accession number SRR14907880 (BioProject PRJNA68669).
Transcriptomic data (PRJNA741387) identifies that eight of the newly identified novel genes (novel_gene_495_5ed78ef1, novel_gene_946_5ed78ef1, novel_gene_431_5ed78ef1, novel_gene_491_5ed78ef1, novel_gene_506_5ed78ef1, novel_gene_82_5ed78ef1, novel_gene_117_5ed78ef1, and novel_gene_303_5ed78ef1) are significantly differentially expressed in at least one RNA-seq treatment; three of these are non-overlapping with current models. All of them are homologous to previously known hypothetical proteins from P. chlamydosporia, Metarhizium anisopliae or Ustilaginoidea virens except novel_gene_491_5ed78ef1, for which no significant sequence similarities can be found (Table 2). None of the split or the merged genes are significantly differentially expressed in any comparison.
Chitosan stimulates the expression of an isoform of Pc123 cns3 gene
Chitosan induces alternative transcript usage in Pc123 (Fig. 1). We have found at least 20 alternatively spliced transcripts significantly expressed in Pc123 treated with chitosan (Table 3). Pc123 locus-tag I1G_00010429 or csn3 (GenBank: RZR62940.1), which encodes chitosanase 3 , shows alternative transcript usage when the fungus is treated with 0.1 mg·mL− 1 chitosan for 4 days. We have found four isoforms for csn3 (Fig. 1A): rna-gnl_WGS:AOSW:I1G_00010429-RA_mrna, rna-gnl_WGS:AOSW:I1G_00010429-RA_mrna.1.5ed7a624, rna-gnl_WGS:AOSW:I1G_00010429-RA_mrna.1.5ed7a624.2.5ee8159f and rna-gnl_WGS:AOSW:I1G_00010429-RA_mrna.1.5ed7a624.3.5ee8159f. We named these isoforms as csn3_v1, csn3_v2, csn3_v3 and csn3_v4, respectively.
In the chitosan-treated fungus Pc123 growing with chitosan, the total expression (all isoforms) of csn3 is 10-fold higher than the total expression in the control treatment without chitosan (Fig. 1B). Adding chitosan to the medium increases csn3_v1 isoform expression by almost 10-fold (Fig. 1C). Relative expression levels of csn3 isoforms (Fig. 1D) in controls are ca. 55% for csn3_v3, ca. 32% for csn3_v2 and ca. 12% for csn3_v4, while for csn3_v1 it is nil. In the chitosan-treated fungus, csn3_v1 isoform represents ca. 80% of gene expression, csn3_v2 ca. 15%, and csn3_v3 and csn3_v4 ca. 2–3%. This means that constitutive expression of isoform csn3_v1 is low or nil in Pc123. However, it becomes the most expressed isoform when chitosan is present in the growth medium of Pc123, suggesting that this isoform is the most efficient at degrading chitosan to chitooligosaccharides, offering a testable hypothesis for future studies.
In this work, we have found that chitosan induces alternative splicing events in csn3 from Pochonia chlamydosporia 123. Alternative splicing occurs naturally in fungi [9, 26]. This is an important emerging issue in the regulation of fungal gene expression. Previous work demonstrates that alternative splicing events are present during root colonization by arbuscular mycorrhizal fungi  and plant infection by Sclerotinia sclerotiorum . Therefore, fungi activate alternative splicing processes under different conditions to adapt to a changing environment. This suggests alternative splicing events may be related to epigenetics  and it may be the environment that determines relative levels of transcript expression. It has been shown that fungi, as well as animals and plants, are highly dependent on epigenetics . Therefore, studying fungal alternative splicing is a promising avenue for future studies related to the effect of environment on gene expression. Understanding the mechanism by which a fungus generates a series of transcripts from a single DNA molecule could help to unravel how it responds to a stimulus. Thus, not taking into account transcriptional variants in RNA-seq analyses may lead to loss of knowledge of alternative transcript usage in key genes. This could lead to incomplete conclusions. RNA-seq analyses should therefore consider splicing variants . Previous studies of alternative splicing in chitosanases demonstrate that differential expression, as in the case of csn3, is a common event . Based on that work, we believe that csn3 isoforms could have different functions or even locations . On the other hand, Pochonia chlamydosporia 123 (Pc123) is known to have high chitosanolytic activity  due to its high content in chitosanases . Future work will extend alternative transcript usage analyses to the rest of Pc123 chitosanases encoded in its genome in order to determine whether they also undergo alternative splicing events. Aranda-Martinez et al.  show that csn3 is induced six-fold during Pc123 RKN parasitism. Besides, expression value of csn3 when Pc123 infects RKN eggs in a medium amended with chitosan almost doubles the value respect to only-chitosan treatment . This suggests csn3 is one of the key genes that take part in RKN parasitism and chitosan enhances its expression. The chitosan-promoting isoform may be related to the degradation of this polymer in bulk form, while other isoforms may be related to the degradation of chitosan in the fungal or RKN egg wall. This could open new insights into understanding fungal resistance to chitosan and RKN egg infection processes.
Availability of data and materials
PacBio and RNA-seq data have been submitted to NCBI SRA database (BioProjects PRJNA68669 and PRJNA741387, respectively). Both can be accessed following these links:
Aktuganov GE, Safina VR, Galimzianova NF, Kuz’mina LY, Gilvanova EA, Boyko TF, et al. Chitosan resistance of Bacteria and Micromycetes differing in ability to produce extracellular Chitinases and Chitosanases. Microbiology. 2018;87:716–24.
Almagro Armenteros JJ, Tsirigos KD, Sønderby CK, Petersen TN, Winther O, Brunak S, et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37:420–3.
Alonge M, Soyk S, Ramakrishnan S, Wang X, Goodwin S, Sedlazeck FJ, et al. RaGOO: fast and accurate reference-guided scaffolding of draft genomes. Genome Biol. 2019;20:224.
Aranda-Martinez A, Lenfant N, Escudero N, Zavala-Gonzalez EA, Henrissat B, Lopez-Llorca LV. CAZyme content of Pochonia chlamydosporia reflects that chitin and chitosan modification are involved in nematode parasitism: CAZome of Pochonia chlamydosporia. Environ Microbiol. 2016;18:4200–15.
Berlato C, Chan KV, Price AM, Canosa M, Scibetta AG, Hurst HC. Alternative TFAP2A isoforms have distinct activities in breast cancer. Breast Cancer Res. 2011;13:R23.
Escudero N, Lopez-Moya F, Ghahremani Z, Zavala-Gonzalez EA, Alaguero-Cordovilla A, Ros-Ibañez C, et al. Chitosan increases tomato root colonization by Pochonia chlamydosporia and their combination reduces root-knot nematode damage. Front Plant Sci. 2017;8:1415.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Grenier J, Asselin A. Some pathogenesis-related proteins are Chitosanases with lytic activity against fungal spores. MPMI. 1990;3:401.
Grutzmann K, Szafranski K, Pohl M, Voigt K, Petzold A, Schuster S. Fungal alternative splicing is associated with multicellular complexity and virulence: a genome-wide multi-species study. DNA Res. 2014;21:27–39.
Haas BJ. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003;31:5654–66.
Hirano K, Arayaveerasid S, Seki K, Adams DJ, Mitsutomi M. Characterization of a Chitosanase from Aspergillus fumigatus ATCC13073. Biosci Biotechnol Biochem. 2012;76:1523–8.
Ibrahim HMM, Kusch S, Didelon M, Raffaele S. Genome-wide alternative splicing profiling in the fungal plant pathogen Sclerotinia sclerotiorum during the colonization of diverse host families. Mol Plant Pathol. 2021;22:31–47.
Kamerling JP, Boons G-J. Comprehensive glycoscience: from chemistry to systems biology. Amsterdam; Boston: Elsevier; 2007.
Kang Y-J, Yang D-C, Kong L, Hou M, Meng Y-Q, Wei L, et al. CPC2: a fast and accurate coding potential calculator based on sequence intrinsic features. Nucleic Acids Res. 2017;45:W12–6.
Klausen MS, Jespersen MC, Nielsen H, Jensen KK, Jurtz VI, Sønderby CK, et al. NetSurfP-2.0: improved prediction of protein structural features by integrated deep learning. Proteins. 2019;87:520–7.
Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. Canu: scalable and accurate long-read assembly via adaptive k -mer weighting and repeat separation. Genome Res. 2017;27:722–36.
Larriba E, Jaime MDLA, Carbonell-Caballero J, Conesa A, Dopazo J, Nislow C, et al. Sequencing and functional analysis of the genome of a nematode egg-parasitic fungus, Pochonia chlamydosporia. Fungal Genet Biol. 2014;65:69–80.
Lin R, Qin F, Shen B, Shi Q, Liu C, Zhang X, et al. Genome and secretome analysis of Pochonia chlamydosporia provide new insight into egg-parasitic mechanisms. Sci Rep. 2018;8:1123.
Liu Y, Li Y, Tong S, Yuan M, Wang X, Wang J, et al. Expression of a Beauveria bassiana chitosanase (BbCSN-1) in Pichia pastoris and enzymatic analysis of the recombinant protein. Protein Expr Purif. 2020;166:105519.
Luco RF, Allo M, Schor IE, Kornblihtt AR, Misteli T. Epigenetics in alternative pre-mRNA splicing. Cell. 2011;144:16–26.
Madhani HD. Unbelievable but true: epigenetics and chromatin in Fungi. Trends Genet. 2021;37:12–20.
Mistry J, Chuguransky S, Williams L, Qureshi M, Salazar GA, Sonnhammer ELL, et al. Pfam: the protein families database in 2021. Nucleic Acids Res. 2021;49:D412–9.
Olivares CM, López-Llorca LV. Fungal egg-parasites of plant-parasitic nematodes from Spanish soils. Rev Iberoam Micol. 2002;19:104–10.
Palma-Guerrero J, Gómez-Vidal S, Tikhonov VE, Salinas J, Jansson H-B, Lopez-Llorca LV. Comparative analysis of extracellular proteins from Pochonia chlamydosporia grown with chitosan or chitin as main carbon and nitrogen sources. Enzym Microb Technol. 2010;46:568–74.
Park JW, Tokheim C, Shen S, Xing Y. Identifying differential alternative splicing events from RNA sequencing data using RNASeq-MATS. In: Shomron N, editor. Deep sequencing data analysis. Totowa: Humana Press; 2013. p. 171–9. https://doi.org/10.1007/978-1-62703-514-9_10.
Sieber P, Voigt K, Kämmer P, Brunke S, Schuster S, Linde J. Comparative study on alternative splicing in human fungal pathogens suggests its involvement during host invasion. Front Microbiol. 2018;9:2313.
Soneson C, Matthes KL, Nowicka M, Law CW, Robinson MD. Isoform prefiltering improves performance of count-based methods for analysis of differential transcript usage. Genome Biol. 2016;17:12.
Stephan-Otto Attolini C, Peña V, Rossell D. Designing alternative splicing RNA-seq studies. Beyond generic guidelines. Bioinformatics. 2015;31:3631–7.
Suarez-Fernandez M, Sambles C, Lopez-Moya F, Nueda MJ, Studholme DJ, Lopez-Llorca LV. Chitosan modulates Pochonia chlamydosporia gene expression during nematode egg parasitism. Environ Microbiol. 2021;23:4980–97. https://doi.org/10.1111/1462-2920.15408.
Thadathil N, Velappan SP. Recent developments in chitosanase research and its biotechnological applications: a review. Food Chem. 2014;150:392–9.
Tobin JD, Haydock PPJ, Hare MC, Woods SR, Crump DH. Effect of the fungus Pochonia chlamydosporia and fosthiazate on the multiplication rate of potato cyst nematodes (Globodera pallida and G. rostochiensis) in potato crops grown under UK field conditions. Biol Control. 2008;46:194–201.
Vitting-Seerup K, Sandelin A. The landscape of isoform switches in human cancers. Mol Cancer Res. 2017;15:1206–20.
Vitting-Seerup K, Sandelin A. IsoformSwitchAnalyzeR: analysis of changes in genome-wide patterns of alternative splicing and its functional consequences. Bioinformatics. 2019;35:4469–71.
Wang S-L, Peng J-H, Liang T-W, Liu K-C. Purification and characterization of a chitosanase from Serratia marcescens TKU011. Carbohydr Res. 2008;343:1316–23.
Yamada T, Hiramatsu S, Songsri P, Fujie M. Alternative expression of a Chitosanase gene produces two different proteins in cells infected with Chlorella virus CVK2. Virology. 1997;230:361–8.
Yang J-I, Loffredo A, Borneman J, Becker JO. Biocontrol efficacy among strains of Pochonia chlamydosporia obtained from a root-knot nematode suppressive soil. J Nematol. 2012;44:67–71.
Zorin EA, Afonin AM, Kulaeva OA, Gribchenko ES, Shtark OY, Zhukov VA. Transcriptome analysis of alternative splicing events induced by Arbuscular Mycorrhizal Fungi (Rhizophagus irregularis) in pea (Pisum sativum L.) roots. Plants. 2020;9:1700.
We would like to thank the Plant Pathology Laboratory members of the University of Alicante for their help and support. We acknowledge the use of the Isca high-performance computing facility at the University of Exeter.
This research was funded by European Project H2020 MUSA no. 727624.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Sambles, C., Suarez-Fernandez, M., Lopez-Moya, F. et al. Chitosan induces differential transcript usage of chitosanase 3 encoding gene (csn3) in the biocontrol fungus Pochonia chlamydosporia 123. BMC Genomics 23, 101 (2022). https://doi.org/10.1186/s12864-021-08232-7