Parasitic infection by larval helminths in Antarctic fishes: pathological changes and impact on the host body condition index.

We examined pathological changes and relationship between body condition index (BCI) and parasitic infection in 5 species of fish, including 42 icefish Chionodraco hamatus (Channichtyidae), 2 dragonfish Cygnodraco mawsoni (Bathydraconidae), 30 emerald rock cod Trematomus bernacchii, 46 striped rock cod T. hansoni and 9 dusty rock cod T. newnesi (Nototheniidae) from the Ross Sea, Antarctica. All parasites were identified by a combination of morphology and mtDNA cytochrome-oxidase-2 sequence (mtDNA cox2) analysis, except Contracaecum osculatum s.l., for which only the latter was used. Five larval taxa were associated with pathological changes including 2 sibling species (D and E) of the C. osculatum species complex and 3 cestodes including plerocercoids of a diphyllobothridean, and 2 tetraphyllidean forms including cercoids with monolocular and bilocular bothridia. The most heavily infected hosts were C. hamatus and C. mawsoni, with C. hamatus most often infected by C. osculatum sp. D and sp. E and diphyllobothrideans, while C. mawsoni was most often infected with tetraphyllidean forms. Histologically, all fish showed varying severity of chronic inflammation associated with larval forms of helminths. Diphyllobothrideans and C. osculatum spp. were located in gastric muscularis or liver and were associated with necrosis and mild to marked fibrosis. Moderate multifocal rectal mucosal chronic inflammation was associated with attached tetraphyllidean scolices. C. hamatus showed a strong negative correlation between BCI and parasite burden.

Different patterns of introgression in a three species hybrid zone among European cave salamanders.

Hybrid zones occur where genetically distinct populations meet, mate and produce offspring with mixed ancestry. In Plethodontid salamanders, introgressive hybridization is a common phenomenon, where hybrids backcross with parental populations leading to the spread of new alleles into the parental genomes. Whereas many hybrid zones have been reported in American Plethodontid salamanders, only a single hybrid zone has been documented in European plethodontids so far, which is located at the Apuan Alps in the Italian Peninsula. Here, we describe a previously unreported hybrid zone in the Northern Apennines involving all the three Plethodontid salamander species inhabiting the Italian Peninsula. We found 21 new Speleomantes sites of occurrence, from a hitherto unexplored area located at the boundaries between three Speleomantes species ranges. Using mitochondrial (Cytb and ND2 genes) and nuclear markers (two diagnostic SNPs at the NCX1 gene), we revealed a three-way contact zone where all the three mainland species hybridize: S. strinatii, S. ambrosii and S. italicus. We observed a strong mitonuclear discordance, with mitochondrial markers showing a conspicuous geographic pattern, while diagnostic nuclear SNPs coexisted in both the same populations and individuals, providing evidence of hybridization in many possible combinations. The introgression is asymmetric, with S. italicus mitogenome usually associated with S. a. ambrosii and, to a lesser extent, to S. strinatii nuclear alleles. This finding confirms that Plethodontid are a group of choice to investigate hybridization mechanisms and suggests that behavioural, genetic and ecological components may concur in determining the direction and extent of introgression.

The use of cyprinodont fish, Aphanius fasciatus, as a sentinel organism to detect complex genotoxic mixtures in the coastal lagoon ecosystem.

In the present work we aimed to standardise the alkaline comet assay with erythrocytes of the cyprinodont, Mediterranean Killifish, Aphanius fasciatus. The aims of the study were to explore the suitability of this fish to assess biomarkers of genotoxic effects and as a sentinel organism to detect complex genotoxic mixtures in coastal lagoon ecosystems. Following proper optimisation, the application and effectiveness of the comet assay in erythrocytes of A. fasciatus were tested by measuring the tail DNA (%) induced by (a) in vivo exposure of individual fish to X-rays (dose, 3Gy) and (b) following in vitro challenge of erythrocytes with restriction endonucleases Fok-I and Eco-RI, which selectively induce double-strand breaks with cohesive and blunt termini, respectively. Furthermore, in order to evaluate whether circulating fish blood contained actively proliferating cells that could influence the extent of DNA damage in control (untreated) fish, we measured the number of "comets" positive for 5-bromo-2'-deoxyuridine (BrdU) by the use of anti-BrdU antibody and immuno-histochemical methods. Both treatments (i.e. with X-rays and restriction endonucleases) induced statistically significant increases in tail DNA (%) values compared with the relevant untreated controls, indicating the effectiveness of the comet assay in the erythrocytes of A. fasciatus to detect different types of DNA lesions. Results from anti-BrdU antibody labelling of erythrocytes indicated a very low percentage (5%) of "comets" positive for BrdU. Following optimisation and validation of the assay under laboratory conditions, fish were collected in the Orbetello lagoon (Tuscany, Italy), considered to be a significantly polluted site. The results showed statistically significant increases for tail DNA (%) compared with corresponding values observed in erythrocytes of fish caught in the unpolluted reference site "Saline di Tarquinia". The effects of physico-chemical parameters of the water (i.e., salinity, pH and oxygen content) did not significantly influence the induction of DNA damage. These results indicate that the comet assay provides a reliable parameter and that A. fasciatus is a promising "sentinel organism" to detect the genotoxic impact of complex mixtures in coastal lagoon ecosystems.

Genetic diversity of the killifish Aphanius fasciatus paralleling the environmental changes of Tarquinia salterns habitat.

The habitat in the Natural Reserve of the Tarquinia salterns, located on the Tyrrhenian coast of central Italy, has undergone dramatic alterations over the last 10 years. After salt production was terminated in 1997 the site was abandoned until 2002, with consequent degradation of habitat quality and stiffening of the environmental conditions. From 2003 to 2006 ecological rehabilitation of the site was carried out, restoring water circulation to its previous equilibrium. The genetic variation in the killifish Aphanius fasciatus inhabiting the salterns was monitored using allozymes from 1998. The results showed that the genetic variability of the killifish strongly reduced through time: a high number of rare alleles were lost and both heterozygosity and allele richness were significantly decreased. The most recent samples, taken after the ecological restoration, showed that to date the genetic erosion of A. fasciatus gene pool has slowed down, since no significant differences have been detected for any genetic variability parameter. Concerning the mechanisms leading to the impoverishment of the genetic variability, the strong loss of rare alleles suggests a role of genetic drift, which accords with the fluctuation of the effective population size recorded over the period of study and with the low gene flow typical of this species. The low levels of gene flow reported for this species imply that once lost, the genetic variability can rarely be restored through immigration from highly variable populations.

Food preference determines the best suitable digestion protocol for analysing microplastic ingestion by fish.

Microplastic presence in the marine environment has generated considerable concern. Many procedures for microplastics detection in fish gastrointestinal tract have been recently developed. In this study, we compared efficiencies of two common procedures applied for the digestion of organic matter (10% KOH; 15% H2O2) with a new proposal (mixture of 5% HNO3 and 15% H2O2). We considered ecological diversity among species and differences in their diet compositions as factors that could affect the efficiency and feasibility of analytical approaches. Our aim was to understand whether either one of the three protocols might be suitable for all species or it might be more advisable to select a method according to the gut content determined by different food preferences. The results showed that the trophic level and feeding habits should be considered for protocol selection. Finally, we applied the best protocols on samples from the Tyrrhenian sea.

A network approach to identify bioregions in the distribution of Mediterranean amphipods associated with Posidonia oceanica meadows.

Although amphipods are key components of the macro-fauna associated with Posidonia oceanica meadows, to date no studies focused on the structure and diversity of their assemblages across the whole Mediterranean Sea. Here, we applied a network approach based on modularity on a dataset mined from literature to identify biogeographic modules and to assess the biogeographic roles of associated localities. We also correlated the patterns evidenced with the biogeographic distribution of amphipod groups by means of a multivariate analysis. Modularity analysis highlighted four biogeographic modules bounded by the main Mediterranean biogeographic divides and evidenced a decrease in species diversity along a NW-SE gradient. Assemblages associated with Central-Western Mediterranean and, to a lesser extent, Tunisian modules showed the highest species richness and were identified as hubs, characterized by species with regional distributions that behave as source in a biogeographic context. The paleogeographic history of the host seagrass and the ecology of associated amphipods, both suggest the joint effect of species persistence and post-Last Glacial Maximum expansion in explaining the pattern of amphipod distribution in the Mediterranean Sea.

Molecular and Structural Characterization of MHC Class II ß Genes Reveals High Diversity in the Cold-Adapted Icefish Chionodraco hamatus.

This study reports the presence of two distinct MHC class II ß genes in the Antarctic icefish Chionodraco hamatus, belonging to the classical (ChhaDAB) and nonclassical (ChhaDBB) evolutionary lineages. By the application of targeted sequencing approach, a remarkable molecular diversity in the exon 2 sequence of the highly expressed gene ChhaDAB has been observed, resulting in an estimate of 92 different variants translated in 87 different peptides from 54 analysed icefish individuals. A highly conservative estimate, based on a 95% sequence identity threshold clustering, translate this variability in 41 different peptide clusters belonging to four different clades and showing the signature of different kinds of selection. In stark contrast, the poorly expressed ChhaDBB gene displayed a very low level of molecular diversity within exon 2, in agreement with expectations for a nonclassical MHC class II ß gene.

Population genetic structure and diversity of the Apennine endemic stream frog, Rana italica--insights on the Pleistocene evolutionary history of the Italian peninsular biota.

For most species in the Western Palaearctic region, southern Mediterranean peninsulas have been identified as major Quaternary refugia and hotspots of intraspecific diversity, and thus, as areas of particular relevance for the conservation of the evolutionary potential. We analysed the patterns of geographical variation among 26 populations of the Italian stream frog, using both nuclear (allozymes) and mitochondrial (partial cytochrome b sequences) markers. Phylogenetic, phylogeographical and population genetic analyses suggested that the species survived the last glacial-interglacial cycles in two distinct refugia, one restricted to the tip of the Calabrian peninsula, at the extreme south of the species' range, the other spanning from central Calabria to central Apennines and showing evidences for further population subdivision therein. Historical demographic tests suggested a significant population expansion from the latter, which most likely began around the last pleniglacial. This expansion would have led to the rapid colonization of the northern Apennines to the north, and to a secondary contact and population admixture with the population from the southern refugium in southern central Calabria. A comparison of the evolutionary history inferred for the Italian stream frog with the data emerging for other codistributed species suggests: (i) the generality of a multiple-refugia scenario for the Italian peninsula, (ii) the possible occurrence of at least one suture zone in southern Italy, and (iii) that for most species, this Pleistocene refugium is not only a hotspot, but also a melting pot of intraspecific genetic diversity. Finally, the conservation implications of these results are also briefly highlighted.

The role of habitat choice in micro-evolutionary dynamics: An experimental study on the Mediterranean killifish Aphanius fasciatus (Cyprinodontidae).

Habitat choice is defined as a nonrandom distribution of genotypes in different microhabitats. Therefore, it could exert a great impact on the genetic variance of natural populations by promoting genetic divergence, local adaptation, and may even lead to sympatric speciation. Despite this potential role in micro- and macro-evolutionary processes, there is little empirical evidence that the various genotypes within a population may differ in habitat choice-related behaviors. Here, we tested whether habitat choice may have contributed to genetic divergence within a local population of the Mediterranean killifish Aphanius fasciatus, which emerged between groups inhabiting microhabitats with different oxygen concentrations during previous field studies. In a first experiment, we studied the distribution of individuals in conditions of hypoxia and normoxia to test whether they had a different ability to shy away from a hypoxic environment; in a second experiment, we analyzed the individual behavior of fish separately in the two conditions, to verify whether they showed peculiar behavioral responses linked to a possible differential distribution. We then analyzed the six allozyme loci, whose allelic and genotypic frequencies were significantly divergent in the previous studies. In the first test, we found that the distribution of the two homozygote genotypes of the glucose-6-phosphate isomerase-1 locus (GPI-1) was significantly different between the hypoxic and the normoxic conditions. During the second test, all individuals were more active in hypoxic conditions, but the two GPI-1 homozygotes showed a significant difference in time spent performing surface breathing, which was consistent with their distribution observed in the first experiment. These results provide evidence that individual behavioral traits, related to genetic features, may lead to a nonrandom distribution of genotypes in heterogeneous although contiguous microhabitats and, consequently, that habitat choice can play a significant role in driving the micro-evolutionary dynamics of this species.

Persistence, isolation and diversification of a naturally fragmented species in local refugia: the case of Hydromantes strinatii.

The study of the European plethodontid salamander Hydromantes strinatii using allozyme and mitochondrial markers showed a strong geographical genetic structure. This was likely the outcome of different evolutionary mechanisms leaving their signature despite the effects of the genetic drift due to the low population size typical of this species. Two highly divergent clades were identified in the eastern and central-western part of the range, with further geographic sub-structure. Nuclear and mitochondrial markers substantially recovered the same population groups but were conflicting in reconstructing their relationships. This apparent incongruence highlighted the action of different mechanisms such as secondary contacts and incomplete lineage sorting in originating the observed genetic variation. The troglophilic habit of this species provided the opportunity to show the importance of caves as local refugia in maintaining the genetic diversity through the persistence of local populations. Accordingly, high nucleotide and haplotype diversity, strong geographic genetic structuring and lack of expansion were evidenced. This signature was found in the populations from the Ligurian and Maritime Alps, in agreement with the complex orography and paleoclimatic history of this Mediterranean hotspot.

Metazoan parasite infection in the swordfish, Xiphias gladius, from the Mediterranean Sea and comparison with Atlantic populations: implications for its stock characterization.

Thirteen parasite taxa were identified in the Mediterranean swordfish by morphological and genetic/molecular methods. The comparison of the identified parasite taxa and parasitic infection values observed in the Mediterranean swordfish showed statistically significant differences with respect to those reported for its Atlantic populations. A stepwise Linear Discriminant Analysis of the individual fish examined showed a separation among three groups: one including fish from the Mediterranean Sea (CTS, STS, and IOS); one consisting of fish from the Central South (CS), Eastern Tropical (ET), and Equatorial (TEQ) Atlantic; and a third comprising the fish sampled from the North-West Atlantic (NW); the CN Atlantic sample was more similar to the first group rather than to the other Atlantic ones. The nematodes Hysterothylacium petteri and Anisakis pegreffii were the species that contributed most to the characterization of the Mediterranean swordfish samples with respect to these Atlantic ones. Anisakis brevispiculata, A. physeteris, A. paggiae, Anisakis sp. 2, Hysterothylacium incurvum, Hepatoxylon trichiuri, Sphyriocephalus viridis, and their high infection levels were associated with the swordfish from the Central and the Southern Atlantic areas. Finally, H. corrugatum, A. simplex (s.s.), Rhadinorhynchus pristis, and Bolbosoma vasculosum were related to the fish from the North-West (NW) Atlantic area. These results indicate that some parasites, particularly Anisakis spp. larvae identified by genetic markers, could be used as "biological tags" and support the existence of a Mediterranean swordfish stock.

Parasite communities of icefish (Chionodraco hamatus) in the Ross Sea (Antarctica): influence of the host sex on the helminth infracommunity structure.

Parasite communities of Chionodraco hamatus were investigated from Terra Nova Bay (Ross Sea, Antarctica) during host spawning time. Special attention was given to helminth infracommunities and effect of host sex on its structure. A total of 21 taxa including 5 ecto-parasites and 16 endo-parasites were identified. The number of ecto and endo-parasite species per individual host ranged from 1 to 3 and 3 to 10, respectively, while the mean numbers of parasite specimens per individual host were 4.7 and 1309.7, respectively. The rich abundance of infection suggests a rich concentration of helminth intermediate/paratenic hosts in the coastal waters of Terra Nova Bay. Chionodraco hamatus serves as a definitive host for 10 helminth taxa, while it acts as an intermediate/paratenic host for 6 helminth taxa. Larvae of 6 helminth taxa for which C. hamatus serves as intermediate/paratenic host represented 98.7% of all specimens found. Of these, the tetraphyllidean and diphyllobothridean cestodes and the nematode Contracaecum osculatum s.l. were the most prevalent and abundant. 'Larval' infracommunities had significantly higher species richness, total abundance and diversity than 'adult' infracommunities, suggesting the important role of C. hamatus in supporting the life cycles of those parasites in the study area as a paratenic/intermediate host. Significant differences in the pattern of helminth infracommunities of larval forms between male and female fish were found. These differences could be caused by physiological, and most probably by behavioral differences between sexes suggesting that sex is an important factor influencing parasite burden in C. hamatus during reproductive season.

Molecular identification, morphological characterization and new insights into the ecology of larval Pseudoterranova cattani in fishes from the Argentine coast with its differentiation from the Antarctic species, P. decipiens sp. E (Nematoda: Anisakidae).

Larvae of the genus Pseudoterranova constitute a risk for human health when ingested through raw or undercooked fish. They can provoke pseudoterranovosis in humans, a fish-borne zoonotic disease whose pathogenicity varies with the species involved, making their correct specific identification a necessary step in the knowledge of this zoonosis. Larvae of Pseudoterranova decipiens s.l. have been reported in several fish species from off the Argentine coasts; however, there are no studies dealing with their specific identification in this region. Here, a genetic identification and morphological characterization of larval Pseudoterranova spp. from three fish species sampled from Argentine waters and from Notothenia coriiceps from Antarctic waters was carried out. Larvae were sequenced for their genetic/molecular identification, including the mitochondrial cytochrome c oxidase subunit II (mtDNA cox2), the first (ITS-1) and the second (ITS-2) internal transcribed spacers of the nuclear ribosomal DNA, and compared with all species of the P. decipiens (sensu lato) species complex (sequences available in GenBank). Further, adults of Pseudoterranova spp. from the definitive host, the southern sea lion, Otaria flavescens, from Argentine and Chilean coasts were sequenced at the same genes. The sequences obtained at the ITS-1 and ITS-2 genes from all the larvae examined from fish of Argentine waters, as well as the adult worms, matched 100% the sequences for the species P. cattani. The sequences obtained at mtDNA cox2 gene for Antarctic larvae matched 99% those available in GenBank for the sibling P. decipiens sp. E. Both MP and BI phylogenetic trees strongly supported P. cattani and P. decipiens sp. E as two distinct phylogenetic lineages and depicted the species P. decipiens sp. E as sister taxon to the remaining taxa of the P. decipiens complex. Larval morphometry was similar between specimens of P. cattani from Argentina, but significantly different from those of P. decipiens sp. E, indicating that larval forms can be distinguished based on their morphology. Pseudoterranova cattani is common and abundant in a variety of fish species from Chile, whereas few host species harbour these larvae in Argentina where they show low levels of parasitism. This pattern could arise from a combination of factors, including environmental conditions, density and dietary preferences of definitive hosts and life-cycle pathways of the parasite. Finally, this study revealed that the life-cycle of P. cattani involves mainly demersal and benthic organisms, with a marked preference by large-sized benthophagous fish.

Terrestrial gastropods (Helix spp) as sentinels of primary DNA damage for biomonitoring purposes: a validation study.

We validated the alkaline comet assay in two species of land snail (Helix aspersa and Helix vermiculata) to test their suitability as sentinels for primary DNA damage in polluted environments. The study was conducted under the framework of a biomonitoring program for a power station in Central Italy that had recently been converted from oil to coal-fired plant. After optimizing test conditions, the comet assay was used to measure the % Tail DNA induced by in vitro exposure of hemocytes to different concentrations of a reactive oxygen species (H2 O2 ). The treatment induced significant increases in this parameter with a concentration effect, indicating the effectiveness of the assay in snail hemocytes. After evaluating possible differences between the two species, we sampled them in three field sites at different distances from the power station, and in two reference sites assumed to have low or no levels of pollution. No species differences emerged. Percent Tail DNA values in snails from the sites near the power station were higher than those from control sites. An inverse correlation emerged between % Tail DNA and distance from the power station, suggesting that the primary DNA damage decreased as distance increased away from the pollution source. Detection of a gradient of heavy metal concentration in snail tissues suggests that these pollutants are a potential cause of the observed pattern. The comet assay appears to be a suitable assay and Helix spp. populations suitable sentinels to detect the genotoxic impact of pollutants.

Infaunal macrobenthic community dynamics in a manipulated hyperhaline ecosystem: a long-term study.

BACKGROUND: Understanding the responses of ecological communities to human-induced perturbations is crucial for establishing conservation goals. Ecological communities are dynamic entities undergoing fluctuations due to their intrinsic characteristics as well as anthropogenic pressures varying over time. In this respect, long-term studies, based on large spatial and temporal datasets, may provide useful information in understanding patterns and processes influencing the communities' structure. Theoretical evidence suggests that a role of biodiversity is acting as a compensatory buffer against environmental variability by decreasing the temporal variance in ecosystem functioning and by raising the level of community response to perturbations through the selection of better performing species. Therefore, the spatial and temporal changes in the specialization of the community components may be used as an effective tool to monitor the effects of natural and anthropogenic alterations of the environment in dynamic systems. We examined the temporal dynamics of macroinvertebrate community structure in the hyperhaline habitat of Tarquinia Saltworks (central Italy). We aimed at: (i) investigating the relationships between the level of community specialization and the alterations of the environment across fourteen years; (ii) comparing the ability of aggregate community parameters such as the average abundance vs. species specialization in describing patterns of community composition. RESULTS: We arranged the data in three sub-sets according to three periods, each characterized by different environmental conditions. The mean abundance of sampled macroinvertebrates showed a significant change (p < 0.01) only in the community inhabiting the saltwork basin closely connected to the sea, characterized by the highest environmental variation (i.e. the coefficient of variation, CV, of the aggregate environmental variability over the study period, CVrange = 0.010 - 0.2). Here we found marine species like Modiolus adriaticus (Lamarck, 1819), Neanthes irrorata (Malmgren, 1867), and Amphiglena mediterranea (Leydig, 1851), which inhabited the saltworks during the halt period but disappeared during the subsequent eutrophication phase. Conversely, species specialization showed a significant decrease for each sampled community in the presence of habitat degradation and a recovery after ecological restoration. The widest fluctuations of specialization were recorded for the community inhabiting the saltwork basin with the highest long-term environmental variability. CONCLUSIONS: Recent advances have shown how the increased temporal and spatial variability of species' abundance within the communities may be a signal of habitat disturbance, even in the absence of an apparent decline. Such approach could also be used as a sensitive monitoring tool, able to detect whether or not communities are subjected to increasing biotic homogenization. Also, the increased functional similarity triggered by habitat degradation may impact on species at higher trophic levels, such as the waterbirds wintering in the area or using it as a stopover during migration.

Phylogeography and historical demography of the Italian treefrog, Hyla intermedia, reveals multiple refugia, population expansions and secondary contacts within peninsular Italy.

We investigated the geographical patterns of genetic diversity in the Italian treefrog through sequence analysis of a mitochondrial cytochrome b gene fragment. Three main mitochondrial lineages were identified, distributed in northern, central and southern Italy, respectively. Their divergence appears indicative of a split time largely predating Late Pleistocene climatic oscillations, and syntopy between them was only observed in the geographically intermediate populations. The historical demographic reconstructions suggest that in both northern and central Italy, an expansion occurred during the last major glacial phase, when a vast widening of the lowland habitats followed the glaciation-induced fall of the sea level. Instead, in southern Italy an expansion event likely followed the end of the last glaciation, although the inference of expansion appears less reliable for the southern clade than for the others. Within this geographical area, a sharp phylogeographic discontinuity separated peninsular from Sicilian populations, and the overall pattern of diversity suggests that the latter derived from a recent colonization of the island, probably through a Late Pleistocene land bridge. Phylogenetic, phylogeographic and historical demographic analyses thus concur in delineating a scenario of multiple refugia, with four groups of populations which survived the last glacial-interglacial cycles in at least three distinct refugia arranged along peninsular Italy, and have recently come into contact following range expansions. Therefore, these results support the hypothesis that a plethora of microevolutionary processes, rather than the prolonged stability of populations, were mainly responsible for shaping the patterns of diversity within this major biodiversity hotspot.

Genetic diversity and phylogeography of the Apennine yellow-bellied toad Bombina pachypus, with implications for conservation.

Genetic variation was investigated in 17 populations of the Italian endemic Apennine yellow-bellied toad using both mitochondrial (598 bp of the cytochrome b gene) and nuclear (21 allozyme loci) markers. Populations from central Calabria (southern Italy) showed the highest levels of intrapopulation genetic variation, whereas samples located north of this region were nearly lacking in variation. This appears to be a typical pattern of 'southern richness and northern purity', usually attributed to the prolonged population stability within southern refugia coupled with the loss of variation during postglacial northward expansion. However, the overall pattern of genetic variation observed has a strong geographical component, suggesting two Calabrian plains, Catanzaro and Crati-Sibari, as historical barriers to dispersal separating three population groups. These findings cannot be explained by the prolonged stability of southern populations alone, and suggest that the southern richness has been at least in part shaped by allopatric differentiation within the refugial range, followed by intermixing of previously differentiated lineages. From a conservation standpoint, Calabria is the major genetic diversity reservoir for this species, thus deserving particular conservation efforts. Furthermore, although the low intrapopulation genetic variation outside Calabria appears to be of clear historical origin, evidence of a current reduction of gene flow suggests that human disturbance has also played a part, particularly in the anthropogenic impacted Volturno river drainage basin.