The principles of tissue engineering (TE) are widely used for bone regeneration concepts. Three-dimensional (3D) cultivation of autologous human mesenchymal stromal cells (MSCs) on porous scaffolds is the basic prerequisite to generate newly formed bone tissue. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a specific and sensitive analytical tool for the measurement of mRNA-levels in cells or tissues. For an accurate quantification of gene expression levels, stably expressed reference genes (RGs) are essential to obtain reliable results. Since the 3D environment can affect a cell's morphology, proliferation, and gene expression profile compared with two-dimensional (2D) cultivation, there is a need to identify robust RGs for the quantification of gene expression. So far, this issue has not been adequately investigated. The aim of this study was to identify the most stably expressed RGs for gene expression analysis of 3D-cultivated human bone marrow-derived MSCs (BM-MSCs). For this, we analyzed the gene expression levels of n=31 RGs in 3D-cultivated human BM-MSCs from six different donors compared with conventional 2D cultivation using qRT-PCR. MSCs isolated from bone marrow aspirates were cultivated on human cancellous bone cube scaffolds for 14 days. Osteogenic differentiation was assessed by cell-specific alkaline phosphatase (ALP) activity and expression of osteogenic marker genes. Expression levels of potential reference and target genes were quantified using commercially available TaqMan® assays. mRNA expression stability of RGs was determined by calculating the coefficient of variation (CV) and using the algorithms of geNorm and NormFinder. Using both algorithms, we identified TATA box binding protein (TBP), transferrin receptor (p90, CD71) (TFRC), and hypoxanthine phosphoribosyltransferase 1 (HPRT1) as the most stably expressed RGs in 3D-cultivated BM-MSCs. Notably, genes that are routinely used as RGs, for example, beta actin (ACTB) and ribosomal protein L37a (RPL37A), were among the least stable genes. We recommend the combined use of TBP, TFRC, and HPRT1 for the accurate and robust normalization of qRT-PCR data of 3D-cultivated human BM-MSCs.