design: Sequencing performed by GENEWIZ, Inc. via Illumina NovaSeq/HiSeq 2x150 bp. Sequenced reads were trimmed for adaptor sequence using Trim Galore! v0.6.3, and HiSat2 v2.1.0 was used to map the raw reads to the reference genome hg38. Expression for each gene were evaluated using featureCounts v1.6.4 based on annotation from GENCODE release 30. The raw counts were then normalized and differential gene expression was determined using DESeq2 v2.11.40.17