The raphidophyte Heterosigma akashiwo is a globally distributed harmful alga that has been associated with fish kills in coastal waters. To understand the mechanisms of H. akashiwo bloom formation, gene expression analysis is often required. To accurately characterize the expression levels of a gene of interest, proper reference genes are essential. In this study, we assessed ten of the previously reported algal candidate genes (rpL17-2, rpL23, cox2, cal, tua, tub, ef1, 18S, gapdh, and mdh) for their suitability as reference genes in this species. We used qRT-PCR to quantify the expression levels of these genes in H. akashiwo grown under different temperatures, light intensities, nutrient concentrations, and time points over a diel cycle. The expression stability of these genes was evaluated using geNorm and NormFinder algorithms. Although none of these genes exhibited invariable expression levels, cal, tub, rpL17-2 and rpL23 expression levels were the most stable across the different conditions tested. For further validation, these selected genes were used to normalize the expression levels of ribulose-1, 5-bisphosphate carboxylase/oxygenase large unite (HrbcL) over a diel cycle. Results showed that the expression of HrbcL normalized against each of these reference genes was the highest at midday and lowest at midnight, similar to the diel patterns typically documented for this gene in algae. While the validated reference genes will be useful for future gene expression studies on H. akashiwo, we expect that the procedure used in this study may be helpful to future efforts to screen reference genes for other algae.