For the gene term enrichment analysis, we used the GOATOOLS Python package (v.1.2.3)92, using hypergeometric tests with the GO database. Like ‘Differential expression analysis: gene term enrichment and hypergeometric analysis’, we defined the gene set universe as all unique genes tested from differential expression analysis.