Adapters were trimmed by TrimGalore. The alignment was performed with the STAR aligner (v.2.5.2a) against the hg19 human genome, and gene features were quantified using the GENCODE v.19 database. Raw counts were normalized into CPM values by edgeR in Bioconductor, and differential expression was calculated using a negative binomial model, with an FDR cutoff <0.05 used to define statistical significance.
