Changes in the transcriptome of morula-stage bovine embryos caused by heat shock: relationship to developmental acquisition of thermotolerance
© Sakatani et al.; licensee BioMed Central Ltd. 2013
Received: 20 November 2012
Accepted: 11 January 2013
Published: 15 January 2013
While initially sensitive to heat shock, the bovine embryo gains thermal resistance as it progresses through development so that physiological heat shock has little effect on development to the blastocyst stage by Day 5 after insemination. Here, experiments using 3’ tag digital gene expression (3’DGE) and real-time PCR were conducted to determine changes in the transcriptome of morula-stage bovine embryos in response to heat shock (40 degrees C for 8 h) that could be associated with thermotolerance.
Using 3’DGE, expression of 173 genes were modified by heat shock, with 94 genes upregulated by heat shock and 79 genes downregulated by heat shock. A total of 38 differentially-regulated genes were associated with the ubiquitin protein, UBC. Heat shock increased expression of one heat shock protein gene, HSPB11, and one heat shock protein binding protein, HSPBP1, tended to increase expression of HSPA1A and HSPB1, but did not affect expression of 64 other genes encoding heat shock proteins, heat shock transcription factors or proteins interacting with heat shock proteins. Moreover, heat shock increased expression of five genes associated with oxidative stress (AKR7A2, CBR1, GGH, GSTA4, and MAP2K5), decreased expression of HIF3A, but did not affect expression of 42 other genes related to free radical metabolism. Heat shock also had little effect on genes involved in embryonic development. Effects of heat shock for 2, 4 and 8 h on selected heat shock protein and antioxidant genes were also evaluated by real-time PCR. Heat shock increased steady-state amounts of mRNA for HSPA1A (P<0.05) and tended to increase expression of HSP90AA1 (P<0.07) but had no effect on expression of SOD1 or CAT.
Changes in the transcriptome of the heat-shocked bovine morula indicate that the embryo is largely resistant to effects of heat shock. As a result, transcription of genes involved in thermal protection is muted and there is little disruption of gene networks involved in embryonic development. It is likely that the increased resistance of morula-stage embryos to heat shock as compared to embryos at earlier stages of development is due in part to developmental acquisition of mechanisms to prevent accumulation of denatured proteins and free radical damage.
The newly-formed preimplantation embryo is very sensitive to disruption of its developmental program by heat shock. In the mouse  and cow , exposure to a temperature that is only 1.5-2.0°C above normal body temperature can disrupt development of 1-cell embryos. While initially sensitive to heat shock, the embryo gains thermal resistance as it progresses through development. Heat shock at the pronuclear stage reduced development of mouse embryos more than heat shock at the 2-cell stage . In the cow, heat shock was more detrimental to development when applied at the 8-cell stage or earlier than when applied later in development [3–6]. Development of bovine zygotes was reduced by a heat shock of 40°C . Moreover, adverse effects of exposure of females to heat stress on embryonic survival become reduced after a few days post-fertilization in sheep , pigs  and cattle .
Two phenomena appear important for developmental acquisition of thermotolerance in the preimplantation embryo. The first is a change in oxidative stress. Free radical production is an important cause of reduced development in heat-shocked embryos; effects of heat shock in vitro [10–13] and heat stress in vivo [14, 15] can be reduced in some cases (although not all) [16, 17] by administration of antioxidants. In the cow, free radical production increases in response to heat shock when embryos were heat-shocked at Day 0 and 2 relative to fertilization but not when heat shock was imposed at Day 4 or 6 . Amounts of the antioxidant glutathione are also low during early development and do not increase until the 9–16 cell stage . Another reason for the developmental gain in resistance of embryos to heat shock is the activation of the embryonic genome to allow stimulation of cellular pathways that protect cells from stress. In both the mouse and bovine, embryos become more resistant to heat shock at about the same time embryonic genome activation takes place (1–2 cell stage in the mouse  and 8–16 cell stage in the cow ). The phenomenon of induced thermotolerance, whereby exposure to a mild heat shock induces changes in cellular function to make cells more resistant to a subsequent, more severe heat shock, can occur in bovine blastocysts  but not in 2-cell embryos . Much of the change in cellular function leading to thermotolerance involves transcriptional activation of genes for heat shock proteins and other proteins .
The purpose of this study was to evaluate changes in the transcriptome of morula-stage bovine embryos in response to heat shock. The goal was to identify genes and gene networks affected by heat shock that would explain the molecular basis for thermal resistance at this stage of development. The transcriptome was assessed using a next-generation sequencing method called 3’ tag digital gene expression (3’DGE). In this procedure, high throughput sequencing is used to read 20–21 bp amplicons of 3’ regions of mRNAs [24, 25]. The procedure is particularly useful for cases where mRNA is limited because more gene counts are generated than for procedures in which the whole transcriptome is analyzed .
The fertilization medium SOF-FERT was made by modifying a custom formulation of Synthetic Oviduct Fluid manufactured by Specialty Media (Millipore, Billerica, MA, USA) to create the following formulation: 107.7 mM NaCl, 7.16 mM KCl, 1.19 mM KH2PO4, 1.17 mM CaCl2·H2O, 5.3 mM sodium lactate, 0.20 mM sodium pyruvate, 25.07 mM NaHCO3, 0.49 mM MgCl2·6H2O, 6 mg/ml essentially fatty acid free BSA (Sigma, St. Louis, MO, USA), 5 μg/ml gentamicin, 10 μg/ml heparin, and 1.0 mM caffeine. The embryo culture medium was SOF-BE1 and was prepared as described elsewhere . Unless otherwise mentioned, all other chemicals and reagents were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Production of embryos in vitro
Bovine embryos were produced from slaughterhouse-derived immature oocytes using procedures for in vitro oocyte maturation and fertilization as described previously . Oocytes were inseminated with a pool of sperm from three Bos taurus bulls. A different pool was used for each replicate. After 8 h of insemination, cumulus cells were removed from putative zygotes by vortexing and the zygotes were cultured in groups of 30 in 50 μl microdrops of SOF-BE1 covered in mineral oil at 38.5°C and 5% CO2 in humidified air. Embryos were cultured until 116 h after insemination (hpi), when they were used for experiments.
Changes in the transcriptome caused by heat shock as determined by 3’DGE
Treatment and RNA purification
At 116 hpi, embryos were either maintained at 38.5°C (control) or moved to an incubator at 40°C (heat shock) and 5% CO2 in humidified air for 8 h. Morulae (defined here as embryos > 16 cells) were collected and the zona pellucidae removed by incubation with 0.1% (w/v) protease from Streptomyces griseus in Dulbecco’s phosphate buffered saline (DPBS) containing 0.2% (w/v) polyvinylpyrrolidone (PVP). Embryos were then washed with DPBS-PVP, transferred in groups of 50 embryos to a 1.5 mL tube containing 50 μL extraction buffer from the PicoPure RNA isolation kit (Applied Biosystems, Carlsbad, CA, USA), incubated at 42°C for 30 min, and centrifuged at 3,000 g for 2 min. The supernatant fractions were used to prepare total RNA using the PicoPure RNA isolation kit. DNA was digested using the RNase-Free DNase Set (Qiagen, Valencia, CA, USA). The concentration and quality of total RNA were assessed by a Nanodrop ND-1000 (ThermoScientific, Wilmington, DE, USA) and Agilent 2000 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Samples were only used for 3’DGE if RNA integrity was > 8.0. A total of three replicates of 50 morulae each that met this criterion were obtained for each treatment.
RNA amplification and 3’-DGE sequencing
Total RNA was processed for cDNA synthesis using the Ovation 3’-DGE system (NuGEN Technologies, Inc., San Carlos, CA, USA). The 3’-DGE libraries were constructed from the resulting double stranded cDNA using the TruSeq™ DNA library preparation kit according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). In brief, 500 ng cDNA was sheared and pooled with 500 ng of un-sheared cDNA and end-repaired by enzymatic polishing with T4 DNA polymerase and E.coli DNA polymerase 1 Klenow fragment. A single ‘A’ base was added (A tailing) to the 3’end of the repaired fragments. Illumina pair-end adaptors, essentially consisting of the sequencing primer-annealing sequences, were then ligated to the A-tailed fragments via a 3’ thymine overhang followed by purification using Agentcourt AmpureXP beads (Beckman Coulter, Inc.). The purified ligated DNA was subjected to 11 cycles of PCR amplification to enrich the adaptor modified cDNA libraries using primers complementary to the ends of the adaptors [PCR primers PE 1.0 and PCR primer PE 2.0 (Illumina)]. The purified ligated population was resolved on a 2% (w/v) TAE-agarose gel and 300–500 bp fragments were excised and purified using Qiaquick gel extraction kit. The amplified libraries were quantified using the Agilent DNA high-sensitivity kit on an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara CA, USA). Based on the calculated values, the barcoded libraries were pooled in equal concentration for cluster generation and sequencing onto two lanes of Illumina GAIIx Genome Analyzer (Illumina, San Diego, CA, USA). Cluster formation, primer hybridization and sequencing runs were done according to manufacturer’s protocols. We ran one lane for single-end 101 cycles and a second lane for paired-end 2x101 cycles. Image analysis and base calling were performed using the Illumina Pipeline, where sequence tags were obtained after purity filtering.
For each treatment, there were a total of three replicates or 6 sequencing data files in total. A technical replication was performed on one control sample and one heat shock sample. Technical replicates were merged as one replicate for further analysis.
Raw sequencing reads were initially processed with GenomeQuest tools . Ambiguous residues were trimmed off from both sides of the sequence. Bases with Phred quality below 12 from the 3’ end of the sequence were removed. Reads that were shorter than 40 bases or that contained more than 10 bases with quality below 12 were discarded. Also excluded were reads consisting of repetitive single bases for more than 60% of the length at the 3’ end. A total of 70-84% of reads were retained after the cleanup step.
For mapping of reads, B. taurus genomic sequences bosTau6 (repeat masked) were downloaded from the UCSC genome browser . Sequencing reads of each sample were mapped independently to the reference sequences using Tophat 1.4.1 (University of Maryland) . Tophat splits reads to segments and joins segment alignments. A maximum of one mismatch in each of the 25 bp segments was allowed. The unmapped reads were collected and mapped to the reference using Bowtie 0.12.8  allowing two mismatches. In the third step, the unmapped reads were further mapped to the mRNA sequences using novoalign 2.08.02  allowing for three mismatches for each read. The mRNA sequences of B. taurus were downloaded from Ensemble . A total of 22,915 genes were obtained.
The number of mapped reads for each individual gene was counted by using HTSeq tool  with intersection-nonempty mode. HTSeq takes two input files - the bam or sam files of mapped reads and a gene model file. Ensemble gene annotation file in GTF format was downloaded from Ensemble genome browser. Three criteria were used to identify a gene as being regulated by heat shock: average count of > 2 in at least one treatment, treatment effect of P<0.05 (paired Student t-test), and a fold change (heat shock/control) of > 1.5 or < 0.67.
Analysis of gene networks and functions
Differentially expressed genes were analyzed for functional and regulatory relationships using Ingenuity Pathways Analysis (IPA; Build 172788; release date Aug 11, 2012). The reference set was the Ingenuity Knowledge Base (genes only) and both direct and indirect relationships that were experimentally observed were included. A total of 151 of the 173 differentially expressed genes was present in the reference database. A P value of 0.05 or less was considered significant when determining pathways and functions in which differentially expressed genes were overexpressed.
The list of differentially expressed genes was also annotated by the Database for Annotation, Visualization and Integrated Discovery (DAVID; (DAVID Bioinformatics Resources 6.7, National Institute of Allergy and Infectious Diseases ). A total of 116 genes were included in the analysis of gene functional classification, ontology and pathway analysis.
Temporal changes in expression of selected genes in response to heat shock
At 116 hpi, embryos were either maintained at 38.5°C or cultured at 40°C and 5% CO2 in humidified air. Morula-stage embryos (n=10) were harvested from microdrops at 2, 4, or 8 h (different drops used for each time point). The experiment was replicated eight times so that there were eight pools of 10 embryos each for each time-treatment combination.
RNA was extracted as follow. The zona pellucida was removed as described previously. Embryos were washed with DPBS-PVP and RNA purified as described earlier except the source of DNAse was New England Biolabs (Ipswich, MA, USA).
RNA was reverse transcribed using a commercial kit (High Capacity cDNA Reverse Transcription Kit, Applied Biosystems). Conditions were 25°C for 10 min, 37°C for 120 min and 85°C for 5 min. The cDNA was subjected to real-time PCR to determine expression of HSPA1A, HSP90AA1, SOD1, and CAT using a CFX96 Real-Time System (Bio-Rad, Hercules, CA), the SsoFast Eva-green Supermix reagent (Bio-Rad) and previously-described primers . Briefly, the reaction mixture consisted of 10 μL of SsoFast EvaGreen Supermix, 7.2 μL of DEPC water with template cDNA (2 μL, 0.5 embryo equivalents) and 0.4 μL each of forward and reverse primer (final concentration of each = 0.2 μM). Amplification reactions were performed in replicates of 5 at 95°C for 30 sec followed by 40 cycles of 95°C for 5 sec and 60°C for 10 sec. A negative control for each sample was included that consisted of amplification without reverse transcriptase. Ct values of each gene were calculated by Bio-Rad CFX Manager ver. 1.6 (Bio-Rad). Gene expression values were normalized by the Ct value of the reference gene GAPDH. Fold changes in expression were calculated relative to values for embryos cultured at 38.5°C for 2 h using the ΔΔCt method.
Gene expression data were analyzed by least-squares analysis of variance using the GLM procedure of the Statistical Analysis System (version 9.2 TA Level 2M2, SAS Institute, Cary, NC, USA). Replicate was considered a random effect while temperature and time were considered as fixed. Tests of significance were performed using appropriate error terms. Data analyzed were ΔΔCt values and results are graphed as fold-change differences. Results are shown as least squares means ± SEM.
Disruption in development of morulae by heat shock
At 116 hpi, microdrops of embryos were either maintained at 38.5°C or were moved to a new incubator and exposed to heat shock at 40°C and 5% CO2 in humidified air for 8 h. Afterwards, heat-shocked embryos were returned to 38.5°C and 5% CO2 in humidified air and cultured through day 8 (day 0 = insemination) when blastocyst formation was determined and blastocysts were harvested for determination of numbers of inner cell mass (ICM) and trophectoderm (TE) cells by differential staining . The experiment was replicated five times. Data were analyzed by least-squares analysis of variance using the GLM procedure of SAS. Replicate was considered a random effect, temperature was considered as fixed and the temperature x replicate interaction was used as the error term for temperature. Results are shown as least squares means ± SEM.
Changes in global gene expression caused by heat shock
A total of 173 genes were modified by heat shock, with 94 genes upregulated by heat shock and 79 downregulated. A list of the differentially expressed genes is presented in Additional file 1: Table S1.
Using the IPA procedure of Ingenuity, differentially expressed genes were associated significantly with 22 molecular and cellular functions, with the most significant associations being with cell morphology (19 molecules), cellular assembly and organization (28 molecules), cell death and survival (20 molecules), cell-to-cell signaling and interaction (17 molecules) and lipid metabolism (18 genes).
Genes involved in cellular stress and survival
Further analysis of 42 genes related to removal of free radicals (peroxidase genes, glutathione metabolism genes, thioredoxin, genes containing thioredoxin domains or involved in thioredoxin metabolism, glutaredoxin genes, superoxide dismutase genes, and catalase) indicated that none were significantly affected by heat shock. For example, the fold-change due to heat shock was 1.3 for CATA, 1.2 for GPX2, 1.1 for GLRX1, 1.0 for GLRX3, 1.0 for SOD1, 1.0 for SOD2, 1.0 for TXN, and 1.1 for TXN2. Hypoxia-induced factors and mTOR are also regulated by free radicals . There was no significant effect of heat shock on expression of HIF1A (0.6 fold-change) or components of the mTOR complex [for example, MTOR, 0.8-fold; RICTOR, 0.80-fold, and RPTOR (1.8-fold)] but HIF3A was significantly reduced by heat shock (0.3 fold; P=0.0007).
Effects of heat shock on expression of selected genes involved in the heat shock protein response
Expression level, control
Expression level, heat shock
Heat shock 70 kDa protein 1A
Heat shock protein family B (small), member 11
Heat shock 27kDa protein 1
Heat shock 105 kDa/110kDa protein 1
heat shock 70kDa binding protein, cytoplasmic cochaperone 1
DnaJ (Hsp40) homolog, subfamily C, member 28
Heat shock protein 90kDa alpha (cytosolic), class A member 1
heat shock 70kDa protein 8 (hsc70)
Heat shock protein 90kDa alpha (cytosolic), class B member 1
Heat shock 10kDa protein 1 (chaperonin 10)
heat shock 60kDa protein 1 (chaperonin)
heat shock protein, alpha-crystallin-related, B6
Heat shock factor protein 1
DnaJ (Hsp40) homolog, subfamily B, member 12
DnaJ (Hsp40) homolog, subfamily C , member 25
Heat shock factor protein 2
Genes important for embryonic development
There were no canonical pathways classified by IPA as being involved in organismal growth and development or classified as involved in cellular growth, proliferation, or development that contained significant overrepresentation of differentially expressed genes.
In addition to pathways already described under cell stress, there were five cell signaling pathways containing differentially expressed genes identified by IPA (Figure 2). The most significant pathway was for aldosterone signaling. Using DAVID, the WNT signaling pathway was also identified as a pathway containing an overrepresentation of differentially expressed genes. In particular, heat shock increased expression of AXIN1 and LEF1 and decreased expression of CUL1 and PPP2A.
An additional eight canonical pathways categorized as metabolic pathways by IPA contained significant representation of differentially expressed genes (Figure 2). Three of these pathways involve synthesis of mucopolysaccharides (heparan sulfate, dematan sulfate, and chondroitin sulfate). Another pathway, CMP-N-acetylneuraminate biosynthesis I, contained a gene downregulated by heat shock, GNE, that is involved in synthesis of a sialic acid precursor. Two pathways involved in synthesis of phosphatidylinositol 4,5-bisphosphate (D-myo-inositol (1,4,5)-trisphosphate biosynthesis and 3-phosphoinositide biosynthesis) contained two genes upregulated by heat shock (PIP5K1C, PI4KA). Glutamate removal from folates was represented because GGH was upregulated by heat shock. Salvage pathway of pyrimidine deoxyribonucleotides was represented because heat shock upregulated expression of AICDA, which participates in RNA editing by causing deamination of cytidine.
PUBMED was also queried to identify differentially regulated genes that have been implicated in preimplantation development. Only two such genes were identified: LEF1, involved in TE differentiation , which was increased by heat shock 1.6 fold, and CUL1, which is necessary for development in the mouse  and was decreased by heat shock (0.58 fold).
Putative transcription factors involved in heat-shock induced changes in gene expression
Temporal response of gene expression to heat shock as determined by real-time PCR
Development of heat-shocked embryos to the blastocyst stage
Effect of exposure of embryos at 116 h after insemination to 40°C for 8 h on percent of inseminated oocytes that cleaved and that developed to the blastocyst stage and on blastocyst cell number [inner cell mass (ICM), trophectoderm (TE) and total]
Number of oocytes (replicates)
Percent developed to blastocyst
ICM cell number
TE cell number
74.0 ± 3.5
30.4 ± 1.1
28.3 ± 1.4
75.2 ± 4.1
103.5 ± 5.0
74.5 ± 5.3
28.8 ± 5.3
28.1 ± 1.1
79.3 ± 3.9
107.4 ± 4.5
Within the period of a few days, the newly-formed bovine embryo is transformed from an organism that is unable to successfully survive exposure to elevated temperature to one that is capable of continuing development after exposure to physiological heat shock. Thermal tolerance becomes acquired at about the 8–16 cell stage [2–6], which is coincident with the period of development when the embryonic genome becomes fully activated . Here we hypothesized that the morula-stage embryo is resistant to heat shock because it can undergo transcriptional changes in genes related to cellular survival that stabilize cellular function at elevated temperature.
In fact, however, the transcriptional response to heat shock was muted. Only a few genes involved in the heat shock protein response or antioxidant function were increased by heat shock. Thus, it is likely that heat shock at this stage of development does not cause a large increase in the signals leading to transcription of heat shock protein genes (denatured protein, ref. ) or antioxidant genes (free radicals, ref [40–42]). Put differently, there was not a large-scale transcriptional response of heat shock protein genes or antioxidant genes in response to heat shock because the increased need for the proteins encoded by these genes is minimal. Embryonic resistance of the morula to physiological heat shock as compared to embryos at earlier, more thermosensitive stages is more likely to be due to possession of mechanisms to prevent accumulation of denatured proteins and free radical damage rather than due to differences in transcriptional regulation of heat shock protein and antioxidant genes. Indeed, the two-cell embryo, which is very sensitive to heat shock is capable of increasing transcription of HSPA1A and synthesis of HSP70 [44, 45] in response to heat shock. In addition, the steady-state amount of mRNA for HSPA1A in the two-cell embryo not exposed to heat shock is higher than amounts in the Day 5 embryo after heat shock .
Heat shock did induce expression of some heat shock protein and antioxidant genes. As determined in the 3’DGE analysis, heat shock increased expression of HSPB11 significantly and tended to increase expression of HSPA1A and HSPB1. Also, four genes regulated by NRF2, a transcription factor activated by oxidative stress , were increased by heat shock (AKR7A2, CBR1, GSTA4, and MAP2K5) as was another gene involved in glutathione metabolism (GGH). The experiment using real-time PCR was performed to determine if effects of heat shock on expression of selected heat shock protein and antioxidant genes depended on the duration of heat shock. Results indicated that responses to heat shock were similar after 2, 4 and 8 h of exposure to elevated temperature, Moreover, the apparent increase in expression of HSPA1A detected by 3’DGE was repeatable and significant, as would be expected from previous papers [2, 43], and there was a tendency for expression of HSP90AA1 to increase slightly after heat shock. The increase was small and of a magnitude similar to that seen in the 3’DGE experiment (although the difference was non-significant). In an earlier experiment, there was no increase in steady-state amounts of HSP90AA1 after 24 h of heat shock in one-cell or Day 5 embryos . There was also no effect of heat shock of any duration on expression of CAT or SOD1, in agreement with the results of 3’DGE and an earlier experiment in which embryos were exposed to 24 h of heat shock .
However, heat shock did not alter expression of 64 other genes that encoded heat shock proteins, heat shock transcription factors or heat shock protein binding proteins or expression of 42 other genes involved in antioxidant defense. Analysis of the literature and the results from the present study provide clues as to why the induction of expression of heat shock protein and antioxidant genes was so muted. The lack of a large-scale heat shock protein response could reflect efficient clearance of denatured protein. The most striking change in response to heat shock in the current study was a change in expression of genes encoding for proteins that bind to UBC. This protein is part of the ubiquitination pathway that is critical for removal of damaged protein including in cells damaged by heat shock . Changes in expression of genes involved in the ubiquitin pathway could be reflective of perturbation of the ubiquitin-proteasome pathway by heat shock. Removal of damaged proteins by this pathway could conceivably reduce the signal for activation of HSF1, the transcription factor controlling expression of heat shock protein genes, because this protein is activated as a result of increased accumulation of denatured protein . The idea that embryos that have developed thermotolerance are less likely to accumulate denatured proteins than embryos at earlier stages is consistent with earlier results in mouse embryos, where a more severe heat shock was required to induce HSP70 synthesis in blastocysts than in eight-cell embryos .
The lack of broad increase in genes involved in antioxidant defense probably reflects the inhibition of free radical accumulation in response to heat shock for embryos at this stage of development. Previous results indicate that heat shock did not increase free radical production at Day 4 or 6 after fertilization (16-cell through morula stage of development) although it did when embryos were at Days 0–2 of development (one-cell to early eight-cell stage) . The most likely cause is increased accumulation of intracellular antioxidants. Indeed, one such antioxidant, glutathione, increases in amount by the 9–16 cell stage .
One mechanism contributing to the thermotolerance of the preimplantation embryo is the acquisition of the capacity for apoptosis at the 8–16 cell stage . At physiological temperatures, heat shock causes apoptosis in about 10-20% of the blastomeres [49–51]. Inhibition of apoptosis through administration of caspase 3 inhibitors exacerbates the effects of heat shock on development [51, 52]. In the present study, effects of heat shock on genes involved in apoptosis were inconsistent. Three genes that are associated with activation of apoptosis were increased by heat shock: BOK, NGFR and STK4, which activated apoptosis through a p53-dependent mechanism . However, another gene that inhibits apoptosis, SERINC3 was increased by heat shock and the pro-apoptotic gene, PDIA3, was decreased by heat shock. The possible involvement of the ubiquitin system in clearance of proteins damaged by heat shock, as indicated by the large number of genes regulated by heat shock that bind UBC, may also have implications for apoptosis since ubiquitination plays a role in activation of apoptosis responses .
There were few genes involved in developmental processes whose expression was affected by heat shock. This is not surprising because, as seen previously [2, 5, 6], the development of the morula-stage embryo to the blastocyst stage was unaffected by heat shock. It is likely that adjustments to cellular function that limited the heat shock protein and antioxidant gene response to heat shock also allowed development to continue normally after exposure to 40°C. It may be, however, that some changes in gene expression induced by heat shock compromise ability of the blastocyst to continue development in later embryonic or fetal life. In one study in which Day 4 embryos were exposed to a more severe heat shock than used here (41°C for 12 h), pregnancy rate of the blastocysts formed after heat shock was lower after transfer into recipient females than for blastocysts formed from control embryos . Among the pathways affected by heat shock that could potentially compromise subsequent development was the WNT signaling pathway, where heat shock increased expression of AXIN1 and LEF1 and decreased expression of CUL1 and PPP2A. WNT signaling regulates cell proliferation, cell fate decision and stem cell maintenance . LEF1 is involved in trophectoderm differentiation , and CUL1 is required for development in the mouse . Heat shock also caused changes in expression of genes regulating mucopolysaccharide synthesis. One mucopolysaccharide, hyaluronan, has been reported to increase competence of morula-stage and early blastocyst-stage bovine embryos to establish pregnancy after transfer into recipients . Heparan sulfate proteoglycans also regulate WNT availability by binding to secreted WNTs .
In conclusion, changes in the transcriptome of the heat-shocked bovine morula indicate that the embryo is largely resistant to effects of heat shock. As a result, activation of transcription of genes involved in thermal protection is muted and there is little indication of disruption in gene networks involved in embryonic development. It is likely that the increased resistance of morula-stage embryos to heat shock as compared to embryos at earlier stages of development is due in part to development of mechanisms to prevent accumulation of denatured proteins and free radical damage.
3’ tag digital gene expression
Database for Annotation, Visualization and Integrated Discovery
Dulbecco’s phosphate-buffered salaine
Inner cell mass
Ingenuity Pathways Analysis
Polymerase chain reaction
The authors thank William Rembert, Silvia Carambula, Jim Moss, Sarah Fields, Justin Fear, and Ricardo Nantes for technical help, Central Beef Packing Co. (Center Hill, FL), for donation of ovaries, Scott A. Randell of Southeastern Semen Services (Wellborn, FL), for donation of semen, and Maureen Long, University of Florida, for assistance with IPA. This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2010-85122-20623 from the USDA National Institute of Food and Agriculture and by the National Agriculture Food Research Organization of Japan.
- Gwazdauskas FC, McCaffrey C, McEvoy TG, Sreenan JM: In vitro preimplantation mouse embryo development with incubation temperatures of 37 and 39oC. J Assist Reprod Genet. 1992, 9: 149-154. 10.1007/BF01203755.View ArticlePubMedGoogle Scholar
- Sakatani M, Alvarez NV, Takahashi M, Hansen PJ: Consequences of physiological heat shock beginning at the zygote stage on embryonic development and expression of stress response genes in cattle. J Dairy Sci. 2012, 95: 3080-3191. 10.3168/jds.2011-4986.View ArticlePubMedGoogle Scholar
- Edwards JL, Hansen PJ: Differential responses of bovine oocytes and preimplantation embryos to heat shock. Mol Reprod Dev. 1997, 46: 138-145. 10.1002/(SICI)1098-2795(199702)46:2<138::AID-MRD4>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Ju JC, Parks JE, Yang X: Thermotolerance of IVM-derived bovine oocytes and embryos after short-term heat shock. Mol Reprod Dev. 1999, 53: 336-340. 10.1002/(SICI)1098-2795(199907)53:3<336::AID-MRD9>3.0.CO;2-M.View ArticlePubMedGoogle Scholar
- Sakatani M, Kobayashi S, Takahashi M: Effects of heat shock on in vitro development and intracellular oxidative state of bovine preimplantation embryos. Mol Reprod Dev. 2004, 67: 77-82. 10.1002/mrd.20014.View ArticlePubMedGoogle Scholar
- Bonilla AQS, Oliveira LJ, Ozawa M, Newsom EM, Lucy MC, Hansen PJ: Developmental changes in thermoprotective actions of insulin-like growth factor-1 on the preimplantation bovine embryo. Mol Cell Endocrinol. 2011, 332: 170-179. 10.1016/j.mce.2010.10.009.View ArticlePubMedGoogle Scholar
- Dutt RH: Detrimental effects of high ambient temperature on fertility and early embryo survival in sheep. Int J Biometeorol. 1964, 8: 47-56. 10.1007/BF02186927.View ArticlePubMedGoogle Scholar
- Tompkins EC, Heidenreich CJ, Stob M: Effect of post-breeding thermal stress on embryonic mortality in swine. J Anim Sci. 1967, 26: 377-380.Google Scholar
- Ealy AD, Drost M, Hansen PJ: Developmental changes in embryonic resistance to adverse effects of maternal heat stress in cows. J Dairy Sci. 1993, 76: 2899-2905. 10.3168/jds.S0022-0302(93)77629-8.View ArticlePubMedGoogle Scholar
- Aréchiga CF, Ealy AD, Hansen PJ: Evidence that glutathione is involved in thermotolerance of preimplantation murine embryos. Biol Reprod. 1995, 52: 1296-1301. 10.1095/biolreprod52.6.1296.View ArticlePubMedGoogle Scholar
- Ealy AD, Drost M, Barros CM, Hansen PJ: Thermoprotection of preimplantation bovine embryos from heat shock by glutathione and taurine. Cell Biol Int Rept. 1992, 16: 125-131. 10.1016/S0309-1651(06)80106-2.View ArticleGoogle Scholar
- Sakatani M, Suda I, Oki T, Kobayashi S, Kobayashi S, Takahashi M: Effects of purple sweet potato anthocyanins on development and intracellular redox status of bovine preimplantation embryos exposed to heat shock. J Reprod Dev. 2007, 53: 605-614. 10.1262/jrd.18124.View ArticlePubMedGoogle Scholar
- Sakatani M, Yamanaka K, Kobayashi S, Takahashi M: Heat shock-derived reactive oxygen species induce embryonic mortality in in vitro early stage bovine embryos. J Reprod Dev. 2008, 54: 496-501. 10.1262/jrd.20017.View ArticlePubMedGoogle Scholar
- Matsuzuka T, Sakamoto N, Ozawa M, Ushitani A, Hirabayashi M, Kanai Y: Alleviation of maternal hyperthermia-induced early embryonic death by administration of melatonin to mice. J Pineal Res. 2005, 39: 217-223. 10.1111/j.1600-079X.2005.00260.x.View ArticlePubMedGoogle Scholar
- Sakamoto N, Ozawa M, Yokotani-Tomita K, Morimoto A, Matsuzuka T, Ijiri D, Hirabayashi M, Ushitani A, Kanai Y: DL-α-tocopherol acetate mitigates maternal hyperthermia-induced pre-implantation embryonic death accompanied by a reduction of physiological oxidative stress in mice. Reproduction. 2008, 135: 489-496. 10.1530/REP-07-0379.View ArticlePubMedGoogle Scholar
- Paula-Lopes FF, Al-Katanani YM, Majewski AC, McDowell LR, Hansen PJ: Manipulation of antioxidant status fails to improve fertility of lactating cows or survival of heat-shocked embryos. J Dairy Sci. 2003, 86: 2343-2351. 10.3168/jds.S0022-0302(03)73827-2.View ArticlePubMedGoogle Scholar
- Paula LA DCe, Hansen PJ: Modification of actions of heat shock on development and apoptosis of cultured preimplantation bovine embryos by oxygen concentration and dithiothreitol. Mol Reprod Dev. 2008, 75: 1338-1350. 10.1002/mrd.20866.View ArticleGoogle Scholar
- Lim JM, Liou SS, Hansel W: Intracytoplasmic glutathione concentration and the role of β-mercaptoethanol in preimplantation development of bovine embryos. Theriogenology. 1996, 46: 429-439. 10.1016/0093-691X(96)00165-3.View ArticlePubMedGoogle Scholar
- Li L, Zheng P, Dean J: Maternal control of early mouse development. Development. 2010, 137: 859-870. 10.1242/dev.039487.PubMed CentralView ArticlePubMedGoogle Scholar
- Memli E, First NL: Zygotic and embryonic gene expression in cow: a review of timing and mechanisms of early gene expression as compared with other species. Zygote. 2000, 8: 87-96. 10.1017/S0967199400000861.View ArticleGoogle Scholar
- Ealy AD, Hansen PJ: Induced thermotolerance during early development of murine and bovine embryos. J Cell Physiol. 1994, 160: 463-468. 10.1002/jcp.1041600309.View ArticlePubMedGoogle Scholar
- Al-Katanani YM, Hansen PJ: Induced thermotolerance in bovine two-cell embryos and the role of heat shock protein 70 in embryonic development. Mol Reprod Dev. 2002, 62: 174-180. 10.1002/mrd.10122.View ArticlePubMedGoogle Scholar
- Lepock JR: How do cells respond to their thermal environment?. Int J Hyperthermia. 2005, 21: 681-687. 10.1080/02656730500307298.View ArticlePubMedGoogle Scholar
- Saha S, Sparks AB, Rago C, Akmaev V, Wang CJ, Vogelstein B, Kinzler KW, Velculescu VE: Using the transcriptome to annotate the genome. Nature Biotech. 2002, 19: 508-512.View ArticleGoogle Scholar
- Morrissy S, Zhao Y, Delaney A, Asano J, Dhalla N, Li I, McDonald H, Pandoh P, Prabhu AL, Tam A, Hirst M, Marra M: Digital gene expression by tag sequencing on the illumina genome analyzer. Curr Protoc Hum Genet. 2010, 11: 1-36.Google Scholar
- Matsumura H, Yoshida K, Luo S, Kimura E, Fujibe T, Albertyn Z, Barrero RA, Krüger DH, Kahl G, Schroth GP, Terauchi R: High-throughput superSAGE for digital gene expression analysis of multiple samples using next generation sequencing. PLoS One. 2010, 5: 12010-12017. 10.1371/journal.pone.0012010.View ArticleGoogle Scholar
- Fields SD, Hansen PJ, Ealy AD: Fibroblast growth factor requirements for in vitro development of bovine embryos. Theriogenology. 2011, 75: 1466-1475. 10.1016/j.theriogenology.2010.12.007.View ArticlePubMedGoogle Scholar
- Zhang J, Chiodini R, Badr A, Zhang G: The impact of next-generation sequencing on genomics. J Genet Genomics. 2011, 38: 95-109. 10.1016/j.jgg.2011.02.003.PubMed CentralView ArticlePubMedGoogle Scholar
- UCSC Genome Bioinformatics. http://hgdownload.cse.ucsc.edu/downloads.html]
- Trapnell C, Pachter L, Salzberg SL: TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25: 1105-1111. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Homer N: A survey of sequence alignment algorithms for next-generation sequencing. Brief Bioinform. 2010, 11: 473-483. 10.1093/bib/bbq015.PubMed CentralView ArticlePubMedGoogle Scholar
- Ensemble ftp. ftp://ftp.ensembl.org/pub/release-67/fasta/bos_taurus/cdna/Bos_taurus.UMD3.1.67.cdna.all.fa.gz,
- HTSeq: Analysing high-throughput sequencing data with Python. http://www-huber.embl.de/users/anders/HTSeq/doc/overview.html#,
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Dewhirst MW: Relationships between cycling hypoxia, HIF-1, angiogenesis and oxidative stress. Radiat Res. 2009, 172: 653-665. 10.1667/RR1926.1.PubMed CentralView ArticlePubMedGoogle Scholar
- He S, Pant D, Schiffmacher A, Meece A, Keefer CL: Lymphoid enhancer factor 1-mediated Wnt signaling promotes the initiation of trophoblast lineage differentiation in mouse embryonic stem cells. Stem Cells. 2008, 26: 842-849. 10.1634/stemcells.2007-0356.View ArticlePubMedGoogle Scholar
- Wang Y, Penfold S, Tang X, Hattori N, Riley P, Harper JW, Cross JC, Tyers M: Deletion of the Cul1 gene in mice causes arrest in early embryogenesis and accumulation of cyclin E. Curr Biol. 1999, 9: 1191-1194. 10.1016/S0960-9822(00)80024-X.View ArticlePubMedGoogle Scholar
- Calderwood SK, Gong J: Molecular chaperones in mammary cancer growth and breast tumor therapy. J Cell Biochem. 2012, 113: 1096-1103. 10.1002/jcb.23461.PubMed CentralView ArticlePubMedGoogle Scholar
- Valko M, Rhodes CJ, Moncol J, Izakovic M, Mazur M: Free radicals, metals and antioxidants in oxidative stress-induced cancer. Chem Biol Interact. 2006, 160: 1-40. 10.1016/j.cbi.2005.12.009.View ArticlePubMedGoogle Scholar
- Shi Q, Gibson GE: Oxidative stress and transcriptional regulation in Alzheimer disease. Alzheimer Dis Assoc Disord. 2007, 21: 276-291. 10.1097/WAD.0b013e31815721c3.PubMed CentralView ArticlePubMedGoogle Scholar
- Nguyen T, Nioi P, Pickett CB: The Nrf2-antioxidant response element signaling pathway and its activation by oxidative stress. J Biol Chem. 2009, 284: 13291-13295. 10.1074/jbc.R900010200.PubMed CentralView ArticlePubMedGoogle Scholar
- Chandolia R, Peltier MR, Tian W, Hansen PJ: Transcriptional control of development, protein synthesis, and heat-induced heat shock protein 70 synthesis in 2-cell bovine embryos. Biol Reprod. 1999, 61: 1644-1648. 10.1095/biolreprod61.6.1644.View ArticlePubMedGoogle Scholar
- Edwards JL, Hansen PJ: Elevated temperature increases heat shock protein 70 synthesis in bovine two-cell embryos and compromises function of maturing oocytes. Biol Reprod. 1996, 55: 340-346.Google Scholar
- Edwards JL, Ealy AD, Monterroso VH, Hansen PJ: Ontogeny of temperature-regulated heat shock protein 70 synthesis in preimplantation bovine embryos. Mol Reprod Dev. 1997, 48: 25-33. 10.1002/(SICI)1098-2795(199709)48:1<25::AID-MRD4>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Riezman H: Why do cells require heat shock proteins to survive heat stress?. Cell Cycle. 2004, 3: 61-63.View ArticlePubMedGoogle Scholar
- Edwards JL, Ealy AD, Hansen PJ: Regulation of heat shock protein 70 synthesis by heat shock in the preimplantation murine embryo. Theriogenology. 1995, 44: 329-337. 10.1016/0093-691X(95)00188-E.View ArticlePubMedGoogle Scholar
- Hansen PJ, Fear JM: Cheating death at the dawn of life: developmental control of apoptotic repression in the preimplantation embryo. Biochem Biophys Res Commun. 2011, 413: 155-158. 10.1016/j.bbrc.2011.08.064.View ArticlePubMedGoogle Scholar
- Paula-Lopes FF, Hansen PJ: Heat shock-induced apoptosis in preimplantation bovine embryos is a developmentally regulated phenomenon. Biol Reprod. 2002, 66: 1169-1177.PubMedGoogle Scholar
- Brad AM, Hendricks KE, Hansen PJ: The block to apoptosis in bovine two-cell embryos involves inhibition of caspase-9 activation and caspase-mediated DNA damage. Reproduction. 2007, 134: 789-797. 10.1530/REP-07-0146.View ArticlePubMedGoogle Scholar
- Jousan FD, Hansen PJ: Insulin-like growth factor-I promotes resistance of bovine preimplantation embryos to heat shock through actions independent of its anti-apoptotic actions requiring PI3K signaling. Mol Reprod Dev. 2007, 74: 189-196. 10.1002/mrd.20527.View ArticlePubMedGoogle Scholar
- Paula-Lopes FF, Hansen PJ: Apoptosis is an adaptive response in bovine preimplantation embryos that facilitates survival after heat shock. Biochem Biophys Res Commun. 2002, 295: 37-42. 10.1016/S0006-291X(02)00619-8.View ArticlePubMedGoogle Scholar
- Rodriguez JM, Glozak MA, Ma Y, Cress WD: Bok, Bcl-2 related ovarian killer, is cell cycle-regulated and sensitizes to stress-induced apoptosis. J Biol Chem. 2006, 281: 22729-22735. 10.1074/jbc.M604705200.PubMed CentralView ArticlePubMedGoogle Scholar
- Mukai J, Hachiya T, Shoji-Hoshino S, Kimura MT, Nadano D, Suvanto P, Hanaoka T, Li Y, Irie S, Greene LA, Sato T: NADE, a p75NTR-associated cell death executor, is involved in signal transduction mediated by the common neurotrophin receptor p75NTR. J Biol Chem. 2000, 275: 17566-17570. 10.1074/jbc.C000140200.View ArticlePubMedGoogle Scholar
- Yuan F, Xie Q, Wu J, Bai Y, Mao B, Dong Y, Bi W, Ji G, Tao W, Wang Y, Yuan Z: MST1 promotes apoptosis through regulating Sirt1-dependent p53 deacetylation. J Biol Chem. 2011, 286: 6940-6945. 10.1074/jbc.M110.182543.PubMed CentralView ArticlePubMedGoogle Scholar
- Bossolasco M, Veillette F, Bertrand R, Mes-Masson AM: Human TDE1, a TDE1/TMS family member, inhibits apoptosis in vitro and stimulates in vivo tumorigenesis. Oncogene. 2006, 25: 4549-4558. 10.1038/sj.onc.1209488.View ArticlePubMedGoogle Scholar
- Pressinotti NC, Klocker H, Schäfer G, Luu VD, Ruschhaupt M, Kuner R, Steiner E, Poustka A, Bartsch G, Sültmann H: Differential expression of apoptotic genes PDIA3 and MAP3K5 distinguishes between low- and high-risk prostate cancer. Mol Cancer. 2009, 8: 130-10.1186/1476-4598-8-130.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Y, Yu X: Regulation of apoptosis: the ubiquitous way. FASEB J. 2008, 17: 790-799.View ArticleGoogle Scholar
- Barros CM, Pegorer MF, Vasconcelos JL, Eberhardt BG, Monteiro FM: Importance of sperm genotype (indicus versus taurus) for fertility and embryonic development at elevated temperatures. Theriogenology. 2006, 65: 210-218. 10.1016/j.theriogenology.2005.09.024.View ArticlePubMedGoogle Scholar
- Van Amerongen R, Nusse R: Towards an integrated view of Wnt signaling in development. Development. 2009, 136: 3205-3214. 10.1242/dev.033910.View ArticlePubMedGoogle Scholar
- Block J, Bonilla L, Hansen PJ: Effect of addition of hyaluronan to embryo culture medium on survival of bovine embryos in vitro following vitrification and establishment of pregnancy after transfer to recipients. Theriogenology. 2009, 71: 1063-1071. 10.1016/j.theriogenology.2008.11.007.View ArticlePubMedGoogle Scholar
- Fuerer C, Habib SJ, Nusse R: A study on the interactions between heparin sulfate proteoglycans and Wnt proteins. Dev Dynam. 2010, 239: 184-190.Google Scholar
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.