- Research article
- Open Access
Effect of experimental treatment on GAPDH mRNA expression as a housekeeping gene in human diploid fibroblasts
BMC Molecular Biology volume 11, Article number: 59 (2010)
Several genes have been used as housekeeping genes and choosing an appropriate reference gene is important for accurate quantitative RNA expression in real time RT-PCR technique. The expression levels of reference genes should remain constant between the cells of different tissues and under different experimental conditions. The purpose of this study was to determine the effect of different experimental treatments on the expression of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) mRNA so that the reliability of GAPDH as reference gene for quantitative real time RT-PCR in human diploid fibroblasts (HDFs) can be validated. HDFs in 4 different treatment groups viz; young (passage 4), senescent (passage 30), H2O2-induced oxidative stress and γ-tocotrienol (GTT)-treated groups were harvested for total RNA extraction. Total RNA concentration and purity were determined prior to GAPDH mRNA quantification. Standard curve of GAPDH expression in serial diluted total RNA, melting curve analysis and agarose gel electrophoresis were used to determine the reliability of GAPDH as reference gene.
HDFs with different experimental treatments exhibited diverse cell morphology with different expression of senescence-associated β-galactosidase (SA β-gal) activity. However the expression level of GAPDH was consistent in all treatment groups.
The study demonstrated that GAPDH is reliable as reference gene for quantitative gene expression analysis in HDFs. Therefore it can be used as housekeeping gene for quantitative real time RT-PCR technique in human diploid fibroblasts particularly in studying cellular senescence.
Cellular senescence has been most widely studied in fibroblast cells in vitro . In many human cells, cellular senescence is characterized by several molecular and cytological markers, such as a large flat morphology, expression of a senescence-associated β-galactosidase activity (SA β-gal) , and altered gene expression . Oxidative stress is a condition that arises when the production of reactive oxygen species (ROS) overwhelms the cellular antioxidant defences. Human diploid fibroblasts (HDFs) can develop a phenotype resembling senescence in response to oxidative stress. It has been reported that treating cells exogenously with certain hydrogen peroxide (H2O2) concentrations can trigger entry into a senescent-like state which is termed 'stress-induced premature senescence' (SIPS) .
The morphology of SIPS cells was found resembled of senescent cells; with gross enlargement and accumulation of granular cytoplasmic inclusions after 2 weeks exposure to low dose of H2O2 . If oxidative stress or free radicals are at least partly responsible for lifespan and aging, it follows that antioxidant should prolong life and retard aging. Tocotrienol is one subclass of vitamin E that can be found abundantly in palm oil, rice bran oil, barley, corn, oats and wheat . It acts effectively as an antioxidant because its hydrogen atom from the hydroxyl group on its chromanol ring can readily be donated to reduce free radicals, each has its own biological activity .
RNA differential displays  and the serial assessment of gene expression  have been applied to explore senescence-associated genes to gain an insight into the molecular mechanisms underlying senescence . The study of the molecular mechanisms underlying senescence has shed light on central aspects of tumor development and has contributed to the research on organismal aging . Quantification of transcription levels of genes plays a central role in the understanding of gene function . Therefore, quantitative real time RT-PCR has become a popular means to assess mRNA expression level; due to its sensitivity, accuracy and ability to amplify mRNA signal . The method allows detection of amplicon accumulation since it is performed using fluorogenic probes or intercalating dyes such as SYBR Green I, rather than by conventional end-point analysis . However, it is essential to control for error between samples when measuring RNA expression. Therefore, in order to control for experimental variations in the amount of RNA used in each quantitative RT-PCR and batch-to-batch variations in PCR reagents, coincident measurement of so-called 'housekeeping' genes has been used for the normalization of target gene expression data .
Housekeeping genes or reference genes are essential endogenous regulatory genes that are involved in various processes in the cell, such as metabolism, cell structure, gene transcription, and homeostasis and are therefore constitutively expressed . Choosing an appropriate reference gene is important for accurate quantitative RNA expression in real time RT-PCR technique. In using the relative quantitative RT-PCR method, the cycle thresholds of the genes of interest were compared to the housekeeping genes to determine relative changes in expression . The expression levels of reference genes should remain constant between the cells of different tissues and under different experimental conditions . If these requirements are not fulfilled then normalization to varying internal references can lead to increased 'noise' or erroneous results . If the chosen housekeeping gene fluctuates randomly between samples, then small differences between genes of interest will be missed. Appropriate validation of internal references is therefore crucial to avoid misinterpretations of study findings .
Several genes have been used as housekeeping genes, including β-actin, β2-microglobulin, cyclooxygenase 1, glyceraldehyde 3-phosphate dehydrogenase (GAPDH), hypoxanthine phosphoribosyl transferase, porphobilinogen deaminase, and the transferring receptor . The RNA encoding GAPDH is universally expressed. GAPDH catalyzes the oxidative phosphorylation of glyceraldehyde 3-phosphate to 1,3-bisphosphoglycerate during glycolysis as well as the reverse reaction in tissues involved in gluconeogenesis. GAPDH has also been implicated in other ubiquitous processes such as DNA replication and repair and apoptosis . This gene has been used for real time comparative gene expression studies. However, recent research has demonstrated that the expression level of housekeeping genes may be altered due to differences in experimental treatment . Therefore, it is important to validate the stability and elucidate the changes of the housekeeping genes between different experimental conditions. In order to determine the suitability of GAPDH as reference gene in senescent and antioxidant studies, total RNA from HDFs after different experimental treatments were used to rectify the consistency of GAPDH expression in different reaction conditions of human skin fibroblast cells senescent model.
Morphology of human diploid fibroblast cells
Young HDFs displayed spindle shape of adherent cells (Figure 1A). When the cells were at the end of replicative lifespan at passage 30, its morphology changed to enlarged cellular with flattening shape and all cells had higher ratio of cytoplasm:nucleus content (Figure 1B). The gradual loss of replicative potential results in reduced in cell yield. Similar senescence morphology was observed in HDFs induced with oxidative stress (Figure 1C). However, no morphological changes was observed in cells treated with γ-tocotrienol when compared to young HDFs (Figure 1D).
Senescence-associated (SA) β-galactosidase expression
Positive senescence-associated (SA) β-galactosidase staining was markedly increased in senescent cells at passage 30 (Figure 2A) and cells with H2O2-induced oxidative stress (Figure 2B). However, no positive staining was observed in young fibroblast cells (Figure 2C) and cells treated with γ-tocotrienol (Figure 2D).
GAPDH mRNA expression
Fluorescence was measured during each PCR cycle and the amount of fluorescence was proportional to the amount of the PCR product. The amplification graph showed that the Ct value for GAPDH of all treatment groups was within the range of 13-16 (Figure 3A) and was inversely proportional to the total RNA concentration isolated from HDFs with different experimental treatments (Table 1) indicating that GAPDH expression level was consistent with the concentration of total RNA (50-100 ng/μl). The melting curve showed the specificity of the GAPDH primer that has been designed, with each PCR product demonstrated one specific melting temperature (Figure 3B). Subsequently, the agarose gel electrophoresis which was performed to confirm the size of the PCR product in various treatment groups showed that the bands were specific, according to 217 bp (Figure 3C). The standard curve drawn with Ct value versus total RNA concentration showed the PCR reaction was at 121.4% of efficiency and correlation coefficient of 0.992 (Figure 3D).
Most type of primary normal cells did not proliferate indefinitely in culture. Instead, after a period of rapid proliferation, their division rate slows down and ultimately ceased altogether. Such cells became unresponsive to mitogenic stimuli yet can remain viable for extended periods of time. Upon entering the state of senescence, cells underwent a dramatic change in morphology; whereby their volume increases and they loss their original shape, acquiring a flattened cytoplasm . These morphological changes were shown in our senescent HDFs and cells with H2O2-induced oxidative stress. This shift was accompanied by changes in nuclear structure, gene expression, protein processing, and metabolism . Although there is limited study on the effects of tocotrienol and human life span, but it was found that tocotrienol administration reduces oxidative protein damage and extends the mean life span of Caenorhabditis elegans, one of the well established aging model . These results strongly suggest that oxidative stress plays a role in bringing the changes in cellular function that occurs during aging.
In the present study, the RNA encoding GAPDH was ubiquitously expressed in all treatment groups. GAPDH is frequently used as an endogenous control for quantitative RT-PCR analysis because its expression is consistent at different time points and various experimental manipulations . With the growing popularity of quantitative RT-PCR, one needs to justify whether the use of GAPDH as the specific internal standard RNA is appropriate . As pointed out by Bustin (2002) , the ease and rapidity of which data was acquired by quantitative RT-PCR can easily create a false sense of objectivity. This assay used small amounts of RNA and therefore, is more prone to errors due to variation in RNA input .
The finding of this study indicated that GAPDH expression level was stable with total RNA concentration in the range of 50-100 ng/μl. This data suggested that GAPDH showed the suitability criteria as a housekeeping gene for human diploid fibroblasts particularly for studying cellular senescence. Both the H2O2 stress-induction and tocotrienol treatment did not essentially change the expression level of GAPDH. Moreover the real-time profile designed in this study with the highly specific primers used caused the reaction efficiency to be optimum with high correlation of the reaction. This finding was similar to Touchberry et al. (2006)  which reported that the average expression level of GAPDH among other genes in human skeletal muscle had the smallest different between old and young groups. Their findings also demonstrated that GAPDH was a stable housekeeping gene and appeared to be an effective gene for group comparisons between young and old populations. However, they also highlighted the limitation of GAPDH as housekeeping gene for comparison within groups due to its large variation from the mean. The high standard deviation scored within groups could be due to mRNA isolated from different individual with high biological variation. In comparison, our study used the same fibroblasts donors for all experiment groups. Thus, the biological variation that may affect GAPDH expression has been minimised.
The suitability of GAPDH as housekeeping gene has also been demonstrated in articular chondrocytes gene expression assessment in hypoxic condition . GAPDH expression was not modulated in the chondrocytes by changes in oxygen tension. Studies on human glioblastoma and cancer cells of different origins for its molecular regulation under hypoxic condition also revealed GAPDH expression was not involved in the gene regulation [29, 30]. In contrast, Zhong & Simons (1999)  showed GAPDH was increased in cancer cells under hypoxic conditions. GAPDH expression was correlated with the upregulation of hypoxia inducible factor-1 alpha. In another study, GAPDH was suggested as the target gene in the evaluation of amino-bisphosphophates effects on prostate and breast cancer cell lines . The expression of GAPDH was significantly decreased in a dose-dependent manner following amino-bisphosphophates treatment to the cancer cells. GAPDH was also not suitable as housekeeping gene for cell lines under mitogens stimulation . Therefore, it was highly recommended that the suitability of GAPDH as housekeeping gene be properly validated for each experiment to confirm that its expression was unaffected by different experimental conditions .
This finding recommended that GAPDH should be the reference gene for future quantitative mRNA analysis involving human diploid fibroblasts particularly in studying cellular senescence.
This research has been approved by the National University of Malaysia Ethical Committee (Approval Project Code: FF-104-2007). Written consents were obtained from all subjects before tissue collection. Primary HDFs were derived from circumcision foreskins of 9-12 year-old boys.
Skin samples were aseptically collected and rinsed several times with 75% alcohol and phosphate buffered saline containing 1% antibiotic-antimycotic (PAA, Austria) solution. After removing the epidermis, the dermis was cut into small pieces and transferred to 50 ml falcon tube (BD Bioscience, USA) containing 0.03% collagenase type I digestive buffer (Worthington Biochemical Corporation, USA). Dermis was digested in an incubator shaker at 37°C for 6 hours. The isolated HDFs were then centrifuged and rinsed with PBS twice. HDFs were cultured in Dulbecco's Modified Eagle Medium (DMEM, Flowlab, Australia) in T25 culture flasks containing 10% fetal bovine serum (FBS, PAA, Austria) and 1% antibiotic-antimycotic at 37°C in 5% CO2 humidified incubator. After 5 to 6 days, the cultured fibroblasts were harvested by trypsinization using 2.5% trypsin containing 0.03% of EDTA (Flowlab, Australia) and culture-expanded in T25 culture flasks. When the subcultures reached 80-90% confluency, serial passaging was done by trypsinization. The HDFs senescent model was consist of 4 different treatment groups; i) young (passage 4), ii) senescent (passage 30), iii) H2O2-induced oxidative stress (10 μM of H2O2 exposure for two weeks) and iv) γ-tocotrienol-treated group (20 μM treatment of γ-tocotrienol prior to H2O2 exposure). The morphology of the cells was assessed daily using inverted microscope.
Morphological analysis and senescence-associated (SA) β-galactosidase staining
The molecular marker of cell-aging in vitro for HDFs (SA β-gal activity) was determined by senescent cells staining kit (Sigma, USA) according to the manufacturer's instruction. Senescent cells were stained blue after 4 h of incubation with β-galactosidase staining solution containing 5-bromo-4-chloro-3-indolyl-β-D-galactosidase (X-gal) at 37°C.
Primers for human GAPDH were designed from listed NIH GenBank database using Primer 3 software and blasted against GenBank database sequences for specificity confirmation. The sequence of GAPDH primers is shown in Table 2.
Total RNA from cultured HDFs in different treatment groups was extracted using TRI Reagent (Molecular Research Center, USA) according to the manufacturer's instruction. Polyacryl Carrier (Molecular Research Center, USA) was added in each extraction to precipitate the total RNA. Extracted RNA pellet was then washed with 75% ethanol and dried before dissolved in RNase and DNase free distilled water . The yield and purity of the extracted total RNA were determined by Nanodrop (Thermo Scientific, USA) in triplicates. Total RNA concentration was then adjusted to 50-100 ng/μl before stored at -80°C immediately after extraction.
Real time RT-PCR
Quantitative real-time RT-PCR reaction was carried out using 1 μl (50-100 ng/μl) total RNA as template; 1 ul of forward and reverse primers of GAPDH and iScript One-Step RT-PCR reagent with SYBR Green (Bio-Rad, USA). Reactions were conducted using Bio-Rad iCycler with reaction profile as follows; cDNA synthesis for 30 min at 50°C; pre-denaturation for 2 min at 94°C; PCR amplification for 38 cycles with 10 sec at 94°C and 30 sec at 61°C. This was followed by a melt curve analysis to determine the reaction specificity. Agarose gel electrophoresis was performed for confirmation of the size of PCR product.
Total RNA extracted from various samples was serially diluted in nuclease-free water by 2 fold. The diluted total RNA was used as template for real-time RT-PCR and the experiment was conducted in duplicate for each sample. Master mix without total RNA was prepared for all reactions and 24 μl was aliquoted into each reaction tube and 1 μl of the diluted total RNA was then added in individually. Mean of Ct value was determined after the reaction in order to determine the linearity of the GAPDH expression level.
Wright WE, Shay JW: Historical claims and current interpretations of replicative aging. Nature Biotechnology. 2002, 20: 682-688. 10.1038/nbt0702-682.
Peter DA: Remodeling of chromatin structure in senescent cells and its potential impact on tumor suppression and aging. Gene. 2007, 397 (1-2): 84-93. 10.1016/j.gene.2007.04.020.
Smith JR, Pereira-Smith OM: Replicative senescence: implications for in vivo aging and tumor suppression. Science. 1996, 273: 63-67. 10.1126/science.273.5271.63.
Lu T, Finkel T: Free radicals and senescence. Experiment Cell Research. 2008, 314: 1918-1922. 10.1016/j.yexcr.2008.01.011.
Duan J, Duan J, Zhang Z, Tong T: Irreversible cellular senescence induced by prolonged exposure to H2O2 involves DNA-damage-and-repair genes and telomere shortening. International Journal of Biochemistry and Cell Biology. 2005, 37: 1407-1420. 10.1016/j.biocel.2005.01.010.
Packer L, Weber SU, Rimbach G: Molecular aspects of alpha tocotrienol antioxidant action and cell signalling. The Journal of Nutrition. 2001, 131: 369S-73.
Xu WL, Liu JR, Liu HK, Qi GY, Sun WG, Chen BQ: Inhibition of proliferation and induction of apoptosis by γ-tocotrienol in human colon carcinoma HT-29 cells. The Journal of Nutrition. 2009, 25: 555-566. 10.1016/j.nut.2008.10.019.
Feng J, Andrews WH, Enlow BE, Saati SM, Tonkin LA, Villeponteau B: Cataloging altered gene expression in young and senescent cells using enhanced differential display. Nucleic Acids Research. 1995, 23: 3244-3251. 10.1093/nar/23.16.3244.
Welle S, Brooks A, Thornton CA: Senescence-related changes in gene expression in muscle: similarities and differences between mice and men. Physiological Genomics. 2001, 5: 67-73.
Yoon IK, Kim HK, Kim YK, Song IH, Kim W, Kim S, Baeka SH, Kim JH, Kim JR: Exploration of replicative senescence-associated genes in human dermal fibroblasts by cDNA microarray technology. Experimental Gerontology. 2004, 39: 1369-1378. 10.1016/j.exger.2004.07.002.
Ben-Porath I, Weinberg RA: When cells get stressed: an integrative view of cellular senescence. The Journal of Clinical Investigation. 2004, 113: 8-13.
Overbergh L, Giulietti A, Valckx D, Decallonne B, Bouillon R, Mathieu C: The use of Real-Time Reverse Transcriptase PCR for the Quantification of Cytokine Gene Expression. Journal of Biomolecular Techniques. 2003, 14: 33-43.
Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. Journal of Biomolecular Techniques. 2004, 15: 155-166.
Peter IR, Helps CR, Hall EL, Day MJ: Real time RT PCR: considerations for efficient and sensitive assay design. Journal of Immunological Methods. 2004, 286: 203-217. 10.1016/j.jim.2004.01.003.
Robert DB, Dan WH, Robert AC, Brian JC: GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiological Genomics. 2005, 21: 389-395. 10.1152/physiolgenomics.00025.2005.
Touchberry CD, Wacker MJ, Richmond SR, Whitman SA, Godard MP: Age-Related Changes in Relative Expression of Real-Time PCR Housekeeping Genes in Human Skeletal Muscle. Journal of Biomolecular Techniques. 2006, 17: 157-162.
Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Heinen E: Housekeeping genes as internal standards: use and limits. Journal of Biotechnology. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7.
Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of Molecular Endocrinology. 2000, 2: 169-193. 10.1677/jme.0.0250169.
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.
Lee PD, Sladek R, Greenwood CMT, Hudson TJ: Control genes and variability: absence of ubiquitous reference transcripts in diverse mammalian expression studies. Genome Research. 2001, 12: 292-297. 10.1101/gr.217802.
Corbin IR, Gong Y, Zhang M, Minuk GY: Proliferative and nutritional dependent regulation of glyceraldehyde-3-phosphate dehydrogenase expression in the rat liver. Cell Proliferation. 2002, 35: 173-182. 10.1046/j.1365-2184.2002.00236.x.
Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: Validation by real time, quantitative RT-PCR. Journal of Biochemical and Biophysical Methods. 2000, 46: 69-81. 10.1016/S0165-022X(00)00129-9.
Campisi J: Cancer, aging and cellular senescence. In Vivo. 2000, 14: 183-188.
Adachi H, Ishii N: Effects of tocotrienols on life span and protein carbonylation in Caenorhabditis elegans. The Journals of Gerontology Series A: Biological Sciences and Medical Sciences. 2000, 55: B280-B285.
Winer J, Jung CK, Shackel I, Williams PM: Development and validation of real-time quantitative reverse transcriptase-polymerase chain reaction for monitoring gene expression in cardiac myocytes in vitro. Analytical Biochemistry. 1999, 270: 41-49. 10.1006/abio.1999.4085.
Butler ST, Marr AL, Pelton SH, Radcliff RP, Lucy MC, Butler WR: Insulin restores GH responsiveness during lactation-induced negative energy balance in dairy cattle: effects on expression of IGF-I and GH receptor 1A. Journal of Molecular Endocrinology. 2003, 176: 205-217.
Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. Journal of Molecular Endocrinology. 2002, 29: 23-39. 10.1677/jme.0.0290023.
Grimshaw MJ, Mason RM: Modulation of bovine articular chondrocyte gene expression in vitro by oxygen tension. Osteoarthritis Cartilage. 2001, 9: 357-64. 10.1053/joca.2000.0396.
Said HM, Hagemann C, Stojic J, Schoemig B, Vince GH, Flentje M, Roosen K, Vordermark D: GAPDH is not regulated in human glioblastoma in vivo or in vitro under hypoxic conditions. BMC Molecular Biology. 2007, 8: 55-10.1186/1471-2199-8-55.
Said HM, Polat B, Hagemann C, Anacker J, Flentje M, Vordermark D: Absence of GAPDH regulation in tumor-cells of different origin under hypoxic conditions in - vitro. BMC Research Notes. 2009, 2: 8-10.1186/1756-0500-2-8.
Zhong H, Simons JW: Direct comparison of GAPDH, beta-actin, cyclophilin, and 28 S rRNA as internal standards for quantifying RNA levels under hypoxia. Biochemistry Biophysical Research Communication. 1999, 259 (3): 523-526. 10.1006/bbrc.1999.0815.
Valenti MT, Bertoldo F, Dalle Carbonare L, Azzarello G, Zenari S, Zanatta M, Balducci E, Vinante O, Lo Cascio V: The effect of bisphosphonates on gene expression: GAPDH as a housekeeping or a new target gene?. BMC Cancer. 2006, 6: 49-10.1186/1471-2407-6-49.
Radonić A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochemistry Biophysical Research Communication. 2004, 313 (4): 856-862. 10.1016/j.bbrc.2003.11.177.
Kenneth JL, Thomas DS: Analysis of relative gene expression data using real time quantitative PCR and the 2-ΔΔCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
Chua KH, Aminuddin BS, Fuzina NH, Ruszymah BHI: Insulin-transferrin-selenium prevent human chondrocyte dedifferentiation and promote the formation of high quality tissue engineered human hyaline cartilage. European Cells and Materials. 2005, 9: 58-67.
This study was funded by the Ministry of Science, Technology and Innovation under the E-Science Fund 02.01.02. SF0027 and National University of Malaysia.
AZ carried out the lab work and drafted the manuscript. KHC and NAR were involved in optimising the HDFs primary culture and revising the manuscript. SM was the Principal Investigator who designed the study and revised the manuscript. All authors read and approved the final manuscript.
Azalina Zainuddin, Kien Hui Chua, Norhazira Abdul Rahim and Suzana Makpol contributed equally to this work.
About this article
Cite this article
Zainuddin, A., Chua, K.H., Rahim, N.A. et al. Effect of experimental treatment on GAPDH mRNA expression as a housekeeping gene in human diploid fibroblasts. BMC Molecular Biol 11, 59 (2010). https://doi.org/10.1186/1471-2199-11-59
- Cellular Senescence
- Senescent Cell
- GAPDH Expression
- Human Diploid Fibroblast
- Hypoxanthine Phosphoribosyl Transferase