Skip to main content

Impact of heat shock transcription factor 1 on global gene expression profiles in cells which induce either cytoprotective or pro-apoptotic response following hyperthermia



Elevated temperatures induce activation of the heat shock transcription factor 1 (HSF1) which in somatic cells leads to heat shock proteins synthesis and cytoprotection. However, in the male germ cells (spermatocytes) caspase-3 dependent apoptosis is induced upon HSF1 activation and spermatogenic cells are actively eliminated.


To elucidate a mechanism of such diverse HSF1 activity we carried out genome-wide transcriptional analysis in control and heat-shocked cells, either spermatocytes or hepatocytes. Additionally, to identify direct molecular targets of active HSF1 we used chromatin immunoprecipitation assay (ChIP) combined with promoter microarrays (ChIP on chip). Genes that are differently regulated after HSF1 binding during hyperthermia in both types of cells have been identified. Despite HSF1 binding to promoter sequences in both types of cells, strong up-regulation of Hsps and other genes typically activated by the heat shock was observed only in hepatocytes. In spermatocytes HSF1 binding correlates with transcriptional repression on a large scale. HSF1-bound and negatively regulated genes encode mainly for proteins required for cell division, involved in RNA processing and piRNA biogenesis.


Observed suppression of the transcription could lead to genomic instability caused by meiotic recombination disturbances, which in turn might induce apoptosis of spermatogenic cells. We propose that HSF1-dependent induction of cell death is caused by the simultaneous repression of many genes required for spermatogenesis, which guarantees the elimination of cells damaged during heat shock. Such activity of HSF1 prevents transmission of damaged genetic material to the next generation.


Somatic cells protect themselves from damage during cellular stress inducing so called heat shock response, which results in heat shock proteins (HSPs) synthesis. There are many different HSPs classified to families based on molecular weight and sequence homology. In mammalian cells there are five major classes of HSP families, namely the HSPH (HSP110), HSPC (HSP90), HSPA (HSP70), DNAJ (HSP40) and HSPB (small HSP exemplified by HSP27), and two chaperonin families: HSPD/E (HSP60/HSP10) and CCT (TRiC) [1]. Generally, HSPs prevent inappropriate protein aggregation and mediate transport of immature proteins to the target organelles for final packaging, degradation or repair. Although some members of HSPs are expressed constitutively in the absence of stress or during gametogenesis and embryogenesis [2, 3], many of them accumulate to a high level in cells subjected to different types of stress stimuli and remain elevated for a prolonged period [4]. Heat shock factor 1 (HSF1) is the primary transcription factor responsible for the activation of Hsp genes following stress. In physiological conditions HSF1 exists as an inactive monomer. Activation of HSF1 in response to cellular stress is connected with its trimerization, phosphorylation and binding to DNA in the promoter regions containing the heat shock elements (HSEs), which are present mainly in heat shock genes [5]. In addition to the regulation of Hsp genes, HSF1 is involved in the transcription of numerous other genes, both in the absence or presence of heat shock. In Saccharomyces or Drosophila the direct transcriptional targets of HSF represent nearly 3% of genes [6, 7]. These genes encode for proteins involved in diverse cellular processes such as RNA splicing, apoptosis, ubiquitinylation and protein degradation, detoxification, energy generation, carbohydrate metabolism, small molecule transport, cell signaling and maintenance of cell integrity [68].

Despite the high degree of conservation of the heat shock response, different cells vary in their ability to induce HSPs synthesis and consequently in sensitivity to damaging agents. HSPs overexpression in various human cancers diminishes the success of anti-cancer treatment by increasing the resistance of cancer cells to therapy [9]. On the other hand, some neurons, pre-ovulatory oocytes, spermatocytes and some stages of embryonic development [1013], as well as certain tumor cell lines (especially of myeloid origin, e.g. lymphomas; [14]), are hypersensitive to elevated temperatures. It has been shown, at least for spermatocytes, that inducible HSP70 expression is blocked in such cells [15, 16]. Opposite to most somatic cells, in which HSF1 is a part of the cytoprotective system, in spermatocytes it acts as a pro-apoptotic factor [17, 18]. Moreover, the testis-specific variant of HSP70 is depleted in cells undergoing HSF1-induced apoptosis [19]. Activation of HSF1 in male germ cells induces massive degeneration of a seminiferous epithelium, which leads to male infertility [1618]. In fact, primary spermatocytes are germ cells the most sensitive to heat stress [13]. Spermatocytes are very unique cells. They originate from spermatogonia and divide by meiosis giving haploid spermatids (that finally differentiate to spermatozoa) [20]. During the whole process of spermatogenesis dramatic changes in patterns of gene expression and chromatin structure are observed. In particular, the first meiotic division (occurring in primary spermatocytes) involves many cell-specific gene products. They are needed for correct processing of chromosome condensation, pairing of homologous chromosomes, formation of the synaptonemal complexes, and genetic recombination. These unique processes should be highly orchestrated and any disturbances at that stage of spermatogenesis could lead to fertility problems [21].

In the majority of mammals, the male gonads are located outside the main body cavity to provide the lower testicular temperature required for correct spermatogenesis and fertility. Increasing of the testis temperature up to the body temperature (or above it) leads to the activation of HSF1 [22]. Active HSF1 acts as a cell-survival factor only in pre-meiotic germ cells, but not in meiotic and post-meiotic ones [15]. Both mitochondria-dependent and death receptor-dependent pathways appear to be involved in the HSF1-induced apoptosis of spermatogenic cells: the levels of BCL-2 family proteins is increased, p53 protein accumulates and expression levels of caspase-8 and death-receptor-interacting proteins (including FADD and TRADD) are elevated [16]. The apoptosis of spermatocytes induced by HSF1 can be an important mechanism involved in the removal of aberrant germ cells. This can prevent maturation of damaged male germ cells and transmission of incorrect genetic information to the next generation.

Transcriptional activity of HSF1 is apparently indispensable in its pro-apoptotic functions within spermatocytes. However, details of molecular mechanisms involved in activation of HSF1-triggered apoptosis of male germ cells are not known at present. We performed studies that aimed to disclose the initiating events leading to apoptosis of male germ cells after temperature elevation. For this purpose, the genome-wide transcriptional analysis was performed both in control and heat-shocked cells, either isolated mouse spermatocytes or hepatocytes representing somatic cells, using the Affymetrix GeneChip system. Genes that are differently expressed after hyperthermia in both types of cells have been identified. To find out genes directly regulated by HSF1, the chromatin immunoprecipitation assay combined with DNA microarray (ChIP on chip) was performed. This approach enabled the identification of genes targeted by this transcription factor, in either somatic or male germ cells.


Differences in transcriptional response to hyperthermia between somatic and spermatogenic cells

The physiological temperature for mouse testes (normally located outside the body cavity) is approximately 32-33°C. Increasing the testicular temperature up to 37°C and above does not induce cytoprotective mechanisms, as it is observed in somatic cells during hyperthermia, but leads to the degeneration of seminiferous epithelium [13]. Aiming to decipher the differences in molecular mechanisms induced by heat shock in mouse somatic and spermatogenic cells we applied the strategy which is schematically illustrated in Figure 1. In the first step, we performed the global gene expression profiling using Affymetrix microarrays in control and heat shocked cells. As model somatic cells we used hepatocytes that respond to hyperthermia in a classical way by induction of heat shock genes transcription (Figure 2). As spermatogenic cells we used a fraction of cells enriched with spermatocytes (isolated by unit gravity sedimentation in bovine serum albumin, BSA, gradient), which are the most sensitive to damage at elevated temperatures [13]. Using isolated spermatocytes we avoided the influence of the somatic testicular component in the final results. We took advantage of Hsp genes as reference genes to monitor response to heat shock, and observed strong induction of Hspa1, Hspa8 and Hsph1 genes transcription two to six hours after hyperthermia in the liver and in whole testes, but not in the fraction of isolated spermatocytes (Figure 2). This result is consistent with previous findings that mRNAs of major heat shock genes are hardly induced in vivo in response to heat shock in male germ cells [15, 23]. On Affymetrix gene chip arrays we analyzed RNA isolated from untreated (control) hepatocytes and after 2 h of recovery at the physiological temperature from heat shock performed at 43°C. In the case of spermatocytes, RNA was isolated from untreated cells and after 2 h of recovery at the physiological temperature from heat shock performed at 38°C or 43°C.

Figure 1
figure 1

General scheme of the study.

Figure 2
figure 2

Monitoring of the heat shock response in isolated spermatocytes, testes and liver. Transcripts of selected Hsp genes up to 24 h after hyperthermia were detected by RT-PCR. Loading control reactions were performed with primers specific to 18S rRNA transcript; N, PCR negative control without template.

To exclude from analyses genes with very low levels of expression thresholds of signals (in arbitrary units in log2) were calculated (Additional file 1: Table S1 [24]). Expression levels of 15 364 genes with signals above the noise threshold level registered in at least one experimental point are presented in Additional file 2[24]. Changes in the level of genes expression in heat-shocked (HS) samples versus control, untreated (C), are presented as a signal log ratio (SLR; the average value of HS/average value of C). For each experimental group significant SLR was calculated, and genes affected by hyperthermia in spermatocytes and hepatocytes were found (Table 1). The analysis revealed that global changes in the levels of expression in spermatocytes were smaller (not exceeding 8-fold changes up or down) than in hepatocytes (exceeding 64-fold up). However, in spermatocytes subjected to hyperthermia at 38°C many more genes were affected than in hepatocytes (Table 1). Importantly, while in hepatocytes the number of up- and down-regulated genes was comparable, in spermatocytes significantly more genes were inhibited than up-regulated.

Table 1 Representation of genes changing their expression following heat shock

The ten top genes with the highest induction of the transcription after hyperthermia are shown in Additional file 3: Table S2 [24]. Approximately 30% of genes activated in spermatocytes subjected to heat shock at 38°C (SLR ≥ 0.82) were expressed in hepatocytes at a very low level (below the noise threshold). In turn, almost 40% encode for proteins with an unknown function. Cluster analysis and the Gene to GO BP (Gene Ontology Biological Process) test showed that, among activated genes, those involved in negative regulation of different biological processes are over-represented (Additional file 4; clusters: 87, 79 [24]). Importantly, activation of gene expression was rather marginal in spermatocytes subjected to heat shock at 43°C (Table 1 and Additional file 5[24]). Genes that were induced in hepatocytes following hyperthermia (SLR ≥ 0.75) were usually not up-regulated in spermatocytes, and for 20% of them the level of expression in spermatocytes was below the noise threshold. The noted exception was activation of Hspa1 gene in both hepatocytes and spermatocytes (after heat shock at 43°C), which suggests that the isolated fraction of spermatogenic cells might be slightly contaminated by somatic cells. As expected, the Gene to GO BP test for over-representation revealed that genes grouped in clusters with the highest degree of induction in hepatocytes are involved mainly in processes corresponding to the stress response and protein folding. Strongly activated genes are also associated with inflammatory response, cytokine production and T cell activation (Additional file 6; clusters: 1, 6, 21, 65, 16, 25, 47 [24]).

The ten top genes with the highest inhibition of the transcription following hyperthermia are shown in Additional file 7: Table S3 [24]. Among genes down-regulated in spermatocytes following heat shock at 38°C (SLR ≤ −0.82) a large portion (6.4%) was not expressed in hepatocytes, and uncharacterized genes represented about 10%. Cluster analysis and the Gene to GO BP test showed that among genes down-regulated in spermatocytes at 38°C strongly over-represented are those involved in RNA metabolism, cellular process involved in reproduction, regulation of gene expression and signal transduction (Additional file 4; clusters: 22, 55, 11, 82, 16, 83, 7, 15 [24]). In spermatocytes subjected to hyperthermia at 43°C inhibition of the transcription was more prominent than activation, although it was weaker than in cells treated at 38°C. Down regulation of genes involved in response to heat and response to oxidative stress was the most characteristic (Additional file 5; clusters: 34, 50 [24]). Genes down-regulated in heat shocked hepatocytes (SLR ≤ −0.75) are involved in lipid, steroid and isoprenoid metabolism, catabolic processes and oxidation-reduction processes (Additional file 6; clusters: 48, 30, 97, 6 [24]); 30% of such genes showed very low level of expression in spermatocytes (below the noise threshold). It should be also noted that for genes with the highest repression following the hyperthermia of hepatocytes the big differences in the level of expression between samples were observed (Additional file 7: Table S3, part C [24]), which makes this group less conclusive.

Our expression data indicated that responses of spermatocytes to hyperthermia at 38°C and at 43°C were quantitatively and qualitatively different. It was previously suggested that heat shock at 43°C resulted rather in rapid changes of the protein localization than their levels (e.g. BAX and cytochrome c redistribution was observed in spermatocytes shortly after treatment; [25]). Hence, comparison of changes induced in hepatocytes treated at 43°C with changes induced in spermatocytes treated at 38°C appeared more relevant (in both cases the temperature of the hyperthermia was 5-6°C above physiological).

Both hepatocytes and spermatocytes showed significant differences in the transcriptional response to stress, and only a few genes responded similarly (Additional file 8: Table S4 [24]). The majority of genes affected by hyperthermia were differentially regulated in hepatocytes and spermatocytes. As expected, Hsp genes and some other genes coding for proteins that enable recovery from stress were strongly induced in hepatocytes, but not in spermatocytes. Furthermore, some Hsp genes were strongly down-regulated in spermatocytes, although there were also examples of chaperones’ expression activation (Table 2).

Table 2 The most significant changes in expression of genes from the HSPs family (and some other chaperones)

Expression of several genes encoding proteins involved (or potentially involved) in cell death was differentially regulated in both types of cells (Table 3). Genes coding for proteins with proved anti-apoptotic properties (e.g. Bag3, Phlda1) were strongly activated only in hepatocytes. In spermatocytes some genes involved in regulation of apoptosis were repressed (e.g. inhibitors of apoptosis: Api5 and Syvn1). It is noteworthy that expression of Pyhin1, which possibly promote ubiquitination and subsequent degradation of MDM2 leading to TP53 stabilization [26], was observed only in spermatocytes. This is consistent with the observation that p53 was accumulated in the testes of transgenic mice expressing active HSF1 [16].

Table 3 The most significant changes in expression of genes involved in cell death regulation

During the heat shock response expression of genes coding for some transcriptional regulators involved in the growth, proliferation and differentiation was strongly activated in hepatocytes, but not in spermatocytes (especially Atf3 and genes coding for components of AP1 transcription factor: JUN and FOS; Table 4). In spermatocytes, a different set of transcriptional regulators was activated (which were not activated in hepatocytes, except for Zfp36l1). In spermatocytes numerous transcriptional regulators were repressed, e.g. some components of general transcription factors machinery (Table 4). Some of the transcription factors, whose mRNAs were elevated in hepatocytes during the heat shock response, play a role in the inflammatory responses (e.g. NFKBIZ, NFIL3). Consistently, some other genes involved in the inflammatory and immune responses were significantly up-regulated in hepatocytes (e.g. Socs3 and genes coding for chemokine ligands; Table 5).

Table 4 Selected transcriptional regulators differentially expressed after hyperthermia
Table 5 The most significant changes in expression of genes involved in the inflammatory and immune responses

Data obtained from microarrays analyses were further verified and visualized by semi-quantitative RT-PCR on independent material. Expression of selected genes was analyzed in isolated spermatocytes, and in whole testes and liver, after 2 h, 4 h, 6 h and 24 h of recovery from heat shock. These analyses revealed that the level of transcripts strongly elevated after hyperthermia at 43°C in somatic cells (both in liver and testes) returned to the steady-state level 24 h after treatment (Figure 2 and Additional file 9: Figure S1 [24]). Analyses confirmed also a lack of activation of Hsp genes (and other genes strongly activated in hepatocytes) in isolated spermatocytes up to 24 h after treatment. RT-PCR analyses generally confirmed data obtained from microarrays studies.

Differences in HSF1 binding to DNA between somatic and spermatogenic cells

HSF1 is activated at elevated temperatures very rapidly, and the strongest HSF1 binding to DNA was observed after 2–20 minutes of heat shock treatment, depending on experimental model [27]. To find out the optimal conditions of HSF1 binding in our experimental system we first analyzed by ChIP-PCR the ability of HSF1 to bind to promoters of Hsph1, Hspa1, Dnaja1 and Phlda1 (genes already known to be regulated by HSF1). Binding of HSF1 was analyzed after 5, 10, 20 and 40 minutes of incubation of cells at elevated temperatures. Isolated hepatocytes were incubated at 43°C. Isolated spermatocytes were incubated at 38°C or 43°C. Observed HSF1 binding varied depending on the gene and conditions of heat shock (time and temperature; Figure 3). In spermatocytes during heat shock at 38°C the HSF1 binding was generally much weaker than at 43°C (both in spermatocytes and hepatocytes), and a prolonged period of hyperthermia (40 min) did not straighten HSF1 binding. Interestingly, in the case of Hspa1 genes, the most activated genes during heat shock in somatic cells, we observed only moderate HSF1 binding to the promoters in hepatocytes (and no binding in spermatocytes). In further in-depth ChIP on chip studies our goal was not study of the kinetics of HSF1 binding but finding the broadest spectrum of its targets. Thus chromatin immunoprecipitation was done on cells heat shocked for 5, 10 and 20 minutes, than combined in one sample. ChIP on chip analyses were performed on Affymetrix mouse promoter tiling arrays to find out the molecular targets of HSF1 in all promoter regions. Untreated cells (control) and cells subjected to hyperthermia (hepatocytes at 43°C and spermatocytes at 38°C or 43°C) were analyzed, and the mean signals (in log2 scale) for each probeset on the array have been calculated in samples with antibody (denoted as AB1) and in samples without antibody (denoted as AB0). Binding of HSF1 was assessed as AB1-AB0 value. An example of analysis showing the graphical representation of AB1-AB0 signals from each probeset distributed on promoters is shown in Figure 4. Finally, numerical value estimating HSF1 binding was calculated in arbitrary units for each promoter (mean signal from all probesets annotated to the promoter). Results of such analyzes obtained for each promoter are presented in Additional file 10[24] (containing additional information about expression level of the corresponding transcript).

Figure 3
figure 3

Kinetics of HSF1 activation. ChIP analysis of HSF1 binding to promoters of selected genes in hepatocytes and spermatocytes subjected to heat shock in different conditions (from 5 to 40 minutes, and temperatures 38°C or 43°C; hepatocytes were heat shocked only at 43°C). N, PCR-negative control without template.

Figure 4
figure 4

Graphical representation of HSF1 binding to the promoters of Hsph1 , Dnaja1 and Mov10l1 genes following hyperthermia in hepatocytes and spermatocytes. Results of ChIP on chip analysis are shown. Each line represents the signal from one probeset and is shown with reference to the position on the chromosome (numbered according to NCBIv36 of the mouse genome).

Searching for HSF1 binding at the normal physiological temperature in spermatocytes we found only a few genes with statistically significant signals (false discovery rate, FDR < 0.05), while in hepatocytes there was no such genes (Table 6). ChIP on chip analyzes revealed that following hyperthermia promoters of several Hsp genes and some other genes already known to be regulated by HSF1 [2830], interact with HSF1 in both hepatocytes and spermatocytes (Figure 5). In this group of genes, binding of HSF1 at the physiological temperature was detected only to the Hsph1 and Hsp90ab1 promoters, although the significance of the signal was only marginal (FDR > 0.05). In general, heat shock strongly induced HSF1 binding (yet binding at 38°C in spermatocytes was usually weaker). Interestingly, HSF1 interacted with more genes in spermatocytes than in hepatocytes, although in hepatocytes the binding was stronger (Table 6). More than 50% of genes interacting with HSF1 in hepatocytes, were bound by HSF1 also in spermatocytes (Figure 6A). Interestingly, in hepatocytes (but not in spermatocytes) the uppermost HSF1 binding following heat shock was observed for sequences located on the Chromosome Y (Additional file 11: Table S5 [24]). Genes showing the uppermost HSF1 binding in spermatocytes also showed strong binding in hepatocytes (Additional file 11: Table S5 [24]). The Gene to GO BP test showed that genes bound by HSF1 in both hepatocytes and spermatocytes have the strongest representation in following classes: RNA processing and transport, DNA replication, repair, recombination and chromatin modification, protein folding, translation and ubiquitination, Golgi organization, zinc ion transport and microtubule cytoskeleton. As could be expected, HSF1-bound genes involved in spermatogenesis, meiosis and mitosis (and related to cell cycle and cell division) have much higher representation in spermatocytes than in hepatocytes (Additional file 12[24]).

Figure 5
figure 5

Binding of HSF1 estimated in ChIP on chip analysis to promoters of selected genes. (A) Binding at physiological and (B) elevated temperatures. HSF1 binding is calculated as AB1-AB0 mean signal in arbitrary units. *FDR > 0.05.

Table 6 Statistics of HSF1 binding
Figure 6
figure 6

Overlap of genes affected by hyperthermia in hepatocytes and spermatocytes. Numbers of genes showing binding of HSF1 (A), modulated expression (B), and both (C). For HSF1 binding values AB1-AB0 ≥ 15 and FDR < 0.05 were taken into consideration; in the case of changes in expression level – range of SLR presented in Table 1.

Correlation between HSF1 binding and modulation of gene expression during heat shock

In the next step of analysis, we selected genes that both were bound by HSF1 and whose expression was changed in cells subjected to hyperthermia. We found significant changes of expression only in a fraction of genes that interacted with HSF1 (Figure 6C). In general, the binding of HSF1 was correlated with an activation of the transcription in hepatocytes and with a repression in spermatocytes (Table 7). Top HSF1-bound genes induced or repressed in spermatocytes and hepatocytes following heat shock are shown in Table 8 (all results are presented in Additional file 13[24]). In the selected genes HSF1 binding was usually observed both in hepatocytes and in spermatocytes. Nevertheless, similar HSF1 binding did not result in similar changes in expression when both types of cells were compared. Even in spermatocytes subjected to hyperthermia at 38°C and 43°C, HSF1-bound genes were differentially regulated after heat shock. This finding suggests that factors responsible for HSF1 activation are not only cell-type specific, but also dependent on the temperature. We also searched for a general statistical correlation between hyperthermia-induced changes in gene expression and HSF1 binding. All genes where in case of expression SLR ≠ 0 and in case of HSF1 binding AB1 > AB0 were included into analysis; results are presented in Additional file 14: Table S6 [24]. The correlation was found between expression and actual HSF1 binding at elevated temperatures but not at physiological temperature. It is noteworthy that a negative statistical correlation was observed for both up- and down-regulated genes in spermatocytes, i.e. enhanced HSF1 binding was associated with a decreased expression. In contrast, a positive statistical correlation was observed for up-regulated genes in hepatocytes, i.e. enhanced HSF1 binding was associated with an increased expression.

Table 7 Correlation of HSF1 binding following heat shock and subsequent changes in the transcription level
Table 8 Top HSF1-bound genes induced or repressed in spermatocytes and hepatocytes following heat shock

Finally, we searched for HSF1-regulated genes putatively involved in cell type- specific response to stress conditions. Among genes that bind HSF1 and are activated in spermatocytes treated at 38°C but not in heat shocked hepatocytes the most interesting seems to be Spen (SPEN homolog, transcriptional regulator (Drosophila)), which is a transcriptional repressor. Among genes that bind HSF1 and are negatively regulated in spermatocytes the biggest group are those involved in DNA replication and cell division, regulation of transcription and RNA processing, protein folding and degradation, and intracellular transport (listed and characterized in Additional file 15: Table S7 [24]). Interestingly, several genes that are bound by HSF1 and are negatively regulated in spermatocytes treated with hyperthermia at 38°C are essential for spermatogenesis (Table 9). It has been shown that their knock-out leads to defects in spermatogenesis. The most interesting in this group are genes involved in the repression of transposable elements: Tdrd1, Piwil2 and Mov10l1. We verified expression of some of these genes using quantitative RT-PCR analysis on independent material and confirmed significant down-regulation of Tdrd1 and Mov10l1 (Additional file 16: Figure S2 [24]).

Table 9 Characteristics of genes necessary for spermatogenesis whose repression following heat shock in spermatocytes correlates with HSF1 binding


The changes in the global gene expression profiles following heat shock were already studied in different organisms (from bacteria to mammals), and in different types of cells or tissues [8, 4046], even in mouse liver [45] and testes [4749]. Although some differences in heat shock response between somatic cells exist, in all of them pro-survival signaling pathways are activated. Their activation was also observed in mouse testes when changes in whole tissue were analyzed [47, 49]. The strategy used in the present study allowed us to find differences in the heat shock response between somatic and spermatogenic germ cells. By isolating the spermatogenic cells we avoided the overlapping of signals from different types of cells. We found that genes selected by Li et al. [49] as activated in whole testes during heat shock at 43°C (e.g. Hspa1a, Egr1, Fos, Socs3, Hspa8, Jun and Dnaja4) are not induced in spermatogenic cells, although they are strongly induced in hepatocytes. In fact, the lack of strong activation of Hsp genes in isolated spermatogenic cells documented in the present study is consistent with some previous findings [15, 16]. Hence, our data confirms basic differences between spermatogenic cells and somatic cells (exemplified by hepatocytes) in the transcriptional response to stress.

Spermatogenic cells do not activate “universal” pro-survival pathways in response to stress. However, because a large fraction of genes affected in spermatocytes during hyperthermia remains uncharacterized, identification of the signaling pathways involved in stress response is difficult. Additionally, considerably more genes are repressed than activated, and primarily pathways connected with RNA processing are strongly repressed under stress conditions. Posttranscriptional regulation is of critical importance during mammalian spermatogenesis. The testis has the greatest enrichment of tissue-specific splicing [50]. Also alternative promoters and polyA sites, antisense transcripts, necessity of the mRNA storage for a longer time before translation, and a different class of regulatory RNAs contribute to transcriptome diversification in testis [51]. Thus, deregulated RNA processing during heat shock could result in serious disturbances of spermatogenesis.

HSF1 binding to DNA was studied on a genome-wide scale in Saccharomyces and Drosophila cells or embryos [6, 7, 52, 53], as well as in mammalian cells [8, 41, 54], including mouse testis at a physiological temperature [55]. All these studies suggest that HSF1 may regulate diverse cellular processes that extend far beyond protein folding, the general role of HSPs that are the major targets of this factor. Here we found that several genes involved in processes other than protein folding (mainly in RNA processing and DNA replication, recombination and repair) were bound by HSF1, and that many genes were similarly bound by HSF1 in both hepatocytes and spermatocytes despite a completely different nature in the resulting transcriptional response. Recently, it has also been noted that HSF1 has important cell-specific functions at a normal physiological temperature. It regulates cancer-specific genes that support oncogenic processes in highly tumorigenic cells, in which it is more active than in less malignant cells [54]. In testis under physiologically normal conditions, HSF1 is required for transcriptional regulation of sex chromosomal multicopy genes [55]. Here we have observed that only a few genes were strongly bound by HSF1 in purified spermatocytes not subjected to hyperthermia.

Although HSF1 binding in hepatocytes and spermatocytes shows a relatively similar pattern, transcription of HSF1-bound genes is regulated in an apparently different way in these cells. Additionally, in spermatocytes heat shocked at 38°C and 43°C the changes in transcription are markedly different despite a rather similar pattern of HSF1 binding. Even though HSF1 binds more effectively to DNA in cells treated at 43°C, changes in the transcription level of HSF1-bound genes are much stronger in spermatocytes treated at 38°C. We have concluded that regulation of transcription by HSF1 binding is not only a cell type-specific mechanism, but it is also affected by temperature. HSF1 binding to the DNA in most cases does not correlate with significant changes in gene expression (on the other hand, HSF1 controls only a subset of the genes altered by heat shock). It is not surprising because the DNA-binding ability and the transactivation competence of HSF1 are regulated independently. HSF1 activity is modulated not only by interactions with chaperones but also by numerous posttranslational modifications including phosphorylation, acetylation, and sumoylation [56]. Additionally, independently of its transactivating activities, HSF1 can control the chromatin organization in response to heat shock [57]. HSF1 mediates a genome-wide and massive histone deacetylation interacting with HDAC1 and HDAC2. Moreover, the transcriptional response to stress may be regulated by microRNAs. Recently, Wilmink et al. [58] identified a group of miRNAs (mostly not annotated) differentially expressed following hyperthermia in dermal fibroblasts. One should not exclude the cell-type specificity of miRNA-related mechanisms in the regulation of heat shock response, which introduces an additional level of complexity in its regulation. In line with this possibility, expression of an extra class of small RNAs molecules involved in gene silencing, so called piRNA, is expressed specifically in spermatogenic cells [59].

The most unexpected finding of our study is that binding of HSF1 could result in inhibition of the transcription on a wide scale. HSF1 is well known as a transcriptional activator, and only a few reports have shown that HSF1 can act as a transcriptional repressor of specific genes involved in the acute phase response and inflammation (like Il1b, Tnf, Fos or Csf1; [60, 61]. However, in that case HSF1 antagonizes the activity of other transcription factors (by direct protein-protein interactions) and blocks the activation of the genes mentioned above, after stimulation (e.g. by lipopolysaccharide). Recently, it was suggested that HSF1 can inhibit transcription on a wider scale by interactions with Alu elements containing HSF binding sites [62]. Such sequences are frequently present in the corresponding genomic regions of the down-regulated transcripts in antisense orientation. Therefore, the antisense-mediated mechanism of inhibition of the transcription following heat shock was proposed. This study additionally implies that the direction of regulation (positive versus negative) is influenced by the location of the HSF1-binding site, what was also suggested by Mendillo et al. [54]. They stated that genes positively regulated by HSF1 were more likely bound at the promoter, whereas negatively regulated were more frequently bound in distal regions. Promoter tiling microarrays used in our study limited analyzes to promoter regions. Hence, it cannot be excluded that analyzing HSF1 binding by ChIP-seq could enable the identification of even more genes repressed following HSF1 binding. Furthermore, other transcription factors (that are differentially expressed in spermatocytes and hepatocytes) could interfere with changes in gene expression following hyperthermia, in either HSF1-dependent or independent mechanisms.

Among HSF1-bound genes negatively regulated in spermatocytes in response to hyperthermia, the biggest group is involved in DNA replication and cell division. They encode for proteins required for centrosome duplication, meiotic recombination and spindle organization, among others. Thus, disturbances of their expression can directly influence the progression of meiosis. Other HSF1-bound genes specifically down-regulated in spermatocytes are involved in RNA processing (synthesis, splicing and transport) and transcriptional regulation, which also may have a huge impact on spermatogenesis. HSF1-dependent pathways of protein folding and degradation were also differentially expressed in spermatocytes and hepatocytes. The greatest difference was noted in the case of Dnaja1 expression. DNAJA1 is generally a widely expressed co-chaperone that together with multiple HSP70 proteins takes part in protein folding [63]. In Dnaja1 null mice, defects in fertility were observed. However, the defects resulted from a lack of DNAJA1 in Sertoli cells rather than in spermatogenic cells. Although a lack of DNAJA1 in Sertoli cells leads to aberrant androgen signaling, disruption of Sertoli–germ cell contact and severe defects in spermiogenesis [31], it is not known what is a consequence of its lack in spermatogenic cells.

Functionally essential for spermatogenesis and male fertility are some other HSF1-bound genes negatively regulated in spermatocytes following heat shock (listed in Table 9). Arrest of spermiogenesis and increased apoptosis was observed in testes from Celf1 knock-out males [32]. CELF1 regulates the alternative splicing, which is of a greater importance to spermatogenic cells than to other cells. However it is broadly expressed and belongs to a bigger family of RNA binding proteins, that could possibly perform overlapping functions. Thus, is seems that heat shock-induced down-regulation of Celf1 (as well as Dnaja1) could be important, but not critical for induction of apoptosis in spermatocytes.

Disruption of spermatogenesis before the end of meiotic division was observed in Spo11, Piwil2, Tdrd1 and Mov10l1 knockout mice. All these genes, which we report to be repressed in HSF1-dependent manner, are expressed at a high level and almost exclusively in testes (according to BioGPS database; [64]). SPO11 is a type II like topoisomerase that generates double-strand breaks during meiosis. Consequently, mouse null Spo11 spermatocytes fail to synapse chromosomes and progress beyond the zygotene stage of meiosis [33, 34]. PIWIL2, TDRD1 and MOV10L1 are all involved in mobile elements silencing, which is carried out in the primordial mouse testis. Lack of these proteins results in de-repression of retrotransposons (mainly LINE1) [38, 65, 66]. PIWIL2 (MILI) belongs to the PIWI (P-element-induced wimpy testis) proteins family which also contains PIWIL1 (MIWI) and PIWIL4 (MIWI2). Following heat shock in spermatocytes we observed down regulation of Piwil1 and Piwil2, while Piwil4 expression was absent (it is expressed primarily in primordial germ cells; [67]). Although mutations for any of the three murine PIWI proteins confer male-specific sterility [35, 68, 69], the exact function of PIWIL1 and PIWIL2 proteins in adult testis is not known. They are a subset of the Argonaute proteins and associate specifically with small RNAs called PIWI-interacting RNAs (piRNAs) [70]. In the absence of PIWIL1 and PIWIL2 proteins, spermatocytes and spermatids are devoid of all piRNAs and spermatogenesis is terminally arrested during prophase I of meiosis [71].

PIWI proteins and piRNAs are enriched in the male germ-cell specific structures: the dense (sex) body and the chromatoid body. The dense body is associated with synapsis and the formation of the XY body during meiosis, however both processes are not affected in the absence of PIWI proteins [71]. The X and Y chromosomes in the dense body are transcriptionally silenced during meiosis at the pachytene stage [72], yet the sex chromosomal multicopy genes can escape the postmeiotic repression [73, 74]. Because PIWIL2 probably participates in chromatin relaxation during DNA repair [75], one could speculate that it might also relax chromatin in sex chromosomes to enable transcription of such genes. It has been shown that expression of multicopy genes located in the dense body is regulated by HSF1 and HSF2 [55, 76]. Because either activation of HSF1 or simultaneous lack of both HSF1 and HSF2 are connected with the diminished fertility of males [16, 17, 77], one could speculate that disturbances in the transcription of the sex chromosomal multicopy genes could be at least one reason for HSF1-related fertility problems in males. The chromatoid body contains abundant RNA binding proteins and has been implicated in the storage, metabolism, and cytoplasm-to-nuclear transport of mRNA, small RNA, and related proteins [78]. Multiple members of the tudor domain containing proteins are also components of the chromatoid body and can associate with PIWI proteins. TDRD1 specifically interacts with PIWIL2 [34, 36, 79]. In addition, PIWIL2 positively regulates Tdrd1 expression at the mRNA level [36]. It has been proposed that TDRD1 can act as a scaffold protein for complex assembly in the piRNA pathway [80]. MOV10L1 (a germ cell–specific putative RNA helicase) also associates with the PIWIL2-piRNA complex in adult testis [38, 39, 65]. Genetic disruption of the MOV10L1 RNA helicase domain in mice is connected with a lack of piRNAs and male sterility. Thus, down-regulation of Piwil2, Tdrd1 and Mov10l1 following heat shock is apparently connected with deregulated piRNA biogenesis.

Our study revealed that transcription of many genes important for spermatogenesis is negatively regulated following hyperthermia independently of HSF1 binding. It has been shown that lack of CDK2 [81], MLH3 [82], SYCP3 [83], MYBL1 [84], GJA1 [85], CLDN11 [86], BSG [87], POU5F2 [88], NRAS [89], CTSL [90], and possibly many others is connected with male sterility. We propose that repression of these genes (and other HSF1-independent genes) in cells subjected to hyperthermia results from activity of SPEN, which is specifically induced in heat shocked spermatocytes. SPEN is a large multidomain protein with ability to bind to DNA and RNA, and is a putative general transcriptional repressor . Importantly, SPEN can interact directly with SMRT (Silencing Mediator for Retinoid and Thyroid receptors; NCOR2) and members of the NuRD (nucleosome remodeling and deacetylating) complex including HDAC1 and HDAC2, MTA2, MBD3 and RbAp48 [91]. It was suggested that SPEN serves as a nuclear matrix platform that organizes and integrates transcriptional responses: recruits proteins involved in histone deacetylation, or acts by sequestration of transcriptional activators. Here we show that heat shock-induced up-regulation of Spen expression correlates with HSF1 binding to its promoter, which indicates its direct HSF1 dependence.


Although in somatic and spermatogenic cells after activation HSF1 binds to many promoters in a similar way, its impact on genes’ expression is completely different. This results in a completely different heat shock response in these two types of cells. In hepatocytes, which are able to survive following heat shock, HSF1 activation is connected rather with up-regulated gene expression, whereas in spermatogenic cells, which induce apoptosis, HSF1 activation is connected with repression of the transcription. HSF1-dependent induction of apoptosis in spermatocytes seems to be caused by the simultaneous repression of many genes essential for spermatogenesis. Some of these genes could be down-regulated due to direct HSF1 binding, while repression of other genes is putatively an indirect effect of HSF1, executed by transcription repressors regulated by this factor. Hence, a broad range of other transcription factors, some of them induced by HSF1, could also contribute to cell-type specific mechanisms regulating transcriptional response to hyperthermia. The obtained results suggest that deregulated RNA (especially piRNA) metabolism could be the most important factor triggering spermatogenic cell death following heat shock.


Experimental animals and ethical statement

Adult (10 to 16-week-old), traditional inbred FVB/N male mice (mean weight: 27.6 g ± 2.1 g) were used for analyses (source: Animal Facility in the Cancer Center and Institute of Oncology, adherent to a high standard of veterinary care). The animals (typically 2–8 animals per cage; minimum floor space per animal – 100 cm2) were maintained under controlled environmental conditions with a 12 h light:12 h darkness cycle and were provided with food and tap water. The animal experiments were carried out according to Polish legislation, and were approved by the Local Committee of Ethics and Animal Experimentation at the Medical University of Silesia in Katowice, Poland (Decision No 82/2009) and by the institutional animal care policy of the Cancer Center and Institute of Oncology (Gliwice, Poland).

Isolation of spermatocytes

The mice were killed by CO2 asphyxiation. The testes (twenty males for one isolation) were decapsulated and cell suspension was obtained by treatment with collagenase type IA (Sigma), trypsin and DNase I as described previously [92, 93]. Population of cells enriched in spermatocytes (up to 80%) was obtained by unit gravity sedimentation in a Sta-Put chamber (SP-120, ProScience Inc, Ontario, Canada) containing a 1-3% (w/v) linear BSA gradient in PBS [94]. The isolated fraction of spermatocytes was contaminated mainly with round spermatids due to their tendency to form multinucleated cells. Over 95% of the cells were viable as assessed by the exclusion of trypan blue. The yield of spermatocytes prepared in one isolation was usually ~5×107. Finally cells were suspended in CO2 saturated RPMI medium (1 × 106 cells per 1 ml) supplemented with 10% (v/v) fetal bovine serum, 0.004% (v/v) gentamycin (KRKA) and 6mM sodium lactate (Sigma). After each isolation cells were equally divided for control group and heat shocked at 38°C or 43°C groups. Totally, 12 isolations were performed.

Isolation of hepatocytes

Mouse hepatocytes were obtained from 10 to 16-week-old FVB/N males by the collagenase perfusion method [95]. Briefly, the liver was cannulated via the portal vein and perfused with 20 ml of warm (37°C) Krebs-Ringer buffer (5 mM HEPES, pH 7.4; 4,8 mM KCl; 1.2 mM MgSO4; 1.2 mM KH2PO4; 120 mM NaCl; 24 mM NaHCO3; 20 mM D-glucose) containing 0.5 mM EGTA at a rate of 5 ml/min with the perfusate exiting through the thoracic anterior vena cava. The liver was then perfused with a 20 ml of Krebs-Ringer buffer without EGTA, and then with a ~40 ml buffer containing 1.37 mM calcium chloride and 60–100 mg collagenase type IA (Sigma) for up to 10 min. After dispersion and washing in cold PBS (137 mM NaCl; 2.7 mM KCl; 8.1 mM Na2HPO4; 1.5 mM KH2PO4, pH 7.2), isolated hepatocytes were filtered through a 100-μm cell strainer. Cells were assessed for viability by trypan blue exclusion. A yield of 0.8 to 1×108 was routinely obtained with over 93% viability. Isolated hepatocytes were resuspended in DMEM medium supplemented with 10% (v/v) fetal bovine serum, 0.004% (v/v) gentamycin (KRKA) and 1 μM insulin. Totally, 12 isolations were performed. For ChIP experiments, after each isolation cells were equally divided for control and heat shocked at 43°C group.

Heat shock treatment

For global gene expression analyses spermatocytes from three independent isolations were heat shocked for one hour at 38°C or 43°C in a CO2 incubator and allowed to recover at the physiological temperature (32°C). Due to difficulties with the isolation of high quality intact RNA from hepatocytes after heat shock done in vitro, whole-body hyperthermia was performed in vivo. Animals were anesthetized using avertin (15–17 μl of a 2.5% solution/g body weight) injected intraperitoneally. The lower half of the torso of each animal was submerged in a water bath at 43°C for 30 min, after which time the animals were dried and returned to their cages. Hepatocytes were collected by collagenase perfusion to obtain cells ready for harvesting after two hours from heat shock. Cells were washed in cold PBS and frozen as a pellet at -70°C until RNA isolation. For RT-PCR analyses, organs (testes and livers) were collected either from control, untreated, or in vivo heat-shocked animals, and frozen on dry ice. For ChIP experiments, an equal volume (10 ml) of CO2 saturated, pre-heated media (to 55°C and 60°C for hepatocytes and spermatocytes, respectively) were added to the cell suspension, which immediately raised the temperature of the media to 42-43°C. The tubes were sealed with parafilm and submerged in a water bath at 43°C for an additional 5 to 40 min. Spermatocytes were heat shocked also at 38°C in a similar way except that the medium was pre-heated to 53°C. Immediately after heat shock cells were fixed by adding formaldehyde according to ChIP protocols. Simultaneously the medium was quickly cooled to room temperature.

RNA isolation, cDNA synthesis and RT-PCR

For global gene expression analyses, the total RNA was purified in triplicate from isolated spermatocytes and hepatocytes using RNeasy® Mini Kit spin columns (Qiagen, Valencia, CA, USA). Validation of the RNA quality was performed using the RNA 6000 Nano Assay on an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). For RT-PCR, the total RNA was isolated from independent samples (isolated spermatocytes, testes, livers) at least twice using the GeneMATRIX Universal RNA Purification Kit (Eurx, Gdańsk, Poland), and then 25–50 μg of RNA was incubated with 2.5 U of DNase I (Worthington Biochemical Corporation, NY, USA) for 15–30 min at 37°C in a buffer containing 40 mM Tris pH 7.9; 10 mM NaCl; 6 mM MgCl2; 0.1 mM CaCl2 to remove genomic DNA contamination. Reaction was stopped by heating for 5 min at 75°C and RNA was precipitated, washed, dissolved in RNase free water and tested for DNA contamination by PCR (35 cycles) with Gapdh primers. For cDNA synthesis 1 μg of RNA was primed with random pentadecamer primers (200 ng) [96] and 18-mer oligo-dT primers (250 ng) and the reaction was carried out at 37°C for 50 min in the presence of 5μM dNTPs and 200 U of reverse transcriptase (Gibco, Invitrogen) in a total volume of 20 μl. Reaction was terminated by heating at 70°C for 15 min, samples were diluted using water to 400 μl, aliquoted and stored at -20°C until use. For RT-PCR 1–3 μl of cDNA template was used and 14–35 cycles were applied depending on the primers set. For each primers set an optimal annealing temperature was established. Quantitative RT-PCR was performed using a Bio-Rad CFX 96™ Real-Time PCR Detection System. A total of 5 pmoles of forward and reverse primer, cDNA template, and Real-Time 2× PCR Master Mix SYBR A (A&A Biotechnology, Gdynia, Poland) were used. Reactions were incubated at 95°C for 3 min, followed by 40 cycles of 95°C for 25 sec, 61°C for 30 sec and 72°C for 30 sec. Analyzed genes were normalized against Gapdh. Primers used in analyses are listed in Additional file 17: Table S8 [24].

Gene expression profiling by DNA microarrays

The global gene expression analysis was performed using the Affymetrix GeneChip® Mouse Genome 430 2.0 Array and GeneChip® Expression 3′Amplification reagents (Affymetrix, Santa Clara, CA) according to manufacturer’s protocols if not stated otherwise. The total RNA was extracted (from three independent samples from each cell type and each experimental point) as described above. cDNA was synthesized from 8 μg of total RNA using the One-Cycle cDNA Synthesis Kit with SuperScriptII polymerase then reverse-transcribed with a T7 RNA polymerase incorporating biotin-labeled nucleotides using the GeneChip IVT Labeling Kit. Biotin-labeled cRNA probe samples were purified, fragmented and hybridized (hybridization cocktail was prepared without DMSO) at 45°C overnight to the array. The arrays were washed, stained and scanned using a Fluidics Station and microarray scanner, which are components of the Affymetrix GeneChip Instrument System. Analyses were done in triplicate. Data was analyzed with R/Bioconductor and Matlab softwares, RMA (Robust Multi-array Analysis) algorithm with Dai (ver. 15.1) annotation were used for data normalization [97, 98]. The level of genes expression was shown in arbitrary units in the logarithmic scale (log2). Noise threshold was estimated by modeling the signal distribution as Gaussian mixture (GM) and defining the first component as a model of noise. The remaining GM components were treated as the models of low and high signals (Additional file 1: Table S1 [24]). Changes in the level of genes expression in heat-shocked samples versus control, untreated, were presented as a signal log ratio (SLR, the average value of HS/average value of C). The GM model applied to SLR distribution allowed for identification of significantly up and down regulated genes (Table 1). Gene clustering based on temperature-related changes was performed with the use of multidimensional GM modeling [99] and the obtained group of genes were tested for overrepresentation of GO terms with the use of conditional hypergeometric test. Characteristics of individual genes biological functions were made on the basis of information available on the websites: [100, 101].

Chromatin Immunoprecipitation (ChIP) Assays and ChIP on chip

For analyses of HSF1 binding kinetics, the ChIP assay was carried out using a ChIP kit (Upstate Biotechnology, Lake Placid, NY) following the manufacturer’s protocol. To obtain recommended shared DNA size fragments (between 200–1,000 bp) chromatin from heat shocked samples had to be sonicated using two times more pulses than chromatin from control samples. For 1 × 106 cells, 0.3 to 3 μg of rabbit anti-HSF1 polyclonal antibody (SPA-901, Stressgen, Victoria, Canada) was tested with similar qualitative results. Immunoprecipitated DNA was analyzed by PCR (ChIP-PCR). Primers characteristics used in analyses are presented in Additional file 18: Table S9 [24]. ChIP on chip analyses were done according to the Affymetrix protocol [102], with one modification: in the second round of immunoprecipitated DNA amplification, 2 mM magnesium chloride was applied instead of the recommended 0.75 mM. For one immunoprecipitation (IP), 5 × 107 cells and 15 μg of anti-HSF1 antibody (SPA-901, Stressgen) were used (sample marked as AB1). A negative control was performed using the same number of cells with no antibody (mock IP, sample marked as AB0). Immunoprecipitated DNA was analyzed on the GeneChip® Mouse Promoter 1.0R Array. ChIP on chip analyses were done in triplicate. Data was quantile normalized with the use of dCHIP software and Mm_PromPR_v02-1_NCBIv36.NR.bpmap file. 5th percentile was used for background subtraction (PM-only). The further analysis was done in Matlab environment, the subsets of probes related to the regions defined by gene Accession Number were analyzed individually. The statistical tests for signal enhancement within the region and nonuniformity of the signal AB1-AB0 difference distribution were performed. Ontology of HSF1 bound genes was analyzed by NucleoAnnot application created within the confines of the GENEPI LowRT project (available on [103]).

Availability of supporting data

Expression microarray and ChIP on chip data discussed in this article have been deposited in NCBI’s Gene Expression Omnibus [104] and are accessible through GEO Series accession number GSE40248 and GSE40390. The data sets supporting the results of this article are available in the LabArchives repository [24].



ChIP samples precipitated with/without antibody


Bovine serum albumin



ChIP on chip:

Chromatin immunoprecipitation assay combined with promoter microarrays


Chromatin immunoprecipitation followed by PCR


ChIP with massively parallel DNA sequencing


False discovery rate


Gene ontology


Gaussian mixture




Heat shock proteins


Heat shock factor


Reverse transcription polymerase chain reaction


Robust multi-array analysis


Signal log ratio.


  1. Kampinga HH, Hageman J, Vos MJ, Kubota H, Tanguay RM, Bruford EA, Cheetham ME, Chen B, Hightower LE: Guidelines for the nomenclature of the human heat shock proteins. Cell Stress Chaperones. 2009, 14: 105-111. 10.1007/s12192-008-0068-7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Christians ES, Zhou Q, Renard J, Benjamin IJ: Heat shock proteins in mammalian development. Semin Cell Dev Biol. 2003, 14: 283-290. 10.1016/j.semcdb.2003.09.021.

    Article  CAS  PubMed  Google Scholar 

  3. Rupik W, Jasik K, Bembenek J, Widłak W: The expression patterns of heat shock genes and proteins and their role during vertebrate’s development. Comp Biochem Physiol A Mol Integr Physiol. 2011, 159: 349-366. 10.1016/j.cbpa.2011.04.002.

    Article  PubMed  Google Scholar 

  4. Lindquist S, Craig EA: The heat-shock proteins. Annu Rev Genet. 1988, 22: 631-677. 10.1146/

    Article  CAS  PubMed  Google Scholar 

  5. Voellmy R: On mechanisms that control heat shock transcription factor activity in metazoan cells. Cell Stress Chaperones. 2004, 9: 122-133. 10.1379/CSC-14R.1.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Hahn JS, Hu Z, Thiele DJ, Iyer VR: Genome-wide analysis of the biology of stress responses through heat shock transcription factor. Mol Cell Biol. 2004, 24: 5249-5256. 10.1128/MCB.24.12.5249-5256.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Birch-Machin I, Gao S, Huen D, McGirr R, White RA, Russell S: Genomic analysis of heat-shock factor targets in Drosophila. Genome Biol. 2005, 6: R63-10.1186/gb-2005-6-7-r63.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Page TJ, Sikder D, Yang L, Pluta L, Wolfinger RD, Kodadek T, Thomas RS: Genome-wide analysis of human HSF1 signaling reveals a transcriptional program linked to cellular adaptation and survival. Mol Biosyst. 2006, 2: 627-639. 10.1039/b606129j.

    Article  CAS  PubMed  Google Scholar 

  9. Toma A, Widlak W, Vydra N: The role of heat shock transcription factor 1 in carcinogenesis. Postepy Biol Komorki. 2012, 39: 269-288.

    CAS  Google Scholar 

  10. Curci A, Bevilacqua A, Mangia F: Lack of heat-shock response in preovulatory mouse oocytes. Dev Biol. 1987, 123: 154-160. 10.1016/0012-1606(87)90437-4.

    Article  CAS  PubMed  Google Scholar 

  11. Fisher BR, Heredia DJ, Brown KM: Induction of hsp72 in heat-treated rat embryos: a tissue-specific response. Teratology. 1995, 52: 90-100. 10.1002/tera.1420520205.

    Article  CAS  PubMed  Google Scholar 

  12. Marcuccilli CJ, Mathur SK, Morimoto RI, Miller RJ: Regulatory differences in the stress response of hippocampal neurons and glial cells after heat shock. J Neurosci. 1996, 16: 478-485.

    CAS  PubMed  Google Scholar 

  13. Yin Y, Hawkins KL, DeWolf WC, Morgentaler A: Heat stress causes testicular germ cell apoptosis in adult mice. J Androl. 1997, 18: 159-165.

    CAS  PubMed  Google Scholar 

  14. Davidson S, Hoj P, Gabriele T, Anderson RL: In vivo growth of a murine lymphoma cell line alters regulation of expression of HSP72. Mol Cell Biol. 1995, 15: 1071-1078.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Izu H, Inouye S, Fujimoto M, Shiraishi K, Naito K, Nakai A: Heat shock transcription factor 1 is involved in quality-control mechanisms in male germ cells. Biol Reprod. 2004, 70: 18-24.

    Article  CAS  PubMed  Google Scholar 

  16. Vydra N, Malusecka E, Jarzab M, Lisowska K, Glowala-Kosinska M, Benedyk K, Widlak P, Krawczyk Z, Widlak W: Spermatocyte-specific expression of constitutively active heat shock factor 1 induces HSP70i-resistant apoptosis in male germ cells. Cell Death Differ. 2006, 13: 212-222. 10.1038/sj.cdd.4401758.

    Article  CAS  PubMed  Google Scholar 

  17. Nakai A, Suzuki M, Tanabe M: Arrest of spermatogenesis in mice expressing an active heat shock transcription factor 1. Embo J. 2000, 19: 1545-1554. 10.1093/emboj/19.7.1545.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Widlak W, Benedyk K, Vydra N, Glowala M, Scieglinska D, Malusecka E, Nakai A, Krawczyk Z: Expression of a constitutively active mutant of heat shock factor 1 under the control of testis-specific hst70 gene promoter in transgenic mice induces degeneration of seminiferous epithelium. Acta Biochim Pol. 2003, 50: 535-541.

    CAS  PubMed  Google Scholar 

  19. Widlak W, Vydra N, Malusecka E, Dudaladava V, Winiarski B, Scieglinska D, Widlak P: Heat shock transcription factor 1 down-regulates spermatocyte-specific 70 kDa heat shock protein expression prior to the induction of apoptosis in mouse testes. Genes Cells. 2007, 12: 487-499. 10.1111/j.1365-2443.2007.01069.x.

    Article  CAS  PubMed  Google Scholar 

  20. Hess RA: Spermatogenesis, Overview. Encyclopedia of Reproduction. Edited by: Knobil E, Neill JD. 1998, Academic Press: San Diego, 539-545.

    Google Scholar 

  21. Ivell R: Lifestyle impact and the biology of the human scrotum. Reprod Biol Endocrinol. 2007, 5: 15-10.1186/1477-7827-5-15.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Sarge KD: Male germ cell-specific alteration in temperature set point of the cellular stress response. J Biol Chem. 1995, 270: 18745-18748.

    Article  CAS  PubMed  Google Scholar 

  23. Widlak W, Winiarski B, Krawczyk A, Vydra N, Malusecka E, Krawczyk Z: Inducible 70 kDa heat shock protein does not protect spermatogenic cells from damage induced by cryptorchidism. Int J Androl. 2007, 30: 80-87. 10.1111/j.1365-2605.2006.00713.x.

    Article  CAS  PubMed  Google Scholar 

  24. LabArchives repository.,

  25. Hikim AP, Lue Y, Yamamoto CM, Vera Y, Rodriguez S, Yen PH, Soeng K, Wang C, Swerdloff RS: Key apoptotic pathways for heat-induced programmed germ cell death in the testis. Endocrinology. 2003, 144: 3167-3175. 10.1210/en.2003-0175.

    Article  CAS  PubMed  Google Scholar 

  26. Ding Y, Lee JF, Lu H, Lee MH, Yan DH: Interferon-inducible protein IFIXalpha1 functions as a negative regulator of HDM2. Mol Cell Biol. 2006, 26: 1979-1996. 10.1128/MCB.26.5.1979-1996.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Yao J, Munson KM, Webb WW, Lis JT: Dynamics of heat shock factor association with native gene loci in living cells. Nature. 2006, 442: 1050-1053. 10.1038/nature05025.

    Article  CAS  PubMed  Google Scholar 

  28. Inouye S, Izu H, Takaki E, Suzuki H, Shirai M, Yokota Y, Ichikawa H, Fujimoto M, Nakai A: Impaired IgG production in mice deficient for heat shock transcription factor 1. J Biol Chem. 2004, 279: 38701-38709. 10.1074/jbc.M405986200.

    Article  CAS  PubMed  Google Scholar 

  29. Trinklein ND, Chen WC, Kingston RE, Myers RM: Transcriptional regulation and binding of heat shock factor 1 and heat shock factor 2 to 32 human heat shock genes during thermal stress and differentiation. Cell Stress Chaperones. 2004, 9: 21-28.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Rossi A, Trotta E, Brandi R, Arisi I, Coccia M, Santoro MG: AIRAP, a new human heat shock gene regulated by heat shock factor 1. J Biol Chem. 2010, 285: 13607-13615. 10.1074/jbc.M109.082693.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Terada K, Yomogida K, Imai T, Kiyonari H, Takeda N, Kadomatsu T, Yano M, Aizawa S, Mori M: A type I DnaJ homolog, DjA1, regulates androgen receptor signaling and spermatogenesis. Embo J. 2005, 24: 611-622. 10.1038/sj.emboj.7600549.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Kress C, Gautier-Courteille C, Osborne HB, Babinet C, Paillard L: Inactivation of CUG-BP1/CELF1 causes growth, viability, and spermatogenesis defects in mice. Mol Cell Biol. 2007, 27: 1146-1157. 10.1128/MCB.01009-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Baudat F, Manova K, Yuen JP, Jasin M, Keeney S: Chromosome synapsis defects and sexually dimorphic meiotic progression in mice lacking Spo11. Mol Cell. 2000, 6: 989-998. 10.1016/S1097-2765(00)00098-8.

    Article  CAS  PubMed  Google Scholar 

  34. Romanienko PJ, Camerini-Otero RD: The mouse Spo11 gene is required for meiotic chromosome synapsis. Mol Cell. 2000, 6: 975-987.

    Article  CAS  PubMed  Google Scholar 

  35. Kuramochi-Miyagawa S, Kimura T, Ijiri TW, Isobe T, Asada N, Fujita Y, Ikawa M, Iwai N, Okabe M, Deng W, Lin H, Matsuda Y, Nakano T: Mili, a mammalian member of piwi family gene, is essential for spermatogenesis. Development. 2004, 131: 839-849. 10.1242/dev.00973.

    Article  CAS  PubMed  Google Scholar 

  36. Wang J, Saxe JP, Tanaka T, Chuma S, Lin H: Mili interacts with tudor domain-containing protein 1 in regulating spermatogenesis. Curr Biol. 2009, 19: 640-644. 10.1016/j.cub.2009.02.061.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Chuma S, Hosokawa M, Kitamura K, Kasai S, Fujioka M, Hiyoshi M, Takamune K, Noce T, Nakatsuji N: Tdrd1/Mtr-1, a tudor-related gene, is essential for male germ-cell differentiation and nuage/germinal granule formation in mice. Proc Natl Acad Sci USA. 2006, 103: 15894-15899. 10.1073/pnas.0601878103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Zheng K, Xiol J, Reuter M, Eckardt S, Leu NA, McLaughlin KJ, Stark A, Sachidanandam R, Pillai RS, Wang PJ: Mouse MOV10L1 associates with Piwi proteins and is an essential component of the Piwi-interacting RNA (piRNA) pathway. Proc Natl Acad Sci USA. 2010, 107: 11841-11846. 10.1073/pnas.1003953107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Frost RJ, Hamra FK, Richardson JA, Qi X, Bassel-Duby R, Olson EN: MOV10L1 is necessary for protection of spermatocytes against retrotransposons by Piwi-interacting RNAs. Proc Natl Acad Sci USA. 2010, 107: 11847-11852. 10.1073/pnas.1007158107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Murray JI, Whitfield ML, Trinklein ND, Myers RM, Brown PO, Botstein D: Diverse and specific gene expression responses to stresses in cultured human cells. Mol Biol Cell. 2004, 15: 2361-2374. 10.1091/mbc.E03-11-0799.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Trinklein ND, Murray JI, Hartman SJ, Botstein D, Myers RM: The role of heat shock transcription factor 1 in the genome-wide regulation of the mammalian heat shock response. Mol Biol Cell. 2004, 15: 1254-1261.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Rezzonico E, Lariani S, Barretto C, Cuanoud G, Giliberti G, Delley M, Arigoni F, Pessi G: Global transcriptome analysis of the heat shock response of Bifidobacterium longum. FEMS Microbiol Lett. 2007, 271: 136-145. 10.1111/j.1574-6968.2007.00704.x.

    Article  CAS  PubMed  Google Scholar 

  43. Laramie JM, Chung TP, Brownstein B, Stormo GD, Cobb JP: Transcriptional profiles of human epithelial cells in response to heat: computational evidence for novel heat shock proteins. Shock. 2008, 29: 623-630.

    CAS  PubMed  Google Scholar 

  44. Beckham JT, Wilmink GJ, Opalenik SR, Mackanos MA, Abraham AA, Takahashi K, Contag CH, Takahashi T, Jansen ED: Microarray analysis of cellular thermotolerance. Lasers Surg Med. 2010, 42: 752-765. 10.1002/lsm.20963.

    Article  PubMed  Google Scholar 

  45. Vallanat B, Anderson SP, Brown-Borg HM, Ren H, Kersten S, Jonnalagadda S, Srinivasan R, Corton JC: Analysis of the heat shock response in mouse liver reveals transcriptional dependence on the nuclear receptor peroxisome proliferator-activated receptor alpha (PPARalpha). BMC Genomics. 2010, 11: 16-10.1186/1471-2164-11-16.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Fiszer-Kierzkowska A, Vydra N, Wysocka-Wycisk A, Kronekova Z, Jarzab M, Lisowska KM, Krawczyk Z: Liposome-based DNA carriers may induce cellular stress response and change gene expression pattern in transfected cells. BMC Mol Biol. 2011, 12: 27-10.1186/1471-2199-12-27.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Rockett JC, Mapp FL, Garges JB, Luft JC, Mori C, Dix DJ: Effects of hyperthermia on spermatogenesis, apoptosis, gene expression, and fertility in adult male mice. Biol Reprod. 2001, 65: 229-239. 10.1095/biolreprod65.1.229.

    Article  CAS  PubMed  Google Scholar 

  48. Li YC, Hu XQ, Xiao LJ, Hu ZY, Guo J, Zhang KY, Song XX, Liu YX: An oligonucleotide microarray study on gene expression profile in mouse testis of experimental cryptorchidism. Front Biosci. 2006, 11: 2465-2482. 10.2741/1983.

    Article  CAS  PubMed  Google Scholar 

  49. Li Y, Zhou Q, Hively R, Yang L, Small C, Griswold MD: Differential gene expression in the testes of different murine strains under normal and hyperthermic conditions. J Androl. 2009, 30: 325-337.

    Article  PubMed Central  PubMed  Google Scholar 

  50. Xu Q, Modrek B, Lee C: Genome-wide detection of tissue-specific alternative splicing in the human transcriptome. Nucleic Acids Res. 2002, 30: 3754-3766. 10.1093/nar/gkf492.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Lee TL, Pang AL, Rennert OM, Chan WY: Genomic landscape of developing male germ cells. Birth Defects Res C Embryo Today. 2009, 87: 43-63. 10.1002/bdrc.20147.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Guertin MJ, Lis JT: Chromatin landscape dictates HSF binding to target DNA elements. PLoS Genet. 2010, 6:

    Google Scholar 

  53. Gonsalves SE, Moses AM, Razak Z, Robert F, Westwood JT: Whole-genome analysis reveals that active heat shock factor binding sites are mostly associated with non-heat shock genes in Drosophila melanogaster. PLoS One. 2011, 6: e15934-10.1371/journal.pone.0015934.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Mendillo ML, Santagata S, Koeva M, Bell GW, Hu R, Tamimi RM, Fraenkel E, Ince TA, Whitesell L, Lindquist S: HSF1 drives a transcriptional program distinct from heat shock to support highly malignant human cancers. Cell. 2012, 150: 549-562. 10.1016/j.cell.2012.06.031.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Akerfelt M, Vihervaara A, Laiho A, Conter A, Christians ES, Sistonen L, Henriksson E: Heat shock transcription factor 1 localizes to sex chromatin during meiotic repression. J Biol Chem. 2010, 285: 34469-34476. 10.1074/jbc.M110.157552.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Anckar J, Sistonen L: Regulation of HSF1 function in the heat stress response: implications in aging and disease. Annu Rev Biochem. 2011, 80: 1089-1115. 10.1146/annurev-biochem-060809-095203.

    Article  CAS  PubMed  Google Scholar 

  57. Fritah S, Col E, Boyault C, Govin J, Sadoul K, Chiocca S, Christians E, Khochbin S, Jolly C, Vourc’h C: Heat-shock factor 1 controls genome-wide acetylation in heat-shocked cells. Mol Biol Cell. 2009, 20: 4976-4984. 10.1091/mbc.E09-04-0295.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Wilmink GJ, Roth CL, Ibey BL, Ketchum N, Bernhard J, Cerna CZ, Roach WP: Identification of microRNAs associated with hyperthermia-induced cellular stress response. Cell Stress Chaperones. 2010, 15: 1027-1038. 10.1007/s12192-010-0189-7.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Siomi MC, Sato K, Pezic D, Aravin AA: PIWI-interacting small RNAs: the vanguard of genome defence. Nat Rev Mol Cell Biol. 2011, 12: 246-258. 10.1038/nrm3089.

    Article  CAS  PubMed  Google Scholar 

  60. Xie Y, Chen C, Stevenson MA, Auron PE, Calderwood SK: Heat shock factor 1 represses transcription of the IL-1beta gene through physical interaction with the nuclear factor of interleukin 6. J Biol Chem. 2002, 277: 11802-11810. 10.1074/jbc.M109296200.

    Article  CAS  PubMed  Google Scholar 

  61. Xie Y, Zhong R, Chen C, Calderwood SK: Heat shock factor 1 contains two functional domains that mediate transcriptional repression of the c-fos and c-fms genes. J Biol Chem. 2003, 278: 4687-4698. 10.1074/jbc.M210189200.

    Article  CAS  PubMed  Google Scholar 

  62. Pandey R, Mandal AK, Jha V, Mukerji M: Heat shock factor binding in Alu repeats expands its involvement in stress through an antisense mechanism. Genome Biol. 2011, 12: R117-10.1186/gb-2011-12-11-r117.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  63. Terada K, Oike Y: Multiple molecules of Hsc70 and a dimer of DjA1 independently bind to an unfolded protein. J Biol Chem. 2010, 285: 16789-16797. 10.1074/jbc.M110.101501.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  64. BioGPS database.,

  65. Vagin VV, Wohlschlegel J, Qu J, Jonsson Z, Huang X, Chuma S, Girard A, Sachidanandam R, Hannon GJ, Aravin AA: Proteomic analysis of murine Piwi proteins reveals a role for arginine methylation in specifying interaction with Tudor family members. Genes Dev. 2009, 23: 1749-1762. 10.1101/gad.1814809.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. De Fazio S, Bartonicek N, Di Giacomo M, Abreu-Goodger C, Sankar A, Funaya C, Antony C, Moreira PN, Enright AJ, O’Carroll D: The endonuclease activity of Mili fuels piRNA amplification that silences LINE1 elements. Nature. 2011, 480: 259-263. 10.1038/nature10547.

    Article  CAS  PubMed  Google Scholar 

  67. Aravin AA, Sachidanandam R, Bourc’his D, Schaefer C, Pezic D, Toth KF, Bestor T, Hannon GJ: A piRNA pathway primed by individual transposons is linked to de novo DNA methylation in mice. Mol Cell. 2008, 31: 785-799. 10.1016/j.molcel.2008.09.003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. Deng W, Lin H: miwi, a murine homolog of piwi, encodes a cytoplasmic protein essential for spermatogenesis. Dev Cell. 2002, 2: 819-830. 10.1016/S1534-5807(02)00165-X.

    Article  CAS  PubMed  Google Scholar 

  69. Carmell MA, Girard A, van de Kant HJ, Bourc’his D, Bestor TH, de Rooij DG, Hannon GJ: MIWI2 is essential for spermatogenesis and repression of transposons in the mouse male germline. Dev Cell. 2007, 12: 503-514. 10.1016/j.devcel.2007.03.001.

    Article  CAS  PubMed  Google Scholar 

  70. Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8: 272-285.

    Article  CAS  PubMed  Google Scholar 

  71. Beyret E, Lin H: Pinpointing the expression of piRNAs and function of the PIWI protein subfamily during spermatogenesis in the mouse. Dev Biol. 2011, 355: 215-226. 10.1016/j.ydbio.2011.04.021.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Turner JM: Meiotic sex chromosome inactivation. Development. 2007, 134: 1823-1831. 10.1242/dev.000018.

    Article  CAS  PubMed  Google Scholar 

  73. Toure A, Clemente EJ, Ellis P, Mahadevaiah SK, Ojarikre OA, Ball PA, Reynard L, Loveland KL, Burgoyne PS, Affara NA: Identification of novel Y chromosome encoded transcripts by testis transcriptome analysis of mice with deletions of the Y chromosome long arm. Genome Biol. 2005, 6: R102-10.1186/gb-2005-6-12-r102.

    Article  PubMed Central  PubMed  Google Scholar 

  74. Mueller JL, Mahadevaiah SK, Park PJ, Warburton PE, Page DC, Turner JM: The mouse X chromosome is enriched for multicopy testis genes showing postmeiotic expression. Nat Genet. 2008, 40: 794-799. 10.1038/ng.126.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. Yin DT, Wang Q, Chen L, Liu MY, Han C, Yan Q, Shen R, He G, Duan W, Li JJ, Wani A, Gao JX: Germline stem cell gene PIWIL2 mediates DNA repair through relaxation of chromatin. PLoS One. 2011, 6: e27154-10.1371/journal.pone.0027154.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Akerfelt M, Henriksson E, Laiho A, Vihervaara A, Rautoma K, Kotaja N, Sistonen L: Promoter ChIP-chip analysis in mouse testis reveals Y chromosome occupancy by HSF2. Proc Natl Acad Sci USA. 2008, 105: 11224-11229. 10.1073/pnas.0800620105.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  77. Wang G, Ying Z, Jin X, Tu N, Zhang Y, Phillips M, Moskophidis D, Mivechi NF: Essential requirement for both hsf1 and hsf2 transcriptional activity in spermatogenesis and male fertility. Genesis. 2004, 38: 66-80. 10.1002/gene.20005.

    Article  PubMed  Google Scholar 

  78. Meikar O, Da Ros M, Korhonen H, Kotaja N: Chromatoid body and small RNAs in male germ cells. Reproduction. 2011, 142: 195-209. 10.1530/REP-11-0057.

    Article  CAS  PubMed  Google Scholar 

  79. Reuter M, Chuma S, Tanaka T, Franz T, Stark A, Pillai RS: Loss of the Mili-interacting Tudor domain-containing protein-1 activates transposons and alters the Mili-associated small RNA profile. Nat Struct Mol Biol. 2009, 16: 639-646. 10.1038/nsmb.1615.

    Article  CAS  PubMed  Google Scholar 

  80. Mathioudakis N, Palencia A, Kadlec J, Round A, Tripsianes K, Sattler M, Pillai RS, Cusack S: The multiple Tudor domain-containing protein TDRD1 is a molecular scaffold for mouse Piwi proteins and piRNA biogenesis factors. Rna. 2012, 18: 2056-2072. 10.1261/rna.034181.112.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  81. Berthet C, Aleem E, Coppola V, Tessarollo L, Kaldis P: Cdk2 knockout mice are viable. Curr Biol. 2003, 13: 1775-1785. 10.1016/j.cub.2003.09.024.

    Article  CAS  PubMed  Google Scholar 

  82. Lipkin SM, Moens PB, Wang V, Lenzi M, Shanmugarajah D, Gilgeous A, Thomas J, Cheng J, Touchman JW, Green ED, Schwartzberg P, Collins FS, Cohen PE: Meiotic arrest and aneuploidy in MLH3-deficient mice. Nat Genet. 2002, 31: 385-390.

    CAS  PubMed  Google Scholar 

  83. Yuan L, Liu JG, Zhao J, Brundell E, Daneholt B, Hoog C: The murine SCP3 gene is required for synaptonemal complex assembly, chromosome synapsis, and male fertility. Mol Cell. 2000, 5: 73-83. 10.1016/S1097-2765(00)80404-9.

    Article  CAS  PubMed  Google Scholar 

  84. Toscani A, Mettus RV, Coupland R, Simpkins H, Litvin J, Orth J, Hatton KS, Reddy EP: Arrest of spermatogenesis and defective breast development in mice lacking A-myb. Nature. 1997, 386: 713-717. 10.1038/386713a0.

    Article  CAS  PubMed  Google Scholar 

  85. Roscoe WA, Barr KJ, Mhawi AA, Pomerantz DK, Kidder GM: Failure of spermatogenesis in mice lacking connexin43. Biol Reprod. 2001, 65: 829-838. 10.1095/biolreprod65.3.829.

    Article  CAS  PubMed  Google Scholar 

  86. Mazaud-Guittot S, Meugnier E, Pesenti S, Wu X, Vidal H, Gow A, Le Magueresse-Battistoni B: Claudin 11 deficiency in mice results in loss of the Sertoli cell epithelial phenotype in the testis. Biol Reprod. 2011, 82: 202-213.

    Article  Google Scholar 

  87. Igakura T, Kadomatsu K, Kaname T, Muramatsu H, Fan QW, Miyauchi T, Toyama Y, Kuno N, Yuasa S, Takahashi M, Senda T, Taguchi O, Yamamura K, Arimura K, Muramatsu T: A null mutation in basigin, an immunoglobulin superfamily member, indicates its important roles in peri-implantation development and spermatogenesis. Dev Biol. 1998, 194: 152-165. 10.1006/dbio.1997.8819.

    Article  CAS  PubMed  Google Scholar 

  88. Pearse RV, Drolet DW, Kalla KA, Hooshmand F, Bermingham JR, Rosenfeld MG: Reduced fertility in mice deficient for the POU protein sperm-1. Proc Natl Acad Sci USA. 1997, 94: 7555-7560. 10.1073/pnas.94.14.7555.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  89. Mangues R, Seidman I, Pellicer A, Gordon JW: Tumorigenesis and male sterility in transgenic mice expressing a MMTV/N-ras oncogene. Oncogene. 1990, 5: 1491-1497.

    CAS  PubMed  Google Scholar 

  90. Wright WW, Smith L, Kerr C, Charron M: Mice that express enzymatically inactive cathepsin L exhibit abnormal spermatogenesis. Biol Reprod. 2003, 68: 680-687.

    Article  CAS  PubMed  Google Scholar 

  91. Shi Y, Downes M, Xie W, Kao HY, Ordentlich P, Tsai CC, Hon M, Evans RM: Sharp, an inducible cofactor that integrates nuclear receptor repression and activation. Genes Dev. 2001, 15: 1140-1151. 10.1101/gad.871201.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  92. Bellve AR, Cavicchia JC, Millette CF, O’Brien DA, Bhatnagar YM, Dym M: Spermatogenic cells of the prepuberal mouse. Isolation and morphological characterization. J Cell Biol. 1977, 74: 68-85. 10.1083/jcb.74.1.68.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  93. Abou-Haila A, Tulsiani DR: Acid Glycohydrolases in Rat Spermatocytes, Spermatids and Spermatozoa: Enzyme Activities, Biosynthesis and Immunolocalization. Biol Proced Online. 2001, 3: 35-42. 10.1251/bpo21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  94. Millette CF, Moulding CT: Cell surface marker proteins during mouse spermatogenesis: two-dimensional electrophoretic analysis. J Cell Sci. 1981, 48: 367-382.

    CAS  PubMed  Google Scholar 

  95. Rokita H, Szuba K: Regulation of acute phase reaction by transforming growth factor beta in cultured murine hepatocytes. Acta Biochim Pol. 1991, 38: 241-249.

    CAS  PubMed  Google Scholar 

  96. Stangegaard M, Dufva IH, Dufva M: Reverse transcription using random pentadecamer primers increases yield and quality of resulting cDNA. Biotechniques. 2006, 40: 649-657. 10.2144/000112153.

    Article  CAS  PubMed  Google Scholar 

  97. Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, Bunney WE, Myers RM, Speed TP, Akil H, Watson SJ, Meng F: Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005, 33: e175-10.1093/nar/gni179.

    Article  PubMed Central  PubMed  Google Scholar 

  98. Jaksik R, Polanska J, Herok R, Rzeszowska-Wolny J: Calculation of reliable transcript levels of annotated genes on the basis of multiple probe-sets in Affymetrix microarrays. Acta Biochim Pol. 2009, 56: 271-277.

    CAS  PubMed  Google Scholar 

  99. Polanska J, Widlak P, Rzeszowska-Wolny J, Kimmel M, Polanski A: Gaussian mixture decomposition of time-course DNA microarray data, in Mathematical Modeling of Biological Systems. Mathematical Modeling of Biological Systems, Cellular Biophysics, Regulatory Networks, Development, Biomedicine, and Data Analysis. Edited by: Deutsch A, Brusch L, Byrne H, de Vries G, Herzel H-P. 2007, Boston: Birkhäuser Basel, 351-359. [Series: Modeling and Simulation in Science, Engineering and Technology, vol 1]

    Google Scholar 

  100. GeneCards database.,

  101. NCBI database.,

  102. Affymetrix Chromatin Immunoprecipitation Assay Protocol.,

  103. Silesian Bioinformatic Platform.,

  104. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


The authors thank Dr. Hanna Rokita for working protocol of liver perfusion by collagenase, Dr. Piotr Widlak for comments and help with preparation of the manuscript, and Dr Adam Makuchowski for implementing the algoritms for qPCR data analysis. This work was supported by the Polish Ministry of Science and Higher Education (grants numbers N N301 002439 and N N401 031837) and by the European Community from the European Social Fund within the RFSD 2 project and within the INTERKADRA project, no. UDA - POKL-04.01.01-00-014/10-00 (to JK and AT).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Wiesława Widłak.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MK-L carried out most the molecular biology experiments. JP performed bioinformatics analyses of expression and promoter microarrays, interpretation including annotation, GO analysis. JK participated in verification of microarrays data by RT-PCR and qPCR. MO hybridized samples on microarrays. NV participated in the ChIP experiments. AT participated in verification of microarrays data. WW conceived, designed, supervised, carried out analysis and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1: Noise threshold in expression microarrays. Signals are in arbitrary units in log2 scale. Available at: (DOCX 12 KB)


Additional file 2: Data from the global gene expression analyses. Excel document contains all data from expression microarrays studies after normalization according to Dai et al. [97]. Legends for abbreviations used in headings are given in comments fields. Available at: (XLSX 2 MB)


Additional file 3: Table S2: Top ten genes identified in spermatocytes (SC) and hepatocytes (HEP) as the most induced (red bold) by heat shock (at 38°C or 43°C) versus control (C). The level of expression is given in arbitrary units in logarithmic scale (log2). Changes in gene expression are shown as SLR. Available at: (DOCX 27 KB)


Additional file 4: Cluster analysis in spermatocytes subjected to heat shock at 38°C and biological processes linked to genes grouped in clusters. Compressed file contains: Excel document with “cluster analysis_summary” sheet (clusters with genes positively regulated following heat shock are marked in red, with genes negatively regulated – in green), and “cluster analysis_all” sheet with all genes shown separately. Legends for abbreviations used in headings are given in comments fields; a folder with results of Gene to GO BP Conditional test for over-representation analyzed in each cluster (html files). Available at: (RAR 1 MB)


Additional file 5: Cluster analysis in spermatocytes subjected to heat shock at 43°C and biological processes linked to the genes grouped in clusters. File contents – as in file 4. Available at: (RAR 1 MB)


Additional file 6: Cluster analysis in hepatocytes subjected to the heat shock at 43°C and biological processes linked to the genes grouped in clusters. File contents – as in file 4. Available at: (RAR 1 MB)


Additional file 7: Table S3: Top ten genes identified in spermatocytes (SC) and hepatocytes (HEP) as the most repressed (green bold) by heat shock (at 38°C or 43°C) versus control (C). The level of expression is given in arbitrary units in logarithmic scale (log2). Changes in gene expression are shown as SLR. Available at: (DOCX 25 KB)


Additional file 8: Table S4: Genes whose expression was changed in the same direction following hyperthermia in spermatocytes (SC) and in hepatocytes (HEP) (at 38°C or 43°C, respectively). The level of expression is given in arbitrary units in logarithmic scale (log2). Changes in gene expression are shown as SLR (up-regulation: red bold, down-regulation: green bold). Available at: (DOCX 15 KB)


Additional file 9: Figure S1: Detection of transcripts of selected genes in isolated spermatocytes, testes and liver up to 24h after hyperthermia, by RT-PCR. (A) Expression of selected genes connected with cell death shown in the main text in Table 3, (B) transcription factors shown in Table 4, (C) genes involved in inflammatory and immune responses shown in Table 5, (D) two uncharacterized genes induced in spermatocytes at 38°C, and (E) reference genes. N, PCR negative control without template. Available at: (TIFF 1 MB)


Additional file 10: Data from ChIP on chip analyses combined with the results of gene expression analyses. Legends for abbreviations used in headings are given in the comments fields. Available at: (XLSX 6 MB)


Additional file 11: Table S5: Top genes with the uppermost HSF1 binding in spermatocytes and hepatocytes at control temperature and following heat shock Binding of HSF1 is expressed as a mean AB1-AB0 value in arbitrary units. Genes marked in blue in group A are present also in groups B and C, marked in blue in group B are present also in group C. Available at: (DOCX 18 KB)


Additional file 12: Gene ontology analysis of genes bound by HSF1 in spermatocytes and hepatocytes following heat shock. Excel documents contain: -sheet 1: gene to GO terms association; -sheet 2: statistics of biological processes (BP) linked to HSF1-bound genes; -sheet 3: statistics of cellular components (CC) linked to HSF1-bound genes; -sheet 4: statistics of molecular functions (MF) linked to HSF1-bound genes. Terms with the biggest representations were marked in red. Available at: (RAR 1020 KB)


Additional file 13: Correlation of HSF1 binding with changes of gene expression following heat shock. Excel document contains lists of genes bound by HSF1 and: -sheet 1: up-regulated in spermatocytes at 38°C;-sheet 2: down-regulated in spermatocytes at 38°C;-sheet 3: up-regulated in spermatocytes at 43°C;-sheet 4: down-regulated in spermatocytes at 43°C;-sheet 5: up-regulated in hepatocytes at 43°C;-sheet 6: down-regulated in hepatocytes at 43°C. For HSF1 binding values AB1-AB0 ≥ 15 and FDR < 0.125 were taken into consideration (values 0.05 < FDR < 0.125 are shown in grey); in the case of changes in expression level – range of SLR presented in Table 1 in the main text. Legends for abbreviations used in headings are given in the comments fields. Available at: (XLSX 104 KB)


Additional file 14: Table S6: Correlation between hyperthermia-induced changes in gene expression and HSF1 binding to promoters in control (C) or heat-shocked (HS) cells. All genes with expression above the noise threshold and SLR ≠ 0, and with binding of HSF1 where AB1 > AB0 were included into analyzes (N); correlations were estimated using the Spearman’s rank correlation coefficient (rho). Available at: (DOCX 18 KB)


Additional file 15: Table S7: Characteristics of genes whose repression following heat shock in spermatocytes correlates with HSF1 binding. Genes necessary for spermatogenesis listed in the main text in Table 9 are marked here in red. Available at: (DOCX 25 KB)


Additional file 16: Figure S2: Changes of expression of selected genes in spermatocytes following heat shock at 38°C (one hour and two hours of recovery) assessed by quantitative RT-PCR. Values are shown in arbitrary units and are calculated against the level of expression at a physiological temperature which is 1.0. Expression was normalized against the level of Gapdh. *p-value > 0.05. Available at: (TIFF 45 KB)


Additional file 17: Table S8: Characteristics of primers used in RT-PCR analyses. Available at: (DOCX 16 KB)


Additional file 18: Table S9: Characteristics of primers used in ChIP-PCR analyses. Available at: (DOCX 13 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kus-Liśkiewicz, M., Polańska, J., Korfanty, J. et al. Impact of heat shock transcription factor 1 on global gene expression profiles in cells which induce either cytoprotective or pro-apoptotic response following hyperthermia. BMC Genomics 14, 456 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: