library_construction_protocol: Messenger RNA was isolated using the NEBNext Poly(A) mRNA Magnetic Isolation Module. RNA quality was determined using an Agilent 2100 BioAnalyzer. Sequencing libraries were constructed using the NEBNext Ultra RNA Library Prep Kit (Illumina Inc., San Diego, CA, USA). The cDNA library was pair-end sequenced on Illumina HiSeq 2500 at 2×150 bp read length. Adapters and low-quality sequence reads were removed using Trimmomatic (version 0.36).