Phylogenetic evidence for the invasion of a commercialized European Phasmarhabditis hermaphrodita lineage into North America and New Zealand

Abstract

Biological control (biocontrol) as a component of pest management strategies reduces reliance on synthetic chemicals, and seemingly offers a natural approach that minimizes environmental impact. However, introducing a new organism to new environments as a classical biocontrol agent can have broad and unanticipated biodiversity effects and conservation consequences. Nematodes are currently used in a variety of commercial biocontrol applications, including the use of Phasmarhabditis hermaphrodita as an agent targeting pest slug and snail species. This species was originally discovered in Germany, and is generally thought to have European origins. P. hermaphrodita is sold under the trade name Nemaslug®, and is available only in European markets. However, this nematode species was discovered in New Zealand and the western United States, though its specific origins remained unclear. In this study, we analyzed 45 nematode strains representing eight different Phasmarhabditis species, collected from nine countries around the world. A segment of nematode mitochondrial DNA (mtDNA) was sequenced and subjected to phylogenetic analyses. Our mtDNA phylogenies were overall consistent with previous analyses based on nuclear ribosomal RNA (rRNA) loci. The recently discovered P. hermaphrodita strains in New Zealand and the United States had mtDNA haplotypes nearly identical to that of Nemaslug®, and these were placed together in an intraspecific monophyletic clade with high support in maximum likelihood and Bayesian analyses. We also examined bacteria that co-cultured with the nematode strains isolated in Oregon, USA, by analyzing 16S rRNA sequences. Eight different bacterial genera were found to associate with these nematodes, though Moraxella osloensis, the bacteria species used in the Nemaslug® formulation, was not detected. This study provided evidence that nematodes deriving from the Nemaslug® biocontrol product have invaded countries where its use is prohibited by regulatory agencies and not commercially available.

About the Authors:

Dana K. Howe

Roles Formal analysis, Investigation, Methodology, Project administration, Resources, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Integrative Biology, Oregon State University, Corvallis, Oregon, United States of America

ORCID logo http://orcid.org/0000-0002-2131-6114

Anh D. Ha

Roles Formal analysis, Investigation, Writing – original draft

Affiliation: Department of Integrative Biology, Oregon State University, Corvallis, Oregon, United States of America

Andrew Colton

Roles Investigation, Resources, Writing – review & editing

Affiliation: Department of Crop and Soil Science, Oregon State University, Corvallis, Oregon, United States of America

Irma Tandingan De Ley

Roles Resources, Writing – review & editing

Affiliation: Department of Nematology, University of California-Riverside, Riverside, California, United States of America

Robbie G. Rae

Roles Resources, Writing – review & editing

Affiliation: School of Biological and Environmental Sciences, Liverpool John Moores University, Liverpool, United Kingdom

Jenna Ross

Roles Resources, Writing – review & editing

Affiliations Department of Conservation Ecology and Entomology, Faculty of AgriSciences, Stellenbosch University, Matieland, South Africa, Institute of Biological and Environmental Sciences, University of Aberdeen, Aberdeen, United Kingdom, Crop Health and Protection (CHAP), National Agri-Food Innovation Campus, Sand Hutton, York, United Kingdom

Michael Wilson

Roles Resources, Writing – review & editing

Affiliation: Independent Researcher/Consultant, Hamilton, New Zealand

Jiří Nermut

Roles Resources, Writing – review & editing

Affiliation: Biology Centre CAS, Institute of Entomology, Branišovská, České Budějovice, Czech Republic

Zhongying Zhao

Roles Resources, Writing – review & editing

Affiliation: Department of Biology, Hong Kong Baptist University, Hong Kong SAR, China

Rory J. Mc Donnell

Roles Conceptualization, Writing – review & editing

Affiliation: Department of Crop and Soil Science, Oregon State University, Corvallis, Oregon, United States of America

ORCID logo http://orcid.org/0000-0003-3415-7295

Dee R. Denver

Roles Conceptualization, Writing – review & editing

Affiliation: Department of Integrative Biology, Oregon State University, Corvallis, Oregon, United States of America

Introduction

Biological control, ‘biocontrol’, is a human-mediated strategy used to suppress populations of pests and other undesirable species that involves the active deployment or conservation of natural enemies of the target species. Biocontrol strategies are common components of integrated pest management (IPM) programs, and are commonly deemed more favorable than chemical-based and other pest control strategies where environmental damage and potential human health impacts are matters of concern [1–3]. The commercialization of biocontrol organisms started in the 1920s with development of Encarsia formosa, a parasitoid wasp species, as an agent used to combat greenhouse whiteflies [4]. Biocontrol products have diversified and increased over recent decades, expanding the overall abundance of these organisms and the geographic range across which they occur.

Although biocontrol strategies are often viewed favorably because they offer more natural approaches as compared to those involving chemicals, biodiversity and conservation concerns are also prominent components of the biocontrol dialogue. In particular, biocontrol approaches that involve the importation of natural enemies to locales where they did not previously exist can have profound unintended ecosystem effects. Introduced species intended for biocontrol applications can negatively impact native, non-target species; most of the well-known cases involve vertebrates. One classic example involved the misguided introduction of the small Indian mongoose, Herpestes javanicus, to Hawaii for the purpose of controlling rats. The mongoose, however, were discovered to be opportunistic feeders that preyed on the eggs and hatchlings of native ground nesting birds, and endangered sea turtles of Hawaii [5,6]. Impacts of introduced species can also be less direct, though still devastating to affected ecosystems. The Hawaiian cane toad, Rhinella marina, was intentionally introduced to Australia in 1935 to control greyback cane beetles, Dermolipida albohirtum, that were a pest in sugarcane production systems. The introduced toads were unable to reach the cane beetles, which resided high on cane stalks beyond the reach of the toads, so instead fed on other insects and rapidly spread across much of northeastern Australia during the subsequent decades. The cane toads took over habitats used by native amphibian species and brought new diseases to the native toads and frogs of Australia that dramatically reduced their populations [7]. Further, the poison produced by the cane toads negatively impacted populations of native predators such as goannas and tiger snakes. The predatory snail Euglandina rosea has been introduced to Pacific and Indian Ocean islands for the purpose of controling giant African land snails (Lissachatina fulica), but has turned out to be ineffective against the target invasive species and damaging to island-endemic gastropods instead [8]. These classic ‘biocontrol gone awry’ and similar stories have motivated the development of rigorous regulatory policies surrounding efforts to develop new biocontrol initiatives, requiring careful environmental impact analyses and other forms of scientific scrutiny prior to biocontrol agent deployment [9–12].

Invertebrates such as parasitic nematode species have also been deployed as biocontrol agents and components of IPM programs, usually targeting insect pests. For example, entomopathogenic nematodes in two genera, Steinernema and Heterorhabditis, parasitize broad ranges of insect hosts and are commonly used to target scarab beetles and weevils [13,14]. These nematodes rely on obligate symbioses with mutualistic bacteria (Xenorhabdus spp. bacteria associate with Steinernema spp. nematodes; Photorhabdus spp. bacteria associate with Heterorhabditis spp. nematodes) to feed on their insect hosts. Deladenus siricidicola offers a second example of a nematode species used in insect biocontrol; this nematode infects Sirex noctilio, a species of woodwasp which causes significant tunnelling damage to the trunks of pine trees [15,16]. D. siricidicola infection results in sterilized female wasps, which also vector the nematodes to nearby trees where they can infect new hosts.

Phasmarhabditis hermaphrodita is a nematode species commercialized for biocontrol applications in Europe [17–19]. Unlike the previously discussed biocontrol nematodes, P. hermaphrodita parasitizes a variety of pest slug and snail species; the most common biocontrol target is the grey field slug, Deroceras reticulatum. P. hermaphrodita is grown in mass production with a bacteria, Moraxella osloensis used as a food source, and sold under the tradenames Nemaslug® and SlugTech®. The relative roles of the nematode and bacteria in causing slug host mortality remain unclear, though the ability of P. hermaphrodita to grow and reproduce by feeding on many different types of bacteria in the lab and other recent findings [20–22] suggest that any putative bacterial symbioses in this system must be facultative in nature. Other species in this nematode genus, such as Phasmarhabditis papillosa and Phasmarhabditis californica, have been proposed as potential additional or alternative biocontrol agents [23,24], though P. hermaphrodita offers the only currently available nematode biocontrol product targeting pest slugs and snails.

P. hermaphrodita was originally discovered in Germany [25] and then later in France [26] and is often assumed to be of European origin. The P. hermaphrodita strain that became Nemaslug® was isolated in Bristol, England [18]. With increased interest in this potential biocontrol there has been increased surveying efforts; P. hermaphrodita has also been discovered in Chile [27], Iran [28], Egypt [29], and more recently in New Zealand [30,31] and the west coast of North America [32,33]. It remains unclear, however, whether these reported observations of P. hermaphrodita in non-European locations constitute previously unknown natural populations of the species, or instead resulted from the introduction of Nemaslug®. The latter possibility would pose a matter of strong conservation and regulatory concern. Nemaslug® is not licensed for sale and has never been sold in the United States or New Zealand, and the potential impacts of P. hermaphrodita on native gastropod species remains unknown.

In this study, we developed and applied a molecular genetic strategy to investigate the origins and diversity of P. hermaphrodita lineages obtained from different parts of the world. Our approach, based on rapidly evolving mitochondrial DNA (mtDNA) sequences, was further applied to additional Phasmarhabditis spp. to provide broader insights into evolutionary relationships in this group of nematodes, and guide tree rooting decisions for phylogenetic analyses focused on P. hermaphrodita. We also examined the bacterial associates of three P. hermaphrodita isolates found in Oregon to add knowledge on the potential relationships between this nematode species and putative bacterial symbionts. Our study demonstrated that mtDNA is a powerful tool for investigating within- and between-species evolutionary patterns in this nematode group, found variable bacteria associated with the three different P. hermaphrodita isolates from Oregon, and revealed strong evidence that recently discovered P. hermaphrodita isolates from California, Oregon, and New Zealand all likely derived from Nemaslug®.

Materials and methods

Nematode populations

Nematodes were isolated from slugs collected from field and nursery sites in Oregon. Briefly, collected slugs were maintained in a laboratory growth chamber at 16.5°C and 12:12 light in 16oz plastic containers (9cm diameter), corresponding to the site of origin. A maximum of 20 slugs were placed in each container to avoid stress or cross contamination. Each container was kept damp with deionized water-soaked towels, and the slugs were fed organic carrot. Dead slugs, and slugs showing signs of illness, were removed and placed in a petridish (5cm diameter) for further monitoring. Upon death, all slugs were visually inspected using a Leica microscope for the presence of nematodes. If detected, nematodes were washed off the slug carcass using M9 buffer and transferred to standard NGM plates for culturing. After 7–10 days, nematodes were collected from the plates for DNA extraction and sequencing. Additional Phasmarhabditis species were isolated from various locations around the globe (Table 1). These nematodes were preserved in DESS (a solution of dimethyl sulphoxide, disodium EDTA, and saturated NaCl) [34] and shipped to Oregon for DNA extraction and sequencing. Most locations did not require permits for slug collection, however, the slugs collected in California fell under California Department of Food and Agriculture (CDFA) permit 2942 (October 2012 to August 2014).

[Figure omitted. See PDF.]

Table 1. Phasmarhabditis spp. strains included in this study.

https://doi.org/10.1371/journal.pone.0237249.t001

PCR and DNA sequencing

Nematodes were collected in a proteinase-K based lysis buffer for DNA extraction as previously described [35]. New nematode isolates’ species were identified using 18S ribosomal DNA amplification and direct sequencing, and collaborator nematode samples were verified using 18S rRNA amplification and direct sequencing [36]. All 18S rRNA sequences were compared to GenBank’s nr database using blastn. When the BLAST match was a Phasmarhabditis species (percent identity of 99–100), the population was incorporated into this mtDNA study.

PCR amplification of ~5 kb of mtDNA was achieved using the Expand Long Range PCR kit (Roche) as previously described [37]. The mtDNA was targeted with ribosomal primers designed to match across a wide range of rhabditid nematodes. This segment of mtDNA included the genes: rrnL, nad3, nad5, tRNA-Ala, tRNA-Pro, tRNA-Val, nad6, nad4L, tRNA-Trp, tRNA-Glu, and rrnS (S1 Fig). Direct end sequencing of the large mtDNA products were generated from all Phasmarhabditis species. For P. hermaphrodita and Phasmarhebditis neopapillosa samples, subsequent sequencing of the large mtDNA products was performed using newly designed genus-specific primers (S1 Table). All newly generated DNA sequences for this study were submitted to GenBank under accession numbers MT179848 –MT179854 for 18S sequences and MT472183 –MT472270 for mtDNA sequences.

Phylogenetic analysis

Two phylogenetic analyses were performed using two different mtDNA sequence data sets. The first analysis targeted the relationships within the broader genus Phasmarhabditis and included 45 nematode strains from eight species, along with three nematode outgroup strains: Pristionchus pacificus NC_015245 [38], Oscheius chongmingensis KP257594 and Caenorhabditis elegans NC_001328 [39]. The second analysis explored the relationship within P. hermaphrodita more closely using a longer region of mtDNA sequences from 24 nematode strains collected from four countries and four P. neopapillosa strains, the nearest outgroup. MEGA7 was used for sequence alignment and phylogenetic model testing of all data sets [40].

For the first analysis, the best fit model for the comparison of 1,121 bp of mtDNA from multiple Phasmarhabditis species was the Hasegawa-Kishino-Yano model, with Gamma distribution and invariable sites (HKY+G+I) [41]. This model was used to generate a maximum-likelihood molecular evolution phylogeny with 1,000 bootstrap replicates in MEGA7. A Bayesian inference phylogeny was also generated with model HKY85+G+I using the MrBayes [42] plugin in Geneious Prime 2019.1.3 (https://www.geneious.com).

For the second analysis, the phylogeny of 2,372 bp of P. hermaphrodita mtDNA, the best fit model was the Hasegawa-Kishino-Yano model, with Gamma distribution. A HKY+G maximum-likelihood evolution phylogeny with 1,000 bootstrap replicates was generated in MEGA7. A Bayesian inference phylogeny was also generated with model HKY85+G using the MrBayes plugin in Geneious Prime.

Isolation and genetic characterization of bacteria

Nematodes of three Oregon P. hermaphrodita strains (DL300, DL307, and DL309) were cultured on NGM plates until a bacterial lawn established. Bacteria were then streaked onto a series of sterile LB nutrient agar plates to obtain single colonies. Twelve colonies from each nematode strain were collected in molecular biology-grade water (product number). DNA was isolated by incubating at 100°C in 10 minutes then immediately cooling on ice for 10 minutes.

Identification of bacteria was performed using 16S ribosomal DNA amplification and direct sequencing. Sequences of length ~1,460 bp in the 16S gene were amplified using the universal primers 27F and 1492R [43,44] (sequences: 27F (5’-AGAGTTTGATCCTGGCTCAG-3’) and 1492R (5’-GGTTACCTTGTTACGACTT-3’)). PCR products were then purified using ChargeSwitch PCR Cleanup Kit (Thermo Fisher Scientific). DNA sequences were determined by Sanger method at the Center for Genome Research and Biocomputing at Oregon State University. Sequences were compared to the NCBI nr database using BLASTN, and the first match with percent identity of ≥ 97% were reported as sequence ID.

Results

Phylogenetics of mtDNA in Phasmarhabditis spp

For the first analysis, a maximum likelihood with bootstrap analysis approach was used to phylogenetically analyze 1,121 bp of mtDNA from 48 nematode strains from eight different Phasmarhabditis species (Fig 1). The Bayesian inference analysis yielded a tree with an identical topology as maximum likelihood, therefore the results were combined. All Phasmarhabditis species formed a monophyletic group with strong support (99/1.00). The mtDNA sequences provide strong support, in both analysis methods, for the previously defined species. For six of the species, we analyzed multiple nematode mtDNA samples, and in every case these species were represented as unique monophyletic groups. The unnamed species, first discovered in earthworms and recently now found in slugs, is a strongly supported internal clade of Phasmarhabditis. Using mtDNA sequence, P. hermaphrodita is most closely related to P. neopapillosa and most distantly related to Phasmarhabditis huizhouensis.

[Figure omitted. See PDF.]

Fig 1. Maximum likelihood phylogeny of ~1,120 bp of mtDNA from Phasmarhabditis species.

A total of 48 nematode mtDNA sequences, representing 8 Phasmarhabditis species and 3 outgroups, were analyzed. Additional information about the nematode populations is found in Table 1. Maximum likelihood and Bayesian phylogenies yielded identical topologies; node support lists maximum likelihood values on left (ML) and Bayesian posterior probabilities are on the right (Bay). The scale bar represents 0.050 substitutions per site.

https://doi.org/10.1371/journal.pone.0237249.g001

Phylogenetics of mtDNA within P. hermaphrodita

For the second analysis, additional mtDNA sequencing, for a total of 2,372 bp, was used in a maximum likelihood with bootstrap analysis approach to phylogenetically analyze 24 populations of P. hermaphrodita, with four populations of P. neopapillosa as an outgroup (Fig 2). A Bayesian analysis yielded an identical topology as that of maximum likelihood, therefore the results were combined. Within the P. hermaphrodita clade, there was strong support for two major subclades, defined here as Ph_mtI (95/0.84) and Ph_mtII (100/1.00). Subclade Ph_mtI contained the Nemaslug® strain, originally isolated in Europe, and all nematodes recently isolated from New Zealand and the United States of America; with one and nine strains, respectively. Ph_mtII contains only strains isolated from Europe, both England and the Czech Republic; nine and four strains respectively.

[Figure omitted. See PDF.]

Fig 2. Maximum likelihood phylogeny of ~2,370 bp of mtDNA from P. hermaphrodita and P. neopapillosa.

A total of 28 nematode mtDNA sequences, originating from five countries, were analyzed. Country designations are CzR = Czech Republic; Eng = England; NZ = New Zealand; Sct = Scotland; USO & USC = United States of America. P. hermaphrodita populations are split into two subclades, Ph_mtI and Ph_mtII. Maximum likelihood and Bayesian phylogenies yielded identical topologies; node support lists maximum likekihood values on left (ML) and Bayesian posterior probabilities are on the right (Bay). The scale bar represents 0.020 substitutions per site.

https://doi.org/10.1371/journal.pone.0237249.g002

Bacterial associates of P. hermaphrodita from Oregon

Twelve bacterial colonies found in association with each of the three Oregon P. hermaphrodita were examined (Fig 3). Out of 12 colonies from the strain DL300, four were identified as bacteria of the genus Sphingobacterium (33%), and the other eight matched to Pseudomonas (77%). Five bacterial genera were found associated with the strain DL307, namely Acinetobacter, Brucella, Microbacterium, Pseudomonas, and Ochrobactrum. For the strain DL309, the three genera Acinetobacter, Pseudomonas, and Stenotromophonas were detected.

[Figure omitted. See PDF.]

Fig 3. Bacterial associates of Oregon P. hermaphrodita.

The numbers of bacteria identified from different specific genera are indicated by different colors in the bar chart. All three P. hermaphrodita strains analyzed were discovered in Oregon (Table 1).

https://doi.org/10.1371/journal.pone.0237249.g003

Pseudomonas bacteria were observed in association with all three strains of nematodes, accounting for 77%, 33%, and 33% of all colonies analyzed in DL300, DL307, and DL309, respectively. Bacteria of the genus Moraxella were not found in any of the colonies examined; however, Acinetobacter, another genus of the same family as Moraxella (i.e. Moraxellaceae) was found associated with both strains DL307 and DL309 (8.3% and 50% of all colonies, respectively). Additional information about the sequence BLASTN matches is provided in S2 Table.

Discussion

The phylogenetic analyses presented in this study revealed strong support for two major intraspecific mtDNA clades of P. hermaphrodita (Figs 1, 2). One clade, Ph_mtI, included all nematode strains isolated from slug hosts in North America (California, Oregon) and New Zealand; this clade also included the Nemaslug® strain (derived from a nematode strain isolated in Bristol, England [18]). The other clade, Ph_mtII, included P. hermaphrodita strains isolated from slugs in Liverpool, England and Benesov u Prahy and Ceske Budejovice, Czech Republic. These patterns, particularly in Ph_mtI, support the hypothesis that recent discoveries of P. hermaphrodita in North America and New Zealand derived from the Nemaslug® strain. Previous sampling efforts (performed in 2007) failed to reveal the presence of P. hermaphrodita in California, Oregon and other U.S. states [45], suggesting that the arrival of this biocontrol species in North America might be very recent [46,47]. But, we cannot rule out an alternative hypothesis that P. hermaphrodita, prior to the development of the Nemaslug®, was a globally distributed and panmictic species, and that under sampling in the previous and present studies has given rise to an artifactual and misleading result. Natural populations of the famous model nematode Caenorhabditis elegans demonstrate worldwide population-genetic patterns consistent with global panmixia, though this pattern is commonly attributed to the effects of human activities [48,49]. The strong and consistent pattern of phylogeographic mtDNA groupings within P. hermaphrodita, coupled with very low levels of sequence variation within Ph_mtI (Fig 2), better support the earlier hypothesis that the North American and New Zealand strains originated from the same strain that was developed into Nemaslug®.

Although the three P. hermaphrodita strains discovered in Oregon shared a common mtDNA haplotype that was nearly identical to the Nemaslug® haplotype, different genera of bacteria were found to co-culture with each of these three strains. In our screen, we found only one bacterial genus common to all three nematode samples (Pseudomonas), and one other genus that was found in two out of three nematode samples (Acinetobacter). For the other six bacterial genera identified, all were unique to a single Oregon nematode strain. Although only limited microbiome inferences are possible from examining small numbers of bacterial colonies able to grow on nutrient agar plates, the patterns suggest that each of the three Oregon P. hermaphrodita strains were associated with multiple and distinct bacterial genera at the time of nematode collection. This pattern is consistent with previous research demonstrating that this species associates with diverse bacterial partners [20–22]. These preliminary characterizations of the Oregon nematode microbiomes must be interpreted with caution. An analysis using high-throughput DNA sequencing-based strategies for microbiome characterization will be required to thoroughly understand the nature of microbial communities associated with P. hermaphrodita.

It is also noteworthy that we did not identify Moraxella osloensis, the bacterial species used in the Nemaslug® commercial formulation. Unlike entomopathogenic nematodes (EPNs) that have a strict symbiotic relationship with bacteria (Steinernema sp. with Xenorhabdus sp. and Heterorhabditis sp. with Photorhabdus sp.) that they rely on to kill insects [50], the relationship P. hermaphrodita has with bacteria is more complex. Initial studies focused on finding a bacterium that P. hermaphrodita could be grown on and produce pathogenic nematodes. Bacteria were isolated from infected slugs and from P. hermaphrodita emerging from dead slugs [20,21]. Many different bacterial species were isolated and tested including: Acinetobacter calcoaceticus, Aeromonas hydrophila, Aeromonas sp., Bacillus cereus, Flavobacterium breve, FIavobacterium odoratum, Moraxella osloensis, Providencia rettgeri, Pseudomonas fluorescens (isolate no. 1a), P. fluorescens (isolate no. 140), P. fluorescens (isolate no. 141), P. fluorescens (pSG), P. paucimobilis, Serratia proteamaculans, Sphingobacterium spiritocorum and Xenorhabdus bovienii. Moraxella osloensis was chosen as it produced consistently high yields of pathogenic nematodes [20,21]. It should be stressed that this bacterium was chosen for commercial production and does not reflect the natural tritrophic interactions that may be occurring between slugs, P. hermaphrodita and bacteria in the wild. A study by Rae et al. showed that P. hermaphrodita, when grown on rotting slugs or emerging after parasitising slugs (D. reticulatum), had no evidence of M. osloensis being present [22]. Similarly, Nermut’ et al. found that P. hermaphrodita strain (DMG0001) lost M. osloensis after repeated culturing [51]. Therefore, it is perhaps not too surprising that wild isolated P. hermaphrodita from the US do not retain M. osloensis as there is little evidence to show this bacterium is vertically transmitted.

Our approach demonstrated the effectiveness of mtDNA in revealing informative phylogenetic patterns within and between species in the genus Phasmarhabditis, similar to past successes in other nematode genera [37,52,53]. The mtDNA phylogenetic framework presented in this study for Phasmarhabditis spp. (Fig 1) revealed numerous topological consistencies when compared to other recent phylogenetic analyses of this genus that relied on nuclear ribosomal RNA (rRNA) genes [54–60]. For example, our analysis identified P. huizhouensis as an early-diverging lineage that is genetically distinct from the other seven species from the genus analyzed here, a pattern consistent with previous analyses of 18S and 28S rRNA loci [57]. Our mtDNA phylogeny revealed an affinity between P. hermaphrodita and P. neopapillosa, consistent with previous nuclear rRNA studies as well [54]. P. californica, a species recently discovered in North America [55], grouped with Phasmarhabditis bonaquaense, consistent with previous nuclear rRNA results [54]. P. papillosa grouped with Phasmarhabditis apuliae, offering a third example of consistency between our mtDNA phylogeny and nuclear rRNA-based analyses [56]. We also note here the discovery of a nematode strain in Oregon (DL308) isolated from Deroceras reticulatum that strongly grouped with two other nematode strains (DF5056, EM434) that were isolated from earthworms [61]. This observation suggests that this unnamed nematode species is capable of associating with different invertebrate hosts.

The Nemaslug® product was originally formulated and commercialized in the United Kingdom, and is available in many European countries. More recently, Dudutech commercialised a similar product under the tradename SlugTech®. However, these biocontrol agents are not commercially available in the United States or New Zealand. Their sale is prohibited by regulatory agencies, primarily due to the unknown impact that P. hermaphrodita might have on native gastropods in these countries. Our findings suggest that, despite these regulatory restrictions, nematodes of the same strain as that of the Nemaslug® product have arrived in California, Oregon, and New Zealand. All discoveries of P. hermaphrodita in North America and New Zealand, to date, derived from nematodes infecting Deroceras reticulatum and other invasive, pest slug species (Table 1). In order to evaluate the potential ecological impacts of these biocontrol-commercialized nematodes, future research is necessary to determine the susceptibilities of New Zealand and North America-native gastropod species to P. hermaphrodita. Broader surveys for the nematode, focused on potential native gastropod hosts, are required to assess whether or not P. hermaphrodita is currently infecting native gastropod species in North America, New Zealand, and other parts of the world.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. Schematic of amplified mitochondrial DNA region for sequence analysis (~5,000 bp).

Linear representation of the fragment of mtDNA amplified using Expand Long Range PCR kit (Roche) with protein-coding genes, ribosomal RNA genes and tRNA genes noted. PCR primers, shown as large arrows, and internal species-specific primers were used for sequencing (see S1 Table). Dashed black lines indicate sequences included in Fig 1 analysis; solid black lines indicate sequences included in Fig 2 analysis.

https://doi.org/10.1371/journal.pone.0237249.s001

(TIF)

S1 Table. PCR and sequencing primers used to target Phasmarhabditis hermaphrodita mtDNA.

https://doi.org/10.1371/journal.pone.0237249.s002

(DOCX)

S2 Table. Bacteria associated with Oregon P. hermaphrodita.

https://doi.org/10.1371/journal.pone.0237249.s003

(DOCX)

Acknowledgments

Thanks to the OSU Center for Genome Research and Biocomputing for DNA sequencing service. Thanks to the following people for assisting with nematode isolation and culturing: Marisa Lutz, Anton Alvarez, Yifei Qiu, Rene Huang, James Cutler, Antoinette Malan, Annika Pieterse and Louwrens Tiedt.

Citation: Howe DK, Ha AD, Colton A, De Ley IT, Rae RG, Ross J, et al. (2020) Phylogenetic evidence for the invasion of a commercialized European Phasmarhabditis hermaphrodita lineage into North America and New Zealand. PLoS ONE 15(8): e0237249. https://doi.org/10.1371/journal.pone.0237249

References

1. Hajek AE, Hurley BP, Kenis M, Garnas JR, Bush SJ, Wingfield MJ, et al. Exotic biological control agents: A solution or contribution to arthropod invasions? Biol Invasions. 2016;18: 953–969.

2. De Clercq P, Mason PG, Babendreier D. Benefits and risks of exotic biological control agents. BioControl. 2011;56: 681–698.

3. Simberloff D. Risks of biological control for conservation purposes. BioControl. 2012: 57: 263. https://doi.org/10.1007/s10526-011-9392-4.

4. Hoddle MS, Van Driesche RG, Sanderson JP. Biology and use of the whitefly parasitoid Encarsia formosa. Annu Rev Entomol. 1998;43(1): 645–669. pmid:15012401

5. Simberloff D, Dayan T, Jones C, Ogura G. Character displacement and release in the small Indian mongoose, Herpestes javanicus. Ecology. 2000;81: 2086–2099.

6. Hays WST, Conant S. Biology and impacts of Pacific Island invasive species. 1. A worldwide review of effects of the small Indian mongoose, Herpestes javanicus (Carnivora: Herpestidae). Pacific Science. 2007;61:3–16.

7. Crossland MR, Brown GP, Anstis M, Shilton CM, Shine R. Mass mortality of native anuran tadpoles in tropical Australia due to the invasive cane toad (Bufo marinus) Biol Conserv. 2008;141: 2387–2394.

8. Wallace MM III, Cowie RH. Feeding Preferences of Two Predatory Snails Introduced to Hawaii and Their Conservation Implications Malacologia, 2010;53(1):135–144.

9. Simberloff D, Stiling P. How Risky is Biological Control? Ecology. 1996;77(7): 1965–1974.

10. Dlugosch KM, Parker IM. Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol. 2007;17(1): 431–449. pmid:17908213

11. Luong LT, Mathot KJ. Facultative parasites as evolutionary stepping-stones towards parasitic lifestyles. Biol Lett. 2019;15(4). http://doi.org/10.1098/rsbl.2019.0058.

12. Roderick GK, Hufbauer R, Navajas M. Evolution and biological control. Evol Appl. 2012; 5(5): 419–423. pmid:22949918

13. Kaya HK, Gaugler R. Entomopathogenic nematodes. Annu. Rev. Entornol. 1993;38:181–206.

14. Dillman AR, Chaston JM, Adams BJ, Ciche TA, Goodrich-Blair H, Stock SP, et al. An Entomopathogenic Nematode by Any Other Name. PLoS Pathog. 2012;8(3): e1002527. pmid:22396642

15. Bedding RA. Biological Control of Sirex noctilio using the nematode Deladenus siricidicola. In: Bedding RA, Akhurst RJ, Kaya HK, editors. Nematodes and the Biological Control of Insect Pests. CSIRO Publications; 1993. pp. 11–20.

16. Fitza KNE, Garnas JR, Lombardero MJ, Ayres MP, Krivak-Tetley FE, Ahumada R, et al. The global diversity of Deladenus siricidicola in native and non-native populations. Biol Control. 2019;132: 57–65. https://doi.org/10.1016/j.biocontrol.2019.01.014.

17. Wilson MJ, Glen DM, George SK, Butler RC. Mass cultivation and storage of the rhabditid nematode Phasmarhabditis hermaphrodita, a biocontrol agent for slugs. Biocontrol Sci Technol. 1993;3(4): 513–521.

18. Wilson MJ, Glen DM, George SK. The rhabditid nematode Phasmarhabditis hermaphrodita as a potential biological control agent for slugs. Biocontrol Sci Technol. 1993;3(4): 503–511.

19. Wilson M, Rae R. Phasmarhabditis hermaphrodita as a Control Agent for Slugs. In: Campos-Herrera Reditor. Nematode Pathogenesis of Insects and Other Pests. Sustainability in Plant and Crop Protection. Springer, Cham; 2015. pp. 509–521.

20. Wilson MJ, Glen DM, Pearce JD, Rodgers PB. Monoxenic culture of the slug parasite Phasmarhabditis hermaphrodita (Nematoda: Rhabditdae) with different bacteria in liquid and solid phase. Fundam. Appl. Nematol. 1995;18(2): 159–166.

21. Wilson MJ, Glen DM, George SK, Pearce JD. Selection of a bacterium for the mass production of Phasmarhabditis hermaphrodita (Nematoda: Rhabditidae) as a biocontrol agent for slugs. Fundam. Appl. Nematol. 1995;18(5): 419–425.

22. Rae RG, Tourna M, Wilson MJ. The slug parasitic nematode Phasmarhabditis hermaphrodita associates with complex and variable bacterial assemblages that do not affect its virulence. J Invertebr Pathol. 2010;104(3): 222–226. pmid:20433847

23. Pieterse A, Tiedt LR, Malan AP, Ross JL. First record of Phasmarhabditis papillosa (Nematoda: Rhabditidae) in South Africa, and its virulence against the invasive slug, Deroceras panormitanum. Nematology. 2017;19(9): 1035–1050. https://doi.org/10.1163/15685411-00003105.

24. Tandingan De Ley I, Holovachov O, Mc Donnell RJ, Bert W, Paine TD, De Ley P. Description of Phasmarhabditis californica n. sp. and first report of P. papillosa (Nematoda: Rhabditidae) from invasive slugs in the USA. Nematology. 2016;18(2): 175–193. https://doi.org/10.1163/15685411-00002952.

25. Schneider A. Über eine Nematodenlarvae und gewisse Verschiedenheiten in den Geschlechtsorganen der Nematoden. Zeitschrift für wissenschaftliche Zoologie. 1859;10, 176–178.

26. Maupas E. Modes et formes de reproduction des nématodes. Archives de Zoologie Expérimentale et Générale. 1900;8, 463–624.

27. France A, Gerding M. Discovery of Phasmarhabditis hermaphrodita in Chile and its pathological differences with the UK isolate in slug control. J Nematol. 2000;32, 430.

28. Karimi J, Kharazi-Pakadel A, Robert SJ. Report of pathogenic nematodes from slugs, Phasmarhabditis hermaphrodita (Nematoda: Rhabditida) in Iran. Journal of Entomological Society of Iran. 2003;22, 77–78.

29. Genena MAM, Mostafa FAM, Fouly AH, Yousef AA. First record of the slug parasitic nematode, Phasmarhabditis hermaphrodita (Schneider) in Egypt. Archives of Phytopathology and Plant Protection. 2011;44, 340–345.

30. Wilson M, Burch G, Tourna M, Aalders L, Barker G. The potential of a New Zealand strain of Phasmarhabditis hermaphrodita for biological control of slugs. N Z Plant Prot. 2012;65: 161–165. https://doi.org/10.30843/nzpp.2012.65.5388.

31. Glare TR O'Callaghan M. Microbial biopesticides for control of invertebrates: Progress from New Zealand. J Invertebr Pathol. 2019;165: Pages 82–88. pmid:29196233

32. Tandingan De Ley I, McDonnell R, Lopez S, Paine TD, De Ley P. Phasmarhabditis hermaphrodita (Nematoda: Rhabditidae), a potential biocontrol agent isolated for the first time from invasive slugs in North America. Nematology. 2014;16: 1129–38.

33. Mc Donnell RJ, Lutz MS, Howe DK, Denver DR. First report of the gastropod-killing nematode, Phasmarhabditis hermaphrodita, in Oregon, USA. J Nematol. 2018;50(1): 77–78.

34. Yoder M, Tandingan De Ley I, King I, Mundo-Ocampo M, Mann J, Blaxter M, et al. DESS: A versatile solution for preserving morphology and extractable DNA of nematodes. Nematology. 2006;8: 367–376.

35. Williams BD, Schrank B, Huynh C, Shownkeen R, Waterston RH. A genetic mapping system in Caenorhabditis elegans based on polymorphic sequence tagged sites. Genetics. 1992;131: 15.

36. Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, et al. A molecular evolutionary framework for the phylum Nematoda. Nature. 1998;392(6671): 71–5. pmid:9510248

37. Howe DK, Denver DR. Muller’s ratchet and compensatory mutation in Caenorhabditis briggsae mitochondrial genome evolution. BMC Evol Biol. 2008;8: 62. pmid:18302772

38. Molnar RI, Bartelmes G, Dinkelacker I, Witte H, Sommer RJ. Mutation rates and intraspecific divergence of the mitochondrial genome of Pristionchus pacificus. Mol Biol Evol. 2011;28 (8): 2317–2326. pmid:21368317

39. Okimoto R, Macfarlane JL, Clary DO, Wolstenholme DR. The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992;130(3): 471–498. pmid:1551572

40. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7): 1870–4. pmid:27004904

41. Hasegawa M, Kishino H, Yano T. Dating the human-ape split by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22: 160–174. pmid:3934395

42. Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenic trees. Bioinformatics. 2001;17: 754–755. pmid:11524383

43. Lane DJ. 16S/23S rRNA sequencing. In: Stackebrandt E, Goodfellow M, editors. Nucleic acid techniques in bacterial systematics. Chichester: John Wiley & Sons; 1991. pp. 115–175.

44. de Lillo A, Ashley FP, Palmer RM, Munson MA, Kyriacou L, Weightman AJ, et al. Novel subgingival bacterial phylotypes detected using multiple universal polymerase chain reaction primer sets. Oral Microbiol Immunol. 2006;21: 61–68. pmid:16390343

45. Ross JL, Ivanova ES, Severns PM, Wilson MJ. The role of parasite release in invasion of the USA by European slugs. Biol Invasions. 2010;12: 603–610.

46. Grewal SK, Grewal PS, Brown I, Tan L, Hammond RB, Gaugler R. First North American survey for the recovery of nematodes associated with mollusks. J. Nematol. 2000;32: 432.

47. Kaya HK, Mitani DR. Molluscicidal nematodes for the biological control of pest slugs. Slosson Report. 2000;14.

48. Sivasundar A, Hey J. Population Genetics of Caenorhabditis elegans: The Paradox of Low Polymorphism in a Widespread Species. Genetics. 2003;163(1): 147–157. pmid:12586703

49. Andersen E, Gerke J, Shapiro J, Crissman JR, Ghosh R, Bloom JS, et al. Chromosome-scale selective sweeps shape Caenorhabditis elegans genomic diversity. Nat Genet. 2012;44: 285–290. pmid:22286215

50. Forst S, Dowds B, Boemare N, Stackebrandt E. Xenorhabdus and Photorhabdus spp.: bugs that kill bugs. Annu Rev Microbiol. 1997;51, 47–72. pmid:9343343

51. Nermut’ J, Půža V, Mráček Z. The effect of different growing substrates on the development and quality of Phasmarhabditis hermaphrodita (Nematoda: Rhabditidae). Biocontrol Science and Technology. 2014;24(9): 1026–1038.

52. Lewis SC, Dyal LA, Hilburn CF, Weitz S, Liau W-S, LaMunyon CW, et al. Molecular evolution in Panagrolaimus nematodes: origins of parthenogenesis, hermaphroditism and the Antarctic species P. davidi. BMC Evol Biology. 2009;9:15.

53. Howe DK, Smith M, Tom DM, Brown AMV, Peetz AB, Zasada IA, et al. Analysis of nematode-endosymbiont coevolution in the Xiphinema americanum species complex using molecular markers of variable evolutionary rates. Nematology. 2018;21(5): 533–546. https://doi.org/10.1163/15685411-00003233.

54. Nermuť J, Půža V, Mekete T, Mráček Z. Phasmarhabditis bonaquaense n. sp. (Nematoda: Rhabditidae), a new slug-parasitic nematode from the Czech Republic. Zootaxa. 2016;4179(3): 530–546. pmid:27811685

55. Tandingan De Ley I, Holovachov O, Mc Donnell RJ, Bert W, Paine TD, De Ley P. Description of Phasmarhabditis californica n. sp. and first report of P. papillosa (Nematoda: Rhabditidae) from invasive slugs in the USA. Nematology. 2016;18(2): 175–193.

56. Nermuť J, Půža V, Mráček Z. Phasmarhabditis apuliae n. sp. (Nematoda: Rhabditidae), a new rhabditid nematode from milacid slugs. Nematology. 2016;18(9): 1095–1112. https://doi.org/10.1163/15685411-00003017.

57. Huang RE, Ye W, Ren X, Zhao Z. Morphological and molecular characterization of Phasmarhabditis huizhouensis sp. nov. (Nematoda: Rhabditidae), a new rhabditid nematode from South China. PLOS ONE. 2015;10(12): e0144386. pmid:26674768

58. Ross JL, Ivanova ES, Spiridonov SE, Waeyenberge L, Moens M, Nicol GW, et al. Molecular phylogeny of slug-parasitic nematodes inferred from 18S rRNA gene sequences. Mol Phylogenet Evol. 2010;55(2): 738–743. pmid:20132899

59. Pieterse A, Tiedt LR, Malan AP, Ross JL. First record of Phasmarhabditis papillosa (Nematoda: Rhabditidae) in South Africa, and its virulence against the invasive slug, Deroceras panormitanum. Nematology. 2017;19(9): 1035–1050. https://doi.org/10.1163/15685411-00003105.

60. Ross JL, Pieterse A, Malan AP, Ivanova E. Phasmarhabditis safricana n. sp. (Nematoda: Rhabditidae), a parasite of the slug Deroceras reticulatum from South Africa. Zootaxa. 2018;4420(3): 391–404. pmid:30313534

61. Kiontke K, Barriere A, Kolotuev I, Podbilewicz B, Sommer R, Fitch DHA, et al. Trends, stasis, and drift in the evolution of nematode vulva development. Curr Biol. 2007;17(22): 1925–1937. pmid:18024125