Email updates

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

Open Access Research article

Expression profiles of the pluripotency marker gene POU5F1 and validation of reference genes in rabbit oocytes and preimplantation stage embryos

Solomon Mamo14*, Arpad Baji Gal1, Zsuzsanna Polgar2 and Andras Dinnyes13

Author Affiliations

1 Genetic Reprogramming Group, Agricultural Biotechnology Centre, Szent Gyorgyi Albert u. 4, H-2100 Gödöllő, Hungary

2 Faculty of Natural Sciences, Constantine the Philosopher University, Slovakia

3 Molecular Animal Biotechnology Laboratory, Szent Istvan University, Pater K. u. 1, H-2103 Gödöllő, Hungary

4 University College of Dublin, Lyons Research Farm, Newcastle Co. Dublin, Ireland

For all author emails, please log on.

BMC Molecular Biology 2008, 9:67  doi:10.1186/1471-2199-9-67

The electronic version of this article is the complete one and can be found online at:

Received:22 April 2008
Accepted:28 July 2008
Published:28 July 2008

© 2008 Mamo et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



The surge in the number of gene expression studies and tendencies to increase the quality of analysis have necessitated the identification of stable reference genes. Although rabbits are classical experimental model animals, stable reference genes have not been identified for normalization. The aims of this study were to compare the expression profiles of the widely used reference genes in rabbit oocytes and preimplantation stage embryos, and to select and validate stable ones to use as reference.


Quantitative real time PCR method was used to evaluate 13 commonly used references (Actb, Gapdh, Hprt1, H2afz, Ubc, Ppia, Eef1e1, Polr2a, Tbp, G6pdx, B2m, Pgk1, and Ywhaz) and POU5F1 (Oct4) genes. Expressions of these genes were examined in multiple individual embryos of seven different preimplantation developmental stages and embryo types (in vivo and in vitro). Initial analysis identified three genes (Ubc, Tbp, and B2m) close to the detection limit with irregular expression between the different stages. As variability impedes the selection of stable genes, these were excluded from further analysis. The expression levels of the remaining ten genes, varied according to developmental stage and embryo types. These genes were ranked using the geNorm software and finally the three most stable references (H2afz, Hprt1, and Ywhaz) were selected. Normalization factor was calculated (from the geometric averages of the three selected genes) and used to normalize the expressions of POU5F1 gene. The results showed the expected expression patterns of the POU5F1 during development.


Compared to the earlier studies with similar objectives, the comparison of large number of genes, the use of multiple individual embryos as compared to pools, and simultaneous analyses of in vitro and in vivo derived embryo samples were unique approaches in our study. Based on quantification, pattern and geNorm analyses, we found the three genes (H2afz, Hprt1, and Ywhaz) to be the most stable across developmental stages and embryo types, and the geometric averages of these genes can be used for appropriate normalization.


Preimplantation embryo development is a dynamic developmental process recognized by changes in the transcript type and quantity [1-3], cell number, total and poly (A) RNA contents [4,5]. These entire phenomena are known to be the basis for the changes in the shape, physiology and functions of the embryos leading to compaction, differentiation, implantation and further development. During this preimplantation period, however, about 15–50% of zygotes die, largely as a result of unknown factors [6]. Therefore, it is imperative to understand the dynamics and search for factors contributing to the losses. Moreover, elucidation of preimplantation development is critical for the management of infertility and refinement of assisted reproductive technology [7]. Gene expression profiles would yield insight into the complex molecular pathways controlling early development. However, the use of classical techniques [8-10] for embryo gene expression analyses was constrained by technical limitations and a dearth of starting materials. To overcome the problem, the production of embryos in large quantities is restricted by associated ethical and cost factors. We (and others) have described in detail the problems associated with the use of common RNA detection and analysis techniques for the application of preimplantation stage embryo studies [11,12].

Real time PCR is a sensitive quantitative method of choice that overcomes the potential sensitivity problems in the earlier classical analysis techniques [13], and detects more subtle changes in gene expression [14]. Its sensitivity allows working with a minimal amount of starting material, while still achieving an accurate quantification of poorly transcribed mRNAs [15]. It is a fast and reliable technique, provided care is taken in all the procedures [16]. The precautions associated with the use of this technique were described elsewhere [12] and, unless properly addressed, could lead to severe misinterpretation of the results [17,18]. To address the issue and account for variations, different normalization strategies were used, and the details were well reviewed [19]. Some researchers use exogenously added references in the form of synthetic RNA or globin RNAs [20,21]. However, this approach has been challenged for competing with endogenous sequences for enzyme and nucleotides during PCR, its demand for extra procedures and associated cost [9,19]. As a result, it has not been adopted widely. The use of endogenous reference genes (commonly known as housekeeping genes) to normalize the expression of gene(s) of interest is the most commonly used approach. Housekeeping genes are constitutively expressed to maintain cellular function, and they are presumed to produce the minimally essential transcripts necessary for normal physiology [22,23]. Despite their ease of use and the wider adoption of this approach, inconsistency in the type and number of reference genes used has made cross-study comparisons difficult (even for closely related studies) and interpretation of the results questionable. As a result, this area has been a hot topic and the focus of many journals publication in recent times.

Rabbits are classical experimental models with close similarity in developmental biology to large animals and humans [24]. They are also preferentially used in pulmonary, cardiovascular, metabolic disorder studies, and for antibody production and drug screening [25]. As an experimental animal for various models, different gene expression studies were carried out using rabbit embryos. Most of these studies [26,27] used a single normalizer gene that differed between studies. The uses of some of these genes (Actb, Gapdh, and ribosomal RNAs) were evaluated in different studies and criticized, due to the observed variations between treatment groups [28-31]. We have earlier published a similar mouse study [12]. However, stable genes have not been identified for use in rabbit preimplantation stage embryo gene expression studies. The aims of the current study were to compare the expression profiles of 13 widely used reference genes in rabbit oocytes and preimplantation stage embryos that were derived in vivo and produced in vitro, to select the most stable ones as reference and finally to validate them by using to normalize the expression of POU5F1 gene. The gene POU5F1 belongs to the POU domain family transcription factors that regulate the transcription of their target genes [32,33]. The rabbit POU5F1 gene was recently cloned [34], but its expression during rabbit embryo preimplantation stage has not been well characterized. Thus, we analyzed the expression of this gene in the oocyte and in vitro produced preimplantation stage rabbit embryos. To our knowledge, this is the first study in rabbit to compare the wide selection of reference genes both in the in vitro and in vivo derived embryo samples with a final validation of the selected genes.


Sequence analyses and product confirmation

For the 14 genes used in this study, the sequences were referred from the databases (Table 1), and primers were designed from these sequences. After optimizing PCR conditions and amplification, the products were sequenced for confirmation. Based on the sequence analysis, all primers amplified the expected amplicon sizes (Table 2).

Table 1. Reference genes selected for the study

Table 2. Primer sequences and optimum PCR conditions for amplifying the products

To amplify beta actin (Actb) gene, we designed primers from 3 different sequences (AY598932, AF404278 and X60733), optimized and compared them. However, the current primer (design from X60733 sequence) gave the best amplification results. It has also the best beta actin similarity hits with orthologs. Although the source sequence was identified as gamma non-muscle actin in the NCBI, and as actin, alpha skeletal in the Ensembl databanks, we assumed this to be a wrong labelling. The feedback we got from the NCBI staff also supports our view. Therefore we used it to design primer and further comparison.

Quality control and primer screening

Quality of the analyzed samples plays a major role in the correct interpretation of the results. Based on the morphological observations, in both in vivo derived and in vitro produced embryo samples, presumably good quality oocytes and preimplantation stage embryos were collected. In the subsequent procedures, in addition to our established quality RNA isolation procedures, the minus RT (reverse transcriptase) reaction during cDNA synthesis, and the design of most primers at the exon-exon junction enabled us to control the absence of contaminating genomic DNA. During the initial screening assays, similar cDNA dilutions from the pooled embryos were used and reaction conditions were optimised for each primer pair separately. Optimum primer concentrations were selected based on the absence of dimer and signal detection recognized by earlier CT values. The CT is defined as the number of cycles needed for the fluorescence to reach a specific threshold level of detection and is inversely correlated with the amount of template nucleic acid present in the reaction [35]. The majority of the selected candidate reference genes were detected in most preimplantation stages, but with various signal intensities. Using a similar low concentration template for all, the ten references (Eef1e1, Polr2a, Ywhaz, Ppia, H2afz, Hprt1, Pgk1, G6pdx, Gapdh, and Actb) and POU5f1 genes were detected at CT values below 33. The other three genes (B2m, Ubc, and Tbp) were detected at CT values above 33 with irregularities in amplification (at some of the stages where the transcripts of that particular gene were at lower concentration). As embryo materials are scarce at the preimplantation stages, we excluded the three genes detected with the latest CT values. Then, the ten candidate references and POU5F1 genes with the earlier signals (earlier CT values) were selected for further comparisons.

Standards and PCR efficiency analyses

For all the selected candidate genes, melt curve analyses were performed at the end of PCR reactions. The specificity and integrity of the PCR products were confirmed by the presence of a single peak (See additional files 1, 2 and 3 as representatives). For the selected genes, the standard curves were deduced from four-fold serial dilutions of the five-pooled embryo cDNA preparations measured at five points. To ensure the comparability of PCR assays, three serial dilutions were made independently that enabled us to determine the CT values and PCR efficiencies of the individual assay and calculate the correlation between them. The assays for the selected candidate reference genes were characterized by a linear correlation coefficient (R2) that varied from 0.963 to 0.997 (average 0.986) and PCR efficiencies between 92.3% and 110.7% (average 97.5%), where as the assay for POU5F1 gene has 82.4% efficiency and 99.8% correlation (Table 3). Based on these results, the assays can be trusted and valid for the quantification of transcripts and further comparisons.

Additional file 1. Supplementary figure 1. Representative melt curve and cycle threshold analyses of the gene Ywhaz.

Format: DOC Size: 66KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 2. Supplementary figure 2. Representative melt curve and cycle threshold analyses of the gene H2afz.

Format: DOC Size: 47KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 3. Supplementary figure 3. Representative melt curve and cycle threshold analyses of the gene Hprt.

Format: DOC Size: 64KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Table 3. Standard curve parameters for the candidate reference genes

Reference genes profile analyses at different developmental stages

In order to select the best stable genes for normalization, expressions of the selected ten candidate reference genes were compared in different preimplantation stage embryos. For gene transcript quantification, five individual embryo cDNA preparations per stage were used in identical experimental procedures but optimised conditions for each gene amplification.

The ten candidate genes were detected in all developmental stages examined with various signal intensities. The oocyte stage was taken as a calibrator for all the genes and the relative expression levels at different developmental stages are shown in Figures 1 and 2.

thumbnailFigure 1. Relative expression levels of selected candidate reference genes in the in vivo derived rabbit embryos. The expression level at the oocyte stage was taken as a reference to calculate the relative levels of the other stages. Stages with different letters are significantly (P ≤ 0.05) different for the levels of the gene.

thumbnailFigure 2. Relative expression levels of selected candidate reference genes in the in vitro produced rabbit embryos. The expression level at the oocyte stage was taken as a reference to calculate the relative levels of the other stages. Stages with different letters are significantly (P ≤ 0.05) different for the levels of the gene. An asterisk (*) indicates significant differences {(* = P ≤ 0.05) and (** = P ≤ 0.01)} with the in vivo levels of the same stage.

Compared to the oocyte stage, the relative levels of most genes (Eef1e1, Ppia, H2afz, Hprt, Pgk1 and G6pdx) reduced at the 1-cell stage, while the relative levels of the remaining genes (Polr2a, Ywhaz, Gapdh, and Actb) increased in various intensities at this stage. The relative levels of some genes (Eef1e1, Polr2a, Ppia, Gapdh, Actb (in vitro) were transiently reduced to the minimum levels at the 8-cell stages, and increased thereafter. Despite differences in the embryo type (in vivo or in vitro), the above expression patterns remained similar (Figure 1 and Figure 2). In both in vivo derived and in vitro produced embryos, the relative levels of H2afz and Hprt1 genes showed a persistent increase in levels after fertilization. The gene Ywhaz (except 2-cell in vivo) also has the same profile.

Comparative analyses in the in vivo derived and in vitro produced embryos

To examine the effects of embryo type (in vivo and in vitro) on gene expression stability, expression profiles of the selected candidate genes were compared in the in vivo derived and in vitro produced embryo samples.

Major pattern differences were not observed due to change in the embryo type (in vivo or in vitro) (Figures 1 and 2). However, the stage-by-stage comparisons revealed differential transcript levels between the two embryo types. At the 2-cell stage, in vivo produced embryos showed higher transcript levels for majority of the genes (except Pgk1 and Actb) compared to their in vitro counterparts, and differences were significant (p < 0.05) for some genes (Eef1e1, Ywhaz, Ppia and G6pdx). However, at the 8-cell stage, the in vitro samples showed higher transcript levels for majority of the genes, with significant difference (p < 0.05) for some (Eef1e1, Ywhaz, Gapdh and Actb). At the morula and blastocyst stages, the relative transcript levels of in vivo derived embryo samples were higher for most genes compared to the in vitro samples. These differences were significant (p < 0.01 for Ywhaz and Ppia, and p < 0.05 for H2afz, Gapdh and Actb) at the morula stage and (p < 0.05 for Polr2a and Hprt1) at the blastocyst stage (Figure 2).

Gene expression stability analyses

The profiles of candidate genes were analysed and the expression stability measure values (M) were calculated using the geNorm software [36]. Following the procedures of the software, the least stable genes were determined (by higher M values) and continuously excluded to recalculate the M values for the rest. Finally all genes were ranked based on the calculated M values and the three most stable genes (with the lowest M values) were selected. Accordingly, for the in vivo derived samples, the genes H2afz and Hprt1 were found to be the most stable, followed by the genes Ywhaz, Actb, Ppia, G6pdx, Gapdh, Eef1e1, Polr2a and Pgk1 in their order of appearance. Similar analyses for the in vitro samples selected the genes H2afz and Ywhaz as the most stable genes, followed by Hprt1, G6pdx, Ppia, Gapdh, Eef1e1, Polr2a, Actb and Pgk1 in their order of appearance. Generally, differences in the embryo type (in vivo or in vitro) had a minor rearrangement effect on the stability order for majority of the genes examined. The gene Actb showed wider stability range between the embryo types with better stability in the in vivo samples compared to the in vitro (Figure 3).

thumbnailFigure 3. Average gene expression stability values of the candidate reference genes in the in vivo derived and in vitro produced embryo samples as calculated by geNorm software and ranking made based on the relative stability values.

Expression profiles of POU5F1 gene in rabbit preimplantation stage embryos

To evaluate the performances of the newly selected reference genes, the normalization factor (calculated from the geometric averages of the three selected genes) was calculated. The expression of POU5F1 gene was quantified in the in vitro produced rabbit oocytes and embryo samples and the results were normalized by the newly calculated factor. The expression at the oocyte stage was used as a calibrator to calculate the relative expression levels in the different developmental stages (Figure 4). Based on analyses, the POU5F1 gene was expressed in all developmental stages with various signal intensities. Comparatively higher expression levels were observed at the oocyte and zygote stages, and declined gradually to the lowest levels at the 8-cell stage. However, starting from the morula stage, the levels increased continuously to the blastocyst stage.

thumbnailFigure 4. Relative expression levels of the POU5F1 gene in the in vitro produced preimplantation stage rabbit embryos. The expression at the oocyte stage was taken as a reference to calculate the relative levels of the other stages. Stages with different letters are significantly (P ≤ 0.05) different for the levels of the gene.


The number of studies revealing the importance of normalization has increased recently and, as a result, the search for study-specific appropriate reference genes has gained momentum. This can be seen in the rising number of relevant publications. In this study, we selected a large number of reference genes and compared their expression in the rabbit oocytes and different preimplantation stage embryos that were derived in vivo or produced in vitro. After detailed analyses, here we show the selection of the three most stable reference genes in rabbit oocytes and embryos that can be used for normalization.

Rabbits have many advantages [24,27], including a convenient reproductive pattern (non-seasonal breeding, induced ovulation and short gestation period of 31 days), the possibility of keeping them in conventional housing in an indoor facility, the lower cost of procurement and handling compared to large animals. Moreover, they have sizeable milk production that allows their use as test animals for therapeutic protein expression in milk. A combination of all these factors has led to the use of rabbits as experimental models and their use in gene expression studies. Moreover, the studies of genetic reprogramming in rabbit embryos [37-41,27], in which our laboratory is also taking part, are surging.

As gene expression results depend on the type of normalizer gene used [42,17], finding appropriate reference genes is timely, and will help toward the correct interpretation of the experimental results. Finding appropriate reference genes implies careful selection of stable genes evaluated for expression stability. Differences arising from the quality and quantity of input RNA, efficiencies of reverse transcription reaction and handling errors can be accounted by normalizing the expression of a target gene to the reference. This allows the direct comparison of normalized expression values between samples [43].

Earlier suggestions to use RNA mass quantity and adding exogenous templates for normalization has been challenged for the inherent technical problems hindering its wider adoption [36,9,19,12]. Moreover, adding exogenous template will compete with endogenous sequences for primers and nucleotides during the PCR reaction [9]. The use of endogenous reference genes for normalization is the widely used approach in most applications. It has gained acceptance over time due to its biological relevance, and its consideration of various errors during the process. Moreover, compared to adding exogenous templates to the samples, it has no extra procedures for its application [19]. To address this issue, the use of individual embryos compared to the pooled ones was in line with earlier findings [44-46] that indicated the significance of using individual samples. It can enable accurate statistical analyses, identify more biological variations and subtle gene expression changes.

Most of the genes (Eef1e1, Polr2a, Ppia, Gapdh and Actb in vitro) had the respective lowest levels at the 8-cell stage. This transcript depression at the 8-cell stage coincides with the lag in development (developmental block) occurring during the maternal to zygote transition (MZT) in rabbit embryos [47-49]. The morphological and molecular variations between embryos of different mammalian species have been discussed [3]. The effects of culture and treatment conditions on the rabbit preimplantation stage embryos have also been described earlier [50-52]. In this study, the comparative transcript levels of the same gene in different embryo types (in vivo vs. in vitro) were described. The profile shows transcript level variations and substantiates the fact that in vitro conditions, in general, are sub-optimal and influence gene expression levels. In our earlier mouse study [12], we have indicated the effects of embryo type on the stability and selection of reference genes. As far as we know, such comparisons and selections were not made for rabbit gene expression studies. Therefore, we made comparison of large number of commonly used reference genes and selected the most stable ones for normalization. Thus, this study is timely and the recommendations derived from the study can be widely applicable for rabbit embryo studies elsewhere, as most rabbit embryo studies deal with both or either of the in vivo/in vitro models.

Despite the traditional ways of using a single reference gene for normalization, the approach has been frequently criticized in a number of studies [28,30,43,17]. The implications of using an inappropriate reference gene have already been discussed [17]. Another elegant study [36] also demonstrated the error related to using a single reference gene and proposed to calculate a normalization factor based on the geometric averages of at least three carefully selected stable reference genes. In this study, detailed comparisons of the expression levels in different developmental stages and embryo types (in vivo vs. in vitro) indicated the comparative stable expression of the genes H2afz, Ywhaz and Hprt1. This was also further confirmed by the geNorm expression stability analyses (Figure 3). The selected constitutive genes were stable in both in vivo and in vitro conditions with slight variations in the order of stability (Figure 3). This variation might be attributed to the differences in expression levels with change in embryo type. Our earlier mouse study [12] has also indicated a shift in the stability order with the changes in embryo type.

In the current study, variations of gene stability values between the in vivo and in vitro comparisons were narrower than in the earlier mouse study. This may indicate better in vitro culture conditions of rabbit embryos than the in vitro mouse embryo culture. Moreover, the comparative stability of H2afz and Hprt genes were also shown in our previous mouse study. Although the gene Ppia was a preferred reference gene in the earlier studies [12,53], it is not among the best in this study. This might be due to species differences or the presence of other better performing genes. In earlier studies, some genes (including Actb and Gapdh) were assumed as universal reference genes, without further evaluation, and used individually for normalization of gene expression data [19]. The current result in rabbit embryos and our earlier mouse study [12] clearly indicated the inappropriateness of some of the widely used reference genes for normalization, at least under the conditions examined, and reinforce the recommendation to evaluate study-specific reference genes before using them for normalization.

To evaluate the performances of the newly selected reference genes, the expression values of POU5F1 gene was normalized by the calculated new normalization factor. After normalization, the observed pattern closely resembled the profiles of POU5F1 orthologue genes. Since POU5F1 orthologue genes, including human, bovine and porcine, are highly conserved, a similar role has been suggested for POU5F1 in all mammals [54]. In line with our finding of the POU5F1 profile, earlier studies in human [55] and bovine [56] embryos detected the expression of POU5F1 throughout all preimplantation stages. The decline of POU5F1 levels during the MZT stage and continuous increase after compaction was also observed in the earlier mouse [57] and human [54] embryo studies. Therefore, the normalized POU5F1 gene profile revealed the conserved pattern and further confirmed the suitability of the selected reference genes.


Our study, for the first time, revealed a detailed reference gene validation and selection for rabbit preimplantation stage embryo studies. The outcomes indicated the possibility of using the same selected genes for both in vivo derived and in vitro produced embryo gene expression studies. Although transcript level variations were observed between individual embryos analyzed for the same gene, the expression patterns were almost similar. Based on detailed analyses of the results, including the pattern and ranking with geNorm, the genes H2afz, Ywhaz and Hprt1 were found to be the most stable. We also believe that the number of reference genes used for normalization depends on several factors [36]. However, using the geometric averages of the above three genes is preferred for accurate gene expression results in rabbit oocytes and preimplantation stage embryo studies. The appropriateness of these genes irrespective of the embryo type, and the conserved patterns of the POU5F1 gene after normalization, further confirm the suitability of these genes.


All chemicals, unless stated otherwise, were purchased from Sigma-Aldrich Chemical Inc. (St. Louis, USA).

The animal experiments were executed in full compliance with European and Hungarian laws and regulations, and were approved by the Agricultural Biotechnology Center, Gödöllö, Animal Experimentation Committee.

Oocyte collection

Hycole hybrid female rabbits were induced to superovulate by administration of 120 IU pregnant mare serum gonadotropin (PMSG, Folligon® Intervet, The Netherlands) i.m. and 72 hours later, these animals were injected i.v. with 170 IU human chorionic gonadotropin (hCG, Choragon® Richter Gedeon Rt., Hungary). Donor females (does) were slaughtered at 16 hours post hCG administration and the oocytes were flushed from the ampullae of the oviducts with M2 medium. The cumulus layer was then removed using 0.1% hyaluronidase in M2 medium, and the denuded metaphase 2 stage oocytes (confirmed by the presence of a single polar body) were collected individually for RNA isolation.

Zygote collection

Hycole hybrid female rabbits were induced to superovulate by administration of PMSG, and 72 hrs later, and shortly before mating, with hCG hormones as described earlier. Each injected female was mated with a male of proven fertility (buck) from the same breed. Donor females were slaughtered at 20 hours post hCG administration and the zygotes were flushed from the oviduct with M2 medium. The flushed zygotes were evaluated morphologically under the microscope and those with two pronuclei, two polar bodies and compact cytoplasm were selected for the experiment.

In vitro embryo production

The in vitro embryos were produced by further culturing the zygotes until the required developmental stages. For this, the zygotes were cultured in 50-μl EBSS (Earle's Balanced Salt Solution) complete drops [56] under mineral oil and incubated at 38.5°C in humidified atmosphere containing 5% CO2 in air. The embryos were cultured until the required developmental stages, and five embryos each at the 2-cell (26 hrs), 8-cell early (44 hrs), 8-cell late (54 hrs), morula (68 hrs), and blastocyst (103 hrs) stages were collected individually.

In vivo embryo production

The procedures of superovulation and mating were as described earlier. For collection of in vivo derived embryos, donor females were slaughtered at respective times (hours post hCG administration) and the embryos at the 2-cell (26 hrs), 8-cell early (41 hrs), 8-cell late (47 hrs), morula (62 hrs), and blastocyst (98 hrs) stages were collected individually.

During sample collection, the denuded oocytes and different developmental stage embryos were washed three times with and collected individually in 2-μl RNase-free water for storage at -80°C until RNA extraction.

RNA isolation and cDNA synthesis

The procedures of RNA isolation and cDNA synthesis were as described in detail earlier [15]. Briefly, messenger RNA was extracted individually from five embryos per developmental stage and embryo type (in vivo or in vitro), using Dynabeads® mRNA DIRECT™ Micro Kit (Dynal A.S, Oslo, Norway), following the manufacturer's instructions. The individually frozen embryos were lysed and incubated with pre-washed magnetic Dynabeads that can base pair with poly (A) tails of mRNA molecules. After hybridisation and subsequent repeated washes with buffers, the RNA was eluted in RNase-free water and reverse transcribed into cDNA, using M-MLV RT kit (Invitrogen, Carlsbad, USA) in a final 20-μl reaction volume. Minus RT reactions were performed to check the absence of contaminating residual DNA. The cDNA synthesis reactions were carried at 42°C for 1 hour, followed by heat inactivation of the enzyme at 75°C for 10 minutes.

Primer design and sequence analyses

A total of 14 genes, most commonly used and considered as reference for normalization by the earlier studies, were selected for evaluation throughout the different embryo developmental stages and embryo types (in vivo and in vitro). Moreover, POUF1 gene was also used for evaluating the selected reference genes. Primers were designed and optimised prior to initial screening and quantitation experiments. Comparison and screening of primers were carried out using optimised protocol for each primer. Template cDNA from the same source was used, and finally those with earlier signals were selected for further quantification. For these genes (Table 1), the expected sizes of the products were confirmed by gel electrophoresis on a 2% agarose gel. The PCR products were also cloned (TOPO TA cloning kit, Invitrogen), and sequenced for confirmation.

Real time PCR reaction conditions and analyses

The details of real time PCR conditions for each gene were as described in Table 2. During quantification of the transcripts, the assay for each gene consisted of five replicates per stage, both in vivo derived and in vitro produced oocytes and preimplantation stage embryo samples, negative and positive controls. All genes were compared from the same stock to avoid inter-assay template variations, and all quantifications were performed consecutively without interruption. Each sample in a run consisted of 0.08 embryo equivalent cDNA template, 200–500 nM of each primer (Table 2), and 50% SYBR® Green JumpStart™ Taq ReadyMix™ in 15-μl reaction volume. The reaction conditions were template denaturation and polymerase activation at 95°C for 2 min followed by 45 cycles of 95°C denaturation for 15 sec, 57°C to 68°C annealing and extension for 45 sec. All reactions were carried out using the Rotor-Gene™ 3000 real time PCR machine (Corbett Research, Mortlake, Australia), and the results were analysed with the integrated Rotor-Gene software (version 6.1). At the end of PCR reactions, melt curve analyses were performed for all the genes. For calculating PCR efficiencies, standard curves were generated from assays made with four fold serial dilutions of 5-pooled blastocyst cDNA preparations. To ensure the compatibility of PCR assays, three independent serial dilutions were made that enabled us to determine the CT values and PCR efficiencies of the individual assay and calculate the correlation between them. PCR efficiencies (E) were calculated with the equation

E = (10 [-1/slope]-1) × 100.

geNorm and expression stability analyses

Analyses of the gene expression stability over the different embryonic stages and types (in vivo vs. in vitro) were performed using the geNorm software [36]. The analyses relies on the principle that the expression ratio of two ideal internal control genes is identical in all samples, regardless of the experimental condition or cell type, and determined as the standard deviation of the logarithmically transformed expression ratios [36]. The internal control gene stability measure value (M) was calculated as the average pair-wise variation of a particular gene with respect to the rest of the genes, and ranking was made based on these values. The lower the M value, the more stable the expression of the gene under consideration. The most stable reference genes were identified by stepwise exclusions of the least stable gene and recalculating the M values with the rest genes.

Authors' contributions

SM conceived the experiment, performed the experimental design, all the molecular biology analyses, including experiment execution, data analyses and interpretation, and was the primary author of the manuscript. ABG participated in the experimental design, performed primer design, data analyses, statistical tests and participated in manuscript preparation. ZP prepared in vivo and in vitro embryos for the experiment. AD supervised the study design, execution, analysis, and approved the final version. All authors read and approved the manuscript.


We are grateful to Mrs Hajnalka Csákány Tolnainé for her excellent technical assistance during the experiment. We also extend our gratitude to the reviewers, particularly reviewer 1 for the comment on beta-actin gene. The work was supported by Wellcome Trust (Grant No.070246), EU FP6 (MEXT-CT-2003-509582), EU FP6 (LSHG-CT-2006-518240), EU FP6 (MRTN-CT-2006-035468).


  1. Rothstein JL, Johnson D, DeLoia JA, Skowronski J, Solter D, Knowles B: Gene expression during preimplantation mouse development.

    Genes Dev 1992, 6:1190-1201. PubMed Abstract | Publisher Full Text OpenURL

  2. Goossens K, Poucke MV, Van Soom A, Vandesompele J, Zeveren AV, Peelman LJ: Selection of reference genes for quantitative real-time PCR in bovine preimplantation embryos.

    BMC Dev Biol 2005, 5:27. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Kuijk EW, du Puy L, van Tol HT, Haagsman HP, Colenbrander B, Roelen BAJ: Validation of reference genes for quantitative real-time PCR studies in porcine oocytes and preimplantation embryos.

    BMC Dev Biol 2007, 7:58. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Piko L, Clegg KB: Quantitative changes in total RNA, total Poly (A), and ribosomes in early mouse embryos.

    Dev Biol 1982, 89:362-378. PubMed Abstract | Publisher Full Text OpenURL

  5. Bachvarova R, De Leon V, Johnson A, Kaplan G, Paynton BV: Changes in total RNA, polyadenylated RNA, and actin mRNA during meiotic maturation of mouse oocytes.

    Dev Biol 1985, 108:325-331. PubMed Abstract | Publisher Full Text OpenURL

  6. Warner CM, Exley GE, McElhinny AS, Tang C: Genetic regulation of preimplantation mouse embryo survival.

    J Exp Zool 1998, 282:272-279. PubMed Abstract | Publisher Full Text OpenURL

  7. Steuerwald N, Cohen J, Herrera RJ, Brenner CA: Analysis of gene expression in single oocytes and embryos by real-time rapid cycle fluorescence monitored RT-PCR.

    Mol Hum Reprod 1999, 5:1034-1039. PubMed Abstract | Publisher Full Text OpenURL

  8. Jeong YJ, Choi HW, Shin HS, Cui XS, Kim NH, Gerton GL, Jun JH: Optimization of real time RT-PCR methods for the analysis of gene expression in mouse eggs and preimplantation embryos.

    Mol Rep Dev 2005, 71:284-289. PubMed Abstract | Publisher Full Text OpenURL

  9. Robert C, McGraw S, Massicotte L, Pravetoni M, Gandolfi F, Sirard M-A: Quantification of housekeeping transcript levels during the development of bovine preimplantation embryos.

    Biol Reprod 2002, 67:1465-1472. PubMed Abstract | Publisher Full Text OpenURL

  10. Bettegowda A, Patel O, Ireland JJ, Smith GW: Quantitative analysis of messenger RNA abundance for ribosomal protein L-15, cyclophilin-A, Phosphoglycerokinase, β-glucuronidase, glyceraldehydes 3-phosphate dehydrogenase. β-actin, and histone H2A during bovine oocyte maturation and early embryogenesis in vitro.

    Mol Reprod Dev 2006, 73:267-278. PubMed Abstract | Publisher Full Text OpenURL

  11. Mamo S, Sargent CA, Affara NA, Tesfaye D, El-Halawany N, Wimmers K, Gilles M, Schellander K, Ponsuksili S: Transcript profiles of some developmentally important genes detected in bovine oocytes and in vitro produced blastocysts using RNA amplification and cDNA microarrays.

    Reprod Dom Anim 2006, 41:527-534. PubMed Abstract | Publisher Full Text OpenURL

  12. 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 7:14. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Heid CA, Stevens J, Livak KJ, Williams PM: Real time quantitative PCR.

    Genome Res 1996, 6:986-994. PubMed Abstract | Publisher Full Text OpenURL

  14. Nailis H, Coenye T, Nieuwerburgh FV, Deforce D, Nelis HJ: Development and evaluation of different normalization strategies for gene expression studies in Candida albicans biofilms by real-time PCR.

    BMC Mol Biol 2006, 7:25. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Spinsanti G, Panti C, Lazzeri E, Marsili L, Casino S, Frati F, Fossi CM: Selection of reference genes for quantitative RT-PCR studies in striped dolphin (Stenella coeruleoalba) skin biopsies.

    BMC Mol Biol 2006, 7:32. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

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

  17. Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GAW, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization.

    Anal Biochem 2005, 344:141-143. PubMed Abstract | Publisher Full Text OpenURL

  18. Zhang QJ, Chadderton A, Clark RL, Augustine-Rauch KA: Selection of normalizer genes in conducting relative gene expression analysis of embryos.

    Birth Defects Res 2003, 67:533-544. PubMed Abstract | Publisher Full Text OpenURL

  19. Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalization; strategies and considerations.

    Genes Immun 2005, 6:279-284. PubMed Abstract | Publisher Full Text OpenURL

  20. De Sousa PA, Watson AJ, Schultz RM: Transient expression of translation initiation factor is conservatively associated with embryonic gene activation in murine and bovine embryos.

    Biol Reprod 1998, 59:969-977. PubMed Abstract | Publisher Full Text OpenURL

  21. McGraw S, Morin G, Vigneault C, Leclerc P, Sirard M-A: Investigation of MYST4 histone acetyl transferase and its involvement in mammalian gametogenesis.

    BMC Dev Biol 2007, 7:123. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Butte AJ, Dzau VJ, Glueck SB: Further defining housekeeping, or "maintenance," genes: Focus on "a compendium of gene expression in normal human tissues".

    Physiol Genomics 2001, 7:95-96. PubMed Abstract OpenURL

  23. Willems E, Mateizel I, Kemp C, Cauffman G, Sermon K, Leyns L: Selection of reference genes in mouse embryos and in differentiating human and mouse ES cells.

    Int J Dev Biol 2006, 50:627-635. PubMed Abstract | Publisher Full Text OpenURL

  24. Dinnyes A, Tian XC, Yang X: Cloning of rabbits. In Principles of cloning. Edited by Cibelli J, Lanza RP, Campbell KHS, West MD. California: Academic Press, Elsevier Science; 2002:343-366. OpenURL

  25. Graves KH, Moreadith RW: Derivation and characterization of putative pluripotential embryonic stem cells from preimplantation rabbit embryos.

    Mol Reprod Dev 1993, 36:424-433. PubMed Abstract OpenURL

  26. Sobajima S, Shimer AL, Chadderdon RC, Kompel JF, Kim JS, Gilbertson LG, Kang JD: Quantitative analysis of gene expression in a rabbit model of intervertebral disc degeneration by real-time polymerase chain reaction.

    Spine J 2005, 5:14-23. PubMed Abstract | Publisher Full Text OpenURL

  27. Fang ZF, Gai H, Huang YZ, Li SG, Chen XJ, Shi JJ, Wu L, Liu A, Xu P, Sheng HZ: Rabbit embryonic stem cell lines derived from fertilized, parthenogenetic or somatic cell nuclear transfer embryos.

    J Exp Cell Res 2006, 312:3669-3682. OpenURL

  28. 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. PubMed Abstract OpenURL

  29. Suzuki T, Higgins PJ, Crawford DR: Control selection for RNA quantitation.

    BioTechniques 2000, 29:332-337. PubMed Abstract OpenURL

  30. Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies.

    Anal Biochem 2002, 309:293-300. PubMed Abstract | Publisher Full Text OpenURL

  31. Boonkusol D, Gal AB, Bodo S, Gorhony B, Kitiyanant Y, Dinnyes A: Gene expression profiles and in vitro development following vitrification of pronuclear and 8-cell stage mouse embryos.

    Mol Reprod Dev 2006, 73:700-708. PubMed Abstract | Publisher Full Text OpenURL

  32. Ovitt CE, Schöler HR: The molecular biology of Oct-4 in the early mouse embryo.

    Mol Hum Reprod 1998, 4:1021-1031. PubMed Abstract | Publisher Full Text OpenURL

  33. Hansis C, Grifo JA, Krey LC: Oct-4 expression in inner cell mass and trophectoderm of human blastocysts.

    Mol Hum Reprod 2000, 6:999-1004. PubMed Abstract | Publisher Full Text OpenURL

  34. Shi JJ, Cai DH, Chen XJ, Sheng HZ: Cloning and characterization of rabbit POU5f1 gene.

    DNA Seq 2007, in press. PubMed Abstract | Publisher Full Text OpenURL

  35. Walker NJ: Tech Sight. A technique whose time has come.

    Science 2002, 296:557-559. PubMed Abstract | Publisher Full Text OpenURL

  36. 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. [] webcite

    Genome Biol 2002, 3:research 0034. OpenURL

  37. Stice SL, Robl JM: reprogramming in nuclear transplant rabbit embryos.

    Biol Reprod 1988, 39:657-664. PubMed Abstract | Publisher Full Text OpenURL

  38. Dinnyes A, Dai Y, Barber M, Liu L, Xu J, Zhou P, Yang X: Development of cloned embryos from adult rabbit fibroblasts: effect of activation treatment and donor cell preparation.

    Biol Reprod 2001, 64:257-263. PubMed Abstract | Publisher Full Text OpenURL

  39. Chesne P, Adenot PG, Viglietta C, Baratte M, Boulanger L, Renard JP: Cloned rabbits produced by nuclear transfer from adult somatic cells.

    Nature Biotechnol 2002, 20:366-369. PubMed Abstract | Publisher Full Text OpenURL

  40. Li WY, Yu WD, Qi B, He X, Ge Wang NY, Liu GS, Du M, Chen QX: Analysis of gene expression in rabbit nuclear transfer embryos: Use of single-embryo mRNA differential display.

    Develop Growth Differ 2003, 45:543-551. OpenURL

  41. Li S, Chen X, Fang Z, Shi J, Sheng HZ: Rabbits generated from fibroblasts through nuclear transfer.

    Reproduction 2006, 131:1085-1090. PubMed Abstract | Publisher Full Text OpenURL

  42. Glare EM, Divjak M, Bailey MJ, Walters EH: β-actin and Gapdh housekeeping gene expression in asthmatic airways is variable and not suitable for normalizing mRNA levels.

    Thorax 2002, 57:765-770. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Bower NI, Moser RJ, Hill JR, Lehnert SA: Universal reference method for real-time PCR gene expression analysis of preimplantation embryos.

    BioTechniques 2007, 42:199-206. PubMed Abstract OpenURL

  44. Jolly RA, Goldstein KM, Wei T, Gao H, Chen P, Huang S, Colet JM, Ryan TP, Thomas CE, Estrem ST: Pooling samples within microarray studies: a comparative analysis of rat liver transcription response to prototypical toxicants.

    Physiol Genomics 2005, 22:346-355. PubMed Abstract | Publisher Full Text OpenURL

  45. Glass A, Henning J, Karopka T, Scheel T, Bansemer S, Koczan D, Gierl L, Rolfs A, Gimsa U: Representation of individual gene expression in completely pooled mRNA samples.

    Biosci Biotechnol Biochem 2005, 69:1098-1103. PubMed Abstract | Publisher Full Text OpenURL

  46. Craig DW, Huentelman MJ, Hu-Lince D, Zismann VL, Kruer MC, Lee AM, Puffenberger EG, Pearson JM, Stephan DA: Identification of disease causing loci using an array-based genotyping approach on pooled DNA.

    BMC Genomics 2005, 6:138. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Manes C: Nucleic acid synthesis in preimplantation rabbit embryos. III. A "dark period" immediately following fertilization, and the early predominance low molecular weight RNA synthesis. Quantitative aspects, relationship to early morphogenesis and protein synthesis.

    J Exp Zool 1969, 172:303-310. PubMed Abstract OpenURL

  48. Manes C: The participation of the embryonic genome during early cleavage in the rabbit.

    Dev Biol 1973, 32:453-459. PubMed Abstract | Publisher Full Text OpenURL

  49. Brunet-Simon A, Henrion G, Renard JP, Duranthon V: Onset of zygotic transcription and maternal transcript legacy in the rabbit embryos.

    Mol Reprod Dev 2001, 58:127-136. PubMed Abstract | Publisher Full Text OpenURL

  50. Makarevich AV, Olexikoval L, Chrenek P, Kubovicova E, Freharova K, Pivko J: The effect of hyperthermia in vitro on stability of rabbit preimplantation embryos.

    Physiol Res, in press.

    2006; Nov 6


  51. Higa K, Shimmura S, Kato N, Kawakita T, Miyashita H, Itabashi Y, Fakuda K, Shimazaki J, Tsubota K: Proliferation and differentiation of transplantable rabbit epithelian sheets engineered with or without an amniotic membrane carrier.

    Invest Opthalmol Vis Sci 2007, 48:597-604. OpenURL

  52. Pierzynski P, Gajda B, Smorag Z, Rasmussen AD, Kuczynski W: Effect of atosiban on rabbit embryo development and human sperm motility.

    Fertil Steril 2007, 87:1147-1152. PubMed Abstract | Publisher Full Text OpenURL

  53. Jung M, Ramankulov A, Roigas J, Johannsen M, Ringsdorf M, Kristiansen G, Jung K: In search of suitable reference genes for gene expression studies of human renal cell carcinoma by real-time PCR.

    BMC Mol Biol 2007, 8:47. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Cauffman G, Velde H, Liebaers I, Steirteghem AV: Oct-4 mRNA and protein expression during human preimplantation development.

    Mol Hum Reprod 2005, 11:173-181. PubMed Abstract | Publisher Full Text OpenURL

  55. Abdel-Rahman B, Fiddler M, Rappolee D, Pergament E: Expression of transcription regulating genes in human preimplantation embryos.

    Hum Reprod 1995, 10:2787-2792. PubMed Abstract | Publisher Full Text OpenURL

  56. van Eijk MJT, van Rooijen MA, Modina S, Scesi L, Folkers G, van Tol HTA, Bevers MM, Fisher SR, Lewin HA, Rakacolli D, Galli C, de Vaureix C, Trounson AO, Mummery CL, Gandolfi F: Molecular cloning, genetic mapping and developmental expression of bovine POU5F1.

    Biol Reprod 1999, 60:1093-1103. PubMed Abstract | Publisher Full Text OpenURL

  57. Pesce M, Wang X, Wolgemuth DJ, Schöler H: Differential expression of the Oct-4 transcription factor during mouse germ cell differentiation.

    Mech Dev 1998, 71:89-98. PubMed Abstract OpenURL

  58. Mitalipov SM, White KL, Farrar VR, Morrey J, Reed WA: Development of nuclear transfer and parthenogenetic rabbit embryos activated with inositol 1, 4, 5-Triphosphate.

    Biol Reprod 1999, 60:821-827. PubMed Abstract | Publisher Full Text OpenURL

  59. Almeida A, Thiery JP, Magdelénat H, Radvanyi F: Gene expression analysis by real-time reverse transcription polymerase chain reaction: influence of tissue handling.

    Anal Biochem 2004, 328:101-108. PubMed Abstract | Publisher Full Text OpenURL