- Open Access
Validation of four reference genes for quantitative mRNA expression studies in a rat model of inflammatory injury
Molecular Painvolume 10, Article number: 55 (2014)
Real-time quantitative PCR (qPCR) is a technique frequently used to measure changes in mRNA expression. To ensure validity of experimental findings, it is important to normalize the qPCR data to reference genes that are stable and unaffected by the experimental treatment to correct for variability among samples. Unlike in some models of neuropathic pain, reference genes for models of inflammatory injury have not been validated. This study examined four candidate reference genes in an effort to identify and validate optimal genes for normalization of transcriptional changes occurring in the dorsal horn of the spinal cord and the rostral ventromedial medulla (RVM) following intraplantar injection of complete Freund’s adjuvant (CFA).
The expression of hypoxanthine phosphoribosyltransferase 1 (Hprt1), beta-actin (Actb), mitogen-activated protein kinase 6 (Mapk6), and beta-2-microglobulin (B2m) was quantified in the dorsal horn and RVM of rats four days or two weeks after intraplantar injection of CFA or saline. The range of expression levels among these four genes differed by as much as 16-fold within the dorsal horn and the RVM. All four of these reference genes were stably expressed in both tissues and did not differ between saline and CFA-treated animals. Analyses using the statistical algorithms in geNorm and NormFinder programs determined that Mapk6 was the most stable gene and recommended the combination of Mapk6 and Actb, or Mapk6 and Hprt1, in such experimental conditions.
This study validated the four genes Hprt1, Actb, Mapk6 or B2m and showed that any one or combination of two of them are good reference genes for normalization of mRNA expression in qPCR experiments in the spinal cord and RVM in the CFA model of inflammatory injury.
A better understanding of the molecular changes that occur during the development and maintenance of pain after inflammatory injury is fundamentally important for the development of effective therapies. Peripheral inflammatory injury activates signaling cascades in the peripheral and central nervous systems as well as in the immune system. These cascades lead to tissue swelling and increased sensitivity to both noxious and non-noxious stimuli that can persist long after the resolution of inflammation [1–4]. The underlying events entail transcriptional or translational changes in numerous genes and proteins including those involved in neurotransmitter release, receptor function and trafficking, subcellular signaling pathways, and regulation of ion channel expression and activity [1–4]. Identifying the nature of these changes can provide insights into the mechanisms of inflammatory injury and may elucidate new approaches for treating inflammatory pain.
Changes in gene expression that occur with peripheral inflammatory injury can be detected and measured in a sensitive and specific way using quantitative real-time polymerase chain reaction (qPCR) assays. Quantitative PCR is a frequently used technique that can be easily adapted to measure mRNA levels for any target protein and is the method of choice for absolute or relative quantification of mRNA expression. However, critics of qPCR maintain that it is often inadequately standardized and frequently inconsistent . Some technical challenges inherent to the technique include isolating high quality RNA, identifying an efficient reverse transcriptase (RT) enzyme to generate cDNA, designing efficient and specific primers to amplify the desired mRNAs, and normalizing results to adequately validated reference genes. In 2009, Bustin et al. published the MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments in an effort to standardize the information needed to ensure the relevance, accuracy, correct interpretation, and repeatability of qPCR experiments .
An essential component of qPCR is normalization of the target mRNA to a reference gene of interest in the same sample to control for variability associated with template input (amount of starting material) as well as RT and qPCR efficiencies. The reference gene mRNA should be stably expressed, and its abundance should show a strong correlation with the total amount of mRNA present in the sample . Control reference genes were initially termed “housekeeping genes” because the genes were historically chosen from cellular maintenance proteins that are ubiquitously expressed and whose mRNA was generally thought to have uniformly unchanging expression in different cells and under different conditions. More recent reports have documented that housekeeping gene expression can vary substantially (e.g. in tumor environments and other pathophysiological states) and that the choice of reference gene can significantly impact the conclusions of a study [7–9]. The development of statistical algorithms (e.g. geNorm, BestKeeper or NormFinder [10–12]) to help researchers determine the most stable reference genes and most appropriate combinations of reference genes has advanced the field.
Many studies of the molecular mechanisms that underlie the development and maintenance of pain after inflammatory injury have used qPCR to examine changes in gene expression in the peripheral and central nervous systems. Yet, to our knowledge, rigorous validation of appropriate reference genes has not been undertaken. This study quantified the expression of four potential reference genes, hypoxanthine phosphoribosyltransferase 1 (Hprt1), beta-actin (Actb), mitogen-activated protein kinase 6 (Mapk6), and beta-2-microglobulin (B2m), in the dorsal horn of the spinal cord and rostral ventromedial medulla (RVM) of the rat. These regions play important roles in the development and maintenance of nociception after peripheral inflammatory injury [1, 4, 13, 14].
Levels of Hprt1, Actb, Mapk6, and B2m mRNA in the dorsal horn and RVM of rats with persistent inflammatory injury
Transcripts for Hprt1, Actb, Mapk6, and B2m were readily detected in the dorsal horn of the spinal cord and the RVM of rats at four days (Figure 1A,C) and two weeks (Figure 1B,D) after intraplantar injection of saline or CFA. The four candidate reference genes were abundantly expressed in both tissues, with Cq values ranging from 21.11 to 21.85 for Hprt1, 17.81 to 18.42 for Actb, 20.56 to 21.54 for Mapk6, and 19.07 to 19.44 for B2m. The relative abundance of these mRNAs in both the dorsal horn and the RVM was Actb > B2m > Mapk6 > Hprt1 (P < 0.001).The levels of transcripts for all four candidate reference genes in CFA-treated rats were comparable to those in the corresponding saline-treated group in the dorsal horn (Figure 1A,B; P > 0.5 for all genes at both time points) and in the RVM (Figure 1C,D; P > 0.3 for all genes at both time points). The expression of each of the four reference genes in the dorsal horn at four days did not differ from that at two weeks in CFA- or saline-treated groups (P > 0.3, each gene). The expression of B2m in the RVM also did not differ at four days and two weeks in CFA- or saline-treated groups (P > 0.2). The expression levels of Hprt1, Actb, and MapK6 in the RVM at four days were statistically different from corresponding values at two weeks in saline- and CFA-treated rats (P < 0.05, each gene). However, the magnitude of the difference was less than 0.2 Cq and unlikely to be biologically significant.
geNorm analysis of reference gene stability
In the geNorm program, an M value under 1.0 is indicative of stable expression of mRNA . In the dorsal horn, all four genes had M values less than 0.4, demonstrating that all four have high expression stability (Figure 2A,B). Comparisons of pair-wise variation in the stability values  of the four genes in the dorsal horn indicated that Mapk6 and Hprt1 were the most stable transcripts at four days and that Mapk6 and Actb were the most stable at two weeks. When the time points were pooled in the interest of following expression changes over time, geNorm analysis indicated that Actb and Mapk6 were the two optimal genes in the dorsal horn (Figure 2C). Similarly, in the RVM, all four genes exhibited high stability of expression at four days and two weeks after injury (Figure 2D,E). Comparisons of pair-wise variation in the stability value indicated that at four days Actb and Hprt1 transcripts were the most stable, whereas at two weeks Hprt1 and Mapk6 were the most stable (Figure 2D,E). Pooling data for both time points indicated that Hprt1 and Mapk6 were the two most stable genes in the RVM (Figure 2F).
NormFinder analysis of reference gene stability
Expression stability for the four candidate reference genes was also calculated using NormFinder (Figure 3). In the dorsal horn, NormFinder analysis indicated that Hprt1 and Mapk6 were the best genes for normalization at four days and two weeks, respectively (Figure 3A,B). When the time points were pooled, NormFinder identified Mapk6 as the most stable gene for normalizing the expression data in the dorsal horn (Figure 3C). In the RVM, NormFinder analysis indicated that Actb and Hprt1 were the best genes for normalization at four days and two weeks, respectively (Figure 3D,E). When the data for both time points were pooled for analysis, NormFinder identified Mapk6 as the best gene for normalization in the RVM (Figure 3F). With respect to the best combination of two genes to normalize gene expression across time points, NormFinder identified Actb and Mapk6 in both the dorsal horn and RVM.NormFinder also includes an analysis of the optimal number of genes for normalization. In both dorsal horn and RVM, NormFinder identified the combination of Mapk6, Actb, and Hprt1 as optimal. However, the addition of one or two more genes caused at best a 0.06 decrease in the accumulated standard deviation, much less than the 0.15 that is considered a significant change (Figure 4). Thus, although the use of all three genes is predicted to be optimal, one gene would be acceptable for normalization of a target gene when resources are restricted.
As qPCR methodology has matured and its use has become more widespread, an increasing number of reports have demonstrated that the expression of classically used reference genes can vary substantially depending on the experimental conditions [7–9, 16]. Moreover, radically different conclusions can be reached depending on the reference gene used for normalization [16–19]. A PubMed literature search (March 20, 2014; search algorithm available from authors) identified 260 papers in English that described qPCR in rat models of nociception. The majority of the papers (83%, n = 216) normalized the expression of the target genes to a single reference gene. Glyceraldehyde 3-phosphate dehydrogenase (Gapdh) was used in 48% (n = 104), Actb in 25% (n = 54), Hprt1 in 4% (n = 8), and other housekeeping genes in the remaining 23% (n = 50) of the papers. Only 5% (n = 12) of the papers used more than one reference gene for normalization. When the search was restricted to studies using the CFA model of inflammatory injury in the rat, only 16 studies were identified. These 16 papers also normalized qPCR data with the same reference genes and in the same relative ratios as observed in the total 260 papers. To our knowledge, the present study is the first to specifically validate the use of reference genes for normalizing qPCR data in the CFA model of inflammatory injury.
Optimally, the transcript of the reference gene(s) selected for any given qPCR analysis will be expressed in an abundance comparable to that of the mRNA of interest . This study identified four different reference genes – Mapk6, Actb, Hprt1, and B2m – that can be used for the normalization of qPCR data in the CFA model of inflammatory injury in rats. These four genes were selected because they are broadly expressed, involved in different homeostatic functions (Table 1) and are not co-regulated . Moreover, two of these genes were also studied for validation as potential reference genes in the spared nerve injury  and spinal nerve ligation  models. The addition of our data allows the assessment of the stability of those genes over different pain models of chronic pain, possibly strengthening their utility as reference genes. The Cq values of the four genes ranged from 17 to 21.85 Cq, reflecting a 24 or a 16-fold difference in expression. Each gene was stably expressed within each time point: four days (subacute injury) and two weeks (chronic injury) after injury in CFA-treated animals compared to corresponding saline-treated animals. Interestingly, the Cq values for these genes in the RVM were very similar to those in the dorsal horn. The NormFinder algorithm indicated that Mapk6 was the reference gene with the least variability and therefore likely the best reference gene of the four in our experimental conditions. It also indicated that the combined use of three reference genes (Mapk6, Actb and Hprt1) was optimal for normalization in both the dorsal horn and the RVM. However, the gain in stability between the use of a single gene and the use of three genes is minimal, as shown by the very small changes in the accumulated SD (< 0.06) between one and three genes, indicating that any one of these three genes would be acceptable as a unique reference gene.
The need to report qPCR data in a standardized format is now widely accepted, and recommendations for the uniformity and reproducibility of qPCR experiments are listed in the MIQE Guidelines . These guidelines emphasize the need to control for sample-to-sample variation by normalization with reference genes. Normalization is necessary for reliable qPCR studies because the starting material, RNA extraction, RT efficiency, and qPCR efficiency can vary among experiments. Moreover, gene expression is highly tissue-specific and often varies based on the pathophysiological status of the organism or experimental treatment. Thus, it is imperative that a pilot study be conducted at the outset of any study to identify the optimal reference gene or combination of genes for that specific experiment. Unfortunately, very few validation studies for the use of reference genes exist in the pain field, and all are restricted to nerve injury models. Several studies have identified suitable qPCR reference genes in dorsal root ganglia in models of neuropathic pain [11, 21, 22]. For example, Mapk6 and Gapdh were identified as the two most stably expressed genes in an analysis of L4 and L5 dorsal root ganglia following L5 spinal nerve ligation . Hprt1 and 18S were validated as stable genes for normalizing expression levels in the dorsal root ganglia in the spared nerve injury model [22, 23], although this conclusion was not shared by Piller et al. . Seven commonly used reference genes, of which Actb was the most stable, were ranked and validated as good reference genes in spinal cord tissue in the spared nerve injury model of neuropathic pain . These reports support the conclusion that the optimal reference gene(s) will be specific to each experimental condition, and likely gender as well.
Gapdh was the single most frequently used reference gene in the literature search but was not included in this analysis for several reasons. Over the last decade, the transcription of Gapdh has been reported to be significantly regulated in different experimental settings and is variable in different tissues [24–26]. Moreover, the rat genome contains 329 Gapdh pseudogenes [27, 28], some of which are transcribed and have the same sequence as the active Gapdh transcript and, as such, can be detected by primers and amplified. For example, amplification using a published Gapdh primer pair : forward (5′-ACCACGAGAAATATGACAACTCCC) and reverse (5′-CCAAAGTTGTCATGGATGACC), designed to amplify a product of 100 base pairs, can give the same size amplicon from numerous transcripts. Performing an in silico PCR test by inserting the primer pair sequences for Gapdh in the UCSC genome website revealed 53 PCR products, each with 100 base pairs and 18 additional matches ranging from 75 to 120 base pairs, derived from genes on 16 different chromosomes. Other published primers for Gapdh yielded similar findings [20, 21]. Even if no genomic DNA contamination exists, and the RNA samples are treated with DNase, a more optimal qPCR assay design is one that does not amplify any pseudogenes. Similarly, Actb has pseudogenes , but the primers used in the present study were carefully selected to amplify Actb cDNA and not any of the pseudogene sequences. The UCSC in silico PCR search identified only one unique product in the rat genome with the Actb primers used here. Hprt1 also has pseudogenes, but the primers selected for this study are unique: one primer spans the splice junction of the gene, and when the primer pair is tested in silico, no genomic product was identified since the primers can only amplify the mRNA.
In summary, this study identified four reference genes that are stably expressed in the dorsal horn and the RVM four days and two weeks after intraplantar injection of CFA or saline. Each assay for the four reference genes was designed with primers that were unique for their target and would not amplify pseudogenes. Mapk6 was the most stable reference gene, although all genes tested met the criteria for a validated reference gene. We suggest that using any one or combination of two of these assays (e.g. Mapk6 and Actb, or Mapk6 and Hprt1) for normalization of data would yield accurate and reproducible results when studying mRNA expression in the spinal cord and RVM in the CFA model of inflammatory injury.
Animals and inflammatory injury model
Male Sprague–Dawley rats (200–350 g; Charles River, Raleigh, NC, USA) were housed in pairs in the University of Iowa Animal Care Facility in rooms with a 12 hr light/dark cycle with water and food provided ad libitum. All studies were conducted in accordance with the Guide for Care and Use of Laboratory Animals published by the National Institutes of Health following the guidelines of the International Association for the Study of Pain. The experiments were approved by the University of Iowa Animal Care and Use Committee (protocol 1107156), and care was taken to minimize the number of animals used and their suffering.
Persistent inflammatory injury was induced by intraplantar injection of CFA. Adult rats were lightly anesthetized with isoflurane, and the thickness of the hindpaw was measured with digital calipers. The plantar surface of the left hindpaw was injected with 150 μl of CFA (150 μg of Mycobacterium butyricum, Calbiochem, La Jolla, CA, USA) or sterile-filtered saline at pH 7.4. The rats were returned to their cages and singly housed for four days or two weeks, depending on the experiment. After the rats were euthanized, the thickness of the hindpaw was measured to verify the presence of inflammation. Four days after injection of CFA, the thickness of the ipsilateral hindpaw had increased from 6.1 ± 0.1 to 9.4 ± 0.4 mm (N = 12; P < 0.01). The thickness of the hindpaw in the saline-treated cohort was unchanged (6.2 ± 0.1 mm; N = 13; P > 0.3). Two weeks after injection of CFA, the thickness of the ipsilateral hindpaw had increased from 6.0 ± 0.1 to 9.0 ± 0.3 mm (N = 10, P < 0.01) whereas the thickness of the hindpaw in the saline-treated cohort was unchanged (6.0 ± 0.1 mm, N = 10, P > 0.4). Measures of nociception were not made to minimize stimulation.
qPCR – Quantitative Real-Time PCR Analysis
On the designated day, the rats were euthanized by CO2 inhalation and a 2-mm transverse slice of the brainstem containing the RVM was immediately isolated on ice and frozen on a platform on dry ice. To obtain the RVM, a 1.5-mm diameter tissue punch (Harris Unicore, Ted Pella Inc., Redding, CA), centered on the midline immediately above the pyramids, was removed from the frozen slice of brainstem tissue. The remainder of the slice was fixed in 10% formalin containing 30% sucrose to allow the verification of the site of the RVM tissue punch. The L4 and L5 portion of the spinal cord was removed, chilled on an ice-cold platform, and the ipsilateral dorsal horn was excised. The tissues were stored at -20°C in RNAlater™ (Ambion, Life Technologies, Carlsbad, CA) until RNA isolation.
DNA oligonucleotide primers were synthesized and purchased from Integrated DNA Technologies (Coralville, IA). The sequence of the forward and reverse primers for each of the four reference genes are listed in Table 1. The primers for Mapk6, B2m, Hprt1, and Actb are described elsewhere [21, 30, 31], and checked with the Primer 3 software (http://biotools.umassmed.edu) . Each qPCR assay consists of primers that hybridize to sequences that lie on different exons or span a splice junction, separated by one or more introns, such that when qPCR is conducted, only the cDNA sequence is amplified. The Ensembl database (http://useast.ensembl.org/) was used to examine the genomic structure of the gene and its transcripts. Special care was used in the selection of primers for Actb and Hprt1, which are known to have pseudogenes [12, 27, 28]. Contaminating genomic DNA, if present, will not amplify any product under the reaction conditions. Thus, only the specific mRNA targets will be measured. When the forward and reverse primer sequences are entered into the UCSC genome bioinformatics website (http://genome.ucsc.edu/) and tested in silico against rat genomic DNA, only a single product is identified for Mapk6 and Actb primers. In the case of the Hprt1 and B2m primers, no genomic match is found because one of primers spans a splice junction and can only hybridize with the correct cDNA sequence. In addition, the PCR amplicons for each qPCR assay were cloned into pSC-A, the PCR cloning vector, according to the manufacturer’s protocol (StrataClone PCR Cloning Kit, Agilent Technologies, Santa Clara, CA), and sequenced at the Iowa Institute of Human Genetics, Genomics Division. Results were aligned with the Genbank sequence for the intended mRNA using the NCBI BLAST program to confirm the specificity of each primer pair. When the primers are used in a qPCR assay, they amplify a unique species with a sharp melting curve, and no primer-dimer products are detected.
RNA isolation, RT, and qPCR
Total RNA was isolated from dorsal horn and RVM tissue according to the manufacturer’s protocol (RNeasy Lipid Tissue Mini Kit, Qiagen). Briefly, each tissue sample was homogenized in 1 ml of QIAzol lysis reagent. The lysate was extracted with chloroform, centrifuged, and the supernatant saved in a clean tube. The supernatant fraction was mixed with an equal volume of 70% ethanol and loaded on the column. Samples were treated with DNase I while on the column for 15 min at room temperature. After washing, the RNA was eluted from the column using RNase-free water. The concentration of total RNA was measured on a Nanodrop spectrophotometer (ND1000 3.8.1, Thermo Scientific, Wilmington, DE). The RNA integrity number was determined for 20% of the RNA samples as further RNA quality control and routinely showed values > 9.0 (out of 10) (Agilent Model 2100 Bioanalyzer, Santa Clara, CA). Reverse transcription was performed according to the manufacturer’s protocol using 100–600 ng of purified RNA and the SuperScript VILO cDNA synthesis kit (Life Technologies, Carlsbad, CA) in a 20 μl reaction volume. No-reverse transcriptase controls were also run for each RNA sample. The qPCR was performed with each well containing the cDNA product generated from 5 ng of input RNA, forward and reverse primers (12 nM), and the fluorogenic DNA-binding dye iQ™ SYBR® Green Supermix (Bio-Rad, Hercules, CA) in 20 μl. Reactions were performed in triplicate on a Bio-Rad CFX96 thermocycler (Bio-Rad, Hercules, CA). The cycle conditions were: 50°C for 2 min, 95°C for 10 min followed by 40 cycles of (95°C for 15 s, 60°C for 1 min and 72°C for 1 min), and then 95°C for 1 min and 55°C for 1 min. A thermal melting curve was generated from 55 to 95°C, at increments of 0.5°C for 10 s. No reverse transcriptase and no template controls for each primer pair were also tested in triplicate and did not amplify any product. Amplification efficiencies, calculated using the Bio-Rad CFX Manager 3.0 software, were similar for all primers and averaged at E = 101.2% ± 2.8, r2 ≥ 0.993 ± 0.002, slope = -3.30 ± 0.07. Cq values for each sample were calculated by the Bio-Rad CFX96 software (Bio-Rad, Hercules, CA).
geNorm and NormFinder
These two programs are available in the GenEX 6 software package (http://genex.gene-quantification.info/). The geNorm program calculates gene expression stability (M) for each reference gene as the average of a pairwise variation for the reference gene relative to the others. The stepwise analysis allows for a ranking of the genes according to the calculated M value. The lower the M value, the more stable the reference gene.
The data for NormFinder are organized in groups: e.g. A) four day RVM, B) two week RVM, C) four day dorsal horn, and D) two week dorsal horn. The program calculates intra-group and inter-group variations in the expression of the reference genes from which a stability value is generated for each reference gene. The candidate gene with the lowest stability value is considered the most stable reference gene. By calculation of the standard deviation of the data, the NormFinder program selects the best number of reference genes for an experiment.
Data are expressed as mean and S.E.M. The rank order of the expression of the four genes was analyzed by a one-way ANOVA followed by Bonferroni’s multiple comparison post-hoc test. A two-way ANOVA (factors: treatment and time) was used to compare the Cq values between saline- and CFA-treated rats, as well as between the four day and two week time points for each gene within each tissue. A P ≤ 0.05 was considered significant.
Complete Freund’s adjuvant
Glyceraldehyde 3-phosphate dehydrogenase
Hypoxanthine phosphoribosyltransferase 1
Mitogen-activated protein kinase 6
Real-time quantitative polymerase chain reaction
Rostral ventromedial medulla.
Latremoliere A, Woolf CJ: Central sensitization: a generator of pain hypersensitivity by central neural plasticity. J Pain 2009, 10: 895–926. 10.1016/j.jpain.2009.06.012
Peruzzotti-Jametti L, Donega M, Giusto E, Mallucci G, Marchetti B, Pluchino S: The role of the immune system in central nervous system plasticity after acute injury. Neuroscience 2014. http://dx.doi.org/10.1016/j.neuroscience.2014.04.036
Rahn EJ, Guzman-Karlsson MC, David Sweatt J: Cellular, molecular, and epigenetic mechanisms in non-associative conditioning: implications for pain and memory. Neurobiol Learn Mem 2013, 105: 133–150.
Xu Q, Yaksh TL: A brief comparison of the pathophysiology of inflammatory versus neuropathic pain. Curr Opin Anaesthesiol 2011, 24: 400–407. 10.1097/ACO.0b013e32834871df
Bustin SA: Why the need for qPCR publication guidelines?–The case for MIQE. Methods 2010, 50: 217–226. 10.1016/j.ymeth.2009.12.006
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 2009, 55: 611–622. 10.1373/clinchem.2008.112797
Bonefeld BE, Elfving B, Wegener G: Reference genes for normalization: a study of rat brain tissue. Synapse 2008, 62: 302–309. 10.1002/syn.20496
Marques TE, De Mendonca LR, Pereira MG, De Andrade TG, Garcia-Cairasco N, Paco-Larson ML, Gitai DL: Validation of suitable reference genes for expression studies in different pilocarpine-induced models of mesial temporal lobe epilepsy. PLoS ONE 2013, 8: e71892. 10.1371/journal.pone.0071892
Nelissen K, Smeets K, Mulder M, Hendriks JJ, Ameloot M: Selection of reference genes for gene expression studies in rat oligodendrocytes using quantitative real time PCR. J Neurosci Methods 2010, 187: 78–83. 10.1016/j.jneumeth.2009.12.018
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 Res 2004, 64: 5245–5250. 10.1158/0008-5472.CAN-04-0496
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.
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: RESEARCH0034.
Heinricher MM, Tavares I, Leith JL, Lumb BM: Descending control of nociception: Specificity, recruitment and plasticity. Brain Res Rev 2009, 60: 214–225. 10.1016/j.brainresrev.2008.12.009
Vanegas H, Schaible HG: Descending control of persistent pain: inhibitory or facilitatory? Brain Res Brain Res Rev 2004, 46: 295–309. 10.1016/j.brainresrev.2004.07.004
Hellemans J, Vandesompele J: Selection of Reliable Reference Genes for RT-qPCR Analysis. Methods Mol Biol 2014, 1160: 19–26. 10.1007/978-1-4939-0733-5_3
Caradec J, Sirab N, Keumeugni C, Moutereau S, Chimingqi M, Matar C, Revaud D, Bah M, Manivet P, Conti M, Loric S: ‘Desperate house genes’: the dramatic example of hypoxia. Br J Cancer 2010, 102: 1037–1043. 10.1038/sj.bjc.6605573
Bas A, Forsberg G, Hammarstrom S, Hammarstrom ML: Utility of the housekeeping genes 18S rRNA, beta-actin and glyceraldehyde-3-phosphate-dehydrogenase for normalization in real-time quantitative reverse transcriptase-polymerase chain reaction analysis of gene expression in human T lymphocytes. Scand J Immunol 2004, 59: 566–573. 10.1111/j.0300-9475.2004.01440.x
Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GA, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal Biochem 2005, 344: 141–143. 10.1016/j.ab.2005.05.022
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
Piller N, Decosterd I, Suter MR: Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search. BMC Res Notes 2013, 6: 266. 10.1186/1756-0500-6-266
Bangaru ML, Park F, Hudmon A, McCallum JB, Hogan QH: Quantification of gene expression after painful nerve injury: validation of optimal reference genes. J Mol Neurosci 2012, 46: 497–504. 10.1007/s12031-011-9628-x
Staaf S, Oerther S, Lucas G, Mattsson JP, Ernfors P: Differential regulation of TRP channels in a rat model of neuropathic pain. Pain 2009, 144: 187–199. 10.1016/j.pain.2009.04.013
Shankarappa SA, Tsui JH, Kim KN, Reznor G, Dohlman JC, Langer R, Kohane DS: Prolonged nerve blockade delays the onset of neuropathic pain. Proc Natl Acad Sci U S A 2012, 109: 17555–17560. 10.1073/pnas.1214634109
Butterfield DA, Hardas SS, Lange ML: Oxidatively modified glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and Alzheimer’s disease: many pathways to neurodegeneration. J Alzheimers Dis 2010, 20: 369–393.
Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun 2004, 313: 856–862. 10.1016/j.bbrc.2003.11.177
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
Liu YJ, Zheng D, Balasubramanian S, Carriero N, Khurana E, Robilotto R, Gerstein MB: Comprehensive analysis of the pseudogenes of glycolytic enzymes in vertebrates: the anomalously high number of GAPDH pseudogenes highlights a recent burst of retrotrans-positional activity. BMC Genomics 2009, 10: 480. 10.1186/1471-2164-10-480
Sun Y, Li Y, Luo D, Liao DJ: Pseudogenes as weaknesses of ACTB (Actb) and GAPDH (Gapdh) used as reference genes in reverse transcription and polymerase chain reactions. PLoS ONE 2012, 7: e41659. 10.1371/journal.pone.0041659
Wan G, Yang K, Lim Q, Zhou L, He BP, Wong HK, Too HP: Identification and validation of reference genes for expression studies in a rat model of neuropathic pain. Biochem Biophys Res Commun 2010, 400: 575–580. 10.1016/j.bbrc.2010.08.106
Hamity MV, Walder RY, Hammond DL: Increased neuronal expression of neurokinin-1 receptor and stimulus-evoked internalization of the receptor in the rostral ventromedial medulla of the rat after peripheral inflammatory injury. J Comp Neurol 2014, 522: 3037–3051. 10.1002/cne.23564
Langnaese K, John R, Schweizer H, Ebmeyer U, Keilhoff G: Selection of reference genes for quantitative real-time PCR in a rat asphyxial cardiac arrest model. BMC Mol Biol 2008, 9: 53. 10.1186/1471-2199-9-53
Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. Totowa, NJ: Humana Press; 2000.
We thank the Iowa Institute of Human Genetics, Genomics Division, for conducting the RIN analyses of RNA samples, for DNA sequencing, and for helpful discussions on normalization and real-time PCR. We thank Xiaomei Gu and Barry Matsumoto for help with the literature search. We thank Frank Jareczek and Chris Sande for helpful discussions. This work was supported by a John J. Bonica fellowship to A.S.W. and by grants R01DA06736 and R01DA23576 to D.L.H.
The authors declare that they have no competing interests.
All authors had full access to all data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. RYW, ASW and DLH originated the study concept and design and performed the statistical analysis. RYW, ASW, and DLH wrote the manuscript. RYW, ASW, MVH and SRW prepared samples. RYW carried out the qPCR experiments. DLH, ASW, MVH, BMP, and RYW analyzed and interpreted the data. All authors discussed and edited the manuscript. All authors read and approved the final manuscript.
Roxanne Y Walder, Anne-Sophie Wattiez contributed equally to this work.