Expression patterns of candidate genes with important functions in animal metabolism can help to identify potential molecular markers for cattle production traits. Reverse transcription followed by polymerase chain reaction is a method for rapid and accurate mRNA quantification. However, for exact comparison of mRNA quantity in various samples or tissues, it is important to choose appropriate reference genes. In cattle, little information is available on the expression stability of housekeeping genes (HKGs). The aim of the present study is to develop a set of reference genes that can be used for normalization of concentrations of mRNAs of genes expressed in the bovine liver, kidney, pituitary and thyroid. The study was performed on 6-, 9-, and 12-month-old bulls of dairy and meat cattle breeds. Six HKGs were investigated: ACTB, GAPDH, HPRTI, SDHA, TBP, and YWHAZ. The most stably expressed potential reference HKGs differed among tissues/organs examined: ACTB, TBP, YWHAZ, GAPDH, HPRTI, and SDHA in the liver; GAPDH and YWHAZ in the kidney; GAPDH and SDHA in the pituitary; and TBP and HPRTI in the thyroid. The results showed that the use of a single gene for normalization may lead to relatively large errors, so it is important to use multiple control genes based on a survey of potential reference genes applied to representative samples from specific experimental conditions.