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

Importance of Windows of Exposure to Maternal High-Fat Diet and Feto-Placental Effects: Discrimination Between Pre-conception and Gestational Periods in a Rabbit Model

Context and Aim: Lipid overnutrition in female rabbits, from prepuberty, leads to impaired metabolism (dyslipidemia and increased adiposity) and follicular atresia, and, when continued during gestation, affects offspring phenotype with intrauterine growth retardation (IUGR) and leads to placental and lipid metabolism abnormalities. Growth retardation is already observed in embryo stage, indicating a possible implication of periconceptional exposure. The objective of this study was to discriminate the effects of preconception and gestational exposures on feto-placental development.

Materials and Methods: Rabbit 1-day zygotes were collected from female donors under control (CD) or high-fat-high-cholesterol (HD) diet and surgically transferred to the left and right uterus, respectively, of each H (n = 6) or C (n = 7) synchronized recipients. Close to term, four combinations, CC (n = 10), CH (n = 13), HC (n = 13), and HH (n = 6), of feto-placental units were collected, for biometry analyses. Fatty acid (FA) profiles were determined in placental labyrinth, decidua, fetal plasma, and fetal liver by gas chromatography and explored further by principal component analysis (PCA). Candidate gene expression was also analyzed by RT-qPCR in the placenta and fetal liver. Data were analyzed by Kruskal–Wallis followed by Dunn’s pairwise comparison test. Combinations of different data sets were combined and explored by multifactorial analysis (MFA).

Results: Compared to controls, HH fetuses were hypotrophic with reduced placental efficiency and altered organogenesis, CH presented heavier placenta but less efficient, whereas HC presented a normal biometry. However, the MFA resulted in a good separation of the four groups, discriminating the effects of each period of exposure. HD during gestation led to reduced gene expression (nutrient transport and metabolism) and big changes in FA profiles in both tissues with increased membrane linoleic acid, lipid storage, and polyunsaturated-to-saturated FA ratios. Pre-conception exposure had a major effect on fetal biometry and organogenesis in HH, with specific changes in FA profiles (increased MUFAs and decreased LCPUFAs).

Conclusion: Embryo origin left traces in end-gestation feto-placental unit; however, maternal diet during gestation played a major role, either negative (HD) or positive (control). Thus, an H embryo developed favorably when transferred to a C recipient (HC) with normal biometry at term, despite disturbed and altered FA profiles.

Screening and stability analysis of reference genes in fasting caecotrophy model in rabbits

BACKGROUND

The selection and validation of stably expressed reference genes is key for accurately quantifying the mRNA abundance of genes under different treatments. In the rabbit model of fasting caecotrophy, reports about the selection of stable reference genes are not available.

METHODS AND RESULTS

This study aims to screen suitable reference genes in different tissues (including uterus, cecum, and liver) of rabbits between control and fasting caecotrophy groups. RT-qPCR was used to analyze the expression levels of eight commonly used reference genes (including GAPDH, 18S rRNA, B2M, CYP, HPRT1, β-actin, H2afz, Ywhaz), and RefFinder (including geNorm, NormFinder, and BestKeeper) was used to analyze the expression stability of these reference genes. Our results showed that the most stable reference genes were different in different tissues and treatments. In the control and fasting caecotrophy groups, CYP, GAPDH and HPRT1 were proven to be the top stable reference genes in the uterus, cecum, and liver tissues, respectively. GAPDH and Ywhaz were proven to be the top two stable reference genes among uterus, cecum, and liver in both control and fasting caecotrophy groups.

CONCLUSIONS

Our results indicated that the combined analysis of three or more reference genes (GAPDH, HPRT1, and Ywhaz) are recommended to be used for RT-qPCR normalization in the rabbit model of fasting caecotrophy, and that GAPDH is a better choice than the other reference genes for normalizing the relative expression of target genes in different tissues of fasting caecotrophy rabbits.

Analysis of the stability of 70 housekeeping genes during iPS reprogramming

Studies on induced pluripotent stem (iPS) cells highly rely on the investigation of their gene expression which requires normalization by housekeeping genes. Whether the housekeeping genes are stable during the iPS reprogramming, a transition of cell state known to be associated with profound changes, has been overlooked. In this study we analyzed the expression patterns of the most comprehensive list to date of housekeeping genes during iPS reprogramming of a mouse neural stem cell line N31. Our results show that housekeeping genes' expression fluctuates significantly during the iPS reprogramming. Clustering analysis shows that ribosomal genes' expression is rising, while the expression of cell-specific genes, such as vimentin (Vim) or elastin (Eln), is decreasing. To ensure the robustness of the obtained data, we performed a correlative analysis of the genes. Overall, all 70 genes analyzed changed the expression more than two-fold during the reprogramming. The scale of this analysis, that takes into account 70 previously known and newly suggested genes, allowed us to choose the most stable of all genes. We highlight the fact of fluctuation of housekeeping genes during iPS reprogramming, and propose that, to ensure robustness of qPCR experiments in iPS cells, housekeeping genes should be used together in combination, and with a prior testing in a specific line used in each study. We suggest that the longest splice variants of Rpl13a, Rplp1 and Rps18 can be used as a starting point for such initial testing as the most stable candidates.

Reference genes selection for real-time quantitative PCR analysis in mouse germinal vesicle oocytes

Reference gene selection in mouse oocytes is an important task required to perform further adequate analysis of target gene expression levels. In the current work we have analyzed expression stability of the seven most commonly used reference genes (Actb, Eef1e1, Gapdh, H2afz, Ppia, Rpl4 and Ubc) in mouse oocytes at the germinal vesicle (GV) stage. We have performed analysis of expression stability of the above-mentioned reference genes with the three most commonly used software tools: geNorm, BestKeeper and NormFinder. Taking into account the results obtained from all of these programmes Gapdh, Rpl4 and H2afz seem to be suitable candidate reference genes in GV oocytes of mouse.

Cryosurvival of rabbit embryos obtained after superovulation with corifollitropin alfa with or without LH

The efficiency of an embryo bank depends on provision of optimal conditions for recovery, cryopreservation and transfer to a breed or strain. In this sense, increasing the number of embryos available using superovulation should improve the cryobank efficiency. However, vagueness of response to conventional protocols to control or increase ovarian response and the quality of oocytes and embryos and their cryotolerance remain a challenge. The aim of our study was to evaluate the effect of corifollitropin alpha (CTP) and a recombinant human FSH (rhFSH), alone or supplemented with rhLH, on embryo cryosurvival by in vitro development and OCT4 and NANOG mRNA abundance at blastocyst stage and offspring rate. In vitro development of vitrified embryos was not significantly affected by superstimulation with or without rhLH supplementation, resulting in similar development rates to those of the control groups (fresh and vitrified embryos from non-superstimulated donor does). Blastocysts developed from vitrified embryos showed higher levels of OCT4 transcript abundance than fresh control, while NANOG transcript abundance was only higher in the blastocysts developed from vitrified embryos after superstimulation treatment in comparison with control groups. The implantation and offspring rates at birth were negatively affected by supplementation with rhLH. Both rhFSH or CTP vitrified embryo groups showed an implantation rate similar to those of the control groups, but an offspring rate lower than control. In conclusion, embryos produced using corifollitropin alpha did not compromise the cryosurvival of vitrified embryos in the rabbit. In addition, this study points out the negative effect of rhLH supplementation in terms of offspring rate on embryo vitrification.

Pre-implantation Development of Domestic Animals

During the first days following fertilization, cells of mammalian embryo gradually lose totipotency, acquiring distinct identity. The first three lineages specified in the mammalian embryo are pluripotent epiblast, which later gives rise to the embryo proper, and two extraembryonic lineages, hypoblast (also known as primitive endoderm) and trophectoderm, which form tissues supporting development of the fetus in utero. Most of our knowledge regarding the mechanisms of early lineage specification in mammals comes from studies in the mouse. However, the growing body of evidence points to both similarities and species-specific differences. Understanding molecular and cellular mechanisms of early embryonic development in nonrodent mammals expands our understanding of basic mechanisms of differentiation and is essential for the development of effective protocols for assisted reproduction in agriculture, veterinary medicine, and for biomedical research. This review summarizes the current state of knowledge on key events in epiblast, hypoblast, and trophoblast differentiation in domestic mammals.

Comparison of culture media indicates a role for autologous serum in enhancing phenotypic preservation of rabbit limbal stem cells in explant culture

In this study, we aimed to compare the effects of six different cell culture media and autologous serum (AS) on the phenotypic characteristics of rabbit limbal epithelial stem cells (LESC) cultivated on porous polyethylene terephthalate (PET) membranes. Limbal explants from rabbit corneas were grown on PET membrane inserts in five different media: DMEM-F12 with fetal bovine serum (FBS) (DMEM-F12-FBS), with pluripotin (DMEM-F12-pluripotin) and with autologous serum (DMEM-F12-AS), Epilife, Keratinocyte Serum Free Medium (KSFM) and Defined-Keratinocyte Serum Free Medium. The effects of different media were evaluated by total cell yield from explants, measuring the expression of proteins by immunofluorescence and gene expression by Real Time PCR. In all five media tested, most of the limbal epithelial cells (LEC) which proliferated from explants were positive for cytokeratin (CK) 14 (85-90%), indicating that all five media support the growth of LESC from explants. The expression of differentiation markers; CK 3 and 12 was highest in DMEM-F12-FBS (56%), was lower in Epilife and KSFM (26 and 19%, respectively), with the lowest values (13%) obtained in DMEM-F12-AS. Gene expression of limbal cultures on PET membrane inserts was compared to fresh limbal tissue. In DMEM-F12-FBS, DMEM-F12-pluripotin, and DMEM-F12-AS, expression of potential LESC markers CXCR4 and polycomb complex protein BMI-1 were similar to limbal tissue. DMEM-F12 with 10% AS maintained a higher percentage of potential stem cell marker genes and lower expression of genes involved in differentiation compared to Epilife or KSFM. Our study shows that rabbit LEC can be cultivated on PET inserts using DMEM-F12 with autologous serum without a requirement for amniotic membrane or feeder cells.

Methods of RNA preparation affect mRNA abundance quantification of reference genes in pig maturing oocytes

To ensure accurate normalization and quantification of target RNA transcripts using reverse transcription quantitative polymerase chain reaction (RT-qPCR), most studies focus on the identification of stably expressed gene(s) as internal reference. However, RNA preparation methods could also be an important factor, especially for test samples of limited quantity (e.g. oocytes). In this study, we aimed to select appropriate reference gene(s), and evaluate the effect of RNA preparation methods on gene expression quantification in porcine oocytes and cumulus cells during in vitro maturation. Expression profiles of seven genes (GAPDH, 18S, YWHAG, BACT, RPL4, HPRT1 and PPIA) were examined, on RNA samples extracted from cumulus cells (RNeasy Kit) and oocytes (RNeasy Kit and Lysis Kit) during in vitro maturation, respectively. Interestingly, different RNA preparation methods were found to potentially affect the quantification of reference gene expression in pig oocytes cultured in vitro. After geNorm analyses, the most suitable genes for normalization were identified, GAPDH/18S for cumulus cells and YWHAG/BACT for oocytes, respectively. Thus, our results provide useful data and information on the selection of better reference genes and RNA preparation method for related functional studies.

Assessment of 'one-step' versus 'sequential' embryo culture conditions through embryonic genome methylation and hydroxymethylation changes

STUDY QUESTION

In comparison to in vivo development, how do different conditions of in vitro culture (‘one step’ versus ‘sequential medium’) impact DNA methylation and hydroxymethylation in preimplantation embryos?

SUMMARY ANSWER

Using rabbit as a model, we show that DNA methylation and hydroxymethylation are both affected by in vitro culture of preimplantation embryos and the effect observed depends on the culture medium used.

WHAT IS KNOWN ALREADY

Correct regulation of DNA methylation is essential for embryonic development and DNA hydroxymethylation appears more and more to be a key player. Modifications of the environment of early embryos are known to have long term effects on adult phenotypes and health; these probably rely on epigenetic alterations.

STUDY DESIGN SIZE, DURATION

The study design we used is both cross sectional (control versus treatment) and longitudinal (time-course). Each individual in vivo experiment used embryos flushed from the donor at the 2-, 4-, 8-, 16- or morula stage. Each stage was analyzed in at least two independent experiments. Each individual in vitro experiment used embryos flushed from donors at the 1-cell stage (19 h post-coïtum) which were then cultured in parallel in the two tested media until the 2-, 4-, 8- 16-cell or morula stages. Each stage was analyzed in at least three independent experiments. In both the in vivo and in vitro experiments, 4-cell stage embryos were always included as an internal reference.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Immunofluorescence with antibodies specific for 5-methylcytosine (5meC) and 5-hydroxymethylcytosine (5hmeC) was used to quantify DNA methylation and hydroxymethylation levels in preimplantation embryos. We assessed the expression of DNA methyltransferases (DNMT), of ten eleven translocation (TET) dioxigenases and of two endogenous retroviral sequences (ERV) using RT-qPCR, since the expression of endogenous retroviral sequences is known to be regulated by DNA methylation. Three repeats were first done for all stages; then three additional repetitions were performed for those stages showing differences or tendencies toward differences between the different conditions in the first round of quantification.

MAIN RESULTS AND THE ROLE OF CHANCE

The kinetics of DNA methylation and hydroxymethylation were modified in in vitro cultured embryos, and the observed differences depended on the type of medium used. These differences were statistically significant. In addition, the expression of TET1 and TET2 was significantly reduced in post-embryonic genome activation (EGA) embryos after in vitro culture in both tested conditions. Finally, the expression of two retroviral sequences was analyzed and found to be significantly affected by in vitro culture.

LIMITATIONS REASONS FOR CAUTION

Our study remains mostly descriptive as no direct link can be established between the epigenetic changes observed and the expression changes in both effectors and targets of the studied epigenetic modifications. The results we obtained suggest that gene expression could be affected on a large scale, but this remains to be confirmed.

WIDER IMPLICATIONS OF THE FINDINGS

Our results are in agreement with the literature, showing that DNA methylation is sensitive to in vitro culture. As we observed an effect of both tested culture conditions on the tested epigenetic marks and on gene expression, we cannot conclude which medium is potentially closest to in vivo conditions. However, as the observed effects are different, additional studies may provide more information and potential recommendations for the use of culture media in assisted reproductive technology.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by an ‘AMP diagnostic prénatal et diagnostic génétique’ 2012 grant from the French Agence de la Biomédecine. This study was performed within the framework of ANR LABEX ‘REVIVE’ (ANR-10-LABX-73). Authors are members of RGB-Net (TD 1101) and Epiconcept (FA 1201) COST actions. The authors declare that there is no competing interest.

Effect of luteinizing hormone on rabbit ovarian superstimulation and embryo developmental potential

Assisted reproduction technologies require ovarian stimulation to increase the number of oocytes and embryos. Currently, superstimulation is achieved by gonadotropin treatment, but the embryo yield and quality are highly variable. Commonly, commercial preparations derived from pituitary and urinary origin are used to superovulate. Hence, ovarian superstimulation protocols have usually included both FSH and LH. The appearance of recombinant gonadotropins manufactured by genetic engineering techniques has ensured high quality and batch-to-batch consistency. Moreover, this enables us to assess the importance of LH in the ovarian stimulation. The main aim of this study was to evaluate the effect of recombinant human LH supplementation (10%) on embryonic development produced by rabbit does superovulated with low or high concentration (18.75 or 37.50 IU) of recombinant human FSH (rhFSH). Females treated with rhFSH increased the ovulation rate, and it was significantly higher when the high FSH dose was supplemented with LH. The superstimulation treatment used did not significantly affect in vitro development rate until the expanded blastocyst stage. The results of this study seem to suggest that, in terms of superovulatory response, when rabbit does are treated with 37.5-IU rhFSH, the use of LH supplementation allows an increase in the number of follicles recruited and the quality of embryos, in terms of ability to develop in vitro until blastocyst, and the expression profile of OCT4, NANOG, and SOX2 genes is not affected.

Cloning and characterization of rabbit POU5F1, SOX2, KLF4, C-MYC and NANOG pluripotency-associated genes

While the rabbit (Oryctolagus cuniculus) is an important research model for aspects of human development and disease that cannot be studied in rodents, the lack of data on the genetic regulation of rabbit preimplantation development is a limitation. To assist in the understanding of this process, our aim was to isolate and characterize genes necessary for the induction and maintenance of cellular pluripotency. We are the first to report the isolation of complete coding regions of rabbit SOX2, KLF4, C-MYC and NANOG, which encode transcription factors that play crucial regulatory roles during early mammalian embryonic development. We determined the exon-intron boundaries and chromosomal localization of these genes using computational analysis. The sequences of mRNA and translated protein of the newly identified genes and those of POU5F1 were aligned to their mammalian orthologs to determine the degree of evolutionary conservation. Furthermore, the expression of these genes in embryonic and adult cells was studied at the mRNA and protein levels. We found the sequences and the expression pattern of these pluripotency-associated genes to be highly conserved between human and rabbit, indicating that the rabbit would be a valuable model for human preimplantation development. Implementing the newly identified genes either as biomarkers or as reprogramming factors might also pave the way towards the creation of stable pluripotent rabbit cell lines.

Vitrification alters rabbit foetal placenta at transcriptomic and proteomic level

Although numerous studies have demonstrated that cryopreservation alters gene expression, less is known about those embryos that implanted successfully and continued in gestation. To raise the question of the neutrality of this technique, we examine the effects of vitrification through gestation in rabbit before and after the implantation. We monitored the distribution of losses of 569 vitrified morulae, observing that embryos which reach the last pre-implantatory stage are able to implant. However, we found that not all implanted embryos had the ability to continue with their gestation. The results reveal that vitrification decreased foetus and maternal placenta weights at mid-gestation, but led to a higher offspring birth weight. A novel finding is that while no differences in gene expression were detected in pre-implantatory embryos at day 6, vitrification affects a gene and protein expression in the placenta at day 14. Our results for first time reveal strong evidence of modifications in implanted embryos subjected to vitrification, suggesting that the crucial step that vitrified embryos must overcome is the placenta formation. On the basis of these findings, our work leaves the question open as to whether the effects we observed that cause vitrification during foetal development could give rise to some type of physiological or metabolic alteration in adulthood.

Candidate gene expression patterns in rabbit preimplantation embryos developed in vivo and in vitro

PURPOSE

The levels and timing of expression of genes like BCLXL, HDAC1 and pluripotency marker genes namely, OCT4, SOX2, NANOG and KLF4 are known to influence preimplantation embryo development. Despite this information, precise understanding of their influence during preimplantation embryo development is lacking. The present study attempts to compare the expression of these genes in the in vivo and in vitro developed preimplantation embryos.

METHODS

The in vivo and in vitro developed rabbit embryos collected at distinct developmental stages namely, pronuclear, 2 cell, 4 cell, 8 cell, 16 cell, Morula and blastocyst were compared at the transcriptional and translational levels using Real Time PCR and immunocytochemical studies respectively.

RESULTS

The study establishes the altered levels of candidate genes at the transcriptional level and translational level with reference to the zygotic genome activation (ZGA) phase of embryo development in the in vivo and in vitro developed embryos. The expression of OCT4, KLF4, NANOG and SOX2 genes were higher in the in vitro developed embryos whereas and HDAC1 was lower. BCLXL expression had its peak at ZGA in in vivo developed embryos. Protein expression of all the candidate genes was observed in the embryos. BCLXL, KLF4 and NANOG exhibited diffused localisation whereas HDAC1, OCT4, and SOX2 exhibited nuclear localisation.

CONCLUSIONS

This study leads to conclude that BCLXL peak expression at the ZGA phase may be a requirement for embryo development. Further expression of all the candidate genes was influenced by ZGA phase of development at the transcript level, but not at the protein level.

Optimisation of reference genes for gene-expression analysis in a rabbit model of left ventricular diastolic dysfunction

Left ventricular diastolic dysfunction (LVDD) is characterized by the disturbance of ventricle’s performance due to its abnormal relaxation or to its increased stiffness during the diastolic phase. The molecular mechanisms underlying LVDD remain unknown. We aimed to identify normalization genes for accurate gene-expression analysis of LVDD using quantitative real-time PCR (RT-PCR) in a new rabbit model of LVDD. Eighteen rabbits were fed with a normal diet (n = 7) or a 0.5% cholesterol-enriched diet supplemented with vitamin D2 (n = 11) for an average of 14.5 weeks. We validated the presence of LVDD in this model using echocardiography for diastolic function assessment. RT-PCR was performed using cDNA derived from left ventricle samples to measure the stability of 10 genes as candidate reference genes (Gapdh, Hprt1, Ppia, Sdha, Rpl5, Actb, Eef1e1, Ywhaz, Pgk1, and G6pd). Using geNorm analysis, we report that Sdha, Gapdh and Hprt1 genes had the highest stability (M <0.2). By contrast, Hprt1 and Rpl5 genes were found to represent the best combination for normalization when using the Normfinder algorithm (stability value of 0.042). Comparison of both normalization strategies highlighted an increase of natriuretic peptides (Bnp and Anp), monocytes chemotactic protein-1 (Mcp-1) and NADPH oxidase subunit (Nox-2) mRNA expressions in ventricle samples of the hypercholesterolemic rabbits compared to controls (P<0.05). This increase correlates with LVDD echocardiographic parameters and most importantly it molecularly validates the presence of the disease in our model. This is the first study emphasizing the selection of stable reference genes for RT-PCR normalization in a rabbit model of LVDD.

Transcript analysis of heat shock protein 72 in vitrified 2-cell mouse embryos and subsequent in vitro development

OBJECTIVE

The aim of the study was to compare the effects of two different concentrations of cryoprotectants by cryotopvitrification on survival, developmental capacity and Heat shock protein 72 (Hsp72) expression of two-cell mouse embryos.

MATERIALS AND METHODS

In this experimental study, transcript analysis of Hsp72 gene was performed on non-vitrified and vitrified 2-cell mouse embryos via a nested quantitative polymerase chain reaction (nqPCR) subsequent to normalization with Hprt1 as the reference gene. The different cryoprotectant combinations were 15% (vit1:7.5% of each ethylene glycol (EG) and dimethyl sulfoxide (DMSO), 30% (vit2:15% EG + 15% DMSO) and control group with no cryoprotectants. Vitrified and fresh 2-cell embryos were cultured to obtain cleavage and blastocyst formation rates. The results were analyzed via one-way analysis of variance and the mean values were compared with least significant difference (LSD) (p< 0.05).

RESULTS

The relative expression of Hsp72 in vit2 (30% v/v) was significantly higher than vit1 (15% v/v). Survival rates were the same for both vitrification treatments and significantly lower than the control group. Cleavage and blastocyst rates in vit1 were significantly higher than vit2 while those in two vitrified groups were significantly lower than the control group.

CONCLUSION

Our developmental data demonstrated that vit1 treatment (7.5% EG and 7.5% DMSO) was more efficient than vit2 (15% EG and 15% DMSO) in mouse embryos. The cryotopvitrification with two concentrations of cryoprotectants caused the relative changes of Hsp72 transcript level, but the stability of the gene in vit1 was significantly higher than vit2 and closer to the fresh 2-cell embryos.

Sexual dimorphism of the feto-placental phenotype in response to a high fat and control maternal diets in a rabbit model

Maternal environment during early developmental stages plays a seminal role in the establishment of adult phenotype. Using a rabbit model, we previously showed that feeding dams with a diet supplemented with 8% fat and 0.2% cholesterol (HH diet) from the prepubertal period and throughout gestation induced metabolic syndrome in adult offspring. Here, we examined the effects of the HH diet on feto-placental phenotype at 28 days post-coïtum (term = 31 days) in relation to earlier effects in the blastocyst (Day 6). At 28 days, both male and female HH fetuses were intrauterine growth retarded and dyslipidemic, with males more affected than females. Lipid droplets accumulated in the HH placentas' trophoblast, consistent with the increased concentrations in cholesteryl esters (3.2-fold), triacylglycerol (2.5-fold) and stored FA (2.12-fold). Stored FA concentrations were significantly higher in female compared to male HH placentas (2.18-fold, p<0.01), whereas triacylglycerol was increased only in HH males. Trophoblastic lipid droplet accumulation was also observed at the blastocyst stage. The expression of numerous genes involved in lipid pathways differed significantly according to diet both in term placenta and at the blastocyst stage. Among them, the expression of LXR-α in HH placentas was reduced in HH males but not females. These data demonstrate that maternal HH diet affects the blastocyst and induces sex-dependent metabolic adaptations in the placenta, which appears to protect female fetuses from developing severe dyslipidemia.

Oviductal and endometrial mRNA expression of implantation candidate biomarkers during early pregnancy in rabbit

Prenatal losses are a complex problem. Pregnancy requires orchestrated communication between the embryo and the uterus that includes secretions from the embryo to signal pregnancy recognition and secretion and remodelling from the uterine epithelium. Most of these losses are characterized by asynchronization between embryo and uterus. To better understand possible causes, an analysis was conducted of gene expression of a set of transcripts related to maternal recognition and establishment of rabbit pregnancy (uteroglobin,SCGB1A1; integrin α1,ITGA1; interferon-γ,IFNG; vascular endothelial growth factor,VEGF) in oviduct and uterine tissue at 16, 72 or 144 h post-ovulation and insemination. In the oviduct tissue, a significant decrease in the level ofSCGB1A1mRNA expression was observed from 144 h post-ovulation. In the case ofITGA1, the transcript abundance was initially lower, but mRNA expression increased significantly at 72 and 144 h post-ovulation. ForIFNG, a huge decrease was observed from 16 to 72 h post-ovulation. Finally, no significant differences were observed in theVEGFtranscript. For the endometrium, the results showed a significant decline in the level ofSCGB1A1mRNA expression from 16 to 144 h post-ovulation induction. The highest levels ofITGA1transcript were detected at 144 h, followed by the 16 h group and lower at 72 h post-ovulation. ForIFNGthere were no significant differences among post-ovulation induction times. Finally, it was possible to observe thatVEGFmRNA abundance was present at low levels at 16 h post-ovulation and remained low at 72 h, but increased at 144 h. The functional significance of these observations may provide new insights into the maternal role in prenatal losses.

Transcriptome profiling of rabbit parthenogenetic blastocysts developed under in vivo conditions

Parthenogenetic embryos are one attractive alternative as a source of embryonic stem cells, although many aspects related to the biology of parthenogenetic embryos and parthenogenetically derived cell lines still need to be elucidated. The present work was conducted to investigate the gene expression profile of rabbit parthenote embryos cultured under in vivo conditions using microarray analysis. Transcriptomic profiles indicate 2541 differentially expressed genes between parthenotes and normal in vivo fertilised blastocysts, of which 76 genes were upregulated and 16 genes downregulated in in vivo cultured parthenote blastocyst, using 3 fold-changes as a cut-off. While differentially upregulated expressed genes are related to transport and protein metabolic process, downregulated expressed genes are related to DNA and RNA binding. Using microarray data, 6 imprinted genes were identified as conserved among rabbits, humans and mice: GRB10, ATP10A, ZNF215, NDN, IMPACT and SFMBT2. We also found that 26 putative genes have at least one member of that gene family imprinted in other species. These data strengthen the view that a large fraction of genes is differentially expressed between parthenogenetic and normal embryos cultured under the same conditions and offer a new approach to the identification of imprinted genes in rabbit.

Investigation of gene expression profiles before and after embryonic genome activation and assessment of functional pathways at the human metaphase II oocyte and blastocyst stage

OBJECTIVE

To compare the oocyte versus the blastocyst transcriptome and provide data on molecular pathways before and after embryonic genome activation.

DESIGN

Prospective laboratory research study.

SETTING

An IVF clinic and a specialist preimplantation genetics laboratory.

PATIENT(S)

Couples undergoing or having completed IVF treatment donating surplus oocytes or cryopreserved blastocysts after patient consent.

INTERVENTION(S)

Sets of pooled metaphase II (MII) oocytes or blastocysts were processed for RNA extraction, RNA amplification, and analysis with the use of the Human Genome Survey Microarrays v2.0 (Applied Biosystems).

MAIN OUTCOME MEASURE(S)

Association of cell type and gene expression profile.

RESULT(S)

Totals of 1,909 and 3,122 genes were uniquely expressed in human MII oocytes and human blastocysts respectively, and 4,910 genes were differentially expressed between the two sample types. Expression levels of 560 housekeeping genes, genes involved in the microRNA processing pathway, as well as hormones and hormone receptors were also investigated.

CONCLUSION(S)

The lists of genes identified may be of use for understanding the processes involved in early embryo development and blastocyst implantation, and for identifying any dysregulation leading to infertility.

Selection of suitable reference genes for normalization of quantitative real-time PCR in cartilage tissue injury and repair in rabbits

When studying the altered expression of genes associated with cartilage regeneration by quantitative real-time RT-PCR (RT-qPCR), reference genes with highly stable expression during different stages of chondrocyte developmental are necessary to normalize gene expression accurately. Until now, no reports evaluating expression changes of commonly used reference genes in rabbit articular cartilage have been published. In this study, defects were made in rabbit articular cartilage, with or without insulin-like growth factor 1 (IGF-1) treatment, to create different chondrocyte living environments. The stability and intensity of the expressions of the candidate reference genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH), 18S Ribosomal RNA (18S rRNA), cyclophilin (CYP), hypoxanthine phosphoribosyl transferase (HPRT1), and beta-2-microglobulin (B2M) were evaluated. The data were analyzed by geNorm and NormFinder. B2M and 18S rRNA were identified to be suitable reference genes for rabbit cartilage tissues.

Up-regulation of insulin-like growth factor I and uteroglobin in in vivo-developed parthenogenetic embryos

Parthenote embryos are being considered as an alternative source of embryonic stem cells. However, as there is still a dearth of knowledge of this kind of embryos, a better understanding of their biology is needed for their application. In this work, we studied the differences and similarities between parthenotes and normal embryos at the blastocyst stage in vivo developed. We analysed the expression of factor OCT-4, vascular endothelial growth factor (VEGF), insulin-like growth factor I (IGF-I) and uteroglobin (UG) by real-time PCR. To do so, oocytes were recovered and after activation procedure were transferred by ventral middle laparoscopy to receptive does to undergo completely in vivo development. Does were slaughtered 6 days post-ovulation induction, and parthenote and normal embryos were recovered for mRNA expression analysis. Our results reported that parthenotes and normal embryos showed similar mRNA expression for OCT-4 and VEGF. However, IGF-I and UG showed to be over-expressed in parthenote embryos. Thus, our study highlights that despite the in vivo development of parthenotes, they still seem to have an altered expression and, therefore, to be different to normal embryos. The altered expression pattern of parthenote embryos suggests that these embryos should be studied carefully before future application.

Spatial and temporal distribution of Oct-4 and acetylated H4K5 in rabbit embryos

Rabbit is a unique species to study human embryology; however, there are limited reports on the key transcription factors and epigenetic events of rabbit embryos. This study examined the Oct-4 and acetylated H4K5 (H4K5ac) patterns in rabbit embryos using immunochemistry staining. The average intensity of the Oct-4 signal in the nuclei of the whole embryo spiked upon fertilization, then decreased until the 8-cell stage and increased afterwards until the compact morula (CM) stage. It decreased thereafter from the CM stage to the early blastocyst (EB) stage, with a minimum at the expanded blastocyst (EXPB) stage and came back to a level similar to that of the CM-stage embryos in the hatching blastocysts (HB). The Oct-4 signal was observed in both the inner cell mass (ICM) and the trophectoderm (TE) cells of blastocysts. The average H4K5ac signal intensity of the whole embryo increased upon fertilization, started to decrease at the 4-cell stage, reached a minimum at the 8-cell stage, increased again at the EXPB stage and peaked at the HB stage. While TE cells maintained similar levels of H4K5ac throughout the blastocyst stages, ICM cells of HB showed higher levels of H4K5ac than those of EB and EXPB. Understanding key genetic and epigenetic events during early embryo development will help to identify factors contributing to embryo losses and consequently improve embryo survival rates. As a preferred laboratory species for many human disease studies such as atherosclerosis, rabbit is also a pioneer species in the development of several embryo biotechnologies, such as IVF, transgenesis, animal cloning, embryo cryopreservation and embryonic stem cells. However, there are limited reports on key transcription factors and epigenetic events of rabbit embryos. In the present study, we documented the temporal and spatial distribution of Oct-4 protein and H4K5 acetylation during early embryo development using the immunostaining approach. We also compared the patterns of these two important biomarkers between the inner cell mass (ICM) and the trophectoderm (TE) cells in blastocyst-stage embryos. Our findings suggest that a combination of Oct-4, H4K5ac and possibly other biomarkers such as Cdx-2 is needed to accurately identify different lineages of cells in morula and blastocyst stage rabbit embryos. Importantly, we revealed a novel wave of Oct-4 intensity change in the ICM cells of rabbit blastocysts. The signal was high at the early blastocyst stage, reached a minimum at the expanded blastocyst stage and returned to a high level at the hatching blastocyst stage. We hypothesize that the signal may have reflected the regulation of Oct-4 through enhancer switching and therefore may be related to cell lineage formation in rabbit embryos. These findings enrich our understanding on key genetic and epigenetic programming events during early embryo development in rabbits.

New challenges in the analysis of gene transcription in bovine blastocysts

In the last years, enormous progress has been made in the analysis of gene transcription at the blastocyst stage. The study of gene expression at this early stage of development is challenging because of the very small amount of starting material, which limits the use of traditional mRNA analysis approaches such as Northern blot. Another problem is the difficulty for data normalization, particularly the identification of the best housekeeping gene with the lowest fluctuation under different developmental conditions. Moreover, the transcriptional analysis of embryo biopsies or individual embryos needs to take into consideration that the blastocyst is a transitional stage of development, which is composed of three different types of cells (trophoblast, epiblast and primitive ectoderm) with different patterns of gene expression, and that there are large differences between male and female blastocysts. In this review, we analyse the different specific and sensitive tools available to compare mRNA expression levels of specific genes at the blastocyst stage, and how the protocol and the analytical method used can influence the results dramatically. Finally, we describe future research challenges to identify candidate genes related to developmental competence of bovine blastocysts, not only in terms of pregnancy rates but also in relation to adverse long-term consequences in the adult animal.

Gene expression profiling of pluripotency and differentiation-related markers in cat oocytes and preimplantation embryos

During mammalian preimplantation development, two successive differentiation events lead to the establishment of three committed lineages with separate fates: the trophectoderm, the primitive endoderm and the pluripotent epiblast. In the mouse embryo, the molecular mechanisms underlying these two cell fate decisions have been studied extensively, leading to the identification of lineage-specific transcription factors. Species-specific differences in expression patterns of key regulatory genes have been reported, raising questions regarding their role in different species. The aim of the present study was to characterise the gene expression patterns of pluripotency (OCT4, SOX2, NANOG) and differentiation (CDX2, GATA6)-related markers during feline early development using reverse transcription–quantitative polymerase chain reaction. In addition, we assessed the impact of in vitro development on gene expression by comparing transcript levels of the genes investigated between in vitro and in vivo blastocysts. To normalise quantitative data within different preimplantation embryo stages, we first validated a set of stable reference genes. Transcript levels of all genes investigated were present and changed over the course of preimplantation development; a highly significant embryo-stage effect on gene expression was observed. Transcript levels of OCT4 were significantly reduced in in vitro blastocysts compared with their in vivo counterparts. None of the other genes investigated showed altered expression under in vitro conditions. The different gene expression patterns of OCT4, SOX2, CDX2 and GATA6 in cat embryos resembled those described in mouse embryos, indicative of a preserved role for these genes during early segregation. However, because of the absence of any upregulation of NANOG transcription levels after embryonic genome activation, it is unlikely that NANOG is a key regular of lineage segregation. Such results support the hypothesis that the behaviour of early lineage markers can be species specific. The present study also revealed a pool of maternal NANOG mRNA transcripts, the role of which remains to be elucidated. Comparing transcription levels of these genes between in vivo and in vitro blastocysts revealed low levels of OCT4 mRNA in the latter, which may contribute to the reduced developmental competence of embryos under suboptimal conditions.

Comparative analysis of nuclear transfer embryo-derived mouse embryonic stem cells. Part II: gene regulation

In a mouse model nuclear transfer embryo-derived embryonic stem cell lines (ntESCs) of various genetic backgrounds and donor cell types were compared with reference ESCs and analyzed comprehensively at molecular level as a second part of a larger study. Expression profiles of ntESCs established by different NT-methods (piezoelectric microinjection or zona-free) were indistinguishable. However, expression profiling analyses identified differentially regulated genes between reference ESCs and ntESCs from different genetic backgrounds. A number of pluripotency and stemness marker genes significantly differed at the mRNA level between the cell lines. However, cluster and lineage analyses revealed that such differences had no effect on cell differentiation and cell fate. Regardless of the donor cell type, gene expression profiles of ntESCs were more similar to each other than to their counterpart fertilized embryo-derived ESCs of the same genotype. Overall, the results indicated that expression profile differences may be related to the genotype rather than to technical variations.

effect of embryonic genotype on reference gene selection for RT-qPCR normalization

To obtain more reliable results with biological significance, it requires data normalization using an appropriate internal control gene. In rabbits, there are classic stable reference genes that have been identified for normalization in oocytes and pre-implantation stage embryos. However, effects of embryonic genotype on reference gene selection have not been elucidated. The aim of this study was to test (i) the stability of mRNA transcription level for histone (H2afz) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) genes in rabbit blastocysts from two lines selected by different criteria (litter size and post-weaning daily weight gain) and (ii) its influence on biological significance examined by means of a set of embryonic transcripts, such as POU5F1 (Oct-4), epidermal growth factor receptor (erbB3), transforming growth factor-beta2, vascular endothelial growth factor and gamma interferon (Ifn-gamma). The geNorm, NormFinder and BestKeeper programs showed similar results, pointing out that H2afz and GAPDH were the most stable reference genes in rabbits selected on litter size at weaning. Moreover, our study revealed that embryonic genotype affected target gene expression when a single reference gene was used to analyse mRNA expression in blastocysts. Results showed that GAPDH gene is better than H2afz for gene expression studies of both embryo genotypes. A normalization factor derived from H2afz and GAPDH is likely to be appropriate when RT-qPCR was performed in rabbit embryos with different genotypes.

Parthenogenic blastocysts cultured under in vivo conditions exhibit proliferation and differentiation expression genes similar to those of normal embryos

Parthenote embryos offer multiple possibilities in biotechnological investigation, such as stem cell research. However, there is still a dearth of knowledge of this kind of embryo. In this study, development and ploidy were analysed in parthenotes under in vitro and in vivo culture conditions. Subsequently, using real-time PCR, the expressions of factor OCT-4, Vascular Endothelial Growth Factor, Epidermal Growth Factor Receptor 3 and Transforming Growth Factor β2 genes were analysed to compare the embryo types at the blastocyst stage. Development and implantation of parthenote embryos were described after transfer at day 10 of pregnancy. Parthenotes showed similar blastocyst development for both culture conditions and most of the parthenotes produced were diploid. However, parthenotes developed under in vivo conditions showed similar mRNA expression of OCT-4, VEGF and TGF-β2 to 5 and 6 day old blastocysts. In contrast, parthenotes developed under in vitro conditions had altered the expression pattern of these genes, except for erbB3 mRNA. Finally, transferred parthenotes had the ability to implant but showed severe growth retardation and lesser size. This is the first demonstration of the influence of culture conditions on parthenote mRNA expression. Our study highlights the importance of culture conditions in subsequent uses of parthenotes, such as the production of stem cell lines.

Effect of different culture systems on mRNA expression in developing rabbit embryos

The rate of zygotes in vitro developed to hatched blastocyst stage was evaluated between two different commercial media (TCM-199 and Ham's F10) and two different culture systems (renewal and non-renewal single medium) to determine the effects of culture conditions on rabbit embryo preimplantation development. The relative transcript abundances of OCT4, vascular endothelial growth factor (VEGF) and epidermal growth factor receptor 3 (erbB3) of resultant blastocysts were also analysed and compared with in vivo developed blastocysts. Results showed an important divergence in mRNA expression between embryos developed under in vivo and in vitro conditions despite there being no significant difference in hatching blastocyst rates between different culture systems and different media. For OCT4, transcript abundance of in vitro culture embryos differs from their in vivo chronological counterparts, but, when the medium is renewed, mRNA expression seemed similar to in vivo developed 4-day-old embryos. In addition, VEGF and erbB3 expression showed marked variation between different in vitro conditions. Therefore, the study of specific transcript abundance in rabbit blastocyst provides a more detailed description of which alterations in gene expression occur due in vitro conditions, and further studies should be carried out to reduce current limitations of long-term culture of rabbit pre-implantation embryos.

mRNA expression in rabbit blastocyst and endometrial tissue of candidate gene involved in gestational losses

Gestation is a complex process that involves different growth factors, cytokines and adhesion proteins related with embryo development, cellular differentiation and proliferation, embryo-endometrium interaction, angiogenesis, maternal-embryonic recognition and growth development of placenta and embryos. In this study, we examine pre-implantational (at 6 days of gestation) and gestational (at 12 days and total from ovulation to birth) losses in two rabbit lines selected by different criteria (post-weaning daily gain and litter size) and the pattern of a set of candidate transcripts, at 6 days of gestation, related with embryo development and implantation process, such as Oct-4, epidermal growth factor receptor 3 (erbB3), Transforming Growth Factor β2, Vascular Endothelial Growth Factor and Interferon γ and related with insulin-like growth factors signalling as insulin growth factors I and II and their receptors in rabbit blastocysts and endometrial tissue. Similar pre-implantational losses were obtained in both lines. However, the gestational losses of the line selected by post-weaning daily gain clearly mirrored an increase in losses by 50% at 12 days and at birth (22.4 vs 9.5 and 50.2 vs 25.4, respectively, between line selected by post-weaning daily gain and line selected by litter size). In blastocysts and endometrial tissue at 6 days of gestation qRT-PCR assays indicated that the mean insulin-like growth factor (IGF)-IIR mRNA expression was down-regulated in line selected by post-weaning daily gain. Dysregulation of the IGF-IIR could be potential reasons for induced gestational losses. We conclude that IGF-IIR gene expression in blastocyst and endometrial tissue at 6th day of gestation tends to decline in line selected by post-weaning daily gain. The functional significance related with gestational losses is uncertain.

Selection of reference genes for normalization of quantitative real-time PCR in organ culture of the rat and rabbit intervertebral disc

Background

The accuracy of quantitative real-time RT-PCR (qRT-PCR) is often influenced by experimental artifacts, resulting in erroneous expression profiles of target genes. The practice of employing normalization using a reference gene significantly improves reliability and its applicability to molecular biology. However, selection of an ideal reference gene(s) is of critical importance to discern meaningful results. The aim of this study was to evaluate the stability of seven potential reference genes (Actb, GAPDH, 18S rRNA, CycA, Hprt1, Ywhaz, and Pgk1) and identify most stable gene(s) for application in tissue culture research using the rat and rabbit intervertebral disc (IVD).

Findings

In vitro, four genes (Hprt1, CycA, GAPDH, and 18S rRNA) in rat IVD tissue and five genes (CycA, Hprt1, Actb, Pgk1, and Ywhaz) in rabbit IVD tissue were determined as most stable for up to 14 days in culture. Pair-wise variation analysis indicated that combination of Hprt1 and CycA in rat and the combination of Hprt1, CycA, and Actb in rabbit may most stable reference gene candidates for IVD tissue culture.

Conclusions

Our results indicate that Hprt1 and CycA are the most stable reference gene candidates for rat and rabbit IVD culture studies. In rabbit IVD, Actb could be an additional gene employed in conjunction with Hprt1 and CycA. Selection of optimal reference gene candidate(s) should be a pertinent exercise before employment of PCR outcome measures for biomedical research.

Establishing reference genes for use in real-time quantitative PCR analysis of early equine embryos

Real-time quantitative PCR (qPCR) is invaluable for investigating changes in gene expression during early development, since it can be performed on the limited quantities of mRNA contained in individual embryos. However, the reliability of this method depends on the use of validated stably expressed reference genes for accurate data normalisation. The aim of the present study was to identify and validate a set of reference genes suitable for studying gene expression during equine embryo development. The stable expression of four carefully selected reference genes and one developmentally regulated gene was examined by qPCR in equine in vivo embryos from morula to expanded blastocyst stage. SRP14, RPL4 and PGK1 were identified by geNorm analysis as stably expressed reference genes suitable for data normalisation. RPL13A expression was less stable and changed significantly during the period of development examined, rendering it unsuitable as a reference gene. As anticipated, CDX2 expression increased significantly during embryo development, supporting its possible role in trophectoderm specification in the horse. In summary, it was demonstrated that evidence-based selection of potential reference genes can reduce the number needed to validate stable expression in an experimental system; this is particularly useful when dealing with tissues that yield small amounts of mRNA. SRP14, RPL4 and PGK1 are stable reference genes suitable for normalising expression for genes of interest during in vivo morula to expanded blastocyst development of horse embryos.

The effects of vitrification on gene expression in mature mouse oocytes by nested quantitative PCR

PURPOSE

this study was conducted on the effects of vitrification cryotop method on gene expression of mature oocytes in Mus musculus.

METHODS

transcript analyses of three mouse genes, namely Mater, Hook1 and Sod1, were performed upon non-vitrified and vitrified oocytes with different concentrations of dimethyl sulfoxide (DMSO) and ethylene glycol (EG),15%: 7.5% DMSO + 7.5% EG, and 30%: 15% DMSO + 15% EG, using cryotop following normalization of transcripts with Hprt1 by nested quantitative PCR.

RESULTS

vitrification caused down-regulation of Mater and Hook1 and up-regulation of Sod1 when lower concentrations of cryoprotectants were used as opposed to the control group. The relative expression of Sod1 in vit(2) (30% v/v) was significantly higher than vit(1) (15% v/v)(.) Quantitative transcript analysis of Mater and Hook1 for the vit(2) condition failed to produce any data. Survival rates were the same for both vitrification treatments and significantly lower than control group.

CONCLUSIONS

although vit(1) treatment had lower survival rate compared to control group, it demonstrated better stability comparing to vit(2) based on the transcript analysis.

Selection of reference genes for quantitative real-time PCR in equine in vivo and fresh and frozen-thawed in vitro blastocysts

BACKGROUND

Application of reverse transcription quantitative real-time polymerase chain reaction is very well suited to reveal differences in gene expression between in vivo and in vitro produced embryos. Ultimately, this may lead to optimized equine assisted reproductive techniques. However, for a correct interpretation of the real-time PCR results, all data must be normalized, which is most reliably achieved by calculating the geometric mean of the most stable reference genes. In this study a set of reliable reference genes was identified for equine in vivo and fresh and frozen-thawed in vitro embryos.

FINDINGS

The expression stability of 8 candidate reference genes (ACTB, GAPDH, H2A/I, HPRT1, RPL32, SDHA, TUBA4A, UBC) was determined in 3 populations of equine blastocysts (fresh in vivo, fresh and frozen-thawed in vitro embryos). Application of geNorm indicated UBC, GAPDH, ACTB and HPRT1 as the most stable genes in the in vivo embryos and UBC, RPL32, GAPDH and ACTB in both in vitro populations. When in vivo and in vitro embryos were combined, UBC, ACTB, RPL32 and GAPDH were found to be the most stable. SDHA and H2A/I appeared to be highly regulated.

CONCLUSIONS

Based on these results, the geometric mean of UBC, ACTB, RPL32 and GAPDH is to be recommended for accurate normalization of quantitative real-time PCR data in equine in vivo and in vitro produced blastocysts.

Promoter analysis of the rabbit POU5F1 gene and its expression in preimplantation stage embryos

BACKGROUND

The POU5F1 gene encodes the octamer-binding transcription factor-4 (Oct4). It is crucial in the regulation of pluripotency during embryonic development and widely used as molecular marker of embryonic stem cells (ESCs). The objective of this study was to identify and to analyse the promoter region of rabbit POU5F1 gene; furthermore to examine its expression pattern in preimplantation stage rabbit embryos.

RESULTS

The upstream region of rabbit POU5F1 was subcloned sequenced and four highly conserved promoter regions (CR1-4) were identified. The highest degree of similarity on sequence level was found among the conserved domains between rabbit and human. Among the enhancers the proximal enhancer region (PE-1A) exhibited the highest degree of homology (96.4%). Furthermore, the CR4 regulator domain containing the distal enhancer (DE-2A) was responsible for stem cell-specific expression. Also, BAC library screen revealed the existence of a processed pseudogene of rabbit POU5F1. The results of quantitative real-time PCR experiments showed that POU5F1 mRNA was abundantly present in oocytes and zygotes, but it was gradually reduced until the activation of the embryonic genome, thereafter a continuous increase in POU5F1 mRNA level was observed until blastocyst stage. By using the XYClone laser system the inner cell mass (ICM) and trophoblast portions of embryos were microdissected and examined separately and POU5F1 mRNA was detected in both cell types.

CONCLUSION

In this study we provide a comparative sequence analysis of the regulatory region of rabbit POU5F1 gene. Our data suggest that the POU5F1 gene is strictly regulated during early mammalian development. We proposed that the well conserved CR4 region containing the DE-2A enhancer is responsible for the highly conserved ESC specific gene expression. Notably, we are the first to report that the rabbit POU5F1 is not restricted to ICM cells only, but it is expressed in trophoblast cells as well. This information may be well applicable to investigate further the possible phylogenetic role and the regulation of POU5F1 gene.