Transcriptomic insights into the effects of CytCo, a novel nematotoxic protein, on the pine wood nematode Bursaphelenchus xylophilus
BMC Genomics volume 22, Article number: 394 (2021)
The pine wood nematode Bursaphelenchus xylophilus is a destructive pest of Pinus trees worldwide and lacks effective control measures. Screening for nematotoxic proteins has been undertaken to develop new strategies for nematode control.
The results of the present study provided initial insights into the responses of B. xylophilus exposed to a nematotoxic cytolytic-like protein (CytCo) based on transcriptome profiling. A large set of differentially expressed genes (DEGs = 1265) was found to be related to nematode development, reproduction, metabolism, motion, and immune system. In response to the toxic protein, B. xylophilus upregulated DEGs encoding cuticle collagens, transporters, and cytochrome P450. In addition, many DEGs related to cell death, lipid metabolism, major sperm proteins, proteinases/peptidases, phosphatases, kinases, virulence factors, and transthyretin-like proteins were downregulated. Gene Ontology enrichment analysis showed that the CytCo treatment substantially affected DEGs involved in muscle contraction, lipid localization, and the mitogen-activated protein kinase cascade. The pathway richness of the Kyoto Encyclopedia of Genes and Genomes showed that the DEGs were concentrated in lysosomes and involved in fatty acid degradation. Weighted co-expression network analysis indicated that the hub genes affected by CytCo were associated with the nematode cuticular collagen.
These results showed that CytCo toxin interferes with gene expression to exert multiple nematotoxic effects, thereby providing insights into its potential use in pine wood nematode control.
The pine wood nematode (PWN), Bursaphelenchus xylophilus (Steiner & Buhrer) Nickle (Tylenchida: Aphenlenchoididae), a serious invasive species and the main cause of pine wilt disease (PWD), is listed as a quarantine pest in the legislation of more than 40 countries [1, 2]. It has caused severe disasters in coniferous forest ecosystems with timber losses of 106 m3 annually [2, 3]. PWNs, spread widely by beetle vectors (Monochamus spp.), invade pine trees by secreting various virulence factors such as expansin-like and venom allergen-like proteins, proliferate explosively, and ultimately kill host plants [2, 4, 5]. Several control measures for PWD are available, including fumigation with methyl bromide for the phytosanitary treatment of log exports, removing and burning dead wood in the infected areas, the trunk injection of nematocidal compounds (e.g., emamectin benzoate and abamectin), monitoring and controlling PWN vectors, and breeding resistant trees [2, 6,7,8].
Pore-forming toxins (PFTs), specifically crystal proteins (Cry) and cytolytic (Cyt) δ-endotoxins, have been widely applied in insect pest control [9,10,11]. Recently, several PFTs, such as Bacillus thuringiensis crystal proteins Cry5B, Cry6A, Cry1E, and Cry55A, were found to have nematotoxic characteristics in bioassays, indicating the potential to develop new strategies for nematode control [12,13,14,15]. For example, the Cry6Aa2 toxin has been found to suppress the hatching of the root-knot nematode Meloidogyne hapla eggs and inhibit its motility and penetration into the host plant . However, the large molecular masses and dimensions of Cry proteins decrease their control efficacy in plant-parasitic nematodes. A Cyt-like protein from the entomopathogenic fungus Conidiobolus obscurus (named CytCo) has been reported to have high nematotoxic effects on PWN, with an inhibitory effect on egg hatching and the reproductive and thrashing behaviors of PWN in bioassays . Structurally, CytCo (less than 22 kDa) has a single domain of β-strands wrapped within a layer of α-helices, which can be easily taken up by PWN, thus presenting potential for nematode control [16, 17].
The modes of action of PFTs are attributed to their cytotoxicity, causing cell lysis by toxin oligomers that assemble pores in cell membranes or elicit lethal reactions in cells through signal transduction, metabolism, and immune responses [10, 18, 19]. Specifically, in Caenorhabditis elegans, Cry6Aa was found to trigger cell death after binding to the receptor of a CUB-like-domain containing protein RBT-1, and Cry5Ba was found to use an invertebrate-specific glycolipid as its receptor for triggering cell lysis [20,21,22]. The different structures of these two Cry toxins might contribute to their different modes of action . The structure of a single α/β domain of most Cyt-like proteins is unlike that of the three-domain Cry toxins [17, 23]. This unique molecular architecture implies that CytCo utilizes a different mode of action on nematode pests. In the present study, we aimed to demonstrate PWN responses to CytCo by transcriptomic profiling to understand its potential mechanism against PWN.
General features of B. xylophilus transcriptome after treatment with CytCo
Following quality assessment and data filtering, the resultant transcriptome contained 358,738,780 clean reads (Table S1). In total, 18,034 genes were predicted by mapping to the PWN reference genome, and 1,265 DEGs (fold change ≥2, 379 upregulated and 886 downregulated) were filtered out from RNA-seq libraries between the CytCo and PBS treatments at 24 h. After consolidating similar sequences, 726 (57.3 %) DEGs were annotated in Wormbase and 541 (42.7 %) DEGs were annotated in Swiss-Prot. In total, 196 (15.5 %) DEGs were associated to GO terms and 163 (12.9 %) DEGs were annotated to 91 KEGG pathways.
The 196 GO-annotated DEGs were divided into 38 classes (level 2 subcategories) in the three ontologies of molecular function (18 classes), cellular component (seven classes), and biological process (13 classes). The largest class of DEGs was single-organism process (50 DEGs upregulated and 63 downregulated), followed by the classes of developmental process (36 upregulated and 52 downregulated) and cellular process (24 upregulated and 57 downregulated) in the biological process ontology (Fig. 1). Most DEGs in the classes of transporter activity, membrane, and membrane part were upregulated, and those in the classes of response to stimulus, multi-organism process, reproduction, metabolic process, and biological regulation were downregulated. The functional enrichment analysis of GO terms further showed that CytCo treatment significantly affected PWN genes related to lipid localization (GO: 0010876), smooth muscle contraction (GO: 0006939), and mitogen-activated protein kinase (MAPK) cascade (GO: 0000165) (Fig. S1). KEGG enrichment analysis showed that the DEGs were concentrated in lysosome, fatty acid degradation, transporters, and drug metabolism by cytochrome P450 (Fig. 2 and Table S2). Eleven DEGs were found to be putatively involved in 14 signaling pathways (Table S3).
Differentially expressed genes reflect nematode responses and CytCo toxicity
Many of the upregulated genes demonstrated potential self-protection of PWNs in response to the nematode toxicity of CytCo (Fig. 3), including 19 DEGs related to nematode cuticular collagen and epidermal growth factor (Table S4), 21 DEGs related to transporters (Table S5), and six DEGs encoding cytochrome P450 (Table S6). In addition, 10 out of 13 DEGs related to programmed cell death were downregulated (Table S7).
More DEGs were found to be downregulated with exposure to CytCo, probably related to the nematotoxicity of the protein (Fig. 3). These included all 16 major sperm protein-related DEGs (Table S8), 17 lipid metabolism-related DEGs (Table S9), 12 out of 13 virulence factor-related DEGs (Table S10), 34 out of 42 proteinase/peptidase-related DEGs (Table S11), 33 out of 35 protein kinase-related DEGs (Table S12), 19 protein phosphatase-related DEGs (Table S13), and 12 out of 14 transthyretin-like-coding DEGs (Table S6).
Co-expression network analysis
To further shed light on the key genes of the PWN responses to CytCo toxin, a weight gene co-expression network analysis (WGCNA) was built based on the pairwise correlation of genes across all samples. Highly interconnected genes were grouped into the same module, and 11 modules were obtained (Fig. 4 a). Among them, the MEturquoise module contained the largest number of genes (Fig. 4b). As is shown in Fig. 4 a, MEyellow had the strongest correlation with the CytCo treatments. Most DEGs in the MEyelllow module were grouped into GO terms related to cuticle development (Fig. 4 c), implying a nematode response to the damage caused by CytCo.
Validation of RNA-seq expression data by RT-qPCR
Nine DEGs in PWNs treated with the CytCo protein for 24 h, as per transcriptomic analysis were selected for further validation: the upregulated cuticle collagen (BXY_1699200) and ATP-binding cassette transporter (BXY_0203900) (Tables S4 and S5), and the downregulated major sperm protein (BXY_0820100), cathepsin (BXY_0408100), cytochrome P450 (BXY_0076600), serine carboxypeptidase (BXY_0963400), arginine kinase (BXY_1237900), elongation of very long chain fatty acids protein (BXY_1705500) and tumor necrosis factor α-induced protein (BXY_0951400) (Tables S6-S12). The relative expression levels of these nine genes in PWNs treated with CytCo protein, PBS, or green fluorescent protein (GFP) for 12 h, 24 h, and 36 h were evaluated by RT-qPCR (Fig. 5). There were no significant differences in the expression of all genes, except the one encoding the serine carboxypeptidase, among the groups (CytCo versus PBS or GFP tratments) at 12 h. However, at 24 h the expression levels of all of the selected DEGs were consistent with those observed in the transcriptome data. Additionally, at 36 h the expression levels became lower in the context of CytCo treatment, compared to those determined 24 h after treatment. Of note, comparing the two control treatments (PBS and GFP) at 12‒36 h, no significant differences in the expression levels of genes were detected, except for the genes encoding for serine carboxypeptidase and cathepsin, suggesting a mild response of PWN to non-toxic proteins (Fig. 5).
The present study showed that the nematotoxic protein, CytCo, exerts diverse effects on PWN development, reproduction, infectivity, motility, and immune defenses based on whole-organism transcriptome profiling, in accordance with the reported bioassay results .
It should be noted that the functions of genes discussed here have been studied in C. elegans and there have been limited studies on PWN. Many DEGs encoding serine carboxypeptidases, aspartyl proteases, cysteine proteases, and zinc metallopeptidases were downregulated in this study, which may be linked to impaired digestive capabilities . Transthyretin-like proteins are involved in nematode innate immune defenses during the interaction of C. elegans with B. thuringiensis ; here, the downregulation of most of these proteins suggests a reduced immune defense during a CytCo attack. Many DEGs annotated as “kinase and phosphatase” were downregulated (Tables S9 and S10), which may have multiple effects on development, reproduction, metabolism, and immune responses. For example, serine/threonine-protein phosphatase 1 (PP1) is essential for sperm meiosis and motility in C. elegans , in agreement with a previous report that showed reduced fecundity of the PWN under CytCo treatment .
In response to protein toxin stress, the first molecular PFT defense pathways identified in nematodes were the MAPK pathways [p38 and c-Jun N-terminal kinase (JNK)-like] in C. elegans in response to Cry5B toxin . The MAPK cascades are central signaling pathways that regulate a wide variety of stimulated cellular processes, including proliferation, differentiation, apoptosis, and stress response . Here, a DEG (BXY_0768000) encoding CRE-HSP-70 was found to be upregulated, which may inhibit apoptosis through a JNK-like MAPK pathway and involve defense against CytCo in PWNs (Fig. S2). In this study, many DEGs related to programmed cell death were downregulated (Table S7), which may be attributed to this activated pathway.
The activation of the necrosis signaling pathway by Cry6Aa has been shown to play an important role in cell death in C. elegans . Necrosis is characterized by the loss of plasma membrane integrity, and the resulting cell death can contribute to inflammation . Two necrosis-related DEGs encoding TFIP8 (tumor necrosis factor α-induced protein 8-like protein) and LITAF (lipopolysaccharide-induced tumor necrosis factor-α factor-like protein) were downregulated in this study. This implies that the Cry and Cyt toxins have differences in their modes of action, and this necessitates the comparison of gene expression patterns under treatment with different nematotoxic proteins. Additionally, the upregulated DEG encoding CREB binding protein isoform X1 (a cyclic AMP-responsive element binding protein) may influence the homeostasis of lipids and proteins in PWNs via the Jak-STAT signaling pathway, which is also involved in the immune system . Downregulated GSK3 (a serine/threonine-protein kinase) may cause an increase in β-catenin and activate Wnt signaling, which is linked to metabolism and stem cell self-renewal [30, 31]. Downregulated ADT2 (an adenine nucleotide translocator) may influence PWN development and body size by modulating the TGF-β signaling activity, which organizes cuticle collagen fibrils as in C. elegans . These genes may favor the development of new molecular targets to control PWN.
In our study, the upregulated DEGs related to sodium/sulfate symporter and potassium channel proteins (Table S5) might be related to PWN response to the pore-forming effects of CytCo on the cell membrane. Moreover, the effects of CytCo as observed in the bioassays were analogous to the adverse effects of the chemical nematicide emamectin benzoate (EB), including reduced fecundity, hatching rate, and thrashing frequency [8, 16]. Substantial transcriptional responses in PWN were observed after 24 h of exposure to EB, and only marginal responses were observed after 12 h; this is similar with the findings of qPCR assays in this study . However, some of the observed DEGs were unique, and even shared DEGs had different expression patterns (Fig. S3). For example, many cuticular collagen-related DEGs were upregulated and programmed cell death-related DEGs were downregulated with the CytCo treatment, but the opposite was reported for the EB treatment . Considering the different functional genes affected, a mixture of protein toxins and chemical agents may have a synergistic effect and could be developed for nematode control. Meanwhile, PWN showed a marginal response to small non-toxic protein GFP at the gene expression level, as per the RT-qPCR result. Considering the lack of information on transcriptomic responses to other nematotoxic proteins in plant parasitic nematodes, it is important to identify genes that are unique or shared in response to different toxic proteins in PWN and elucidate their modes of action in the future.
Materials and methods
Preparation of CytCo
CytCo protein was expressed and purified according to the method described by Zhou et al. . Briefly, Escherichia coli Arctic-Express™ cells (Agilent Technologies, Santa Clara, CA, USA) with the recombinant plasmid (pCzn1-CytCo) was inoculated for heterologous expression. The CytCo-expressing cells were harvested by centrifugation and lysed by sonication in an ice-water bath. CytCo was eluted from affinity chromatography by loading the cleared bacterial lysate onto a 1-mL Ni-IDA-Sepharose Cl-6B affinity column (Novagen, Madison, WI, USA). The protein was extensively dialyzed overnight with PBS (pH 7.4), and the final protein concentration was assessed using the Bradford Protein Assay Kit (Takara Bio Inc., Shiga, Japan) and bovine serum albumin as a standard.
Preparation of nematodes
PWNs (isolate NB-6) were collected from forests with PWD outbreaks in Ningbo City, Zhejiang, China, and fed on 7-d-cultivated Botrytis cinerea Pers. by using potato dextrose agar (PDA) plates at 25 °C. Newly emerged stage larvae (L2) were collected and inoculated on B. cinerea plates in batches. After three days, the larvae developed into adults. The Baermann funnel method was used to separate the nematodes from each PDA plate, and the nematode samples (10,000 nematodes/ml) were collected by centrifugation (4000 g) for 4 min . PWN adults (2000 nematodes/sample) were collected after being treated with 20 µg/mL purified CytCo or phosphate-buffered saline (PBS, pH 7.4) or 20 µg/mL GFP (Sangon Biotech, Shanghai, China) in the dark for 12 h, 24 h, and 36 h at 25 °C, according to the nematotoxic effect of CytCo on PWN, as previously described .
Three biological replicates were used for the total RNA extraction per treatment. The TRIzol Max Bacterial RNA Isolation Kit (Thermo Fisher Scientific, New York, NY, USA) was used according to the manufacturer’s protocol. The RNA concentration and purity were measured using a NanoDrop2000 (Thermo Fisher Scientific), and the integrity was verified by 1 % agarose gel electrophoresis and on an Agilent 2100 Bioanalyzer (Agilent Technologies). The extracted RNA samples (~ 3 µg of total RNA per sample) were stored at − 80 °C and then sent to Woosen Co. (Hangzhou, China) for sequencing or used for RT-qPCR analysis.
Nematode mRNA was enriched from each total RNA sample (24 h treatment) using oligo(dT) magnetic beads. Paired-end RNA-seq libraries of different treatments were prepared following Illumina’s library construction protocol, and the libraries were then sequenced on the Illumina HiSeq 2500 platform (Illumina, San Diego, CA, USA). FASTQ files were then produced and sorted. To obtain high-quality clean reads, raw reads were removed if they contained adapters, had ≥10 % unknown nucleotides, or possessed other low-quality indicators. The leftover clean reads were mapped to a reference genome (PRJEA64437, www.wormbase.org) by using TopHat2 (version 184.108.40.206). The raw sequencing data were deposited in the China National GeneBank Database (CNGBdb, https://db.cngb.org/) with accession number CNP0001233.
The R package RSEM was used to calculate the fragments per kilobase of exon per million fragments mapped (FPKM) value . Differentially expressed genes (DEGs) between libraries were filtered using the R package DEGsEq. DEGs were identified with a fold change ≥2 and a false discovery rate (FDR) <0.001 . To provide further insight into the DEGs involved in the modes of CytCo working on PWNs, the functions of DEGs were predicted by annotation, using several databases. Specifically, the genes were annotated by BLASTx search (E-value < 10−5) against the Wormbase, Swiss-Prot (www.uniprot.org), Gene Ontology (GO, www.geneontology.org), and Kyoto Encyclopedia of Genes and Genomes (KEGG; www.genome.jp/kegg/kegg2.html) databases .
Functional enrichment analysis for GO terms and KEGG pathways was carried out using the R package clusterProfiler and an FDR value of ≤0.05, which indicated a significant difference . Gene numbers were calculated for each GO term or pathway, and significantly enriched GO terms and pathways in DEGs compared to the genome background were defined using a hypergeometric test. The calculated P-values were subjected to FDR correction, considering FDR ≤ 0.05 as the threshold. The GO terms and pathways meeting this criterion were defined as significantly enriched GO terms or pathways in the DEGs.
Co-expression network construction
To screen the distinct hub genes influenced by CytCo toxin, gene co-expression networks were established using the WGCNA (v1.69) package of R . Gene expression values were imported into WGCNA to construct co-expression modules by using the automatic network construction function block-wise modules with default settings. The expression levels of the DEGs were log-transformed using log2 (FPKM + 1). Pearson’s correlation coefficient was used to measure the co-expression relationship between each pair of genes. The WGCNA network was constructed with a soft thresholding power of β = 17, a minimum module size of 30 genes, and the TOM-Type was unsigned, and the merge cut height was 0.25. The module-trait relationship was used to differentiate the hub genes between the CytCo and emamectin benzoate (EB) treatments. Original data of the transcriptome of EB treatments (EB12 and EB24 indicate the treatment of PWNs with EB agent for 12 and 24 h; each treatment contained three samples) and the control treatments (ddH2O, EC12, and EC24; EC, emulsifier control) were downloaded from the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) under accession number GSE135014.
RT-qPCR analysis of selected DEGs
RNA samples were extracted from PWNs exposed to either 20 µg/mL purified CytCo or 20 µg/mL GFP or PBS for 12, 24, and 36 h. RT-qPCR was performed to assess the transcript levels of selected genes, namely those encoding cuticle collagen (BXY_1699200), serine carboxypeptidase (BXY_0963400), cathepsin (BXY_0408100), cytochrome P450 (BXY_007660), major sperm protein (BXY_0820100), arginine kinase (BXY_1237900), ATP-binding cassette transporter (BXY_0203900), elongation of very long chain fatty acid protein (BXY_1705500), and tumor necrosis factor α-induced protein (BXY_0951400). For RT-qPCR analysis, total RNA (1 µg) was reverse-transcribed into cDNA by using a PrimeScript™ RT Reagent Kit with gDNA Eraser (TaKaRa, Tokyo, Japan). qPCR of the cDNA samples was performed using a SYBR Green PCR kit (SYBR Premix Ex Taq™ II; TaKaRa). Paired primers were designed and are listed in Table S1. PCR was performed on a real-time PCR thermal cycler (qTOWER 2.2; Analytik Jena, Jena, Germany) under the following conditions: 95 °C for 2 min, followed by 40 cycles at 95 °C for 5 s, 55 °C for 30 s, and 72 °C for 30 s. The data were analyzed using qPCRsoft 1.1 (Analytik Jena). The transcript quantification of each gene was performed using at least three independent replicates. The relative fold-change in gene expression was normalized to that of elongation factor 1-alpha (ef-1α) (BXY_0569100).
Availability of data and materials
Availability of data and materials, include information: Illumina sequence data have been submitted to CNGBdb (https://db.cngb.org/) under the accession number CNP0001233. NCBI non-redundant protein (Nr) database: http://www.ncbi.nlm.nih.gov; Swiss-Prot protein database: http://www.expasy.ch/sprot; KEGG Database: https://www.kegg.jp/; KEGG PATHWAY.
Database: https://www.kegg.jp/kegg/pathway; Gene Ontology Database: http://geneontology.org/; The DESeq2 and ClusterProfiler package are from R (https://bioconductor.org/). All data generated or analyzed during this study are included in this published article and its supplementary information files.
Zhao L, Jiang P, Humble LM, Sun J. Within-tree distribution and attractant sampling of propagative pinewood nematode, Bursaphelenchus xylophilus: An early diagnosis approach. Forest Ecol Manag. 2009; 258: 1932–1937.
Foit J, Čermák V, Gaar V, Nový KH, Rolincová P. New insights into the life history of Monochamus galloprovincialis can enhance surveillance strategies for the pinewood nematode. J Pest Sci. 2019; 92: 1203–1215.
Futai K. Pine Wood Nematode, Bursaphelenchus xylophilus. Annu Rev Phytopathol. 2013; 51: 61–83.
Kikuchi T, Li H, Karim N, Kennedy M, Moens M, Jones J. Identification of putative expansin-like genes from the pine wood nematode, Bursaphelenchus xylophilus and evolution of the expansin gene family within the Nematoda. Nematology. 2009; 11: 355–364.
Kang JS, Koh YH, Moon YS, Lee SH. Molecular properties of a venom allergen-like protein suggest a parasitic function in the pinewood nematode Bursaphelenchus xylophilus. Int J Parasitol. 2012; 42: 63–70.
Yang ZQ, Wang XY, Zhang YN. Recent advances in biological control of important native and invasive forest pests in China. Biol Control. 2014; 68: 117–128.
Seabright K, Davila-Flores A, Myers S, Taylor A. Efficacy of methyl bromide and alternative fumigants against pinewood nematode in pine wood samples. J Plant Dis Prot. 2020; 127: 393–400.
Lu F, Guo K, Chen AL, Chen SN, Lin H, Zhou X. Transcriptomic profiling of effects of emamectin benzoate on the pine wood nematode Bursaphelenchus xylophilus. Pest Manag Sci. 2020; 76: 747–757.
Li J, Koni PA, Ellar DJ. Structure of the mosquitocidal δ-endotoxin CytB from Bacillus thuringiensis sp. kyushuensis and implications for membrane pore formation. J Mol Biol. 1996; 257: 129–152.
Bravo A, Gill SS, Soberón M. Mode of action of Bacillus thuringiensis Cry and Cyt toxins and their potential for insect control. Toxicon. 2007; 49: 423–435.
Cohen S, Albeck S, Ben-Dov E, Cahan R, Firer M, Zaritsky A, Dym O. Cyt1Aa toxin: crystal structure reveals implications for its membrane-perforating function. J Mol Biol. 2011; 413: 804–814.
Yu Z, Xiong J, Zhou Q, Luo H, Hu S, Xia L, Sun M, Li L, Yu Z. The diverse nematicidal properties and biocontrol efficacy of Bacillus thuringiensis Cry6A against the root-knot nematode Meloidogyne hapla. J Invertebr Pathol. 2015; 125: 73–80.
Zhang F, Peng D, Cheng C, Zhou W, Ju S, Wan D, Yu Z, Shi J, Deng Y, Wang F, Ye X, Hu Z, Lin J, Ruan L, Sun M. Bacillus thuringiensis Crystal Protein Cry6Aa Triggers Caenorhabditis elegans Necrosis Pathway Mediated by Aspartic Protease (ASP-1). Plos Pathog. 2016; 12: e1005389.
Cheng CH, Qin JL, Wu CF, Lei MY, Wang YJ, Zhang LQ. Suppressing a plant-parasitic nematode with fungivorous behavior by fungal transformation of a Bt cry gene. Microb Cell Fact. 2018; 17: 116.
Huang T, Lin Q, Qian X, Zheng Y, Yao J, Wu H, Li M, Jin X, Pan X, Zhang L, Guan X. Nematicidal activity of Cry1Ea11 from Bacillus thuringiensis BRC-XQ12 against the pine wood nematode (Bursaphelenchus xylophilus). Phytopathology. 2018; 108: 45–52.
Zhou X, Chen SN, Lu F, Guo K, Huang LL, Su X, Chen Y. Nematotoxicity of a Cyt-like protein toxin from Conidiobolus obscurus (Entomophthoromycotina) on the pine wood nematode Bursaphelenchus xylophilus. Pest Manag Sci. 2021; 77: 686–692.
Wang Y, Chen SN, Wang JH, Zhou X. Characterization of a cytolytic-like gene from the aphid-obligate fungal pathogen Conidiobolus obscurus. J Invertebr Pathol. 2020; 173: 107366.
López-Diaz JA, Cantón PE, Gill SS, Soberón M, Bravo A. Oligomerization is a key step in Cyt1Aa membrane insertion and toxicity but not necessary to synergize Cry11Aa toxicity in Aedes aegypti larvae. Environ Microbiol. 2013; 15: 3030–3039.
Xu C, Wang BC, Yu Z, Sun M. Structural Insights into Bacillus thuringiensis Cry, Cyt and Parasporin Toxins. Toxins. 2014; 6: 2732–2770.
Griffitts JS, Haslam SM, Yang T, Garczynski SF, Mulloy, B. Glycolipids as receptors for Bacillus thuringiensis crystal toxin. Science. 2005; 307: 922–925.
Hui F, Scheib U, Hu Y, Sommer RJ, Aroian RV, Ghosh P. Structure and glycolipid binding properties of the nematicidal protein Cry5B. Biochemistry. 2012; 51: 9911–9921.
Shi J, Peng D, Zhang F, Ruan L, Sun M. The Caenorhabditis elegans CUB-like-domain containing protein RBT-1 functions as a receptor for Bacillus thuringiensis Cry6Aa toxin. PLoS Pathog. 2020; 16: e1008501.
Soberón M, Lópezdíaz JA, Bravo A. Cyt toxins produced by Bacillus thuringiensis: a protein fold conserved in several pathogenic microorganisms. Peptides. 2013; 41: 87–93.
Treitz C, Cassidy L, Höckendorf A, Leippe M, Tholey A. Quantitative proteome analysis of Caenorhabditis elegans upon exposure to nematicidal Bacillus thuringiensis. J Proteomics. 2015; 113: 337–350.
Wu JC, Go AC, Samson M, Cintra T, Mirsoian S, Wu TF, Jow MM, Routman EJ, Chu DS. Sperm Development and Motility are Regulated by PP1 Phosphatases in Caenorhabditis elegans. Genetics. 2012; 190: 143–157.
Huffman DL, Abrami L, Sasik R, Corbeil J, van der Goot FG, Aroian RV. Mitogen-activated protein kinase pathways defend against bacterial pore-forming toxins. P Natl Acad Sci USA. 2004; 101: 10995–11000.
Tong SM, Feng MG. Insights into regulatory roles of MAPK-cascaded pathways in multiple stress responses and life cycles of insect and nematode mycopathogens. Appl Microbiol Biotechnol. 2019; 103: 577–587.
Chan KM, Luz NF, Moriwaki K. Programmed Necrosis in the Cross Talk of Cell Death and Inflammation. Annu Rev Immunol. 2015; 33: 79–106.
Ivashkiv LB. Jak-STAT signaling pathways in cells of the immune system. Rev Immunogenet. 2000; 2: 220–230.
Clevers H, Nusse R. Wnt/β-catenin signaling and disease. Cell. 2012; 149: 1192–1205.
Gómez-Orte E, Sáenz-Narciso B. Moreno S, Cabello J. Multiple functions of the noncanonical Wnt pathway. Trends Genet. 2013; 29: 545–553.
Fernando T, Flibotte S, Xiong S, Yin J, Yzeiraj E, Moerman DG, Meléndez A, Savage-Dunn C. C. elegans ADAMTS ADT-2 regulates body size by modulating TGFβ signaling and cuticle collagen organization. Dev Biol. 2011; 352: 92–103.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf. 2011; 12: 323–338.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential genes and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012; 7: 562–578.
Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000; 28: 27–30.
Yu G, Wang LG, Han Y, He QY. ClusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. Omics. 2012; 16: 284–287.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowite 2. Nat Methods. 2012; 9: 357.
Thanks for the technical support from Dr. Linlin Pan (Beijing Biomedi Technology Company) and Mr. Junhao Chen (Saint Louis University).
This research was funded by the National Natural Science Foundation of China(31870637、31200487)and jointly funded by Zhejiang Key Research Plan (2016C32016、2019C02024).
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, Y., Zhou, X., Guo, K. et al. Transcriptomic insights into the effects of CytCo, a novel nematotoxic protein, on the pine wood nematode Bursaphelenchus xylophilus. BMC Genomics 22, 394 (2021). https://doi.org/10.1186/s12864-021-07714-y
- Cyt δ-endotoxin
- Plant parasitic nematode
- Transcriptome profiling