For RNA-seq, libraries were statistically analyzed in R to determine differential expression between treatment conditions. Raw P-values are shown, as well as corrected P-values (false discovery rate < 0.05). The freely available Ballgown package (Pertea et al. 2016) afforded generalized linear modeling statistics for differential expression analysis and data visualization tools. Libraries were normalized with Ballgown package default – sum of each library’s log expression measurements below the 75th percentile (Frazee et al. 2015). For intersectional analyses, hypergeometric statistics were performed on pairwise comparisons and considered different from chance at P-value < 0.05.
