Search

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomic resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomic resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, anti-predator strategies, and resilience and adaptive responses. They also serve as essential models for studying broad genomic traits, such as evolutionary genome expansions and contractions, as they exhibit the widest range of genome sizes among all animal taxa and possess multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The emergence of long-read sequencing technologies, combined with advanced molecular and computational techniques that improve scaffolding and reduce computational workloads, is now making it possible to address some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC, https://mvs.unimelb.edu.au/amphibian-genomics-consortium ) in early 2023. This burgeoning community already has more than 282 members from 41 countries. The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and call on the research and conservation communities to unite as part of the AGC to enable amphibian genomics research to "leap" to the next level., (© 2024. The Author(s).)

More efficient methods for extensive biodiversity monitoring are required to support rapid measures to address the biodiversity crisis. While environmental DNA (eDNA) metabarcoding and quantitative PCR (qPCR) methods offer advantages over traditional monitoring approaches, their large-scale application is limited by the time and labour required for developing assays and/or for analysis. CRISPR (clustered regularly interspaced short palindromic repeats) diagnostic technologies (Dx) may overcome some of these limitations, but they have been used solely with species-specific primers, restricting their versatility for biodiversity monitoring. Here, we demonstrate the feasibility of designing species-specific CRISPR-Dx assays in silico within a short metabarcoding fragment using a general primer set, a methodology we term 'ampliscanning', for 18 of the 22 amphibian species in Switzerland. We sub-selected nine species, including three classified as regionally endangered, to test the methodology using eDNA sampled from ponds at nine sites. We compared the ampliscanning detections to data from traditional monitoring at these sites. Ampliscanning was successful at detecting target species with different prevalences across the landscape. With only one visit, we detected more species per site than three traditional monitoring visits (visual and acoustic detections by trained experts), in particular more elusive species and previously undocumented but expected populations. Ampliscanning detected 25 species/site combinations compared to 12 with traditional monitoring. Sensitivity analyses showed that larger numbers of field visits and PCR replicates are more important for reliable detection than many technical replicates at the CRISPR-Dx assay level. Given the reduced sampling and analysis effort, our results highlight the benefits of eDNA and CRISPR-Dx combined with universal primers for large-scale monitoring of multiple endangered species across landscapes to inform conservation measures., (© 2024 The Author(s). Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Green algae usually assigned to the genus Oophila are known to colonize egg capsules of amphibian egg masses across the Nearctic and Palearctic regions. We study the phylogenetic relationships of these algae using a phylotranscriptomic data set of 76 protein-coding single-copy nuclear genes. Our data set includes novel RNAseq data for six amphibian-associated and five free-living green algae, and draft genomes of two of the latter. Within the Oophila clade (nested within Moewusinia), we find samples from two European frogs (Rana dalmatina and R. temporaria) closely related to those of the North American frog R. aurora (Oophila subclade III). An isolate from the North American R. sylvatica (subclade IV) appears to be sister to the Japanese isolate from the salamander Hynobius nigrescens (subclade J1), and subclade I algae from Ambystoma maculatum are sister to all other lineages in the Oophila clade. Two free-living algae (Chlamydomonas nasuta and Cd. pseudogloeogama) are nested within the Oophila clade, and a strain of the type species of Chlorococcum (Cc. infusionum) is related to this assemblage. Our phylotranscriptomic tree suggests that recognition of different species within the Oophila clade ("clade B" of earlier studies) is warranted, and calls for a comprehensive taxonomic revision of Moewusinia., 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 The Author(s). Published by Elsevier Inc. All rights reserved.)

Among tetrapod (terrestrial) vertebrates, amphibians remain more closely tied to an amphibious lifestyle than amniotes, and their visual opsin genes may be adapted to this lifestyle. Previous studies have discussed physiological, morphological, and molecular changes in the evolution of amphibian vision. We predicted the locations of the visual opsin genes, their neighboring genes, and the tuning sites of the visual opsins, in 39 amphibian genomes. We found that all of the examined genomes lacked the Rh2 gene. The caecilian genomes have further lost the SWS1 and SWS2 genes; only the Rh1 and LWS genes were retained. The loss of the SWS1 and SWS2 genes in caecilians may be correlated with their cryptic lifestyles. The opsin gene syntenies were predicted to be highly similar to those of other bony vertebrates. Moreover, dual syntenies were identified in allotetraploid Xenopus laevis and X. borealis. Tuning site analysis showed that only some Caudata species might have UV vision. In addition, the S164A that occurred several times in LWS evolution might either functionally compensate for the Rh2 gene loss or fine-tuning visual adaptation. Our study provides the first genomic evidence for a caecilian LWS gene and a genomic viewpoint of visual opsin genes by reviewing the gains and losses of visual opsin genes, the rearrangement of syntenies, and the alteration of spectral tuning in the course of amphibians' evolution., 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 Elsevier Ltd. All rights reserved.)

Unusually for oceanic islands, the granitic Seychelles host multiple lineages of endemic amphibians. This includes an ancient (likely ca. 60 million years) radiation of eight caecilian species, most of which occur on multiple islands.These caecilians have a complicated taxonomic history and their phylogenetic inter-species relationships have been difficult to resolve. Double-digest RAD sequencing (ddRADseq) has been applied extensively to phylogeography and increasingly to phylogenetics but its utility for resolving ancient divergences is less well established. To address this, we applied ddRADseq to generate a genome-wide SNP panel for phylogenomic analyses of the Seychelles caecilians, whose phylogeny has so far not been satisfactorily resolved with traditional DNA markers. Based on 129,154 SNPs, we resolved deep and shallow splits, with strong support. Our findings demonstrate the capability of genome-wide SNPs for evolutionary inference at multiple taxonomic levels and support the recently proposed synonymy of Grandisonia Taylor, 1968 with Hypogeophis Peters, 1879. We revealed three clades of Hypogeophis (large-, medium- and short-bodied) and identify a single origin of the diminutive, stocky-bodied and pointy-snouted phenotype., 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 The Authors. Published by Elsevier Inc. All rights reserved.)

Natural selection can drive organisms to strikingly similar adaptive solutions, but the underlying molecular mechanisms often remain unknown. Several amphibians have independently evolved highly adhesive skin secretions (glues) that support a highly effective antipredator defence mechanism. Here we demonstrate that the glue of the Madagascan tomato frog, Dyscophus guineti, relies on two interacting proteins: a highly derived member of a widespread glycoprotein family and a galectin. Identification of homologous proteins in other amphibians reveals that these proteins attained a function in skin long before glues evolved. Yet, major elevations in their expression, besides structural changes in the glycoprotein (increasing its structural disorder and glycosylation), caused the independent rise of glues in at least two frog lineages. Besides providing a model for the chemical functioning of animal adhesive secretions, our findings highlight how recruiting ancient molecular templates may facilitate the recurrent evolution of functional innovations., (© 2024. The Author(s).)

Most described species have not been explicitly included in phylogenetic trees-a problem named the Darwinian shortfall-owing to a lack of molecular and/or morphological data, thus hampering the explicit incorporation of evolution into large-scale biodiversity analyses. We investigate potential drivers of the Darwinian shortfall in tetrapods, a group in which at least one-third of described species still lack phylogenetic data, thus necessitating the imputation of their evolutionary relationships in fully sampled phylogenies. We show that the number of preserved specimens in scientific collections is the main driver of phylogenetic knowledge accumulation, highlighting the major role of biological collections in unveiling novel biodiversity data and the importance of continued sampling efforts to reduce knowledge gaps. Additionally, large-bodied and wide-ranged species, as well as terrestrial and aquatic amphibians and reptiles, are phylogenetically better known. Future efforts should prioritize phylogenetic research on organisms that are narrow-ranged, small-bodied and underrepresented in scientific collections, such as fossorial species. Addressing the Darwinian shortfall will be imperative for advancing our understanding of evolutionary drivers shaping biodiversity patterns and implementing comprehensive conservation strategies.

Most vertebrates develop distinct females and males, where sex is determined by repeatedly evolved environmental or genetic triggers. Undifferentiated sex chromosomes and large genomes have caused major knowledge gaps in amphibians. Only a single master sex-determining gene, the dmrt1-paralogue (dm-w) of female-heterogametic clawed frogs (Xenopus; ZW♀/ZZ♂), is known across >8740 species of amphibians. In this study, by combining chromosome-scale female and male genomes of a non-model amphibian, the European green toad, Bufo(tes) viridis, with ddRAD- and whole genome pool-sequencing, we reveal a candidate master locus, governing a male-heterogametic system (XX♀/XY♂). Targeted sequencing across multiple taxa uncovered structural X/Y-variation in the 5'-regulatory region of the gene bod1l, where a Y-specific non-coding RNA (ncRNA-Y), only expressed in males, suggests that this locus initiates sex-specific differentiation. Developmental transcriptomes and RNA in-situ hybridization show timely and spatially relevant sex-specific ncRNA-Y and bod1l-gene expression in primordial gonads. This coincided with differential H3K4me-methylation in pre-granulosa/pre-Sertoli cells, pointing to a specific mechanism of amphibian sex determination., (© 2024. The Author(s).)

Effective population size (N e ) is a particularly useful metric for conservation as it affects genetic drift, inbreeding and adaptive potential within populations. Current guidelines recommend a minimum N e of 50 and 500 to avoid short-term inbreeding and to preserve long-term adaptive potential respectively. However, the extent to which wild populations reach these thresholds globally has not been investigated, nor has the relationship between N e and human activities. Through a quantitative review, we generated a dataset with 4610 georeferenced N e estimates from 3829 populations, extracted from 723 articles. These data show that certain taxonomic groups are less likely to meet 50/500 thresholds and are disproportionately impacted by human activities; plant, mammal and amphibian populations had a <54% probability of reaching N ̂ e = 50 and a <9% probability of reaching N ̂ e = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller median N ̂ e than unlisted populations, and this was consistent across all taxonomic groups. N ̂ e was reduced in areas with a greater Global Human Footprint, especially for amphibians, birds and mammals, however relationships varied between taxa. We also highlight several considerations for future works, including the role that gene flow and subpopulation structure plays in the estimation of N ̂ e in wild populations, and the need for finer-scale taxonomic analyses. Our findings provide guidance for more specific thresholds based on N e and help prioritise assessment of populations from taxa most at risk of failing to meet conservation thresholds., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past., (© 2024 The Author(s). Physiological Reports published by Wiley Periodicals LLC on behalf of The Physiological Society and the American Physiological Society.)

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

Surveying biodiversity has taken a quantum leap with environmental DNA (eDNA) metabarcoding, an immensely powerful approach lauded for its efficiency, sensitivity, and non-invasiveness. This approach emerges as a game-changer for the elusive realm of endangered and rare species-think nocturnal, environmentally elusive amphibians. Here, we have established a framework for constructing a reliable metabarcoding pipeline for amphibians, covering primer design, performance evaluation, laboratory validation, and field validation processes. The Am250 primer, located on the mitochondrial 16S gene, was optimal for the eDNA monitoring of amphibians, which demonstrated higher taxonomic resolution, smaller species amplification bias, and more extraordinary detection ability compared to the other primers tested. Am250 primer exhibit an 83.8% species amplification rate and 75.4% accurate species identification rate for Chinese amphibians in the in silico PCR and successfully amplified all tested species of the standard samples in the in vitro assay. Furthermore, the field-based mesocosm experiment showed that DNA can still be detected by metabarcoding even days to weeks after organisms have been removed from the mesocosm. Moreover, field mesocosm findings indicate that eDNA metabarcoding primers exhibit different read abundances, which can affect the relative biomass of species. Thus, appropriate primers should be screened and evaluated by three experimental approaches: in silico PCR simulation, target DNA amplification, and mesocosm eDNA validation. The selection of a single primer set or multiple primers' combination should be based on the monitoring groups to improve the species detection rate and the credibility of results., (© 2024 John Wiley & Sons Ltd.)

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals., (© 2023 The Authors. Journal of Experimental Zoology Part A: Ecological and Integrative Physiology published by Wiley Periodicals LLC.)

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches., 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 © 2023 Elsevier Inc. All rights reserved.)

In recent years, environmental DNA (eDNA) has received attention from biologists due to its sensitivity, convenience, labor and material efficiency, and lack of damage to organisms. The extensive application of eDNA has opened avenues for the monitoring and biodiversity assessment of amphibians, which are frequently small and difficult to observe in the field, in areas such as biodiversity survey assessment and detection of specific, rare and threatened, or alien invasive species. However, the accuracy of eDNA can be influenced by factors such as ambient temperature, pH, and false positives or false negatives, which makes eDNA an adjunctive tool rather than a replacement for traditional surveys. This review provides a concise overview of the eDNA method and its workflow, summarizes the differences between applying eDNA for detecting amphibians and other organisms, reviews the research progress in eDNA technology for amphibian monitoring, identifies factors influencing detection efficiency, and discusses the challenges and prospects of eDNA. It aims to serve as a reference for future research on the application of eDNA in amphibian detection., (© 2023 Wiley Periodicals LLC.)

Background: Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are two pathogenic fungi that are a significant threat to amphibian communities worldwide. European populations are strongly impacted and the monitoring of the presence and spread of these pathogens is crucial for efficient decision-making in conservation management., Results: Here we proposed an environmental DNA (eDNA) monitoring of these two pathogenic agents through droplet digital PCR (ddPCR) based on water samples from 24 ponds in Luxembourg. In addition, amphibians were swabbed in eight of the targeted ponds in order to compare the two approaches at site-level detection. This study allowed the development of a new method taking below-Limit of Detection (LOD) results into account thanks to the statistical comparison of the frequencies of false positives in no template controls (NTC) and below-LOD results in technical replicates. In the eDNA-based approach, the use of this method led to an increase in Bd and Bsal detection of 28 and 50% respectively. In swabbing, this resulted in 8% more positive results for Bd. In some samples, the use of technical replicates allowed to recover above-LOD signals and increase Bd detection by 35 and 33% respectively for eDNA and swabbing, and Bsal detection by 25% for eDNA., Conclusions: These results confirmed the usefulness of technical replicates to overcome high levels of stochasticity in very low concentration samples even for a highly sensitive technique such as ddPCR. In addition, it showed that below-LOD signals could be consistently recovered and the corresponding amplification events assigned either to positive or negative detection via the method developed here. This methodology might be particularly worth pursuing in pathogenic agents' detection as false negatives could have important adverse consequences. In total, 15 ponds were found positive for Bd and four for Bsal. This study reports the first record of Bsal in Luxembourg., (© 2023. The Author(s).)

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Since Ficetola et al. (2008) alerted ecologists and conservation biologists to the existence of environmental DNA (eDNA), the number of studies using eDNA has exploded, with a rapidly increasing diversity of research, monitoring, and management objectives. Initial applications focused on amphibians and fishes while today's taxonomic targets span the phylogenetic tree. The environmental media that are sampled have expanded from freshwater to saltwater to soils, and, most recently, to air. In this issue of Molecular Ecology Resources, Lynggaard et al. (Molecular Ecology Resources, 2023) use eDNA captured on air filters to census vertebrate biodiversity in a forest. With a three day, six sample period, 143 sample effort in a nature park in a rural area of Zealand, Denmark, their wild species detections comprised about 25% of the terrestrial vertebrates that are known to occur in the area, including about 33% of the mammal, 17% of the bird, and 60% of the amphibian species. This study demonstrates that air sampling for eDNA has the potential to become a powerful standard method for terrestrial biodiversity assessment that is complementary to traditional methods (e.g., trapping, visual and acoustic observation, collection of scat and hair)., (© 2023 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)