Influences of ski-runs, meadow management and climate on the occupancy of reptiles and amphibians in a high-altitude environment of Italy.

Alpine ecosystems harbour a rich and highly specialised biodiversity, which is particularly susceptible to anthropogenic disturbances such as habitat loss and fragmentation as well as to climate change. Combined with other forms of land-use conversion, construction and maintenance of ski resorts can have severe consequences on alpine biodiversity. In this study, we show how one amphibian and two reptile species, namely Rana temporaria, Zootoca vivipara and Vipera berus, respond to such impacts by means of a multi-season occupancy analysis. We found all three species both in and outside ski-runs, showing that these habitats do not necessarily preclude their occurrence. Contrarily, this is influenced more by microhabitat availability, such as ground vegetation, humid areas and rock cover, rather than by macro-characteristics like elevation or habitat type. Moreover, we found a climatic influence on the year-to-year occupancy change of the species, with activity-month conditions being more relevant than overwintering ones. Our results demonstrate how, in the specific case of reptiles and amphibians, ski resorts do not necessarily limit species' occurrence and that a mild series of management actions might secure the species' persistence in the area.

Diverse aging rates in ectothermic tetrapods provide insights for the evolution of aging and longevity.

Comparative studies of mortality in the wild are necessary to understand the evolution of aging; yet, ectothermic tetrapods are underrepresented in this comparative landscape, despite their suitability for testing evolutionary hypotheses. We present a study of aging rates and longevity across wild tetrapod ectotherms, using data from 107 populations (77 species) of nonavian reptiles and amphibians. We test hypotheses of how thermoregulatory mode, environmental temperature, protective phenotypes, and pace of life history contribute to demographic aging. Controlling for phylogeny and body size, ectotherms display a higher diversity of aging rates compared with endotherms and include phylogenetically widespread evidence of negligible aging. Protective phenotypes and life-history strategies further explain macroevolutionary patterns of aging. Analyzing ectothermic tetrapods in a comparative context enhances our understanding of the evolution of aging.

Decision-making of citizen scientists when recording species observations.

Citizen scientists play an increasingly important role in biodiversity monitoring. Most of the data, however, are unstructured-collected by diverse methods that are not documented with the data. Insufficient understanding of the data collection processes presents a major barrier to the use of citizen science data in biodiversity research. We developed a questionnaire to ask citizen scientists about their decision-making before, during and after collecting and reporting species observations, using Germany as a case study. We quantified the greatest sources of variability among respondents and assessed whether motivations and experience related to any aspect of data collection. Our questionnaire was answered by almost 900 people, with varying taxonomic foci and expertise. Respondents were most often motivated by improving species knowledge and supporting conservation, but there were no linkages between motivations and data collection methods. By contrast, variables related to experience and knowledge, such as membership of a natural history society, were linked with a greater propensity to conduct planned searches, during which typically all species were reported. Our findings have implications for how citizen science data are analysed in statistical models; highlight the importance of natural history societies and provide pointers to where citizen science projects might be further developed.

Response of Common and Rare Beetle Species to Tree Species and Vertical Stratification in a Floodplain Forest.

Vertical stratification and host tree species are factors with a high influence on the structure of communities of xylobiont beetles. However, little is known about how this influence varies between common and rare species. Based on estimated species richness, we compared alpha and beta diversity patterns of common and rare species in the canopy of the Leipzig floodplain forest to assess their response to vertical stratification and tree species. We used two measures of rarity: threat level in red lists and abundance based on octaves. The understory displayed a significantly higher number of common species than the canopy strata. Conversely, the canopy strata harbored a higher number of rare species. Turnover was always dominant over richness differences in beta diversity partitions. Using Raup-Crick null models and non-metric multidimensional scaling, we found that the vertical strata accounted for 19% of the overall beta diversity of common species and for 15% of the overall beta diversity of rare species. The tree species accounted for 7% of the overall beta diversity of the common species and 3% of the beta diversity of the rare species. Our results indicate that studies carried out in the understory alone do not allow drawing conclusions regarding the biodiversity in the canopy strata, and thus regarding the overall community structure of xylobiont beetles in the canopy.

Sampling and modelling rare species: Conceptual guidelines for the neglected majority.

Biodiversity conservation faces a methodological conundrum: Biodiversity measurement often relies on species, most of which are rare at various scales, especially prone to extinction under global change, but also the most challenging to sample and model. Predicting the distribution change of rare species using conventional species distribution models is challenging because rare species are hardly captured by most survey systems. When enough data are available, predictions are usually spatially biased towards locations where the species is most likely to occur, violating the assumptions of many modelling frameworks. Workflows to predict and eventually map rare species distributions imply important trade-offs between data quantity, quality, representativeness and model complexity that need to be considered prior to survey and analysis. Our opinion is that study designs need to carefully integrate the different steps, from species sampling to modelling, in accordance with the different types of rarity and available data in order to improve our capacity for sound assessment and prediction of rare species distribution. In this article, we summarize and comment on how different categories of species rarity lead to different types of occurrence and distribution data depending on choices made during the survey process, namely the spatial distribution of samples (where to sample) and the sampling protocol in each selected location (how to sample). We then clarify which species distribution models are suitable depending on the different types of distribution data (how to model). Among others, for most rarity forms, we highlight the insights from systematic species-targeted sampling coupled with hierarchical models that allow correcting for overdispersion and spatial and sampling sources of bias. Our article provides scientists and practitioners with a much-needed guide through the ever-increasing diversity of methodological developments to improve the prediction of rare species distribution depending on rarity type and available data.

Accelerated Evolution of Tissue-Specific Genes Mediates Divergence Amidst Gene Flow in European Green Lizards.

The European green lizards of the Lacerta viridis complex consist of two closely related species, L. viridis and Lacerta bilineata that split less than 7 million years ago in the presence of gene flow. Recently, a third lineage, referred to as the "Adriatic" was described within the L. viridis complex distributed from Slovenia to Greece. However, whether gene flow between the Adriatic lineage and L. viridis or L. bilineata has occurred and the evolutionary processes involved in their diversification are currently unknown. We hypothesized that divergence occurred in the presence of gene flow between multiple lineages and involved tissue-specific gene evolution. In this study, we sequenced the whole genome of an individual of the Adriatic lineage and tested for the presence of gene flow amongst L. viridis, L. bilineata, and Adriatic. Additionally, we sequenced transcriptomes from multiple tissues to understand tissue-specific effects. The species tree supports that the Adriatic lineage is a sister taxon to L. bilineata. We detected gene flow between the Adriatic lineage and L. viridis suggesting that the evolutionary history of the L. viridis complex is likely shaped by gene flow. Interestingly, we observed topological differences between the autosomal and Z-chromosome phylogenies with a few fast evolving genes on the Z-chromosome. Genes highly expressed in the ovaries and strongly co-expressed in the brain experienced accelerated evolution presumably contributing to establishing reproductive isolation in the L. viridis complex.

Patterns of richness across forest beetle communities-A methodological comparison of observed and estimated species numbers.

Species richness is a frequently used measure of biodiversity. The compilation of a complete species list is an often unattainable goal. Estimators of species richness have been developed to overcome this problem. While the use of these estimators is becoming increasingly popular, working with the observed number of species is still common practice.To assess whether patterns of beetle communities based on observed numbers may be compared among each other, we compared patterns from observed and estimated numbers of species for beetle communities in the canopy of the Leipzig floodplain forest. These patterns were species richness and the number of shared species among three tree species and two canopy strata.We tested the applicability of the asymptotic Chao1 estimator and the estimate provided by the nonasymptotic rarefaction-extrapolation method for all tree species and both upper canopy and lower canopy. In the majority of cases, the ranking patterns of species richness for host tree species and strata were the same for the observed and estimated number of species. The ranking patterns of the number of species shared among host tree species and strata, however, were significantly different between observed and estimated values.Our results indicate that the observed number of species under-represents species richness and the number of shared species. However, ranking comparisons of published patterns based on the number of observed species may be acceptable for species richness but likely not reliable for the number of shared species. Further studies are needed to corroborate this conclusion. We encourage to use estimators and to provide open access to data to allow comparative assessments.

A GIS-based policy support tool to determine national responsibilities and priorities for biodiversity conservation.

Efficient biodiversity conservation requires that limited resources be allocated in accordance with national responsibilities and priorities. Without appropriate computational tools, the process of determining these national responsibilities and conservation priorities is time intensive when considering many species across geographic scales. Here, we have developed a computational tool as a module for the ArcGIS geographic information system. The ArcGIS National Responsibility Assessment Tool (NRA-Tool) can be used to create hierarchical lists of national responsibilities and priorities for global species conservation. Our tool will allow conservationists to prioritize conservation efforts and to focus limited resources on relevant species and regions. We showcase our tool with data on 258 bird species and various biophysical regions, including Environmental Zones in 58 Asian countries and regions. Our tool provides a decision support system for conservation policy with attractive and easily interpretable visual outputs illustrating national responsibilities and priorities for species conservation. The graphical output allows for smooth integration into assessment reports, such as the European Article 17 report, the Living Planet Index report, or similar regional and global reports.

Correction: Predicting strike susceptibility and collision patterns of the common buzzard at wind turbine structures in the federal state of Brandenburg, Germany.

[This corrects the article DOI: 10.1371/journal.pone.0227698.].

Differential effects of habitat loss on occupancy patterns of the eastern green lizard Lacerta viridis at the core and periphery of its distribution range.

The effects of habitat loss on the distribution of populations are often linked with species specialization degree. Specialist species can be more affected by changes in landscape structure and local patch characteristics compared to generalist species. Moreover, the spatial scale at which different land covers (eg. habitat, cropland, urban areas) affect specialist species can be smaller. Specialization is usually assumed as a constant trait along the distribution range of species. However, for several taxa, there is evidence of higher specialization degree in peripheral populations compared with populations in the core. Hence, peripheral populations should have a higher sensitivity to habitat loss, and strongest effects should be found at a smaller spatial scale. To test these expectations, we implemented a patch-landscape approach at different spatial scales, and compared effects of landscape structure and patch characteristics on occupancy probability among northern peripheral, more specialized populations (Czech Republic) and core populations (Bulgaria) of the eastern green lizard Lacerta viridis. We found that landscape structure and patch characteristics affect differently the occupancy probability of Lacerta viridis in each region. Strongest effects of habitat loss were found at a spatial scale of 150m around patches in the periphery, but at a scale of 500m in the core. In the periphery occupancy probability of populations was principally affected by landscape composition, and the effect of habitat quality was stronger compared to core populations. In the core, persistence of populations was mainly explained by characteristics of the spatial configuration of habitat patches. We discuss possible ecological mechanisms behind the relationship between sensitivity to habitat loss, populations' specialization degree and position in the distribution range, and suggest conservation measures for L. viridis.

Predicting strike susceptibility and collision patterns of the common buzzard at wind turbine structures in the federal state of Brandenburg, Germany.

With the increase in wind turbines, bird collisions have developed as a potential hazard. In the federal state of Brandenburg, Germany, despite the on-going mitigation efforts of increasing the distances of wind turbines from the breeding areas of the more severely affected populations of red kites (Milvus milvus), the additional detrimental influences on the buzzard populations (Buteo buteo) have added to the challenges for wind power expansion. Using data on the regional distribution of the buzzards, along with their carcass detections around the wind turbines (WTs), we aimed to better understand their collision distribution patterns in relation to their habitat use patterns to predict their exposure to collision risk using boosted regression trees (BRTs). Additionally, we integrated the developed collision potential map with the regional density map of buzzards to identify areas of increased strike susceptibility in turbine installations. Our study showed that the buzzard collisions were primarily concentrated at the turbines situated at sensitive distances from the edges of watercourses (>1000 metres), as well as those along the edges of grasslands (>750 metres), in the green open areas around/areas with minimal settlements (750 metres-1750 metres), and along the edges of bushlands (>1500 metres), together explaining 58% of the variance in their collision distribution. Conclusively, our study is applicable to conservation because it demonstrates the identification of potential collision areas along with the causes of the collisions, in addition to demonstrating the benefits of incorporating a species collision dataset as a proxy for species presence into species distribution models to make informed management decisions to eventually combat biodiversity loss.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation.

BACKGROUND: Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes. FINDINGS: Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated. CONCLUSION: The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Realized niche and microhabitat selection of the eastern green lizard (Lacerta viridis) at the core and periphery of its distribution range.

The available range of habitats and suitable abiotic conditions like temperature and radiation tends to be narrower toward the periphery of the distribution range of species. Peripheral populations of generalist species could then be more specialized and have a smaller and differentiated realized niche (habitat niche in our study) compared to populations at the core. Likewise, patterns of microhabitat selection can differ between periphery and core. In our study, we compared niche size and microhabitat selection among core (Bulgaria) and northern peripheral (Germany, Czech Republic) populations of Lacerta viridis and estimated niche differentiation among regions. We collected data on vegetation structure and abiotic parameters at the microhabitat scale in each region. In order to compare niche size among regions and estimate niche differentiation, we built multidimensional niche hypervolumes. We applied generalized linear mixed models and model averaging, accounting for spatial autocorrelation when necessary, to analyze microhabitat differences among regions and microhabitat selection in each region. Peripheral populations were more specialized, having a smaller niche than core ones, and their niche differed from that in the core (Sørensen overlap in all comparisons <0.3). Microhabitats at the periphery had lower radiation and soil compaction and less structured vegetation. Microhabitat selection at the core depended solely on abiotic parameters, while at the periphery it was defined by only vegetation structure (Czech Republic) or a combination of both, vegetation structure, and abiotic factors (Germany). Thus, peripheral populations seem to compensate for overall harsher climatic conditions by responding to different parameters of the microhabitat compared to core populations. We suggest specific conservation measures for L. virids in each studied region and point out the general implications of a higher specialization degree of peripheral populations in relation to climate change and habitat fragmentation.

Collision sensitive niche profile of the worst affected bird-groups at wind turbine structures in the Federal State of Brandenburg, Germany.

Biodiversity-related impacts at wind energy facilities have increasingly become a cause of conservation concern, central issue being the collision of birds. Utilizing spatial information of their carcass detections at wind turbines (WTs), we quantified the detections in relation to the metric distances of the respective turbines to different land-use types. We used ecological niche factor analysis (ENFA) to identify combinations of land-use distances with respect to the spatial allocation of WTs that led to higher proportions of collisions among the worst affected bird-groups: Buntings, Crows, Larks, Pigeons and Raptors. We also assessed their respective similarities to the collision phenomenon by checking for overlaps amongst their distance combinations. Crows and Larks showed the narrowest "collision sensitive niche"; a part of ecological niche under higher risk of collisions with turbines, followed by that of Buntings and Pigeons. Raptors had the broadest niche showing significant overlaps with the collision sensitive niches of the other groups. This can probably be attributed to their larger home range combined with their hunting affinities to open landscapes. Identification of collision sensitive niches could be a powerful tool for landscape planning; helping avoid regions with higher risks of collisions for turbine allocations and thus protecting sensitive bird populations.

Too hot to die? The effects of vegetation shading on past, present, and future activity budgets of two diurnal skinks from arid Australia.

Behavioral thermoregulation is an important mechanism allowing ectotherms to respond to thermal variations. Its efficiency might become imperative for securing activity budgets under future climate change. For diurnal lizards, thermal microhabitat variability appears to be of high importance, especially in hot deserts where vegetation is highly scattered and sensitive to climatic fluctuations. We investigated the effects of a shading gradient from vegetation on body temperatures and activity timing for two diurnal, terrestrial desert lizards, Ctenotus regius, and Morethia boulengeri, and analyzed their changes under past, present, and future climatic conditions. Both species' body temperatures and activity timing strongly depended on the shading gradient provided by vegetation heterogeneity. At high temperatures, shaded locations provided cooling temperatures and increased diurnal activity. Conversely, bushes also buffered cold temperature by saving heat. According to future climate change scenarios, cooler microhabitats might become beneficial to warm-adapted species, such as C. regius, by increasing the duration of daily activity. Contrarily, warmer microhabitats might become unsuitable for less warm-adapted species such as M. boulengeri for which midsummers might result in a complete restriction of activity irrespective of vegetation. However, total annual activity would still increase provided that individuals would be able to shift their seasonal timing towards spring and autumn. Overall, we highlight the critical importance of thermoregulatory behavior to buffer temperatures and its dependence on vegetation heterogeneity. Whereas studies often neglect ecological processes when anticipating species' responses to future climate change the strongest impact of a changing climate on terrestrial ectotherms in hot deserts is likely to be the loss of shaded microhabitats rather than the rise in temperature itself. We argue that conservation strategies aiming at addressing future climate changes should focus more on the cascading effects of vegetation rather than on shifts of species distributions predicted solely by climatic envelopes.

Genetic diversity and spatial structure of the Rufous-throated Antbird (Gymnopithys rufigula), an Amazonian obligate army-ant follower.

Amazonian understory antbirds are thought to be relatively sedentary and to have limited dispersal ability; they avoid crossing forest gaps, and even narrow roads through a forest may limit their territories. However, most evidence for sedentariness in antbirds comes from field observations and plot-based recapture of adult individuals, which do not provide evidence for lack of genetic dispersal, as this often occurs through juveniles. In this study, we used microsatellite markers and mitochondrial control-region sequences to investigate contemporary and infer historical patterns of genetic diversity and structure of the Rufous-throated Antbird (Gymnopithys rufigula) within and between two large reserves in central Amazonia. Analyses based on microsatellites suggested two genetically distinct populations and asymmetrical gene flow between them. Within a population, we found a lack of genetic spatial autocorrelation, suggesting that genotypes are randomly distributed and that G. rufigula may disperse longer distances than expected for antbirds. Analyses based on mitochondrial sequences did not recover two clear genetic clusters corresponding to the two reserves and indicated the whole population of the Rufous-throated Antbird in the region has been expanding over the last 50,000 years. Historical migration rates were low and symmetrical between the two reserves, but we found evidence for a recent unilateral increase in gene flow. Recent differentiation between individuals of the two reserves and a unilateral increase in gene flow suggest that recent urban expansion and habitat loss may be driving changes and threatening populations of Rufous-throated Antbird in central Amazonia. As ecological traits and behavioral characteristics affect patterns of gene flow, comparative studies of other species with different behavior and ecological requirements will be necessary to better understand patterns of genetic dispersal and effects of urban expansion on Amazonian understory antbirds.

Environmental characteristics drive variation in Amazonian understorey bird assemblages.

Tropical bird assemblages display patterns of high alpha and beta diversity and, as tropical birds exhibit strong habitat specificity, their spatial distributions are generally assumed to be driven primarily by environmental heterogeneity and interspecific interactions. However, spatial distributions of some Amazonian forest birds are also often restricted by large rivers and other large-scale topographic features, suggesting that dispersal limitation may also play a role in driving species' turnover. In this study, we evaluated the effects of environmental characteristics, topographic and spatial variables on variation in local assemblage structure and diversity of birds in an old-growth forest in central Amazonia. Birds were mist-netted in 72 plots distributed systematically across a 10,000 ha reserve in each of three years. Alpha diversity remained stable through time, but species composition changed. Spatial variation in bird-assemblage structure was significantly related to environmental and topographic variables but not strongly related to spatial variables. At a broad scale, we found bird assemblages to be significantly distinct between two watersheds that are divided by a central ridgeline. We did not detect an effect of the ridgeline per se in driving these patterns, indicating that most birds are able to fly across it, and that differences in assemblage structure between watersheds may be due to unmeasured environmental variables or unique combinations of measured variables. Our study indicates that complex geography and landscape features can act together with environmental variables to drive changes in the diversity and composition of tropical bird assemblages at local scales, but highlights that we still know very little about what makes different parts of tropical forest suitable for different species.

Setting temporal baselines for biodiversity: the limits of available monitoring data for capturing the full impact of anthropogenic pressures.

Temporal baselines are needed for biodiversity, in order for the change in biodiversity to be measured over time, the targets for biodiversity conservation to be defined and conservation progress to be evaluated. Limited biodiversity information is widely recognized as a major barrier for identifying temporal baselines, although a comprehensive quantitative assessment of this is lacking. Here, we report on the temporal baselines that could be drawn from biodiversity monitoring schemes in Europe and compare those with the rise of important anthropogenic pressures. Most biodiversity monitoring schemes were initiated late in the 20th century, well after anthropogenic pressures had already reached half of their current magnitude. Setting temporal baselines from biodiversity monitoring data would therefore underestimate the full range of impacts of major anthropogenic pressures. In addition, biases among taxa and organization levels provide a truncated picture of biodiversity over time. These limitations need to be explicitly acknowledged when designing management strategies and policies as they seriously constrain our ability to identify relevant conservation targets aimed at restoring or reversing biodiversity losses. We discuss the need for additional research efforts beyond standard biodiversity monitoring to reconstruct the impacts of major anthropogenic pressures and to identify meaningful temporal baselines for biodiversity.

The complete mitochondrial genome of Lacerta bilineata and comparison with its closely related congener L. Viridis.

We sequenced the mitochondrial genome of the Western green lizard (Lacerta bilineata) using Illumina technology and additional Sanger sequencing. The assembled 17 086 bp mitogenome had a GC content of 40.32% and consisted of 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region (CR), with a gene order identical to the chordate consensus. In addition, we re-sequenced the mitogenome of the closely related Eastern green lizard L. viridis using the same techniques as for L. bilineata. The mitogenomes of L. bilineata and L. viridis showed a sequence identity of 94.4% and 99.9%, respectively, relative to the previously published L. viridis mitogenome. The phylogenetic reconstruction based on 17 Lacertinae mitogenomes using Anolis carolinensis as the outgroup supported L. bilineata and its sister species L. viridis as distinct lineages.