Honey bees are both important pollinators and model insects due to their highly developed sociality and colony management. To better understand the molecular mechanisms underlying honey bee colony management, it is important to investigate the expression of genes putatively involved in colony physiology. Although quantitative real-time PCR (qRT-PCR) can be used to quantify the relative expression of target genes, internal reference genes (which are stably expressed across different conditions) must first be identified to ensure accurate normalisation of target genes. To identify reliable reference genes in honey bee (Apis mellifera) colonies, therefore, we evaluated seven candidate genes (ACT , EIF, EF1, RPN2, RPS5, RPS18 and GAPDH) in samples collected from three honey bee tissue types (head, thorax and abdomen) across all four seasons using three analysis programmes (NormFinder, BestKeeper and geNorm). Subsequently, we validated various normalisation methods using each of the seven reference genes and a combination of multiple genes by calculating the expression of catalase (CAT). Although the genes ranked as the most stable gene were slightly different on conditions and analysis methods, our results suggest that RPS5, RPS18 and GAPDH represent optimal honey bee reference genes for target gene normalisation in qRT-PCR analysis of various honey bee tissue samples collected across seasons. The Western honey bee, Apis mellifera L., plays an important role as a pollinator 1. In addition, the honey bee is considered to be a key model insect due to its relatively complex behaviours, including sociality, labour division and colony management 2. Previous studies have demonstrated that endocrine system status and gene expression are important factors for flexible honey bee colony management, which involves colonies seasonally regulating their labour division and population dynamics 3-5. In order to extend our understanding of the molecular mechanisms that underlie the regulation of honey bee colony physiology, information on the physiological functions of the genes putatively associated with colony management can be determined by analysing their expression profiles among different seasons and honey bee tissues 6,7. In quantitative real-time PCR (qRT-PCR), gene-specific mRNA (or cDNA) is quantified; this method has been used extensively because of its relative speed, sensitivity, replicability and accuracy 8,9. Therefore, qRT-PCR would be an ideal method for analysing the expression patterns of honey bee genes putatively involved in the plasticity of colony molecular physiology in samples collected across different seasons and tissues. However, because qRT-PCR results are highly sensitive to the initial amount of RNA content in the amplification reaction, the interpretation of target gene expression levels among various conditions would result in appreciable errors without the use of a reliable internal standard 7-10. Therefore, prior to analysing target gene expression levels among conditions, reference genes are requir...