Skip to main content
  • Research article
  • Open access
  • Published:

Irregular transcriptome reprogramming probably causes thec developmental failure of embryos produced by interspecies somatic cell nuclear transfer between the Przewalski’s gazelle and the bovine



Interspecies somatic cell nuclear transfer (iSCNT) has been regarded as a potential alternative for rescuing highly endangered species and can be used as a model for studying nuclear–cytoplasmic interactions. However, iSCNT embryos often fail to produce viable offspring. The alterations in normal molecular mechanisms contributing to extremely poor development are for the most part unknown.


Przewalski’s gazelle–bovine iSCNT embryos (PBNT) were produced by transferring Przewalski’s gazelle fibroblast nuclei into enucleated bovine oocytes. The percentages of PBNT embryos that developed to morula/blastocyst stages were extremely low even with the use of various treatments that included different SCNT protocols and treatment of embryos with small molecules. Transcriptional microarray analyses of the cloned embryos showed that the upregulation of reprogramming-associated genes in bovine–bovine SCNT (BBNT) embryos was significantly higher than those observed in PBNT embryos (1527:643). In all, 139 transcripts related to various transcription regulation factors (TFs) were unsuccessfully activated in the iSCNT embryos. Maternal degradation profiles showed that 1515 genes were uniquely downregulated in the BBNT embryos, while 343 genes were downregulated in the PBNT embryos. Incompatibilities between mitochondrial DNA (mtDNA) and nuclear DNA revealed that the TOMM (translocase of outer mitochondrial membrane)/TIMM (translocase of inner mitochondrial membrane) complex-associated genes in BBNT embryos had the highest expression levels, while the PBNT embryos exhibited much lower expression rates.


Improper degradation of maternal transcripts, incomplete activation of TFs and abnormal expression of genes associated with mitochondrial function in PBNT embryos likely contributed to incomplete reprogramming of the donor cell nuclei and therefore led to the developmental failure of these cloned embryos.


Przewalski’s gazelle (Procapra przewalskii) is one of the more critically endangered Eurasian large mammals and is unique to China. Only around 350–400 mature individuals are thought to remain [1]. Their range exists today only in a small area surrounding Qinghai Lake [2, 3]. This gazelle’s fate is considered to be even more precarious than the giant panda [4]. Rescue and conservation programs are a challenge for wildlife biologists and ecologists, although management efforts are underway to provide for a more sustainable population [5]. Somatic cell nuclear transfer (SCNT) has been successfully utilized in the production of many mammal species including laboratory and domestic animals. One potential application of this technology is that it might be useful for the propagation of rare and endangered species. However, the major limitation in using this technology for species rescue is that oocytes and suitable recipients are rare, so intraspecies cloning of endangered species becomes an even more daunting task. Interspecies SCNT (iSCNT) where endangered animal somatic cell nuclei are transferred to domestic oocyte cytoplasts is an approach that might minimize the limitations of SCNT. Many trials of iSCNT have been reported in wildlife species such as the giant panda (Ailuropoda melanoleuca) [6, 7], Tibetan antelope (Pantholops hodgsonii) [8], Banteng (Bos javanicus) [9], yak (Bos grunniens) [10], Siberian tiger (Panthera tigris altaica) [11], and Sei whale (Balaenoptera borealis) [12]. The iSCNT cloned embryos generated in these studies had extremely poor development to the blastocyst stage. The best results of iSCNT in mammals occurred when using subspecies and sibling species that can hybridize naturally, such as among the cloned argali (Ovis ammon) [13] and the river buffalo (Bubalus bubalus arnee) [14].

The current hypotheses for the high developmental failure of iSCNT are that there is genomic incompatibility between the nucleus and the host ooplasm/cytoplasm [15] and between mitochondrial DNA (mtDNA) and nuclear DNA [16]. A major barrier that hinders the developing iSCNT embryo first occurs at the time of zygotic genome activation (ZGA) or embryonic genome activation (EGA) [17]. The zygote genome of the cloned embryo manifests itself independently from the maternal transcripts [18]. Compared with the development of normally fertilized preimplantation embryos, SCNT-derived embryos have to overcome many more challenges in silencing their somatic-specific genes while reactivating all of the embryo-related genes [19]. During the process, it also has to shed its differentiated phenotype and establish a totipotent state [20]. At the present time there are very few reports on embryonic gene function at the time of EGA in iSCNT embryos. Although it has been shown that small molecules such as valproic acid (VPA) [21], trichostatin A (TSA) [22], 5-aza-2′-deoxycytidine [23] can improve SCNT embryo development in vitro and in vivo, such beneficial effects have not been observed in interspecies cloned embryo development.

We transferred Przewalski’s gazelle fibroblast nuclei into enucleated bovine oocytes in this study. The Przewalski’s gazelle–bovine nuclear transfer (PBNT) embryos were treated with VPA and TSA and monitored for embryo development. However, the trials did not result in a significant improvements in iSCNT embryo development. To better understand why the embryos failed to develop and thrive, we used genomic computational methods to analyze the global reprogramming transcriptome between PBNT embryos and bovine–bovine nuclear transfer (BBNT) embryos at the time of the maternal–zygotic genomic transition (MZT). We first identified and then systematically analyzed the different gene expression patterns between PBNT and BBNT embryos in the global transcriptome, maternal mRNA degradation, transcription regulation-related genes and the aberrant expression of genes associated with mtDNA. Result of the quantitative PCR in the SCNT embryo and in the iSCNT embryos revealed excellent agreement with the microarray data. Putative mechanisms affecting developmental potential between SCNT and iSCNT embryos are discussed.


Treatment of PBNT embryos with VPA did not affect embryo development

Fused PBNT embryos were treated with VPA at 0.5, 1.0, 2.0 and 4.0 mM for 24 h. Two- and 8–16 cell embryo development rates were significantly higher in the 0.1 and 1.0 mM and control groups than in the 2.0 and 4.0 mM VPA groups (Table 1). There were no beneficial effects on morula or blastocyst development rates among VPA treatments and the controls. Fewer cloned embryos developed to blastocysts in the 0.5 mM (0.7%), 1.0 mM (0.7%) and control groups (1.5%). Better results were obtained in the 0.5 mM VPA group. When the embryos were treated with VPA at 0.5 mM for 24 h, the initial cleavage and 8–16-cell developmental rates were significantly better than in treatment times of 5, 12 and 48 h (Table 1). This improvement in development was not observed in the morula and blastocyst stage. There was no difference among the treatments and the controls. Longer exposure times to VPA did not improve PBNT embryo development.

Table 1 PBNT embryo development after VPA treatment

The PBNT embryos derived from green fluorescent protein (GFP)-expressing Oct-4-eGFP transgenic cells were further treated with 0.5 mM VPA for 24 h. Cleavage rates, 8–16-cell, morula and blastocyst development rates were no different between treatment groups and controls (Additional file 1: Table S1). Treatment of Oct-4-eGFP-derived PBNT embryos with VPA had no effect on embryo development.

Reverse nuclear transfer (RNT) did not improve PBNT embryo development

The reverse nuclear transfer (RNT) method significantly improved the blastocyst development rate in BBNT (36% in RNT vs 26% in SCNT). Whereas RNT did not improve the rate of PBNT embryo development (Table 2).

Table 2 Development of PBNT embryos derived from the reverse nuclear transfer protocol

Scatterplot comparison of different microarray datasets

Gene array analysis was performed using the Affymetrix gene chip bovine genome array (Santa Clara, CA, USA). A total of 1150 bovine oocytes (BOs), 309 8- to 16-cell BBNT embryos, 527 8- to 16-cell PBNT embryos, Przewalski’s gazelle fibroblasts (PCs) and bovine fibroblasts (BCs) were used in the computational analyses. The developmental stage and the morphology of the iSCNT embryos were with no obvious different from the control intra-species NT embryos (Additional file 2: Figure S1). High reproducibility was obtained between the replicates and datasets. The scatterplot compared the results of the log transformed gene expression levels and the differentially expressed gene profiles between two cell types (Figure 1). All of the treatments were repeated at least three times. High reproducibility was obtained between the replicates and datasets. The scatter plot compared the results of the log transformed gene expression levels and the differentially expressed gene profiles between two cell types (Figure 1). The Spearman correlation coefficients between the profiles of different cell types reflect the degree of change in transcriptomes (Additional file 3: Table S2). The results showed the BBNT oocytes and PBNT embryos had the most similar transcriptional profiles (R2 = 0.97) than other comparisons.

Figure 1
figure 1

Scatter plot compares the results of log transformed gene expression levels and the differentially expressed gene distribution pattern between the two cell types. Green indicates down- and red up-regulation of gene expression.

Reprogramming transcriptome analysis of BBNT and PBNT embryos

Hierarchical clustering of complete transcriptional profiling between the two different chip platforms indicated separations between any two sample groups (Figure 2A). More than 10,000 transcripts were detected in each sample, with similar profile patterns clustered between any two close cell or embryo types based on average Euclidean distance. The BCs and PCs showed highly unique and consistent expression patterns. The clustered coherent profiles revealed that BBNT embryos had more upregulated genes than did PBNT embryos at the same development stages (Additional file 4: Table S3). A Venn diagram illustrating shared and unique genes between BBNT and PBNT embryos, PCs and BCs, bovine oocytes and various regulatory factors is illustrated in Figure 2B.

Figure 2
figure 2

The hierarchical analyses and shared co-expression transcripts for different cell types. (A). Hierarchical clusters of overall gene expression profiles. Green indicates the down-and red the up-regulated gene expression. (B). Venn diagram of shared and unique genes among different transcriptomes.

Upon subtracting the genes expressed in the bovine oocytes and in the PCs and BCs from the transcriptome, there were, 1527 and 643 upregulated genes associated with nuclear reprogramming in the BBNT and PBNT embryos, respectively (Venn diagrams in Figure 3; Additional file 5: Figure S2 and Additional file 6: Table S4), which were 2.4 times higher in BBNT than in PBNT embryos. Three hundred nine reprogramming related genes were identified as being co-expressed in both the BBNT embryo and PBNT embryo samples. These co-expressed genes are associated with several important biological processes such as membrane-enclosed lumen, RNA processing/splicing, RNA biosynthetic processing and transcriptional regulation (Figure 3, Additional file 7: Table S5). These results were consistent with the recently published studies of single-cell RNA sequencing analysis [24, 25].

Figure 3
figure 3

The co-expression of reprogramming genes and the functional classification categories between BBNT and PBNT embryos. Green indicates down-and red up-regulation of gene expression.

The Venn diagrams in Figure 3 show that there were 1218 uniquely upregulated genes in BBNT embryos and 334 upregulated genes in PBNT embryos. The heatmap profile showed that these two gene clusters are significantly diverse (Additional file 8: Figure S3). In the top five upregulated gene clusters, the most significant biological processes in Table 3 were related to nuclear composition, ribosome biogenesis, tRNA/rRNA metabolic processes and mRNA splicing. Significant differences between the mRNA expression profiles were related to transcriptional regulation. It has been demonstrated that transcription regulation is critical for early transcriptional activation during the MZT [26].

Table 3 Functional annotation of specific expressed transcripts between BBNT and PBNT embryos

There were 25 of the transcription regulator-associated genes co-upregulated in both BBNT and PBNT embryos (Figure 4B). One hundred thirty-nine transcripts associated with transcriptional regulation were upregulated in BBNT embryos, whereas these genes had only low expression levels in PBNT embryos (Figure 4A; Additional file 9: Table S6). These genes contribute to multiple biological functions such as general transcription factors, mediator complexes, nuclear receptor subfamily kinase, anchor proteins, RNA polymerases, and zinc finger proteins. As an example, the basic transcription factors for RNA polymerase (e.g., TBP, TFIIB, TAF1D, SP1 and TAF2) were all upregulated in BBNT embryos, but only TBP and TAF2 were upregulated in PBNT embryos (Additional file 10: Figure S4). Unsupervised hierarchical clustering revealed that the expression pattern of these transcription factors in PBNT embryos were more similar to the donor somatic cells, which indicates that the key transcription regulatory pathway in PBNT embryos did not activate as in BBNT embryos (Figure 3).

Figure 4
figure 4

Expression clusters of transcriptional regulated related genes. (A). The 139 unique upregulated transcripts in BBNT embryos that associated with transcriptional regulation. (B). The 25 co-upregulated transcripts of transcription regulation in both BBNT and PBNT embryos. Green represents down- and red up-regulation of gene expression.

Incompatibility between mtDNA and nuclear DNA is regarded as another one of the major problems impairing iSCNT embryo development. Most mitochondrial proteins are encoded in the nucleus, and thus have to be transported into one of the four sub-compartments of the organelle (Figure 5). These processes are mainly mediated by a general translocase in the outer mitochondrial membrane (TOMM: TOMM40), and two distinct translocases in the mitochondrial inner mitochondrial membrane (TIMM: TIMM23 and TIMM22 complexes) [27]. For example, TOMM40 is the channel-forming subunit of the TOMM complex and has an essential role for protein import into the mitochondria. TIMM proteins mediate the import and insertion of hydrophobic membrane proteins into the mitochondrial inner membrane. In this study, the gene transcripts related to the TOMM/TIMM complex in BBNT embryos had the highest expression when compared with the BOs and fibroblasts. Most of the mitochondrial protein import genes in the PBNT embryos such as TOMM40, TIMM10, and TIMM50, had lower expression levels than in the BBNT embryos.Ninety-seven percent of the PBNT embryos did not develop beyond the 8- to 16-cell stage, which is regarded as the crucial time for MZT. The maternal degradation profiles showed that 3822 genes in BBNT embryos and 2650 genes in PBNT embryos were downregulated. There were 2307 genes shared between BBNT and PBNT embryos. Of these genes, 1515 were uniquely downregulated in the BBNT embryos, whereas 343 genes were downregulated in PBNT embryos (Figure 6A). These results suggest that the BBNT embryos were capable of eliciting a more significant degradation of maternal RNA than the PBNT embryos.

Figure 5
figure 5

Cooperation of translocase complexes in mitochondrial protein import and expression comparisons in different cell types.

Figure 6
figure 6

Venn diagram of differentshared and uniquetranscripts in PBNT and BBNT embryos. A. Venn diagrams showing the same down-regulated genes between iSCNT embryosvs bovine oocytesand SCNT embryosvs bovine oocytes (2,307 genes). The blue circle (2,650 genes) indicates the number of down-regulated genes in iSCNT embryosvs the bovine oocytes, and theyellow circle (3,822 genes)indicates the number of down-regulated genes in the SCNT embryosvs the bovine oocytes. B. Venn diagram showing the same non-regulated genes between iSCNT embryosvs Przewalski’s gazelle fibroblasts and SCNT embryosvs bovine fibroblasts (12,513 genes). The blue circle (16,320 genes) indicates the number of non-regulated genes in the iSCNT embryosvs the Przewalski’s gazelle fibroblasts, andyellow circle (15,310 genes)indicates the number of non-regulated genes in the SCNT embryosvs bovine fibroblasts. C. Venn diagram showing common gene expression between iSCNT and SCNT embryos (466 genes). The blue circle denotes transcripts that have similar expression levels between iSCNT embryos vs SCNT embryos (21,185 genes), andthe yellow circle indicates reprogrammed transcriptomes in SCNT embryos (1,524 genes). D. Venn diagram showing somatic gene expression in 8- to 16-cell stage iSCNT embryos (438 genes) compared toSCNT embryos. The blue circle denotes up-regulated genes in iSCNT compared to SCNT embryos (1,100), and yellow circle denotes up-regulated genes in bovine fibroblasts compared to SCNT embryos (4,394).

Theoretically, only donor cell nuclei that completely reprogram in the host cytoplast have a reasonable probability of developing to term following SCNT. After analysis of the donor-specific genes, we found there were more donor-specific transcripts expressed in PBNT embryos than in BBNT embryos (16,320 vs 15,310). The two kinds of embryos shared 12,513 transcripts (Figure 6B). We also observed 438 abnormal fibroblast-specific genes expressed in the PBNT embryos (Figure 6D). The host cytoplast did not silence enough of the donor-specific genes, which is consistent with other reports [28, 29]. For example, the collagen-related protein, Col4A1, was upregulated 2.2 times higher in PBNT than in BBNT embryos. There were 21,855 nonsignificant discrepancy transcripts shared between the two kinds of cloned embryos. When we compared these genes with the SCNT embryonic reprogrammed transcriptome, 466 genes were shared in both datasets (Figure 6C).

Quantitative PCR confirmation of array data

For the validation of microarray datasets results, nine genes showing high levels of significance (TOMM40, TIMM22, TCEB3, ATF1, POLR1C, POLR2B, POLR3A, AMT and NR2F2 were selected and their expression were determined by quantitative PCR (Figure 7). Two polymerase (RNA) II polypeptide genes, POLR2B and POLR3A, showed higher expression in the SCNT embryos with an average of 6.3- and 5.2- fold differences, respectively, when measured with qPCR and 2.8- and 3.2-fold differences, respectively, using microarray measure. The other two genes of mitochondrial translocase, TOMM40 and TIMM22, showed an average of 2.9- and 3.1-fold differences, respectively using qPCR and 2.0- and 1.9-fold difference, respectively, with microarray measure. The activating transcription factor 1 (ATF1) showed highest expression in the SCNT embryos with the average of 10.9- fold difference. Thus, all of the evaluated genes showed similar patterns of mRNA abundance in microarray analysis and qPCR measure.

Figure 7
figure 7

Quantification of relative abundance of nine transcripts between the 8–16 cell stage of the embryos (72h post activation). Przewalski’s gazelle-bovine interspecies SCNT embryos (P-B) and bovine intra-species SCNT embryos(B-B) in expression of nine select genes as determined by quantitative PCR. Data are presented as mean value of triplicate measurements including standard errors.


Various studies have shown that the majority of cloned mammalian embryos fail to development because of incomplete genomic reprogramming [30]. Theoretically, the iSCNT embryos derived from inter-class, inter-order or inter-family donor cells and recipient oocytes should result in much more incomplete nuclear reprogramming than that of intra-species cloning because of greater genetic divergence. To assist in the relaxation of chromatin structure and provide for better cell nuclear reprogramming, histone-deacetylase inhibitors and DNA methyl-transferase have been used in SCNT protocols. TSA, VPA and 5-aza-2′-deoxycytidine have been reported to increase cloned blastocyst development rates significantly [21, 3133]. However treatments with these small molecules do not enhance full-term in vivo development [34]. In iSCNT studies, the use of TSA has not been shown to have any beneficial effect on cloned embryo development in the gaur-bovine [35], human-rabbit [36], and dog-porcine [37]. In our study, the PBNT morula/blastocyst development also did not improve when using TSA, VPA, or Oct 4-transgenic treatments. The coexistence of the donor cell nucleus and oocyte nucleus using the RNT protocol has been reported to facilitate donor cell reprogramming and improve cloned embryo development in mouse [38] and bovine [39]. In our study, PBNT hybrid embryo development was not improved with the RNT protocol.

Genomic incompatibility between the nuclear donor cell and the cytoplast is a major contributing factor causing iSCNT hybrid embryos to arrest at the 8- to 16-cell stage. Wang et al. [17] observed in chimpanzee–bovine iSCNT embryos that the bovine ooplasm partially remodeled the chimpanzee somatic cell nuclei. However, the cloned embryos still arrested at the 8- to 16-cell stages. Lagutina et al. [40] reported that pig–bovine iSCNT embryo development was completely blocked at the 16-cell stage, while the bovine–pig iSCNT embryos stopped earlier at the 4-cell stage. These studies demonstrated that embryo development was completely arrested before ZGA.

Cytoplasmic factors begin the genomic reprogramming process shortly after the somatic cell nucleus is transferred into the oocyte. It is generally agreed that most cloned embryos are not fully activated and that incomplete cellular reprogramming is the major cause of the extremely poor outcome in SCNT and iSCNT experiments. The timing of large-scale activation of the embryonic genome is species-specific. In the bovine embryo it occurs at the 8-cell stage [41]. At the beginning of MZT, the embryo is dependent on a vast store of mRNA and proteins in the recipient cytoplasm. Complete EGA occurs in the later preimplantation development stage. At the MZT stage the embryo should completely reprogram the somatic cell nucleus, silencing the accompanying donor-specific transcripts and then begin its self-sustained cellular directed development [42]. Complete activation of the embryonic genome is required to produce a normal full-term cloned fetus [43]. In inter-specific cloning, the reprogramming process is greatly compromised and becomes an extremely complicated disarray when the two genetically different genotypes are constrained to coexist in the same embryo. We observed 309 genes that showed increased co-expression in the BBNT and PBNT embryos. These normally control important biological processes such as the formation of membrane-enclosed lumens, RNA processing/splicing, RNA biosynthetic processes and transcriptional regulation. There were significant differences in the upregulated genes associated with nuclear reprogramming between PBNT and BBNT embryos. Greater reprogramming occurred in the BBNT embryos than in the PBNT embryos. Our findings are consistent with the most recently published single-cell RNA sequencing analysis, which identified three pathways in human 8-cell stage embryos [24, 25], and they are also consistent with transcriptional activation waves based on a single human embryo study [26].

Communication between the nucleus and cytoplasm is critical for ZGA or EGA, and is essential for normal embryo development [26, 44]. An orderly and complete transition from the maternal to the zygotic genome is crucial for a normal term pregnancy. Irregular or incomplete activation of the embryonic genome will alter the processes required to direct normal embryo development. During MZT, the oocyte’s maternal mRNAs gradually undergo degradation along with gradual activation of the zygotic genes [45]. In this study, maternal degradation profiles differed between PBNT embryos with 2650 genes and BBNT embryos with 3822 genes;, 1515 genes were uniquely downregulated in BBNT embryos whereas there were 343 in PBNT embryos. There were more gazelle cell-specific transcripts detected in PBNT than in BBNT embryos. This observation was consistent with the presence of 438 abnormal fibroblast-specific genes expressed in PBNT embryos. These results suggested that the BBNT embryos were more efficient in degrading maternal RNA than the PBNT embryos, and that the bovine cytoplasts were not very effective in silencing the gazelle cell-specific genes. The inability to silence the reprogramming-associated genes in PBNT embryos might in part arise from various nuclear-cytoplasmic incompatibilities, thus leading to failure of subsequent embryo development [28, 29]. Past studies on SCNT embryos have established a direct relationship between embryonic genome activation and the ability of the embryo to develop [46, 47]. Aberrations in the expression of housekeeping genes and genes dependent on the major embryonic genome activation is a major problem of iSCNT embryos for different species combinations [42]. Human–human, human–bovine, and human–rabbit interspecific clones all had similar rates of development to the 8- to 16-cell stage. Gene expression profiles of human–human embryos had similar levels of upregulation as the control fertilized embryos, which is not the case for interspecies clones [18, 42].

The extremely low developmental rates of PBNT embryos in this study strongly suggest that aberrant gene expression and altered transcriptional regulation are major contributing factors leading to poor embryo development. Even though there were 25 transcription regulators co-upregulated in the BBNT and PBNT embryos, more than 100 genes involved in transcriptional regulation were widely expressed in BBNT embryos but not in the PBNT embryos. The unsupervised hierarchical clustering utilized in this study showed that expression patterns of the transcription factors in PBNT embryos were more closely aligned to the donor somatic cell nucleus, which indicates in all likelihood that the key transcription regulatory pathway of PBNT was not activated as it was in BBNT embryos. In chimpanzee–bovine iSCNT experiment where the transcription factors Oct-4, Stella, Crabp1, CCNE2, CXCL6, PTGER4, H2AFZ, c-MYC, KLF4, and GAPDH were expressed, the chimpanzee nuclei were only partially remodeled and the cloned embryos did not develop beyond the 8- to 16-cell stages [17]. Even when using Oct-4 transgenic cells, blastocyst development was also not improved in our PBNT embryos. It appears that the downregulation of genes associated with transcription and transcription regulation was only partially reprogrammed in PBNT embryos, which was not enough to support normal zygotic genome activation.

Another major obstacle preventing nuclear–cytoplasmic compatibility is the interaction between donor cell nuclear DNA, mtDNA and oocyte mtDNA [4851]. An iSCNT study has indicated that the donor cell mtDNA gradually becomes dominant in the post-implantation embryo and during fetal development [35]. Mitochondria promote a broad range of critical functions and have to provide sufficient energy for reprogramming and embryonic genome activation during embryonic and fetal development [40, 48]. Mitochondrial impairment or insufficiency has been reported to impair the efficacy of SCNT [52]. Our data revealed that the TOMM/TIMM complex-associated genes in BBNT embryos had significantly higher expression levels than in PBNT embryos. It has been estimated that over 1000 proteins involved in mitochondrial function are encoded by nuclear DNA, synthesized in the cytosol, and then targeted and imported into the mitochondria by specific transfer pathways [5254]. Any alterations in the mitochondrial protein import pathway resulting from incompatibilities of heterogeneous mtDNA and genome DNA would likely lead to developmental failure of PBNT embryos. Thus, inefficiency in mitochondrial protein importation mechanisms might be a major contributor to iSCNT embryo failure.


The gene expression network determines the cell’s identity and behavior. The latest transcriptome profiles from single-cell RNA Sequencing provides convincing evidence that the human preimplantation transcriptional organization is highly preserved, that gene activation is sequentially ordered and genetic programming is essential for preimplantation development. Our study provided a comprehensive comparison of intra/inter bovine embryonic transcriptome during preimplantation development. Significant differences were observed in the mRNA expression profiles between the BBNT and PBNT cloned embryos. Comparisons of the reprogrammed transcriptomes identified major differences in expressed genes and transcriptional regulation, some of which are most likely first stage-specific modules that are necessary for embryonic genome activation. A failure of these transcriptional regulatory pathways would hinder MZT. The inefficiency of mitochondrial protein import might also act as another cause for the developmental failure of iSCNT embryos. Further functional analysis studies on the pathways and genes essential for early embryo development are needed if the efficiency of SCNT and iSCNT is to be improved.


Ethics statement

All the bovine Oocytes and embryos were handled and studies were carried out according to the guidelines of The Inner Mongolia University Animal Care and Use Committee. The bovine ovaries used in this study were collected with permission of the Hohhot slaughterhouse. The animal protocol was approved by The Animal Care and Use Committee of Inner Mongolia University. A small piece of ear tissue of an adult Przewalski’s gazelle was collected in Qinghai Wildlife Garden (Xining) with the permission of Qinghai Forestry Bureau.

The production of cloned embryos using classic SCNT and reverse nuclear transfer (RNT) protocols

The culture of Przewalski’s gazelle fibroblasts and the in vitro maturation of bovine oocytes were described in our previous reports [55, 56]. The donor cells were transferred to the perivitelline space of enucleated oocytes by the classic SCNT [55, 57] and RNT protocols was described by Meng et al. [39]. Successful enucleation was confirmed by UV illumination of the Hoechst-stained karyoplasts as described above. The enucleated oocytes were allowed to recover for 30 min in an incubator and then activated with ionomycin and CHX treatments. The resulted embryos were cultured for further development.

The reconstructed couplets were electrically fused and then chemically activated by ionomycin and cycloheximide. The activated cloned embryos were cultured in SOFaa in the presence of VPA or TSA at 38.5°C in a humidified atmosphere containing 5% CO2. Cleavage and embryo development were observed with a light microscope at 100X. A more detail description of the SCNT and iSCNT protocols can be found in the supplemental materials (Additional file 11).

Microarray dataset

The Affymetrix Gene Chip Bovine Genome array contains approximately 23,000 transcripts including assemblies from ~19,000 UniGene clusters. The array images were first quantified using Gene Chip Operating Software (GCOS, Affymetrix). The biological replicates of the datasets have high reproducibility.

Expression pattern analysis

Genes with similar expression patterns are likely to have functional correlations, therefore we performed a cluster analysis of the gene expression patterns using Cluster 3.0 [58] and JavaTreeview [59] software. Expression differences were clustered by the Hierarchical Complete Linkage Clustering method using an uncentered correlation similarity matrix. R packages were used for the Venn diagram and expression analysis of different transcripts.

Gene ontology analysis

Functional annotation was performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resource [60]. Gene ontology terms shown in this study summarized all similar sub-terms into an overarching term, and Benjamani-Hochberg adjusted P values were shown for the representative term.

Statistical analysis

Differences in embryo development between experimental groups were analyzed using Student’s t-test and χ2 analysis, SPSS16.0 was used for statistical analysis. P < 0.05 was considered statistically significant.


  1. Group ISAS: Specialist Group: Procapra przewalskii. IUCN 2008. IUCN Red List of Threatened Species. Retrieved 10 December 2008. Database entry includes justification for why this species is Endangered. 2008

    Google Scholar 

  2. Leslie DM, Groves CP, Abramov AV: Procapra przewalskii (Artiodactyla: Bovidae). Mamm Species. 2010, 42: 124-137.

    Article  Google Scholar 

  3. Li Z, Jiang Z, Beauchamp G: Nonrandom mixing between groups of Przewalski’s gazelle and Tibetan gazelle. J Mammalogy. 2010, 91: 674-680. 10.1644/09-MAMM-A-203.1.

    Article  Google Scholar 

  4. Jiang Z, Lei R, Liu B, Li C: A review on the researches of Przewalski’s gazelle. Chinese J Zool. 2002, 38: 129-132.

    Google Scholar 

  5. Liu B, Jiang Z: Impacts of grassland fencing on plant communities and conservation of a rare gazelle, the Przewalski’s gazelle. Chin Biodiversity. 2001, 10: 326-331.

    Google Scholar 

  6. Chen D, Sun Q, Liu J, Li G, Lian L, Wang M, Han Z, Song X, Li J, Sun Q: The giant panda (Ailuropoda melanoleuca) somatic nucleus can dedifferentiate in rabbit ooplasm and support early development of the reconstructed egg. Sci Chin Ser C: Life Sci. 1999, 42: 346-353. 10.1007/BF02882053.

    Article  CAS  Google Scholar 

  7. Chen D, Han Z, Ziyu Z, Ruicheng W, Zhonghua L, Minkang W, Li L, Jinsong L, Jun D, Pengyan W, Hemin Z: The effect of series nuclear transfer protocol on the embryonic development of the giant panda interpecies cloning. Chin Sci Bulletin. 2001, 46: 1899-1901.

    Google Scholar 

  8. Zhao ZJ, Li RC, Cao HH, Zhang QJ, Jiang MX, Ouyang YC, Nan CL, Lei ZL, Song XF, Sun QY, Chen DY: Interspecies nuclear transfer of Tibetan antelope using caprine oocyte as recipient. Mol Reprod Dev. 2007, 74: 412-419. 10.1002/mrd.20608.

    Article  CAS  PubMed  Google Scholar 

  9. Sansinena MJ, Hylan D, Hebert K, Denniston RS, Godke RA: Banteng (Bos javanicus) embryos and pregnancies produced by interspecies nuclear transfer. Theriogenology. 2005, 63: 1081-1091. 10.1016/j.theriogenology.2004.05.025.

    Article  CAS  PubMed  Google Scholar 

  10. Li Y, Dai Y, Du W, Zhao C, Wang L, Wang H, Liu Y, Li R, Li N: In vitro development of yak (Bos grunniens) embryos generated by interspecies nuclear transfer. Anim Reprod Sci. 2007, 101: 45-59. 10.1016/j.anireprosci.2006.09.018.

    Article  PubMed  Google Scholar 

  11. Song J, Hua S, Song K, Zhang Y: Culture, characteristics and chromosome complement of Siberian tiger fibroblasts for nuclear transfer. In Vitro Cell Dev Biol Anim. 2007, 43: 203-209.

    Article  PubMed  Google Scholar 

  12. Bhuiyan MMU, Suzuki Y, Watanabe H, Lee E, Hirayama H, Matsuoka K, Fujise Y, Ishikawa H, Ohsumi S, Fukui Y: Production of sei whale (Balaenoptera borealis) cloned embryos by inter-and intra-species somatic cell nuclear transfer. J Reprod Dev. 2010, 56: 131-139. 10.1262/jrd.09-158A.

    Article  CAS  PubMed  Google Scholar 

  13. Loi P, Ptak G, Barboni B, Fulka J, Cappai P, Clinton M: Genetic rescue of an endangered mammal by cross-species nuclear transfer using post-mortem somatic cells. Nat Biotech. 2001, 19: 962-964. 10.1038/nbt1001-962.

    Article  CAS  Google Scholar 

  14. Yang CY, Li RC, Pang CY, Yang BZ, Qin GS, Chen MT, Zhang XF, Huang FX, Zheng HY, Huang YJ, Liang XW: Study on the inter-subspecies nuclear transfer of river buffalo somatic cell nuclei into swamp buffalo oocyte cytoplasm. Anim Reprod Sci. 2010, 121: 78-83. 10.1016/j.anireprosci.2010.05.011.

    Article  CAS  PubMed  Google Scholar 

  15. Amarnath D, Choi I, Moawad AR, Wakayama T, Campbell KH: Nuclear-cytoplasmic incompatibility and inefficient development of pig-mouse cytoplasmic hybrid embryos. Reproduction. 2011, 142: 295-307. 10.1530/REP-11-0044.

    Article  CAS  PubMed  Google Scholar 

  16. Bowles EJ, Lee JH, Alberio R, Lloyd RE, Stekel D, Campbell KH, St John JC: Contrasting effects of in vitro fertilization and nuclear transfer on the expression of mtDNA replication factors. Genetics. 2007, 176: 1511-1526. 10.1534/genetics.106.070177.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Wang K, Beyhan Z, Rodriguez RM, Ross PJ, Iager AE, Kaiser GG, Chen Y, Cibelli JB: Bovine ooplasm partially remodels primate somatic nuclei following somatic cell nuclear transfer. Clon Stem Cell. 2009, 11: 187-202. 10.1089/clo.2008.0061.

    Article  CAS  Google Scholar 

  18. Chung Y, Bishop CE, Treff NR, Walker SJ, Sandler VM, Becker S, Klimanskaya I, Wun WS, Dunn R, Hall RM, Su J, Lu SJ, Maserati M, Choi YH, Scott R, Atala A, Dittman R, Lanza R: Reprogramming of human somatic cells using human and animal oocytes. Clon Stem Cell. 2009, 11: 213-223. 10.1089/clo.2009.0004.

    Article  CAS  Google Scholar 

  19. Narbonne P, Simpson DE, Gurdon JB: Deficient induction response in a Xenopus nucleocytoplasmic hybrid. PLoS Biol. 2011, 9: e1001197-10.1371/journal.pbio.1001197.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Esteves TC, Balbach ST, Pfeiffer MJ: Araúzo‒Bravo MJ, Klein DC, Sinn M, Boiani M: Somatic cell nuclear reprogramming of mouse oocytes endures beyond reproductive decline. Aging Cell. 2011, 10: 80-95. 10.1111/j.1474-9726.2010.00644.x.

    Article  CAS  PubMed  Google Scholar 

  21. Xu W, Wang Y, Li Y, Wang L, Xiong X, Su J, Zhang Y: Valproic acid improves the in vitro development competence of bovine somatic cell nuclear transfer embryos. Cell Reprogram. 2012, 14: 138-145.

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Akagi S, Matsukawa K, Mizutani E, Fukunari K, Kaneda M, Watanabe S, Takahashi S: Treatment with a histone deacetylase inhibitor after nuclear transfer improves the preimplantation development of cloned bovine embryos. J Reprod Develop. 2011, 57: 120-126. 10.1262/jrd.10-058A.

    Article  CAS  Google Scholar 

  23. Sun L, Wu KL, Zhang D, Wang HY, Wang Y, Xu ZY, Huang XY, Chen ZJ, Liu HQ: Increased cleavage rate of human nuclear transfer embryos after 5-aza-2’-deoxycytidine treatment. Reprod Biomed Online. 2012, 25: 425-433. 10.1016/j.rbmo.2012.06.018.

    Article  CAS  PubMed  Google Scholar 

  24. Xue Z, Huang K, Cai C, Cai L, Jiang CY, Feng Y, Liu Z, Zeng Q, Cheng L, Sun YE, Liu JY, Horvath S, Fan G: Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013, 500: 593-597. 10.1038/nature12364.

    Article  CAS  PubMed  Google Scholar 

  25. Yan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J, Huang J, Li M, Wu X, Wen L, Lao K, Li R, Qiao J, Tang F: Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol. 2013, 20: 1131-1139. 10.1038/nsmb.2660.

    Article  CAS  PubMed  Google Scholar 

  26. Vassena R, Boue S, Gonzalez-Roca E, Aran B, Auer H, Veiga A, Izpisua Belmonte JC: Waves of early transcriptional activation and pluripotency program initiation during human preimplantation development. Development. 2011, 138: 3699-3709. 10.1242/dev.064741.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Bauer MF, Hofmann S, Neupert W, Brunner M: Protein translocation into mitochondria: the role of TIM complexes. Trends Cell Biol. 2000, 10: 25-31. 10.1016/S0962-8924(99)01684-0.

    Article  CAS  PubMed  Google Scholar 

  28. Wang K, Otu HH, Chen Y, Lee Y, Latham K, Cibelli JB: Reprogrammed transcriptome in rhesus-bovine interspecies somatic cell nuclear transfer embryos. PLoS One. 2011, 6: e22197-10.1371/journal.pone.0022197.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Vassena R, Han Z, Gao S, Baldwin DA, Schultz RM, Latham KE: Tough beginnings: alterations in the transcriptome of cloned embryos during the first two cell cycles. Dev Biol. 2007, 304: 75-89. 10.1016/j.ydbio.2006.12.015.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Narbonne P, Miyamoto K, Gurdon JB: Reprogramming and development in nuclear transfer embryos and in interspecific systems. Curr Opin Genet Dev. 2012, 22: 450-458. 10.1016/j.gde.2012.09.002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Wen DC, Bi CM, Xu Y, Yang CX, Zhu ZY, Sun QY, Chen DY: Hybrid embryos produced by transferring panda or cat somatic nuclei into rabbit MII oocytes can develop to blastocyst in vitro. J Exp Zool A Comp Exp Biol. 2005, 303: 689-697.

    Article  PubMed  Google Scholar 

  32. Wee G, Shim JJ, Koo DB, Chae JI, Lee KK, Han YM: Epigenetic alteration of the donor cells does not recapitulate the reprogramming of DNA methylation in cloned embryos. Reproduction. 2007, 134: 781-787. 10.1530/REP-07-0338.

    Article  CAS  PubMed  Google Scholar 

  33. Enright BP, Kubota C, Yang X, Tian XC: Epigenetic characteristics and development of embryos cloned from donor cells treated by trichostatin A or 5-aza-2’-deoxycytidine. Biol Reprod. 2003, 69: 896-901. 10.1095/biolreprod.103.017954.

    Article  CAS  PubMed  Google Scholar 

  34. Sawai K, Fujii T, Hirayama H, Hashizume T, Minamihashi A: Epigenetic status and full-term development of bovine cloned embryos treated with trichostatin A. J Reprod Dev. 2012, 58: 302-309. 10.1262/jrd.2011-020.

    Article  CAS  PubMed  Google Scholar 

  35. Srirattana K, Matsukawa K, Akagi S, Tasai M, Tagami T, Nirasawa K, Nagai T, Kanai Y, Parnpai R, Takeda K: Constant transmission of mitochondrial DNA in intergeneric cloned embryos reconstructed from swamp buffalo fibroblasts and bovine ooplasm. Anim Sci J. 2011, 82: 236-243. 10.1111/j.1740-0929.2010.00827.x.

    Article  CAS  PubMed  Google Scholar 

  36. Shi LH, Miao YL, Ouyang YC, Huang JC, Lei ZL, Yang JW, Han ZM, Song XF, Sun QY, Chen DY: Trichostatin A (TSA) improves the development of rabbit-rabbit intraspecies cloned embryos, but not rabbit-human interspecies cloned embryos. Dev Dyn. 2008, 237: 640-648. 10.1002/dvdy.21450.

    Article  PubMed  Google Scholar 

  37. Sugimura S, Narita K, Yamashiro H, Sugawara A, Shoji T, Terashita Y, Nishimori K, Konno T, Yoshida M, Sato E: Interspecies somatic cell nucleus transfer with porcine oocytes as recipients: A novel bioassay system for assessing the competence of canine somatic cells to develop into embryos. Theriogenology. 2009, 72: 549-559. 10.1016/j.theriogenology.2009.04.011.

    Article  CAS  PubMed  Google Scholar 

  38. Wakayama S, Cibelli JB, Wakayama T: Effect of timing of the removal of oocyte chromosomes before or after injection of somatic nucleus on development of NT embryos. Clon Stem Cell. 2003, 5: 181-189. 10.1089/153623003769645848.

    Article  CAS  Google Scholar 

  39. Meng Q, Bai C, Liu Y, Wu X, Bunch TD, Li GP: In vitro development and chromosomal configuration of bovine somatic cloned embryos with nonenucleated metaphase II oocytes. Cell Reprogram. 2010, 12: 481-490. 10.1089/cell.2009.0114.

    Article  CAS  PubMed  Google Scholar 

  40. Lagutina I, Fulka H, Brevini TA, Antonini S, Brunetti D, Colleoni S, Gandolfi F, Lazzari G, Fulka J, Galli C: Development, embryonic genome activity and mitochondrial characteristics of bovine-pig inter-family nuclear transfer embryos. Reproduction. 2010, 140: 273-285. 10.1530/REP-09-0578.

    Article  CAS  PubMed  Google Scholar 

  41. Kues W, Sudheer S, Herrmann D, Carnwath J, Havlicek V, Besenfelder U, Lehrach H, Adjaye J, Niemann H: Genome-wide expression profiling reveals distinct clusters of transcriptional regulation during bovine preimplantation development in vivo. PNAS. 2008, 105: 19768-19773. 10.1073/pnas.0805616105.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Lagutina I, Zakhartchenko V, Fulka H, Colleoni S, Wolf E, Fulka J, Lazzari G, Galli C: Formation of nucleoli in interspecies nuclear transfer embryos derived from bovine, porcine, and rabbit oocytes and nuclear donor cells of various species. Reproduction. 2011, 141: 453-465. 10.1530/REP-10-0266.

    Article  CAS  PubMed  Google Scholar 

  43. Alizadeh Z, Kageyama S, Aoki F: Degradation of maternal mRNA in mouse embryos: selective degradation of specific mRNAs after fertilization. Mol Reprod Dev. 2005, 72: 281-290. 10.1002/mrd.20340.

    Article  CAS  PubMed  Google Scholar 

  44. Minami N, Suzuki T, Tsukamoto S: Zygotic gene activation and maternal factors in mammals. J Reprod Dev. 2007, 53: 707-715. 10.1262/jrd.19029.

    Article  CAS  PubMed  Google Scholar 

  45. Shin SW, Shimizu N, Tokoro M, Nishikawa S, Hatanaka Y, Anzai M, Hamazaki J, Kishigami S, Saeki K, Hosoi Y, Iritani A, Murata S, Matsumoto K: Mouse zygote-specific proteasome assembly chaperone important for maternal-to-zygotic transition. Biol Open. 2013, 2: 170-182. 10.1242/bio.20123020.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Maddox-Hyttel P, Svarcova O, Laurincik J: Ribosomal RNA and nucleolar proteins from the oocyte are to some degree used for embryonic nucleolar formation in cattle and pig. Theriogenology. 2007, 68 (Suppl 1): S63-S70.

    Article  CAS  PubMed  Google Scholar 

  47. Laurincik J, Maddox-Hyttel P: Nucleolar remodeling in nuclear transfer embryos. Adv Exp Med Biol. 2007, 591: 84-92. 10.1007/978-0-387-37754-4_6.

    Article  PubMed  Google Scholar 

  48. Jiang Y, Kelly R, Peters A, Fulka H, Dickinson A, Mitchell DA, St John JC: Interspecies somatic cell nuclear transfer is dependent on compatible mitochondrial DNA and reprogramming factors. PLoS One. 2011, 6: e14805-10.1371/journal.pone.0014805.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Hua S, Zhang Y, Song K, Song J, Zhang Z, Zhang L, Zhang C, Cao J, Ma L: Development of bovine-ovine interspecies cloned embryos and mitochondria segregation in blastomeres during preimplantation. Anim Reprod Sci. 2008, 105: 245-257. 10.1016/j.anireprosci.2007.03.002.

    Article  CAS  PubMed  Google Scholar 

  50. Gomez MC, Pope CE, Ricks DM, Lyons J, Dumas C, Dresser BL: Cloning endangered felids using heterospecific donor oocytes and interspecies embryo transfer. Reprod Fertil Dev. 2009, 21: 76-82.

    Article  CAS  PubMed  Google Scholar 

  51. Sansinena MJ, Lynn J, Bondioli KR, Denniston RS, Godke RA: Ooplasm transfer and interspecies somatic cell nuclear transfer: heteroplasmy, pattern of mitochondrial migration and effect on embryo development. Zygote. 2011, 19: 147-156. 10.1017/S0967199410000419.

    Article  CAS  PubMed  Google Scholar 

  52. Hiendleder S, Zakhartchenko V, Wolf E: Mitochondria and the success of somatic cell nuclear transfer cloning: from nuclear-mitochondrial interactions to mitochondrial complementation and mitochondrial DNA recombination. Reprod Fertil Dev. 2005, 17: 69-83.

    Article  CAS  PubMed  Google Scholar 

  53. Taylor SW, Fahy E, Zhang B, Glenn GM, Warnock DE, Wiley S, Murphy AN, Gaucher SP, Capaldi RA, Gibson BW, Ghosh SS: Characterization of the human heart mitochondrial proteome. Nat Biotechnol. 2003, 21: 281-286. 10.1038/nbt793.

    Article  CAS  PubMed  Google Scholar 

  54. Cotter D, Guda P, Fahy E, Subramaniam S: MitoProteome: mitochondrial protein sequence database and annotation system. Nucleic Acids Res. 2004, 32: D463-D467. 10.1093/nar/gkh048.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Li GP, White KL, Aston KI, Meerdo LN, Bunch TD: Conditioned medium increases the polyploid cell composition of bovine somatic cell nuclear-transferred blastocysts. Reproduction. 2004, 127: 221-228. 10.1530/rep.1.00089.

    Article  CAS  PubMed  Google Scholar 

  56. Bai C, Liu H, Liu Y, Wu X, Cheng L, Bou S, Li GP: Diploid oocyte formation and tetraploid embryo development induced by cytochalasin B in bovine. Cell Reprogram. 2011, 13: 37-45. 10.1089/cell.2010.0038.

    Article  CAS  PubMed  Google Scholar 

  57. Li GP, Bunch TD, White KL, Aston KI, Meerdo LN, Pate BJ, Sessions BR: Development, chromosomal composition, and cell allocation of bovine cloned blastocyst derived from chemically assisted enucleation and cultured in conditioned media. Mol Reprod Dev. 2004, 68: 189-197. 10.1002/mrd.20071.

    Article  CAS  PubMed  Google Scholar 

  58. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. PNAS. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Saldanha AJ: Java Treeview–extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.

    Article  CAS  PubMed  Google Scholar 

  60. da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.

    Article  PubMed  Google Scholar 

Download references


We would like to thank Dr. Thomas D. Bunch, Department of Animal, Dairy and Veterinary Sciences, Utah State University for his critical reading and edition of this manuscript. This work was partly supported by the National Basic Research Program of China (973, 2012CB22306), the National Natural Science Foundation of China (31372289), the High-level Scientific Research Foundation Award for introduction of talent, Inner Mongolia University (No.115115), the Specialized Research Fund for the Doctoral Program of Higher Education (20131501120009) and the Natural Science Foundation of Inner Mongolia Autonomous Region (No:2013MS0503); as well as the support of the Program for Chang jiang Scholars and Innovative Research Team in Inner Mongolia University.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Shorgan Bou or Guangpeng Li.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GPL and SB designed the study, YCZ, YG, GHS, ZYW, KL and CLB performed the experiments, YCZ, GPL, GHS and QZL analyzed the data, YCZ and GPL prepared the manuscript. All authors read and approved the manuscript.

Yongchun Zuo, Yu Gao, Guanghua Su contributed equally to this work.

Electronic supplementary material


Additional file 1: Table S1: The Cleavage rates, 8-16-cell, morula and blastocyst development of PBNT embryos derived from Oct-4-eGFP transgenic cells. (DOC 34 KB)


Additional file 2: Figure S1: Photomicrographs of the Przewalski’s gazelle-bovine interspecies SCNT embryos and bovine intra-species SCNT embryos cultured in vitro. The figures of A, B, C, D are different development ages of Przewalski’s gazelle-bovine interspecies SCNT embryos. The figures of A’, B’, C’, D’ are different development ages of bovine intra-species SCNT embryos. A, A’ are 2–4 cell stage of the embryos (36h post activation). B, B’ are 8–16 cell stage of the embryos (72 h post activation). C, C’ are morula stage of the embryos (120 h post activation). D, D’ are blastocyst stage of the embryos (160 h post activation). Each scale bar represents 100 μm. (DOC 2 MB)


Additional file 3: Table S2: The Spearman correlation coefficients between the profiles of different cell types reflect the degree of change in transcriptomes. (XLS 20 KB)


Additional file 4: Table S3: The expression of 10,000 transcripts in different cell types. BO: bovine oocytes, BBNT: bovine SCNT 8- to 16-cell embryos, PBNT: iSCNT 8- to 16-cell embryos, PC: Przewalski’s gazelle, BC: bovine fibroblast cells. (XLS 1 MB)


Additional file 5: Figure S2: Hierarchical clusters of reprogramming related gene expression profiles. Green indicates the down-and red the up-regulated gene expression. Venn diagram of shared and unique genes among different transcriptomes. (JPEG 3 MB)


Additional file 6: Table S4: The expression of reprogramming related genes in the BBNT(1527) and PBNT(643) embryos. (XLS 2 MB)


Additional file 7: Table S5: The GO annotation of 309 reprogramming related genes that expressed both in BBNT and PBNT embryos. (XLS 42 KB)


Additional file 8: Figure S3: The heatmap profile of the reprogramming related gene that uniquely upexpressed in BBNT embryos and PBNT embryos, respectively. There are respectively 1,218 uniquely up-regulated genes occurred in BBNT embryos and 334 up-regulated genes in PBNT embryos. (PDF 403 KB)


Additional file 9: Table S6: The expression of transcription regulation related genes. One hundred thirty-nine transcripts related to transcription regulation were up-regulated in BBNT embryos; whereas, these genes were down-regulated in PBNT embryos. (XLS 218 KB)


Additional file 10: Figure S4: The expression diversity for basal transcription factors of RNA polymerase in two different embryos. (JPEG 1 MB)

Additional file 11: The detail description of the SCNT and iSCNT protocols in this study.(DOC 30 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zuo, Y., Gao, Y., Su, G. et al. Irregular transcriptome reprogramming probably causes thec developmental failure of embryos produced by interspecies somatic cell nuclear transfer between the Przewalski’s gazelle and the bovine. BMC Genomics 15, 1113 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: