Search

The geographic distribution of genetic variation within a species reveals information about its evolutionary history, including responses to historical climate change and dispersal ability across various habitat types. We combine genetic data from salamander species with geographic, climatic, and life history data collected from open-source online repositories to develop a machine learning model designed to identify the traits that are most predictive of unrecognized genetic lineages. We find evidence of hidden diversity distributed throughout the clade Caudata that is largely the result of variation in climatic variables. We highlight some of the difficulties in using machine-learning models on open-source data that are often messy and potentially taxonomically and geographically biased., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Parsons et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Recognizing the influence of pathogen diversity on infection dynamics is crucial for mitigating emerging infectious diseases. Characterising such diversity is often complex, for instance when multiple pathogen variants exist that interact differently with the environment and host. Here, we explore genotypic and phenotypic variation of Batrachochytrium salamandrivorans (Bsal), an emerging fungal pathogen that is driving declines among an increasing number of European amphibian species. For thirteen isolates, spanning most of the known temporal and geographical Bsal range in Europe, we mapped phenotypic diversity through numerous measurements that describe varying reproductive rates in vitro across a range of temperatures. Bsal isolates are revealed to have different thermal optima and tolerances, with phenotypic variation correlating with genomic diversity. Using a mechanistic niche model of the fire salamander (Salamandra salamandra) as an example, we illustrate how host steady-state body temperature and Bsal thermal range variation may influence pathogen growth through space and time across Europe. Our combined findings show how the identity of emergent pathogen variants may strongly influence when and which host populations are most at risk., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Kelly et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Due to traditional classification methods' limitations, some cryptic species remain undiscovered. To better explore the existence of the Schrenck salamander (Salamandrella tridactyla, a cryptic species of Siberian salamander S. keyserlingii) in China, we conducted a molecular phylogenetic analysis to confirm the taxonomic relationship among Salamandrella species and investigate genetic variation. We used complete sequences of the mitochondrial COI gene from 65 specimens collected across a wide range in Northeastern China. Thirty-five haplotypes were obtained from six populations. They showed medium-high haplotype diversity (Hd) and low nucleotide polymorphism (π). The phylogenetic tree and haplotype network analysis revealed that populations from Greater Khingan Ridge (Huma: HM) and Lesser Khingan Ridge (Tieli: TL) belong to S. keyserlingii, while populations from Changbai Mountain (Shangzhi-zhuziying: SZ, Shangzhi-cuijia: SC, Hailin: HL, and Baishan: BS) belong to S. tridactyla. This indicates the monophyly of Salamandrella and each of the two species. There was a substantial level of genetic differentiation between different species and within populations of the same species. This differentiation was significantly related to geographical distance. At last, the mismatch distribution and neutrality analyses indicated that the TL populations have undergone expansion of history. The study supplements the distributional range of Schrenck salamander. And it provides a theoretical basis for species conservation of Salamandrella species., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Liu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The molecular detection of DNA fragments that are shed into the environment (eDNA) has become an increasingly applied tool used to inventory biological communities and to perform targeted species surveys. This method is particularly useful in habitats where it is difficult or not practical to visually detect or trap the target organisms. Central Texas Eurycea salamanders inhabit both surface and subterranean aquatic environments. Subterranean surveys are challenging or infeasible, and the detection of salamander eDNA in water samples is an appealing survey technique for these situations. Here, we develop and validate an eDNA assay using quantitative PCR for E. chisholmensis, E. naufragia, and E. tonkawae. These three species are federally threatened and constitute the Septentriomolge clade that occurs in the northern segment of the Edwards Aquifer. First, we validated the specificity of the assay in silico and with DNA extracted from tissue samples of both target Septentriomolge and non-target amphibians that overlap in distribution. Then, we evaluated the sensitivity of the assay in two controls, one with salamander-positive water and one at field sites known to be occupied by Septentriomolge. For the salamander-positive control, the estimated probability of eDNA occurrence (ψ) was 0.981 (SE = 0.019), and the estimated probability of detecting eDNA in a qPCR replicate (p) was 0.981 (SE = 0.011). For the field control, the estimated probability of eDNA occurring at a site (ψ) was 0.938 (95% CRI: 0.714-0.998). The estimated probability of collecting eDNA in a water sample (θ) was positively correlated with salamander relative density and ranged from 0.371 (95% CRI: 0.201-0.561) to 0.999 (95% CRI: 0.850- > 0.999) among sampled sites. Therefore, sites with low salamander density require more water samples for eDNA evaluation, and we determined that our site with the lowest estimated θ would require seven water samples for the cumulative collection probability to exceed 0.95. The estimated probability of detecting eDNA in a qPCR replicate (p) was 0.882 (95% CRI: 0.807-0.936), and our assay required two qPCR replicates for the cumulative detection probability to exceed 0.95. In complementary visual encounter surveys, the estimated probability of salamanders occurring at a known-occupied site was 0.905 (SE = 0.096), and the estimated probability of detecting salamanders in a visual encounter survey was 0.925 (SE = 0.052). We additionally discuss future research needed to refine this method and understand its limitations before practical application and incorporation into formal survey protocols for these taxa., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Adcock et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In Tennessee, populations of the state endangered Streamside Salamander (Ambystoma barbouri) are in decline as their distribution lies mostly within rapidly developing areas in the Nashville Basin. Information regarding the partitioning of genetic variation among populations of A. barbouri and the taxonomic status of these populations relative to northern populations and their congener, the Small-mouthed Salamander (A. texanum), have important implications for management and conservation of this species. Here we combined mitochondrial sequencing and genome-wide single nucleotide polymorphism (SNP) data generated using Genotyping-by-Sequencing (GBS) to investigate patterns of genetic variation within Tennessee populations of A. barbouri, to assess their relationship to populations in Kentucky, and to examine their phylogenetic relationship to the closely related A. texanum. Results from phylogenetic reconstructions reveal a complex history of Tennessee A. barbouri populations with regards to northern populations, unisexual A. barbouri, and A. texanum. Patterns of mitochondrial sequence variation suggest that A. barbouri may have originated within Tennessee and expanded north multiple times into Kentucky, Ohio, Indiana, and West Virginia. Phylogenetic reconstructions based on genome-wide SNP data contradict results based on mitochondrial DNA and correspond to geographic and taxonomic boundaries. Variation in allele frequencies at SNP genotypes, as identified by multivariate analyses and Bayesian assignment tests, identified three evolutionary significant units (ESUs) for A. barbouri within Tennessee. Collectively, these results emphasize the need for prioritizing conservation needs for Tennessee populations of A. barbouri to ensure the long-term persistence of this species., Competing Interests: The authors have declared no competing interests.

Determining the genetic consequences of both historical and contemporary events can clarify the effects of the environment on population connectivity and inform conservation decisions. Historical events (like glaciations) and contemporary factors (like logging) can disrupt gene flow between populations. This is especially true among species with specialized ecological requirements and low dispersal ability, like amphibians. We test for the genetic consequences of historical and contemporary disturbances in the coastal giant salamander (Dicamptodon tenebrosus) in the Pacific Northwest of the United States. We consider predictions based on the contemporary landscape (habitat connectivity, logging, forest fires, and topography), in addition to relatively ancient post-Pleistocene range expansion (following the last glacial retreat). To assess local versus larger-scale effects, we sampled 318 individuals across 23 sites, which were clustered in five sampling regions. Genetic variation was assessed using five microsatellite markers. We found evidence of (i) historical regional isolation, with decreased genetic diversity among more recently colonized northern sites, as well as (ii) high levels of inbreeding and loss of heterozygosity at local scales, despite relatively low overall population differentiation (FST) or strong evidence for population bottlenecks. Genetic diversity was not associated with contemporary disturbances (logging or fire), and there were no detectable effects on the genetic connectivity of populations based on intervening landscape features (habitat fragmentation and topography). However, lower genetic diversity in more northern regions indicates a lag in recovery of genetic diversity following post-Pleistocene expansion. Additionally, some populations had evidence of having undergone a recent genetic bottleneck or had high inbreeding (FIS) values. Lower genetic diversity in more northern sites means populations may be more vulnerable to future environmental changes, and managing for connectivity alone may not be sufficient given low mobility. Recent apparent reductions in some populations were not clearly linked to anthropogenic disturbances we examined. This suggests the type of disturbances this species is sensitive to may not be well understood., Competing Interests: The authors have declared that no competing interests exist.

Infectious diseases are considered major threats to biodiversity, however strategies to mitigate their impacts in the natural world are scarce and largely unsuccessful. Chytridiomycosis is responsible for the decline of hundreds of amphibian species worldwide, but an effective disease management strategy that could be applied across natural habitats is still lacking. In general amphibian larvae can be easily captured, offering opportunities to ascertain the impact of altering the abundance of hosts, considered to be a key parameter affecting the severity of the disease. Here, we report the results of two experiments to investigate how altering host abundance affects infection intensity in amphibian populations of a montane area of Central Spain suffering from lethal amphibian chytridiomycosis. Our laboratory-based experiment supported the conclusion that varying density had a significant effect on infection intensity when salamander larvae were housed at low densities. Our field experiment showed that reducing the abundance of salamander larvae in the field also had a significant, but weak, impact on infection the following year, but only when removals were extreme. While this suggests adjusting host abundance as a mitigation strategy to reduce infection intensity could be useful, our evidence suggests only heavy culling efforts will succeed, which may run contrary to objectives for conservation., Competing Interests: The authors have declared that no competing interests exist.

Background: Body shape can be both a consequence and cause of a species' evolution and ecology. There are many examples of phenotypes associated with specific ecological niches, likely as a result of specific selective regimes. A classic example of this is the phenotypic change associated with colonization of caves, including body and limb elongation. However, studies explicitly testing for differences in body shape between cave-dwelling and non-cave-dwelling lineages have been limited and so the role of the cave environment in determining morphological characteristics is still not completely understood. Here we examine variation in body shape among 405 individuals representing 20 species in the salamander genus Eurycea (Plethodontidae) and select outgroups exhibiting great diversity in morphology, ecological niche, and life history., Results: After analyzing morphometric data in a phylogenetic context using phylogenetic MANOVA and examination of the phylomorphospace, we found significant differences in body shape among cave-dwelling and non-cave-dwelling species and between aquatic and terrestrial species. Notably, limb elongation and reduced body and tail size characterized cave-dwelling species. Terrestrial species also exhibited elongation of the limbs and digits. We also observed differences in shape variance among paedomorphic and biphasic species. Our results suggest that the functional limitations imposed by habitat and life history played a key role in the evolution of body shape in this group in the context of their phylogenetic history., Competing Interests: The authors have declared that no competing interests exist.

The salamander family Sirenidae is represented by four extant species that are restricted to North America. Sirens are abundant throughout the southern United States and are among the world's largest amphibians, yet the biology, ecology, and phylogeography of this group is poorly-known. In this study we use morphological and genetic evidence to describe a previously unrecognized species from southern Alabama and the Florida panhandle. We name this species the Reticulated Siren, Siren reticulata. Future studies will enable more precise phylogenetic information about S. reticulata and will almost surely reveal additional undescribed species within the family., Competing Interests: The authors have declared that no competing interests exist.

The red-backed salamander (Plethodon cinereus) has long-served as a model system in ecology, evolution, and behavior, and studies surveying molecular variation in this species have become increasingly common over the past decade. However, difficulties are commonly encountered when extending microsatellite markers to populations that are unstudied from a genetic perspective due to high levels of genetic differentiation across this species' range. To ameliorate this issue, we used 454 pyrosequencing to identify hundreds of microsatellite loci. We then screened 40 of our top candidate loci in populations in Virginia, Pennsylvania, and Ohio-including an isolated island population ~ 4.5 km off the shore of Lake Erie (South Bass Island). We identified 25 loci that are polymorphic in a well-studied region of Virginia and 11 of these loci were polymorphic in populations located in the genetically unstudied regions of Ohio and Pennsylvania. Use of these loci to examine patterns of variation within populations revealed that South Bass Island has low diversity in comparison to other sites. However, neither South Bass Island nor isolated populations around Cleveland are inbred. Assessment of variation between populations revealed three well defined genetic clusters corresponding to Virginia, mainland Ohio/Pennsylvania, and South Bass Island. Comparisons of our results to those of others working in various parts of the range are consistent with the idea that differentiation is lower in regions that were once glaciated. However, these comparisons also suggest that well differentiated isolated populations in the formerly glaciated portion of the range are not uncommon. This work provides novel genetic resources that will facilitate population genetic studies in a part of the red-backed salamander's range that has not previously been studied in this manner. Moreover, this work refines our understanding of how neutral variation is distributed in this ecologically important organism.

Environmental DNA (eDNA) has recently been used for detecting the distribution of macroorganisms in various aquatic habitats. In this study, we applied an eDNA method to estimate the distribution of the Japanese clawed salamander, Onychodactylus japonicus, in headwater streams. Additionally, we compared the detection of eDNA and hand-capturing methods used for determining the distribution of O. japonicus. For eDNA detection, we designed a qPCR primer/probe set for O. japonicus using the 12S rRNA region. We detected the eDNA of O. japonicus at all sites (with the exception of one), where we also observed them by hand-capturing. Additionally, we detected eDNA at two sites where we were unable to observe individuals using the hand-capturing method. Moreover, we found that eDNA concentrations and detection rates of the two water sampling areas (stream surface and under stones) were not significantly different, although the eDNA concentration in the water under stones was more varied than that on the surface. We, therefore, conclude that eDNA methods could be used to determine the distribution of macroorganisms inhabiting headwater systems by using samples collected from the surface of the water.

Environmental DNA (eDNA) holds great promise for conservation applications like the monitoring of invasive or imperiled species, yet this emerging technique requires ongoing testing in order to determine the contexts over which it is effective. For example, little research to date has evaluated how seasonality of organism behavior or activity may influence detection probability of eDNA. We applied eDNA to survey for two highly imperiled species endemic to the upper Black Warrior River basin in Alabama, US: the Black Warrior Waterdog (Necturus alabamensis) and the Flattened Musk Turtle (Sternotherus depressus). Importantly, these species have contrasting patterns of seasonal activity, with N. alabamensis more active in the cool season (October-April) and S. depressus more active in the warm season (May-September). We surveyed sites historically occupied by these species across cool and warm seasons over two years with replicated eDNA water samples, which were analyzed in the laboratory using species-specific quantitative PCR (qPCR) assays. We then used occupancy estimation with detection probability modeling to evaluate both the effects of landscape attributes on organism presence and season of sampling on detection probability of eDNA. Importantly, we found that season strongly affected eDNA detection probability for both species, with N. alabamensis having higher eDNA detection probabilities during the cool season and S. depressus have higher eDNA detection probabilities during the warm season. These results illustrate the influence of organismal behavior or activity on eDNA detection in the environment and identify an important role for basic natural history in designing eDNA monitoring programs., Competing Interests: The authors have declared that no competing interests exist.

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Species are a fundamental unit of biodiversity, yet can be challenging to delimit objectively. This is particularly true of species complexes characterized by high levels of population genetic structure, hybridization between genetic groups, isolation by distance, and limited phenotypic variation. Previous work on the Cumberland Plateau Salamander, Plethodon kentucki, suggested that it might constitute a species complex despite occupying a relatively small geographic range. To examine this hypothesis, we sampled 135 individuals from 43 populations, and used four mitochondrial loci and five nuclear loci (5693 base pairs) to quantify phylogeographic structure and probe for cryptic species diversity. Rates of evolution for each locus were inferred using the multidistribute package, and time calibrated gene trees and species trees were inferred using BEAST 2 and *BEAST 2, respectively. Because the parameter space relevant for species delimitation is large and complex, and all methods make simplifying assumptions that may lead them to fail, we conducted an array of analyses. Our assumption was that strongly supported species would be congruent across methods. Putative species were first delimited using a Bayesian implementation of the GMYC model (bGMYC), Geneland, and Brownie. We then validated these species using the genealogical sorting index and BPP. We found substantial phylogeographic diversity using mtDNA, including four divergent clades and an inferred common ancestor at 14.9 myr (95% HPD: 10.8-19.7 myr). By contrast, this diversity was not corroborated by nuclear sequence data, which exhibited low levels of variation and weak phylogeographic structure. Species trees estimated a far younger root than did the mtDNA data, closer to 1.0 myr old. Mutually exclusive putative species were identified by the different approaches. Possible causes of data set discordance, and the problem of species delimitation in complexes with high levels of population structure and introgressive hybridization, are discussed.

We describe an Oligocene newt specimen from western Germany that has gone practically unnoticed in the literature despite having been housed in the Museum für Naturkunde (Berlin) for a century. It is referable to the coeval Chelotriton, but is unusually peramorphic; for many characters it is more peramorphic than all other caudates or even all other lissamphibians. Most noticeable are the position of the jaw joints far caudal to the occiput, the honeycombed sculpture on the maxilla, and the possible presence of a septomaxilla (which would be unique among salamandrids). Referral to a species would require a revision of the genus, but the specimen likely does not belong to the type species. A phylogenetic analysis of nonmolecular characters of Salamandridae, far larger than all predecessors, confirms the referral to Chelotriton. It further loosely associates the Oligocene Archaeotriton and the Miocene Carpathotriton with the extant Lissotriton, though the former may alternatively lie outside Pleurodelinae altogether. The Miocene? I. randeckensis may not belong to the extant Ichthyosaura. The Miocene "Triturus" roehrsi is found neither with the extant Ommatotriton nor with Lissotriton, but inside an Asian/aquatic clade or, when geographic distribution is included as a character, as the sister-group to all other European molgins. The main cause for discrepancies between the results and the molecular consensus is not heterochrony, but adaptations to a life in mountain streams; this is the most likely reason why the Paleocene Koalliella from western Europe forms the sister-group to some or all of the most aquatic extant newts in different analyses. We would like to urge neontologists working on salamandrids to pay renewed attention to the skeleton, not limited to the skull, as a source of diagnostic and phylogenetically informative characters.

The dynamic geologic history of the southeastern United States has played a major role in shaping the geographic distributions of amphibians in the region. In the phylogeographic literature, the predominant pattern of distribution shifts through time of temperate species is one of contraction during glacial maxima and persistence in refugia. However, the diverse biology and ecology of amphibian species suggest that a "one-size-fits-all" model may be inappropriate. Nearly 10% of amphibian species in the region have a current distribution comprised of multiple disjunct, restricted areas that resemble the shape of Pleistocene refugia identified for other temperate taxa in the literature. Here, we apply genetics and spatially explicit climate analyses to test the hypothesis that the disjunct regions of these species ranges are climatic refugia for species that were more broadly distributed during glacial maxima. We use the salamander Plethodon serratus as a model, as its range consists of four disjunct regions in the Southeast. Phylogenetic results show that P. serratus is comprised of multiple genetic lineages, and the four regions are not reciprocally monophyletic. The Appalachian salamanders form a clade sister to all other P. serratus. Niche and paleodistribution modeling results suggest that P. serratus expanded from the Appalachians during the cooler Last Glacial Maximum and has since been restricted to its current disjunct distribution by a warming climate. These data reject the universal applicability of the glacial contraction model to temperate taxa and reiterate the importance of considering the natural history of individual species.

The study of the European plethodontid salamander Hydromantes strinatii using allozyme and mitochondrial markers showed a strong geographical genetic structure. This was likely the outcome of different evolutionary mechanisms leaving their signature despite the effects of the genetic drift due to the low population size typical of this species. Two highly divergent clades were identified in the eastern and central-western part of the range, with further geographic sub-structure. Nuclear and mitochondrial markers substantially recovered the same population groups but were conflicting in reconstructing their relationships. This apparent incongruence highlighted the action of different mechanisms such as secondary contacts and incomplete lineage sorting in originating the observed genetic variation. The troglophilic habit of this species provided the opportunity to show the importance of caves as local refugia in maintaining the genetic diversity through the persistence of local populations. Accordingly, high nucleotide and haplotype diversity, strong geographic genetic structuring and lack of expansion were evidenced. This signature was found in the populations from the Ligurian and Maritime Alps, in agreement with the complex orography and paleoclimatic history of this Mediterranean hotspot.

Osteological variation is recorded among and within four of the most distinctive species of the Mexican salamander genus Chiropterotriton. Analysis of the data is consistent with the monophyletic status of the genus and documents previously unrecorded intraspecific and interspecific variation. Most of the recorded variation involves qualitative and quantitative proportional differences, but four fixed differences constitute autapomorphic states that affirm and diagnose some species (C. dimidiatus, C. magnipes). Osteological variation in 15 characters is analyzed with respect to predictions generated from four hypotheses: 1) phylogeny, 2) adaptation to specific habitats (the four species include cave-dwelling, terrestrial, and arboreal forms), 3) size-free shape, and 4) size. High levels of intraspecific variation suggest that the characters studied are not subject to rigid functional constraints in salamanders, regardless of size. The pattern predicted by the hypothesis based on size differences seen among these four Chiropterotriton species matches most closely the observed pattern of relative skull robustness. Since size change and heterochrony are often associated in plethodontid evolution, it is likely that changes in developmental timing play a role in the morphological transitions among these morphologically diverse taxa. Webbed feet, miniaturization, body shape, and an unusual tarsal arrangement are morphologies exhibited in species of Chiropterotrition that are shown to be homoplastic with other clades of tropical plethodontids. Although extensive homoplasy in salamanders might be seen as a roadblock to unraveling phylogenetic hypotheses, the homologous developmental systems that appear to underlie such homoplasy may reveal common and consistent evolutionary processes at work.

The Chinese giant salamander, Andrias davidianus, is an important species in the course of evolution; however, there is insufficient genomic data in public databases for understanding its immunologic mechanisms. High-throughput transcriptome sequencing is necessary to generate an enormous number of transcript sequences from A. davidianus for gene discovery. In this study, we generated more than 40 million reads from samples of spleen and skin tissue using the Illumina paired-end sequencing technology. De novo assembly yielded 87,297 transcripts with a mean length of 734 base pairs (bp). Based on the sequence similarities, searching with known proteins, 38,916 genes were identified. Gene enrichment analysis determined that 981 transcripts were assigned to the immune system. Tissue-specific expression analysis indicated that 443 of transcripts were specifically expressed in the spleen and skin. Among these transcripts, 147 transcripts were found to be involved in immune responses and inflammatory reactions, such as fucolectin, β-defensins and lymphotoxin beta. Eight tissue-specific genes were selected for validation using real time reverse transcription quantitative PCR (qRT-PCR). The results showed that these genes were significantly more expressed in spleen and skin than in other tissues, suggesting that these genes have vital roles in the immune response. This work provides a comprehensive genomic sequence resource for A. davidianus and lays the foundation for future research on the immunologic and disease resistance mechanisms of A. davidianus and other amphibians.

Bromeliad phytotelmata are frequently used by several Neotropical amphibian taxa, possibly due to their high humidity, microclimatic stability, and role as a refuge from predators. Indeed, the ability of phytotelmata to buffer against adverse environmental conditions may be instrumental in allowing some amphibian species to survive during periods of environmental change or to colonize sub-optimal habitats. Association between bromeliad traits and salamanders has not been studied at a fine scale, despite the intimate association of many salamander species with bromeliads. Here, we identify microhabitat characteristics of epiphytic bromeliads used by two species of the Bolitoglossa morio group (B. morio and B. pacaya) in forest disturbed by volcanic activity in Guatemala. Specifically, we measured multiple variables for bromeliads (height and position in tree, phytotelma water temperature and pH, canopy cover, phytotelma size, leaf size, and tree diameter at breast height), as well as salamander size. We employed a DNA barcoding approach to identify salamanders. We found that B. morio and B. pacaya occurred in microsympatry in bromeliads and that phytotelmata size and temperature of bromeliad microhabitat were the most important factors associated with the presence of salamanders. Moreover, phytotelmata with higher pH contained larger salamanders, suggesting that larger salamanders or aggregated individuals might modify pH. These results show that bromeliad selection is nonrandom with respect to microhabitat characteristics, and provide insight into the relationship between salamanders and this unique arboreal environment.

In response to pervasive sexual selection, protein sex pheromones often display rapid mutation and accelerated evolution of corresponding gene sequences. For proteins, the general dogma is that structure is maintained even as sequence or function may rapidly change. This phenomenon is well exemplified by the three-finger protein (TFP) superfamily: a diverse class of vertebrate proteins co-opted for many biological functions - such as components of snake venoms, regulators of the complement system, and coordinators of amphibian limb regeneration. All of the >200 structurally characterized TFPs adopt the namesake "three-finger" topology. In male red-legged salamanders, the TFP pheromone Plethodontid Modulating Factor (PMF) is a hypervariable protein such that, through extensive gene duplication and pervasive sexual selection, individual male salamanders express more than 30 unique isoforms. However, it remained unclear how this accelerated evolution affected the protein structure of PMF. Using LC/MS-MS and multidimensional NMR, we report the 3D structure of the most abundant PMF isoform, PMF-G. The high resolution structural ensemble revealed a highly modified TFP structure, including a unique disulfide bonding pattern and loss of secondary structure, that define a novel protein topology with greater backbone flexibility in the third peptide finger. Sequence comparison, models of molecular evolution, and homology modeling together support that this flexible third finger is the most rapidly evolving segment of PMF. Combined with PMF sequence hypervariability, this structural flexibility may enhance the plasticity of PMF as a chemical signal by permitting potentially thousands of structural conformers. We propose that the flexible third finger plays a critical role in PMF:receptor interactions. As female receptors co-evolve, this flexibility may allow PMF to still bind its receptor(s) without the immediate need for complementary mutations. Consequently, this unique adaptation may establish new paradigms for how receptor:ligand pairs co-evolve, in particular with respect to sexual conflict.

Antimicrobial peptides have been widely identified from amphibian skins except salamanders. A novel antimicrobial peptide (CFBD) was isolated and characterized from skin secretions of the salamander, Cynops fudingensis. The cDNA encoding CFBD precursor was cloned from the skin cDNA library of C. fudingensis. The precursor was composed of three domains: signal peptide of 17 residues, mature peptide of 41 residues and intervening propeptide of 3 residues. There are six cysteines in the sequence of mature CFBD peptide, which possibly form three disulfide-bridges. CFBD showed antimicrobial activities against Staphylococcus aureus, Bacillus subtilis, Candida albicans and Escherichia coli. This peptide could be classified into family of β-defensin based on its sequence similarity with β-defensins from other vertebrates. Evolution analysis indicated that CFBD was close to fish β-defensin. As far as we know, CFBD is the first β-defensin antimicrobial peptide from salamanders.

Conservation genetics is a powerful tool to assess the population structure of species and provides a framework for informing management of freshwater ecosystems. As lotic habitats become fragmented, the need to assess gene flow for species of conservation management becomes a priority. The eastern hellbender (Cryptobranchus alleganiensis alleganiensis) is a large, fully aquatic paedamorphic salamander. Many populations are experiencing declines throughout their geographic range, yet the genetic ramifications of these declines are currently unknown. To this end, we examined levels of genetic variation and genetic structure at both range-wide and drainage (hierarchical) scales. We collected 1,203 individuals from 77 rivers throughout nine states from June 2007 to August 2011. Levels of genetic diversity were relatively high among all sampling locations. We detected significant genetic structure across populations (Fst values ranged from 0.001 between rivers within a single watershed to 0.218 between states). We identified two genetically differentiated groups at the range-wide scale: 1) the Ohio River drainage and 2) the Tennessee River drainage. An analysis of molecular variance (AMOVA) based on landscape-scale sampling of basins within the Tennessee River drainage revealed the majority of genetic variation (∼94-98%) occurs within rivers. Eastern hellbenders show a strong pattern of isolation by stream distance (IBSD) at the drainage level. Understanding levels of genetic variation and differentiation at multiple spatial and biological scales will enable natural resource managers to make more informed decisions and plan effective conservation strategies for cryptic, lotic species.

The salamander Pachyhynobius shangchengensis (Hynobiidae) is a vulnerable species restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes in the Mount Dabieshan region in southeastern China. However, molecular phylogeography and population genetic structure of P. shangchengensis remain poorly investigated. In this study, we explored the genetic structure and phylogeography of P. shangchengensis based on partial sequences of the mitochondrial DNA (mtDNA) cytochrome b and cytochrome c oxidase subunit I genes. Fifty-one haplotypes and four clades were found among 93 samples. Phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographic regions separated by complicated topography and long distances. The distinct geographic distributions of all lineages and the estimated divergence time suggest spatial and temporal separation coinciding with climatic changes during the Pleistocene. Analysis of molecular variance indicated that most of the observed genetic variation occurred among the four groups, implying long-term interruption of gene flow, and the possible separation of P. shangchengensis into four management units for conservation.

Very little is known about the factors that cause variation in regenerative potential within and between species. Here, we used a genetic approach to identify heritable genetic factors that explain variation in tail regenerative outgrowth. A hybrid ambystomatid salamander (Ambystoma mexicanum x A. andersoni) was crossed to an A. mexicanum and 217 offspring were induced to undergo metamorphosis and attain terrestrial adult morphology using thyroid hormone. Following metamorphosis, each salamander's tail tip was amputated and allowed to regenerate, and then amputated a second time and allowed to regenerate. Also, DNA was isolated from all individuals and genotypes were determined for 187 molecular markers distributed throughout the genome. The area of tissue that regenerated after the first and second amputations was highly positively correlated across males and females. Males presented wider tails and regenerated more tail tissue during both episodes of regeneration. Approximately 66-68% of the variation in regenerative outgrowth was explained by tail width, while tail length and genetic sex did not explain a significant amount of variation. A small effect QTL was identified as having a sex-independent effect on tail regeneration, but this QTL was only identified for the first episode of regeneration. Several molecular markers significantly affected regenerative outgrowth during both episodes of regeneration, but the effect sizes were small (<4%) and correlated with tail width. The results show that ambysex and minor effect QTL explain variation in adult tail morphology and importantly, tail width. In turn, tail width at the amputation plane largely determines the rate of regenerative outgrowth. Because amputations in this study were made at approximately the same position of the tail, our results resolve an outstanding question in regenerative biology: regenerative outgrowth positively co-varies as a function of tail width at the amputation site.

Background: The Siberian salamander (Ranodon sibiricus), distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations., Methodology/principal Findings: We used two genetic markers (mtDNA and microsatellites) to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC) method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines., Conclusion/significance: Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.

Detecting the factors that determine the interruption of gene flow between populations is key to understanding how speciation occurs. In this context, caves are an excellent system for studying processes of colonization, differentiation and speciation, since they represent discrete geographical units often with known geological histories. Here, we asked whether discontinuous calcareous areas and cave systems represent major barriers to gene flow within and among the five species of Sardinian cave salamanders (genus Hydromantes) and whether intraspecific genetic structure parallels geographic distance within and among caves. We generated mitochondrial cytochrome b gene sequences from 184 individuals representing 48 populations, and used a Bayesian phylogeographic approach to infer possible areas of cladogenesis for these species and reconstruct historical and current dispersal routes among distinct populations. Our results show deep genetic divergence within and among all Sardinian cave salamander species, which can mostly be attributed to the effects of mountains and discontinuities in major calcareous areas and cave systems acting as barriers to gene flow. While these salamander species can also occur outside caves, our results indicate that there is a very poor dispersal of these species between separate cave systems.

Comparative genome analysis of non-avian reptiles and amphibians provides important clues about the process of genome evolution in tetrapods. However, there is still only limited information available on the genome structures of these organisms. Consequently, the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes in tetrapods remain poorly understood. We constructed chromosome maps of functional genes for the Chinese soft-shelled turtle (Pelodiscus sinensis), the Siamese crocodile (Crocodylus siamensis), and the Western clawed frog (Xenopus tropicalis) and compared them with genome and/or chromosome maps of other tetrapod species (salamander, lizard, snake, chicken, and human). This is the first report on the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes inferred from comparative genomic analysis of vertebrates, which cover all major non-avian reptilian taxa (Squamata, Crocodilia, Testudines). The eight largest macrochromosomes of the turtle and chicken were equivalent, and 11 linkage groups had also remained intact in the crocodile. Linkage groups of the chicken macrochromosomes were also highly conserved in X. tropicalis, two squamates, and the salamander, but not in human. Chicken microchromosomal linkages were conserved in the squamates, which have fewer microchromosomes than chicken, and also in Xenopus and the salamander, which both lack microchromosomes; in the latter, the chicken microchromosomal segments have been integrated into macrochromosomes. Our present findings open up the possibility that the ancestral amniotes and tetrapods had at least 10 large genetic linkage groups and many microchromosomes, which corresponded to the chicken macro- and microchromosomes, respectively. The turtle and chicken might retain the microchromosomes of the amniote protokaryotype almost intact. The decrease in number and/or disappearance of microchromosomes by repeated chromosomal fusions probably occurred independently in the amphibian, squamate, crocodilian, and mammalian lineages.

Understanding constraints on phenotypic plasticity is central to explaining its evolution and the evolution of phenotypes in general, yet there is an ongoing debate on the classification and relationships among types of constraints. Since plasticity is often a developmental process, studies that consider the ontogeny of traits and their developmental mechanisms are beneficial. We manipulated the timing and reliability of cues perceived by fire salamander larvae for the future desiccation of their ephemeral pools to determine whether flexibility in developmental rates is constrained to early ontogeny. We hypothesized that higher rates of development, and particularly compensation for contradictory cues, would incur greater endogenous costs. We found that larvae respond early in ontogeny to dried conspecifics as a cue for future desiccation, but can fully compensate for this response in case more reliable but contradictory cues are later perceived. Patterns of mortality suggested that endogenous costs may depend on instantaneous rates of development, and revealed asymmetrical costs of compensatory development between false positive and false negative early information. Based on the results, we suggest a simple model of costs of development that implies a tradeoff between production costs of plasticity and phenotype-environment mismatch costs, which may potentially underlie the phenomenon of ontogenetic windows constraining plasticity.

Stream ecosystems harbor many secretive and imperiled species, and studies of vertebrates in these systems face the challenges of relatively low detection rates and high costs. Environmental DNA (eDNA) has recently been confirmed as a sensitive and efficient tool for documenting aquatic vertebrates in wetlands and in a large river and canal system. However, it was unclear whether this tool could be used to detect low-density vertebrates in fast-moving streams where shed cells may travel rapidly away from their source. To evaluate the potential utility of eDNA techniques in stream systems, we designed targeted primers to amplify a short, species-specific DNA fragment for two secretive stream amphibian species in the northwestern region of the United States (Rocky Mountain tailed frogs, Ascaphus montanus, and Idaho giant salamanders, Dicamptodon aterrimus). We tested three DNA extraction and five PCR protocols to determine whether we could detect eDNA of these species in filtered water samples from five streams with varying densities of these species in central Idaho, USA. We successfully amplified and sequenced the targeted DNA regions for both species from stream water filter samples. We detected Idaho giant salamanders in all samples and Rocky Mountain tailed frogs in four of five streams and found some indication that these species are more difficult to detect using eDNA in early spring than in early fall. While the sensitivity of this method across taxa remains to be determined, the use of eDNA could revolutionize surveys for rare and invasive stream species. With this study, the utility of eDNA techniques for detecting aquatic vertebrates has been demonstrated across the majority of freshwater systems, setting the stage for an innovative transformation in approaches for aquatic research.