Mitochondrial ghost lineages blur phylogeography and taxonomy of Natrix helvetica and N. natrix in Italy and Corsica

Genetic identification, morphology and distribution of Natrixhelvetica subspecies in southern and western Switzerland (Reptilia, Squamata, Serpentes)

Most of Switzerland is inhabited by the nominotypical subspecies of the barred grass snake (Natrixhelveticahelvetica), which is characterized by mitochondrial DNA lineage E. Only in the northeast of the country, the common grass snake (N.natrix) occurs and hybridizes with N.h.helvetica in a narrow contact zone. However, we discovered that in southern and western Switzerland barred grass snakes representing another mtDNA lineage (lineage C) are widely distributed. Lineage C is typical for Alpine populations of the southern subspecies N.h.sicula. Our microsatellite analyses of the Swiss samples revealed differences between the two subspecies and also a substructure with two clusters in each subspecies. Furthermore, we discovered a contact and hybrid zone of N.h.helvetica and N.h.sicula along the northern shore of Lake Geneva and also confirm that interbreeding with alien common grass snakes (N.n.moreotica, mtDNA lineage 7) occurs there. This finding is of concern for nature conservation and measures should be taken to prevent further genetic pollution. Using morphometrics, we found no differences between the two subspecies of N.helvetica, while N.natrix was slightly distinct from N.helvetica.

Natrix natrix after dark: citizen science sheds light on the common grass snake's nightlife

Activity patterns in animals are often species-specific, and can be generally categorized as diurnal, crepuscular, or nocturnal. Understanding these patterns provides insight into ecological adaptations and behaviors. The common grass snake (Natrix natrix), one of the most common and widespread European snake species, is traditionally considered diurnal, with scarce evidence of its crepuscular and nocturnal activity. We aimed to document the distribution, environmental conditions, and potential phenotype associations of nighttime activity in N. natrix. We used citizen science data from iNaturalist (1992-2022), Observation.org (2012-2022), together with personal field observations (2010-2023) to collect 127 crepuscular and nocturnal activity records. Most observations occurred between May and August, coinciding with the peak activity period of grass snakes across their distribution range. Statistical analyses revealed no significant difference in mean daily temperatures between crepuscular and nocturnal observations. However, striped individuals displayed nocturnal activity at higher temperatures, consistent with their distribution in warmer regions, but failed to register any difference when tested on a geographic subsample, that accounted for sympatry of the phenotypes. Surprisingly, we found no significant impact of moon presence or moonlight on nighttime activity or age class, contrary to expectations based on other snake species' responses. While our study reveals that nocturnal activity in the common grass snake is geographically widespread, further research is warranted to understand its drivers and ecological implications. This study highlights the value of citizen science platforms for biological and ecological research, offering unparalleled spatial and temporal coverage by their users. In conclusion, our work extends the knowledge of nocturnal behavior in N. natrix and underlines the critical role of citizen science in discovering behavioral aspects of common and widespread species.

DNA Barcoding in Species Delimitation: From Genetic Distances to Integrative Taxonomy

Over the past two decades, DNA barcoding has become the most popular exploration approach in molecular taxonomy, whether for identification, discovery, delimitation, or description of species. The present contribution focuses on the utility of DNA barcoding for taxonomic research activities related to species delimitation, emphasizing the following aspects:(1) To what extent DNA barcoding can be a valuable ally for fundamental taxonomic research, (2) its methodological and theoretical limitations, (3) the conceptual background and practical use of pairwise distances between DNA barcode sequences in taxonomy, and (4) the different ways in which DNA barcoding can be combined with complementary means of investigation within a broader integrative framework. In this chapter, we recall and discuss the key conceptual advances that have led to the so-called renaissance of taxonomy, elaborate a detailed glossary for the terms specific to this discipline (see Glossary in Chap. 35 ), and propose a newly designed step-by-step species delimitation protocol starting from DNA barcode data that includes steps from the preliminary elaboration of an optimal sampling strategy to the final decision-making process which potentially leads to nomenclatural changes.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

The Endangered Sardinian Grass Snake: Distribution Update, Bioclimatic Niche Modelling, Dorsal Pattern Characterisation, and Literature Review

The Sardinian grass snake, Natrix helvetica cetti, is an endangered endemic snake subspecies with a restricted and highly fragmented geographic distribution. Information on its ecology and detailed geographic distribution are scarce and may negatively impact on its conservation status. Therefore, a literature review on its taxonomy, morphology, ecology, and conservation is presented here. Moreover, field records from the authors, citizen science and the existing literature provide an updated geographic distribution highlighting its presence within 13 new and 7 historic 10 × 10 km cells. Bioclimatic niche modelling was then applied to explore patterns of habitat suitability and phenotypic variation within N. h. cetti. The geographic distribution of the species was found to be positively correlated with altitude and precipitation values, whereas temperature showed a negative correlation. Taken together, these outcomes may explain the snake's presence, particularly in eastern Sardinia. In addition, analysis of distribution overlap with the competing viperine snake (N. maura) and the urodeles as possible overlooked trophic resources (Speleomantes spp. and Euproctus platycephalus) showed overlaps of 66% and 79%, respectively. Finally, geographical or bioclimatic correlations did not explain phenotypic variation patterns observed in this highly polymorphic taxon. Perspectives on future research to investigate N. h. cetti's decline and support effective conservation measures are discussed.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.)

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

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

Background

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.

Results

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.

Conclusion

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.

Erythrism in the Eastern Grass Snake, Natrix natrix (Linnaeus, 1758)

We describe the unusual case of erythrism in the Eastern Grass Snake, Natrix natrix. This colour morph is very rare and has not been reported in the literature before. Despite having observed thousands of N. natrix in the field, we personally detected this morph in only three individuals originating from Slovakia, Romania, and mainland Greece, while photos of a fourth individual from a Greek island were provided to us. In addition, a recent study with a large data set from citizen scientists was unable to produce a single reddish Eastern Grass Snake. Such colouration is likewise uncommon in the western members of Grass Snakes (N. helvetica, N. astreptophora), with two examples provided herein. Because the potential biological importance of erythristic colouration is unclear, we encourage other field herpetologists and naturalists to publish their observations of reddish Grass Snakes in the printed literature.

The European Pine Marten Martes martes (Linnaeus, 1758) Is Autochthonous in Sicily and Constitutes a Well-Characterised Major Phylogroup within the Species (Carnivora, Mustelidae)

Simple Summary

The faunal assemblages currently occurring on islands are often a mélange of native and non-native species, and the actual status of some of them is doubtful at present. Since different laws are enforced for native and non-native species, sound knowledge about their status is pivotal for both their management and our understanding of the natural history of the studied areas. In the frame of this work, the Sicilian population of the European pine marten is genetically characterized for the first time, based on mitochondrial DNA sequences. Our results prove that the European pine marten is native on the island, where it arrived during the Pleistocene, and is represented there by a well-differentiated and endemic evolutionary lineage.

Abstract

No molecular data are currently available for the Sicilian populations of the European pine marten Martes martes, thus preventing any sound inference about its native or non-native status on the island, as well as the local phylogeography of the species. In order to investigate these issues, we sequenced two mtDNA markers in road-killed specimens collected in Sicily. Both markers consistently demonstrated the existence of a well-characterised Sicilian clade of the species, which is endemic to the island and constitutes the sister group of a clade including the Mediterranean and Central–North European major phylogroups of the European pine marten. Such evidence supports the autochthony of Martes martes in Sicily and points to a natural Pleistocene colonisation of the island followed by isolation. The occurrence of a, to date undetected, major phylogroup of the species in Sicily calls for the dedicated monitoring of the Sicilian populations of the species in order to preserve this evolutionarily significant unit.

Phylogeography of Sarmarutilus rubilio (Cypriniformes: Leuciscidae): Complex Genetic Structure, Clues to a New Cryptic Species and Further Insights into Roaches Phylogeny

Italy hosts a large number of endemic freshwater fish species due to complex geological events which promoted genetic differentiation and allopatric speciation. Among them, the South European roach Sarmarutilus rubilio inhabits various freshwater environments in three different ichthyogeographic districts. We investigated the genetic diversity of S. rubilio using two different mitochondrial markers (COI and CR), aiming to define its relationship with other similar taxa from the Balkan area and, from a phylogeographic perspective, test the effects of past hydrogeological dynamics of Italian river basins on its genetic structure and demographic history. Our analysis highlighted a marked genetic divergence between S. rubilio and all other roach species and, among Italian samples, revealed the existence of three deeply divergent geographic haplogroups, named A, B and C. Haplogroup C likely corresponds to a new putative cryptic species and is located at the northern border of the South European roach range; haplogroup B is restricted to Southern Italy; and haplogroup A is widespread across the entire range and in some sites it is in co-occurrence with C or B. Their origin is probably related to the tectonic uplifting of the Apuan Alps in the north and of the Colli Albani Volcano in the south during the Pleistocene, which promoted isolation and vicariance followed by secondary contacts.

Biogeography of Italy revisited: genetic lineages confirm major phylogeographic patterns and a pre-Pleistocene origin of its biota

Warm-adapted species survived the cold stages of the past glacial-interglacial cycles in southern European peninsulas and recolonized major parts of Central and Northern Europe in the wake of postglacial warming. However, many of the genetic lineages which differentiated within these refugia predate even the Pleistocene. One of these refugia is the Italian Peninsula with its surrounding islands. In this meta-analysis, we compiled phylogeographic patterns of multiple species across this region. We transformed existing geographic information on 78 animal and plant species (with a total of 471 genetic lineages) within 17 predefined areas into presence/absence matrices. We elaborated three matrices: (i) only old lineages exclusively distinguished by deep splits, (ii) only young lineages distinguished by shallow (i.e. recent) splits, and (iii) presence/absence of the respective species. To infer biogeographic relationships between the predefined areas, we performed bootstrapped neighbour joining cluster analyses on these three matrices. In addition, we reviewed the geological history of Italy to identify causes of the observed biogeographic patterns. We found Sardinia and Corsica to be biogeographically closely linked with each other, and that they diverge strongly from all other regions. Sicily also diverges strongly from all other regions, while the intra-island differentiation was comparatively low. On the Italian mainland, Calabria exhibited the most pronounced biogeographic differentiation, often with several lineages present, resulting from old vicariance events within the region. Furthermore, southern Apulia and the Po Plain with adjoining areas of northern peninsular Italy displayed considerable distinctiveness. Admixture prevailed in the areas between these three regions. The ancient isolation of Sicily, as well as Sardinia plus Corsica, resulted in endemic lineages with only moderate recent exchange with adjacent mainland regions. Pronounced diversification occurs within the Italian Peninsula. The complex tectonic activities, such as shifting (micro)plates, submergence of major parts of peninsular Italy with the genesis of numerous Pliocene islands, in combination with the climatic cycles during the Pleistocene have most likely generated the current biogeographic pattern of species. Immigrations from the Balkan Peninsula into northern Italy partly accounted for the distinctiveness of this region.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

How often do they do it? An in-depth analysis of the hybrid zone of two grass snake species (Natrix astreptophora and Natrix helvetica)

We examined the contact zone of two parapatric species of grass snake (Natrix astreptophora and Natrix helvetica) in southern France. To this end, we used comprehensive sampling, analysed mtDNA sequences and microsatellite loci, and built Species Distribution Models for current and past climatic conditions. The contact zone had established by the mid-Holocene during range expansions from glacial refuges in the Iberian Peninsula (N. astreptophora) and southern or western France (N. helvetica). The contact zone represents a narrow bimodal hybrid zone, with steep genetic transition from one taxon to the other and rare hybridization, supporting species status for N. astreptophora and N. helvetica. Our results suggest that the steepness of the clines is a more robust tool for species delimitation than cline width. In addition, we discovered in western France, beyond the hybrid zone, a remote population of N. helvetica with genetic signatures of hybridization with N. astreptophora, most likely the result of human-mediated long-distance dispersal. For N. helvetica, we identified a southern and a northern population cluster, connected by broad-scale gene flow in a unimodal hybrid zone running across France. This pattern either reflects genetic divergence caused by allopatry in two microrefuges and subsequent secondary contact or introgression of foreign alleles into the southern cluster.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.