Andrias davidianus
RNA sequencing was performed by the Illumina Hiseq4000 system (LC Sciences, Hangzhou, China). Raw reads were subjected to an in-house program, ACGT101-miR (v4.2) to remove adapter dimers, junk, low complexity, common RNA families (rRNA, tRNA, snRNA, snoRNA) and repeats. Subsequently, unique sequences with length in 18~26 nucleotide were mapped to specific species precursors in miRBase 22.1 by BLAST search to identify known miRNAs and novel 3p- and 5p- derived miRNAs