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

Toads release toxic dry secretions from glands in their skin. Toxin possesses a wide range of biological effects, but little is known about its specific gene expression pattern and regulatory mechanisms. The Asiatic toad (Bufo gargarizans) is widely used to produce toxin. Here, we explored the gene expression of 30 tissue samples from three different skin sites (parotoid gland, dorsal skin, and abdomen skin) of B. gargarizans. After de novo assembly, 783,130 unigenes with an average length of 489 bp (N50 = 556 bp) were obtained. A total of 9,248 significant differentially expressed genes (DEGs) were detected. There were 8,819 DEGs between the parotoid gland and abdomen skin and 1,299 DEGs between the dorsal skin and abdomen skin, while only 1,283 DEGs were obtained between the parotoid gland and dorsal skin. Through enrichment analysis, it was found that the detected differential gene expressions corresponded to the different functions of different skin sites. Our key findings were the genetic expression of toxin secretion, the protection function of skin, and the related genes such as HSD3B, Cyp2c, and CAT, LGALS9. In conclusion, we provide useful transcript resources to study the gene expression and gene function of B. gargarizans and other amphibians. The detected DEGs between different sites of the skin provided better insights into the genetic mechanisms of toxin secretion and the protection function of skin for amphibians.