Honey bee has been widely used as a model insect for biological sciences because of its sociality and specialized labor division. For the investigation of the seasonal and labor-dependent expression patterns of genes putatively involved in its sociality, quantitative real-time PCR (qRT-PCR) can be applied to quantify gene expression level and selection of reliable reference gene(s) for normalization is an accurate step. In this study, using three softwares (geNorm, NormFinder and BestKeeper), we evaluated seasonal expression stabilities of four reference genes that have been widely used for qRT-PCR in forager and nurse heads. Among four candidates, two genes, rpS18 and gapdh, were suggested to be the optimal reference genes for qRT-PCR.