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.

Average heterozygosity is frequently used as a proxy for genetic health, and to compare genetic diversity between species and populations. However, this measurement could be misleading if the distribution of heterozygosity across the genome is highly skewed. We investigated this pitfall in methodology using whole-genome sequencing of the adder ( Vipera berus ), a species experiencing dramatic declines in the UK. We find that mean heterozygosity in adders is notably high, exceeding that of other vertebrates typically regarded as genetically diverse. Their genome-wide distribution of heterozygosity, however, approximates a negative exponential distribution, with most genome regions showing extremely low heterozygosity. Modelling approaches show that this pattern is likely to have resulted from a recent, severe bottleneck and fragmentation most likely caused by anthropogenic activity in a previously large, interconnected adder population. Our results highlight that high standing diversity may mask severe genetic erosion when declines are recent and rapid. In such situations, whole-genome sequencing may provide the best option for genetic risk assessment and targeted conservation actions.