RNA-seq analysis provides insight into molecular adaptations of Andrias davidianus

Comparative transcriptomes of nine tissues for the Heilongjiang brown frog (Rana amurensis)

The Heilongjiang brown frog (Rana amurensis) is widely used in traditional Chinese medicine. In particular, the oviduct and skin have been developed into various health products. However, limited numbers of complete genomes of amphibian species have been reported, excluding the Heilongjiang brown frog. Here, the transcriptomes of 45 samples from the liver, spleen, heart, ovaries, thigh muscles, skin, oviduct, stomach and intestine of five Heilongjiang brown frog were reassembled and analyzed. A total of 1,085,532 unigenes with an average length of 676.6 bp and N50 of 722 bp were obtained. Comparative transcriptomics of different tissues detected tissue-specific expression. There were 3248 differentially expressed genes (DEGs) in the ovary, and the number of unique DEGs between the ovary and spleen was the largest. The results of DEGs enrichment showed there were many pathways and items related to protein synthesis and metabolism in the oviduct. The DEGs of the skin were enriched with many bacterial defense items, indicating that there were a large number of antimicrobial peptides in the skin. Thus, these were suitable as biological sources for the development and extraction of antimicrobial peptides. Through the assembly of transcriptome sequencing data and functional annotation of the Heilongjiang brown frog genome, this study provides reference materials for further exploring and utilizing functional gene resources of frogs and lays a foundation for medical research and the development of new products.

Full-length transcriptome assembly of Andrias davidianus (Amphibia: Caudata) skin via hybrid sequencing

The Chinese giant salamander, Andrias davidianus, is the largest amphibian species in the world; it is thus an economically and ecologically important species. The skin of A. davidianus exhibits complex adaptive structural and functional adaptations to facilitate survival in aquatic and terrestrial ecosystems. Here, we report the first full-length amphibian transcriptome from the dorsal skin of A. davidianus, which was assembled using hybrid sequencing and the PacBio and Illumina platforms. A total of 153,038 transcripts were hybrid assembled (mean length of 2039 bp and N50 of 2172 bp), and 133,794 were annotated in at least one database (nr, Swiss-Prot, KEGG, KOGs, GO, and nt). A total of 58,732, 68,742, and 115,876 transcripts were classified into 24 KOG categories, 1903 GO term categories, and 46 KEGG pathways (level 2), respectively. A total of 207,627 protein-coding regions, 785 transcription factors, 27,237 potential long non-coding RNAs, and 8299 simple sequence repeats were also identified. The hybrid-assembled transcriptome recovered more full-length transcripts, had a higher N50 contig length, and a higher annotation rate of unique genes compared with that assembled in previous studies using next-generation sequencing. The high-quality full-length reference gene set generated in this study will help elucidate the genetic characteristics of A. davidianus skin and aid the identification of functional skin proteins.