In lactating dairy cattle, the corpus luteum (CL) is a dynamic endocrine tissue vital for pregnancy maintenance, fertility, and cyclicity. Understanding processes underlying luteal physiology is therefore necessary to increase reproductive efficiency in cattle. A common technique for investigating luteal physiology is reversetranscription quantitative PCR (RT-qPCR), a valuable tool for quantifying gene expression. However, reference-gene-based RT-qPCR quantification methods require utilization of stably expressed genes to accurately assess mRNA expression. Historically, selection of reference genes in cattle has relied on subjective selection of a small pool of reference genes, many of which may have significant expression variation among different tissues or physiologic states. This is particularly concerning in dynamic tissues such as the CL, with its capacity for rapid physiologic changes during luteolysis, and likely in the less characterized period of CL maintenance during pregnancy. Thus, there is a clear need to identify reference genes well suited for the bovine CL over a wide range of physiological states. Whole-transcriptome RNA sequencing stands as an effective method to identify new reference genes by enabling the assessment of the expression profile of the entire pool of mRNA transcripts. We report the identification of 13 novel putative reference genes using RNA sequencing in the bovine CL throughout early pregnancy and luteolysis: RPL4,