Total RNA was isolated and underwent quality assessment with Agilent 2000 bioanalyzer (Agilent Technologies, Netherlands). The library was generated with NEBNext Ultra II Directional RNA-Seq Kit (New England Biolabs, UK). Next, the mRNA was isolated, fragmented and synthesized using the Poly(A) RNA Selection Kit (Lexogen, Austria) according to the manufacturer’s instructions. The Illumina indexes 1-12 were used for indexing. The raw sequencing data underwent quality control by FastQC(https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). HiSeq X 10 (Illumina, USA) was utilized for high-throughput paired-end 100 sequencing, and the adapter and low-quality reads (< Q20) were excluded by FASTX_Trimmer ((http://hannonlab.cshl.edu/fastx_toolkit/) and BBMap ((https://sourceforge.net/projects/bbmap). For analysis, the reads were mapped to the reference genome via TopHat [23 (link)]and the expression levels were calculated based on the Fragments Per Kilobase Million (FPKM) mapped reads from Cufflinks [24 (link)]. Quantile normalization was performed via EdgeR (R Development Core Team, 2016) and ExDEGA (E- biogen, Inc., Korea) was used for mining and graphic visualization of the data. A more detailed description of the method is provided in our previous study [4 (link)].
Free full text: Click here