Transcriptional reprogramming of gene expression in bovine somatic cell chromatin transfer embryos
© Rodriguez-Osorio et al. 2009
Received: 12 January 2008
Accepted: 24 April 2009
Published: 24 April 2009
Skip to main content
© Rodriguez-Osorio et al. 2009
Received: 12 January 2008
Accepted: 24 April 2009
Published: 24 April 2009
Successful reprogramming of a somatic genome to produce a healthy clone by somatic cells nuclear transfer (SCNT) is a rare event and the mechanisms involved in this process are poorly defined. When serial or successive rounds of cloning are performed, blastocyst and full term development rates decline even further with the increasing rounds of cloning. Identifying the "cumulative errors" could reveal the epigenetic reprogramming blocks in animal cloning.
Bovine clones from up to four generations of successive cloning were produced by chromatin transfer (CT). Using Affymetrix bovine microarrays we determined that the transcriptomes of blastocysts derived from the first and the fourth rounds of cloning (CT1 and CT4 respectively) have undergone an extensive reprogramming and were more similar to blastocysts derived from in vitro fertilization (IVF) than to the donor cells used for the first and the fourth rounds of chromatin transfer (DC1 and DC4 respectively). However a set of transcripts in the cloned embryos showed a misregulated pattern when compared to IVF embryos. Among the genes consistently upregulated in both CT groups compared to the IVF embryos were genes involved in regulation of cytoskeleton and cell shape. Among the genes consistently upregulated in IVF embryos compared to both CT groups were genes involved in chromatin remodelling and stress coping.
The present study provides a data set that could contribute in our understanding of epigenetic errors in somatic cell chromatin transfer. Identifying "cumulative errors" after serial cloning could reveal some of the epigenetic reprogramming blocks shedding light on the reprogramming process, important for both basic and applied research.
The process of early embryonic development is determined by activation of the embryonic genome, which for bovine embryos begins as a "minor genome activation" at the 1-cell stage  ascending to a "major genome activation" during the 8-cell to 16-cell stage . In the absence of proper genome activation, the developing embryo will die because it can no longer support its essential developmental functions [3, 4]. In the case of embryos produced by somatic cell nuclear transfer (SCNT) the somatic nucleus has to be reprogrammed in order to restart and continue the developmental process. It is believed that, guided by the ooplasm, the somatic nucleus aborts its own program of somatic gene expression and re-establishes a particular program of embryonic gene expression necessary for normal embryo development .
Embryos produced by SCNT have lower developmental rates than their in vitro and in vivo produced counterparts . Embryos produced by SCNT also have a greater incidence of apoptosis and consequently a lower number of cells . Additionally, SCNT derived embryos have greater rates of embryo and fetal mortality, stillbirths and perinatal deaths, which bring down the overall efficiency of cloning. These alterations may be caused, at least partially, by incomplete epigenetic reprogramming of the somatic nuclei [5, 7]. Somatic cell chromatin transfer (SCCT) attempts to facilitate the reprogramming process by exposing the somatic cells, prior to the transfer, to a mitotic cell extract, which is supposed to induce chromosome condensation and promote the removal and solubilisation of nuclear factors, enhancing nuclear remodelling . Compared with nuclear transfer, SCCT shows greater survival of cloned calves up to at least 1 month and could be a useful tool in understanding the mechanisms of reprogramming . Remarkably, a recent study did not detect any significant differences in the global gene expression profiles of SCCT and SCNT embryos .
Embryos derived from nuclear transfer have an abnormal pattern of DNA methylation, in some cases resembling that of somatic cells [10–12]. This aberrant DNA methylation pattern has been inversely correlated with the developmental potential of the cloned embryos . Treatment of donor cells with DNA demethylation agents, prior the nuclear transfer, may remove epigenetic marks improving the ability of the somatic cells to be fully reprogrammed by the recipient karyoplast . Global alteration of gene expression has been another finding in embryos produced by cloning. The abnormal expression of genes playing important roles in early embryonic development, implantation and fetal development is of particular interest. Conversely, other studies have reported a significant reprogramming for SCNT embryos by the blastocyst stage and similar transcriptome profiles to those of embryos produced in vitro or in vivo, suggesting that defects in gene expression for SCNT embryos may occur later during redifferentiation and organogenesis [15, 16].
Among the abnormally expressed genes reported in bovine cloned embryos are IL6, FGF4, and FGFr2 ; FGF4, DNMT1, Mash2, HSP70, and interferon tau ; Acrogranin, Cdx2, and ERR2 . Cytokeratin 19, Cytokeratin 8, Vimentin, Hsp27, Nidogen2 and MHC-I ; HDAC-1, 2, and 3, DNMT3A, and OCT4 . Lower levels of transcripts involved in the retinoic acid signalling pathway (RARB, CRAB1, HLA-A, THBS2, and SERPINB5) were reported for cloned bovine embryos . There have been conflicting results when it comes to the expression of particular genes in SCNT and IVF embryos. Such is the case of the developmentally important POU5F1 gene, which has been reported as misregulated in cloned embryos compared to IVF derived blastocysts in some studies [21, 23], while being detected at similar concentration in others [17, 24].
SCNT is often used for the production of human proteins in the milk of transgenic animals. For the achievement of some specific transgenic phenotypes, multiple genetic modifications need to be completed through sequential modifications in primary cells prior to nuclear transfer . Since transfection and selection of transgenic cells requires nearly the entire lifespan of a cell, only one genetic modification can be completed in each cell lifespan . Therefore, consecutive rounds of cloning (also referred to "repeated cloning", "serial cloning", "recloning" or "nuclear recycling") are performed. It has been proposed that consecutive rounds of cloning, allow for rejuvenation and selection of transformed cultured cells [27–30] and that it may improve the efficiency of SCNT by increasing the reprogramming potential of the somatic cells [31, 32]. Conversely, other reports suggest that epigenetic errors could accumulate in the embryos as a result of serial cloning and prolonged in vitro culture decreasing cloning efficiency. After serial cloning up to the sixth generation was performed in mice, no signs of telomere shortening or premature ageing were observed. However, cloning efficiency significantly decreased with increasing rounds of cloning . A greatly reduced in vitro and in vivo developmental capaCity was reported for bovine embryos derived after several rounds of serial cloning [34, 35]. It has been suggested that extended culture associated with transfection and selection procedures may induce changes of somatic cells, which decrease the efficiency of nuclear transfer and that these changes cannot be reversed by recloning .
The objective of the present study was to identify the "cumulative errors" on global gene expression, caused by serial rounds of chromatin transfer, by comparing the transcriptome profile of IVF derived blastocysts to that of SCCT derived blastocysts from the first and fourth rounds of cloning (CT1 and CT4) using oligonucleotide microarray analysis (Affymetrix Bovine GeneChips). Donor cells used for first and fourth rounds of cloning (DC1 and DC4) were also the target of the study as we compared the global gene expression of the SCCT embryos with their respective donor cells. Additionally, we analyzed the expression patterns of a panel of selected genes, in fetal fibroblasts obtained from foetuses from zero to fifth rounds of chromatin transfer. Our results show that a substantial reprogramming has taken place in the cloned embryos from both generations of chromatin transfer. However, there was a set of differential expressed genes in both groups of cloned embryos compared to their IVF counterparts. The number and functions of these genes could suggest accumulative misregulations probably caused by the successive rounds of cloning.
The Affymetrix GeneChip® Bovine Genome Array contains 24,129 probe sets representing over 23,998 bovine transcripts, including assemblies from approximately 19,000 UniGene Clusters. In order to assess the influence of the two cycles of linear amplification, on the representation of original transcripts, we compared microarray experiments from one-cycle and two-cycle amplifications using total RNA from DC1. The results showed that amplification of messages using 1 vs. 2 cycles were highly consistent with a correlation coefficient of 0.93 (data not shown). These data confirm the manufacturer's results using 1 and 2 cycles of linear amplification.
Microarray experiments were performed in three biological replicates for all blastocysts (CT1, CT4 and IVF) and donor cells (DC1 and DC4). Images were processed with GCOS and data extracted using MAS 5.0. However, one of the CT1 blastocyst chips did not pass the quality control analysis  and was excluded from the study. The analyses for CT1 are based on the remaining two chips in this group, which showed an appropriate p-value distribution. The GCOS software expression data report showed that 56% of the probe sets were called "Present" (P) for all donor cell chips. This number was lower for all blastocyst chips with 44%, 41%, and 47% for IVF, CT1, and CT4 respectively. Probe sets that were called "Absent" (A) in all the samples were excluded from the analysis. Therefore only 16,521 probe sets were included in the analysis.
Number of differentially expressed transcripts in pairwise comparisons between IVF embryos, CT1 embryos, CT4 embryos, DC1 cells, and DC4 cells (p-value < 0.01 and fold change > 2.0)
Group 1 vs. Group 2
Differentially expressed transcripts
Higher in the first group
Higher in the second group
IVF embryos vs. CT-1 embryos
IVF embryos vs. CT-4 embryos
IVF embryos vs. DC-1 cells
IVF embryos vs. DC-4 cells
CT-1 embryos vs. CT-4 embryos
CT-1 embryos vs. DC-1 cells
CT-1 embryos vs. DC-4 cells
CT-4 embryos vs. DC-1 cells
CT-4 embryos vs. DC-4 cells
DC-1 cells vs. DC-4 cells
Classification of differentially expressed probe sets in pairwise comparisons
Probe set category
IVF vs . CT1
IVF vs . CT4
IVF vs . DC1
IVF vs . DC4
CT1 vs . CT4
CT1 vs . DC1
CT1 vs . DC4
CT4 vs . DC1
CT4 vs . DC4
DC1 vs . DC4
Transcripts with strong similarity to a known gene
Transcripts with moderate similarity to a known gene
Transcripts with weak similarity to a known gene
Top 30 upregulated genes, in IVF blastocysts compared to CT blastocysts sorted by P-value
Probe set ID
Fold change IVF/CT1
Fold change IVF/CT4
Peptidase inhibitor 3, skin-derived (SKALP)
Polo-like kinase 3 (Drosophila)
Zinc finger protein 183
Fatty acid binding protein 7, brain
ATPase, H+ transporting, lysosomal V0 subunit a1
Heat shock 70 kDa protein 1A
Bcl-2 inhibitor of transcription
Zinc finger protein 410
LSM1 homolog, U6 small nuclear RNA associated
Testis expressed 12
Nuclear factor, interleukin 3 regulated
Sodium channel modifier 1
Signal transducer and activator of transcription 3
Serpin peptidase inhibitor, E member 2
Sjogren's syndrome/scleroderma autoantigen
DEAD (Asp-Glu-Ala-Asp) polypeptide 56
General transcription factor IIF, polypeptide
Ubiquitin-conjugating enzyme E2A (RAD6 homolog)
BUB3 budding uninhibited by benzimidazoles 3
Transmembrane 4 superfamily member 1
Intraflagellar transport 20 homolog (Chlamydomonas)
Nucleoside-diphosphate kinase NBR-A
Keratin 10 (epidermolytic hyperkeratosis)
Collaborates/cooperates with ARF protein
High mobility group nucleosomal binding domain 3
TSR2, 20S rRNA accumulation, homolog
Sin3A-associated protein, 18 kDa
Top 30 genes upregulated in CT blastocysts and donor cells compared to IVF blastocysts, sorted by P-value
Probe set ID
Fold change CT1/IVF
Fold change CT4/IVF
Adaptor-related protein complex 3, sigma 2
X-ray repair complementing defective repair in Chinese hamster cells 1
insulin receptor substrate 4
Crystallin, lambda 1
Nuclear casein kinase and cyclin-dependent kinase substrate 1
Arginine/serine-rich coiled-coil 1
thyroid hormone receptor associated protein 5
vacuolar protein sorting 26 homolog A
ARP3 actin-related protein 3 homolog (yeast)
Nuclear autoantigenic sperm protein
phosphatidylinositol binding clathrin assembly protein
Palladin, cytoskeletal associated protein
Phenylalanyl-tRNA synthetase 2, mitochondrial
Mitochondrial ribosomal protein S35
Aldehyde oxidase 1
Eukaryotic translation initiation factor 4H isoform 2
Tumor differentially expressed 2-like
Coagulation factor II receptor-like 1
Casein kinase 1, epsilon
Chromosome 1 open reading frame 35
Interferon induced transmembrane protein 1
Anti-oxidant protein 2 (independent phospholipase A2)
Heat shock 90 kDa protein 1, beta
LAG1 homolog, ceramide synthase 2
Primers used for Real time PCR validation
Primer sequences and positions (5' - 3')
Fragment size (bp)
Genes with putative cummulative downregulation in blastocysts obtained after serial rounds of chromatin tranfer
Probe Set ID
Fold change IVF/CT1
Fold change IVF/CT4
heat shock 70 kD protein 1
neuroguidin, EIF4E binding protein
high mobility group nucleosomal binding domain 3
neuroguidin, EIF4E binding protein
intraflagellar transport protein 20
Bcl-2 inhibitor of transcription
Sjogren's syndrome/scleroderma autoantigen 1
TSR2, 20S rRNA accumulation, homolog (S. cerevisiae)
chromosome 14 open reading frame 10
collaborates/cooperates with ARF (alternate reading frame) protein
sodium channel modifier 1
myosin, heavy polypeptide 7, cardiac muscle, beta
proton-dependent gastrointestinal peptide transporter
protease inhibitor 3, skin-derived (SKALP)
hypertension-related calcium-regulated gene
dispatched homolog 1 (Drosophila)
SERTA domain containing 1
transition protein 1 (during histone to protamine replacement)
microtubule-associated protein, RP/EB family, member 2
oviduct specific glycoprotein
neurofilament, medium polypeptide
glycoprotein (transmembrane) nmb
fucosyltransferase 10 (alpha (1,3) fucosyltransferase)
solute carrier family 6 (neurotransmitter transporter, dopamine), member 3
membrane-associated ring finger (C3HC4) 2
Bernardinelli-Seip congenital lipodystrophy 2
RPGR-interacting protein 1
DNA replication factor
Nucleolar protein family A, member 1
purinergic receptor P2Y G-protein coupled, 2
nucleotide phosphodiesterase, 3'-5'-cyclic
anterior gradient 2 homologue
Genes with putative cummulative upegulation in blastocysts obtained after serial rounds of chromatin tranfer
Probe Set ID
Fold change IVF/CT1
Fold change IVF/CT4
NADH dehydrogenase (ubiquinone) Fe-S protein 2, 49 kDa (NADH-coenzyme Q reductase)
palladin, cytoskeletal associated protein
non-muscle myosin heavy chain
guanidine nucleotide binding protein, (G protein), alpha inhibiting activity polypeptide 2
golgi autoantigen, golgin subfamily a, 7
zinc finger protein 313
chromatin modifying protein 1B
mitochondrial carrier homolog 2
F-box protein 9
methylmalonyl Coenzyme A mutase
leukotriene B4 12-hydroxydehydrogenase
adaptor-related protein complex 3, sigma 2 subunit
Abl-philin 2 isoform 2
thrombospondin repeat containing 1
lysosomal-associated membrane protein 1
adrenergic, beta 3, receptor
myosin, heavy polypeptide 10, non-muscle
trophoblast Kunitz domain protein 1
interferon-stimulated protein, 15 kDa
5,10-methylenetetrahydrofolate reductase (NADPH)
prion protein interacting protein
Guanine nucleotide binding protein (G protein), alpha stimulating activity polypeptide 1
Zinc finger protein 325 (gonadotropin inducible transcription repressor-3)
Guanine nucleotide binding protein (G protein), alpha stimulating activity polypeptide 1
interferon-stimulated protein, 15 kDa
It has been reported that in vitro culture conditions alter gene expression and may lead to developmental aberrations in IVF derived cattle, commonly referred to as the large offspring syndrome [38–40]. In the case of embryos produced by SCNT, besides the alterations due to in vitro culture conditions, gene expression defects may be caused by improper silencing and activation of specific genes, altered chromatin remodelling, and epigenetic alterations . But identifying key genes responsible for the general developmental failure in cloned embryos is not an easy task, since the alterations may be caused by a variety of factors including donor cell type, cell cycle stage, nuclear transfer protocol, source of the oocytes, embryo culture system, embryo transfer procedure, recipient cows management, and operators' skills . Consequently, there is a big variety of alterations that are not always shared by all cloned embryos. Still, the common thread uniting many of the SCNT failures can be traced to epigenetic alterations, specifically failures in chromatin remodelling and DNA and histone methylation [13, 43, 44]. The fetal fibroblast cells used in this study are not fully representative of adult somatic cells. However, these cells were chosen because of their practicality and higher efficiency in SCNT studies.
Microarray analysis has been used to explore the transcriptome profile of cloned embryos relative to that of the donor cells and IVF embryos as a control. However, the appropriate microarray platform is crucial in order to detect changes in particular genes. Smith and colleagues reported similar transcriptome profiles for cloned blastocysts and blastocysts produced by artificial insemination . However, the cDNA microarray used by Smith and colleagues consisted of placenta and spleen cDNA libraries, lacking embryonic genes, which therefore were not analysed. The results from the present study show an extensive reprogramming in cloned embryos by the blastocyst stage. However, the data point to a group of differentially expressed transcripts between IVF and cloned blastocysts.
Serial cloning is often performed for the production of transgenic animals. Although apparently healthy animals can be obtained after serial cloning, the efficiency of cloning decreases from generation to generation despite comparable blastocyst and early pregnancy rates. This increase in pregnancy losses and perinatal deaths could be caused by gene expression defects accumulated throughout the serial cloning procedures, which could be detected in blastocysts, although no phenotypic alterations are observed at this stage. Furthermore, it has been proposed that the extended culture, associated with transfection and selection procedures, may induce changes in the donor cells . The present studies show that serial cloning does not significantly affect transcriptional reprogramming of cloned blastocysts. The global transcriptome profile of blastocysts from four consecutive rounds of cloning did not significantly differ from the one obtained from blastocysts after only one round of cloning. However, for a set of genes, misregulation was significantly greater in the blastocysts obtained from four rounds of cloning (see Tables 6 and 7). However the observed differences between blastocysts from the first and fourth rounds of cloning could be due to the fact that these are different donor cells and not of the same clonal origin.
To our knowledge this is the first study to focus on the influence of serial chromatin transfer on global transcriptome profile of embryos and donor cells. Only a small proportion of the data set generated by the present study corresponded to fully annotated bovine genes (Table 2). The rest of the probe sets were excluded from further analyses due to lack of annotations. Progress in the annotation of the bovine genome will greatly facilitate global gene expression studies in the bovine species.
In the present study, multiple comparisons revealed five distinctive patterns of differential gene expression among all embryos and donor cells. The first pattern corresponded to 1,183 transcripts (30.74% of the data set) that had similar abundance in all five groups. Housekeeping genes like GAPD and Actin showed this pattern of expression. The second pattern corresponded to genes that had similar expression in IVF and CT embryos, but had a very different pattern of expression in both donor cell lines. We hypothesised that these are genes that switched from the "donor cell gene expression mode" to the "embryo gene expression mode". The majority of the genes in the data set (76.49%) showed this pattern, including some imprinted and embryonic specific genes such as the Oct-4 protein coding gene (POU5F1), which has been reported as differentially expressed for cloned embryos in previous studies [21, 23]. Placenta specific 8 (PLAC8) also shows this pattern of expression (Figure 6A). It is possible that some genes, due to their methylation pattern in the somatic cells or to their location in the chromosome, are more likely to be reprogrammed by the oocyte factors.
The third pattern corresponded to genes with a similar pattern of expression for CT embryos and donor cells, and a very different expression pattern in IVF embryos. These were 147 (3.81%) genes with apparently incomplete reprogramming, probably with a somatic cell pattern of expression. The heat shock 70 kD protein 1 (HSPA1A), involved in cell protection from stress and apoptosis was significantly higher in IVF embryos when compared to CT embryos and donor cells (Figure 6B). Important embryonic genes showed this pattern of expression. Desmocollin 3 (DSC3) a transmembrane glycoprotein, involved in cell adhesion that belongs to the cadherin family, was present in IVF embryos but was absent in CT embryos and donor cells. The signal transducer and activator of transcription 3 (STAT3), was significantly upregulated in IVF embryos when compared to both groups of cloned embryos and donor cells. A similar pattern was observed for high mobility group nucleosomal binding domain 3 (HMGN3) a gene involved in chromatin remodelling, a vital process during embryonic genome activation (Figure 6C). The importance of both genes during morulae and blastocyst formation could make them good candidates in understanding the lower developmental rates of cloned embryos.
The fourth group of genes corresponded to only 85 probe sets (2.21%) with a marked differential expression in all cloned embryos compared to the one observed in both IVF embryos and donor cells. The misregulation of these genes could point to a compensation mechanism after chromatin transfer. Genes with this kind of expression pattern included prostaglandin-endoperoxide synthase (PTGS2) and the transcription factor GATA-2. Both genes had a greater microarray signal in all CT embryos, but low expression in IVF and donor cells. The imprinted gene glycine amidinotransferase (GATM), showed significantly greater values in the cloned embryos compared to IVF embryos and donor cells. Two interesting genes in this group were DNMT3A and DNMT3B transcripts, which are responsible for de novo methylation. Both genes were significantly greater in CT-1 and CT-4 embryos compared to IVF blastocysts (Figure 6D and 6E), which is consistent with the hypermethylation often reported in cloned blastocysts. These results do not agree with previous findings, in which DNMT3A was downregulated in NT embryos compared to IVF embryos . Zhou et al., reported similar levels of DNMT3B for embryos produced in vivo, in vitro, and by different nuclear transfer methods, including chromatin transfer . These contrasting results confirm that alterations greatly vary and are not shared by all cloned embryos. One limitation of our study is that we have not used in vivo blastocysts which might have provided more biological means and as the physiological standard against in vitro culture conditions.
A fifth pattern corresponded to genes that had an increasing or a decreasing pattern of expression from IVF embryos through donor cells showing an intermediate pattern of expression in CT embryos. In total, 245 probe sets showed this pattern of expression with 119 (3.08%) increasing, and 126 (3.28%) decreasing from IVF through DC. It could be assumed that these genes have been partially reprogrammed, since their transcript abundance is in between IVF and donor cells. The imprinted gene insulin-like growth factor 2 receptor (IGF2R), one of the most studied genes in the large offspring syndrome, showed similar expression values in IVF and CT1 embryos, but significantly higher signals in CT4 embryos, and very high signals in both donor cells (Figure 6F). These higher mRNA levels in the fourth generation of cloning could indicate a cumulative misregulation of this gene. The Bcl-2 inhibitor of transcription (BIT1) showed the greatest values in IVF embryos, intermediate values in CT embryos and the lowest values in donor cells (Figure 7A). The nuclear transcription factor Y, alpha (NFYA), showed a similar expression pattern in both IVF and CT1 embryos; although it was significantly lower in CT4 embryos and donor cells. Neuroguidin (NGDN), an eukaryotic translation initiation factor with important functions in embryonic development was another gene with a decreasing pattern of expression (Figure 6B). Genes with and increasing pattern of expression included F-box protein 9 (FBXO9), and guanine nucleotide binding protein alpha inhibiting activity polypeptide 2 (GNAI2) represented in Figure 6C and Figure 6D, respectively. Transcripts for the progesterone receptor (PGR) were significantly higher in IVF embryos compared to CT embryos and donor cells (Figure 6E). Among this group of transcripts could be genes that are cumulatively affected be serial cloning.
Based on the difference in gene expression for RARB, CRAB1, THBS, SERPINB5, and HLA-A, Beyhan et al. suggest a possible role for the retinoic acid signalling pathway in the failures observed in cloned bovine embryos . However the bovine GeneChip does not contain a Retinoic Acid Receptor Beta (RARB) probe set. It only contains a probe set that corresponds to a bovine EST with similarity to the rat RARB (Bt.21044.2.A1_at). In the present data, CRAB1 and THBS2 were slightly higher in IVF embryos, although without statistical significance. They also found differential gene expression among several genes in both donor cells (CDKN1C, COPG2, DCN, GATM, MEST, NDN, NNAT, PON3, and SGCE). In the current study GATM was significantly downregulated in donor cells from the fifth successive generation of chromatin transfer (Figure 7).
At the blastocyst stage there is an extensive reprogramming of cloned embryos leading to very similar transcriptomes in IVF and CT blastocysts. However, there were around 200 differentially expressed genes in both CT embryos compared to IVF. For some genes, the differences were significantly greater in CT4 when compared to CT1, suggesting a possible cumulative missregulation caused by serial cloning. Genes involved in transcription, cellular proliferation, embryonic development, cellular death, and response to stress are over represented in IVF embryos; many of these genes are present in the nucleus, which was the cell component overrepresented in IVF embryos. Genes involved in cell morphology, cell development, and metabolism were over expressed in donor cells and in cloned embryos when compared to IVF, suggesting that they were not properly silenced in the donor nucleus. The up regulation of genes involved in metabolism should be further explored as it could be linked to the large size of cloned animals.
As gene expression profile can only show one step in cell phenotype and function control, namely transcriptome regulation, proteomic analysis could complement this study by providing a more complete picture of the regulation of embryonic development. With a more complete bovine genome annotation, more of the differentially expressed transcripts could be analyzed further providing more information for the currently unidentified transcripts, which, in the present study represented around 18% of the dataset. Gene Ontology information for a proportion of the differentially expressed genes is still incomplete. Thus, for some of the genes the cellular component is known, but the biological process and/or its molecular function is not documented. It is interesting that the majority of genes upregulated in CT blastocysts participate in metabolism processes, while the percentage of metabolism genes in IVF blastocyst was lower compared to signalling pathway genes.
The present study provides a data set that could be useful in identifying epigenetic errors in cloning and may facilitate our understanding of the reprogramming process in SCCT. Future studies should involve more of the successive generations of cloned embryos and their respective donor cells to identify cumulative misregulated genes. Gene expression studies from fetal, newborn, and placental tissues could identify genes that are responsible for abnormalities, abortions, stillborns and low birth rate. Functional studies should target particular genes that play key roles in molecular reprogramming and early embryo development and manipulate their mRNA concentrations in SCCT embryos, to mimic that of IVF embryos.
Bovine oocytes were aspirated from 2–8 mm follicles of abattoir-obtained ovaries from Holstein cows and matured in Tissue Culture Medium (TCM-199, Gibco/Invitrogen, Grand Island, NY) supplemented with 0.2 mM pyruvate, 0.5 μg/ml FSH (Sioux Biochemicals, Sioux City, IA), 5 μg/ml LH (Sioux Biochemicals, Sioux City, IA), 10% FCS (Gibco/Invitrogen, Grand Island, NY), 100 U/ml penicillin and 100 μg/ml streptomycin (Gibco/Invitrogen, Grand Island, NY) in 5% CO2 in air at 38.5°C. For fertilization, matured oocytes were transferred to fertilization medium and were fertilized using thawed sperm from a Holstein bull separated by Percoll density gradient and further incubated for 24 hours. Presumptive zygotes were transferred to Gardner's culture medium 1 (G1) for 3 days, followed by 3–4 days culture in Gardner's culture medium 2 (G2). Blastocysts were evaluated and graded according the International Embryo Transfer Society (IETS) guidelines . Grade 1 blastocysts were selected, pooled in groups of 3 blastocysts per tube, frozen (with addition of lyses buffer from RNeasy MicroKit (Qiagen Valencia, CA) in liquid nitrogen and stored in -80°C until RNA isolation.
In vitro-matured oocytes were enucleated at 20 hpm. Bovine fetal fibroblasts after one and four rounds of cloning were trypsinized and washed in Ca/Mg Hank's Balanced Salt Solution (HBSS) and permeabilized by incubation of 50,000 - 100,000 cells in 31.25 units Streptolysin O (SLO-Sigma, St. Louis, MO) in 100 μl for 30 minutes in a 37°C H2O bath. Permeabilized fibroblasts were washed, pelleted and incubated in 40 μl of mitotic extract prepared from MDBK cells containing an ATP-generating system (1 mM ATP, 10 mM creatine phosphate and 25 μg/ml creatine kinase) for 30 min at 38°C. At the end of incubation, the reaction mix was diluted with 500 μl of cell culture media (Alpha MEM with 10% FBS), pelleted and resuspended in TL Hepes. These cells were fused to enucleated oocytes, activated 26 h after maturation with 5 μM calcium ionophore for 4 min followed by 10 μg/ml of cycloheximide and 2.5 μg/ml of cytochalasin D for 5 h. After activation, embryos were washed, and cultured in SOF medium for the first 4 days with 8 mg/ml BSA and the last three days with 10% fetal calf serum at 38.5°C and 5% CO2 in air. Grade 1 blastocysts were pooled (3 per tube) and frozen, with addition of lysis buffer. Embryos were stored in -80°C until RNA isolation.
For subsequent rounds of cloning, CT derived bovine blastocysts from the first generation were transferred into hormonally synchronized cows. At seventy-days, pregnancies were interrupted, and foetuses recovered. Fetal fibroblast cultures were established and used for the next chromatin transfer process. The same procedure was done 3 times to provide a fourth generation of clones. Grade 1 blastocysts from the fourth generation were pooled (3 per tube) and frozen, with addition of lysis buffer. Embryos were stored in -80°C until RNA isolation.
Fetal cell lines were developed according. Seventy-day old male bovine foetuses were recovered and transported to the laboratory in Dulbecco's PBS (DPBS) with 16 ml/ml of antibiotic-antimycotic (Gibco, Grand Island, NY), 4 ml/ml tylosin tartrate (Sigma, St. Louis, MO), and 8 ml/ml fungizone (Gibco). Foetuses were rinsed in DPBS, the head and internal organs were removed, and remaining tissues were finely chopped into pieces with a scalpel blade. The fibroblasts were separated from the tissue pieces using 0.08% trypsin and 0.02% EDTA in PBS (trypsin-EDTA). The cells were seeded onto 100-mm tissue culture plates (Corning, VWR, Chicago, IL) in a minimal essential medium (a-MEM; Gibco) supplemented with 10% fetal calf serum (FCS; Hyclone, Logan, UT), 0.15 g/ml glutamine (Sigma), 0.003% b-mercaptoethanol (Gibco), and antibiotic-antimycotic (Gibco). On the same day of cloning (day 3 of seeding), the cells were harvested using DPBS with trypsin-EDTA solution and were counted. One million cells were frozen in MEM with 10% FCS, dimethyl sulfoxide (Sigma), and lysis buffer.
RNA was isolated from IVF blastocysts, SCCT blastocysts, and donor cells using the RNeasy MicroKit (Qiagen Valencia, CA) according to the manufacturer's specifications. Briefly, embryos and cells frozen at -80°C in lysis buffer were transferred to silica-gel membrane spin columns and washed with RW1 wash buffer and 80% ethanol. Final RNA elution was conducted using 14 μl of RNAse free water provided in the kit. Concentration and purity of isolated RNA were determined using a NanoDrop® ND-1000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE). Integrity and quality were analyzed using a Bioanalyzer 2100 RNA 6000 Picochip kit (Agilent Technologies, Palo Alto, CA).
Microarray hybridizations were performed in triplicate for each of the experimental groups using Affymetrix Bovine DNA Chips as described by the manufacturer (Affymetrix Santa Clara, CA). Briefly, complementary DNA (cDNA) synthesis was performed from 10 ng total RNA using the Two-Cycle cDNA Synthesis Kit (Affymetrix Santa Clara, CA). The MEGAscript® T7 Kit (Ambion, Inc.) was used for the first in vitro transcription (IVT). GeneChip IVT Labeling Kit was used for the second IVT and labelling of RNA. Complementary RNA (cRNA) was fragmented and 10 μg of fragmented cRNA were hybridized to the Genechips in a Hybridization Oven, set to 45°C and rotations of 60 rpm for 16 hours. The chips were then washed and stained with streptavidin/phycoerythrin (SAPE) antibody solution using an Affymetrix FS-450 fluidics station. GeneChips were scanned using the Affymetrix GeneChip scanner 3300.
Images were processed with the Affymetrix GeneChip® Operating Software (GCOS) and expression quantified with MAS 5.0, which also provides information on signal, detection and calculated the detection p-value. Signal information is a numeric value indicating transcript abundance for a particular probe set. Detection information indicates whether the transcript is detected (P, present), undetected (A, absent), or if it is at the limit of detection (M, marginal). Detection p-value indicates the significance of the detection call for a probe set. Only probe sets that were called Present in at least one of the five groups were included in the analysis. A total of 5,599 probe sets were excluded from the analysis as they were called Absent in all groups. The data set for further analysis included 18,396 probe sets.
Metrics like noise, background, Scale factor, and the ratio of intensities of 3' probes to 5' probes for Actin and GAPDH genes were analyzed for chip quality control. Information about the intensities of the spiked in controls (B. subtilis genes lys, phe, thr, and dap), which were mixed with the total RNA at known concentrations at the beginning of the experiment was used to monitor the linear amplification and labelling process independently from the target samples. The performance of the hybridization control genes (E. coli genes BioB, BioC and BioD and P1 Bacteriophage cre) was also used for determining the quality of each chip.
For data visualization, the raw GeneChip signals were uploaded into GeneTraffic UNO (Iobion Informatics LLC), which generated scatter plots of pairwise hybridization comparisons and Heat maps from all hybridizations using hierarchical clustering. Power Atlas, a web-based resource from the University of Alabama at Birmingham, was used to estimating the power of the hybridization given the sample size . HDBStat was used for statistical analysis. Data were quantile-quantile normalized and examined for outliers using Person's correlation. Quality control statistics included a deleted residuals approach [37, 47, 48]. False discovery rates (FDR) for the genes were calculated using t-test . Fold changes were calculated based upon the unadjusted data means in pairwise comparisons. Probe sets in each pairwise comparison with a p < 0.01, and FDR of <20%, and a Fold Change (FC) in excess of 2.0 were considered to be significant and examined further. For multiple comparisons, One-way analysis of variance (ANOVA) from PROC GLM in SAS 9.1 (SAS Institute inc. Carey, NC) was performed on the complete data set. The Least Significant Difference (LSD) test was used to detect significant differences between groups.
The probe sets corresponding to differentially expressed genes were uploaded into the Affymetrix Netaffx Analysis Center (Bovine GeneChip annotation from November 6 2007) to retrieve updated information regarding gene symbol, gene title, Biological Process (BP), Molecular function (MF), and Cellular Component (CC) . To complement the annotation from Netaffx, we used the GOAnna tool  from AgBase, a Mississippi State University curated, web-accessible resource for functional analysis of agricultural plant and animal gene products. For data visualization, all the GO terms associated to each gene were uploaded into GOSlimViewer  another AgBase tool that provides a high level summary of the GO categories found in the dataset allowing a better visualization of the data.
Ingenuity Pathway Analysis 5.0 from Ingenuity Systems was used for data modelling and the analysis of networks related to the generated data sets. Genes upregulated in IVF embryos compared to CT embryos and donor cells (Figure 8) and genes downregulated in IVF embryos compared to CT embryos and donor cells (Figure 9) were uploaded in the Ingenuity Pathway Analysis 5.0. Since Ingenuity Pathway Analysis database is based on human, mouse, and rat genes, some of the bovine names were not recognized by the software, mostly because of different gene symbols. For those genes, we manually identified the human orthologous symbol.
DNA microarray derived gene expression results for genes DNMT3A, DNMT3B, IGF2R, PLAC8, PGR, BIT1, HMGN3, HSPA1A, NGDN, FBXO9, and GNAI2 were confirmed by Real time PCR using GAPDH as the reference gene. Complementary DNA was generated with the First-Strand cDNA Synthesis system for RT-PCR using SuperScript III Platinum® Two-Step qRT-PCR Kit (Invitrogen Life Technologies, Carlsbad, CA) according to the manufacturer's protocol. The samples were incubated for 10 min at 25°C, 50 min at 42°C and at 85°C for 5 min. Then 2 U of E. coli Rnase H was added to each tube and incubated at 37°C for 20 min. The cDNA was used for quantitative real-time PCR amplification with SYBR Green I chemistry (Roche Applied Sciences, Indianapolis, IN). Real-time quantitative PCR was performed using the LightCyclerTM instrument (Roche Applied Sciences, Indianapolis, IN). The real time PCR reactions were carried out in a total volume of 10 μl according to the manufacturer's manuals for DNA Master SYBR Green I mix (Roche Applied Sciences, IN). The primer concentrations were adjusted to 0.5 μM for each gene. Primers were designed using Primer Premier 5 software (Premier Biosoft International, Palo Alto, CA). Primer sequences used for real time PCR are shown in Table 1. The cycling parameters were 30 seconds at 95°C for denaturation, 50 cycles of 2 seconds at 95°C, 10 seconds at 55°C for amplification (quantification was performed at this step), and 12 seconds at 72°C for extension. The specifiCity of all individual amplification reactions was confirmed by melting curve analysis. Real-time expression values were calculated through the relative standard curve method, using 10-fold serial dilutions for both the target and the endogenous reference genes by measuring the cycle number at which exponential amplification occurred in a dilution series of samples. Values were normalized to the relative amounts of the control mRNA, which were obtained from a similar standard curve. In real time PCR reactions, the same initial amounts of target molecules were used, and the Cp values of control mRNA were constant in all samples.
Donor cell lines included in the study were fibroblasts from non-cloned foetuses (DC0), and fetal fibroblasts from first, second, fourth, and fifth rounds of cloning (DC1, DC2, DC4, and DC5). RNA isolation from donor cells and subsequent cDNA synthesis were performed according to the above mentioned protocols. Relative mRNA abundance was determined for paladin (PALLD), nuclear transcription factor Y alpha (NFYA), glycine amidinotransferase (GATM) and Taspase 1 (C20orf13). Quantitative assessment of RNA amplification was detected by SYBR® GreenER™ qPCR SuperMixes for iCycler (Invitrogen Life Technologies, Carlsbad, CA, 11761-100). Real-time PCR reactions were performed using the iCycler iQ Real-Time PCR instrument (BIO-RAD). The cycling parameters were 50°C for 2 min, 95°C for 8 min 30 s for denaturation, 40 cycles of 15 s at 95°C and 30 s at 60°C and 30 s at 72°C for amplification and extension respectively. The melting curve was performed starting at 55°C with a 0.5°C increase for 10 s in 80 cycles. Expression values were calculated using the relative standard curve method. Standard curves were generated using 10-fold serial dilutions for both GAPDH and 18S ribosomal RNA. Standard curves were also generated for all target genes by measuring the cycle number at which exponential amplification occurred.
Results from different groups were analyzed by one-way analysis of variance (ANOVA) by SAS 9.1 (SAS Institute inc. Carey, NC). Differences at p < 0.001 were considered significant. An additional analysis was performed using Relative expression software tool (REST©, 384-beta version May 2005) to compare all samples of each group. The mathematical model used in the REST software is based on the PCR efficiencies (E) and the crossing point deviation between the samples (CP) [53–55].
Bovine Serum Albumin
Fetal Calf Serum
in vitro fertilization
Somatic cell nuclear transfer
Somatic cell chromatin transfer
This study was funded by the Life Sciences and Biotechnology Institute (LSBI), Mississippi Agricultural and Forestry Experiment Station (MAFES) Mississippi State University (Assigned manuscript number: J-11231).
The authors acknowledge: Ms. Juliet Tang at the Life Sciences and Biotechnology Institute for technical assistance in microarray experiments; Dr. Fiona McCarthy of the College of Veterinary Medicine, Mississippi State University for technical assistance in gene ontology annotation with AgBase; Angie Edwards for her contribution to cell culture, shipping materials and scientific discussions, and Rachel Paulson for her contribution in coordinating scientific meetings between Hematech Inc. and the Memili lab.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.