Selection of internal reference genes for SYBR green qRT-PCR studies of rhesus monkey (Macaca mulatta) tissues
© Ahn et al; licensee BioMed Central Ltd. 2008
Received: 18 March 2008
Accepted: 10 September 2008
Published: 10 September 2008
The rhesus monkey (Macaca mulatta) is a valuable and widely used model animal for biomedical research. However, quantitative analyses of rhesus gene expression profiles under diverse experimental conditions are limited by a shortage of suitable internal controls for the normalization of mRNA levels. In this study, we used a systematic approach for the selection of potential reference genes in the rhesus monkey and compared their suitability to that of the corresponding genes in humans.
Eight housekeeping genes (HKGs) (GAPDH, SDHA, ACTB, RPL13A, RPL32, UBA52, PGK1Y, and YWHAZ) from rhesus monkeys and humans were selected to test for normalization of expression levels in six different tissue types (brain, colon, kidney, liver, lung, and stomach). Their stability and suitability as reference genes were validated by geNorm, NormFinder and BestKeeper programs. Intriguingly, RPL13A and RPL32 were selected as ideal reference genes only in rhesus monkeys.
The results clearly indicated the necessity of using different reference genes for normalization of expression levels between rhesus monkeys and humans in various tissues.
The rhesus monkey (Macaca mulatta) has been one of the most valuable experimental model species for biomedical studies in various research fields including microbiology, neuroscience, and biochemistry. Recently, the genome sequence of the rhesus monkey was determined by the Rhesus Macaque Genome Sequencing and Analysis Consortium . The consortium revealed that the average sequence similarity between humans and rhesus monkeys is 93%. Approximately twenty thousand genes were annotated or predicted by the human-based comparative analysis. In order to conduct biomedical research in areas such as infectious disease, cancer research, and drug development using various monkey tissues, quantitative gene expression analysis is essential.
Quantitative real-time PCR (qRT-PCR) is a recently developed and popularized molecular technique that can be used for gene expression studies that require highly sensitive and accurate quantification of mRNA levels from various samples . These PCR-based analyses allow investigators to analyze a limited amount of material of interest while guaranteeing accurate quantification of target mRNAs. However, to achieve optimal results in qRT-PCR analysis, certain steps are indispensable including: (1) quality control of the material of interest and the primer sets, (2) determination of the PCR efficiency, and (3) selection of suitable internal control genes for normalization . The availability of advanced biomaterials including sophisticated PCR instrumentation, kits for RNA isolation from various materials, and highly purified Taq polymerase, has overcome many of these difficult steps with the exception of identifying appropriate internal control genes . Different expression profiling results from the same target genes can be obtained by using different HKGs as internal references . Therefore, selection of proper HKGs is important for the accurate analysis and normalization of data generated by real-time qRT-PCR .
Before the systematic selection of reference genes, researchers interested in quantitative expression analysis have used empirically determined housekeeping genes such as ACTB and GAPDH. However, recent studies have provided solid evidence that transcription levels of housekeeping genes vary between cell types [6–11], developmental stages [12, 13], and experimental conditions [14, 15], individuals [16–20]. Therefore, stable reference gene selection procedures have to be undertaken before quantitative mRNA expression analysis can be conducted.
In this study, three different specific tools (geNorm, Normfinder, and BestKeeper) that can be used for the validation of the stability of selected HKGs (GAPDH, SDHA, ACTB, RPL13A, RPL32, UBA52, PGK1Y, and YWHAZ) were investigated in rhesus monkey and human tissues (brain, colon, kidney, liver, lung, stomach) using real-time qRT-PCR with SYBR green.
Results and discussion
Reference gene and qRT-PCR primer selection
GenBank accession numbers and functions of 8 reference genes in rhesus monkeys and humans
Cytoskeletal structural protein
Ribosomal protein L13A
Member of ribosome protein
60S ribosomal protein L32
Member of ribosome protein
Signal transduction by binding to phosphorilated serine residue on a variaty of signaling molecules
Third enzyme of the heme biosynthetic pathway and catalyzes the head to tail condensation of four porphobilinogen molecules into the linear hydroxymethylbilane.
TATA box binding protein
Transcription initiation from RNA polymerase II promoter
Hypoxanthine phosphoribosyltransferase 1
Metabolic salvage of purines in mammals
Details of primers used for quantitative RT-PCR
Forward Primer Sequence [5' → 3']
Position in cDNA
Reverse Primer Sequence [5' → 3']
Position in cDNA
5, 6th Exon
2, 3th Exon
The expression levels of HKGs were analyzed using the relative quantification (delta-Ct method) . The amplification efficiencies and correlation coefficients (R2) of the eight HKGs were generated using the slopes of the standard curves obtained by serial dilutions. Standard curves with a ten-fold dilution series were used to calculate the amplification efficiency (additional files 1 and 2). Similar amplification efficiencies for the corresponding human and rhesus monkey reference genes are a prerequisite for accurate real-time RT-PCR quantification. The amplification efficiency was calculated by the formula: efficiency (%) = (10(-1/slope) -1)*100. The efficiency range of the qRT-PCR amplifications for all of the genes tested was between 80% and 115%. Overall, similar efficiency levels were observed between human and rhesus samples (additional file 3). The amplification specificity for each qRT-PCR analysis was confirmed by melting curve analysis (additional files 4 and 5).
The normalization step must be conducted very carefully for accurate gene expression level analysis. However, the selection of stable HKGs for normalization is very complicated. geNorm is a useful tool for the selection of the most stable internal control genes using the principle that the expression ratio of two ideal internal control genes is identical in all tested samples .
As the most stable reference genes, RPL32 and RPL13A yielded a V value of 0.200 in the rhesus monkey (Fig. 1b). That value indicated that it would not be necessary to include additional reference genes in qRT-PCR of rhesus monkey tissues. However, in the case of humans, the addition of a third and fourth gene significantly affected the reliability of the selected HKGs (Fig. 1d).
The results of the geNorm tool indicated different results between humans and rhesus monkeys. RPL32 and RPL13A in the rhesus monkey and GAPDH, HMBS, ACTB and RPL32 in humans are the most stable HKGs for the accurate calculation of normalization factors.
Candidate reference genes for normalization and the best combination of two genes listed according to their expression stability calculated by the NormFinder program
Stablility value of rhesus monkey
Stablility value of human
Results from Bestkeeper descriptive statistical analysis
geo Mean [CP]
ar Mean [CP]
std dev [± CP]
CV [% CP]
Based on the Bestkeeper program, the RPL13A gene in the rhesus monkey and the RPL32 gene in the human were selected as the best reference genes.
To validate our results from the rhesus monkey, one unknown rhesus monkey was tested using same reference genes. And results were analyzed for the selection of stable HKGs using three authorized programs (geNORM, NormFinder, and BestKeeper) (additional files 6, 7 and 8). Although, slight differences were observed, overall patterns were coincided with our results of RPL13A and RPL32 as ideal reference genes in rhesus monkeys.
In this study, we investigated for the first time the most reliable HKGs for the normalization of real-time qRT-PCR data in rhesus monkey tissues (brain, colon, kidney, liver, lung, stomach) using three programs (geNorm, Normfinder, and Bestkeeper). The results from the rhesus monkey were compared to those of humans.
Our results indicated that the two most stable genes, RPL13A and RPL32, covering a broad expression range, could be used as reference genes for relative gene quantification and normalization purposes in gene profiling studies of the rhesus monkey as determined by three reference gene identification programs.
In humans, RPL32, HMBS and GAPDH were found to be suitable reference genes. In the case of RPL32, two programs (Bestkeepe r and Normfinder) identified this gene as a reliable reference gene. However, geNorm indicated that RPL32 was an unsuitable reference gene by virtue of having a high M value. Likewise, the GAPDH gene was rejected by the Bestkeeper program as a result of having a high SD [± Ct] value. Finally, HMBS was ranked as the sixth best gene by the Normfinder program.
Intriguingly, the GAPDH and ACTB genes are the most widely used reference genes in human studies. However, GAPDH and ACTB are ranked as the least stable genes in the rhesus monkey. Genomic and environmental differences between humans and rhesus monkeys might result in the selection of different suitable HKGs in the six different tissues tested. Thus, different reference genes (RPL32 and RPL13A) should be used for the analysis of rhesus monkey tissue than those used for human tissue analysis
Materials and methods
Total RNA from rhesus macaque tissues (colon, liver, brain, lung, stomach, and kidney) and human total RNA (colon, liver, brain, lung, stomach, and kidney) were purchased from Clontech and Ambion, respectively. For the validation of our results for reference gene selection, corresponding unknown rhesus macaque tissues (colon, liver, brain, lung, stomach, and kidney) were provided by the National Primate Research Center (NPRC) of Korea. Total RNA was quantified using a NanoDrop® ND-1000 UV-Vis Spectrophotometer.
Rhesus macaque and human tissue cDNA synthesis
To eliminate DNA contamination from the total RNA samples, Turbo DNA-freeTM (Ambion) was used. A no-RT control was also amplified to confirm the absence of DNA contamination. M-MLV (Moloney-Murine Leukemia Virus) reverse transcriptase with an annealing temperature of 42°C was used for the reverse transcription reaction with RNase inhibiter (Promega).
Candidate reference genes and primers for real-time quantitative RT-PCR
Rhesus monkey reference genes were identified by comparison with human HKGs using the NCBI (National Center for Biotechnology Information) database and the UCSC (University of California Santa Cruz) genome browser (Table 1). Their sequences were used to design primers design using Primer3 software . BLAST searches were performed to confirm the total gene specificity of the primer sequences, and the results showed the absence of multi-locus matching at the primer site. All primers except for RPL32 spanned at least one intron to avoid inaccuracy due to genomic DNA contamination in RNA samples.
Real-time quantitative RT-PCR
SYBR green qRT-PCR was performed on a Rotor Gene 3000 (Corbett Research). In each run, 1 ul of cDNA was used as template for amplification per reaction. The sample was added to 19 ul of reaction mixture containing 7 ul H2O, 10 ul QuantiTect® SYBR® Green PCR Master Mix (Qiagen) and 1 ul forward and reverse primers (Table 2). Real-time qRT-PCR amplification of the HKGs was carried out for 50 cycles of 94°C for 10 sec, 58°C for 15 sec, and 72°C for 15 sec. The temperature range for the analysis of melting curves was 55°C to 99°C over 30 sec. Three independent experiments were performed.
Characterization of analysis programs
The geNorm  program provides a measure of gene expression stability (M), which is the mean pair-wise variation between an individual gene and all other tested control genes. This method differs from model-based approaches by comparing genes based on the similarity of their expression profiles. Ct values were converted to scale expression quantities using the delta-Ct method and entered into geNorm, which then ranks the genes based on M values where the gene with the most stable expression has the lowest value. NormFinder  focuses on finding the two genes with the least intra- and inter-group expression variation. A BestKeeper  index was created using the geometric mean of each candidate gene's Ct values. This index was then compared to each individual candidate housekeeping gene by pair-wise correlation analyses, with each combination assigned a value for the Pearson correlation coefficient (r) and the probability (p).
This work was supported by the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korea government (MOST) (No. R01-2007-000-20035-0).
- Rhesus Macaque Genome Sequencing and Analysis Consortium: Evolutionary and Biomedical Insights from the Rhesus Macaque Genome. Science 2007, 316: 222. 10.1126/science.1139247.View 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 Biol 2002, 3: 34. 10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nat Protoc 2006, 1(3):1559-82. 10.1038/nprot.2006.236.View ArticlePubMedGoogle Scholar
- Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative realtime PCR. Biochem Biophys Res Commun 2004, 313: 856-862. 10.1016/j.bbrc.2003.11.177.View ArticlePubMedGoogle Scholar
- Sellars MJ, Vuocolo T, Leeton LA, Coman GJ, Degnan BM, Preston NP: Real-time RT-PCR quantification of Kuruma shrimp transcripts: a comparison of relative and absolute quantification procedures. J Biotechnol 2007, 129(3):391-9. 10.1016/j.jbiotec.2007.01.029.View 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-9.View ArticlePubMedGoogle Scholar
- Mogal A, Abdulkadir SA: Effects of Histone Deacetylase Inhibitor (HDACi); Trichostatin-A (TSA) on the expression of housekeeping genes. Mol Cell Probe 2006, 20: 81-86. 10.1016/j.mcp.2005.09.008.View ArticleGoogle Scholar
- Li YP, Bang DD, Handberg KJ, Jorgensen PH, Zhang MF: Evaluation of the suitability of six host genes as internal control in real-time RT-PCR assays in chicken embryo cell cultures infected with infectious bursal disease virus. Vet Microbiol 2005, 110(3-4):155-165. 10.1016/j.vetmic.2005.06.014.View ArticlePubMedGoogle Scholar
- Ohl F, Jung M, Radoniæ A, Sachs M, Loening SA, Jung K: Identification and Validation of Suitable Endogenous Reference Genes for Gene Expression Studies of Human Bladder Cancer. J Urology 2006, 175: 1915-1920. 10.1016/S0022-5347(05)00919-5.View ArticleGoogle Scholar
- Fischer M, Skowron M, Berthold F: Reliable transcript quantification by real-time reverse transcriptase-polymerase chain reaction in primary neuroblastoma using normalization to averaged expression levels of the control genes HPRT1 and SDHA . J Mol Diagn 2005, 7(1):89-96.PubMed CentralView ArticlePubMedGoogle Scholar
- Goidin D, Mamessier A, Staquet MJ, Schmitt D, Berthier-Vergnes O: Ribosomal 18S RNA prevails over glyceraldehyde-3-phosphate dehydrogenase and beta-actin genes as internal standard for quantitative comparison of mRNA levels in invasive and noninvasive human melanoma cell subpopulations. Anal Biochem 2001, 295(1):17-21. 10.1006/abio.2001.5171.View ArticlePubMedGoogle Scholar
- Goossens K, Van Poucke M, Van Soom A, Vandesompele J, Van Zeveren A, Peelman LJ: Selection of reference genes for quantitative real-time PCR in bovine preimplantation embryos. BMC Dev Biol 2005, 5: 2. 10.1186/1471-213X-5-27.View ArticleGoogle Scholar
- Kuijk EW, du Puy L, van Tol HT, Haagsman HP, Colenbrander B, Roelen BA: Validation of reference genes for quantitative RT-PCR studies in porcine oocytes and preimplantation embryos. BMC Dev Biol 2007, 7: 58. 10.1186/1471-213X-7-58.PubMed CentralView ArticlePubMedGoogle Scholar
- Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Anal Biochem 2002, 309: 293-300. 10.1016/S0003-2697(02)00311-1.View ArticlePubMedGoogle Scholar
- Zhang X, Ding L, Sandford AJ: Selection of reference genes for gene expression studies in human neutrophils by real-time PCR. BMC Mol Biol 2005, 6: 4. 10.1186/1471-2199-6-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Overbergh L, Kyama CM, Valckx D, Debrock S, Mwenda JM, Mathieu C, Hooghe TD: Validation of real-time RT-PCR assays for mRNA quantification in baboons. Cytokine 2005, 31(6):454-458. 10.1016/j.cyto.2005.07.002.View ArticlePubMedGoogle Scholar
- Nygard AB, Jorgensen CB, Cirera S, Fredholm M: Selection of reference genes for gene expression studies in pig tissues using SYBR green qPCR. BMC Mol Biol 2007, 8: 67. 10.1186/1471-2199-8-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Brinkhof B, Spee B, Rothuizen J, Penning LC: Development and evaluation of canine reference genes for accurate quantification of gene expression. Anal Biochem 2006, 356: 36-43. 10.1016/j.ab.2006.06.001.View ArticlePubMedGoogle 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 Mol Biol 2006, 7: 32. 10.1186/1471-2199-7-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Olsvik PA, Lie KK, Jordal AE, Nilsen TO, Hordvik I: Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol Biol 2005, 6: 21. 10.1186/1471-2199-6-21.PubMed CentralView ArticlePubMedGoogle Scholar
- Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC PlantBiol 2004, 4: 14. 10.1186/1471-2229-4-14.View ArticleGoogle Scholar
- geNorm Software[http://medgen.ugent.be/~jvdesomp/genorm/]
- NormFinder Software[http://www.mdl.dk/publicationsnormfinder.htm]
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology Totowa: Humana Press; 2000, 365-386. [http://frodo.wi.mit.edu/]Google 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-33.PubMed CentralView ArticlePubMedGoogle Scholar
- Falk O, Monika J, Chuanliang X, Carsten S, Anja R, Mick B, Andreas N, Glen K, Stefan AL, Aleksandar R, Klaus J: Gene expression studies in prostate cancer tissue: which reference gene should be selected for normalization? J Mol Med 2005, 83: 1014-1024. 10.1007/s00109-005-0703-z.View ArticleGoogle Scholar
- Erkens T, Van Poucke M, Vandesompele J, Goossens K, Van Zeveren A, Peelman LJ: Development of a new set of reference genes for normalization of real-time RT-PCR data of porcine backfat and longissimus dorsi muscle, and evaluation with PPARGC1A . BMC Biotechnol 2006, 6: 41. 10.1186/1472-6750-6-41.PubMed CentralView 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.47.View ArticlePubMedGoogle 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.