Normalization with genes encoding ribosomal proteins but not GAPDH provides an accurate quantification of gene expressions in neuronal differentiation of PC12 cells
© Zhou et al; licensee BioMed Central Ltd. 2010
Received: 24 August 2009
Accepted: 29 January 2010
Published: 29 January 2010
Gene regulation at transcript level can provide a good indication of the complex signaling mechanisms underlying physiological and pathological processes. Transcriptomic methods such as microarray and quantitative real-time PCR require stable reference genes for accurate normalization of gene expression. Some but not all studies have shown that housekeeping genes (HGKs), β-actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which are routinely used for normalization, may vary significantly depending on the cell/tissue type and experimental conditions. It is currently unclear if these genes are stably expressed in cells undergoing drastic morphological changes during neuronal differentiation. Recent meta-analysis of microarray datasets showed that some but not all of the ribosomal protein genes are stably expressed. To test the hypothesis that some ribosomal protein genes can serve as reference genes for neuronal differentiation, a genome-wide analysis was performed and putative reference genes were identified based on stability of expressions. The stabilities of these potential reference genes were then analyzed by reverse transcription quantitative real-time PCR in six differentiation conditions.
Twenty stably expressed genes, including thirteen ribosomal protein genes, were selected from microarray analysis of the gene expression profiles of GDNF and NGF induced differentiation of PC12 cells. The expression levels of these candidate genes as well as ACTB and GAPDH were further analyzed by reverse transcription quantitative real-time PCR in PC12 cells differentiated with a variety of stimuli including NGF, GDNF, Forskolin, KCl and ROCK inhibitor, Y27632. The performances of these candidate genes as stable reference genes were evaluated with two independent statistical approaches, geNorm and NormFinder.
The ribosomal protein genes, RPL19 and RPL29, were identified as suitable reference genes during neuronal differentiation of PC12 cells, regardless of the type of differentiation conditions. The combination of these two novel reference genes, but not the commonly used HKG, GAPDH, allows robust and accurate normalization of differentially expressed genes during PC12 differentiation.
During development, neurons make networks of connections with other neurons by growing axons and dendrites. These neuronal out-growths are regulated by extracellular cues that signal to cells resulting in phenotypic changes. A major challenge is the identification of molecular mechanisms underlying this highly complex and interactive network in terms of the functions of genes and proteins.
Currently, transcriptomic methods are widely used as an initial step in unraveling the complex signaling mechanisms underlying physiological and pathological processes and in neuronal differentiation [2–5]. Gene microarray offers a high throughput platform for the analysis of the entire transcriptome to identify differentially expressed genes. Reverse transcription quantitative real-time PCR (RT-qPCR), with a wider dynamic range of quantification and higher assay sensitivity and precision, is often used to corroborate microarray findings [6, 7]. Regardless of the method used, normalization, a critical process of adjusting the expression measurements between samples to compensate for various sources of variability in the assay, is essential to allow accurate comparisons of the results between different samples and conditions [8, 9]. Normalization with internal reference gene is used to control for technical and biological variations introduced during both sample preparation and detection by RT-qPCR . It has also been shown to be suitable for the normalization of partially degraded RNA samples [11–13].
With nearly all normalization methods, the assumption that one or more reference genes are constitutively expressed at near-constant levels under all experimental conditions is implicit and the expression levels of all other genes in the sample are then scaled to these reference genes accordingly. It is common to use reference genes selected from an assumed list of "housekeeping" genes (HKGs) which typically include transcripts such as GAPDH and ACTB [9, 10, 14]. A number of studies have now shown that the expressions of these genes, in some but not all experimental conditions, are altered significantly[15–18], thus, making the choice of using these HKGs for normalization uncertain without a priori knowledge.
A variety of approaches have been employed to enable better selection of reference genes. One approach is the use of statistical algorithms, for example, geNorm , Best keeper , NormFinder , Global Pattern Recognition , and Equivalence tests , to evaluate the relative expression stabilities of genes from a pool of predefined lists of candidates. While this approach is certainly more robust than using the single gene methods, it too is based on potentially unfounded assumptions about which genes may be stably expressed in the conditions studied. These genes are still required to be pre-selected and incorporated into the experimental designs without any a priori evidence to support their use. An alternative and less biased approach is the meta-analysis of large scale gene expression profiles to identify stably expressed genes [23–26]. A selected number of potential references genes can then be validated experimentally and the stability of expressions analyzed by the above mentioned statistical algorithms in defined experimental settings.
To date, reference genes validated for neuronal differentiation studies have not been reported yet. The present study aims to identify suitable reference genes during chemically induced neuronal differentiation of PC12, a cell-line derived from a pheochromocytoma of the rat adrenal medulla. Because of its unique cellular properties, suitability for genetic and biochemical manipulations, the PC12 cell-line is widely regarded as a convenient alternative to endogenous neuronal cells, and serves as a commonly used model system for studies on neuronal differentiation [27–29]. For example, in response to NGF, PC12 cells stop dividing, elaborate neuronal processes, are electrically excitable, and have the potential to form synapses when co-cultured with muscle cells . Here, we measured the temporal expression of twenty novel candidate reference genes identified from microarray studies and the commonly used HKGs, ACTB and GAPDH, at various stages of PC12 differentiation. Based on two independent statistical approaches ["pairwise comparison"  and "model based variation" , the expressions of ribosomal protein genes RPL19 and RPL29 were found to be highly stable regardless of pharmacological treatments and stages of differentiation. The combination of the two reference genes was sufficient to allow robust and accurate normalization of differentiation related genes.
Selection of candidate reference genes from microarray data
Selection of candidate reference genes from microarray data
Ribosomal protein L29
Ribosomal protein L10A
Vps20-associated 1 homolog
Similar to ribosomal protein L15
Similar to ribosomal protein L7-like 1
Ribosomal protein S15
Acidic ribosomal phosphoprotein P0
Ribosomal protein L14
Eukaryotic translation elongation factor 1 alpha 1
Ribosomal protein S15a
Ribosomal protein L18
RalBP1 associated Eps domain containing protein (predicted)
chromatin modifying protein 2B
CCR4-NOT transcription complex, subunit 8
RNA terminal phosphate cyclase domain 1
Ribosomal protein L19
NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 6,
Ribosomal protein L9
Similar to 40S ribosomal protein S3
Ribosomal protein L3
Real-time PCR validation of novel candidate reference genes
Stabilities of candidate reference genes and common housekeeping genes
Further analysis of candidate gene stabilities in each treatment group (Additional File 2A) or at specific time point (Additional File 2B) revealed that the stability rankings of candidate genes do vary among different subgroups. However, with the exception of RPL29 in KCl treated samples, the two genes RPL19 and RPL29 were consistently ranked among top 5 in all subgroups. In contrast, the stability rankings of GAPDH and ACTB varied considerably among different subgroups and they were ranked among the least stable ones within the group of 22 genes in several subgroups. The data indicated that the two novel candidate genes RPL19 and RPL29 have higher expression stabilities than both GAPDH and ACTB, and may serve as better normalizers for gene expression in neuronal differentiation of PC12 cells.
Comparison of the normalization factors generated by different reference gene(s)
Effect of different reference genes on the interpretation of target gene regulation
Similarly, we investigated the effect of different reference gene(s) on normalized target gene expressions in GDNF, Forskolin, KCl and ROCK inhibitor Y27632 treated samples. Similar to the case of NGF treatment, with GDNF stimulated PC12-GFRa1a/RET9 and PC12-GFRa1a/RET51 cells, normalization by GAPDH resulted in the underestimations of GDNF induced upregulation of EGR1 and ITGA1 expressions; and over-estimated CRYAB down-regulation (Additional File 5A, B). Interestingly, normalization by ACTB was found to overestimate the expression of EGR1 and ITGA1 expressions in PC12-GFRa1a/RET51 but not RET9 cells (Additional File 5B, 72 h), highlighting the subtle differences between GFRa1a/RET9 and GFRa1a/RET51 systems.
For Forskolin and ROCK inhibitor Y23672 differentiated samples, normalization by ACTB consistently led to the over-estimations of target gene expressions (Additional File 5C, D). Depending on the time point analyzed, normalization by GAPDH was shown to result in either underestimations or overestimations of target gene expressions (Additional File 5C, D). In KCl treated samples, no statistical significant difference was observed among NFtop2 (RPL19/REPS1), NFRPL19/RPL29, NFACTB and NFGAPDH normalized target gene expression, which suggested that all four were acceptable reference gene(s) for this particular experimental condition (data not shown).
It is thus evident that the stabilities of the two most commonly used HKGs, GAPDH and ACTB vary across different experimental conditions during neuronal differentiation of PC12 cells. They were acceptable reference genes under some conditions but may significantly under- or over-estimate target gene expression under others. On the contrary, the two novel candidate reference genes RPL19 and RPL29 were stably expressed among all conditions analyzed and allowed accurate normalization of differentially regulated genes during PC12 differentiation. It is worthy to note that at early time points (0.5 h and 6 h), the expressions of EGR1, ITGA1 and CRYAB did not show any significant differences when scaled with either NFtop2, NFRPL19/RPL29, NFACTB or NFGAPDH (data not shown). This is consistent with the observation that the expression of GAPDH did not change significantly at the early time points (Figure 6. and Additional File 4).
Twenty candidate reference genes that showed little variation but high expression in PC12 cells differentiated with NGF and GDNF were first selected from microarray datasets using two independent statistical algorithms. Together with two well studied HKGs, the expression stabilities of these candidate reference genes were further analyzed using RT-qPCR in cells differentiated with other stimuli. From these studies, unexpectedly, RPL19 and RPL29 but not the HKGs, were identified as suitable reference genes that can be used for normalization of gene expression in neuronal differentiation of PC12 induced by a variety of chemical stimuli.
Neuronal differentiation is a process where cells undergo enormous morphological changes, over a period of several days. It is accompanied by substantial biochemical changes including cell cycle exit , changes in metabolism [37, 38] and alteration in structural proteins [39, 40]. Since the commonly used reference genes are mostly structural proteins or enzymes involved in metabolism, it is especially important to validate the stabilities of these genes during the process of differentiation. Many of these studies investigated gene expression changes in PC12 but few have evaluated the suitability of HKGs as normalizers in this model. Our microarray analysis revealed that a group of novel candidate genes was more stably expressed than commonly used HKGs ACTB and GAPDH, suggesting that ACTB and GAPDH may not be ideal reference genes in neuronal differentiation of PC12 cells.
In an effort to gain an insight into the temporal regulation of genes during neuronal differentiation, it is necessary that the reference genes used are stably expressed over a period of days. GAPDH and ACTB have been used for normalization in more than 90% of previous reports , often without proper validation of their stabilities. Numerous publications have reported that such HKGs can be differentially expressed under various experimental paradigms and are therefore inappropriate for normalization [9, 15, 16]. However, there are also recent reports that these HKGs are stably expressed and can serve as reference genes [17, 18, 42]. Most genes, including GAPDH and ACTB, examined in this study were stable at early stages of differentiation (0.5 h or 6 h). However, as differentiation proceeded with dramatic morphological changes and concomitant biochemical changes, the instability of expressions of GAPDH and many of the genes examined was obvious. In the case of GAPDH, this instability issue correlated well with the temporal increase of expression level, which peaked at 28 h and was sustained over a period of 72 h. While GAPDH may still serve as a reference gene for PC12 cells under specific conditions, the validity of using this gene and other less stable ones should be experimentally verified. However, the two RP genes (RPL19 and RPL29) that showed good stability in expression over the period of differentiation provided an optimal pair of reference genes for the entire period of and various experimental conditions for neuronal differentiation.
Among the twenty candidate genes selected, thirteen were ribosomal protein genes, suggesting that the family of ribosomal protein genes may become yet another source of reference genes. Several recent publications have validated and recommended the use of ribosomal protein genes as reference genes [23, 25], while others have reported their tissue-dependent variations . A plausible explanation for such disparity is the large number of ribosomal protein genes present in mammalian systems (80 genes in human, mouse and rat genome), which may be stably or differentially expressed depending on the tissue type or experimental conditions. At present, relatively little is known about these mammalian ribosomal proteins, as compared to their bacterial and archael counterparts . While bacterial ribosomal protein genes exist largely in clusters, the mammalian RP genes are dispersed throughout the genome . Some have suggested that all of these proteins are intimately involved in ribosome production and could be co-regulated. Depletion of a particular ribosomal protein would generally cause a reduction of all other ribosomal proteins in the same ribosome sub-unit . Other reports have shown that some ribosomal protein genes could be regulated independent of others . Recently, extra-ribosomal functions of some of these proteins have been reported [48–51], suggesting that they may be individually regulated. A previous study comparing random ESTs from naïve and NGF-treated PC12 cells, reported an NGF-promoted decrease in the expressions of RPL19 . However, this decrease in RPL19 was not observed in other studies using SAGE  or microarray . Similar to the latter studies, we too did not observe changes in RPL19 transcripts with NGF-treated PC12. Moreover, the SAGE study but not the microarray analysis reported a significant decrease in RPL29 expression. Using both microarray and RT-qPCR, we have also shown that RPL29 was unchanged when the cells were differentiated. The reasons for these discrepancies are unclear and may be due to the differences in methods used. We have shown here by quantitative real-time PCR that some ribosomal protein genes, RPL19 and RPL29, are highly stably expressed and are thus suitable reference genes, whereas others like RPL9 and RPL18 can vary significantly during differentiation.
Unlike some studies that attempted to identify ideal reference genes through meta-analysis of many publically available microarray data which includes a diverse range of tissue types and experimental conditions, this study was designed to specifically identify a set of suitable reference genes for PC12 cells undergoing neuronal differentiation. We have performed both the microarray analysis and RT-qPCR validation on biological samples prepared with the same techniques and reagents, thus minimizing variations introduced by differences in sample preparation methods and assay platforms. We have also systematically evaluated the effect of the use of NFs of inappropriate reference gene(s) on the expression changes of the target genes and the erroneous results they resulted in. Thus, with neuronal differentiation of PC12 cells, scaling with the geometric means of the expressions of RPL19 and RPL29 is recommended for the accurate normalization of gene expression. Whether these two genes are suitable for normalization of neuronal differentiation in other systems remains to be evaluated.
Twenty novel candidate reference genes were identified and their expression stabilities were analyzed and compared to that of commonly used HKGs ACTB and GAPDH. Through this systematic study that included both microarray analysis and RT-qPCR, we have found two ribosomal protein genes RPL19, and RPL29 to be stably expressed during neuronal differentiation of PC12 cells, induced by five different chemical stimuli, over 72 h. The combination of these two novel reference genes allowed robust and accurate normalization of differentially expressed genes, regardless of stimuli and stages of differentiation. In contrast, the use of an inappropriate reference gene like GAPDH led to significant erroneous estimation of differentially expressed genes.
The rat pheochromocytoma cell line PC12 (catalog # CRL-1721; American Type Culture Collection) cells were grown in DMEM supplemented with 10% heat-inactivated fetal bovine serum (FBS; Hyclone, Logan, UT) and 5% Horse Serum (HS), in a humidified atmosphere with 5% CO2 at 37°C. Wild type PC12 cells, that do not endogenously express GFRa or RET (data not shown), were co-infected with murine GFRa1a (NM_010279) and RET9 (NM_001080780) or RET51 (NM_009050) in pQCXIN or pQCXIH vector by retro-viral infection (Clontech, Mountain View, CA) and selected with 0.4 mg/ml G418 and 0.1 mg/ml Hygromycin, over a period of 2 months.
Differentiation, sample collection and assessment of neurite outgrowth
Two million wild type or infected PC12 cells were seeded in 75 cm2 flask (NUNC, Finland) overnight, in DMEM supplemented with 10% FBS and 5% HS, followed by serum depletion for 12 h. PC12 cells were then treated with 50 ng/ml recombinant human GDNF (Peprotech, NJ), 50 ng NGF (Peprotech, NJ), 10 μM Forskolin (Sigma, St. Louis, MO), 5 mM KCl (Sigma, St. Louis, MO) or 25 μM ROCK inhibitor Y27632 (Calbiochem, USA) in DMEM to induce neurite outgrowth. Total RNA was isolated from control and treated cells at 0.5, 6, 24 or 72 h. For neurite outgrowth assessment, cells bearing at least one neurite with the length equivalent to the cell bodies were scored at 72 h by independent observers. More than 400 cells from three different fields were counted per flask.
Control and treated PC12 Cells were fixed with 4% paraformaldehyde in 1xPBS for 15 min at 37°C, and subsequently in methanol at -20°C for an additional 15 min. After three washes with 1xPBS, cells were permeabilized and blocked with normal goat serum (1:10; Dako, Glostrup, Denmark) in 0.5% Triton X-100/PBS for 30 min at 37°C. The cells were then incubated with high-molecular-weight neurofilament protein (NF-200) antibody (Sigma, St. Louis, MO) at 1:80 dilution in 0.1% Triton X-100/0.1% BSA/1xPBS for 1 h at 37°C and washed three times in 1xPBS. Subsequently, the cells were incubated with goat anti-rabbit fluorescent secondary antibody (Alexa Fluor 488; Invitrogen, CA) diluted 1:1000 in 0.1% Triton X-100/0.1% BSA/1xPBS for 1 h. The cells were then washed three times in 1xPBS. Image acquisition was performed using the Zeiss inverted Axovert 25 microscope equipped with fluorescence detection (Oberkochen, Germany).
RNA Purification and cDNA Preparation
Total RNA from PC12 cells was prepared using TRIzol reagent (Invitrogen, CA) according to manufacturer's instruction. Total RNA was collected from samples in quadruplicate at each treatment time point and the integrity of the RNA validated by denaturing agarose gel electrophoresis and using the StdSens analysis chip on the Experion Automated Electrophoresis System (BioRad, CA) according to manufacturer's instructions. The Experion Automated Electrophoresis System assigns a RQI to each RNA electropherogram which ranges from 10 (intact RNA) to 1 (completely degraded RNA). RNA concentration was quantified using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE), and the 260/280 and 260/230 ratios were examined for protein and solvent contamination. Five micrograms of total RNA were reverse transcribed in a total volume of 20 μl containing 400 U of ImpromII and 0.5 μg of random hexamer (Promega, Madison, WI) for 60 min at 42°C according to the manufacturer's instructions. The reaction was terminated by heating at 70°C for 5 min, and the cDNA was diluted 1:20 for quantitative real-time PCR.
PC12 cells were seeded on 25 cm2 flask in complete medium and subsequently incubated for 12 h in serum free DMEM. The cells were then treated with GDNF (50 ng/ml) or NGF (50 ng/ml) for 0.5 h or 72 h in duplicates. Total RNA was isolated, quantified and integrity verified before it was amplified using Ambion Illumina RNA Amplification kit (Ambion, TX, USA). Briefly, total RNA (500 ng) was reverse transcribed by ArrayScript in the presence of T7 Oligo(dT) primer. Second strand of the cDNA was synthesized by DNA polymerase at 16°C for 2 h. The cDNA was purified and in vitro transcribed with T7 RNA polymerase and biotin-NTPs. Biotin-labeled cRNA samples were purified and quantified by ND-1000 spectrophotometer (NanoDrop, Fisher Thermo, DE, USA). Each cRNA (750 ng) was hybridized to RatRef-12 Expression BeadChip (Illumina, San Diego, CA, USA) containing 22,523 probes for a total of 21,910 rat genes selected primarily from the NCBI RefSeq database (Release 16) according to instruction provided by Experienced User Card (11286340 Rev A, Illumina). After hybridization, washing and blocking, the BeadChip was incubated with Streptavidin-Cy3 solution (Amersham Biosciences, Piscataway, NJ, USA). Fluorescent signals were obtained from scans on the high resolution Illumina BeadArray reader, using a two-channel, 0.8 μm resolution confocal laser scanner. The Illumina BeadStudio software (Version 2.0) was used to extract fluorescence intensities and the raw fluorescent data was background subtracted and used for analysis. Background is defined as the average signal intensity estimated from the negative control bead types. Outliers are removed using the median absolute deviation method. Detection p-values produced by the BeadStudio software were corrected for multiple hypothesis testing.
Primer Design and Plasmid Standards
The Genbank accession for each candidate reference gene was retrieved from the Illumina microarray probe set and compared to the NCBI RefSeq database (Release 16; http://www.ncbi.nlm.nih.gov. Transcript splicing sites were retrieved from Ensembl http://www.ensembl.org. Where more than one transcript matched the probe, the sequences were aligned and the primers were designed to amplify the consensus region. Vector NTI Advance 10 was used to design two sets of primers for each target gene. The first set of primers generates an amplicon of ~300 bp and is used as a template for RT-qPCR of the targeted gene. The template was subcloned into pGEMT-easy (Promega) vector as previously described . The second set of primers was used for RT-qPCR and was designed to amplify a ~100 bp region within the ~300 bp template. Both primer sets were exon spanning to avoid amplification from genomic sequences. Where possible, primers for RT-qPCR were designed to target the same exons used in the Illumina Expression BeadChip. All primer sequences were evaluated for possible false priming to known rat sequences using the NCBI BLAST tool http://blast.ncbi.nlm.nih.gov/. All products generated after amplifications were verified by gel-electrophoresis and DNA sequencing.
Quantitative Real-Time PCR
Real-time PCR was performed on Biorad iCycler 4 Real-Time PCR Detection System (Bio-Rad, Hercules, CA) using SYBR Green I. The threshold cycles (Ct) were calculated using the iQ5 Optical system software version 2.0. Real-time PCR was performed after an initial denaturation for 10 min at 95°C, followed by 40 cycles of 30 s denaturation at 95°C, 30 s annealing at 60°C, and 30 s extension at 72°C. Fluorescent detection was performed at the annealing phase. Melt curve analysis was carried out at the end of the cycling to confirm that a single product had been amplified. Primer dimer formation in all the assays showed distinct melt characteristics from the correct amplicons. The reaction was performed in a total volume of 40 μl in 1× XtensaMix-SG (BioWORKS, Singapore), containing 2.5 mM MgCl2, 10 pmol of each primer, and 0.5 U of KlearTaq Hotstart DNA polymerase (KBioscience, UK). All real-time PCR quantification was performed simultaneously with linearized plasmid standards and a non-template control . As PCR is an exponential process, it can be described by the equation, Nn = N0(1 + ε)n, where Nn is the number of target molecules at cycle n, N0 is the initial number of target molecules, ε is the efficiency of amplification and n is the number of cycles. The efficiency of target amplification of an assay was determined from the slope of a plot of Ct (Threshold cycle) versus -log10 concentration of the initial number of target molecules. High efficiency of amplification has a slope approaching the value of 3.32 cycles for every 10-fold dilution of the target. The gene expression levels were interpolated from standard curves and expressed as absolute copy numbers. All Real-time PCR experiments were compliant with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines .
Statistical data analysis
Gene expression stability analysis using two publicly available software tools, geNorm http://medgen.ugent.be/genorm/ and NormFinder http://www.mdl.dk/, were carried out according to authors' instruction.
The authors would like to thank NUS Department of Biochemistry for providing the necessary laboratory equipments and administrative support. The authors would like to acknowledge Ms Chen Peiyi from NUS Department of Statistics and Applied Probability for her help with the statistical analysis.
- Qiu Z, Ghosh A: A brief history of neuronal gene expression: regulatory mechanisms and cellular consequences. Neuron. 2008, 60 (3): 449-455. 10.1016/j.neuron.2008.10.039.PubMedView ArticleGoogle Scholar
- Dijkmans TF, van Hooijdonk LW, Schouten TG, Kamphorst JT, Vellinga AC, Meerman JH, Fitzsimons CP, De Kloet ER, Vreugdenhil E: Temporal and functional dynamics of the transcriptome during nerve growth factor-induced differentiation. Journal of neurochemistry. 2008, 105 (6): 2388-2403. 10.1111/j.1471-4159.2008.05338.x.PubMedView ArticleGoogle Scholar
- Perez-Iratxeta C, Palidwor G, Porter CJ, Sanche NA, Huska MR, Suomela BP, Muro EM, Krzyzanowski PM, Hughes E, Campbell PA: Study of stem cell function using microarray experiments. FEBS letters. 2005, 579 (8): 1795-1801. 10.1016/j.febslet.2005.02.020.PubMedView ArticleGoogle Scholar
- Rhodes DR, Chinnaiyan AM: Integrative analysis of the cancer transcriptome. Nature genetics. 2005, 37 (Suppl): S31-37. 10.1038/ng1570.PubMedView ArticleGoogle Scholar
- Preuss TM, Caceres M, Oldham MC, Geschwind DH: Human brain evolution: insights from microarrays. Nat Rev Genet. 2004, 5 (11): 850-860. 10.1038/nrg1469.PubMedView ArticleGoogle Scholar
- Rajeevan MS, Ranamukhaarachchi DG, Vernon SD, Unger ER: Use of real-time quantitative PCR to validate the results of cDNA array and differential display PCR technologies. Methods (San Diego, Calif). 2001, 25 (4): 443-451.View ArticleGoogle Scholar
- Chuaqui RF, Bonner RF, Best CJ, Gillespie JW, Flaig MJ, Hewitt SM, Phillips JL, Krizman DB, Tangrea MA, Ahram M: Post-analysis follow-up and validation of microarray experiments. Nature genetics. 2002, 32 (Suppl): 509-514. 10.1038/ng1034.PubMedView ArticleGoogle Scholar
- Do JH, Choi DK: Normalization of microarray data: single-labeled and dual-labeled arrays. Molecules and cells. 2006, 22 (3): 254-261.PubMedGoogle Scholar
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes and immunity. 2005, 6 (4): 279-284. 10.1038/sj.gene.6364190.PubMedView ArticleGoogle Scholar
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. Journal of biotechnology. 1999, 75 (2-3): 291-295. 10.1016/S0168-1656(99)00163-7.PubMedView ArticleGoogle Scholar
- Antonov J, Goldstein DR, Oberli A, Baltzer A, Pirotta M, Fleischmann A, Altermatt HJ, Jaggi R: Reliable gene expression measurements from degraded RNA by quantitative real-time PCR depend on short amplicons and a proper normalization. Laboratory investigation; a journal of technical methods and pathology. 2005, 85 (8): 1040-1050. 10.1038/labinvest.3700303.PubMedView ArticleGoogle Scholar
- Fleige S, Pfaffl MW: RNA integrity and the effect on the real-time qRT-PCR performance. Molecular aspects of medicine. 2006, 27 (2-3): 126-139. 10.1016/j.mam.2005.12.003.PubMedView ArticleGoogle Scholar
- Fleige S, Walf V, Huch S, Prgomet C, Sehm J, Pfaffl MW: Comparison of relative mRNA quantification models and the impact of RNA integrity in quantitative real-time RT-PCR. Biotechnology letters. 2006, 28 (19): 1601-1613. 10.1007/s10529-006-9127-2.PubMedView ArticleGoogle 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 biology. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.PubMed CentralPubMedView ArticleGoogle Scholar
- Waxman S, Wurmbach E: De-regulation of common housekeeping genes in hepatocellular carcinoma. BMC genomics. 2007, 8: 243-10.1186/1471-2164-8-243.PubMed CentralPubMedView ArticleGoogle Scholar
- Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, 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-47.PubMed CentralPubMedView ArticleGoogle Scholar
- Rhinn H, Marchand-Leroux C, 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-62.PubMed CentralPubMedView ArticleGoogle Scholar
- Spinsanti G, Panti C, Lazzeri E, Marsili L, Casini S, Frati F, Fossi CM: Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies. BMC molecular biology. 2006, 7: 32-10.1186/1471-2199-7-32.PubMed CentralPubMedView ArticleGoogle 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 (6): 509-515. 10.1023/B:BILE.0000019559.84305.47.PubMedView 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 research. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Akilesh S, Shaffer DJ, Roopenian D: Customized molecular phenotyping by quantitative gene expression and pattern recognition analysis. Genome research. 2003, 13 (7): 1719-1727. 10.1101/gr.533003.PubMed CentralPubMedView ArticleGoogle Scholar
- Haller F, Kulle B, Schwager S, Gunawan B, von Heydebreck A, Sultmann H, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Analytical biochemistry. 2004, 335 (1): 1-9. 10.1016/j.ab.2004.08.024.PubMedView ArticleGoogle Scholar
- Popovici V, Goldstein DR, Antonov J, Jaggi R, Delorenzi M, Wirapati P: Selecting control genes for RT-QPCR using public microarray data. BMC bioinformatics. 2009, 10: 42-10.1186/1471-2105-10-42.PubMed CentralPubMedView ArticleGoogle Scholar
- Lee S, Jo M, Lee J, Koh SS, Kim S: Identification of novel universal housekeeping genes by statistical analysis of microarray data. Journal of biochemistry and molecular biology. 2007, 40 (2): 226-231.PubMedView ArticleGoogle Scholar
- de Jonge HJ, Fehrmann RS, De Bont ES, Hofstra RM, Gerbens F, Kamps WA, De Vries EG, Zee van der AG, te Meerman GJ, ter Elst A: Evidence based selection of housekeeping genes. PLoS ONE. 2007, 2 (9): e898-10.1371/journal.pone.0000898.PubMed CentralPubMedView ArticleGoogle Scholar
- Kwon MJ, Oh E, Lee S, Roh MR, Kim SE, Lee Y, Choi YL, In YH, Park T, Koh SS: Identification of novel reference genes using multiplatform expression data and their validation for quantitative gene expression analysis. PLoS ONE. 2009, 4 (7): e6162-10.1371/journal.pone.0006162.PubMed CentralPubMedView ArticleGoogle Scholar
- Ravni A, Bourgault S, Lebon A, Chan P, Galas L, Fournier A, Vaudry H, Gonzalez B, Eiden LE, Vaudry D: The neurotrophic effects of PACAP in PC12 cells: control by multiple transduction pathways. Journal of neurochemistry. 2006, 98 (2): 321-329. 10.1111/j.1471-4159.2006.03884.x.PubMedView ArticleGoogle Scholar
- Levi A, Biocca S, Cattaneo A, Calissano P: The mode of action of nerve growth factor in PC12 cells. Molecular neurobiology. 1988, 2 (3): 201-226. 10.1007/BF02935346.PubMedView ArticleGoogle Scholar
- Vaudry D, Stork PJ, Lazarovici P, Eiden LE: Signaling pathways for PC12 cell differentiation: making the right connections. Science (New York, NY). 2002, 296 (5573): 1648-1649.View ArticleGoogle Scholar
- Fujita K, Lazarovici P, Guroff G: Regulation of the differentiation of PC12 pheochromocytoma cells. Environmental health perspectives. 1989, 80: 127-142. 10.2307/3430738.PubMed CentralPubMedView ArticleGoogle Scholar
- Sysi-Aho M, Katajamaa M, Yetukuri L, Oresic M: Normalization method for metabolomics data using optimal selection of multiple internal standards. BMC bioinformatics. 2007, 8: 93-10.1186/1471-2105-8-93.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu W, Dave N, Tseng GC, Richards T, Xing EP, Kaminski N: Comparison of normalization methods for CodeLink Bioarray data. BMC bioinformatics. 2005, 6: 309-10.1186/1471-2105-6-309.PubMed CentralPubMedView ArticleGoogle Scholar
- Ravni A, Vaudry D, Gerdin MJ, Eiden MV, Falluel-Morel A, Gonzalez BJ, Vaudry H, Eiden LE: A cAMP-dependent, protein kinase A-independent signaling pathway mediating neuritogenesis through Egr1 in PC12 cells. Molecular pharmacology. 2008, 73 (6): 1688-1708. 10.1124/mol.107.044792.PubMed CentralPubMedView ArticleGoogle Scholar
- Mark MD, Liu Y, Wong ST, Hinds TR, Storm DR: Stimulation of neurite outgrowth in PC12 cells by EGF and KCl depolarization: a Ca(2+)-independent phenomenon. The Journal of cell biology. 1995, 130 (3): 701-710. 10.1083/jcb.130.3.701.PubMedView ArticleGoogle Scholar
- Pacary E, Petit E, Bernaudin M: Concomitant inhibition of prolyl hydroxylases and ROCK initiates differentiation of mesenchymal stem cells and PC12 towards the neuronal lineage. Biochemical and biophysical research communications. 2008, 377 (2): 400-406. 10.1016/j.bbrc.2008.09.145.PubMedView ArticleGoogle Scholar
- Politis PK, Thomaidou D, Matsas R: Coordination of cell cycle exit and differentiation of neuronal progenitors. Cell cycle (Georgetown, Tex). 2008, 7 (6): 691-697.View ArticleGoogle Scholar
- Burkhalter J, Fiumelli H, Allaman I, Chatton JY, Martin JL: Brain-derived neurotrophic factor stimulates energy metabolism in developing cortical neurons. J Neurosci. 2003, 23 (23): 8212-8220.PubMedGoogle Scholar
- Salton SR: Neurotrophins, growth-factor-regulated genes and the control of energy balance. The Mount Sinai journal of medicine, New York. 2003, 70 (2): 93-100.PubMedGoogle Scholar
- da Silva JS, Dotti CG: Breaking the neuronal sphere: regulation of the actin cytoskeleton in neuritogenesis. Nature reviews. 2002, 3 (9): 694-704. 10.1038/nrn918.PubMedView ArticleGoogle Scholar
- Pak CW, Flynn KC, Bamburg JR: Actin-binding proteins take the reins in growth cones. Nature reviews. 2008, 9 (2): 136-147. 10.1038/nrn2236.PubMedView ArticleGoogle Scholar
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. BioTechniques. 2000, 29 (2): 332-337.PubMedGoogle 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 molecular biology. 2006, 7: 33-10.1186/1471-2199-7-33.PubMed CentralPubMedView ArticleGoogle Scholar
- Thorrez L, Van Deun K, Tranchevent LC, Van Lommel L, Engelen K, Marchal K, Moreau Y, Van Mechelen I, Schuit F: Using ribosomal protein genes as reference: a tale of caution. PLoS ONE. 2008, 3 (3): e1854-10.1371/journal.pone.0001854.PubMed CentralPubMedView ArticleGoogle Scholar
- Lafontaine DL, Tollervey D: The function and synthesis of ribosomes. Nat Rev Mol Cell Biol. 2001, 2 (7): 514-520. 10.1038/35080045.PubMedView ArticleGoogle Scholar
- Kenmochi N, Kawaguchi T, Rozen S, Davis E, Goodman N, Hudson TJ, Tanaka T, Page DC: A map of 75 human ribosomal protein genes. Genome research. 1998, 8 (5): 509-523.PubMedGoogle Scholar
- Robledo S, Idol RA, Crimmins DL, Ladenson JH, Mason PJ, Bessler M: The role of human ribosomal proteins in the maturation of rRNA and ribosome production. RNA (New York, NY). 2008, 14 (9): 1918-1929.View ArticleGoogle Scholar
- Faliks D, Meyuhas O: Coordinate regulation of ribosomal protein mRNA level in regenerating rat liver. Study with the corresponding mouse cloned cDNAs. Nucleic acids research. 1982, 10 (3): 789-801. 10.1093/nar/10.3.789.PubMed CentralPubMedView ArticleGoogle Scholar
- Kirn-Safran CB, Dayal S, Martin-DeLeon PA, Carson DD: Cloning, expression, and chromosome mapping of the murine Hip/Rpl29 gene. Genomics. 2000, 68 (2): 210-219. 10.1006/geno.2000.6283.PubMedView ArticleGoogle Scholar
- Wanzel M, Russ AC, Kleine-Kohlbrecher D, Colombo E, Pelicci PG, Eilers M: A ribosomal protein L23-nucleophosmin circuit coordinates Miz1 function with cell growth. Nat Cell Biol. 2008, 10: 1051-1061. 10.1038/ncb1764.PubMedView ArticleGoogle Scholar
- He H, Sun Y: Ribosomal protein S27L is a direct p53 target that regulates apoptosis. Oncogene. 2007, 26 (19): 2707-2716. 10.1038/sj.onc.1210073.PubMedView ArticleGoogle Scholar
- Jeon YJ, Kim IK, Hong SH, Nan H, Kim HJ, Lee HJ, Masuda ES, Meyuhas O, Oh BH, Jung YK: Ribosomal protein S6 is a selective mediator of TRAIL-apoptotic signaling. Oncogene. 2008, 27 (31): 4344-4352. 10.1038/onc.2008.73.PubMedView ArticleGoogle Scholar
- Lee NH, Weinstock KG, Kirkness EF, Earle-Hughes JA, Fuldner RA, Marmaros S, Glodek A, Gocayne JD, Adams MD, Kerlavage AR: Comparative expressed-sequence-tag analysis of differential gene expression profiles in PC-12 cells before and after nerve growth factor treatment. Proceedings of the National Academy of Sciences of the United States of America. 1995, 92 (18): 8303-8307. 10.1073/pnas.92.18.8303.PubMed CentralPubMedView ArticleGoogle Scholar
- Angelastro JM, Klimaschewski L, Tang S, Vitolo OV, Weissman TA, Donlin LT, Shelanski ML, Greene LA: Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling. Proceedings of the National Academy of Sciences of the United States of America. 2000, 97 (19): 10424-10429. 10.1073/pnas.97.19.10424.PubMed CentralPubMedView ArticleGoogle Scholar
- Too HP: Real time PCR quantification of GFRalpha-2 alternatively spliced isoforms in murine brain and peripheral tissues. Brain Res Mol Brain Res. 2003, 114 (2): 146-153. 10.1016/S0169-328X(03)00169-4.PubMedView ArticleGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical chemistry. 2009, 55 (4): 611-622. 10.1373/clinchem.2008.112797.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.