Vocalizations of anuran amphibians have received much attention in studies of behavioral ecology and physiology, but also provide informative characters for identifying and delimiting species. We here review the terminology and variation of frog calls from a perspective of integrative taxonomy, and provide hands-on protocols for recording, analyzing, comparing, interpreting and describing these sounds. Our focus is on advertisement calls, which serve as premating isolation mechanisms and, therefore, convey important taxonomic information. We provide recommendations for terminology of frog vocalizations, with call, note and pulse being the fundamental subunits to be used in descriptions and comparisons. However, due to the complexity and diversity of these signals, an unequivocal application of the terms call and note can be challenging. We therefore provide two coherent concepts that either follow a note-centered approach (defining uninterrupted units of sound as notes, and their entirety as call) or a call-centered approach (defining uninterrupted units as call whenever they are separated by long silent intervals) in terminology. Based on surveys of literature, we show that numerous call traits can be highly variable within and between individuals of one species. Despite idiosyncrasies of species and higher taxa, the duration of calls or notes, pulse rate within notes, and number of pulses per note appear to be more static within individuals and somewhat less affected by temperature. Therefore, these variables might often be preferable as taxonomic characters over call rate or note rate, which are heavily influenced by various factors. Dominant frequency is also comparatively static and only weakly affected by temperature, but depends strongly on body size. As with other taxonomic characters, strong call divergence is typically indicative of species-level differences, whereas call similarities of two populations are no evidence for them being conspecific. Taxonomic conclusions can especially be drawn when the general advertisement call structure of two candidate species is radically different and qualitative call differences are thus observed. On the other hand, quantitative differences in call traits might substantially vary within and among conspecific populations, and require careful evaluation and analysis. We provide guidelines for the taxonomic interpretation of advertisement call differences in sympatric and allopatric situations, and emphasize the need for an integrative use of multiple datasets (bio-acoustics, morphology, genetics), particularly for allopatric scenarios. We show that small-sized frogs often emit calls with frequency components in the ultrasound spectrum, although it is unlikely that these high frequencies are of biological relevance for the majority of them, and we illustrate that detection of upper harmonics depends also on recording distance because higher frequencies are attenuated more strongly. Bioacoustics remains a prime approach in integrative taxonomy of anurans if uncertainty due to possible intraspecific variation and technical artifacts is adequately considered and acknowledged.

Abstract. A precise knowledge of the spatial distribution of taxa is essential for decision-making processes in land management and biodiversity conservation, both for present and under future global change scenarios. This is a key base for several scientific disciplines (e.g. macro-ecology, biogeography, evolutionary biology, spatial planning, or environmental impact assessment) that rely on species distribution maps. An atlas summarizing the distribution of European amphibians and reptiles with 50 × 50 km resolution maps based on ca. 85 000 grid records was published by the Societas Europaea Herpetologica (SEH) in 1997. Since then, more detailed species distribution maps covering large parts of Europe became available, while taxonomic progress has led to a plethora of taxonomic changes including new species descriptions. To account for these progresses, we compiled information from different data sources: published in books and websites, ongoing national atlases, personal data kindly provided to the SEH, the 1997 European Atlas, and the Global Biodiversity Information Facility (GBIF). Databases were homogenised, deleting all information except species names and coordinates, projected to the same coordinate system (WGS84) and transformed into a 50 × 50 km grid. The newly compiled database comprises more than 384 000 grid and locality records distributed across 40 countries. We calculated species richness maps as well as maps of Corrected Weighted Endemism and defined species distribution types (i.e. groups of species with similar distribution patterns) by hierarchical cluster analysis using Jaccard's index as association measure. Our analysis serves as a preliminary step towards an interactive, dynamic and online distributed database system (NA2RE system) of the current spatial distribution of European amphibians and reptiles. The NA2RE system will serve as well to monitor potential temporal changes in their distributions. Grid maps of all species are made available along with this paper as a tool for decision-making and conservation-related studies and actions. We also identify taxonomic and geographic gaps of knowledge that need to be filled, and we highlight the need to add temporal and altitudinal data for all records, to allow tracking potential species distribution changes as well as detailed modelling of the impacts of land use and climate change on European amphibians and reptiles.

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.

Species distributed across vast continental areas and across major biomes provide unique model systems for studies of biotic diversification, yet also constitute daunting financial, logistic and political challenges for data collection across such regions. The tree frog Dendropsophus minutus (Anura: Hylidae) is a nominal species, continentally distributed in South America, that may represent a complex of multiple species, each with a more limited distribution. To understand the spatial pattern of molecular diversity throughout the range of this species complex, we obtained DNA sequence data from two mitochondrial genes, cytochrome oxidase I (COI) and the 16S rhibosomal gene (16S) for 407 samples of D. minutus and closely related species distributed across eleven countries, effectively comprising the entire range of the group. We performed phylogenetic and spatially explicit phylogeographic analyses to assess the genetic structure of lineages and infer ancestral areas. We found 43 statistically supported, deep mitochondrial lineages, several of which may represent currently unrecognized distinct species. One major clade, containing 25 divergent lineages, includes samples from the type locality of D. minutus. We defined that clade as the D. minutus complex. The remaining lineages together with the D. minutus complex constitute the D. minutus species group. Historical analyses support an Amazonian origin for the D. minutus species group with a subsequent dispersal to eastern Brazil where the D. minutus complex originated. According to our dataset, a total of eight mtDNA lineages have ranges >100,000 km2. One of them occupies an area of almost one million km2 encompassing multiple biomes. Our results, at a spatial scale and resolution unprecedented for a Neotropical vertebrate, confirm that widespread amphibian species occur in lowland South America, yet at the same time a large proportion of cryptic diversity still remains to be discovered.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analyzed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions, compared to environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes, and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host's macro-environment mediate microbial diversity.

The dominant hypothesis for the evolutionary origin of snakes from 'lizards' (non-snake squamates) is that stem snakes acquired many snake features while passing through a profound burrowing (fossorial) phase. To investigate this, we examined the visual pigments and their encoding opsin genes in a range of squamate reptiles, focusing on fossorial lizards and snakes. We sequenced opsin transcripts isolated from retinal cDNA and used microspectrophotometry to measure directly the spectral absorbance of the photoreceptor visual pigments in a subset of samples. In snakes, but not lizards, dedicated fossoriality (as in Scolecophidia and the alethinophidian Anilius scytale) corresponds with loss of all visual opsins other than RH1 (k max 490-497 nm); all other snakes (including less dedicated burrowers) also have functional sws1 and lws opsin genes. In contrast, the retinas of all lizards sampled, even highly fossorial amphisbaenians with reduced eyes, express functional lws, sws1, sws2 and rh1 genes, and most also express rh2 (i.e. they express all five of the visual opsin genes present in the ancestral vertebrate). Our evidence of visual pigment complements suggests that the visual system of stem snakes was partly reduced, with two (RH2 and SWS2) of the ancestral vertebrate visual pigments being eliminated, but that this did not extend to the extreme additional loss of SWS1 and LWS that subsequently occurred (probably independently) in highly fossorial extant scolecophidians and A. scytale. We therefore consider it unlikely that the ancestral snake was as fossorial as extant scolecophidians, whether or not the latter are para-or monophyletic.

No abstract

Many tropical organisms show large genetic differences among populations, yet the prevalent drivers of the underlying divergence processes are incompletely understood. We explored the effect of several habitat and natural history features (body size, macrohabitat, microhabitat, reproduction site, climatic heterogeneity, and topography) on population genetic divergence in tropical amphibians, based on a data set of 2680 DNA sequences of the mitochondrial cytochrome b gene in 39 widely distributed frog species from Brazil, Central America, Cuba, and Madagascar. Generalized linear models were implemented in an information-theoretic framework to evaluate the effects of the six predictors on genetic divergence among populations, measured as spatially corrected pairwise distances. Results indicate that topographic complexity and macrohabitat preferences have a strong effect on population divergence with species specialized to forest habitat and/or from topographically complex regions showing higher phylogeographic structure. This relationship changed after accounting for phylogenetic relatedness among Electronic supplementary material The online version of this article (taxa rendering macrohabitat preferences as the most important feature shaping genetic divergence. The remaining predictors showed negligible effects on the observed genetic divergence. A similar analysis performed using the population-scaled mutation rate (H) as response variable showed little effect of the predictors. Our results demonstrate greater evolutionary independence among populations of anurans from forested regions versus species from open habitats. This pattern may result from lower vagility and stringency in reproductive requirements of rainforest species. Conversely, open landscapes may offer ephemeral and unstable breeding sites suitable for vagile generalist species, resulting in reduced intraspecific divergence. Our results predict that, for a given period of time, there should be a higher chance of speciation in tropical anurans living in forests than in species adapted to open habitats.