Total RNA was extracted by using TRIzol/chloroform extraction method and further treated with DNase by using RNA Clean & Concentrator (Zymo Research #R1014). Poly(A)+ RNA was selected by using the NEBNext® Poly(A) mRNA Magnetic Isolation Module (NEB #E7490S). RNA-seq libraries were prepared according to the manual for NEBNext® Ultra™ II RNA Library Prep Kit (NEB #E7770L).