Efficient nutritional assimilation and energy metabolism in the stomachs of yaks contribute to their adaption to harsh environments. Accurate gene expression profile analysis will help further reveal the molecular mechanism of nutrient and energy metabolism in the yak stomach. RT-qPCR is regarded as an accurate and dependable method for analyzing gene expression. The selection of reference genes is essential to obtain meaningful RT-qPCR results, especially in longitudinal gene expression studies of tissues and organs. Our objective was to select and validate optimal reference genes from across the transcriptome as internal controls for longitudinal gene expression studies in the yak stomach. In this study, 15 candidate reference genes (CRGs) were determined according to transcriptome sequencing (RNA-seq) results and the previous literature. The expression levels of these 15 CRGs were quantified using RT-qPCR in the yak stomach, including the rumen, reticulum, omasum and abomasum at five stages: 0 days, 20 days, 60 days, 15 months and three years old (adult). Subsequently, the expression stabilities of these 15 CRGs were evaluated via four algorithms: geNorm, NormFinder, BestKeeper and the comparative C method. Furthermore, RefFinder was employed to obtain a comprehensive ranking of the stability of CRGs. The analysis results indicate that , and are the most stable genes in the yak stomach throughout the growth cycle. In addition, to verify the reliability of the selected CRGs, the relative expression levels of were quantified via RT-qPCR using the three most stable or the three least stable CRGs. Overall, we recommend combining , and as reference genes for the normalization of RT-qPCR data in the yak stomach throughout the growth cycle.
