Identification of valid reference genes during the differentiation of human myoblasts
- Jens Stern-Straeter†1,
- Gabriel A Bonaterra†2Email author,
- Karl Hörmann,
- Ralf Kinscherf†2 and
- Ulrich R Goessler†1
© Stern-Straeter et al; licensee BioMed Central Ltd. 2009
Received: 08 April 2009
Accepted: 02 July 2009
Published: 02 July 2009
Analysis of RNA expression using real-time PCR (qRT-PCR) traditionally includes reference genes (RG) as an internal control. This practice is being questioned as it becomes increasingly clear that RG may vary considerably under certain experimental conditions. Thus, the validity of a particular RG must be determined for each experimental setting. We used qRT-PCR to measure the levels of six RG, which have been reported in the literature to be invariant. The RG were analyzed in human myoblast cultures under differentiation conditions. We examined the expression by qRT-PCR of mRNA encoding Beta-actin (ACTB), Beta-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), peptidylprolyl isomerase A (PPIA), TATA box binding protein (TBP) and ribosomal protein, large, P0 (RPLPO). The mRNA expression of the following genes of interest (GOI) were analyzed: skeletal muscle alpha 1 actin (ACTA1), myogenin/myogenic factor 4 (MYOG), embryonic skeletal muscle myosin heavy chain 3 (MYH3) and the activity of creatine phosphokinase (CK). The geNorm, NormFinder and BestKeeper software programs were used to ascertain the most suitable RG to normalize the RNA input.
Using the geNorm program, RPLPO and TBP were found to be the most stable genes, additionally a suitable normalization factor (NF) was calculated. The NormFinder software showed that RPLPO was the most stable, whereas TBP ranked second. BestKeeper program also revealed that RPLPO and TBP as stable genes, but PPIA was the most stable reference gene, whereas GAPDH and ACTB were the worst ranked.
RNA expression analyses including three independent softwares revealed that RPLPO, TBP as reference genes or NF calculated by geNorm software, are suitable to normalize the mRNA expression in myoblast after culture under differentiation conditions. Significant correlations can be identified between the differentiations markers ACTA1, MYOG, MYH3 and creatine phosphokinase (CK) activity, when the expression is normalized with the NF calculated with RPLPO and TBP.
The real-time polymerase chain reaction (qRT-PCR) has revolutionized the field of gene expression analysis in living organisms. In comparison to classical semi-quantitative reverse transcription-PCR (sqRT-PCR), the main advantages of qRT-PCR are its higher sensitivity, specificity, and broad quantification range [1, 2]. Despite being an extremely powerful technique, qRT-PCR suffers from certain pitfalls, the most important being the normalization with a reliable reference gene, habitually called 'housekeeping gene (HKG)'. The housekeeping term was initially given to genes that are necessary for cell function and being constitutively expressed in each cell type  but obviously it would be more precise to call it 'reference gene (RG)'. RG are often taken from the literature and used across a variety of experimental conditions, some of which may induce differences in the RG own expression under certain conditions . Thus, experimental results are highly dependent on the RG gene chosen . If unrecognized, unexpected changes in RG expression could result in erroneous conclusions about real biological effects such as responses to drugs [4–6]. Such errors could be introduced at a number of stages throughout the experimental protocol (input sample, RNA extraction, reverse transcription, etc.) [1, 2, 7]. Therefore, it is critical to correct any errors between samples when measuring RNA expression. During the selection of an adequate RG, several indications should be taken into consideration: 1) the expression of the RG must remain constant throughout the intervention; 2) the amplification efficiency of the RG should be similar to that of the genes of interest (GOI); and 3) the abundance of the RG should be similar to that of the genes of interest . Therefore, appropriate RG validation of internal references is crucial in order to avoid misinterpretations of study findings . If the selected RG fluctuate randomly among samples, then subtle differences between GOI will be lost. For example, variation in presumably stable RG was shown in studies examining their expression in serum-stimulated fibroblasts [9, 10].
Skeletal muscle tissue engineering aims at the reconstruction of skeletal muscle loss. The field of skeletal muscle tissue engineering has involved to create efficient muscle tissue by using the regenerative potential of stem cells and their potential for proliferation and maturation . The preferred sources of cells for skeletal muscle tissue engineering applications are primary satellite cells . Human satellite cells can be successfully extracted and expanded in vitro and differentiated into myofibers using different stimuli. In this context, previous reports have used conventional RG , but there is no data available about the validity of these RG during serum-dependent differentiation of human myoblasts in vitro. Our main GOI were differentiation markers that have unknown RNA expression under experimental conditions. Consequently, we could anticipate that differences between the treatments might be small. Therefore, it was highly important to find an internal reference with minimal variability. Thus, we used qRT-PCR to investigate the levels of six RG expressed in cultured myoblasts. Here we report the most suitable RG to detect small variations during the differentiation of myoblasts under certain culture conditions.
Differentiation of Myoblasts
Samples indicated various stages of differentiation. We initially verified that our samples truly represented differentiating cells by analyzing the creatine phosphokinase (CK) activity level, a well-known and accepted marker of differentiation .
All RNA samples were examined as to their concentration, purity, and integrity. RNA purity was measured using the NanoDrop® Spectrophotometer. Based on the absorbance ratio at 260 nm/280 nm (mean ± standard deviation [SD], 1.90 ± 0.05), all RNA samples were pure and protein free. The mean (± SEM) A260/230 ratio of 1.90 ± 0.20 (range from 1.65 – 2.10) indicated that the RNA was free of phenol and ethanol. RNA integrity was assessed by the calculation of RIN values using the Agilent 2100 Bioanalyzer. Differentiated and undifferentiated myoblast samples revealed RIN values between 6.0 and 8.0.
Stability of RG within Sample Groups
The stability of the six RG was assessed for each sample group. Gene expression levels were measured by qRT-PCR and the expression stabilities were evaluated by the three most commonly used software-based methods: 1- geNorm , 2- NormFinder  and 3-BestKeeper .
Expression Stability of RG Evaluated by NormFinder Software
N and D
RPLPO and TBP
3- BestKeeper (BK)
Expression Stability of RG Evaluated by BestKeeper Software
n = 6
n = 3
SD [± CP]
CV [% CP]
Expression Stability of GOI Evaluated by BestKeeper Software
SD [± CP]
CV [% CP]
SD [± x-fold]
Evaluation of Selected Candidate RG and Normalization Approach
Correlation Between CK Activity, RG and Differentiation Markers
Myogenesis of skeletal muscle is a highly complex phenomenon controlled and regulated by multiple intracellular signalling mechanisms . In order to investigate optimal culture conditions and characterization of human myoblasts, we analyzed cell behaviour under different culture media and showed that cultures treated with differentiation medium expressed several characteristic features of mature skeletal muscle . Furthermore, in myoblasts cultured with growth medium, certain markers of differentiation were detected, demonstrating that, for an accurate characterization, multiple markers (e.g. genes, enzymes and structure proteins) must be analyzed . During skeletal muscle differentiation, mononucleated proliferating myoblasts stop dividing, activate muscle-specific gene expression, and fuse into multinucleated myotubes ; these events can produce drastic or small variation in the expression of genes and proteins. Markers of differentiation were also detected in cultures treated with non-differentiation medium, but there was no formation of myotubes. In the enzymatic assay of CK, a well-known and accepted marker of differentiation, cultures treated with differentiation medium showed a higher activity, evidencing a higher degree of differentiation . In this context, the strategy used to normalize the mRNA with adequate RG is critical. To identify the best RG as normalisers in gene expression studies, several strategies including computer programs have been recommended [5, 6, 9, 14–20]. The absence of the differential expression or variability of the candidate RG examined is the strongest proof of suitability and stability . Therefore, we proposed to use several software programs to analyze six RG and select the RG or normalization factors to correct the RNA input. The major result of our study – evaluating six candidate RG in myoblasts under normal and differentiation conditions, during 8 and 12 days of treatment – was that the genes RPLPO or TBP were identified by the three independent software programs as suitable RG that did not differ in their expression in differentiated and undifferentiated myoblasts. Consequently, only RPLPO, identified by geNorm, NormFinder and BestKeeper software programs fulfils the criterion of expression stability between samples and can be recommended as an accurate normalizer for relative gene quantification in cultured myoblasts. First, the expression of the candidate RG between the respective conditions was compared using an adequate test such as geNorm, NormFinder and BestKeeper. It can be assumed that genes with significantly different expressions are not suited for target gene normalization, since they are affected by the study condition of interest. Thus, those genes should be excluded as normalizers. In the present study, the expression of all genes in the differentiated and undifferentiated myoblast samples showed different grades of variability but the most stable were RPLPO and TBP. We used the BestKeeper software to confirm the results obtained by geNorm and NormFinder. The raw CPs analyzed by BestKeeper seem to be good estimators of the expression levels as they are (in most cases) normally distributed and a parametric test can thus be performed . This also gives the CP datasets the Gaussian distribution, justifying usage of parametric methods . Low numbers of expressed genes where CPs were obtained somewhere surely show different variances compared to highly expressed genes with CPs. These higher levels of CV and SD invalidate the use of BestKeeper as a reference to normalizing the expression of GOI . The correlation analysis of the normalized gene with NF showed a bad correlation between the RG RPLPO, TBP or GAPDH and the different differentiation markers (CK, ACTA1, MYOG and MYH3). In this context, the absence of correlation confirmed that the RG selected by geNorm, NormFinder and BestKeeper are suitable as reference for the normalization. In addition, we found a significant correlation between B2M or ACTB and the GOI which indicates a variation associated with the differentiation status of the myoblast. Thus, they cannot be used as reference genes.
Two or three genes represent a realistic calculation basis in a common laboratory and the minimum necessary number for a good analysis . RPLPO STD dev CP [± SD] is identical to BestKeeper calculated with RPLPO, TBP and PPIA. In this context, the RG, RPLPO and TBP, detected by BestKeeper agree with the best RG detected by geNorm and NormFinder.
The normalization with NF, RPLPO or TBP allows measurements of differences in the expression of the differentiation gene-markers between myoblasts treated with differentiation medium or normal medium (control).
In conclusion, the RG RPLPO, TBP as well as a NF calculated by geNorm software are recommended as references for relative gene quantification in gene profiling studies of single genes in myoblasts under serum-differentiation conditions. The confirmation of RG by the three software programs is a suitable method to perform an adequate normalization of the mRNA input in quantitative PCR experiments.
Human skeletal muscle biopsies were obtained from 15 patients during head and neck surgery. The median age of these patients was 58, ranging from 41 to 72. The study was approved by the Ethics Committee of the Mannheim Faculty of Medicine, University of Heidelberg, Germany, and the patients gave informed consent. Satellite cells were dissociated from the minced muscles by digestion with collagenase B (Roche, Mannheim, Germany) for 60 min and 0.05% trypsin-0.02% EDTA (Promo Cell, Heidelberg, Germany) for 45 min at 37°C, filtered through a sterile 70-μm cell strainer (Becton Dickinson, Franklin Lakes, NJ, USA) and purified with the pre-plating technique as recently described . Purity of myoblast cultures (>80%) was evaluated by anti-desmin immunostainings. Cells were grown on 0.2% gelatine-coated culture flasks (Sigma, Deisenhofen, Germany) and in Ham's F10 growth medium, containing 1% penicillin/streptomycin/fungizone-solution (PSF), 2 mM L-glutamine (all from PromoCell, Heidelberg, Germany) and 10% foetal bovine serum (FBS) (PAA Laboratories, Linz, Austria) (GM). Cells were maintained at 37°C in a humidified atmosphere of 5% CO2 and 95% air, and the medium was changed every 72 h.
Differentiation of skeletal muscle myoblasts
Skeletal myoblasts were cultured to ~60% confluence, washed with phosphate-buffered saline (PBS) and induced to differentiate by a change in medium (D) consisting of minimal essential medium (MEM) (PromoCell) supplemented with 2% horse serum (PAA Laboratories), 2 mM L-glutamine and PSF (both from PromoCell).
Total RNA was isolated with the RNA Mini Kit (Qiagen, Hilden, Germany), according to the manufacturer's instructions. RNA concentration, purity, and integrity was determined by A260 and A280 (A260/A280 = 1.7–2.0) measurements using a NanoDrop 8000 Spectrophotometer (Thermo Scientific, Schwerte, Germany) and Agilent 2100 Bioanalyzer (Agilent Technologies, Waldbronn, Germany).
cDNA Synthesis and Real-time PCR
An aliquot of 0.5 μg total RNA was treated with 1 unit DNAse (Fermentas, St. Leon-Rot, Germany) 30 min at 37°C. Reverse transcription of RNA (0.5 μg) was performed with oligo (dT)12–18 primer and 200 units of SUPERSCRIPT II (Invitrogen, Karlsruhe, Germany) and 24 units of Ribo LockTM RNAse inhibitor (Fermentas) for 1 h at 42°C. The cDNA was used for PCR analysis. All cDNA probes were analyzed for: ACTA1, (NM_001100), amplicon length 85 bp; MYOG, (NM_002479), amplicon length 113 bp; MYH3, (NM_002470), amplicon length 84 bp and the RG: ACTB (NM_001101), amplicon length 104 bp; B2M, (NM_004048), amplicon length 98 bp; GAPDH, (NM_002046), amplicon length 119 bp; cyclophilin A/PPIA, (NM_203430), amplicon length 121 bp; RPLPO, (NM_001002), amplicon length 170 bp; TBP, (NM_003194), amplicon length 132 bp. The QuantiTect/PrimerAssays were purchased from QIAGEN GmbH (Hilden, Germany). cDNAs were amplified with Brilliant® II SYBR® Green QRT-PCR Master Mix (Stratagene-Agilent Technologies, Waldbronn, Germany). The thermal profile consisted of 1 cycle at 50°C for 2 minutes followed by 1 cycle at 95°C (2 min), 45 cycles at 95°C (15 sec), 60°C (1 min). Amplification was performed using the Mx3005P™ QPCR System (Stratagene). For relative quantification, a standard curve was generated in every individual run. Shortly, total RNA was pooled from muscle biopsies of healthy human volunteers, reverse transcription was performed and a serial dilution of the cDNA was used to perform the calibration curve. The data were analyzed using the relative standard curve method. For each unknown sample, the relative amount is calculated using linear regression analysis from their respective standard curves. Data were analyzed using the Mx3005P analysis software (Stratagene-Agilent Technologies, Waldbronn, Germany). The efficiencies of all GOI and RG were calculated in every individual run (Table 1).
Creatine Phosphokinase (CK) Analysis
Differentiation of myoblasts was identified by measuring the creatine phosphokinase (CK) activity of developing myofibers using the CK assay (Sigma-Aldrich Chemie, Taufkirchen, Germany). Cells were lysed with a buffer containing NP-40, protease inhibitors and were stored at -80°C. Samples were assayed according to the manufacturer's protocol. Enzymatic activity was normalized over the total protein content determined by the Lowry protein assay (Bio-Rad, Munich, Germany).
Analysis of Expression Stability
For stability comparisons of candidate RG within sample groups, the software geNorm, version 3.4  (Visual Basic application tool for Microsoft Excel), NormFinder  and BestKeeper , were used according to developer's recommendations.
GeNorm uses a gene-stability measure M, which is defined as the average pair wise variation between a particular gene and all other control genes. It calculates the optimal number of genes necessary for normalization of a target gene and combines them into a normalization factor (NF). The NormFinder program uses a model-based approach for estimation of expression and enables the identification of the single best genes as well as giving a ranking order. The BestKeeper software determines the best suited reference genes and combines them into an index (BK). The index can be compared with further target genes to decide whether they are differentially expressed under certain conditions. In addition, the option in NormFinder to define groups was applied to compare the effect of the treatment of myoblasts with differentiation medium.
The software Sigma Plot was used to carry out statistical analysis, by paired Student's t-test or Mann-Whitney-U test, as well as Pearson's product moment correlation test of gene expression among experimental groups. Gene expression results were expressed as the mean. P < 0.05 was considered statistically significant.
The authors thank Mrs. Ulrike Traut and Petra Prohaska for the great technical assistance and Michael Collins for his help with the manuscript.
- Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. BioTechniques. 2004, 37: 112-119.PubMedGoogle Scholar
- Hendriks-Balk MC, Michel MC, Alewijnse AE: Pitfalls in the normalization of real-time polymerase chain reaction data. Basic Res Cardiol. 2007, 102: 195-197. 10.1007/s00395-007-0649-0PubMed CentralView ArticlePubMedGoogle Scholar
- Rhinn H, Leroux CM, Croci N, Plotkine M, Scherman D, Escriou V: Housekeeping while brain's storming Validation of normalizing factors for gene expression studies in a murine model of traumatic brain injury. BMC Molecular Biology. 2008, 9: 62- 10.1186/1471-2199-9-62PubMed CentralView ArticlePubMedGoogle Scholar
- Dheda K, Huggett JF, Chang JS, Kima LU, Bustin SA, Johnson MA, Rook GAW, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Analytical Biochemistry. 2005, 344: 141-143. 10.1016/j.ab.2005.05.022View ArticlePubMedGoogle Scholar
- Szabo A, Perou CM, Karaca M, Perreard L, Quackenbush JF, Bernard PS: Statistical modeling for selecting housekeeper genes. Genome Biology. 2004, 5: R59- 10.1186/gb-2004-5-8-r59PubMed CentralView ArticlePubMedGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190View ArticlePubMedGoogle Scholar
- Ståhlberg A, Håkansson J, Xian X, Semb H, Kubista M: Properties of the Reverse Transcription Reaction in mRNA Quantification. Clinical Chemistry. 2004, 50: 509-515. 10.1373/clinchem.2003.026161View ArticlePubMedGoogle Scholar
- Mahoney DJ, Carey K, Fu M-H, Snow R, Smith DC, Parise G, Tarnopolsky MA: Real-time RT-PCR analysis of housekeeping genes in human skeletal muscle following acute exercise. Physiol Genomics. 2004, 18: 226-231. 10.1152/physiolgenomics.00067.2004View ArticlePubMedGoogle Scholar
- Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol. 2000, 25: 169-193. 10.1677/jme.0.0250169View ArticlePubMedGoogle Scholar
- Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods. 2000, 46: 69-81. 10.1016/S0165-022X(00)00129-9View ArticlePubMedGoogle Scholar
- Stern-Straeter J, Bran G, Riedel F, Sauter A, Hörmann K, Goessler UR: Characterization of human myoblast cultures for tissue engineering. Int J Mol Med. 2008, 21: 49-56.PubMedGoogle Scholar
- Machida S, Spangenburg EE, Booth FW: Primary rat muscle progenitor cells have decreased proliferation and myotube formation during passages. Cell Prolif. 2004, 37: 267-277. 10.1111/j.1365-2184.2004.00311.xView ArticlePubMedGoogle Scholar
- 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: research 0034.1-0034.11. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- 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: 5245-5250. 10.1158/0008-5472.CAN-04-0496View ArticlePubMedGoogle Scholar
- 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. Biotechnology Letters. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47View ArticlePubMedGoogle Scholar
- Brand-Saberi B: Genetic and epigenetic control of skeletal muscle development. Ann Anat. 2005, 187: 199-207. 10.1016/j.aanat.2004.12.018View ArticlePubMedGoogle Scholar
- Haller F, Kulle B, Schwager S, Gunawan B, Sultmann von HA, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Anal Biochem. 2004, 335: 1-9. 10.1016/j.ab.2004.08.024View ArticlePubMedGoogle Scholar
- 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: 509-515. 10.1023/B:BILE.0000019559.84305.47View ArticlePubMedGoogle Scholar
- Silver N, Best S, Jiang J, Thein SL: Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol Biol. 2006, 7: 33- 10.1186/1471-2199-7-33PubMed CentralView ArticlePubMedGoogle Scholar
- Jung M, Ramankulov A, Roiga J, Johannsen M, Ringsdorf MG, Kristiansen G, Jung K: In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR. BMC Molecular Biology. 2007, 8: 47- 10.1186/1471-2199-8-47PubMed CentralView ArticlePubMedGoogle Scholar