show Abstracthide AbstractTotal RNA extraction was performed by using the hot acidic phenol method. Then genomic DNA was removed from RNA sample using DNase I. Integrity and size distribution was checked on a Bioanalyzer 2100 with RNA 6000 Nano Labchips (Agilent technologies, Santa Clara, CA, USA). The samples was standardized to 500 ng/ul. Two libraries was constructed by TruSeq RNA Sample Preparation Kit according to the product instruction (Illumina), and sequenced on Illumina HiSeq 2500 (Illumina, San Diego, USA) with 150 bp paired-end reads.