- Research article
- Open Access
Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies
© Spinsanti et al; licensee BioMed Central Ltd. 2006
- Received: 06 June 2006
- Accepted: 19 September 2006
- Published: 19 September 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.
- Ribosomal Protein Gene
- Striped Dolphin
- Succinate Dehydrogenase Complex
- Optimal Control Gene
- BestKeeper Index
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.
GenBank Accession Numbers of partially sequenced genes in Stenella coeruleoalba.
Cytoskeletal structural protein
Cytoskeletal protein involved in cell locomotion
hypoxantine phosphoribosyltransferase 1
Metabolic salvage of purines in mammals
phosphoglycerate kinase 1
ribosomal Protein L4
Member of ribosome proteins
ribosomal Protein S18
Member of ribosome proteins
succinate dehydrogenase complex
Electron transporter in the TCA cycle and respiratory chain
FP, SDH2, SDHF
TFR; p90; CD71; TFR1; Trfr; Mtvr1; Mtvr-1
Signal transduction by binding to phosphorilated serine residues on a variety of signaling molecules
Details of primers and amplicons for each of the 10 evaluated genes.
Forward Primer Sequence [5'→3']
Position in cDNA
Reverse Primer Sequence [5'→3']
Position in cDNA
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
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.
Average M value
Candidate reference genes for normalization of qRT-PCR listed according to their expression stability calculated by the NormFinder VBA applet.
Results from BestKeeper descriptive statistical analysis.
SD [± Ct]
CV [% Ct]
Results from BestKeeper correlation analysis.
coeff. of corr. [r]
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.
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.
- 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-research0034PubMed CentralView ArticlePubMedGoogle Scholar
- 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-8View ArticlePubMedGoogle Scholar
- 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.6364190View ArticlePubMedGoogle Scholar
- 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.177View ArticlePubMedGoogle Scholar
- Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation. Biotechniques. 2000, 29 (2): 332-337.PubMedGoogle Scholar
- 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.47View ArticlePubMedGoogle Scholar
- 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.022View ArticlePubMedGoogle Scholar
- 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.081View ArticlePubMedGoogle Scholar
- 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.0376View ArticlePubMedGoogle Scholar
- 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.21215View ArticlePubMedGoogle Scholar
- 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-zView ArticlePubMedGoogle Scholar
- 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. 2006Google Scholar
- 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.20333View ArticlePubMedGoogle Scholar
- 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.xView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.009View ArticlePubMedGoogle Scholar
- 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/eri285View ArticlePubMedGoogle Scholar
- 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-3PubMed CentralView ArticlePubMedGoogle Scholar
- 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-21PubMed CentralView ArticlePubMedGoogle Scholar
- 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-0496View ArticlePubMedGoogle Scholar
- 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-14PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.10075View ArticlePubMedGoogle Scholar
- 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-3View ArticlePubMedGoogle Scholar
- 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:1026298032009View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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-9View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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-6View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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-7View ArticlePubMedGoogle Scholar
- geNorm Software [http://medgen.ugent.be/~jvdesomp/genorm/].Google Scholar
- NormFinder Software [http://www.mdl.dk/publicationsnormfinder.htm].Google Scholar
- BestKeeper Software [http://www.gene-quantification.de/bestkeeper.html].Google Scholar
- Winfrey MR, Rott MA, Wortman AT: Unraveling DNA: Molecular Biology for the Laboratory. 1997, Upper Saddle River, NJ , Prentice-HallGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- ClustalW [http://www.ebi.ac.uk/clustalw/].Google Scholar
- BioEdit. Biological Sequence Alignment Editor for Windows 95/98/NT/XP. [http://www.mbio.ncsu.edu/BioEdit/bioedit.html].Google Scholar
- Basic Local Alignment Search Tool (BLAST). [http://www.ncbi.nlm.nih.gov/BLAST/].Google Scholar
- 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.1232297100PubMed CentralView ArticlePubMedGoogle Scholar
- 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.118View ArticlePubMedGoogle Scholar
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.