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.

Bat dynamics modelling as a tool for conservation management in subterranean environments.

Bat species inhabit subterranean environments (e.g., caves and mines) in small areas with specific microclimatic conditions, during various periods of their life cycle. Bats can be negatively influenced by microclimatic changes within their roosts if optimal habitat patches become unavailable. Therefore, proper management solutions must be applied for the conservation of vulnerable bat populations, especially in show caves. We have pursued an ensemble species distribution modelling approach in subterranean environments to identify sensible patches for bats. Using multi-annual temperature monitoring and bat distribution surveys performed within ten caves and mines, including show caves, we modelled relevant habitat patches for five bat species. The temperature-based variables generated from this approach proved to be effective when processed via species distribution models, which generated optimal validation results, even for bats that were heavily clustered in colonies. Management measures are proposed for each show cave to help in long-time conservation of hibernation and maternity colonies. These measures include creating suitable microclimatic patches within the caves by ecological reconstruction measures, tourist management practices in relation to bats, and show cave fitting recommendations. This approach has never been performed at this scale due to the complex geostatistical challenges involving subterranean environment mapping and can be further used as best practice guidelines for future conservation projects.

Ecological corridors for the amphibians and reptiles in the Natura 2000 sites of Romania.

Landscape heterogeneity and fragmentation are key challenges for biodiversity conservation. As Earth's landscape is increasingly dominated by anthropogenic land use, it is clear that broad-scale systems of nature reserves connected by corridors are needed to enable the dispersal of flora and fauna. The European Union currently supports a continent-wide network of protected areas, the Natura 2000 program, but this program lacks the necessary connectivity component. To examine whether a comprehensive network could be built in order to protect amphibians and reptiles, two taxonomic groups sensitive to environmental changes due to their physiological constrains and low dispersal capacity, we used species' distribution maps, the sites of community interest (SCIs) in Romania, and landscape resistance rasters. Except Vipera ursinii rakosiensis, all amphibians and reptiles had corridors mapped that, when assembled, provided linkages for up to 27 species. Natura 2000 species were not good candidates for umbrella species as these linkages covered only 17% of the corridors for all species. Important Areas for Connectivity were identified in the Carpathian Mountains and along the Danube River, further confirming these regions as hot spots for biodiversity in Europe, where successful linkages are most likely. In the end, while such corridors may not be created just for amphibians and reptiles, they can easily be incorporated into more complex linkages with corridors for more charismatic species, therefore enhancing the corridors' value in terms of quality and structure.

Temperature shapes the costs, benefits and geographic diversification of sexual coloration in a dragonfly.

The environment shapes the evolution of secondary sexual traits by determining how their costs and benefits vary across the landscape. Given the thermal properties of dark coloration generally, temperature should crucially influence the costs, benefits and geographic diversification of many secondary sexual colour patterns. We tested this hypothesis using sexually selected wing coloration in a dragonfly. We find that greater wing coloration heats males - the magnitude of which improves flight performance under cool conditions but dramatically reduces it under warm conditions. In a colder region of the species' range, behavioural observations of a wild population show that these thermal effects translate into greater territorial acquisition on thermally variable days. Finally, geo-referenced photographs taken by citizen scientists reveal that this sexually selected wing coloration is dramatically reduced in the hottest portions of the species' range. Collectively, our results underscore temperature's capacity to promote and constrain the evolution of sexual coloration.

Resolving complex phylogeographic patterns in the Balkan Peninsula using closely related wall-lizard species as a model system.

The Balkan Peninsula constitutes a biodiversity hotspot with high levels of species richness and endemism. The complex geological history of the Balkans in conjunction with the climate evolution are hypothesized as the main drivers generating this biodiversity. We investigated the phylogeography, historical demography, and population structure of closely related wall-lizard species from the Balkan Peninsula and southeastern Europe to better understand diversification processes of species with limited dispersal ability, from Late Miocene to the Holocene. We used several analytical methods integrating genome-wide SNPs (ddRADseq), microsatellites, mitochondrial and nuclear DNA data, as well as species distribution modelling. Phylogenomic analysis resulted in a completely resolved species level phylogeny, population level analyses confirmed the existence of at least two cryptic evolutionary lineages and extensive within species genetic structuring. Divergence time estimations indicated that the Messinian Salinity Crisis played a key role in shaping patterns of species divergence, whereas intraspecific genetic structuring was mainly driven by Pliocene tectonic events and Quaternary climatic oscillations. The present work highlights the effectiveness of utilizing multiple methods and data types coupled with extensive geographic sampling to uncover the evolutionary processes that shaped the species over space and time.

Hidden diversity in the Podarcis tauricus (Sauria, Lacertidae) species subgroup in the light of multilocus phylogeny and species delimitation.

The monophyletic species subgroup of Podarcis tauricus is distributed in the western and southern parts of the Balkans, and includes four species with unresolved and unstudied inter- and intra-specific phylogenetic relationships. Using sequence data from two mitochondrial and three nuclear genes and applying several phylogenetic methods and species delimitation approaches to an extensive dataset, we have reconstructed the phylogeny of the Podarcis wall lizards in the Balkans, and re-investigated the taxonomic status of the P. tauricus species subgroup. Multilocus analyses revealed that the aforementioned subgroup consists of five major clades, with P. melisellensis as its most basal taxon. Monophyly of P. tauricus sensu stricto is not supported, with one of the subspecies (P. t. ionicus) displaying great genetic diversity (hidden diversity or cryptic species). It comprises five, geographically distinct, subclades with genetic distances on the species level. Species delimitation approaches revealed nine species within the P. tauricus species subgroup (P. melisellensis, P. gaigeae, P. milensis, and six in the P. tauricus complex), underlining the necessity of taxonomic re-evaluation. We thus synonymize some previously recognized subspecies in this subgroup, elevate P. t. tauricus and P. g. gaigeae to the species level and suggest a distinct Albanian-Greek clade, provisionally named as the P. ionicus species complex. The latter clade comprises five unconfirmed candidate species that call for comprehensive studies in the future.

A reinvestigation of phylogeny and divergence times of the Ablepharus kitaibelii species complex (Sauria, Scincidae) based on mtDNA and nuDNA genes.

Morphological and DNA data support that the East Mediterranean snake-eyed skink Ablepharus kitaibelii represents a species complex that includes four species A. kitaibelii, A. budaki, A. chernovi, and A. rueppellii, highlighting the need of its taxonomic reevaluation. Here, we used Bayesian and Maximum Likelihood methods to estimate the phylogenetic relationships of all members of the complex based on two mitochondrial (cyt b, 16S rRNA) and two nuclear markers (MC1R, and NKTR) and using Chalcides, Eumeces, and Eutropis as outgroups. The biogeographic history of the complex was also investigated through the application of several phylogeographic (BEAST) and biogeographic (BBM) analyses. Paleogeographic and paleoclimatic data were used to support the inferred phylogeographic patterns. The A. kitaibelli species complex exhibits high genetic diversity, revealing cases of hidden diversity and cases of non-monophyletic species such as A. kitaibelii and A. budaki. Our results indicate that A. pannonicus branches off first and a group that comprises specimens of A. kitaibelli and A. budaki from Kastelorizo Island group (southeast Greece) and southwest Turkey, respectively is differentiated from the rest A. kitaibelli and A. budaki populations and may represent a new species. The estimated divergence times place the origin of the complex in the Middle Miocene (∼16Mya) and the divergence of most currently recognized species in the Late Miocene. The inferred ancestral distribution suggests that the complex originated in Anatolia, supposing that several vicariance and dispersal events that are related with the formation of the Mid-Aegean Trench, the Anatolian Diagonal and the orogenesis of the mountain chains in southern and eastern Anatolia have led to current distribution pattern of A. kitaibelii species complex in the Balkans and Middle East.

A revision of the distribution of sea kraits (Reptilia, Laticauda) with an updated occurrence dataset for ecological and conservation research.

The genus Laticauda (Reptilia: Elapidae), commonly known as sea kraits, comprises eight species of marine amphibious snakes distributed along the shores of the Western Pacific Ocean and the Eastern Indian Ocean. We review the information available on the geographic range of sea kraits and analyze their distribution patterns. Generally, we found that south and south-west of Japan, Philippines Archipelago, parts of Indonesia, and Vanuatu have the highest diversity of sea krait species. Further, we compiled the information available on sea kraits' occurrences from a variety of sources, including museum records, field surveys, and the scientific literature. The final database comprises 694 occurrence records, with Laticauda colubrina having the highest number of records and Laticauda schistorhyncha the lowest. The occurrence records were georeferenced and compiled as a database for each sea krait species. This database can be freely used for future studies.

Going west - A subtropical lineage (Vincetoxicum, Apocynaceae: Asclepiadoideae) expanding into Europe.

Vincetoxicum sensu lato is a tropical lineage comprising two clades that have reached high northern latitudes. Of the temperate clades, one is restricted to the Far East, the other one (Vincetoxicum s. str. Clade) extends into Europe, but their ranges overlap in Central China and Japan. Three species invasive in North America, V. hirundinaria, V. nigrum and V. rossicum, are members of the Vincetoxicum s. str. Clade. We explore the prerequisites for the range expansion in the Vincetoxicum s. str. Clade performing Bayesian and Maximum Likelihood phylogenetic analyses on sequences of the nuclear internal transcribed spacer (ITS) region, the nuclear external transcribed spacer region (ETS), and five plastid markers. The resulting phylogeny is used to conduct biogeographic analysis using BioGeoBEARS to reconstruct ancestral species ranges. Moreover, we map the known occurrences of two rare characters in Asclepiadoideae, the possession of phenanthroindolizidine alkaloids and reported cases of autogamy onto our phylogeny. Finally, we have conducted ecological niche modelling using Maxent on a total of 220 spatially unique occurrences of nine Vincetoxicum s. str. species spanning more than 4,000 km along the east-west gradient to learn about the climatic conditions along the presumed migration route. Our results indicate a north-westward migration in Vincetoxicum s. str. along the Asian mountain chains to Europe. Climatic preferences of the nine species sampled are dissimilar, except for the common exposure to at least one month of subfreezing temperatures, indicating a rather wide climatic tolerance for the clade as a whole. The three species invasive in North America belong to the northern Eurasian subclade and show the rare combination of phenanthroindolizidine alkaloids and autogamy.

Refining climate change projections for organisms with low dispersal abilities: a case study of the Caspian whip snake.

Climate warming is one of the most important threats to biodiversity. Ectothermic organisms such as amphibians and reptiles are especially vulnerable as climatic conditions affect them directly. Ecological niche models (ENMs) are increasingly popular in ecological studies, but several drawbacks exist, including the limited ability to account for the dispersal potential of the species. In this study, we use ENMs to explore the impact of global climate change on the Caspian whip snake (Dolichophis caspius) as model for organisms with low dispersal abilities and to quantify dispersal to novel areas using GIS techniques. Models generated using Maxent 3.3.3 k and GARP for current distribution were projected on future climatic scenarios. A cost-distance analysis was run in ArcGIS 10 using geomorphological features, ecological conditions, and human footprint as "costs" to dispersal of the species to obtain a Maximum Dispersal Range (MDR) estimate. All models developed were statistically significant (P<0.05) and recovered the currently known distribution of D. caspius. Models projected on future climatic conditions using Maxent predicted a doubling of suitable climatic area, while GARP predicted a more conservative expansion. Both models agreed on an expansion of suitable area northwards, with minor decreases at the southern distribution limit. The MDR area calculated using the Maxent model represented a third of the total area of the projected model. The MDR based on GARP models recovered only about 20% of the total area of the projected model. Thus, incorporating measures of species' dispersal abilities greatly reduced estimated area of potential future distributions.