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.

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

Field-based molecular detection of Batrachochytrium dendrobatidis in critically endangered Atelopus toads and aquatic habitats in Ecuador.

Batrachochytrium dendrobatidis (Bd) is a lethal fungal species that parasitizes vertebrates and is associated with the worldwide decline of amphibian populations. The development of sensitive, rapid detection methods, particularly DNA-based techniques, is critical for effective management strategies. This study evaluates the efficacy of DNA extraction and a portable PCR device in a mountable field laboratory setup for detecting Bd near the habitats of three critically endangered Atelopus toad species in Ecuador. We collected skin swabs from Atelopus balios, A. nanay, and A. bomolochos, and environmental DNA (eDNA) samples from streams in Andean and coastal regions of Ecuador. For eDNA, a comparison was made with duplicates of the samples that were processed in the field and in a standard university laboratory. Our findings revealed Bd detection in eDNA and swabs from 6 of 12 water samples and 10 of 12 amphibian swab samples. The eDNA results obtained in the field laboratory were concordant with those obtained under campus laboratory conditions. These findings highlight the potential of field DNA-based monitoring techniques for detecting Bd in amphibian populations and their aquatic habitats, particularly in remote areas. Furthermore, this research aligns with the National Action Plan for the Conservation of Ecuadorian Amphibians and contributes to the global effort to control this invasive and deadly fungus.

Using genome-wide data to ascertain taxonomic status and assess population genetic structure for Houston toads (Bufo [= Anaxyrus] houstonensis).

The Houston toad (Bufo [= Anaxyrus] houstonensis) is an endangered amphibian with a small geographic range. Land-use changes have primarily driven decline in B. houstonensis with population supplementation predominant among efforts to reduce its current extinction risk. However, there has been historic uncertainty regarding the evolutionary and conservation significance of B. houstonensis. To this end, we used 1170 genome-wide nuclear DNA markers to examine phylogenetic relationships between our focal taxon, representatives of the Nearctic B. americanus group, and B. nebulifer, a sympatric Middle American species. Phylogenetic analyses indicate B. houstonensis is a taxon that is distinct from B. americanus. We corroborated such genetic distinctiveness with an admixture analysis that provided support for recent reproductive isolation between B. americanus and B. houstonensis. However, ABBA-BABA tests for ancient admixture indicated historic gene flow between Nearctic species while no signal of historic gene flow was detected between Nearctic and Middle-American species. We used an admixture analysis to recognize four Management Units (MU) based on observed genetic differentiation within B. houstonensis and recommend captive propagation, population supplementation, and habitat restoration efforts specific to each MU. Our results re-affirm the evolutionary novelty of an endangered relict.

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo).

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

Comparative molecular and conventional cytogenetic analyses of three species of Rhinella (Anura; Bufonidae).

The genus Rhinella corresponds to a group of anurans characterized by numerous taxonomic and systemic challenges, leading to their organization into species complexes. Cytogenetic data for this genus thus far are limited to the diploid number and chromosome morphology, which remain highly conserved among the species. In this study, we analyse the karyotypes of three species of the genus Rhinella (Rhinella granulosa, Rhinella margaritifera, and Rhinella marina) using both classical (conventional staining and C-banding) and molecular (FISH-fluorescence in situ hybridization with 18S rDNA, telomeric sequences, and microsatellite probes) cytogenetic approaches. The aim of this study is to provide data that can reveal variations in the distribution of repetitive sequences that can contribute to understanding karyotypic diversification in these species. The results revealed a conserved karyotype across the species, with 2n = 22 and FN = 44, with metacentric and submetacentric chromosomes. C-banding revealed heterochromatic blocks in the pericentromeric region for all species, with a proximal block on the long arms of pairs 3 and 6 in R. marina and on the short arms of pairs 4 and 6 in R. margaritifera. Additionally, 18S rDNA probes hybridized to pair 5 in R. granulosa, to pair 7 in R. marina, and to pair 10 in R. margaritifera. Telomeric sequence probes displayed signals exclusively in the distal region of the chromosomes, while microsatellite DNA probes showed species-specific patterns. These findings indicate that despite a conserved karyotypical macrostructure, chromosomal differences exist among the species due to the accumulation of repetitive sequences. This variation may be attributed to chromosome rearrangements or differential accumulation of these sequences, highlighting the dynamic role of repetitive sequences in the chromosomal evolution of Rhinella species. Ultimately, this study emphasizes the importance of the role of repetitive DNAs in chromosomal rearrangements to elucidate the evolutionary mechanisms leading to independent diversification in the distinct phylogenetic groups of Rhinella.