In the investigation of the expression levels of target genes, reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is the most accurate and widely used method. However, a normalization step is a prerequisite to obtain accurate quantification results from RT-qPCR data. Therefore, many studies regarding the selection of reference genes have been carried out. Recently, these studies have involved large-scale gene analysis methods such as microarray and next generation sequencing. In our previous studies, we analyzed large amounts of transcriptome data from the cynomolgus monkey. Using a modification of this large-scale transcriptome sequencing dataset, we selected and compared 12 novel candidate reference genes (ARFGAP2, ARL1, BMI1, CASC3, DDX3X, MRFAP1, ORMDL1, RSL24D1, SAR1A, USP22, ZC3H11A, and ZRANB2) and 4 traditionally used reference genes (ACTB, GAPDH, RPS19, and YWHAZ) in 13 different whole-body tissues by the 3 well-known programs geNorm, NormFinder, and BestKeeper. Combined analysis by these 3 programs showed that ADP-ribosylation factor GTPase activating protein 2 (ARFGAP2), morf4 family associated protein 1 (MRFAP1), and ADP-ribosylation factor-like 1 (ARL1) are the most appropriate reference genes for accurate normalization. Interestingly, 4 traditionally used reference genes were the least stably expressed in this study. For this reason, selection of appropriate reference genes is vitally important, and large-scale analysis is a good method for finding new candidate reference genes. Our results could provide reliable reference gene lists for future studies on the expression of various target genes in the cynomolgus monkey.
References
[1]
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: RESEARCH0034.
[2]
Park SJ, Huh JW, Kim YH, Lee SR, Kim SH, et al. (2012) Selection of Internal Reference Genes for Normalization of Quantitative Reverse Transcription Polymerase Chain Reaction (qRT-PCR) Analysis in the Canine Brain and Other Organs. Mol Biotechnol. In press.
[3]
Yperman J, De Visscher G, Holvoet P, Flameng W (2004) Beta-actin cannot be used as a control for gene expression in ovine interstitial cells derived from heart valves. J Heart Valve Dis 13: 848–853.
[4]
Maccoux LJ, Clements DN, Salway F, Day PJ (2007) Identification of new reference genes for the normalisation of canine osteoarthritic joint tissue transcripts from microarray data. BMC Mol Biol 8: 62.
[5]
Cheng WC, Chang CW, Chen CR, Tsai ML, Shu WY, et al. (2011) Identification of reference genes across physiological states for qRT-PCR through microarray meta-analysis. PLoS One 6: e17347.
[6]
Zhou K, Zhou L, Lim Q, Zou R, Stephanopoulos G, et al. (2011) Novel reference genes for quantifying transcriptional responses of Escherichia coli to protein overexpression by quantitative PCR. BMC Mol Biol 12: 18.
[7]
Noriega NC, Kohama SG, Urbanski HF (2010) Microarray analysis of relative gene expression stability for selection of internal reference genes in the rhesus macaque brain. BMC Mol Biol 11: 47.
[8]
Chang Z, Ling C, Yamashita M, Welham NV (2010) Microarray-driven validation of reference genes for quantitative real-time polymerase chain reaction in a rat vocal fold model of mucosal injury. Anal Biochem 406: 214–221.
[9]
Demidenko NV, Logacheva MD, Penin AA (2011) Selection and validation of reference genes for quantitative real-time PCR in buckwheat (Fagopyrum esculentum) based on transcriptome sequence data. PLoS One 6: e19434.
[10]
De Pablos V, Barcia C, Martinez S, Gomez A, Ros-Bernal F, et al. (2009) MPTP administration increases plasma levels of acute phase proteins in non-human primates (Macaca fascicularis). Neurosci Lett 463: 37–39.
[11]
Register TC (2009) Primate models in women's health: inflammation and atherogenesis in female cynomolgus macaques (Macaca fascicularis). Am J Primatol 71: 766–775.
[12]
Race B, Meade-White KD, Miller MW, Barbian KD, Rubenstein R, et al. (2009) Susceptibilities of nonhuman primates to chronic wasting disease. Emerg Infect Dis 15: 1366–1376.
[13]
Higgs S, Ziegler SA (2010) A nonhuman primate model of chikungunya disease. J Clin Invest 120: 657–660.
[14]
Uchida A, Sasaguri H, Kimura N, Tajiri M, Ohkubo T, et al. (2012) Non-human primate model of amyotrophic lateral sclerosis with cytoplasmic mislocalization of TDP-43. Brain 135: 833–846.
[15]
Heo JH, Lee SR, Lee ST, Lee KM, Oh JH, et al. (2011) Spatial distribution of glucose hypometabolism induced by intracerebroventricular streptozotocin in monkeys. J Alzheimers Dis 25: 517–523.
[16]
Carlsson HE, Schapiro SJ, Farah I, Hau J (2004) Use of primates in research: a global overview. Am J Primatol 63: 225–237.
[17]
Gibbs RA, Rogers J, Katze MG, Bumgarner R, Weinstock GM, et al. (2007) Evolutionary and biomedical insights from the rhesus macaque genome. Science 316: 222–234.
[18]
Yan G, Zhang G, Fang X, Zhang Y, Li C, et al. (2011) Genome sequencing and comparison of two nonhuman primate animal models, the cynomolgus and Chinese rhesus macaques. Nat Biotechnol 29: 1019–1023.
[19]
Taylor K (2010) Clinical veterinarian's perspective of non-human primate (NHP) use in drug safety studies. J Immunotoxicol 7: 114–119.
[20]
Huh JW, Kim YH, Park SJ, Kim DS, Lee SR, et al. (2012) Large-scale transcriptome sequencing and gene analyses in the crab-eating macaque (macaca fascicularis) for biomedical research. BMC Genomics 13: 163.
[21]
Uno Y, Suzuki Y, Wakaguri H, Sakamoto Y, Sano H, et al. (2008) Expressed sequence tags from cynomolgus monkey (Macaca fascicularis) liver: a systematic identification of drug-metabolizing enzymes. FEBS Lett 582: 351–358.
[22]
Osada N, Hida M, Kusuda J, Tanuma R, Iseki K, et al. (2001) Assignment of 118 novel cDNAs of cynomolgus monkey brain to human chromosomes. Gene 275: 31–37.
[23]
Osada N, Hida M, Kusuda J, Tanuma R, Hirata M, et al. (2002) Cynomolgus monkey testicular cDNAs for discovery of novel human genes in the human genome sequence. BMC Genomics 3: 36.
[24]
Osada N, Hirata M, Tanuma R, Kusuda J, Hida M, et al. (2005) Substitution rate and structural divergence of 5'UTR evolution: comparative analysis between human and cynomolgus monkey cDNAs. Mol Biol Evol 22: 1976–1982.
[25]
Osada N, Hashimoto K, Kameoka Y, Hirata M, Tanuma R, et al. (2008) Large-scale analysis of Macaca fascicularis transcripts and inference of genetic divergence between M. fascicularis and M. mulatta. BMC Genomics 9: 90.
[26]
Osada N, Hirata M, Tanuma R, Suzuki Y, Sugano S, et al. (2009) Collection of Macaca fascicularis cDNAs derived from bone marrow, kidney, liver, pancreas, spleen, and thymus. BMC Res Notes 2: 199.
[27]
Andersen CL, Jensen JL, Orntoft TF (2004) Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 64: 5245–5250.
[28]
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP (2004) Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett 26: 509–515.
[29]
Rozen S, Skaletsky H (2000) Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 132: 365–386.
[30]
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, et al. (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55: 611–622.
[31]
Ahn K, Huh JW, Park SJ, Kim DS, Ha HS, et al. (2008) Selection of internal reference genes for SYBR green qRT-PCR studies of rhesus monkey (Macaca mulatta) tissues. BMC Mol Biol 9: 78.
