All analyses apart from QC and read mapping were performed using R version 3.6.1 [20]. An overview of the analysis workflow can be found in Fig. 1. The Benjamini–Hochberg (FDR) [21] procedure was used to correct for multiple testing. Differentially expressed genes were considered statistically significant at FDR < 0.05. All downstream analyses were performed using genes significantly differentially expressed at FDR < 0.25.
