Email updates

Keep up to date with the latest news and content from BMC Molecular Biology and BioMed Central.

Open Access Research article

Evaluation of potential reference genes for real time RT-PCR studies in Atlantic halibut (Hippoglossus Hippoglossus L.); during development, in tissues of healthy and NNV-injected fish, and in anterior kidney leucocytes

Aina-Cathrine Øvergård12*, Audun Helge Nerland12 and Sonal Patel1

Author Affiliations

1 Institute of Marine Research, PO Box 1870 Nordnes, N-5817 Bergen, Norway

2 The Gade Institute, University of Bergen, PO Box 7804, N-5021 Bergen, Norway

For all author emails, please log on.

BMC Molecular Biology 2010, 11:36  doi:10.1186/1471-2199-11-36


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2199/11/36


Received:11 December 2009
Accepted:11 May 2010
Published:11 May 2010

© 2010 Øvergård et al; licensee 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.

Abstract

Background

Real time RT-PCR has become an important tool for analyzing gene expression in fish. Although several housekeeping genes have been evaluated in Atlantic halibut (Hippoglossus Hippoglossus L.), appropriate reference genes for low copy mRNA transcripts at the earliest developmental stages have not been identified. No attempts have been reported to identify suitable reference genes in halibut infected with NNV or in stimulated halibut leucocytes. In this study, β-actin1 (ACTB1), elongation factor 1 alpha (EF1A1), hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1), ribosomal protein L7 (RPL7), tubulin beta 2C (Tubb2C), and ubiquitin-conjugating enzyme (UbcE) were evaluated as reference genes for normalization of real time RT-PCR data during Atlantic halibut development, in tissue of healthy and NNV-infected fish, and in in vivo and in vitro stimulated anterior kidney leucocytes.

Results

The expression of all six genes was relatively stable from the unfertilized egg until 12 day degrees post fertilization (ddpf). However, none of the selected genes were found to be stably expressed throughout halibut development. The mRNA levels of the six genes increased from 18 ddpf, when zygotic transcription is likely to be activated, and stabilized at different time points. The Excel-based software programs BestKeeper, geNorm, and NormFinder ranked EF1A1 and UbcE as the best candidate reference genes before activation of zygotic transcription, and RPL7 and EF1A1 as the best candidates after hatching. EF1A1 and RPL7 were also listed as the best reference genes when exploring the expression levels of the six genes in various halibut organs, both in non-injected fish and in mock- and NNV-injected fish. None of the reference genes were found optimal for normalization of real time RT-PCR data from in vitro stimulated anterior kidney leucocytes.

Conclusion

Generally, it was found that EF1A1 and RPL7 were the genes that showed least variation, with HPRT1 and UbcE as intermediate genes, and ACTB1 and Tubb2C as the least stable ones. None of the six reference genes can be recommended as reference gene candidates in ConA-PMA stimulated leucocytes. However, UbcE can be a good candidate in other experimental setups. This study emphasizes the need for reference gene evaluation, as universal reference genes have not been identified.

Background

Real time reverse transcriptase polymerase chain reaction (real time RT-PCR) has become a widely used method for gene expression analysis, and it is a useful method for studying immune related genes and host-pathogen interactions. It is more accurate and sensitive than traditional methods like RT-PCR and northern blotting [1], but normalization of the assay is critically important as differences in loading amounts of total RNA in the RT reaction, variations in RT efficiency and RNA integrity, instrumental errors, and the presence of PCR inhibitors have to be accounted for [2]. Housekeeping genes are often used as internal reference genes. Ideally, genes chosen should have stable gene expression among individuals, organs and cells, during different developmental stages, and various experimental treatments. The housekeeping genes chosen should thus be validated for each new experimental setup. Also, the use of a single housekeeping gene has been found to be insufficient [3]. Thus, it is important to evaluate and establish a two-gene normalization strategy for normalization of real time RT-PCR data. While establishing such a strategy one should bear in mind not to use genes involved in the same biological process to avoid co-regulation.

Larvae hatching at a primitive state, followed by a long developmental period has made the farming of the marine flatfish Atlantic halibut (Hippoglossus Hippoglossus L.) challenging [4,5]. Several microorganisms have been associated with high mortality of halibut eggs and larvae at this stage of life when the halibut immune system is poorly developed [6]. One of the most important pathogen in economical terms affecting halibut during larval and early juvenile stages is the nervous necrosis virus (NNV). NNV is the causative agent of Viral Encephalopathy and Retinopathy (VER), and the major site for virus replication is within the central nervous system [6]. Much work has been done to characterize various NNV strains and in vaccine development [7-11]. However, analyzing halibut immune related genes in response to NNV-infection has not been optimal as suitable reference genes for such experimental setups have not been evaluated.

Several commonly used reference genes have been applied in real time RT-PCR studies of Atlantic halibut gene expression, including β-actin (ACTB), 18S rRNA, elongation factor 1 alpha (EF1A), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) [12-15]. Recently, several housekeeping genes have been evaluated during halibut development where GAPDH was found to be unsuitable as a reference gene in halibut egg and larvae [16,17]. Moreover, Fernandes et al. [16] found 18S rRNA to be rather stable from the two cell stage at about 1 day degree post fertilization (ddpf) to the first feeding stage. However, other genes tested, including ACTB appeared to be developmentally regulated at early stages [16], as seen in Japanese flounder (Paralichthys olivaceus) [18]. Also, the time points sampled before activation of zygotic transcription were not sufficient to evaluate potential reference gene candidates, and appropriate reference genes for low copy mRNA transcripts at the earliest developmental stages remains to be identified. Some potential reference genes have been analyzed in several tissues of Atlantic halibut, including ACTB and some ribosomal proteins. ACTB was found to be inappropriate for normalization of muscle samples during fasting and re-feeding of halibut males [19], but the ACTB gene expression was not tested in other tissues. Paralogous genes of ribosomal protein L19 (RPL19) and RPL22 were tested in halibut larva and in several organs of three year old halibut with promising results [20]. However, no statistical evaluation of the gene expression stability was conducted in this study, and further evaluation of ribosomal proteins as reference gene candidates is needed.

Based on earlier gene expression studies in halibut [16,17], Japanese flounder [18], medaka (Oryzias latipes) [21], and Atlantic salmon (Salmo salar) [22,23], six housekeeping genes, ACTB1, EF1A1, hypoxanthine-guanine phosphoribosyltransferase 1 (HPRT1), RPL7, tubulin beta 2C (Tubb2C), and ubiquitin-conjugating enzyme E2L3 (UbcE) were selected. In the present study, the selected genes were evaluated as reference genes for normalization of real time RT-PCR studies during halibut development from the unfertilized egg to juveniles, as well as in different organs of healthy fish. Evaluation of these genes was also carried out in tissues sampled from fish infected with NNV and in anterior kidney leucocytes isolated from in vivo NNV-stimulated fish and further in vitro stimulated with ConA-PMA. This study will give valuable information regarding reference genes for real time RT-PCR studies during halibut development, as it will give a more detailed evaluation of the earliest developmental stages before activation of zygotic transcription. Also, results may serve as a good foundation for further search of reference genes in experimental setups including infection and leucocyte stimulation, both in halibut and other fish species.

Methods

Fish stocks and sample collection

Halibut eggs and larvae were collected in two different experiments conducted a year apart, and reared at the Austevoll Aquaculture Reseach Station, Norway and at the Institute of Marine Research (IMR, Bergen), Norway following normal rearing conditions [4,24,25]. The larvae were fed enriched Artemia until weaning to commercially available dry pellets at 1150 day degrees post hatching (ddph). In the first experiment, samples were taken from unfertilized eggs at stripping, and at 0.04, 0.08, 0.17, 0.25, 0.9, 6, 12, 18, 24, 30, 36, 42, 48, and 72 ddpf, and further of larvae at 0, 18, 30, 42 and 60 ddph. In the second experiment, samples were taken from unfertilized egg at stripping, and from fertilized eggs at 0, 30, and 72 day degrees post fertilization (ddpf), and at 0, 18, 42, 60, 84, 102, 126, 144, 168, 186, 210, 228, 252, 270, 318, 354, 438, 522, 606, 690, 774, 858, 954, 1026, 1110, 1194, 1278, 1458, and 1638 ddph. For developmental stages corresponding to the different time-points, see Pittman et al. [26] and Patel el al. [27]. For total RNA purification, five eggs or larvae were pooled at each time point until 210 ddph. From 228 ddph to 1638 ddph, either one individual, or a pool of 2-5 larvae or juveniles depending on the size and biomass, were employed after dividing them in two just behind the intra-peritoneal cavity, and the anterior part was used for total RNA isolation. For each time point four parallels were purified.

Generally, fish 6 months or older were acclimatized upon arrival at IMR (Bergen), Norway, and reared in seawater (34 ppm salinity) at 9°C, and fed commercial pellets twice a day. If injected, fish were anesthetized with benzocain (The Norwegian medicine depot) at a concentration of 60 mg/l seawater, while for tissue sampling an overdose of benzocain was employed.

Four individuals, approximately 1 year old and weighing between 70 - 150 g, were obtained from Austevoll Aquaculture Reseach station, Norway. For total RNA isolation, samples were collected from thymus, spleen, skin, heart, anterior and posterior kidney, pectoral fin, gills, brain, liver, anterior and posterior gut, red and white muscle, stomach, and eye.

Fish approximately 6 months old and weighing approximately 30 g, were obtained from Aga Marine, Bømlo, Norway. They were injected intra peritoneal (i.p) with either 200 μl of L-15 medium (Sigma) or 200 μl of 1 × 108.5 TCID50 NNV. The NNV used was isolated from natural outbreak affecting halibut [28], and propagated in cell culture as described previously [29]. Samples from brain, retinal nerve, thymus, anterior kidney, spleen, and liver were taken from six individuals from both NNV-injected and mock-injected group (L-15 medium injected group) at 6, 12, 24, 72 hours (hrs) and 1, 2, 3, and 4 weeks post-injection. Anterior kidney leucocytes were isolated as described previously [30] from four fish from each group 10 weeks post-injection. The isolated leucocytes were further divided into two aliquots and plated out in 24-well culture plates. One aliquot served as a non-stimulated control group and the other group was stimulated with 10 μg/ml with 5 ng/ml phorbol myristate acetate (Calbiochem) (ConA-PMA) as previously described [30]. Cells were harvested at 0 hrs (non-stimulated) and from both non-stimulated control and ConA-PMA stimulated group after 4, 12, 18, and 24 hrs of incubation, pelleted, and frozen at -80°C until use.

All eggs, larva, and tissue samples were snap-frozen in liquid nitrogen immediately after sampling, and stored at -80°C until use. All fish experiments were performed according to national legislation and approved by the Norwegian animal experimentation board.

Isolation of total RNA and cDNA synthesis

Total RNA from developmental stages and from healthy fish was isolated using TRI reagent (Sigma) according to the Trizol reagent protocol described by Invitrogen, with a few modifications as described previously [14]. The tissue samples from injected fish were purified using the iPrep™ TRIzol® Plus RNA Kit (Invitrogen). Total RNA from leucocytes was purified by a combined Trizol (Invitrogen) and RNeasy (Qiagen) method. Briefly, the aqueous phase from the chloroform phase separation (Trizol method) were added 1 volume of 70% ethanol, transferred to an RNeasy spin column, and total RNA was further purified according to the RNeasy® Mini Handbook (Qiagen).

The concentration and the purity of the total RNA were assessed with a NanoDrop Spectrophotometer (NanoDrop Technologies), and the quality of random samples was analyzed with an Agilent 2100 Bioanalyzer (Agilent Technologies). A 260/280 nm absorbance ratio above 1.8 was accepted as pure, and RNA integrity numbers (RIN) above 5 were considered as good quality total RNA. Total RNA was reversed transcribed using a Reverse Transcriptase Core Kit (Eurogentec) and random nonamers as primers. All total RNA was reversed transcribed in 30 μl leucocytes where 300 ng of total RNA was used due to low total RNA concentrations. cDNA was stored at -20°C until use.

Cloning and DNA sequencing

Expressed sequence tags (EST) representing EF1A1 [GeneBank: GE628160], ACTB1 [GeneBank: GE628929, GE629238, GE629769, GE631177, GE631405, GE631443], Tubb2C [GeneBank: GE630028, GE631294, GE631482], and HPRT1 [GeneBank: GE627320, GE631255] were identified in a sequence database generated on the basis of Atlantic halibut cDNA libraries [30].The RPL7 EST [GeneBank: EB035355] and the UbcE EST [GeneBank: CF931657] were retrieved from NCBI http://www.ncbi.nlm.nih.gov/sites/entrez webcite[31,32]. The EF1A1 clone was amplified in XL 10-Gold® ultracompetent E. coli, and plasmid DNA was purified using the QIAprep® Spin Miniprep Kit for Microcentrifuge (Qiagen). The plasmid insert was sequenced by primer walking using BigDye® Terminator v3.1 Cycle Sequencing kit and run on an ABI prism 7700 automated sequencing apparatus (Applied Biosystems). The PCR cycling was carried out as follows: 96°C for 5 min followed by 25 cycles of 96°C for 10 sec, 50°C for 5 sec, and 60°C for 4 min. The full length sequence [GeneBank: EU561357] was retrieved by RLM-RACE (Invitrogen) as described previously [14], and the PCR products were purified using QIAquick PCR purification kit (Qiagen) and sequenced as described above.

To confirm the EST sequences of the ACTB1, HPRT1, RPL7, Tubb2C, and UbcE gene, primers were designed based on EST alignments. Amplification of cDNA was performed in 25 μl reactions GoTaq® Flexi DNA Polymerase with (Promega), using 200 nM of each PCR primers, and 0.5 μl cDNA. The PCR cycling was carried out as follows: 95°C for 5 min followed by 30 cycles of 95°C for 30 s, 55°C for 30 s, 72°C for 1 min per kb, and finally 72°C for 7 min. The PCR products were enzymatically cleaned using Shrimp Alkaline Phosphatase in combination with Exonuclease I (USB Corporation), and sequenced as described above. To confirm the exon-exon boundaries, genomic sequences [ACTB1; GeneBank: GQ465758, EF1A1; GeneBank: EU561358, HPRT1; GeneBank: GQ465759, RPL7; GeneBank: FJ008911, Tubb2C; GeneBank: GQ465760, and Ubc; GeneBank: FJ008912] were amplified by PCR as described above, using 100 ng of halibut genomic DNA purified from whole blood erythrocytes by the phenol/chloroform method, protocol I for frozen blood [33] as template. The PCR products were cloned using the TOPO TA Cloning® Kit for Sequencing (Invitrogen). The plasmids were purified and the inserts were sequenced as described above.

Primer and probe design

Primers and probes were designed with Primer express software 3.0 (Applied Biosystems), according to the manufacturer's guidelines. Either the probe or one of the primers was designed such that they spanned an exon-exon boundary, to avoid amplification of genomic DNA. Both primers and probes were screened for homo- or cross-dimers and hairpin structures that could affect the efficiency of the PCR reaction. Five-point standard curves of 4-fold dilution series (1:1 - 1:256) from pooled cDNA were used for calculation of the PCR efficiency, given by the equation E% = (101/slope - 1) × 100 [34], and for revealing PCR poisoning. The slope was calculated from the linear regression model fitted from the log-transformed cDNA concentrations plotted against the Ct values. Generally, a PCR efficiency between 90-110% is considered acceptable (Applied Biosystems). The primer and probe sequences with corresponding PCR efficiencies are listed in Table 1. Retinal nerve samples were analyzed for the presence of viral RNA using a real time RT-PCR assay previously described [35].

Table 1. Primers and probes used for real time RT-PCR analysis.

Real time RT-PCR assay and data analysis

The PCR reaction mix contained 1× TaqMan Fast PCR Master Mix (Applied Biosystems), 900 nM of each primer, 200 nM TaqMan probe The PCR cycling was carried out as follows: 95°C for 20 sec, 40 cycles of 95°C for 1 sec followed by 60°C for 20 sec. Two technical replicates were run for each sample on the 7900 HT Fast Real-Time PCR System (Applied Biosystems), and if the percent deviation was below 5% the mean Ct values for each sample were used for further analysis. Non template controls and samples negative for the reverse transcriptase enzyme were included. All assays were tested for amplification of genomic DNA, using 200 ng of genomic DNA as template in the real time RT-PCR reaction. The data collected was analyzed by the Excel-based programs BestKeeper [36], geNorm [3], and NormFinder [37]. Three programs were used as none of these are accepted methods on their own for evaluation and ranking of reference genes, thus by using three independent methods a more reliable result was expected.

Results and Discussion

For real time RT-PCR analysis, it is essential to use good quality total RNA with high degree of purity [38]. High concentrations of contaminants may give PCR poisoning and decrease reverse transcriptase activity, while degraded total RNA can give inaccurate measurement of target gene RNA if the region of interest is degraded before analysis. In the present study, RNA quality and purity measured by Bioanalyzer and NanoDrop respectively, showed that high quality total RNA was extracted with RIN values above 6, normally close to 10, and ratios of 260/280 nm absorbance between 1.8 and 2.1. When using the 2-ΔΔCt method it is crucial to have a primer and probe design that gives PCR efficiency close to 100% when the correlation coefficient (R2) of the linear regression line is close to one [39]. The PCR efficiency of the six assays established in this study was between 90 and 98% (Table 1), and PCR poisoning could not be seen (R2 close to 1). No signal was detected when samples lacking the reverse transcriptase enzyme or genomic DNA were analyzed, indicating that the assays did not detect any genomic DNA contaminants.

When the collected data was analyzed by BestKeeper, geNorm, and NormFinder, often the three software programs deviated in the ranking order of the analyzed genes reflecting the differences in the estimation approach by the three programs. Both BestKeeper and geNorm uses a pair-wise comparison approach, and are highly dependent on the assumption that none of the genes being analyzed are co-regulated [3,36]. BestKeeper ranks reference genes according to a correlation with a BestKeeper index, calculated based on the geometric mean of the Ct values of the candidate reference genes with a standard deviation (SD) below 1. The geNorm program calculates an M value that corresponds to the average pair-wise variation of a single reference gene to all other genes, and allows for a repetitive procedure where the least stably expressed candidate is removed and new M values are calculated. The gene with the lowest M value should be the most stably expressed gene, and the recalculation of the M values can alter the ranking of the reference genes and decrease the M values (Table 2). NormFinder on the other hand uses a model-based approach and ranks according to a minimal combined inter- and intra-group expression variation [37], and should be a more robust approach. However, in all the three programs used it seemed like the ranking of a given reference gene was highly dependent on the set of candidate reference genes included in the analysis. Not only was the reliability dependent on the assumption of minimal co-regulation, but also that the set of reference genes included in the analysis were relatively stably expressed. In the present study, the chosen reference genes were highly deviating in the expression stability. Often, two of the selected genes had a high degree of variation within the given sample set, one or two of the genes were intermediate, and two or three genes were relatively stable in expression. Seemingly, the software programs often chose the intermediate reference gene as the best candidate, as the expression pattern of this gene was more similar to the general expression pattern than to the most stably expressed gene. Therefore, conclusions drawn from this study were not based on the ranking order made by the three software programs alone.

Table 2. Average expression stability (M values) of the six potential reference genes according to geNorm.

Expression of candidate reference genes during halibut development

As there is a delay between fertilization and the activation of zygotic transcription, many gene products are deposited in the egg during oogenesis to execute certain basic cellular functions [40]. The development of Atlantic halibut eggs and larvae progresses through major changes throughout development [26], evidently the activation of zygotic transcription is a major factor affecting the expression of several halibut housekeeping genes [16]. All six reference genes analyzed here were detected likely as maternal transcripts in the halibut eggs, as seen in other studies of gene expression during early teleost development [16,18,41-43]. The two sampling series, which were taken one year apart and with different brood-stock, revealed concurrence in developmental regulation of the genes with stage-specific expression stability, and the halibut eggs showed a relatively high abundance in mRNA storage (Table 3). All genes were found to have a relatively stable expression level from the unfertilized egg until 12 ddpf (Figure 1). After 12 ddpf the mRNA level of all genes except ACTB1 increased, as the zygotic transcription was likely to be activated, and stabilized at different time points. A minor drop in the ACTB1 mRNA level at 18 ddpf could be seen, and an increase in the mRNA level was not present until 24 ddpf. ACTB1 expression was previously reported during early halibut development [16], correlating with the expression pattern reported here. However, a more accurate estimate of the onset of zygotic transcription was shown in the present study.

Table 3. Average Ct values ± SD of the six reference genes in the different experimental setups.

thumbnailFigure 1. Experiment 1. Expression profiles of the six potential reference genes during Atlantic halibut early development. Data is shown from unfertilized egg until 60 ddph as average Ct values ± SD where n = 4. Abbreviations: Unf - unfertilized, f - day degrees post fertilization, h - day degrees post hatching.

Likely, the stage specific gene expression of the selected genes would result in an invalid data analysis using BestKeeper, geNorm, and NormFinder, as the three software programs assume that the candidate reference genes chosen are not co-regulated [3,36,37]. The data collected throughout the developmental stages was therefore not analyzed by the three software programs. However, UbcE showed the smallest level of regulation during halibut development (Figure 2A), as it decreased from an average Ct value of 23.7 before 12 ddpf to 22.9 after hatching (Table 3). This is supported by findings in Japanese flounder where UbcE, during embryogenesis, was found to be the second best reference gene after 18S rRNA [18]. Two studies have found 18S rRNA to be rather stable during Atlantic halibut development with no developmental regulation [16,17]. But due to the lack of introns [44], high expression level [45], and as the mRNA fraction of total RNA is not always represented properly by the rRNAs expression level [46], 18S rRNA is not an optimal candidate for normalization of low copy mRNA target genes. However, it can be questioned whether the use of maternal and zygotic transcripts as references within the same experiment would be optimal, and thus one should be critical in the use of a reference gene throughout development.

thumbnailFigure 2. Experiment 2. Variation of the potential reference genes during Atlantic halibut development presented as fold change for A) developmental stages from the unfertilized egg until 1638 ddpf, B) developmental stages from the unfertilized egg until 12 ddph, and C) developmental stages from hatching (0 ddph) until 1638 ddph. Expression levels of selected genes are shown as medians (square), 25th percentile to 75th percentile (boxes), non-outlier ranges (whiskers), and outliers (asterisks).

All genes were relatively stable from the unfertilized egg until 12 ddpf, although, EF1A1, UbcE, and RPL7 showed least variation in expression (Figure 2B). After hatching, all genes appeared to be stabilized at different stability levels (Figure 2C). To assess the most stably expressed candidate reference gene at the different expression levels before and after zygotic activation, the data collected was analyzed by BestKeeper, geNorm, and NormFinder. Before zygotic activation, from the unfertilized egg until 12 ddpf, EF1A1 was found to be the best candidate by all programs (Table 4). UbcE was ranked as the second best by geNorm and BestKeeper, while NormFinder identified Tubb2C as the second best gene. However, the UbcE and RPL7 expression were apparently more stable than the Tubb2C expression (Figure 2B), and the NormFinder ranking could thus be questioned. Also, BestKeeper showed a relatively good correlation between EF1A1 and UbcE (r = 0,789), indicating that they could be suitable for a two-gene normalization strategy. After hatching, all programs ranked RPL7 as the best gene, while Tubb2C and ACTB1 were identified as the most unstably expressed genes. The software programs differed in the ranking of the intermediate reference genes after hatching, as NormFinder listed UbcE > EF1A1 > HPRT1, geNorm listed EF1A1 > HPRT1 > UbcE, while BestKeeper listed EF1A1 > UbcE > HPRT1. The expression level of EF1A1 was seemingly more stable than UbcE (Figure 2C), and thus RPL7 and EF1A1 could be recommended for normalization after hatching, consistent with previous reports [17,21]. However, they cannot be recommended for a two-gene normalization strategy, as both genes are involved in protein biosynthesis, and some co-regulation may be expected. If a two-way normalization strategy is to be chosen, UbcE in combination with RPL7 or EF1A1 is likely to be the best alternative amongst the genes analyzed, as BestKeeper estimated a higher inter-gene correlation for RPL7/UbcE (r = 0.829) and EF1A1/UbcE (r = 0,769) than for RPL7/HPRT1 (r = 0.828) and EF1A1/HPRT1 (r = 0.666).

Table 4. Ranking output of the six reference genes according to their expression stability.

Expression of candidate reference genes in organs of one year old halibut

When the expression levels of the six candidate reference genes were explored in various organs of one year old halibut (Figure 3), the highest mRNA level was seen for EF1A1, followed by RPL7, Tubb2C, ACTB1, HPRT1, and UbcE (Table 3). This is in accordance with the mRNA levels of ACTB1, EF1A1, HPRT1, Tubb2C, and UbcE seen during developmental stages after activation of zygotic transcription. However, the average mRNA level of RPL7 in 1 year old halibut increased by 2 Ct, indicating that RPL7 has an decrease in gene expression level at some stage after the age of 6 months. Highest variation in Ct values amongst various organs was seen for ACTB1, while EF1A1 showed least variation (Figure 3). Considering individual variations seen within organs, it could be noted that the variation was generally high in anterior and posterior gut for all genes, and also in muscle and eye for ACTB1 expression, and in stomach, liver, and eye for Tubb2C expression. However, it should be mentioned that the conclusion drawn is based on four individuals only.

thumbnailFigure 3. Expression profiles of the six potential reference genes in various organs of adult Atlantic halibut. Data is shown as average Ct values ± SD (n = 4). Abbreviations: T - thymus, S - spleen, AK - anterior kidney, PK - posterior kidney, G - gills, F - pectoral fin, AG - anterior gut, PG - posterior gut, Sk - skin, M - white muscle, RM - red muscle, H - heart, St - stomach, L - liver, E - eye, and B - brain.

Tubb2C and ACTB1 were listed as the most unstably expressed genes by all three software programs when the real time RT-PCR data from various organs was analyzed with geNorm, BestKeeper, and NormFinder. EF1A1 was ranked as the best candidate by geNorm and BestKeeper, followed by RPL7, HPRT1, and UbcE, while NormFinder proposed UbcE as the best candidate gene, followed by RPL7, EF1A1, and HPRT1 (Table 4). The different expression pattern of the six genes chosen are likely to affect the NormFinder ranking, as UbcE is seemingly not the most stably expressed gene (Figure 3). EF1A1 and RPL7 can thereby be recommended as the best candidate reference genes in one year old halibut, supported by findings in medaka [21], Atlantic salmon [23,47], and zebrafish (Danio reio) [43]. As mentioned previously, EF1A1 and RPL7 may be co-regulated, and HPRT1 in combination with EF1A1 (r = 0.656) or RPL7 (r = 0.633) was estimated by the BestKeeper program to be the best choice for a two-way normalization strategy in organs of adult halibut.

Expression of candidate reference genes in organs of NNV-injected fish

Generally, altered mRNA levels of host genes can be expected following viral replication and production of virus antigens, and the importance of testing reference genes under such conditions is high. A previously established assay detecting NNV RNA2 was applied to retinal nerve samples, as to ensure viral replication in the NNV-injected fish. All mock-injected fish were negative for NNV, while 99% of the injected fish were positive (Figure 4). As early as 6 hrs post-injection NNV could be detected in the retinal nerve samples, likely due to circulating virus particles. Four weeks post NNV-injection the Ct values were increased, presumably as the virus had started to propagate in the fish retina.

thumbnailFigure 4. Detection of NNV in halibut retinal nerve samples. The data is represented as average Ct values ± SD (n = 6). Control samples that were undetected for NNV were set to 40 Ct. Abbreviations: h - hrs, w - weeks.

When the mRNA level of the six candidate reference genes was explored in various organs of mock-injected and NNV-injected fish, minor changes in the average Ct values (Table 3) and the general expression levels (Figure 5) were found between the two groups, presumably as the viral replication at these time points was not sufficient to reduce the level of housekeeping gene mRNA. Tubb2C was the gene showing the highest average increase in mRNA level in the NNV-infected group, with a decrease of 0.3 Ct in spleen. As expected, the highest mRNA level was seen for EF1A1, followed by RPL7, Tubb2C, ACTB1, HPRT1, and UbcE mRNA, respectively (Table 3). ACTB1 was the gene showing the highest degree of variation between organs with a SD of 2.5, with a difference in average Ct of 5.5 between anterior kidney and brain. The average expression level of RPL7 was at the same level as larvae and juveniles, increased compared to one year old fish, also when comparing the specific organs. As the fish used in the infection experiment were approximately 6 months, this supports the proposed decrease of expression between 6 months and one year. The use of RPL7 as a reference gene in this period should therefore be closely monitored, as to reveal a possible decrease in mRNA level during this developmental period.

thumbnailFigure 5. Variation of the potential reference genes during a NNV-infection. Data is shown as fold change of the six reference genes in all organs of mock- and NNV-injected fish. Fold change is shown as medians (square), 25th percentile to 75th percentile (boxes), non-outlier ranges (whiskers), and outliers (asterisks).

NormFinder, geNorm, and BestKeeper showed some concurrence in the ranking of reference genes when analyzing the inter-organ variation, listing EF1A1 as the best reference gene in contrast to data published on European seabass (Dicentrarchus labrax) [42], and Tubb2C and ACTB1 as the worst (Table 4). However, some disagreement in the ranking of the intermediate genes HPRT1, RPL7, and UbcE was seen, as NormFinder ranked UbcE > RPL7 > HPRT1 in mock-injected group and RPL7 > UbcE > HPRT1 in NNV-injected group, while geNorm and BestKeeper ranked RPL7 > HPRT1 > UbcE in both groups. However, HPRT1 and RPL7 were seemingly lesser regulated than UbcE (Figure 5), indicating that the BestKeeper and geNorm ranking could be more reliable.

When analyzing the intra-organ variation the ranking of the reference genes was generally quite stable between the treatments, though with variable ranking patterns given by the different software programs (Table 4). Again, this reflects the differences in the software programs used, and the varying expression pattern between the six genes selected. As seen in tissue samples of healthy fish, ACTB1 and Tubb2C were ranked as the most unstable genes in the different tissues tested after mock- and NNV-injection, with some exceptions. The generally high variation in the ACTB1 and Tubb2C expression compared to the other genes tested, gave a basis of comparison that was highly variable between the samples. The BestKeeper software allowed for the removal of genes with a SD above 1 [36]. However, liver was the only organ where the two most unstably expressed genes, ACTB1 and Tubb2C, could be removed. The BestKeeper ranking in liver was thereby the most reasonable compared to the general expression pattern of the six genes (Figure 6), ranking RPL7 as the best gene, followed by EF1A1, UbcE, HPRT1, Tubb2C, and ACTB1 as the most unstable gene. Also, NormFinder seemingly failed to identify the most stably expressed reference gene, and listed the more intermediate gene as the best candidate instead. This gave UbcE the best ranking in brain, spleen, anterior kidney, and liver, while HPRT1 was listed as the best candidate in eye and thymus. The geNorm ranking was found to be the most reasonable compared to the general expression pattern of the reference genes (Figure 6), as the program compare the expression ratio of different genes by eliminating the worst ranked reference gene in repetitive steps. However, some apparently poor ranking was seen, illustrated by the ranking of UbcE as the best candidate reference gene in thymus, when RPL7 was seemingly the most stably expressed gene (Figure 6).

thumbnailFigure 6. Variation of the potential reference genes in different organs of NNV-injected fish. Data is shown as fold change of the six reference genes in the different organs of mock- and NNV-injected fish. Fold change is shown as medians (square), 25th percentile to 75th percentile (boxes), non-outlier ranges (whiskers), and outliers (asterisks).

Expression of candidate reference genes in in vivo and in vitro stimulated leucocytes

Anterior kidney leucocytes were isolated from both mock-injected group and NNV-injected group 10 weeks post-injection and further in vitro stimulated with ConA-PMA, mitogens that have previously shown to stimulate fish leucocytes [48,49]. The mRNA levels of the six genes were all decreased compared to those in tissue samples (Table 3), likely due to the lesser amount of total RNA used in the reverse transcriptase reaction. A general decrease in the mRNA level of the six genes was seen during the first 18 hrs (Figure 7), probably as many leucocytes were going into apoptosis. Apoptosis is a tightly regulated process involving activation and inhibition of several genes at the transcriptional level and prelytic DNA fragmentation [50,51]. The up-regulation of such apoptosis-regulated genes can lead to difficulties in the analysis of real time RT-PCR data, as the stability of the reference genes may be affected by the apoptotic cells. Also, during the first 18 hrs post-stimulation the mRNA levels of the six reference genes were generally lower in the anterior kidney leucocytes stimulated with ConA-PMA compared to the non-stimulated control cells. Various apoptosis-related genes and immune-response genes were isolated in a cDNA library based on anterior kidney leucocytes stimulated with ConA-PMA from Japanese flounder [52,53]. This indicates an up-regulation of such genes in response to ConA-PMA stimuli, and could be diluting or regulate the expression of reference gene mRNA. Interestingly, at 24 hrs the expression of ACTB1, HPRT1, and Tubb2C compared to the levels at 18 hrs were increased in the non-stimulated cells that were isolated from in vivo NNV stimulated fish, and in the cells stimulated in vitro with ConA-PMA. However, due to high variation, the small sample size (n = 4), and as the experimental design was limited regarding the effect of tanks, no conclusion could be drawn.

thumbnailFigure 7. Expression profiles of the six potential reference genes in anterior kidney leucocytes after in vivo stimulation with NNV and in vitro stimulation with ConA-PMA. Data is shown as average Ct values ± SD (n = 4). Abbreviations: Ctr-Mock - non-stimulated control cells isolated from mock-injected fish, Ctr-NNV - non-stimulated control cells isolated from NNV-injected fish, Stim-Mock - ConA-PMA stimulated cells isolated from mock-injected fish, Stim-NNV - ConA-PMA stimulated cells isolated from NNV-injected fish.

When analyzing the data from anterior kidney leucocytes by NormFinder, geNorm, and BestKeeper, the software programs deviated in the ranking of the six genes (Table 4). UbcE and RPL7 were generally listed as good candidates, followed by EF1A1 and HPRT1, and with ACTB1 and Tubb2C at the bottom having the lowest rankings. Despite the good ranking of UbcE and RPL7, the difference in reference gene expression through time and between the groups renders them not optimal for normalization of real time RT-PCR data from ConA-PMA stimulated anterior kidney leucocytes. This is in contrast to the findings in Atlantic salmon where EF1A and ACTB1 were found suitable as reference genes in LPS stimulated anterior kidney leucocytes 72 hrs post stimulation [47]. However, it should be noted that in this study only one time point was sampled, thus the variation seen over time in this study cannot be directly compared. In three-spined stickleback (Gasterosteus aculeatus), hardly any of the 9 genes tested in anterior kidney cell cultures, including ACTB1, EF1A1, and HPRT1, were expressed with SD lower than 1 [54]. UbcE having a SD of 0.5 in non-stimulated control cells (Table 3) could be a good candidate for normalization of anterior kidney leucocytes in other experimental setups not using ConA-PMA as stimuli.

Conclusion

This study reports the expression level of six commonly used reference genes during halibut development, in different tissues of healthy and NNV-infected halibut, and in anterior kidney leucocytes both in vivo and in vitro stimulated. The study emphasizes the need for such pilot studies, as universal reference genes have not been identified. Generally, it was found that EF1A1 and RPL7 were the genes that showed the least variation, with HPRT1 and UbcE as intermediate genes, and ACTB1 and Tubb2C as the least stable ones. However, EF1A1 and RPL7 cannot be recommended for a two-gene normalization strategy, as both genes are involved in protein biosynthesis and some co-regulation may be expected. During development, UbcE was found to be the most stable reference gene during activation of zygotic transcription. But, some stage specific expression is seen, and it can be questioned if the use of maternal and zygotic transcripts as references within the same experiment is optimal. EF1A1 and RPL7 were found to be the best candidate reference genes for a normalization approach in halibut larvae and juveniles, and in tissue samples from one year old halibut. This is the first report exploring reference gene expression during a NNV infection in halibut, and amongst the genes tested especially RPL7 was shown to be a good candidate, but also EF1A1 and HPRT1. Based on this data set, none of the six housekeeping genes were optimal as reference gene candidates in ConA-PMA stimulated leucocytes, but UbcE could be a good candidate in other experimental setups showing relatively stable expression in the non-stimulated control groups. This study will facilitate further work on developmental processes and will help in studying immune responses during host-pathogen interactions in Atlantic halibut, both in vivo and in vitro.

Authors' contributions

ACØ was involved in the design and sampling of adult fish and NNV-infection experiment, responsible for the sequencing, real time RT-PCR, data analysis, and drafted the manuscript. AHN participated as a supervisor in the study design and analysis, and helped in drafting the manuscript. SP was involved in study design and sampling in all the studies, and helped in drafting the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was financially supported by The Research Council of Norway. The authors thank the staff at Austevoll Aquaculture Reseach Station for halibut rearing and maintenance, especially Svanild Gokstad and Sissel Waage Kalvenes for assistance during sampling of the halibut larvae and juveniles. Thanks to Rolf Hetlelid Olsen, Ina Nepstad, and Elin Sørhus for excellent technical assistance. And thanks to Rasmus Skern for the combined Trizol and RNAeasy method, and Craig Morton for critical reading.

References

  1. Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays.

    J Mol Endocrino 2000, 25(2):169-193. Publisher Full Text OpenURL

  2. Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR--a perspective.

    J Mol Endocrino 2005, 34(3):597-601. Publisher Full Text OpenURL

  3. 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.0031-0034.0011. BioMed Central Full Text OpenURL

  4. Harboe T, Tuene S, Mangor-Jensen A, Rabben H, Huse I: Design and operation of an incubator for yolk-sac larvae of Atlantic halibut.

    Prog Fish-Cult 1994, 56(3):188-193. Publisher Full Text OpenURL

  5. Harboe T, Mangor-Jensen A: Time of first feeding of Atlantic halibut, Hippoglossus hippoglossus L., larvae.

    Aquaculture Res 1998, 29(12):913-918. Publisher Full Text OpenURL

  6. Bergh O, Nilsen F, Samuelsen OB: Diseases, prophylaxis and treatment of the Atlantic halibut Hippoglossus hippoglossus: a review.

    Dis Aquat Organ 2001, 48(1):57-74. PubMed Abstract | Publisher Full Text OpenURL

  7. Grotmol S, Bergh O, Totland GK: Transmission of viral encephalopathy and retinopathy (VER) to yolk-sac larvae of the Atlantic halibut Hippoglossus hippoglossus: occurrence of nodavirus in various organs and a possible route of infection.

    Dis Aquat Organ 1999, 36(2):95-106. PubMed Abstract | Publisher Full Text OpenURL

  8. Johansen R, Grove S, Svendsen AK, Modahl I, Dannevig B: A sequential study of pathological findings in Atlantic halibut, Hippoglossus hippoglossus (L), throughout one year after an acute outbreak of viral encephalopathy and retinopathy.

    J Fish Dis 2004, 27(6):327-341. PubMed Abstract | Publisher Full Text OpenURL

  9. Sommerset I, Lorenzen E, Lorenzen N, Bleie H, Nerland AH: A DNA vaccine directed against a rainbow trout rhabdovirus induces early protection against a nodavirus challenge in turbot.

    Vaccine 2003, 21(32):4661-4667. PubMed Abstract | Publisher Full Text OpenURL

  10. Sommerset I, Nerland AH: Complete sequence of RNA1 and subgenomic RNA3 of Atlantic halibut nodavirus (AHNV).

    Dis Aquat Organ 2004, 58(2-3):117-125. PubMed Abstract | Publisher Full Text OpenURL

  11. Sommerset I, Skern R, Biering E, Bleie H, Fiksdal IU, Grove S, Nerland AH: Protection against Atlantic halibut nodavirus in turbot is induced by recombinant capsid protein vaccination but not following DNA vaccination.

    Fish Shellfish Immunol 2005, 18(1):13-29. PubMed Abstract | Publisher Full Text OpenURL

  12. Hildahl J, Galay-Burgos M, Sweeney G, Einarsdóttir IE, Björnsson BT: Identification of two isoforms of Atlantic halibut insulin-like growth factor-I receptor genes and quantitative gene expression during metamorphosis.

    Comp Biochem Physiol B Biochem Mol Biol 2007, 147(3):395-401. PubMed Abstract | Publisher Full Text OpenURL

  13. Matsuoka MP, van Nes S, Andersen O, Benfey TJ, Reith M: Real-time PCR analysis of ovary- and brain-type aromatase gene expression during Atlantic halibut (Hippoglossus hippoglossus) development.

    Comp Biochem Physiol B Biochem Mol Biol 2006, 144(1):128-135. PubMed Abstract | Publisher Full Text OpenURL

  14. Patel S, Overgard AC, Nerland AH: CD8α and CD8β in Atlantic Hippoglossus hippoglossus: Cloning, characterization and gene expression during viral and bacterial infection.

    Fish Shellfish Immunol 2008, 25(5):570-580. PubMed Abstract | Publisher Full Text OpenURL

  15. van Nes S, Andersen O: Temperature effects on sex determination and ontogenetic gene expression of the aromatases cyp19a and cyp19b, and the estrogen receptors esr1 and esr2 in Atlantic halibut (Hippoglossus hippoglossus).

    Mol Reprod Dev 2006, 73(12):1481-1490. PubMed Abstract | Publisher Full Text OpenURL

  16. 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. PubMed Abstract | Publisher Full Text OpenURL

  17. Infante C, Matsuoka MP, Asensio E, Canavate JP, Reith M, Manchado M: Selection of housekeeping genes for gene expression studies in larvae from flatfish using real-time PCR.

    BMC Mol Biol 2008, 9(1):28. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Zhong Q, Zhang Q, Wang Z, Qi J, Chen Y, Li S, Sun Y, Li C, Lan X: Expression Profiling and Validation of Potential Reference Genes During Paralichthys olivaceus Embryogenesis.

    Mar Biotechnol 2008, 3:310-318. Publisher Full Text OpenURL

  19. Hagen O, Fernandes JM, Solberg C, Johnston IA: Expression of growth-related genes in muscle during fasting and refeeding of juvenile Atlantic halibut, Hippoglossus hippoglossus L.

    Comp Biochem Physiol B Biochem Mol Biol 2009, 152(1):47-53. PubMed Abstract | Publisher Full Text OpenURL

  20. Matsuoka MP, Infante C, Reith M, Canavate JP, Douglas SE, Manchado M: Translational machinery of senegalese sole (Solea senegalensis Kaup) and Atlantic halibut (Hippoglossus hippoglossus L.): comparative sequence analysis of the complete set of 60s ribosomal proteins and their expression.

    Mar Biotechnol 2008, 10(6):676-691. PubMed Abstract | Publisher Full Text OpenURL

  21. Zhang Z, Hu J: Development and validation of endogenous reference genes for expression profiling of medaka (Oryzias latipes) exposed to endocrine disrupting chemicals by quantitative real-time RT-PCR.

    Toxicol Sci 2007, 95(2):356-368. PubMed Abstract | Publisher Full Text OpenURL

  22. Jørgensen SM, Kleveland EJ, Grimholt U, Gjøen T: Validation of reference genes for real-time polymerase chain reaction studies in Atlantic salmon.

    Mar Biotechnol 2006, 8(4):398-408. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Harboe T, Mangor-Jensen A, Naas KE, Næss T: A tank design for first feeding of Atlantic halibut, Hippoglossus hippoglossus L., larvae.

    Aquaculture Res 1998, 29(12):919-923. Publisher Full Text OpenURL

  25. Mangor-Jensen A, Harboe T, Hennoe JS, Troland R: Design and operation of Atlantic halibut, Hippoglossus hippoglossus L., egg incubators.

    Aquaculture Res 1998, 29(12):887-892. Publisher Full Text OpenURL

  26. Pittman K, Bergh Ø, Opstad I, Skiftesvik AB, Skjolddal L, Strand H: Development of egg and yolk sac larvae of halibut (Hippoglossus hippoglossus L.).

    J Appl Ichthyol 1990, 6(3):142-160. Publisher Full Text OpenURL

  27. Patel S, Sorhus E, Fiksdal IU, Espedal PG, Bergh O, Rodseth OM, Morton HC, Nerland AH: Ontogeny of lymphoid organs and development of IgM-bearing cells in Atlantic halibut (Hippoglossus hippoglossus L.).

    Fish Shellfish Immunol 2009, 26(3):385-395. PubMed Abstract | Publisher Full Text OpenURL

  28. Grotmol S, Nerland AH, Biering E, Totland GK, Nishizawa T: Characterization of the capsid protein gene from a nodavirus strain affecting the Atlantic halibut Hippoglossus hippoglossus and design of an optimal reverse-transcriptase polymerase chain reaction (RT-PCR) detection assay.

    Dis Aquat Organ 2000, 39(2):79-88. PubMed Abstract | Publisher Full Text OpenURL

  29. Dannevig BH, Nilsen R, Modahl I, Jankowska M, Taksdal T, Press CM: Isolation in cell culture of nodavirus from farmed Atlantic halibut Hippoglossus hippoglossus in Norway.

    Dis Aquat Organ 2000, 43(3):183-189. PubMed Abstract | Publisher Full Text OpenURL

  30. Patel S, Malde K, Lanzen A, Olsen RH, Nerland AH: Identification of immune related genes in Atlantic halibut (Hippoglossus hippoglossus L.) following in vivo antigenic and in vitro mitogenic stimulation.

    Fish Shellfish Immunol 2009, 27(6):729-738. PubMed Abstract | Publisher Full Text OpenURL

  31. Park KC, Osborne JA, Tsoi SC, Brown LL, Johnson SC: Expressed sequence tags analysis of Atlantic halibut (Hippoglossus hippoglossus) liver, kidney and spleen tissues following vaccination against Vibrio anguillarum and Aeromonas salmonicida.

    Fish Shellfish Immunol 2005, 18(5):393-415. PubMed Abstract | Publisher Full Text OpenURL

  32. Douglas SE, Knickle LC, Kimball J, Reith ME: Comprehensive EST analysis of Atlantic halibut (Hippoglossus hippoglossus), a commercially relevant aquaculture species.

    BMC genomics 2007, 8:144. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  33. Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: a laboratory manual. 2nd edition. Cold Spring Harbor Laboratory Press; 1987.

  34. 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. PubMed Abstract | Publisher Full Text OpenURL

  35. Korsnes K, Devold M, Nerland AH, Nylund A: Viral encephalopathy and retinopathy (VER) in Atlantic salmon Salmo salar after intraperitoneal challenge with a nodavirus from Atlantic halibut Hippoglossus hippoglossus.

    Dis Aquat Organ 2005, 68(1):7-15. PubMed Abstract | Publisher Full Text OpenURL

  36. 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. PubMed Abstract | Publisher Full Text OpenURL

  37. 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. PubMed Abstract | Publisher Full Text OpenURL

  38. Fleige S, Pfaffl MW: RNA integrity and the effect on the real-time qRT-PCR performance.

    Mol Aspects Med 2006, 27(2-3):126-139. PubMed Abstract | Publisher Full Text OpenURL

  39. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt Method.

    Methods 2001, 25(4):402-408. PubMed Abstract | Publisher Full Text OpenURL

  40. Pelegri F: Maternal factors in zebrafish development.

    Dev Dyn 2003, 228(3):535-554. PubMed Abstract | Publisher Full Text OpenURL

  41. Kane DA, Kimmel CB: The zebrafish midblastula transition.

    Development 1993, 119(2):447-456. PubMed Abstract | Publisher Full Text OpenURL

  42. Mitter K, Kotoulas G, Magoulas A, Mulero V, Sepulcre P, Figueras A, Novoa B, Sarropoulou E: Evaluation of candidate reference genes for QPCR during ontogenesis and of immune-relevant tissues of European seabass (Dicentrarchus labrax).

    Comp Biochem Physiol B Biochem Mol Biol 2009, 153(4):340-347. PubMed Abstract | Publisher Full Text OpenURL

  43. Tang R, Dodd A, Lai D, McNabb WC, Love DR: Validation of zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization.

    Acta Biochim Biophys Sin 2007, 39(5):384-390. PubMed Abstract | Publisher Full Text OpenURL

  44. Mishra RK, Eliceiri GL: Three small nucleolar RNAs that are involved in ribosomal RNA precursor processing.

    Proc Natl Acad Sci 1997, 94(10):4972-4977. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Hildahl J, Sweeney G, Galay-Burgos M, Einarsdóttir IE, Björnsson BT: Cloning of Atlantic halibut growth hormone receptor genes and quantitative gene expression during metamorphosis.

    Gen Comp Endocrinol 2007, 151(2):143-152. PubMed Abstract | Publisher Full Text OpenURL

  46. Solanas M, Moral R, Escrich E: Unsuitability of using ribosomal RNA as loading control for Northern blot analyses related to the imbalance between messenger and ribosomal RNA content in rat mammary tumors.

    Anal Biochem 2001, 288(1):99-102. PubMed Abstract | Publisher Full Text OpenURL

  47. Ingerslev HC, Pettersen EF, Jakobsen RA, Petersen CB, Wergeland HI: Expression profiling and validation of reference gene candidates in immune relevant tissues and cells from Atlantic salmon (Salmo salar L.).

    Mol immunol 2006, 43(8):1194-1201. PubMed Abstract | Publisher Full Text OpenURL

  48. Secombes CJ: Lymphokine-release from rainbow trout leucocytes stimulated with concanavalin A. Effects upon macrophage spreading and adherence.

    Dev Comp Immunol 1987, 11(3):513-520. PubMed Abstract | Publisher Full Text OpenURL

  49. Lin GL, Ellsaesser CF, Clem LW, Miller NW: Phorbol ester/calcium ionophore activate fish leukocytes and induce long-term cultures.

    Dev Comp Immunol 1992, 16(2-3):153-163. PubMed Abstract | Publisher Full Text OpenURL

  50. Xu D, Woodfield SE, Lee TV, Fan Y, Antonio C, Bergmann A: Genetic control of programmed cell death (apoptosis) in Drosophila.

    Fly 2009, 3(1):78-90. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Krumschnabel G, Podrabsky JE: Fish as model systems for the study of vertebrate apoptosis.

    Apoptosis 2009, 14(1):1-21. PubMed Abstract | Publisher Full Text OpenURL

  52. Nam BH, Hirono I, Aoki T: Bulk isolation of immune response-related genes by expressed sequenced tags of Japanese flounder Paralichthys olivaceus leucocytes stimulated with Con A/PMA.

    Fish Shellfish Immunol 2003, 14(5):467-476. PubMed Abstract | Publisher Full Text OpenURL

  53. Arma NR, Hirono I, Aoki T: Characterization of expressed genes in kidney cells of Japanese flounder Paralichthys olivaceus following treatment with ConA/PMA and LPS.

    Fish Pathology 2004, 39(4):189-196. OpenURL

  54. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL