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

Abstract Phoretic nematodes associated with two mass-occurring populations of the millipede Parafontaria laminata were examined, focusing on Pristionchus spp. The nematodes that propagated on dissected millipedes were genotyped using the D2-D3 expansion segments of the 28S ribosomal RNA gene. Four Pristionchus spp. were detected: P. degawai, P. laevicollis, P. fukushimae, and P. entomophagus. Of the four, P. degawai dominated and it was isolated from more than 90% of the millipedes examined. The haplotypes of partial sequences of mitochondrial cytochrome oxidase subunit I examined for Pristionchus spp. and P. degawai showed high haplotype diversity.

The genus Pristionchus (Kreis, 1932) is a satellite model system in many different fields of biology . The flagship species P. pacificus (Sommer et al., 1996) is used to study phenotypic plasticity, kin recognition, and chemical biology (e.g., Ragsdale et al., 2013aRagsdale et al., , 2013bLightfoot et al., 2019), and its congeners show high biological and physiological diversity (e.g., fig-associated species) and provide comparative information on genome-level diversification and speciation (e.g., Rödelsperger et al., 2014;Susoy et al., 2016). Therefore, dense taxon sampling of the genus has been conducted (e.g., Herrmann et al., 2006;Mayer et al., 2007;Kanzaki et al., 2012a), but the diversity of the genus is far from saturated . The genus is distributed widely, i.e., at least one species including undescribed ones have been isolated from all major continents Wang et al., 2015;Susoy et al., 2016), and in the previous studies, the genus has been isolated from relatively nutrient-rich substrates in European countries, partially because the taxonomy of Pristionchus has not been conducted in other areas of the world (e.g., Herrmann et al., 2015;Ragsdale et al., 2015). However, recently, their close association with wide-ranged invertebrates, mostly insects, has been recognized (e.g., Herrmann et al., 2006;Kanzaki et al., 2012aKanzaki et al., , 2012bKanzaki et al., , 2013, and the isolation from these hosts, especially Scarabaeidae , are increasing. In addition, recent surveys have found that many nominal and undescribed Pristionchus species are also associated with soil arthropods, such as millipedes Kanzaki unpubl. obs.).
Biologically, Pristionchus is phoretic and necromenic nematodes. The nematodes are isolated from their host/carrier insects as dauer ( = dormant and dispersal) stage, which can be reared on artificial media, suggesting phoretic association, but the worms also propagate on the carcass of their host/carrier insects utilizing the carcass as the substrate of their food bacteria (e.g., Rae et al., 2008;Cinkornpumin et al., 2014). Thus, the insects which can be utilized in multiple ways are very important for the nematodes, and the genus possibly provide an insight for analyzing evolution of feeding preference and host/carrier usage.
To isolate Pristionchus spp., phoretic/necromenic nematodes were examined in two mass-occurring populations of Parafontaria laminata in Nagano, Japan. The mass-occurring population of P. laminata is known as 'train millipede' covers the ground and even disturb the train service covering the railroad (Toyota et al., 2006), and it was previously considered an independent subspecies, 'P. l. armigera (Verhoeff)', but is now synonymized to the original subspecies (Tanabe, 2002). Parafontaria laminata is the only species undergoing mass-occurrence every 8 years in the collection area, and easily identified based on its mass-occurrence and general morphology (Fig. 1). Four species of Pristionchus, an Oscheius sp., and a hind gut parasite (Thelastomatidae) were isolated, and two Pristionchus spp. that were undescribed at the time were described taxonomically as P. degawai  and P. laevicollis  before this study .
The present study provides extensive information on the isolation of these four Pristionchus species.

Collection of millipedes
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.
Nematode isolation, culture, and genotyping 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 .
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).

Species diversity
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). Bootstrap values exceeding 50% are given on the appropriate clades. Some P. fukushimae strains did not amplify with the universal mtCOI primers (Kanzaki and Futai, 2002), probably because of a mutation in the primer region, and these strains were treated as a genotype (they do not appear in the tree  ; and X which was not amplified with universal primer set is P. fukushimae (Mat110). In addition to these strains, haplotypes A and B (P. fukushimae strains established in the first isolation) were included the phylogenetic analysis.
In addition, an unidentified Oscheius sp. was isolated from all examined individual millipedes, although further study was not conducted for the species. Previous and present isolation records of Pristionchus spp. from millipedes are summarized in Figure 2. Pristionchus entomophagus is widely distributed in Europe 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(Ragsdale et al., , 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  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 , and dominated the Pristionchus spp. associated with the millipede in this study.
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 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 habitat, but can utilize wide-ranged arthropod as both carrier and substrate.

Intraspecific variation in P. degawai
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.

Additional remarks
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 Morgan et al., 2012Morgan et al., , 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.
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.