1. Introduction
Entomopathogenic nematodes (EPNs) have come to be widely accepted as some of the most optimal biological control agents to replace chemical control methods in insect pest management [
1,
2]. Two well-known EPNs include two independently originated families of Heterorhabditidae and Steinernematidae
, which exhibit similar parasitic modes of life cycles through a convergent evolution [
3]. They are typically found in the soil and search for their hosts based on characteristic behaviors wherein they act as ambushers or cruisers [
4]. Ambusher nematodes implement a strategy where they sit and wait in soil in search of hosts whereas cruiser nematodes use active searching to seek out potential target insects using the host-emitted volatile organic compounds [
5]. Once they are damaged by herbivores, plant roots release terpenes such as α-pinene, linalool, d-limonene, and pregeijerene, which are detected by these EPNs [
6].
Once the infective juveniles (IJs) enter the host insect hemocoel through natural openings such as the mouth, anus, or spiracles, or through direct cuticle penetration, they release their symbiotic bacteria [
7]. The symbiosis between the nematodes and their bacteria is species-specific, where
Steinernema and
Heterorhabditis nematodes respectively have
Photorhabdus and
Xenorhabdus bacteria in their gut [
8]. The bacteria are located in a specific receptacle in the intestine near the esophagus of the nematodes, and they are released by intestinal contraction in response to insect hemolymph factor(s) [
9]. In the complexes of
S. carpocapsae and
X. nematophila, nematode intestine localization (nil) factors A, B, and C are identified as molecular components that explain the specificity between the nematode and bacteria [
10]. NilB and NilC are an outer membrane beta barrel protein and periplasmic lipoprotein, respectively [
11,
12], suggesting their roles in the molecular interactions with host gut epithelium.
The bacteria released in the hemocoel suppress the insect immune responses to protect the bacteria from the fatal cellular and humoral responses, with the intention of killing the insect by septicemia [
13]. To suppress the immune responses, these bacteria release several secondary metabolites to inhibit phospholipase A
2 (PLA
2) [
14]. PLA
2 catalyzes phospholipids to release free fatty acids, with are typically arachidonic acid or its biosynthetic precursor in insects [
15]. AA is then oxygenated to produce prostaglandins, leukotrienes, or eicosatrienoic acid, which are collectively called eicosanoids [
16]. Eicosanoids mediate various immune responses in insects [
17]. In this way, the inhibitory activity by the bacterial metabolites leads to host immunosuppression. In the insect cadaver, the IJs are then developed to the reproductive stage to produce progeny. After consuming all nutrients in the host insect, the nematodes develop into IJs through re-association with their specific bacteria and then emerge from the insects to look for new hosts.
Nematodes can be effective in controlling coleopteran larvae inhabiting the rhizosphere. For example, one study showed that, in peanut fields, the white grub,
Holotrichia parallela, was effectively controlled by the nematodes,
S. longicaudum and
H. bacteriophora at 10,000 and 5,000 IJs per plant, respectively [
18]. Insect pests inhabiting the aquatic environment are also targets of control by nematodes. For example, mosquitoes and tephritid flies have previously been controlled by
Steinernema and
Heterorhabditis nematodes [
19,
20].
In addition to
Steinernema and
Heterorhabditis, other rhabditid nematodes have been evaluated as efficient EPNs, with notable examples including
Oscheius tipulae,
O. carolinensis [
21],
O. chongmingensis [
22], and
O. gingeri [
23]. In particular,
O. tipulae has been proposed to control tephritid flies due to its high virulence in addition to its previously documented pathogenicity to a lepidopteran insect [
24]. This relatively wide host range has led to the nematode,
O. tipulae, being proposed an effective biological control agent against insect pests. The present study identified a strain of
O. tipulae in Korea and evaluated its pathogenic spectrum and its potential to control insect pests.
2. Materials and Methods
2.1. Test Insects
Altogether, this study used three different lepidopterans (
Spodoptera exigua,
Maruca vitrata, and
Plutella xylostella), three different dipterans (
Delia platura,
Drosophila melanogaster, and
Aedes albopictus), and one coleopteran (
Tenebrio molitor) species. The rearing conditions in the laboratory were a temperature of 25 ± 2°C with a relative humidity of 60 ± 5% and a photoperiod of 16:8 h (L:D). The larvae of
S. exigua, which were collected from a field in which they were infesting Welsh onion (
Allium fistulosum) in Andong, Korea, were reared on an artificial diet [
25] and underwent five instars (L1-L5). For bioassay, L4 larvae were used. The larvae of
P. xylostella collected from Andong, Korea were reared on fresh cabbage (
Brassica rapa) and underwent four instars, in which L4 larvae were used for bioassay. The larvae of
M. vitrata, which were collected from an Adzuki bean (
Vigna angularis) field in Suwon, Korea, were reared on an artificial diet [
26] and underwent five instars. For bioassay, L4 larvae were used. The larvae of
T. molitor, which were supplied by Bio Utility (Andong, Korea), were reared on a wheat bran-based diet [
27] and underwent ~20 instars [
28]. For bioassay, two groups of larvae were selected based on their body length: ~10 mm for young larvae and ~25 mm for old larvae.
D. melanogaster adults were purchased from Biological Nara (Seoul, Korea) and reared on an artificial diet (Hansol Tech, Yangsan, Korea). The larvae of
D. platura were collected from a Welsh onion field in Andong, Korea and reared on bulb onion. The larvae of
A. albopictus were donated by Professor Donghun Kim in Kyungpook National University (Sangju, Korea). First instar larvae were reared in a plastic box (12 × 8 × 5 cm) containing 1 L distilled water and provided with a fish food mini pellet (Ilsung, Daegu, Korea), while adults were maintained on 10% sugar solution using a cotton plug. Female mosquitoes were fed mouse blood as described by Choi et al. [
29]. For laboratory bioassays of the dipteran insects, the penultimate larval instar (L3) was used.
2.2. Culturing Nematodes of Two Steinernema carpocapsae Strains
Two
S. carpocapsae strains, Sc-AD and Sc-DR, were respectively collected from Andong and Dongrae in Korea [
30]. Sc-AD was multiplied using
S. exigua host while Sc-DR was cultured in
T. molitor. The infection and collection of the nematodes were conducted following the method described by Park et al. [
30]. Briefly, approximately 200 IJs in 500 μL distilled water were applied to Petri dishes (9 cm diameter, 3 cm height) that each contained an individual insect larva. Treated larvae were reared with artificial diet under the laboratory rearing conditions. At 5 days after the nematode infection, the emerging IJs were collected from the dead larvae using a White trap and then transferred to a 6 cm culture flask (Bibby Sterilin, Staffordshire, UK). IJs were kept in 10 ± 1°C in distilled water under darkness until being used.
2.3. Chemicals
Arachidonic acid (AA: 5,8,11,14-eicosatetraenoic acid) and dexamethasone (DEX: (11β,16α)-9-fluoro-11,17,21-trihydroxy-16-methylpregna-1,4-diene-3) were obtained from Sigma-Aldrich Korea (Seoul, Korea) and dissolved in dimethyl sulfoxide (DMSO). Chelex 100 Resin was purchased from Bio-Rad (Hercules, CA, USA) and used for DNA extraction. Luria-Bertani (LB) medium was purchased from MBcell (Seoul, Korea). Phosphate-buffered saline (PBS) solution was prepared with 100 mM phosphoric acid and adjusted to pH 7.4 with 1 N NaOH. A commercial version (Captin®, Kyungnong, Seoul, Korea) of fluxamethamide, an insecticide, was purchased from a local distributor.
2.4. Soil Samplings and EPN Collection
Soils were collected from three different field sites (SC1, SC2, and YS) located at Andong, Korea (
Figure S1). The sampling sites were host to different plants, including
Gingko biloba,
Quercus dentata,
Cucumis sativus, and
Juglans regia. At each site, approximately 200-400 g of each soil sample was gathered using a hand shovel and stored in a plastic bag once a week from March through September.
Live nematodes were collected from these gathered soil samples using a modified Baermann funnel technique, as described by Lee et al. [
31]. The nematode suspensions were topically applied on the larvae of
T. molitor and
S. exigua. Subsequent dead larvae were collected and placed onto a White trap that was kept at room temperature to facilitate the emergence of IJs [
32]. The IJs were multiplied using
S. exigua that were collected using the method described above.
2.5. Nematode Identification with Morphological Characteristics
For morphological characterization, 20 IJs were randomly selected from T. molitor cadavers. The nematodes were fixed using warm water (60°C) in 4% formaldehyde solution and mounted on slide glass with the addition of glycerin. A glass rod was used to prevent sample flattening while the cover glass was overlaid. Observations were made using a stereoscopic microscope (Stemi SV11, Zeiss, Jena, Germany) and a light microscope (Eclipse 80i, Nikon, Tokyo, Japan) at 50× magnification.
2.6. Nematode Identification with Molecular Characters
To begin, a sample of about 0.5 g of IJs was homogenized in 100 μL of 5% Chelex-100 (Sigma-Aldrich, Seoul, Korea) and incubated in water bath at 100°C for 10 min. After centrifugation at 12,500 ×
g for 5 min, the extracted DNA was transferred to a new tube and stored at -20°C. A rRNA gene covering two internal transcribed spaces (ITS-1 and ITS-2) was amplified using universal primers: 5ʹ-TTGATTACGTCCCTTGCCCTT-3ʹ and 5ʹ-TTTCACTCGCCGTTACTAAGG-3ʹ [
33]. The PCR mixture consisted of 5 μL of 10× PCR buffer, 15 μL of deionized water, 1 μL of dNTP, 1 μL of Taq polymerase, 1 μL of forward primer, 1 μL of reverse primer, and 1 μL of template DNA (~100 ng). PCR comprised of 35 amplification cycles, where each cycle was programmed to consist of 1 min at 94°C for denaturation, 1 min at 46°C for annealing, and 1 min at 72°C for extension. After cloning into pCR2.1 cloning vector, the PCR product was bidirectionally sequenced with M13 forward and reverse universal primers by Macrogen sequencing company (Seoul, Korea). The obtained nucleotide sequences were analyzed using the BlastN program from the National Center for Biotechnology Information (
www.ncbi.nlm.nih.gov) to match the most similar sequences of specific species. Phylogenetic tree analysis was conducted using a Neighbor- Joining method with MEGA6 [
34]. Bootstrap values on the branches were estimated with 1,000 repetitions.
2.7. Nematode Virulence Test
The insecticidal activity of the nematode isolate was determined by applying different concentrations (0, 100, 200, 400 IJs/larva) to filter paper in a 9-cm diameter Petri dish containing 10 larvae. The experimental unit (i.e., each Petri dish) was replicated three times per nematode concentration. For mosquito bioassay, L3 larvae of A. albopictus were bathed in 4 mL of distilled water and then placed in 8 mL plate well (4 cm in diameter). Four concentrations of nematodes were tested (0, 50, 100, and 200 IJs/mL). Distilled water was used for a control. For this bioassay, three lepidopteran species (S. exigua, M. vitrata, and P. xylostella) at L4 instar and D. melanogaster at L3 instar were used. Further, the larvae of a coleopteran insect (T. molitor) with approximately 25 mm body length were used in the assay. Treated larvae were provided with diets and kept in the previously described rearing conditions. Mortality was assessed at 4 days after treatment (DAT).
2.8. Nematode Life Cycle Assessment
The life cycle of the nematode isolate was examined using M. vitrata larvae as hosts. Each L4 larva was infected with 50 IJs per larva as described above. The treated insects were dissected daily for 10 days. Each time point was replicated with three different infected larvae and each nematode life stage was counted. The time to emergence, i.e., the number of days after infection to the moment of IJ release from the cadaver, was determined on the White traps through daily observation until the first infective juveniles started to come out of the cadavers.
2.9. Isolation of the Nematode-Associated Bacteria
The bacterial extraction from the nematodes was conducted following the method described by Loulou et al. [
35]. Briefly, approximately 200 IJs were washed three times with sterile PBS and incubated in 1% sodium hypochlorite (Yakuri Pure Chemicals, Kyoto, Japan) in PBS for 5 min with gentle agitation. After removing the incubation solution, the nematode sample was incubated again with 1% sodium hypochlorite. After several washes with sterile PBS, the nematodes were sonicated with a power strength of 100 for 1 min using a ultrasonicator (MS73, Bandelin Sonoplus, Heinrichstraße, Berlin, Germany). After centrifugation at 650 ×
g for 10 min at 4°C, the supernatant was used for bacterial culture or metagenomics analysis. For bacterial culture, the extracts were plated on LB agar plates for 24-48 h at 28°C. The resulting colonies were then classified into three types based on the colony shape: fuzzy, small, and big. These colonies were used for subsequent virulence testing.
2.10. Bacterial Virulence Assessment
To assess the insecticidal activities of the isolated bacteria, freshly overnight cultured bacterial cells were washed three times with sterilized PBS through centrifugation at 8,000 rpm for 2 min at 4°C, and bacterial pellets were resuspended in PBS. After counting the bacterial cells, the different bacterial concentrations (0, 101, 102, 103, 104, and 105 colony-forming unit (cfu)/larva) were prepared through a serial dilution and injected into the L5 larvae of M. vitrata. One microliter of bacterial solution was injected into the larvae hemocoel using a 10 μL micro syringe (Hamilton, Reno, NV, USA). Ten surface-sterilized L5 larvae were injected per bacterial concentration. The control larvae were injected with PBS buffer. The larvae were placed in sterile plastic boxes and kept in the dark at 25°C; insect mortality was then determined at 72 h post-injection. All assays were replicated three times.
2.11. Metagenomics Analysis of the Nematode-Associated Bacteria
Three different nematodes were used in this metagenomics analysis: Sc-DR, Sc-AD, and
O. tipulae (Ot). The bacteria were extracted while following the method described above. The DNA was amplified for the V3-V4 region of the 16S rRNA gene with the Herculase II Fusion DNA Polymerase kit using 27-F (5ʹ-TTGATTACGTCCCTTGCCCTT-3ʹ and 1492-R (5ʹ-TTTCACTCGCCGTTACTAAGG-3ʹ). The library preparation was done using the DNA Polymerase Nextera XT kit, which is suitable for working with small genomes. Indexed libraries were then joined at equimolar concentrations, and the product was sequenced with a MiSeq device, and it was set to sequence the PCR product 300 base pairs from both ends. The MiSeq benchtop sequencer output is over 100,000 reads with a length of 300 and a total of more than 30 million bases. Upon completion of this sequencing, Illumina MiSeq raw data were sorted by sample utilizing the index sequences, and paired-end FASTQ files were generated for each sample. According to Martin [
36], Cutadapt v3.2 was used to eliminate sequencing adapter and primer sequences. The DADA2 package in R v4.0.3 was used to rectify errors in the amplicon sequencing process [
37]. Briefly, paired-end reads with expected errors exceeding 2 were excluded. Post-preprocessing, an error model was constructed for each batch to eliminate sample noise. The corrected paired-end sequences were merged into a single sequence, and chimera sequences were filtered out using DADA2’s consensus method, which ultimately resulted in the formation of Amplicon Sequence Variants (ASVs). ASVs shorter than 350 bp were removed. Moreover, for microbial community comparison analysis, subsampling was utilized for normalization via QIIME v1.9 [
38,
39] based on the minimum read count across all samples. Each ASV sequence was taxonomically annotated by BLAST+ v2.9.0 against the reference database (NCBI 16S Microbial DB), where taxonomy information was assigned based on the organism with the highest similarity. Utilizing the abundance and taxonomy data of ASVs, various comparative analyses of microbial communities were conducted using QIIME, and species diversity and evenness within microbial communities were evaluated using Shannon and Simpson index [
40].
2.12. Statistical Analysis
The means and variances of the treatments were analyzed by one-way ANOVA by PROC GLM of the SAS program [
41]. All the data were plotted in terms of mean ± standard error using Sigma Plot (Systat Software, Inc., Point Richmond, CA, USA) and GraphPad Prism software (GraphPad Software, Inc., La Jolla, CA, USA). The least squared difference (LSD) test was used to compare the means with a Type I error of 0.05.
4. Discussion
Entomopathogenic nematodes have been successfully used to control insect pests. Although there is substantial nematode species diversity, which has been estimated at ~30,000 described species [
44], two families of EPNs have been the main objects of focus since the first proposal [
45]: Steinernematidae and Heterorhabditidae. This study isolated a non-model EPN species and analyzed its application in the control of insect pests.
An EPN isolate from a soil sample killed lepidopteran and coleopteran insects and was found to multiply in insect cadavers. It was identified as
Oscheius tipulae based on analyses of both morphological and molecular characters. The species was named as such because it was first isolated from a dipteran insect,
Tipula paludosa [
46]. Our phylogenetic analysis using ITS sequences showed that
O. tipulae was clustered with other
Oscheius spp. and
Heterorhabditis spp. and separated from
Steinernema spp. This supports the nematode classification [
47] that classified two separate suborders by clustering
Oscheius and
Heterorhabditis into
Rhabditina while
Steinernema was classified into
Tylenchina in the Order Rhabditida.
The infective juveniles (= 3rd instar) developed into adults in 4-5 days.
O. tipulae has previously been shown to have a fast developmental rate in in vitro culture, where it develops 3-4 days at room temperature after four juvenile stages [
48]. Some delay in the parasitic life cycle might consist of elapsed time taken during the infection stage, such as penetration to insect hemocoel and suppression of host immune defense.
The insect-killing activity of
O. tipulae was unusual because nematodes have been known to be bacterivorous and used as animals in genetic research due to their close relationship with a well-known
Caenorhabditis elegans and its easy culture [
49,
50]. The insecticidal activity of
O. tipulae was assessed in a Mediterranean fly,
Ceratitis capitata, in which the nematodes killed all the immature stages of the flies, especially with high susceptibility of pupae, even at a few IJs per pupa [
24]. This study suggested that
O. tipulae could be used as a biological control agent against
C. capitata flies. In
Oscheius, about 30 species are recorded, of which more than 13 are known to kill insects [
51].
O. tipulae appear to occur in a wide range of soil environments. A survey of EPNs using White trap baited with
Galleria mellonella collected
O. tipulae from agricultural soils in Brazil in addition to
Heterorhabditis and
Steinernema [
52]. Our current study showed that
O. tipulae was obtained from a non-agricultural area and showed a wide host range, including lepidopteran, coleopteran, and dipteran insects, for its parasitic life.
Suppression of host immunity was required for the insecticidal activity of
O. tipulae. Eicosanoids mediate various immune responses in insects against virus, bacterial, and fungal infections [
17]. PLA
2 catalyzes the biosynthetic pathway of the eicosanoids by releasing arachidonic acid (AA) as the biosynthetic substrate [
53]. Dexamethasone treatment significantly increased the nematode insecticidal activity whereas AA treatment suppressed the activity. This suggests that eicosanoids play a crucial role in defending against nematode infection. For EPNs such as
Heterorhabditis and
Steinernema, their symbiotic bacteria suppress the target insect immune responses by inhibiting PLA
2 [
14]. The enzyme activity is inhibited by the secondary metabolites released by the nematode-specific symbiotic bacteria,
Photorhabdus and
Xenorhabdus [
54]. These findings led us to examine the bacteria associated with
O. tipulae: Three different bacteria were isolated from the surface-sterilized IJs of
O. tipulae and showed insecticidal activities. These bacteria were identified to be
S. marsescens,
S. nematodiphila, and
K. pneumonia from 16S rRNA sequences. They were also predicted from the metagenomics analysis of
O. tipulae. Surprisingly, the metagenomics analysis also suggested that
X. nematophila, which is specific to
S. carpocapsae, is associated with
O. tipulae. Our previous work showed that
X. nematophila produces at least eight different secondary metabolites, inhibiting PLA
2 [
55]. These findings suggest that the insecticidal activity of
O. tipulae may be caused by its associated bacteria. However, it is too premature to conclude that
X. nematophila is associated with
O. tipulae until the bacteria can be isolated. Our current metagenomics analysis showed that two strains of the same species,
S. carpocapsae, exhibited a bacterial fauna that was presumably attributable to their insect hosts. Even though there were differences in the two strains of
S. carpocapsae,
O. tipulae had lower species diversity than
S. carpocapsae. The relatively low bacterial species of
O. tipulae occurred in both entomopathogenic and non-entomopathogenic bacteria from the metagenomics analysis.
The use of EPNs has been regarded as one of the highly promising biological control strategies due to their capacity to efficiently control different insect pests with a well-developed mass production technique [
56]. Our current study applied the nematode isolate identified as
O. tipulae to control an onion fly,
Delia platura, of which larvae damages the seedlings of Welsh onions, particularly at the rhizosphere [
57].
O. tipulae gave significant insecticidal activities to other dipteran insects, including a mosquito,
A. albopictus. This mosquito showed high resistance to the other EPN,
S. carpocapsae, which did not show more than 40% mosquitocidal activity. The resistance of
A. albopictus to
S. carpocapsae might be a cellular immune defense of the mosquito in exhibiting hemocytic encapsulation around the IJs infected to the mosquito hemocoel [
58]. This then leads to the question of why
O. tipulae can evade or suppress the mosquito immune responses. It could be attributable to be the immunosuppressive activity of the nematode-associated bacteria by releasing secondary metabolites. This needs to be explored in a subsequent study. Our current study also showed that
O. tipulae demonstrated outstanding control efficacy to Drosophila melanogaster. These results suggest that
O. tipulae exhibits a wide range of pathogenic spectrum against various insect pests, which is a useful characteristic for its application to biological control. Our current study also showed the potential for the long-term storage of
O. tipulae using 10°C in distilled water. The storability may be enhanced by formulating IJs in alginate gel to increase survival without impairing infectivity [
59].
Author Contributions
Conceptualization, Y.K., E.A., M.E.; Methodology, Y.K., E, A. M.E., J.K., G.J.; formal analysis, E.A., M.E., J.K., G.J.; Data curation and software, E.A., M.E., J.K., G.J.; Investigation, Y.K. E.A., M.E., J.K., G.J.; Funding acquisition, Y.K.; supervision, Y.K.; writing—original draft, Y.K., E.A; writing—review & editing, Y.K. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Identification of a nematode isolate killing insects using morphological characteristics. (A) Release of infective juveniles (IJs) from the dead T. molitor larva. The right panel shows an enlarged picture of the square region in the left panel as well as IJs. (B) Morphological traits of male and female specimens: TL for tail length, RL for rectum length, ovary, BDM for body diameter, BL for body length, and SL for stylet length. (C) Genus identification of the nematode isolate using the morphological characters. Each measurement used 10 nematode samples and was reported as mean ± SD. (D) Species identification of the nematode isolate using the morphological characteristics. ‘+’ and ‘-’ in the parenthesis indicate being identical and non-identical with the nematode isolate with respect to each morphological characteristic.
Figure 1.
Identification of a nematode isolate killing insects using morphological characteristics. (A) Release of infective juveniles (IJs) from the dead T. molitor larva. The right panel shows an enlarged picture of the square region in the left panel as well as IJs. (B) Morphological traits of male and female specimens: TL for tail length, RL for rectum length, ovary, BDM for body diameter, BL for body length, and SL for stylet length. (C) Genus identification of the nematode isolate using the morphological characters. Each measurement used 10 nematode samples and was reported as mean ± SD. (D) Species identification of the nematode isolate using the morphological characteristics. ‘+’ and ‘-’ in the parenthesis indicate being identical and non-identical with the nematode isolate with respect to each morphological characteristic.
Figure 2.
Molecular identification of a nematode isolate killing insects using ITS sequence. (A) Nucleotide sequence of internal transcribed spacer (‘ITS’) region of the nematode isolate. The sequence was submitted with the GenBank access number of PP358559. With sequence comparison of the ITS sequences of Heterorhabditis spp. (Adams et al., 1998), the following putative rDNA regions were identified: partial 18S rDNA (1-161), ITS-1 (162-470), 5.8S rDNA (470-600), ITS-2 (600-843), and partial 28S rDNA (843-929). (B) Phylogenetic analysis of the ITS sequence of the nematode isolate with those of the closely related species. The tree was constructed with MEGA.11 based on Neighbor-Joining analysis, in which the bootstrap values were obtained with 1,000 repetitions.
Figure 2.
Molecular identification of a nematode isolate killing insects using ITS sequence. (A) Nucleotide sequence of internal transcribed spacer (‘ITS’) region of the nematode isolate. The sequence was submitted with the GenBank access number of PP358559. With sequence comparison of the ITS sequences of Heterorhabditis spp. (Adams et al., 1998), the following putative rDNA regions were identified: partial 18S rDNA (1-161), ITS-1 (162-470), 5.8S rDNA (470-600), ITS-2 (600-843), and partial 28S rDNA (843-929). (B) Phylogenetic analysis of the ITS sequence of the nematode isolate with those of the closely related species. The tree was constructed with MEGA.11 based on Neighbor-Joining analysis, in which the bootstrap values were obtained with 1,000 repetitions.
Figure 3.
Biological characteristics of O. tipulae. (A) Parasitic life cycle of O. tipulae in M. vitrata larvae. Parasitized host larvae with O. tipulae IJs (= L3) were kept at room temperature and dissected every day to assess the nematode development of different larval instars (L1-L4). Each measurement consists of three replications (= three parasitized larvae). Relative proportions were calculated by the ratios of individual numbers among different developmental stages. (B) Comparison of IJ yield between M. vitrata (‘Mv’) and T. molitor (‘Tm’) larvae. IJs were collected using a White trap for 3 days after the first IJ emergence from the cadaver. Each measurement consisted of three replications, and the error bar indicates the standard deviation. The different letters above the standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). (C) Survival analysis of IJs at different temperatures and preservatives for a storage period of 30 days. The different letters following the survival means at 30 days indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 3.
Biological characteristics of O. tipulae. (A) Parasitic life cycle of O. tipulae in M. vitrata larvae. Parasitized host larvae with O. tipulae IJs (= L3) were kept at room temperature and dissected every day to assess the nematode development of different larval instars (L1-L4). Each measurement consists of three replications (= three parasitized larvae). Relative proportions were calculated by the ratios of individual numbers among different developmental stages. (B) Comparison of IJ yield between M. vitrata (‘Mv’) and T. molitor (‘Tm’) larvae. IJs were collected using a White trap for 3 days after the first IJ emergence from the cadaver. Each measurement consisted of three replications, and the error bar indicates the standard deviation. The different letters above the standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). (C) Survival analysis of IJs at different temperatures and preservatives for a storage period of 30 days. The different letters following the survival means at 30 days indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 4.
Insecticidal spectrum of the infective juveniles (IJs) of O. tipulae. (A) Virulence of O. tipulae against larvae of four different insects: T. molitor (‘Tm’), S. exigua (‘Se’), P. xylostella (‘Px’), and M. vitrata (‘Mv’). The experimental unit was a Petri dish containing 10 larvae, and each one was replicated three times. Mortality was recorded at 4 days after treatment (‘DAT’). (B) Medium lethal time (LT50) at a concentration of 400 IJs per larva. (C) Pathogenicity of O. tipulae against different developmental stages of T. molitor at 400 IJs per larva. Different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 4.
Insecticidal spectrum of the infective juveniles (IJs) of O. tipulae. (A) Virulence of O. tipulae against larvae of four different insects: T. molitor (‘Tm’), S. exigua (‘Se’), P. xylostella (‘Px’), and M. vitrata (‘Mv’). The experimental unit was a Petri dish containing 10 larvae, and each one was replicated three times. Mortality was recorded at 4 days after treatment (‘DAT’). (B) Medium lethal time (LT50) at a concentration of 400 IJs per larva. (C) Pathogenicity of O. tipulae against different developmental stages of T. molitor at 400 IJs per larva. Different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 5.
Entomopathogenic bacteria isolated from O. tipulae (‘Ot’). (A) Virulence of three different colonies of small (‘S’), big (‘B’), or fuzzy (‘F’) shapes against M. vitrata larvae. The bacteria (103 cfu/larva) were injected into the larvae of M. vitrata. The experimental unit was a Petri dish containing 10 larvae and each unit was replicated three times. Mortality was recorded at 48 h after the treatment. (B) Identification of the bacterial colonies using 16S rRNA sequence. (C) Influence of eicosanoids on the virulence of the nematode against two insect hosts: M. vitrata (‘Mv’) and T. molitor (‘Tm’). Ot was treated at 400 IJs per larva to the two insects. Dexamethasone (‘DEX’) or arachidonic acid (‘AA’) was injected into the larvae at 1 μg per larva. The experimental unit was a Petri dish containing 10 larvae, and each unit was replicated three times. Mortality was recorded at 4 days after the nematode treatment (‘DAT’). The different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 5.
Entomopathogenic bacteria isolated from O. tipulae (‘Ot’). (A) Virulence of three different colonies of small (‘S’), big (‘B’), or fuzzy (‘F’) shapes against M. vitrata larvae. The bacteria (103 cfu/larva) were injected into the larvae of M. vitrata. The experimental unit was a Petri dish containing 10 larvae and each unit was replicated three times. Mortality was recorded at 48 h after the treatment. (B) Identification of the bacterial colonies using 16S rRNA sequence. (C) Influence of eicosanoids on the virulence of the nematode against two insect hosts: M. vitrata (‘Mv’) and T. molitor (‘Tm’). Ot was treated at 400 IJs per larva to the two insects. Dexamethasone (‘DEX’) or arachidonic acid (‘AA’) was injected into the larvae at 1 μg per larva. The experimental unit was a Petri dish containing 10 larvae, and each unit was replicated three times. Mortality was recorded at 4 days after the nematode treatment (‘DAT’). The different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 6.
Metagenomics analysis of the gut bacteria isolated from the nematodes of O. tipulae (‘Ot’) and two S. carpocapsae (‘Sc-AD’ and ‘Sc-DR’). (A) Comparison of bacterial diversity based on Amplicon Sequence Variants (ASVs). The figures within the bars indicate the total numbers of the bacterial species based on ASVs. (B) Relative bacterial abundance in taxonomic composition at phylum level. (C) Comparative analysis of bacterial diversity on two community indices: Shannon index and Gini-Simpson index. (D) Comparison of the entomopathogenic bacteria identified from the nematodes. Presence and absence are respectively denoted by + and -.
Figure 6.
Metagenomics analysis of the gut bacteria isolated from the nematodes of O. tipulae (‘Ot’) and two S. carpocapsae (‘Sc-AD’ and ‘Sc-DR’). (A) Comparison of bacterial diversity based on Amplicon Sequence Variants (ASVs). The figures within the bars indicate the total numbers of the bacterial species based on ASVs. (B) Relative bacterial abundance in taxonomic composition at phylum level. (C) Comparative analysis of bacterial diversity on two community indices: Shannon index and Gini-Simpson index. (D) Comparison of the entomopathogenic bacteria identified from the nematodes. Presence and absence are respectively denoted by + and -.
Figure 7.
High virulence of O. tipulae against dipteran insect pests of D. platura (‘Dp’), D. melanogaster (‘Dm’), and A. albopictus (‘Aa’). (A) Comparative virulence analysis of the nematodes of O. tipulae (‘Ot’) and two S. carpocapsae (‘Sc-AD’ and ‘Sc-DR’) against different insect hosts including lepidopteran (‘Lep’), coleopteran (‘Col’), and dipteran (‘Dip’) species. The Lep and Col insects were M. vitrata (‘Mv’) and T. molitor (‘Tm’) larvae, respectively. All larvae were treated at 400 IJs per larva. The experimental unit was a Petri dish (each well of a 24-well plate) containing 10 larvae, and each unit was replicated three times. Mortality was recorded at 4 days after the nematode treatment (‘DAT’). Asterisks above the standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). ‘ns’ represents no significant difference. (B) Ot virulence against three Dip insects. Each treatment used 10 larvae and was replicated three times. Different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). Photos indicate the emerging IJs from the cadavers.
Figure 7.
High virulence of O. tipulae against dipteran insect pests of D. platura (‘Dp’), D. melanogaster (‘Dm’), and A. albopictus (‘Aa’). (A) Comparative virulence analysis of the nematodes of O. tipulae (‘Ot’) and two S. carpocapsae (‘Sc-AD’ and ‘Sc-DR’) against different insect hosts including lepidopteran (‘Lep’), coleopteran (‘Col’), and dipteran (‘Dip’) species. The Lep and Col insects were M. vitrata (‘Mv’) and T. molitor (‘Tm’) larvae, respectively. All larvae were treated at 400 IJs per larva. The experimental unit was a Petri dish (each well of a 24-well plate) containing 10 larvae, and each unit was replicated three times. Mortality was recorded at 4 days after the nematode treatment (‘DAT’). Asterisks above the standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). ‘ns’ represents no significant difference. (B) Ot virulence against three Dip insects. Each treatment used 10 larvae and was replicated three times. Different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). Photos indicate the emerging IJs from the cadavers.
Figure 8.
Control efficacies of O. tipulae (‘Ot’) against D. platura in pot (A) and field (B) assays. In pot assay, each Welsh onion was grown in a pot and inoculated with 30 nematodes. Different concentrations of the nematodes were applied to the pot. The nematode mortality was recorded at six days after treatment (‘DAT’). In the field assay, each onion plant infested by D. platura was drenched with the nematode suspension containing 1,000 IJs. Fluxametamide (Flx, 500 ppm) was treated as a positive control in the experiment. Control efficacy was estimated by comparing survival IJ numbers in untreated and treated plants. The different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Figure 8.
Control efficacies of O. tipulae (‘Ot’) against D. platura in pot (A) and field (B) assays. In pot assay, each Welsh onion was grown in a pot and inoculated with 30 nematodes. Different concentrations of the nematodes were applied to the pot. The nematode mortality was recorded at six days after treatment (‘DAT’). In the field assay, each onion plant infested by D. platura was drenched with the nematode suspension containing 1,000 IJs. Fluxametamide (Flx, 500 ppm) was treated as a positive control in the experiment. Control efficacy was estimated by comparing survival IJ numbers in untreated and treated plants. The different letters above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test).
Table 1.
Virulence of the three bacterial colonies isolated from O. tipulae against L5 larvae of M. vitrata. The virulence was estimated using the median lethal dose (LD50). Each bioassay used 10 larvae at the L5 stage and was replicated three times.
Table 1.
Virulence of the three bacterial colonies isolated from O. tipulae against L5 larvae of M. vitrata. The virulence was estimated using the median lethal dose (LD50). Each bioassay used 10 larvae at the L5 stage and was replicated three times.
Colony shape |
Number of larvae used for bioassay |
LD50 (95% CI) cfu/larva |
Slope ± SE |
X2 |
Fuzzy |
50 |
249.3 (89.1-698.6) |
0.49 ± 0.1 |
8.37 |
Small |
50 |
1,649.7 (667-4,350.4) |
0.47 ± 0.1 |
8.94 |
Big |
50 |
8,313.9 (3,564.4-2,4477.0) |
0.53 ± 0.2 |
7.96 |