Technical note: Selection of suitable reference genes for studying gene expression in milk somatic cell of yak (Bos grunniens) during the lactation cycle.

Quantitative real-time PCR is the most sensitive technique for gene expression analysis. Data normalization is essential to correct for potential errors incurred in all steps from RNA isolation to PCR amplification. The commonly accepted approach for normalization is the use of reference gene. Until now, no suitable reference genes have been available for data normalization of gene expression in milk somatic cells of lactating yaks across lactation. In the present study, we evaluated the transcriptional stability of 10 candidate reference genes in milk somatic cells of lactating yak, including ACTB, B2M, GAPDH, GTP, MRPL39, PPP1R11, RPS9, RPS15, UXT, and RN18S1. Four genes, RPS9, PPP1R11, UXT, and MRPL39, were identified as being the most stable genes in milk somatic cells of lactating yak. Using the combination of RPS9, PPP1R11, UXT, and MRPL39 as reference genes, we further assessed the relative expression of 4 genes of interest in milk somatic cells of yak across lactation, including ELF5, ABCG2, SREBF2, and DGAT1. Compared with expression in colostrum, the overall transcription levels of ELF5, ABCG2, and SREBF2 in milk were found to be significantly upregulated in early, peak, and late lactation, and significantly downregulated thereafter, before the dry period. A similar pattern was observed in the relative expression of DGAT1, but no significant difference was revealed in its expression in milk from late lactation compared with colostrum. Based on these results, we suggest that the geometric mean of RPS9, PPP1R11, UXT, and MRPL39 can be used for normalization of real-time PCR data in milk somatic cells of lactating yak, if similar experiments are performed.