- Research article
- Open Access
The genome-wide role of HSF-1 in the regulation of gene expression in Caenorhabditis elegans
© The Author(s). 2016
- Received: 15 February 2016
- Accepted: 15 June 2016
- Published: 5 August 2016
The heat shock response, induced by cytoplasmic proteotoxic stress, is one of the most highly conserved transcriptional responses. This response, driven by the heat shock transcription factor HSF1, restores proteostasis through the induction of molecular chaperones and other genes. In addition to stress-dependent functions, HSF1 has also been implicated in various stress-independent functions. In C. elegans, the HSF1 homolog HSF-1 is an essential protein that is required to mount a stress-dependent response, as well as to coordinate various stress-independent processes including development, metabolism, and the regulation of lifespan. In this work, we have performed RNA-sequencing for C. elegans cultured in the presence and absence of hsf-1 RNAi followed by treatment with or without heat shock. This experimental design thus allows for the determination of both heat shock-dependent and -independent biological targets of HSF-1 on a genome-wide level.
Our results confirm that C. elegans HSF-1 can regulate gene expression in both a stress-dependent and -independent fashion. Almost all genes regulated by HS require HSF-1, reinforcing the central role of this transcription factor in the response to heat stress. As expected, major categories of HSF-1-regulated genes include cytoprotection, development, metabolism, and aging. Within both the heat stress-dependent and -independent gene groups, significant numbers of genes are upregulated as well as downregulated, demonstrating that HSF-1 can both activate and repress gene expression either directly or indirectly. Surprisingly, the cellular process most highly regulated by HSF-1, both with and without heat stress, is cuticle structure. Via network analyses, we identify a nuclear hormone receptor as a common link between genes that are regulated by HSF-1 in a HS-dependent manner, and an epidermal growth factor receptor as a common link between genes that are regulated by HSF-1 in a HS-independent manner. HSF-1 therefore coordinates various physiological processes in C. elegans, and HSF-1 activity may be coordinated across tissues by nuclear hormone receptor and epidermal growth factor receptor signaling.
This work provides genome-wide HSF-1 regulatory networks in C. elegans that are both heat stress-dependent and -independent. We show that HSF-1 is responsible for regulating many genes outside of classical heat stress-responsive genes, including genes involved in development, metabolism, and aging. The findings that a nuclear hormone receptor may coordinate the HS-induced HSF-1 transcriptional response, while an epidermal growth factor receptor may coordinate the HS-independent response, indicate that these factors could promote cell non-autonomous signaling that occurs through HSF-1. Finally, this work highlights the genes involved in cuticle structure as important HSF-1 targets that may play roles in promoting both cytoprotection as well as longevity.
- Heat shock response
- C. elegans
- Transcript analysis
When organisms are exposed to protein-denaturing stressors such as heat, the heat shock response (HSR) is engaged to manage protein damage and restore proteostasis . The HSR is highly conserved across species and is regulated by the transcription factor heat shock factor 1 (HSF1). During basal conditions, HSF1 exists as a monomer in the cytoplasm and nucleus, and during stress conditions undergoes trimerization and accumulation in the nucleus, where it binds to heat shock elements in the promoters of heat shock protein (hsp) genes . HSPs primarily act as molecular chaperones which refold the misfolded proteins that accumulate during stress, but they can also have essential functions in protein synthesis, processing, and degradation [3, 4]. Thus the HSR, and HSPs, play a large role in maintaining organismal proteostasis.
The soil-dwelling, free-living, nematode Caenorhabditis elegans is a powerful model organism that has provided insights into the regulation of a number of stress response pathways, including the HSR. HSF-1, the C. elegans homolog to mammalian HSF1, contains conserved N-terminal DNA-binding and trimerization domains, as well as a putative transactivation domain at the C-terminus . It has recently been shown that the same activity steps required for mammalian HSF1 activation, including trimerization, hyperphosphorylation, and induction of DNA-binding, are also required for worm HSF-1 activation [6, 7].
Studies in C. elegans show that HSF-1 plays a central role not only in the HSR, but also in contributing to organismal physiology. HSF-1 is essential to worm viability, as a truncated mutant that lacks the C-terminal putative activation domain is defective in chaperone induction and egg laying, and also has a decreased lifespan . In addition, this strain has a temperature-sensitive developmental arrest phenotype, with arrest occurring at the L2-L3 transition . Various experiments using hsf-1 RNA interference (RNAi) have shown that HSF-1 regulates the expression of specific hsp genes upon heat shock (HS), and have also implicated a non-stress-induced role for HSF-1 in processes including development, metabolism, and longevity [5, 8–14]. Interestingly, studies in C. elegans have identified the HSR as a cell non-autonomous process that requires thermosensory neurons for hsp induction . Upon the completion of sequencing of the C. elegans genome, over 40 % of the predicted protein products were found to be significantly conserved in other organisms , and many signaling pathways are conserved . C. elegans is thus an excellent model system for studying the role of HSF-1 in stress responses and other physiological processes in a simple multicellular organism.
In this study, we have performed RNA-sequencing (RNA-seq) with synchronous larval stage L4 wild-type C. elegans fed empty vector (EV) control RNAi or hsf-1 RNAi treated with or without HS. We show that significant numbers of genes are upregulated as well as downregulated by HSF-1 under both conditions. In addition to hsp genes, HSF-1 is required for the regulation of genes involved in a wide variety of cellular processes including cytoprotection, development, metabolism, and aging. Network analysis points to possible routes by which HSF-1 signaling may be coordinated across tissues. A nuclear hormone receptor may coordinate the HS-induced HSF-1 transcriptional response, while an epidermal growth factor receptor may coordinate the HS-independent response. Surprisingly, the top HSF-1-regulated gene category, both with and without heat stress, is cuticle structure. This result, together with other recent studies, thus links regulation of the extracellular matrix to HSF-1, cytoprotection, and longevity.
Experimental set-up for genome-wide analysis of regulation of gene expression by HSF-1
Previous experiments have shown that HSF-1 regulates the expression of specific hsp genes upon HS, and have also implicated a non-stress-induced role for HSF-1 in development, metabolism, and longevity [5, 8–14]. To examine HS-dependent vs. -independent gene regulation by HSF-1 on a genome-wide level, we used whole transcriptome RNA-sequencing. We treated synchronous L1 larval stage nematodes with RNAi against hsf-1 [indicated as hsf-1(-)] or with an EV control plasmid [indicated as hsf-1(+)] until the L4 larval stage. At the L4 stage, we then treated nematodes from both groups with or without a 30 min 33 °C HS, as diagrammed (see Additional file 1: Figure S1a). Experiments were performed in biological duplicates. The L4 stage was chosen for our studies as this is a time when the response to HS is strong, prior to a sharp decline that occurs shortly after the transition to adulthood [18, 19]. These treatment conditions, optimized for our studies, resulted in a ~9- log2-fold induction of the hsp-70 gene C12C8.1, a classical HSF-1 target gene, in hsf-1(+) animals treated with HS (see Additional file 1: Figure S1b, black bars). As expected, RNAi against hsf-1 blunted hsp-70 induction by HS (see Additional file 1: Figure S1b, purple bars). The efficiency of our RNAi feeding strategy was assessed by testing the effects of hsf-1 RNAi on transcription driven by two heat shock protein promoter- GFP reporter worm strains (see Additional file 1: Figure S1c-d). HS increases GFP expression, and this effect is dependent on HSF-1 as demonstrated with hsf-1 RNAi. Using an HSF-1::GFP overexpression worm strain, we also show that HSF-1 protein levels are reduced 80 % in response to hsf-1 RNAi treatment (see Additional file 1: Figure S1e-f). Overall, these data validate our HS treatment conditions and RNAi feeding strategy.
Cluster analysis performed on biological replicate RNA-seq samples revealed conserved patterns of expression induced by each treatment condition (see Additional file 1: Figure S2). We normalized each condition to the hsf-1(+);-HS control in order to determine fold changes in relative RNA abundance (see Additional file 1: Figure S3). A complete list of the significant genes altered in response to each condition, after normalization to the control, is also provided (see Additional file 2: Table S1). Volcano plot analyses show that while a limited group of genes for each comparison have a log2-fold change of 6 or higher, the majority of genes have a log2 -fold change of approximately 4 or less (see Additional file 1: Figure S4a-c). As growth temperature and HSF-1 expression levels can affect the rate of development in the worm , we verified that our observed gene expression changes were not simply due to a change in the rate of development between each treatment condition. To do this, we analyzed several genes that are known to be differentially expressed during development and molting, including abu-11, wrt-2, his-24, lin-29, abu-10, abu6, pqn-47, pin-42, ptr-3, abu-8, abu-7, and sdz-37 , and detected no significant expression differences in these genes across treatment groups (see Additional file 1: Figure S5). Overall, our data indicate that the worms in our four treatment conditions are developmentally synchronous to one another, and that the biological replicates for each condition share a similar expression profile, thus validating our experimental conditions.
Next, to visualize total HS-dependent vs. -independent gene expression regulated by HSF-1, we constructed a Venn diagram with the differentially expressed genes for each condition which were determined to be statistically significant as compared to the hsf-1(+);-HS control (see Additional file 1: Figure S6). The shaded areas of the Venn diagram correspond to HS-dependent and -independent processes regulated by HSF-1 (as indicated by the red and pink shaded areas in Figure S6, respectively), and these are the transcripts that we have focused our subsequent analyses on. Altogether, we found that 942 genes are significantly regulated by HSF-1 during HS (Figure S6, red shaded area), and that 2,436 genes are significantly regulated by HSF-1 independently of HS (Figure S6, pink shaded area), highlighting that HSF-1 regulates both HS-dependent and -independent transcriptional processes. Interestingly, only 4 genes are significantly regulated by HS independently of HSF-1. HSF-1 is thus not only an essential transcriptional regulator for a majority of the genes altered by HS, but global RNA expression analysis supports a gene-regulatory role for HSF-1 under both heat stress and non-stress conditions. A greater number of gene changes that depend on HSF-1 are independent of heat stress, thus highlighting the important role of HSF-1 in regulating gene expression under non-stress conditions.
Genes that are regulated by HSF-1 in response to HS
Genes that are normally upregulated by HSF-1 in response to HS
Top 15 genes normally upregulated by HSF-1 in response to HS
Fold Change (log2) hsf-1(+);+HS vs. control
F44E5.5 encodes a member of the Hsp70 family of heat shock proteins
F44E5.4 encodes a member of the Hsp70 family of heat shock proteins
hsp-16.2 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
hsp-16.41 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
hsp-16.11 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
hsp-70 encodes a heat-shock protein that is a member of the HSP70 family of molecular chaperones
hsp-16.1 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
hsp-16.48 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
nspe-1 is a nematode-specific peptide that has an unknown function
hsp-16.49 encodes a 16-kD heat shock protein (HSP) that is a member of the hsp16/hsp20/alphaB-crystallin (HSP16) family of heat shock proteins
nurf-1 encodes the C. elegans ortholog of Drosophila NURF301, a component of the NURF chromatin remodeling complex
scp-1 encodes PTC-related protein that contains a sterol-sensing domain related to human Sterol regulatory element binding protein (SREBP) cleavage activating protein
col-149 is predicted to be a structural constituent of the cuticle
We then further examined the induction characteristics of the top 15 HSF-1-dependent genes induced by HS. The log2-fold changes for a subset of these genes are plotted (see Additional file 1: Figure S7a, black bars) and compared to the expression of the same genes in the presence of HS but in the absence of hsf-1 (see Additional file 1: Figure S7a, purple bars). The fact that hsf-1 RNAi completely eliminates HS-inducibility of these genes highlights their dependency on HSF-1. Independent quantitative RT-PCR (qRT-PCR) for the same subset of highly induced genes (see Additional file 1: Figure S7b) validates our RNA-seq data.
In order to better visualize the global patterns of transcript upregulation by HSF-1 in response to HS, we constructed a heat map of the 654 significantly upregulated genes (Fig. 1b, lane 1). Interestingly, and consistent with the data for a set of the top 15 upregulated genes (see Additional file 1: Figure S7a-b), genes that are normally upregulated by HSF-1 upon HS are either unchanged or downregulated under HS conditions upon hsf-1 knockdown on a global level (Fig. 1b, lane 2). The fact that many of the 654 HS-induced genes that require HSF-1 are downregulated upon hsf-1 knockdown implies that HSF-1 may play a role in the basal regulation of these genes, which is then enhanced upon HS. Together, these data demonstrate the HSF-1 dependency of most HS-induced genes.
To identify the various functional processes normally upregulated by HSF-1 during HS, we used the Database for Annotation, Visualization, and Integrated Discovery (DAVID) classification tool to define the top 5 gene ontology terms for all of the 654 genes found to be significantly upregulated (Fig. 1c). The complete output from DAVID is also provided (see Additional file 4: Table S3). Surprisingly, the top functional category, with an enrichment score of 97, contains genes involved in cuticle structure. The next four categories, all with enrichment scores under 12, include genes involved in translation, the response to stress, the regulation of growth, and amine catabolic processes. We thus find that the largest functional category of genes regulated by HS is not the expected heat stress-responsive gene-set, but instead genes associated with forming cuticle structure.
We next tested the effects of a cuticle collagen gene, col-123, on induction of the HSR by measuring hsp-70 promoter activity in phsp-70::GFP worms (see Additional file 1: Figure S8a-c). We see that col-123 RNAi decreases HS-induced hsp-70 promoter activity, and may control tissue-specific regulation of the HSR. Testing the effects of other collagen genes on regulation of the HSR may provide insight into a signaling role for collagens in coordinating stress responses.
Genes that are normally downregulated by HSF-1 in response to HS
Top 15 genes normally downregulated by HSF-1 in response to HS
Fold Change (log2) hsf-1(+);+HS vs. control
acs-2 encodes an acyl-CoA synthetase; by homology, ACS-2 is predicted to catalyze conversion of a fatty acid to Acyl-CoA for subsequent beta oxidation
dct-1 encodes a protein with similarity to the mammalian BNIP3 proteins that interact with Bcl-2 and the Adenovirus E1B proteins which have been shown to have pro-apoptotic activity
tsp-1 is a part of the tetraspanin family that encodes an ortholog of a human CD151 molecule
fbxa-66 encodes an FboxA protein that has an unknown function
fbxa-21 encodes an FboxA protein that has an unknown function
nep-26 is a thermolysin-like zinc metallopeptidase found on the surface of cells that negatively regulates small signaling peptides
glc-1 is the alpha subunit of a glutamate-gated chloride channel
delm-2 encodes an ortholog of human acid-sensing ion channel family member 4, and is predicted to have sodium channel activity
We then further studied the top 15 genes that are normally downregulated by HSF-1 upon HS. The log2-fold changes of a subset of the top 15 HSF-1-dependent genes repressed by HS are plotted (see Additional file 1: Figure S9a, black bars), and these data are compared to the expression of the same genes in the presence of HS but in the absence of hsf-1 (see Additional file 1: Figure S9a, purple bars). Interestingly, we found that all of the genes that are downregulated by HS in the presence of HSF-1 are upregulated by HS in the absence of HSF-1. One way this could occur is if HSF-1 normally suppresses the expression of these genes, and HS-activated HSF-1 suppresses them even further. To verify our RNA-seq data, we performed independent qRT-PCR for the same subset of highly downregulated genes, and found that the qRT-PCR data was consistent with our RNA-seq data (see Additional file 1: Figure S9b).
To visualize the patterns of transcripts that are normally downregulated by HSF-1 upon HS, we constructed a heat map to visualize the log2-fold changes of the 288 significantly downregulated genes (Fig. 2b, lane 1). As a comparison, the expression of the same transcripts under HS conditions upon hsf-1 knockdown is shown (Fig. 2b, lane 2). As with the data for the top downregulated genes (see Additional file 1: Figure S9a-b), many of the 288 genes that are normally downregulated by HSF-1 under HS conditions are conversely upregulated by HS in the absence of HSF-1. Thus, HS can have completely opposite effects on gene expression depending on the presence or absence of HSF-1.
Upon DAVID analysis of gene ontology terms for the repressed genes, the top 5 functional categories all had enrichment scores of 4 or lower. These functional categories include genes that encode proteins with kinase activity, ion binding activity, transcription, ATP-binding activity, and reproductive development (Fig. 2c). The complete output from DAVID is also provided (see Additional file 4: Table S3). Overall, these results show that a diverse group of genes are normally suppressed by HSF-1 during HS.
Genes that are regulated by HSF-1 independently of HS
Genes that are normally upregulated by HSF-1 independently of HS
Top 15 genes normally upregulated by HSF-1 independently of HS
Fold Change (log2) control vs. hsf-1(-);-HS
acdh-1 encodes a short-chain acyl-CoA dehydrogenase. ACDH-1 is predicted to be a mitochondrial enzyme that catalyzes the first step of fatty acid beta-oxidation, and thus plays a key role in energy production
vit-3 encodes a vitellogenin, a precursor of the lipid-binding protein related to vertebrate vitellogenins and mammalian ApoB-100, a core LDL particle constituent
vit-4 is involved in embryo development and is predicted to have lipid transporter activity
vit-5 encodes a vitellogenin, a lipid-binding protein precursor related to vertebrate vitellogenins and mammalian ApoB-100, a core LDL particle constituent
vit-1 is predicted to have lipid transporter activity
K11G9.3 encodes an ortholog of human butyrylcholinesterase
folt-2 encodes a putative folate transporter and is orthologous to the human folate transporters SLC19A1, SLC19A2, and SLC19A3
ilys-5 is involved in embryo development and is predicted to have lysozyme activity
ZC266.1 is predicted to have G-protein coupled receptor activity
fat-7 encodes an essential delta-9 fatty acid desaturase that is required for the synthesis of monounsaturated fatty acids
K10B2.2 encodes an ortholog of human cathepsin A and is predicted to have carboxypeptidase activity
Y52E8A.4 encodes an ortholog of human major facilitator superfamily domain containing 11
ugt-22 encodes an ortholog of human UDP glucuronosyltransferase 1 family polypeptide, and is predicted to have transferase activity
We then further characterized the induction characteristics of the top 15 genes that are normally upregulated by HSF-1 independently of heat stress. The log2-fold changes of a subset of these top 15 genes are plotted (see Additional file 1: Figure S10a, orange bars), and are compared to the expression of the same genes in the presence of HS and absence of hsf-1 (see Additional file 1: Figure S10a, purple bars). This data shows that HS does not affect the expression of these genes in the absence of hsf-1. The expression of these mRNAs was also verified with qRT-PCR, and the results are consistent with the RNA-seq data (see Additional file 1: Figure S10b).
To investigate a HS-independent role for HSF-1 in the induction of gene expression on a global level, we generated a heat map of all 1,353 genes found via Venn diagram to be normally upregulated by HSF-1 in the absence of HS (Fig. 3b, lane 1). As a comparison, genes that are normally upregulated by HSF-1 in the presence of HS are plotted (Fig. 3b, lane 2). The transcripts in this group that are normally upregulated by HSF-1 in the absence of HS remain upregulated or unchanged in the presence of HS, verifying the heat stress-independent induction of this gene group.
We next used DAVID to determine the top 5 functional processes that are normally upregulated by HSF-1 independently of HS. We found that cuticle structure was again the gene category with the highest enrichment score (9.7), as was also the case for the genes that are upregulated by HSF-1 during HS (Fig. 3c). This indicates that HSF-1 may regulate the basal expression of genes involved in cuticle structure, and that these genes are then further induced upon HS. Other functional processes, with enrichment scores less than 6, include genes that encode proteins involved in the mitochondrial envelope, acyl-CoA dehydrogenase activity, peptidase activity, and oxidoreductase activity. The complete output from DAVID is also available (see Additional file 5: Table S4). The mitochondrial envelope, acyl-CoA dehydrogenase activity, peptidase activity, and oxidoreductase activity are all processes that can be linked to metabolism, further substantiating a functional role for HSF-1 in regulating this process.
Genes that are normally downregulated by HSF-1 independently of HS
Top 15 genes normally downregulated by HSF-1 independently of HS
Fold Change (log2) control vs. hsf-1(-);-HS
T22F3.11 encodes an ortholog of human solute carrier family 17
eol-1 is required in the URX sensory neurons for inhibition of olfactory learning
B0348.2 encodes an ortholog of human lipopolysaccharide-induced TNF factor
col-158 is predicted to be a structural constituent of the cuticle
fbxa-163 encodes a protein containing an F-box, a motif predicted to mediate protein-protein interactions
clec-174 is predicted to have carbohydrate binding activity
srg-31 is involved in embryo development and is predicted to have transmembrane signaling receptor activity
clec-60 appears to play a role in the innate immune response to some bacterial pathogens
B0507.8 encodes an ortholog of human cingulin-like 1
clec-13 is predicted to have carbohydrate binding activity
T08E11.1 encodes a protein containing an F-box motif predicted to mediate protein-protein interactions
F22G12.1 encodes an ortholog of human GRB10 interacting GYF protein 2
The induction characteristics of the top 15 genes in this category were then analyzed. The log2-fold changes from the RNA-seq data for a subset these genes are plotted (see Additional file 1: Figure S10c, orange bars), and are compared to the expression of the same genes in the presence of HS in the absence of hsf-1 (see Additional file 1: Figure S10c, purple bars). The expression of these mRNAs was also verified with qRT-PCR, and the results are consistent with the RNA-seq data (see Additional file 1: Figure S10d). Altogether, these data confirm that HS does not affect the expression of these genes in the absence of hsf-1.
To investigate a HS-independent role for HSF-1 in the suppression of gene expression, we constructed a heat map of all 1,083 genes found via Venn diagram to be normally suppressed by hsf-1 in the absence of HS (Fig. 4b, lane 1). As a comparison, the expression patterns of the same transcripts suppressed by hsf-1 in the presence of HS are shown (Fig. 4b, lane 2). The transcripts in this group that are significantly downregulated by HSF-1 in the absence of HS remain downregulated or unchanged by HSF-1 in the presence of HS, further verifying that regulation of this subset of genes by HSF-1 is independent of HS.
We next used DAVID to determine the functional processes that are normally downregulated by HSF-1 independently of HS. Genes involved in cell cycle processes were most abundant, with an enrichment score of 19, followed by genes involved in epithelium development, regulation of growth, genitalia development, and cell migration, all with enrichment scores under 6 (Fig. 4c). The complete output from DAVID is also available (see Additional file 5: Table S4). Cell cycle processes, epithelium development, regulation of growth, genitalia development, and cell migration are all processes that can be linked to development, thus confirming a HS-independent role for HSF-1 in development.
Regulation of gene expression by HSF-1
C. elegans is a useful model organism for identifying regulatory processes that are shared between species. While the transcription factor HSF-1 has classically been studied as a factor that is responsive to cytoplasmic proteotoxic stress, it is becoming increasingly evident that this transcription factor also has major non-stress-induced roles in coordinating gene expression. With the RNA-sequencing experiments performed here, we confirm that HSF-1 can regulate gene expression under both heat stress and non-stress conditions. In addition to genes that are classically stress-responsive, such as chaperones, our work shows that HSF-1 also regulates sets of genes involved in a variety of cellular processes including metabolism, development, and longevity.
Cuticle structure genes are normally upregulated by HSF-1 via HS-dependent and -independent mechanisms
A surprising result of our study is that genes controlling cuticle structure comprise the largest gene ontology group that is upregulated by HSF-1 in both a HS-dependent and -independent manner. The C. elegans cuticle is an exoskeletal structure that creates a barrier between the animal and the environment, provides body shape, and allows movement via attachment to muscle. Many of the genes in the cuticle structure category are collagens, structural proteins that form an extracellular matrix composing the exoskeleton, or cuticle, of the nematode. There are ~154 distinct collagen genes in C. elegans, and they are expressed in a tissue-specific fashion and at distinct temporal times . Cuticle structure is controlled by enzymes involved in collagen processing, and the polymerization pattern is dictated by actin filaments that are organized in specific patterns around the body of the worm. In humans, collagens comprise about one-third of all expressed protein . Aside from the structural role of collagens, these genes can also participate in signal transduction [37–39]. Collagen genes were recently found to be upregulated by SKN-1, the C. elegans oxidative stress-responsive transcription factor . In future work, it will be interesting to test whether collagen can act to relay signals to stress-specific transcription factors including SKN-1 and HSF-1.
Roles for HSF-1 in regulating metabolism and development in a HS-independent manner
Although HSF-1 has classically been studied as a transcription factor that responds to HS and cytoplasmic proteotoxic stress, HSF-1 has recently been gaining importance as a transcription factor that is involved in non-stress processes including development and metabolism [41, 42]. In mice, HSF family members have been documented to be involved in diverse developmental processes, including oogenesis, spermatogenesis, and corticogenesis . HSF-1 in the worm has also been shown to be regulated by insulin/IGF-1, TGF-β, and cGMP signaling to control development . The finding that mammalian HSF1 can be regulated by SIRT1, a deacetylase that is under metabolic control, provides evidence that HSF-1 and metabolism are linked . This finding is also true for C. elegans HSF-1, as the SIRT1 homolog SIR-2.1 regulates the C. elegans HSR . Additionally, the insulin-like signaling regulators DDL-1/2 have been linked to HSF-1 regulation . Here, we show that under non-stress conditions, HSF-1 regulates a number of genes involved in developmental and metabolic processes. Therefore, our work further highlights the links between HSF-1 and these non-stress processes.
Network analysis identifies a nuclear hormone receptor as a common link between processes regulated by HSF-1 upon HS
To determine how the genes regulated by HSF-1 during HS may interact with each other, we performed network analysis using genes associated with the top GO-terms as determined by DAVID (Fig. 5). We used the MiMI plugin to integrate data from protein interaction databases (including gene ontology databases, MeSH, and PubMed) to allow for the creation of interaction networks using the network-building software Cytoscape. This analysis enabled us to identify interacting partners shared by at least two genes regulated by HSF-1 during HS. The transcripts induced by HSF-1 during HS are shown in red, while the transcripts downregulated by HSF-1 during HS are shown in blue, with the intensity of color correlating to the fold change. Genes that are not colored were not found to be affected by HSF-1 during HS in our dataset, but are neighbors shared by at least two genes that were affected.
Network analysis identifies a tyrosine kinase as a common link between various developmental processes regulated by HSF-1 independently of HS
To uncover C. elegans interaction networks associated with processes regulated by HSF-1 independently of HS, we performed network analysis using genes associated with the top GO-terms as determined by DAVID (Fig. 5b). Network linkages were found for genes involved in the mitochondrial envelope, peptidase activity, oxidoreductase activity, cell cycle processes, and epithelium development. We found that let-23, an epidermal growth factor receptor tyrosine kinase (EGF-RTK) , is predicted to interact with many of these transcripts. The EGF pathway in C. elegans has been linked to multiple developmental pathways . This RTK may thus allow for signaling to HSF-1 during non-stress conditions to modulate developmental gene expression. Interestingly, HSF1 null mouse embryonic fibroblasts are defective in both basal and EGF-induced cell migration , so the link between HSF-1 and the EGF signaling pathways may be conserved across species. It will thus be worthwhile in future work to test for the involvement of let-23 in the regulation of HSF-1 activity in C. elegans.
HSF-1 impacts aging-regulated gene expression
Next, a comparison between aging-associated genes and those we found to be regulated by HSF-1 independently of HS was done. A Venn diagram shows that 275 aging-associated genes overlapped with our dataset (Fig. 6d). A complete list of genes shared between both data-sets is also provided (see Additional file 7: Table S6). The functional processes regulated by this gene set was then determined via DAVID analysis. Cuticle structure was again the largest category, with an enrichment score of 9. The other categories, with enrichment scores of 5 or lower, include peptidase activity, C-type lectin, vitamin binding, and lifespan-associated processes (Fig. 6e). The transcriptional impact of these HSF-1 and aging-regulated genes that are independent of HS was then determined by network generation (Fig. 6f). Interestingly, the interaction network of aging-associated genes regulated by HSF-1 independently of HS is four-fold larger than the network of aging genes that depend on HS. We thus conclude that HSF-1 may have a role in impacting longevity that can be separated from its role in stress responses.
HSF-1 regulates collagen genes which may affect the aging process
It is interesting that cuticle structure genes constitute the largest overlap with aging-related genes. In humans, mutations in collagens lead to a large number of heritable human diseases such as osteoporosis and musculoskeletal diseases . Collagens are long-lived proteins known to accumulate damage during aging, leading to a decline in tissue health . Also, type I collagens become resistant to proteolysis upon age [55, 56], affecting their turnover. Interestingly, mice expressing cleavage-resistant type I collagen go through an accelerated aging process . Thus, cellular aging can be affected by the state of the extracellular matrix in mammals.
Recently, collagen production and extracellular matrix remodeling were determined to be essential for longevity in C. elegans. Collagen may directly affect signaling processes associated with longevity in C. elegans, including signaling through SKN-1 [40, 58]. We note that HSF-1 was also recently shown to regulate cytoskeletal integrity in a process that can influence stress resistance and longevity in C. elegans . Thus, the linkage of both the extracellular matrix and the cytoskeleton to HSF-1 may provide a mechanism by which HSF-1 promotes longevity.
Next generation sequencing has allowed us to uncover highly varied roles for C. elegans HSF-1 in both HS-dependent and -independent mechanisms, including roles in the regulation of development, cytoprotection, metabolism, and aging (for a model, see Additional file 1: Figure S11). Network analyses show that under HS conditions, the nuclear hormone receptor NHR-111 may allow coordination of the HSF-1 response across tissues, while under basal conditions, the EGF receptor LET-23 may regulate a similar coordination. These findings warrant further studies in order to further understand the methods of cell non-autonomous signaling across tissues. A striking result shown here is that multiple genes involved in cuticle structure, including collagen genes, are enriched as HSF-1 targets in HS-dependent and -independent manners. As recent studies link collagen to cytoprotection and longevity, the regulation of collagen expression may be one method by which HSF-1 enhances lifespan. Harnessing the ability of HSF-1 to regulate collagen could thus have broad appeal in the treatment of diseases of aging.
C. elegans strains and maintenance
The wild-type N2 strain, phsp-70(C12C8.1)::GFP , phsp-16.2(Y46H3A.3)::GFP , and EQ73 (HSF-1::GFP)  strains were used in this study. Worms were maintained at 23 °C on standard NGM plates seeded with Escherichia coli OP50 [61–63]. A synchronous population of nematodes was obtained by standard 20 % hypochlorite treatment, and a 24 h rotation at 220 rpm in M9 buffer without food.
RNA interference and heat shock conditions
Approximately 4,000 wild-type nematodes were synchronized and placed at the L1 larval stage onto standard NGM plates supplemented with 50 μg/mL ampicillin and 1 mM isopropyl-beta- -thiogalactopyranoside seeded with either HT115 bacteria containing an empty plasmid (L4440, control), or sequence-verified gene-specific RNAi isolated from the Ahringer RNAi library . RNAi bacteria were allowed to induce on the plates overnight at room-temperature. Synchronized animals developed on RNAi plates before being treated at the L4 stage with a 30 min 33 °C HS by submerging the plates into a water bath. The time and duration of HS was optimized for this experiment (see Additional file 1: Figure S1b-d). Worms were then collected for RNA extraction.
Immunoblotting and quantification
Animals were harvested in Buffer C (20 mM HEPES pH 7.9, 25 % Glycerol, 0.42 M NaCl, 1.5 mM MgCl2, 0.2 mM EDTA, and 0.5 mM DTT) with the addition of Halt™ protease inhibitors (Pierce, cat# 78430). Protein was extracted by sonication with a Diagenode Bioruptor 300 for 10 min with 30 s pulses. Protein was quantified by Bradford assay, resolved on a 10 % SDS-PAGE gel, and transferred to a PVDF membrane. The blot was incubated with an α-GFP polyclonal antibody (Abcam, cat# ab290) at a 1:2500 dilution and with α-Actin (Amersham, cat# JLA20-C) at a 1:750 dilution. Quantification of band intensity was performed using ImageJ Software (v. 1.44; http://imagej.nih.gov/ij/).
RNA preparation for RNA-seq
Total RNA was prepared using TRIzol® reagent (Ambion®, cat# 15596-026) by standard protocols, and then cleaned up on RNeasy columns (QIAgen, cat# 74104). RNA integrity analysis, sample preparation, and RNA-sequencing was performed at the Yale Center for Genome Analysis using the Illumina HiSeq 2000 sequencing system.
RNA-seq data analysis
A quality-control analysis of raw RNA-seq reads was performed using the FastQC program . Short reads were aligned to the C. elegans reference genome (ws200 release) using Bowtie software . The program TopHat was used to discover transcript splicing junctions . The program Cufflinks was chosen to assemble the aligned reads, estimate their abundance, and calculate the fragments per kilobase of exon per million fragments mapped (FPKM) values . Transcripts that were differentially expressed in different conditions, compared to the hsf-1(+);-HS control, were determined with CuffDiff, which uses the Benjamini-Hochberg correction for multiple testing to obtain the q-value (the FDR-adjusted the p-value) . The results were visualized with a dendogram using the program CummeRbund . The RNA-seq data has been deposited in NCBI SRA database (Access ID: SRP078295).
Volcano plot analysis
Volcano plots were made using GraphPad Prism Software (GraphPad Software, La Jolla California USA, http://www.graphpad.com), where the Y-axis represents the q-value (FDR-corrected p-value) after being adjusted to reflect a -log10 value, and the X-axis represents the log2-fold change of each mRNA after comparison to the hsf-1(+);-HS control.
Venn diagram analysis
Venny 2.0  was used to construct Venn diagrams with the significantly altered mRNAs for each condition (q-value < 0.05) as compared to the hsf-1(+);-HS control.
qRT-PCR was performed to validate the top hits from our RNA-seq data. An aliquot of the RNA samples that were used for sequencing were reverse transcribed into cDNA using a High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, cat# 4368814) according to the manufacturer’s instructions. cDNA was diluted to 50 ng/μl to be used as a template for qRT-PCR which was performed with the Step One Plus Real-time PCR system (Applied Biosystems) using iTaq™ Universal SYBR® Green Supermix (BioRad, cat# 172-5121) according to manufacturer’s instructions. Data analysis was performed according to standard calculations using the comparative Ct method . Relative mRNA levels were normalized to gapdh, and calculated from two biological replicates and technical triplicates. Primer sequences are available upon request. Statistical analyses were carried out with GraphPad Prism Software (GraphPad Software, La Jolla California USA, http://www.graphpad.com) using ANOVA followed by the Bonferroni post-hoc test. Error bars are representative of standard deviation between independent biological replicates.
Animals were anesthetized with 10 mM Levamisole and photographed using an EVOS fluorescence microscope. Image processing was accomplished using Adobe Photoshop© (Adobe Systems Incorporated, San Jose, CA).
Heat map generation
The heat maps were organized using Cluster 3, by organizing the genes into 3 clusters, using K-means and 100 runs, and the Euclidean distance similarity metric .
Gene ontology analysis via DAVID
The Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used to identify over-represented gene ontology terms using the Functional Annotation Clustering tool and a high classification stringency . The enrichment score provided by DAVID takes into account the probability that the members of a gene cluster are present randomly in the gene list. The enrichment score determines biologically significant functional groups by using the p-values for a cluster of genes to determine the geometric mean of that cluster (in negative log scale), where if the geometric mean of the p-values = 1e−10, then the enrichment score would be 10.
Network analysis with the Cytoscape platform
Network analysis was performed using the MiMI plugin for the Cytoscape platform . The MiMI plugin integrates data from protein interaction databases including gene ontology databases, MeSH, and PubMed to allow the creation of interaction networks using the network-building software Cytoscape. Interacting partners shared by at least two mRNAs were identified and used to construct interaction pathways.
DAVID, Database for Annotation, Visualization, and Integrated Discovery; EGF-RTK, Epidermal Growth Factor Receptor Tyrosine Kinase; EV, Empty Vector; HS, Heat shock; HSF-1, Heat Shock Factor 1; HSP(s), Heat Shock Protein(s); HSR, Heat Shock Response; qRT-PCR, quantitative RT-PCR; RNAi, RNA interference; RNA-seq, RNA-sequencing.
The authors would like to thank Dr. R.I. Morimoto for providing the phsp-70(C12C8.1)::GFP strain, Dr. A. Hsu for providing the HSF-1::GFP (EQ73) strain, and the Caenorhabditis Genetics Center, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440), for providing the N2 strain. We would also like to thank Dr. Keith Choe, Dr. Nancy Denslow, Dr. Chris Martyniuk, and Dr. Rachel Raynes for helpful discussions and advice, and the Yale Center for Genome Analysis for performing the RNA-sequencing analysis.
This work was supported by a departmental start-up grant from the Department of Cell Biology, Microbiology, and Molecular Biology at the University of South Florida to Sandy D. Westerheide.
Availability of data and materials
The data-sets supporting the results of this article are available in the NCBI BioProject repository (Accession Number PRJNA311958, http://www.ncbi.nlm.nih.gov/bioproject/PRJNA311958) and in the SRA database (Accession Number SRP078295).
JB carried out the sequencing experiment, data analysis, and Fig. generation. SM carried out qRT-PCR experiments. YL and FC contributed to data analysis, normalization of the sequencing data, and dendogram clustering of samples. JB, FC, and SW participated in the design of the study. JB and SW wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gidalevitz T, Prahlad V, Morimoto RI: The stress of protein misfolding: from single cells to multicellular organisms. Cold Spring Harb Perspect Biol 2011, 3(6). doi: 10.1101/cshperspect.a009704.
- Anckar J, Sistonen L. Heat shock factor 1 as a coordinator of stress and developmental pathways. Adv Exp Med Biol. 2007;594:78–88.View ArticlePubMedGoogle Scholar
- Parsell DA, Lindquist S. The function of heat-shock proteins in stress tolerance: degradation and reactivation of damaged proteins. Annu Rev Genet. 1993;27:437–96.View ArticlePubMedGoogle Scholar
- Young JC, Agashe VR, Siegers K, Hartl FU. Pathways of chaperone-mediated protein folding in the cytosol. Nat Rev Mol Cell Biol. 2004;5(10):781–91.View ArticlePubMedGoogle Scholar
- Hajdu-Cronin YM, Chen WJ, Sternberg PW. The L-type cyclin CYL-1 and the heat-shock-factor HSF-1 are required for heat-shock-induced protein expression in Caenorhabditis elegans. Genetics. 2004;168(4):1937–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Chiang WC, Ching TT, Lee HC, Mousigian C, Hsu AL. HSF-1 regulators DDL-1/2 link insulin-like signaling to heat-shock responses and modulation of longevity. Cell. 2012;148(1-2):322–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Morton EA, Lamitina T. Caenorhabditis elegans HSF-1 is an essential nuclear protein that forms stress granule-like structures following heat shock. Aging Cell. 2013;12(1):112–20.View ArticlePubMedGoogle Scholar
- Garigan D, Hsu AL, Fraser AG, Kamath RS, Ahringer J, Kenyon C. Genetic analysis of tissue aging in Caenorhabditis elegans: a role for heat-shock factor and bacterial proliferation. Genetics. 2002;161(3):1101–12.PubMedPubMed CentralGoogle Scholar
- Hsu AL, Murphy CT, Kenyon C. Regulation of aging and age-related disease by DAF-16 and heat-shock factor. Science. 2003;300(5622):1142–5.View ArticlePubMedGoogle Scholar
- Morley JF, Morimoto RI. Regulation of longevity in Caenorhabditis elegans by heat shock factor and molecular chaperones. Mol Biol Cell. 2004;15(2):657–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Walker GA, Thompson FJ, Brawley A, Scanlon T, Devaney E. Heat shock factor functions at the convergence of the stress response and developmental pathways in Caenorhabditis elegans. FASEB J. 2003;17(13):1960–2.PubMedGoogle Scholar
- Minsky N, Roeder RG. Direct link between metabolic regulation and the heat-shock response through the transcriptional regulator PGC-1alpha. Proc Natl Acad Sci U S A. 2015;112(42):E5669–5678.View ArticlePubMedPubMed CentralGoogle Scholar
- Peng ZY, Serkova NJ, Kominsky DJ, Brown JL, Wischmeyer PE. Glutamine-mediated attenuation of cellular metabolic dysfunction and cell death after injury is dependent on heat shock factor-1 expression. JPEN J Parenter Enteral Nutr. 2006;30(5):373–8. discussion 379.View ArticlePubMedGoogle Scholar
- Westerheide SD, Anckar J, Stevens Jr SM, Sistonen L, Morimoto RI. Stress-inducible regulation of heat shock factor 1 by the deacetylase SIRT1. Science. 2009;323(5917):1063–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Prahlad V, Cornelius T, Morimoto RI. Regulation of the cellular heat shock response in Caenorhabditis elegans by thermosensory neurons. Science. 2008;sas320(5877):811–4.View ArticleGoogle Scholar
- Consortium CeS. Genome sequence of the nematode C. elegans: a platform for investigating biology. Science. 1998;282(5396):2012–8.View ArticleGoogle Scholar
- Kaletta T, Hengartner MO. Finding function in novel targets: C. elegans as a model organism. Nat Rev Drug Discov. 2006;5(5):387–98.View ArticlePubMedGoogle Scholar
- Ben-Zvi A, Miller EA, Morimoto RI. Collapse of proteostasis represents an early molecular event in Caenorhabditis elegans aging. Proc Natl Acad Sci U S A. 2009;106(35):14914–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Labbadia J, Morimoto RI. Repression of the Heat Shock Response Is a Programmed Event at the Onset of Reproduction. Mol Cell. 2015;59(4):639–50.View ArticlePubMedGoogle Scholar
- Klass MR. Aging in the nematode Caenorhabditis elegans: major biological and environmental factors influencing life span. Mech Ageing Dev. 1977;6(6):413–29.View ArticlePubMedGoogle Scholar
- Hendriks GJ, Gaidatzis D, Aeschimann F, Grosshans H. Extensive oscillatory gene expression during C. elegans larval development. Mol Cell. 2014;53(3):380–92.View ArticlePubMedGoogle Scholar
- Andersen EC, Lu X, Horvitz HR. C. elegans ISWI and NURF301 antagonize an Rb-like pathway in the determination of multiple cell fates. Development. 2006;133(14):2695–704. View ArticlePubMedGoogle Scholar
- Van Gilst MR, Hadjivassiliou H, Jolly A, Yamamoto KR. Nuclear hormone receptor NHR-49 controls fat consumption and fatty acid composition in C. elegans. PLoS Biol. 2005;3(2):e53.View ArticlePubMedPubMed CentralGoogle Scholar
- Pinkston-Gosse J, Kenyon C. DAF-16/FOXO targets genes that regulate tumor growth in Caenorhabditis elegans. Nat Genet. 2007;39(11):1403–9.View ArticlePubMedGoogle Scholar
- Moribe H, Yochem J, Yamada H, Tabuse Y, Fujimoto T, Mekada E. Tetraspanin protein (TSP-15) is required for epidermal integrity in Caenorhabditis elegans. J Cell Sci. 2004;117(Pt 22):5209–20.View ArticlePubMedGoogle Scholar
- Turner AJ, Isaac RE, Coates D. The neprilysin (NEP) family of zinc metalloendopeptidases: genomics and function. Bioessays. 2001;23(3):261–9.View ArticlePubMedGoogle Scholar
- Ghosh R, Andersen EC, Shapiro JA, Gerke JP, Kruglyak L. Natural variation in a chloride channel subunit confers avermectin resistance in C. elegans. Science. 2012;335(6068):574–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Han L, Wang Y, Sangaletti R, D'Urso G, Lu Y, Shaham S, Bianchi L. Two novel DEG/ENaC channel subunits expressed in glia are needed for nose-touch sensitivity in Caenorhabditis elegans. J Neurosci. 2013;33(3):936–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Spieth J, Denison K, Kirtland S, Cane J, Blumenthal T. The C. elegans vitellogenin genes: short sequence repeats in the promoter regions and homology to the vertebrate genes. Nucleic Acids Res. 1985;13(14):5283–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Murphy CT, McCarroll SA, Bargmann CI, Fraser A, Kamath RS, Ahringer J, Li H, Kenyon C. Genes that act downstream of DAF-16 to influence the lifespan of Caenorhabditis elegans. Nature. 2003;424(6946):277–83.View ArticlePubMedGoogle Scholar
- Van Gilst MR, Hadjivassiliou H, Yamamoto KR. A Caenorhabditis elegans nutrient response system partially dependent on nuclear receptor NHR-49. Proc Natl Acad Sci U S A. 2005;102(38):13496–501.View ArticlePubMedPubMed CentralGoogle Scholar
- Balamurugan K, Ashokkumar B, Moussaif M, Sze JY, Said HM. Cloning and functional characterization of a folate transporter from the nematode Caenorhabditis elegans. Am J Physiol Cell Physiol. 2007;293(2):C670–681.View ArticlePubMedGoogle Scholar
- Watts JL, Browse J. A palmitoyl-CoA-specific delta9 fatty acid desaturase from Caenorhabditis elegans. Biochem Biophys Res Commun. 2000;272(1):263–9.View ArticlePubMedGoogle Scholar
- Shen Y, Zhang J, Calarco JA, Zhang Y. EOL-1, the homolog of the mammalian Dom3Z, regulates olfactory learning in C. elegans. J Neurosci. 2014;34(40):13364–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnstone IL. Cuticle collagen genes. Expression in Caenorhabditis elegans. Trends Genet. 2000;16(1):21–7.View ArticlePubMedGoogle Scholar
- Di Lullo GA, Sweeney SM, Korkko J, Ala-Kokko L, San Antonio JD. Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen. J Biol Chem. 2002;277(6):4223–31.View ArticlePubMedGoogle Scholar
- Fu HL, Valiathan RR, Arkwright R, Sohail A, Mihai C, Kumarasiri M, Mahasenan KV, Mobashery S, Huang P, Agarwal G et al. Discoidin domain receptors: unique receptor tyrosine kinases in collagen-mediated signaling. J Biol Chem. 2013;288(11):7430–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Munger JS, Sheppard D. Cross talk among TGF-beta signaling pathways, integrins, and the extracellular matrix. Cold Spring Harb Perspect Biol. 2011;3(11):a005017.View ArticlePubMedPubMed CentralGoogle Scholar
- Seeger-Nukpezah T, Golemis EA. The extracellular matrix and ciliary signaling. Curr Opin Cell Biol. 2012;24(5):652–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewald CY, Landis JN, Abate JP, Murphy CT, Blackwell TK: Dauer-independent insulin/IGF-1-signalling implicates collagen remodelling in longevity. Nature 2014.Google Scholar
- Seo K, Choi E, Lee D, Jeong DE, Jang SK, Lee SJ. Heat shock factor 1 mediates the longevity conferred by inhibition of TOR and insulin/IGF-1 signaling pathways in C. elegans. Aging Cell. 2013;12(6):1073–81.View ArticlePubMedGoogle Scholar
- Steinkraus KA, Smith ED, Davis C, Carr D, Pendergrass WR, Sutphin GL, Kennedy BK, Kaeberlein M. Dietary restriction suppresses proteotoxicity and enhances longevity by an hsf-1-dependent mechanism in Caenorhabditis elegans. Aging Cell. 2008;7(3):394–404.View ArticlePubMedPubMed CentralGoogle Scholar
- Akerfelt M, Morimoto RI, Sistonen L. Heat shock factors: integrators of cell stress, development and lifespan. Nat Rev Mol Cell Biol. 2010;11(8):545–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Barna J, Princz A, Kosztelnik M, Hargitai B, Takacs-Vellai K, Vellai T. Heat shock factor-1 intertwines insulin/IGF-1, TGF-beta and cGMP signaling to control development and aging. BMC Dev Biol. 2012;12:32.View ArticlePubMedPubMed CentralGoogle Scholar
- Raynes R, Leckey Jr BD, Nguyen K, Westerheide SD. Heat Shock and Caloric Restriction have a Synergistic Effect on the Heat Shock Response in a sir2.1-dependent Manner in Caenorhabditis elegans. J Biol Chem. 2012;287(234):29045–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Taubert S, Ward JD, Yamamoto KR. Nuclear hormone receptors in nematodes: evolution and function. Mol Cell Endocrinol. 2011;334(1-2):49–55.View ArticlePubMedGoogle Scholar
- Weber KP, Alvaro CG, Baer GM, Reinert K, Cheng G, Clever S, Wightman B. Analysis of C. elegans NR2E nuclear receptors defines three conserved clades and ligand-independent functions. BMC Evol Biol. 2012;12:81.View ArticlePubMedPubMed CentralGoogle Scholar
- Aroian RV, Koga M, Mendel JE, Ohshima Y, Sternberg PW. The let-23 gene necessary for Caenorhabditis elegans vulval induction encodes a tyrosine kinase of the EGF receptor subfamily. Nature. 1990;348(6303):693–9.View ArticlePubMedGoogle Scholar
- Yu S, Driscoll M. EGF signaling comes of age: promotion of healthy aging in C. elegans. Exp Gerontol. 2011;46(2-3):129–34.View ArticlePubMedGoogle Scholar
- O'Callaghan-Sunol C, Sherman MY. Heat shock transcription factor (HSF1) plays a critical role in cell migration via maintaining MAP kinase signaling. Cell Cycle. 2006;5(13):1431–7.View ArticlePubMedGoogle Scholar
- Anckar J, Sistonen L. Regulation of HSF1 function in the heat stress response: implications in aging and disease. Annu Rev Biochem. 2011;80:1089–115.View ArticlePubMedGoogle Scholar
- Budovskaya YV, Wu K, Southworth LK, Jiang M, Tedesco P, Johnson TE, Kim SK. An elt-3/elt-5/elt-6 GATA transcription circuit guides aging in C. elegans. Cell. 2008;134(2):291–303.View ArticlePubMedPubMed CentralGoogle Scholar
- Myllyharju J, Kivirikko KI. Collagens, modifying enzymes and their mutations in humans, flies and worms. Trends Genet. 2004;20(1):33–43.View ArticlePubMedGoogle Scholar
- Toyama BH, Hetzer MW. Protein homeostasis: live long, won't prosper. Nat Rev Mol Cell Biol. 2013;14(1):55–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Hamlin CR, Kohn RR. Evidence for progressive, age-related structural changes in post-mature human collagen. Biochim Biophys Acta. 1971;236(2):458–67.View ArticlePubMedGoogle Scholar
- Hamlin CR, Kohn RR. Determination of human chronological age by study of a collagen sample. Exp Gerontol. 1972;7(6):377–9.View ArticlePubMedGoogle Scholar
- Vafaie F, Yin H, O'Neil C, Nong Z, Watson A, Arpino JM, Chu MW, Wayne Holdsworth D, Gros R, Pickering JG. Collagenase-resistant collagen promotes mouse aging and vascular cell senescence. Aging Cell. 2014;13(1):121–30.View ArticlePubMedGoogle Scholar
- Shin H, Lee H, Fejes AP, Baillie DL, Koo HS, Jones SJ. Gene expression profiling of oxidative stress response of C. elegans aging defective AMPK mutants using massively parallel transcriptome sequencing. BMC Res Notes. 2011;4:34.View ArticlePubMedPubMed CentralGoogle Scholar
- Baird NA, Douglas PM, Simic MS, Grant AR, Moresco JJ, Wolff SC, Yates JR, 3rd, Manning G, Dillin A. HSF-1-mediated cytoskeletal integrity determines thermotolerance and life span. Science. 2014;346(6207):360–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Link CD, Cypser JR, Johnson CJ, Johnson TE. Direct observation of stress response in Caenorhabditis elegans using a reporter transgene. Cell Stress Chaperones. 1999;4(4):235–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77(1):71–94.PubMedPubMed CentralGoogle Scholar
- Stiernagle T: Maintenance of C. elegans. WormBook : the online review of C elegans biology 2006:1-11.Google Scholar
- Hirsh D, Oppenheim D, Klass M. Development of the reproductive system of Caenorhabditis elegans. Dev Biol. 1976;49(1):200–19.View ArticlePubMedGoogle Scholar
- Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M et al. Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003;421(6920):231–7.View ArticlePubMedGoogle Scholar
- FastQC: a quality control tool for high throughput sequence data [http://www.bioinformatics.babraham.ac.uk/projects/fastqc]
- Langmead B: Aligning short sequencing reads with Bowtie. Current protocols in bioinformatics / editoral board, Andreas D Baxevanis [et al] 2010, Chapter 11:Unit 11 17.Google Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25(9):1105–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.View ArticlePubMedPubMed CentralGoogle Scholar
- An interactive tool for comparing lists with Venn's diagrams [http://bioinfogp.cnb.csic.es/tools/venny/index.html]
- Bookout AL, Mangelsdorf DJ. Quantitative real-time PCR protocol for analysis of nuclear receptor signaling pathways. Nucl Recept Signal. 2003;1:e012.View ArticlePubMedPubMed CentralGoogle Scholar
- de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20(9):1453–4.View ArticlePubMedGoogle Scholar
- Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4(5):3.View ArticleGoogle Scholar
- Gao J, Ade AS, Tarcea VG, Weymouth TE, Mirel BR, Jagadish HV, States DJ. Integrating and annotating the interactome using the MiMI plugin for cytoscape. Bioinformatics. 2009;25(1):137–8.View ArticlePubMedGoogle Scholar