Single-end deep sequencing was performed on all cDNA libraries with a
75-nucleotide (nt) read length using HiSeq-2500 Genome Analyzer (Illumina, San
Diego, CA). Adaptors and low quality regions were trimmed from raw sequences
using btrim with options “−3 -P -l 15” [35 (link)]. The trimmed sequences were mapped to human
genome (hg38) with BWA [36 (link)]. For microRNA
annotation, miRBase v21 was used [37 (link)].
Differential gene expression analysis was performed using R package
“DESeq2” [38 (link)]. MicroRNAs
with base mean expression (BME) < 10 were excluded, and differential
expression was defined as absolute log2(fold-change) > 1.0.