RNA was extracted using the Qiagen RNeasy mRNA extraction kit and sent to the Genome Analysis Core for sequencing as previously described (17 (link); 18 (link)). Libraries were prepared using the Illumina Stranded mRNA Prep and sequenced on an Illumina HiSeq 4000. Reads were aligned using STAR RNA-seq aligner to hg38 genome build. The gene and exon counts were generated by FeatureCounts. Differential expression analysis was carried out using a negative binomial generalized log-linear model in the edgeR R package. Sequencing data are available in the Gene Expression Omnibus database (GEO number GSE185957).