Connectivity of local amphibian populations: modelling the migratory capacity of radio-tracked natterjack toads

Landscape resistance index aiming at functional forest connectivity

Resistance models may quantify the ability of the landscape to impede species' movement and represent suitable habitats. Moreover, the performance of resistance models parameterized by land-use/land cover attributes evidence that the integrity of the environments subject to urban sprawl is poorly understood. In this sense, the study assumed we could identify the forest functional connectivity in a landscape considering the disparity in the landscape mosaic. In this context, we sought to develop a landscape resistance index through structural equation modeling (SEM), supported by the criteria of heat emission, biomass, and anthropogenic barriers, obtained by remote sensing, called observed variables. The landscape studied in the Green Belt Biosphere Reserve of São Paulo has significant remnants of the Atlantic Forest, a biodiversity hotspot. However, our results indicated criteria variability in the landscape modeled through the SEM, obtaining a significant adjustment of the landscape resistance index, with comparative fit index (CFI) of 1.00 and root mean square error of approximation (RMSEA) of 0.00. The index reflects the resistance levels of the land use/land cover, expressed by the class interval, ranging from 0% (1.73) to 100% (493.88), with the highest values associated with the anthropized uses and forest isolation. Thus, our index based on environmental attributes reflects the structure of functional forest connectivity and offers a framework to design forest corridors across landscapes.

Integrating dispersal, breeding and abundance data with graph theory for the characterization and management of functional connectivity in amphibian pondscapes

CONTEXT

Robust assessment of functional connectivity in amphibian population networks is essential to address their global decline. The potential of graph theory to characterize connectivity among amphibian populations has already been confirmed, but the movement data on which modelled graphs rely are often scarce and inaccurate. While probabilistic methods that account for intraspecific variability in dispersal better reflect the biological reality of functional connectivity, they must be informed by systematically recorded individual movement data, which are difficult to obtain for secretive taxa like amphibians.

OBJECTIVES

Our aim is to assess the applied potential of probabilistic graph theory to characterize overall connectivity across amphibian pondscapes using fine-scale capture-recapture data, and to inform conservation management based on the role of ponds on functional connectivity.

METHODS

We monitored an amphibian community in a pondscape located in a Spanish "dehesa" for 2 years. Photoidentification was used to build capture histories for individuals of six species, from which dispersal kernels and population sizes were estimated to model probabilistic graphs.

RESULTS

We obtained kernels of variable robustness for six species. Node importance for connectivity varied between species, but with common patterns such as shared road crossing areas and the presence of coincident interconnected pond clusters.

CONCLUSIONS

The combination of photoidentification, capture-recapture data and graph theory allowed us to characterize functional connectivity across the pondscape of study accounting for dispersal variability and identify areas where conservation actions could be most efficient. Our results highlight the need to account for interspecific differences in the study and management of amphibian pondscapes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10980-022-01520-x.

Dispersal without drivers: Intrinsic and extrinsic variables have no impact on movement distances in a terrestrial amphibian

Dispersive movements are often thought to be multicausal and driven by individual body size, sex, conspecific density, environmental variation, personality, and/or other variables. Yet such variables often do not account for most of the variation among dispersive movements in nature, leaving open the possibility that dispersion may be indeterministic. We assessed the amount of variation in 24 h movement distances that could be accounted for by potential drivers of displacement with a large empirical dataset of movement distances performed by Fowler's Toads (Anaxyrus fowleri) on the northern shore of Lake Erie at Long Point, Ontario (2002–2021, incl.). These toads are easy to sample repeatedly, can be identified individually and move parallel to the shoreline as they forage at night, potentially dispersing to new refuge sites. Using a linear mixed‐effect model that incorporated random effect terms to account for sampling variance and inter‐annual variation, we found that all potential intrinsic and extrinsic drivers of movement accounted for virtually none of the variation observed among 24 h distances moved by these animals, whether over short or large spatial scales. We examined the idea of movement personality by testing variance per individual toad and found no evidence of individuality in movement distances. We conclude that deterministic variables, whether intrinsic or extrinsic, neither can be shown to nor are necessary to drive movements in this population over all spatial scales. Stochastic, short time‐scale movements, such as daily foraging movements, can instead accumulate over time to produce large spatial‐scale movements that are dispersive in nature.

Did decades of glyphosate use have selected for resistant amphibians in agricultural habitats?

Glyphosate-based herbicides are used worldwide, and glyphosate's primary metabolite (aminomethylphosphonic acid: AMPA), is globally retrieved in surface waters. AMPA induces various adverse effects on aquatic wildlife, including selective mortality, which suggests that glyphosate exposure may have selected for AMPA-resistant individuals. We tested this hypothesis using spined toads (Bufo spinosus), an amphibian found in a variety of habitats, from AMPA-exposed agricultural lands to AMPA-free forested areas. We predicted that the offspring of individuals originating from agricultural habitats would develop AMPA-resistance - and be less prone to develop adverse effects from- AMPA exposure. To investigate this question, we performed a common garden brood-rearing experiment. The embryos and larvae of 40 spined toad pairs captured in agricultural and forest ponds were exposed either to an environmental relevant concentration of AMPA (0.4 μg L^-1) or to control conditions (n = 8160 embryos, n = 240 tadpoles). We monitored development durations, developmental abnormalities and morphology, measured across key developmental stages. Although we observed significant effects of AMPA on fitness parameters in each group, these effects were not exacerbated in individuals from AMPA-free habitats. We suggest that temporal and/or spatial dynamics of contamination, as well as gene flow between exposed and preserved populations, may hinder adaptive divergence between populations. Yet, we show strong adverse effects of AMPA exposure at early developmental stages. AMPA could therefore be one of the numerous causes of declining wild amphibian populations.

Breeding behaviour of ectotherms at high latitudes: the case of the natterjack toad Epidalea calamita at its northern range limit

We explored the breeding behaviour of a threatened amphibian, the natterjack toad, at its northern range limit in Estonia, to determine the extent to which reproduction is affected by harsh and unstable climatic conditions. Using photo identification of specimens, we found that in optimal weather conditions males formed three breeding cohorts, while in adverse conditions only a single cohort occurred and under extreme conditions reproduction was skipped entirely. During the extended breeding season, larger males participated in reproduction throughout the breeding period, while smaller males appeared in later cohorts. Breeding success was related to the calling effort of a male, where larger males had greater mating success than smaller ones. We found that the natterjack toad males exhibit significant plasticity in reproductive behaviour at the northern range limit, which, given the energetic cost of reproduction and the increased risk of predation, allows them to increase their fitness at high latitudes.

Modeling habitat connectivity in support of multiobjective species movement: An application to amphibian habitat systems

Reasoning about the factors underlying habitat connectivity and the inter-habitat movement of species is essential to many areas of biological inquiry. In order to better describe and understand the ways in which the landscape may support species movement, an increasing amount of research has focused on identification of paths or corridors that may be important in providing connectivity among habitat. The least-cost path problem has proven to be an instrumental analytical tool in this sense. A complicating aspect of such path identification methods is how to best reconcile and integrate the array of criteria or objectives that species may consider in traversal of a landscape. In cases where habitat connectivity is thought to be influenced or guided by multiple objectives, numerous solutions to least-cost path problems can exist, representing tradeoffs between the objectives. In practice though, identification of these solutions can be very challenging and as such, only a small proportion of them are typically examined leading to a weak characterization of habitat connectivity. To address this computational challenge, a multiobjective optimization framework is proposed. A generalizable multiobjective least-cost path model is first detailed. A non-inferior set estimation (MONISE) algorithm for identifying supported efficient solutions to the multiobjective least-cost path model is then described. However, it is well known that unsupported efficient solutions (which are equally important) can also exist, but are typically ignored given that they are more difficult to identify. Thus, to enable the identification of the full set of efficient solutions (supported and unsupported) to the multiobjective model, a multi-criteria labeling algorithm is then proposed. The developed framework is applied to assess different conceptualizations of habitat connectivity supporting amphibian movement in a wetland system. The results highlight the range of tradeoffs in characterizations of connectivity that can exist when multiple objectives are thought to contribute to movement decisions and that the number of unsupported efficient solutions (which are typically ignored) can vastly outweigh that of the supported efficient solutions.

A Multispecies Assessment to Identify the Functional Connectivity of Amphibians in a Human-Dominated Landscape

Preserving functional connectivity is a key goal of conservation management. However, the spatially confined conservation areas may not allow for dispersal and gene flow for the intended long-term persistence of populations in fragmented landscapes. We provide a regional multi-species assessment to quantify functional connectivity for five amphibian species in a human dominated landscape in the Swiss lowlands. A set of resistance maps were derived based on expert opinion and a sensitivity analysis was conducted to compare the effect of each resistance scenario on modelled connectivity. Deriving multi-species corridors is a robust way to identify movement hotspots that provide valuable baseline information to reinforce protective measures and green infrastructure.

Reproductive behavior drives female space use in a sedentary Neotropical frog

Longer-range movements of anuran amphibians such as mass migrations and habitat invasion have received a lot of attention, but fine-scale spatial behavior remains largely understudied. This gap is especially striking for species that show long-term site fidelity and display their whole behavioral repertoire in a small area. Studying fine-scale movement with conventional capture-mark-recapture techniques is difficult in inconspicuous amphibians: individuals are hard to find, repeated captures might affect their behavior and the number of data points is too low to allow a detailed interpretation of individual space use and time budgeting. In this study, we overcame these limitations by equipping females of the Brilliant-Thighed Poison Frog (Allobates femoralis) with a tag allowing frequent monitoring of their location and behavior. Neotropical poison frogs are well known for their complex behavior and diverse reproductive and parental care strategies. Although the ecology and behavior of the polygamous leaf-litter frog Allobates femoralis is well studied, little is known about the fine-scale space use of the non-territorial females who do not engage in acoustic and visual displays. We tracked 17 females for 6 to 17 days using a harmonic direction finder to provide the first precise analysis of female space use in this species. Females moved on average 1 m per hour and the fastest movement, over 20 m per hour, was related to a subsequent mating event. Traveled distances and activity patterns on days of courtship and mating differed considerably from days without reproduction. Frogs moved more on days with lower temperature and more precipitation, but mating seemed to be the main trigger for female movement. We observed 21 courtships of 12 tagged females. For seven females, we observed two consecutive mating events. Estimated home ranges after 14 days varied considerably between individuals and courtship and mating associated space use made up for ∼30% of the home range. Allobates femoralis females spent large parts of their time in one to three small centers of use. Females did not adjust their time or space use to the density of males in their surroundings and did not show wide-ranging exploratory behavior. Our study demonstrates how tracking combined with detailed behavioral observations can reveal the patterns and drivers of fine-scale spatial behavior in sedentary species.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

When new human-modified habitats favour the expansion of an amphibian pioneer species: Evolutionary history of the natterjack toad (Bufo calamita) in a coal basin

Human activities affect microevolutionary dynamics by inducing environmental changes. In particular, land cover conversion and loss of native habitats decrease genetic diversity and jeopardize the adaptive ability of populations. Nonetheless, new anthropogenic habitats can also promote the successful establishment of emblematic pioneer species. We investigated this issue by examining the population genetic features and evolutionary history of the natterjack toad (Bufo [Epidalea] calamita) in northern France, where populations can be found in native coastal habitats and coalfield habitats shaped by European industrial history, along with an additional set of European populations located outside this focal area. We predicted contrasting patterns of genetic structure, with newly settled coalfield populations departing from migration-drift equilibrium. As expected, coalfield populations showed a mosaic of genetically divergent populations with short-range patterns of gene flow, and native coastal populations indicated an equilibrium state with an isolation-by-distance pattern suggestive of postglacial range expansion. However, coalfield populations exhibited (i) high levels of genetic diversity, (ii) no evidence of local inbreeding or reduced effective population size and (iii) multiple maternal mitochondrial lineages, a genetic footprint depicting independent colonization events. Furthermore, approximate Bayesian computations suggested several evolutionary trajectories from ancient isolation in glacial refugia during the Pleistocene, with biogeographical signatures of recent expansion probably confounded by human-mediated mixing of different lineages. From an evolutionary and conservation perspective, this study highlights the ecological value of industrial areas, provided that ongoing regional gene flow is ensured within the existing lineage boundaries.

Movement ecology of amphibians: from individual migratory behaviour to spatially structured populations in heterogeneous landscapes,

Both genetic cohesion among local populations of animals and range expansion depend on the frequency of dispersers moving at an interpatch scale. Animal movement has an individual component that reflects behaviour and an ecological component that reflects the spatial organization of populations. The total movement capacity of an individual describes maximum movement distance theoretically achievable during a lifetime, whereas its variation among the members of a local population determines the magnitude of interpatch movements and thus of gene flow between neighbouring patches within metapopulation or patchy population systems. Here, I review information on dispersal and migration as components of the movement capacity of juvenile and adult pond-breeding amphibians and discuss how these components inform the spatial structure of populations. Amphibians disperse as juveniles and adults, but movement distances detected in tracking or capture–mark–recapture studies are usually far below the corresponding estimates based on molecular gene-flow data. This discrepancy reflects the constraints of available tracking methods for free-ranging individuals leading to inappropriate surrogates of annual movement capacity, but can be resolved using probabilistic approaches based on dispersal functions. There is remarkable capacity for and plasticity in movements in amphibians. Annual within-patch movements (migrations) of individuals can be large and likely represent an underestimated capacity for movement at the interpatch scale. Landscape resistance may influence the paths of dispersing amphibians, but rarely impedes interpatch movements. Juveniles emigrating unpredictably far from the natal pond and adults switching from within-patch migrations to dispersal to another patch demonstrate the plasticity of individual movement behaviour. Three basic conclusions can be drawn with respect to the linkage of individual movement behaviour and spatial or genetic structure of local amphibian populations embedded in a heterogeneous landscape: (1) individual movements or consecutive short-term series of movements are misleading surrogate measures of total movement capacity; (2) probabilistic modelling of movement capacity is the best available behavioural predictor of interpatch gene flow; (3) connectivity of local populations in heterogeneous landscapes is less affected by landscape resistance than previously expected.

Genetic variability in geographic populations of the natterjack toad (Bufo calamita)

Across altitudinal and latitudinal gradients, the proportion of suitable habitats varies, influencing the individual dispersal that ultimately can produce differentiation among populations. The natterjack toad (Bufo calamita) is distributed across a wide geographic range that qualifies the species as interesting for a geographic analysis of its genetic variability. Five populations of B. calamita in the Sierra de Gredos (Spain) were studied in an altitudinal gradient ranging from 750 to 2270 m using microsatellite markers. In addition, we analyzed the latitudinal genetic variation in B. calamita within a global European distribution using genetic diversity parameters (mean number of alleles per locus [M_a] and expected heterozygosity [H_E]) obtained from our results and those published in the literature. The low level of genetic differentiation found between populations of B. calamita (F_st ranging from 0.0115 to 0.1018) and the decreases in genetic diversity with altitude (M_a from 13.6 to 8.3, H_E from 0.82 to 0.74) can be interpreted by the combined effects of discontinuous habitat, produced mainly by the high slopes barriers and geographic distance. In the latitudinal gradient, genetic diversity decreases from south to north as a consequence of the colonization of the species from the Pleistocene refugium. We conclude that the genetic variability in B. calamita along its wide altitudinal and latitudinal geographic distribution mainly reflects the colonization history of the species after the last glacial period.