BACKGROUND
dugesia sicula lepori,  <dig> is a freshwater planarian  typically found in ponds, streams and springs close to the mediterranean coast.
the first record of d. sicula described a population comprising both sexual and asexual  individuals in catania, sicily  <cit> .
fissiparity is an asexual form of reproduction that involves transversal division of the individuals into two fragments and subsequent regeneration of the absent structures.
fissiparous planarians do not develop a copulatory apparatus, which precludes the assignment of these populations to any species since this structure provides most of the defining characters for species description in planarians.
however, the morphological characters of the copulatory apparatus of the sexual specimens allowed the species to be described.
d. sicula was subsequently found on the islands of elba  <cit>  and mallorca  <cit> .
none of the originally described sexual populations can be currently found  <cit> ; but, some sexual populations in northern africa  have been recently reported  <cit> .
moreover, a molecular-based study  <cit>  showed that most of the numerous fissiparous populations of genus dugesia present in the western mediterranean belong to d. sicula.
under certain conditions  triploid fissiparous planarians can develop hyperplasic ovaries and a copulatory apparatus  <cit> , these individuals are known as ex-fissiparous.
in d. sicula, it has been observed that some ex-fissiparous individuals laid cocoons but remained sterile  <cit> , which indicated that asexual populations have lost the capability to sexually reproduce.
the literature suggests that d. sicula is widely distributed in the mediterranean area from morocco to greece  <cit> , and our own observations show that this species mostly comprise fissiparous populations, most likely sexually sterile, distributed throughout the mediterranean, as well as the canary islands, being the only species of this genus in the mediterranean with a such a wide distribution.
contrasting with this wide distribution the species presents extremely low variability for the coi gene  <cit>  among populations located hundreds of kilometres apart.
these two features  could be a consequence of a recent expansion promoted by human activities, as has been shown for other organisms  <cit> .
alternatively, the low levels of variability could be exclusively a consequence of their fissiparous reproduction, as described in schmidtea mediterranea <cit>  and other organisms  <cit> , and its wide distribution may be a consequence of its active spread over the years.
in both cases, the hypotheses contrast with the poor dispersion capability assumed for freshwater planarians  <cit> .
freshwater planarians do not exhibit larval dispersal stages or forms resistant to desiccation; these individuals thus require contiguous freshwater bodies to survive and disperse  <cit> .
however, introductions due to human activities have been documented in other freshwater planarians, such as girardia tigrina <cit>  and schmidtea polychroa <cit>  and more recently have also been proposed for other dugesia species  <cit> , but in all these cases the introduced animals either have a restricted distribution  or its introduction and progression in the new areas has been followed by scientist.
fissiparous reproduction can have effects on the distribution of genetic variability within the populations and even on individuals that could perhaps help disentangle this situation.
regeneration after fission is driven by neoblasts  <cit> , which are pluripotent stem cells.
a new mutation in a neoblast would expand in the individual when this neoblast generates new tissues after the fission of the animal, thus increasing frequency of the mutation in the individual.
following replication via several fission cycles, the new mutation could disappear or spread in the population of cells within the individual and the organism’s population by genetic drift .
in this manner, we would expect to observe cells with different sequences within one individual.
this could be considered a special case of heteroplasmy that differs from cases in which heteroplasmy is transmitted sexually or parthenogenetically through the mitochondria present in the oocyte  where all cells from a new individual will present the same haplotypes.
heteroplasmy has been detected in many species  <cit> , and some studies in humans have shown it increases with age in somatic cells and is more frequent in certain tissues, such as muscle tissue, most likely because of increases in the number of mutations caused by the presence of oxidative radicals  <cit> .
heteroplasmy has also been detected in coi sequences of fissiparous populations of d. japonica <cit> .
due to their characteristics, in fissiparous planarians the somatic heteroplasmy would be maintained in populations and increased over time.
if this situation holds in fissiparous populations of d. sicula, the analysis of their genetic variability could help in discriminating between a recent or ancient origin of the species populations and hence provide clues to the origin of its distribution.
the absence or low frequency of heteroplasmic specimens in a specific locality of fissiparous reproducing individuals could be interpreted as a consequence of a recent bottleneck, potentially because of recent colonisation by a few individuals or recovery from a recent population-wide decrease.
the alternate condition, the presence of many heteroplasmic individuals, particularly those exhibiting a high number of different low-frequency haplotypes, could then potentially be related to aged individuals and populations, although a recent colonisation by highly heteroplasmic individuals could also explain this outcome.
following the same reasoning, one would not expect to find heteroplasmy, or different copies of mitochondrial genes within an individual in sexual organisms.
in the present study, we analysed the genetic diversity and structure  of a broad sampling of d. sicula from the mediterranean, and the potential species’ distribution was modelled to determine  whether intrapopulation and intraindividual nucleotide and haplotype diversity varies within the distribution;  whether the species is currently expanding.
with this data we expect to be able to answer questions on where the centre of origin of d. sicula’s present distribution is located, if the species has expanded recently assisted by humans and whether it is expected to continue its expansion.
methods
dna extraction and sequence amplification
fifty-eight populations were sampled from the mediterranean coast .
individuals from the boussadia and soukra populations  in tunisia reproduce sexually  <cit> .
individuals collected from crete , greece, also possess a copulatory apparatus.
approximately  <dig> individuals  were analysed from each location.
high-molecular-weight dna was purified from live or ethanol-fixed specimens using dnazol reagent .
a 711-bp fragment of the mitochondrial gene coi was amplified and sequenced with specific primers .
in some individuals, the sequence had to be obtained in two non-overlapping fragments since we were not able to obtain amplifications in a single fragment.
sequencing was performed using big dye  and an automated sequencer, abi prism  <dig> applied biosystems/hitachi  or at macrogen inc. .
the dna sequences were aligned in bioedit  <cit>  by sight based on the amino acid sequence and submitted to genbank .
all populations were identified as d. sicula by their coi sequence except populations marked with * that were previously identified by karyotype;  <dig> cohabited with d. liguriensis;  <dig> sexual populations with karyotype 2n = 18;  <dig> collected in the same pond where d. maghrebiana has been reported;  <dig> individuals with copulatory apparatus;  <dig> cohabited with d. elegans;  <dig> exfissiparous.
1
f forward, r reverse, *only for amplification, **only for sequencing  <dig> annealing temperature.
*only for amplification, **only for sequencing.
intraindividual polymorphism
to analyse the presence of multiple copies of the mitochondrial genome within individuals, some of the pcr products were cloned using an htp topo ta cloning kit for sequencing  following the manufacturers’ protocol.
we cloned the pcr products from sexual  and fissiparous individuals, within the latter group we included both specimens who showed double peaks in the sequences obtained directly from the pcr product  and some that did not show this pattern .
approximately ten colonies from each individual were amplified and sequenced using the t <dig> and t <dig> primers .
alignment was performed by sight based on the amino acid sequence, and the sequences were submitted to genbank .
we tested whether the number of substitutions observed in the cloned sequences agreed with the expected random accumulation of mutations promoted by the polymerase errors to determine whether the variation observed in the sequences was caused by errors in the pcr process.
therefore, we assumed that errors in the pcr process followed a poisson distribution, and the expected mean was computed using the polymerase error rate  obtained by savolainen et al.
<cit>  and following their procedure.
we rejected the random accumulation of substitutions due to the occurrence of polymerase errors when the observed number of mutations fell outside the 99% confidence interval computed for the poisson distribution .
in a functional protein coding sequence, more changes are expected in the 1st and 3rd codon positions.
to test the functionality of the sequences obtained in this study we calculated the expected number of changes per codon site in the case of a non-functional sequence  and tested whether the sequences obtained from the cloned pcr products fit the expectation, using a chi-square test.
nucleotide and haplotype diversity analyses
the dnasp v <dig>  program  <cit>  was used to determine the haplotype and nucleotide diversity for each population and the species using only individuals that were not polymorphic, since we could not know the exact sequence  of haplotypes for all the individuals presenting multiple polymorphic sites.
those populations for which only one sequence was available were not included.
haplotype and nucleotide diversity were also calculated for each of the cloned individuals.
two haplotype networks were constructed to study the relationships among haplotypes using tcs v <dig>   <cit> :  a haplotype network including only sequences from non-polymorphic individuals obtained after direct sequencing the pcr products and  a network containing only the cloned sequences from the polymorphic specimens.
distribution modelling
to study the possible expansion of d. sicula species, their potential distribution was modelled with maxent software v <dig> .3 k  <cit>  using the default settings, a random test of 25%, the threshold rule of  <dig> percentile training presence, and  <dig> replicates.
the independence of  <dig> climatic variables were tested based on the r <dig> statistic: annual mean temperature, mean diurnal range, isothermality, temperature seasonality, maximum temperature of warmest month, minimum temperature of coldest month, temperature annual range, mean temperature of wettest quarter, mean temperature of driest quarter, mean temperature of warmest quarter, mean temperature of coldest quarter, annual precipitation, precipitation of wettest month, precipitation of driest month, precipitation seasonality, precipitation of wettest quarter, precipitation of driest quarter, precipitation of warmest quarter and precipitation of coldest quarter.
finally four independent climatic variables were used: isothermality, mean temperature of the wettest quarter, mean temperature of the driest quarter and precipitation seasonality.
from  <dig> known localities ,  <dig> observation points were used, and  <dig> locations were discarded because of their geographic proximity .
RESULTS
population analyses of non-polymorphic individuals
in total, for the  <dig> locations sampled ,  <dig> sequences were obtained, and  <dig>  of these sequences showed double peaks in the chromatogram and were discarded for the population analyses.
because some populations had only one sequenced specimen  or all of their individuals exhibited polymorphic positions, only  <dig> individuals from  <dig> populations could be analysed.
only  <dig> populations showed more than one haplotype, and the majority of these showed two haplotypes.
the haplotype network was based on an alignment of  <dig> sequences and a length of  <dig> base pairs.
the total number of haplotypes for the coi gene in the populations was  <dig> .
the most frequent haplotype  was observed in  <dig> individuals from  <dig> populations throughout the distribution, and presented  <dig> differences regarding haplotype b .
haplotype b was present in  <dig> individuals from  <dig> populations, the majority of which were found northwest of the mediterranean.
haplotype c  was only observed in six individuals from southwestern populations, and presented six differences regarding haplotype a. the other  <dig> haplotypes were unique or observed in only two individuals.
haplotypes a and b were observed simultaneously in  <dig> populations from the centre of the mediterranean.
intra-individual polymorphism
between one and ten polymorphic positions  were observed in  <dig> individuals from  <dig> populations  when sequencing the pcr product directly.
twenty-one pcr products  from  <dig> populations were cloned, and approximately  <dig> clones per individual  were sequenced to determine the haplotypes of those individuals.
the final alignment had  <dig> sequences and a length of 537 bp.
versus
heteroplasmy test
s <dig>
1
s <dig>
1
s <dig>
1
expected number of point mutations for each individual based on a polymerase error rate of  <dig>  × 10- <dig> <cit> , point mutations found in cloned sequences, and probability that observed point mutations exceed the expected based on a poisson distribution.
sequence length was 537 bp, except in some cases marked with 1; n bp, total length amplified ; *significant differences based on 99% criterion.
in all cases we found multiple haplotypes within each individual even for those not showing polymorphic sites in the direct pcr sequence.
however, some of the point mutations detected in the cloned sequences could have been an artefact of the pcr and/or cloning processes.
savolainen et al.
<cit>  estimated this error rate to be  <dig>  × 10- <dig>  we compared the number of expected point mutations based on this rate with the number observed in each cloned specimen, taking into account the entire length of the amplified and sequenced dna .
all of the individuals with polymorphic sites exhibited significant differences between the two values , which indicated that the detected point mutations could not be accounted for by the addition of erroneous nucleotides by the polymerases during the pcr reactions.
none of the sexual individuals analysed had a significant result.
moreover, the cloning results showed that haplotypes a, b and c were present simultaneously in some of the specimens .
however, three of the specimens in which we had not detected polymorphic sites when sequencing directly from pcr, also showed significant differences between the expected and observed point mutations .
this result suggests that, in these individuals, despite that some point mutations could be caused by polymerase error, there are most likely low-frequency haplotypes that we cannot detect in the direct sequences using pcr.
in the subsequent analyses, we considered all of the individuals that showed a significant result to be heteroplasmic and the remaining individuals to be non-heteroplasmic.
for the haplotype network reconstruction , the single haplotypes sequenced from the non-heteroplasmic individuals were discarded since those must most probably be the result of a polymerase error as shown in the previous tests.
the final network contained  <dig> sequences and  <dig> different haplotypes were observed; haplotypes a, b and c were the most frequent and appeared  <dig>   <dig> and  <dig> times, respectively.
the remaining haplotypes were unique in the majority of the cases and showed a single difference from the a, b or c haplotypes.
all of the analysed individuals, except specimen s <dig>  , exhibited at least one of the most frequent haplotypes , and all of the individuals showed at least one unique haplotype .
the haplotype and nucleotide diversity, number of polymorphic sites and number of haplotypes were higher in heteroplasmic than non-heteroplasmic populations .
the values of π for heteroplasmic individuals  is slightly higher than that for all populations considering only non-polymorphic individuals , and in general intra-individual π values are as high or even more than the π value found for the whole species.
for h
d,
there is an increase when we only consider the heteroplasmic individuals due to the presence of multiple haplotypes at similar frequencies.
predicted distribution
the maxent results show the potential distribution of d. sicula from the mediterranean coasts, extending throughout the western region: northern tunisia, algeria and morocco, southern iberian peninsula, the italic peninsula, and islands of the balearic archipelago, sicily and sardinia.
all of the populations were observed within the predicted area, except for population s <dig>  and s <dig> , .
discussion
polymorphic sequences within individuals
heteroplasmy, the presence of different copies of the mitochondrial genome within the cells of an individual, has been reported for several taxa  <cit>  and was an expected outcome in the case of fissiparous flatworms, where somatic heteroplasmy would be inherited through clonal reproduction .
however, there are multiple hypotheses other than heteroplasmy that can explain polymorphic sequences; these explanations include the presence of numts or errors introduced by polymerases during the pcr reaction.
in the present case, the statistical analysis  indicated that the number of mutations observed in some individuals cannot be explained by errors introduced by the polymerases during pcr.
moreover, although the full evidence for a mitochondrial sequence to be functional and not present in the nucleus  can be reached only by sequencing mrnas  from the individuals  <cit> , in our case a variety of evidence makes us think that the diverse sequences found within single planarians are not numts.
numts typically have short sequences  and their mutation percentages are between  <dig> and 25%  <cit> .
the length of the shortest amplified fragments in this study was 228 bp, and the polymorphic positions observed in these fragments were also present when we amplified the larger fragment  .
the percentage of variation for our polymorphic sequences was below  <dig> %.
most of the substitutions were synonymous, and no stop codons were observed.
moreover, it is expected that the distribution of substitutions for a functional sequence shows a bias with more differences observed at the third and first position than at the second.
we tested the distribution of changes in the cloned sequences of both those individuals in which the previous test could not reject that the substitutions seen are due to taq error  and in those where the enzyme cannot explain the number of changes observed .
the result shows that in the latter there is a significant deviation from the 1:1: <dig> expected ratio .
however, some cases of recent numts have been described  <cit>  that still presented the characteristics of a functional gene .
there is still one more thing in our case that provides support for all the copies being functional, in the cloned individuals presenting heteroplasmy, we found copies of at least two or the three most frequent haplotypes  that are also found alone in some individuals and in whole populations, and that do not present stop codons, no traces of being non-functional .
although fission could easily explain the coexistence of low-frequency haplotypes that differ in one or a few substitutions in a given specimen , the origin of heteroplasmic individuals that simultaneously present haplotypes a, b and c is difficult to envision.
they can be a consequence of past sexual reproduction episodes, or a result of the appearance of the third haplotype in individuals that already presented the other two.
in any case, the finding of more unique haplotypes , differing mostly by one nucleotide from haplotypes a, b or c, in the heteroplasmic individuals supports the prediction that asexuals will tend to accumulate rare variants in all loci, which will be more accentuated in older asexuals .
on the other hand, as expected, heteroplasmy is not observed in sexual individuals.
although sexual individuals can also reproduce by fission, gametogenesis and zygote formation could represent a limitation for heteroplasmy propagation, which is a situation also observed in plants that reproduce both sexually and asexually  <cit> .
however, some specimens from the s <dig> population from crete possessed a copulatory apparatus and were consequently considered sexual individuals.
the observation that  <dig> out of the  <dig> sequenced individuals from this population were heteroplasmic could indicate that the individuals with a copulatory apparatus were actually ex-fissiparous or that the population is composed of a mix of sexual and fissiparous individuals.
origin and expansion
a recent study  <cit>  emphasised that d. sicula distribution could be a consequence of a relatively recent dispersion driven by human activities in accordance with a previous report  <cit> .
in the present study we found that in d. sicula populations, only two haplotypes showed high frequencies, differing at  <dig> sites over 624 bp, and most of the populations had a single haplotype .
differing from what was found in s. mediterranea <cit> , a non-overlapping geographical distribution for the different haplotypes, the most frequent d. sicula haplotypes are found practically throughout the distribution of the species.
this lack of structure seems to reinforce the idea of a recent expansion.
moreover, the haplotype and nucleotide diversity values are extremely low at the intra-population level and when the coi values for the entire set of populations were compared to those observed in other planarian species , such as the freshwater schmidtea mediterranea <cit>  or the terrestrial cephaloflexa bergi and imbira marcusi <cit> .
however, the presence of numerous heteroplasmic individuals bearing rare haplotypes in many fissiparous populations found in the present study  suggests the ancient age of these lineages and likely for the populations where we find multiple heteroplasmic specimens.
h
d
π
imbira marcusi
a
a this species has recently been renamed , in Álvarez-presas et al.
<cit>  appears as geoplana goetschi marcus,  <dig>
similar to other long-lived asexual lineages, d. sicula follows the expected accumulation of neutral or slightly deleterious substitutions over time and independent of one another, which results in high heterozygosity in their populations  <cit> .
however, in the case of d. sicula this accumulation is found within individuals, which does not allow us to exactly quantify the levels of intrapopulation variability .
however, taking into account our predictions of the expected levels of heteroplasmy and of heteroplasmic individuals in populations as related to their age, we can draw some conclusions on the origin and distribution of the species.
the high frequency and allocation of haplotypes a and b appears to indicate their antiquity in the species, particularly haplotype a, which is observed from morocco to greece .
in the northwestern populations , all  <dig> sequenced individuals showed haplotype b .
the low genetic variability in this region could indicate a relatively recent colonisation of this area by individuals deriving from geographically distant populations situated in the central or eastern regions .
haplotype c was observed only in western populations and is most likely derived from haplotype a based on the haplotype network .
haplotype c is observed in non-heteroplasmic individuals in morocco and the southern iberian peninsula and may have spread from there to the canary islands, where c was the only haplotype observed.
haplotype c was not present in any sexual specimen, but it was observed in heteroplasmic individuals from the western mediterranean islands and northern italy .
based on these data, we cannot determine with certainty where the haplotype originated from, but we can assume it had a relatively more recent origin than the other two haplotypes  and dispersed only through the western region.
however, we can speculate that the presence of haplotype c in heteroplasmic individuals in the central regions could indicate its origin is there, and its prevalence in the s <dig>  s <dig> and s <dig> populations could be the result of a founder effect during colonisation.
insufficient information provided by the few most common haplotype distributions and the fact that the presence of multiple single haplotypes in a population not necessarily are a consequence of an old age for the population at that locality, make it difficult to determine where colonisation began in the mediterranean.
however, some facts can help:  the evolutionary sister group of this species, d. aethiopica <cit>  is distributed exclusively in africa  <cit> , and their split is dated around  <dig> my   <cit> , well after the mediterranean basin was formed and all islands had reached their present location  <cit> ,  most of the populations  in which haplotypes a and b cohabit are in the central region of the mediterranean basin ,  known sexual populations originate only in either sicily  or tunisia and algeria.
collectively, these data could indicate that the central mediterranean is the region from where dispersal took place.
however, the question as to how this species attained such a wide distribution with such a low dispersal capability  <cit>  remains.
various factors could have contributed, including their fissiparous reproduction and their capability to resist high water temperatures.
it has been shown  <cit>  that asexual reproduction in planarians allows a more rapid increase in population size under ecological conditions in which food is limited.
moreover, a study by vila-farré  <cit>  have shown that d. sicula is sensitive to low temperatures and to abrupt temperature changes, while charni et al.
<cit>  found that the diverse tunisian populations lived in water with temperatures ranging from  <dig> to 22°c.
it has also been observed that during spring and autumn, higher rates of fissioning occur  <cit> .
therefore, they can resist high temperatures of water, and can reproduce and increase their populations very fast.
but could the species by itself reach its present distribution area?
given the maximum age of the species situated around  <dig> my   <cit> , when the mediterranean basin was already formed and its islands in its present situation, it seems difficult to explain how the animals could have moved between islands and the continent  even in one direction or the other.
if the origin of the species preceded the messinian salinity crisis or else as a consequence of the sea level changes during the pleistocene glaciations, we could maybe hypothesize the contact among some fluvial basins in the central mediterranean that would allow its expansion there, but still this would not explain its presence in places situated as far away as the greek islands, the canary islands or how haplotype b could have reached the french coast.
these latter cases provide support for a mixed history of self-dispersal and human introductions.
de vries  <cit>  already proposed that given the antiquity of human trade in the mediterranean basin there may have been multiple opportunities for these animals to be transported around, including modern commercial trade of fish, plants or timber, as has been proposed for other planarian species  <cit> , which also could explain this expansion.
our hypothesis is then that d. sicula most likely originated in northern africa and dispersed through the central mediterranean region.
after one or multiple populations became triploid and fissiparous, they colonized the whole mediterranean basin either by their own means or assisted by human activities .
species distribution and future expansion
under the hypothesis that d. sicula is a relatively recent coloniser, it was advisable to test whether this species had reached its current range of fitted distribution or was still expanding.
the potential distribution predicted with the climatic models for d. sicula , fits well with what has been observed, mainly occupying coastal regions.
it has been reported that d. sicula does not cope well with abrupt changes in temperature  <cit>  and is typically observed near the coast, where the water temperature is highest in the winter and the thermal changes are less abrupt than in the high mountains.
all known populations, except for s <dig> and s <dig>  have been observed within the predicted area; however, population s <dig> inhabited an oasis in the tunisian desert, the special conditions in this reduced area could be ideal for the species although precision of climate layers used in the analysis do not allow for detecting them.
population s <dig> in the canary islands is found within garajonay national park, in a area of laurel forest, thought to be the remnant of subtropical woods that covered the mediterranean area during the tertiary period, which likely provides d. sicula the conditions needed to survive, although the models, again as a consequence of the scale, may have not been able to detect it.
some regions that are potentially suitable for d. sicula have not yet been widely sampled, such as northern africa , the southwestern iberian peninsula and parts of the italic peninsula.
however, restriction pattern analyses showed similarities among some populations north of the italic peninsula and d. sicula <cit> , although the italian populations could not be identified; moreover, benazzi and deri  <cit>  reported populations in tuscany , ponza island , rome  and calabria that could be attributed to d. sicula based on morphological features, and some populations with a chromosomal complement of 27+ 2b chromosomes have been reported in calabria by deri et al.
<cit> .
in other areas, d. sicula has never been reported  or is relegated to the periphery , although these regions appear to fulfil the climatic requirements of the species.
however, several populations of d. benazzii have been observed in these regions.
the competition between d. sicula and d. benazzii may explain why d. sicula is absent in corsica and has a low presence in sardinia.
however, d. sicula is present together with other planarian species in localities such as lebna in tunisia , which points to the possibility that historical reasons could also explain its absence in some regions .
for the eastern mediterranean, some populations described in israel as type  <dig> by bromley  <cit>  showed karyotypes and morphological characteristics similar to those of d. sicula.
although type  <dig> was related to d. biblica in a study by bromley  <cit> , these two species cannot be definitively differentiated  <cit> .
an accurate study of these two species is necessary to clarify their status and confirm whether d. sicula is also present in that region or whether the presence of a closely related species potentially precluded its colonisation.
recent exhaustive samplings performed in the iberian peninsula in search of dugesia subtentaculata , allows us to confidently state that d. sicula is only observed on the east and south coast, as the prediction analysis suggests.
in this area, d. sicula cohabits with other freshwater planarian species such as s. mediterranea.
in recent years, populations of s. mediterranea <cit>  appeared to have vanished, whereas populations of d. sicula became more frequent.
additionally, in the greek population from rhodes , d. sicula cohabits with the autochthonous d. elegans, and only a few animals belong to this latter species.
this result could indicate that d. sicula can have an adverse effect on the native triclad species, which contrasts with what appears to occur in sardinia.
a similar situation has been reported in wales, where the presence of two immigrant species, planaria torva and girardia tigrina has been followed via three surveys over 50 years  <cit> .
the most recent report  <cit>  indicates that the invasion of g. tigrina has been more successful, which is most likely because of its asexual reproduction and opportunistic feeding characteristics, although other factors, such as its capability to adhere to surfaces better than other triclads, have also been considered.
moreover, g. tigrina has been shown to have a negative effect on native species of triclads, which is most likely explained by food competition.
CONCLUSIONS
in conclusion, it appears that d. sicula has reached a large proportion of the area of its potentially favourable distribution in the mediterranean basin, being a remarkable case of a broad colonisation, in extreme contrast with the rest of mediterranean dugesia species, with all of them being endemic or with very restricted distributions.
d. sicula expansion is now limited to spreading to new freshwater basins within the areas it currently inhabits.
however, future changes increasing the temperature, such as those predicted by climate change hypothesis, could expand its fitted area to more northern and interior areas.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mr and eml designed the study.
eml performed the molecular work and obtained the sequence data.
eml and mr conducted all of the analyses and wrote the manuscript.
both authors read and approved the final manuscript.
supplementary material
additional file 1: figure s1
schematic representation of the spreading of a new mutation in fissiparous organisms.
the yellow and orange colours represent the tissues bearing the new mutations.
the arrows represent fission  and regeneration  cycles.
figure s <dig>  the haplotype network based on cloned sequences.
each circle represents a haplotype observed in the sample.
the black dots represent intermediate  haplotypes, the lines connecting the haplotypes  represent one nucleotide change, and the size of each circle is proportional to the haplotype frequency in the sample.
the colours in the network correspond to those on the map and represent the studied locations.
the crossed-out names of the individuals represent non-heteroplasmic specimens, and their sequences were not included in the network.
click here for file
additional file 2: table s1
genetic diversity in the populations.
table s <dig>  genetic diversity in the cloned individuals.
table s <dig>  chi-square test for the distribution of substitutions along the three codon positions .
click here for file
acknowledgements
this research was supported by dgi-mec grant cgl2005-00371/bos, cgl2008-00378/bos and cgl2011- <dig> to m.r.
we thank p. aguilera, n. bonada, m. charni, k. gritzalis, a. h. harrath, c. hernando, l. leria, c. maritur, e. mateos, m. pala, s. pons, i. ribera, r. sluys, e. solà, j. solà, g. a. stocchino, m. vila-farré and m. yacoubi for their collaboration on sample collection.
g. a. stocchino also for sharing information on the likely presence of d. sicula in itlay.
j. rozas for his advice on the analytical methodologies.
we also thank the serveis cientificotècnics .
jim hesson of academicenglishsolutions.com proofread the english.

 dugesia sicula lepori, is a freshwater planarian typically found in ponds, streams and springs close to the mediterranean coast. the first record of d. sicula described a population comprising both sexual and asexual individuals in catania, sicily. fissiparity is an asexual form of reproduction that involves transversal division of the individuals into two fragments and subsequent regeneration of the absent structures. fissiparous planarians do not develop a copulatory apparatus, which precludes the assignment of these populations to any species since this structure provides most of the defining characters for species description in planarians. however, the morphological characters of the copulatory apparatus of the sexual specimens allowed the species to be described. d. sicula was subsequently found on the islands of elba and mallorca. none of the originally described sexual populations can be currently found; but, some sexual populations in northern africa have been recently reported. moreover, a molecular-based study showed that most of the numerous fissiparous populations of genus dugesia present in the western mediterranean belong to d. sicula. under certain conditions triploid fissiparous planarians can develop hyperplasic ovaries and a copulatory apparatus, these individuals are known as ex-fissiparous. in d. sicula, it has been observed that some ex-fissiparous individuals laid cocoons but remained sterile, which indicated that asexual populations have lost the capability to sexually reproduce. the literature suggests that d. sicula is widely distributed in the mediterranean area from morocco to greece, and our own observations show that this species mostly comprise fissiparous populations, most likely sexually sterile, distributed throughout the mediterranean, as well as the canary islands, being the only species of this genus in the mediterranean with a such a wide distribution. contrasting with this wide distribution the species presents extremely low variability for the coi gene among populations located hundreds of kilometres apart. these two features could be a consequence of a recent expansion promoted by human activities, as has been shown for other organisms. alternatively, the low levels of variability could be exclusively a consequence of their fissiparous reproduction, as described in schmidtea mediterranea and other organisms, and its wide distribution may be a consequence of its active spread over the years. in both cases, the hypotheses contrast with the poor dispersion capability assumed for freshwater planarians. freshwater planarians do not exhibit larval dispersal stages or forms resistant to desiccation; these individuals thus require contiguous freshwater bodies to survive and disperse. however, introductions due to human activities have been documented in other freshwater planarians, such as girardia tigrina and schmidtea polychroa and more recently have also been proposed for other dugesia species, but in all these cases the introduced animals either have a restricted distribution or its introduction and progression in the new areas has been followed by scientist. fissiparous reproduction can have effects on the distribution of genetic variability within the populations and even on individuals that could perhaps help disentangle this situation. regeneration after fission is driven by neoblasts, which are pluripotent stem cells. a new mutation in a neoblast would expand in the individual when this neoblast generates new tissues after the fission of the animal, thus increasing frequency of the mutation in the individual. following replication via several fission cycles, the new mutation could disappear or spread in the population of cells within the individual and the organism’s population by genetic drift. in this manner, we would expect to observe cells with different sequences within one individual. this could be considered a special case of heteroplasmy that differs from cases in which heteroplasmy is transmitted sexually or parthenogenetically through the mitochondria present in the oocyte where all cells from a new individual will present the same haplotypes. heteroplasmy has been detected in many species, and some studies in humans have shown it increases with age in somatic cells and is more frequent in certain tissues, such as muscle tissue, most likely because of increases in the number of mutations caused by the presence of oxidative radicals. heteroplasmy has also been detected in coi sequences of fissiparous populations of d. japonica. due to their characteristics, in fissiparous planarians the somatic heteroplasmy would be maintained in populations and increased over time. if this situation holds in fissiparous populations of d. sicula, the analysis of their genetic variability could help in discriminating between a recent or ancient origin of the species populations and hence provide clues to the origin of its distribution. the absence or low frequency of heteroplasmic specimens in a specific locality of fissiparous reproducing individuals could be interpreted as a consequence of a recent bottleneck, potentially because of recent colonisation by a few individuals or recovery from a recent population-wide decrease. the alternate condition, the presence of many heteroplasmic individuals, particularly those exhibiting a high number of different low-frequency haplotypes, could then potentially be related to aged individuals and populations, although a recent colonisation by highly heteroplasmic individuals could also explain this outcome. following the same reasoning, one would not expect to find heteroplasmy, or different copies of mitochondrial genes within an individual in sexual organisms. in the present study, we analysed the genetic diversity and structure of a broad sampling of d. sicula from the mediterranean, and the potential species’ distribution was modelled to determine whether intrapopulation and intraindividual nucleotide and haplotype diversity varies within the distribution; whether the species is currently expanding. with this data we expect to be able to answer questions on where the centre of origin of d. sicula’s present distribution is located, if the species has expanded recently assisted by humans and whether it is expected to continue its expansion. methods dna extraction and sequence amplification fifty-eight populations were sampled from the mediterranean coast. individuals from the boussadia and soukra populations in tunisia reproduce sexually. individuals collected from crete, greece, also possess a copulatory apparatus. approximately individuals were analysed from each location. high-molecular-weight dna was purified from live or ethanol-fixed specimens using dnazol reagent. a 711-bp fragment of the mitochondrial gene coi was amplified and sequenced with specific primers. in some individuals, the sequence had to be obtained in two non-overlapping fragments since we were not able to obtain amplifications in a single fragment. sequencing was performed using big dye and an automated sequencer, abi prism applied biosystems/hitachi or at macrogen inc.. the dna sequences were aligned in bioedit by sight based on the amino acid sequence and submitted to genbank. all populations were identified as d. sicula by their coi sequence except populations marked with * that were previously identified by karyotype; cohabited with d. liguriensis; sexual populations with karyotype 2n = 18; collected in the same pond where d. maghrebiana has been reported; individuals with copulatory apparatus; cohabited with d. elegans; exfissiparous. 1 f forward, r reverse, *only for amplification, **only for sequencing annealing temperature. *only for amplification, **only for sequencing. intraindividual polymorphism to analyse the presence of multiple copies of the mitochondrial genome within individuals, some of the pcr products were cloned using an htp topo ta cloning kit for sequencing following the manufacturers’ protocol. we cloned the pcr products from sexual and fissiparous individuals, within the latter group we included both specimens who showed double peaks in the sequences obtained directly from the pcr product and some that did not show this pattern. approximately ten colonies from each individual were amplified and sequenced using the t and t primers. alignment was performed by sight based on the amino acid sequence, and the sequences were submitted to genbank. we tested whether the number of substitutions observed in the cloned sequences agreed with the expected random accumulation of mutations promoted by the polymerase errors to determine whether the variation observed in the sequences was caused by errors in the pcr process. therefore, we assumed that errors in the pcr process followed a poisson distribution, and the expected mean was computed using the polymerase error rate obtained by savolainen et al. and following their procedure. we rejected the random accumulation of substitutions due to the occurrence of polymerase errors when the observed number of mutations fell outside the 99% confidence interval computed for the poisson distribution. in a functional protein coding sequence, more changes are expected in the 1st and 3rd codon positions. to test the functionality of the sequences obtained in this study we calculated the expected number of changes per codon site in the case of a non-functional sequence and tested whether the sequences obtained from the cloned pcr products fit the expectation, using a chi-square test. nucleotide and haplotype diversity analyses the dnasp v program was used to determine the haplotype and nucleotide diversity for each population and the species using only individuals that were not polymorphic, since we could not know the exact sequence of haplotypes for all the individuals presenting multiple polymorphic sites. those populations for which only one sequence was available were not included. haplotype and nucleotide diversity were also calculated for each of the cloned individuals. two haplotype networks were constructed to study the relationships among haplotypes using tcs v: a haplotype network including only sequences from non-polymorphic individuals obtained after direct sequencing the pcr products and a network containing only the cloned sequences from the polymorphic specimens. distribution modelling to study the possible expansion of d. sicula species, their potential distribution was modelled with maxent software v.3 k using the default settings, a random test of 25%, the threshold rule of percentile training presence, and replicates. the independence of climatic variables were tested based on the r statistic: annual mean temperature, mean diurnal range, isothermality, temperature seasonality, maximum temperature of warmest month, minimum temperature of coldest month, temperature annual range, mean temperature of wettest quarter, mean temperature of driest quarter, mean temperature of warmest quarter, mean temperature of coldest quarter, annual precipitation, precipitation of wettest month, precipitation of driest month, precipitation seasonality, precipitation of wettest quarter, precipitation of driest quarter, precipitation of warmest quarter and precipitation of coldest quarter. finally four independent climatic variables were used: isothermality, mean temperature of the wettest quarter, mean temperature of the driest quarter and precipitation seasonality. from known localities, observation points were used, and locations were discarded because of their geographic proximity.  population analyses of non-polymorphic individuals in total, for the locations sampled, sequences were obtained, and of these sequences showed double peaks in the chromatogram and were discarded for the population analyses. because some populations had only one sequenced specimen or all of their individuals exhibited polymorphic positions, only individuals from populations could be analysed. only populations showed more than one haplotype, and the majority of these showed two haplotypes. the haplotype network was based on an alignment of sequences and a length of base pairs. the total number of haplotypes for the coi gene in the populations was. the most frequent haplotype was observed in individuals from populations throughout the distribution, and presented differences regarding haplotype b. haplotype b was present in individuals from populations, the majority of which were found northwest of the mediterranean. haplotype c was only observed in six individuals from southwestern populations, and presented six differences regarding haplotype a. the other haplotypes were unique or observed in only two individuals. haplotypes a and b were observed simultaneously in populations from the centre of the mediterranean. intra-individual polymorphism between one and ten polymorphic positions were observed in individuals from populations when sequencing the pcr product directly. twenty-one pcr products from populations were cloned, and approximately clones per individual were sequenced to determine the haplotypes of those individuals. the final alignment had sequences and a length of 537 bp. versus heteroplasmy test s 1 s 1 s 1 expected number of point mutations for each individual based on a polymerase error rate of × 10-, point mutations found in cloned sequences, and probability that observed point mutations exceed the expected based on a poisson distribution. sequence length was 537 bp, except in some cases marked with 1; n bp, total length amplified; *significant differences based on 99% criterion. in all cases we found multiple haplotypes within each individual even for those not showing polymorphic sites in the direct pcr sequence. however, some of the point mutations detected in the cloned sequences could have been an artefact of the pcr and/or cloning processes. savolainen et al. estimated this error rate to be × 10- we compared the number of expected point mutations based on this rate with the number observed in each cloned specimen, taking into account the entire length of the amplified and sequenced dna. all of the individuals with polymorphic sites exhibited significant differences between the two values, which indicated that the detected point mutations could not be accounted for by the addition of erroneous nucleotides by the polymerases during the pcr reactions. none of the sexual individuals analysed had a significant result. moreover, the cloning results showed that haplotypes a, b and c were present simultaneously in some of the specimens. however, three of the specimens in which we had not detected polymorphic sites when sequencing directly from pcr, also showed significant differences between the expected and observed point mutations. this result suggests that, in these individuals, despite that some point mutations could be caused by polymerase error, there are most likely low-frequency haplotypes that we cannot detect in the direct sequences using pcr. in the subsequent analyses, we considered all of the individuals that showed a significant result to be heteroplasmic and the remaining individuals to be non-heteroplasmic. for the haplotype network reconstruction, the single haplotypes sequenced from the non-heteroplasmic individuals were discarded since those must most probably be the result of a polymerase error as shown in the previous tests. the final network contained sequences and different haplotypes were observed; haplotypes a, b and c were the most frequent and appeared and times, respectively. the remaining haplotypes were unique in the majority of the cases and showed a single difference from the a, b or c haplotypes. all of the analysed individuals, except specimen s, exhibited at least one of the most frequent haplotypes, and all of the individuals showed at least one unique haplotype. the haplotype and nucleotide diversity, number of polymorphic sites and number of haplotypes were higher in heteroplasmic than non-heteroplasmic populations. the values of π for heteroplasmic individuals is slightly higher than that for all populations considering only non-polymorphic individuals, and in general intra-individual π values are as high or even more than the π value found for the whole species. for h d, there is an increase when we only consider the heteroplasmic individuals due to the presence of multiple haplotypes at similar frequencies. predicted distribution the maxent results show the potential distribution of d. sicula from the mediterranean coasts, extending throughout the western region: northern tunisia, algeria and morocco, southern iberian peninsula, the italic peninsula, and islands of the balearic archipelago, sicily and sardinia. all of the populations were observed within the predicted area, except for population s and s,. discussion polymorphic sequences within individuals heteroplasmy, the presence of different copies of the mitochondrial genome within the cells of an individual, has been reported for several taxa and was an expected outcome in the case of fissiparous flatworms, where somatic heteroplasmy would be inherited through clonal reproduction. however, there are multiple hypotheses other than heteroplasmy that can explain polymorphic sequences; these explanations include the presence of numts or errors introduced by polymerases during the pcr reaction. in the present case, the statistical analysis indicated that the number of mutations observed in some individuals cannot be explained by errors introduced by the polymerases during pcr. moreover, although the full evidence for a mitochondrial sequence to be functional and not present in the nucleus can be reached only by sequencing mrnas from the individuals, in our case a variety of evidence makes us think that the diverse sequences found within single planarians are not numts. numts typically have short sequences and their mutation percentages are between and 25%. the length of the shortest amplified fragments in this study was 228 bp, and the polymorphic positions observed in these fragments were also present when we amplified the larger fragment. the percentage of variation for our polymorphic sequences was below %. most of the substitutions were synonymous, and no stop codons were observed. moreover, it is expected that the distribution of substitutions for a functional sequence shows a bias with more differences observed at the third and first position than at the second. we tested the distribution of changes in the cloned sequences of both those individuals in which the previous test could not reject that the substitutions seen are due to taq error and in those where the enzyme cannot explain the number of changes observed. the result shows that in the latter there is a significant deviation from the 1:1: expected ratio. however, some cases of recent numts have been described that still presented the characteristics of a functional gene. there is still one more thing in our case that provides support for all the copies being functional, in the cloned individuals presenting heteroplasmy, we found copies of at least two or the three most frequent haplotypes that are also found alone in some individuals and in whole populations, and that do not present stop codons, no traces of being non-functional. although fission could easily explain the coexistence of low-frequency haplotypes that differ in one or a few substitutions in a given specimen, the origin of heteroplasmic individuals that simultaneously present haplotypes a, b and c is difficult to envision. they can be a consequence of past sexual reproduction episodes, or a result of the appearance of the third haplotype in individuals that already presented the other two. in any case, the finding of more unique haplotypes, differing mostly by one nucleotide from haplotypes a, b or c, in the heteroplasmic individuals supports the prediction that asexuals will tend to accumulate rare variants in all loci, which will be more accentuated in older asexuals. on the other hand, as expected, heteroplasmy is not observed in sexual individuals. although sexual individuals can also reproduce by fission, gametogenesis and zygote formation could represent a limitation for heteroplasmy propagation, which is a situation also observed in plants that reproduce both sexually and asexually. however, some specimens from the s population from crete possessed a copulatory apparatus and were consequently considered sexual individuals. the observation that out of the sequenced individuals from this population were heteroplasmic could indicate that the individuals with a copulatory apparatus were actually ex-fissiparous or that the population is composed of a mix of sexual and fissiparous individuals. origin and expansion a recent study emphasised that d. sicula distribution could be a consequence of a relatively recent dispersion driven by human activities in accordance with a previous report. in the present study we found that in d. sicula populations, only two haplotypes showed high frequencies, differing at sites over 624 bp, and most of the populations had a single haplotype. differing from what was found in s. mediterranea, a non-overlapping geographical distribution for the different haplotypes, the most frequent d. sicula haplotypes are found practically throughout the distribution of the species. this lack of structure seems to reinforce the idea of a recent expansion. moreover, the haplotype and nucleotide diversity values are extremely low at the intra-population level and when the coi values for the entire set of populations were compared to those observed in other planarian species, such as the freshwater schmidtea mediterranea or the terrestrial cephaloflexa bergi and imbira marcusi. however, the presence of numerous heteroplasmic individuals bearing rare haplotypes in many fissiparous populations found in the present study suggests the ancient age of these lineages and likely for the populations where we find multiple heteroplasmic specimens. h d π imbira marcusi a a this species has recently been renamed, in Álvarez-presas et al. appears as geoplana goetschi marcus, similar to other long-lived asexual lineages, d. sicula follows the expected accumulation of neutral or slightly deleterious substitutions over time and independent of one another, which results in high heterozygosity in their populations. however, in the case of d. sicula this accumulation is found within individuals, which does not allow us to exactly quantify the levels of intrapopulation variability. however, taking into account our predictions of the expected levels of heteroplasmy and of heteroplasmic individuals in populations as related to their age, we can draw some conclusions on the origin and distribution of the species. the high frequency and allocation of haplotypes a and b appears to indicate their antiquity in the species, particularly haplotype a, which is observed from morocco to greece. in the northwestern populations, all sequenced individuals showed haplotype b. the low genetic variability in this region could indicate a relatively recent colonisation of this area by individuals deriving from geographically distant populations situated in the central or eastern regions. haplotype c was observed only in western populations and is most likely derived from haplotype a based on the haplotype network. haplotype c is observed in non-heteroplasmic individuals in morocco and the southern iberian peninsula and may have spread from there to the canary islands, where c was the only haplotype observed. haplotype c was not present in any sexual specimen, but it was observed in heteroplasmic individuals from the western mediterranean islands and northern italy. based on these data, we cannot determine with certainty where the haplotype originated from, but we can assume it had a relatively more recent origin than the other two haplotypes and dispersed only through the western region. however, we can speculate that the presence of haplotype c in heteroplasmic individuals in the central regions could indicate its origin is there, and its prevalence in the s s and s populations could be the result of a founder effect during colonisation. insufficient information provided by the few most common haplotype distributions and the fact that the presence of multiple single haplotypes in a population not necessarily are a consequence of an old age for the population at that locality, make it difficult to determine where colonisation began in the mediterranean. however, some facts can help: the evolutionary sister group of this species, d. aethiopica is distributed exclusively in africa, and their split is dated around my, well after the mediterranean basin was formed and all islands had reached their present location, most of the populations in which haplotypes a and b cohabit are in the central region of the mediterranean basin, known sexual populations originate only in either sicily or tunisia and algeria. collectively, these data could indicate that the central mediterranean is the region from where dispersal took place. however, the question as to how this species attained such a wide distribution with such a low dispersal capability remains. various factors could have contributed, including their fissiparous reproduction and their capability to resist high water temperatures. it has been shown that asexual reproduction in planarians allows a more rapid increase in population size under ecological conditions in which food is limited. moreover, a study by vila-farré have shown that d. sicula is sensitive to low temperatures and to abrupt temperature changes, while charni et al. found that the diverse tunisian populations lived in water with temperatures ranging from to 22°c. it has also been observed that during spring and autumn, higher rates of fissioning occur. therefore, they can resist high temperatures of water, and can reproduce and increase their populations very fast. but could the species by itself reach its present distribution area? given the maximum age of the species situated around my, when the mediterranean basin was already formed and its islands in its present situation, it seems difficult to explain how the animals could have moved between islands and the continent even in one direction or the other. if the origin of the species preceded the messinian salinity crisis or else as a consequence of the sea level changes during the pleistocene glaciations, we could maybe hypothesize the contact among some fluvial basins in the central mediterranean that would allow its expansion there, but still this would not explain its presence in places situated as far away as the greek islands, the canary islands or how haplotype b could have reached the french coast. these latter cases provide support for a mixed history of self-dispersal and human introductions. de vries already proposed that given the antiquity of human trade in the mediterranean basin there may have been multiple opportunities for these animals to be transported around, including modern commercial trade of fish, plants or timber, as has been proposed for other planarian species, which also could explain this expansion. our hypothesis is then that d. sicula most likely originated in northern africa and dispersed through the central mediterranean region. after one or multiple populations became triploid and fissiparous, they colonized the whole mediterranean basin either by their own means or assisted by human activities. species distribution and future expansion under the hypothesis that d. sicula is a relatively recent coloniser, it was advisable to test whether this species had reached its current range of fitted distribution or was still expanding. the potential distribution predicted with the climatic models for d. sicula, fits well with what has been observed, mainly occupying coastal regions. it has been reported that d. sicula does not cope well with abrupt changes in temperature and is typically observed near the coast, where the water temperature is highest in the winter and the thermal changes are less abrupt than in the high mountains. all known populations, except for s and s have been observed within the predicted area; however, population s inhabited an oasis in the tunisian desert, the special conditions in this reduced area could be ideal for the species although precision of climate layers used in the analysis do not allow for detecting them. population s in the canary islands is found within garajonay national park, in a area of laurel forest, thought to be the remnant of subtropical woods that covered the mediterranean area during the tertiary period, which likely provides d. sicula the conditions needed to survive, although the models, again as a consequence of the scale, may have not been able to detect it. some regions that are potentially suitable for d. sicula have not yet been widely sampled, such as northern africa, the southwestern iberian peninsula and parts of the italic peninsula. however, restriction pattern analyses showed similarities among some populations north of the italic peninsula and d. sicula, although the italian populations could not be identified; moreover, benazzi and deri reported populations in tuscany, ponza island, rome and calabria that could be attributed to d. sicula based on morphological features, and some populations with a chromosomal complement of 27+ 2b chromosomes have been reported in calabria by deri et al.. in other areas, d. sicula has never been reported or is relegated to the periphery, although these regions appear to fulfil the climatic requirements of the species. however, several populations of d. benazzii have been observed in these regions. the competition between d. sicula and d. benazzii may explain why d. sicula is absent in corsica and has a low presence in sardinia. however, d. sicula is present together with other planarian species in localities such as lebna in tunisia, which points to the possibility that historical reasons could also explain its absence in some regions. for the eastern mediterranean, some populations described in israel as type by bromley showed karyotypes and morphological characteristics similar to those of d. sicula. although type was related to d. biblica in a study by bromley, these two species cannot be definitively differentiated. an accurate study of these two species is necessary to clarify their status and confirm whether d. sicula is also present in that region or whether the presence of a closely related species potentially precluded its colonisation. recent exhaustive samplings performed in the iberian peninsula in search of dugesia subtentaculata, allows us to confidently state that d. sicula is only observed on the east and south coast, as the prediction analysis suggests. in this area, d. sicula cohabits with other freshwater planarian species such as s. mediterranea. in recent years, populations of s. mediterranea appeared to have vanished, whereas populations of d. sicula became more frequent. additionally, in the greek population from rhodes, d. sicula cohabits with the autochthonous d. elegans, and only a few animals belong to this latter species. this result could indicate that d. sicula can have an adverse effect on the native triclad species, which contrasts with what appears to occur in sardinia. a similar situation has been reported in wales, where the presence of two immigrant species, planaria torva and girardia tigrina has been followed via three surveys over 50 years. the most recent report indicates that the invasion of g. tigrina has been more successful, which is most likely because of its asexual reproduction and opportunistic feeding characteristics, although other factors, such as its capability to adhere to surfaces better than other triclads, have also been considered. moreover, g. tigrina has been shown to have a negative effect on native species of triclads, which is most likely explained by food competition.  in conclusion, it appears that d. sicula has reached a large proportion of the area of its potentially favourable distribution in the mediterranean basin, being a remarkable case of a broad colonisation, in extreme contrast with the rest of mediterranean dugesia species, with all of them being endemic or with very restricted distributions. d. sicula expansion is now limited to spreading to new freshwater basins within the areas it currently inhabits. however, future changes increasing the temperature, such as those predicted by climate change hypothesis, could expand its fitted area to more northern and interior areas. competing interests the authors declare that they have no competing interests. authors’ contributions mr and eml designed the study. eml performed the molecular work and obtained the sequence data. eml and mr conducted all of the analyses and wrote the manuscript. both authors read and approved the final manuscript. supplementary material additional file 1: figure s1 schematic representation of the spreading of a new mutation in fissiparous organisms. the yellow and orange colours represent the tissues bearing the new mutations. the arrows represent fission and regeneration cycles. figure s the haplotype network based on cloned sequences. each circle represents a haplotype observed in the sample. the black dots represent intermediate haplotypes, the lines connecting the haplotypes represent one nucleotide change, and the size of each circle is proportional to the haplotype frequency in the sample. the colours in the network correspond to those on the map and represent the studied locations. the crossed-out names of the individuals represent non-heteroplasmic specimens, and their sequences were not included in the network. click here for file additional file 2: table s1 genetic diversity in the populations. table s genetic diversity in the cloned individuals. table s chi-square test for the distribution of substitutions along the three codon positions. click here for file acknowledgements this research was supported by dgi-mec grant cgl2005-00371/bos, cgl2008-00378/bos and cgl2011- to m.r. we thank p. aguilera, n. bonada, m. charni, k. gritzalis, a. h. harrath, c. hernando, l. leria, c. maritur, e. mateos, m. pala, s. pons, i. ribera, r. sluys, e. solà, j. solà, g. a. stocchino, m. vila-farré and m. yacoubi for their collaboration on sample collection. g. a. stocchino also for sharing information on the likely presence of d. sicula in itlay. j. rozas for his advice on the analytical methodologies. we also thank the serveis cientificotècnics. jim hesson of academicenglishsolutions.com proofread the english.
the distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: the species most probably originated in africa and dispersed through the central mediterranean. after one or multiple populations became triploid and fissiparous, the species colonized the mediterranean basin, likely both by its own means and helped by human activities. it mostly comprises fissiparous populations, most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the mediterranean area. the species presents most frequent haplotypes, differing in a maximum of base pairs. we have estimated haplotype and nucleotide diversity of cytochrome oxidase i within and among individuals along the species distribution of a broad sample of d. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change. the distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central mediterranean region. the potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data. its present distribution practically fulfils its potential distribution as modelled with climatic variables. however, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones.
 dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the mediterranean area. it mostly comprises fissiparous populations, most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. however, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones. we have estimated haplotype and nucleotide diversity of cytochrome oxidase i within and among individuals along the species distribution of a broad sample of d. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. our aim was to determine the centre of colonisation origin, whether the populations are recent, and whether the species is expanding.  the species presents most frequent haplotypes, differing in a maximum of base pairs. as expected from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple heteroplasmic haplotypes. the distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central mediterranean region. the potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data.  the distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: the species most probably originated in africa and dispersed through the central mediterranean. after one or multiple populations became triploid and fissiparous, the species colonized the mediterranean basin, likely both by its own means and helped by human activities. its present distribution practically fulfils its potential distribution as modelled with climatic variables. its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change. heteroplasmyfissiparityasexualrecent colonizationmediterranean basintrade routesniche modelling 