Background:Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is commonly considered as the best-established technique for gene expression assay. However, appropriate reference gene set selection remains the critical and challenging subject for the proper understanding of gene expression pattern. Mixed opinions pertain in how to choose optimal reference gene set according to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guideline. Therefore, the purpose of this study was to investigate which schemes was the most feasible for the identification of reference genes in bone and cartilage bioengineering experiments. In this study, rat bone mesenchymal stem cells (rBMSCs), skeletal muscle tissue and adipose tissue were utilized, undergoing either chondrogenic or osteogenic induction, to investigate the optimal reference gene set identification scheme ensuring the stable and accurate interpretation of gene expression in bone and cartilage bioengineering. Results:The stability and pairwise variance of eight candidate reference genes were analyzed using geNorm. The V0.15- vs. Vmin-based normalization scheme in rBMSCs had no significant effect on the eventual normalization of target genes. In terms of the muscle tissue, the results of the correlation of NF values between the V0.15 and Vmin schemes and the variance of target genes expression levels generated by these two schemes showed different schemes may have a significant effect on the eventual normalization of target genes. Three selection schemes were adopted in terms of the adipose tissue, including the three optimal reference genes (Opt3), V0.20 and Vmin schemes, and the analysis of NF values and eventual normalization of target genes showed that the different selection schemes have significant effect on the eventual normalization of target genes.Conclusions:Based on the results, the proposed cut-off value of Vn/n+1 under 0.15, according to geNorm algorithm, should be considered with caution, especially in certain tissue types such as skeletal muscle and adipose tissue. Instead, the minimum Vn/n+1 should be used as the cut-off for choosing the optimal reference gene set especially when the stability and variation of candidate reference genes in a specific study are unclear.