- Research article
- Open Access
Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies
BMC Molecular Biology volume 7, Article number: 32 (2006)
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 . Several works [8–10] 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 , 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 [11–19] 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) . 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 [21–25].
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 . EDCs mime sex steroid hormones, both androgens and estrogens, binding to hormone receptors and influencing cellular pathways . 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 . In addition, a geographical contamination trend was accurately described in the Mediterranean basin for the striped dolphin , 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 . 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 [26–32].
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.
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).
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 . 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 . 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).
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 . 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.
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 . 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 . 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 . 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).
Calculation of the InVar parameter (as described in ) 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 . 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 . 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 . 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.
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 . 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 , 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 . 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.
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
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
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
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
Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29 (2): 332-337.
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
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
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
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
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
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
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
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
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
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.
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
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
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
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
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
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
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.
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
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
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
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.
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
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.
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
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.
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.
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
geNorm Software [http://medgen.ugent.be/~jvdesomp/genorm/].
NormFinder Software [http://www.mdl.dk/publicationsnormfinder.htm].
BestKeeper Software [http://www.gene-quantification.de/bestkeeper.html].
Winfrey MR, Rott MA, Wortman AT: Unraveling DNA: Molecular Biology for the Laboratory. 1997, Upper Saddle River, NJ , Prentice-Hall
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.
BioEdit. Biological Sequence Alignment Editor for Windows 95/98/NT/XP. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html].
Basic Local Alignment Search Tool (BLAST). [http://www.ncbi.nlm.nih.gov/BLAST/].
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
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
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.
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)
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). https://doi.org/10.1186/1471-2199-7-32
- Ribosomal Protein Gene
- Striped Dolphin
- Succinate Dehydrogenase Complex
- Optimal Control Gene
- BestKeeper Index