Suppressing a plant-parasitic nematode with fungivorous behavior by fungal transformation of a Bt cry gene
Microbial Cell Factories volume 17, Article number: 116 (2018)
Pine wilt disease, caused by the pinewood nematode Bursaphelenchus xylophilus (PWN), is an important destructive disease of pine forests worldwide. In addition to behaving as a plant-parasitic nematode that feeds on epithelial cells of pines, this pest relies on fungal associates for completing its life cycle inside pine trees. Manipulating microbial symbionts to block pest transmission has exhibited an exciting prospect in recent years; however, transforming the fungal mutualists to toxin delivery agents for suppressing PWN growth has received little attention.
In the present study, a nematicidal gene cry5Ba3, originally from a soil Bacillus thuringiensis (Bt) strain, was codon-preferred as cry5Ba3Φ and integrated into the genome of a fungus eaten by PWN, Botrytis cinerea, using Agrobacterium tumefaciens-mediated transformation. Supplementing wild-type B. cinerea extract with that from the cry5Ba3Φ transformant significantly suppressed PWN growth; moreover, the nematodes lost fitness significantly when feeding on the mycelia of the cry5Ba3Φ transformant. N-terminal deletion of Cry5Ba3Φ protein weakened the nematicidal activity more dramatically than did the C-terminal deletion, indicating that domain I (endotoxin-N) plays a more important role in its nematicidal function than domain III (endotoxin-C), which is similar to certain insecticidal Cry proteins.
Transformation of Bt nematicidal cry genes in fungi can alter the fungivorous performance of B. xylophilus and reduce nematode fitness. This finding provides a new prospect of developing strategies for breaking the life cycle of this pest in pines and controlling pine wilt disease.
Plant-parasitic nematodes (PPNs) are major pathogenic factors in many cash crops, including potato, soybean, and tomato [1, 2], and woody species, such as olive tree and pines [3, 4]. Destructive PPNs are very difficult to control because most of them are endoparasites. These nematodes spend most of their lives in the plant tissues, which can protect them from routine control strategies. Traditional control using highly toxic synthetic nematicides caused severe environmental problems and induced the production of chemical-resistant PPN strains. Biological agents for the control of PPNs have received greater attention in recent years, because they appear to be better solutions for crop protection against these devastating parasites. Bioactive compounds from plants and microorganisms [5, 6], novel nanoparticle delivery systems for biopesticides , nematophagous fungi and bacteria [8, 9], as well as advances in in planta transgenic or RNA interference technology [10,11,12] have significantly extended the outlook for controlling PPNs. However, relatively little attention has been paid to the microbial mutualists of nematodes that could be explored as promising targets for achieving efficient control.
The plant pathogenic pinewood nematode Bursaphelenchus xylophilus (PWN) is in ecological balance with native pine species in North America , but it has become an invasive alien species in Japan, South Korea, and China, and spread into Portugal and Spain [4, 14]. This nematode is the major causal agent of pine wilt disease, which has devastated more than one million hectares of pine forests in China . In recent years, compounds originally produced by Streptomyces species (such as avermectin, emamectin, milbemycin, and their derivatives) [16,17,18], nematicidal constituents screened from plants [19, 20], and endoparasitic fungi of nematode , were found applicable to be biopesticides for controlling PWN.
Crystal (Cry) proteins, produced by the soil bacterium Bacillus thuringiensis (Bt), are effective to control insects that destroy crops, as are the Cry proteins expressed in transgenic plants . Moreover, increasing evidence has shown that Bt Cry proteins kill a wide range of nematodes [23, 24] and their nematicidal activities can be effectively delivered to crops for controlling PPNs through transgenic modification [25,26,27]. This suggested that using Bt cry gene-modified pine trees may be a sustainable and effective strategy for the conifer-parasitic nematode, PWN. Nevertheless, genetic engineering of tree genomes is highly challenging owing to the large costs and long-term evaluation of transgenic efficacy.
Although PWN mainly feeds on xylem parenchyma cells of pines during initial infection, the nematode eats blue stain fungi, which flourish later after the pine host is killed . Spores of blue stain fungi could cling to the body surface of adult Monochamus alternatus and be transmitted to the twigs of healthy pine trees for the next cycle of PWN infection . Ophiostomatoid fungi such as Ophiostoma and Sporothrix have been reported to be associated with PWN and the insect vector in different geographic regions [29,30,31,32]. PWN seems to benefit from the proliferation of blue stain fungi around the insect pupal chambers, because nematode reproduction increases when feeding on the blue stain fungi. Moreover, the abundance of mutualistic fungi correlates with the severity of pinewood disease . The nutritional symbiotic partnership between PWN and its fungal associates implies that expressing toxic Bt protein from the fungus eaten by PWN might be a favorable alternative method for breaking the multispecies interactions among insect vectors, fungi, and the nematodes.
From the soil of the Tianmu Mountain, Bacillus thuringiensis zjfc85 was isolated and found to produce an approximately 130 kDa crystal protein Cry5Ba3, which acts as a strong nematicide against PWN, causing abnormal morphology within 48 h . However, the characteristics of Cry5Ba3 and its potential applicability in biocontrol of PWN has not been well understood. In the present study, a codon-preferred cry5Ba3Φ gene was transformed into the filamentous fungus Botrytis cinerea, which is a good diet fungus for the laboratory population of PWN, via Agrobacterium tumefaciens-mediated transformation (ATMT). We demonstrate that fungal transformation with Bt expression provides significant nematicidal activity against PWN, suggesting a prospective strategy for delivering toxins by fungus to sites where the nematode forages. This method could be used to test the efficacies of the Cry protein family against parasitic nematodes with facultative fungivorous behavior.
PWN, Botrytis cinerea, and the cry5Ba3 gene from Bacillus thuringiensis zjfc85
The PWN Bursaphelenchus xylophilus was first isolated from dead Masson pine trees in Zhejiang Province by Baermann funnels in 2010, and then cultured in the laboratory with the fungus B. cinerea (CGMCC NO.: 3.18906) on potato dextrose agar (PDA; BD Difco, Detroit, MI, USA).
Bacillus thuringiensis zjfc85 was one of 467 isolates in soils that were collected from Tianmu Mountain, Zhejiang Province . Isolate zjfc85 exhibited significantly higher nematicidal activity than did the other collected strains, because of the single crystal protein Cry5Ba3 it harbors, which is approximately 130 kDa in mass (Additional file 1: Figure S1) .
Construction of the plasmid pTFCM-cry5Ba3Φ
According to the codon usage bias described previously , a new coding gene was designed to enhance Cry protein expression in filamentous fungi and named cry5Ba3Φ (GenBank Accession Number MG737676, Additional file 1: Figure S2), which encodes a protein consisting of 698 amino acids, with molecular mass approximately 78.5 kDa. The designed cry5Ba3Φ was then synthesized, combined with trpC promoter/terminator and sticky ends XhoI/SpeI, and linked with a pUC57 plasmid to obtain pUC57-cry5Ba3Φ (Genscript Co. Ltd., Nanjing, Jiangsu Province, China) (Fig. 1a). The plasmid pTFCM containing the T-DNA border repeat sequence and the hph gene, with Aspergillus nidulans PtrpC/TtrpC (Fig. 1a), was maintained in Escherichia coli DH5α and kept in National Joint Engineering Laboratory of Biopesticide Preparation, Zhejiang A&F University.
Plasmids (pUC57-cry5Ba3Φ and pTFCM) were extracted using Plasmid Mini Prepare (Axygen Biosciences, Union City, CA, USA) according to manufacturing instructions. These plasmids were then digested at the XhoI and SpeI sites, followed by dephosphorylation of the pTFCM plasmid and purifying the digested DNA segment containing cry5Ba3Φ by AxyPrep DNA gel recycle kit (Axygen Biosciences). The cry5Ba3Φ gene with PtrpC and TtrpC was ligated to the plasmid pTFCM by the action of T4 ligase at 16 °C overnight to obtain the plasmid pTFCM-cry5Ba3Φ. The ligated product pTFCM-cry5Ba3Φ was then determined by XhoI/SpeI and another set of restriction enzyme XhoI/SacI.
Agrobacterium tumefaciens-mediated transformation
Agrobacterium tumefaciens competent cells were prepared using 10% glycerol and the plasmid pTFCM-cry5Ba3Φ was transformed into A. tumefaciens AGL-1 strain according to an electroporation method . Cells of A. tumefaciens AGL-1 strain (100 μl) carrying the plasmid pTFCM-cry5Ba3Φ were grown in 5 ml Luria–Bertani (LB) broth supplemented with 50 μg/ml kanamycin and streptomycin at 28 °C and 150 rpm for 12 h. After transferring 200 µl bacterial cells into 10 ml of induction medium (IM)  containing 200 µM acetosyringone (AS), they were grown for 6 h at 28 °C and 180 rpm. A. tumefaciens cells were finally diluted to achieve an optical density at 600 nm (OD600) of 0.15–0.3. The bacterial cells were mixed with an equal volume of B. cinerea (1 × 106 conidia/ml); thereafter, 200 µl of the mixture was spread onto sterilized cellulose membranes (cellulose nitrate) with a pore size of 0.45 µm, overlaid on co-cultivation medium (IM + AS, 10 mmol/l glucose). After co-cultivation at 28 °C for 2 days, the membranes were transferred to PDA medium amended with 100 μg/ml hygromycin B and 50 μg/ml cefotaxime to select fungal transformants and to kill A. tumefaciens cells. After incubation at 28 °C for 3–5 days, individual colonies were cultured on potato dextrose broth (PDB) containing hygromycin B (100 μg/ml) and cefotaxime (50 μg/ml) at 28 °C for another 5 days.
Genomic DNA isolation, polymerase chain reaction confirmation, and Southern blotting
Transformed and wild-type B. cinerea strains (100 μl of 1 × 106 conidia/ml) were grown in 50 ml of PDB at 25 °C and 180 rpm for 7 days. The mycelia were filtered with autoclaved gauze, washed with sterile ddH2O three times, and then ground to powder in liquid nitrogen. Genomic DNA of B. cinerea was extracted using the Cetyltrimethylammonium Bromide (CTAB) method. Polymerase chain reaction (PCR) was used to confirm the presence of cry5Ba3Φ gene by amplifying a 1.2 kb region comprising portions of hph and cry5Ba3Φ genes. Primers cry-F (5ʹ-ACTACCTCAGACCACCACA-3ʹ) and cry-R (5ʹ-TCTCAAGCCTACAGGACAC-3ʹ) were internal in cry5Ba3Φ and hph, respectively. Each PCR consisted of 1 µg of genomic DNA, 2.5 U Taq DNA polymerase, 10 µl 10 × polymerase buffer, 1.5 µM MgCl2, 200 µM dNTP, and 0.5 µM each primer, adding ddH2O to adjust the volume of the reaction mixture to 50 µl. The cycling conditions were as follows: an initial denaturation of 5 min at 94 °C, followed by 35 cycles of 45 s for denaturation at 92 °C, 1 min for annealing at 58 °C, and 1 min for polymerization at 72 °C, with a final extension at 72 °C for 10 min.
Southern blotting was done to analyze the frequency of cry5Ba3Φ gene insertion in B. cinerea. Genomic DNA of the transformant and wild-type strains was extracted by DNeasy Mini Kit (Qiagen, Valencia, CA, USA). One hundred microgram of DNA from the strains and plasmid pTFCM-cry5Ba3Φ were digested with HindIII for 24 h, followed by electrophoresis on a 0.8% agarose gel at 1 V/cm for over 10 h in 0.5% Tris–acetate-ethylenediaminetetraacetic acid (EDTA) buffer. DNA fragments were transferred to a Nylon membrane using 10× SSC (saline sodium citrate; 1.5 M NaCl plus 0.15 M sodium citrate). The probe was prepared by PCR amplification using the primers of cry-Sf (5ʹ-TGCTGAAGCTGCTGTTC-3ʹ) and cry-Sr (5ʹ-CCTTGATGGTAGTTATGGGT-3ʹ) with the DIG-High prime DNA Labeling Kit I (Roche, Mannheim, Germany), according to the manufacturer’s instructions. Hybridizations were performed at 42 °C for 12–16 h and detected using the DIG-High Detection Kit I (Roche). After hybridization, the blots were washed twice in 2× SSC plus 0.1% sodium dodecyl sulfate (SDS) for 5 min at 25 °C. Thereafter, they were washed twice again in 0.5× SSC plus 0.1% SDS for 15 min at 42 °C.
Real-time PCR analysis
Mycelia of the B. cinerea transformant and wild-type strains were collected, and the total RNA was extracted using AxyPrep™ multisource total RNA miniprep Kit (Axygen), followed by reverse transcription of mRNA with PrimeScript™ RT reagent Kit (TaKaRa, Tokyo, Japan). The PCR primers used were cry-rtF (5ʹ-CTCCCACTCACCCAACTC-3ʹ) and cry-rtR (5ʹ-TCACCCTTGGAAGCGTAT-3ʹ) and quantification of cry5Ba3Φ expression was performed using an Applied Biosystems 7300 Real-time PCR system (ABI, Foster City, CA, USA) with SYBR® Premix Ex Taq™ II (TaKaRa). The housekeeping gene β-actin (GenBank Accession No. AJ000335) was chosen as control using the primers Bcactin-F (5ʹ-AAGTGTGATGTTGATGTCC-3ʹ) and Bcactin-R (5ʹ-CTGTTGGAAAGTAGACAAAG-3ʹ).
Mitotic stability of Botrytis cinerea transformant
To determine the mitotic stability of the B. cinerea transformants with cry5Ba3Φ, they were cultured on PDA without hygromycin B for 5 days. After nine successive transfers, the colonies were tested for growth on PDA amended with 100 µg/ml hygromycin B. PCR amplification was further applied to confirm cry5Ba3Φ insertion using the primers cry-F and cry-R following the procedures mentioned above.
Nematicidal activity of Botrytis cinerea transformant
Botrytis cinerea wild-type and transformant strains were cultured on PDB for 5 days at 25 °C and 180 rpm. Culture extracts from both the strains were filtered using 0.22 µm sterile filters. The solutions prepared for this test comprised 1 ml of extracts with different wild-type:transformant ratios: 1:0, 0.9:0.1, 0.5:0.5, 0.1:0.9, and 0:1. PWNs with mixed juvenile stages were washed out with sterile ddH2O from B. cinerea wild-type plates that have been fed by nematodes for 5 days. Thereafter, 50 µl of the suspensions (containing 600 nematodes on an average) were added into the extracts. Following 24 and 48 h of contact with extracts, live or dead nematodes were recorded. Five replicates were used in all the treatments.
In another experiment, 50 µl of the nematode suspensions were added to PDA freshly grown with wild-type and cry5Ba3Φ-transgenic B. cinerea strains, respectively. After 5 days, the nematodes were flushed with 3 ml sterile distilled H2O and the number of live or dead nematodes were recorded. Six replicates were used for each treatment.
Botrytis cinerea transformants with different lengths of cry5Ba3Φ
National Center for Biotechnology Information (NCBI) Conserved Domain search (CD-search) analysis predicted that Cry5Ba3Φ contains an endotoxin-N (aa 91–327) and an endotoxin-C (aa 562–695). According to the structure of homologous crystal protein Cry5Ba (GenBank no. Q45712, Additional file 1: Figure S1), endotoxin-N of Cry5Ba3Φ may contain several helical regions (aa 115–161, 175–177, 187–207, 208–211, 222–256, 262–290, 299–315, 321–326) and endotoxin-C may consist of two helical regions (aa 571–573 and 584–586) and several beta-strand regions. Based on this information, aa 115, 202, 560, and 572 were selected to break down endotoxin-N and endotoxin-C, respectively. To determine whether integrity of segments anterior to endotoxin-N affects Cry5Ba3Φ nematicidal activity, another locus, aa 74, was selected.
Six pairs of primers were prepared to amplify fragments of cry5Ba3Φ aa 74–698 (1875 bp), 115–698 (1758 bp), 202–698 (1497 bp), 1–572 (1719 bp), 1–560 (1683 bp), and 74–572 (1503 bp). The pTFCM-TRP vector backbone was amplified using pTFCM-cry5Ba3Φ as template and pTFCM-phiF/-phiR as primers (pTFCM-phiF: 5ʹ-TACCTATTCTACCCAAGCATCCAAGATATCAGTAGATGCCGACCGGGA-3ʹ; pTFCM-phiR: 5ʹ-TTGGATGCTTGGGTAGAATAGGT-3ʹ). The steps mentioned above have been described in detail in Additional file 1: Methods. Vector (300 ng) and DNA fragment (1.20 µg) were ligated using GIBSON Assembly Cloning Kit (New England Biolabs, Ipswich, MA, USA). The ligated products were transformed into E. coli XL10-Gold for multiplication followed by A. tumefaciens-mediated transformation as described above, which resulted in the six B. cinerea transformant strains with different lengths of cry5Ba3Φ genes. During the above process, primers hph-F (5ʹ-TTCGATGTAGGAGGGCGTGGAT-3ʹ) and hph-R (5ʹ-CATTGCAGATGAGCTGTATCTGG-3ʹ) were used to determine whether truncated cry5Ba3Φ genes were successfully transformed into AGL-1.
Confirmation of Botrytis cinerea with truncated cry5Ba3Φ by PCR and SDS-PAGE
Each of the six B. cinerea transformant strains were cultivated in 5 ml of PDB supplemented with the antibiotics hygromycin B and cefotaxime at 28 °C for 5 days. Genomic DNA was then extracted for PCR amplification using primers IDF (5ʹ-ACTAGTCATTGCAGATGAGCTG-3ʹ) and IDR (5ʹ-ACTAGTCATTGCAGATGAGCTGTATCTGGA-3ʹ) to certify the presence of truncated cry5Ba3Φ genes with different lengths. The PCR program has been presented in Additional file 1: Methods.
A colony of each of the successfully transformed B. cinerea strains was transferred to 50 ml of PDB containing hygromycin B and cefotaxime at 28 °C for 5 days. After centrifugation, the pellet of each strain was resuspended in 0.5% NaOH for repeated freeze–thaw cycles. Soluble proteins, including the Cry proteins, were detected by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) followed by Coomassie blue staining.
Nematicidal activity of Botrytis cinerea transformants with truncated cry5Ba3Φ
Botrytis cinerea transformed with pTFCM and the six cry5Ba3Φ-mutant strains were grown on PDB at 25 °C for 3 days to obtain fungal extracts. Nematodes cultivated on wild-type B. cinerea plates (containing approximately 600 nematodes in 50 µl ddH2O) were added into each of the extracts. Following 24 and 48 h of contact with the extracts, live or dead nematodes were recorded. Five replicates were used for all the strains studied.
Wild-type and cry5Ba3Φ-transgenic B. cinerea strains were cultured to collect fungal extracts. The extract of the cry5Ba3Φ-transgenic B. cinerea strain was supplemented either with 100 µl of 20 mg/ml elastase (Sangon Biotech, Shanghai, China) or with 100 µl sterile ddH2O. Forty microliter of PWNs (containing approximately 600 nematodes) was added into each of the extracts (1 ml) mentioned above. Following 48 h of contact with the extracts, live or dead nematodes were recorded. Five replicates were used in all the treatments.
Differences in cry5Ba3Φ gene expression and effects on nematode fitness between transformant and wild-type strains were compared by independent samples t test. Differences in number of live nematodes and ratio of live:dead nematodes among treatments were performed by one-way analysis of variance (ANOVA). When data were not normally distributed or had no variance homogeneity, the data were rank transformed, log transformed, or were applied for non-parametric Kruskal–Wallis test. The Student–Newman–Keuls (SNK) method or Mann–Whitney U test (with adjusted α-level) was used for pairwise comparisons.
Transformation and expression of cry5Ba3Φ in Botrytis cinerea, the diet fungus of PWN in laboratory
The plasmid pUC57-cry5Ba3Φ, prepared by whole-sequence synthesis, was validated by digestion with the restriction enzymes XhoI and SpeI, resulting in an approximately 2600 bp DNA fragment (PtrpC, cry5Ba3Φ, and TtrpC), which was absent from the gel in case of the plasmid pUC57 (Fig. 1b). The ligated product pTFCM-cry5Ba3Φ, combining plasmid pUC57-cry5Ba3Φ with plasmid pTFCM, was confirmed by the restriction enzymes Xho I and Sac I, which cut off an approximately 5600 bp DNA fragment containing expression elements of both cry5Ba3Φ and the hygromycin B resistant (hph) gene (Fig. 1c). PCR amplification using primers cry-F and cry-R showed that plasmid pTFCM-cry5Ba3Φ was successfully carried by A. tumefaciens AGL-1 (Fig. 1d).
After co-cultivation of AGL-1 (carrying pTFCM-cry5Ba3Φ) and B. cinerea conidia on IM with acetosyringone, five B. cinerea transformants were picked out from hygromycin B-containing PDA medium. To make sure that the cry5Ba3Φ-transgenic fungus received the correct gene for transcription, one of the five fungal transformants was randomly selected for PCR amplification, Southern blotting, and qRT-PCR analyses. Using primers cry-F and cry-R, internal regions of cry5Ba3Φ and hph were amplified. As shown in Fig. 1e, the expected 1.2 kb fragments were observed at the correct size for cry5Ba3Φ-containing B. cinerea and plasmid pTFCM-cry5Ba3Φ, but not for the wild-type strain of B. cinerea. The PCR products were sequenced for confirming the presence of cry5Ba3Φ and hph genes further. Southern hybridization using the cry5Ba3Φ gene fragment as a probe showed that a single copy of T-DNA was integrated into the genome of the B. cinerea transformant (Fig. 1f). The qRT-PCR further determined that cry5Ba3Φ was expressed only in B. cinerea transformant, with an expression quantity approximately six times higher than that of the housekeeping gene actin (Fig. 1g).
Assessment of the genetic stability of cry5Ba3Φ-containing B. cinerea showed that the transformant maintained hygromycin B resistance after being cultured on hygromycin B-free PDA for nine successive generations (5 days for one generation; Additional file 1: Figure S3). The ninth generation was confirmed to contain cry5Ba3Φ gene using PCR detection (data not shown).
Cry5Ba3Φ confers nematicidal activity to Botrytis cinerea
The cry5Ba3Φ-transgenic B. cinerea secreted soluble proteins, including the Cry5Ba3Φ protein, which was detected by SDS-PAGE (Fig. 1h). As compared with fungal extracts of the wild-type B. cinerea, those supplemented with over 10% extracts of cry5Ba3Φ-containing B. cinerea caused significantly lower number of live PWNs both after 24 h (Fig. 2a) and 48 h (Fig. 2b). After 24 h, lower live:dead ratio of the nematode was found along with higher proportion of cry5Ba3Φ-transgenic B. cinerea extract (Fig. 2c). After 48 h, no significant difference in live:dead ratio was found among solutions containing different proportions of transformant extract, but the live:dead ratios of them were all significantly lower than that of the wild-type extract (Fig. 2d). After 5 days, the number of live nematodes feeding on cry5Ba3Φ-transgenic B. cinerea mycelia was significantly lower than that of live nematodes feeding on the wild-type strain (Fig. 2e). Compared to PWNs on the wild-type B. cinerea strain, dead:live ratio of the nematodes was much higher on the transgenic strain (Fig. 2f).
Differential effects of truncated Cry5Ba3Φ proteins on nematicidal activity
Amplified DNA fragments of cry5Ba3Φ mutants (Fig. 3a) and pTFCM-TRP vector backbone were assembled and then transformed into E. coli XL10-Gold. Putative recombinant plasmids, which showed slower mobility rates than the control on gel, were confirmed by PCR amplification (Fig. 3b, c) and sequencing. The pTFCM plasmids with different lengths of cry5Ba3Φ were then extracted for transformation into A. tumefaciens AGL-1 (Fig. 3d). Because of the co-occurrence of truncated cry5Ba3Φ and hph genes, PCR amplification using primers hph-F and hph-R indicated that cry5Ba3Φ-truncating genes on the plasmid pTFCM were successfully carried by A. tumefaciens AGL-1 (Fig. 3e). Transformant strains of B. cinerea containing truncated cry5Ba3Φ genes were constructed successfully (Additional file 1: Figure S4), because the cry5Ba3Φ gene fragments were amplified by the primers IDF and IDR, which targeted the cry5Ba3Φ expression cassette (Fig. 3f). Furthermore, expressed proteins, including truncated Cry5Ba3Φ, were detected from corresponding mutants of B. cinerea by SDS-PAGE (Fig. 3g).
As compared with the treatment using the B. cinerea transformed with pTFCM, treatments with all the cry5Ba3Φ-truncating mutants resulted in a significantly lower number of live PWNs both after 24 h (Fig. 4a) and 48 h (Fig. 4b). Significant differences were found in dead:live ratio of PWNs among the B. cinerea mutant strains. After 24 h, N-terminal deletion of the first 114 amino acids (mutant 2) significantly weakened the toxicity against nematodes, which caused lower dead:live ratio of PWNs than did the C-terminal truncation (mutant 4 and mutant 5; Fig. 4c). After 48 h, N-terminal truncated strains (mutants 1, 2, and 3) showed lower dead:live ratio of PWNs than did the C-terminal truncated strains (mutants 4, 5) (Fig. 4d). Supplementing the extract of the cry5Ba3Φ-transgenic strain of B. cinerea with elastase significantly alleviated Cry5Ba3Φ toxicity against PWN. After 48 h, as compared with the PWNs in the cry5Ba3Φ-transgenic strain extract without elastase, the nematodes in the extract containing elastase had higher number of individuals in the population (Fig. 5a) and lower dead:live ratio (Fig. 5b).
Bt crystal proteins are toxic to a variety of nematodes, particularly those from Cry5 subfamily [36,37,38]. Multiple free-living nematode species exhibited susceptibilities to Cry5B, Cry14A, and Cry21A . Cry5A, Cry5B, and Cry13 proteins were identified with significant inhibition against the free-living larval stages of nematode parasites of livestock . Cry5B was isolated from supernatant of B. thuringiensis AB88 culture . This protein was further found to be highly active in vivo against the human nematode parasite Ancylostoma ceylanicum  and was shown to have excellent potential as a control agent of the root-knot nematode Meloidogyne incognita , suggesting its wide-spectrum virulence on various parasitic nematodes. Phylogenetic analysis indicated high similarity of amino acid sequences between Cry5Ba3 and Cry5B (Additional file 1: Figure S1), supporting the toxicity of Cry5Ba3 or Cry5Ba3Φ protein against the PWNs (Fig. 2). Other toxins, belonging to the Cry5 subfamily, such as Cry5Ca1, Cry5Da1, and Cry14 were recently found to act on PPNs, such as M. incognita and M. javanica [42, 43]. The nematicidal cry5 subfamily genes, however, appeared to occur in low frequency in natural environments . This finding was consistent with our previous study, which showed that strikingly low frequency of cry genes in this subfamily can be detected in the soil samples collected .
Insecticidal Cry proteins generally contain three domains, and the function of each domain has been studied. In previous studies, researchers argued that domain I of Cry proteins was involved in the perforation of the intestinal tract during the insecticidal process [45, 46]. Domain II determined the insecticidal specificity of the Cry proteins [47, 48], binding to the gut cadherin receptors of target insects [47, 49,50,51]. Domain III appeared to prevent Cry proteins from being excessively degraded by proteases in the intestine of insects  and to modulate the permeability of gut epithelial cell channels . The nematicidal Cry5B has been crystallized for structure determination. It showed a familiar three-domain arrangement seen in insecticidal Cry proteins, but with a more structurally divergent domain II, which was implicated in interaction with glycolipid receptors of nematodes [53, 54].
Considering the high similarity of amino acid sequences between Cry5Ba3 and Cry5B, the putative domain II of Cry5Ba3Φ was retained without any excision, whereas domain I (in endotoxin-N) and III (in endotoxin-C) were designed to be broken down to detect weakening in nematicidal activity after the protein expression in B. cinerea. As compared with the wild-type B. cinerea strain, extracts from all the cry5Ba3Φ-mutant strains displayed stronger toxicities against PWNs (Fig. 4a, b). In contrast, differential virulence was found among the mutant strains of B. cinerea, particularly for their effect on dead:live ratio of the nematodes (Fig. 4c, d). The virulence of N-terminal mutants to nematodes was significantly impaired, while the mutant strains expressing C-terminal truncated Cry5Ba3Φ proteins had greater toxicities than other mutant types. This implied that N-terminal portion (domain I) was of greater significance for the toxicity of Cry5Ba3Φ protein than was the C-terminal portion (domain III). For the closely related protein Cry5B, cleavage by elastase, which disintegrated the N-terminal portion of this protein, yielded two residues of 112–170 and 173–698 amino acids . In the present study, supplementing the extract of the cry5Ba3Φ-transgenic strain of B. cinerea with elastase caused a significant decrease in its toxicity against PWN (Fig. 5), which further supported the crucial role of domain I of Cry5Ba3Φ in nematicidal ability. In the future, variations in the glycolipid-binding activities of truncated Cry5Ba3Φ proteins and their structures should be determined. Interesting, the loss of segments anterior to endotoxin-N also affected Cry5Ba3Φ nematicidal activity, which needs to be further demonstrated by experiments.
Mutagenesis using ATMT has been widely applied into filamentous fungi, including Aspergillus awamori, Penicillium digitatum, and Umbilicaria muehlenbergii [55,56,57,58]. In this study, ATMT was used to generate insertional mutation in the fungus, B. cinerea. The strong nematicidal cry5Ba3Φ gene was successfully carried and expressed by this diet fungus of PWNs. PCRs and Southern hybridizations demonstrated that cry5Ba3Φ gene was integrated into B. cinerea as a single copy, with expression level approximately six times higher than that of the housekeeping actin gene. These studies and our previous research might be among the first comprehensive attempts of searching for a novel Bt-based control strategy against PWN. In this process, B. thuringiensis strains were isolated from soils, a Cry protein was identified to have efficient nematicidal effects on PWN, and the cry gene was transformed into B. cinerea. The cry-transformed B. cinerea was found to be toxic to PWN, and finally, various cry mutant strains of B. cinerea were constructed to characterize the structure–function relationship of the Cry protein.
Botrytis cinerea may be applied as a prospective attractant for PWN in the field, albeit in usual as a diet fungus for laboratory PWN population. B. cinerea, Pestalotia, and Microzyme have been compared in their capabilities to attract PWN. Botrytis cinerea was most attractive to PWN among these fungi, possibly through secretion of extracellular active substances . This implied that B. cinerea is a suitable expression receptor for nematicidal cry genes, which could be developed as a “sweet toxin” biocontrol agent for PWNs. Recently, the recombination of symbiotic microbes of pests for producing toxic molecules has attracted attention, and some of these achieved efficient control. For instance, with the purpose of controlling human malaria parasite Plasmodium falciparum, symbiotic bacteria of its vector (Anopheles mosquito) were genetically engineered for secretion of anti-Plasmodium effector proteins to interfere with the development of P. falciparum in mosquitoes [60, 61]. Therefore, the platform we constructed using B. cinerea would inspire explorations of fungal mutualists, naturally associated with the PWN-vector Monochamus complex, which are genetically engineered to secrete nematicidal Cry proteins. A recent study reported that a native fungal symbiont, Sporothrix sp. 1, was dominant in the sites where tree infestation was higher, and its presence significantly improved the fitness of both PWN and the vector M. alternatus . This finding may provide an opportunity for further genetic engineering to produce recombinant Sporothrix sp. 1 carrying cry genes in the future. Symbioses affect the pest status of many groups of insects and nematodes, signifying a wide prospect for the use of genetically modified symbiotic microbes as powerful tools for combating forest and agricultural animal pests.
The data presented in this study demonstrate that nematicidal Cry5Ba3 found in a soil B. thuringiensis strain can be successfully expressed in B. cinerea, a diet fungus of PWN in the laboratory, via ATMT technology. The mortality of the nematodes exposed to extracts from cry5Ba3Φ-containing B. cinerea is comparable with that in the presence of purified Cry5Ba3, as previously reported. In addition, our results indicate that impairing Cry5Ba3Φ at different loci result in distinct levels of effects on its toxicity against PWN. These findings will not only inspire explorations into Bt-transgenic fungal mutualists of PWN, but also help in establishing a platform suitable for characterizing the structure–function relationships of various candidate Cry proteins against this nematode pest.
Agrobacterium tumefaciens-mediated transformation
potato dextrose agar
polymerase chain reaction
saline sodium citrate
sodium dodecyl sulfate
National Center for Biotechnology Information
conserved domain search
sodium dodecyl sulphate polyacrylamide gel electrophoresis
analysis of variance
Jung C, Wyss U. New approaches to control plant parasitic nematodes. Appl Microbiol Biotechnol. 1999;51:439–46.
Williamson VM, Gleason CA. Plant–nematode interactions. Curr Opin Plant Biol. 2003;6:327–33.
Ali N, Chapuis E, Tavoillot J, Mateille T. Plant-parasitic nematodes associated with olive tree (Olea europaea L.) with a focus on the Mediterranean Basin: a review. CR Biol. 2014;337:423–42.
Futai K. Pine wood nematode, Bursaphelenchus xylophilus. Annu Rev Phytopathol. 2013;51:61–83.
Dong L, Li X, Li H, Ying G, Zhong L, Zheng Y, Zuo Y. Lauric acid in crown daisy root exudate potently regulates root-knot nematode chemotaxis and disrupts Mi-flp-18 expression to block infection. J Exp Bot. 2014;65:131–41.
Liarzi O, Bucki P, Miyara SB, Ezra D. Bioactive volatiles from an endophytic Daldinia cf. concentrica isolate affect the viability of the plant parasitic nematode Meloidogyne javanica. PLoS ONE. 2016;11:e0168437.
Cao J, Guenther RH, Sit TL, Lommel SA, Opperman CH, Willoughby JA. Development of abamectin loaded plant virus nanoparticles for efficacious plant parasitic nematode control. ACS Appl Mater Interfaces. 2015;7:9546–53.
Li J, Zou C, Xu J, Ji X, Niu X, Yang J, Huang X, Zhang K-Q. Molecular mechanisms of nematode-nematophagous microbe interactions: basis for biological control of plant-parasitic nematodes. Annu Rev Phytopathol. 2015;53:67–95.
Moosavi MR, Zare R. Fungi as biological control agents of plant-parasitic nematodes. In: Mérillon JM, Ramawat KG, editors. Plant defence: biological control, vol. 12. Berlin: Springer; 2012. p. 67–107 (Progress in biological control).
Banerjee S, Banerjee A, Gill SS, Gupta OP, Dahuja A, Jain PK, Sirohi A. RNA interference: a novel source of resistance to combat plant parasitic nematodes. Front Plant Sci. 2017;8:834.
Davies LJ, Elling AA. Resistance genes against plant-parasitic nematodes: a durable control strategy? Nematology. 2015;17:249–63.
Dutta TK, Papolu PK, Banakar P, Choudhary D, Sirohi A, Rao U. Tomato transgenic plants expressing hairpin construct of a nematode protease gene conferred enhanced resistance to root-knot nematodes. Front Microbiol. 2015;6:260.
Wingfield MJ, Blanchette RA, Nicholls TH. Is the pinewood nematode an important pathogen in the United States? J Forest. 1984;82:232–5.
Mota MM, Vieira PR, editors. Pine wilt disease: a worldwide threat to forest ecosystems. 1st ed. Dordrecht: Springer; 2008.
Zhao L, Mota M, Vieira P, Butcher RA, Sun J. Interspecific communication between pinewood nematode, its insect vector, and associated microbes. Trends Parasitol. 2014;30:299–308.
Bi Z, Gong Y, Huang X, Yu H, Bai L, Hu J. Efficacy of four nematicides against the reproduction and development of pinewood nematode, Bursaphelenchus xylophilus. J Nematol. 2015;47:126–32.
Li J-S, Du M-N, Zhang H, Zhang J, Zhang S-Y, Wang H-Y, Chen A-L, Wang J-D, Xiang W-S. New milbemycin metabolites from the genetically engineered strain Streptomyces bingchenggensis BCJ60. Nat Prod Res. 2017;31:780–4.
Sousa E, Naves P, Vieira M. Prevention of pine wilt disease induced by Bursaphelenchus xylophilus and Monochamus galloprovincialis by trunk injection of emamectin benzoate. Phytoparasitica. 2013;41:143–8.
Barbosa P, Faria JMS, Mendes MD, Dias LS, Tinoco MT, Barroso JG, Pedro LG, Figueiredo AC, Mota M. Bioassays against pinewood nematode: assessment of a suitable dilution agent and screening for bioactive essential oils. Molecules. 2012;17:12312–29.
Liu G, Lai D, Liu Q, Zhou L, Liu Z. Identification of nematicidal constituents of Notopterygium incisum rhizomes against Bursaphelenchus xylophilus and Meloidogyne incognita. Molecules. 2016;21:1276–85.
Wang CY, Fang ZM, Wang Z, Zhang DL, Gu LJ, Lee MR, Liu L, Sung CK. Biological control of the pinewood nematode Bursaphelenchus xylophilus by application of the endoparasitic fungus Esteya vermicola. Biocontrol. 2011;56:91–100.
Knox OGG, Gupta VVSR, Nehl DB, Stiller WN. Constitutive expression of Cry proteins in roots and border cells of transgenic cotton. Euphytica. 2006;154:83–97.
Chen Q, Yang B, Liu X, Chen F, Ge F. Long-term cultivation of Bt rice expressing the Cry1Ab/1Ac gene reduced phytoparasitic nematode abundance but did not affect other nematode parameters in paddy fields. Sci Total Environ. 2017;607–608:463–74.
Höss S, Reiff N, Nguyen HT, Jehle JA, Hermes H, Traunspurger W. Small-scale microcosms to detect chemical induced changes in soil nematode communities—effects of crystal proteins and Bt-maize plant material. Sci Total Environ. 2014;472:662–71.
Ali MA, Azeem F, Abbas A, Joyia FA, Li H, Dababat AA. Transgenic strategies for enhancement of nematode resistance in plants. Front Plant Sci. 2017;8:750.
Li X-Q, Tan A, Voegtline M, Bekele S, Chen C-S, Aroian RV. Expression of Cry5B protein from Bacillus thuringiensis in plant roots confers resistance to root-knot nematode. Biol Control. 2008;47:97–102.
Li X-Q, Wei J-Z, Tan A, Aroian RV. Resistance to root-knot nematode in tomato roots expressing a nematicidal Bacillus thuringiensis crystal protein. Plant Biotechnol J. 2007;5:455–64.
Maehara N. Interactions of pine wood nematodes, wood-inhabiting fungi, and vector beetles. In: Zhao BG, Futai K, Sutherland JR, Takeuchi Y, editors. Pine wilt disease. Tokyo: Springer; 2008. p. 286–98.
Suh DY, Hyun MW, Kim JJ, Son SY, Kim SH. Ophiostoma ips from pinewood nematode vector, Japanese pine sawyer beetle (Monochamus alternatus), in Korea. Mycobiology. 2013;41:59–62.
Hyun MW, Kim JH, Suh DY, Lee SK, Kim SH. Fungi isolated from pine wood nematode, its vector Japanese pine sawyer, and the nematode-infected Japanese black pine wood in Korea. Mycobiology. 2007;35:159–61.
Wingfield MJ. Fungi associated with the pine wood nematode, Bursaphelenchus xylophilus, and cerambycid beetles in Wisconsin. Mycologia. 1987;79:325–8.
Zhao L, Lu M, Niu H, Fang G, Zhang S, Sun J. A native fungal symbiont facilitates the prevalence and development of an invasive pathogen-native vector symbiosis. Ecology. 2013;94:2817–26.
Wang Y, Mei L, Wu J, Zhang L. Detection and characterisation of a Bacillus thuringiensis crystal protein with nematicidal activity against the pinewood nematode Bursaphelenchus xylophilus. Biocontrol Sci Technol. 2012;22:1143–53.
den Dulk-Ras A, Hooykaas PJJ. Electroporation of Agrobacterium tumefaciens. In: Nickoloff JA, Totowa NJ, editors. Plant cell electroporation and electrofusion protocols. New York: Springer; 1995. p. 63–72.
Bundock P, den Dulk-Ras A, Beijersbergen A, Hooykaas PJ. Trans-kingdom T-DNA transfer from Agrobacterium tumefaciens to Saccharomyces cerevisiae. EMBO J. 1995;14:3206–14.
Rae R, Riebesell M, Dinkelacker I, Wang Q, Herrmann M, Weller AM, Dieterich C, Sommer RJ. Isolation of naturally associated bacteria of necromenic Pristionchus nematodes and fitness consequences. J Exp Biol. 2008;211:1927–36.
Schulenburg H, Müller S. Natural variation in the response of Caenorhabditis elegans towards Bacillus thuringiensis. Parasitology. 2004;128:433–43.
Wei J-Z, Hale K, Carta L, Platzer E, Wong C, Fang S-C, Aroian RV. Bacillus thuringiensis crystal proteins that target nematodes. Proc Natl Acad Sci USA. 2003;100:2760–5.
Kotze AC, O’Grady J, Gough JM, Pearson R, Bagnall NH, Kemp DH, Akhurst RJ. Toxicity of Bacillus thuringiensis to parasitic and free-living life-stages of nematode parasites of livestock. Int J Parasitol. 2005;35:1013–22.
Kostichka K, Warren GW, Mullins M, Mullins AD, Palekar NV, Craig JA, Koziel MG, Estruch JJ. Cloning of a cryV-type insecticidal protein gene from Bacillus thuringiensis: the cryV-encoded protein is expressed early in stationary phase. J Bacteriol. 1996;178:2141–4.
Cappello M, Bungiro RD, Harrison LM, Bischof LJ, Griffitts JS, Barrows BD, Aroian RV. A purified Bacillus thuringiensis crystal protein with therapeutic activity against the hookworm parasite Ancylostoma ceylanicum. Proc Natl Acad Sci USA. 2006;103:15154–9.
Baghaee Ravari S, Mahdikhani Moghaddam E. Efficacy of Bacillus thuringiensis Cry14 toxin against root knot nematode, Meloidogyne javanica. Plant Prot Sci. 2015;51:46–51.
Geng C, Liu Y, Li M, Tang Z, Muhammad S, Zheng J, Wan D, Peng D, Ruan L, Sun M. Dissimilar crystal proteins Cry5Ca1 and Cry5Da1 synergistically act against Meloidogyne incognita and delay Cry5Ba-based nematode resistance. Appl Environ Microbiol. 2017;83:e03505-16.
Bravo A, Sarabia S, Lopez L, Ontiveros H, Abarca C, Ortiz A, Ortiz M, Lina L, Villalobos FJ, Peña G, et al. Characterization of cry genes in a Mexican Bacillus thuringiensis strain collection. Appl Environ Microbiol. 1998;64:4965–72.
Bravo A, Gómez I, Porta H, García-Gómez BI, Rodriguez-Almazan C, Pardo L, Soberón M. Evolution of Bacillus thuringiensis Cry toxins insecticidal activity. Microb Biotechnol. 2013;6:17–26.
Lin X, Parthasarathy K, Surya W, Zhang T, Mu Y, Torres J. A conserved tetrameric interaction of cry toxin helix α3 suggests a functional role for toxin oligomerization. Biochem Biophys Acta. 2014;1838:1777–84.
Boonserm P, Davis P, Ellar DJ, Li J. Crystal structure of the mosquito-larvicidal toxin Cry4Ba and its biological implications. J Mol Biol. 2005;348:363–82.
Boonserm P, Mo M, Angsuthanasombat C, Lescar J. Structure of the functional form of the mosquito larvicidal Cry4Aa toxin from Bacillus thuringiensis at a 2.8-Angstrom resolution. J Bacteriol. 2006;188:3391–401.
Atsumi S, Inoue Y, Ishizaka T, Mizuno E, Yoshizawa Y, Kitami M, Sato R. Location of the Bombyx mori 175 kDa cadherin-like protein-binding site on Bacillus thuringiensis Cry1Aa toxin. FEBS J. 2008;275:4913–26.
Gómez I, Arenas I, Benitez I, Miranda-Ríos J, Becerril B, Grande R, Almagro JC, Bravo A, Soberón M. Specific epitopes of domains II and III of Bacillus thuringiensis Cry1Ab toxin involved in the sequential interaction with cadherin and aminopeptidase-N receptors in Manduca sexta. J Biol Chem. 2006;281:34032–9.
Xie R, Zhuang M, Ross LS, Gomez I, Oltean DI, Bravo A, Soberon M, Gill SS. Single amino acid mutations in the cadherin receptor from Heliothis virescens affect its toxin binding ability to Cry1A toxins. J Biol Chem. 2005;280:8416–25.
Li J, Carroll J, Ellar DJ. Crystal structure of insecticidal δ-endotoxin from Bacillus thuringiensis at 2.5 Å resolution. Nature. 1991;353:815–21.
Griffitts JS, Haslam SM, Yang T, Garczynski SF, Mulloy B, Morris H, Cremer PS, Dell A, Adang MJ, Aroian RV. Glycolipids as receptors for Bacillus thuringiensis crystal toxin. Science. 2005;307:922–5.
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–21.
de Groot MJ, Bundock P, Hooykaas PJ, Beijersbergen AG. Agrobacterium tumefaciens-mediated transformation of filamentous fungi. Nat Biotechnol. 1998;16:839–42.
Michielse CB, Hooykaas PJJ, van den Hondel CAMJJ, Ram AFJ. Agrobacterium-mediated transformation of the filamentous fungus Aspergillus awamori. Nat Protoc. 2008;3:1671–8.
Park S-Y, Jeong M-H, Wang H-Y, Kim JA, Yu N-H, Kim S, Cheong YH, Kang S, Lee Y-H, Hur J-S. Agrobacterium tumefaciens-mediated transformation of the lichen fungus, Umbilicaria muehlenbergii. PLoS ONE. 2013;8:e83896.
Wang J-Y, Li H-Y. Agrobacterium tumefaciens-mediated genetic transformation of the phytopathogenic fungus Penicillium digitatum. J Zhejiang Univ Sci B. 2008;9:823–8.
Long R. Studies on migration of Bursaphelenchus xylophilus and the extracts from medium cultured with Botrytis cinerea and its inducing activity to B. xylophilus. Xiamen: Xiamen University; 2007.
Wang S, Dos-Santos ALA, Huang W, Liu KC, Oshaghi MA, Wei G, Agre P, Jacobs-Lorena M. Driving mosquito refractoriness to Plasmodium falciparum with engineered symbiotic bacteria. Science. 2017;357:1399–402.
Wang S, Ghosh AK, Bongio N, Stebbings KA, Lampe DJ, Jacobs-Lorena M. Fighting malaria with engineered symbiotic bacteria from vector mosquitoes. Proc Natl Acad Sci USA. 2012;109:12734–9.
LZ, YW, and CC designed the study; CC and JQ contributed equally to the manuscript; CC, JQ, CW, and ML performed the experimental works; CC and JQ analyzed the data; CC and JQ wrote the paper. All authors read and approved the final manuscript.
This study was funded by the National Natural Science Foundation of China (Grant Numbers 31702018 and 31550004) and Huzhou Municipal Natural Science Foundation (Grant Number 2017YZ07). We thank doctoral student Yang LI (Nanjing Forestry University) and technician Wangying WANG (Huzhou University) for their assistances in this research.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets supporting the conclusions of this study are available upon request from the corresponding author.
Consent for publication
Ethics approval and consent to participate
All applicable international, national, and/or institutional guidelines for the care and use of animals were followed. This article does not contain any studies with human participants performed by any of the authors.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This file includes: Methods. Figure S1. Phylogenetic analysis of Cry5Ba3 with other homologous cry5 subfamily proteins. Figure S2. Codon modification of cry5Ba3 to cry5Ba3Φ. Figure S3. Fungal colony morphologies. Figure S4. Fungal colony morphologies of Botrytis cinerea transformants with different lengths of cry5Ba3Φ.
About this article
Cite this article
Cheng, C., Qin, J., Wu, C. et al. Suppressing a plant-parasitic nematode with fungivorous behavior by fungal transformation of a Bt cry gene. Microb Cell Fact 17, 116 (2018). https://doi.org/10.1186/s12934-018-0960-5