Dugesia naiadis, 2013

Sluys, Ronald, Solà, Eduard, Gritzalis, Konstantinos, Vila-Farré, Miquel, Mateos, Eduardo & Riutort, Marta, 2013, Integrative delineation of species of Mediterranean freshwater planarians (Platyhelminthes: Tricladida: Dugesiidae), Zoological Journal of the Linnean Society 169 (3), pp. 523-547 : 526

publication ID	https://doi.org/ 10.1111/zoj.12077
persistent identifier	https://treatment.plazi.org/id/03DB7724-E201-836C-FC32-FDD5FD5FFEEC
treatment provided by	Marcus
scientific name	Dugesia naiadis
status

Treatment

DUGESIA View in CoL POPULATIONS

We performed a GMYC approach ( Pons et al., 2006; Fontaneto et al., 2007) to compare the units delimited by this method with those identified in the morphological analysis and to detect possible cryptic species. We used the partial COI sequences of 155 individuals of Dugesia from 34 localities ( Table S1). GMYC detects the change from population processes (coalescence of alleles) to speciation and extinction processes through analysis of branching rate patterns, setting a threshold between the inter- and intraspecific relationships. To obtain the ultrametric tree necessary for this approach, we conducted a phylogenetic analysis in BEAST v1.7.3 ( Drummond & Rambaut, 2007), using a fragment of COI (745 bp) from 2–5 individuals per sampling locality ( Table S1). A lognormal relaxed clock with a substitution rate of 0.017 substitutions per lineage and per million years was applied (cf. Solà et al., 2013). The analysis was run under a GTR + I + Γ evolutionary model. Three monophyletic clades were forced: (1) Dugesia species , without D. sicula and D. naiadis (used as outgroup); (2) Dugesia species , without D. sicula , D. naiadis and Dugesia from Central Europe; (3) Dugesia species , without D. sicula Lepori, 1948 , D. naiadis Sluys sp. nov., Dugesia from Central Europe and D. cretica ( Meixner, 1928; Solà et al., 2013). Monte Carlo Markov chains were run for 150 000 000 generations, sampling every 15 000 trees. The parameters were checked to have reached an effective sampling size (ESS) value of over 100 after a 10% burn-in with Tracer v.1.5 ( Rambaut & Drummond, 2007).

The BEAST tree obtained was submitted to the SPLITS (SPecies LImits by Threshold Statistics; Ezard, Fujisawa & Barraclough, 2009) package for R (available at http://r-forge.r-project.org/projects/ splits/), which implements the GMYC approach. The program also performs likelihood ratio tests (LRTs) between (a) the null and GMYC models to test whether one or multiple species are involved, and (b) single and multiple threshold options.

Abbreviations used in Figures 3–18 View Figure 3 View Figures 4, 5 View Figure 6 View Figures 7, 8 View Figure 9 View Figure 10 View Figures 11, 12 View Figure 13 View Figure 14 View Figures 15, 16 View Figure 17 View Figure 18 : bc, bursal canal; cb, copulatory bursa; cg, cement glands; cod, common oviduct; cs, cyanophilic secretion; dpf, dorsal penial fold; ed, ejaculatory duct; fl, flap; go, gonopore; in, intestine; od, oviduct; pg, penial glands; ph, pharynx; pp, penis papilla; sg, shell gland; spf, spermatphore; sv, seminal vesicle; te, testis; vd, vas deferens.