The analyses were performed using the topGO and kegga function from the edgeR R Bioconductor packages. The statistical metrics presented are as for S2 and S3 Tables respectively.