Skip to main content

Housekeeping while brain's storming Validation of normalizing factors for gene expression studies in a murine model of traumatic brain injury

Abstract

Background

Traumatic brain injury models are widely studied, especially through gene expression, either to further understand implied biological mechanisms or to assess the efficiency of potential therapies. A large number of biological pathways are affected in brain trauma models, whose elucidation might greatly benefit from transcriptomic studies. However the suitability of reference genes needed for quantitative RT-PCR experiments is missing for these models.

Results

We have compared five potential reference genes as well as total cDNA level monitored using Oligreen reagent in order to determine the best normalizing factors for quantitative RT-PCR expression studies in the early phase (0–48 h post-trauma (PT)) of a murine model of diffuse brain injury. The levels of 18S rRNA, and of transcripts of β-actin, glyceraldehyde-3P-dehydrogenase (GAPDH), β-microtubulin and S100β were determined in the injured brain region of traumatized mice sacrificed at 30 min, 3 h, 6 h, 12 h, 24 h and 48 h post-trauma.

The stability of the reference genes candidates and of total cDNA was evaluated by three different methods, leading to the following rankings as normalization factors, from the most suitable to the less: by using geNorm VBA applet, we obtained the following sequence: cDNA(Oligreen); GAPDH > 18S rRNA > S100β > β-microtubulin > β-actin; by using NormFinder Excel Spreadsheet, we obtained the following sequence: GAPDH > cDNA(Oligreen) > S100β > 18S rRNA > β-actin > β-microtubulin; by using a Confidence-Interval calculation, we obtained the following sequence: cDNA(Oligreen) > 18S rRNA; GAPDH > S100β > β-microtubulin > β-actin.

Conclusion

This work suggests that Oligreen cDNA measurements, 18S rRNA and GAPDH or a combination of them may be used to efficiently normalize qRT-PCR gene expression in mouse brain trauma injury, and that β-actin and β-microtubulin should be avoided.

The potential of total cDNA as measured by Oligreen as a first-intention normalizing factor with a broad field of applications is highlighted. Pros and cons of the three methods of normalization factors selection are discussed. A generic time- and cost-effective procedure for normalization factor validation is proposed.

Background

Real-time RT-PCR, which allows to measure any chosen RNA with great accuracy over a large dynamic range, has become the gold-standard for nucleic acid quantification. It has also opened new investigations fields, since very small amount of RNA is needed, allowing transcripts from low-expressed genes or from very small samples to be quantified.

If constant developments in both reagents and data analysis make real-time PCR measurements more and more accurate and reliable, many considerations have to be taken to convert this technical precision into biologically relevant data. Factually, real-time RT-PCR gives access to the number of copies of a chosen sequence in a cDNA solution, which is obtained from RNA extracted from a known quantity of tissue. The biologically relevant information that has to be ultimately obtained is the global expression level of the chosen gene in the tested sample, at least relatively to another sample.

The quantification of a target gene in a given sample needs three majors steps: RNA/mRNA extraction, reverse transcription of the extracted RNA, and qPCR (quantitative PCR) processing of the synthesised cDNA. A control normalization may be performed at each step to level out dissimilarities between samples [1].

The first possible normalization is by equalizing samples size, such as cell number or tissue weight. This is the easiest and the most intuitive measure. However, its position upstream of the reactions sequence does not allow to correct for the distortions generated by downstream manipulations, especially by RNA extraction, whose efficiency may broadly vary from one sample to another.

The second method consists in normalizing samples according to RNA content after its extraction. This however does not take into account the reverse transcription efficiency, known to vary from one sample to the other [2].

Thus, a downstream normalization method appears to be the most effective. This may be performed by measuring the expression level of a gene transcript expressed in the sample, as an endogenous control for the different reaction steps. The housekeeping term, which is often applied to these genes, was initially given to genes that are necessary for the function of each cell. As a matter of fact, they have to be expressed in each cell type. The most typical case is β-actin, a cornerstone of the inner architecture of the cell. This makes housekeeping genes suitable for organism-wide positive controls for many cDNA-based techniques, but does not ensure their expression levels to be equals. The expression levels of usual housekeeping genes has however been shown to vary in some conditions [1].

If the belief in the existence of perfect reference genes, whose levels would remain unchanged in each cell whatever the tissue or the experimental conditions, is known to be more idealistic than real, reference genes have to be chosen specifically for a given experiment, on the basis of the stability of their expression in the subset of studied tissues and experimental conditions one is interested in. In consequence, the use of internal controls implies a proper validation for each experimental condition, as the use of unappropriate normalizing factors, with unconstant expression levels, would generate discrepancies in normalized expression results [3]. Several methods have been described for that purpose [47].

In addition to the use of reference genes, alternative normalization procedures have been proposed for RT-PCR: the addition of known quantity of artificial RNA molecules to extracted RNA prior to RT reaction [8], or the use of an oligo-dT linked artificial target sequence quantifiable by quantitative PCR (qPCR) [9] or the total cDNA quantification by fluorescent dyes [10].

In the specific context of brain trauma, sample size measurement is especially doubtful, since the oedema resulting from the trauma [11] may enhance brain water content, hence lowering the cellular density of a given weight of tissue. This effect may also be strengthened, as necrosis and apoptosis are known to be major post-trauma events [12, 13].

We have reviewed the use of PCR to quantify mRNA levels in mouse and rat traumatic mechanical brain injury models for the five last years. We have identified 22 qPCR and 17 semi quantitative PCR studies. About 85% of them used a reference gene for normalization, while the others only used RNA quantification for this purpose.

GAPDH has been up to now the most commonly used reference gene in qPCR [1424] and the second one [2528] for sqPCR (semi-quantitative PCR). The use of β-actin has been predominant for sqPCR [2936], but less frequent for qPCR, to the benefit of 18S rRNA [3741]. We have found no correlation between the choice of a reference gene and the pathological model (Controlled Cortical Injury (CCI), Closed Head Injury (CHI) or Fluid Percussion Injury (FPI)) or species. Strikingly, a single reference gene has been used in all the experiments, without reference to a prior validation, whereas the use of more than one reference gene is widely encouraged [4], and the need for validation highlighted [1, 3]. Indeed, the choice of reference genes is especially crucial for brain trauma studies, where many biological pathways are implicated (inflammation, hypoxia, apoptosis, neovascularisation...), resulting in numerous gene expression changes, which are thus likely to include or affect potential reference genes.

The present study aims at the validation of optimal nomalizing factors(NF) for RT-qPCR-based transcriptional studies of an early phase of a murine model of traumatic diffuse brain injury.

RT-qPCR was used to determine the expression levels of five potential reference genes. Total cDNA level was also measured using Oligreen reagent and considered as an other potential NF. These six candidates were compared using different NF-selection methods.

Results

Reference genes and total cDNA level

After induction of trauma, 10 mice were sacrificed at 30 min, 3 h, 6 h, 12 h, 24 h and 48 h, and RNA was extracted from the lesionnal zone of injured brains or from the equivalent zone of uninjuried control mice brains, then reverse transcribed. The expression level of 18S rRNA, β-microtubulin, S100β, β-actin and GAPDH were measured by real-time PCR for each individual cDNA.

Oligreen reagent is a dye that emits fluorescence at 500/520 nm when bound to single-stranded DNA, but is insensitive to free nucleotides and to very short oligonucleotides. However, Oligreen emits fluorescence in the presence of RNA at room temperature; its use for cDNA quantification therefore implies to work at 80°C to specifically quantify reverse-transcribed cDNA [43]. Levels of reference genes at different time post-trauma, as measured by qPCR and of total cDNA as measured using Oligreen reagent, are represented on Fig. 1.

Figure 1
figure 1

Level of normalizing factors at various time post-trauma. The levels of 18S rRNA, β-microtubulin, S100β, β-actin and GAPDH as measured by real-time PCR, and cDNA level as measured with Oligreen were determined at the indicated time point post-trauma. Each point corresponds to the measured level for one individual sample, normalized to the mean of all individual values at all times Horizontal bars represent the median of the measures for the group, vertical crosses represent the arithmetic mean. (C): control.

ANOVA

A statistical test was applied to look for significant differences between two experimental conditions for each NF level. As the variance equality hypothesis was verified for each of them (no significant variance differences detected with a Hartley test for a confidence level of 95%), a one-factor ANOVA test was used, with a Fischer's test to detect significant differences (p < 0.05) between two groups, for a given NF expression level.

The expression level fluctuations between time groups were found to be significant for β-actin, S100β and β-microtubulin. No significant differences were detected from one group to another for GAPDH, 18S rRNA and cDNA (Oligreen).

Significant differences between experimental groups clearly make the tested NF unsuitable to normalize sample of the tested experimental set. The expression level of β-actin, S100β and β-microtubulin should thus be avoided as a NF for the present experimental model.

To further study the suitability of the studied normalization factors, three published selection methods were applied: geNorm analysis [4], Normfinder analysis [6] and a confidence interval based method [7].

GeNorm analysis

Gene stability analysis was performed using the geNorm VBA applet as described in [4]. Briefly, the basic assumption of this method is that the ratio of two perfect reference genes should be constant throughout the different experimental conditions. The inter-condition variability of this ratio is thus evaluated for each experimental condition and for each couple of reference gene, and a gene-stability-measure M is calculated for each candidate. This reflects the average pair-wise variability between one putative reference gene and all the others. The less stable candidate (i.e. with the highest M value) is then excluded as the least suitable NF, and a new step of M-values calculation has then to be performed because the former inclusion of the less stable gene had influenced other M values. By sequential exclusions, the two most-stably expressed genes are selected. The ranking of the studied NF according to their M value, as calculated by the GeNorm software was, from the most stable to the least:

GAPDH – cDNA(Oligreen)>18S rRNA > S100β > β-actin > β-microtubulin (Fig. 2A).

Figure 2
figure 2

Genorm output. A: Expression stability values (M) of the candidate normalization factors, during the stepwise exclusion of least stable normalization factors. Candidate normalization factors are ranked from the least stable to the most stable (left to right). B: Determination of the optimal number of normalization factors to be used, based on the analysis of pairwise variation (V) between normalization factors. A low value of V means that the inclusion of an additional normalization factor would not lead to a increase in stability.

As proposed by Vandersompele et al. [4], the use of a geometric mean of multiple NF among the best ranked should ensure a more accurate normalization. The choice of NF to be included starts from the most stable NF, by a step wise inclusion of the next-ranking NF and the influence of the inclusion on the overall stability is plotted by the GeNorm software (Fig. 2B). A strong instability elevation implies that the inclusion is detrimental, and that the last included NF should thus be discarded. In our case, GeNorm method leads to the use of the geometric mean of total cDNA and GAPDH levels for the best normalization.

- Normfinder analysis

A model-based approach, was used as an another technique to rank the best potential reference genes, using the Excel Spreadsheet Normfinder [6]. The criteria for sample sizes (n>8) and candidate RG (n>3) were met. Briefly, for one given sample and one given candidate reference gene, the log-transformed measured level is formulated as the sum of three terms: the general expression level of the gene in the experimental group to which the sample belongs, the amount of mRNA in the sample, and the random variation caused by biological and experimental factors. The latter term is used to calculate an intra-group variation, while the first term permits to evaluate an inter-group variability for the considered NF. Both variations estimations are combined in a stability value used to rank the NF.

The ranking of each NF according to its stability value as calculated by the Normfinder software was, from the most stable to the least:

GAPDH>cDNA(Oligreen)> S100β >18S rRNA > β-actin > β-microtubulin (Fig. 3).

Figure 3
figure 3

Normfinder stability values output.

- Confidence-Interval method

A third approach has been used, based on the one presented in Haller et al. [7]. Equivalence testing is similar to classical statistical tests, but while most of them are based on the rejection of a relative closeness hypothesis to prove a significant difference between two means, equivalence testing relies on the rejection of an hypothese of a relative difference between two means to prove their significant closeness.

This can be tested by classical t-test, but, for a given level of confidence alpha, one can proceed more visually through a confidence interval calculation, which is the most-used method for statistical assessment of clinical trials.

Briefly, two samples Ti and Tj, with expression levels μi and μj, are considered to be equivalent for a confidence level alpha if their confidence interval Iij, centered on μi - μj and whose length is based on a Student t-distribution depending of the distributions of the two samples (detailed in Haller et al.) is such as I ij [-ε, ε], epsilon being an arbitrarly chosen treshold. Moreover, if 0 I ij , the two groups are significantly different for the confidence level.

When dealing with expression levels, intrinsic meaning has to be found in ratio rather than in absolute difference. We thus worked with log-transformed expression levels for equivalence testing. In this way, the limit epsilon that has to be chosen is equivalent to a fold-change. We chose to set it to 1, which means an expression ratio between the two compared groups inferior to two.

The confidence intervals for log-transformed normalized NF levels were calculated for each couple of experimental groups and for all NF, with 95% confidence levels (Fig. 4). Two groups were said to be equivalent with 95% confidence if their confidence level is included in [-1;1].

Figure 4
figure 4

95% Confidence intervals for equivalence between experimental groups. TG/TG' corresponds to the interval calculated between group G and group G'. With T0: control group and T1, T2, T3, T4, T5 and T6 trauma groups sacrificed at 30 min, 3 h, 6 h, 12 h, 24 h and 48 h respectively. For symmetrical reasons, since Iij = Iji, only (Ii,j)i≤j are plotted. For each of the normalization factor: If 0 I ij , its level is significantly different between experimental conditions i and j. If I ij [-ε; ε], its level is equivalent between experimental conditions i and j. A: Oligreen measurements; B: 18S rRNA; C: S100β; D: GAPDH; E: β-actin; F: β-microtubulin.

Measured levels were found to be equivalent between all groups for cDNA(Oligreen).

For 18S rRNA and GAPDH, the expression levels were found not to be equivalent for two couples of groups out of twenty-one. It was the case for 8 combinations for S100β, 18 for β-microtubulin and 19 for β-actin (Fig. 5) An equivalence-based ranking thus appears to be:

Figure 5
figure 5

Schematic results for confidence interval calculations. Confidence interval are calcultated using a cut-off fold change = 2 and 95% confidence level. White squares: the two groups are equivalent for the considered NF. Grey squares: the two groups are not equivalent for the considered NF. Dashed squares: the two groups are significantly different for the considered NF. With T0: control group and T1, T2, T3, T4, T5 and T6 trauma groups sacrificed at 30 min, 3 h, 6 h, 12 h, 24 h and 48 h respectively.

cDNA(Oligreen)> 18S rRNA-GAPDH> S100β > β-microtubulin > β-actin

Discussion

We have used the Oligreen reagent to quantify the cDNA effectively synthesised during the RT step. Whatever the procedure for NF choice, cDNA as measured by Oligreen ranked among the two best. As its use appears to be validated for the present model, it should also be very attractive for other cases, since small variability is unlikely to depend on experimental conditions, in contrast to other reference genes.

To ensure that the overall high stability of Oligreen-assayed cDNA is not an artefact, we compared for each sample the mean of all studied reference genes levels to the measured Oligreen level. The aim is to ensure that Oligreen reflects the level of cDNA in each sample. The application of the previously used selection procedures for RNA input – equal by definition to 1 μ g for each sample – would indeed have led to the systemic selection of RNA quantification as the best normalizing factor, with both intra- and inter- group variation equal to zero. RNA level is however obviously not an appropriate normalizer, as it does not encompasses RT efficiency. As the mean of all the gene expression levels measured is likely to be proportionnal to the cDNA content, studying the correlation between this average expression level and the Oligreen level for each sample appears as a good mean to evaluate whether Oligreen measurement is proportionnal to the cDNA content or not. A significant correlation was found using a Pearson test (p-value < 0.0001, Pearson coefficient: 0.64) and the correlation plot (Fig. 6) shows that Oligreen measurements effectively reflect cDNA levels in the samples, as it correlates well with the the mean of different genes expression levels.

Figure 6
figure 6

Correlation between Oligreen measurements and other reference genes. Each point corresponds to one sample, as the mean of normalized measured levels for all tested reference genes in function of the normalized measured levels of Oligreen. A 95% confidence ellipse is drawn.

Oligreen thus appears to be a good candidate for a generic NF in quantitative RT-PCR experiments. Its use may be furthermore refined with alternative RT protocols, such as an oligo dT-primed cDNA synthesis. In such a case, Oligreen would reflect the level of effectively reverse-transcribed mRNA, in contrast to total RNA in the present study.

The use of Oligreen combined to a standard curve of known DNA concentration could moreover allows a precise cDNA quantification. In parallel to an absolute quantification of target genes, one could thus formulate the expression level of a studied gene, as the number of copies per microgram of cDNA, which would have a stronger intrinsic meaning and could open the gate for easier intra-experiments meta-analysis.

The question remains open whether the fact of combining Oligreen measurements to the expression level of internal control (18S rRNA and GAPDH in our case) brings more accuracy.

Following a confidence-interval based methodology, Oligreen appeared as the only NF for which all groups were found to be equivalent to each other, but one has to keep in mind the arbitrary choice of fold-change cut-off level and confidence levels: 18S rRNA and GAPDH appear to be very close to equivalence and thus appear as acceptable normalizing factors.

When working with confidence levels, an arbitrary level of fold change has to be set. A possible criterion for assessing the pertinence of a chosen cut-off fold-change may come from the observation of confidence intervals Iii calculated for a group with itself: as a group is trivially equivalent to itself, the calculated confidence interval is nothing but the reflection of the group's standard deviation. It would thus make little sense to choose a cut-off fold change that would not lead to the conclusion that a group is equivalent to itself. This consideration excluded the initially-considered choice of 1.5 as a cut-off fold-change in our case, and led to the use of 2 as a cut-off fold-change for the study. A standardized confidence interval approach, which remains to be developped, should thus set the cut-off fold change in function of the intra-groups confidence intervals. One has however to keep in mind that alpha level is arbitrary too, and that its value has a direct effect (through t-distribution) on the interval size: lowering cut-off fold change is in some extent equivalent to setting a higher confidence level alpha. These two values have thus to be chosen in a concerted manner to make sense.

Among the different procedures of normalizing factors that were applied, two criteria appear to be determinant. On one hand, one criterion is whether the method takes into account the belonging of a sample to an experimental group (either a tissue, a differentially treated sample, or a given point of a time series). This should ultimately be a requirement, as the final goal of many, if not all, gene expression experiments consists in levels comparison between different groups. This is handled by model-based and confidence intervals based methodology, but not by pair-wise approach. However, if the Normfinder software evaluates inter-group variability, its effect is combined to intra-groups variability to calculate a candidate's stability. This implies that a NF with a low variability in each group, but with an expression level clearly different in one single group, will be considered as the optimal NF, when compared to candidate with greater variability in each group but with almost no differences between groups (data not shown). This is also the case for the GeNorm software, and is probably the worst type of error that may occur while selecting a NF, as it will generate systematic errors in normalized gene expression results.

On the other hand, a second determinant criterion for eliciting a NF is the absolute or relative character of its evaluation as a potential NF. This means whether its sole expression values, as measured for every sample, are sufficient or if it has to be compared to other potential normalizing factor. The relative comparison is clearly the basis of the pair-wise method, but is also implicitly present in the model-base approach for the calculation of both intra- and inter- group variations (through the use of γ ^ MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafq4SdCMbaKaaaaa@2D8C@ , variance of α ig in the calculation of the stability value ρ i (Eq. C and D, [6])). The main drawback of these circular methodologies is the risk of selecting co-regulated genes. One would thus need to examinate many carefully chosen potential NF (i.e. from biological pathways as distinct as possible) in order to ensure that a majority of them is not co-regulated, the quality of the NF selection being strongly dependant of the number of studied candidates.

The equivalence interval approach enables to evaluate each NF candidate separately. This permits to imagine a sequential approach to assess reference gene: rather than gathering an as large as possible pool of potential NF candidates, then measuring their expression level and sorting them to determine the most appropriate among them, one may choose a first-intention potential NF and test it ab nihilo. If the test concords with previously fixed criteria, then the candidate is validated as a NF, otherwise, an additional candidate has to be tested. As it seems best to normalize with more than one reference gene, one may stop the selection procedure when two or three candidates are positively evaluated.

An additional interesting feature of an absolute evaluation is that it enables researchers to proceed to a post-hoc validation of their previously chosen NF on the sole basis of already-existing data analysis, with no need to run additional experiments with additional NF. This appears to be especially interesting when the biological material was totally used for previous reactions.

Despite its less restricted use, probably due to the lack of an easy-to-handle Excel spreadsheet as it is the case with GeNorm and Normfinder, the confidence interval based method appears to be the most powerful and most adequate for NF selection. A scoring procedure based on this approach, allowing a ranking of different NF, along with an automated easy-to-handle automated application, may be an interesting tool to develop.

Conclusion

In conclusion, we present here the first study aiming at the identification of optimal normalizing factors in a model of traumatic brain injury. If the use of β-actin and β-microtubulin appears to be avoided, the combination of different methods leads us to suggest the use of a geometric mean of 18S rRNA, GAPDH and total cDNA as measured with Oligreen. The use of cDNA measurement with Oligreen has been validated in the present case and is encouraged as first-intention generic normalization factor. The present study also highlights the interest of using confidence interval for normalizing factor accuracy assessing, and opens the way for an iterative, time and cost-effective normalizing factor selection procedure.

Methods

Closed head injury model

Male CD-1 mice (IFFA Credo, France), weighting 21–24 g, were housed at room temperature under controlled light conditions (12 h light: 12 h dark) with food and water ad libitum. The murine model of diffuse head injury was used as previously described [11]. Mice were briefly anesthetized under 2% halothane balanced with air and oxygen. Closed head trauma was induced by a 50 g weight dropped from 32-cm height along a stainless steel string. This model typically results in a mortality of 30% within the first 10 min following the impact, with no observed delayed mortality. Animal care was in compliance with French regulations on protection of animals used for experimental research and with the EC regulations (Official Journal of European Community L35812/18/1986).

cDNA preparation

Mice were sacrificed at different times (30 min, 3 h, 6 h, 12 h, 24 h and 48 h post-trauma). Brains were immediately removed. Column-shaped samples of approximately 10 mg were taken vertically around the lesion site using a 4 mm punch and put in RNAlater solution (QIAGEN), and were kept at 4°C for further extraction of total RNA using the Rneasy kit (QIAGEN). RNA concentration was assessed using a NanoDrop ND-1000 Spectrophotometer (Nanodrop). One microgram of total RNA was reverse transcribed in a final volume of 20 μ l containing 4 μ L of 5× RT buffer, 20 units of RNasin RNase inhibitor (Promega), 10 mM DDT, 100 units of Superscript II RNase H-reverse transcriptase (Invitrogen), 3 mM random hexamers (Invitrogen). The samples were incubated at 20°C for 10 min, 42°C for 30 min and 99°C for 5 min.

Real-time PCR

PCR primers for tested reference gene were chosen in published articles, for their common use as reference genes and their belonging to different biological pathways (Fig. 7). Real-time PCR reactions were carried out using ABI PRISM 7900HT Sequence Detection System (Applied Biosystems) in a 384-well, clear optical reaction plate with optical adhesive covers (Applied Biosystems). Reactions were run in a 5 μ l volume in duplicate, with 2 μ L of cDNA solution and 3 μ L of a homemade target-specific mix composed of 5/6 2× Power SYBR Green Master Mix (Applied Biosystems) and 1/6 of 100 mM primers solution. The PCR program was: 95°C for 10 min, followed by 45 cycles of (15 seconds at 95°C; 1 min at 60°C). Product specificity was assessed by 3% agarose gel electrophoresis followed by ethidium bromide staining.

Figure 7
figure 7

Chosen primers.

For each PCR well, raw fluorescence data were exported and a curve-fitting method was applied to evaluate individual efficiencies [42]. As for each target, calculated efficiencies appeared to be homogeneous throughout samples, the efficiency for each target was considered as shared for all wells. This efficiency was set as the mean of the fitted individual efficiencies. Individual Ct values exported from the SDS software (Applied Biosystems) were used to calculate the relative expression level of target t in sample s: L t,s = e f f t C t ( s , t ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaeeitaW0aaSbaaSqaaiabbsha0jabbYcaSiabbohaZbqabaGccqGH9aqpcqWGLbqzcqWGMbGzcqWGMbGzdaqhaaWcbaGaemiDaqhabaGaeyOeI0Iaem4qamKaemiDaqNaeiikaGIaem4CamNaeiilaWIaemiDaqNaeiykaKcaaaaa@4065@

As dealing with numerous normalizing factors, these relative expression levels were normalized for each target, for easier calculation and comparison.

Oligreen measurements

Fluorescence measurements was carried out using ABI PRISM 7900HT Sequence Detection System (Applied Biosystems) in a 384-well, clear optical reaction plate with optical adhesive covers (Applied Biosystems). The reactions were run in a 5 μ l volume in duplicate, with 2 μ L of cDNA solution, 3 μ L of Oligreen stock solution (Invitrogen) 200× diluted in TE buffer. Fluorescence was continuously read at 80°C for 1 min, a condition which ensures exclusive cDNA measurement [43]. Measurements linearity was assessed by a cDNA standard curve. After background substraction, fluorescence values were normalized throughout all samples to be treated in the same way as PCR-based NF values.

Data analysis

Statistical analysis

Statistical calculations were run with Microsoft Excel and XLSTAT 2006 add-on (Addinsoft).

Confidence interval calculation

For each reference gene, groups are tested two-by-two for their equivalences. For two groups (Ti;Tj) with the same number of elements (N), with means ( X ^ i ; X ^ j MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGafmiwaGLbaKaadaWgaaWcbaGaemyAaKgabeaakiabcUda7iqbdIfayzaajaWaaSbaaSqaaiabdQgaQbqabaaaaa@327F@ ) and standard deviations (S i ;S j ); the confidence interval is calculated on a Microsoft Excel spreadsheet using the formula:

I i j = [ X ^ i X ^ j t 2 N 2 ; 1 α S i 2 + S j 2 N ; X ^ i X ^ j + t 2 N 2 ; 1 α S i 2 + S j 2 N ] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemysaK0aaSbaaSqaaiabdMgaPjabdQgaQbqabaGccqGH9aqpdaWadaqaaiqbdIfayzaajaWaaSbaaSqaaiabdMgaPbqabaGccqGHsislcuWGybawgaqcamaaBaaaleaacqWGQbGAaeqaaOGaeyOeI0IaemiDaq3aaSbaaSqaaiabikdaYiabd6eaojabgkHiTiabikdaYiabcUda7iabigdaXiabgkHiTiabeg7aHbqabaGcdaGcaaqaaKqbaoaalaaabaGaem4uam1aa0baaeaacqWGPbqAaeaacqaIYaGmaaGaey4kaSIaem4uam1aa0baaeaacqWGQbGAaeaacqaIYaGmaaaabaGaemOta4eaaaWcbeaakiabcUda7iqbdIfayzaajaWaaSbaaSqaaiabdMgaPbqabaGccqGHsislcuWGybawgaqcamaaBaaaleaacqWGQbGAaeqaaOGaey4kaSIaemiDaq3aaSbaaSqaaiabikdaYiabd6eaojabgkHiTiabikdaYiabcUda7iabigdaXiabgkHiTiabeg7aHbqabaGcdaGcaaqcfayaamaalaaabaGaem4uam1aa0baaeaacqWGPbqAaeaacqaIYaGmaaGaey4kaSIaem4uam1aa0baaeaacqWGQbGAaeaacqaIYaGmaaaabaGaemOta4eaaaWcbeaaaOGaay5waiaaw2faaaaa@6B75@

t 2N-2;1-α:1-α quantile of the two-tails Student t-distribution with 2N-2 degrees of freedom.

Abbreviations

NF:

Normalizing Factor

Ct:

Cycle Treshold

PCR:

Polymerase Chain Reaction

qPCR:

quantitative PCR

sqPCR:

semi-quantitative PCR

RT:

Reverse Transcription

GAPDH:

glyceraldehyde-3P-dehydrogenase

CCI:

Controlled Cortical Injury

CHI:

Closed Head Injury

FPI:

Fluid Percussion Injury.

References

  1. Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun 2005, 6: 279-284. 10.1038/sj.gene.6364190.

    Article  CAS  PubMed  Google Scholar 

  2. Stahlberg A, Kubista M, Pfaffl M: Comparison of reverse transcriptases in gene expression analysis. Clin Chem 2004, 50: 1678-1680. 10.1373/clinchem.2004.035469.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    Article  CAS  PubMed  Google Scholar 

  4. 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-research0034.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol Lett 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47.

    Article  CAS  PubMed  Google Scholar 

  6. 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.

    Article  CAS  PubMed  Google Scholar 

  7. Haller F, Kulle B, Schwager S, Gunawan B, von Heydebreck A, Sultmann H, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Anal Biochem 2004, 335: 1-9. 10.1016/j.ab.2004.08.024.

    Article  CAS  PubMed  Google Scholar 

  8. Smith RD, Brown B, Ikonomi P, Schechter AN: Exogenous reference RNA for normalization of real-time quantitative PCR. Biotechniques 2003, 34: 88-91.

    CAS  PubMed  Google Scholar 

  9. Argyropoulos D, Psallida C, Spyropoulos CG: Generic normalization method for real-time PCR. Application for the analysis of the mannanase gene expressed in germinating tomato seed. Febs J 2006, 273: 770-777. 10.1111/j.1742-4658.2006.05109.x.

    Article  CAS  PubMed  Google Scholar 

  10. Lundby C, Nordsborg N, Kusuhara K, Kristensen KM, Neufer PD, Pilegaard H: Gene expression in human skeletal muscle: alternative normalization method and effect of repeated biopsies. Eur J Appl Physiol 2005, 95: 351-360. 10.1007/s00421-005-0022-7.

    Article  CAS  PubMed  Google Scholar 

  11. Hellal F, Bonnefont-Rousselot D, Croci N, Palmier B, Plotkine M, Marchand-Verrecchia C: Pattern of cerebral edema and hemorrhage in a mice model of diffuse brain injury. Neurosci Lett 2004, 357: 21-24. 10.1016/j.neulet.2003.12.036.

    Article  CAS  PubMed  Google Scholar 

  12. Raghupathi R, Graham DI, McIntosh TK: Apoptosis after traumatic brain injury. J Neurotrauma 2000, 17: 927-938.

    Article  CAS  PubMed  Google Scholar 

  13. Friedrich DH, Plesnila N, Eriskat J, Stoffel M, Baethmann A: Cerebral blood flow and the secondary growth of brain tissue necrosis after trauma. Acta Neurochir Suppl 2000, 76: 257-259.

    CAS  PubMed  Google Scholar 

  14. Ringger NC, Tolentino PJ, McKinsey DM, Pike BR, Wang KK, Hayes RL: Effects of injury severity on regional and temporal mRNA expression levels of calpains and caspases after traumatic brain injury in rats. J Neurotrauma 2004, 21: 829-841. 10.1089/0897715041526177.

    Article  CAS  PubMed  Google Scholar 

  15. Li HH, Lee SM, Cai Y, Sutton RL, Hovda DA: Differential gene expression in hippocampus following experimental brain trauma reveals distinct features of moderate and severe injuries. J Neurotrauma 2004, 21: 1141-1153. 10.1089/0897715041953777.

    Article  PubMed  Google Scholar 

  16. Larner SF, Hayes RL, McKinsey DM, Pike BR, Wang KK: Increased expression and processing of caspase-12 after traumatic brain injury in rats. J Neurochem 2004, 88: 78-90.

    Article  CAS  PubMed  Google Scholar 

  17. Leinhase I, Schmidt OI, Thurman JM, Hossini AM, Rozanski M, Taha ME, Scheffler A, John T, Smith WR, Holers VM, Stahel PF: Pharmacological complement inhibition at the C3 convertase level promotes neuronal survival, neuroprotective intracerebral gene expression, and neurological outcome after traumatic brain injury. Exp Neurol 2006, 199: 454-464. 10.1016/j.expneurol.2006.01.033.

    Article  CAS  PubMed  Google Scholar 

  18. Giza CC, Maria NS, Hovda DA: N-methyl-D-aspartate receptor subunit changes after traumatic injury to the developing brain. J Neurotrauma 2006, 23: 950-961. 10.1089/neu.2006.23.950.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Lu KT, Wang YW, Wo YY, Yang YL: Extracellular signal-regulated kinase-mediated IL-1-induced cortical neuron damage during traumatic brain injury. Neurosci Lett 2005, 386: 40-45. 10.1016/j.neulet.2005.05.057.

    Article  CAS  PubMed  Google Scholar 

  20. Lu KT, Wang YW, Yang JT, Yang YL, Chen HI: Effect of interleukin-1 on traumatic brain injury-induced damage to hippocampal neurons. J Neurotrauma 2005, 22: 885-895. 10.1089/neu.2005.22.885.

    Article  PubMed  Google Scholar 

  21. Larner SF, McKinsey DM, Hayes RL, KK WW: Caspase 7: increased expression and activation after traumatic brain injury in rats. J Neurochem 2005, 94: 97-108. 10.1111/j.1471-4159.2005.03172.x.

    Article  CAS  PubMed  Google Scholar 

  22. Johnson EA, Svetlov SI, Pike BR, Tolentino PJ, Shaw G, Wang KK, Hayes RL, Pineda JA: Cell-specific upregulation of survivin after experimental traumatic brain injury in rats. J Neurotrauma 2004, 21: 1183-1195. 10.1089/neu.2004.21.1183.

    Article  PubMed  Google Scholar 

  23. Iino M, Nakatome M, Ogura Y, Fujimura H, Kuroki H, Inoue H, Ino Y, Fujii T, Terao T, Matoba R: Real-time PCR quantitation of FE65 a beta-amyloid precursor protein-binding protein after traumatic brain injury in rats. Int J Legal Med 2003, 117: 153-159.

    PubMed  Google Scholar 

  24. Tolentino PJ, DeFord SM, Notterpek L, Glenn CC, Pike BR, Wang KK, Hayes RL: Up-regulation of tissue-type transglutaminase after traumatic brain injury. J Neurochem 2002, 80: 579-588. 10.1046/j.0022-3042.2001.00726.x.

    Article  CAS  PubMed  Google Scholar 

  25. Shimamura M, Garcia JM, Prough DS, Dewitt DS, Uchida T, Shah SA, Avila MA, Hellmich HL: Analysis of long-term gene expression in neurons of the hippocampal subfields following traumatic brain injury in rats. Neuroscience 2005, 131: 87-97. 10.1016/j.neuroscience.2004.10.041.

    Article  CAS  PubMed  Google Scholar 

  26. Hellmich HL, Capra B, Eidson K, Garcia J, Kennedy D, Uchida T, Parsley M, Cowart J, DeWitt DS, Prough DS: Dose-dependent neuronal injury after traumatic brain injury. Brain Res 2005, 1044: 144-154. 10.1016/j.brainres.2005.02.054.

    Article  CAS  PubMed  Google Scholar 

  27. Ke C, Poon WS, Ng HK, Lai FM, Tang NL, Pang JC: Impact of experimental acute hyponatremia on severe traumatic brain injury in rats: influences on injuries, permeability of blood-brain barrier, ultrastructural features, and aquaporin-4 expression. Exp Neurol 2002, 178: 194-206. 10.1006/exnr.2002.8037.

    Article  CAS  PubMed  Google Scholar 

  28. Yakovlev AG, Ota K, Wang G, Movsesyan V, Bao WL, Yoshihara K, Faden AI: Differential expression of apoptotic protease-activating factor-1 and caspase-3 genes and susceptibility to apoptosis during brain development and after traumatic brain injury. J Neurosci 2001, 21: 7439-7446.

    CAS  PubMed  Google Scholar 

  29. Yoshiya K, Tanaka H, Kasai K, Irisawa T, Shiozaki T, Sugimoto H: Profile of gene expression in the subventricular zone after traumatic brain injury. J Neurotrauma 2003, 20: 1147-1162. 10.1089/089771503770802844.

    Article  PubMed  Google Scholar 

  30. Kamm K, Vanderkolk W, Lawrence C, Jonker M, Davis AT: The effect of traumatic brain injury upon the concentration and expression of interleukin-1beta and interleukin-10 in the rat. J Trauma 2006, 60: 152-157.

    Article  PubMed  Google Scholar 

  31. Truettner JS, Suzuki T, Dietrich WD: The effect of therapeutic hypothermia on the expression of inflammatory response genes following moderate traumatic brain injury in the rat. Brain Res Mol Brain Res 2005, 138: 124-134. 10.1016/j.molbrainres.2005.04.006.

    Article  CAS  PubMed  Google Scholar 

  32. Nicole O, Goldshmidt A, Hamill CE, Sorensen SD, Sastre A, Lyuboslavsky P, Hepler JR, McKeon RJ, Traynelis SF: Activation of protease-activated receptor-1 triggers astrogliosis after brain injury. J Neurosci 2005, 25: 4319-4329. 10.1523/JNEUROSCI.5200-04.2005.

    Article  CAS  PubMed  Google Scholar 

  33. Vitarbo EA, Chatzipanteli K, Kinoshita K, Truettner JS, Alonso OF, Dietrich WD: Tumor necrosis factor alpha expression and protein levels after fluid percussion injury in rats: the effect of injury severity and brain temperature. Neurosurgery 2004, 55: 416-24; discussion 424-5. 10.1227/01.NEU.0000130036.52521.2C.

    Article  PubMed  Google Scholar 

  34. Long Y, Zou L, Liu H, Lu H, Yuan X, Robertson CS, Yang K: Altered expression of randomly selected genes in mouse hippocampus after traumatic brain injury. J Neurosci Res 2003, 71: 710-720. 10.1002/jnr.10524.

    Article  CAS  PubMed  Google Scholar 

  35. Tehranian R, Andell-Jonsson S, Beni SM, Yatsiv I, Shohami E, Bartfai T, Lundkvist J, Iverfeldt K: Improved recovery and delayed cytokine induction after closed head injury in mice with central overexpression of the secreted isoform of the interleukin-1 receptor antagonist. J Neurotrauma 2002, 19: 939-951. 10.1089/089771502320317096.

    Article  PubMed  Google Scholar 

  36. Kinoshita K, Chatzipanteli K, Vitarbo E, Truettner JS, Alonso OF, Dietrich WD: Interleukin-1beta messenger ribonucleic acid and protein levels after fluid-percussion brain injury in rats: importance of injury severity and brain temperature. Neurosurgery 2002, 51: 195-203; discussion 203. 10.1097/00006123-200207000-00027.

    Article  PubMed  Google Scholar 

  37. Sang Q, Kim MH, Kumar S, Bye N, Morganti-Kossman MC, Gunnersen J, Fuller S, Howitt J, Hyde L, Beissbarth T, Scott HS, Silke J, Tan SS: Nedd4-WW domain-binding protein 5 (Ndfip1) is associated with neuronal survival after acute cortical brain injury. J Neurosci 2006, 26: 7234-7244. 10.1523/JNEUROSCI.1398-06.2006.

    Article  CAS  PubMed  Google Scholar 

  38. Yao XL, Liu J, Lee E, Ling GS, McCabe JT: Progesterone differentially regulates pro- and anti-apoptotic gene expression in cerebral cortex following traumatic brain injury in rats. J Neurotrauma 2005, 22: 656-668. 10.1089/neu.2005.22.656.

    Article  PubMed  Google Scholar 

  39. He J, Evans CO, Hoffman SW, Oyesiku NM, Stein DG: Progesterone and allopregnanolone reduce inflammatory cytokines after traumatic brain injury. Exp Neurol 2004, 189: 404-412. 10.1016/j.expneurol.2004.06.008.

    Article  CAS  PubMed  Google Scholar 

  40. Burbach GJ, Dehn D, Del Turco D, Staufenbiel M, Deller T: Laser microdissection reveals regional and cellular differences in GFAP mRNA upregulation following brain injury, axonal denervation, and amyloid plaque deposition. Glia 2004, 48: 76-84. 10.1002/glia.20057.

    Article  PubMed  Google Scholar 

  41. Raghavendra Rao VL, Dhodda VK, Song G, Bowen KK, Dempsey RJ: Traumatic brain injury-induced acute gene expression changes in rat cerebral cortex identified by GeneChip analysis. J Neurosci Res 2003, 71: 208-219. 10.1002/jnr.10486.

    Article  PubMed  Google Scholar 

  42. Rhinn H, Scherman D, Escriou V: One-step quantification of single-stranded DNA in the presence of RNA using Oligreen in a real-time polymerase chain reaction thermocycler. Anal Biochem 2008, 372: 116-118. 10.1016/j.ab.2007.08.023.

    Article  CAS  PubMed  Google Scholar 

  43. Rutledge RG: Sigmoidal curve-fitting redefines quantitative real-time PCR with the prospective of developing automated high-throughput applications. Nucleic Acids Res 2004, 32: e178. 10.1093/nar/gnh177.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Hervé Rhinn, Catherine Marchand-Leroux, Nicole Croci, Michel Plotkine, Daniel Scherman or Virginie Escriou.

Additional information

Authors' contributions

HR performed all the in vitro experimental procedures and the data analysis, and was primary author of the manuscript. VE and DS supervised the study design and contributed to writing the manuscript. MP and C M-L provided the experimental model of traumatic brain injury. NC performed the experimental trauma.

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Rhinn, H., Marchand-Leroux, C., Croci, N. et al. Housekeeping while brain's storming Validation of normalizing factors for gene expression studies in a murine model of traumatic brain injury. BMC Molecular Biol 9, 62 (2008). https://doi.org/10.1186/1471-2199-9-62

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2199-9-62

Keywords