library_construction_protocol: Total RNA was extracted using TAKARA RNAiso Plus#9109 following the manufacturer's instructions and checked for a RIN number to inspect RNA integrity by an Agilent Bioanalyzer 2100 (Agilent technologies, Santa Clara, CA, US). RNA degradation and contamination was monitored on 1% agarose gels. RNA purity was checked using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA) . RNA concentration was measured using Qubit® RNA Assay Kit in Qubit®2.0 Flurometer (Life Technologies, CA, USA). Quality RNA samples were converted into cDNA libraries using VAHTSTM mRNA-seq V2 Library Prep Kit for Illumina® (NR601, Vazyme, NanJing, China). VAHTSTM DNA Clean Beads (Vazyme #N411) in this study was used to purify the fragments during the process of library constructing. The purified products were enriched with 12-15 cycles of PCR to create the final cDNA library. Finally, libraries were sequenced on the Illumina Hiseq X Ten according the manufacturer's protocols.