Quantitative reverse transcription polymerase chain reaction (qRT-PCR) in thyroid tumors require accurate data normalization, however, there are no sufficient studies addressing the suitable reference genes for gene expression analysis in malignant and normal thyroid tissue specimens. The purpose of this study was to identify valid internal control genes for normalization of relative qRT-PCR studies in human papillary thyroid carcinoma tissue samples. The expression characteristics of 12 candidate reference genes (GAPDH, ACTB, HPRT1, TBP, B2M, PPIA, 18SrRNA, HMBS, GUSB, PGK1, RPLP0, and PGM1) were assessed by qRT-PCR in 45 thyroid tissue samples (15 papillary thyroid carcinoma, 15 paired normal tissues and 15 multinodular goiters). These twelve candidate reference genes were selected by a systematic literature search. GeNorm, NormFinder, and BestKeeper statistical algorithms were applied to determine the most stable reference genes. The three algorithms were in agreement in identifying GUSB and HPRT1 as the most stably expressed genes in all thyroid tumors investigated. According to the NormFinder software, the pair of genes including ‘GUSB and HPRT1’ or ‘GUSB and HMBS’ or ‘GUSB and PGM1’ were the best combinations for selection of pair reference genes. The optimal number of genes required for reliable normalization of qPCR data in thyroid tissues would be three according to calculations made by GeNorm algorithm. These results suggest that GUSB and HPRT1 are promising reference genes for normalization of relative qRT-PCR studies in papillary thyroid carcinoma.