Validation of housekeeping genes for quantitative real-time PCR in in-vivo and in-vitro models of cerebral ischaemia

  • Carme Gubern1,

    Affiliated with

    • Olivia Hurtado2,

      Affiliated with

      • Rocío Rodríguez1,

        Affiliated with

        • Jesús R Morales2,

          Affiliated with

          • Víctor G Romera1,

            Affiliated with

            • María A Moro2,

              Affiliated with

              • Ignacio Lizasoain2,

                Affiliated with

                • Joaquín Serena1 and

                  Affiliated with

                  • Judith Mallolas1Email author

                    Affiliated with

                    BMC Molecular Biology200910:57

                    DOI: 10.1186/1471-2199-10-57

                    Received: 03 November 2008

                    Accepted: 16 June 2009

                    Published: 16 June 2009



                    Studies of gene expression in experimental cerebral ischaemia models can contribute to understanding the pathophysiology of brain ischaemia and to identifying prognostic markers and potential therapeutic targets. The normalization of relative qRT-PCR data using a suitable reference gene is a crucial prerequisite for obtaining reliable conclusions. No validated housekeeping genes have been reported for the relative quantification of the mRNA expression profile activated in in-vitro ischaemic conditions, whereas for the in-vivo model different reference genes have been used.

                    The present study aims to determine the expression stability of ten housekeeping genes (Gapdh, β2m, Hprt, Ppia, Rpl13a, Oaz1, 18S rRNA, Gusb, Ywhaz and Sdha) to establish their suitability as control genes for in-vitro and in-vivo cerebral ischaemia models.


                    The expression stability of the candidate reference genes was evaluated using the 2-ΔC'T method and ANOVA followed by Dunnett's test. For the in-vitro model using primary cultures of rat astrocytes, all genes analysed except for Rpl13a and Sdha were found to have significantly different levels of mRNA expression. These different levels were also found in the case of the in-vivo model of pMCAO in rats except for Hprt, Sdha and Ywhaz mRNA, where the expression did not vary. Sdha and Ywhaz were identified by geNorm and NormFinder as the two most stable genes.


                    We have validated endogenous control genes for qRT-PCR analysis of gene expression in in-vitro and in-vivo cerebral ischaemia models. For normalization purposes, Rpl13a and Sdha are found to be the most suitable genes for the in-vitro model and Sdha and Ywhaz for the in-vivo model. Genes previously used as housekeeping genes for the in-vivo model in the literature were not validated as good control genes in the present study, showing the need for careful evaluation for each new experimental setup.


                    Oxygen-glucose deprivation (OGD) and middle cerebral artery occlusion (MCAO) are accepted as in-vitro and in-vivo cerebral ischaemia models [14]. The permanent MCAO (pMCAO) focal ischaemia model is widely used in stroke research and is helpful in elucidating physiopathological causes and identifying prognostic markers and potential therapeutic targets [59].

                    Real-time polymerase chain reaction (RT-PCR) is a powerful technique for gene expression studies due to its high sensitivity, specificity and broad quantification range for high throughput and accurate expression profiling of selected genes [10]. The data obtained by quantitative RT-PCR (qRT-PCR) is typically normalised with an internal control, often referred to as a housekeeping gene. Up to now, several mathematical methods, such as NormFinder [11], geNorm [12] and 2-ΔC'T method [13], have been developed to analyse the variability of the expression of candidate housekeeping genes. The ideal housekeeping gene for qRT-PCR would be one whose mRNA is consistently expressed at the same level in all samples under investigation, regardless of tissue type, disease state, medication or experimental conditions, and would have expression levels comparable to that of the target [14, 15]. However, no ideal reference gene has yet been discovered and it is reasonable to suppose that such a gene may not exist. Although different genes have been used, several previous studies have shown that the suitability of a particular reference gene depends on the insult, lesion and experimental conditions applied to the cell culture or experimental animal [16, 17]. Therefore, careful evaluation and validation of control genes is required prior to using them in each individual case to avoid possible inaccuracies stemming from the use of an unsuitable reference gene [14, 16, 1820].

                    In addition, no validated housekeeping genes have been reported for the relative quantification of the mRNA expression profile activated in in-vitro ischaemic conditions. With regards to the in-vivo model (pMCAO), Glyceraldehyde-3-phosphate dehydrogenase (Gapdh) [16, 21, 22], Peptidyl-prolyl cis-trans isomerasa (Ppia) [2123], Ribosomal Protein L13A (Rpl13a)[23] and Hypoxanthine phosphoribosyltransferase (Hprt) [16] have previously been accepted as control genes. The aim of this study was to confirm the validity of these reference genes in our pMCAO samples and to test their usefulness as endogenous controls in the in-vitro model using primary cultures of rat astrocytes. Furthermore, we aimed to analyse the expression stability of other reference genes that have been used, Beta-2 microglobulin (β2m) and 18S rRNA [24], Glucoronidase B (Gusb) [25], Succinate dehydrogenase complex, subunit A (Sdha) [14, 25, 26], Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide (Ywhaz) [14, 26] and a new candidate gene for the in-vivo model, Ornithine decarboxylase antizyme 1 (Oaz1), which has been hypothesised as being suitable in experimental conditions involving high protein turn-over [10].


                    Evaluation and validation of selected candidate housekeeping genes in an in-vitro ischaemia model

                    The results of gene expression of Gapdh, Hprt, Ppia, 18S rRNA, β2m, Oaz1, Gusb, Ywhaz, Sdha and Rpl13a from astrocytes exposed to OGD during 3 hours and collected at 0 minutes, 30 minutes, 2, 6, 24 and 48 hours after reperfusion are shown in Figure 1.
                    Figure 1

                    Effects of OGD in astrocyte expression of housekeeping genes: (A) Gapdh, (B) Hprt, (C) Ppia, (D) 18S rRNA, (E) β2m, (F) Oaz1, (G) Gusb, (H) Ywhaz, (I) Sdha and (J) Rpl13a. Fold change in gene expression analysed by the 2-ΔC'T (see methods for details). Data are mean ± S.E.M., n = 3: * p < 0.05, ** p < 0.01, # p < 0.0001 vs. control.

                    The expression of Gapdh was found to have significant fluctuations in the levels of mRNA (p < 0.0001) increasing significantly at 0 h after reperfusion (p < 0.05) whereas no significant differences were observed at other points (Figure 1A). The expression of Hprt, in comparison to the control group, was found to fluctuate significantly in all the groups as a whole (p < 0.0001). This fluctuation was a significant increase at 0 min (p < 0.0001), 30 min (p < 0.01) and 6 h after reperfusion (p < 0.0001) (Figure 1B). Significant differences in the expression of Ppia were found between time groups (p < 0.0001). Levels increased at 0 and 30 min (p < 0.0001 in both cases) and at 2 h and 6 h (p < 0.05 in both cases) after reperfusion in comparison with the control group (Figure 1C). The 18S rRNA mRNA showed significant level fluctuations (p < 0.0001). Levels decreased at 24 h (p < 0.01) and 48 h (p < 0.0001) after reperfusion (Figure 1D). Changes in β2m mRNA levels were found to be significant between time groups (p < 0.0001). β2m expression decreased at 0 min (p < 0.01), 30 min and 6 h (p < 0.0001 in both cases) and 48 h (p < 0.05) in comparison with the control group (Figure 1E). Oaz1 also showed significant differences in expression between groups (p < 0.0001). The expression was stable at the first four analysis times but decreased significantly at 24 h (p < 0.01) and 48 h (p < 0.05) (Figure 1F). Significant fluctuations in the mRNA expression of Gusb were found (p < 0.0001) decreasing at 2 h after reperfusion (p < 0.0001) (Figure 1G). The expression of Ywhaz was found to have significant fluctuations in mRNA levels (p < 0.0001) (Figure 1H). Analysis of the same samples with a specific gene expression assay for Sdha and Rpl13a did not reveal any significant differences in expressions at the different time points when analysis was made (p = 0.065 and p = 0.244, respectively) (Figure 1I and 1J).

                    Evaluation and validation of selected candidate housekeeping genes in an in-vivo cerebral ischaemia model

                    The gene expression of candidate housekeeping genes, β2m, Rpl13a, Oaz1, Gapdh, Ppia, Gusb, 18SrRNA, Sdha, Ywhaz and Hprt, in cerebral cortex samples collected at 30 minutes and 2, 6, 24, 48 and 72 hours after pMCAO is shown in Figure 2.
                    Figure 2

                    Effects of pMCAO in expression of housekeeping genes: (A) β2m, (B) Rpl13a, (C) Oaz1, (D) Ppia, (E) Gusb, (F) 18SrRNA, (G) Gapdh, (H) Sdha, (I) Hprt and (J) Ywhaz. Fold change in gene expression analysed by the 2-ΔC'T (see methods for details). Data are mean ± S.E.M., n = 6: * p < 0.05, ** p < 0.01, # p < 0.0001 vs. control.

                    β2m expression presented a significant increase at 48 and 72 h (p < 0.0001 in both cases) and was significantly different between groups (p < 0.0001) (Figure 2A). The expression of Rpl13a increased 72 h after pMCAO (p < 0.05) with a significant fluctuation between groups (p < 0.01) (Figure 2B). The level of mRNA variation of Oaz1 was found to be significant for all the time groups taken as a whole (p < 0.01) (Figure 2C). Significant differences were found between groups for the expression of Ppia (p < 0.0001) with mRNA decreasing significantly at 2 and 6 h (p < 0.01 in both cases) after pMCAO (Figure 2D). The expression of Gusb was found to fluctuate significantly between the different test points (p < 0.0001). At 48 h and 72 h, levels of Gusb mRNA had increased (p < 0.01 and p < 0.0001 respectively) (Figure 2E). The level of the mRNA of 18S rRNA showed significant fluctuation (p < 0.05) decreasing significantly at 48 h (p < 0.05) and 72 h (p < 0.01) (Figure 2F). Changes in Gapdh mRNA levels were found to be significant between time groups (p < 0.0001). Gapdh expression decreased at 30 min (p < 0.01), 2 h (p < 0.05) and 6 h (p < 0.01) in comparison with the control group (Figure 2G). mRNA expression was invariable under in-vivo ischaemic brain conditions for the Sdha (p = 0.067), Hprt (p = 0.236) and Ywhaz (p = 0.244) genes (Figure 2H–2J).

                    The stability of these three genes was evaluated using geNorm and NormFinder algorithms. GeNorm identified Ywhaz and Sdha as the most stable pairwise combination of reference genes for the in-vivo ischaemia model (Figure 3A). Figure 3B shows the pairwise variation value (V) between two sequential normalization factors containing an increasing number of genes. As we only analysed three genes, the figure shows the pairwise variation value of V2/3 (0.000). According to Vandesompele et al. [12], the ideal pairwise variation value is less than 0.15 and therefore two genes, in our case Ywhaz and Sdha, are enough for normalization. NormFinder identified Ywhaz as the most stable gene with a stability value of 0.023 and, as with geNorm, the best two genes for normalization were Ywhaz and Sdha with a stability value of 0.022 (Table 1).
                    Figure 3

                    Gene expression stability and determination of the optimum number of genes for normalization of the candidate reference genes for the in-vivo ischaemia model using geNorm analysis. (A) Expression stability plot showing average expression stability values M. (B) Pairwise variation analysis to determine the optimal number of reference genes for use in RT qPCR data normalization. The use of the two most stable genes is sufficient for accurate normalization (cutoff 0.15 according to [12])

                    Table 1

                    Stability values of housekeeping genes obtained by NormFinder


                    Gene name

                    Stability value










                    Best combination of two genes

                    Sdha and Ywhaz



                    The selection of the housekeeping gene is crucial to the normalization of quantitative gene expression results. This is even more important in those settings in which gene expression is severely affected, as after cerebral ischaemia.

                    No validated housekeeping genes have been reported for the relative quantification of the mRNA expression profile activated in in-vitro ischaemic conditions. Our present study discards the use of Gapdh, Ppia, β2m, 18S rRNA, Hprt, Oaz1, Gusb and Ywhaz as reference genes and finds Rpl13a and Sdha to be the best control genes for the in-vitro model in the conditions we have tested.

                    With regards to the pMCAO model, four studies have reported the suitability of several reference genes for this ischaemia model. In Table 2 we summarise the methodologies employed and results obtained both for previously published and the present studies. Earlier investigations only tested from two to four candidate control genes finding invariable expressions for Rpl13a, Hprt, Ppia and Gapdh genes. In the case of Gapdh, results were contradictory: three studies validated the gene [16, 21, 22] whereas another rejected it [23]. Our results show that Gapdh, Rpl13a and Ppia together with other candidate genes, namely β2m, Gusb, 18S rRNA and Oaz1, are not appropriate as housekeeping genes for the pMCAO model. As previously reported by Meldgaard et al. when using the pMCAO model in mice [16], we have found Hprt to be a suitable housekeeping gene for the in-vivo ischaemia model. Moreover we have identified two new housekeeping genes for this model, Sdha and Ywhaz, which have not been previously analysed. According to the results obtained with the NormFinder and geNorm algorithms, these two genes are more stable and so the best suited for normalization of relative quantification gene expression studies. Several factors have an influence on the variability of control genes and may lead to differences between studies. In the same way that gene control expression may vary depending on whether the model is permanent (pMCAO) or transient (tMCAO) [27], different methodologies employed to perform the pMCAO, electrocautery or intraluminal filament, may have an influence on the gene expression. Sample collection and the control sample used are other factors which must be taken into consideration. Three different cerebral regions have been collected in each previous study: the whole of the infarcted cortex [21, 22], the whole of the infarcted hemisphere [16] and the ischaemic penumbra cortex [23]. In our study we have collected the infarcted and peri-infarcted cortical area to analyse changes in gene expression in the affected region. With regards to the control sample there are also different possibilities: brain samples from sham or naive rats, or healthy contralateral hemispheres. We consider sham rats to be the best control because this control minimises bias due to surgery. In our opinion, the use of the healthy contralateral hemisphere may lead to erroneuos conclusions because studies have reported structural changes as well as changes in gene expression in the contralateral hemisphere after ischaemia [2830]. Finally, the time course studied is also a relevant factor. We have validated Sdha, Ywhaz and Hprt as a suitable housekeeping genes for the pMCAO model conducting an analysis of the control situation and after ischaemia at frequent time intervals up to 72 hours, which vary from those used in earlier studies [16, 2123]. Considering all the variables that may influence gene expression, the validation of candidate reference genes in each individual experimental condition is crucial.
                    Table 2

                    Housekeeping gene validation in the pMCAO model


                    Harrison DC 2000

                    Medhurst AD 2000

                    Meldgaard M 2006

                    Tian YF 2007

                    Gubern C

                    Focal ischaemia model

                    pMCAO intraluminal filament

                    pMCAO intraluminal filament

                    pMCAO electrocautery

                    pMCAO intraluminal filament

                    pMCAO electrocautery


                    Male Sprague-Dawley rats

                    Male Sprague-Dawley rats

                    Male SJL mice

                    Male Sprague-Dawley rats

                    Male Fisher rats

                    Collected sample

                    Whole infarcted cortex

                    Whole infarcted cortex

                    Whole infarcted hemisphere

                    Ischaemic penumbra cortex

                    Infarcted and peri-infarcted cortex region

                    Time points

                    Naive, 3, 6, 12 and 24 h (pMCAO and sham)

                    Contralateral hemisphere, 6 and 24 h

                    Sham, 30 min, 1, 2, 4, 6, 12, 24 and 48 h, 5 and 10d

                    Sham, 2 h, 24 h, 3, 7 and 21d

                    Sham, 30 min, 2, 6, 24, 48 and 72 h

                    Real-Time PCR

                    Custom TaqMan assays

                    Custom TaqMan assays

                    Custom TaqMan assays


                    Inventoried TaqMan assays


                    Standard curve

                    Standard curve

                    %mRNA = 100 × (1+E)ΔCT



                    Statistical analyses


                    Student's unpaired t-test

                    Kruskal-Wallis + Dunn's

                    ANOVA + Tukey

                    ANOVA + Dunnett's

                    Suitable controls

                    Ppia & Gapdh

                    Ppia & Gapdh

                    Hprt & Gapdh

                    Rpl13a & Ppia

                    Hprt, Sdha & Ywhaz

                    Unsuitable controls



                    Gapdh, β-actin

                    Oaz1, Gapdh, Ppia, Rpl13a, B2m, Gusb, 18SrRNA

                    The geNorm [12] and NormFinder [11] algorithms are now widely used to determine the most stable reference genes from a set of candidates genes with invariable expression [26, 3134]. We have only used algorithms on genes identified as stable housekeeping genes in the in-vivo model in our analysis as the software requires a minimum of three genes for the analysis.

                    In order to check the value of the validation of endogenous controls, we have used different housekeeping genes analysed in this study to normalise the expression of cyclooxygenase 2 (Cox2) and Bcl2, genes which are known to be involved in cerebral ischaemia (data not shown). Cox2 is an inflammatory mediator gene which undergoes a significant increase due to ischaemia [35, 36]. Bcl2 is a protooncogene that codifies an antiapoptotic protein, which has been described as participating in neuroprotective mechanisms in ischaemic stroke although changes in the expression are limited to an increase in the neurons within the peri-infarct region [3739]. In the two cases we have observed differences in the results obtained when suitable and unsuitable control genes are used. These are even more critical when the differences in expression of the analysed gene are minor, as is the case with Bcl2, which has a different expression pattern depending on the housekeeping gene that is used.


                    In designing any study into relative gene expression quantification it is necessary to perform a validation of the housekeeping genes used which takes into account the specific experimental conditions to be encountered and the time course to be employed.

                    We report a detailed investigation of the effects of experimental ischaemia on the expression of potential housekeeping genes including genes typically used and new candidates as internal controls in qRT-PCR assays. This analysis was facilitated by the sensitivity and efficiency of real-time quantitative PCR. We determined that Gapdh, Rpl13a and Ppia, which have previously been validated as housekeeping genes in quantitative expression studies with the pMCAO model, had significant differences in their expression in our samples collected at different time points. Our study confirms the suitability of Hprt as a housekeeping gene for this in-vivo model and identifies two new housekeeping genes, Sdha and Ywhaz, which are more stable than Hprt. We have also identified, to the best of our knowledge for the first time, reference genes (Rpl13a and Sdha) for the OGD in-vitro model using primary astrocyte cultures.

                    The study also highlights the importance of validating control genes for each individual study as gene expression is influenced by many different factors which can lead to mistaken conclusions.


                    All procedures were approved by the Committee of Animal Care of the Universidad Complutense de Madrid, in accordance with the regulations of the European Union (86/609/CEE) and Spanish legislation (RD223/88).

                    Reference gene selection

                    Candidate reference genes were selected from those most commonly used in the literature [14, 2426], a new candidate [10], and those previously described as suitable for pMCAO studies [16, 2123] (see Table 3).
                    Table 3

                    Candidate reference genes

                    Gene Symbol

                    Gene Name

                    TaqMan Assay Number*

                    Genebank accession


                    Hypoxanthine phosphoribosyltransferase




                    Glyceraldehyde-3-phosphate dehydrogenase




                    Ornithine decarboxylase antizyme 1




                    Beta-2 microglobulin




                    Peptidyl-prolyl cis-trans isomerasa




                    Ribosomal Protein L13A




                    Eukaryotic 18S rRNA




                    Tyrosine 3-monooxygenase




                    Succinate dehydrogenase complex, subunit A




                    Glucuronidase, beta



                    *"_m" indicates an assay whose probe spans an exon junction and will not detect genomic DNA. "_s" indicates an assay whose probes and primers are designed within a single exon and, hence, will detect genomic DNA. "_g" indicates an assay that may detect genomic DNA as the assay probe and primers may also be within a single exon.

                    In-vitro ischaemia

                    Rat astrocyte cultures

                    Primary astrocyte cultures were prepared from neonatal (P0) Wistar rat cortex, as previously described [40, 41].

                    Exposure of astrocyte culture to oxygen-glucose deprivation (OGD)

                    Oxygen-glucose deprivation was performed as previously described [42]. Culture medium was replaced by a solution containing (in mM) NaCl (130), KCl (5.4), CaCl2 (1.8), NaHCO3 (26), MgCl2 (0.8), NaH2PO4 (1.18) and 2% FS bubbled with 95% N2/5% CO2 for OGD cells (OGD solution). OGD cells were transferred to an anaerobic chamber (Forma Scientific, Hucoa Erloss, Spain) containing a gas mixture of 95% N2/5% CO2 and humidified at 37°C, and maintained at a constant pressure of 0.15 bar. Time of exposure to OGD was 180 min. OGD was terminated by replacing the exposure medium with oxygenated MEM containing 0.6% glucose, 0.029% glutamine, 50 I.U./ml penicillin, 50 Ag/ml streptomycin, 10% FS (reperfusion medium) and returned to the normoxic incubator. Control cultures in a solution identical to OGD solution except that it contained glucose (33 mM; control solution) were kept in the normoxic incubator for the same time period as the OGD, and then incubation solution was replaced with the reperfusion buffer. Cultures were returned to the normoxic incubator until the end of the experiment.

                    Experimental groups

                    Astrocyte samples were collected from controls and from the OGD group at 0 minutes, 30 minutes, 2 hours, 6 hours, 24 hours and 48 hours after reperfusion (n = 5 in each group). Astrocytes were collected in RNA Protect Cell Reagent (Qiagen, Barcelona, Spain) and frozen at -80°C until RNA isolation was undertaken.

                    In-vivo cerebral ischaemia

                    Permanent middle cerebral artery occlusion (pMCAO) in rats

                    Experiments were performed on male Fischer rats weighing 250–300 g. Rats were anaesthetised with 2.5% halothane in a mixture of 70% nitrogen/30% oxygen. Permanent focal cerebral ischemia was induced by ligature of the left common carotid artery (CCA) and occlusion of the ipsilateral distal middle cerebral artery (MCA) as previously described [41, 42]. Briefly, for the CCA ligature, a midline ventral cervical incision was made, and the CCA was isolated and permanently occluded with a silk ligature. For the MCA occlusion, a 1-cm incision perpendicular to the line connecting the lateral canthus of the left eye and the external auditory canal was made to expose and retract the temporalis muscle. A 2-mm burr hole was drilled and the MCA was exposed by cutting and retracting the dura. The MCA was elevated and cauterised. Rats in which the MCA was exposed but not occluded served as sham-operated controls (SHAM). Following surgery, subjects were returned to their cages and allowed free access to water and food. The body temperature of animals was monitored throughout the experiment and was maintained at 37.5 ± 0.5°C using a heating pad.

                    Experimental groups

                    Brain were removed and cortical infarcted and peri-infarcted tissue was collected from controls and the pMCAO group at 30 minutes, 2 hours, 6 hours, 24 hours, 48 hours and 72 hours after artery occlusion (n = 6 in each group). Samples were quickly transferred to 1.5 mL tubes and frozen at -80°C until RNA isolation was undertaken.

                    RNA extraction, quantification and retrotranscription

                    Cerebral cortex frozen samples were immediately transferred to QIAzol lysis reagent (Qiagen) and then homogenised using an Ultra-Turrax T25 homogeniser. RNA was isolated using the RNeasy® Lipid Tissue Mini Kit (Qiagen) in accordance with the manufacturer's instructions and treated with DNase (RNase-free DNase Set, Qiagen) in order to remove any trace of genomic DNA. RNA was quantified spectrophotometrically (GeneQuant, Biochrom) at 260 nm (A260) and purity was estimated by an A260/A280 ratio > 1.8. The integrity, purity, and amount of RNA were verified by visualization of rRNAs after electrophoresis on agarose gel.

                    Total RNA was isolated from primary glial cultures using the RNeasy® Plus Mini Kit (Qiagen) and total RNA concentration was measured by fluorescence using a Quant-iT RNA Assay Kit (Invitrogen, Barcelona, Spain). Manufacturers' instructions were followed in both cases.

                    2 μg of RNA from each cortex samples and 1 μg of RNA from astrocyte samples were reverse transcribed using a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Barcelona, Spain) according to the manufacturer's instructions. Two independent RT reactions were performed for all samples.

                    Real-time PCR

                    RT-Q-PCR reactions were carried out for all genes of interest in each sample using specific TaqMan® Gene Expression Assays (Table 3) and an ABI PRISM® 7000 Instrument (Applied Biosystems). Each reaction was performed in a final volume of 25 μL containing 11.25 μL cDNA diluted with H2O, 1.25 μL TaqMan® Gene Expression Assay and 12.5 μL of TaqMan® Universal PCR Master Mix (Applied Biosystems). Amplifications were performed starting with a 2 min activation step for Amperase UNG at 50°C, 10 min template denaturation step at 95°C, followed by 40 cycles of 95°C for 15s and 60°C for 1 min. As we used TaqMan® Gene Expression Assays validated with PCR efficiency near 100% for all of them, the expression of each housekeeping gene was obtained using the 2-ΔC'T method, a modification of the 2-ΔΔCt method described by K. Livak et al., where ΔCt (CtTime X - CtControl) and Control is the calibrator sample representing the 1× expression of each gene [13, 24]. All samples were run in triplicate and average values were calculated. Two independent reverse transcriptions were tested for each gene.

                    As some of the Taqman Gene Expression assays employed are not cDNA specific, we confirmed no genomic DNA contamination of our samples using no-reverse transcription controls. A no-reverse transcription control is a reaction that has been prepared for reverse transcription with RNA and all the components of the retrotranscription kit except for reverse transcriptase.

                    Statistical analysis

                    Data are mean ± standard error of the mean (SEM). A statistical test was applied to look for significant differences between experimental conditions for each candidate control gene in each experimental model. The variance equality hypothesis was confirmed for each gene in the two experimental models as Bartlett's test failed to reveal any significant variance with a confidence level of 95% (data not shown). A one-way analysis of variance (ANOVA) test was therefore conducted, followed by Dunnett's post-hoc analysis to compare each time set against the control group mean. A p value < 0.05 was considered statistically significant. Calculations were performed using the 11.0.0 version of SPSS software (SPSS Inc).

                    Determination of reference gene expression stability

                    After determining which housekeeping genes had an invariable expression in the in-vivo model using 2-ΔC'T method and ANOVA + Dunnett's post hoc analysis, we used two publicly available software tools, geNorm [12] and NormFinder [11], to analyse their gene expression stability. GeNorm calculates a gene-stability measure M, which is the average pairwise variation of a particular gene compared with all other candidate reference genes. A lower value of M indicates the greater stability of the reference gene. NormFinder estimates the overall expression variation of the candidate normalization genes and the variation between sample subgroups of the sample set using a model-based approach.



                    This work was partially supported by the Spanish Ministry of Science and Innovation and the European Regional Development Fund (ERDF) PROFIT grant nos. CIT 090100-2005-15, CIT 090100-2007-43, CIT 090000-2008-11 (J. Mallolas and I. Lizasoain), Spanish Ministry of Science and Innovation, Institut de Salud Carlos III (ISCIII) grant nos. FIS-PI070322 (J. Mallolas), ISCIII retics-RENEVAS-RD06/0026/0005 (MA. Moro) and RD07/0026/2002 (J. Serena) and the Centre d'Innovació i Desenvolupament Empresarial (CIDEM) of the Generalitat de Catalunya (grant no.RDITSCON07-1-0006). C. Gubern was recipient of a BE grant no.2007 BE-100293 from the Agència de Gestió d'Ajuts Universitaris i de Recerca (AGAUR).

                    Authors’ Affiliations

                    Servei de Neurologia, Fundació Privada Institut d'Investigació Biomèdica de Girona Dr. Josep Trueta (IdIBGi), Hospital Universitari de Girona Dr. Josep Trueta
                    Departamento de Farmacología, Facultad de Medicina, Universidad Complutense de Madrid (UCM)


                    1. Goldberg MP, Choi DW: Combined oxygen and glucose deprivation in cortical cell culture: calcium-dependent and calcium-independent mechanisms of neuronal injury. J Neurosci. 1993, 13: 3510-3524.PubMed
                    2. Brint S, Jacewicz M, Kiessling M, Tanabe J, Pulsinelli W: Focal brain ischemia in the rat: methods for reproducible neocortical infarction using tandem occlusion of the distal middle cerebral and ipsilateral common carotid arteries. J Cereb Blood Flow Metab. 1988, 8: 474-485.View ArticlePubMed
                    3. Tamura A, Graham DI, McCulloch J, Teasdale GM: Focal cerebral ischaemia in the rat: 1. Description of technique and early neuropathological consequences following middle cerebral artery occlusion. J Cereb Blood Flow Metab. 1981, 1: 53-60.View ArticlePubMed
                    4. Tamura A, Graham DI, McCulloch J, Teasdale GM: Focal cerebral ischaemia in the rat: 2. Regional cerebral blood flow determined by [14C]iodoantipyrine autoradiography following middle cerebral artery occlusion. J Cereb Blood Flow Metab. 1981, 1: 61-69.View ArticlePubMed
                    5. Krieger DW, Yenari MA: Therapeutic hypothermia for acute ischemic stroke: what do laboratory studies teach us?. Stroke. 2004, 35: 1482-1489. 10.1161/01.STR.0000126118.44249.5cView ArticlePubMed
                    6. Schwarz DA, Barry G, Mackay KB, Manu F, Naeve GS, Vana AM, Verge G, Conlon PJ, Foster AC, Maki RA: Identification of differentially expressed genes induced by transient ischemic stroke. Brain Res Mol Brain Res. 2002, 101: 12-22. 10.1016/S0169-328X(02)00135-3View ArticlePubMed
                    7. Caso JR, Pradillo JM, Hurtado O, Lorenzo P, Moro MA, Lizasoain I: Toll-like receptor 4 is involved in brain damage and inflammation after experimental stroke. Circulation. 2007, 115: 1599-1608. 10.1161/CIRCULATIONAHA.106.603431View ArticlePubMed
                    8. Hurtado O, Pradillo JM, Fernandez-Lopez D, Morales JR, Sobrino T, Castillo J, Alborch E, Moro MA, Lizasoain I: Delayed post-ischemic administration of CDP-choline increases EAAT2 association to lipid rafts and affords neuroprotection in experimental stroke. Neurobiol Dis. 2008, 29: 123-131. 10.1016/j.nbd.2007.08.004View ArticlePubMed
                    9. Morales JR, Ballesteros I, Deniz JM, Hurtado O, Vivancos J, Nombela F, Lizasoain I, Castrillo A, Moro MA: Activation of liver x receptors promotes neuroprotection and reduces brain inflammation in experimental stroke. Circulation. 2008, 118: 1450-1459. 10.1161/CIRCULATIONAHA.108.782300View ArticlePubMed
                    10. de Jonge HJ, Fehrmann RS, de Bont ES, Hofstra RM, Gerbens F, Kamps WA, de Vries EG, Zee van der AG, te Meerman GJ, ter Elst A: Evidence based selection of housekeeping genes. PLoS ONE. 2007, 2: e898- 10.1371/journal.pone.0000898PubMed CentralView ArticlePubMed
                    11. 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-0496View ArticlePubMed
                    12. 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- 10.1186/gb-2002-3-7-research0034PubMed CentralView ArticlePubMed
                    13. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262View ArticlePubMed
                    14. Coulson DT, Brockbank S, Quinn JG, Murphy S, Ravid R, Irvine GB, Johnston JA: Identification of valid reference genes for the normalization of RT qPCR gene expression data in human brain tissue. BMC Mol Biol. 2008, 9: 46- 10.1186/1471-2199-9-46PubMed CentralView ArticlePubMed
                    15. Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29: 332-337.PubMed
                    16. Meldgaard M, Fenger C, Lambertsen KL, Pedersen MD, Ladeby R, Finsen B: Validation of two reference genes for mRNA level studies of murine disease models in neurobiology. J Neurosci Methods. 2006, 156: 101-110. 10.1016/j.jneumeth.2006.02.008View ArticlePubMed
                    17. Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7View ArticlePubMed
                    18. Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques. 2004, 37: 112-119.PubMed
                    19. 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.022View ArticlePubMed
                    20. Derks NM, Muller M, Gaszner B, Tilburg-Ouwens DT, Roubos EW, Kozicz LT: Housekeeping genes revisited: Different expressions depending on gender, brain area and stressor. Neuroscience. 2008, 156: 305-309. 10.1016/j.neuroscience.2008.07.047View ArticlePubMed
                    21. Medhurst AD, Harrison DC, Read SJ, Campbell CA, Robbins MJ, Pangalos MN: The use of TaqMan RT-PCR assays for semiquantitative analysis of gene expression in CNS tissues and disease models. J Neurosci Methods. 2000, 98: 9-20. 10.1016/S0165-0270(00)00178-3View ArticlePubMed
                    22. Harrison DC, Medhurst AD, Bond BC, Campbell CA, Davis RP, Philpott KL: The use of quantitative RT-PCR to measure mRNA expression in a rat model of focal ischemia – caspase-3 as a case study. Brain Res Mol Brain Res. 2000, 75: 143-149. 10.1016/S0169-328X(99)00305-8View ArticlePubMed
                    23. Tian YF, Zhang PB, Xiao XL, Zhang JS, Zhao JJ, Kang QY, Chen XL, Qiu F, Liu Y: The quantification of ADAMTS expression in an animal model of cerebral ischemia using real-time PCR. Acta Anaesthesiol Scand. 2007, 51: 158-164. 10.1111/j.1399-6576.2006.01161.xView ArticlePubMed
                    24. 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-9View ArticlePubMed
                    25. Cook NL, Vink R, Donkin JJ, van den HC: Validation of reference genes for normalization of real-time quantitative RT-PCR data in traumatic brain injury. J Neurosci Res. 2009, 87: 34-41. 10.1002/jnr.21846View ArticlePubMed
                    26. 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-53PubMed CentralView ArticlePubMed
                    27. Bond BC, Virley DJ, Cairns NJ, Hunter AJ, Moore GB, Moss SJ, Mudge AW, Walsh FS, Jazin E, Preece P: The quantification of gene expression in an animal model of brain ischaemia using TaqMan real-time RT-PCR. Brain Res Mol Brain Res. 2002, 106: 101-116. 10.1016/S0169-328X(02)00417-5View ArticlePubMed
                    28. Buga AM, Sascau M, Pisoschi C, Herndon JG, Kessler C, Popa-Wagner A: The genomic response of the ipsilateral and contralateral cortex to stroke in aged rats. J Cell Mol Med. 2008, 12: 2731-2753. 10.1111/j.1582-4934.2008.00252.xPubMed CentralView ArticlePubMed
                    29. Spiegler M, Villapol S, Biran V, Goyenvalle C, Mariani J, Renolleau S, Charriaut-Marlangue C: Bilateral changes after neonatal ischemia in the P7 rat brain. J Neuropathol Exp Neurol. 2007, 66: 481-490. 10.1097/01.jnen.0000263875.22306.3cView ArticlePubMed
                    30. Keyvani K, Witte OW, Paulus W: Gene expression profiling in perilesional and contralateral areas after ischemia in rat brain. J Cereb Blood Flow Metab. 2002, 22: 153-160. 10.1097/00004647-200202000-00003View ArticlePubMed
                    31. Rhinn H, Marchand-Leroux C, Croci N, Plotkine M, Scherman D, Escriou V: Housekeeping while brain's storming Validation of normalizing factors for gene expression studies in a murine model of traumatic brain injury. BMC Mol Biol. 2008, 9: 62- 10.1186/1471-2199-9-62PubMed CentralView ArticlePubMed
                    32. McNeill RE, Miller N, Kerin MJ: Evaluation and validation of candidate endogenous control genes for real-time quantitative PCR studies of breast cancer. BMC Mol Biol. 2007, 8: 107- 10.1186/1471-2199-8-107PubMed CentralView ArticlePubMed
                    33. Bonefeld BE, Elfving B, Wegener G: Reference genes for normalization: a study of rat brain tissue. Synapse. 2008, 62: 302-309. 10.1002/syn.20496View ArticlePubMed
                    34. Tanic N, Perovic M, Mladenovic A, Ruzdijic S, Kanazir S: Effects of aging, dietary restriction and glucocorticoid treatment on housekeeping gene expression in rat cortex and hippocampus-evaluation by real time RT-PCR. J Mol Neurosci. 2007, 32: 38-46. 10.1007/s12031-007-0006-7View ArticlePubMed
                    35. Sobrado M, Pereira MP, Ballesteros I, Hurtado O, Fernandez-Lopez D, Pradillo JM, Caso JR, Vivancos J, Nombela F, Serena J, et a;.: Synthesis of lipoxin A4 by 5-lipoxygenase mediates PPARgamma-dependent, neuroprotective effects of rosiglitazone in experimental stroke. J Neurosci. 2009, 29: 3875-3884. 10.1523/JNEUROSCI.5529-08.2009View ArticlePubMed
                    36. Candelario-Jalil E, Fiebich BL: Cyclooxygenase inhibition in ischemic brain injury. Curr Pharm Des. 2008, 14: 1401-1418. 10.2174/138161208784480216View ArticlePubMed
                    37. Dirnagl U, Simon RP, Hallenbeck JM: Ischemic tolerance and endogenous neuroprotection. Trends Neurosci. 2003, 26: 248-254. 10.1016/S0166-2236(03)00071-7View ArticlePubMed
                    38. Matsushita K, Matsuyama T, Kitagawa K, Matsumoto M, Yanagihara T, Sugita M: Alterations of Bcl-2 family proteins precede cytoskeletal proteolysis in the penumbra, but not in infarct centres following focal cerebral ischemia in mice. Neuroscience. 1998, 83: 439-448. 10.1016/S0306-4522(97)00391-6View ArticlePubMed
                    39. Chen J, Graham SH, Chan PH, Lan J, Zhou RL, Simon RP: bcl-2 is expressed in neurons that survive focal ischemia in the rat. Neuroreport. 1995, 6: 394-398.View ArticlePubMed
                    40. McCarthy KD, de Vellis J: Preparation of separate astroglial and oligodendroglial cell cultures from rat cerebral tissue. J Cell Biol. 1980, 85: 890-902. 10.1083/jcb.85.3.890View ArticlePubMed
                    41. Mallolas J, Hurtado O, Castellanos M, Blanco M, Sobrino T, Serena J, Vivancos J, Castillo J, Lizasoain I, Moro MA, et al.: A polymorphism in the EAAT2 promoter is associated with higher glutamate concentrations and higher frequency of progressing stroke. J Exp Med. 2006, 203: 711-717. 10.1084/jem.20051979PubMed CentralView ArticlePubMed
                    42. Hurtado O, Moro MA, Cardenas A, Sanchez V, Fernandez-Tome P, Leza JC, Lorenzo P, Secades JJ, Lozano R, Davalos A, et al.: Neuroprotection afforded by prior citicoline administration in experimental brain ischemia: effects on glutamate transport. Neurobiol Dis. 2005, 18: 336-345. 10.1016/j.nbd.2004.10.006View ArticlePubMed


                    © Gubern et al; licensee BioMed Central Ltd. 2009

                    This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.