Libraries were prepared using the standard protocol for messenger RNA sequencing (Illumina, TruSeq Stranded mRNA) and sequenced on HiSeq2000 instrument (Illumina) according to the manufacturer’s protocol. The resulting reads were aligned to the GRCh37/hg19 assembly of the human genome using TopHat (42 (link)) with the default parameters. Finally, the gene expression level was calculated as the number of reads per gene, computed using HTSeq (43 (link)) and gene features as defined in the GRCh37.75 release of the human genome (gtf file) and expressed as RPKMs. Differential gene expression analysis was performed with the Bioconductor (44 (link)) edgeR package (45 (link)) for the R statistical Software (http://www.R-project.org/). The raw reads and processed data derived from RNA-seq experiments are available at NCBI’s Gene Expression Omnibus (46 (link),47 (link)) and are accessible through the following GEO Series accession numbers: GSE89831, GSE89838, GSE89839, GSE89840. All computations were performed using R software package (http://www.R-project.org/).