- Research note
- Open Access
Selection and validation of reference gene for RT-qPCR studies in co-culture system of mouse cementoblasts and periodontal ligament cells
BMC Research Notes volume 15, Article number: 57 (2022)
RT-qPCR is a reliable method for gene expression analysis, but the accuracy of the quantitative data depends on the appropriate selection of reference genes. A Co-culture system consisting of periodontal ligament cells (SV-PDL) and cementoblasts (OCCM-30) to investigate the crosstalk between these two cell lines under orthodontic condition is essential for experimental orthodontic setups in-vitro. Therefore, we aimed to identify a set of reliable reference genes suitable for RT-qPCR studies for prospective co-culture systems of OCCM-30 and SV-PDL cells.
The results demonstrated that PPIB, GUSB and RPLP0 turned out to be the three most stable reference genes for OCCM-30 in the co-culture system, while PPIB, POLR2A and RPLP0 have the three highest rankings for SV-PDL cells in the co-culture system. The most stable gene combination were PPIB and POLR2A in the co-culture system. In conclusion, PPIB is overall the most stably expressed reference gene for OCCM-30 or SV-PDL cell line in the system. The combination of PPIB and POLR2A as reference genes are indicated to be the potential and mandatory to obtain accurate quantification results for normalizing RT-qPCR data in genes of interest expression in these two cell lines co-culture systems.
Cementoblasts are located on the cementum covered root surface and have the lifelong capability to produce cementum . Periodontal ligament cells are fibroblast-like cells characterized by collagen production .
Reverse Transcriptase-quantitative polymerase chain reaction (RT-qPCR)  is a versatile molecular technique for quantification of the expression of genes of interest due to the method’s merits concerning its high sensitivity, simplicity and specificity as well as its accuracy [4, 5]. Reference genes are considered to be consistently expressed in various tissues and treatments , which guarantee precise gene expression quantification by accurate and valid data normalization .
Our laboratory will deliver a co-culture system of periodontal ligament cells with cementoblasts which in-vitro is used to mimic the biological conditions to explore the interaction between these two-cell lines [8,9,10,11]. Therefore, a commonly stable reference gene selection is of vital role for the co-culture system RT-qPCR experiment set.
Together from previous studies [12,13,14,15], the ribosomal 60S protein L22 (RPL22), Peptidylprolyl isomerase B (PPIB), polymerase RNA II polypeptide A (POLR2A), Ribosomal protein, large, P0 (RPLP0), glucuronidase, beta (GUSB), Actin-beta (β-actin), TATA-binding protein (TBP), ubiquitin C (UBC), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Tyrosine 3-monooxygenase/tryptophan 5-mono-oxygenase activate protein, zeta (YWHAZ), eukaryotic translation elongation factor 1 alpha 1 (EEF1A1), Ribosomal protein L18 (RPL) and beta-5 class I (TUBB) were identified as reliable reference genes.
Since co-culture system could be a promising method to analyze the biological effect response to mimic orthodontically induced tooth movement in-vitro. The identification of the most reliable reference genes for RT-qPCR analysis on OCCM-30 and SV-PDL, is an essential step to facilitate further research in this area.
Materials and methods
Mono cell culture
Immortalized murine mouse cementoblast (OCCM-30) cell line  and immortalized mouse murine periodontal ligament (SV-PDL) cells  were provided by Prof. Martha J. Somerman (Laboratory of Oral Connective Tissue Biology, NIH, Bethesda, USA). The optimal density of OCCM-30  and SV-PDL cells  were previously determined, thus 60—100% confluence status of cells were used for cell confluence experiments .
Both cells were cultured at a density of 1*106 cell/well in D-MEM (31885-023, Gibco) containing 10% FBS (10270-106, Gibco), 1% Penicillin/Streptomycin (15140-122, Gibco) and 1% HEPES (15630-056, Gibco) in a humidified atmosphere of 5% CO2 at 37 °C.
Direct cell–cell contact culture
The direct cell–cell contact system was established by seeding each cell line by same number of 1*106 cells/well in the same 6-well plate. After incubation overnight to allow for firm adherence to the bottom, the cell–cell contact system was established . For the experiment, control group (as 0 h) was set when cells reached approximately 60% confluence. The mRNA was harvested at 0, 12 and 24 h at the same day.
The OCCM-30/SV-PDL co-culture system was established through 6-well plate and ThinCert@ Cell Culture Inserts (pore size 0.4 µm, porosity/transparent membrane) (657,641, Greiner Bio-One) enabling the cells to exchange soluble factors  as previous described . Briefly, the SV-PDL cells (1*106 cell/well) were seeded into 6-well plate and OCCM-30 cells (1*106 cell/well) were seeded into the ThinCert@ inserts. After 6 h cultivation to allow for adherence, the inserts containing OCCM-30 cells are placed on top in the 6-well plate containing SV-PDL cells (Fig. 2A) on bottom. Then, the co-culture system was established after co-incubation for another 10 h. For the experiment, control group (as 0 h) was set when cells reached approximately 60% confluence. The mRNA was harvested at 0, 12, 24 h at the same day.
Cells were harvested with 350 μL buffer RLT (Qiagen, Germany). Afterwards, RNA was isolated with RNase Mini Kit (Qiagen, Germany) following an on-column DNA digestion (RNase-Free DNase, Qiagen, Germany) including DNase step for removal of genomic DNA. After isolation, the eluted RNA purity and quantity of each sample was verified photometrically by OD readings of the A260/280 nm ratio (Nanodrop 2000, Thermo Fisher Scientific, USA).
For every RT-qPCR 20 μL volume reaction, we used 8 μL DNase-free water (Sigma-Aldrich), 10 μL SsoAdvanced™ Universal SYBR Green Supermix (Bio-Rad), 1.0 μL cDNA and 1.0 μL primer . Primers are designed by Bio-Rad (Additional file 1: Table S1). For analysis, a Cq cut-off of 40 was applied.
Stability assessment and statistical analysis
Statistical analysis regarding reference gens stability was performed by using four different mathematical procedures: geNorm (qBase + , Biogazelle) , NormFinder (version 0.953) , BestKeeper (version 1)  and Comparative ΔCq method . The values of cycle threshold (Cq) were inputted to all programs (Additional file 2: Table S2). For evaluation, the selected reference genes were listed based on their stability values (geNorm: M value; NormFinder: Vn/Vn+1; BestKeeper: Pearson’s r value; and comparative ΔCq: mean of SD value). Graphic were produced by GraphPad software (version 8.0). Descriptive statistics are shown as arithmetic mean values ± standard deviation (SD). The ranking sum for each gene is calculated by the summation of four respective rankings (Table 1).
Expression levels of candidate reference genes in mono cell culture
For OCCM-30 cells, as showed in Fig. 1A, the β-actin, GAPDH, EEF1A1 and RPLP0 are the candidate reference genes most abundantly expressed with Cq values below 2. The genes TBP, RPL22, PPIB, YWHAZ, POLR2A, GUSB, UBC and RPL are all moderately expressed with Cq values ranging from 20 to 30. Due to the cut-off applied, 10 out of 12 measurements for TUBB (Cq values 37.17) in the mono-cultured OCCM-30 dataset are removed from Fig. 1A.
For SV-PDL cells, Fig. 1B shows that the β-actin, GAPDH, EEF1A1 and RPLP0 are the candidate reference genes most abundantly expressed with Cq values below 20. The genes TBP, RPL22, PPIB, YWHAZ, POLR2A, TUBB, GUSB, UBC and RPL are all moderately expressed with Cq values ranging from 20 to 30 in Fig. 1B.
Stability analysis of candidate reference genes in mono cell culture system
For studies with monocultured cementoblasts, total ranking results in Table 1 show that RPL22 is the least regulated reference gene in the preselected panel on OCCM-30 cells. Similarly, although PPIB is not as stable as RPL22, it ranks higher than POLR2A in all calculation in cementoblasts as showed in Fig. 1C.
It is revealed that GUSB reached the best stability values on monocultured SV-PDL cells. The geNorm and the ΔCq method in Fig. 1E show the same results. In this case, GUSB ranked highest in the comparison, but GAPDH was less stable compared to RPLP0.
The geNorm analysis revealed that the use of two reference genes in this case GUSB and GAPDH for normalization in RT-qPCR is adequate for studies in monoculture of OCCM-30 cells (Fig. 1D) and SV-PDL cells (Fig. 1F). Notably, the output results of geNorm in the selection study showed that the M values of TUBB and RPL in OCCM-30 cells and RPL in SV-PDL cells were missing, indicating their exclusion for further analysis. The ranking order and the stability values calculated with the geNorm and NormFinder programs did not change when TUBB was excluded from the dataset.
In concordance with the above given results, PPIB, GUSB, RPLP0, POLR2A and RPL22 were selected as the most five stable reference genes based on ranking sum for further analysis in the direct cell–cell contact and co-culture system.
Stability analysis of 5 chosen reference genes in direct cell – cell contact culture and co - culture system
According to the single cell-culture results, the three highest ranking genes were selected for each cell line, respectively. These were GUSB, POLR2A, RPLP0, RPL22 and PPIB. The entire ranking shows that GUSB, POLR2A and RPLP0 were the least regulated reference genes when both cell lines were cultivated with direct contact. The suitable number of reference targets in the direct cell–cell contact experimental situation was 2. As such, the suitable normalization factor can be calculated as the geometric mean  of reference targets GUSB and POLR2A for the direct cell–cell contact culture system (Table 1, Fig. 2B, C).
As analyzed using different incubation times during the same day in the indirect co-culture systems, results showed that across three different time period, PPIB, GUSB and RPLP0 have the best stability values for co-cultured OCCM-30 cells. The PPIB, POLR2A and RPLP0 reached the best stability values for co-cultured SV-PDL cells. geNorm analysis shows that the suitable number of reference targets in this experimental situation was 2 and can be calculated as the geometric mean of reference targets PPIB and POLR2A (Table 1; Fig. 2D–G). Altogether, PPIB was the most stable reference gene for the co-culture system.
The in-vitro co-culture model is based on the location of OCCM-30 and SV-PDL cells which was subjected to mimicked specific conditions during OTM. The porosity membrane allows the cells to exchange soluble factors in the co-culture setup (Fig. 2A). Thus, the co-culture system seems reasonable to investigate the intercellular communication between these two cell lines of the periodontal compartment which become closer when orthodontic force is applied [10, 25]. However, reference gene selection is depending on the exact research questions and thus applies to both the experimental condition and the corresponding control, in this sense, the suitable reference genes in hypoxic- or orthodontic force-induced conditions need to be further investigated .
In the present work, from our analysis PPIB achieved the most stable results in all four algorithms methods for the use as a reference gene according to the comparison of all potential reference genes in co-culture system at different time point in the same day. This matches with previous research’s publication analyzing the combined dental, periodontal and alveolar bone tissue of rat which showed that PPIB and YWHAZ were the most stabile reference genes for RT-qPCR analysis in untreated rats with additional periodontitis . PPIB is also reported to have the highest expression stability values and reliability on hPDLF subjected to static mechanical strain . Besides PPIB, GUSB and RPLP0 were ranked as the most stable reference genes for co-cultured OCCM-30 cells at different time periods. This is in concordance with the ranking in the control group, indicating that these three highest ranking genes are stably for the co-cultured OCCM-30 cells. Similarly, PPIB, POLR2A and RPLP0 were recommended as the three most stable reference genes for co-cultured SV-PDL cells at different time points within the same day.
Direct cell–cell contact culture is more closely to the in-vivo condition that both cell types are cultivated together. However, it would be difficult to compare the gene expression in the direct cell–cell contact culture compared to the mono-cultured cells. Thus, for the purpose, discrimination of two types of cells using special surface markers assessed by flow cytometry would be more accurate to provide separate results of the different reference genes. Besides, for the indirect co-culture system, the Cq values expression of the reference genes may be changed dependent upon the specific placement of cells within the experimental setup. Although this study did not specifically refer to this issue, one might speculate that each cell type exerts different on the other when the position of two types of cells changed as shown in Fig. 2A. Furthermore, in the present proposed setup, the cells of two types in the insert and on the bottom of the wells are < 1 mm apart. Therefore, it might be necessary to investigate if the magnitude of gravity interferes with diffusion.
We concluded that PPIB, GUSB and RPLP0 are the most stable reference genes for normalization in RT-qPCR studies using OCCM-30 cells in a co-culture system. PPIB, POLR2A and RPLP0 were demonstrated to be the most reliable normalizers for SV-PDL cells used for RT-qPCR gene expression analysis in the co-culture system. The PPIB is an ideal reference and combination of PPIB and POLR2A for RT-qPCR experiments can improve the normalization in co-culture systems.
In direct cell–cell contact culture, the gene expression results should be considered as a mean of both cell types, thus it’s necessary to separate these two types of cells using FACS by their special membrane marker to provide separate results of reference genes. Different orthodontic induced conditions such as hypoxia and mechanical forces are not included in the stimulation, which needs to be further investigated.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
- Cq :
Fluorescence-activated cell sorting
Human periodontal ligament fibroblasts
Orthodontic tooth movement
Reverse transcriptase-quantitative polymerase chain reaction
Yong J, von Bremen J, Ruiz-Heiland G, Ruf S. Adiponectin interacts in-vitro with cementoblasts influencing cell migration, proliferation and cementogenesis partly through the MAPK signaling pathway. Front Pharmacol. 2020;11:585346.
Jonsson D, Nebel D, Bratthall G, Nilsson BO. The human periodontal ligament cell: a fibroblast-like cell acting as an immune cell. J Periodontal Res. 2011;46(2):153–7.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.
Adamski MG, Gumann P, Baird AE. A method for quantitative analysis of standard and high-throughput qPCR expression data based on input sample quantity. PLoS ONE. 2014;9(8):e103917.
Bustin SA. Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002;29(1):23–39.
Ma R, Xu S, Zhao Y, Xia B, Wang R. Selection and validation of appropriate reference genes for quantitative real-time PCR analysis of gene expression in Lycoris aurea. Front Plant Sci. 2016;7:536.
Pfaffl MW, Hageleit M. Validities of mRNA quantification using recombinant RNA and recombinant DNA external calibration curves in real-time RT-PCR. Biotech Lett. 2001;23(4):275–82.
Kats A, Gerasimcik N, Nareoja T, Nederberg J, Grenlov S, Lagnohed E, Desai S, Andersson G, Yucel-Lindberg T. Aminothiazoles inhibit osteoclastogenesis and PGE2 production in LPS-stimulated co-cultures of periodontal ligament and RAW 264.7 cells, and RANKL-mediated osteoclastogenesis and bone resorption in PBMCs. J Cell Mol Med. 2019;23(2):1152–63.
Shi J, Baumert U, Folwaczny M, Wichelhaus A. Influence of static forces on the expression of selected parameters of inflammation in periodontal ligament cells and alveolar bone cells in a co-culture in vitro model. Clin Oral Investig. 2019;23(6):2617–28.
Shi J, Folwaczny M, Wichelhaus A, Baumert U. Differences in RUNX2 and P2RX7 gene expression between mono- and coculture of human periodontal ligament cells and human osteoblasts under compressive force application. Orthod Craniofac Res. 2019;22(3):168–76.
Steller D, Scheibert A, Sturmheit T, Hakim SG. Establishment and validation of an in vitro co-culture model for oral cell lines using human PBMC-derived osteoclasts, osteoblasts, fibroblasts and keratinocytes. Sci Rep. 2020;10(1):16861.
Yong J, von Bremen J, Groeger S, Ruiz-Heiland G, Ruf S. Hypoxia-inducible factor 1-alpha acts as a bridge factor for crosstalk between ERK1/2 and caspases in hypoxia-induced apoptosis of cementoblasts. J Cell Mol Med. 2021;25:9710–23.
Yong J, von Bremen J, Ruiz-Heiland G, Ruf S. Adiponectin as well as compressive forces regulate in vitro beta-catenin expression on cementoblasts via mitogen-activated protein kinase signaling activation. Front Cell Dev Biol. 2021;9:645005.
Kirschneck C, Batschkus S, Proff P, Kostler J, Spanier G, Schroder A. Valid gene expression normalization by RT-qPCR in studies on hPDL fibroblasts with focus on orthodontic tooth movement and periodontitis. Sci Rep. 2017;7(1):14751.
Niederau C, Craveiro RB, Azraq I, Brockhaus J, Bastian A, Kirschneck C, Wolf M. Selection and validation of reference genes by RT-qPCR for murine cementoblasts in mechanical loading experiments simulating orthodontic forces in vitro. Sci Rep. 2020;10(1):10893.
D’Errico JA, Berry JE, Ouyang H, Strayhorn CL, Windle JJ, Somerman MJ. Employing a transgenic animal model to obtain cementoblasts in vitro. J Periodontol. 2000;71(1):63–72.
Hakki SS, Wang D, Franceschi RT, Somerman MJ. Bone sialoprotein gene transfer to periodontal ligament cells may not be sufficient to promote mineralization in vitro or in vivo. J Periodontol. 2006;77(2):167–73.
Mercado KP, Helguera M, Hocking DC, Dalecki D. Estimating cell concentration in three-dimensional engineered tissues using high frequency quantitative ultrasound. Ann Biomed Eng. 2014;42(6):1292–304.
Bloemen V, Schoenmaker T, de Vries TJ, Everts V. Direct cell-cell contact between periodontal ligament fibroblasts and osteoclast precursors synergistically increases the expression of genes related to osteoclastogenesis. J Cell Physiol. 2010;222(3):565–73.
Bustin SA, Beaulieu JF, Huggett J, Jaggi R, Kibenge FS, Olsvik PA, Penning LC, Toegel S. MIQE precis: practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Mol Biol. 2010;11:74.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):RESEARCH0034.
Andersen CL, Jensen JL, Orntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64(15):5245–50.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26(6):509–15.
Silver N, Cotroneo E, Proctor G, Osailan S, Paterson KL, Carpenter GH. Selection of housekeeping genes for gene expression studies in the adult rat submandibular gland under normal, inflamed, atrophic and regenerative states. BMC Mol Biol. 2008;9:64.
Yong J, Elisabeth Groeger S, Ruf S, Ruiz-Heiland G. Influence of leptin and compression in GAS-6 mediated homeostasis of periodontal ligament cell. Oral Dis. 2021. https://doi.org/10.1111/odi.14092.
Kirschneck C, Proff P, Fanghanel J, Wolf M, Roldan JC, Romer P. Reference genes for valid gene expression studies on rat dental, periodontal and alveolar bone tissue by means of RT-qPCR with a focus on orthodontic tooth movement and periodontitis. Ann Anat. 2016;204:93–105.
Nazet U, Schroder A, Spanier G, Wolf M, Proff P, Kirschneck C. Simplified method for applying static isotropic tensile strain in cell culture experiments with identification of valid RT-qPCR reference genes for PDL fibroblasts. Eur J Orthod. 2020;42(4):359–70.
We sincerely thank Prof. J. Deschner and Dr. M. Nokhbehsaim (Department of Periodontology, University of Bonn, Germany) to facilitate the sending of OCCM-30 cells and Prof. M. Somerman, (Laboratory of Oral Connective Tissue Biology, NIH, Bethesda, MD, USA) to kindly provide the SV-PDL cells. We also thank the China Scholarship Council (CSC) for Ph.D. financial support.
Open Access funding enabled and organized by Projekt DEAL.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yong, J., Groeger, S., Ruiz-Heiland, G. et al. Selection and validation of reference gene for RT-qPCR studies in co-culture system of mouse cementoblasts and periodontal ligament cells. BMC Res Notes 15, 57 (2022). https://doi.org/10.1186/s13104-022-05948-x
- Reference genes
- Periodontal ligament cells