Functional annotation clustering of differentially expressed genes in mated queen spermatheca vs. virgin queen spermatheca using the online tool DAVID [50,51].

Recognized gene clusters were identified using the lowest stringency of the software's Functional Annotation Clustering settings. The reported genes are those that were clustered and had an FDR value < 0.05. Genes are reported using their BeeBase gene identifiers and their protein descriptions. Genes in bold had their expression additionally examined by RT-qPCR. See S4 Table for a complete list of all the genes that were differentially expressed between mated queen spermatheca and virgin queen spermatheca. (XLSX)