Four Pristionchus species associated with two mass-occurring Parafontaria laminata populations

Millipedes were collected manually in Koumi near Matsubara Lake (Mat) and Nobeyama (Nob), Minamimaki Village, Nagano, Japan in 2016. These sites are relatively cool mountain areas (> 1,100 m a.s.l.) in central Japan.

Collected millipedes were brought to the laboratory, and kept at 10 to 15°C until dissection.

First, millipedes (33 individuals from Mat and 38 from Nob) were individually dissected on water agar (2.0% agar in a φ = 90 mm Petri dish) and kept at ca. 20°C for 1 week. After the propagation of Pristionchus spp. was confirmed, 3 to 5 individual nematodes were hand-picked from each plate, transferred individually to nematode digestion buffer (Kikuchi et al., 2009; Tanaka et al., 2012), and genotyped based on the D2-D3 LSU sequence, as described in Ye et al. (2007). Three species of Pristionchus were recognized: P. fukushimae from Mat, P. laevicollis from Nob, and P. degawai from both localities; all examined millipedes harbored one or two Pristionchus species. Pristionchus degawai dominated in both millipede populations (Table 1). Here, two strains of P. fukushimae, one strain of P. laevicollis, and one strain of P. degawai were established. The last two were used as the type strains of each species (Kanzaki et al., 2018).

To examine more detailed associations, the partial sequence of the mitochondrial cytochrome oxidase subunit I (mtCOI) gene was determined for the established cultures. Millipedes (20 per population) were dissected on the agar plate, and kept at ca. 20°C, as described above. The agar plates were examined under a dissecting microscope (S8 Apo, Leica) daily for 2 weeks.

When any Pristionchus species was recognized, 5 to 10 first-found gravid females were transferred separately into Nematode Growth Medium (φ = 40 mm Petri dish) previously inoculated with the Escherichia coli OP50 strain. The propagated nematodes were subcultured and kept as laboratory strains. Using this procedure, 215 temporal strains were established. An individual nematode was hand-picked from each culture and genotyped based on the D2-D3 LSU sequence (as described above). In addition, all temporal strains and two strains of P. fukushimae established in the first isolation were genotyped based on the mtCOI sequence according to the methods in Kanzaki and Futai (2002), and the sequences were analyzed phylogenetically with the online version of PhyML (http://www.atgc-montpellier.fr/phyml/) in which the analytical parameters were selected automatically (Guindon et al., 2010).

The 217 strains were separated into four species: P. entomophagus (one strain from Mat; one mtCOI haplotype), P. fukushimae (11 strains from Mat; three mtCOI haplotypes), P. laevicollis (three strains from Nob; three mtCOI haplotypes), and P. degawai (97 strains from Mat and 105 from Nob; 47 mtCOI haplotypes) (Fig. 2A and Table 2). The D2-D3 LSU sequences of these four species were identical to those deposited in GenBank (https://www.ncbi.nlm.nih.gov/genbank/?). The newly determined mtCOI sequences were deposited in GenBank with accession numbers LC589007-LC589060 (Table 3). In addition, an unidentified Oscheius sp. was isolated from all examined individual millipedes, although further study was not conducted for the species.

Figure 2:

Previous and present isolation records of Pristionchus spp. from millipedes are summarized in Figure 2. Pristionchus entomophagus is widely distributed in Europe (Kanzaki et al., 2014; Ragsdale et al., 2015), but this is the first report of the species in East Asia. Pristionchus fukushimae was originally described from stag beetle collected from northeast Japan (Fukushima Prefecture) (Ragsdale et al., 2013a, 2013b) and subsequently found from stag beetles and a soil sample in several relatively cool areas of Japan (Kanzaki, unpubl. obs.). Therefore, the species is considered widespread in the cool areas of Japan (and possibly other East Asian countries) mostly associated with decomposed plant materials and their related arthropods. Pristionchus laevicollis was previously isolated from Aegus laevicollis subnitidus in Nagoya (Kanzaki et al., 2018) and Kyoto (Kanzaki, unpubl. obs.), Japan, a relatively warm area, and might be distributed widely in Japan, although this remains unknown. So far, P. degawai has been isolated only from P. laminata (Kanzaki et al., 2018), and dominated the Pristionchus spp. associated with the millipede in this study.

Figure 3:

Considering millipedes’ bionomics, i.e., feeding on decomposing plant materials and inhabiting soil (Toyota et al., 2006), their habitat partially overlap stag beetles (Tanahashi et al., 2009), and it is not surprising that these hosts share the same nematode species. In addition, relatively wide host (carrier) range of Pristionchus spp. (P. fukushimae, P. laevicollis, and P. entomophagus) was confirmed. Pristionchus spp., except for fig associates, is not only phoretic, but also necromenic species. Thus, considering their host as a substrate, host specificity is not necessary to become strict. Contrastingly, many species inhabiting newly-dead wood, e.g., Micoletzkya (diplogastrid bacteria feeder) and Bursaphelenchus (aphelenchoidid fungal feeder), have relatively high host/carrier specificity (Susoy and Herrmann, 2014; Kanzaki and Giblin-Davis, 2018). This could be partially because of the stability of habitat, i.e., newly dead wood environment, especially inner bark, does not last long, and the nematodes need to be transferred to appropriate habitat by specialized carrier. While the soil/litter habitat are rather consistent, i.e., leaf litter is constantly supplied. Thus, the nematodes are not necessary to be transported long distance to specific habitat, but can utilize wide-ranged arthropod as both carrier and substrate.

Pristionchus degawai is highly divergent in the mtCOI sequence, separated into 47 haplotypes, with 24, 15, and 8 types isolated from Mat, Nob, and both localities, respectively. The average number of haplotypes isolated from an individual millipede was 5.0 (SD 1.5; range 3-8) for Mat and 3.0 (1.2; 1-5) for Nob (Table 2), suggesting that the Mat population is genetically more divergent than the Nob population.

The haplotypes were separated into four groups according to the bootstrap values (Fig. 2B). Groups 1 and 2 were mostly found from Nob, group 3 was found from Mat, and group 4 consisted of strains from both localities. Considering the phylogenetic relationships (Fig. 2B), the Mat population could be more derived. However, to clarify the population structure, more materials from different localities are necessary. Compared with the first isolation, in which the nematodes were examined 1 week after the dissection, the species and genotypes were more divergent in the extensive survey, where the first-found individuals were cultured. In addition, P. degawai, which was found from all millipedes in the extensive survey, was not found in several millipedes in the first isolation (Tables 1 and 2). This suggests that one species/genotype dominates quickly (during 1 week of culture), and the nematodes should be examined in earlier cultures to evaluate diversity.

The primary purpose of this study was to establish many Pristionchus cultures to find undescribed species; so the sample size for the haplotype analysis was limited. However, regardless of the limitation, the haplotype diversity of P. degawai, 47 types, was much higher than expected. So far, population genetic studies of Pristionchus species have been conducted only on P. pacificus (Herrmann et al., 2010; Morgan et al., 2012, 2014), and no other species has been examined. P. pacificus is a hermaphroditic and cosmopolitan species. Thus, as a comparative system, similar analysis of more locally distributed gonochoristic species could be valuable. Further analysis of P. degawai and other gonochoristic species focusing on haplotype diversity (without establishing strains) will give more detailed population genetics information.

To date, several Pristionchus spp. have been isolated from millipedes, including P. laevicollis, P. degawai, P. riukiariae (Kanzaki et al., 2018), P. arcanus (Kanzaki et al., 2012a, 2012b), P. pacificus (Kanzaki et al., 2016, 2018), and several other species (Kanzaki, unpubl. data). Pristionchus pacificus is a hermaphroditic species with a cosmopolitan distribution that dominates in the lowland temperate zone in Japan (Herrmann et al., 2007; Kanzaki et al., 2016; Kanzaki, unpubl. obs.). Nevertheless, none of our more than 200 tentative strains was P. pacificus, suggesting that Pristionchus spp. segregate by altitude, temperature, or both.

Although many thelastomatid and rhigonematid gut parasites have been described from millipedes, phoretic/necromenic species have not been examined systematically, i.e., in addition to above species descriptions, Carta et al., (2018) isolated Oscheius rugaoensis, Oscheius necromenus (Sudhaus and Schulte, 1989) and Mononchoides sp. from Chamberlinius hualienensis collected at Hachijojima Island, Japan and Kanzaki et al. (2016) reported P. pacificus, P. arcanus, and Oscheius spp. from Riukiaria spp. collected in three localities in Japan (Fig. 3). In the present study, although detailed analysis was not conducted, an unidentified Oscheius sp. was associated with all examined millipedes. Considering the isolation records provided by Kanzaki et al. (2016) and Carta et al. (2018), the genus seems commonly associated with millipedes. Oscheius is a typical soil dwelling bacteria feeder, and all member of the genus are hermaphroditic species (Sudhaus, 2011), and presumed to be a competitor of Pristionchus spp. for the food (substrate for food bacteria). Further analysis of its diversity and ecology, e.g., whether they are competing or segregating, will give new insight to understand the biological interaction among Pristionchus spp. and other millipede associates.

In addition, P. laminata has other populations that do not undergo typical mass-occurrence and has several congeners in Japan (e.g., Tanabe, 2002). Further surveys of phoretic nematodes of millipedes, especially in cool mountain areas, will reveal further nematode diversity, including the satellite model group, Pristionchus spp. as well as other phoretic/necromenic species.