To facilitate genome annotation, RNA sequencing (RNA-Seq) was conducted. For Illumina RNA-Seq, the total RNA independently extracted from the six tissues of the female sample was first qualified using a Nanodrop spectrophotometer, Qubit fluorometer, and Agilent 2100 Bioanalyzer to evaluate purity (A260/A280 ratio raging from 1.8 and 2.0), concentration, and integrity (RNA Integrity Number, RIN  7.5). Subsequently, six tissue-specific cDNA libraries were constructed independently from the aforementioned tissues and sequenced on the Illumina NovaSeq 6000 platform with a PE 150bp mode. Raw reads were filtered to remove adapter-contaminated reads, reads containing more than 10% unknown nucleotides (Ns) and low-quality bases (with more than 20% bases having a Phred quality score < 5)