Selection of suitable internal controls for gene expression normalization in rats with spinal cord injury.

There is a lack of systematic research on the expression of internal control genes used for gene expression normalization in real-time reverse transcription polymerase chain reaction in spinal cord injury research. In this study, we used rat models of spinal cord hemisection to analyze the expression stability of 13 commonly applied reference genes: Actb, Ankrd27, CypA, Gapdh, Hprt1, Mrpl10, Pgk1, Rictor, Rn18s, Tbp, Ubc, Ubxn11, and Ywhaz. Our results show that the expression of Ankrd27, Ubc, and Tbp were stable after spinal cord injury, while Actb was the most unstable internal control gene. Ankrd27, Ubc, Tbp, and Actb were consequently used to investigate the effects of internal control genes with differing stabilities on the normalization of target gene expression. Target gene expression levels and changes over time were similar when Ankrd27, Ubc, and Tbp were used as internal controls but different when Actb was used as an internal control. We recommend that Ankrd27, Ubc, and Tbp are used as internal control genes for real-time reverse transcription polymerase chain reaction in spinal cord injury research. This study was approved by the Administration Committee of Experimental Animals, Jiangsu Province, China (approval No. 20180304-008) on March 4, 2018.

Chinese Library Classification No. R446; R74; R318

Introduction

Spinal cord injury impairs sensation and muscle functions and often leads to paraplegia or quadriplegia. Spinal cord injury has complex pathophysiology involving primary injury and subsequent secondary injury (Silva et al., 2014). Primary spinal cord injury refers to the initial mechanical injury and occurs within 2 hours, whereas secondary injury refers to a series of pathological processes, including apoptosis, hypoxia, edema, and inflammation, and can last weeks or even months (Anwar et al., 2016; Ahuja et al., 2017a). Secondary injury aggravates cell death and tissue damage and expands the zone and severity of tissue damage (Ahuja et al., 2017b; Zhang et al., 2020). These multifaceted injury events involve many types of cells, including neurons, glial cells (such as astrocytes, oligodendrocytes, and microglia), and inflammatory cells (such as macrophages, neutrophils, and T cells) (O’Shea et al., 2017; Huang et al., 2019). Recently, with the development of high-throughput analysis, gene profiles after spinal cord injury have been revealed, facilitating the understanding of the underlying molecular basis of and the discovery of potential targets for the treatment of spinal cord injury (Chen et al., 2013; Yu et al., 2019; Gong et al., 2020). However, for the precise determination of gene abundance, appropriate internal control genes should be used to normalize gene expression during quantification.

In quantitative real-time reverse transcription polymerase chain reaction (RT-PCR) experiments, internal controls are used to account for variations in gene expression determination caused by variations in experimental procedures and RNA sample quantity and quality (Vandesompele et al., 2002). Many reference genes have been identified as suitable internal controls in diverse species, tissues, and physiological and pathological processes (Benak et al., 2019; Herath et al., 2020; Liu et al., 2020; Xia et al., 2021). Nevertheless, there are few consensuses on proper internal controls for studies of spinal cord injury and repair. Therefore, we investigated the expression stability of 13 reference genes in normal or injured spinal cord with the aim to select optimal internal control genes following spinal cord injury. The 13 reference genes were Actb (beta-actin), Ankrd27 (ankyrin repeat domain 27), CypA (cyclophilin A), Gapdh (glyceraldehyde-3-phosphate dehydrogenase), Hprt1 (hypoxanthine phosphoribosyltransferase 1), Mrpl10 (mitochondrial ribosomal protein L10), Pgk1 (phosphoglycerate kinase 1), Rictor (RPTOR independent companion of MTOR, complex 2), Rn18s (18S ribosomal RNA), Tbp (TATA box binding protein), Ubc (ubiquitin C), Ubxn11 (UBX domain protein 11), and Ywhaz (tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta).

Materials and Methods

Animals

The animal experiments were ethically approved by the Administration Committee of Experimental Animals, Jiangsu Province, China (approval No. 20180304-008) on March 4, 2018. Surgical procedures were performed following the Institutional Animal Care Guidelines of Nantong University, Nantong, China. All experiments were designed and reported according to the Animal Research: Reporting of In Vivo Experiments (ARRIVE) guidelines.

A total of 45 adult female Sprague-Dawley rats (specific-pathogen-free; 2 months old; 180–220 g) were purchased from the Animal Center of Nantong University (license no. SCXK [Su] 2019-0001 and SYXK [Su] 2017-0046). The rats were randomized into experiment and sham groups. Rats were anesthetized by intraperitoneal injection of mixed narcotics (MilliporeSigma, Burlington, MA, USA; Merck KGaA, Darmstadt, Germany) containing 85 mg/kg trichloroacetaldehyde monohydrate, 42 mg/kg magnesium sulfate, and 17 mg/kg sodium pentobarbital (Wang et al., 2018, 2019). A longitudinal incision was made in the dorsal midline skin, and the spinal cord was exposed by a laminectomy at vertebral segments T9–10. For the experiment group, a hemisection was performed at vertebral segment T9 with an ophthalmic iris knife (Beaver-Visitec International, Waltham, MA, USA). For the sham group, the hemisection was not performed. All rats were housed in temperature- and humidity-controlled home cages with free access to water and food after surgery.

Immunostaining observation

Rats were anesthetized with mixed narcotics at 1, 4, 7, 14, and 21 days after spinal cord injury and were transcardially perfused with 4% paraformaldehyde. Spinal cord segments at the caudal region of the injured site were harvested, dehydrated in 30% sucrose, embedded in Tissue-Tek O.C.T. compound (Sakura Finetek USA, Inc., Torrance, CA, USA), and coronally cryo-sectioned with a cryostat microtome (Leica, Wetzlar, Germany). The 20-μm spinal cord slices were mounted on microscope slides and incubated in primary antibody overnight at 4°C. The primary antibodies were rabbit anti-glial fibrillary acidic protein (GFAP; 1:400, Cat# Z0334, RRID:AB_10013382; Agilent Technologies, Santa Clara, CA, USA) or rabbit anti-ionized calcium binding adapter molecule 1 (IBA1; 1:200, Cat# 019-19741, RRID:AB_839504; FUJIFILM Wako Shibayagi Corporation, Osaka, Japan). This was followed by incubation with the secondary antibody Alexa Fluor 488-conjugated anti-rabbit IgG (1:400, Cat# ab150077; Abcam, Cambridge, UK) for 2 hours at room temperature. Spinal cord slices were photographed under fluorescence microscopy (Axio Imager M2, Carl Zeiss Meditec AG, Jena, Germany).

Quantitative real-time RT-PCR

At 1, 4, 7, and 14 days after surgery, 5-mm-long rat spinal cord segments at the caudal region were collected after laminectomy or hemisection. Total RNA samples were prepared from spinal cord segments with TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, USA). RNA integrity was confirmed by agarose gel electrophoresis. RNA samples (1 μg) were reverse transcribed to complementary DNA using HiScript III RT SuperMix for qPCR (+gDNA wiper) (Vazyme Biotech Co., Ltd., Nanjing, China). Sequence-specific primers were designed according to previous publications (Peinnequin et al., 2004; Bonefeld et al., 2008; Martínez-Beamonte et al., 2011; Bangaru et al., 2012; Gambarotta et al., 2014) or with Primer Express software (v3.0.1; Thermo Fisher Scientific) and synthesized by Sangon Biotech (Shanghai) Co., Ltd. (Shanghai, China). The primer sequences of the candidate reference genes tested are in . Primer sequences of NONRATT020851.2 were (sense) 5′-TTG GCA ATG GCA GTG TAT GTC-3′ and (antisense) 5′-TCA AGA GCA ACG GAA CAA GAA C-3′, and primer sequences of NONRATT007410.2 were (sense) 5′-GCA CGG TCA GAG GTC AGT T-3′ and (antisense) 5′-GAT GGC GGC TCA GTG GTA A-3′. Quantitative RT-PCR was performed with ChamQ Universal SYBR qPCR Master Mix (Vazyme) and specific primers in a StepOne Real-Time PCR System (Thermo Fisher Scientific). The thermocycler program was one initial denaturation cycle at 95°C for 5 minutes followed by 40 cycles of 95°C for 15 seconds and 60°C for 30 seconds. Each sample was tested in triplicate and analyzed using the 2–??CT method (Livak and Schmittgen, 2001) to determine the relative gene expression.

Bioinformatic analysis

Raw Ct values of candidate internal control genes were analyzed using geNorm and NormFinder bioinformatic software. geNorm analysis (v3.4; https://genorm.cmgg.be/) was applied using the Microsoft Excel 2007 Visual Basic plug-in (Microsoft Corporation, Redmond, WA, USA) to calculate the gene expression stability value (M-value), which is negatively correlated with gene expression stability. Additionally, we used geNorm software to calculate the pairwise variation (Vn/n+1), a normalization factor. Vn/n+1 less than 0.15 indicates that no further (n + 1) internal control genes are needed (Vandesompele et al., 2002). NormFinder analysis (v19; https://moma.dk/normfinder-software) was applied using the Microsoft Excel 2007 Visual Basic plug-in (Microsoft Corporation) to calculate intra-group and inter-group variations and estimate expression stability values (Andersen et al., 2004).

Statistical analysis

Statistical analyses were performed using GraphPad Prism 6.0 (GraphPad Software, San Diego, CA, USA). Differences between the experiment and control groups were analyzed by one-way analysis of variance followed by the Holm-Šídák multiple comparisons test. Significant differences were set at P < 0.05.

Results

Spinal cord hemisection elicits significant cellular changes

Rat spinal cord segments at the caudal region of the injured site were collected at 1, 4, 7, 14, and 21 days after spinal cord hemisection at vertebral segment T9 and immunostained to observe injury-induced morphological changes. The expressions of two main glial cell types, astrocytes and microglia, were visualized by immunofluorescence assay, as glial cells make up a large proportion of cells in mammals and play essential roles after spinal cord injury (Bahney and von Bartheld, 2018).

Few astrocytes were present at 1 day after spinal cord injury by immunostaining with the astrocyte marker GFAP (Zhang et al., 2017). GFAP immunopositivity was elevated at later time points, reaching its peak at 7 days after injury and gradually reduced toward the basal value at 14 and 21 days after spinal cord injury (). Immunostaining from the astrocyte marker GFAP was mainly localized in the gray matter, whereas immunostaining from the microglia marker IBA1 (Waller et al., 2019) was observed in the gray and white matter. The immunopositivity of IBA1 was elevated at 4 days after spinal cord injury and recovered to the basal level at 21 days after injury (), similar to GFAP results. These observations showed that there are robust morphological and molecular changes after spinal cord hemisection. These remarkable cellular changes indicate that many genes, including reference genes, might be differentially expressed after spinal cord injury.

Immunostaining of rat spinal cord segments at the caudal region at 1, 4, 7, 14, and 21 days following spinal cord injury.

(A) GFAP immunofluorescence images (green, Alexa Fluor 488) and immunopositivity of GFAP in injured rat spinal cord. The immunopositivity of GFAP gradually elevated and reached a peak value at 7 days after injury. (B) IBA1 immunofluorescence images (green, Alexa Fluor 488) and immunopositivity of IBA1 in injured rat spinal cord. The immunopositivity of IBA1 increased at 4 days after spinal cord injury and recovered to the basal level at 21 days after injury. Scale bars: 500 μm. GFAP: Glial fibrillary acidic protein; IBA1: ionized calcium binding adapter molecule 1.

Changes in candidate reference gene expression in spinal cord

To select potential internal controls for the normalization of target genes, the commonly used reference genes Actb, Ankrd27, CypA, Gapdh, Hprt1, Mrpl10, Pgk1, Rictor, Rn18s, Tbp, Ubc, Ubxn11, and Ywhaz were analyzed by RT-PCR. All these reference genes had linear correlation coefficients (R) close to 1 and amplification efficiencies between 90% to 110%, which suggests high efficiency and quality of the synthesized primers of these reference genes ().

Next, the raw Ct values of these reference genes were determined (). The raw Ct values of Rn18s were less than 7, which are the lowest values among the tested reference genes and are consistent with the high abundance of the Rn18s ribosomal RNA gene. Ubxn11 had the highest raw Ct values (~28), and other reference genes had raw Ct values ranging from 16 to 25. Notably, some reference genes, such as Rn18s and Actb, showed curved trends in their raw Ct values, especially in samples from spinal cord after hemisection. Conversely, the raw Ct values of some reference genes, such as Ankrd27, Gapdh, Mrpl10, Tbp, Ubc, and Ywhaz, did not exhibit obvious changes.

Raw Ct values of candidate reference genes at 1, 4, 7, and 14 days following spinal cord injury.

Raw Ct values of Actb(A), Ankrd27(B), CypA(C), Gapdh(D), Hprt1(E), Mrpl10(F), Pgk1(G), Rictor(H), Rn18s(I), Tbp(J), Ubc(K), Ubxn11(L), and Ywhaz(M) after laminectomy (Sha) or hemisection (Exp) surgery of rat spinal cord. Data are expressed as mean ± SD (n = 3 rats/group) and were analyzed by one-way analysis of variance followed by the Holm-Šídák multiple comparisons test. Actb: Beta-actin; Ankrd27: Ankyrin repeat domain 27; CypA: cyclophilin A; Gapdh: glyceraldehydes-3-phosphate dehydrogenase; Hprt1: hypoxanthine phosphoribosyltransferase 1; Mrpl10: mitochondrial ribosomal protein L10; Pgk1: phosphoglycerate kinase 1; Rictor: RPTOR independent companion of mTOR, complex 2; Rn18s: 18S ribosomal RNA; Tbp: TATA box binding protein; Ubc: ubiquitin C; Ubxn11: UBX domain protein 11; Ywhaz: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta.

Determination of expression stability of candidate reference genes

The raw Ct values of the 13 reference genes were analyzed using geNorm to determine gene expression stability, and the M-value, which negatively correlates with gene stability, was calculated. Actb had the highest M-value, which suggests that the stability of the expression of Actb was relatively low. Conversely, Ankrd27 and Ubc had the lowest M-values and the highest gene expression stabilities (). Calculation of pairwise variation (Vn/n+1), a normalization factor that advises the ideal number of internal control genes, showed that the pairwise variation of V2/3 was low (). Bioinformatic analysis using NormFinder indicated that Tbp, Ankrd27, and Ubc were the most stable genes, while Rn18s and Actb were the least stable genes in spinal cord segments following rat spinal cord injury (). Therefore, the reference genes Tbp, Ankrd27, and Ubc can be used as internal controls for the determination of gene expressions in injured spinal cord segments.

geNorm-determined expression stability of candidate reference genes.

(A) Average expression stability of the candidate reference genes Actb, Ankrd27, CypA, Gapdh, Hprt1, Mrpl10, Pgk1, Rictor, Rn18s, Tbp, Ubc, Ubxn11, and Ywhaz. (B) Pairwise variations of reference genes. Actb: Beta-actin; Ankrd27: Ankyrin repeat domain 27; CypA: cyclophilin A; Gapdh: glyceraldehydes-3-phosphate dehydrogenase; Hprt1: hypoxanthine phosphoribosyltransferase 1; Mrpl10: mitochondrial ribosomal protein L10; Pgk1: phosphoglycerate kinase 1; Rictor: RPTOR independent companion of mTOR, complex 2; Rn18s: 18S ribosomal RNA; Tbp: TATA box binding protein; Ubc: ubiquitin C; Ubxn11: UBX domain protein 11; Ywhaz: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta. Vn/n+1below 0.15 indicates that no further (n + 1) internal control genes are needed. For example, when V2/3= 0.044, Ankrd27 and Ubc are sufficient for gene normalization, and Tbp is not needed.

NormFinder-determined gene expression stability of candidate reference genes.

Stability values of the candidate reference genes Actb, Ankrd27, CypA, Gapdh, Hprt1, Mrpl10, Pgk1, Rictor, Rn18s, Tbp, Ubc, Ubxn11, and Ywhaz. Actb: Beta-actin; Ankrd27: Ankyrin repeat domain 27; CypA: cyclophilin A; Gapdh: glyceraldehydes-3-phosphate dehydrogenase; Hprt1: hypoxanthine phosphoribosyltransferase 1; Mrpl10: mitochondrial ribosomal protein L10; Pgk1: phosphoglycerate kinase 1; Rictor: RPTOR independent companion of mTOR, complex 2; Rn18s: 18S ribosomal RNA; Tbp: TATA box binding protein; Ubc: ubiquitin C; Ubxn11: UBX domain protein 11; Ywhaz: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta.

Application of stable or unstable candidate reference genes as internal control genes

The candidate reference genes were used to determine expression patterns of the genes NONRATT020851.2 and NONRATT007410.2, which were randomly selected from previously obtained sequencing data of rat spinal cord segments at 1, 4, and 7 days after injury (stored in the National Omics Data Encyclopedia [NODE], https://www.biosino.org/node/, accession number OEP000369) (Yu et al., 2019). Sequencing data showed that the expression levels of NONRATT020851.2 were elevated at 1, 4, and 7 days after spinal cord injury compared with the control (P < 0.05), reaching a peak at 4 days after injury (). The expression patterns of NONRATT020851.2 were different when using Actb, a reference gene with low stability, as the internal control for normalization. The relative abundance of NONRATT020851.2 decreased at 1 day after injury and increased to peak levels at 7 days after injury when using the unstable reference gene Actb to calculate ΔΔCt (). Conversely, the use of the stable reference genes Tbp or Ankrd27+Ubc as internal controls achieved similar trends of NONRATT020851.2 expression as sequencing data (Figure and ). Similarly, the expression patterns of NONRATT007410.2 were more comparable to sequencing observations when using Tbp or Ankrd27+Ubc as internal controls instead of Actb (Figure ).

Relative expression of target genes .

(A) Gene expression patterns of NONRATT020851.2 from sequencing data. (B–D) Gene expression of NONRATT020851.2 was determined by qRT-PCR using Actb(B), Tbp(C), and Ankrd27+Ubc(D) as internal controls. (E) Gene expression patterns of NONRATT007410.2 from sequencing data. (F–H) Gene expression of NONRATT007410.2 was determined by qRT-PCR using Actb(F), Tbp(G), and Ankrd27+Ubc(H) as internal controls. Actb is an unstable internal control, and Tbp and Ankrd27+Ubc are stable internal controls. Data are expressed as mean ± SD (n = 3 rats/group). *P < 0.05, vs. control group (Con) using one-way analysis of variance followed by the Holm-Šídák multiple comparisons test. Actb: Beta-Actin; Ankrd27: Ankyrin repeat domain 27; Tbp: TATA box binding protein; Ubc: ubiquitin C.

Discussion

qPCR is the gold standard for the determination of gene abundance and the validation of large-scale transcript data, such as from microarrays or sequencing (Derveaux et al., 2010). qPCR has advantages, such as high sensitivity and simplicity, and disadvantages, including stringent sample quality requirements, operator variability, and the inefficiency of reverse transcription (Bustin and Nolan, 2004). The use of internal controls largely reduces variations caused by technological defects or operational errors and increases the reliability and accuracy of RT-PCR results (Kadegowda et al., 2009).

Reference genes are essential genes with relatively constant abundance in all cells for various biological activities (Eisenberg and Levanon, 2013). Therefore, reference genes are commonly used as internal controls for calibration of RT-PCR. However, emerging studies have demonstrated that reference gene expression is tissue- and condition-dependent. For instance, Actb and Gapdh, two widely used reference genes, had dissimilar histological expression profiles in the liver, the vascular system, and the digestive tract (Lin and Redies, 2012). Our previous sequencing analyses also revealed that the expression level of Actb was upregulated in rat sciatic nerve segments after peripheral nerve injury, which suggests that the expression of reference genes may be affected by pathophysiological changes (Wang et al., 2018). Recently, Košuth et al. (2019) analyzed the expression stability of several reference genes in the rat spinal cord using two different injury models, minimal and contusion spinal cord injury. They found Gapdh and Hprt1 were suitable for contusion spinal cord injury experiments, while Gapdh was suitable for minimal spinal cord injury experiments. The expression stability of reference genes depended on the severity of the injury. Therefore, it is of great importance to select reference genes as internal controls for normalization that are stably expressed in the target cells, tissues, and biological conditions of interest.

Previously, we examined the stability of many frequently used reference genes in different tissues after rat peripheral nerve injury. We discovered that Mrpl10, Tbp, and Hprt1 are suitable internal controls for sciatic nerve segments, dorsal root ganglion samples, and gastrocnemius muscles, respectively (Wang et al., 2017, 2019). In the current study, we focused on spinal cord injury—injury to the central nervous system that induces even more serious consequences—and measured the expression of reference genes in spinal cord segments after central nervous system injury. For all tested reference genes, geNorm-calculated gene expression stability M-values were less than the default limit of 1.5 and pairwise variation V-values were less than the default limit of 0.15 (Bangaru et al., 2012). Therefore, the tested reference genes met the basic requirements for use as internal controls. Raw Ct values and geNorm and NormFinder outcomes showed that Tbp, Ankrd27, and Ubc were the most stably expressed reference genes in spinal cords injured by hemisection. Moreover, a comparison study of the gene expression of NONRATT020851.2 and NONRATT007410.2 using the stable reference genes Tbp, Ankrd27, and Ubc and an unstable reference gene Actb further validated and supported the expression stability of Tbp, Ankrd27, and Ubc. The outcome of our current study is suitable for the spinal cord hemisection model, but we still need to find suitable internal control genes for use in other injury models, such as contusion or full transection. We used female rats in our study, as their use is preferred in spinal cord injury research (Onifer et al., 2007; Stewart et al., 2020). Many sex-related genes contribute to morphological, physiological, and behavioral differences between males and females (Shi et al., 2019), so it is worth exploring in the future whether sex differences affect the selection of internal control genes.

In sum, the identification of suitable internal controls provides guidance for RT-PCR experiments and facilitates the exposure and understanding of the underlying genetic infrastructure of spinal cord injury and regeneration.

Table 1

Primer sequences of candidate reference genes

Symbol	Accession number	Primer sequences (5’–3’)	Amplicon size (bp)	Reference
Actb	NM_031144.2	Sense: TGT CAC CAA CTG GGA CGA TA	165	Bonefeld et al., 2008
		Antisense: GGG GTG TTG AAG GTC TCA AA
Ankrd27	NM_001271264.1	Sense: CCC AGG ATC CGA GAG GTG CTG TC	95	Gambarotta et al., 2014
		Antisense: CAG AGC CAT ATG GAC TTC AGG GGG
CypA	NM_017101.1	Sense: CCA AAC ACA AAT GGT TCC CAG T	135	Bangaru et al., 2012
		Antisense: ATT CCT GGA CCC AAA ACG CT
Gapdh	NM_017008.4	Sense: CCA CCA ACT GCT TAG CCC CC	91	Gambarotta et al., 2014
		Antisense: GCA GTG ATG GCA TGG ACT GTG G
Hprt1	NM_012583.2	Sense: CTC ATG GAC TGA TTA TGG ACA GGA C	123	Peinnequin et al., 2004
		Antisense: GCA GGT CAG CAA AGA ACT TAT AGC C
Mrpl10	NM_001109620.1	Sense: CTC CTC CCA AGC CCC CCA AG	97	Gambarotta et al., 2014
		Antisense: CAG ACA GCT ATC ATT CGG TTG TCC C
Pgk1	NM_053291.3	Sense: GCA GAT TGT TTG GAA CGG TCC	113	Martínez-Beamonte et al., 2011
		Antisense: TAG TGA TGC AGC CCC TAG ACG T
Rictor	XM_001055633.3	Sense: TCC GAA TAC GAG GGC GGA A	142	*
		Antisense: AGA TGG CCC AGC TTT CTC ATA
Rn18s	X01117.1	Sense: ACT CAA CAC GGG AAA CCT CA	114	Bangaru et al., 2012
		Antisense: AAT CGC TCC ACC AAC TAA GA
Tbp	NM_001004198.1	Sense: TAA TCC CAA GCG GTT TGC TG	111	Martínez-Beamonte et al., 2011
		Antisense: TTC TTC ACT CTT GGC TCC TGT G
Ubc	NM_017314.1	Sense: TCG TAC CTT TCT CAC CAC AGT ATC TAG	82	Gambarotta et al., 2014
		Antisense: GAA AAC TAA GAC ACC TCC CCA TCA
Ubxn11	NM_138853.2	Sense: GCG AGA CTG GAT GAA GGC CAA G	120	Gambarotta et al., 2014
		Antisense: CCC TCC ACC ACC AGC TCA CTC
Ywhaz	NM_013011.3	Sense: TTG AGC AGA AGA CGG AAG GT	136	Bonefeld et al., 2008
		Antisense: GAA GCA TTG GGG ATC AAG AA

* Designed by Primer Express® software. Actb: Beta-actin; Ankrd27: Ankyrin repeat domain 27; CypA: cyclophilin A; Gapdh: glyceraldehydes-3-phosphate dehydrogenase; Hprt1: hypoxanthine phosphoribosyltransferase 1; Mrpl10: mitochondrial ribosomal protein L10; Pgk1: phosphoglycerate kinase 1; Rictor: RPTOR independent companion of mTOR, complex 2; Rn18s: 18S ribosomal RNA; Tbp: TATA box binding protein; Ubc: ubiquitin C; Ubxn11: UBX domain protein 11; Ywhaz: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta.

Table 2

The amplification efficiencies of candidate reference genes

Symbol	Slope	R 2	Efficiency (%)
Actb	–3.192	0.9714	106
Ankrd27	–3.174	0.9996	107
CypA	–3.543	0.9985	92
Gapdh	–3.177	0.9871	106
Hprt1	–3.341	0.986	99
Mrpl10	–3.316	0.9971	100
Pgk1	–3.146	0.9915	108
Rictor	–3.483	0.9961	94
Rn18s	–3.257	0.9655	103
Tbp	–3.105	0.9585	110
Ubc	–3.276	0.9948	102
Ubxn11	–3.526	0.9867	92
Ywhaz	–3.378	0.9999	98

Actb: Beta-actin; Ankrd27: Ankyrin repeat domain 27; CypA: cyclophilin A; Gapdh: glyceraldehydes-3-phosphate dehydrogenase; Hprt1: hypoxanthine phosphoribosyltransferase 1; Mrpl10: mitochondrial ribosomal protein L10; Pgk1: phosphoglycerate kinase 1; Rictor: RPTOR independent companion of mTOR, complex 2; Rn18s: 18S ribosomal RNA; Tbp: TATA box binding protein; Ubc: ubiquitin C; Ubxn11: UBX domain protein 11; Ywhaz: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta.