Molecular phylogeny of tubificid oligochaetes with special emphasis on Tubifex tubifex (Tubificidae)

Tubifex tubifex is a cosmopolitan freshwater oligochaete whose presence has been studied as a health indicator of the aquatic environment and as a host for several myxozoan parasites of fish. Unfortunately, current morphological criteria used to distinguish Tubifex spp. (Tubificidae) are inadequate. We therefore developed mitochondrial 16S ribosomal DNA markers to examine phylogenetic relationships among aquatic oligochaetes and to distinguish species of Tubifex that might serve as hosts for a particular myxozoan parasite, Myxobolus cerebralis. Our phylogenetic analyses of oligochaetes based on a 378-bp segment yielded one most parsimonious tree with three major groups that corresponded to the families Lumbricidae, Sparganophilidae, and Tubificidae. T. tubifex and T. ignotus formed a monophyletic assemblage, and a sister relationship between the genera Tubifex and Limnodrilus was strongly supported. A second analysis of the relationship within the genus Tubifex identified six genetically distinct lineages of T. tubifex from North America and Europe that were separated by genetic distances comparable to those found for "well-defined" species of Limnodrilus. Therefore, the existence of several morphologically indistinguishable, thus cryptic, species of Tubifex in North America and Europe is suggested.     

The world in a river? A preliminary analysis of the 16S rDNA variability of Tubifex species (Clitellata: Tubificidae) from the Lambro River

Tubifex tubifex Müller, 1774 is a cosmopolitan freshwater tubificid widely used as a model in ecotoxicology, population dynamics and developmental biology. It is traditionally recognized as a polytypic species and in Lambro River (Milano, Northern Italy) it occurs in two of the three recognized forms, named "tubifex" and "blanchardi", alternatively considered as ecological forms or distinct species. To investigate the genetic differentiation of the populations occurring in the Lambro River we sequenced a fragment of the 16S rDNA mitochondrial gene. T. blanchardi, characterized by a low genetic diversity, was genetically segregated from the other sympatric T. tubifex. The ancestral state reconstruction was used to define the morphological traits that support its distinctness. On the contrary, the other T. tubifex from the Lambro community, although morphologically indistinguishable, revealed an astonishing degree of genetic variability, both between and within the three identified clades that proved to be genetically isolated. Using samples from the mixed Lambro River community and from other countries around the world we present an overview of the species complex' 16S rDNA variability. Our results show that the genetic variability did not sensibly increase widening the data set, suggesting that the Lambro River populations meet the species' worldwide genetic variability.     

Prevalence and susceptibility of infection to Myxobolus cerebralis,and genetic differences among populations of Tubifex tubifex

The prevalence of infection and susceptibility of the aquatic oligochaete Tubifex tubifex to Myxobolus cerebralis, was examined in 2 studies on the upper Colorado River, Colorado, USA, where whirling disease occurs in wild trout populations. In the first study, the prevalence of infection ranged from 0.4 to 1.5%, as determined by counting the number of T. tubifex releasing triactinomyxons of M. cerebralis directly following their collection from the field. The susceptibility of those T. tubifex not releasing triactinomyxons was assessed by the number of these oligochaetes releasing triactinomyxons 3 mo following experimental exposures to spores of M. cerebralis. The prevalence of infection following experimental exposures of these T. tubifex ranged from 4.2 to 14.1%. In a second study, all T. tubifex collected at 2 different times directly from the 2 field sites in Colorado were exposed to spores of M. cerebralis. Individual oligochaetes representing those groups of T. tubifex releasing and those groups not releasing triactinomyxons at 3 mo were screened with molecular genetic markers. T. tubifex populations found at the 2 study sites consisted of 4 genetically distinct lineages that varied with respect to their susceptibility to experimental exposure to M. cerebralis. Lineages I and III contained the most oligochaetes susceptible to M. cerebralis and were the most prominent lineages at Windy Gap Reservoir, a site of high infectivity for wild rainbow trout on the upper Colorado River. In contrast, at the Breeze Bridge site which is below Windy Gap Reservoir and where M. cerebralis infections are less severe in wild trout, oligochaetes in lineages V and VI that are resistant to M. cerebralis were more prominent. These results suggest that certain habitats, such as Windy Gap Reservoir, are conducive to large and more homogenous populations of susceptible T. tubifex lineages that may serve as point sources of infection for M. cerebralis. Although not a direct objective of this study, there was no evidence of M. cerebralis infections among any oligochaetes other than those that would be classified as T. tubifex by standard morphological characteristics.     

The severity of whirling disease among wild trout corresponds to the differences in genetic composition of Tubifex tubifex populations in central Colorado

We analyzed the geographic distribution of Tubifex tubifex from various river drainages in central Colorado by genetic screening with specific mitochondrial 16S ribosomal DNA (mt 16S rDNA) markers. Four distinct mt 16S rDNA lineages are evident. The sites varied with respect to land- and water-use practices. All sites represented habitats presumed capable of supporting oligochaetes. At the locations where whirling disease has had the greatest impact on resident rainbow trout, T. tubifex, representing lineages I and III (genotypes known to be susceptible to Mxyobolus cerebralis), were most commonly found. In contrast, at sites less affected by whirling disease, T. tubifex of lineages V and VI that are more resistant to M. cerebralis infections were more abundant. The predominance of resistant T. tubifex worms (lineages V and VI) at low-impact sites supports the conclusion that when these genotypes are in greater abundance, the potential for more severe effects of whirling disease on wild rainbow trout populations may be diminished.     

Prevalence of Myxobolus cerebralis infections among genetic lineages of Tubifex tubifex at three locations in the Madison River, Montana

Host biodiversity can impact disease risk and influence the transmission of parasitic disease. Stream sediment-dwelling worms, Tubifex tubifex (Clitellata: Oligochaeta), are the definitive host of the parasite Myxobolus cerebralis (Myxozoa: Myxosporea), which causes whirling disease in salmonid fishes. Genetic diversity of T. tubifex is correlated with host susceptibility to M. cerebralis , and mitochondrial Lineage III is generally shown to be more likely to be infected and produce the triactinomyxon (TAM) spores than other lineages. We determined the mitochondrial lineage, relative abundance, and prevalence of infection of T. tubifex collected at 3 sites in the Madison River, Montana, where previous study had shown variation in whirling disease prevalence and severity in caged trout fry. We also compared visual identification of TAMs released from cultured worms with a molecular genetic assay (diagnostic polymerase chain reaction [PCR]) for parasite detection of both infected and uninfected worms. We estimated that mitochondrial Lineage III was most abundant at the site previously shown to have high fish disease and was also most likely to be infected. The 2 techniques for detecting parasite infection did not always agree, and the likelihood of PCR (+) and spore (-) was not significantly different from PCR (-) and spore (+). Differences in the relative infection prevalence for these 2 lineages may explain the wide range of infection in natural streams.     

Tubifex tubifex from Alaska and their susceptibility to Myxobolus cerebralis

Although widespread throughout the continental United States, Myxobolus cerebralis, the myxozoan parasite that causes whirling disease in salmonids, has not been reported from the state of Alaska. As part of a risk assessment for the introduction and establishment of M. cerebralis into Alaska, the distribution of the invertebrate host Tubifex tubifex was surveyed, and its genetic composition and susceptibility to the parasite were determined. Many oligochaetes, but no T. tubifex, were collected from southeastern Alaska; however, 4 lineages of T. tubifex (I, III, IV, and VI) were identified from south-central Alaska. Lineage IV had not been previously described in North America, and its susceptibility to M. cerebralis was unknown. When lineage IV T. tubifex and 3 mixed lineage (I, III, IV, and VI) groups were exposed to M. cerebralis, only lineage III became infected under our experimental conditions. Infection occurred in this lineage even when it made up just 3% of the population. Implications of the presence of nonsusceptible lineages of T. tubifex on Alaskan salmonids would be significant in areas where these lineages dominate T. tubifex populations.     

Nuclear DNA recapitulates the cryptic mitochondrial lineages of Lumbricus rubellus and suggests the existence of cryptic species in an ecotoxological soil sentinel

Mitochondrial DNA analysis has revealed two distinct phylogenetic lineages within the ecotoxological sentinel earthworm model Lumbricus rubellus Hoffmeister, 1843. The existence of these lineages could complicate ecotoxicological studies that use the species as a sentinel for soil contamination testing, as they may respond differently to contamination; however, as mitochondrial haplotypes are not always expected to segregate in the same way as chromosomal DNA in natural populations, we further investigated this issue by using nuclear DNA markers (microsatellites) to measure genetic diversity, differentiation, and gene flow in sympatric populations of the two L. rubellus lineages at two sites in South Wales. Our results show that sympatric populations of the two lineages are more genetically differentiated than geographically distant populations of the same lineage, and Bayesian clustering analysis revealed no evidence of gene flow between the lineages at either site. Additionally, DNA sequencing of these microsatellite loci uncovered substantial differentiation between lineages at homologous flanking regions. Overall our findings indicate a high degree of nuclear genetic differentiation between the two lineages of L. rubellus, implying reproductive isolation at the two study sites and therefore the potential existence of cryptic species. The existence of two cryptic taxa has major implications for the application of L. rubellus as an ecotoxicological sentinel. It may therefore be necessary to consider the lineages as separate taxa during future ecotoxicological studies. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 780–795.     

Spatial mixing of mitochondrial lineages and greater genetic diversity in some invasive populations of the American mink (<Emphasis Type="Italic">Neovison vison</Emphasis>) compared to native populations

The genetic characteristics of introduced populations have a relevant impact on their ability to establish and spread. The American mink (Neovison vison), native to North America, is an important invasive species in the Iberian Peninsula. Here, we used mitochondrial DNA sequences data to investigate the genetic diversity and phylogeographic structure of invasive versus native populations of this species. We also evaluated whether genetic diversity in invasive populations could be explained by the genetic characteristics of the native sources from which they derived. Phylogenetic analysis revealed two major lineages in the native range, which indicated a clear separation between western and eastern populations. On the contrary, we found no evidence of genetic structure in the invasive range. This was probably the result of the diverse origins of the released specimens and the rapid expansion and encounters of the introduced populations. We detected spatial mixing of both North American lineages in several sampling localities of the north central area of the Iberian Peninsula, giving rise to high levels of genetic diversity in some areas compared to North American populations. This could potentially lead to higher fitness of these individuals and thus increase the population viability and invasiveness of this species. These results point to the need to better study the populations in which lineages mix and, if necessary, intensify control efforts in them.     

Patterns of morphological variation among populations of the widespread annual killifish Nothobranchius orthonotus are independent of genetic divergence and biogeography

Populations of annual killifish of the genus Nothobranchius occur in patchily distributed temporary pools in the East African savannah. Their fragmented distribution and low dispersal ability result in highly structured genetic clustering of their populations. In this study, we examined body shape variation in a widely distributed species, Nothobranchius orthonotus with known phylogeographic structure. We tested whether genetic divergence of major mitochondrial lineages forming two candidate species is congruent with phenotypic diversification, using linear and geometric morphometry analyses of body shape in 23 wild populations. We also conducted a common‐garden experiment with two wild‐derived populations to control for the effect of local environmental conditions on body shape. We identified different allometric trajectories for different mitochondrial lineages and candidate species in both sexes. However, in a principal components analysis of population‐level body shape, the separation among mitochondrial lineages was incomplete. Higher similarity of mitochondrial lineages belonging to different candidate species than that of same candidate species prevented distinction of the two candidate species on the basis of body shape. Analysis at the individual level demonstrated that N. orthonotus express high intrapopulation variability, with major overlap among individuals from all populations. In conclusion, we suggest that N. orthonotus be considered as a single species with an extensive geographic range, strong population genetic structure and high morphological variability.     

Host-associated genetic differentiation in Thrips tabaci (Insecta; Thysanoptera), as determined from mtDNA sequence data

We tested for host-associated genetic differentiation in 22 populations of Thrips tabaci collected from tobacco and leek, respectively. Clustering analyses and haplotype networks based on sequence variation at a fragment of the mitochondrial cytochrome oxidase I gene yielded three major evolutionary lineages; two were clearly associated with leek and the third with tobacco. These genetic findings corroborated recent experimental observations on the heterogeneity of T. tabaci populations with regard to host-plant preference and their capacity to be vectors for tomato spotted wilt virus. Estimated divergence times suggested an ancient divergence of these lineages dating back to the Miocene 28-21 million years ago. F(ST) values between these lineages ranged between 0.824 and 0.954 (P<0.001 for all comparisons), and sequence divergences ranged between 4 and 11%. Given these findings and by the standards of genetic and ecological differentiation in other published species groups, T. tabaci must be considered a complex of cryptic (sub)species.     

Genetic diversity and invasion history of the European subterranean termite <Emphasis Type="Italic">Reticulitermes urbis</Emphasis> (Blattodea,Termitoidae, Rhinotermitidae)

Biological invasions are among key factors of ecological changes, and social insects appear as highly successful invasive animals. Subterranean termites of the holarctic genus Reticulitermes are present in Europe with six native and one invasive (the nearctic R. flavipes) species. The species R. urbis shows a disjunct distribution in the Western Balkans, Eastern Italy and Southern France. Previous molecular and population genetics data suggested that the taxon originated from the Balkans, and that Italian and French populations are invasive, but it is still unknown how many introduction events occurred and from which Balkan source populations. To address these questions, a population genetics analysis was performed on a larger sampling than previous studies, using mitochondrial cytochrome oxidase II and 6 microsatellite markers on 47 colonies collected across the whole distribution area. Mitochondrial analysis confirmed the presence of two major lineages where colonies from Balkans, Italy, and France intermingle. Similarly, microsatellite loci analysis indicated the presence of two genetic clusters, though not corresponding to the two mitochondrial clades, each including colonies from the three sampled areas and with individuals showing mixed cluster membership. Overall, French and Italian populations showed indications of bottleneck (reduced genetic diversity and change of allele frequencies) and do not appear genetically differentiated from the Balkan population. Results presented here support a history of multiple introductions in Italy and France, in a scenario consistent with continuous exchanges between native and invasive areas, as expected along human trades routes.     

Variability in triactinomyxon production from Tubifex tubifex populations from the same mitochondrial DNA lineage infected with Myxobolus cerebralis, the causative agent of whirling disease in salmonids

Myxobolus cerebralis, the causative agent of whirling disease, infects both salmonid fish and an aquatic oligochaete, Tubifex tubifex. Although M. cerebralis has been detected in river drainages throughout the United States, disease severity among wild fish populations has been highly variable. Tubifex tubifex populations have been genetically characterized using sequences from the 16S mitochondrial DNA (mtDNA) gene, the 18S ribosomal RNA gene, the internal transcribed spacer region 1 (ITS1), and randomly amplified polymorphic DNA (RAPD). Our earlier work indicated that large differences in compatibility between the parasite and populations of T. tubifex may play a substantial role in the distribution of whirling disease and resulting mortality in different watersheds. In the present study, we examined 4 laboratory populations of T. tubifex belonging to 16S mtDNA lineage III and 1 population belonging to 16S mtDNA lineage I for triactinomyxon (TAM) production after infection with M. cerebralis myxospores. All 4 16S mtDNA lineage III populations produced TAMs, but statistically significant differences in TAM production were observed. Most individuals in the 16S mtDNA lineage III-infected populations produced TAMs. The 16S mtDNA lineage I population produced few TAMs. Further genetic characterization of the 16S mtDNA lineage III populations with RAPD markers indicated that populations producing similar levels of TAMs had more genetic similarity.     

Global major histocompatibility Class II β (mh-IIβ)-polymorphism in Arctic charr Salvelinus alpinus

This study explored the use of the gene encoding the β subunit of the major histocompatibility (MH) receptor as a population marker in Arctic charr Salvelinus alpinus. The use of this polymorphic marker allowed differentiation of the S. alpinus lineages previously defined using mitochondrial DNA (mtDNA) but also allowed differentiation between the populations studied within those lineages. The majority of the variation observed here occurred prior to the last glaciation event. Nevertheless, all S. alpinus populations were differentiated using both MH Class II β (mh-IIβ) sequences and allelic frequencies. The fact that all the populations studied presented high rates of non-synonymous: synonymous substitutions and high levels of interpopulation variation, suggested mh-IIβ as an ideal marker to assess differentiation among S. alpinus populations in ways that may represent divergence both by genetic drift and natural adaptation to the local environment.     

Phylogeography, genetic structure and diversity in the endangered bearded vulture (Gypaetus barbatus, L) as revealed by mitochondrial DNA

Bearded vulture populations in the Western Palearctic have experienced a severe decline during the last two centuries that has led to the near extinction of the species in Europe. In this study we analyse the sequence variation at the mitochondrial control region throughout the species range to infer its recent evolutionary history and to evaluate the current genetic status of the species. This study became possible through the extensive use of museum specimens to study populations now extinct. Phylogenetic analysis revealed the existence of two divergent mitochondrial lineages, lineage A occurring mainly in Western European populations and lineage B in African, Eastern European and Central Asian populations. The relative frequencies of haplotypes belonging to each lineage in the different populations show a steep East-West clinal distribution with maximal mixture of the two lineages in the Alps and Greece populations. A genealogical signature for population growth was found for lineage B, but not for lineage A; futhermore the Clade B haplotypes in western populations and clade A haplo-types in eastern populations are recently derived, as revealed by their peripheral location in median-joining haplotype networks. This phylogeographical pattern suggests allopatric differentiation of the two lineages in separate Mediterranean and African or Asian glacial refugia, followed by range expansion from the latter leading to two secondary contact suture zones in Central Europe and North Africa. High levels of among-population differentiation were observed, although these were not correlated with geographical distance. Due to the marked genetic structure, extinction of Central European populations in the last century re-sulted in the loss of a major portion of the genetic diversity of the species. We also found direct evidence for the effect of drift altering the genetic composition of the remnant Pyrenean population after the demographic bottleneck of the last century. Our results argue for the management of the species as a single population, given the apparent ecological exchangeability of extant stocks, and support the ongoing reintroduction of mixed ancestry birds in the Alps and planned reintroductions in Southern Spain.     

Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species

While Quaternary climatic changes are considered by some to have been a major factor promoting speciation within the neotropics, others suggest that much of the neotropical species diversity originated before the Pleistocene. Using mitochondrial and nuclear sequence data, we evaluate the relative importance of Pleistocene and pre-Pleistocene events within the evolutionary history of the Mexican iguana Ctenosaura pectinata , and related species. Results support the existence of cryptic lineages with strong mitochondrial divergence (> 4%) among them. Some of these lineages form zones of secondary contact, with one of them hybridizing with C. hemilopha . Evolutionary network analyses reveal the oldest populations of C. pectinata to be those of the northern and southern Mexican coastal regions. Inland and mid-latitudinal coastal populations are younger in age as a consequence of a history of local extinction within these regions followed by re-colonization. Estimated divergence times suggest that C. pectinata originated during the Pliocene, whereas geographically distinct mitochondrial DNA lineages first started to diverge during the Pliocene, with subsequent divergence continuing through the Pleistocene. Our results highlight the influence of both Pliocene and Pleistocene events in shaping the geographical distribution of genetic variation within neotropical lowland organisms. Areas of high genetic diversity in southern Mexico were detected, this finding plus the high levels of genetic diversity within C. pectinata , have implications for the conservation of this threatened species.     

Admixture of genetic lineages of different glacial origin: a case study of <Emphasis Type="Italic">Abies alba</Emphasis> Mill. in the Carpathians

Genetic exchange between divergent lineages of silver fir (Abies alba Mill.) was studied in the Ukrainian Carpathians where two expanding populations originating from different glacial refugia meet. The study included 33 silver fir populations from Ukraine, Romania and Slovakia which were investigated using the maternally inherited mitochondrial nad5-4 marker and biparentally inherited nuclear microsatellites. The boundary between mitochondrial lineages is very sharp; only six populations containing a mixture of different haplotypes were found. Bayesian analysis of population structure based on seven nSSR loci revealed the existence of two clusters which coincided fairly well with mitochondrial lineages. Both haplotype frequencies and proportions of clusters identified by the Bayesian analysis exhibited a clinal transition over the contact zone, with cline widths of 17.6 km for mitochondrial haplotype frequencies (reflecting gene flow by seeds) and 119.6 km for Bayesian clusters based on nSSR (reflecting gene flow by pollen). Allelic richness and gene diversity differ significantly between mitochondrial lineages, the Balkan group being more variable, but an increase in gene diversity towards the boundary between lineages was observed only within the Balkan lineage. The observed patterns are suggested to reflect the postglacial colonization and historical gene flow. They demonstrated that in silver fir as a wind-dispersed species, the colonization front is quite continuous, and the survival of migrant seeds in already established populations is low. The relative contribution of pollen-mediated gene flow to genetic exchange between divergent lineages associated with glacial refugia is much higher than in the case of seeds, but pollen dispersal distance is lower than suggested by earlier studies.     

Opening a can of worms: unprecedented sympatric cryptic diversity within British lumbricid earthworms

Earthworms play a major role in many aspects of soil fertility, food web ecology and ecosystem functioning, and hence are frequently the subjects of, for example, ecological and toxicological research. Our aim was to examine the genetic structure of common earthworm species, to identify cryptic lineages or species that may be distinct ecotypes or biotypes (and hence confound current research based upon morphotypes) and to try to explain the massive cryptic diversity that eventually emerged. We demonstrated that species such as Allolobophora chlorotica, Aporrectodea longa, Aporrectodea rosea and Lumbricus rubellus all comprise highly divergent lineages with species-level divergence at the mitochondrial cytochrome oxidase I (COI) gene. In Allo. chlorotica alone, we found 55 haplotypes for COI, with 35 of these being found in pink and 20 in green morph worms. There were no cases of the two colour morphs sharing COI haplotypes. Phylogenetic analyses of mitochondrial COI and 16S genes showed the presence of five highly divergent lineages, suggesting the presence of multiple cryptic species within Allo. chlorotica. There was no clear geographical pattern to lineage distribution and many populations were polymorphic for both mitochondrial DNA lineage and colour morph. Amplified fragment length polymorphism results, based on two primer combinations, were broadly congruent with mitochondrial DNA results with one significant exception. Despite showing over 14% divergence at COI, amplified fragment length polymorphism markers showed that the two green morph lineages may be interbreeding and therefore represent a single taxon. The cryptic diversity revealed by these results has profound consequences for all areas of earthworm research.     

Genetic structure and distinctness of Apis mellifera L. populations from the Canary Islands

The genetic structure of Apis mellifera populations from the Canary Islands has been assessed by mitochondrial (restriction fragment length polymorphisms of the intergenic transfer RNAleu-COII region) and nuclear (microsatellites) studies. These populations show a low level of genetic variation in terms of average number of alleles and degree of heterozygosity. Significant differences in the distribution of alleles were found in both data sets, confirming the genetic differentiation among some of the islands but not within them. Two mitochondrial haplotypes characteristic of the Canary Islands are found at high frequencies, although populations are introgressed by imported honeybees of eastern European C lineage. This introgression is rather high on Tenerife and El Hierro and low on Gran Canaria and La Gomera, whereas on La Palma it has not been recorded. The finding of microsatellite alleles characteristic of the eastern European lineage corroborates the genetic introgression. Phylogenetic analyses indicate that the Canarian honeybees are differentiated from other lineages and provide genetic evidence of their African origin.     

Sister species within the Triops cancriformis lineage (Crustacea, Notostraca)

We investigated the phylogenetic relationships among the three presently recognized subspecies of the tadpole shrimp, Triops cancriformis , using mitochondrial 16S and 12S rDNA sequences. Our results indicate that the taxon is divided into two distinct lineages. One lineage is formed of T. c. cancriformis populations and samples from northern Spain that had been classified as T. c. simplex in the most recent literature. The second lineage comprises all populations of T. c. mauritanicus and northern African populations of T. c. simplex . These two main lineages separated 2.3 to 8.9 million years ago, based on the range of inferred molecular clocks recognized for crustacean mtDNA sequence divergence. Percentages of divergence are in the range reported for recognized species in other notostracan lineages and we therefore propose to recognize them as two species, Triops cancriformis and Triops mauritanicus . The latter would comprise two subspecies in northern Africa, one consisting of the Moroccan populations of the former T. c. mauritanicus , the other comprising the African populations of the former T. c. simplex . It also includes three as-yet unnamed lineages. A comparison of morphological characters with the molecular data revealed that the former T. c. simplex cannot be reliably separated from T. c. cancriformis , using morphological characters that have hitherto been used to distinguish among subspecies of T. cancriformis . Our investigation is the first to demonstrate the presence of T. c. cancriformis in Africa (Tunisia). The genetic haplotypes of these populations are identical with haplotypes also occurring in Central and Western Europe, as well as in Sicily. Therefore, we hypothesize that the African populations of T. c. cancriformis represent a result of repeated long-distance dispersal across the Mediterranean Sea.