As presented above and detailed in Table 5, 47 genes with score ≥ 2 had promising q-value improvement, and were used as our high priority list for follow up bioinformatics analysis. These genes also had evidence from at least 2 different species (Figure 1). We first performed functional enrichment analysis of Gene Ontology (GO) terms implemented in the WebGestalt tool. In this tool, each gene set is tested its functional enrichment with GO annotations based on the hypergeometric test. As shown in Table 6 the most significantly enriched functional terms belonged to the groupings of neurotransmitter receptor activity, ion binding, and synaptic structure. The most significantly enriched functional terms were behavior (pBH = 5.08 × 10-5), gamma-aminobutyric acid (GABA) signaling pathway (pBH = 8.51 × 10-5) and cell communication (pBH = 0.0001) in the GO category of "Biological Process"; GABA receptor activity (pBH = 3.35 × 10-9), GABA-A receptor activity (pBH = 7.38 × 10-8), and neurotransmitter binding (pBH = 7.38 × 10-8) in the category of "Molecular Function"; and postsynaptic membrane (pBH = 1.86 × 10-7), chloride channel complex
