<title>Background</title> <p>Preimplantation embryo development is a dynamic developmental process recognized by changes in the transcript type and quantity [##REF##1628826##1##, ####REF##16324220##2##, ##REF##17540017##3####17540017##3##], cell number, total and poly (A) RNA contents [##REF##6173273##4##,##REF##2416609##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 [##REF##9723184##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 [##REF##10541565##7##]. Gene expression profiles would yield insight into the complex molecular pathways controlling early development. However, the use of classical techniques [##UREF##0##8##, ####REF##12390877##9##, ##REF##16261607##10####16261607##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 [##REF##17107512##11##,##REF##17341302##12##].</p> <p>Real time PCR is a sensitive quantitative method of choice that overcomes the potential sensitivity problems in the earlier classical analysis techniques [##REF##8908518##13##], and detects more subtle changes in gene expression [##REF##16889665##14##]. Its sensitivity allows working with a minimal amount of starting material, while still achieving an accurate quantification of poorly transcribed mRNAs [##REF##16984641##15##]. It is a fast and reliable technique, provided care is taken in all the procedures [##REF##14706621##16##]. The precautions associated with the use of this technique were described elsewhere [##REF##17341302##12##] and, unless properly addressed, could lead to severe misinterpretation of the results [##REF##16054107##17##,##REF##14632301##18##]. To address the issue and account for variations, different normalization strategies were used, and the details were well reviewed [##REF##15815687##19##]. Some researchers use exogenously added references in the form of synthetic RNA or globin RNAs [##REF##9746750##20##,##REF##17980037##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 [##REF##12390877##9##,##REF##15815687##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 [##REF##11773595##22##,##REF##16892176##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.</p> <p>Rabbits are classical experimental models with close similarity in developmental biology to large animals and humans [##UREF##1##24##]. They are also preferentially used in pulmonary, cardiovascular, metabolic disorder studies, and for antibody production and drug screening [##REF##8305204##25##]. As an experimental animal for various models, different gene expression studies were carried out using rabbit embryos. Most of these studies [##REF##15653081##26##,##UREF##2##27##] used a single normalizer gene that differed between studies. The uses of some of these genes (<italic>Actb</italic>, <italic>Gapdh</italic>, and ribosomal RNAs) were evaluated in different studies and criticized, due to the observed variations between treatment groups [##REF##10617337##28##, ####REF##10948434##29##, ##REF##12413463##30##, ##REF##16541460##31####16541460##31##]. We have earlier published a similar mouse study [##REF##17341302##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 <italic>in vivo </italic>and produced <italic>in vitro</italic>, to select the most stable ones as reference and finally to validate them by using to normalize the expression of <italic>POU5F1 </italic>gene. The gene <italic>POU5F1 </italic>belongs to the POU domain family transcription factors that regulate the transcription of their target genes [##REF##9835353##32##,##REF##11044462##33##]. The rabbit <italic>POU5F1 </italic>gene was recently cloned [##REF##17852331##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 <italic>in vitro </italic>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 <italic>in vitro </italic>and <italic>in vivo </italic>derived embryo samples with a final validation of the selected genes.</p>
<title>Methods</title> <p>All chemicals, unless stated otherwise, were purchased from Sigma-Aldrich Chemical Inc. (St. Louis, USA).</p> <p>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.</p> <title>Oocyte collection</title> <p>Hycole hybrid female rabbits were induced to superovulate by administration of 120 IU pregnant mare serum gonadotropin (PMSG, Folligon<sup>® </sup>Intervet, The Netherlands) <italic>i.m</italic>. and 72 hours later, these animals were injected <italic>i.v</italic>. with 170 IU human chorionic gonadotropin (hCG, Choragon<sup>® </sup>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.</p> <title>Zygote collection</title> <p>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.</p> <title><italic>In vitro </italic>embryo production</title> <p>The <italic>in vitro </italic>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 [##REF##10208969##56##] under mineral oil and incubated at 38.5°C in humidified atmosphere containing 5% CO<sub>2 </sub>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.</p> <title><italic>In vivo </italic>embryo production</title> <p>The procedures of superovulation and mating were as described earlier. For collection of <italic>in vivo </italic>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.</p> <p>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.</p> <title>RNA isolation and cDNA synthesis</title> <p>The procedures of RNA isolation and cDNA synthesis were as described in detail earlier [##REF##16984641##15##]. Briefly, messenger RNA was extracted individually from five embryos per developmental stage and embryo type (<italic>in vivo </italic>or <italic>in vitro</italic>), using Dynabeads<sup>® </sup>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.</p> <title>Primer design and sequence analyses</title> <p>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 (<italic>in vivo </italic>and <italic>in vitro</italic>). Moreover, <italic>POUF1 </italic>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 ##TAB##0##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.</p> <title>Real time PCR reaction conditions and analyses</title> <p>The details of real time PCR conditions for each gene were as described in Table ##TAB##1##2##. During quantification of the transcripts, the assay for each gene consisted of five replicates per stage, both <italic>in vivo </italic>derived and <italic>in vitro </italic>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 ##TAB##1##2##), and 50% SYBR<sup>® </sup>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 C<sub>T </sub>values and PCR efficiencies of the individual assay and calculate the correlation between them. PCR efficiencies (E) were calculated with the equation</p> <p></p> <title>geNorm and expression stability analyses</title> <p>Analyses of the gene expression stability over the different embryonic stages and types (<italic>in vivo </italic>vs. <italic>in vitro</italic>) were performed using the geNorm software [##REF##12184808##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 [##REF##12184808##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.</p>
<title>Results</title> <title>Sequence analyses and product confirmation</title> <p>For the 14 genes used in this study, the sequences were referred from the databases (Table ##TAB##0##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 ##TAB##1##2##).</p> <p>To amplify beta actin (<italic>Actb) </italic>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.</p> <title>Quality control and primer screening</title> <p>Quality of the analyzed samples plays a major role in the correct interpretation of the results. Based on the morphological observations, in both <italic>in vivo </italic>derived and <italic>in vitro </italic>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 C<sub>T </sub>values. The C<sub>T </sub>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 [##REF##11964485##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 (<italic>Eef1e1, Polr2a, Ywhaz, Ppia, H2afz, Hprt1, Pgk1, G6pdx, Gapdh</italic>, and <italic>Actb</italic>) and <italic>POU5f1 </italic>genes were detected at C<sub>T </sub>values below 33. The other three genes (<italic>B2m, Ubc, and Tbp</italic>) were detected at C<sub>T </sub>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 C<sub>T </sub>values. Then, the ten candidate references and <italic>POU5F1 </italic>genes with the earlier signals (earlier C<sub>T </sub>values) were selected for further comparisons.</p> <title>Standards and PCR efficiency analyses</title> <p>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 ##SUPPL##0##1##, ##SUPPL##1##2## and ##SUPPL##2##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 C<sub>T </sub>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 (R<sup>2</sup>) 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 <italic>POU5F1 </italic>gene has 82.4% efficiency and 99.8% correlation (Table ##TAB##2##3##). Based on these results, the assays can be trusted and valid for the quantification of transcripts and further comparisons.</p> <title>Reference genes profile analyses at different developmental stages</title> <p>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.</p> <p>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 ##FIG##0##1## and ##FIG##1##2##.</p> <p>Compared to the oocyte stage, the relative levels of most genes (<italic>Eef1e1, Ppia, H2afz, Hprt, Pgk1 and G6pdx</italic>) reduced at the 1-cell stage, while the relative levels of the remaining genes (<italic>Polr2a, Ywhaz, Gapdh</italic>, and <italic>Actb</italic>) increased in various intensities at this stage. The relative levels of some genes (<italic>Eef1e1, Polr2a, Ppia, Gapdh, Actb (in vitro</italic>) were transiently reduced to the minimum levels at the 8-cell stages, and increased thereafter. Despite differences in the embryo type (<italic>in vivo </italic>or <italic>in vitro</italic>), the above expression patterns remained similar (Figure ##FIG##0##1## and Figure ##FIG##1##2##). In both <italic>in vivo </italic>derived and <italic>in vitro </italic>produced embryos, the relative levels of <italic>H2afz </italic>and <italic>Hprt1 </italic>genes showed a persistent increase in levels after fertilization. The gene <italic>Ywhaz </italic>(except 2-cell in vivo) also has the same profile.</p> <title>Comparative analyses in the <italic>in vivo </italic>derived and <italic>in vitro </italic>produced embryos</title> <p>To examine the effects of embryo type (<italic>in vivo </italic>and <italic>in vitro</italic>) on gene expression stability, expression profiles of the selected candidate genes were compared in the <italic>in vivo </italic>derived and <italic>in vitro </italic>produced embryo samples.</p> <p>Major pattern differences were not observed due to change in the embryo type (<italic>in vivo </italic>or <italic>in vitro) </italic>(Figures ##FIG##0##1## and ##FIG##1##2##). However, the stage-by-stage comparisons revealed differential transcript levels between the two embryo types. At the 2-cell stage, <italic>in vivo </italic>produced embryos showed higher transcript levels for majority of the genes (except <italic>Pgk1 </italic>and <italic>Actb</italic>) compared to their <italic>in vitro </italic>counterparts, and differences were significant (p &lt; 0.05) for some genes (<italic>Eef1e1, Ywhaz, Ppia </italic>and <italic>G6pdx</italic>). However, at the 8-cell stage, the <italic>in vitro </italic>samples showed higher transcript levels for majority of the genes, with significant difference (p &lt; 0.05) for some (<italic>Eef1e1, Ywhaz</italic>, <italic>Gapdh </italic>and <italic>Actb</italic>). At the morula and blastocyst stages, the relative transcript levels of <italic>in vivo </italic>derived embryo samples were higher for most genes compared to the <italic>in vitro </italic>samples. These differences were significant (p &lt; 0.01 for <italic>Ywhaz </italic>and <italic>Ppia</italic>, and p &lt; 0.05 for <italic>H2afz</italic>, <italic>Gapdh </italic>and <italic>Actb</italic>) at the morula stage and (p &lt; 0.05 for <italic>Polr2a </italic>and <italic>Hprt1</italic>) at the blastocyst stage (Figure ##FIG##1##2##).</p> <title>Gene expression stability analyses</title> <p>The profiles of candidate genes were analysed and the expression stability measure values (M) were calculated using the geNorm software [##REF##12184808##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 <italic>in vivo </italic>derived samples, the genes <italic>H2afz </italic>and <italic>Hprt1 </italic>were found to be the most stable, followed by the genes <italic>Ywhaz, Actb, Ppia, G6pdx, Gapdh, Eef1e1</italic>, <italic>Polr2a </italic>and <italic>Pgk1 </italic>in their order of appearance. Similar analyses for the <italic>in vitro </italic>samples selected the genes <italic>H2afz </italic>and <italic>Ywhaz </italic>as the most stable genes, followed by <italic>Hprt1, G6pdx, Ppia, Gapdh, Eef1e1, Polr2a, Actb </italic>and <italic>Pgk1 </italic>in their order of appearance. Generally, differences in the embryo type (<italic>in vivo </italic>or <italic>in vitro</italic>) had a minor rearrangement effect on the stability order for majority of the genes examined. The gene <italic>Actb </italic>showed wider stability range between the embryo types with better stability in the <italic>in vivo </italic>samples compared to the <italic>in vitro </italic>(Figure ##FIG##2##3##).</p> <title>Expression profiles of <italic>POU5F1 </italic>gene in rabbit preimplantation stage embryos</title> <p>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 <italic>POU5F1 </italic>gene was quantified in the <italic>in vitro </italic>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 ##FIG##3##4##). Based on analyses, the <italic>POU5F1 </italic>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.</p>
<title>Discussion</title> <p>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 <italic>in vivo </italic>or produced <italic>in vitro</italic>. 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.</p> <p>Rabbits have many advantages [##UREF##1##24##,##UREF##2##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 [##REF##3196797##37##, ####REF##11133682##38##, ##REF##11923842##39##, ##UREF##3##40##, ##REF##16735547##41####16735547##41##,##UREF##2##27##], in which our laboratory is also taking part, are surging.</p> <p>As gene expression results depend on the type of normalizer gene used [##REF##12200519##42##,##REF##16054107##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 [##REF##17373485##43##].</p> <p>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 [##REF##12184808##36##,##REF##12390877##9##,##REF##15815687##19##,##REF##17341302##12##]. Moreover, adding exogenous template will compete with endogenous sequences for primers and nucleotides during the PCR reaction [##REF##12390877##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 [##REF##15815687##19##]. To address this issue, the use of individual embryos compared to the pooled ones was in line with earlier findings [##REF##15914576##44##, ####REF##15973040##45##, ##REF##16197552##46####16197552##46##] that indicated the significance of using individual samples. It can enable accurate statistical analyses, identify more biological variations and subtle gene expression changes.</p> <p>Most of the genes (<italic>Eef1e1, Polr2a, Ppia, Gapdh </italic>and <italic>Actb in vitro) </italic>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 [##REF##5395775##47##, ####REF##4789702##48##, ##REF##11139224##49####11139224##49##]. The morphological and molecular variations between embryos of different mammalian species have been discussed [##REF##17540017##3##]. The effects of culture and treatment conditions on the rabbit preimplantation stage embryos have also been described earlier [##UREF##4##50##, ####UREF##5##51##, ##REF##17224148##52####17224148##52##]. In this study, the comparative transcript levels of the same gene in different embryo types (<italic>in vivo </italic>vs. <italic>in vitro</italic>) were described. The profile shows transcript level variations and substantiates the fact that <italic>in vitro </italic>conditions, in general, are sub-optimal and influence gene expression levels. In our earlier mouse study [##REF##17341302##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 <italic>in vivo/in vitro </italic>models.</p> <p>Despite the traditional ways of using a single reference gene for normalization, the approach has been frequently criticized in a number of studies [##REF##10617337##28##,##REF##12413463##30##,##REF##17373485##43##,##REF##16054107##17##]. The implications of using an inappropriate reference gene have already been discussed [##REF##16054107##17##]. Another elegant study [##REF##12184808##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 (<italic>in vivo </italic>vs. <italic>in vitro</italic>) indicated the comparative stable expression of the genes <italic>H2afz, Ywhaz </italic>and <italic>Hprt1</italic>. This was also further confirmed by the geNorm expression stability analyses (Figure ##FIG##2##3##). The selected constitutive genes were stable in both <italic>in vivo </italic>and <italic>in vitro </italic>conditions with slight variations in the order of stability (Figure ##FIG##2##3##). This variation might be attributed to the differences in expression levels with change in embryo type. Our earlier mouse study [##REF##17341302##12##] has also indicated a shift in the stability order with the changes in embryo type.</p> <p>In the current study, variations of gene stability values between the <italic>in vivo </italic>and <italic>in vitro </italic>comparisons were narrower than in the earlier mouse study. This may indicate better <italic>in vitro </italic>culture conditions of rabbit embryos than the <italic>in vitro </italic>mouse embryo culture. Moreover, the comparative stability of <italic>H2afz </italic>and <italic>Hprt </italic>genes were also shown in our previous mouse study. Although the gene <italic>Ppia </italic>was a preferred reference gene in the earlier studies [##REF##17341302##12##,##REF##17559644##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 <italic>Actb </italic>and <italic>Gapdh</italic>) were assumed as universal reference genes, without further evaluation, and used individually for normalization of gene expression data [##REF##15815687##19##]. The current result in rabbit embryos and our earlier mouse study [##REF##17341302##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.</p> <p>To evaluate the performances of the newly selected reference genes, the expression values of <italic>POU5F1 </italic>gene was normalized by the calculated new normalization factor. After normalization, the observed pattern closely resembled the profiles of <italic>POU5F1 </italic>orthologue genes. Since <italic>POU5F1 </italic>orthologue genes, including human, bovine and porcine, are highly conserved, a similar role has been suggested for <italic>POU5F1 </italic>in all mammals [##REF##15695770##54##]. In line with our finding of the <italic>POU5F1 </italic>profile, earlier studies in human [##REF##8567814##55##] and bovine [##REF##10208969##56##] embryos detected the expression of <italic>POU5F1 </italic>throughout all preimplantation stages. The decline of <italic>POU5F1 </italic>levels during the MZT stage and continuous increase after compaction was also observed in the earlier mouse [##REF##9507072##57##] and human [##REF##15695770##54##] embryo studies. Therefore, the normalized <italic>POU5F1 </italic>gene profile revealed the conserved pattern and further confirmed the suitability of the selected reference genes.</p>
<title>Conclusion</title> <p>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 <italic>in vivo </italic>derived and <italic>in vitro </italic>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 <italic>H2afz, Ywhaz </italic>and <italic>Hprt1 </italic>were found to be the most stable. We also believe that the number of reference genes used for normalization depends on several factors [##REF##12184808##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 <italic>POU5F1 </italic>gene after normalization, further confirm the suitability of these genes.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>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.</p> <title>Results</title> <p>Quantitative real time PCR method was used to evaluate 13 commonly used references (<italic>Actb, Gapdh, Hprt1, H2afz, Ubc, Ppia, Eef1e1, Polr2a, Tbp, G6pdx, B2m, Pgk1</italic>, and <italic>Ywhaz</italic>) and <italic>POU5F1 (Oct4) </italic>genes. Expressions of these genes were examined in multiple individual embryos of seven different preimplantation developmental stages and embryo types (<italic>in vivo </italic>and <italic>in vitro</italic>). Initial analysis identified three genes (<italic>Ubc, Tbp</italic>, and <italic>B2m) </italic>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 (<italic>H2afz, Hprt1</italic>, and <italic>Ywhaz</italic>) were selected. Normalization factor was calculated (from the geometric averages of the three selected genes) and used to normalize the expressions of <italic>POU5F1 </italic>gene. The results showed the expected expression patterns of the POU5F1 during development.</p> <title>Conclusion</title> <p>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 <italic>in vitro </italic>and <italic>in vivo </italic>derived embryo samples were unique approaches in our study. Based on quantification, pattern and geNorm analyses, we found the three genes (<italic>H2afz, Hprt1</italic>, and <italic>Ywhaz</italic>) to be the most stable across developmental stages and embryo types, and the geometric averages of these genes can be used for appropriate normalization.</p>
<title>Authors' contributions</title> <p>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 <italic>in vivo </italic>and <italic>in vitro </italic>embryos for the experiment. AD supervised the study design, execution, analysis, and approved the final version. All authors read and approved the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>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).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Relative expression levels of selected candidate reference genes in the <italic>in vivo </italic>derived rabbit embryos. </bold>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.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Relative expression levels of selected candidate reference genes in the <italic>in vitro </italic>produced rabbit embryos.</bold> 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 <italic>in vivo </italic>levels of the same stage.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Average gene expression stability values of the candidate reference genes in the <italic>in vivo </italic>derived and <italic>in vitro </italic>produced embryo samples as calculated by geNorm software and ranking made based on the relative stability values.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Relative expression levels of the <italic>POU5F1 </italic>gene in the <italic>in vitro </italic>produced preimplantation stage rabbit embryos.</bold> 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.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Reference genes selected for the study </p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Symbol</bold></td><td align="left"><bold>Gene name</bold></td><td align="left"><bold>Source sequences</bold></td><td align="left"><bold>Some of the References</bold></td></tr></thead><tbody><tr><td align="left"><italic>Actb</italic></td><td align="left">Actin, beta, cytoplasmic</td><td align="left">X60733</td><td align="left">2, 16, 17, 35</td></tr><tr><td align="left"><italic>Gapdh</italic></td><td align="left">Glyceraldehyde-3-phosphate dehydrogenase</td><td align="left">NM_001082253</td><td align="left">2, 16, 17, 35</td></tr><tr><td align="left"><italic>Hprt1</italic></td><td align="left">Hypoxanthine guanine phosphoribosyl transferase 1</td><td align="left">ENSOCUG00000003186</td><td align="left">2, 16, 17, 20</td></tr><tr><td align="left"><italic>H2afz</italic></td><td align="left">H2A histone family, member Z</td><td align="left">ENSOCUG00000001888</td><td align="left">8, 12</td></tr><tr><td align="left"><italic>Ubc</italic></td><td align="left">Ubiquitin</td><td align="left">ENSOCUG00000001288</td><td align="left">20, 34</td></tr><tr><td align="left"><italic>Ppia</italic></td><td align="left">Peptidylprolyl isomerase A</td><td align="left">ENSOCUG00000016662</td><td align="left">10, 16, 17</td></tr><tr><td align="left"><italic>Eef1e1</italic></td><td align="left">Eukaryotic translation elongation factor 1 epsilon 1</td><td align="left">ENSOCUG00000012587</td><td align="left">57, 12</td></tr><tr><td align="left"><italic>Polr2a</italic></td><td align="left">Polymerase (RNA) II (DNA directed) polypeptide A</td><td align="left">ENSOCUG00000017929</td><td align="left">16</td></tr><tr><td align="left"><italic>Tbp</italic></td><td align="left">TATA box binding protein</td><td align="left">ENSOCUG00000007979</td><td align="left">2, 16, 17, 20</td></tr><tr><td align="left"><italic>G6pdx</italic></td><td align="left">Glucose-6-phosphate dehydrogenase X-linked</td><td align="left">ENSOCUG00000007866</td><td align="left">16, 21</td></tr><tr><td align="left"><italic>B2m</italic></td><td align="left">Beta-2-microglobulin</td><td align="left">ENSOCUG00000017117</td><td align="left">10, 16, 17, 20</td></tr><tr><td align="left"><italic>Pgk1</italic></td><td align="left">Phosphoglycerate kinase 1</td><td align="left">ENSOCUG00000014726</td><td align="left">17, 20</td></tr><tr><td align="left"><italic>Ywhaz</italic></td><td align="left">Tyrosin 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide</td><td align="left">ENSOCUG00000000734</td><td align="left">2, 34</td></tr><tr><td align="left"><italic>POU5F1</italic></td><td align="left">POU domain, class 5, transcription factor 1</td><td align="left">NM_001099957</td><td align="left">32</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Primer sequences and optimum PCR conditions for amplifying the products</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Gene</td><td align="left">Forward (5' to 3')<break/>Reverse (5' to 3')</td><td align="left">Product size (bp)</td><td align="left">Primer Used (nM)</td><td align="left">Annealing Temp (T°)</td></tr></thead><tbody><tr><td align="left"><italic>Actb</italic></td><td align="left">TCCGCCGCCGGCCCACACT<break/>AGTCCTTCTGGCCCATGC</td><td align="center">188</td><td align="center">0.40<break/>0.40</td><td align="center">60</td></tr><tr><td align="left"><italic>Gapdh</italic></td><td align="left">CAAGTTCCACGGCACGGTCA<break/>CTCGGCACCAGCATCACCC</td><td align="center">118</td><td align="center">0.40<break/>0.40</td><td align="center">68</td></tr><tr><td align="left"><italic>Hprt1</italic></td><td align="left">ACGTCGAGGACTTGGAAAGGGTGTT<break/>GGCCTCCCATCTCCTTCATCACATC</td><td align="center">96</td><td align="center">0.40<break/>0.40</td><td align="center">68</td></tr><tr><td align="left"><italic>H2afz</italic></td><td align="left">AGAGCCGGCTGCCAGTTCC<break/>CAGTCGCGCCCACACGTCC</td><td align="center">85</td><td align="center">0.20<break/>0.33</td><td align="center">59</td></tr><tr><td align="left"><italic>Ubc</italic></td><td align="left">GTGACACCATCGAGAAT<break/>GACACCTCCCCTCAGAC</td><td align="center">173</td><td align="center">0.40<break/>0.30</td><td align="center">60</td></tr><tr><td align="left"><italic>Ppia</italic></td><td align="left">TCCAGGGTTTATGTGCCAGGGTG<break/>CGTTTGCCATGGACAAGATGCC</td><td align="center">137</td><td align="center">0.30<break/>0.30</td><td align="center">68</td></tr><tr><td align="left"><italic>Eef1e1</italic></td><td align="left">ACCGCAGAAGAGAAAGCCATAG<break/>AGCGATGTAGCCCATAGTAGAGGA</td><td align="center">190</td><td align="center">0.50<break/>0.50</td><td align="center">68</td></tr><tr><td align="left"><italic>G6pdx</italic></td><td align="left">AGCCCGCCTCCACTGACTCC<break/>ACCACGTTGTCCGCCTGCAC</td><td align="center">88</td><td align="center">0.40<break/>0.30</td><td align="center">60</td></tr><tr><td align="left"><italic>Tbp</italic></td><td align="left">GCTGAATATAATCCCAAGCGGTTTGC<break/>AAATCAGCGCTGTGGTTCGTGGCTCTC</td><td align="center">73</td><td align="center">0.30<break/>0.30</td><td align="center">68</td></tr><tr><td align="left"><italic>Polr2a</italic></td><td align="left">AGACTTCTCGGCCCGCACTG<break/>ACTTGGCGCCTGGGTACTGG</td><td align="center">176</td><td align="center">0.40<break/>0.40</td><td align="center">68</td></tr><tr><td align="left"><italic>B2m</italic></td><td align="left">GCTCCGTCTTGGGCTTG<break/>CGGATGAAACCCAGATACATAG</td><td align="center">152</td><td align="center">0.33<break/>0.33</td><td align="center">60</td></tr><tr><td align="left"><italic>Pgk1</italic></td><td align="left">TGTTGGTCGGGCGAAGCAG<break/>CAGTGTCTCCACCGCCGATG</td><td align="center">149</td><td align="center">0.50<break/>0.20</td><td align="center">60</td></tr><tr><td align="left"><italic>Ywhaz</italic></td><td align="left">GGTCTGGCCCTTAACTTCTCTGTGTTCTA<break/>GCGTGCTGTCTTTGTATGATTCTTCACTT</td><td align="center">142</td><td align="center">0.50<break/>0.50</td><td align="center">68</td></tr><tr><td align="left"><italic>POU5F1</italic></td><td align="left">CGAGTGAGAGGCAACTTGG<break/>CGGTTACAGAACCACACACG</td><td align="center">125</td><td align="center">0.45<break/>0.20</td><td align="center">57</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Standard curve parameters for the candidate reference genes</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Genes</td><td align="left">Slope (m)</td><td align="left">Intercept (b)</td><td align="left">Efficiency (E)</td><td align="left">Correlation (R<sup>2</sup>)</td></tr></thead><tbody><tr><td align="left"><italic>Actb</italic></td><td align="left">-3.37</td><td align="left">21.90</td><td align="left">98.0</td><td align="left">97.9</td></tr><tr><td align="left"><italic>Gapdh</italic></td><td align="left">-3.43</td><td align="left">18.58</td><td align="left">97.7</td><td align="left">99.4</td></tr><tr><td align="left"><italic>Hprt1</italic></td><td align="left">-3.25</td><td align="left">22.56</td><td align="left">103.1</td><td align="left">98.9</td></tr><tr><td align="left"><italic>H2afz</italic></td><td align="left">-3.52</td><td align="left">18.85</td><td align="left">92.3</td><td align="left">99.4</td></tr><tr><td align="left"><italic>Ppia</italic></td><td align="left">-3.51</td><td align="left">16.07</td><td align="left">92.6</td><td align="left">99.7</td></tr><tr><td align="left"><italic>Eef1e1</italic></td><td align="left">-3.48</td><td align="left">20.50</td><td align="left">94.0</td><td align="left">97.2</td></tr><tr><td align="left"><italic>G6pdx</italic></td><td align="left">-3.09</td><td align="left">23.69</td><td align="left">110.7</td><td align="left">99.2</td></tr><tr><td align="left"><italic>Polr2a</italic></td><td align="left">-3.38</td><td align="left">22.10</td><td align="left">97.7</td><td align="left">96.3</td></tr><tr><td align="left"><italic>Pgk1</italic></td><td align="left">-3.51</td><td align="left">16.89</td><td align="left">92.8</td><td align="left">99.4</td></tr><tr><td align="left"><italic>Ywhaz</italic></td><td align="left">-3.42</td><td align="left">20.09</td><td align="left">96.0</td><td align="left">98.9</td></tr><tr><td align="left"><italic>POU5F1</italic></td><td align="left">-3.83</td><td align="left">17.32</td><td align="left">82.4</td><td align="left">99.8</td></tr></tbody></table></table-wrap>
<disp-formula>E = (10 [-1/slope]-1) × 100.</disp-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Supplementary figure 1. Representative melt curve and cycle threshold analyses of the gene <italic>Ywhaz</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Supplementary figure 2. Representative melt curve and cycle threshold analyses of the gene <italic>H2afz</italic>.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Supplementary figure 3. Representative melt curve and cycle threshold analyses of the gene <italic>Hprt</italic>.</p></caption></supplementary-material>
<graphic xlink:href="1471-2199-9-67-1"/> <graphic xlink:href="1471-2199-9-67-2"/> <graphic xlink:href="1471-2199-9-67-3"/> <graphic xlink:href="1471-2199-9-67-4"/>
<media xlink:href="1471-2199-9-67-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2199-9-67-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2199-9-67-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Jeong", "Choi", "Shin", "Cui", "Kim", "Gerton", "Jun"], "given-names": ["YJ", "HW", "HS", "XS", "NH", "GL", "JH"], "article-title": ["Optimization of real time RT-PCR methods for the analysis of gene expression in mouse eggs and preimplantation embryos"], "source": ["Mol Rep Dev"], "year": ["2005"], "volume": ["71"], "fpage": ["284"], "lpage": ["289"]}, {"surname": ["Dinnyes", "Tian", "Yang", "Cibelli J, Lanza RP, Campbell KHS, West MD"], "given-names": ["A", "XC", "X"], "article-title": ["Cloning of rabbits"], "source": ["Principles of cloning"], "year": ["2002"], "publisher-name": ["California: Academic Press, Elsevier Science"], "fpage": ["343"], "lpage": ["366"]}, {"surname": ["Fang", "Gai", "Huang", "Li", "Chen", "Shi", "Wu", "Liu", "Xu", "Sheng"], "given-names": ["ZF", "H", "YZ", "SG", "XJ", "JJ", "L", "A", "P", "HZ"], "article-title": ["Rabbit embryonic stem cell lines derived from fertilized, parthenogenetic or somatic cell nuclear transfer embryos"], "source": ["J Exp Cell Res"], "year": ["2006"], "volume": ["312"], "fpage": ["3669"], "lpage": ["3682"]}, {"surname": ["Li", "Yu", "Qi", "He", "Ge Wang", "Liu", "Du", "Chen"], "given-names": ["WY", "WD", "B", "X", "NY", "GS", "M", "QX"], "article-title": ["Analysis of gene expression in rabbit nuclear transfer embryos: Use of single-embryo mRNA differential display"], "source": ["Develop Growth Differ"], "year": ["2003"], "volume": ["45"], "fpage": ["543"], "lpage": ["551"]}, {"surname": ["Makarevich", "Olexikoval", "Chrenek", "Kubovicova", "Freharova", "Pivko"], "given-names": ["AV", "L", "P", "E", "K", "J"], "article-title": ["The effect of hyperthermia in vitro on stability of rabbit preimplantation embryos"], "source": ["Physiol Res"], "comment": ["2006; Nov 6"]}, {"surname": ["Higa", "Shimmura", "Kato", "Kawakita", "Miyashita", "Itabashi", "Fakuda", "Shimazaki", "Tsubota"], "given-names": ["K", "S", "N", "T", "H", "Y", "K", "J", "K"], "article-title": ["Proliferation and differentiation of transplantable rabbit epithelian sheets engineered with or without an amniotic membrane carrier"], "source": ["Invest Opthalmol Vis Sci"], "year": ["2007"], "volume": ["48"], "fpage": ["597"], "lpage": ["604"]}]
BMC Mol Biol. 2008 Jul 28; 9:67