Libraries were prepared using the NeoPrep™ system according to
the TruSeq Stranded mRNA Library Prep protocol (Illumina) and sequenced using
the Illumina HiSeq 2500 or 3000. Reads were aligned using the STAR v.2.5 25, the
abundance of each gene was quantified using StringTie v.1.3.3.26 and the
differential gene expression was performed using DESeq2 v.1.20 (
Benjamini-Hochberg correction was used to adjust for multiple comparisons and a
corrected p-value (q-value) of 0.05 or less was considered statistically
significant.