design: For next-generation RNA-sequencing (RNA-Seq), transcriptome libraries were constructed using Ion Total RNA-Seq Kit v2 (Thermo Fisher Scientific, Massachusetts, USA). Briefly, an appropriate amount of polyA RNA was enzymatically fragmented for 20 seconds and purified. The yield and size distribution were assessed with Ag-ilent 2100 Bioanalyzer Instrument and Agilent RNA 6000 Pico Kit (Agilent Tech-nologies, Santa Clara, California, USA). This was followed by hybridisation and liga-tion of adapters to the fragmented mRNA and reverse transcription. Prepared cDNA fragments were purified using magnetic beads, followed by barcoding of each sample using Ion Xpress RNASeq Barcode 116 Kit (Thermo Fisher Scientific,) and PCR-amplification. Barcoded cDNA libraries were purified, and yield and size distri-bution of libraries were assess using Agilent 2100 Bioanalyzer instrument with the DNA 1000 Kit or the Agilent High Sensitivity DNA Kit. Prepared cDNA libraries were diluted to the same molar concentrations and pooled in equal volumes. After-wards, pooled sample was prepared for sequencing using an Ion OneTouch 2 system and Ion PI Hi-Q OT2 200 Kit (Thermo Fisher Scientific) according to accompanying manufacturers instructions. Sequencing was performed on an Ion Proton system, us-ing the Ion PI Hi-Q sequencing 200 Kit (Thermo Fisher Scientific).