Spatio-temporal analysis of toxigenic genes expression in the growing Bufo gargarizans based on RNA sequencing data.

BACKGROUND: Bufo gargarizans Cantor, a widely distributed amphibian species in Asia, produces and releases toxins through its retroauricular and granular glands. Although various tissues have been sequenced, the molecular mechanisms underlying the toxin production remain unclear. To elucidate these mechanisms, abdominal skin (non-toxic secretory glands) and retroauricular gland (toxic secreting glands) samples were collected at different time points (3, 6, 12, 24, and 36 months) for RNA sequencing (RNA-seq) and analysis. RESULTS: In comparison to the S group during the same period, a total of 3053, 3026, 1516, 1028, and 2061 differentially expressed genes (DEGs) were identified across five developmental stages. Gene Ontology (GO) analysis revealed that DEGs were primarily enriched in biological processes including cellular processes, single-organism processes, metabolic processes, and biological regulation. In terms of cellular components, the DEGs were predominantly localized in the cell and cell parts, whereas molecular function indicated significant enrichment in binding and catalytic activity. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that the metabolism and synthesis of various substances, such as lipid metabolism, cofactor and vitamin metabolism, tryptophan metabolism, steroid biosynthesis, and primary bile acid biosynthesis, were accompanied by the development of toads. Additionally, using trend analysis, we discovered candidate genes that were upregulated in the retroauricular glands during development, and the abundance of these genes in the abdominal skin was extremely low. Finally, we identified 26 genes that are likely to be involved in toxin production and that are likely to be involved in toxin anabolism. CONCLUSION: Overall, these results provide new insights into the genes involved in toxin production in B. gargarizans, which will improve our understanding of the molecular mechanisms underlying toxigenic gene expression.

Context-dependent effects of thermal acclimation on physiological correlates of animal personality in Asiatic toads.

Persistent individual variation in behaviour, or 'personality', is a widespread phenomenon in animals, and understanding the evolution of animal personality is a key task of current biology. Natural selection has been proposed to promote the integration of personality with animal 'intrinsic states', such as metabolic or endocrine traits, and this integration varies with ecological conditions. However, these external ecological modulatory effects have rarely been examined. Here, we investigate the effects of thermal acclimation on between-individual covariations between physiology and behaviour in Asiatic toads (Bufo gargarizans) along an altitudinal gradient. Our results reveal that the thermal modulatory effects on the covariations depend on the altitudinal population. Specifically, at low altitudes, between-individual covariations are highly plastic, with risk-taking behaviour covarying with baseline glucocorticoids (GCs) under warm acclimation, but risk-taking and exploration behaviour covarying with resting metabolic rate (RMR) under cold acclimation. In contrast, between-individual covariations are relatively fixed at high altitudes, with risk-taking behaviour consistently covarying with baseline GCs. Furthermore, at low altitudes, changes in covariations between RMR and personality are associated with adjustment of energy management models. Evidently, animal physiological states that determine or covary with personality can adapt according to the seasonal thermal environment and the thermal evolutionary background of populations. Our findings highlight the importance of a multi-system physiological approach to understand the evolution of animal personality.

Archival mitogenomes identify invasion by the Batrachochytrium dendrobatidis CAPE lineage caused an African amphibian extinction in the wild.

Outbreaks of emerging infectious diseases are influenced by local biotic and abiotic factors, with host declines occurring when conditions favour the pathogen. Deterioration in the population of the micro-endemic Tanzanian Kihansi spray toad (Nectophrynoides asperginis) occurred after the construction of a hydropower dam, implicating habitat modification in this species decline. Population recovery followed habitat augmentation; however, a subsequent outbreak of chytridiomycosis caused by Batrachochytrium dendrobatidis (Bd) led to the spray toad's extinction in the wild. We show using spatiotemporal surveillance and mitogenome assembly of Bd from archived toad mortalities that the outbreak was caused by invasion of the BdCAPE lineage and not the panzootic lineage BdGPL. Molecular dating reveals an emergence of BdCAPE across southern Africa overlapping with the timing of the spray toad's extinction. That our post-outbreak surveillance of co-occurring amphibian species in the Udzungwa Mountains shows widespread infection by BdCAPE yet no signs of ill-health or decline suggests these other species can tolerate Bd when environments are stable. We conclude that, despite transient success in mitigating the impact caused by dams' construction, invasion by BdCAPE caused the ultimate die-off that led to the extinction of the Kihansi spray toad.

Insights into the evolution of photoreceptor oil droplets in frogs and toads.

Photoreceptor oil droplets (ODs) are spherical organelles placed most commonly within the inner segment of the cone photoreceptors. Comprising neutral lipids, ODs can be either non-pigmented or pigmented and have been considered optically functional in various studies. Among living amphibians, ODs were only reported to occur in frogs and toads (Anura), while they are absent in salamanders and caecilians. Nonetheless, the limited understanding of their taxonomic distribution in anurans impedes a comprehensive assessment of their evolution and relationship with visual ecology. We studied the retinae of 134 anuran species, extending the knowledge of the distribution of ODs to 46 of the 58 currently recognized families, and providing a new perspective on this group that complements the available information from other vertebrates. The occurrence of ODs in anurans shows a strong phylogenetic signal, and our findings revealed that ODs evolved at least six times during the evolutionary history of the group, independently from other vertebrates. Although no evident correlation was found between OD occurrence, adult habits and diel activity, it is inferred that each independent origin involves distinct scenarios in the evolution of ODs concerning photic habits. Furthermore, our results revealed significant differences in the size of the ODs between nocturnal and arrhythmic anurans relative to the length of the cones' outer segment.

Harlequin frog tadpoles-comparative buccopharyngeal morphology in the gastromyzophorous tadpoles of the genus Atelopus (Amphibia, Anura, Bufonidae), with discussion on the phylogenetic and evolutionary implication of characters.

The Neotropical genus Atelopus is the most diverse genus of bufonids comprising 99 species. Tadpoles of these frogs are readily distinguished based on the presence of a belly sucker, used by them to stay attached to rocks in fast-flowing streams. Despite their intriguing biology, information about their anatomy is scarce and many morphological systems are unknown. We describe the buccopharyngeal cavity of five Atelopus species. The Atelopus buccopharyngeal cavity is characterized by (1) presence of a pendulum-like papillae in the prenarial arena, (2) presence of a glandular zone in the prenarial arena, (3) narial vacuities, (4) conical median ridge, (5) absence of buccal roof arena papillae, (6) absence of buccal roof pustulations, (7) single pair of infralabial papillae, (8) absence of lingual papillae, and (9) absence of pustulations in the buccal floor. We propose that characters 1, 2, and 3 are new synapomorphies for the genus. We also propose that the presence of a single pair of infralabial papillae is a synapomorphy for bufonid. Finally, we discuss the convergent evolution of gastromyzophorous and suctorial tadpoles withing anurans.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians.

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Baseline and stress-induced steroid plasma levels and immune function vary annually and are associated with vocal activity in male toads (Rhinella icterica).

Theoretical models predict that elevated androgen and glucocorticoid levels in males during the reproductive season promote immunosuppression. However, some studies report decreased stress response during this season. This study investigated annual variation in plasma corticosterone and testosterone levels, plasma bacterial killing ability (BKA), and neutrophil to lymphocyte ratio (NLR) in free-living male toads (Rhinella icterica). Toads were sampled in the field (baseline) and 1 h-post restraint over five months, and we considered the occurrence of vocal activity. Baseline corticosterone, testosterone, and BKA showed higher values during the reproductive period, specifically in calling male toads. The NLR was similar throughout the year, but higher values were observed in calling toads. Moreover, baseline NLR and BKA were positively correlated with both testosterone and corticosterone, suggesting higher steroid levels during reproduction are associated with enhanced cellular and humoral immunity. Despite fluctuation of baseline values, post-restraint corticosterone levels remained uniform over the year, indicating that toads reached similar maximum values throughout the year. Testosterone levels decreased following restraint before one specific reproductive period but increased in response to restraint during and after this period. Meanwhile, BKA decreased due to restraint only after the reproductive period, indicating immune protection and resilience to immunosuppression by stressors associated with steroid hormones during reproduction. Our results show that baseline and stress-induced hormonal and immune regulation varies throughout the year and are associated with vocal activity in R. icterica males, indicating a possible compromise between steroids and immune function in anuran males.

Thermal conditions predict intraspecific variation in senescence rate in frogs and toads.

Variation in temperature is known to influence mortality patterns in ectotherms. Even though a few experimental studies on model organisms have reported a positive relationship between temperature and actuarial senescence (i.e., the increase in mortality risk with age), how variation in climate influences the senescence rate across the range of a species is still poorly understood in free-ranging animals. We filled this knowledge gap by investigating the relationships linking senescence rate, adult lifespan, and climatic conditions using long-term capture-recapture data from multiple amphibian populations. We considered two pairs of related anuran species from the Ranidae (Rana luteiventris and Rana temporaria) and Bufonidae (Anaxyrus boreas and Bufo bufo) families, which diverged more than 100 Mya and are broadly distributed in North America and Europe. Senescence rates were positively associated with mean annual temperature in all species. In addition, lifespan was negatively correlated with mean annual temperature in all species except A. boreas In both R. luteiventris and A. boreas, mean annual precipitation and human environmental footprint both had negligible effects on senescence rates or lifespans. Overall, our findings demonstrate the critical influence of thermal conditions on mortality patterns across anuran species from temperate regions. In the current context of further global temperature increases predicted by Intergovernmental Panel on Climate Change scenarios, a widespread acceleration of aging in amphibians is expected to occur in the decades to come, which might threaten even more seriously the viability of populations and exacerbate global decline.

Immune and endocrine alterations at the early stage of inflammatory assemblage in toads after stimulation with heat-killed bacteria (Aeromonas hydrophila).

The red-leg syndrome in amphibians is a condition commonly associated with the bacteria Aeromonas hydrophila and has led to population declines. However, there is little information concerning the inflammatory assemblage in infected anurans. We evaluated immune and endocrine alterations induced by stimulation with heat-killed A. hydrophila injected in Rhinella diptycha toads. Control animals were not manipulated, while the others were separated into groups that received intraperitoneal injection of 300 µl of saline or heat-killed bacteria: groups A1 (3 × 107 cells), A2 (3 × 108 cells), and A3 (3 × 109 cells). Animals were bled and euthanized six hours post-injection. We evaluated neutrophil: lymphocyte ratio (NLR), plasma bacterial killing ability (BKA), testosterone (T), melatonin (MEL), and corticosterone (CORT) plasma levels. Heat-killed A. hydrophila increased CORT and NLR, and decreased MEL, especially at higher concentrations. There was no effect of treatment on T and BKA. We then selected the saline and A3 groups to conduct mRNA expression of several genes including glucocorticoid receptor (GR), toll-like receptor-4 (TLR-4), interferon-γ (IFN-γ), interleukin (IL)-1ß, IL-6, and IL-10. We found higher expression of IL-6, IL-1ß, IL-10, and IFN-γ in group A3 compared to the saline group. These results indicate the beginning of an inflammatory assemblage, notably at the two highest concentrations of bacteria, and give a better understanding of how anurans respond to an infection within an integrated perspective, evaluating different physiological aspects. Future studies should investigate later phases of the immune response to elucidate more about the inflammation in amphibians challenged with A. hydrophila.

Toxicity assessment of a tebuconazole-based fungicide on the embryo-larval development of the common south American toad Rhinella arenarum.

Agrochemicals cause diverse effects on aquatic communities, and amphibian species are particularly threatened due the high susceptibility to contamination. Present study evaluates the toxicity of a widely used fungicide tebuconazole (Trigal®) by the assessment of mortality and developmental alterations at acute, subchronic, and chronic exposure during the embryo-larval development of the South American toad Rhinella arenarum. Also, the sensitivity of the different embryonic stages was evaluated with 24-h pulse exposure treatments. The results demonstrated that larvae were more sensitive than embryos at acute exposure (LC50-24 and 96 h = 74.62, 31.92 mg/L and 24.27, 16.81 mg/L for embryos and larvae, respectively). Nevertheless, embryos toxicity increased significantly achieving a sensitivity very similar to larvae at chronic exposure (LC50-168 and 504 h = 13.31, 4.35 mg/L and 14.47, 6.83 mg/L for embryos and larvae, respectively). Embryos exhibited several sublethal effects from 5 mg/L at 96 h onwards, such as delayed development, reduce body size, edemas, tail/axial flexures, weakness, and absence of movements. The teratogenic index at 96 h was 10.13, indicating the severe teratogenic potential of the fungicide. 24-h pulse exposure treatments showed an increased sensitivity in intermediate stages as S.11, S.18, S20, and S.23 (NOEC-96 h = 100, 200, 75, and 20 mg/L, respectively), while stage S.25 was the most sensitive to the fungicide (NOEC-96 h = 5 mg/L). About metamorphic process, tebuconazole caused an acceleration of metamorphosis at the lowest concentration (0.001 mg/L), but also an increase in mortality and in addition, significant differences in the weight in all treatments. The results obtained throughout this work indicate that tebuconazole cause several adverse effects in Rhinella arenarum embryo-larval development.

Elevational variation in metabolic rate, feeding capacity and their associations in the Asiatic toad Bufo gargarizans.

Foraging behavior is known to place demands on the metabolic characteristics of anurans. Active foragers feeding on sedentary prey typically have high aerobic capacity and low anaerobic capacity, whereas sit-and-wait foragers feeding on active and mobile prey have the opposite pattern. Thus, the energetic demands of foraging may influence their metabolic adaptations to harsh environments, such as high elevations. Anurans that engage in active foraging have been found to increase maximum metabolic rate (MMR) and aerobic scope (AS, the difference between MMR and resting metabolic rate, RMR) at high elevations. However, data are lacking in amphibian ambush foragers. In this study, we examined the RMR, MMR, AS, and feeding capacity of a sit-and-wait forager ─the Asiatic toad (Bufo gargarizans), from two populations that are in close geographic proximity but differ by 1350 m in elevation. Our results show that there is no elevational variation in RMR and feeding capacity in either males or females. However, there are sex-specific variations in MMR and AS along an elevational gradient; females from high elevations have lower MMR and smaller net AS than their counterparts from low elevations while males maintain similar MMR and net AS across elevations. Furthermore, aerobic performances do not appear to be associated with feeding capacity at either the individual or population level. Our results support the hypothesis that sit-and-wait foragers may not increase their aerobic capacity as a strategy in hypoxic and low food availability environments and the role of sex in these adaptive adjustments should not be overlooked.

Changes in thermal sensitivity of Rhinella arenarum tadpoles (Anura: Bufonidae) exposed to sublethal concentrations of different pesticide fractions (Lorsban® 75WG).

The intensive use of agrochemicals and the rapid increase of global temperatures have modified the thermal conditions of aquatic environments, thus increasing amphibians' vulnerability to global warming and positioning them at great risk. Commercial formulations of chlorpyrifos (COM) are the pesticides most widely used in agricultural activities, with a high toxic potential on amphibians. However, little is known about the separate effects of the active ingredient (CPF) and adjuvants (AD). We studied the thermal sensitivity at different concentrations and pesticide fractions in Rhinella arenarum tadpoles, on thermal tolerance limits (CTmax = Critical thermal maximum and CTmin = Critical thermal minimum), swimming speed (Ss), Optimum temperature (Top), and Thermal breadth 50 (B50). Our results demonstrate that the pesticide active ingredient, the adjuvants, and the commercial formulation of chlorpyrifos differentially impair the thermal sensitivity of R. arenarum tadpoles. The pesticide fractions affected the heat and the cold tolerance (CTmax and CTmin), depending on the concentrations they were exposed to. The locomotor performance (Ss, Top, and B50) of tadpoles also varied among fractions, treatments, and environmental temperatures. In the context of climate change, the outcomes presented are particularly relevant, as mean temperatures are increasing at unprecedented rates, which suggests that tadpoles inhabiting warming and polluted ponds are currently experiencing deleterious conditions. Considering that larval stages of amphibians are the most susceptible to changing environmental conditions and the alarming predictions about environmental temperatures in the future, it is likely that the synergism between high temperatures and pesticide exposure raise the threat of population deletions in the coming years.

Response Surface Methodology (RSM) for Optimizing Protein Extraction from Housefly (Musca domestica) Larvae Fed with Toad and Its Structural Characterization.

With the global population on the rise, an escalating interest exists in environmentally sustainable and friendly protein sources. Insects have emerged as multifaceted resources, viewed not only as potential food items, but also as sources of traditional medicines and proteins. This study utilized response surface methodology (RSM) to ascertain the optimal extraction conditions for proteins from Musca domestica used in toad feeding, denoted as MDPs-T. The yield of MDPs-T was elevated to 18.3% ± 0.2% under these optimized conditions. Subsequently, the particle size, ζ-potentials, and structures of MDPs-T were analyzed and compared with the proteins derived from Musca domestica fed on a normal diet (MDPs-ND). This comparative analysis utilized a range of advanced techniques, involving UV spectroscopy, Fourier transform infrared spectroscopy (FT-IR), X-ray diffraction (XRD), high-performance gel permeation chromatography (HPGPC), and scanning electron microscopy (SEM). The outcomes have revealed a marginal disparity in the physical and chemical properties between MDPs-T and MDPs-ND. Derosination led to a reduction in the particle size of the MDPs by 10.98% to 62.81%. MDPs-T exhibited a higher proportion of low-molecular-weight components relative to MDPs-ND. Additionally, in a comparative analysis of amino acids, MDPs-T displayed a greater abundance of essential and total amino acids relative to MDPs-ND. Consequently, MDPs-T holds potential as a valuable food supplement for human consumption or as a nutrient-rich feed supplement for animals.

Isolation, Identification and Chemical Modification of Bufadienolides from Bufo melanostictus Schneider and Their Cytotoxic Activities against Prostate Cancer Cells.

The traditional Chinese medicine toad venom (Venenum bufonis) has been extensively used to treat various diseases, including cancers, in China and other Southeast Asian countries. The major constituents of toad venom, e.g., bufadienolides and alkaloids, exhibit broad-spectrum pharmacological effects in cancers. Herein, two new bufadienolides (1 and 2), along with eleven known compounds (3-13) were successfully isolated from Bufo melanostictus Schneider. Their structures were elucidated by extensive spectroscopic data and X-ray diffraction analysis. Furthermore, four lactam derivatives were synthesized through the transformation of bufadienolides lactones. The inhibitory effects of these compounds against human prostate cancer cell lines PC-3 and DU145 were evaluated. The outcomes indicated a notable trend, with a substantial subset displaying nanomolar range IC50 values against PC-3 and DU145 cells, underscoring their pronounced cytotoxicity. Moreover, a noteworthy distinction surfaces, wherein lactones consistently outperformed their lactam counterparts, further validating their heightened potency for the treatment of prostate cancer. This study contributes significant preclinical evidence substantiating the therapeutic viability of bufadienolides and toad venom as intervention strategies for prostate cancer.

Diversity of Acanthocephala parasites in Neotropical amphibians.

Acanthocephalans constitute a small taxonomic group related to rotifers and specialized in a parasitic lifestyle. Anurans act as paratenic and definitive hosts and infections always occur trophically. Our objective is to describe and summarize the richness of acanthocephalans in Neotropical anurans. We conducted a literature review in the main research databases, compiling data published until August 2021. We identified 66 articles with records of acanthocephalan-anuran association, 53.03% were carried out in Brazil. We detected 108 species of anurans from 11 families parasitized by acanthocephalans. With the exception of Bufonidae, Hylidae and Leptodactylidae, which are relatively well-studied families, interaction with acanthocephalans remains largely unexplored for most anuran species. We found six families of acanthocephalans: Centrorhynchidae, Echinorhynchidae, Oligacanthorhynchidae, Cavisomidae, Neoechinorhynchidae and Plagiorhynchidae. Centrorhynchidae and Echinorhynchidae presented the largest number of taxa associated with anurans. The largest number of records corresponded to acanthocephalans in the larval stage (cystacanths), for which anurans act as paratenic hosts. We observed a lack of specific taxonomic resolution in the identifications of most reports, because a large part of the records in the larval stage make morphological identification difficult. Brazil, Mexico, Paraguay, Argentina, Ecuador and Peru are the countries with the most records, while Costa Rica, Venezuela, Colombia, Chile and Uruguay exhibited the lowest publication numbers, resulting in gaps in the distribution of acanthocephalans. We expanded the known number of anuran species parasitized by acanthocephalans, compared to the last published review. Overall, we aim to contribute to the understanding of diversity within this intriguing but understudied group.

[Synthesis and tissue distribution of bufadienolides in Bufo bufo gargarizans].

This study explored the biosynthesis of bufadienolides(BDs) in Bufo bufo gargarizans to solve the dilemma of the decreasing resources of B. bufo gargarizans and provide a theoretical basis for the sustainable utilization of the resources. Ultra-high performance liquid chromatography-Orbitrap-mass spectrometry(UHPLC-Orbitrap-MS) was employed to detect the synthesis sites of BDs in B. bufo gargarizans, and the results were verified by desorption electrospray ionization-mass spectrometry imaging(DESI-MSI) and homogenate incubation experiments. BDs in B. bufo gargarizans had the highest content in the liver and the highest concentration in the gallbladder, in addition to the parotid gland and skin, which suggested that the liver could synthesize BDs. The results of DESI-MSI also showed that BDs were mainly enriched in the liver rather than the immature parotid gland. The incubation experiment of liver homogenates demonstrated the liver of B. bufo gargarizans had the ability to synthesize BDs. This study showed that the liver was a major organ for the synthesis of BDs in B. bufo gargarizans during metamorphosis, development, and growth, which provided strong theoretical support for the biosynthesis of BDs and the sustainable utilization of B. bufo gargarizans resources.

[Molecular identification of animal species of Bufonis Venenum from secretions].

The animal species is one of the key factors affecting the quality of Bufonis Venenum. The quality of Bufonis Venenum derived from Bufo bufo gargarizans is significantly higher than that from B. melanostictus. Since Bufonis Venenum is from secretions, the conventional identification methods are difficult to identify the animal species due to the lack of the appearance and morphology of the animals. The rapid development of molecular identification technology has provided new methods for the identification of Bufonis Venenum. However, because of the low content and serve degradation of residual DNA in secretions, the research on the molecular identification of Chinese medicinal materials from secretions remains to be carried out. To understand the animal species of Bufonis Venenum, this study collected 83 samples of Bufonis Venenum, including 7 commercially available samples, 5 reference medicinal materials, and 71 animal samples from which Bufonis Venenum was prepared according to the method in the 2020 edition of the Chinese Pharmacopoeia. Different DNA extraction methods were used and compared, and the mitochondrial 16S rRNA gene fragments were amplified, on the basis of which the phylogenetic trees were built. Finally, molecular identification of the animal species of the samples was performed. The results showed that the DNA extracted from Bufonis Venenum by the reagent kit had good quality, and 16S rRNA sequences were successfully amplified from 80 out of the 83 samples. In addition, 71 16S rRNA sequences of the animal species of Bufonis Venenum were downloaded from GenBank. The phylogenetic trees constructed based on the neighbor-joining(NJ) method and the Bayesian inference(BI) method showed that the samples derived from B. bufo gargarizans and B. melanostictus were clustered into separate monophyletic clades, with the support of 100%(NJ) and 1.00(BI), respectively. The animal species of both commercially available samples and reference medicinal materials were B. bufo gargarizans. In conclusion, DNA can be extracted from Bufonis Venenum derived from secretions, and the 16S rRNA gene sequences can be amplified, which can be used for molecular identification of the animal species of Bufonis Venenum. The findings provide a reference for the quality control of Bufonis Venenum and the identification of animal species of medicinal materials derived from secretions.

Mining and Characterization of Indolethylamine N-Methyltransferases in Amphibian Toad Bufo gargarizans.

Strong, psychedelic indolethylamines (IAAs) are typically present in trace amounts in the majority of species, but they build up significantly in the skin of amphibian toads, especially N-methylated 5-hydroxytryptamine (5-HT) analogues. However, there is no pertinent research on the investigation of indoleamine N-methyltransferase (INMT) in amphibians, nor is there any adequate information on the key amino acids that influence the activity of known INMTs from other species. Herein, we focused on Bufo toad INMT (BINMT) for the first time and preliminarily identified BINMT 1 from the transcriptomes of Bufo gargarizans active on tryptamine, 5-HT, and N-methyl-5-HT. We established the enzyme kinetic characteristics of BINMT 1 and identified the essential amino acids influencing its activity via molecular docking and site-directed mutagenesis. Subsequently, we carried out sequence alignment and phylogenetic tree analysis on 43 homologous proteins found in the genome of B. gargarizans with BINMT 1 as the probe and selected seven of them for protein expression and activity assays. It was found that only three proteins possessing the highest similarity to BINMT 1 had INMT activity. Our research unveils the binding residues of BINMT for 5-HT analogues for the first time and initiates the study of INMTs in amphibian toads, serving as a tentative reference for further study of BINMT and providing insight into the comprehension of BINMT's catalytic mechanism and its role in the biosynthesis of 5-HT analogues in Bufo toads. It also contributes to the expansion of the INMT library to help explore and explain interspecies evolution in the future.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease.

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

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

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