Total RNA was extracted using a GenElute Mammalian Total RNA Miniprep Kit (Sigma) according to the manufacturer’s protocol. RNA samples were treated with RNase-free DNase (Sigma), and quality control was determined using a Qubit assay and an Agilent 2100 Bioanalyzer, as per each manufacturer’s specifications. Sample amounts were normalized and 1000 ng was used for library preparation using the NEB Ultra RNA Library Preparation Kit as per the manufacturer’s instructions. Library QC and quantitation were performed on all individual libraries using the Qubit assay and the Agilent 2100 Bioanalyzer. Libraries were normalized and pooled via Qubit measurement. The final pool was quantitated via qPCR. Sequencing was performed on the Illumina HiSeq 2500 Rapid Run Mode on one flowcell (two lanes) as per the system manufacturer. Raw RNASeq data were processed, normalized, and mapped to the human reference genome (hg19) using CLC Genomics Workbench 8 (Qiagen). Differentially expressed genes were identified using DESeq2 (35 (link)) with a significance cutoff of 0.001. Hierarchical clustering was performed using Cluster 3.0/Java TreeView.