- Research article
- Open Access
Specific gene-regulation networks during the pre-implantation development of the pig embryo as revealed by deep sequencing
BMC Genomics volume 15, Article number: 4 (2014)
Because few studies exist to describe the unique molecular network regulation behind pig pre-implantation embryonic development (PED), genetic engineering in the pig embryo is limited. Also, this lack of research has hindered derivation and application of porcine embryonic stem cells and porcine induced pluripotent stem cells (iPSCs).
We identified and analyzed the genome wide transcriptomes of pig in vivo-derived and somatic cell nuclear transferred (SCNT) as well as mouse in vivo-derived pre-implantation embryos at different stages using mRNA deep sequencing. Comparison of the pig embryonic transcriptomes with those of mouse and human pre-implantation embryos revealed unique gene expression patterns during pig PED. Pig zygotic genome activation was confirmed to occur at the 4-cell stage via genome-wide gene expression analysis. This activation was delayed to the 8-cell stage in SCNT embryos. Specific gene expression analysis of the putative inner cell mass (ICM) and the trophectoderm (TE) revealed that pig and mouse pre-implantation embryos share regulatory networks during the first lineage segregation and primitive endoderm differentiation, but not during ectoderm commitment. Also, fatty acid metabolism appears to be a unique characteristic of pig pre-implantation embryonic development. In addition, the global gene expression patterns in the pig SCNT embryos were different from those in in vivo-derived pig embryos.
Our results provide a resource for pluripotent stem cell engineering and for understanding pig development.
Embryonic stem cells (ESCs) can differentiate into any cell type of the three germ layers as well as into male and female germ cells [1, 2]. As such, they offer great potential for regenerative medicine and animal breeding. ESCs have been derived from the mouse, rat and humans, but only mouse ESCs have been tested in tetraploid complementation assays [3–6]. The derivation of such cell lines later helped researchers understand the molecular mechanisms governing pluripotency and early embryonic cell fate commitment. However, because of the short lifespan of the mouse, mouse models are insufficient for evaluating the long-term effects of cell replacement or cell therapy. On the other hand, the in vivo developmental potential of ESCs cannot be directly tested in humans. Pigs are an ideal model for preclinical development and design of therapeutic approaches because their organs are morphologically and functionally similar to humans [7, 8]. For this reason, pig ES cell lines must be generated with the same in vivo developmental potential as mouse ES cells. Since the 1990s, attempts have been made to derive pig ESCs and the generation of porcine iPSCs has been recently reported. However, during this time, no stable porcine cell lines have been capable of germ-line transmission; thus, they are not authentically pluripotent [9–12]. These problems may be due to inadequacies in the currently used in vitro culture conditions that cannot support pluripotency maintenance. This is partly a result of a lack of information regarding unique molecular mechanisms of early pig embryonic development [8, 13].
The progression from fertilization to implantation among mammals is highly conserved, and the morphologic stages are similar . Few interspecies differences do occur, such as time spent at each stage, timing of zygotic genome activation (ZGA) and cell lineage commitment initiation. Gene-regulation networks of mouse PEDs have been extensively studied and reported [15–17], but scarce information regarding molecular mechanism of pig early embryonic development as well as other large domestic animals has limited our knowledge of developmental biology and aspects of engineering their stem cells.
Genome-wide transcriptome analysis may reveal unique gene regulation networks during PED that would be useful in the biological studies of undifferentiated ESCs and pre-implantation embryos [18, 19]. Transcriptome analysis of early pig embryos may also elucidate differentiation characteristics of putative porcine ESCs and iPSCs to optimize in vitro culture conditions for the generation of true pig ESCs. Therefore, we compared early pig in vivo fertilized-derived and mouse in vivo fertilized-derived embryo transcriptomes and mapped a putative gene regulation network during pig PED. This work represents a significant step towards characterizing normal and cloned pig early embryos using genome-wide gene expression patterns.
Dynamic gene expression landscapes of PED
To identify the gene regulation networks that act during in vivo pig PED, we isolated mRNA for deep sequencing from in vivo porcine (Pnm) and mouse embryos (Ms) and porcine SCNT embryos (Pnt) harvested at different stages during PED (Additional file 1: Table S1 and Figure 1A). At the blastocyst stage, we physically separated the blastocyst into two parts, one containing pure trophectoderm (TE) and the other (ICMTE) containing the inner cell mass (ICM) and part of the TE (Additional file 2: Figure S1). The protocol for small cell number mRNA-sequencing was optimized to identify and analyze pig and mouse genome-wide transcriptomes at various stages of early development . To verify the sequencing data reproducibility, we collected two sets of pig in vivo embryos at each stage (Additional file 1: Table S1). The Pearson correlation coefficient for replicates calculated by log10 RPKM ranged from 0.865 to 0.985, indicating reliable sequencing data. Representative results observed in oocyte samples are shown in Figure 1B. To further validate the data and analysis methods, we performed three correlation measurements to estimate the transcriptome similarity between mouse and pig embryos derived in vivo and between in vivo pig embryos and pig SCNT embryos (Figure 1C). All three correlation coefficient calculations revealed similar patterns at different stages with no obvious differences between mouse and pig embryos from 4-cell stage. However, a major difference was found at the 4-cell stage between the pig SCNT-derived and in vivo fertilized-derived pig embryos. Also a prominent difference existed between TE cells at the blastocyst stage (Figure 1C and Additional file 3: Table S2).
To explore the relationships between different developmental stages, unsupervised hierarchical clustering analysis was used to evaluate similarities in global gene expression patterns. Data showed that the biggest difference between consecutive time points during mouse PED occurred from the 1- to 2-cell stage, and the second biggest difference occurred from the 8-cell to morula stages. A similar pattern was observed in normal pig embryos but the major difference occurred between the 2-cell and 4-cell stages (Figure 1D and Additional file 3: Table S2). In addition, the morula clustered together with the ICMTE (Figure 1D) in mice, while the ICMTE clustered with the TE in pigs. In the human, a different clustering pattern exists (Additional file 2: Figure S2A) . Thus pig, mouse and human PED have different patterns of gene expression.
An obvious difference was observed in the hierarchical order between pig in vivo derived embryos and pig SCNT embryos (Figure 1D; Additional file 2: Figure S2B). In SCNT embryos, the largest difference existed from the 4- to 8-cell stage and the next largest difference was noted from the morula to blastocyst stage (Figure 1D, Additional file 3: Table S2). Thus the events of reprogramming that occurred after reconstruction of SCNT embryos changed the gene expression patterns during their pre-implantation development.
Comparative analysis of the gene regulation networks important to maternal deposition and zygotic gene activation between pigs and mice
After fertilization, there is a transition from maternal to zygotic developmental control which requires both degradation of maternal RNAs and zygotic genome activation (ZGA) . We first identified maternal transcripts from oocyte transcriptomes and zygotic activated transcripts in mice and pigs based on transcription trends (Figure 2; Additional file 4: Table S3) . In vivo fertilized-derived pig embryos shared 81.3% of their maternal mRNA with in vitro SCNT-derived embryos, and 47.3% of transcripts were shared with mouse embryos (Figure 2A left). In vivo pig embryos shared 17.7% of their zygotic activation transcripts with those in mice embryos, and 50.0% were shared in SCNT-derived embryos (Figure 2A right). These results suggest that porcine maternal transcripts are more conserved than pig zygotic activated transcripts, a finding that is consistent with previous reports on humans and cattle embryos . Maternal transcript specialization may explain the low efficiency of interspecies animal cloning [23–27]. These data also suggest that transcriptional factors in the ooplasm can more effectively identify donor-cell-specific DNA domains of the same species.
To evaluate the differences in zygotic activation between mice and pigs, we counted the number of up- and down-regulated genes at different stages of early cleavage development. Major bursts of transcription occurred at the 2-cell stage in mice, at the 4-cell stage in normal pig embryos, and at the 8-cell stage in pig SCNT embryos (Figure 2B). These results may correspond to the zygotic activation stage . The timing of ZGA varies among species and occurs between 4- and 8-cell stages in humans and reportedly between 1- and 2-cell stages in mice [16, 28]. Based on our transcriptome analysis, porcine ZGA occurs at 4-cell stage, consistent with previous reports . ZGA appeared to be delayed by one cell cycle in pig SCNT embryos, compared with in vivo normal embryos (Figure 2B). It is associated with the transcripts during SCNT embryonic development that might result from artificial micromanipulation . For example, reconstructed embryos undergo only active demethylation at a slow and gradual pace, whereas normal embryos exhibit both active and passive demethylation at faster rates . We found that genes associated with demethylation, such as Dnmt3b and Dnmt1, changed their expression level after demethylation in both normal and SCNT embryos (Additional file 2: Figure S3).
Transition from maternal to zygotic transcripts is accompanied by expression of specific genes whose products associate with protein binding. Genes of known ontology (GO:0005515, protein binding) were used to construct heat-maps of hierarchical clustering (Additional file 2: Figure S4A), showing that the protein binding-associated transcripts were enriched in mouse 2-cell embryos, pig normal 4-cell stage and morula embryos, and pig 8-cell SCNT embryos. This supports the concept that ZGA timing events differ in mice and pigs (Additional file 2: Figure S4B). Also, mouse and pig embryos showed differences in transcripts associated with ATP-synthesis-coupled-proton transportation. The column dendrogram was reordered to show the ATP synthesis genes are differentially expressed in the two species. ATP-synthesis-coupled-proton transports were highly expressed at the 4-cell stage in normal pig embryos, the 8-cell stage in pig SCNT embryos, and the 2-cell stage in mouse embryos, which supports the previous findings that ZGA events occur at these stages (Additional file 2: Figure S4C). Genes and GO analysis data regarding maternal deposit and ZGA are shown in Additional file 5: Table S4.
Comparative analysis of gene regulation in pig and mouse during the three committed lineages of embryonic development
Two lineage segregation events occur in mammalian PED. In mice, the first event occurs at E3.5, when ICM and TE formation are mediated by Pou5f1 (also known as Oct4). The second event is regulated by Nanog and Gata6, resulting in the segregation of the primitive endoderm and epiblast lineages from the ICM at the blastocyst stage . To evaluate differences in the regulation of lineage segregation between pig and mouse embryos, we first examined gene expression in pig ICM cells.
Profile identification of putative pig ICM transcripts
We split the porcine blastocyst into two parts using an ultra-sharp blade. One half contained the TE only and the other comprised both ICM and TE cells (ICMTE). We then identified specific genes expressed in porcine ICM cells by comparing transcriptomes of pure TE and ICMTE in pig normal embryos. To validate this method, we identified 1,531 ICM-specific genes through differential expression analysis between ICMTE and TE mRNA sequencing data in mouse embryos (Figure 3A). 74.9% of the ICM-specific genes (Additional file 2: Figure S5A) and 71.5% of the ICMTE-specific genes (Additional file 2: Figure S5B) were included in ICM datasets that have been previously reported . This suggests that our data is reliable (Figure 3A; Additional file 5: Table S4). With this method, we identified 2,201 putative porcine ICM-specific genes in in vivo-derived embryos and 581 genes in SCNT embryos from the first sample set (Figure 3A). Similar results were obtained with the second set of normal pig embryos. With these, 2,559 human ICM-specific genes were identified . We found that 30.5% of the putative mouse ICM-specific genes and 35.3% of the human ICM-specific genes are shared with pig in vivo-derived embryos (Figure 3B and Additional file 2: Figure S5C), and that these contained ES cell markers such as Pou5f1, Tbx3 and Gata6. (Additional file 5: Table S4). GO analysis results suggested that mouse and pig ICM cells share pathways including those within the cell cycle, cell division, in utero embryonic development, the TGF-beta receptor signaling pathway, and positive regulation of transcription. However, pig ICM-specific genes were clustered with unique pathways, such as fatty acid metabolism, heat shock protein binding, and fatty acid beta-oxidation, which may be important in pig early embryonic development (Additional file 2: Figure S5B; Additional file 5: Table S4). Many ICM-specific genes in the mouse embryos were associated with mouse ES cell pluripotency. Typical mES cell markers such as Pou5f1, Sox2, Esrrb, Klf4, Mest and Tbx3 were expressed in the putative ICM, and these are believed to act within the pathways involved with stem cell maintenance and response to retinoic acid (Additional file 5: Table S4).
To determine whether the same key regulation factors act in the same signaling pathways in both pig and mouse embryos during PED, we analyzed signaling pathway patterns relevant to development or pluripotency: the TGF-beta, MAPK, Jak-Stat and Wnt pathway. Interestingly, different genes within the same pathways were activated in mouse and in pig embryos. Acvr2b, Id1 and Amhr2 were highly expressed in mouse ICM in the TGF-beta pathway, whereas Smad1, Smurf1 and Id4 were highly expressed in pig ICM. Likewise, differences were also observed for the Jak-Stat, MAPK and Wnt signaling pathways (Additional file 2: Figure S7). These data indicate that unique regulatory signaling pathways may be associated with porcine ICM development, and that these are different during lineage segregation between mouse and pig embryos.
Analysis of molecular markers for lineage commitment during PED
Next, we analyzed the expression of genes that govern early lineage segregation, specifically Oct4, Cdx2, Nanog, Sox2 and Gata6, in putative ICM and TE cells. Consistent with previous reports, we found that Oct4 expression was restricted to the ICM at the blastocyst stage during mouse PED (Additional file 2: Figure S8A). This was consistent with analysis of deep sequencing data (Figure 3C). The pattern of Oct4 expression in in vivo pig embryos was similar to that of mouse embryos (Figure 3C,D). These results suggest that Oct4, which regulates early development in mice, may have a similar function for lineage segregation in pigs. In mice, the expression of Cdx2 increased gradually from the 8-cell stage, and was predominant in TE at the blastocyst stage (Figure 3C). We observed a similar tendency of Cdx2 expression during pig PED, but the expression level was much lower than in mice (Figure 3C; Additional file 2: Figure S8B). Nanog expression was very low and could not be detected in the morula or blastocyst by immunostaining during pig PED (Figure 3C; Additional file 2: Figure S8C). In addition, the pattern of Gata6 and Sox2 expression in pig PED was different to that observed in mice (Figure 3C; Additional file 2: Figure S8D,E). This suggests that the regulation of second lineage segregation events during PED might differ between mice and pigs.
Comparison of gene regulation in lineage segregation during pig and mouse PED
To understand the mechanisms underlying lineage segregation in pig PED, we analyzed co-expressed genes in the morula and ICM and genes shared by the morula and TE in mouse and pig (Figure 4A). Both gene groups may be relevant to ICM development and TE formation independently because both ICM and TE differentiate from the morula.
Among the co-expressed genes in the morula and ICM, we observed that 295 genes (27.7% in ms; 17.3% in Pnm) were shared between Pnm and Ms, including Oct4, Tbx3, Gata6 and Smad4. These genes may be involved in regulating second lineage segregation (Figure 4A left; Additional file 6: Table S5). Using GO analysis of the shared genes co-expressed by the morula and ICM between the mouse and pig, we found that most gene expression is associated with embryonic biology, including in utero embryonic and endoderm development. Mouse and pig embryos also differed in co-expressed genes related to fatty acid metabolic processes, lipid metabolic processes, and biological aspects of the cytoplasm, nucleus, mitochondria, and protein binding (Additional file 6: Table S5).
A total of 1,217 genes were commonly expressed in the morula and TE cells in Ms and in Pnm (53.6% in mice and 54.8% in pigs) (Figure 4A right). Genes related to lineage commitment, such as Tcf3, Tcfap2c, Cnot7, Grb2 and Smad2, were found in both pigs and mice. Mouse morula and TE co-expressed genes included Cdx2, Fgfr1op, Tcf15, Trim11, Trim 27, Eed, Gata3, Fgfr1, Cnot8, Grb7, and Tcf25. Pnm morula and TE co-expressed genes included Trim15, Trim37, Trim35, Gata2, Fgf7 and Grb10 but not Cdx2. (Additional file 6: Table S5). These data indicate that TE lineage differentiation is regulated by different signaling pathways or by different levels of gene expression in pig and mouse embryos. Comparative analysis of the expression of specific markers of ICM and TE in the different stages of embryos was performed (Figure 4B and Additional file 2: Figure S6). A significant difference in the heat map patterns for these marker genes suggests that regulation of ICM development among pigs, mice and the human were more conserved than those of TE development.
Expression of enzymes relevant to fatty acid metabolism in pig early embryonic development
Comparing gene expression in mouse and pig embryos and their relevant pathways during PED, we observed an enrichment in the pathways relevant to fatty acid metabolism during ZGA and lineage segregation in pig embryos. We then evaluated the pathways of fatty acid biosynthesis, fatty acid elongation, and fatty acid metabolism during PED. Genes for enzymes that regulate fatty acid biosynthesis and elongation were highly expressed at the 4-cell and the morula stage in normal pig embryos, compared with mouse and pig SCNT embryos. In contrast, a slight increase in expression of these enzymes was observed in mouse blastocysts and in pig 8-cell SCNT embryos (Figure 5; Additional file 7: Table S6). Thus, fatty acid biosynthesis and elongation in pre-implantation embryos may be important to pig embryonic and fetal development. In addition, enzymes related to fatty acids are activated during different stages of development and at different levels across species. Under current culture conditions, oocytes for pig SCNT embryos were allowed to mature in vitro, and the reconstructed embryos were cultured in fatty-acid-free medium. This may explain the low efficiency of pig cloning. Further studies are warranted to understand the mechanisms by which fatty acid metabolism is regulated during pig PED and in in vitro culture of pig pre-implantation embryos, ES cells and iPSCs.
Nuclear transfer micromanipulation and global gene expression patterns during reconstructed embryo development
The ooplasm contains factors that erase somatic epigenetic imprints, rendering the somatic nucleus totipotent. Gene expression patterns that occur after nuclear transfer may be essential to later development. Identification and comparison of the unique gene expression patterns of normal and cloned pig embryos may help determine the unique requirements of cloned embryos and inform us about specific factors in pig iPSCs reprogramming.
To explore mechanisms underlying reprogramming of SCNT, we analyzed serial gene expression of normal and cloned embryos at different development stages. The correlation coefficient between normal ICMTE and SCNT embryonic ICMTE was found to be 0.93 (Additional file 3: Table S2). We observed fewer putative ICM-specific genes in the cloned embryos than in the normal embryos (Figure 3B). Perhaps cloned embryos contain fewer ICM cells than normal embryos [34, 35]. GO analysis of genes expressed in cloned embryonic putative ICM showed that they participated in pathways involving protein binding and regulation in the cytoplasm, mitochondria and nucleus, presumably because of interactions between the ooplasm and the transferred somatic nucleus. However, genes expressed in normal embryos were related to signaling pathways involving the cytoplasm, nucleus, mitochondria, ATP binding, nucleotide binding and fatty acid metabolic processes (Additional file 2: Figure S5B; Additional file 5: Table S4). Fatty acid metabolism appears to be unique to pig early embryonic development, and its function is unclear. In addition, pig normal and SCNT embryos differed in signaling pathways relevant to development including FGF, MAPK, Jak-Stat, and Wnt pathways (Additional file 2: Figure S9).
In mouse SCNT, functional defects in the trophoblast cell lineage are the main cause of low cloning success rates . Pnm and Pnt comparative analysis of the transcriptomes of embryos at different stages revealed that the correlation of TE from blastocysts is lower than that of ICMTE (Figure 1C). GO analysis revealed differences in biological processes in the nucleus, cytoplasm, and cell organs between normal and cloned embryonic TE (Additional file 8: Table S7). This suggests that the functional defects in the trophoblastic cell lineage might also contribute to low efficiency in pig cloning.
Putative gene regulation networks during pig PED
From our data, we assumed that most of the two lineage segregations during PED were conserved across the three-layer networks of mouse and pig embryos collected in vivo. Marker genes (§ Experimental procedures) for the three layers were selected based on previous findings [18, 37]. Co-expression networks of those marker genes were constructed based on the Pearson correlation coefficient (Additional file 2: Figures S10A and B), and singleton sub-networks (common targeted genes from two distinct marker genes that are less than 5% we defined as a singleton network) were removed from the final network.
Network evaluation revealed that pig and mouse pre-implantation embryos share some hub genes (about 50%) with regulatory networks for lineage segregations. In sub-network analysis, ICM-lineage segregation and primitive endoderm differentiation markers such as Esrrb, Gata6, Pou5f1, Tcl1, Sall4 and Tbx3, and TE commitment makers such as Hand1, Ccdc3, Ets2, Wnt7a and Eomes, were found to be common to pigs and mice (Additional file 2: Figures S10A and B). KEGG pathway analysis of all of the networks involved with these genes revealed that the regulation process in pig embryos is mostly associated metabolic pathways whereas that in mouse embryos, is mostly associated with signal transduction (Additional file 2: Figure S10C and Additional file 9: Table S8). These results also suggest that metabolic pathways, including fatty acid metabolism, are important in lineage segregation during pig PED, although details of the underlying mechanism remain to be further explored.
We developed a platform for mRNA sequencing of porcine pre-implantation embryos and a method for evaluating putative ICM gene expression by comparing the ICMTE and TE transcriptome. Genome-wide transcriptome analysis revealed unique molecular gene regulation networks which regulate ZGA, lineage segregation and embryonic metabolism during pig PED.
Some highly expressed genes in pig putative ICM, most of which are related to ES cell pluripotency, differed from those highly expressed in mouse ICM . Even with shared pathways that occur at the same stage of PED, the specific members of special signaling pathways that are highly expressed were different between mice and pigs. This indicates that pig PED is regulated by a species-specific network of genetic factors. Identification of the activated genes in ES cell-specific pathways such as the TGF-beta, MAPK, Jak-Stat and Wnt pathways support this conclusion. SCNT is a reliable platform for studying somatic cell reprogramming. Differences between the Pnm and Pnt embryos, especially between their ICMs, provide insights into the mechanisms involved in somatic reprogramming, including those involved in the generation of iPSCs. For example, we found that fatty acid metabolic processes to be significantly more active in Pnm than in mouse or SCNT pig pre-implantation embryos. These results indicate that fatty acid metabolism may be important for pig PED and lineage segregation events. It may also be important for the maintenance of pluripotency, which may be a useful nutritional element during derivation of true pig ES cell lines.
The molecular regulators of early lineage segregation in mouse embryos are well known but those in pig embryos are not. Transcriptome analysis of genes co-expressed by the morula and ICM, and co-expressed by the morula and TE in the pig and comparison with genes common to those in mice revealed that some markers of early embryonic segregation, such as Cdx2 and Nanog, may not be restricted to the ICM or the TE in the pig blastocyst. This indicates that either a different mechanism regulates the differentiation of the pig pre-implantation embryonic ectoderm lineage or these events occur at a later stage in pig PED. Analysis of genes co-expressed by the morula and ICM together with morula and TE revealed that regulation of ICM differentiation is more conserved. However, the regulatory mechanism of TE differentiation may be subtly different in mouse and pig embryos because Pnm ICM expresses Eed and Gata3, which are restricted to the TE in mice. It is unclear whether Oct4 expression may limit Cdx2 expression in the pig TE as it does in mouse embryos. It is also unclear if Oct4 expression in pig TE causes a decline in the control of TE differentiation or if there are other factors that are expressed in the Pnm TE that regulate TE differentiation.
The events that occur shortly after fertilization in normal embryos include segregation of maternal chromosomes, breakdown of the sperm’s nuclear envelope, repackaging of the sperm’s chromatin, and the formation of two pronuclei. In cloned embryos, the sub-cellular changes include the breakdown of the somatic cell nuclear envelope, the condensation of the somatic cell chromosomes, and the formation of pseudo-pronucleus . The different processes that occur during this transient period may cause gene expression patterns to differ between normal and cloned embryos, which may then affect the regulation of later development. Aberrant global gene expression and delay of the ZGA in SCNT pre-implantation embryos may be the reasons for the incomplete reprogramming observed in SCNT embryos.
Evidence from previous studies indicates that individual fatty acids may affect oocyte maturation and embryonic development [38, 39]. Endogenous lipids are more abundant in pig oocytes and pre-implantation embryos than in mouse, cattle and sheep. This can cause significant damage during cryopreservation . Lipids have important functions in membrane receptor biology, signal transduction, and growth regulation [41–46]. However, little is known about the fatty acid metabolism network and its regulation in pre-implantation embryos. Information concerning lipid uptake and utilization by oocytes and embryos may therefore be crucial to improving cryopreservation practices, in vitro culture systems for the derivation of pig ESCs, and the combination of the factors used for reprogramming into pig iPSCs.
In summary, our study represents a significant step in the characterization of pig pre-implantation embryos and provides insight into the dynamic molecular regulation of pig PED. We provide genome-level evidence and gene expression patterns for events such as ZGA and lineage segregation during PED. GO and KEGG analyses of each trancriptome suggests that pathways which regulate epiblast and primitive endoderm lineage commitment may be more conserved than those that regulate ectoderm segregation. Our data also provide a resource for pig pluripotent stem cell engineering and for understanding porcine development.
Embryo collection and RNA isolation
Young adult female C57 mice (Vital River Laboratories, Beijing, China) and young adult female Nong Da Xiang pigs (China Agricultural University pig farm, Zhuo Zhou, China) were kept in a 12 h light/dark cycle and given water and food ad libitum. All mice and pigs were handled and studies were carried out according to the guidelines of The State Key Laboratory Animal Care and Use Committee.
Oocytes were collected from the oviduct of C57 mice 14 h after injection with human chorionic gonadotropin (hCG). Then pre-implantation embryos were collected at various points after hCG injection and mating as depicted in Table 1. Pure TE was separated physically using an ultra-sharp splitting blade (Bioniche, Animal Health US, Inc) under a stereomicroscope (Additional file 2: Figure S1).
Normal pig embryos were washed from the oviduct or uterus using PBS with 5% FBS at the indicated time points (Table 1) after estrus and natural mating. Donors of normal embryos were mini pigs (Nong Da Xiang, a local strain), and the same strain was used to harvest fibroblasts that were used as donor cells for SCNT. SCNT embryos were collected after activation of reconstructed embryos at time points given in Additional file 1: Table S1.
Next, 4–10 embryos at the same stage were pooled together for each sample and transferred into extraction buffer from the PicoPure RNA Isolation Kit (Arcturus, KIT0204, Life Technologies, US) at 42°C for 30 min. Samples were either stored at −80°C for up to one month or used immediately for analysis. RNA from each sample was extracted according to the manufacturer’s instructions and eluted into 10-μL elution buffer.
We performed mRNA sequencing using the Applied Biosystems SOLiD 4 System as follows: RNAs isolated from pre-implantation embryos were used for double-stranded cDNA synthesis and PCR amplification. The procedure for fragment library preparations was performed according to the Library Preparation Protocol for whole transcriptome analysis of a single cell (http://www.appliedbiosystems.com) and the Applied Biosystems SOLiD 4 System Library Preparation Guide (http://www.appliedbiosystems.com). mRNA sequencing data are available from the Lab Archive (http://www.ncbi.nlm.nih.gov/sra) under accession number SRA076823.
Transcriptome analysis tools included in the BioScope software package (Applied Biosystems) were used to map the corresponding sequenced reads against the mouse mm9 genome and the pig susScr2 genome. Mouse (NCBIM37.65) and pig (Sscrofa9.65) gene and transcript annotation files were downloaded from the Ensembl database (http://asia.ensembl.org/info/data/ftp/index.html).
Prior to differential gene expression analysis, read count tables were generated from binary sequence alignment/map (BAM) files using HTseq software (http://pypi.python.org/pypi/HTSeq). The value for reads per kilobase of coding sequence per million mapped reads (RPKM) was calculated to estimate gene expression under each set of conditions . For each sequenced library, the read counts were adjusted using the edgeR software package through a one-scaling normalized factor [48, 49]. The DEGseq software package was used to calculate differences in gene expression between the two assigned groups. A gene was considered significant if the Benjamini and Hochberg corrected P-value was less than 0.05 and the fold-change was greater than 2 . Differentially expressed isoforms were estimated using Cufflinks v1.30 while treating early pre-implantation embryo samples as a time-series input.
We also downloaded human pre-implantation embryonic datasets from GSE36552, deposited in GEO databases (http://www.ncbi.nlm.nih.gov/geo/) . For each developmental embryonic stage, pooled read counts were used for recalculating RPKM values. EPI, PE and TE single cells datasets in blastocysts were analyzed as described in this paper. All genes with RPKM > 0.1 were defined as expressed in individual cells at same stage. Based on our results, we then compared pig and mouse in vivo datasets with human PED datasets.
Gene ontology (GO) was implemented using the GOseq software package, in which gene length bias was adjusted . The BioMart system was used to abstract orthologous relationships between pig and mouse (one-to-one) gene pairs and to convert mapping information regarding gene identities between mouse Ensembl and reference sequences [52, 53].
ICM segregation markers included Pou5f1, Sox2, Klf2, Nanog, Rex1, Utf1, Zfx, Esrrb, Tbx3, Tcl1 and Klf4. Primitive endoderm differentiation markers included Gata6, Gata4, Gdf1, Gdf3, Hnf4a (Nr2a1), Mixl1, Sall4, Sox7 and Sox17. TE commitment markers included Cdx2, Eomes, Hand1, Fgfr2, Ets2, Tcfap2, Elf5, Etv4, Furin, Ccdc3, Pace4, Casq1, Wnt7a, Fgf5, Pax6, and Tead4 [18, 37]. Only genes with orthologous relationships between pigs and mice were kept for analysis (Klf2, Rex1, Gata4, Mixl1, Sox7, and Pace4 were removed). Three types of markers were chosen for co-expression sub-network construction analysis, which were displayed graphically using Cytoscape 2.8 . The absolute value of the correlation between expression profiles includes all points in time during embryonic development. The correlation between markers and predicted targets must be larger than 0.9.
Embryos were fixed in 4% paraformaldehyde in washing solution (PBS with 0.1% Tween 20 and 0.01% Triton X-100) for 30 min at room temperature. Then embryos were permeabilized in 1% Triton X-100 in PBS for 4 h at 4°C, blocked with 1% BSA in washing solution (blocking solution), and incubated with the following antibodies: Oct4 (Santacruz Sc-8628), Cdx2 (Biogenex MU392A-UC), Nanog (Abnova PAB6837), Sox2 (Santacruz Sc-17320) and Gata6 (Abcam Ab22600) in blocking solution for 1 h at room temperature. After incubation with secondary antibody for 1 h at room temperature, embryos were counterstained with 10 μg/mL Hoechst 3342 (Sigma B2261) in washing solution for 10 min at room temperature.
Hubner K, Fuhrmann G, Christenson LK, Kehler J, Reinbold R, De La Fuente R, Wood J, Strauss JF, Boiani M, Scholer HR: Derivation of oocytes from mouse embryonic stem cells. Science. 2003, 300 (5623): 1251-1256. 10.1126/science.1083452.
Geijsen N, Horoschak M, Kim K, Gribnau J, Eggan K, Daley GQ: Derivation of embryonic germ cells and male gametes from embryonic stem cells. Nature. 2004, 427 (6970): 148-154. 10.1038/nature02247.
Buehr M, Meek S, Blair K, Yang J, Ure J, Silva J, McLay R, Hall J, Ying QL, Smith A: Capture of authentic embryonic stem cells from rat blastocysts. Cell. 2008, 135 (7): 1287-1298. 10.1016/j.cell.2008.12.007.
Nagy A, Gocza E, Diaz EM, Prideaux VR, Ivanyi E, Markkula M, Rossant J: Embryonic stem cells alone are able to support fetal development in the mouse. Development. 1990, 110 (3): 815-821.
Evans MJ, Kaufman MH: Establishment in culture of pluripotential cells from mouse embryos. Nature. 1981, 292 (5819): 154-156. 10.1038/292154a0.
Thomson JA, Itskovitz-Eldor J, Shapiro SS, Waknitz MA, Swiergiel JJ, Marshall VS, Jones JM: Embryonic stem cell lines derived from human blastocysts. Science. 1998, 282 (5391): 1145-1147.
Hall V: Porcine embryonic stem cells: a possible source for cell replacement therapy. Stem Cell Rev. 2008, 4 (4): 275-282. 10.1007/s12015-008-9040-2.
Brevini TA, Antonini S, Cillo F, Crestan M, Gandolfi F: Porcine embryonic stem cells: facts, challenges and hopes. Theriogenology. 2007, 68 (Suppl 1): S206-S213.
Ezashi T, Telugu BP, Alexenko AP, Sachdev S, Sinha S, Roberts RM: Derivation of induced pluripotent stem cells from pig somatic cells. Proc Natl Acad Sci USA. 2009, 106 (27): 10993-10998. 10.1073/pnas.0905284106.
Wu Z, Chen J, Ren J, Bao L, Liao J, Cui C, Rao L, Li H, Gu Y, Dai H, et al: Generation of pig induced pluripotent stem cells with a drug-inducible system. J Mol Cell Biol. 2009, 1 (1): 46-54. 10.1093/jmcb/mjp003.
Esteban MA, Xu J, Yang J, Peng M, Qin D, Li W, Jiang Z, Chen J, Deng K, Zhong M, et al: Generation of induced pluripotent stem cell lines from Tibetan miniature pig. J Biol Chem. 2009, 284 (26): 17634-17640. 10.1074/jbc.M109.008938.
West FD, Uhl EW, Liu Y, Stowe H, Lu Y, Yu P, Gallegos-Cardenas A, Pratt SL, Stice SL: Brief report: chimeric pigs produced from induced pluripotent stem cells demonstrate germline transmission and no evidence of tumor formation in young pigs. Stem Cells. 2011, 29 (10): 1640-1643. 10.1002/stem.713.
Hall VJ, Christensen J, Gao Y, Schmidt MH, Hyttel P: Porcine pluripotency cell signaling develops from the inner cell mass to the epiblast during early development. Dev Dyn. 2009, 238 (8): 2014-2024. 10.1002/dvdy.22027.
Xie D, Chen CC, Ptaszek LM, Xiao S, Cao X, Fang F, Ng HH, Lewin HA, Cowan C, Zhong S: Rewirable gene regulatory networks in the preimplantation embryonic development of three mammalian species. Genome Res. 2010, 20 (6): 804-815. 10.1101/gr.100594.109.
Duncan FE, Schultz RM: Gene expression profiling of mouse oocytes and preimplantation embryos. Methods Enzymol. 2010, 477: 457-480.
Hamatani T, Carter MG, Sharov AA, Ko MS: Dynamics of global gene expression changes during mouse preimplantation development. Dev Cell. 2004, 6 (1): 117-131. 10.1016/S1534-5807(03)00373-3.
Sharov AA, Piao Y, Matoba R, Dudekula DB, Qian Y, VanBuren V, Falco G, Martin PR, Stagg CA, Bassey UC, et al: Transcriptome analysis of mouse stem cells and early embryos. PLoS Biol. 2003, 1 (3): E74-
Tang F, Barbacioru C, Bao S, Lee C, Nordman E, Wang X, Lao K, Surani MA: Tracing the derivation of embryonic stem cells from the inner cell mass by single-cell RNA-Seq analysis. Cell Stem Cell. 2010, 6 (5): 468-478. 10.1016/j.stem.2010.03.015.
Fang H, Yang Y, Li C, Fu S, Yang Z, Jin G, Wang K, Zhang J, Jin Y: Transcriptome analysis of early organogenesis in human embryos. Dev Cell. 2010, 19 (1): 174-184. 10.1016/j.devcel.2010.06.014.
Tang F, Barbacioru C, Nordman E, Li B, Xu N, Bashkirov VI, Lao K, Surani MA: RNA-Seq analysis to capture the transcriptome landscape of a single cell. Nat Protoc. 2010, 5 (3): 516-535. 10.1038/nprot.2009.236.
Yan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J, et al: Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol. 2013, 20 (9): 1131-1139. 10.1038/nsmb.2660.
Schier AF: The maternal-zygotic transition: death and birth of RNAs. Science. 2007, 316 (5823): 406-407. 10.1126/science.1140693.
Wen DC, Yang CX, Cheng Y, Li JS, Liu ZH, Sun QY, Zhang JX, Lei L, Wu YQ, Kou ZH, et al: Comparison of developmental capacity for intra- and interspecies cloned cat (Felis catus) embryos. Mol Reprod Dev. 2003, 66 (1): 38-45. 10.1002/mrd.10333.
Chung Y, Bishop CE, Treff NR, Walker SJ, Sandler VM, Becker S, Klimanskaya I, Wun WS, Dunn R, Hall RM, et al: Reprogramming of human somatic cells using human and animal oocytes. Cloning Stem Cells. 2009, 11 (2): 213-223. 10.1089/clo.2009.0004.
Lorthongpanich C, Laowtammathron C, Chan AW, Ketudat-Cairns M, Parnpai R: Development of interspecies cloned monkey embryos reconstructed with bovine enucleated oocytes. J Reprod Dev. 2008, 54 (5): 306-313. 10.1262/jrd.20049.
Beyhan Z, Iager AE, Cibelli JB: Interspecies nuclear transfer: implications for embryonic stem cell biology. Cell Stem Cell. 2007, 1 (5): 502-512. 10.1016/j.stem.2007.10.009.
Chen T, Zhang YL, Jiang Y, Liu JH, Schatten H, Chen DY, Sun QY: Interspecies nuclear transfer reveals that demethylation of specific repetitive sequences is determined by recipient ooplasm but not by donor intrinsic property in cloned embryos. Mol Reprod Dev. 2006, 73 (3): 313-317. 10.1002/mrd.20421.
Braude P, Bolton V, Moore S: Human gene expression first occurs between the four- and eight-cell stages of preimplantation development. Nature. 1988, 332 (6163): 459-461. 10.1038/332459a0.
Jarrell VL, Day BN, Prather RS: The transition from maternal to zygotic control of development occurs during the 4-cell stage in the domestic pig, Sus scrofa: quantitative and qualitative aspects of protein synthesis. Biol Reprod. 1991, 44 (1): 62-68. 10.1095/biolreprod44.1.62.
Latham KE: Early and delayed aspects of nuclear reprogramming during cloning. Biol Cell. 2005, 97 (2): 119-132. 10.1042/BC20040068.
Mann MR, Bartolomei MS: Epigenetic reprogramming in the mammalian embryo: struggle of the clones. Genome Biol. 2002, 3 (2): REVIEWS1003-
Arnold SJ, Robertson EJ: Making a commitment: cell lineage allocation and axis patterning in the early mouse embryo. Nat Rev Mol Cell Biol. 2009, 10 (2): 91-103. 10.1038/nrm2618.
Zhang JD, Wiemann S: KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor. Bioinformatics. 2009, 25 (11): 1470-1471. 10.1093/bioinformatics/btp167.
Terashita Y, Sugimura S, Kudo Y, Amano R, Hiradate Y, Sato E: Improving the quality of miniature pig somatic cell nuclear transfer blastocysts: aggregation of SCNT embryos at the four-cell stage. Reprod Domest Anim. 2011, 46 (2): 189-196. 10.1111/j.1439-0531.2010.01614.x.
Niemann H, Rath D: Progress in reproductive biotechnology in swine. Theriogenology. 2001, 56 (8): 1291-1304. 10.1016/S0093-691X(01)00630-6.
Lin J, Shi L, Zhang M, Yang H, Qin Y, Zhang J, Gong D, Zhang X, Li D, Li J: Defects in trophoblast cell lineage account for the impaired in vivo development of cloned embryos generated by somatic nuclear transfer. Cell Stem Cell. 2011, 8 (4): 371-375. 10.1016/j.stem.2011.02.007.
Loh YH, Wu Q, Chew JL, Vega VB, Zhang W, Chen X, Bourque G, George J, Leong B, Liu J, et al: The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nat Genet. 2006, 38 (4): 431-440. 10.1038/ng1760.
Waterman RA, Wall RJ: Lipid interactions with in vitro development of mammalian zygotes. Gamete Res. 1988, 21 (3): 243-254. 10.1002/mrd.1120210306.
Nonogaki T, Noda Y, Goto Y, Kishi J, Mori T: Developmental blockage of mouse embryos caused by fatty acids. J Assist Reprod Genet. 1994, 11 (9): 482-488. 10.1007/BF02215713.
McEvoy TG, Coull GD, Broadbent PJ, Hutchinson JS, Speake BK: Fatty acid composition of lipids in immature cattle, pig and sheep oocytes with intact zona pellucida. J Reprod Fertil. 2000, 118 (1): 163-170.
Loh HH, Law PY: The role of membrane lipids in receptor mechanisms. Annu Rev Pharmacol Toxicol. 1980, 20: 201-234. 10.1146/annurev.pa.20.040180.001221.
Watkins EB, Miller CE, Majewski J, Kuhl TL: Membrane texture induced by specific protein binding and receptor clustering: active roles for lipids in cellular function. Proc Natl Acad Sci USA. 2011, 108 (17): 6975-6980. 10.1073/pnas.1014579108.
Harder T: Lipid raft domains and protein networks in T-cell receptor signal transduction. Curr Opin Immunol. 2004, 16 (3): 353-359. 10.1016/j.coi.2004.03.013.
Bader B, Kuhn K, Owen DJ, Waldmann H, Wittinghofer A, Kuhlmann J: Bioorganic synthesis of lipid-modified proteins for the study of signal transduction. Nature. 2000, 403 (6766): 223-226. 10.1038/35003249.
Rondeau I, Picard S, Bah TM, Roy L, Godbout R, Rousseau G: Effects of different dietary omega-6/3 polyunsaturated fatty acids ratios on infarct size and the limbic system after myocardial infarction. Can J Physiol Pharmacol. 2011, 89 (3): 169-176. 10.1139/Y11-007.
Sellmayer A, Danesch U, Weber PC: Effects of different polyunsaturated fatty acids on growth-related early gene expression and cell growth. Lipids. 1996, 31 (Suppl): S37-S40.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.
Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010, 26 (1): 139-140. 10.1093/bioinformatics/btp616.
Robinson MD, Oshlack A: A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11 (3): R25-10.1186/gb-2010-11-3-r25.
Wang L, Feng Z, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138. 10.1093/bioinformatics/btp612.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515. 10.1038/nbt.1621.
Vilella AJ, Severin J, Ureta-Vidal A, Heng L, Durbin R, Birney E: EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 2009, 19 (2): 327-335.
Smedley D, Haider S, Ballester B, Holland R, London D, Thorisson G, Kasprzyk A: BioMart biological queries made easy. BMC Genomics. 2009, 10 (1): 22-10.1186/1471-2164-10-22.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.
This study was supported by grants from the China National Basic Research Program (2011CBA01001, 2009CB941003, 2011CBA01102, 2010CB945404), National Thousand Talents Program of China, the Program for New Century Excellent Talents in University (NCET-11-0482), 111 Project (B12008) and the project from State Key Laboratories for Agrobiotechnology (201103). We thank Professor Lin Liu from Nankai University, Professor Bruce Whitelaw from The Roslin Institute, United Kingdom, and Dr. See Ting Leong from the Genome Institute of Singapore for their assistance and suggestions.
The authors declare that they have no competing interests.
JH and NL designed the project. SC, XW and XR carried out the library preparation and the RNA sequencing. JW, SC and JH carried out the bioinformatics analysis. SL, ZL and QL performed embryonic immunostaining experiments. Both JH and BL drew conclusions and JH drafted the manuscript. All authors have read and approved the final manuscript.
Suying Cao, Jianyong Han, Jun Wu contributed equally to this work.
Electronic supplementary material
Additional file 2:Supplement figures. Figure S1. Splitting of pig blastocyst with an ultra-sharp splitting blade. Figure S2. (A) Unsupervised hierarchical clustering of the expression profiles . Figure S3. Dnmt3b and Dnmt1 expression. Figure S4. Heat map of protein binding-associated transcripts under different conditions. Figure S5. Analysis of the ICM-specific genes in mouse and pig. Figure S6. Heat map of ICM and TE marker gene clustering in human pre-implantation embryos. Figure S7. Heat map clusters of ICM-specific genes in TGF-beta, MAPK, Jak-Stat and Wnt signaling pathways expressed during normal mouse PED and normal pig PED. Figure S8. Immunostaining of the pluripotent markers in pig and mouse pre-implantation embryos. Figure S9. Heat map clusters of ICM-specific genes in TGF-beta, MAPK, Jak-Stat and Wnt signaling pathways expressed during normal pig in vivo and SCNT pig PED. Figure S10. Putative regulatory networks for lineage segregation during PED. (PDF 1 MB)
About this article
Cite this article
Cao, S., Han, J., Wu, J. et al. Specific gene-regulation networks during the pre-implantation development of the pig embryo as revealed by deep sequencing. BMC Genomics 15, 4 (2014). https://doi.org/10.1186/1471-2164-15-4
- Fatty Acid Metabolism
- Inner Cell Mass
- Oct4 Expression
- Zygotic Genome Activation
- Reconstructed Embryo