Enrichment analysis was performed using tools available at WEB-based GEne SeT AnaLysis Toolkit (WebGestalt) [24, 25]. Entrez gene IDs were used as unique gene identifiers. Target gene lists for main effects and interactions and a background gene set consisting of all 9530 genes were used to identify enrichment in GO, KEGG, Transcription Factor Target, Microarray Target, Protein Interaction Network Module, and Phenotype Analysis in WebGestalt using Fisher’s exact test and BH correction for multiple testing.