Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans)

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.

Transcriptomic and gene-family dynamic analyses reveal gene expression pattern and evolution in toxin-producing tissues of Asiatic toad (Bufo gargarizans)

Comprising a major clade of Anura, toads produce and secrete numerous toxins from both the parotoid glands behind their eyes and their dorsal skin. These toxins, made of various proteins and compounds, possess pharmacological potential to be repurposed to benefit human health. However, the detailed genetic regulation of toad toxin production is still poorly understood. A recent publication uncovering the genome of the representative Asiatic toad (Bufo gargarizans) provides a good reference to resolve this issue. In the present study, we sequenced the transcriptomes of parotoid gland, dorsal skin and liver from the Asiatic toad. Combining our data with 35 previously published transcriptomes across eight different tissues from the same species but from different locations, we constructed a comprehensive gene co-expression network of the Asiatic toad with the assistance of the reference genome assembly. We identified 2,701 co-expressed genes in the toxin-producing tissues (including parotoid gland and dorsal skin). By comparative genomic analysis, we identified 599 expanded gene families with 2,720 genes. Through overlapping these co-expressed genes in the toad toxin-producing tissues, we observed that three cytochrome P450 (Cyp) family members (Cyp27a1, Cyp2c29, and Cyp2c39) were significantly enriched in pathways related to cholesterol metabolism. Cholesterol is a critical precursor to steroids, and the known main steroidal toxins of bufadienolides are considered as the major bioactive components in the parotoid glands of Asiatic toad. We found 3-hydroxy-methylglutaryl CoA reductase (hmgcr), encoding the major rate-limiting enzyme for cholesterol biosynthesis, appears with multiple copies in both Asiatic toad and common toad, possibly originating from a tandem duplication event. The five copies of hmgcr genes consistently displayed higher transcription levels in the parotoid gland when compared with the abdominal skin, suggesting it as a vital candidate gene in the involvement of toad toxin production. Taken together, our current study uncovers transcriptomic and gene-family dynamic evidence to reveal the vital role of both expanded gene copies and gene expression changes for production of toad toxins.