“…This was not surprising, as nucleotide substitution rates in the mitochondrial genome in sponges and other lower metazoans are up to 100 times slower than other animal lineages (Hellberg, 2006; Shearer et al, 2002). In this sense, although COI has been previously found to be useful to explain speciation and phylogeographic patterns in some haplosclerid sponges (see, for example, López-Legentil and Pawlik, 2009; DeBiasse and Hellberg, 2015), it seems to be unsuitable for resolving phylogenetic relationships in other sponges (see, for example, Duran et al , 2004; Dailianis et al , 2011). In light of these results, two possible explanations arise: the first one is that they are separate but closely related species that have undergone secondary contact favoring introgression, and a second explanation may be that I. variabilis and I. fasciculata are the same species but are currently undergoing incipient speciation without reproductive isolation yet occurring, and therefore a certain level of mixing still exists.…”