The DE (differentially expressed) was determined using edgeR, with the option of ‘tagwise’ dispersion; the threshold for significance was set at a false discovery rate (FDR) < 0.05. For gene DV (differentially variable), we utilized the ‘var.test’ procedure in the R ‘stats’ package. For exon DV, we computed pairwise distances between all samples (Canberra metric) and then utilized the ‘mrpp’ function in the ‘vegan’ R package. The mrpp function is sensitive to differences in spread/dispersion of pairwise distances, as well as within/between group distance differences.
