Search

Transposable elements (TEs) compose a substantial proportion of the largest eukaryotic genomes. TE diversity has been hypothesized to be negatively correlated with genome size, yet empirical demonstrations of such a relationship in a phylogenetic context are largely lacking. Furthermore, the most abundant type of TEs in genomes varies across groups, and it is not clear if there are patterns of TE activity consistent with genome size among different taxa with large genome sizes. We use low-coverage sequencing of 16 species of Neotropical salamanders, which vary ∼7-fold in genome size, to estimate TE relative abundance and diversity for each species. We also compare the divergence of copies of each TE superfamily to estimate patterns of TE activity in each species. We find a negative relationship between TE diversity and genome size, which is consistent with the hypothesis that either competition among TEs or reduced selection against ectopic recombination may result in lower diversity in the largest genomes. We also find divergent activity patterns in the largest versus the smallest genomes, suggesting that the history of TE activity may explain differences in genome size. Our results suggest that both TE diversity and relative abundance may be predictable, at least within taxonomic groups., Competing Interests: Conflict of Interest: The authors declare that they have no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Background: Gram-negative bacteria are the main bacterial pathogens infecting Chinese giant salamanders (Andrias davidianus; CGS) in captivity and the wild, causing substantial economic losses in the CGS industry. However, the molecular mechanisms underlying pathogenesis following infection remain unclear., Results: Spleen-derived macrophages from healthy CGS were isolated, cultured, and identified using density gradient centrifugation and immunofluorescence. A macrophage transcriptome database was established 0, 6, and 12 h post lipopolysaccharide stimulation using RNA-sequencing. In the final database 76,743 unigenes and 4,698 differentially expressed genes (DEGs) were functionally annotated. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment results showed that DEGs were concentrated in toll-like receptor-nuclear factor kappa B-related immune pathways. Ten DEGs were validated 12 h after lipopolysaccharide (LPS) stimulation. Although the common LPS recognition receptor toll-like receptor 4 was not activated and the key adaptor protein MyD88 showed no significant response, we observed significant up-regulation of the following adaptors: toll/interleukin-1 receptor domain-containing adaptor inducing interferon-β, tumour necrosis factor receptor-associated factor 6, and transforming growth factor-β activated kinase 1, which are located downstream of the non-classical MyD88 pathway., Conclusions: In contrast to that in other species, macrophage activation in CGS could depend on the non-classical MyD88 pathway in response to bacterial infection. Our study provides insights into the molecular mechanisms regulating CGS antibacterial responses, with implications for disease prevention and understanding immune evolution in amphibians., (© 2024. The Author(s).)

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.)

High-density genotyping methods have revolutionized the field of population and conservation genetics in the past decade. To exploit the technological and analytical advances in the field, access to high-quality genetic material is a key component. However, access to such samples in endangered and rare animals is often challenging or even impossible. Here, we used a minimally invasive sampling method (MIS) in the endangered cave salamander Proteus anguinus, the olm, to generate thousands of genetic markers using ddRADseq for population and conservation genomic analyses. Using tail clips and MIS skin swabs taken from the same individual, we investigated genotyping data properties of the two different sampling types. We found that sufficient DNA can be extracted from swab samples to generate up to 200,000 polymorphic SNPs in divergent Proteus lineages. Swab and tissue samples were highly reproducible exhibiting low SNP genotyping error rates. We found that SNPs were most frequently (~50%) located within genic regions, while the rest mapped to mostly flanking regions of repetitive DNA. The vast majority of DNA recovered from swabbing was host DNA. However, a fraction of DNA recovered from swabs contained additional ecological information on the species, including eDNA from the surrounding environment and bacterial skin fauna. Most exogenous DNA recovered from swabs were bacteria (~80%), followed by vertebrates (~20%). Our results demonstrate that MIS can be used to (i) generate tens of thousands of ddRADseq markers for conservation and population genomic analyses and (ii) inform on the species health status and ecology from exogenous DNA., (© 2024 The Author(s). Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Taiwan harbors five endemic species of salamanders (Hynobius spp.) that inhabit distinct alpine regions, contributing to population fragmentation across isolated "sky islands". With an evolutionary history spanning multiple glacial-interglacial cycles, these species represent an exceptional paradigm for exploring biogeography and speciation. However, a lack of suitable genetic markers applicable across species has limited research efforts. Thus, developing cross-amplifying markers is imperative. Expressed sequence-tag simple-sequence repeats (EST-SSRs) that amplify across divergent lineages are ideal for species identification in instances where phenotypic differentiation is challenging. Here, we report a suite of cross-amplifying EST-SSRs from the transcriptomes of the five Hynobius species that exhibit an interspecies transferability rate of 67.67%. To identify individual markers exhibiting cross-species polymorphism and to assess interspecies genetic diversity, we assayed 140 individuals from the five species across 84 sampling sites. A set of EST-SSRs with a high interspecies polymorphic information content (PIC = 0.63) effectively classified these individuals into five distinct clusters, as supported by discriminant analysis of principal components (DAPC), STRUCTURE assignment tests, and Neighbor-joining trees. Moreover, pair-wise F ST values > 0.15 indicate notable between-cluster genetic divergence. Our set of 20 polymorphic EST-SSRs is suitable for assessing population structure within and among Hynobius species, as well as for long-term monitoring of their genetic composition., (© 2024. The Author(s).)

This paper offers an historical introduction to Pere Alberch's evolutionary thought and his contributions to Evo-Devo, based on his unique approach to experimental teratology. We will take as our point of reference the teratogenic experiments developed by Alberch and Emily A. Gale during the 1980s, aimed at producing monstrous variants of frogs and salamanders. We will analyze his interpretation of the results of these experiments within the framework of the emergence of evolutionary developmental biology (or "Evo-Devo"). The aim is understand how Alberch interpreted teratological anomalies as highly revealing objects of study for understanding the development of organic form, not only in an ontogenetic sense-throughout embryonic development-but also phylogenetically-throughout the evolution of species. Alberch's interpretation of monsters reflects the influence of a long tradition of non-Darwinian evolutionary thought, which began in the nineteenth century and was continued in the twentieth century by people such as Richard Goldschmidt, Conrad H. Waddington, and Stephen Jay Gould. They all proposed various non-gradualist models of evolution, in which embryonic development played a central role. Following this tradition, Alberch argued that, in order to attain a correct understanding of the role of embryological development in evolution, it was necessary to renounce the gradualist paradigm associated with the Darwinian interpretation of evolution, which understood nature as a continuum. According to Alberch, the study of monstrous abnormalities was of great value in understanding how certain epigenetic restrictions in development could give rise to discontinuities and directionality in morphological transformations throughout evolution., (© 2024. The Author(s).)

Global climate change is expected to have a profound effect on species distribution. Due to the temperature constraints, some narrow niche species could shift their narrow range to higher altitudes or latitudes. In this study, we explored the correlation between species traits, genetic structure, and geographical range size. More specifically, we analyzed how these variables are affected by differences in fundamental niche breadth or dispersal ability in the members of two sympatrically distributed stream-dwelling amphibian species (frog, Quasipaa yei; salamander, Pachyhynobius shangchengensis), in Dabie Mountains, East China. Both species showed relatively high genetic diversity in most geographical populations and similar genetic diversity patterns (JTX, low; BYM, high) correlation with habitat changes and population demography. Multiple clustering analyses were used to disclose differentiation among the geographical populations of these two amphibian species. Q. yei disclosed the relatively shallow genetic differentiation, while P. shangchengensis showed an opposite pattern. Under different historical climatic conditions, all ecological niche modeling disclosed a larger suitable habitat area for Q. yei than for P. shangchengensis; these results indicated a wider environment tolerance or wider niche width of Q. yei than P. shangchengensis. Our findings suggest that the synergistic effects of environmental niche variation and dispersal ability may help shape genetic structure across geographical topology, particularly for species with extremely narrow distribution., (© 2023 International Society of Zoological Sciences, Institute of Zoology/Chinese Academy of Sciences and John Wiley & Sons Australia, Ltd.)

Chinese giant salamander (Andrias davidianus) is the largest extant urodela species and has unique evolutionary position. Studying the immune system of Chinese giant salamander contributes to understanding the evolution of immune systems of vertebrates. The NLR-related protein 3 (NLRP3) inflammasome comprised of NLRP3, ASC and caspase-1 play important roles in the host innate immunity. However, little is know about the NLRP3 inflammasome components in Chinese giant salamander. In this study, the NLRP3, apoptosis-associated speck-like protein (ASC) and caspase-1 (adaNLRP3, adaASC and adaCaspase-1) were characterized from Chinese giant salamander. The proteins of these three genes shared similar motifs and structures with their mammalian counterparts, with a PYD motif, a nucleotide-binding domain (NACHT) motif, and four leucine-rich repeat domain (LRR) motifs identified in adaNLRP3, a pyrin domain (PYD) motif and a caspase recruitment domain (CARD) motif in adaASC, and a CARD motif and a CASc motif in adaCaspase-1. These three genes were constitutively expressed in the skin, heart, lung, kidney, muscle, brain, spleen, and liver of Chinese giant salamander. Following Aeromonas hydrophia infection, all the three genes were up-regulated in various tissues. Molecular docking analysis revealed that the key residues involved in forming the adaNLRP3/adaASC complex were located in the PYD motifs, and that involved in forming the adaASC/adaCaspase-1 complex were located in the CARD motifs. Further analysis revealed that the hydrogen bonds and salt bridges had crucial roles in the formation of adaNLRP3/acaASC and adaASC/adaCaspase-1 complexes. To the best of our knowledge, this is the first report on the NLRP3 inflammasome components in Chinese giant salamander which will be helpful in further understanding the function of the NLRP3 inflammasome and in elucidating its role in the immune response to microbes., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Genome size variation in eukaryotes has myriad effects on organismal biology from the genomic to whole-organism level. Large genome size may be associated with lower selection efficiency because lower effective population sizes allow fixation of deleterious mutations via genetic drift, increasing genome size and decreasing selection efficiency. Because of a hypothesized negative relationship between genome size and recombination rate per base pair, increased genome size could also increase the effect of linked selection in the genome, decreasing the efficiency with which natural selection can fix or remove mutations. We used a transcriptomic dataset of 15 and a subset of six Neotropical salamander species ranging in genome size from 12 to 87 pg to study the relationship between genome size and efficiency of selection. We estimated dN/dS of salamanders with small and large genomes and tested for relaxation of selection in the larger genomes. Contrary to our expectations, we did not find a significant relationship between genome size and selection efficiency or strong evidence for higher dN/dS values in species with larger genomes for either species set. We also found little evidence for relaxation of selection in species with larger genomes. A positive correlation between genome size and range size (a proxy of population size) in this group disagrees with predictions of stronger drift in species with larger genomes. Our results highlight the complex interactions between the many forces shaping genomic variation in organisms with genomic gigantism., (© 2024. The Author(s).)

Bifenthrin (BF) is a broad-spectrum insecticide that has gained widespread use due to its high effectiveness. However, there is limited research on the potential toxic effects of bifenthrin pollution on amphibians. This study aimed to investigate the 50 % lethal concentration (LC 50 ) and safety concentration of Chinese giant salamanders (CGS) exposed to BF (at 0, 6.25,12.5,25 and 50 μg/L BF) for 96 h. Subsequently, CGS were exposed to BF (at 0, 0.04, and 4 μg/L BF) for one week to investigate its toxic effects. Clinical poisoning symptoms, liver pathology, oxidative stress factors, DNA damage, and transcriptome differences were observed and analyzed. The results indicate that exposure to BF at 4 μg/L significantly decreased the adenosine-triphosphate (ATP), superoxide dismutase (SOD), glutathione (GSH), and catalase (CAT) contents in the brain, liver, and kidney of CGS. Additionally, the study found that the malondialdehyde (MDA), reactive oxygen species (ROS), and 8-hydroxydeoxyguanosine (8-OHdG) contents were increased. The liver tissue exhibited significant inflammatory reactions and structural malformations. RNA-seq analysis of the liver showed that BF caused abnormal antioxidant indices of CGS. This affected molecular function genes such as catalytic activity, ATP-dependent activity, metabolic processes, signaling and immune system processes, behavior, and detoxification, which were significantly upregulated, resulting in the differential genes significantly enriched in the calcium signaling pathway, PPARα signaling pathway and NF-kB signaling pathway. The results suggest that BF induces the abnormal production of free radicals, which overwhelms the body's self-defense system, leading to varying degrees of oxidative stress. This can result in oxidative damage, DNA damage, abnormal lipid metabolism, autoimmune diseases, clinical poisoning symptoms, and tissue inflammation. This work provides a theoretical basis for the rational application of bifenthrin and environmental risk assessment, as well as scientific guidance for the conservation of amphibian populations., Competing Interests: Declaration of competing interest All authors disclosed no relevant relationships., (Copyright © 2024 Elsevier B.V. All rights reserved.)

The Chinese giant salamander (CGS) Andrias davidianus is the largest extant amphibian and has recently become an important species for aquaculture with high economic value. Meanwhile, its wild populations and diversity are in urgent need of protection. Exploring the mechanism of its early gonadal differentiation will contribute to the development of CGS aquaculture and the recovery of its wild population. In this study, transcriptomic and phenotypic research was conducted on the critical time points of early gonadal differentiation of CGS. The results indicate that around 210 days post-hatching (dph) is the critical window for female CGS's gonadal differentiation, while 270 dph is that of male CGS. Besides, the TRPM1 gene may be the crucial gene among many candidates determining the sex of CGS. More importantly, in our study, key genes involved in CGS's gonadal differentiation and development are identified and their potential pathways and regulatory models at early stage are outlined. This is an initial exploration of the molecular mechanisms of CGS's early gonadal differentiation at multiple time points, providing essential theoretical foundations for its captive breeding and offering unique insights into the conservation of genetic diversity in wild populations from the perspective of sex development., Competing Interests: Declaration of competing interest There are no declarations of interest between authors., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Assessing direct fitness effects of individual genetic diversity is challenging due to the intensive and long-term data needed to quantify survival and reproduction in the wild. But resolving these effects is necessary to determine how inbreeding and outbreeding influence eco-evolutionary processes. We used 8 years of capture-recapture data and single nucleotide polymorphism genotypes for 1906 individuals to test for effects of individual heterozygosity on stage-specific survival probabilities in the salamander Gyrinophilus porphyriticus. The life cycle of G. porphyriticus includes an aquatic larval stage followed by metamorphosis into a semi-aquatic adult stage. In our study populations, the larval stage lasts 6-10 years, metamorphosis takes several months, and lifespan can reach 20 years. Previous studies showed that metamorphosis is a sensitive life stage, leading us to predict that fitness effects of individual heterozygosity would occur during metamorphosis. Consistent with this prediction, monthly probability of survival during metamorphosis declined with multi-locus heterozygosity (MLH), from 0.38 at the lowest MLH (0.10) to 0.06 at the highest MLH (0.38), a reduction of 84%. Body condition of larvae also declined significantly with increasing MLH. These relationships were consistent in the three study streams. With evidence of localised inbreeding within streams, these results suggest that outbreeding disrupts adaptations in pre-metamorphic and metamorphic individuals to environmental gradients along streams, adding to evidence that headwater streams are hotspots of microgeographic adaptation. Our results also underscore the importance of incorporating life history in analyses of the fitness effects of individual genetic diversity and suggest that metamorphosis and similar discrete life stage transitions may be critical periods of viability selection., (© 2024 John Wiley & Sons Ltd.)

Excessive antibiotic use has led to the spread of antibiotic resistance genes (ARGs), impacting gut microbiota and host health. However, the effects of antibiotics on amphibian populations remain unclear. We investigated the impact of oxytetracycline (OTC) and ciprofloxacin (CIP) on Chinese giant salamanders (Andrias davidianus), focusing on gut microbiota, ARGs, and gene expression by performing metagenome and transcriptome sequencing. A. davidianus were given OTC (20 or 40 mg/kg) or CIP (50 or 100 mg/kg) orally for 7 days. The results revealed that oral administration of OTC and CIP led to distinct changes in microbial composition and functional potential, with CIP treatment having a greater impact than OTC. Antibiotic treatment also influenced the abundance of ARGs, with an increase in fluoroquinolone and multi-drug resistance genes observed post-treatment. The construction of metagenome-assembled genomes (MAGs) accurately validated that CIP intervention enriched fish-associated potential pathogens Aeromonas hydrophila carrying an increased number of ARGs. Additionally, mobile genetic elements (MGEs), such as phages and plasmids, were implicated in the dissemination of ARGs. Transcriptomic analysis of the gut revealed significant alterations in gene expression, particularly in immune-related pathways, with differential effects observed between OTC and CIP treatments. Integration of metagenomic and transcriptomic data highlighted potential correlations between gut gene expression and microbial composition, suggesting complex interactions between the host gut and its gut microbiota in response to antibiotic exposure. These findings underscore the importance of understanding the impact of antibiotic intervention on the gut microbiome and host health in amphibians, particularly in the context of antibiotic resistance and immune function., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge." Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for > 100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical. [Amphibian; hybridization; introgression; Plethodontidae; stream capture.]., (© The Author(s) 2023. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Applications of environmental DNA (eDNA) analysis methods for biomonitoring have grown exponentially over the last decade and provide a wealth of new information on the distribution of species. However, eDNA methods have limited application for estimating population-level metrics. Environmental RNA (eRNA) has the potential to address ecological questions by gathering population demographic information from environmental media but may be challenging to detect and analyze. We developed gene-specific eRNA assays targeting keratin-associated genes in two focal species, American bullfrogs (Lithobates catesbeianus) and tiger salamanders (Ambystoma mavortium) to answer an important question in amphibian management: whether species detections represent breeding populations versus transitory adults. We performed an extensive laboratory validation with amphibians housed across development stages, where we collected 95 and 127 environmental samples for bullfrogs and salamanders, respectively. Both assays were highly specific to the larval stage and amplified with high sensitivity (90% in bullfrog and 88.4% in tiger salamander samples). We then applied our validated assays to multiple natural systems. When larvae were present, we found 74.1% overall detection in bullfrog field samples and 70.8% and 48.5% overall detection in field samples from ponds with A. macrodactylum and A. californiense larvae, correlating with eDNA detection rates. When only adults were present, we did not detect larvae-specific eRNA in A. macrodactylum ponds, despite high eDNA detection rates. Although much work is ahead for optimizing assay design, sampling and filtering methods, we demonstrate that eRNA can successfully be used to discern life stages with direct application for ecology and conservation management., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

The Chinese giant salamander ( Andrias davidianus ), listed as an endangered species under "secondary protection" in China, faces significant threats due to ecological deterioration and the expansion of human activity. Extensive field investigations are crucial to ascertain the current status in the wild and to implement effective habitat protection measures to safeguard this species and support its population development. Traditional survey methods often fall short due to the elusive nature of the A. davidianus , presenting challenges that are time-consuming and generally ineffective. To overcome these obstacles, this study developed a real-time monitoring method that uses environmental DNA (eDNA) coupled with recombinase polymerase amplification and lateral flow strip (RPA-LFD). We designed five sets of species-specific primers and probes based on mitochondrial genome sequence alignments of A. davidianus and its close relatives. Our results indicated that four of these primer/probe sets accurately identified A. davidianus, distinguishing it from other tested caudata species using both extracted DNA samples and water samples from a tank housing an individual. This method enables the specific detection of A. davidianus genomic DNA at concentrations as low as 0.1 ng/mL within 50 min, without requiring extensive laboratory equipment. Applied in a field survey across four sites in Huangshan City, Anhui Province, where A. davidianus is known to be distributed, the method successfully detected the species at three of the four sites. The development of these primer/probe sets offers a practical tool for field surveying and monitoring, facilitating efforts in population recovery and resource conservation for A. davidianus .

Knowledge of the phylogeographic history of organisms is valuable for understanding their evolutionary processes. To the best of our knowledge, the phylogeographic structure of Hokuriku salamander, Hynobius takedai , an endangered species, remains unclear. This study aimed to elucidate the phylogeographic history of H. takedai , which is expected to be strongly influenced by paleogeographic events. Phylogenetic analysis based on partial sequences of the mitochondrial DNA cytochrome b gene confirmed the genetic independence of H. takedai , and the divergence time with closely related species was estimated to be from the Late Pliocene to the Early Pleistocene. In the phylogenetic tree, two clades were identified within H. takedai , and their haplotypes were found in samples collected from the west and east of the distribution range. These intraspecific divergences were strongly influenced by geohistorical subdivisions of the current major distribution areas in the Middle Pleistocene. One clade was further subdivided and its formation may have been influenced by sea level changes in the Late Pleistocene.

In comparison to biodiversity on Earth's surface, subterranean biodiversity has largely remained concealed. The olm (Proteus anguinus) is one of the most enigmatic extant cave inhabitants, and until now little was known regarding its genetic structure and evolutionary history. Olms inhabit subterranean waters throughout the Dinaric Karst of the western Balkans, with a seemingly uniform phenotypic appearance of cave-specialized traits: an elongate body, snout and limbs, degenerated eyes and loss of pigmentation ("white olm"). Only a single small region in southeastern Slovenia harbours olms with a phenotype typical of surface animals: pigmented skin, eyes, a blunt snout and short limbs ("black olm"). We used a combination of mitochondrial DNA and genome-wide single nucleotide polymorphism data to investigate the molecular diversity, evolutionary history and biogeography of olms along the Dinaric Karst. We found nine deeply divergent species-level lineages that separated between 17 and 4 million years ago, while molecular diversity within lineages was low. We detected no signal of recent admixture between lineages and only limited historical gene flow. Biogeographically, the contemporaneous distribution of lineages mostly mirrors hydrologically separated subterranean environments, while the historical separation of olm lineages follows microtectonic and climatic changes in the area. The reconstructed phylogeny suggests at least four independent transitions to the cave phenotype. Two of the species-level lineages have miniscule ranges and may represent Europe's rarest amphibians. Their rarity and the decline in other lineages call for protection of their subterranean habitats., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

When invasive and endangered native taxa hybridize, the resulting admixture introduces novel conservation challenges. Across a large region of central California, a hybrid swarm consisting of admixed endangered California tiger salamanders (CTS) (Ambystoma californiense) and introduced barred tiger salamanders (BTS) (Ambystoma mavortium) has replaced native populations, threatening the genetic integrity of CTS and the vernal pool systems they inhabit. We employed a large-scale, genomically informed field experiment to test whether shortening breeding pond hydroperiod would favor native CTS genotypes. We constructed 14 large, seminatural ponds to evaluate the effect of hydroperiod duration on larval survival and mass at metamorphosis. We tracked changes in non-native allele frequencies with a 5237-gene exon capture array and employed a combination of custom Bayesian and generalized linear models to quantify the effect of pond duration on salamander fitness. Earlier work on this system showed hybrid superiority under many conditions and suggested that hybrids are favored in human-modified ponds with artificially long hydroperiods. Consistent with these earlier studies, we found overwhelming evidence for hybrid superiority. Very short hydroperiods substantially reduced the mass (1.1-1.5 fold) and survival probability (10-13 fold) of both native and hybrid larvae, confirming that hydroperiod likely exerts a strong selective pressure in the wild. We identified 86 genes, representing 1.8% of 4723 screened loci, that significantly responded to this hydroperiod-driven selection. In contrast to earlier work, under our more natural experimental conditions, native CTS survival and size at metamorphosis were always less than hybrids, suggesting that hydroperiod management alone will not shift selection to favor native larval genotypes. However, shortening pond hydroperiod may limit productivity of hybrid ponds, complementing other strategies to remove hybrids while maintaining vernal pool ecosystems. This study confirms and expands on previous work that highlights the importance of hydroperiod management to control invasive aquatic species., (© 2023 The Authors. Conservation Biology published by Wiley Periodicals LLC on behalf of Society for Conservation Biology.)

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 South China giant salamander, Andrias sligoi, is one of the largest extant amphibian species worldwide. It was recently distinguished from another Chinese species, the Chinese giant salamander, Andrias davidianus, which is considered Critically Endangered according to the International Union for Conservation of Nature (IUCN) Red List. It appears too late to save this extremely rare and large amphibian in situ. Another extant species of the same genus, Andrias japonicus, inhabits Japan. However, the introduction of Chinese giant salamanders into some areas of Japan has resulted in hybridization between the Japanese and Chinese species. During our genetic screening of giant salamanders in Japan, we unexpectedly discovered four individuals of the South China giant salamander: two were adult males in captivity, and one had recently died. The last individual was a preserved specimen. In this study, we report these extremely rare individuals of A. sligoi in Japan and discuss the taxonomic and conservational implications of these introduced individuals., (© 2024. The Author(s).)

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species., (© 2023 John Wiley & Sons Ltd.)

For a long time, it has been debated whether the two giant salamanders, Andrias japonicus from Japan and A. davidianus from China, are conspecific or heterospecific. Morphological information about their diagnostic characteristics has been limited, without considering sexual dimorphism and/or body size variation. Recently, A. davidianus, which was introduced into Japan sometime in the past, has been found to hybridize with A. japonicus in situ. Taxonomic identification of individuals involved in this unusual breeding is made based on mitochondrial and nuclear DNA analyses. This identification method is time-consuming and costly. Thus, developing easier methods of identification, such as utilizing external morphological characteristics, is urgently needed. In this study, we verify previous descriptions showing that A. davidianus has a longer relative tail length than A. japonicus, and the tubercles on the lower jaw and throat were present in both sexes of A. davidianus. In addition, many head characteristics were found to be relatively larger in A. davidianus than in A. japonicus, which were new distinguishing characters. These morphological differences help support the idea that these are heterospecific lineages. In hybrids, relative values of head width and tail length were larger than those of A. japonicus, and the tubercles on their lower jaw and throat were present as in A. davidianus, suggesting that the hybrids and A. davidianus are distinguishable from A. japonicus.

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.)

Thioredoxin-like protein-1 (TXNL1) is the member of thioredoxin superfamily, a family of thiol oxidoreductases. TXNL1 plays an important role in scavenging ROS and the maintenance of cellular redox balance. However, its physiological functions in Andrias davidianus have not been well understood. In the present study, the full-length cDNA encoding thioredoxin-like protein-1 (AdTXNL1) of A. davidianus was cloned, the mRNA tissue distribution was analyzed, and the function was characterized. The Adtxnl1 cDNA contained an open reading frame (ORF) of 870 bp encoding a polypeptide of 289 amino acids with the N-terminal TRX domain, a Cys 34 -Ala 35 -Pro 36 -Cys 37 (CAPC) motif, and the C-terminal proteasome-interacting thioredoxin domain (PITH). The mRNA of AdTXNL1 was expressed in a wide range of tissues, with the highest level in the liver. The transcript level of AdTXNL1 was significantly up-regulated post Aeromonas hydrophila challenge in liver tissue. Moreover, the recombinant AdTXNL1 protein was produced and purified, and used to investigate the antioxidant activity. In the insulin disulfide reduction assay, rAdTXNL1 exhibited strong antioxidant capability. Altogether, the thioredoxin-like protein-1 may be involved in reduction/oxidation (redox) balance and as an important immunological gene in A. davidianus., Competing Interests: Declaration of competing interest The authors declare that they have no conflict of interest., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

The key role of symbiotic skin bacteria communities in amphibian resistance to emerging pathogens is well recognized, but factors leading to their dysbiosis are not fully understood. In particular, the potential effects of population translocations on the composition and diversity of hosts' skin microbiota have received little attention, although such transfers are widely carried out as a strategy for amphibian conservation. To characterize the potential reorganization of the microbiota over such a sudden environmental change, we conducted a common-garden experiment simulating reciprocal translocations of yellow-spotted salamander larvae across three lakes. We sequenced skin microbiota samples collected before and 15 days after the transfer. Using a database of antifungal isolates, we identified symbionts with known function against the pathogen Batrachochytrium dendrobatidis, a major driver of amphibian declines. Our results indicate an important reorganization of bacterial assemblages throughout ontogeny, with strong changes in composition, diversity and structure of the skin microbiota in both control and translocated individuals over the 15 days of monitoring. Unexpectedly, the diversity and community structure of the microbiota were not significantly affected by the translocation event, thus suggesting a strong resilience of skin bacterial communities to environmental change-at least across the time-window studied here. A few phylotypes were more abundant in the microbiota of translocated larvae, but no differences were found among pathogen-inhibiting symbionts. Taken together, our results support amphibian translocations as a promising strategy for this endangered animal class, with limited impact on their skin microbiota., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Significant advances have been made in species delimitation and numerous methods can test precisely defined models of speciation, though the synthesis of phylogeography and taxonomy is still sometimes incomplete. Emerging consensus treats distinct genealogical clusters in genome-scale data as strong initial evidence of speciation in most cases, a hypothesis that must therefore be falsified under an explicit evolutionary model. We can now test speciation hypotheses linking trait differentiation to specific mechanisms of divergence with increasingly large data sets. Integrative taxonomy can, therefore, reflect an understanding of how each axis of variation relates to underlying speciation processes, with nomenclature for distinct evolutionary lineages. We illustrate this approach here with Seal Salamanders (Desmognathus monticola) and introduce a new unsupervised machine-learning approach for species delimitation. Plethodontid salamanders are renowned for their morphological conservatism despite extensive phylogeographic divergence. We discover 2 geographic genetic clusters, for which demographic and spatial models of ecology and gene flow provide robust support for ecogeographic speciation despite limited phenotypic divergence. These data are integrated under evolutionary mechanisms (e.g., spatially localized gene flow with reduced migration) and reflected in emergent properties expected under models of reinforcement (e.g., ethological isolation and selection against hybrids). Their genetic divergence is prima facie evidence for species-level distinctiveness, supported by speciation models and divergence along axes such as behavior, geography, and climate that suggest an ecological basis with subsequent reinforcement through prezygotic isolation. As data sets grow more comprehensive, species-delimitation models can be tested, rejected, or corroborated as explicit speciation hypotheses, providing for reciprocal illumination of evolutionary processes and integrative taxonomies. [Desmognathus; integrative taxonomy; machine learning; species delimitation.]., (© The Author(s) 2022. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Hynobius akiensis sensu lato has recently been split into three species based on short sequence analyses of cyt-b gene of mtDNA and without data of nuclear DNA, and strange sympatric distribution in some areas has been indicated in two species. We analyzed nuclear DNA marker (SNPs) and complete sequence of cyt-b in H. akiensis sensu lato to reassess species delimitation and genetic introgression among species. As a result, we found two lineages with discordant mitochondrial and nuclear DNA in some areas. Of H. akiensis sensu lato, each of the two contains the type locality of two species recently reported (H. sumidai and H. geiyoensis), and the use of these names has been previously advocated. However, their sympatric distribution was rejected based on nuclear DNA data, which we consider is more reliable than mtDNA. We thus clarify geographic boundary of these two species and revise the species delimitations.

Loss of biodiversity is a serious concern, and amphibians are particularly threatened. Most small salamanders in Japan are endangered. Distributional information is fundamental to the conservation of these rare species; however, small salamanders are generally difficult to locate or catch. Environmental DNA analysis is an effective survey method for monitoring such rare species. The conventional polymerase chain reaction (PCR) method, which combines PCR amplification with subsequent electrophoresis, and the real-time PCR method, which uses fluorescent material, are commonly used for this purpose. In this study, a comparison of these two detection methods was conducted using a rare salamander species, Hynobius boulengeri, as a model case. We compared three points: (i) detection sensitivity, (ii) influence of environmental factors related to detection, and (iii) time and financial costs of the two methods. To perform this comparison, we developed a real-time PCR detection assay, conducted field surveys, and compared the time and financial costs of conventional and real-time PCR methods. The comparison showed no statistical difference in the detection sensitivity from field samples, and the effects of environmental factors tended to be similar. In addition, the financial cost was lower for the conventional PCR method while the time cost was lower for the real-time PCR method. Therefore, selecting eDNA detection methods based on objectives, time, and financial costs will promote efficient monitoring and contribute to the conservation of rare species., (© 2023. The Author(s), under exclusive licence to The Japan Society for Analytical Chemistry.)

How does amphibian MHC diversity fit in the landscape of vertebrate evolution? Mimnias et al. (2022) address this gap in the field of MHC evolution by focusing on the lesser described MHC class I in salamanders. These findings contribute to understanding MHC diversity and the susceptibility of amphibians to pathogens, which could lead to future research on a major threat to amphibian biodiversity, chytrid fungi., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Building a blastema from the stump is a key step of salamander limb regeneration. Stump-derived cells temporarily suspend their identity as they contribute to the blastema by a process generally referred to as dedifferentiation. Here, we provide evidence for a mechanism that involves an active inhibition of protein synthesis during blastema formation and growth. Relieving this inhibition results in a higher number of cycling cells and enhances the pace of limb regeneration. By small RNA profiling and fate mapping of skeletal muscle progeny as a cellular model for dedifferentiation, we find that the downregulation of miR-10b-5p is critical for rebooting the translation machinery. miR-10b-5p targets ribosomal mRNAs, and its artificial upregulation causes decreased blastema cell proliferation, reduction in transcripts that encode ribosomal subunits, diminished nascent protein synthesis, and retardation of limb regeneration. Taken together, our data identify a link between miRNA regulation, ribosome biogenesis, and protein synthesis during newt limb regeneration., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Ranavirus is a large nucleocytoplasmic DNA virus. Chinese giant salamander iridovirus (CGSIV) belongs to the ranavirus genus, and its replication involves a series of essential viral genes. Viral PCNA is a gene closely associated with viral replication. CGSIV-025L also encodes PCNA-like genes. We have described the function of CGSIV-025L in virus replication. The promoter of CGSIV-025L is activated during viral infection, and it is an early (E) gene that can be effectively transcribed after viral infection. CGSIV-025L overexpression promoted viral replication and viral DNA replication. siRNA interfered with CGSIV-025L expression and attenuated viral replication and viral DNA replication. The Δ025L-CGSIV strain with the deletion of CGSIV-025L could not replicate normally and could be rescued by the replenishment of 025L. CGSIV-025L was proven to be an essential gene for CGSIV by overexpression, interference, and deletion mutation experiments. CGSIV-025L was found to interact with CGSIV-062L by yeast two-hybrid, CoIP, and GST pulldown. Thus, the current study demonstrated that CGSIV-025L is an essential gene of CGSIV, which may be involved in viral infection by participating in viral DNA replication and interacting with replication-related proteins.

Whole-mount in situ hybridization (WISH) is widely used to visualize transcribed gene sequences (mRNA) in developing embryos, larvae, and other nucleotide probe permeable tissue samples. This methodology involves the hybridization of an antisense nucleotide probe to the target mRNA, followed by chromogen or fluorescence-based detection. Here we describe a protocol for the spatiotemporal analysis of mRNA transcripts in axolotl embryos/larvae using digoxigenin-labeled riboprobes, anti-digoxigenin alkaline phosphatase, Fab fragments antibody, and NBT/BCIP chromogen detection., (© 2023. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Polystomatidae is a monogenean family whose representatives infect mainly (semi)-aquatic tetrapods. Species of Sphyranura Wright, 1879 exhibit ectoparasitism on salamander hosts, with molecular work supporting their inclusion within Polystomatidae, at an early diverging, yet unresolved, position in the clade of otherwise endoparasitic polystomatid parasites of batrachian hosts. Records of representatives of Sphyranura are scarce with genetic data only available for S. oligorchis Alvey, 1933. Based on detailed morphological examination and comparison with type material, we identified worms belonging to Sphyranura infecting Oklahoma salamander (Eurycea tynerensis) as S. euryceae Hughes & Moore, 1943. Along with an amended diagnosis of Sphyranura, we provide the first molecular data for S. euryceae in the form of a mitochondrial genome and nuclear (18S, 28S rRNA) markers. Close morphological similarity between the two species of Sphyranura is reflected in low genetic divergence. Mitochondrial level comparison reveals instances of tRNA gene rearrangements in polystomatids. Although the phylogenetic reconstruction supports Sphyranura as early branching in the lineage of polystomatid monogeneans infecting batrachians, certain nodes remain unresolved., (© S.J. Leeming et al., published by EDP Sciences, 2023.)

The DNAJ family, a class of chaperone proteins involved in protein folding, assembly, and transport, plays an essential role in viral infections. However, the role of DNAJA4 (DnaJ Heat Shock Protein Family (Hsp40) Member A4) in the ranavirus infection has not been reported. This study demonstrates the function of the epithelial papilloma of carp (EPC) DNAJA4 in Chinese giant salamander (Andrias davidianus) iridovirus (CGSIV) replication. DNAJA4 consists of 1479 base pairs and encodes a 492 amino acid polypeptide. Sequence analysis has shown that EPC DNAJA4 contains a conserved J domain and shares 84% homology with Danio rerio DNAJA4 and 68% homology with Homo sapiens DNAJA4. EPC DNAJA4 was localized in the cytoplasm, and its expression was significantly upregulated after CGSIV infection. Overexpression of EPC DNAJA4 promotes CGSIV replication and CGSIV DNA replication. siRNA knockdown of DNAJA4 expression attenuates CGSIV replication and viral DNA replication. Overexpression and interference experiments have proved that EPC DNAJA4 is a pro-viral factor. Co-IP, GST-pulldown, and immunofluorescence confirmed the interaction between EPC DNAJA4 and CGSIV proliferating cell nuclear antigen (PCNA). Our results demonstrate for the first time that EPC DNAJA4 is involved in viral infection by promoting viral DNA replication and interacting with proteins associated with viral replication.

AbstractAvoiding inbreeding is considered a key driver of dispersal evolution, and dispersal distances should be especially important in mediating inbreeding risk because the likelihood of mating with relatives decreases with dispersal distance. However, a lack of direct data on dispersal distances has limited empirical tests of this prediction, particularly in the context of the multiple selective forces that can influence dispersal. Using the headwater stream salamander Gyrinophilus porphyriticus , we tested whether spatial variation in environmental conditions leads to differences in dispersal distances, resulting in spatial variation in the effect of dispersal on inbreeding risk. Using capture-recapture and population genomic data from five streams, we found that dispersal distances were greater in downstream reaches than upstream reaches. Inbreeding risk trended lower for dispersers than nondispersers in downstream reaches but not in upstream reaches. Furthermore, stream reaches did not differ in spatial patterns of individual relatedness, indicating that variation in inbreeding risk was in fact due to differences in dispersal distances. These results demonstrate that environmentally associated variation in dispersal distances can cause the inbreeding consequences of dispersal to vary at fine spatial scales. They also show that selective pressures other than inbreeding avoidance maintain phenotypic variation in dispersal, underscoring the importance of addressing alternative hypotheses in dispersal research.

Amphibians have undergone important evolutionary transitions in reproductive modes and life-cycles. We compare large-scale macroevolutionary patterns in these transitions across the three major amphibian clades: frogs, salamanders, and caecilians. We analyse matching reproductive and phylogenetic data for 4025 species. We find that having aquatic larvae is ancestral for all three groups and is retained by many extant species (33-44%). The most frequent transitions in each group are to relatively uncommon states: live-bearing in caecilians, paedomorphosis in salamanders, and semi-terrestriality in frogs. All three groups show transitions to more terrestrial reproductive modes, but only in caecilians have these evolved sequentially from most-to-least aquatic. Diversification rates are largely independent of reproductive modes. However, in salamanders direct development accelerates diversification whereas paedomorphosis decreases it. Overall, we find a widespread retention of ancestral modes, decoupling of trait transition rates from patterns of species richness, and the general independence of reproductive modes and diversification., (© 2022. The Author(s).)

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution., (© 2022 The Authors. Evolution © 2022 The Society for the Study of Evolution.)

Transposable elements (TEs) are repetitive sequences of DNA that replicate and proliferate throughout genomes. Taken together, all the TEs in a genome form a diverse community of sequences, which can be studied to draw conclusions about genome evolution. TE diversity can be measured using models for ecological community diversity that consider species richness and evenness. Several models predict TE diversity decreasing as genomes expand because of selection against ectopic recombination and/or competition among TEs to garner host replicative machinery and evade host silencing mechanisms. Salamanders have some of the largest vertebrate genomes and highest TE loads. Salamanders of the genus Plethodon, in particular, have genomes that range in size from 20 to 70 Gb. Here, we use Oxford Nanopore sequencing to generate low-coverage genomic sequences for four species of Plethodon that encompass two independent genome expansion events, one in the eastern clade (Plethodon cinereus, 29.3 Gb vs. Plethodon glutinosus, 38.9 Gb) and one in the western clade (Plethodon vehiculum, 46.4 Gb vs Plethodon idahoensis, 67.0 Gb). We classified the TEs in these genomes and found > 40 TE superfamilies, accounting for 22-27% of the genomes. We calculated Simpson's and Shannon's diversity indices to quantify overall TE diversity. In both pairwise comparisons, the diversity index values for the smaller and larger genome were almost identical. This result indicates that, when genomes reach extremely large sizes, they maintain high levels of TE diversity at the superfamily level, in contrast to predictions made by previous studies on smaller genomes., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

The "mountain" dusky is a charismatic ecomorph of Dusky Salamander (Desmognathus): smaller (~20-50mm SVL), more terrestrial and montane woodland species, typically with round tails and extremely variable, often colorful patterns. Originally considered a single species (D. ochrophaeus), populations across the Appalachians are now divided into D. abditus, D. carolinensis, D. ochrophaeus, D. ocoee, and D. orestes, which do not form a clade but are instead scattered across the phylogeny of Desmognathus. A Coastal Plain species with a mountain dusky phenotype (D. apalachicolae) is also nested within a lineage related to D. ocoee. Within the current concept of D. ocoee, there are up to 8 genetically distinct geographic lineages which display significantly different but substantially overlapping morphologies. Recent studies have suggested that as many as six species could reasonably be delimited based on genetic data. We evaluate the strength of that evidence here and find support for five species in the D. apalachicolae + D. ocoee group. The first is D. apalachicolae in essentially its originally described form from the Coastal Plain. The second and third are D. adatsihi sp. nov. and D. balsameus sp. nov., endemic to the Great Smoky/Plott Balsam and Great Balsam Mountains, respectively. The fourth corresponds to the resurrected D. perlapsus for populations from the Alarka and Cowee Mountains in North Carolina, through the Chattahoochee River drainage to the Fall Line in the Piedmont of Georgia and Alabama. The fifth is D. ocoee, here restricted to populations in the Nantahala and Unicoi Mountains and southernmost Blue Ridge escarpment, along with previously recognized populations from the Cumberland Escarpment and Cumberland Plateau of south-central Tennessee and northeastern Alabama. This taxon concept of D. ocoee includes three deeply divergent geographic genetic segments that hybridize across two contact zones, while D. apalachicolae may be nested within some southern Blue Ridge populations of D. ocoee. The resurrected form of D. perlapsus also exhibits some admixture with the newly restricted D. ocoee. While this taxonomic arrangement is robust and stably derived from mitochondrial, morphological, and nuclear data, more sophisticated future analyses sampling additional populations and loci to estimate relationships may reach more subtle conclusions regarding the identity of some lineages within D. ocoee.

Most swamp-dwelling dusky salamanders of the genus Desmognathus from the Coastal Plain were long treated as a single species (Desmognathus auriculatus) ranging from east Texas to southeastern Virginia. This taxon concept was based on the name Salamandra auriculata Holbrook, 1838 with type locality Riceboro, Liberty County, Georgia and a type series that could not be located by later authors. Recent workers have been unable to locate or verify swamp-dwelling populations from east Texas and western Louisiana, which appear to be extirpated and may not have represented a distinct taxon from co-occurring lineages of D. conanti. Recent molecular phylogenies have supported at least four distinct species-level taxa within D. auriculatus. Populations from the Gulf Coastal Plain in eastern Louisiana, Mississippi, and southwestern Alabama were recently described as D. valentinei Means, Lamb, and Bernardo, 2017 and D. pascagoula Pyron, O'Connell, Lamb, and Beamer, 2022. This leaves two remaining species-level lineages with uncertain taxonomy and nomenclature: D. auriculatus A (Alabama, Florida, and Georgia), and D. auriculatus B/C (Georgia, South Carolina, and North Carolina), both of which occur near the type locality. We recently located a specimen at the Muséum national d'Histoire naturelle in Paris (MNHN 0.4675) that we concluded is one of Holbrook's syntypes and designated it as the lectotype, but without allocation. Here, we use linear morphometrics to confidently allocate it to D. auriculatus A, bolstered by examination of three historical topotypic collections. This requires a new name for D. auriculatus B/C, which we describe as D. valtos sp. nov. (suggested common name: Carolina Swamp Dusky Salamander) from Otter Creek, Craven County, North Carolina. Other related and sympatric species of Desmognathus remain to be described from the Atlantic Coastal Plain and adjacent Piedmont of the southeastern United States.

Background: Hynobius hirosei is a lotic-breeding salamander endemic to Shikoku Island in western Japan. Significant allozymic and morphological differences have been found among the populations of this species; however, the degree and pattern of intraspecific variation have not been surveyed using a sufficient number of samples., Methods: For the taxonomic revision of H. hirosei , we conducted genetic and morphological surveys using samples collected throughout the distribution. Phylogenetic analysis using the cytochrome b region of mitochondrial DNA and population structure analysis using single nucleotide polymorphisms were conducted to evaluate the population structure within the species and the degree of genetic differentiation. Subsequently, a morphological survey based on multivariate and univariate analyses was performed to assess the morphological variation., Results: Genetic analyses revealed three genetic groups (Tsurugi, Central, and Nanyo) within H. hirosei , with the Nanyo group distributed allopatrically from the others, and the Tsurugi and Central groups distributed parapatrically with the formation of a hybrid zone between them. The Nanyo group was morphologically distinguishable from the remaining samples, including the topotype of H. hirosei , based on a smaller body size and several ratio values of characters to snout-vent length, longer axilla-groin distance, shorter tail length, shorter internarial distance, longer upper eyelid length, and larger medial tail width. These results support the notion that the Nanyo group is an undescribed species. However, the remaining genetically differentiated groups could not be divided in the present study. Herein, we described the Nanyo group as a new species., Competing Interests: The authors declare that they have no competing interests., (© 2022 Kanamori et al.)

One or more members of four living amphibian clades have independently dispensed with pulmonary respiration and lack lungs, but little is known of the developmental basis of lung loss in any taxon. We use morphological, molecular, and experimental approaches to examine the Plethodontidae, a dominant family of salamanders, all of which are lungless as adults. We confirm an early anecdotal report that plethodontids complete early stages of lung morphogenesis: Transient embryonic lung primordia form but regress by apoptosis before hatching. Initiation of pulmonary development coincides with expression of the lung-specification gene Wnt2b in adjacent mesoderm, and the lung rudiment expresses pulmonary markers Nkx2.1 and Sox9 . Lung developmental-genetic pathways are at least partially conserved despite the absence of functional adult lungs for at least 25 and possibly exceeding 60 million years. Adult lung loss appears associated with altered expression of signaling molecules that mediate later stages of tracheal and pulmonary development.

Metamorphosis is a critical stage in the adaptive development of amphibians from aquatic to terrestrial animals. Metamorphosis of the Chinese giant salamander is mainly manifested by the loss of external gills with consequent changes in the respiratory pattern. The loss of the external gill is regulated by the pathway of apoptosis in which caspase genes are the key factors. This study cloned and expressed the caspase 3/7/8/9 genes of the Chinese giant salamander. The main results were as follows: the complete open reading frames (ORFs) were 885 bp, 960 bp, 1461 bp and 1279 bp, respectively; caspase 3/7/8/9 genes all contained the CASc domain, and most of the motifs were located in CASc domain; and caspase 8 possessed two DED structural domains and caspase 9 possessed a CARD structural domain. Furthermore, results from the tissue distribution analysis indicated that caspase 3/7/8/9 genes were all significantly expressed in the external gill, and at 9 and 10 months of age (MOA), which is the peak time for the loss, the EXPRESSION level of caspase 3/7/8/9 genes was obviously high, which was consistent with the histological result. Moreover, the loss of external gills of the Chinese giant salamander may result from activation of both the apoptosis-related death receptor pathway and the mitochondrial pathway. Finally, it was discovered that thyroid hormone (TH) treatment could both advance the time point at which the external gills of the Chinese giant salamander began to degenerate and shorten this process. Interestingly, at the peak of its metamorphosis (9 MOA), the Chinese giant salamander further accelerated the metamorphosis rate of TH treatment, which suggested a promotive effect on the loss of external gills via the superimposition of the exogenous TH and caspase genes. The study of caspase genes in this experiment was conducive to understanding the mechanism of external gill loss in the Chinese giant salamander, as well as improving our understanding of the metamorphosis development of some Caudata species.

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.

Effective conservation of threatened biota relies on accurate assessments and scientific guidance. As an unfortunate example, Chinese giant salamanders ( Andrias , CGS) remain critically endangered in nature. Misguided conservation efforts, e.g., commercial propagation and releasing of millions of likely non-indigenous or interspecific hybrids, have further compromised conservation initiatives. Limited information on wild populations of CGS poses a significant conservation challenge. Following 18-month long field monitoring, we now report the discovery of a wild population of CGS in a closed nature reserve in Jiangxi Province, China. Genomic assessments reveal its genetic distinctiveness and do not detect genetic admixture with other species. Based on morphological and molecular evidences, we describe this CGS as a new species Andrias jiangxiensis sp. nov. This is the only known species of CGS today with a genetically pure, reproducing, in situ population. This discovery emphasizes the important role that closed nature reserves play in protecting species, and the necessity of integrating long-term field monitoring and genetic assessments. It sets a new pathway for discovering and conserving endangered species, especially for those biotas that are similarly being extirpated by anthropogenic translocations and overexploitation.

We describe a new, narrowly endemic species of swamp-dwelling dusky salamander (Plethodontidae: Desmognathus pascagoula sp. nov.) from the Gulf Coastal Plain of southeastern Mississippi and southwestern Alabama based on linear morphometrics, mitochondrial DNA, and single nucleotide polymorphisms from 881 loci produced using genotype-by-sequencing. Some populations of the new species were historically referred to as D. auriculatus, a polyphyletic assemblage of at least three species in the Atlantic and Gulf Coastal Plain from Texas to North Carolina. Populations of D. auriculatus from the Gulf Coastal Plain in Louisiana and Mississippi were recently described as D. valentinei. The new species includes populations that were tentatively referred to D. valentinei, but we find it is morphologically, genetically, and geographically distinct. It is smaller, has a more defined dorsal color pattern, more irregular whitish portholes in up to three rows on the lateral surfaces of the body and tail, and a brighter orange or yellowish orange postocular stripe. At present, the new species is known from only six extant populations in the lower Pascagoula, Escatawpa, and Mobile drainages. The latter represents a distinct phylogeographic lineage. We also refer a historical collection from the northeastern side of the Mobile-Tensaw River Delta to this species, suggesting a much broader range in the past. We suspect that more populations remain to be discovered in the area, and their potential species-level distinctiveness should be tested further. This discovery increases knowledge of the biodiversity in the southeastern United States Coastal Plain, a candidate region meeting the global criteria for a biodiversity hotspot, and underscores the amount of cryptic diversity likely remaining to be discovered and described in Nearctic salamanders.

Transposable elements (TEs) are sequences that replicate and move throughout genomes, and they can be silenced through methylation of cytosines at CpG dinucleotides. TE abundance contributes to genome size, but TE silencing variation across genomes of different sizes remains underexplored. Salamanders include most of the largest C-values - 9 to 120 Gb. We measured CpG methylation levels in salamanders with genomes ranging from 2N = ∼58 Gb to 4N = ∼116 Gb. We compared these levels to results from endo- and ectothermic vertebrates with more typical genomes. Salamander methylation levels are approximately 90%, higher than all endotherms. However, salamander methylation does not differ from other ectotherms, despite an approximately 100-fold difference in nuclear DNA content. Because methylation affects the nucleotide compositional landscape through 5-methylcytosine deamination to thymine, we quantified salamander CpG dinucleotide levels and compared them to other vertebrates. Salamanders and other ectotherms have comparable CpG levels, and ectotherm levels are higher than endotherms. These data show no shift in global methylation at the base of salamanders, despite a dramatic increase in TE load and genome size. This result is reconcilable with previous studies that considered endothermy and ectothermy, which may be more important drivers of methylation in vertebrates than genome size., (© 2022 The Authors. Evolution © 2022 The Society for the Study of Evolution.)