Amphibian phylogeography: a model for understanding historical aspects of species distributions

Apennine-Pyrenees disjunct distribution: an unusual biogeographic pattern revealed in flea beetles of the Longitarsus candidulus species-group (Coleoptera, Chrysomelidae)

European mountain systems have played a crucial role in shaping the distribution of species and of their genetic diversity during the Quaternary climatic changes, with the establishment of allopatric patterns across main mountain ranges. Here we investigated the evolutionary history of flea beetles of the Longitarsus candidulus species-group showing an uncommon disjunct biogeographic pattern across the Apennine and the Pyrenees. We applied a multilocus molecular approach and multispecies coalescent models to establish a phylogenetic and systematic framework for this morphologically homogeneous species-group and to estimate the time of main cladogenetic events underlying the origin of the Apennine-Pyrenees pattern. We found strong support for the monophyly of the candidulus group with a sister relationship between Longitarsus laureolae and L. leonardii endemic to the Apennine and the Pyrenees mountains respectively. The timing of speciation events in the candidulus species-group coincides with 2 major climatic transitions during the Early and Middle Pleistocene which resulted in significant environmental changes in Europe and suggest a scenario of allopatric isolation and divergence on distinct mountain ranges. The split between the thermophilic species L. candidulus and the ancestor of the temperate species L. laureolae and L. leonardii is estimated at ~3 Ma during the transition from Pliocene to Pleistocene and was probably triggered by their segregation in xerophilous and temperate habitats. The speciation between L. laureolae and L. leonardii, estimated at ~1 Ma during the Mid-Pleistocene Transition, can be explained by the establishment of unfavorable conditions in West Alps and Central Massif underlying the onset of the Apennine-Pyrenees disjunct pattern. Finally, the strict association between members of the candidulus group and distinct Thymelaeaceae plants suggests further studies to address the hypothesis that speciation in these flea beetles might have been also associated with Pleistocene range changes of their host plants.

Mitonuclear and phenotypic discordance in an Atlantic Forest frog hybrid zone

Discordance between mitochondrial and nuclear DNA is common among animals and can be the result of a number of evolutionary processes, including incomplete lineage sorting and introgression. Particularly relevant in contact zones, mitonuclear discordance is expected because the mitochondrial genome is haploid and primarily uniparentally inherited, whereas nuclear loci are evolving at slower rates. In addition, when closely related taxa come together in hybrid zones, the distribution of diagnostic phenotypic characters and their concordance with the mitochondrial or nuclear lineages can also inform on historical and ongoing dynamics within hybrid zones. Overall, genetic and phenotypic discordances provide evidence for evolutionary divergence and processes that maintain boundaries among sister species or lineages. In this study, we characterized patterns of genetic and phenotypic variation in a contact zone between Cycloramphus dubius and Cycloramphus boraceiensis, two sister species of frogs endemic to the Atlantic Coastal Forest of Brazil. We examined genomic-scale nuclear diversification across 19 populations, encompassing the two parental forms and a contact zone between them. We compared the distribution of genomic DNA variability with that of a mitochondrial locus (16S) and two morphological traits (dorsal tubercles and body size). Our results reveal multiple divergent lineages with ongoing admixture. We detected discordance in patterns of introgression across the three data types. Cycloramphus dubius males are significantly larger than C. boraceiensis males, and we posit that competition among males in the hybrid zone, coupled with mate choice by females, may be one mechanism leading to patterns of introgression observed between the species.

Differential Mitochondrial Genome Expression of Four Hylid Frog Species under Low-Temperature Stress and Its Relationship with Amphibian Temperature Adaptation

Extreme weather poses huge challenges for animals that must adapt to wide variations in environmental temperature and, in many cases, it can lead to the local extirpation of populations or even the extinction of an entire species. Previous studies have found that one element of amphibian adaptation to environmental stress involves changes in mitochondrial gene expression at low temperatures. However, to date, comparative studies of gene expression in organisms living at extreme temperatures have focused mainly on nuclear genes. This study sequenced the complete mitochondrial genomes of five Asian hylid frog species: Dryophytes japonicus, D. immaculata, Hyla annectans, H. chinensis and H. zhaopingensis. It compared the phylogenetic relationships within the Hylidae family and explored the association between mitochondrial gene expression and evolutionary adaptations to cold stress. The present results showed that in D. immaculata, transcript levels of 12 out of 13 mitochondria genes were significantly reduced under cold exposure (p < 0.05); hence, we put forward the conjecture that D. immaculata adapts by entering a hibernation state at low temperature. In H. annectans, the transcripts of 10 genes (ND1, ND2, ND3, ND4, ND4L, ND5, ND6, COX1, COX2 and ATP8) were significantly reduced in response to cold exposure, and five mitochondrial genes in H. chinensis (ND1, ND2, ND3, ND4L and ATP6) also showed significantly reduced expression and transcript levels under cold conditions. By contrast, transcript levels of ND2 and ATP6 in H. zhaopingensis were significantly increased at low temperatures, possibly related to the narrow distribution of this species primarily at low latitudes. Indeed, H. zhaopingensis has little ability to adapt to low temperature (4 °C), or maybe to enter into hibernation, and it shows metabolic disorder in the cold. The present study demonstrates that the regulatory trend of mitochondrial gene expression in amphibians is correlated with their ability to adapt to variable climates in extreme environments. These results can predict which species are more likely to undergo extirpation or extinction with climate change and, thereby, provide new ideas for the study of species extinction in highly variable winter climates.

Phylogeographic Substructuring in the Southernmost Refugium of the European Common Frog Rana temporaria

Rana temporaria is one of the most widespread Palearctic brown frogs. We aimed to clarify distribution pattern of two main genetic clades in the understudied Balkan peninsula by using 16SrRNA and MT-CYTB sequences, already widely applied in analyses of populations from other parts of Europe, while focusing on the broad area along the Morava river (central Balkans) as a known gap in the species distribution. Additionally, we were interested in revealing the extent of haplotype diversity within the main genetic clades in the Balkans, particularly around the supposed suture zone. The results revealed a suture zone between the Western and Eastern Clades in the central part of the Balkan Peninsula. This indicated the existence of a historical barrier between the Balkan Mountain Belt and geographically close mountains surrounding the Vlasina Plateau (Rhodope/Serbian-Macedonian Massif). The overall observed haplotype diversity in populations of R. temporaria from the Balkan Peninsula seems high. Harboring both main genetic clades of R. temporaria qualifies the Balkan Peninsula as another important center of species' genetic diversity, as well as rich in unique haplotypes. This points out the necessity of applying conservation measures focused on the common European frog populations and habitats in this part of the species' distribution area.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Characterization of Two Mitogenomes of Hyla sanchiangensis (Anura: Hylidae), with Phylogenetic Relationships and Selection Pressure Analyses of Hylidae

Hyla sanchiangensis (Anura: Hylidae) is endemic to China and is distributed across Anhui, Zhejiang, Fujian, Guangdong, Guangxi, Hunan, and Guizhou provinces. The mitogenomes of H. sanchiangensis from two different sites (Jinxiu, Guangxi, and Wencheng, Zhejiang) were sequenced. Phylogenetic analyses were conducted, including 38 mitogenomes of Hylidae from the NCBI database, and assessed the phylogenetic relationship of H. sanchiangensis within the analyzed dataset. Two mitogenomes of H. sanchiangensis showed the typical mitochondrial gene arrangement with 13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA (tRNA) genes, and one non-coding control region (D-loop). The lengths of the 12S rRNA and 16S rRNA genes from both samples (Jinxiu and Wencheng) were 933 bp and 1604 bp, respectively. The genetic distance (p-distance transformed into percent) on the basis of the mitogenomes (excluding the control region) of the two samples was calculated as 4.4%. Hyla sanchiangensis showed a close phylogenetic relationship with the clade of (H. annectans + H. tsinlingensis), which was supported by ML and BI analyses. In the branch-site model, five positive selection sites were found in the clade of Hyla and Dryophytes: Cytb protein (at position 316), ND3 protein (at position 85), and ND5 protein (at position 400) have one site, respectively, and two sites in ND4 protein (at positions 47 and 200). Based on the results, we hypothesized that the positive selection of Hyla and Dryophytes was due to their experience of cold stress in historical events, but more evidence is needed to support this conclusion.

Morphological variation and its correlation with bioclimatic factors in Odorrana graminea sensu stricto

The large green cascade frog (Odorrana graminea sensu stricto) shows significant genetic differentiation in China, forming western, southern, and eastern clades. However, the morphological differentiation among the three clades is unclear, and the influence of bioclimatic factors on morphological variation among clades is unknown. Based on 20 morphological traits of 309 specimens from 28 localities, the present study explored the morphological differentiation and variation among clades and their correlation with bioclimatic factors through the multivariate statistical analysis. The results of the present study showed that O. graminea sensu stricto was divided into western, southern, and eastern morphological groups, and the gene flow between neighboring populations had caused an individual misidentification. With the three-step terrain and population distribution latitude and humidity, the annual mean temperature (Bio1) was significantly different between the southern and eastern–western clades; the maximum temperature of the hottest month (Bio5) was significantly different between the southern and western clades, and the mean temperature of the wettest quarter (Bio8) and the precipitation seasonality (Bio15) were significantly different between the eastern and western–southern clades. The southern clade that was affected by a high temperature had a smaller body size and larger sensory organs, and the eastern clade distributed in highly humid areas had a larger body size and smaller sensory organs. Moreover, the annual mean temperature range (Bio7) was the dominant factor in the morphological variation of O. graminea sensu stricto, and it had significant negative correlations with seven traits of male frogs and four traits of female frogs. The effect of precipitation factors on the morphological differentiation of each clade remained unclear. The local adaptation caused by climatic differences was the main reason for the morphological differentiation among clades. These findings will help us to understand amphibians’ abilities to adapt to environmental variation.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from İyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Şavşat town of Artvin to Ardeşen town in Rize. In addition, the two species coexist in Pazar, Hemşin, and Çamlıhemşin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Unraveling global impacts of climate change on amphibians distributions: A life-history and biogeographic-based approach

Climate change can affect species distribution patterns in three different ways: pushing them to disperse to new suitable areas, forcing them to adapt to novel climatic conditions, or driving them to extinction. However, the biological and geographical traits that lead to these different responses remain poorly explored. Here, we evaluated how ecological and biogeographic traits influence amphibians’ response to climate change. We performed a systematic review searching for studies that evaluated the effects of future climate change on amphibian’s distribution. Our research returned 31 articles that projected the distribution of 331 amphibians. Our results demonstrate that species inhabiting an elevation above 515 m will lose a significant portion of their climatically suitable area. We also found that as isothermality increases, the amount of area suitable in response to climate change also increases. Another important discovery was that as the size of the baseline area increases, the greater must be the loss of climatically suitable areas. On the other hand, species with very small areas tend to keep their current climatically suitable area in the future. Furthermore, our results indicate that species that inhabit dry habitats tend to expand their suitable area in response to climate change. This result can be explained by the environmental characteristics of these habitats, which tend to present extreme seasonal climates with well-defined periods of drought and rain. We also found that anurans that inhabit exclusively forests are projected to lose a greater portion of their suitable areas, when compared to species that inhabit both forest and open areas, wetlands, and dry and rupestrian environments. The biogeographical realm also influenced anuran’s range shifts, with Afrotropic and Nearctic species projected to expand their geographical ranges. The assessment of climate change effects on amphibian distribution has been the focus of a growing number of studies. Despite this, some regions and species remain underrepresented. Current literature evaluates about 4% of the 7,477 species of Anura and 8% of the 773 species of Caudata and some regions rich in amphibian species remain severely underrepresented, such as Madagascar. Thus, future studies should focus on regions and taxas that remain underrepresented.

Historical biogeography highlights the role of Miocene landscape changes on the diversification of a clade of Amazonian tree frogs

The diversification processes underlying why Amazonia hosts the most species-rich vertebrate fauna on earth remain poorly understood. We studied the spatio-temporal diversification of a tree frog clade distributed throughout Amazonia (Anura: Hylidae: Osteocephalus, Tepuihyla, and Dryaderces) and tested the hypothesis that Miocene mega wetlands located in western and central Amazonia impacted connectivity among major biogeographic areas during extensive periods. We assessed the group’s diversity through DNA-based (16S rRNA) species delimitation to identify Operational Taxonomic Units (OTUs) from 557 individuals. We then selected one terminal for each OTU (n = 50) and assembled a mitogenomic matrix (~14,100 bp; complete for 17 terminals) to reconstruct a Bayesian, time-calibrated phylogeny encompassing nearly all described species. Ancestral area reconstruction indicates that each genus was restricted to one of the major Amazonian biogeographic areas (western Amazonia, Guiana Shield and Brazilian Shield, respectively) between ~10 and 20 Mya, suggesting that they diverged and diversified in isolation during this period around the Pebas mega wetland. After 10 Mya and the transition to the modern configuration of the Amazon River watershed, most speciation within each genus continued to occur within each area. In Osteocephalus, only three species expanded widely across Amazonia (< 6 Mya), and all were pond-breeders. Species with other breeding modes remained mostly restricted to narrow ranges. The spectacular radiation of Osteocephalus was probably driven by climatic stability, habitat diversity and the acquisition of new reproductive modes along the Andean foothills and western Amazonia. Our findings add evidence to the importance of major hydrological changes during the Miocene on biotic diversification in Amazonia.

Unravelling the cryptic diversity and evolution of the dwarf swamp frog Pseudopaludicola mystacalis (Anura, Leptodactylidae) in open habitats of South America

Many studies on the diversity of the South American biota support the role of ecological and geological events as main drivers of species diversification. For many groups, geomorphological events are the key drivers of diversification, while the influence of Pleistocene climate oscillations is prominent for others. To precisely indicate which events were key for the development of the astonishing biodiversity in South America, studies on widely distributed species are paramount. One such species, the dwarf swamp frog Pseudopaludicola mystacalis (Leptodactylidae, Leiuperinae), is widely distributed in open habitats of South America and we herein investigate population differentiation and diversification in this species across its geographic range. We sequenced a 1374 bp mtDNA fragment from 64 specimens across 25 localities. We used population assignment and species delimitation methods to assess genetic structure and lineage limits across the species distribution. We estimated, for each lineage, intraspecific diversity, divergence times, and demographic histories. Our results recovered ten lineages with up to 5% of genetic divergence among them. Diversification occurred mainly during the Tertiary, suggesting that Miocene-Pliocene topographic events had a major influence on the diversification of P. mystacalis. Pleistocene climatic oscillations also played a role on evolutionary history of P. mystacalis, causing demographic changes in one lineage.

Genetic diversity and demography of Bufo japonicus and B. torrenticola (Amphibia: Anura: Bufonidae) influenced by the Quaternary climate

The Quaternary climate affected the present species richness and geographic distribution patterns of amphibians by limiting their activities during the glacial period. The present study examined the phylogenetic relationships of Japanese toads (Bufo japonicus and B. torrenticola) and the demography of each lineage from the past to the present based on mitochondrial sequences and ecological niche models. Japanese toads are a monophyletic group with two main clades (clades A and B). Clade A represents B. j. formosus, including three clades (clades A1, A2, and A3). Clade B contains three clades, two of which corresponded to B. j. japonicus (clades B1 and B2) and the other to B. torrenticola. Clade B2 and B. torrenticola made a sister group, and, thus, B. j. japonicus is paraphyletic. Clades A and B diverged in the late Miocene 5.7 million years ago (Mya) during the period when the Japanese archipelago was constructed. The earliest divergence between the three clades of clade A was estimated at 1.8 Mya. Clades A1 and A2 may have diverged at 0.8 Mya, resulting from the isolation in the multiple different refugia; however, the effects of the glacial climate on the divergence events of clade A3 are unclear. Divergences within clade B occurred from the late Pliocene to the early Pleistocene (3.2-2.2 Mya). Niche similarity between the parapatric clade in clade B (clades B1 and B2) indicated their allopatric divergence. It was suggested that niche segregation between B. japonicus and B. torrenticola contributed to a rapid adaptation of B. torrenticola for lotic breeding. All clade of Japanese toads retreated to each refugium at a low elevation in the glacial period, and effective population sizes increased to construct the current populations after the Last Glacial Maximum. Furthermore, we highlight the areas of climate stability from the last glacial maximum to the present that have served as the refugia of Japanese toads and, thus, affected their present distribution patterns.

Phenotypic and genomic diversification with isolation by environment along elevational gradients in a neotropical treefrog

Understanding how geographic and environmental heterogeneity drive local patterns of genetic variation is a major goal of ecological genomics and a key question in evolutionary biology. The tropical Andes and inter-Andean valleys are shaped by markedly heterogeneous landscapes, where species experience strong selective processes. We examined genome-wide SNP data together with behavioural and ecological traits (mating calls and body size) known to contribute to genetic isolation in anurans in the banana tree-dwelling frog, Boana platanera, distributed across an environmental gradient in Central Colombia (northern South America). Here, we analysed the relationships between environmentally (temperature and precipitation) associated genetic and phenotypic differentiation and the potential drivers of isolation by environment along an elevation gradient. We identified candidate SNPs associated with temperature and body size, which follow a clinal pattern of genome-wide differentiation tightly coupled with phenotypic variation: as elevation increases, B. platanera exhibits larger body size and longer call duration with more pulses but lower pulse rate and frequency. Thus, the environmental landscape has rendered a scenario where isolation by environment and candidate loci show concordance with phenotypic divergence in this tropical frog along an elevation gradient in the Colombian Andes. Our study sets the basis for evaluating the role of temperature in the genetic structure and local adaptation in tropical treefrogs and its putative effect on life cycle (embryos, tadpoles, adults) along elevation gradients.

Updating the bionomy and geographical distribution of Anopheles (Nyssorhynchus) albitarsis F: A vector of malaria parasites in northern South America

Anopheles albitarsis F is a putative species belonging to the Albitarsis Complex, recognized by rDNA, mtDNA, partial white gene, and microsatellites sequences. It has been reported from the island of Trinidad, Venezuela and Colombia, and incriminated as a vector of malaria parasites in the latter. This study examined mitochondrially encoded cytochrome c oxidase I (MT-CO1) sequences of An. albitarsis F from malaria-endemic areas in Colombia and Venezuela to understand its relations with other members of the Complex, revised and update the geographical distribution and bionomics of An. albitarsis F and explore hypotheses to explain its phylogenetic relationships and geographical expansion. Forty-five MT-CO1 sequences obtained in this study were analyzed to estimate genetic diversity and possible evolutionary relationships. Sequences generated 37 haplotypes clustered in a group where the genetic divergence of Venezuelan populations did not exceed 1.6% with respect to Colombian samples. Anopheles albitarsis F (π = 0.013) represented the most recent cluster located closer to An. albitarsis I (π = 0.009). Barcode gap was detected according to Albitarsis Complex lineages previously reported (threshold 0.014–0.021). Anopheles albitarsis F has a wide distribution in northern South America and might play an important role in the transmission dynamics of malaria due to its high expansion capacity. Future studies are required to establish the southern distribution of An. albitarsis F in Venezuela, and its occurrence in Guyana and Ecuador.

Paso doble: A two-step Late Pleistocene range expansion in the Tyrrhenian tree frog Hyla sarda

The Tyrrhenian tree frog, Hyla sarda, is an amphibian endemic to the Tyrrhenian islands (Western Mediterranean). Previous investigations of its Pleistocene evolutionary history suggested that it colonised the northern portion of its current range, through a spatial diffusion process from the Sardinia island, during the last glaciation. However, southern and northern portions of the species' range experienced markedly different climatic conditions during the Late Pleistocene, suggesting the possibility of an unusual two-step process of demographic expansion. Here, we use Bayesian phylogeographic approaches to locate the ancestral area in Sardinia and to characterise better the demographic component of this expansion event. These analyses located the ancestral area for H. sarda populations along the central-eastern coast of the Sardinia island, within an area previously shown to host suitable bioclimatic conditions for H. sarda populations throughout the Late Pleistocene. Historical demographic reconstructions clearly showed that a two-step process of demographic growth fits well the data, with northern populations expanding later than Sardinia populations. The harsher climatic conditions occurred in northern islands during the glacial epoch, as compared to Sardinia, likely delayed tree frog colonisation of northern territories, and the associated demographic growth.

Conservation genetics of yellow-bellied toads (Bombina variegata): a matter of geographical scale and isolation

Amphibian populations world-wide are threatened by declines and extinctions mainly due to habitat loss and fragmentation. Habitat fragmentation threatens the yellow-bellied toad Bombina variegata in the northern and western regions of its distribution where it is strictly protected. We studied the genetic structure and diversity of populations at three geographical scales using microsatellite loci to detect potential threats for population persistence. At the local scale, we sampled four neighbouring localities at 1–2.6 km distance to detect effects of short-term (decades) fragmentation on connectivity. At the regional scale, five additional localities in the mountains of the Westerwald (Rhineland-Palatinate, Germany) were studied at up to 50.1 km distance to analyse genetic diversity and population structure. At the continental scale, we included data from regions in the northern distribution with fragmented populations (Hesse and Lower Saxony, Germany) and more continuous populations in the South (Alsace, France; Geneva, Switzerland; Trentino, Italy) to evaluate variation of genetic diversity. At the local scale, short-term fragmentation caused significant genetic differentiation between breeding assemblages only 1.4 km apart from each other. At the regional scale, we found notable genetic distance among localities. At the continental scale, we identified Alsace, Trentino and Geneva in the South as regions with low genetic structuring and high allelic richness, and the northern remaining regions in Germany as deeply structured with reduced allelic richness. We suggest that reduced genetic diversity and habitat fragmentation in northern regions makes these populations particularly vulnerable to decline. In conclusion, informed conservation management of B. variegata should focus on measures maintaining or improving connectivity among neighbouring populations.

The Influence of Environmental Variation on the Genetic Structure of a Poison Frog Distributed Across Continuous Amazonian Rainforest

Biogeographic barriers such as rivers have been shown to shape spatial patterns of biodiversity in the Amazon basin, yet relatively little is known about the distribution of genetic variation across continuous rainforest. Here, we characterize the genetic structure of the brilliant-thighed poison frog (Allobates femoralis) across an 880-km-long transect along the Purus-Madeira interfluve south of the Amazon river, based on 64 individuals genotyped at 7609 single-nucleotide polymorphism (SNP) loci. A population tree and clustering analyses revealed 4 distinct genetic groups, one of which was strongly divergent. These genetic groups were concomitant with femoral spot coloration differences, which was intermediate within a zone of admixture between two of the groups. The location of these genetic groups did not consistently correspond to current ecological transitions between major forest types. A multimodel approach to quantify the relative influence of isolation-by-geographic distance (IBD) and isolation-by-environmental resistance (IBR) nevertheless revealed that, in addition to a strong signal of IBD, spatial genetic differentiation was explained by IBR primarily linked to dry season intensity (r2 = 8.4%) and canopy cover (r2 = 6.4%). We show significant phylogenetic divergence in the absence of obvious biogeographical barriers and that finer-scaled measures of genetic structure are associated with environmental variables also known to predict the density of A. femoralis.

The roles of climate, geography and natural selection as drivers of genetic and phenotypic differentiation in a widespread amphibian Hyla annectans (Anura: Hylidae)

The role of geological events and Pleistocene climatic fluctuations as drivers of current patterns of genetic variation in extant species has been a topic of continued interest among evolutionary biologists. Nevertheless, comprehensive studies of widely distributed species are still rare, especially from Asia. Using geographically extensive sampling of many individuals and a large number of nuclear single nucleotide polymorphisms (SNPs), we studied the phylogeography and historical demography of Hyla annectans populations in southern China. Thirty-five sampled populations were grouped into seven clearly defined genetic clusters that closely match phenotype-based subspecies classification. These lineages diverged 2.32-5.23 million years ago (Ma), a timing that closely aligns with the rapid and drastic uplifting of the Qinghai-Tibet Plateau and adjacent southwest China. Demographic analyses and species distribution models indicate that different populations of this species have responded differently to past climatic changes. In the Hengduan Mountains, most populations experienced a bottleneck, whereas the populations located outside of the Hengduan Mountains have gradually declined in size since the end of the last glaciation. In addition, the levels of phenotypic and genetic divergence were strongly correlated across major clades. These results highlight the combined effects of geological events and past climatic fluctuations, as well as natural selection, as drivers of contemporary patterns of genetic and phenotypic variation in a widely distributed anuran in Asia.

Evolutionary principles guiding amphibian conservation

The Anthropocene has witnessed catastrophic amphibian declines across the globe. A multitude of new, primarily human-induced drivers of decline may lead to extinction, but can also push species onto novel evolutionary trajectories. If these are recognized by amphibian biologists, they can be engaged in conservation actions. Here, we summarize how principles stemming from evolutionary concepts have been applied for conservation purposes, and address emerging ideas at the vanguard of amphibian conservation science. In particular, we examine the consequences of increased drift and inbreeding in small populations and their implications for practical conservation. We then review studies of connectivity between populations at the landscape level, which have emphasized the limiting influence of anthropogenic structures and degraded habitat on genetic cohesion. The rapid pace of environmental changes leads to the central question of whether amphibian populations can cope either by adapting to new conditions or by shifting their ranges. We gloomily conclude that extinction seems far more likely than adaptation or range shifts for most species. That said, conservation strategies employing evolutionary principles, such as selective breeding, introduction of adaptive variants through translocations, ecosystem interventions aimed at decreasing phenotype-environment mismatch, or genetic engineering, may effectively counter amphibian decline in some areas or for some species. The spread of invasive species and infectious diseases has often had disastrous consequences, but has also provided some premier examples of rapid evolution with conservation implications. Much can be done in terms of setting aside valuable amphibian habitat that should encompass both natural and agricultural areas, as well as designing protected areas to maximize the phylogenetic and functional diversity of the amphibian community. We conclude that an explicit consideration and application of evolutionary principles, although certainly not a silver bullet, should increase effectiveness of amphibian conservation in both the short and long term.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

RNA-seq analysis provides insight into molecular adaptations of Andrias davidianus

The Chinese giant salamander Andrias davidianus is regarded as an ideal model for studying local adaptations, such as longevity, tolerance to starvation, and cutaneous respiration. Transcriptome analysis is useful for studying the large and complex genomes of amphibians. Based on the coding gene set of adult A. davidianus, dozens of A. davidianus-specific genes were identified and three signaling pathway (JAK-STAT, HIF-1, and FoxO) genes were expanded as compared with other amphibians. The results of the pathway analysis of A. davidianus-specific genes indicated that the molecular adaptation of A. davidianus may have required a more rapid evolution of the immune system. Additionally, for the first time, the gene expressions in different parts of the skin tissue were compared. The results of the comparison analysis demonstrated that lateral skin could be more focused on mucus secretion, dorsal skin on immunity and melanogenesis, and abdominal skin on water and salt metabolism. This study provides the first insight into studying longevity and starvation tolerance in A. davidianus, and offers a basis for further investigation of the molecular mechanisms of adaptations in amphibians.

Genetic diversity and structure of two endangered mole salamander species of the Trans-Mexican Volcanic Belt

The most important factor leading to amphibian population declines and extinctions is habitat degradation and destruction. To help prevent further extinctions, studies are needed to make appropriate conservation decisions in small and fragmented populations. The goal of this study was to provide data from the population genetics of two micro-endemic mole salamanders from the Trans-Mexican Volcanic Belt. Nine microsatellite markers were used to study the population genetics of 152 individuals from twoAmbystomaspecies. We sampled 38 individuals in two localities forA. altamiraniandA. rivualre. We found medium to high levels of genetic diversity expressed as heterozygosity in the populations. However, all the populations presented few alleles per locus and genotypes. We found strong genetic structure between populations for each species. Effective population size was small but similar to that of the studies from other mole salamanders with restricted distributions or with recently fragmented habitats. Despite the medium to high levels of genetic diversity expressed as heterozygosity, we found few alleles, evidence of a genetic bottleneck and that the effective population size is small in all populations. Therefore, this study is important to propose better management plans and conservation efforts for these species.

Mitochondrial Introgression, Color Pattern Variation, and Severe Demographic Bottlenecks in Three Species of Malagasy Poison Frogs, Genus Mantella

Madagascar is a biodiversity hotspot particularly rich in amphibian diversity and only a few charismatic Malagasy amphibians have been investigated for their population-level differentiation. The Mantellamadagascariensis group is composed of two rainforest and three swamp forest species of poison frogs. We first confirm the monophyly of this clade using DNA sequences of three nuclear and four mitochondrial genes, and subsequently investigate the population genetic differentiation and demography of the swamp forest species using one mitochondrial, two nuclear and a set of nine microsatellite markers. Our results confirm the occurrence of two main mitochondrial lineages, one dominated by Mantellaaurantiaca (a grouping supported also by our microsatellite-based tree) and the other by Mantellacrocea + Mantellamilotympanum. These two main lineages probably reflect an older divergence in swamp Mantella. Widespread mitochondrial introgression suggests a fairly common occurrence of inter-lineage gene flow. However, nuclear admixture seems to play only a limited role in this group, and the analyses of the RAG-1 marker points to a predominant incomplete lineage sorting scenario between all five species of the group, which probably diverged relatively recently. Our demographic analyses show a common, severe and recent demographic contraction, inferred to be in temporal coincidence with the massive deforestation events that took place in the past 1000 years. Current data do not allow to conclusively delimit independent evolutionary units in these frogs, and we therefore refrain to suggest any taxonomic changes.

Phylogeographic patterns and conservation implications of the endangered Chinese giant salamander

Understanding genetic diversity patterns of endangered species is an important premise for biodiversity conservation. The critically endangered salamander Andrias davidianus, endemic to central and southern mainland in China, has suffered from sharp range and population size declines over the past three decades. However, the levels and patterns of genetic diversity of A. davidianus populations in wild remain poorly understood. Herein, we explore the levels and phylogeographic patterns of genetic diversity of wild-caught A. davidianus using larvae and adult collection with the aid of sequence variation in (a) the mitochondrial DNA (mtDNA) fragments (n = 320 individuals; 33 localities), (b) 19 whole mtDNA genomes, and (c) nuclear recombinase activating gene 2 (RAG2; n = 88 individuals; 19 localities). Phylogenetic analyses based on mtDNA datasets uncovered seven divergent mitochondrial clades (A-G), which likely originated in association with the uplifting of mountains during the Late Miocene, specific habitat requirements, barriers including mountains and drainages and lower dispersal ability. The distributions of clades were geographic partitioned and confined in neighboring regions. Furthermore, we discovered some mountains, rivers, and provinces harbored more than one clades. RAG2 analyses revealed no obvious geographic patterns among the five alleles detected. Our study depicts a relatively intact distribution map of A. davidianus clades in natural species range and provides important knowledge that can be used to improve monitoring programs and develop a conservation strategy for this critically endangered organism.

Landscape genetics reveal broad and fine-scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east-west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east-west split occurred ~13-18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5-11 kya during a period of extreme drought cycles. A range-wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid-Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.

Phylogeny and species delimitation of near Eastern Neurergus newts (Salamandridae) based on genome-wide RADseq data analysis

We reconstruct the molecular phylogeny of Near Eastern mountain brook newts of the genus Neurergus (family Salamandridae) based on newly determined RADseq data, and compare the outcomes of concatenation-based phylogenetic reconstruction with species-tree inference. Furthermore, we test the current taxonomy of Neurergus (with four species: Neurergus strauchii, N. crocatus, N. kaiseri, and N. derjugini) against coalescent-based species-delimitation approaches of our genome-wide genetic data set. While the position of N. strauchii as sister species to all other Neurergus species was consistent in all of our analyses, the phylogenetic relationships between the three remaining species changed depending on the applied method. The concatenation approach, as well as quartet-based species-tree inference, supported a topology with N. kaiseri as the closest relative to N. derjugini, while full-coalescent species-tree inference approaches supported N. crocatus as sister species of N. derjugini. Investigating the individual signal of gene trees highlighted an extensive variation among gene histories, most likely resulting from incomplete lineage sorting. Coalescent-based species-delimitation models suggest that the current taxonomy might underestimate the species richness within Neurergus and supports seven species. Based on the current sampling, our analysis suggests that N. strauchii, N. derjugini and N. kaiseri might each be subdivided into further species. However, as amphibian species are known to be composed of deep conspecific lineages that do not always warrant species status, these results need to be cautiously interpreted in an integrative taxonomic framework. We hypothesize that the rather shallow divergences detected within N. kaiseri and N. derjugini likely reflect an ongoing speciation process and thus require further investigation. On the contrary, the much deeper genetic divergence found between the two morphologically and geographically differentiated subspecies of N. strauchii leads us to propose that N. s. barani should be considered a distinct species, Neurergus barani Öz, 1994.

Environment and evolutionary history shape phylogenetic turnover in European tetrapods

Phylogenetic turnover quantifies the evolutionary distance among species assemblages and is central to understanding the main drivers shaping biodiversity. It is affected both by geographic and environmental distance between sites. Therefore, analyzing phylogenetic turnover in environmental space requires removing the effect of geographic distance. Here, we apply a novel approach by deciphering phylogenetic turnover of European tetrapods in environmental space after removing geographic land distance effects. We demonstrate that phylogenetic turnover is strongly structured in environmental space, particularly in ectothermic tetrapods, and is well explained by macroecological characteristics such as niche size, species richness and relative phylogenetic diversity. In ectotherms, rather recent evolutionary processes were important in structuring phylogenetic turnover along environmental gradients. In contrast, early evolutionary processes had already shaped the current structure of phylogenetic turnover in endotherms. Our approach enables the disentangling of the idiosyncrasies of evolutionary processes such as the degree of niche conservatism and diversification rates in structuring biodiversity.

Anuran diversity indicates that Caatinga relictual Neotropical forests are more related to the Atlantic Forest than to the Amazon

The relationships among the morphoclimatic domains of South America have been a major biogeographical issue of recent years. Palynological, geological and phytogeographical data suggest that the Amazon Forest and the Atlantic Forest were connected during part of the Tertiary and Quaternary periods. This study uses amphibians as model organisms to investigate whether relict northeastern forests are a transition between the Amazon Forest and the Atlantic Forest. We compiled matrices of species composition for four different phytogeographic formations and "Brejos de Altitude," and analyzed them using clustering methods and Cladistic Analysis of Distributions and Endemism. Our results indicate that the anurofauna of these northeastern forest relicts is most similar in composition to the areas of the Atlantic Forest included in this study, and most dissimilar to the Amazon Forest, which leads us to affirm that events of biotic exchange were more frequent within the Atlantic Forest areas.

Phylogeography of Ticks (Acari: Ixodida)

Improved understanding of tick phylogeny has allowed testing of some biogeographical patterns. On the basis of both literature data and a meta-analysis of available sequence data, there is strong support for a Gondwanan origin of Ixodidae, and probably Ixodida. A particularly strong pattern is observed for the genus Amblyomma, which appears to have originated in Antarctica/southern South America, with subsequent dispersal to Australia. The endemic Australian lineages of Ixodidae (no other continent has such a pattern) appear to result from separate dispersal events, probably from Antarctica. Minimum ages for a number of divergences are determined as part of an updated temporal framework for tick evolution. Alternative hypotheses for tick evolution, such as a very old Pangean group, a Northern hemisphere origin, or an Australian origin, fit less well with observed phylogeographic patterns.

Past vicariance promoting deep genetic divergence in an endemic frog species of the Espinhaço Range in Brazil: The historical biogeography of Bokermannohyla saxicola (Hylidae)

The highland endemic frog Bokermannohyla saxicola occurs within the Espinhaço Range, the most extensive and continuous orogenic belt of the Brazilian territory, located in southeastern Brazil. We used mitochondrial DNA markers to test for spatial structure, to investigate the likely influence of past vicariant events, to evaluate demographic dynamics along the species range, and to understand the role of habitat discontinuities in promoting connectivity and diversity along the range. We found four major monophyletic lineages, each one associated with distinct mountain tops. The divergence time found between the four main clades clearly pre-dated the Pleistocene, except for the most recent separation. We observed no signs of population expansion for most of the sampling sites along the range, and a higher genetic diversity in the most continuous and central highland plateau, compared to smaller marginal regions. The Espinhaço Range harbors four deeply divergent lineages of B. saxicola within areas restricted by barriers for millions of years. These relatively isolated populations were kept apart by discontinuities represented by lowland habitats between mountain tops. Most of the lineage divergences occurred earlier than the Pleistocene, thus they cannot be solely explained by climatic oscillations of this epoch. However, within-lineage divergence times were all dated from the Pleistocene, suggesting an important effect in population dynamics. We also suggest that some marginal populations like those from Serra Negra and Serra de Itacambira can be the result of recent colonization events. Finally, in the southern Espinhaço region, the most continuous central highland area shows greater genetic diversity than the marginal discontinuous areas, where we have also observed a significant correlation between genetic and geographic distances. Bokermannohyla saxicola is a good model to study the biogeography of the Espinhaço Range because its high genetic structure reflects ancient as well as recent geological/climatic events, with important implications for conservation.

Connections and containers: Using genetic data to understand how watershed evolution and human activities influence cutthroat trout biogeography

Species with large geographic distributions often exhibit complex patterns of diversity that can be further complicated by human activities. Cutthroat trout (Oncorhynchus clarkii) are one of the most widely distributed freshwater fish species in western North America exhibiting substantial phenotypic and genetic variability; however, fish stocking practices have translocated populations outside of their native range and may have obscured intraspecific boundaries. This study focuses on cutthroat trout populations representing three distinct evolutionary clades that are found intermixed within a contact zone between the Bonneville and upper Snake River watersheds in the western United States. We used mitochondrial and microsatellite genetic data, as well as historical stocking records, to evaluate whether populations of cutthroat trout in the contact zone are native or are introduced. We found significant genetic differentiation and fine-scale genetic population structure that was organized primarily by watershed boundaries. While we detected increased genetic diversity in some areas in close proximity to the greatest number of stocking events, the highly organized population structure both within and between areas of the contact zone indicates that the populations are native to the watersheds. Intermixing of distinct evolutionary lineages of cutthroat trout appears to be the result of historical connections between paleodrainages. Our analyses provide a context for understanding how genetic data can be used to assess the status of populations as native or introduced.

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.

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency-dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color-divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome-wide markers from a double-digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.