Changes in the transcriptome of morula-stage bovine embryos caused by heat shock: relationship to developmental acquisition of thermotolerance
Reproductive Biology and Endocrinology volume 11, Article number: 3 (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
Gene networks and gene relationships identified using IPA were initially scrutinized for effects of heat shock on genes and gene networks involved in cell death and survival. A total of 38 differentially-regulated genes (25% of the differentially expressed genes available for analysis) were associated with the ubiquitin protein, UBC (Figure 1). These included molecules that bind to UBC (11 genes upregulated by heat shock and 18 downregulated by heat shock) as well as molecules that interact with molecules that bind UBC (6 upregulated and 3 downregulated genes). There were also three canonical pathways related to cell stress and survival identified as having significant representation of differentially expressed genes (Figure 2). These were glutathione-mediated detoxification (two genes upregulated by heat shock), protein ubiquitination pathway (3 upregulated and 2 downregulated) and NRF2-mediated oxidative stress response (4 upregulated genes). Using DAVID, the positive regulation of apoptosis ontology was also found to contain overrepresentation of differentially expressed genes, with BOK, NGFR, STK4, and SERINC3 increased by heat shock and PDIA3 decreased by heat shock.
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).
Expression of heat shock protein genes was not generally affected by heat shock. We looked at expression of 68 genes involved in heat shock protein responses, including heat shock proteins, heat shock transcription factors and heat shock protein binding proteins. Results from a selection of these genes are shown in Table 1. Heat shock increased expression of HSPB11 (P<0.05) and tended to increase expression of HSPA1A (P=0.06) and HSPB1 (P=0.09). Heat shock also increased (P<0.05) expression of the heat shock binding protein gene, HSPBP1. None of the other heat shock protein genes, heat shock transcription factors, or heat shock protein binding proteins was affected by heat shock.
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
Ingenuity was used to identify transcription factors that regulate genes that were regulated in response to heat shock. A total of 14 transcription factors were significant. Together, these control expression of 15 of the differentially-regulated genes (Figure 3). HSF1 or HSF2, which control expression of heat shock protein genes, and NRF2, which controls regulation of antioxidant genes, were not identified.
Temporal response of gene expression to heat shock as determined by real-time PCR
Results from 3’DGE indicated that heat shock caused little change in expression of genes for heat shock proteins or antioxidant enzymes, with a tendency for a slight increase in expression of HSPA1A and no effect on most heat shock protein genes or antioxidant genes. To confirm these observations and to determine whether responses to heat shock would be different after other times of exposure, the effects of heat shock on steady-state amounts of mRNA for HSPA1A, HSP90AA, SOD1 and CAT after 2, 4 or 8 h of exposure were determined (Figure 4). Heat shock increased steady-state amounts of mRNA for HSPA1A (P<0.05) and tended to increase amounts of HSP90AA1 (P<0.07) but had no effect on expression of SOD1 or CAT. As indicated by the lack of a time x temperature interaction, the increase in mRNA for HSPA1A and HSP90AA1 caused by heat shock occurred after 2 h of heat shock and was similar in magnitude at later times examined.
Development of heat-shocked embryos to the blastocyst stage
An experiment was conducted to confirm that exposure of morula-stage embryos to a heat shock of the magnitude used for the gene expression studies (40°C for 8 h) would not compromise development to the blastocyst stage. There was no effect of heat shock on the proportion of inseminated oocytes that became blastocysts or on blastocyst cell number (ICM, TE or total) (Table 2).
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
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.
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.
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.
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.
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.
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.
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.
Tompkins EC, Heidenreich CJ, Stob M: Effect of post-breeding thermal stress on embryonic mortality in swine. J Anim Sci. 1967, 26: 377-380.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Li L, Zheng P, Dean J: Maternal control of early mouse development. Development. 2010, 137: 859-870. 10.1242/dev.039487.
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.
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.
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.
Lepock JR: How do cells respond to their thermal environment?. Int J Hyperthermia. 2005, 21: 681-687. 10.1080/02656730500307298.
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.
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.
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.
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.
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.
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.
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.
Li H, Homer N: A survey of sequence alignment algorithms for next-generation sequencing. Brief Bioinform. 2010, 11: 473-483. 10.1093/bib/bbq015.
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.
Dewhirst MW: Relationships between cycling hypoxia, HIF-1, angiogenesis and oxidative stress. Radiat Res. 2009, 172: 653-665. 10.1667/RR1926.1.
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.
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.
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.
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.
Shi Q, Gibson GE: Oxidative stress and transcriptional regulation in Alzheimer disease. Alzheimer Dis Assoc Disord. 2007, 21: 276-291. 10.1097/WAD.0b013e31815721c3.
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.
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.
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.
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.
Riezman H: Why do cells require heat shock proteins to survive heat stress?. Cell Cycle. 2004, 3: 61-63.
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.
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.
Paula-Lopes FF, Hansen PJ: Heat shock-induced apoptosis in preimplantation bovine embryos is a developmentally regulated phenomenon. Biol Reprod. 2002, 66: 1169-1177.
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.
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.
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.
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.
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.
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.
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.
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.
Yang Y, Yu X: Regulation of apoptosis: the ubiquitous way. FASEB J. 2008, 17: 790-799.
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.
Van Amerongen R, Nusse R: Towards an integrated view of Wnt signaling in development. Development. 2009, 136: 3205-3214. 10.1242/dev.033910.
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.
Fuerer C, Habib SJ, Nusse R: A study on the interactions between heparin sulfate proteoglycans and Wnt proteins. Dev Dynam. 2010, 239: 184-190.
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.
The authors state that they have no competing interests.
Conceived and designed experiments: MS, PJH. Carried out experiments: MS, KBD, JB, MO, SS. Analyzed data: MS, JY. Wrote the initial drafts of the paper: MS, SS, JY, PJH. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Genes whose expression was differentially expressed by heat shock in bovine morula. Shown (from left to right) are the ENSEMBL gene identification number, the number of counts for each of three control samples (C5, C6, C14) and three heat-shock samples (H5, H6, H14), the average count, the mean counts for control (MeanC) and heat shock (MeanH), the fold change (heat shock/control), the standard error of the difference, the p value, gene name and gene description. Genes that are upregulated by heat shock are shaded light red while those that are downregulated are shaded green. (XLS 80 KB)
About this article
Cite this article
Sakatani, M., Bonilla, L., Dobbs, K.B. et al. Changes in the transcriptome of morula-stage bovine embryos caused by heat shock: relationship to developmental acquisition of thermotolerance. Reprod Biol Endocrinol 11, 3 (2013). https://doi.org/10.1186/1477-7827-11-3