Accurate measures of species abundance are essential to identify conservation strategies. N-mixture models are increasingly used to estimate abundance on the basis of species counts. In this study we tested whether abundance estimates obtained using N-mixture models provide consistent results with more traditional approaches requiring capture (capture-mark recapture and removal sampling). We focused on endemic, threatened species of amphibians and reptiles in Italy, for which accurate abundance data are needed for conservation assessments: the Lanza’s Alpine salamander Salamandra lanzai, the Ambrosi’s cave salamander Hydromantes ambrosii and the Aeolian wall lizard Podarcis raffonei. In visual counts, detection probability was variable among species, ranging between 0.14 (Alpine salamanders) and 0.60 (cave salamanders). For all the species, abundance estimates obtained using N-mixture models showed limited differences with the ones obtained through capture-mark-recapture or removal sampling. The match was particularly accurate for cave salamanders in sites with limited abundance and for lizards, nevertheless non-incorporating heterogeneity of detection probability increased bias. N-mixture models provide reliable abundance estimates that are comparable with the ones of more traditional approaches, and offer additional advantages such as a smaller sampling effort and no need of manipulating individuals, which in turn reduces the risk of harming animals and spreading diseases.

Many organisms live in networks of local populations connected by dispersing individuals, called spatially structured populations (SSPs), where the long-term persistence of the entire network is determined by the balance between two processes acting at the scale of local populations: extinction and colonization. When multiple threats act on an SSP, a comparison of the different factors determining local extinctions and colonizations is essential to plan sound conservation actions. Here we assessed the drivers of long-term population dynamics of multiple amphibian species at the regional scale. We used dynamic occupancy models within a Bayesian framework to identify the factors determining persistence and colonization of local populations. Since connectivity among patches is fundamental for SSPs dynamics, we considered two measures of connectivity acting on each focal patch: incidence of the focal species and incidence of invasive crayfish. We used meta-analysis to summarize the effect of different drivers at the community level. Persistence and colonization of local populations were jointly determined by factors acting at different scales. Persistence probability was positively related to the area and the permanence of wetlands, while it showed a negative relationship with the occurrence of fish. Colonization probability was highest in semipermanent wetlands and in sites with a high incidence of the focal species in nearby sites, while it showed a negative relationship with the incidence of invasive crayfish in the landscape. By analyzing long-term data on amphibian population dynamics, we found a strong effect of some classic features commonly used in SSP studies, such as patch area and focal species incidence. The presence of an invasive alien species at the landscape-scale emerged as one of the strongest drivers of colonization dynamics, suggesting that studies on SSPs should consider different connectivity measures more frequently, such as the incidence of predators, especially when dealing with biological invasions.

Invasive predators can strongly affect native populations. If alien predator pressure is strong enough, it can induce anti-predator responses, including phenotypic plasticity of exposed individuals and local adaptations of impacted populations. Furthermore, maternal investment is an additional pathway that could provide resources and improve performance in the presence of alien predators. We investigated the potential responses to an alien predator crayfish (Procambarus clarkii) in a threatened frog (Rana latastei) by combining field observations with laboratory measurements of embryo development rate, to assess the importance of parental investment, origin and exposure to the crayfish cues. We detected a strong variation in parental investment amongst frog populations, but this variation was not related to the invasion status of the site of origin, suggesting that mothers did not modulate parental investment in relation to the presence of alien predators. However, cues of the invasive crayfish elicited plastic responses in clutches and tadpoles development: embryos developed faster when exposed to the predator. Furthermore, embryos from invaded sites reached Gosner’s development stage 25 faster than those from non-invaded sites. This ontogenetic shift can be interpreted as a local adaptation to the alien predator and suggests that frogs are able to recognise the predatory risk. If these plastic responses and local adaptation are effective escape strategies against the invasive predator, they may improve the persistence of native frog populations.

Aim Understanding which factors determine the variation in population size across space and time is crucial to plan sound conservation interventions. Amphibians are often characterized by large demographic changes, therefore a better understanding of factors driving these changes could help mitigate their global crisis. We investigated drivers of abundance dynamics of two similar frog species to understand why they met different fates in the same study area. Location Northern Italy. Methods In seven different years between 2004 and 2020, we performed repeated counts of egg masses of two similar frog species, the agile frog (Rana dalmatina) and the Italian agile frog (Rana latastei), in 31 wetlands, and used Bayesian models to estimate the relationships between frog abundance and candidate drivers acting at the scale of (i) site: wetland surface, shading percentage and the presence of the invasive red swamp crayfish (Procambarus clarkii); (ii) landscape: forest cover around wetlands; (iii) climate: yearly precipitation and (iv) spatially structured population (SSP): clutch incidence of the focal species and crayfish incidence in the surrounding landscape. Results The two species showed sharp differences in population size and trends: R. dalmatina was abundant and showed a stable trend throughout the entire study period; R. latastei showed low abundance in the first years and then almost disappeared. The abundance of R. dalmatina was positively related to forest cover, shading, wetland area and precipitation, while negatively related to the occurrence of invasive crayfish at both local and SSP level. R. latastei abundance increased with wetland area and precipitation, while models were not able to detect relationships with other factors. Main conclusions The high sensitivity to drought and demographic stochasticity could have contributed to the quasi‐extinction of R. latastei in the study area, highlighting that similar species can meet different fates even under the same environmental conditions.

Amphibians are an exemplary case of the current biodiversity crisis, being among the vertebrates suffering the fastest decline. Population dynamics of amphibians can result from processes acting at different scales. Both the local characteristics of breeding wetlands and the features of the surrounding landscape can strongly affect the temporal dynamics of amphibian populations. European newts are particularly threatened by land‐use change and invasive alien species. While it is known that newts are declining across Europe, few studies have performed broad‐scale assessments of their decline, either because abundance dynamics are more complex to analyse than presence/absence data or because they require a high sampling effort and long‐term monitoring. In this study, we show that long‐term distribution data can be combined with demographic models to quantify the decline in abundance of newt species at the regional scale, and to assess the importance of multiple factors in determining abundance dynamics. We performed multiple surveys between 1996 and 2020 and used N‐mixture models in a Bayesian framework. We then calculated abundance changes between the first and the last sampling season, which were performed with an average timespan of 13 years across all wetlands. Both Italian crested newts and smooth newts showed large declines, with an average estimated abundance loss between the first and last sampling season of 57% and 63%, respectively. Local characteristics of the wetlands were the main determinants of abundance dynamics: the abundance of both species showed a positive relationship with the area and the permanence of the wetland and a negative relationship with the presence of fish. Additionally, the abundance of Italian crested newts was negatively related to the presence of invasive crayfish. No relationship was detected between abundance and terrestrial habitat availability or connectivity measures. Despite uncertainties in the absolute values of estimated abundance, the striking regional‐scale decline of newts is evident. Among the major determinants of population dynamics, fish and crayfish presence increased their prevalence in the study area, while other factors remained more stable. Management actions aimed at eradicating or controlling invasive fish and crayfish might halt abundance loss and even revert this declining trend. The application of N‐mixture models to long‐term data from representative sites permits the analysis of temporal trends of species at the regional scale even when data come from complex monitoring schemes. We found large declines in abundance of two newt species, suggesting that European newts may be more threatened than previously thought.

Most animal species are detected imperfectly and overlooking individuals can result in a biased inference of the abundance patterns and underlying processes. Several techniques can incorporate the imperfect detection process for a more accurate estimation of abundance, but most of them require repeated surveys, i.e., more sampling effort compared to single counts. In this study, we used the dependent double-observer approach to estimate the detection probability of the egg clutches of two brown frog species, Rana dalmatina and R. latastei. We then simulated the data of a declining population at different levels of detection probability in order to assess under which conditions the double counts provided better estimates of population trends compared to naïve egg counts, given the detectability of frog clutches. Both species showed a very high detection probability, with average values of 93% for Rana dalmatina and 97% for R. latastei. Simulations showed that not considering imperfect detection reduces the power of detecting population trends if detection probability is low. However, at high detection probability (>80%), ignoring the imperfect detection does not bias the estimates of population trends. This suggests that, for species laying large and easily identifiable egg clutches, a single count can provide useful estimates if surveys are correctly timed.

Hybridization is a widespread phenomenon in animals, and hybrid heterosis/breakdown could be key processes determining the evolutionary dynamics of hybrids. Indeed, hybrids are not consistently disadvantaged compared to the parental lineages, as was historically assumed. Multiple processes could lead to performance differences between parental lineages and their hybrids. Despite many studies evaluated the performance of hybrids, a quantitative synthesis is required to assess the general pattern. Here we used meta-analytic and meta-regression approaches to quantify the fitness differences between parental lineages and their hybrids, and to identify possible processes that could lead to these differences. Specifically, we tested biological and methodological parameters that could determine differences in performance between hybrids and parental lineages. Hybrid performance was extremely variable across studies, being often significantly higher or lower compared to the mean performance of their parents. Nevertheless, the averaged hybrid performance was similar to the fitness of parental lineages, with differences across studies related to how performance was assessed. Genetic divergence between parental lineages, and the approach used to identify hybrids were the parameters most strongly related to variation in hybrid performance. Performance was lower for hybrids between distantly related lineages. Furthermore, study settings and the use of imprecise approaches for hybrid identification (e.g. morphology-based) can bias assessments of performance. Studies performed on wild populations and using genetic approaches for hybrid identification detected more often a decreased hybrid performance, compared to laboratory studies. We highlight the importance of appropriate settings for a realistic understanding of the evolutionary impacts of hybridization.

The evolution of sexual dimorphism (SD) is driven by intricate interplays between sexual and natural selection. When it comes to SD variation within populations, however, environmental factors play a major role. Sexually selected traits are expected to be strongly dependent on individual body condition, which is influenced by the local environment that individuals experience. As a consequence, the degree of SD may also depend on resource availability. Here, we investigated the potential drivers of SD expression at two sexually dimorphic morphometric traits, body size (snout vent length) and head shape (head geometric morphometrics), in the Italian wall lizard (Podarcis siculus). We assessed the existence of condition- and context-dependent SD across ten islands of the Aeolian archipelago (southern Italy), at within- and among-population scales. We observed strong geographical variation of SD among islands, and tested three potential SD predictors related to resource availability (individual body condition, ecosystem productivity, temperature). Body condition and ecosystem productivity were the main drivers of body size SD variation, and body condition was also the main driver for head shape SD. Our results highlight that the expression of SD in the Italian wall lizard is both condition- and context-dependent. These results are congruent at within- and among-populations scales highlighting that spatial multi-scale analysis represents a useful approach to understand patterns of SD expression.