The RNA samples were prepared by  using either the RNeasy Mini Kit (QIAGEN) or the NucleoSpin RNA Mini Kit (MACHEREY-NAGEL). 50ng of total RNA was used for preparing sequencing library. RNA libraries for RNA seq were prepared using a bulk version of the MARS-Seq procedure. Briefly, 50 ng of total RNA from each sample was barcoded during reverse transcription (RT) and then pooled. After cleanup with Agencourt AMPure XP beads (Beckman Coulter), the pooled samples underwent second-strand synthesis and were linearly amplified by T7 in-vitro transcription. The resulting RNA was fragmented and converted into a sequencing-ready library by adding Illumina sequences during ligation, RT, and PCR. Libraries were quantified using Qubit and TapeStation, as well as by qPCR for the ACTB gene.