Accurate normalization of gene expression with qRT-PCR depends on the use of appropriate reference genes (RGs) for the species under a given set of experimental conditions. Multiple RGs for gene expression analysis of soybean exposed to heavy metal stress treatment have not been reported in the literature. In this study, we evaluated the expression stability of ten candidate RGs in leaves, roots and stems of two soybean cultivars exposed to cadmium (Cd). Based on the geNorm and NormFinder analysis, ACT3, PP2A, ELF1B and F-box were the most stable RGs in these gene expression studies. In contrast, G6PD, UBC2, TUB, and ELF1A were the most variable ones and should not be used as RGs in these experimental conditions.