Skip to main content

Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann’s tortoise Testudo hermanni



Hotspots of intraspecific genetic diversity represent invaluable resources for species to cope with environmental changes, and their identification is increasingly recognized as a major goal of conservation ecology research. However, even for iconic and endangered species, conservation strategies are often planned without thorough information on the geographic patterns of genetic variation. Here, we investigated the spatial patterns of genetic variation of the endangered Hermann’s tortoise Testudo hermanni in the Italian Peninsula by genotyping 174 individuals at 7 microsatellite loci, with the aim to contribute to planning effective conservation strategies.


Ordination-based and Bayesian clustering analyses consistently identified three main genetic clusters, one spread in the central and northern part of the peninsula, and two restricted to southern Italy and Sicily, respectively. The highest levels of genetic diversity were found in populations of the southern cluster and, in particular, at the northern edges of its distribution (He > 0.6, Ar > 2.8 ), that correspond to areas of putative secondary contact and admixture between distinct lineages. Our results clearly identify a hotspot of genetic diversity for the Hermann’s tortoise in southern Italy.


We inferred the evolutionary history and the spatial patterns of genetic variation of the Hermann’s tortoise in the Italian Peninsula. We identified three main genetic clusters along the peninsula and a hotspot of intraspecific diversity in southern Italy. Our results underline the urgent need for conservation actions to warrant the long-term persistence of viable tortoise populations in this area. Furthrmore, these data add further evidence to the role of southern Italy as a biodiversity hotspot for temperate fauna, claiming for higher consideration of this area in large scale conservation programs.

Peer Review reports


Hotspots of genetic diversity are geographic areas harbouring exceptionally high levels of intraspecific genetic variation and represent appealing resources to investigate the mechanisms shaping the geographic structure of biodiversity [1,2,3]. Still, as genetic variation provides populations with the potential to adapt to environmental changes [4, 5], hotspots of genetic diversity also represent invaluable resources for species to cope with the global change [3, 6]. Thus, the identification of hotspots is a primary goal of evolutionary research, as well as a crucial step for designing effective conservation strategies for the long-term persistence of species and populations [7,8,9,10].

In the western Palaearctic region, Mediterranean peninsulas have been identified as major hotspots of intraspecific genetic variation for temperate species [11], and references therein]. The concentration of hotspots in these areas has been linked to the outcomes of Pleistocene glacial cycles, which induced populations to repeated cycles of fragmentation and allopatric differentiation, followed by secondary contact and admixture [2, 11,12,13,14]. However, during the last decade a new generation of phylogeographic studies has outlined more complex patterns of intraspecific genetic structure within Mediterranean peninsulas. In particular, two prominent patterns have been recurrently observed: (i) the occurrence of multiple, distinct evolutionary lineages within each single peninsula [15,16,17,18,19,20,21], and (ii) remarkable variation in the levels of genetic diversity among populations [14, 22,23,24,25]. Still, several studies have shown that the spatial resolution of a hotspot is smaller than previously believed, as in these areas the species genetic structure can be very scattered [10, 20]. Therefore, the single peninsula cannot be treated as a unique, homogenous hotspot, especially in a conservation genetic perspective, where the proper definition of evolutionary and management units is mandatory to successfully define conservation plannings [26]. This evidence stresses the need for more detailed data on the geographic structure of genetic variation of threatened/endangered taxa, in order to set proper conservation programs accounting for population evolutionary history and genetic diversity levels.

In this study, we focus on the geographic patterns of genetic variation of an endangered tortoise species, the Hermann’s tortoise Testudo hermanni Gmelin, 1789. The Hermann’s tortoise is a land tortoise inhabiting coastal and sub-mountain regions of the Italian and Balkan peninsulas, some Mediterranean islands, and two small areas of eastern Spain and southern France [27]. Hermann’s tortoise populations have declined markedly in the last 40 years, mainly due to coastal habitat destruction (including fires) and overharvesting for pet-trades [28,29,30]. Genetic investigations on mitochondrial markers identified two main genetic clusters, corresponding to the two recognized subspecies: T. h. hermanni, which includes populations from the Italian Peninsula, the western Mediterranean islands, and the isolated Spanish and French populations; T. h. boettgeri, which includes populations from the Balkan Peninsula and northern Italy [31]. Further genetic sub-structuring has been recognized between the eastern and the western Balkan populations, between the Italian and the western Mediterranean populations - although western Mediterranean populations were considered as introduced by humans during the Neolithic - and among populations within the Italian Peninsula [32]. However, previous studies did not analyse population structure and genetic diversity in southern Italy, a region usually harbouring multiple evolutionary lineages and non-negligible levels of diversity in many terrestrial vertebrates [14, 18, 20]. Interestingly, a recent study aimed at identifying the geographic origin of confiscated tortoises using microsatellite loci revealed some level of genetic differentiation in Hermann’s tortoise individuals from southern Italy [30, 33].

Here we employed a set of seven microsatellite loci to analyse the fine-scale population genetic structure of the Hermann’s tortoise in the Italian Peninsula, with a particular effort on the southernmost part of its range. We aimed to (i) providing further insights on the Hermann’s tortoise genetic structure and evolutionary history within the Italian Peninsula, (ii) identifying its hotspot of genetic diversity, (iii) mapping the geographic distribution of its evolutionary lineages. Considering the ongoing reduction of the Hermann’s tortoise habitat, this information will provide fundamental new insights to define effective strategies for the conservation of its populations.


The dataset consisted of a multi-locus genotype for 174 wild Hermann’s tortoise individuals (see Table 1) at seven microsatellite loci, with 10.3% of missing data (available on ZENODO repository following the link The inspection of null allele occurrences with Micro-Checker [34] revealed the possible null alleles at locus Gal236 in populations 8, 11, 13 and 14, and at Test10 in populations 3. However, no significant deviation from the Hardy-Weinberg and linkage equilibria was found after the Bonferroni correction was applied in FSTAT [35].

Table 1 Sample number, collecting locality, geographic coordinates, sample size (analysed specimens), allelic richness (Ar) and expected heterozygosity (He)

Allelic richness and the unbiased mean expected heterozygosity estimated for each population using the adegenet R package [36] are shown in Table 1. Population 20 (Roggiano Gravina), 23 (Santa Rania) and 24 (Roccabernarda) from the north-central Calabrian region showed the highest values of genetic diversity (both heterozygosity and allelic richness), whereas the lowest values of genetic diversity were observed in population 10 (Baia Domizia, central Italy).

The spatial Principal Component Analysis (sPCA) performed by the adegenet R package [36] revealed a significant geographic structure in the analysed populations. The global permutation test on the eigenvalues showed significant global structure (max(t) = 0.038, P < 0.001), but did not show any significant local structure (max(t) = 0.0189, P = 0.131). The scree plots of the eigenvectors suggest that the first 2 global axes were the most informative axes and were thus retained. For each sample, the scores of the first two principal components were plotted on a distribution map (Fig. 1). The first PC identify two groups, one ranging from Sicily to the north of Calabrian region, and one ranging from the north of Calabrian region to the northern peninsula. The second PC separated the Sicilian samples from the Calabrian samples. Last, we plotted the lagged scores of the first two axes together in a single colour plot, so that individuals that are closely related in the multivariate space are more similar in colour. With this approach, the three main clusters are clearly distinguished on map (Fig. 1c).

Fig. 1
figure 1

Results from the spatial principal component analysis (sPCA) on the Hermann’s tortoise Testudo hermanni in the Italian Peninsula. (a-b) projection of the individual scores from the first two principal components on a distribution map; the color of the boxes (black or white) corresponds to the sign of the score and the area is proportional to the absolute value of the score; therefore large circles indicate large differentiation with smaller circles indicating smaller differences; (c) lagged scores of the first two axes together in a single colour plot, so that individuals that are closely related in the multivariate space are more similar in colour; (d) barplot showing eigenvalues, positive eigenvalues (red, on the left) correspond to global structures, while negative eigenvalues (blues, on the right) indicate local patterns

The Bayesian clustering analyses carried out with TESS 2.3.1 [37, 38] showed a geographic structuring of genetic variation consistent with results from the sPCA. The plots of the DIC values (i.e. deviance information criterion) versus K values (i.e. the number of clusters) reached a plateau at K = 3 and only a minor decrease in the DIC values was observed at higher K values. Furthermore, the inspection of the plotted membership coefficients for higher K values (see Additional File 1) showed that only three meaningful clusters were represented [39]. Bar-plots showing the individual admixture proportions and pie-charts showing the average proportion of each cluster within each sampled population are given in Fig. 2. The spatial distribution of the three clusters showed a clear geographical structure: one cluster extends in Sicily, one cluster ranges from the south to the north of Calabria, and the third cluster ranges from the north of Calabria to the rest of the peninsula. Large genetic admixture is observed in individuals from the north of Calabria (pop. 19, 20 and 21), as well as in those from the south of Calabrian region (pop. 27).

Fig. 2
figure 2

Genetic structure of Italian populations of Testudo hermanni at 7 microsatellite loci estimated using TESS. The bar plot on the left shows the admixture proportions of each individual for the three genetic clusters recovered; the pie diagrams on the maps show the frequency distributions of each cluster among the populations; circle size is proportional to the sample size; the line chart shows the mean values of the DIC statistics (averaged over 100 runs) for the number of genetic clusters (K) ranging from 2 to 9. The map was drawn using the software Canvas 11 (ACD Systems of America, Inc.)

Finally, results from the standard Principal Component Analysis (PCA) performed on individual genotypes showed a shallow but clear individual clustering in three main groups, which are consistent with the results from previous analyses (see Fig. 3). Indeed, the PC1 allows to easily separate northern from southern individuals, while PC2 show a shallow demarcation between Sicilian and Calabrian individuals.

Fig. 3
figure 3

Principal Component Analysis (PCA) on individual genotypes for 7 microsatellites loci. Dots represent individuals, lines connect individuals from a same group. Colors represent the different clusters identified by the TESS analysis (cf. Figure 2)

The analysis of molecular variance (AMOVA) in ARLEQUIN [40] was performed on the groupings resulting from the previous clustering analyses: Apennine group, pop. 1–19; Calabrian group, pop. 20–27; Sicilian group, pop. 28–31. With this grouping, 14.36% of variation had been attributed to the among-group level (FCT: 0.14), 6.57% to the among-population within groups level (FSC: 0.08), whereas the most of variation was attributed to the within population level (79.05%, FST: 0.20); all variance components and fixation indices were statistically significant (P < 0.001). The genetic differentiation among the three clusters, expressed as pairwise Fst estimation, was of the same magnitude, spanning from 0.153 to 0.155 (population pairwise Fst estimates are resumed in the Additional File 2). Estimates of genetic diversity for the three genetic clusters were as follow: Apennine cluster, Ar: 8.9, He: 0.56 (sd 0.21), Ho: 0.46 (sd 0.16); Calabrian cluster: Ar: 8.0, He: 0.59 (sd. 0.24), Ho: 0.54 (0.25); Sicilian cluster: Ar: 5.4, He 0.56 (sd 0.31) Ho 0.51 (sd 0.32).


The analysis of genetic variation of Testudo hermanni in the Italian Peninsula and Sicily showed the existence of three main genetic clusters. Previous studies identified some genetic differentiation between populations from the Italian Peninsula and Sicily, and among populations within the peninsula [30, 32]. However, the analysis of the geographic patterns of genetic variation applied here allowed us to identify marked sub-structuring within the southernmost region of the Italian Peninsula, and to define the geographic distributions of the distinct genetic lineages. In fact, we identified one genetic cluster spread in Sicily, one cluster spread from the Aspromonte massif to the Sila mountain chain, and one cluster spread throughout the rest of the Italian Peninsula. Furthermore, we identified two restricted areas of genetic admixture, one in the south of the Aspromonte massif, and one in the north of Calabria, corresponding to the Pollino massif. Overall, our results unveil a hotspot of genetic diversity for Testudo hermanni in southern Italy, in an area spanning from the Pollino massif to the Aspromonte massif, and suggest that the interplay between high topographic complexity and Pleistocene climate changes in this region triggered the formation of this hotspot.

The occurrence of genetic sub-structuring and distinct genetic clusters within the Italian Peninsula has been observed in several temperate species including amphibians [16, 17, 22, 23, 25], mammals [14, 18, 20, 41] and reptiles [21], among others [12]. In most of these species, the Calabrian region was identified as a hotspot of genetic lineages [16, 20,21,22,23]. Within this region, major mountain areas are arranged along the north-south axis and are separated by lowland fluvial valleys. This topographic structure led glacio-eustatic sea-level oscillations to turn mountain massifs into paleo-islands [42,43,44,45,46,47]. In particular, these dynamics repeatedly insularized the Sila and Aspromonte massifs, heavily affecting the population structure of terrestrial animal species inhabiting these areas, and leaving detectable imprints in their current genetic structure [14, 16, 22, 23]. This scenario is concordant also with the genetic structure that we identified for the Hermann’s tortoise in southern Italy. Indeed, looking at the distribution of the three genetic clusters, it is possible to identify at least three putative areas acting as Pleistocene refugia for Hermann’s tortoise populations: one located in the southern part of the Calabrian region, one located somewhere north of this area, and one located in Sicily. In the absence of molecular dating analyses, our hypothesis on the Pleistocene history Hermann’s tortoise populations in southern Italy should be taken with caution. However, support for this hypothesis comes also from the fossil record, which identified sites where populations survived during the Late Pleistocene to be located mainly in the southern part of the peninsula, between the Campania and the Calabrian regions [48]. Further investigations involving a higher number of markers (e.g. SNPs and/or nuclear sequence markers) could shed light on the level of divergence among populations as well as on the demographic and evolutionary histories of the three lineages, opening for studies on local adaptation.

Our data clearly identified two areas of genetic admixture, one located in the northern and one in the southern edge of the Calabrian region. These areas, likely originating from secondary contacts among distinct lineages, closely match with areas of secondary contact and admixture observed in several other taxa [14, 22, 23, 49]. Within these areas, gene flow between differentiated lineages boosted the level of population genetic diversity, leading to the comparatively high values of both heterozygosity and allelic richness observed (see Table 1). As a consequence, the whole Calabrian region emerge as a structured hotspot of intraspecific genetic variation for the Hermann’s tortoise, where both unique lineages and high levels of population genetic diversity are found.

Our results have remarkable implications for the management of the Hermann’s tortoise populations. Because of the widespread population decline, T. hermanni is considered as Near Threatened by the IUCN red list of threatened species at a global scale [50] and as Endangered in Italy [51]. However, we identified three unique evolutionarily significant units [52], two of them with narrow and endemic ranges. Assessments of their demographic consistence, as well as of the current threats to their populations have to be planned in the near future, in order to integrate the genetic information into the regional strategies for biodiversity conservation. Furthermore, identifying a hotspot of intraspecific genetic variation in a previously under-investigated region claims for a more detailed investigation on the status of populations inhabiting the hotspot, which might represent a valuable resource for the conservation and management of this species [5]. Indeed, intraspecific genetic variation provides populations with the potential to adapt to the ongoing changes in their biotic and abiotic environment [53, 54]. At the same time, because of the link between genetic diversity and effective population size, these populations are less likely to be affected by the detrimental consequences of genetic drift and inbreeding depression [5, 55,56,57]. Thus, Hermann’s tortoise populations from the southern part of the Italian Peninsula clearly represent a conservation priority for this species. Finally, the sharp genetic structure identified here, which define the proper geographic distribution of the distinct management units, provides valid support for more informed relocation programs of confiscated animals in the wild [30].


This study emphasizes the importance of integrating a multi-marker approach with a thorough sampling scheme in investigating the geographic structure of intraspecific genetic variation, in order to unveil hidden patterns of fine-scale genetic structuring. By applying this approach, we identified sharp genetic structure in the Hermann’s tortoise populations from southern Italy. Still, our results add further evidence to the role of Calabrian region as a hotspot of biodiversity, and claim for an assessment of the population genetic structure in other supposed well-known taxa inhabiting this region, with special attention on threatened taxa.


We investigated the geographic structure of genetic variation of Testudo h. hermanni within the Italian Peninsula and Sicily, by genotyping 174 wild individuals at seven microsatellite loci (Test10, Test56, Test71, Test76, Gal136, Gal75, and Gal263), following the protocols described in Biello et al. [30].

Collecting sites and sample sizes are given in Table 1. Samples were collected from 2010 to 2018, mainly from March to June and from September to October, i.e. the period of maximum species activity. Considering the low dispersal capacity of this species [58], and the fragmentation of the species habitat in Italy [59], each collection site can be considered a distinct breeding population. Genomic DNA was obtained from blood samples collected from the nape or coccygeal vein. About 75 µL of blood were spotted on FTA® Classic Cards (Whatman™, GE Healthcare) and stored at room temperature. Alternatively, whole blood samples (100 µL – 1 mL) were treated with K3-EDTA and stored at -20 °C. DNA was extracted from both FTA-Cards and whole blood samples using a suspension of 5% Chelex® 100 Resin. Field works, collection of tissues, and the experimental protocols were performed in accordance with the relevant guidelines and regulations (including ethics guidelines and regulations), and were approved by the Italian Ministry of Environment, Land and Sea Protection (permit codes: 0044068–4/12/2012-PNM-II; 0001805/PNM − 4/2/2015; ISPRA 68,754/T-A31–28/11/2016).

Fragment analysis of PCR products was performed by Macrogen Inc. on an ABI 3730xl Genetic Analyser (Applied Biosystems) with a 400HD size standard. Electropherograms checking and allele calling was performed by using GENEMAPPER® 4.1. Micro-Checker 2.2.3 [34] was used to test for null alleles and large-allele dropout influences. Allelic frequencies were then computed by using GENETIX 4.05 [60], while FSTAT [35] was used to test for deviations from the expected Hardy-Weinberg and linkage equilibria, as well as to estimate population pairwise Fst values. We also estimated the mean allelic richness and the mean observed and expected (unbiased) heterozygosity using the adegenet package in R environment [36]; allelic richness was computed using the rarefaction method [61].

The extent of spatial patterns of genetic variation was investigated by a spatial Principal Components Analysis (sPCA), as implemented in the adegenet R package [36]. The analysis takes into account the spatial autocorrelation (Moran’s I) of genetic data to estimate coefficients of similarity among individuals. Jittered geographic coordinates of sampling location for each individual were used to build two different networks, a Delaunay triangulation and a Neighbourhood by distance network. Since no significant differences were observed in preliminary results, we only retained the Neighbourhood by distance network, which best represents realistic connections among the sampled populations (i.e. closer populations are more likely to be connected than farther populations). The significance of global (neighbouring individuals are more similar than expected) and local (neighbouring individuals are more dissimilar than expected) structures of genetic variation was assessed with a Monte Carlo-based test (9999 iterations). Results from the first two principal components were mapped in geographic space and examined to cluster samples into genetically and geographically distinct populations.

The genetic clustering of individuals across the study area was also investigated using the Bayesian clustering algorithm implemented in TESS 2.3.1 and the geographical location of individuals as prior information [37, 38]. The analysis was performed by modelling admixture using a conditional autoregressive model (CAR). Preliminary analyses were carried out to assess model performance, with 20 000 steps (the first 5 000 were discarded as burn-in) and 10 replicates for each K value (i.e. the number of clusters) between 2 and 31. The final analysis contained 100 replicates for each K value, with K = 2–9; each run consisted of 80 000 steps, with the first 30 000 discarded as burn-in. The spatial interaction parameter was initially kept at the default value (0.6), and the updating option was activated. The model that best fitted the data was selected using the deviance information criterion (DIC). DIC values were averaged over the 100 replicates for each K value, and the most probable K value was selected as the one at which the average DIC reached a plateau. For the selected K value, the estimated admixture proportions of the 10 runs with the lowest DIC were averaged using CLUMPP 1.1.2 [62].

Finally, the genetic clustering of individuals was also evaluated without accounting for spatial information by means of a standard PCA of genotypic data, as implemented in the adegenet R package [36].

In order to estimate the amount of variation attributable to differences among the clusters identified by the clustering methods, among populations within clusters, and within populations, a locus by locus analysis of molecular variance (AMOVA) was performed using ARLEQUIN [40]. Populations resulting as admixed were attributed to the most represented cluster. The significance of variance components and fixation indices was tested using 1092 permutations. Using this grouping, we estimated the pairwise Fst among group, the mean allelic richness, and the mean observed and expected (unbiased) heterozygosity within the main clusters using the adegenet R package.

Data availability

The datasets generated and/or analysed during the current study are available in the ZENODO repository at the following link:


  1. Petit RJ, Aguinagalde I, de Beaulieu J-L, Bittkau C, Brewer S, Cheddadi R, et al. Glacial refugia: hotspots but not melting pots of genetic diversity. Science. 2003;300:1563–5.

    Article  CAS  PubMed  Google Scholar 

  2. Hewitt GM. The structure of biodiversity – insights from molecular phylogeography. Front Zool. 2004;1:4.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Zachos FE, Habel JC. Biodiversity hotspots: distribution and protection of conservation priority areas. Springer Science & Business Media; 2011.

  4. Lande R, Shannon S. The role of genetic variation in adaptation and population persistence in a changing environment. Evolution. 1996;50:434.

    Article  PubMed  Google Scholar 

  5. Kardos M, Armstrong EE, Fitzpatrick SW, Hauser S, Hedrick PW, Miller JM, et al. The crucial role of genome-wide genetic variation in conservation. Proc Natl Acad Sci USA. 2021;118:e2104642118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Hampe A, Petit RJ. Conserving biodiversity under climate change: the rear edge matters: Rear edges and climate change. Ecol Lett. 2005;8:461–7.

    Article  PubMed  Google Scholar 

  7. Willis KJ, Birks HJB. What is natural? The need for a long-term perspective in biodiversity conservation. Science. 2006;314:1261–5.

    Article  CAS  PubMed  Google Scholar 

  8. Vandergast AG, Bohonak AJ, Hathaway SA, Boys J, Fisher RN. Are hotspots of evolutionary potential adequately protected in southern California? Biol Conserv. 2008;141:1648–64.

    Article  Google Scholar 

  9. Cheng E, Hodges KE, Melo-Ferreira J, Alves PC, Mills LS. Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare. Mol Ecol. 2014;23:2929–42.

    Article  CAS  PubMed  Google Scholar 

  10. Zampiglia M, Bisconti R, Maiorano L, Aloise G, Siclari A, Pellegrino F, et al. Drilling down hotspots of intraspecific diversity to bring them into on-ground conservation of threatened species. Front Ecol Evol. 2019;7:205.

    Article  Google Scholar 

  11. Hewitt G. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13.

    Article  CAS  PubMed  Google Scholar 

  12. Hewitt GM. Mediterranean peninsulas: the evolution of hotspots. In: Zachos FE, Habel JC, editors. Biodiversity Hotspots: Distribution and Protection of Conservation Priority Areas. Berlin: Springer; 2011. pp. 123–47.

    Chapter  Google Scholar 

  13. Gómez A, Lunt DH. Refugia within refugia: patterns of phylogeographic concordance in the iberian peninsula. In: Weiss S, Ferrand N, editors. Phylogeography of Southern European Refugia: Evolutionary perspectives on the origins and conservation of European biodiversity. Dordrecht: Springer Netherlands; 2007. pp. 155–88.

    Chapter  Google Scholar 

  14. Canestrelli D, Aloise G, Cecchetti S, Nascetti G. Birth of a hotspot of intraspecific genetic diversity: notes from the underground. Mol Ecol. 2010;19:5432–51.

    Article  PubMed  Google Scholar 

  15. Miraldo A, Hewitt GM, Paulo OS, Emerson BC. Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones. BMC Evol Biol. 2011;11:170.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Canestrelli D, Sacco F, Nascetti G. On glacial refugia, genetic diversity, and microevolutionary processes: Deep phylogeographical structure in the endemic newt Lissotriton italicus. Biol J Linn Soc 2012: 105(1); 42–55.

  17. Canestrelli D, Salvi D, Maura M, Bologna MA, Nascetti G. One species three Pleistocene evolutionary histories: phylogeography of the Italian crested newt, Triturus carnifex. PLoS ONE. 2012;7:e41754.

  18. Bisconti R, Aloise G, Siclari A, Fava V, Provenzano M, Arduino P, et al. Forest Dormouse (Dryomys nitedula) populations in southern Italy belong to a deeply divergent evolutionary lineage: implications for taxonomy and conservation. Hystrix. 2018;29(1):75–9. doi:

    Article  Google Scholar 

  19. Querejeta M, Castresana J. Evolutionary history of the endemic water shrew Neomys anomalus: Recurrent phylogeographic patterns in semi-aquatic mammals of the Iberian Peninsula. Ecol Evol. 2018;8:10138–46.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Chiocchio A, Colangelo P, Aloise G, Amori G, Bertolino S, Bisconti R, et al. Population genetic structure of the bank vole Myodes glareolus within its glacial refugium in peninsular Italy. J Zool Syst Evol Res. 2019;57:959–69.

    Article  Google Scholar 

  21. Schultze N, Spitzweg C, Corti C, Delaugerre M, Di Nicola MR, Geniez P, et al. Mitochondrial ghost lineages blur phylogeography and taxonomy of Natrix helvetica and N. natrix in Italy and Corsica. Zool Scr. 2020;49:395–411.

    Article  Google Scholar 

  22. Canestrelli D, Cimmaruta R, Costantini V, Nascetti G. Genetic diversity and phylogeography of the Apennine yellow-bellied toad Bombina pachypus, with implications for conservation. Mol Ecol. 2006;15:3741–54.

    Article  CAS  PubMed  Google Scholar 

  23. Canestrelli D, Cimmaruta R, Nascetti G. Population genetic structure and diversity of the Apennine endemic stream frog, Rana italica - insights on the Pleistocene evolutionary history of the Italian peninsular biota. Mol Ecol. 2008;17:3856–72.

    Article  PubMed  Google Scholar 

  24. Chiocchio A, Bisconti R, Zampiglia M, Nascetti G, Canestrelli D. Quaternary history, population genetic structure and diversity of the cold-adapted Alpine newt Ichthyosaura alpestris in peninsular Italy. Sci Rep. 2017;7:2955.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Chiocchio A, Arntzen JW, Martínez-Solano I, de Vries W, Bisconti R, Pezzarossa A, et al. Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo). Sci Rep. 2021;11:260.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Moritz CC, Potter S. The importance of an evolutionary perspective in conservation policy planning. Mol Ecol. 2013;22:5969–71.

    Article  PubMed  Google Scholar 

  27. Sillero N, Campos J, Bonardi A, Corti C, Creemers R, Crochet P-A, et al. Updated distribution and biogeography of amphibians and reptiles of Europe. Amphibia-Reptilia. 2014;35:1–31.

    Article  Google Scholar 

  28. Stubbs D, Swingland IR. The ecology of a Mediterranean tortoise (Testudo hermanni): a declining population. Can J Zool. 2011.

    Article  Google Scholar 

  29. Couturier T, Besnard A, Bertolero A, Bosc V, Astruc G, Cheylan M. Factors determining the abundance and occurrence of Hermann’s tortoise Testudo hermanni in France and Spain: Fire regime and landscape changes as the main drivers. Biol Conserv. 2014;170:177–87.

    Article  Google Scholar 

  30. Biello R, Zampiglia M, Corti C, Deli G, Biaggini M, Crestanello B, et al. Mapping the geographic origin of captive and confiscated Hermann’s tortoises: A genetic toolkit for conservation and forensic analyses. Forensic Sci Int-Gen. 2021;51:102447.

    Article  CAS  Google Scholar 

  31. Fritz U, Auer M, Bertolero A, Cheylan M, Fattizzo T, Hundsdorfer AK, et al. A rangewide phylogeography of Hermann’s tortoise, Testudo hermanni (Reptilia: Testudines: Testudinidae): implications for taxonomy. Zool Scripta. 2006;35:531–43.

    Article  Google Scholar 

  32. Perez M, Livoreil B, Mantovani S, Boisselier M-C, Crestanello B, Abdelkrim J, et al. Genetic Variation and Population Structure in the Endangered Hermann’s Tortoise: The Roles of Geography and Human-Mediated Processes. J Hered. 2014;105:70–81.

    Article  PubMed  Google Scholar 

  33. AA.VV. Linee guida per il ricollocamento in natura di individui del genere Testudo spp. confiscati. Ministero dell’Ambiente e della Tutela del Territorio e del Mare. Museo di Storia Naturale dell’ Università degli Studi di Firenze. Università degli Studi di Perugia. 2019.

  34. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

    Article  Google Scholar 

  35. Goudet J. FSTAT (Version 1.2): A Computer Program to Calculate F-Statistics. J Hered. 1995;86:485–6.

    Article  Google Scholar 

  36. Jombart T. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5.

    Article  CAS  PubMed  Google Scholar 

  37. Chen C, Durand E, Forbes F, François O. Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study. Mol Ecol Notes. 2007;7:747–56.

    Article  Google Scholar 

  38. François O, Durand E. Spatially explicit Bayesian clustering models in population genetics: SPATIAL CLUSTERING MODELS. Mol Ecol Res. 2010;10:773–84.

    Article  Google Scholar 

  39. Durand E, Jay F, Gaggiotti OE, Francois O. Spatial Inference of Admixture Proportions and Secondary Contact Zones. Mol Biol Evol. 2009;26:1963–73.

    Article  CAS  PubMed  Google Scholar 

  40. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010;10:564–7.

    Article  Google Scholar 

  41. Vega R, Amori G, Aloise G, Cellini S, Loy A, Searle JB. Genetic and morphological variation in a Mediterranean glacial refugium: evidence from Italian pygmy shrews, Sorex minutus (Mammalia: Soricomorpha). Biol J Linn Soc. 2010;100:774–87.

    Article  Google Scholar 

  42. Ghisetti F. L’evoluzione strutturale del bacino pliopleistocenico di Reggio Calabria nel quadro geodinamico dell’arco calabro. Bollettino della Società Geologica Italiana. 1981;100:433–66.

    Google Scholar 

  43. Caloi L, Malatesta A, Palombo MR. Biogeografia della Calabria meridionale durante il quaternario. Accad Peloritana dei Pericolanti Classe FFMMNN. 1989;67:261–78.

    Google Scholar 

  44. Tortorici L, Monaco C, Tansi C, Cocina O. Recent and active tectonics in the Calabrian arc (Southern Italy). Tectonophysics. 1995;243:37–55.

    Article  Google Scholar 

  45. Bonfiglio L, Mangano G, Marra AC, Masini F, Pavia M, Petruso D. Pleistocene Calabrian and Sicilian bioprovinces. Geobios. 2002;35:29–39.

    Article  Google Scholar 

  46. Cucci L. Raised marine terraces in the Northern Calabrian Arc (Southern Italy): a ~ 600 kyr-long geological record of regional uplift. 2004.

  47. Tansi C, Muto F, Critelli S, Iovine G. Neogene-Quaternary strike-slip tectonics in the central Calabrian Arc (southern Italy). J Geodyn. 2007;43:393–414.

    Article  Google Scholar 

  48. Morales Pérez JV, Serra AS. The Quaternary fossil record of the genus Testudo in the Iberian Peninsula. Archaeological implications and diachronic distribution in the western Mediterranean. J Archaeol Sci. 2009;36:1152–62.

    Article  Google Scholar 

  49. Scalercio S, Cini A, Menchetti M, Vodă R, Bonelli S, Bordoni A, et al. How long is 3 km for a butterfly? Ecological constraints and functional traits explain high mitochondrial genetic diversity between Sicily and the Italian Peninsula. J Anim Ecol. 2020;89:2013–26.

    Article  PubMed  Google Scholar 

  50. Amori G, Hutterer R, Kryštufek B, Yigit N, Mitsain G, Palomo LJ, et al. Testudo hermanni. The IUCN Red List of Threatened Species 2008: e.T4973A11104409.

  51. Rondinini C, Battistoni A, Peronace V, Teofli C. Lista Rossa IUCN dei Vertebrati Italiani. Roma: Comitato Italiano IUCN e Ministero dell’Ambiente e della Tutela del Territorio e del Mare; 2013.

    Google Scholar 

  52. Moritz C. Defining ‘Evolutionarily Significant Units’ for conservation. Trends Ecol Evol. 1994;9:373–5.

    Article  CAS  PubMed  Google Scholar 

  53. Hoffmann AA, Sgrò CM. Climate change and evolutionary adaptation. Nature. 2011;470:479–85.

    Article  CAS  PubMed  Google Scholar 

  54. Alberto FJ, Aitken SN, Alía R, González-Martínez SC, Hänninen H, Kremer A, et al. Potential for evolutionary responses to climate change – evidence from tree populations. Glob Change Biol. 2013;19:1645–61.

    Article  Google Scholar 

  55. Allendorf FW, Funk WC, Aitken SN, Byrne M, Luikart G. Conservation and the Genomics of Populations. Oxford University Press; 2022.

  56. Frankham R. Challenges and opportunities of genetic approaches to biological conservation. Biol Conserv. 2010;143:1919–27.

    Article  Google Scholar 

  57. Lanfear R, Kokko H, Eyre-Walker A. Population size and the rate of evolution. TREE. 2014;29:33–41.

    PubMed  Google Scholar 

  58. Fasola M, Mazzotti S. & Pisapia A Activity and home range of Testudo hermanni in Northern Italy. Amph Rept. 2002;23(3):305–12.

    Article  Google Scholar 

  59. Corti C, Bassu L, Biaggini M, Bressi N, Capula M, Di Cerbo AR, et al Updated distribution of Testudo hermanni hermanni in Italy. In: Caron, S, editor, Proceedings of the International workshop on the management and restoration of Hermann’s tortoise habitats and populations. Atelier international sur la gestion et la restauration des populations et habitats de la Tortue d’Hermann. International workshop on the management and restoration of Hermann’s tortoise populations and habitats. Soptom, Gonfaron, France, 2013 (pp. 28–33).

  60. Belkhir K. GENETIX 4.05, logiciel sous Windows TM pour la genetique des populations. 2004.

  61. Petit RJ, El Mousadik A, Pons O. Identifying Populations for Conservation on the Basis of Genetic Markers. Conserv Biol. 2008;12:844–55.

    Article  Google Scholar 

  62. Jakobsson M, Rosenberg NA. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007;23:1801–6.

    Article  CAS  PubMed  Google Scholar 

Download references


We are grateful to Francesco Pellegrino, Gianluca Deli, Antonio Romano and Sandro Tripepi for their support on fieldwork.


This work was supported by the Italian Ministry for Environment, Land and Sea Protection (Direzione generale per la Protezione della Natura e del Mare).

Author information

Authors and Affiliations



DC and CC conceived and designed the study. CC, MB, AC, DC, DTL, LFL, ES, OO, MTM performed fieldwork and provided samples. MZ, AC and RB performed laboratories experiments. AC and MZ analysed the data. AC drafted the manuscript, with contributions from DC. All the Authors read and approved the manuscript.

Corresponding author

Correspondence to Andrea Chiocchio.

Ethics declarations

Ethics approval and consent to participate

All experimental protocols were performed in accordance with the relevant guidelines and regulations, including ethics guidelines and regulations, and were approved by the Italian Ministry of Environment, Land and Sea Protection (permit codes: 0044068–4/12/2012-PNM-II; 0001805/PNM − 4/2/2015; ISPRA 68754/T-A31–28/11/2016). All methods are reported in accordance with ARRIVE guidelines ( for the reporting of animal experiments.

Consent for publication

Not applicable.

Competing interests

The Authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chiocchio, A., Zampiglia, M., Biaggini, M. et al. Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann’s tortoise Testudo hermanni. BMC Ecol Evo 22, 131 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Biodiversity hotspots
  • Conservation genetics
  • Italian Peninsula
  • Genetic structure
  • Phylogeography
  • Threatened species