Selection of internal reference genes for SYBR green qRT-PCR studies of rhesus monkey (Macaca mulatta) tissues
BMC Molecular Biology volume 9, Article number: 78 (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
For the selection of suitable reference genes from the human and rhesus monkey genomes, the eight most commonly used and available HKGs were selected: Hydroxymethylbilane synthase (HMBS), Beta-actin (ACTB), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Ribosomal protein 13 (RPL13A), Ribosomal protein 32 (RPL32), Tyrosine 3-monooxygenase (YWHAZ), Hypoxanthine phosphoribosyltransferase 1(HPRT1), and TATA box-binding protein (TBP). Rhesus monkey HKGs are identified by comparison of orthologous loci of the human HKGs. Summarized sequence information for the human and rhesus monkey genes is shown in Table 1. This gene information was also used to design qRT-PCR primers (Table 2).
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 .
Using the geNorm program, the eight HKGs selected for normalization were ranked according to their M values. The M value is the average pair-wise variation of a particular gene with all other control genes. The M values for RPL13A, RPL32, TBP, HMBS, ACTB, YWHAZ, HPRT1, and GAPDH in rhesus were lower than 1.5.RPL32 and RPL13A displayed the most stability of all eight genes according to geNorm analysis. Both genes encode components of the 60S ribosomal subunit that is involved in protein biosynthesis (Fig. 1a). In contrast to rhesus monkeys, human tissue analysis indicated that GAPDH and HMBS were the most stable genes (Fig. 1c). The M values of all tested human genes except for HPRT1 were lower than 1.5. To investigate the optimal number of reference genes, it was determined whether the stepwise inclusion of less stable genes significantly affected the normalization factors, considering the term V (Fig. 1).
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.
NormFinder is another reference tool for identifying the ideal normalization genes among various candidate HKGs . NormFinder ranks the various candidate reference genes according to their expression variation between inter and intra groups. This analysis identified HMBS and RPL32 as the most stable genes in rhesus monkeys and humans, respectively (Table 3). These rankings appear to be generally different from the results of geNorm. However, RPL32, GAPDH and ACTB in humans and HMBS and RPL32 in rhesus monkeys are still ranked as the most stable reference genes. Although geNorm and NormFinder did not identify the exact same reference genes, GAPDH in humans and HMBS in rhesus monkeys are still considered candidate reference genes.
Intriguingly, human HMBS showed opposite results between geNorm and Normfinder program. In the geNorm software, the most important factor is gene expression stability (M). The principle of geNorm is that the expression ratio of two ideal internal control genes is identical in all tested samples . Thus, the geNorm program did not consider the other HKGs. However, Normfinder program focuses on the two genes with the least intra-and inter-group expression variation. In other words, the Normfinder program considered other HKGs for the selection of good HKGs. For example, high rate of expression variation of HMBS among HKGs could result in the discordance between Normfinder and geNorm (Fig. 2).
The BestKeeper program index was created using the geometric mean of each candidate gene's Ct values . Initial analysis of the data, based on the inspection of raw Ct values, calculated variations for all the housekeeping genes in the six samples (SD [± Ct] and CV [%Ct]), and showed an overall stability in gene expression. Genes with SD values higher than 1 were considered inconsistent. ACTB, HMBS, RPL13A, RPL32, TBP, YWHAZ showed a SD value lower than 1 in the rhesus monkey. Among them, RPL13A (SD = 0.22) had the lowest Ct value variation. GAPDH had the highest Ct value variation. In human, HMBS and RPL32 exhibit SD values lower than 1. RPL32 had the lowest Ct value variation and HPRT1 had the highest Ct value variation (Table 4). In the case of HMBS and RPL32, BestKeeper also showed a different result compared to Normfinder and geNorm program, respectively. Because the variation of Ct value of HMBS and RPL32 is ranked very low, the BestKeeper program indicated HMBS and RPL32 as good candidate HKGs. However, compared to BestKeeper, geNorm is considered a pair-wise variation between the gens as the most important factor. Thus geNorm did not indicate one gene as a best reference gene but indicated a gene pair.
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).
Rhesus Macaque Genome Sequencing and Analysis Consortium: Evolutionary and Biomedical Insights from the Rhesus Macaque Genome. Science 2007, 316: 222. 10.1126/science.1139247.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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/]
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.
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.
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.
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.
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).
Electronic supplementary material
Additional file 1: Standard curves for calculation of PCR efficiency and quantification according to reference genes in the rhesus monkey. Amplification of 10-fold serial dilutions of the plasmid standard ranging from 100 to 105 copies per reaction was carried out in duplicate. (PDF 66 KB)
Additional file 2: Standard curves for calculation of PCR efficiency and quantification according to reference genes in humans. Amplification of 10-fold serial dilutions of the plasmid standard ranging from 100 to 105 copies per reaction was carried out in duplicate. (PDF 54 KB)
Additional file 4: Melting curve analysis. Melting curve analysis of 8 different primer sets for reference genes in the rhesus monkey. (PDF 155 KB)
Additional file 5: Melting curve analysis. Melting curve analysis of 8 different primer sets for reference genes in humans. (PDF 176 KB)
Additional file 6: Selection of the most suitable reference genes for normalization using geNorm analysis using the unknown rhesus monkey samples for validation. (PDF 23 KB)
Additional file 7: Candidate reference genes for normalization and the best combination of two genes listed according to their expression stability calculated by the NormFinder program using the unknown rhesus monkey samples for validation. (XLS 18 KB)
Additional file 8: Results from Bestkeeper descriptive statistical analysis using the unknown rhesus monkey samples for validation. (XLS 17 KB)
About this article
Cite this article
Ahn, K., Huh, JW., Park, SJ. et al. Selection of internal reference genes for SYBR green qRT-PCR studies of rhesus monkey (Macaca mulatta) tissues. BMC Molecular Biol 9, 78 (2008). https://doi.org/10.1186/1471-2199-9-78