Skip to main content


Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies

Article metrics



Odontocete cetaceans occupy the top position of the marine food-web and are particularly sensitive to the bioaccumulation of lipophilic contaminants. The effects of environmental pollution on these species are highly debated and various ecotoxicological studies have addressed the impact of xenobiotic compounds on marine mammals, raising conservational concerns. Despite its sensitivity, quantitative real-time PCR (qRT-PCR) has never been used to quantify gene induction caused by exposure of cetaceans to contaminants. A limitation for the application of qRT-PCR is the need for appropriate reference genes which allow the correct quantification of gene expression. A systematic evaluation of potential reference genes in cetacean skin biopsies is presented, in order to validate future qRT-PCR studies aiming at using the expression of selected genes as non-lethal biomarkers.


Ten commonly used housekeeping genes (HKGs) were partially sequenced in the striped dolphin (Stenella coeruleoalba) and, for each gene, PCR primer pairs were specifically designed and tested in qRT-PCR assays. The expression of these potential control genes was examined in 30 striped dolphin skin biopsy samples, obtained from specimens sampled in the north-western Mediterranean Sea. The stability of selected control genes was determined using three different specific VBA applets (geNorm, NormFinder and BestKeeper) which produce highly comparable results. Glyceraldehyde-3P-dehydrogenase (GAPDH) and tyrosine 3-monooxygenase (YWHAZ) always rank as the two most stably expressed HKGs according to the analysis with geNorm and Normfinder, and are defined as optimal control genes by BestKepeer. Ribosomal protein L4 (RPL4) and S18 (RPS18) also exhibit a remarkable stability of their expression levels. On the other hand, transferrin receptor (TFRC), phosphoglycerate kinase 1 (PGK1), hypoxanthine ribosyltransferase (HPRT1) and β-2-microglobin (B2M) show variable expression among the studied samples and appear as less suitable reference genes for data normalization.


In this work, we have provided essential background information for the selection of control genes in qRT-PCR studies of cetacean skin biopsies, as a molecular technique to investigate ecotoxicological hazard in marine mammals. Of 10 HKGs tested, those encoding for YWHAZ and GAPDH appear as the most reliable control genes for the normalization of qRT-PCR data in the analysis of striped dolphin skin biopsies. Potentially useful reference genes are also those encoding for ribosomal proteins L4 and S18.


Quantitative real-time PCR (qRT-PCR) represents a rapid and reliable method for the detection and quantification of mRNA transcription levels of a selected gene of interest (GOI). In ecotoxicological studies, qRT-PCR may be effectively used to evaluate the levels of expression of biomarker genes under the induction of xenobiotic contaminants. The sensitivity of qRT-PCR allows to work with a minimal amount of starting material, still achieving an accurate quantification of poorly transcribed mRNAs. Problems associated with the use of this assay are linked to the variability associated with the various steps of the experimental procedure, and could lead to severe misinterpretation of the results: different amounts and quality of starting material, variable enzymatic efficiencies (i.e. efficiency of retrotranscription from RNA to cDNA, and PCR efficiency) or differences between tissues or cells in overall transcriptional activity [1, 2]. Among several strategies proposed [2, 3], house-keeping genes (HKGs) are commonly accepted and frequently used to normalize qRT-PCR and to reduce possible errors generated in the quantification of gene expression. In this normalization strategy, internal controls are subjected to the same conditions as the RNA of interest and their expression measured by qRT-PCR. The success of this procedure is highly dependent on the choice of the appropriate control genes. Although many studies using qRT-PCR have relied upon only one endogenous control [4, 5], to date the use of a single HKG appears to be insufficient, and normalization by multiple controls is required [1, 6]. An ideal HKG, exposed to the same experimental protocol of the gene of interest (GOI), should present stable expression levels. If the expression of the reference gene is altered by the experimental conditions or by external factors, such as contamination, and is affected by a large variation, the noise of the assay is increased and detection of small changes becomes unfeasible, producing results that may be entirely incorrect [7]. Several works [810] prove how some of the most commonly used HKGs can not always be considered as reliable controls and/or they show different behaviour in various tissues [11], emphasizing the importance of preliminary evaluation studies, aimed at identifying the most stable HKGs in different organisms. The number of articles concerning the evaluation and selection of the best HKGs for each single experiment [1119] is rapidly increasing, together with the number of softwares which use statistical methods to evaluate stability of selected HKGs [1, 6, 20]. Some HKGs (those encoding for Act-B, GAPDH, HPRT1 and 18S ribosomal RNA) have been used as reference for many years in Northern blots, RNase protection tests and conventional quantitative PCR (qPCR) [3]. With the introduction of qRT-PCR, other ordinarily used HKGs, involved with basic and ubiquitous cellular functions and belonging to different functional classes, have been introduced [2125].

Top predators, such as odontocete cetaceans, are known to accumulate high concentrations of persistent organic pollutants (POPs), including endocrine disrupting chemicals (EDCs), thereby incurring in high toxicological hazard [26]. EDCs mime sex steroid hormones, both androgens and estrogens, binding to hormone receptors and influencing cellular pathways [27]. Xenobiotic compounds exhibit lipophilic properties and tend to accumulate preferentially in the fat tissue of top predators. Striped dolphins (Stenella coeruleoalba) inhabiting the Mediterranean basin appear to be particularly exposed to the bioaccumulation of OCs if compared to other worldwide distributed cetacean populations [27]. In addition, a geographical contamination trend was accurately described in the Mediterranean basin for the striped dolphin [28], where levels of OCs in the body of these animals decrease from north-west to south-east. It is evident that cetaceans inhabiting the Ligurian Sea appear to be much more exposed to OCs if compared to those sampled in the Eolian Archipelago and the Ionian Sea [27]. Specimens of the striped dolphin from the Ligurian Sea have been chosen as a test population to validate qRT-PCR studies aimed at a better understanding of the role of POPs in the induction of potential biomarker genes.

Previous ecotoxicological studies have assessed that cetacean skin biopsies constitute a suitable, non-lethal, biological material for the assessment of toxicological hazard in Mediterranean marine mammals [29, 30]. They enable the evaluation of levels of OCs with endocrine disrupting capacity and of Cytochrome P450 (CYP1A1) induction, as an early warning sign of exposure to organochlorines [2632].

Despite its sensitivity and reliability, qRT-PCR has never been used in cetacean skin biopsies to investigate gene induction. In this context, our work constitutes the first effort towards the routine application of this molecular technique for the evaluation of toxicological hazard in cetacean species. Ten among the most common HKGs have been described and ranked, according to the stability of their expression, in skin biopsies of the striped dolphin. The animals were sampled in the "Pelagos Sanctuary for Mediterranean Marine Mammals", a special marine protected area extending for about 90.000 km2 in the north-western Mediterranean Sea between Italy, France and the Sardinia Island, encompassing Corsica Island and the Tuscan Archipelago.

Results and discussion

We obtained the partial sequence of ten of the most commonly used housekeeping genes: β-Actin (Act-B), glyceraldehyde-3P-dehydrogenase (GAPDH), hypoxanthine ribosyltransferase (HPRT1), β-2-microglobin (B2M), phosphoglycerate kinase 1 (PGK1), succinate dehydrogenase complex, subunit A (SDHA), transferrin receptor (TFRC), tyrosine 3-monooxygenase (YWHAZ), L4 (RPL4) and S18 (RPS18) ribosomal proteins. Sequence information, deposited in GenBank under the Accession Numbers shown in Table 1, was used to design RT-PCR primers (Table 2). Target amplicons, never exceeding 171 bp, were sequenced to confirm specificity of the products. The efficiency of real-time amplifications for each single primer pair was determined in RT-PCR using 1:5 cDNA serial dilutions retrotranscribed from 1 μg total RNA isolated in striped dolphin skin biopsies (Table 2). Expression patterns of the ten HKGs have been subsequently evaluated in skin biopsies (n = 30) obtained using non-lethal techniques, in a Mediterranean sub-population of the striped dolphin from the Ligurian Sea.

Table 1 GenBank Accession Numbers of partially sequenced genes in Stenella coeruleoalba.
Table 2 Details of primers and amplicons for each of the 10 evaluated genes.

RNA quality

All RNA samples were examined as to their concentration, purity and integrity. Based on absorbance ratio at 260/280 nm and at 230/260 nm, all samples were pure, free from protein and organic pollutants derived from RNA extraction. Overall sample integrity was confirmed by denaturing formaldehyde agarose gel electrophoresis, showing sharp and intense 18S and 28S ribosomal RNA bands with a total absence of smears.

Expression levels of candidate reference genes

The 10 studied HKGs displayed a relatively wide range of expression level, from the lowest mean Ct values (16.1) in RPL4, to the highest (24.6) in SDHA. We grouped the expression levels in two arbitrary categories as shown in Figure 1. Highly expressed genes have mean Ct values below 20 cycles (RPL4, RPS18, GAPDH, Act-B and YWHAZ), and lowly expressed genes have mean Ct values above 20 cycles (PGK1, B2M, HPRT1, TFRC and SDHA). Individual HKGs have different expression levels across all studied samples, as shown in Figure 1. RPL4 and GAPDH have the smallest gene expression variation (below 3 cycles) while, PGK1, Act-B, SDHA and B2M are the genes with the most variable levels of expression (over 4 cycles).

Figure 1

Expression levels of candidate HKGs. Expression levels of candidate control genes in the tested striped dolphin skin biopsies (n = 30). Values are given as qRT-PCR cycle threshold numbers (Ct values). Boxes represent mean Ct values and whiskers the range of Ct values in the 30 samples. The arbitrary line at Ct 20 distinguishes two groups of similarly expressed housekeeping genes.

Data analysis

The data obtained for each skin biopsy and each HKG were analysed using three different VBA applets, geNorm [1, 33], NormFinder [20, 34] and BestKeeper [6, 35].

geNorm provides a ranking of the tested genes, based on their expression stability, determining the two most stable HKGs or a combination of multiple stable genes for normalization. Selected HKGs were ranked according to the determined control gene-stability measure (M, average pair-wise variation of a particular gene with all other control genes), from the most stable (lowest M values) to the least stable (highest M values): GAPDH/YWHAZ – RPS18 – RPL4 – SDHA – Act-B – TFRC – PGK1 – HPRT1 – B2M (Table 3; Figure 2a). All studied genes reach a high expression stability with low M values, less than 1, below the default limit of M = 1.5 [1]. Interestingly, the expression of the gene encoding for Act-B, widely used as a HKG in many studies, appears to be less stable than other genes in skin samples where Act-B should be constitutively expressed at high levels, presenting a large expression variation as calculated by Ct values. It is also remarkable that the two ribosomal proteins (S18 and L4) exhibit a similar behaviour (ranking as the 3rd and 4th most stably expressed genes, respectively) and are characterized by highly comparable average M values (Table 3). Additionally, the assessment of the normalization factor allows the identification of the optimal number of control genes. The geNorm software suggests that an accurate normalization factor of qRT-PCR data can be calculated by using these two most stably expressed genes (YWHAZ and GAPDH). As shown in Figure 2b, the addition of further HKGs will not significantly affect the reliability of the determined normalization factor, yielding a V2/3 value (pair-wise variation between two sequential normalization factors) of 0.141, lower than the default cut-off value of 0.15 [1]. According to the geNorm stability rank of the HKGs studied, the third gene to include in the calculation of a reliable normalization factor should be chosen between the two ribosomal protein genes S18 and L4 (Table 3; Figure 2a).

Table 3 Candidate reference genes for normalization of qRT-PCR ranked according to their expression stability (calculated as the average M value after stepwise exclusion of worst scoring genes) by the geNorm VBA applet.
Figure 2

geNorm output charts. (a) average expression stability measure (M) of control genes during stepwise exclusion of the least stable control genes; (b) determination of the optimal number of control genes for normalization calculated on the basis of the pair-wise variation (V) analysis; V values under 0.15 threshold line indicate no need to include further HKG for calculation of a reliable normalization factor.

NormFinder is another VBA applet based on an algorithm for identifying the optimal normalization gene(s) among a set of candidates. It ranks the set of candidate genes according to their expression stability value in a given sample set and a given experimental design [20]. The results of the NormFinder analysis applied to our data are shown in Table 4. This ranking appears to be identical to the one previously determined using geNorm, except for the position of the two ribosomal protein genes S18 and L4, which are inverted. YWHAZ, GAPDH, RPL4 and RPS18 still occupy the highest positions, while TFRC, PGK1, HPRT1 and B2M are equally defined as the least reliable controls. As one would expect given their identical cellular function, the genes encoding for the ribosomal proteins S18 and L4 occupy the same position as in the geNorm rank.

Table 4 Candidate reference genes for normalization of qRT-PCR listed according to their expression stability calculated by the NormFinder VBA applet.

BestKeeper is an Excel-based tool determining the "optimal" HKGs by using a pair-wise correlation analysis of all pairs of candidate genes, and calculating the geometric mean of the "best" suited ones [6]. The first analysis of the data, based on the inspection of raw Ct values calculated variations for all the HKGs in the 30 biopsies samples (SD [± Ct] and CV [%Ct]), shows an overall stability in gene expression (Table 5). Since all tested HKGs exhibit a SD value lower than 1, none of them can be clearly considered inconsistent, and therefore excluded from the calculation of the BestKeeper index [6]. For this reason, all control genes have been retained in the calculation of the BestKeeper index, which finally exhibit a moderate SD variation of 0.51. Most important, BestKeeper allows a comparative analysis across HKGs, by estimating correlations in the expression levels between all the possible candidates. Highly correlated control genes are combined into an index. Afterwards, the pair-wise correlation between genes and the correlation between each gene and the index are calculated, describing the relation between the index and the contributing HKG [6]. The 10 control genes tested in our analysis correlate well one with another, and also if compared with the BestKeeper index (Table 6). The best correlation between the HKGs and the BestKeeper index is obtained for YWHAZ (r = 0.938), followed by Act-B and GAPDH (Table 6). The statistically significant correlation shown by Act-B (r = 0.935) with the BestKeeper index (Table 6) appears inconsistent with the performance of this gene as assessed by geNorm and NormFinder. YWHAZ and GAPDH are still ranked as two of the most reliable control genes, showing respectively the best (0.938) and the third best (0.927) correlation values (Table 6). It is also remarkable how TFRC, PGK1, HPRT1 and B2M are again classified as the least reliable HKGs, showing the worst correlations with the determined BestKeeper index (Table 6). As previously mentioned, the genes encoding for the ribosomal proteins L4 and S18, which share identical cellular functions, appear to be strictly correlated and show highly comparable SD and r values (Table 5 and 6).

Table 5 Results from BestKeeper descriptive statistical analysis.
Table 6 Results from BestKeeper correlation analysis.

Calculation of the InVar parameter (as described in [6]) shows an overall good sample integrity.


This work constitutes the first effort for the selection of optimal control genes in qRT-PCR studies designed for the assessment of ecotoxicological hazard in cetacean populations. The three softwares tested (geNorm, NormFinder and BestKeeper), based on different algorithms and analytical procedures, produced highly comparable results.

Our study let us conclude that GAPDH and YWHAZ are the most reliable HKGs of this set and we therefore strongly recommend their use in future qRT-PCR studies on the striped dolphin. However, it is common opinion that the reliability of the normalization factor increases together with the number of stably expressed HKGs included in its calculation and the use of at least three HKGs for the correct normalization of qRT-PCR data has been proposed [1]. Our analysis suggests that the ribosomal protein genes RPL4 and RPS18 constitute additional useful HKGs.

On the other hand, TFRC, PGK1, HPRT1 and B2M show unstable expression patterns and are always classified as the least reliable control genes of this group.

It appears, from Figure 1, that the best reference genes selected from this study are also those with higher levels of expression (lower Ct values). While highly expressed genes might in principle exhibit lower levels of inter-individual variation, it should be pointed out that the softwares tested here do not assess stability as a function of variation in Ct values. In fact, the reliability of a given HKG is measured as a function of the variation of the ratios between the level of expression of pairs of HKGs across all samples. It is, therefore, a relative measure (to all other HKGs tested), not an absolute measure of variation of expression of a single control gene across samples. The gene encoding for Act-B is indeed one of the most commonly used HKGs, and it has found application in expression studies of various animal taxa. In our study, however, this gene is outperformed by several other genes (mostly GAPDH, YWHAZ, RPL4 and RPS18) as an optimal HKG. It is worth noting that Act-B shows a significant correlation (Table 6) with the BestKeeper index, comparable to that of GAPDH and YWHAZ, an outcome that would suggest some potential usefulness for this gene.

Ideally, in an experimental application, one could choose to use preferentially HKGs that have the same expression levels as the target gene, in order to enhance the uniformity of the analysis. However, two considerations must be made in this context. First, the use of target genes with extremely low physiological levels of expression should always be avoided, because of the considerable variance due to experimental fluctuations. Second, by definition, the expression levels of a good target gene are supposed to be under a strong influence by any contaminant whose effects are being studied. It is therefore likely that these levels vary across a wide range of Ct values. Nevertheless, if a target gene with very low level of expression is chosen, we may recommend the inclusion, among the control genes, of SDHA, which has demonstrated to possess a sub-optimal suitability as a reference gene.


Skin biopsies sampling

Sub-samples of skin biopsies (around 30 mg in weight) of 30 different specimens of the striped dolphin were obtained during several sampling trips in the Ligurian Sea, in the summer of 2005 (International CITES permit: IT 007). Dolphins were sampled from the prow of the boat using biopsy tips mounted on an aluminium pole [29]. Biopsy samples were immediately stored in RNA later (Ambion) at 4°C overnight and, subsequently, at -80°C.

Total RNA and genomic DNA isolation and characterization

Skin biopsies were homogenized with a TissueLyser (Qiagen, Hilden, Germany).

Total RNA was subsequently extracted using the RNeasy Fibrous Tissue Mini Kit (Qiagen), according to the manufacturer's instructions. Genomic DNA was eliminated by a DNase-on-column treatment supplied with the kit. The quantity of extracted RNA was assessed with the NanoDrop® ND-1000 UV-Vis Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA); absorbance ratio at 260/280 nm and 260/230 nm was used to assess purity of the RNA samples [36]. The integrity of the samples was determined by denaturing agarose gel (1.2%) electrophoresis and ethidium bromide staining.

Genomic DNA was isolated from homogenized skin biopsies using the Wizard® SV Genomic DNA Purification System (Promega, Madison, Wisconsin, USA). Final elution was in 80 μl Nuclease Free Water (Promega) and 1 μl was subsequently used in PCR amplification reactions.

Sequencing HKGs

Ten potential HKGs have been evaluated in this study: β-Actin (Act-B), glyceraldehyde-3P-dehydrogenase (GAPDH), hypoxanthine ribosyltransferase (HPRT1), β-2-microglobin (B2M), phosphoglycerate kinase 1 (PGK1), succinate dehydrogenase complex, subunit A (SDHA), transferrin receptor (TFRC), tyrosine 3-monooxygenase (YWHAZ), and the ribosomal proteins L4 (RPL4) and S18 (RPS18). To reduce possible errors due to co-regulation of different genes, we decided to include in our analysis HKGs involved in different cellular functions (Table 1). The only potential co-regulated genes are those encoding for the two ribosomal proteins S18 and L4. Partial nucleotide sequences of all selected HKGs have been determined in the striped dolphin. PCR primers were designed in conserved regions by comparing sequences of all genes available from GenBank for several mammal species. Alignments were performed with ClustalW [37, 38] and BioEdit [39]. Both alignments and primer sequences are available from the authors upon request. Primers were used in PCR reactions on cDNA retrotranscribed from striped dolphin skin biopsies. Retrotranscriptions were carried out with 1 μg total RNA using iScript cDNA synthesis kit (Bio-Rad, Hercules, California, USA), according to the manufacturer's instructions. 1 μl cDNA was used as template in a standard 25 μl PCR reaction, using 10 μM of each primer and 0.05 U Taq polimerase (Promega). In order to extend the sequence information for some genes, 3'-RACE experiments were carried out for Act-B, RPS18 and GAPDH genes, using the 3'-RACE System for Rapid Amplification of cDNA Ends (Invitrogen, Carlsbad, California, USA). A 5'-RACE experiment was also performed with the 5'-RACE System (Invitrogen) on the cDNA of the B2M gene.

Amplification products were first loaded and checked on a 1% agarose gel, and sequenced on a automatic sequencer CEQ8000XL (Beckman Coulter, Fullerton, California, USA) using DTCS system, purified with Wizard SV Gel and PCR Clean-Up System (Promega). Electropherograms and sequences were compared and manually corrected using Sequencher 4.5 (GeneCodes, Ann Arbor, Michigan, USA). Sequences, checked for their specificity using Blast [40], were deposited in GenBank (see Table 1 for Accession Numbers).

Design of RT-PCR primers

PCR primers for real-time assays were designed on the determined nucleotide sequences of each HKG using Beacon Designer 2.06 (Premier Biosoft International, Palo Alto, California, USA). Special attention was given to primer length, annealing temperature, base composition and 3'-end stability. To ensure optimal DNA polymerization efficiency, amplicons length ranged between 84 and 171 bp. Exon and intron boundaries were determined aligning sequenced genes with corresponding human genomic sequences downloaded from GenBank. Position of introns in the HKGs, deduced from the alignments with Homo sapiens, was subsequently confirmed by PCR on striped dolphin genomic DNA. As expected [41, 42], intron position was highly conserved between cetaceans and human. On the basis of this knowledge, except for GAPDH, we designed RT-PCR primers located on different exons or directly spanning exon-exon junction of each cDNA (Table 2). After preliminary real-time assays, the optimal primer concentration appeared to be 300 mM, which generates the lowest Ct value and a sharp peak, with no amplification of non-specific products or primer-dimer artifacts (see additional files 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10). For each different pair of primers, efficiency of RT-PCR (E), slope values, and correlation coefficients (R2) were determined, using serial 1:5 dilutions of template cDNA, on a iQ5 machine (Bio-Rad) (Table 2). Products were subsequently run on 2% agarose gel to check for size specificity, and, eventually, sequenced.

Real-time quantitative PCR (qRT-PCR)

To avoid any genomic DNA contamination during qRT-PCR, retrotranscription of 1 μg total RNA was performed with the QuantiTect Reverse Transcription Kit (Qiagen), according to the manufacturer's instructions. This kit ensures complete digestion of genomic DNA by a brief incubation of the sample at 42°C with a specific Wipeout buffer before retrotranscription. For each of the thirty samples, 1 μl of treated RNA was collected before proceeding with the retrotranscription step, and stored at -80°C. Treated RNA samples have been subsequently used in RT-PCR, as internal controls to check for potential genomic DNA contamination.

The second part of this two-step RT-PCR protocol was characterized by real-time amplifications, using SYBR Green detection chemistry, run in triplicate on 96-wells reaction plates with the iQ5 machine (Bio-Rad). Reactions were prepared in a total volume of 20 μl containing: 1 μl cDNA, 0.8 μl of each 10 μM primer (300 mM each; Invitrogen), 10 μl of iQ™ SYBR® Green Supermix (Bio-Rad) and 8 μl RNase,/DNase-free sterile water (Qiagen). Blank controls were run in triplicate for each master mix. The cycle conditions were set as follows: initial template denaturation at 95°C for 1 min, followed by 40 cycles of denaturation at 95°C for 10 s, and combined primer annealing/elongation at 60°C for 30 s, as described in [4]. This cycle was followed by a melting curve analysis, ranging from 56°C to 95°C, with temperature increasing steps of 0.5°C every 10 s.

Baseline and threshold values were automatically determined for all plates using the Bio-Rad iQ5 Software 1.0. In order to ensure comparability between data obtained from different experimental plates, the threshold value has been subsequently manually set to the value corresponding to the arithmetic mean between the automatically determined thresholds annotated previously; then all data have been reanalyzed. Raw Ct values were transformed to quantities using an Excel spreadsheet generated by the authors, based on the comparative Ct method. The data obtained have been converted into correct input files, according to the requirements of the software, and analysed using geNorm (version 3.4), NormFinder (version 0.953) and BestKeeper (version 1)VBA applets.


  1. 1.

    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 (7): RESEARCH0034- 10.1186/gb-2002-3-7-research0034

  2. 2.

    Ginzinger DG: Gene quantification using real-time quantitative PCR: an emerging technology hits the mainstream. Exp Hematol. 2002, 30 (6): 503-512. 10.1016/S0301-472X(02)00806-8

  3. 3.

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

  4. 4.

    Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004, 313 (4): 856-862. 10.1016/j.bbrc.2003.11.177

  5. 5.

    Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29 (2): 332-337.

  6. 6.

    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 (6): 509-515. 10.1023/B:BILE.0000019559.84305.47

  7. 7.

    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 (1): 141-143. 10.1016/j.ab.2005.05.022

  8. 8.

    Bemeur C, Ste-Marie L, Desjardins P, Hazell AS, Vachon L, Butterworth R, Montgomery J: Decreased beta-actin mRNA expression in hyperglycemic focal cerebral ischemia in the rat. Neurosci Lett. 2004, 357 (3): 211-214. 10.1016/j.neulet.2003.12.081

  9. 9.

    Selvey S, Thompson EW, Matthaei K, Lea RA, Irving MG, Griffiths LR: Beta-actin--an unsuitable internal control for RT-PCR. Mol Cell Probes. 2001, 15 (5): 307-311. 10.1006/mcpr.2001.0376

  10. 10.

    Lee JH, Fitzgerald JB, Dimicco MA, Grodzinsky AJ: Mechanical injury of cartilage explants causes specific time-dependent changes in chondrocyte gene expression. Arthritis Rheum. 2005, 52 (8): 2386-2395. 10.1002/art.21215

  11. 11.

    Ohl F, Jung M, Xu C, Stephan C, Rabien A, Burkhardt M, Nitsche A, Kristiansen G, Loening SA, Radonic A, Jung K: Gene expression studies in prostate cancer tissue: which reference gene should be selected for normalization?. J Mol Med. 2005, 83 (12): 1014-1024. 10.1007/s00109-005-0703-z

  12. 12.

    Dydensborg AB, Herring E, Auclair J, Tremblay E, Beaulieu JF: Normalizing Genes for Quantitative RT-PCR in Differentiating Human Intestinal Epithelial Cells and Adenocarcinomas of the Colon. Am J Physiol Gastrointest Liver Physiol. 2006

  13. 13.

    Bettegowda A, Patel OV, Ireland JJ, Smith GW: Quantitative analysis of messenger RNA abundance for ribosomal protein L-15, cyclophilin-A, phosphoglycerokinase, beta-glucuronidase, glyceraldehyde 3-phosphate dehydrogenase, beta-actin, and histone H2A during bovine oocyte maturation and early embryogenesis in vitro. Mol Reprod Dev. 2006, 73 (3): 267-278. 10.1002/mrd.20333

  14. 14.

    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 (4): 770-777. 10.1111/j.1742-4658.2006.05109.x

  15. 15.

    Chen J, Rider DA, Ruan R: Identification of valid housekeeping genes and antioxidant enzyme gene expression change in the aging rat liver. J Gerontol A Biol Sci Med Sci. 2006, 61 (1): 20-27.

  16. 16.

    Meller M, Vadachkoria S, Luthy DA, Williams MA: Evaluation of housekeeping genes in placental comparative expression studies. Placenta. 2005, 26 (8-9): 601-607. 10.1016/j.placenta.2004.09.009

  17. 17.

    Nicot N, Hausman JF, Hoffmann L, Evers D: Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot. 2005, 56 (421): 2907-2914. 10.1093/jxb/eri285

  18. 18.

    Garcia-Crespo D, Juste RA, Hurtado A: Selection of ovine housekeeping genes for normalisation by real-time RT-PCR; analysis of PrP gene expression and genetic susceptibility to scrapie. BMC Vet Res. 2005, 1: 3- 10.1186/1746-6148-1-3

  19. 19.

    Olsvik PA, Lie KK, Jordal AE, Nilsen TO, Hordvik I: Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol Biol. 2005, 6: 21- 10.1186/1471-2199-6-21

  20. 20.

    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 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496

  21. 21.

    Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol. 2004, 4: 14- 10.1186/1471-2229-4-14

  22. 22.

    Aerts JL, Gonzales MI, Topalian SL: Selection of appropriate control genes to assess expression of tumor antigens using real-time RT-PCR. Biotechniques. 2004, 36 (1): 84-6, 88, 90-1.

  23. 23.

    Zhang QJ, Chadderton A, Clark RL, Augustine-Rauch KA: Selection of normalizer genes in conducting relative gene expression analysis of embryos. Birth Defects Res A Clin Mol Teratol. 2003, 67 (8): 533-544. 10.1002/bdra.10075

  24. 24.

    Loseke S, Grage-Griebenow E, Wagner A, Gehlhar K, Bufe A: Differential expression of IFN-alpha subtypes in human PBMC: evaluation of novel real-time PCR assays. J Immunol Methods. 2003, 276 (1-2): 207-222. 10.1016/S0022-1759(03)00072-3

  25. 25.

    Kim BR, Nam HY, Kim SU, Kim SI, Chang YJ: Normalization of reverse transcription quantitative-PCR with housekeeping genes in rice. Biotechnol Lett. 2003, 25 (21): 1869-1872. 10.1023/A:1026298032009

  26. 26.

    Marsili L, Casini C, Marini L, Regoli A, Focardi S: Age, growth and organochlorines (HCB, DDTs and PCBs) in Mediterranean striped dolphins Stenella coeruleoalba stranded in 1988-1994 on the coasts of Italy. Marine Ecology Progress Series. 1997, 151: 273-282.

  27. 27.

    Fossi MC, Marsili L, Neri G, Natoli A, Politi E, Panigada S: The use of a non-lethal tool for evaluating toxicological hazard of organochlorine contaminants in Mediterranean cetaceans: new data 10 years after the first paper published in MPB. Mar Pollut Bull. 2003, 46 (8): 972-982. 10.1016/S0025-326X(03)00113-9

  28. 28.

    Marsili L: Lipophilic contaminants in marine mammals: review of the results of ten years work at the Department of Environmental Biology, Siena University (Italy). International Journal of the Environment and Pollution. 2000, 416-452. 10.1504/IJEP.2000.002329.

  29. 29.

    Fossi MC, Marsili L, Neri G, Casini S, Bearzi G, Politi E, Zanardelli M, Panigada S: Skin biopsy of Mediterranean cetaceans for the investigation of interspecies susceptibility to xenobiotic contaminants. Mar Environ Res. 2000, 50 (1-5): 517-521. 10.1016/S0141-1136(00)00127-6

  30. 30.

    Fossi MC, Marsili L, Leonzio C, Notarbartolo Di Sciara G, Zanardelli M, Focardi S: The use of non-destructive biomarker in Mediterranean cetaceans: Preliminary data on MFO activity in skin biopsy. Marine Pollution Bulletin. 1992, 24 (9): 459-461. 10.1016/0025-326X(92)90346-8.

  31. 31.

    Stegeman JJ: Cytochrome P450 gene diversity and function in marine animals: past, present, and future. Marine Environmental Research. 2000, 50 (1-5): 61-62. 10.1016/S0141-1136(00)00136-7.

  32. 32.

    Monaci F, Borrel A, Leonzio C, Marsili L, Calzada N: Trace elements in striped dolphins (Stenella coeruleoalba) from the western Mediterranean. Environ Pollut. 1998, 99 (1): 61-68. 10.1016/S0269-7491(97)00174-7

  33. 33.

    geNorm Software [].

  34. 34.

    NormFinder Software [].

  35. 35.

    BestKeeper Software [].

  36. 36.

    Winfrey MR, Rott MA, Wortman AT: Unraveling DNA: Molecular Biology for the Laboratory. 1997, Upper Saddle River, NJ , Prentice-Hall

  37. 37.

    Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680.

  38. 38.

    ClustalW [].

  39. 39.

    BioEdit. Biological Sequence Alignment Editor for Windows 95/98/NT/XP. [].

  40. 40.

    Basic Local Alignment Search Tool (BLAST). [].

  41. 41.

    Roy SW, Fedorov A, Gilbert W: Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain. Proc Natl Acad Sci U S A. 2003, 100 (12): 7158-7162. 10.1073/pnas.1232297100

  42. 42.

    Rogozin IB, Sverdlov AV, Babenko VN, Koonin EV: Analysis of evolution of exon-intron structure of eukaryotic genes. Brief Bioinform. 2005, 6 (2): 118-134. 10.1093/bib/6.2.118

Download references


This project is partially supported by grant from the Italian Ministry for the Environment (Convention 12/12/05 MATT).

We would like to thank Dr. Tommaso Mancuso and Dr. Marco Giorgio Bianchi (Bio-Rad System Specialists – Gene eXpr. Division) for excellent technical assistance. We also thank Dr. Jo Vandesompele, Dr. Michael W. Pfaffl, Dr. Claus Lindbjerg Andersen and Dr. Pål Olsvik for their valuable suggestions and comments.

Author information

Correspondence to Giacomo Spinsanti.

Additional information

Authors' contributions

GS performed all the experimental procedure and was primary author of the manuscript. CP and EL participated in the experimental process and data analysis. LM and MCF provided the samples. MCF, FF, LM and SC supervised the study design and contributed to writing the manuscript.

Electronic supplementary material

Additional File 1: Melting curve analyses obtained for the Act-B gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the Act-B gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 161 KB)

Additional File 2: Melting curve analyses obtained for the GAPDH gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the GAPDH gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 182 KB)

Additional File 3: Melting curve analyses obtained for the HPRT1 gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the HPRT1 gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 165 KB)

Additional File 4: Melting curve analyses obtained for the B2M gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the B2M gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 183 KB)

Additional File 5: Melting curve analyses obtained for the PGK1 gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the PGK1 gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 170 KB)

Additional File 6: Melting curve analyses obtained for the SDHA gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the SDHA gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 174 KB)

Additional File 7: Melting curve analyses obtained for the TFRC gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the TFRC gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 174 KB)

Additional File 8: Melting curve analyses obtained for the YWHAZ gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the YWHAZ gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 169 KB)

Additional File 9: Melting curve analyses obtained for the RPL4 gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the RPL4 gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 168 KB)

Additional File 10: Melting curve analyses obtained for the RPS18 gene. Melting curve analyses image (jpg format) collected using the iQ5 Optical System Software 1.0 (Bio-Rad) during calibration experiments of the selected primer pair for the RPS18 gene. Data were obtained using 1:5 dilutions of template cDNA (retrotranscribed from striped dolphin skin biopsy isolated total RNA) on a iQ5 machine (Bio-Rad). (JPEG 191 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

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 (, 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

Spinsanti, G., Panti, C., Lazzeri, E. et al. Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies. BMC Molecular Biol 7, 32 (2006) doi:10.1186/1471-2199-7-32

Download citation


  • Ribosomal Protein Gene
  • Striped Dolphin
  • Succinate Dehydrogenase Complex
  • Optimal Control Gene
  • BestKeeper Index