Selection and validation of a set of reliable reference genes for quantitative RT-PCR studies in the brain of the Cephalopod Mollusc Octopus vulgaris
- Maria Sirakov†1, 5,
- Ilaria Zarrella†1,
- Marco Borra2,
- Francesca Rizzo3,
- Elio Biffali2,
- Maria Ina Arnone3 and
- Graziano Fiorito1, 4Email author
© Sirakov et al; licensee BioMed Central Ltd. 2009
Received: 27 August 2008
Accepted: 14 July 2009
Published: 14 July 2009
Quantitative real-time polymerase chain reaction (RT-qPCR) is valuable for studying the molecular events underlying physiological and behavioral phenomena. Normalization of real-time PCR data is critical for a reliable mRNA quantification. Here we identify reference genes to be utilized in RT-qPCR experiments to normalize and monitor the expression of target genes in the brain of the cephalopod mollusc Octopus vulgaris, an invertebrate. Such an approach is novel for this taxon and of advantage in future experiments given the complexity of the behavioral repertoire of this species when compared with its relatively simple neural organization.
We chose 16S, and 18S rRNA, actB, EEF1A, tubA and ubi as candidate reference genes (housekeeping genes, HKG). The expression of 16S and 18S was highly variable and did not meet the requirements of candidate HKG. The expression of the other genes was almost stable and uniform among samples. We analyzed the expression of HKG into two different set of animals using tissues taken from the central nervous system (brain parts) and mantle (here considered as control tissue) by BestKeeper, geNorm and NormFinder. We found that HKG expressions differed considerably with respect to brain area and octopus samples in an HKG-specific manner. However, when the mantle is treated as control tissue and the entire central nervous system is considered, NormFinder revealed tubA and ubi as the most suitable HKG pair. These two genes were utilized to evaluate the relative expression of the genes FoxP, creb, dat and TH in O. vulgaris.
We analyzed the expression profiles of some genes here identified for O. vulgaris by applying RT-qPCR analysis for the first time in cephalopods. We validated candidate reference genes and found the expression of ubi and tubA to be the most appropriate to evaluate the expression of target genes in the brain of different octopuses. Our results also underline the importance of choosing a proper normalization strategy when analyzing gene expression by qPCR taking into appropriate account the experimental setting and variability of the sample of animals (and tissues), thus providing a set of HGK which expression appears to be unaffected by the experimental factor(s).
Relative quantification of mRNA transcripts requires endogenous normalizers, i.e. reference genes. The normalization of the expression of target genes favors the elimination of unspecific variation caused by differences in starting material, RNA extraction, enzymatic efficiencies, transcriptional activity and presence of inhibitors in the samples. A suitable reference gene has to i. be adequately expressed in the tissue of interest, ii. not be co-regulated, and iii. show comparable expression levels with target genes. In addition, it should show reduced variability in expression levels among samples and experimental conditions [1–4].
A brief review of the literature shows that a limited number of genes is commonly utilized in RT-qPCR experiments as reference. However, strategies have been developed to validate the use of reference genes experimentally [2, 4–10] and to apply computational approaches that allow to monitor the stability of genes in respect to others in all samples [11, 12], thus allowing proper normalization of RT-qPCR data.
The aim of the present study is to identify candidate genes in the cephalopod mollusc Octopus vulgaris that could be used in RT-qPCR experiments as internal reference (housekeeping genes, HKGs) to normalize the expression of target genes. Such an approach is novel for this taxon and may result of great advantage in future experiments given the complexity of the behavioral repertoire of this animal when compared with its relatively simple neural organization [13, 14].
The brain (supra- and sub-oesophageal masses and optic lobes ) and, as control tissue, the mantle of naïve O. vulgaris were utilized as samples in this study. Two different sets of individuals were used: a broad collection (for both size and age) caught throughout the year (set a) and a restricted sample (in body size) collected exclusively during the summer (set b).
Expression stability of candidate reference genes
Gene products, accession number (AN) and gene ontology accession number (GO) of genes considered in this work.
Primer sequence 5' – 3'
Forkhead box protein P
Ubiquitin/ribosomal protein S27a
The analysis of samples taken from a more homogeneous set of animals (set b, highlighted row) provided, as expected, a consensus between the three algorithms. In fact, each software indicated the pair tubA-ubi as the best combination of reference genes (ubi is the best gene according to NormFinder).
A different pattern appeared when both sets of samples were considered (Figure 2). Again, geNorm was not capable to provide a combination of reference genes because the variability in expression resulted in M values above the suggested level of acceptance. Similarly, BestKeeper indicated tubA as the best reference gene, but suggested that the levels of ubi mRNA in the different samples were too variable to be recommended as reference. NormFinder found ubi as the best HKG.
NormFinder also allows to define experimental conditions/samples in a data set . For our purposes, here we considered the mantle as control tissue when studying the expression of different genes in the brain. In these circumstances, NormFinder revealed ubi and tubA as the best combination of reference genes for evaluating the expression of target genes in animals belonging to different sets.
Expression level of target genes
Discussion and conclusion
In this study, we analyzed the expression profiles of some genes here identified for Octopus vulgaris for the first time. We applied RT-qPCR analysis, which is unprecedented in cephalopods, and, as a prerequisite, we also validated candidate reference genes. Normalization strategies are required to correct sample-to-sample variability in order to reveal gene-specific variation among experimental conditions and/or tissues. Similarly to what occurs in other animal models [e.g. [18, 19]], normalization of mRNA levels to the cell number is not possible for O. vulgaris tissues.
Our experiments showed that the expression of several so called housekeeping genes vary among different conditions (possibly body size or seasonality) and/or tissues (i.e. mantle vs brain), but also when the different parts in the brain of the octopus (masses) are considered. However, in such circumstances an algorithm that takes into account different conditions in a given experimental design (NormFinder ) results to be beneficial for the identification of suitable reference genes to be used in RT-qPCR experiments. This was our case. NormFinder allowed us to identify ubi and tubA as the appropriate reference genes for most of the analyzed samples. Moreover, this gene pair resulted to be chosen as the best combination when a tissue (mantle in our condition) is used as control to study the expression of different genes in the brain. We applied this reference gene pair to analyze the level of expression of our target genes (creb, dat, TH, FoxP) in different regions of Octopus brain applying a geometric approach that take into account primer efficiencies for both reference and target genes.
In our experimental conditions, 16S and 18S rRNAs, EEF1A and actB resulted not reliable as reference genes in terms of stability and relative levels of expression, in analogy to what is reported for other species [e.g. [6, 9, 10, 20]].
In particular, the expression of 16S rRNA allowed to identify two discrete groups of animals (60-fold differences in Ct; data not shown) in our sample. This did not correlate with sex, seasonality or the body size of the animals and was independent from the tissues considered. The fact that 16S rRNA is a mitochondrial gene suggests that this may be related with physiological (and/or metabolic) status of the animals. Similarly, we cannot exclude that 18S rRNA expression may be influenced by physiological conditions. Moreover, we found that actB reached higher levels of expression in the mantle when compared to the brain; it also resulted expressed with high interquartile ranges among samples, similarly to what resulted in other experimental settings . Surprisingly, the three algorithms did not indicate EEF1A as a suitable reference gene in our experimental conditions (except for Bestkeeper: mantle set a and SUB set b). However, its expression is comparable to the other putative reference genes (Figure 1). This resembles what recently found in other systems .
In order to account for the largest inter-individual variability in the octopus, we analyzed some target genes separately in two different groups of animals, one representing the whole population (set a) and one more restricted to sub-adults of the summer season (set b). The expression levels detected for each animal group suggest a correlation between data set heterogeneity and gene expression variability as deduced by the elevated Ct value variability obtained from samples belonging to octopuses of set a.
Moreover, octopuses of our samples showed low expression levels for transcription factors when compared with mRNA coding for proteins highly utilized in the metabolic pathway (i.e. FoxP and creb vs dat and TH; data not shown). It is interesting to note that dopamine-related genes were abundantly expressed in the optic lobes of O. vulgaris, confirming previously published data [review in ].
In conclusion, we were able to identify reference genes to be utilized for normalization in particular conditions in RT-qPCR experiments designed to test gene expression in different tissues of O. vulgaris and to describe inter-individual variability in gene expression in naïve octopuses.
All the octopuses were collected from the wild and maintained in seawater until sacrifice. The animals were anesthetized in a sea water solution using MgCl2 . A portion of the mantle (without the skin) and the brain was dissected out the animal. The brain was divided into its parts: supra-oesophageal mass (SEM), sub-oesophageal mass (SUB) and optic lobes (OL). Samples were transferred in plastic tubes (1 ml of Eurozol, Euroclone), immediately frozen in liquid nitrogen, and stored at -80°C until further processed.
The octopuses (Octopus vulgaris, N = 10) were collected from different locations of the Bay of Napoli (Italy) throughout the year (2006). Their weight ranged from 30 to 2100 g, spanning across a broad size/age range.
The octopuses (O. vulgaris, N = 15) were caught in the Bay of Napoli in the same season (June – July, 2006); their size ranged between 200–550 g.
RNA extraction and quantification
Total RNA was extracted using Eurozol (Euroclone) according to the manufacturer's instructions. Contaminating DNA was degraded by treating each sample with Turbo DNase Kit (Ambion) according to the instruction manual.
For all RNA samples, the absence of DNA contamination was tested by performing PCR with β-actin primers. The quantity and purity of total RNA extracted was estimated monitoring both the absorbance at 260 nm and ratios 260/280 nm and 260/230 nm by Nanodrop (ND-1000 UV-Vis Spectrophotometer; NanoDrop Technologies). The quality of RNA was evaluated by gel electrophoresis. Intact rRNA subunits (28S and 18S) were observed on the gel indicating minimal degradation of the RNA.
For each sample, 2 μg of total RNA extracted was retrotranscribed with SuperScript™ First-Strand Synthesis System for RT-PCR (Invitrogen) following the manufacturer's instructions. The cDNA was stored at -20°C until further use. cDNA was diluted 1:100 with H2O prior to use in RT-qPCR experiments, and 1:1500 when using 18S primers.
Isolation of target and reference genes
Target and reference genes were obtained from O. vulgaris using primers designed on conserved regions by means of bioinformatic analysis comparing homologous sequences available in Gene Bank for the different taxa (see Additional file 1).
To isolate FoxP, retrotranscription reactions were performed using SuperScript™ One-Step RT-PCR Systems (Invitrogen) on total RNA extracted from the brain of adult octopuses using TRIZOL® (Invitrogen) with inosinate degenerate primers (RPPfw 5'-AGACCGCCITTYACITAYGC-3' and AVWrev 5'-TCRTCIACVGTCCAIACIGC-3' for forward and reverse, respectively). Primers WKNfw 5'-TGGAAGAATGCCGTGCGCCA CA-3' and WKNrev 5'-TGTGGCGCACGGCATTCTTCCA-3' were used to obtain the 5' and 3' ends respectively using GeneRacer™ (Invitrogen) following the manifacturer's protocol. At the same time a PCR screening of cDNA libraries from octopus brain were performed to confirm the transcripts obtained (1111 bp).
For all cases (18S rRNA, ubi, dopamine transporter: dat, and tyrosine hydroxylase: TH) cDNAs were synthesized using SuperScript™ First-Strand Synthesis System for RT-PCR (Invitrogen) from total RNA extracted from brains of adult octopuses using Eurozol (Invitrogen) according to the manufacturer's instructions. The cDNA obtained was used as template in PCR reactions to amplify genes of interest using primer pairs designed by Primer3 software .
For 18S rRNA, a 442 bp fragment was amplified with primers designed on the basis of Eledone cirrhosa 18S rRNA sequence (GeneBank accession number: AY557467): 18Sfw 5'-CGTTTTCCTCGATCAAGAGC-3' and 18Srev 5'-CGAACTCGCGAAAGAAGAAG-3'.
For ubiquitin (ubi), a 324 bp cDNA sequence coding for ubi/ribosomal S27a protein was identified using the following primers: UBfw 5'-TGTCAAGGCAAAGATTCAAGA-3' and Ubrev 5'-GGCCATAAACACACCAGCTC-3'. These primers were designed on the basis of sequences present in Octopus eye EST library (http://cib.nig.ac.jp/dda/database/octopus.htm, last visited: Aug, 2008).
For dopamine transporter (dat) we analyzed the ortholog sequence alignment of vertebrates and invertebrates. Degenerated oligonucleotides were designed to amplify conserved regions. A 506 bp fragment coding for dat was amplified using the following primers: DATfor1 5'-TCKGGIAARGTDGTBTGGTT-3' and DATrev3 5'-ATIGCYCIGADCCNCCRAA-3'.
For tyrosine hydroxylase (TH), the most conserved regions of vertebrate and invertebrate ortholog sequences were identified and degenerated oligos were designed for PCR reactions. The primers Thfor1 5'-RTSTTYCAGWSYACICAGTA-3' and Threv2 5'-AAYTCVACRGTGAACCAGTA-3' were used to amplify a fragment of 539 bp.
cDNA amplified fragments were purified from gel agarose using QIAquick® Gel extraction kit (Qiagen) and cloned into pCRII-TOPO vector (Invitrogen) according to the instruction manual. The resulting plasmids were sequenced using M13 reverse primer and T7. Sequences of the PCR products obtained were analyzed by BLASTX and BLASTP programs.
The c-AMP response element binding protein (creb) cDNA sequence was found after screening O. vulgaris cDNA library constructed from mRNA of the supra-oesophageal mass. The cDNA sequence coding for Aplysia creb1α was used to screen the library (gift from Dr E.R. Kandel Laboratory). The sequence identified was 4313 bp long and coded for a protein of 296 amminoacids.
The resulting alignments of O. vulgaris translated sequences with invertebrate and vertebrate orthologues are summarized in Additional file 1.
Primer design for RT-qPCR
All RT-qPCR primers were designed by Primer 3 software . Table 1 lists the primers' sequences together with amplicon size. Target genes' sequences amplified by the primer pairs were evaluated with MFOLD software  in order to check for secondary structures at the site of primer binding. Specificity of PCR products was checked by melting curve analysis followed by gel electrophoresis and DNA sequencing.
Real Time PCR
The efficiency of each primer pair (Table 1) was calculated according to standard methods curves using the equation E = 10-1/slope. Five serial dilutions were set up to determine Ct values and reaction efficiencies for all primer pairs. Standard curves were generated for each oligonucleotide pairs using the Ct value versus the logarithm of each dilution factor.
Diluted cDNA was used as template in a reaction containing a final concentration of 0.3 μM for each primer and 1× FastStart SYBR Green master mix (total volume of 25 μl). PCR amplifications were performed in a Chromo4™ Real-Time Detector (Biorad) thermal cycler using the following thermal profile: 95°C for 10 min, one cycle for cDNA denaturation; 95°C for 15 sec and 60°C for 1 min, 40 cycles for amplification; 72°C for 5 min, one cycle for final elongation; one cycle for melting curve analysis (from 60°C to 95°C) to verify the presence of a single product. Each assay included a no-template control for each primer pair. To capture intra-assay variability all RT-qPCR reactions were carried out in triplicate (set a) or duplicate (set b). Fluorescence was measured using Opticon Monitor 3.1 (Biorad).
Stability analysis of reference genes
The distribution of the Ct was first calculated for each tissue considering all the samples and reported as box plot. The Ct values were obtained from the average of each animal tissue.
Three different gene normalization algorithms were utilized in this work: BestKeeper , geNorm  and NormFinder . For each one a step-wise exclusion method have been applied in order to identify the best candidate reference genes.
Expression analysis of target genes
trg: target gene
Raw Ct data were exported to a worksheet for further analysis. SPSS 16 was used for statistical analysis (One-way ANOVA followed Dunnett post-hoc test). Significance level was set at 5%.
We are indebted to Drs M. Branno and A. Arcucci who greatly contributed to creb identification and isolation in O. vulgaris. We are grateful to Luciana Borrelli for critical reading of the manuscript and suggestions. This work was funded by Istituto Banco di Napoli, Fondazione (Italy) and the network of Excellence "Marine Genomics Europe" (GOCE-04-505403).
- Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.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: 856-862. 10.1016/j.bbrc.2003.11.177View ArticlePubMedGoogle Scholar
- Skern R, Frost P, Nilsen F: Relative transcript quantification by quantitative PCR: roughly right or precisely wrong?. BMC Mol Biol. 2005 Apr 26;6(1):10. 2005, 6 (1): 10- 10.1186/1471-2199-6-10Google Scholar
- Thellin O, Zorzi W, Lakaye B, De Borman B, Coumans B, Hennen G, Grisar T, Igout A, Heinen E: Housekeeping genes as internal standards: use and limits. J Biotechnol. 1999, 75: 291-295. 10.1016/S0168-1656(99)00163-7View ArticlePubMedGoogle Scholar
- Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques. 2004, 37: 112-114.PubMedGoogle 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: 141-143. 10.1016/j.ab.2005.05.022View ArticlePubMedGoogle Scholar
- Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Commun. 2006, 345: 646-651. 10.1016/j.bbrc.2006.04.140View ArticlePubMedGoogle Scholar
- Liu DW, Chen ST, Liu HP: Choice of endogenous control for gene expression in nonsmall cell lung cancer. Eur Respir J. 2005, 26: 1002-1008. 10.1183/09031936.05.00050205View ArticlePubMedGoogle Scholar
- Mamo S, Gal AB, Bodo S, Dinnyes A: Quantitative evaluation and selection of reference genes in mouse oocytes and embryos cultured in vivo and in vitro. BMC Dev Biol. 2007, 7: 14- 10.1186/1471-213X-7-14PubMed CentralView ArticlePubMedGoogle Scholar
- Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods. 2000, 46: 69-81. 10.1016/S0165-022X(00)00129-9View ArticlePubMedGoogle Scholar
- Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30 (9 ): e36- 10.1093/nar/30.9.e36PubMed CentralView ArticlePubMedGoogle Scholar
- 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
- Borrelli L, Fiorito G: Behavioral Analysis of Learning and Memory in Cephalopods. Learning Theory and Behavior. Vol. 1 (Menzel, R., Volume Editor) of Learning and Memory: A Comprehensive Reference, 4 vols. (Byrne, J., Editor). 2008, 605-628.Google Scholar
- Hochner B, Shomrat T, Fiorito G: The octopus: a model for a comparative analysis of the evolution of learning and memory mechanisms. Biol Bull. 2006, 210: 308-317. 10.2307/4134567View ArticlePubMedGoogle Scholar
- Young JZ: The Anatomy of the Central Nervous System of Octopus vulgaris. 1971, Claredon Press, Oxford, Google 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. Biotechnology Letters. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47View 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: 5245-5250. 10.1158/0008-5472.CAN-04-0496View ArticlePubMedGoogle Scholar
- Fernandes JM, Mommens M, Hagen O, Babiak I, Solberg C: Selection of suitable reference genes for real-time PCR studies of Atlantic halibut development. Comp Biochem Physiol B Biochem Mol Biol. 2008, 150: 23-32. 10.1016/j.cbpb.2008.01.003View ArticlePubMedGoogle Scholar
- Olsvik PA, Lie KK, Jordal AEO, Nilsen TO, Hordvik I: Evaluation of potential reference genes in real-time RT-PCR studies of Atlantic salmon. BMC Mol Biol. 2005 Nov 17;6:21. 2005, 6: 21-Google Scholar
- Zhong Q, Zhang Q, Wang Z, Qi J, Chen Y, Li S, Sun Y, Li C, Lan X: Expression profiling and validatio of potential reference genes during Paralichthylus olivaceus embryogenesis. Mar Biotechnol (NY). 2008, 10 (3): 310-8. 10.1007/s10126-007-9064-7View ArticleGoogle Scholar
- Toegel S, Huang WW, Piana C, Unger FM, Wirth M, Goldring MB, Gabor F, Viernstein H: Selection of reliable reference genes for qPCR studies on chondroprotective action. BMC Mol Biol. 2007, 8: 13- 10.1186/1471-2199-8-13PubMed CentralView ArticlePubMedGoogle Scholar
- Hibbeler S, Scharsack JP, Becker S: Housekeeping genes for quantitative expression studies in the three-spined stickleback Gasterosteus aculeatus. BMC Mol Biol. 2008, 9: 18- 10.1186/1471-2199-9-18PubMed CentralView ArticlePubMedGoogle Scholar
- Messenger JB: Neurotransmitters of cephalopods. Invertebrate Neuroscience. 1996, 2: 95-114. 10.1007/BF02214113.View ArticleGoogle Scholar
- Messenger JB, Nixon M, Ryan KP: Magnesium chloride as an anaesthetic for cephalopods. Comp Biochem Physiol C Comp Pharmacol. 1985, 82: 203-205. 10.1016/0742-8413(85)90230-0.View ArticleGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, 365-386. Totowa, NJ: Humana PressGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595PubMed CentralView 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.