The total RNA was extracted according to the instruction manual of the TRlzol Reagent (Life technologies, California, USA). Sequencing libraries were generated using NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample.