New insights into phylogenetic relationships of Rhabdocoela (Platyhelminthes) including members of Mariplanellida
BMC Zoology volume 8, Article number: 9 (2023)
Previous flatworm phylogenetic research has been carried out analysing 18S and 28S DNA markers. Through this methodology, Mariplanellinae subfamily has been recently re-classified as Mariplanellida status novus. This new classification implied that 3 genera belonged to Mariplanellida: Mariplanella, Lonchoplanella and Poseidoplanella. In this study, we aim to clarify some of the relationships within Rhabdocoela analysing 18S and 28S DNA markers of a total of 91 species through Maximum Likelihood and Bayesian Inference methodologies. A total of 11 species and genera, including Lonchoplanella, from the island of Sylt are included and had not previously been involved in any molecular phylogenetic analyses.
Our phylogenetic results support Mariplanellida as an independent group within Rhabdocoela and its status as an infraorder. Our study suggests that Lonchoplanella axi belongs to Mariplanellida. Within Rhabdocoela, Haloplanella longatuba is nested within Thalassotyphloplanida, instead of Limnotyphloplanida. Within Kalyptorhynchia, the taxon Eukalyptorhynchia turned out to be paraphyletic including members of Schizorhynchia. These results also support the position of the genus Toia separate from Cicerinidae.
Lonchoplanella axi belongs to Mariplanellida, whose status as infraorder is herein confirmed. The genus Toia belongs separate from Cicerinidae. Further research is needed to clarify the phylogenetic relationships of Hoploplanella. Most of the species, genera and families included in this study with more than one terminal are monophyletic and well supported. Adding gene markers and complementary morphological studies will help to clarify those relationships that remain uncertain.
Flatworms are a large group in terms of diversity, with more than 26,500 described species [1, 2]. Most of them are parasitic, while around 6500 species of them are “free-living Platyhelminthes”. The parasitic flatworms, Neodermata (comprising Trematoda, Monogenea, and Cestoda), are a well-defined and supported clade characterized by a syncytial, nonciliated epidermis whose nuclei-bearing parts lie sunken below the musculature. The rest of platyhelminthes are mostly free-living but also symbiotic flatworms, also known as turbellarians (non-cladistic group); most of them with a ciliated, cellular epidermis .
The free living Rhabdocoela and Proseriata are the two most diverse microturbellarian orders, which, in recent phylogenetic hypotheses, are basally branching within Euneoophora  (Sup. Figure 1). Proseriata has above 400 described species , and Rhabdocoela, has about 1530 described species . Nevertheless, these numbers represent a scarce amount of the estimated total number of microturbellarian species present on Earth, which was estimated to be around 44.000 by Armonies . Further research is needed in order to record their diversity more accurately .
Proseriata, is a monophyletic order ; based on morphological traits Sopott-Ehlers  divided it into Unguiphora (taxa without cirrus and a statocyst but with pigment in the mantle cells of rhabdomeric receptors) and Lithophora (statocyst present and no pigment in the mantle cells). Later the monophyly of both clades was confirmed [9, 11]), though with some dispute in the Unguiphora as it was underrepresented in the studies supporting its monophyly and appears to be paraphyletic in others [11, 12]). Within Proseriata, the position of the genus Ciliopharyngiella Ax, 1952 has been recently debated. Curini-Galletti et al.  suggested that it belongs inside Proseriata clustered with Unguiphora, while Van Steenkiste & Leander  suggested Ciliopharyngiella being the sister lineage of Proseriata.
Rhabdocoela is a monophyletic order  of free-living Platyhelminthes, traditionally subdivided into two major groups: Kalyptorhynchia and Dalytyphloplanida, with and without an anterior proboscis, respectively. Within Rhabdocoela, the family Mariplanellidae was initially included in Dalytyphloplanida though not fully fitting with the diagnostic characters. Recently, Steenkiste & Leander  re-classified the subfamily Mariplanellinae to Mariplanellida status novus, representing the monophyletic group sister to a large clade comprising Kalyptorhynchia and Dalytyphloplanida. According to this classification, Mariplanellidae (the only family of Mariplanellida) currently comprises three genera: Mariplanella Ax & Heller, 1970, Lonchoplanella Ehlers, 1974, and Poseidoplanella Willems et al., 2005. Morphologically, these genera share the characters single ovary and a double connection in the female reproductive system. However, the phylogenetic analyses by Van Steenkiste and Leander  only included two species of Mariplanella, M. piscadera Van Steenkiste & Leander, 2022 and M. frisia Ax & Heller, 1970, while members of Lonchoplanella and Poseidoplanella have not been included in any molecular phylogenetic analysis up to date.
Trees obtained from independent Maximum Likelihood (ML) analyses of single markers (18S, 28S), as well as the one obtained from the analyses of the concatenated data matrix (18S + 28S) showed congruent results (SFig. 2, SFig. 3 and Fig. 1, respectively). The concatenated data matrix (18S + 28S) run through Bayesian Inference (BI) yielded highly congruent results (SFig. 4) with practically the same tree topology and high support values. The trees were rooted in Ciliopharyngiella constricta Martens & Schockaert, 1981, excepting the tree of 28S which was rooted in Proseriata.
In the ML concatenated analysis (Fig. 1), the taxon Proseriata is not well supported (68 B, where B represents bootstrap support values). Within Proseriata, Unguiphora (excepting Nematoplana sp.) and Lithophora are monophyletic and well supported (99B and 80B, respectively). Within Unguiphora, both Nematoplana coelogynoporoides Meixner, 1938 are in a monophyletic group; however, Nematoplana sp. is sister to Lithophora. The genus Polystyphora Ax, 1958, represented by three species, is monophyletic. Lithophora is subdivided into Otoplanidae and Monocelididae (100 B), and Coelogynoporidae and Calviriidae (100 B), respectively, as sister groups.
Rhabdocoela is well supported (100 B). Within Rhabdocoela, the three clades, Mariplanellida, Dalytyphloplanida and Kalyptorhynchia, are strongly supported (100 B). Mariplanellida is sister to Dalytyphloplanida and Kalyptorhynchia, which are joined in a not highly supported clade (69 B). In the concatenated analysis (Fig. 1), the species Lonchoplanella axi is located within Mariplanellida, sister to M. axi (86 B), showing the genus Mariplanella, as currently delineated, as paraphyletic. The 18S tree (SFig. 2) shows Mariplanella as monophyletic but is not supported (61 B) with Lonchoplanella axi as its sister group. Within Dalytyphloplanida, Neodalyellida, Limnotyphloplanida (excepting Haloplanella longatuba Ax & Heller, 1970) and Thalassotyphloplanida are monophyletic and, excepting the latter, well supported (100 B, 100 B, 96 B, and 84 B, respectively). Limnotyphloplanida and Thalassotyphloplanida are sister to each other (100B). The species Haloplanella longatuba is located within Thalassotyphloplanida, out of the Limnotyphloplanida.
Within Kalyptorhynchia, Eukalyptorhynchia is paraphyletic, including Schizorhynchia. The genus Toia Markus, 1952 is located as the sister group of a large clade (100 B) containing the rest of Eukalyptorhynchia (83 B) and Schizorhynchia, which is not supported (48B). These two latter groups are better supported in the BI analisis (SFig. 4) with a posterior probability value of 1. The family Cheliplanidae appears to be paraphyletic (100 B) with Cheliplanilla caudata Meixner 1938 being the sister taxon of a group compromised by Carcharodorhynchus Meixner, 1938 and Schizorhynchidae, represented by more than one species, and also well supported (100 B). The genera Carchadorhynchus, Schizochilus Boaden, 1963 and Proschizorhynchus Meixner, 1928, and Psammorhynchus tubulipenis Meixner, 1938, including more than one terminal, are respectively monophyletic, though herein represented by few species.
Our results show Mariplanellida as an independent group within Rhabdocoela as found by Van Steenkiste & Leander , and previously suggested by other molecular analyses . Its higher-level status as an additional infraorder proposed by Van Steenkiste & Leander  within Rhabdocoela is supported herein as well. Our study included Lonchoplanella axi, a genus and species not incorporated before in a molecular phylogeny. Lonchoplanella axi shares a number of conspicuous characters with Mariplanella frisia, including two types of adenal rhabdites (needle-shaped and elongate viscous) and a muscular copulatory bursa with a sclerotized basal membrane  and had been described as a member of Mariplanellidae . The topology of the concatenated tree confirms Lonchoplanella axi as belonging to the Mariplanellidae family and Mariplanellida infraorder. More information (terminals and markers) is needed to clarify the monophyly of both genera, and further analyses including Poseidoplanella halleti Willems et al., 2005 (presumably the only missing taxon of Mariplanellidae) will clarify the relationships between the three genera.
Our study included six new terminals and five species of Kalyptorhynchia in addition to those included by Van Steenkiste & Leander . Results corroborate Eukalyptorhynchia as non-monophyletic supporting Willems et al.  and Tessens et al.  results. Our results also support the position of the genus Toia separate from Cicerinidae, as suggested by Tessens, et al. . The rest of Cicerinidae are monophyletic and the sister group of Schizorhynchia. Cicerinidae (except Toia) were located in a polytomy with Schizorhynchia and most members of Eukalyptorhynchia in the analysis of Van Steenkiste & Leander .
Most of the species, genera and families included in this study with more than one terminal are monophyletic and well supported. However, the terminals Nematoplana sp. and Haloplanella longatuba need further consideration. Within Proseriata, Nematoplana sp. (downloaded from GenBank) is not closely related to N. coelogynoporoides; its sister relationship to Lithophora might indicate that this terminal could have been misidentified or the result of contamination in the 18S sequence. Additionally, the length of branches of Nematoplana coelogynoporoides from Sylt (included herein) and N. coelogynoporoides from Roscoff (downloaded from GenBank) might indicate they are not be the same species.
In Rhabdocoela, Haloplanella longatuba is one of the species we incorporate in a phylogenetic analysis for the first time with molecular information. The relationships of brackish and marine water Typhloplanidae species, such as Haloplanella longatuba have been previously discussed [6, 16, 17]. Rieger  describes the resemblance between certain genera in the family Typhloplanidae (Limnotyphloplanida parvorder) and various genera in the Thalassotyphloplanida parvorder, thus encouraging future reorganizations within these taxonomic groups. Hochberg & Cannon  remarks on the presence of an unusual character in some genera of the Typhloplanidae family, such as Haloplanella Luther, 1946 and Pratoplana Ax 1960, where a stylet is present in the copulatory apparatus instead of a cirrus, which is one of the family's ground pattern traits. Moreover, several similarities between this species and members of Thalassotyphloplanida have been found, such as the female genital canal and the proboscis structure. Van Steenkiste et al.  already suggested that several brackish water and marine Typhloplanidae taxa might be closely related to Byrsophlebidae (Thalassotyphloplanida). Our results show Haloplanella longatuba nested within Thalassotyphloplanida, which supports that the position of this taxon has to be taken into further consideration, and revised in future studies.
In this study, the selected gene markers, 18S and 28S were used because they are already available for a large number of species and have been previously found useful to discern phylogenetic relationships within Proseriata and Rhabdocoela . Nevertheless, Next Generation Sequencing (NGS) techniques will be undoubtedly useful to discern and clarify the evolution of these groups.
Our results support the infraorder Mariplanellida as an independent group within Rhabdocoela and confirms, by the first time with molecular data, that Lonchoplanella axi belongs to Mariplanellida. Eukalyptorhynchia is paraphyletic including members of Schizorhynchia, and the genus Toia separate from Cicerinidae. Haloplanella longatuba is herein nested within Thalassotyphloplanida, and not in Typhloplanidae (Limnotyphloplanida), which suggests that further studies are needed to clarify its phylogenetic relationships.
More terminals and information from morphological studies, as well as new markers of NGS techniques may clarify the gaps and the still doubtful relationships.
Material and methods
In this study we aim to provide a more robust phylogenetic hypothesis of Rhabdocoela relationships. For this purpose, we introduced the species Lonchoplanella axi to phylogenetic analyses to test whether or not it belongs to Mariplanellida. With respect to Kalyptorhynchia and Dalytyphloplanida we added species of Cystiplana Karling, 1964, Haloplanella, Marirhynchus Schilke, 1970, Lonchoplanella, and Schizorhynchoides Meixner, 1928 (Table 1).
Sampling and species identification
A total of 37 samples were taken during two days sampling from intertidal sand in the island of Sylt, Germany. Samples were obtained by digging on the substrate with a 10 cm long shovel. The substrate was kept in zip bags and stored in fridges at 4ºC.
Sampling sites were the beach besides List Harbour (55.015337N, 8.435999E) and the beach in front of Alfred-Wegener-Institute building (55.023745, 8.439049), always during low tide. The collected sediment samples were all coarse sand enriched with variable amounts of organic material. Meiofauna was separated from the sediment using the MgCl2 decantation method .
Flatworms were morphologically identified under Leica S APO stereomicroscope and Leica DM 2500 microscope and photographed (stylets) under a portable Leica MC 190 HD attached camera (Fig. 2).
DNA extraction, amplification and sequencing
For the DNA extraction the DNeasy® Blood & Tissue Kit (QIAGEN) was used. Manufacturer’s instructions were followed, with the exception that DNA was eluted in 60µL of preheated AE elution buffer (60 °C). For samples with low concentration, this protocol was followed by the Amplification of purified genomic DNA protocol from QIAGEN REPLI-g® kit. Thermocycling conditions and primers from Van Steenkiste & Leander  were used to sequence markers 18S and 28S (see supplementary material, S1). Sequencing was carried out by Eurofins Genomics (Konstanz, Germany). All new sequences were deposited in GenBank, and sequence accession numbers are provided in Table 1.
Once the sequences were obtained, they were blasted using blastn through ncbi-blast + v.2.12.0 to confirm that platyhelminthes’ DNA was amplified during the PCRs. The rest of the sequences were obtained from GenBank attempting to gather a broad representation of the different families and subfamilies. Those terminals for which both markers (18S and 28S) were available were selected for this study. Several terminals of Proseriata, including Ciliopharyngiella constricta were also included to root the tree.
Sequences were visually checked in Geneious v10.2.3 and aligned using MAFFT v.7.305b  using the iterative refinement method E-INSI. Single genes (18S, 28S) were concatenated using FASconCAT-G [20, 21]. The maximum likelihood (ML) analysis of the single markers, as well as the concatenated matrix was performed through IQtree v.188.8.131.52 [22, 23], with best fitting models selected by Modelfinder  (18S: GTR + F + I + G4; 28S: GTR + F + I + G4). In all analyses, each partition was allowed to have its own set of branch lengths. (-spp option). Support values were estimated based on 1000 bootstrap pseudo replicates (B). iTol v.6. and Adobe Illustrator (2020) were used to edit the phylogenetic trees. The concatenated matrix (18S + 28S) was also analysed through Bayesian inference (BI). For BI analyses, two independent runs of 1,342,000 generations and four chains, each (one cold, three heated) were run in MrBayes 3.2.7 . The most similar models available in MrBayes (-mset option) to those selected by Modelfinder for each partition were applied. All parameters were unlinked, rates were allowed to vary freely over partitions andtrees were sampled every 1000 generations. The runs were stopped when the standard deviation reached the value of 0,007. After discarding 25% first trees as burn-in, trees from the stationary phase were combined to obtain a majority rule consensus and posterior node probabilities .
Availability of data and materials
The data generated during and/or analyzed during the current study are available in GenBank. Accession numbers are displayed in Table 1.
Caira JN, Littlewood DTJ. Worms, Platyhelminthes. Encyclopedia of Biodiversity. 2013:437–69. https://doi.org/10.1016/b978-0-12-384719-5.00166-0.
Tyler S, Schilling S, Hooge M, Bush LF. Turbellarian taxonomic database. Version 2.0. 2006–2022. http://turbellaria.umaine.edu. Accessed 07 Jan 2022.
Tyler S, Hooge M. Comparative morphology of the body wall in flatworms (Platyhelminthes). Can J Zool. 2004;82(2):194–210. https://doi.org/10.1139/z03-222.
Laumer CE, Giribet G. Inclusive taxon sampling suggests a single, stepwise origin of ectolecithality in Platyhelminthes. Biol J Lin Soc. 2014;111(3):570–88. https://doi.org/10.1111/bij.12236.
Van Steenkiste N, Tessens B, Willems W, Backeljau T, Jondelius U, Artois T. A Comprehensive Molecular Phylogeny of Dalytyphloplanida (Platyhelminthes: Rhabdocoela) Reveals Multiple Escapes from the Marine Environment and Origins of Symbiotic Relationships. PLoS ONE. 2013;8(3). https://doi.org/10.1371/journal.pone.0059917.
Armonies W. Uncharted biodiversity in the marine benthos: The void of the smallish with description of ten new Platyhelminth taxa from the well-studied North Sea. Helgoland Marine Research. 2018;72(1). https://doi.org/10.1186/s10152-018-0520-8.
Schockaert ER, Hooge M, Sluys R, Schilling S, Tyler S, Artois T. Global diversity of free living flatworms (Platyhelminthes, “Turbellaria”) in freshwater. Hydrobiologia. 2008;595(1):41–8. https://doi.org/10.1007/s10750-007-9002-8.
Willems WR, Wallberg A, Jondelius U, Littlewood DTJ, Backeljau T,Schockaert ER, Artois TJ. Filling a gap in the phylogeny offlatworms: relationships within the Rhabdocoela (Platyhelminthes), inferred from 18S ribosomal DNA sequences. Zoologica Scripta. 2006;35:1–17.
Sopott-Ehlers B. The phylogenetic relationships within the Seriata (Platyhelminthes). In: Conway Morris S, George DJ, Gibson R, Platt HM, editors. The origin and relationships of lower invertebrates. Oxford: Oxford University Press; 1985. p. 159–67.
Curini-Galletti M, Webster BL, Huyse T, Casu M, Schockaert ER, Artois TJ, Littlewood DTJ. New insights on the phylogenetic relationships of the Proseriata (Platyhelminthes), with proposal of a new genus of the family Coelogynoporidae. Zootaxa. 2010;2537:1–18. https://doi.org/10.11646/zootaxa.2537.1.1.
Laumer CE, Giribet G, Curini-Galletti M. Prosogynopora riseri, gen. et sp. nov., a phylogenetically problematic lithophoran proseriate (Platyhelminthes:Rhabditophora) with inverted genital pores from the New England coast. Invertebrate Systematics. 2014;28(3):309–325. https://doi.org/10.1071/IS13056.
Van Steenkiste NWL, Leander BS. The molecular phylogenetic position of Mariplanella piscadera sp. nov. reveals a new major group of rhabdocoel flatworms: Mariplanellida status novus (Platyhelminthes: Rhabdocoela). Org Divers Evol. 2022:0123456789. https://doi.org/10.1007/s13127-022-00542-2.
Ehlers U. Interstitielle Typhloplanoida (Turbellaria) aus dem Litoral der Nordseeinsel Sylt. Mikrofauna des Meeresbodens, 49. vol. 102. Mainz: Akademie der Wissenschaften und der Literatur; 1974. p. 427–526. ISBN 3-515-02037-3.
Tessens B, Janssen T, Artois T. Molecular phylogeny of Kalyptorhynchia (Rhabdocoela, Platyhelminthes) inferred from ribosomal sequence data. Zoolog Scr. 2014;43(5):519–30. https://doi.org/10.1111/zsc.12066.
Rieger MR. A new group of Turbellaria-Typhloplanoida with a proboscis and its relationship to Kalyptorhynchia. In: Morse MP, editor. Riser NW. Mc- Graw-Hill Company: Biology of the Turbellaria. Libbie H Hyman Memorial Volume. New York; 1974. p. 23–62.
Hochberg R, Cannon LRG. Two new freshwater rhabdocoels, Austrodalyellia gen. nov. and Haplodidymos gen. nov. (Platyhelminthes), from Queensland, Australia. Zootaxa. 2002;44(1):1. https://doi.org/10.11646/zootaxa.44.1.1.
Schockaert ER. The Importance of Turbellarians in Ecosystems. In: Hall GS, editor. Methods for the Examination of Organismal Diversity in Soils and Sediments. Wallingford: CAB International; 1996. p. 211–25.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6. https://doi.org/10.1093/bib/bbx108.
Kück P, Longo GC. FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Front Zool. 2014;11:81. https://doi.org/10.1186/s12983-014-0081-x.
Kück P, Meusemann K. FASconCAT: Convenient handling of data matrices. Mol Phylogenet Evol. 2010;56(3):1115–8. https://doi.org/10.1016/j.ympev.2010.04.024.
Chernomor O, von Haeseler A, Minh BQ. Terrace aware data structure for phylogenomic inference from supermatrices. Syst Biol. 2016;65:997–1008. https://doi.org/10.1093/sysbio/syw037.
Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol Biol Evol. 2015;32:268–74. https://doi.org/10.1093/molbev/msu300.
Kalyaanamoorthy S, Minh BQ, Wong T, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9. https://doi.org/10.1038/nmeth.4285.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42. https://doi.org/10.1093/sysbio/sys029.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5. https://doi.org/10.1093/bioinformatics/17.8.754.
We would like to acknowledge Conrad Helm, Leonard Spohr and Christoph Bleidorn (Georg August University, Göttingen) for their help during the field trips, sorting samples and support to conduct this study. We are also grateful to the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Wattenmeerstation Sylt for hosting us during our sampling trips.
Open Access funding enabled and organized by Projekt DEAL. This study was financed by the Biodiversitätsmuseum Göttingen (PI:MTA).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary table 1
Phylogenetic relationships between the mayor clades of Platyhelminthes phylum (retrieved from Laumer & Giribet, 2014).
ML phylogenetic tree inferred from 18S gene. Species provided by this study in red. Bootstrap support values under /beside nodes. Values below 70% not represented.
ML phylogenetic tree inferred from 28S gene. Species provided by this study in red. Bootstrap support values under /beside nodes. Values below 70% not represented.
Majority rule consensus tree from BI analysis obtained from the concatenated data set (18S+28S). Posterior probability support values close to each node.
Supplementary figure legends
About this article
Cite this article
Vicente-Hernández, Í., Armonies, W., Henze, K. et al. New insights into phylogenetic relationships of Rhabdocoela (Platyhelminthes) including members of Mariplanellida. BMC Zool 8, 9 (2023). https://doi.org/10.1186/s40850-023-00171-y
- Free-living Platyhelminthes
- Maximum likelihood
- Bayesian Inference