First, we mapped differentially expressed genes to each term of the GO database (http://www.geneontology.org/) and calculated the gene number of each term, and thus we got a GO function gene list and gene number statistics. Then, we used hypergeometric inspection, compared with the entire genome background, to find out significantly enriched GO entries in the differentially expressed genes.