Skip to main content

Symbiosis of the millipede parasitic nematodes Rhigonematoidea and Thelastomatoidea with evolutionary different origins

Abstract

Background

How various host–parasite combinations have been established is an important question in evolutionary biology. We have previously described two nematode species, Rhigonema naylae and Travassosinema claudiae, which are parasites of the xystodesmid millipede Parafontaria laminata in Aichi Prefecture, Japan. Rhigonema naylae belongs to the superfamily Rhigonematoidea, which exclusively consists of parasites of millipedes. T. claudiae belongs to the superfamily Thelastomatoidea, which includes a wide variety of species that parasitize many invertebrates. These nematodes were isolated together with a high prevalence; however, the phylogenetic, evolutionary, and ecological relationships between these two parasitic nematodes and between hosts and parasites are not well known.

Results

We collected nine species (11 isolates) of xystodesmid millipedes from seven locations in Japan, and found that all species were co-infected with the parasitic nematodes Rhigonematoidea spp. and Thelastomatoidea spp. We found that the infection prevalence and population densities of Rhigonematoidea spp. were higher than those of Thelastomatoidea spp. However, the population densities of Rhigonematoidea spp. were not negatively affected by co-infection with Thelastomatoidea spp., suggesting that these parasites are not competitive. We also found a positive correlation between the prevalence of parasitic nematodes and host body size. In Rhigonematoidea spp., combinations of parasitic nematode groups and host genera seem to be fixed, suggesting the evolution of a more specialized interaction between Rhigonematoidea spp. and their host. On the other hand, host preference of Thelastomatoidea spp. was not specific to any millipede species, indicating a non-intimate interaction between these parasites and their hosts.

Conclusions

The two nematode superfamilies, Rhigonematoidea and Thelastomatoidea, have phylogenetically distinct origins, and might have acquired xystodesmid millipede parasitism independently. Currently, the two nematodes co-parasitize millipedes without any clear negative impact on each other or the host millipedes. Our study provides an example of balanced complex symbioses among parasitic nematodes and between parasitic nematodes and host millipedes, which have been established over a long evolutionary history.

Background

How parasite diversity has evolved is an important question to understand host–parasite interactions, which form the basis of sustainable development. Nematoda is a large and ancient phylum that exhibits significant ecological variety. Free-living nematodes occupy almost all terrestrial, marine, and freshwater habitats, ranging from tropical to polar environments, with the number of individuals exceeding all other organisms [20, 38, 69]. Parasitic or phoretic nematodes live on or inside other organisms, and their diversity represents the fact that each host is associated with at least one unique nematode [5, 57]. Some parasitic nematodes live together with hosts as commensals [12, 24, 25, 54] or mutualistic partners [39], but others are harmful to the hosts [4, 7, 17, 56]. These four interactions—phoresis, commensalism, parasitism, and mutualism—are all recognized as subcategories of symbiosis [15], and are commonly observed in the phylum Nematoda. Their original symbiotic interactions with native hosts are sometimes disrupted by human activities, which can trigger devastating epidemics. Examples include pine wilt disease in Eastern Asia and Europe caused by the invasive pine-wood nematode Bursaphelenchus xylophilus [22], and crop diseases caused by a variety of plant parasitic/pathogenic nematodes in agricultural fields globally [61].

According to phylogenetic analysis, the phylum Nematoda is divided into five clades [5]. Like all living organisms, the Nematoda are believed to be of marine origin. The most primitive subclass is Enoplia (Clade II), most of which consists of free-living marine species. However, some nematode groups in the subclass Enoplia have evolved in terrestrial ecosystems and have developed plant and animal parasitism [6]. The next subclass to branch from Enoplia was Dorylaimia (Clade I), which contains free-living freshwater and terrestrial nematodes, and plant- and animal-parasitic nematodes. For example, known animal and plant pests, such as Trichinella and Dorylaimida, occur in the subclass Dorylaimia [6]. Another subclass that branched off, Chromadoria, is divided into three suborders: Rhabditina (Clade V), Tylenchina (Clade IV), and Spirurina (Clade III) [5, 6, 70]. Animal and plant parasitism have evolved independently in these five clades, and all nematodes belonging to the suborder Spirurina are animal parasites [5, 6, 70].

Since millipedes (class Diplopoda) are believed to be the first animals to inhabit terrestrial environments during the Silurian period (ca. 420 mya) [23, 41], it is a good model for understanding the evolutionary and ecological relationships between these parasitic nematodes and between hosts and parasites. We previously described the two nematode species, Rhigonema naylae (Rhigonematoidea) and Travassosinema claudiae (Thelastomatoidea), which are parasites of the xystodesmid millipede species Parafontaria laminata (Polydesmida: Xystodesmidae) in Aichi Prefecture, Japan [45, 46]. Both nematode species are in the suborder Spirurina, but in different infraorders. Rhigonema naylae belongs to the infraorder Rhigonematomorpha, which includes the millipede parasites Ransomnematoidea [40, 42] and Rhigonematidae [32, 45, 46], and the amphibian and reptile parasite Cosmocercoidea [8, 9]. Travassosinema claudiae is a member of the Oxyuridomorpha, which includes a wide variety of parasites, such as vertebrate-parasitic Oxyuroidea [49], invertebrate-parasite Coronostomatoidea [55], and Thelastomatoidea [53, 47, 48]. The family Thelastomatidae mainly includes cockroach parasites and is a highly diversified group characterized by low host specificity [2, 35, 54].

Here, we reveal that the two parasitic nematodes, Rhigonematoidea spp. and Thelastomatoidea spp., have different evolutionary origins but co-infect the xystodesmid millipedes without clear negative impacts between each parasitic nematode and host millipedes.

Results

Two parasitic nematodes R. naylae and T. claudiae isolated from the two millipedes

Parafontaria laminata CU and Parafontaria tonominea species complex CU

We isolated two parasitic nematodes, R. naylae and T. claudiae, from the xystodesmid millipede, P. laminata, in the Chubu University campus in Aichi Prefecture, Japan (Table 1, Fig. 1). Hereafter, this millipede species is referred to as P. laminata CU. Rhigonema naylae has a large body size, short pharynx, round head, and short tail. The body size of T. claudiae is smaller than that of R. naylae; the females have a typical umbrella-like head and long tail, and the males have a tiny body with a spicule (Fig. 2) [45, 46]. No clear sexual dimorphism was detected in either nematode species during the juvenile stage; however, it was possible to distinguish the two species at any developmental stage.

Table 1 Summary of collection data of millipedes and parasitic nematodes
Fig. 1
figure 1

Geographical location of millipede collection sites. The zoomed maps A and B from the red squared on the map of Japan are created using Google Earth. Location numbers are, 1 Chubu University, Kasugai City, Aichi Prefecture. 2 Mt. Kinka, Gifu City, Gifu Prefecture. 3 Hyakunen Park, Seki City, Gifu Prefecture. 4 Embara, Yamagata City, Gifu Prefecture. 5 Mt. Shimono, Yamaga City, Kumamoto Prefecture. 6 Miyanoura, Kagoshima City, Kagoshima Prefecture. 7 Shiroyama, Kagoshima City, Kagoshima Prefecture

Fig. 2
figure 2

Two parasitic nematodes, Rhigonema naylae and Travassosinema claudiae parasites of the xystodesmid millipede Parafontaria laminata CU and P. tonominea species complex CU. Nomarski differential interference construct images of adult female and male. Body sizes (Average ± S.D.) are adapted from [45, 46]

A total of 113 P. laminata CU millipedes were captured from April 2018 to July 2019, and parasitic nematode infections were examined. The total infection prevalence (counted as present when at least one individual was detected from any developmental stage) of R. naylae and T. claudiae was 31.0% and 27.4%, respectively, and co-infection prevalence (infected by the two nematode species together) was 10.6% (Table 2). The density (number of males, females, and juveniles of nematodes in individual hosts) and infection percentage at each stage are also presented in Table 2. We found 11.5% of males and 12.4% of females of R. naylae, and 1.8% of males and 3.5% of females of T. claudiae (Table 2). All adult male and female nematodes appeared healthy, and all females matured, and many eggs were observed in the uteri.

Table 2 Population of the two parasitic nematodes in the two Parafontaria millipedes

Another xystodesmid millipede in the same area was also infected with the two parasitic nematodes. From the body shape and male genitalia, this millipede was determined to be a member of the P. tonominea species complex (Additional file 1: Figure S1) [67]. Hereafter, this millipede species is referred to as P. tonominea species complex CU. A total of 82 P. tonominea species complex CU millipedes were captured from May 2018 to July 2019, and parasitic nematode infections were examined. The infection prevalence of both parasitic nematodes in the P. tonominea species complex CU was higher than that in P. laminata. The total infection prevalence (counted as present when at least one individual was detected from any developmental stage) of R. naylae and T. claudiae was 96% and 72%, respectively, and the co-infection prevalence was 72% (Table 2). A total of 89% of males and 87% of females of R. naylae, and 52% of males and 17% of females of T. claudiae were detected (Table 2). All females of both nematode species in the P. tonominea species complex CU were mature, and many eggs were present in the uteri. From these results, we deduced that both parasitic nematodes, R. naylae and T. claudiae, were able to use the two millipedes as hosts; however, the infection prevalence and density of T. claudiae were lower than those of R. naylae.

Population of R. naylae was not negatively affected by co-infection with T. claudiae

We analyzed whether the two parasitic nematodes, R. naylae and T. claudiae, affected each other when they co-infected the host millipede P. tonominea species complex CU. There were three infection patterns: (1) co-infection with the two parasitic nematodes (N = 59), (2) infection with only R. naylae (N = 20), and (3) no nematode infection (N = 3). We did not observe any millipede host P. tonominea species complex CU infected with T. claudiae alone. We studied the infection prevalence % and mean density (mean number of nematodes per infected host) of male, female, and juvenile R. naylae in: (1) all the hosts (N = 82), (2) hosts infected with only R. naylae (N = 20), and (3) hosts co-infected with R. naylae and T. claudiae (N = 59). The infection prevalence of R. naylae juveniles in hosts co-infected with T. claudiae was significantly higher than that in hosts infected with only R. naylae, but there were no differences in the prevalence of R. naylae adults (Fig. 3A). In addition, the mean densities of R. naylae females and juveniles were significantly higher when co-infected with T. claudiae (Fig. 3A). The effects of host body size on these infection conditions was not significantly different (Additional file 2: Figure S2). Thus, co-infection with T. claudiae did not affect the population of R. naylae in the alimentary tract of the P. tonominea species complex CU.

Fig. 3
figure 3

Infection prevalence (%) (orange dot, right Y-axis) and mean density (mean number of nematodes per infected host) (blue dot, left Y-axis) of male, female, and juvenile of the parasitic nematodes in Parafontaria tonominea species complex CU. A Infection prevalence and density of R. naylae in all the hosts (N = 82), hosts infected with only Rhigonema naylae (N = 20), and the hosts co-infected with R. naylae and Travassosinema claudiae (N = 59). B Prevalence and density of R. naylae in all the hosts collected during spring (N = 33) and summer (N = 49). C Infection prevalence and density of T. claudiae in all hosts collected during spring (N = 33) and summer (N = 49). Error bars indicate confidence interval with 95% confidence limit. *p < 0.05, **p < 0.005, Fisher’s exact test for comparing prevalence, and Bootstrap 2-sample t-test for comparing mean densities

We also checked the effects of seasonal differences (spring vs. summer) on the prevalence and densities of the two parasitic nematodes in the host millipede P. tonominea species complex CU. The infection prevalence of R. naylae at all stages was high, and there was no significant difference between spring (91%, N = 33) and summer (96%, N = 49) (Fig. 3B). The density of R. naylae juveniles was significantly higher during summer (Fig. 3B). As indicated previously, the infection prevalence of T. claudiae was lower than that of R. naylae; however, female prevalence in summer was higher (86%, N = 42), and those of the females and juveniles of T. claudiae significantly increased during summer (Fig. 3C). The density of T. claudiae females was also higher during summer (Fig. 3C). The effects of host body size on these seasonal conditions were not significantly different (Additional file 2: Figure S2). Thus, while the density and prevalence of R. naylae remained high year-round, those of T. claudiae were strongly affected by seasonal changes.

Two nematodes commonly co-infected in Parafontaria millipedes

We also collected five Parafontaria species from three other locations in Gifu Prefecture (Fig. 1, Table 1). From morphological observations, one was P. laminata Kinka (collected from Mt Kinka, Gifu City), and belonged to the same species as P. laminata CU; another was P. longa Embara (collected from Embara, Yamagata City); and the three other millipedes were different species that were tentatively classified as members of the P. tonominea species complex (Additional file 1: Figure S1). Similar to the P. tonominea species complex CU and P. laminata CU, these five Parafontaria millipedes tended to be infected with the two parasitic nematodes, R. naylae and T. claudiae. All D2D3 LSU sequences within the species were almost identical (Additional files 3, 4: Tables S1 and S2). Infection prevalence and mean density of parasitic nematode males, females, and juveniles are shown in Fig. 4 and Additional file 5: Table S3. The characteristics of the two nematodes in the Parafontaria species were similar, and the infection prevalence of R. naylae was higher than that of T. claudiae; however, all adult male and female nematodes matured and reproduced.

Fig. 4
figure 4

Infection prevalence (%) of the two parasitic nematodes males (M), females (F), and juveniles (J) in each Parafontaria millipede. Details of the host information are in Table 1

In addition to the two parasitic nematodes, unknown species were isolated from the four host species. Two individuals of P. tonominea species complex Kinka (2/43, 4.5%), two of P. laminata Kinka (2/35, 5.7%), eight of P. tonominea species complex Embara (8/20, 40%), and one of seven P. longa Embara (1/7, 14%) were infected. These nematodes were morphologically distinguishable from R. naylae and T. claudiae. D2D3 LSU sequencing data from the nematode isolated from the P. tonominea species complex Embara showed that it had a similar sequence to Cephalobellus brevicaudatus [14], MF668725.1 with 92% identity (671/727 bp). Combined with morphological observations, this nematode is an undescribed species, and we temporarily named it Cephalobellus sp.1 (Table 1). We could not identify the other unknown nematodes because of the small number of specimens.

Parasitic nematodes in Riukiaria millipedes

To determine the infection patterns of parasitic nematodes in other xystodesmid millipedes, we collected Riukiaria spp. from Kyushu Island. Three species (four isolates) were collected from three locations in Kyushu (Table 1, Additional file 1: Figure S1). All of the species tended to be infected with the three parasitic nematodes, Rhigonematoidea sp. 1, T. claudiae, and Thelastomatidae spp. (Table 1). From the D2D3 LSU sequencing data, Thelastomatidae sp. 1 and sp. 2 and Rhigonematoidea sp.1 were undescribed species (Additional files 3, 6: Tables S1 and S4). Surprisingly, T. claudiae, the species that was isolated from Parafontaria millipedes, was also found in all Riukiaria millipedes in our samples (Table 1, Additional file 4: Table S2). Thelastomatidae sp. 1 and Thelastomatidae sp. 2 were combined and referred to as Thelastomatidae spp.

Similar to R. naylae in Parafontaria millipedes, the infection prevalence of male, female, and juvenile Rhigonematoidea sp.1 was high (67–100%), whereas the infection prevalence of T. claudiae and Thelastomatidae spp. was lower (0–87%) (Additional file 7: Table S5). In addition, males and juveniles of T. claudiae and Thelastomatidae spp. in Riukiaria hosts were detected only in a few cases. Infection data of only nematode females were selected and are shown in Fig. 5. Since all of the females were mature, and many eggs in their uteri were visible, we concluded that all of these parasitic nematodes used Riukiaria millipedes as hosts and could reproduce.

Fig. 5
figure 5

Infection prevalence (%) (orange dot, right Y-axis) and mean density (the mean number of nematodes per infected host) (blue dot, left Y-axis) of the females of the parasitic nematodes. Error bars indicate confidence intervals with 95% confidence limit

We then analyzed whether the two parasitic nematodes, T. claudiae and Thelastomatidae spp., affect each other when co-infecting host millipedes. The mean densities of T. claudiae adult females in all Riukiaria spp. were combined and calculated for two infection patterns: (1) infected with only T. claudiae, and (2) co-infected with T. claudiae and Thelastomatidae spp. In addition, the mean densities of Thelastomatidae spp. adult females in Riukiaria spp. were calculated for two infection patterns: (3) infected with only Thelastomatidae spp., and (4) co-infected with T. claudiae and Thelastomatidae spp. (Fig. 6). Since 99% of the Riukiaria millipedes (N = 97 of 98) were infected with Rhigonematoidea sp.1, we could not examine the effect of Rhigonematoidea sp.1 on the densities of Thelastomatidae spp. and T. claudiae. The densities of both parasitic nematodes in co-infected hosts were significantly higher than those in single infections (Fig. 6).

Fig. 6
figure 6

Mean density of Travassosinema claudiae adult female in all Riukiaria spp. (1) infected with only T. claudiae and (2) co-infected with Thelastomatidae spp., and mean density of Thelastomatidae spp. adult female in all Riukiaria spp. (3) infected with only Thelastomatidae spp. and (4) co-infected with T. claudiae. Error bars indicate confidence interval with 95% confidence limit. *p < 0.05, Bootstrap 2-sample t-test for comparing mean densities

Furthermore, we found that the mean densities of Rhigonematoidea sp.1 adult females were significantly higher when the two parasitic nematodes Thelastomatidae spp. and T. claudiae were co-infected in the host at the same time, compared with those with only Thelastomatidae spp. or T. claudiae (Fig. 7A). Since the mean body size of Riukiaria spp. was larger when co-infected with the two parasitic nematodes Thelastomatidae spp. and T. claudiae (Fig. 7B), the prevalence of these two parasites was positively correlated with host body size. Thus, the number of parasites in millipede hosts increased with increasing host size, suggesting that these nematodes did not compete and eliminate each other.

Fig. 7
figure 7

A Mean density of Rhigonematoidea sp. 1 adult female, in Riukiaria spp. infected with only Travassosinema claudiae, only Thelastomatidae spp., and both parasites. Error bars indicate confidence intervals with 95% confidence limit. *p < 0.05, Bootstrap 2-sample t-test for comparing mean densities. B Average ± S.D. of body length of host millipede Riukiaria spp., infected with only T. claudiae, only Thelastomatidae spp., or both parasites. **p < 0.005, Statistical differences were analyzed by Tukey’s multiple comparison test followed by Bonferroni correction

Parasitic nematodes in xystodesmid millipedes with phylogenetically different origins

We used D2D3 LSU rDNA genes selected from representatives of the four infraorders Oxyuridomorpha, Rhigonematomorpha, Gnathostomatomorpha, Spirurinomorpha, and the two families Dracunculoidea and Camallanoidea (Additional file 8: Table S6). A phylogenetic tree covering the suborder Spirurina was constructed that supported the taxonomic relationship previously reported [37, 43], these infraorders were clearly separated and clustered (Fig. 8A). Rhigonematoidea was clustered with Cosmocercoidea and Ransomnematoidea. As Travassosinematidae appeared to be close to Thelastomatidae, and both are classified in Thelastomatoidea (Fig. 8A), the xystodesmid millipedes examined in this study tended to be infected with parasitic nematodes belonging to these two infraorders.

Fig. 8
figure 8figure 8figure 8

Phylogenetic trees of parasitic nematodes inferred from the D2D3 LSU rRNA. Bootstrap values with 1000 replicates are shown next to the branches. A Phylogenetic tree of the suborder Spirurina, constructed using the neighbor-joining (NJ) method. B Phylogenetic tree of Rhigonematomorpha, constructed using the maximum likelihood (ML) method. Newly obtained sequences are in red. *The position of Rhigonematoidea sp.1. is emphasized. C Phylogenetic tree of the Thelastomatoidea constructed using the ML method. Newly obtained sequences are in red, and their hosts are in blue. **The position of Cephalobellus sp.1 is emphasized

In the maximum likelihood (ML) tree, Rhigonematoidea sp. 1 was clearly separated from the Rhigonematidae cluster but could be within Rhigonematoidea (Fig. 8B). All parasitic nematodes in Rhigonematomorpha are believed to be millipede parasites, but Cosmocercoidea (reptile and amphibian parasites) were included in Rhigonematomorpha in our phylogenetic analysis (Fig. 8B). The parasitic nematodes T. claudiae, Cephalobellus sp. 1, Thelastomatidae sp. 1, and Thelastomatidae sp. 2 in our experiments were positioned in Thelastomatoidea (Fig. 8C). The closest group to Cephalobellus sp. 1 was Travassosinematidae, and it was clearly separated from Thelastomatidae. Both Thelastomatidae sp. 1 and sp. 2 were positioned close to the cockroach-parasitic nematode Thelastoma bulhoesi (Fig. 8C).

Discussion

Two nematode species, R. naylae and Rhigonematoidea sp. 1 in the superfamily Rhigonematoidea, were isolated from millipede hosts Parafontaria spp. and Riukiaria spp., respectively, with high infection prevalence (Table 1). Phylogenetic analysis with rRNA SSU, rRNA LSU, and mitochondrial DNA indicated that Rhigonematomorpha is not monophyletic, but is nested within Ascaridomorpha [36, 43]. The superfamily Cosmocercoidea was also positioned within the Rhigonematomorpha in our analysis (Fig. 8B). This superfamily contains reptile and amphibian parasitic nematodes, and is classified within Ascaridomorpha [8, 9, 50]. Thus, further studies are required on the taxonomy of these groups; however, these infraorders are clearly closely related phylogenetically. In addition, the infraorder Rhigonematomorpha is composed of Rhigonematoidea and Ransomnematoidea and contains exclusively millipede parasites [19, 33]. We also showed that the combination of parasitic nematode groups and host genera seems to be specifically fixed. For example, Riukiaria hosts were exclusively infected with Rhigonematoidea sp. 1, whereas Parafontaria hosts were exclusively infected with R. naylae (Fig. 8B). Combining these data, we predicted the appearance of Rhigonematoidea, which might have divided from a common ancestor with Ascaridomorpha and acquired millipede parasitism at an earlier period before millipedes began to inhabit terrestrial environments.

The infection cycle of Rhigonematoidea is not well known; however, nematode eggs laid by adult females are probably deposited within their host feces, similar to parasitic nematodes in Oxyuridomorpha. When eggs were collected from the adult females of R. naylae and incubated in phosphate buffered saline (PBS) buffer, they began to develop, and larvae hatched (data not shown). These parasites might have a special mechanism that specifically infects only millipedes. For example, larvae are only attracted to specific hosts. In addition, autoinfection might occur similar to that in human parasitic nematode Strongyloides stercoralis [73], resulting in repeated generations occurring in the same host individual, leading to high infection prevalence throughout the season (Fig. 3B).

Four nematode species in Thelastomatoidea, T. claudiae, Cephalobellus sp. 1, Thelastomatidae sp. 1, and Thelastomatidae sp. 2 were isolated from xystodesmid millipedes (Table 1). Cephalobellus sp. 1 was only found in the P. tonominea species complex Embara, and Thelastomatidae sp. 1 and Thelastomatidae sp. 2 were only found in Riukiaria spp. Interestingly, R. naylae infected all of the millipedes. Cephalobellus is a member of the family Thelastomatidae and is found in many invertebrate hosts, including mole crickets, crane flies, white grubs, cockroaches, and millipedes [14, 29, 60]. In addition, we isolated nematode genera from the mole cricket Gryllotalpa orientalis, white grub Protaetia orientalis submarumorea, and pill millipedes Hyleoglomeris japonica (data not shown). Travassosinema is currently one of three genera (together with Indiana and Pulchrocephala) in the family Travassosinematidae [1, 28]. Although reported mainly in millipedes, it has also been found in other invertebrate hosts, including the larvae of a scarabaeid beetle, wood-burrowing cockroach, and earthworms [1, 30, 31, 34, 46, 58, 63]. Because almost all studies thus far have been nematode isolation records and taxonomical descriptions, the host–parasite relationship for each species has not yet been clarified, but the broad host range is a clear characteristic of these genera.

The family Thelastomatidae has also been isolated from many invertebrates; however, overwhelmingly, a large number of species has been found in Blattodea hosts [2, 52, 53, 64]. From our field surveys and previous studies, the host specificity of nematodes belonging to the family Thelastomatidae appeared high [13, 47, 48, 52, 53, 64]. Yet, the same nematode species have been isolated from different cockroach species [35, 54, 62]. We showed, experimentally, that Leidynema appendiculatum, a thelastomatid parasitic nematode, was capable of infecting five cockroach species from three families in two suborders [54]. While ecological interactions between host and parasitic nematodes determine the host range, a broad host range might still be possible in Thelastomatoidea.

Parasitic nematodes in Thelastomatoidea are thought to share a simple infection cycle; nematode eggs are laid by adult females, deposited within their host’s feces, and then released outside the host’s body; ingestion of the eggs by new host individuals leads to infection [3, 18, 53]. When eggs were collected from adult females of T. claudiae and incubated in a PBS buffer, they developed until the infection stage but never hatched (data not shown). Such characteristics of embryogenesis are shared with the cockroach-parasitic nematode family Thelastomatidae. Because it is necessary for the eggs of T. claudiae to exit the host body once to generate the next developmental stage, infection prevalence might have been lower (Figs. 3A, 4, 5) and affected by the season (Fig. 3C) when compared with those of Rhigonematoidea. In addition, T. claudiae has a population structure similar to that of the cockroach-parasitic nematodes, with one adult male and a few adult females being present in the host [52, 54, 74]. These properties are also frequently observed in Thelastomatoidea.

Since infection of the two parasitic nematodes Rhigonematoidea and Thelastomatoidea in xystodesmid millipedes in Japan is universal, as far as we have investigated, the relationship between host millipedes and parasitic nematodes appears not to be pathogenic; rather, it is commensalism (obligatory for nematode, host is not affected) or parasitism (obligatory for nematode, host is inhibited; [15]. Only a few studies have shown the effects of these parasitic nematodes on their invertebrate hosts (e.g., [71, 72]), and they have generally been found to be harmless [3, 52, 54]. The composition of cockroach gut microbiomes is influenced by parasitic nematode species [71, 72], suggesting direct or indirect effects of nematodes on the host. As the parasitic nematodes in our study are all obligatory parasites, and culturing methods of their host millipedes have not yet been established, laboratory experiments controlling infection conditions are not possible. In addition, Thelastomatoidea was not eliminated by Rhigonematoidea when co-infected (Figs. 3, 7). From our ecological, parasitological, and phylogenetic studies, we hypothesized that some ancestors of cockroach-parasitic nematodes might have changed their host to the Passalidae beetle at an early stage, and become the family Hystrignathidae. Similarly, another ancestor might have changed its host to the millipede and become the family Travassosinematidae and Cephalobellus (Fig. 8C). Their properties might allow Thelastomatoidea nematodes to use the millipede as a host and also to tolerate co-infection with Rhigonematoidea. The exchange of parasitic nematode between hosts living in the same ecological niches could occur (for example, from cockroaches to millipedes), and new parasitic nematodes may have evolved in millipedes after reproductive isolation. This parasitic nematode might also have been isolated from other invertebrate hosts with similar ecology. Speciation of nematodes could also result from genetic isolation, which is largely influenced by their associated animals and plants [44].

Interestingly, all of the millipedes studied in this work tended to be infested with the same species, T. claudiae. These millipedes were mainly sampled from forests where the planted Japanese cedar Cryptomeria japonica was dominant (Table 1). This species might have spread as a result of artificial planting and, because of its low host specificity, T. claudiae could have established relationships with local millipedes. This hypothesis could be clarified using molecular markers, such as the whole mitochondrial genome or ITS regions, which reflect intraspecific polymorphisms.

Conclusions

In this study, we found parasitic nematodes of the two superfamilies Rhigonematoidea and Thelastomatoidea, commonly co-infecting xystodesmid millipedes. Both superfamilies were in the suborder Spirurina, but in apparently different infraorders. We found that the infection prevalence and densitiy of Thelastomatoidea were lower than those of Rhigonematoidea, which reflects the difference in infection mechanism. However, the two nematode superfamilies Rhigonematoidea and Thelastomatoidea were not in a competitive relationship, and co-infected all millipede hosts. Our field study shows an example of the complex symbiosis among parasitic nematodes and between hosts and parasites, established throughout evolution.

Methods

Field collection of millipedes

Adult xystodesmid millipedes were collected from seven locations in Japan (Table 1). All millipedes were manually collected and maintained at low temperatures until dissection. Before dissection, all millipedes were first confirmed to be alive, and then the species and sexes were identified following the relevant illustration references [65, 67, 68].

Millipede dissection and parasitic nematode observation

After body length and sex (male or female) were recorded, millipedes were dissected to determine the presence or absence of parasitic nematodes in the alimentary tract. However, once the host millipede died, the parasitic nematode also died and degraded immediately, which made their identification difficult. All the millipedes used in this experiment were alive when dissected, and the parasitic nematodes were active after being released from the host gut. Nematode species, developmental stages (adult or larvae), sex (male or female), and numbers were recorded (Additional file 9: Table S7). The prevalence % of infected millipedes among all the millipedes examined, mean density (the mean number of nematodes per infected host), and bootstrap confidence interval of each nematode species and stage (adult male, adult female, and larva) were calculated using Quantitative Parasitology 3.0 [10, 59].

Morphological characteristics were captured using Nomarski DIC optics (Eclipse E600, Nikon, Japan) equipped with a CCD camera (VTCH1.4ICE; Visualix, Japan) as described by Ozawa and Hasegawa [54]. Molecular characterization was also performed as described by Ozawa and Hasegawa [54]. Briefly, genomic DNA was extracted from individual female nematode, and the D2D3 expansion segment of the 28S ribosomal RNA gene (D2D3 LSU) was amplified using the universal primers D2A (5-ACA AGT ACC GTG AGG GAA AGT TG-3) and D3B (5-TCG GAA GGA ACC AGC TAC TA-3) [51]. Samples were submitted to Eurofins Genomics (Tokyo, Japan) for sequencing from both strands, using the same PCR primers. Sequences were deposited in GenBank NCBI (http://www.ncbi.nlm.nih.gov/genbank/).

Pairwise comparison of percent differences (D) within each nematode group (Rhigonematoidea spp., 730 bp; Travassosinematidae spp., 704 bp; Thelastomatidae spp., 676 bp) were performed using the formula D = (M/L) × 100 [16], where M is the number of alignment positions at which the two sequences have a base in common, and L is the total number of alignment positions.

Population analysis of the two parasitic nematodes in P. tonominea species complex CU

To study the effect of co-infection with R. naylae and T. claudiae on the prevalence and densities of the host, P. tonominea species complex CU was classified into three infection patterns: (1) all the hosts, (2) the hosts only infected with R. naylae, or (3) the hosts co-infected with R. naylae and T. claudiae. As no host infected with only T. claudiae was observed in this experiment, we compared the prevalence and density of R. naylae in hosts infected with R. naylae and those co-infected with R. naylae and T. claudiae.

In addition, to study the effects of seasonal differences on the prevalence and density of the two parasitic nematodes in the P. tonominea species complex CU, hosts were classified into two seasonal patterns, collected during spring (from March to June) or summer (July–October) and the prevalence and density of R. naylae and T. claudiae between the two seasons were compared.

Significant differences in nematode infection prevalence and mean densities were calculated using Fisher’s exact test and the bootstrap 2-sample t-test, respectively [59]. Differences in host body size in these infection conditions or seasonal conditions were analyzed using Tukey’s multiple comparison test followed by Bonferroni correction.

Population analysis of the three parasitic nematodes in Riukiaria spp.

To study the effect of co-infection with the three parasitic nematodes, Rhigonematoidea sp. 1, T. claudiae, and Thelastomatidae spp. on the prevalence and density, all Riukiaria hosts were combined and categorized into the three infection patterns: (1) the host infected with only T. claudiae, (2) the host co-infected with T. claudiae and Thelastomatidae spp., or (3) the host infected with only Thelastomatidae spp. Instead of the high infection prevalence of all stages and sexes of Rhigonematoidea sp. 1, the prevalence of adult males and juveniles of T. claudiae and Thelastomatidae spp. was low, and we used data only from adult females. We compared the densities of T. claudiae and Thelastomatidae spp. in hosts infected with a single or two parasite species.

Furthermore, the densities of Rhigonematoidea sp. 1 females were compared between the three host conditions: (1) hosts infected with only T. claudiae, (2) hosts co-infected with T. claudiae and Thelastomatidae spp., or (3) hosts infected with only Thelastomatidae spp. Significant differences in nematode mean densities were calculated using a bootstrap 2-sample t-test [59]. Differences in host body size under these infection conditions were analyzed using Tukey’s multiple comparison test followed by Bonferroni correction. Differences in host body size under these infection conditions were analyzed using Tukey’s multiple comparison test followed by Bonferroni correction.

Phylogenetic analysis

For the phylogenetic analysis, D2D3 LSU sequences obtained in this experiment (Table 1) and the data uploaded to the NCBI database and published in the paper were used (Additional file 8: Table S6). ClustalW multiple alignment was conducted using BioEdit version 7.2.6 [26] and sequence alignments were trimmed automatically by trimAI with default settings [11]. A phylogenetic tree of the suborder Spirurina was constructed from evolutionary distances by the neighbor-joining method using Mega 6.0 software and the Tamura-Nei model [66]. Phylogenetic trees of the Rhigonematomorpha and Thelastomatoidea were constructed from evolutionary distances by the maximum likelihood (ML) method using Mega 6.0 software [66] and the Hasegawa-Kishino-Yano model [27]. Phylogenetic robustness was inferred by bootstrap analysis using 1000 iterations [21].

Availability of data and materials

Sequence data produced in this study are available in the NCBI (https://www.ncbi.nlm.nih.gov/) with Accession No.: MT988313-MT988378.

Abbreviations

CU:

Chubu University as a collection site

D2D3 LSU:

D2D3 expansion segment of the large subunit ribosomal RNA gene

SSU:

Small subunit ribosomal RNA gene

References

  1. Adamson ML. Oxyuridan (Nematoda) parasites of Scaphiostreptus seychellarum with comments on the families Pulchrocephalidae Kloss, 1959 and Travassosinematidae Rao, 1958. Can J Zool. 1987;65:2747–54. https://doi.org/10.1139/z87-416.

    Article  Google Scholar 

  2. Adamson ML, van Waerebeke D. Revision of the Thelastomatidae, Oxyurida of invertebrate hosts I. Thelastomatidae. Syst Parasitol. 1992;21:21–63.

    Article  Google Scholar 

  3. Adamson ML. Evolutionary patterns in life histories of Oxyurida. Int J Parasitol. 1994;24:1167–77.

    Article  CAS  Google Scholar 

  4. Baniya A, Huguet-Tapia JC, DiGennaro P. A draft genome of Steinernema diaprepesi. J Nematol. 2020;52:1–4. https://doi.org/10.21307/jofnem-2020-069.

    Article  CAS  PubMed  Google Scholar 

  5. Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, Vanfleteren JR, Mackey LY, Dorris M, Frisse LM, Vida JT, Thomas WK. A molecular evolutionary framework for the phylum Nematoda. Nature. 1998;392:71–5. https://doi.org/10.1038/32160.

    Article  CAS  PubMed  Google Scholar 

  6. Blaxter ML, Koutsovoulos G. The evolution of parasitism in Nematoda. Parasitol. 2015;142:S26-39. https://doi.org/10.1017/S0031182014000791.

    Article  Google Scholar 

  7. Bungiro R, Cappello M. Hookworm infection: new developments and prospects for control. Curr Opin Infect Dis. 2004;17:421–6. https://doi.org/10.1097/00001432-200410000-00006.

    Article  PubMed  Google Scholar 

  8. Bursey CR, Goldberg SR. Helminths in Mesaspis monticola (Squamata: Anguidae) from Costa Rica, with the description of a new species of Entomelas (Nematoda: Rhabdiasidae) and a new species of Skrjabinodon (Nematoda: Pharyngodonidae). Parasite. 2006;13:183–91. https://doi.org/10.1051/parasite/2006133183.

    Article  CAS  PubMed  Google Scholar 

  9. Bursey CR, Goldberg SR, Telford SR. Gastrointestinal helminths of 14 species of lizards from Panama with descriptions of five new species. Comp Parasitol. 2007;74:108–40. https://doi.org/10.1654/4228.1.

    Article  Google Scholar 

  10. Bush AO, Lafferty KD, Lotz JM, Shostak AW. Parasitology meets ecology on its own terms: Margolis et al. revisited. J Parasitol. 1997;83:575–83.

    Article  CAS  Google Scholar 

  11. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAI: a tool for automated alignment trimming in large-scale phylogenetic analysis. Bioinformatics. 2009;25:1972–3.

    Article  Google Scholar 

  12. Carreno RA. The systematics and evolution of pinworms (Nematoda: Oxyurida: Thelastomatoidea) from invertebrates. J Parasitol. 2014;100:553–60. https://doi.org/10.1645/14-529.1.

    Article  PubMed  Google Scholar 

  13. Carreno RA. New species of Hammerschmidtiella Chitwood, 1932, and Blattophila Cobb, 1920, and new geographical records for Severianoia annamensis Van Luc & Spiridonov, 1993 (Nematoda: Oxyurida: Thelastomatoidea) from Cockroaches (Insecta: Blattaria) in Ohio and Florida, U.S.A. Zootaxa. 2017;4426:429–41. https://doi.org/10.11646/zootaxa.4226.3.6.

    Article  Google Scholar 

  14. Carreno RA, Kiebler L, Tuhela L. First Record of Cephalobellus brevicaudatus (Leidy, 1851) Christie, 1933 (Nematoda: Oxyurida: Thelastomatoidea), from Cranefly Larvae (Diptera: Tipulidae) in Ohio, USA. Comp Parasitol. 2018;85:133–40. https://doi.org/10.1654/1525-2647-85.2.133.

    Article  Google Scholar 

  15. Chen TC. Is parasitism symbiosis? A definition of terms and the evolution of concepts. In: Toft CA, Aeschlimann A, Bolis L, editors. Parasite-host associations. Coexistence or conflict? Oxford: Oxford University Press; 1991. p. 15–36.

    Google Scholar 

  16. Chilton NB, Gasser RB, Beveridge I. Differences in a ribosomal DNA sequence of morphologically indistinguishable species within the Hypodontus macropi complex (Nematoda: Strongyloidea). Int J Parasitol. 1995;25:647–51.

    Article  CAS  Google Scholar 

  17. Ciche T. The biology and genome of Heterorhabditis bacteriophora, WormBook, ed. The C. elegans Research Community, WormBook. 2007.

  18. Cook GC. Enterobius vermicularis infection. Gut. 1994;35:1159–62. https://doi.org/10.1136/gut.35.9.1159.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. De Ley P, Blaxter ML. Systematic position and phylogeny. In: Lee DL, editor. The biology of nematodes. London: Taylor and Francis; 2002. p. 1–30.

    Chapter  Google Scholar 

  20. De Ley P. A quick tour of nematode diversity and the backbone of nematode phylogeny (January 25, 2006), WormBook, ed. The C. elegans Research Community, WormBook. 2006.

  21. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.

    Article  Google Scholar 

  22. Futai K. Pine wood nematode, Bursaphelenchus xylophilus. Annu Rev Phytopathol. 2013;51:61–83. https://doi.org/10.1146/annurev-phyto-081211-172910.

    Article  CAS  PubMed  Google Scholar 

  23. Garwood RJ, Edgecombe GD. Early terrestrial animals, evolution, and uncertainty. Evol Edu Outreach. 2011;4:489–501. https://doi.org/10.1007/s12052-011-0357-y.

    Article  Google Scholar 

  24. Giblin-Davis RM, Norden BB, Batra SWT, Eickwort GC. Commensal nematodes in the glands, genitalia, and brood cells of bees (Apoidea). J Nematol. 1990;22:150–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Giblin-Davis RM, Davis KA, Taylor GS, Thomas WK. Entomophilic nematode models for studying biodiversity and cospeciation. In: Chen ZX, Chen SY, Dickson DW, editors. Nematology, advances and perspectives, vol. 1. Wallingford: CABI Publishing; 2003. p. 493–540.

    Google Scholar 

  26. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.

    CAS  Google Scholar 

  27. Hasegawa M, Kishino H, Yano T. Dating the human-ape split by a molecular clock of mitochondrial DNA. J Mol Evol. 1985;22:160–74. https://doi.org/10.1007/BF02101694.

    Article  CAS  PubMed  Google Scholar 

  28. Hunt DJ. On Artigasia horridospina n. sp., Longior semialata n. sp., Mentecle magnifica n. sp., Paraxyo ensicrinatus n. sp. (Oxyurida: Hystrignathidae) and Pulchrocephala ? pulchrocephala Travassos, 1925 (Oxyurida: Pulchrocephalidae). Syst Parasitol. 1981;3:33–52. https://doi.org/10.1007/BF00012238.

    Article  Google Scholar 

  29. Hunt DJ. Observations on the genus Indiana Chakravarty, 1943 (Oxyurida: Pulchrocephalidae). Syst Parasitol. 1983;5:215–21.

    Article  Google Scholar 

  30. Hunt DJ. Two new species of Travassosinema Rao, 1958 (Nematoda: Travassosinematidae) from diplopods in Sulawesi and Papua New Guinea. Afro-Asian J Nematol. 1993;3:196–200.

    Google Scholar 

  31. Hunt DJ. Travassosinema thyropygi sp. n. (Nematoda: Travassosinematidae) from a spirobolid millipede from Vietnam with SEM observations on Heth imias Spiridonov, 1989 (Nematoda: Hethidae). Fundam Appl Nematol. 1996;19:7–14.

    Google Scholar 

  32. Hunt DJ, Luc PV, Spiridonov SE. Bizarre Ichthyocephalidae (Nematoda: Rhigonematida) from Vietnam with description of Ichthyocephaloides comatus sp. n., Xystrognathus priscus gen. n., sp. n. and X. phrissus gen. n., sp. n. Nematology. 2002;4:829–43.

    Article  Google Scholar 

  33. Hunt DJ, Moore D. Rhigonematida from New Britain diplopods. 1. The genus Carnoya Gilson, 1898 (Ransomnematoidea: Carnoyidae) with descriptions of three new species. Fundam Appl Nematol. 1998;21:281–97.

    Google Scholar 

  34. Jex AR, Schneider MA, Rose HA, Cribb TH. The Thelastomatoidea (Nematoda: Oxyurida) of two sympatric Panesthiinae (Insecta: Blattodea) from southeastern Queensland, Australia: taxonomy, species richness and host specificity. Nematology. 2005;7:543–75. https://doi.org/10.1163/156854105774384741.

    Article  Google Scholar 

  35. Jex AR, Schneider TH, Cribb TH. The importance of host ecology in thelastomatid (Nematoda: Oxyurida) host specificity. Parasitol Int. 2006;55:169–74. https://doi.org/10.1016/j.parint.2006.03.001.

    Article  PubMed  Google Scholar 

  36. Kim T, Kim J, Cho S, Min GS, Park C, Carreno RA, Nadler SA, Park JK. Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome of Rhigonema thysanophora (Nematoda: Chromadorea). Zool Scripta. 2014;43:289–303. https://doi.org/10.1111/zsc.12047.

    Article  Google Scholar 

  37. Laetsch DR, Heitlinger EG, Taraschewski H, Nadler SA, Blaxter ML. The phylogenetics of Anguillicolidae (Nematoda: Anguillicolidea), swim bladder parasites of eels. BMC Evol Biol. 2012;12:60. https://doi.org/10.1186/1471-2148-12-60.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Lambshead PJD. Recent developments in marine benthic biodiversity. Oceanis. 1993;19:5–24.

    Google Scholar 

  39. Ledón-Rettig CC, Moczek AP, Ragsdale EJ. Diplogastrellus nematodes are sexually transmitted mutualists that alter the bacterial and fungal communities of their beetle host. PNAS. 2018;115:10696–701. https://doi.org/10.1073/pnas.1809606115.

    Article  CAS  PubMed  Google Scholar 

  40. Malysheva SV, Sumaya NHN. Two new species of Heth Cobb, 1898 (Nematoda: Ransomnematoidea) from Northern Mindanao, the Philippines. Nematology. 2017;19:1003–15. https://doi.org/10.1163/15685411-00003102.

    Article  Google Scholar 

  41. Means JC, Hennen DA, Tanabe T, Marek PE. Phylogenetic systematics of the millipede family Xystodesmidae. Insect Syst Diversity. 2021;5:1–26. https://doi.org/10.1093/isd/ixab003.

    Article  Google Scholar 

  42. Mejia-Madrid HH. Three new species of Heth Cobb, 1898 (Ransomnematoidea: Hethidae) from spiroboloid millipedes from Mexico. Nematology. 2014;16:83–98. https://doi.org/10.1163/15685411-00002747.

    Article  Google Scholar 

  43. Mejia-Madrid HH. A molecular phylogeny of the Rhigonematomorpha De Ley & Blaxter, 2002 as inferred from SSU and LSU rDNA sequences. Nematology. 2018;20:547–65. https://doi.org/10.1163/15685411-00003161.

    Article  CAS  Google Scholar 

  44. Mitter C, Brooks DR. Phylogenetic aspects of coevolution. In: Futuyma DJ, Slatkin M, editors. Coevolution: Sinauer Associates, Inc; 1983. p. 65–98.

  45. Morffe J, Hasegawa K. Rhigonema naylae n. sp. (Rhigonematomorpha: Rhigonematidae) a new parasitic nematode from a Japanese polydesmid millipede (Polydesmida: Xystodesmidae). Zootaxa. 2017;4269:277–86. https://doi.org/10.11646/zootaxa.4269.2.6.

    Article  PubMed  Google Scholar 

  46. Morffe J, Hasegawa K. Morphological and molecular characterization of Travassosinema claudiae n. sp. (Oxyuridomorpha: Travassosinematidae) from the Japanese millipede Parafontaria laminata (Attems, 1909) (Polydesmida: Xystodesmidae). Zootaxa. 2017;4282:166–78. https://doi.org/10.11646/zootaxa.4282.1.10.

    Article  Google Scholar 

  47. Morffe J, García N, Hasegawa K, Carreno RA. A new species of Aoruroides Travassos & Kloss, 1958 (Nematoda: Oxyuridomorpha: Thelastomatidae) parasite of the wood-burrowing cockroach Panesthia angustipennis spadica (Shiraki, 1906) (Blattodea: Blaberidae: Panesthiinae) from Japan with comments on the validity of the genus Aoruroides. Zootaxa. 2019;4712:365–76. https://doi.org/10.11646/zootaxa.4712.3.3.

    Article  Google Scholar 

  48. Morffe J, Carreno RA, Ozawa S, Sriwati R, Hasegawa K. On the status of Suifunema peregrinatum (Carreno, 2017) n comb (Nematoda: Oxyuridomorpha: Thelastomatidae) and first record of the species for Japan. Zootaxa. 2019;4651:351–65. https://doi.org/10.11646/zootaxa.4651.2.9.

    Article  Google Scholar 

  49. Nadler SA, Carreno RA, Mejía-Madrid H, Ullberg J, Pagan C, Houston R, Hugot JP. Molecular phylogeny of clade III nematodes reveals multiple origins of tissue parasitism. Parasitol. 2007;134:1421–42. https://doi.org/10.1017/S0031182007002880.

    Article  CAS  Google Scholar 

  50. Ni XF, Barton DP, Chen HX, Li L. A new species of Cosmocerca (Nematoda, Ascaridomorpha) from the marine toad Rhinella marina (Linnaeus) (Anura, Bufonidae) in Australia. Zookeys. 2020;931:11–20. https://doi.org/10.3897/zookeys.931.50478.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Nunn GB. Nematode molecular evolution. An investigation of evolutionary patterns among nematodes based upon DNA sequences. PhD Thesis, University of Nottingham; 1992.

  52. Ozawa S, Vicente CSL, Sato K, Yoshiga T, Kanzaki N, Hasegawa K. First report of the nematode Leidynema appendiculata from Periplaneta fuliginosa. Acta Parasitol. 2014;59:219–28. https://doi.org/10.2478/s11686-014-0230-6.

    Article  PubMed  Google Scholar 

  53. Ozawa S, Morffe J, Vicente CSL, Ikeda K, Shinya R, Hasegawa K. Morphological, molecular and developmental characterization of the thelastomatid nematode Thelastoma bulhoesi (de Magalhães, 1900) (Oxyuridomorpha: Thelastomatidae) parasite of Periplaneta americana (Linnaeus, 1758) (Blattodea: Blattidae) in Japan. Acta Parasitol. 2016;61:241–54. https://doi.org/10.1515/ap-2016-0034.

    Article  PubMed  Google Scholar 

  54. Ozawa S, Hasegawa K. Broad infectivity of Leidynema appendiculatum (Nematoda: Oxyurida: Thelastomatidae) parasite of the smokybrown cockroach Periplaneta fuliginosa (Blattodea: Blattidae). Ecol Evol. 2018;8:3908–18. https://doi.org/10.1002/ece3.3948.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Phillips G, Bernard EC, Pivar RJ, Moulton JK, Shelley RM. Coronostoma claireae n. sp. (Nematoda: Rhabditida: Oxyuridomorpha: Coronostomatidae) from the Indigenous Milliped Narceus gordanus (Chamberlain, 1943) (Diplopoda: Spirobolida) in Ocala National Forest, Florida. J Nematol. 2016;48:159–69. https://doi.org/10.21307/jofnem-2017-023.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Pilotte N, Unnasch TR, Williams SA. The current status of molecular xenomonitoring for Lymphatic filariasis and onchocerciasis. Trends Parasitol. 2017;33:788–98. https://doi.org/10.1016/j.pt.2017.06.008.

    Article  PubMed  Google Scholar 

  57. Poinar GO, Ed. Entomogenous nematodes: a manual and host list of insect-nematode associations. Brill, Leiden, 1975. pp. 317.

  58. Rao PN. Studies on the nematode parasites of insects and other arthropods. Arquivos do Museo Nacional, Rio de Janeiro. 1958;46:33–84.

    Google Scholar 

  59. Rózsa L, Reiczigel J, Majoros G. Quantifying parasites in samples of hosts. J Parasitol. 2000;86:228–32. https://doi.org/10.1645/0022-3395(2000)086[0228:QPISOH]2.0.CO;2.

    Article  PubMed  Google Scholar 

  60. Rusconi JM, Camino NB, Achinelly MF. Nematodes (Mermithidae) parasitizing grasshoppers (Orthoptera: Acrididae) in the Pampean region, Argentina. Braz J Biol. 2017;77:12–5. https://doi.org/10.1590/1519-6984.06015.

    Article  CAS  PubMed  Google Scholar 

  61. Singh SK, Hodda M, Ash GJ, Banks NC. Plant-parasitic nematodes as invasive species: characteristics, uncertainty and biosecurity implications. Ann Appl Biol. 2013;163:319–514. https://doi.org/10.1111/aab.12065.

    Article  Google Scholar 

  62. Sinnott D, Carreno RA, Herrera H. Distribution of thelastomatid nematode (Nematode: Oxyurida) in endemic and introduced cockroaches on the Galápagos island Archipelago. Ecuador J Parasitol. 2015;101:445–57. https://doi.org/10.1645/15-721.1.

    Article  PubMed  Google Scholar 

  63. Spiridonov SE, Cribb TH. Travassosinema dalei sp. n. (Thelastomatoidea; Nematoda) from the hind gut of an Australian scarabaeid larva. New Zealand J Zool. 2012;39:291–300. https://doi.org/10.1080/03014223.2012.668855.

    Article  Google Scholar 

  64. Sriwati R, Ozawa S, Morffe J, Hasegawa K. First record of Hammerschmidtiella diesingi (Hammerschmidt, 1838) (Oxyuridomorpha: Thelastomatidae) parasite of Periplaneta Americana (Linnaeus, 1758) (Blattodea: Blattidae) in Japan, morphological and molecular characterization. Acta Parasitol. 2016;61:720–8. https://doi.org/10.1515/ap-2016-0101.

    Article  PubMed  Google Scholar 

  65. Takakuwa Y. Rysodesmus Arten aus Japan. Trans Natl History Soc Formosa. 1941;31:413–5.

    Google Scholar 

  66. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9. https://doi.org/10.1093/molbev/mst197.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Tanabe T. Revision of the milipede genus Parafontaria Verhoeff 1936 (Diplopoda, Xystodesmidae). J Nat History. 2002;36:2139–83.

    Article  Google Scholar 

  68. Tanabe T. Two new species of the Genus Riukiaria from Kyushu and Is. Yaku-shima, Japan (Diplopoda: Polydesmida: Xystodesmidae). Acta Arachnol. 1988;37:37–45.

    Article  Google Scholar 

  69. van den Hoogen J, Geisen S, Crowther TW. Soil nematode abundance and functional group composition at a global scale. Nature. 2019;572:194–8. https://doi.org/10.1038/s41586-019-1418-6.

    Article  CAS  PubMed  Google Scholar 

  70. van Megen H, van den Elsen S, Holterman M, Karssen G, Mooyman P, Bongers T, Holovachov O, Bakker J, Helder J. A phylogenetic tree of nematodes based on about 1200 full-length small subunit ribosomal DNA sequences. Nematology. 2009;11:927–50. https://doi.org/10.1163/156854109X456862.

    Article  CAS  Google Scholar 

  71. Vicente CSL, Ozawa S, Hasegawa K. Composition of the cockroach gut microbiome in the presence of parasitic nematodes. Microbes Environ. 2016;31:314–20. https://doi.org/10.1264/jsme2.ME16088.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Vicente CSL, Ozawa S, Hasegawa K. The composition of hindgut microbiota of Periplaneta japonica in the presence of thelastomatid parasitic nematodes. Nematol Res. 2018;48:19–26. https://doi.org/10.3725/jjn.48.19.

    Article  Google Scholar 

  73. Viney ME, Lok JB. The biology of Strongyloides spp., WormBook, ed. The C. elegans Research Community, WormBook, 2015.

  74. Zervos S. Evidence for population self-regulation, reproductive competition and arrhenotoky in a thelastomatid nematode of cockroaches. Parasitol. 1988;96:369–79. https://doi.org/10.1017/S0031182000058352.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Yuto Koike, Chernysh Yelizaveta, Fumi Takasuka, and Anna Yasuda (Hasegawa lab, Chubu University), for their help with sample collection. We also thank Jans Morffe (Instituto de Ecología y Sistemática, Cuba) for his advice on nematode identification.

Funding

This work was supported by research funds from the Research Institute for Biological Function, and the College of Bioscience and Biotechnology, Chubu University.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, KH; Funding Acquisition, KH; Host collection and identification, SN, TT, KH; Nematode species identification, SN, KH; Parasite population analysis, SN, KS, KH; Molecular phylogenetic analysis, SN, KH; Supervision, KH; Writing, KH, KS All authors have read and approved the manuscript.

Corresponding author

Correspondence to Koichi Hasegawa.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Photographs of xystodesmid millipedes studied in this experiment.

Additional file 2: Figure S2.

Average ± S.D. of host body size (mm) in infection condition (infected with only R. naylae or co-infected with R. naylae and T. claudiae) or in different seasonal conditions (spring or summer).

Additional file 3: Table S1.

Pairwise differences % in the D2D3 sequence (730 bp) between 15 samples of Rhigonematoidea spp.

Additional file 4: Table S2.

Pairwise differences % in the D2D3 sequence (704 bp) between 13 samples of Travassosinematidae spp.

Additional file 5: Table S3.

Population of the two parasitic nematodes in Parafontaria millipedes.

Additional file 6: Table S4.

Pairwise differences % in the D2D3 sequence (676 bp) between eight samples of Thelastomatidae spp.

Additional file 7: Table S5.

Population of the parasitic nematodes in Riukiaria millipedes.

Additional file 8: Table S6.

List of DNA sequence used for molecular phylogenetic analysis.

Additional file 9: Table S7.

Datasheet of parasitic nematodes isolated from xystodesmid millipedes.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nagae, S., Sato, K., Tanabe, T. et al. Symbiosis of the millipede parasitic nematodes Rhigonematoidea and Thelastomatoidea with evolutionary different origins. BMC Ecol Evo 21, 120 (2021). https://doi.org/10.1186/s12862-021-01851-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-021-01851-4

Keywords