Gene ontology analyses were performed using the clusterProfiler package20 as implemented in R. All differentially expressed genes that passed the threshold of significance at 25% FDR were included in the analysis. Results for the enrichment analysis were extracted and plotted using the ggplot2 package in R.