Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2011; 6(9): e25476.
Published online 2011 September 29. doi:  10.1371/journal.pone.0025476
PMCID: PMC3183046

HIF-1 and SKN-1 Coordinate the Transcriptional Response to Hydrogen Sulfide in Caenorhabditis elegans

Todd Lamitina, Editor


Hydrogen sulfide (H2S) has dramatic physiological effects on animals that are associated with improved survival. C. elegans grown in H2S are long-lived and thermotolerant. To identify mechanisms by which adaptation to H2S effects physiological functions, we have measured transcriptional responses to H2S exposure. Using microarray analysis we observe rapid changes in the abundance of specific mRNAs. The number and magnitude of transcriptional changes increased with the duration of H2S exposure. Functional annotation suggests that genes associated with protein homeostasis are upregulated upon prolonged exposure to H2S. Previous work has shown that the hypoxia-inducible transcription factor, HIF-1, is required for survival in H2S. In fact, we show that hif-1 is required for most, if not all, early transcriptional changes in H2S. Moreover, our data demonstrate that SKN-1, the C. elegans homologue of NRF2, also contributes to H2S-dependent changes in transcription. We show that these results are functionally important, as skn-1 is essential to survive exposure to H2S. Our results suggest a model in which HIF-1 and SKN-1 coordinate a broad transcriptional response to H2S that culminates in a global reorganization of protein homeostasis networks.


Exogenous H2S has dramatic effects on mammalian physiology that can improve survival in changing environmental conditions. Mice exposed to H2S enter into a hibernation-like state that allows them to endure periods of low metabolic rate and decreased core body temperature without apparent ill effects [1]. The H2S-induced state enables mice to survive exposure to otherwise lethal hypoxic conditions [2], and improves outcome in rodent models of severe blood loss [3], myocardial infarction [4], aortic occlusion [5] and hepatic ischemia/reperfusion [6].

C. elegans grown in H2S have increased thermotolerance and lifespan [7]. Increased lifespan and thermotolerance require the conserved sirtuin homologue sir-2.1, though mutant animals with deletions in sir-2.1 grow normally in H2S. In contrast, the hif-1 transcription factor is required to survive exposure to H2S [8]. hif-1 is a highly conserved bHLH transcription factor that is well-known for its role in coordinating the transcriptional response to hypoxia in all animals, including C. elegans [9], [10]. Sirtuins have been shown to modulate lifespan in yeast, worms, flies and mice [11]. Recent work has demonstrated that HIF-1 activity can influence lifespan in C. elegans [12], [13], [14]. Thus, the response to H2S involves at least two genes, hif-1 and sir-2.1, which influence lifespan.

In this study, we investigated the transcriptional response to H2S in C. elegans. Using an unbiased microarray approach, we show that there are rapid and progressive changes in mRNA abundance associated with exposure to H2S. Functional genomic analysis suggests that adaptation to H2S results in significant changes to protein homeostasis pathways. We found that hif-1 is required for nearly all of the early changes observed, though there is little overlap between genes regulated in response to H2S and those that have been reported to change in hypoxia. Moreover, our data show that other factors contribute to coordinate the response to H2S, as we found that some H2S-induced transcriptional changes require the skn-1 transcription factor. We demonstrate that, like hif-1, skn-1 is required to survive in low concentrations of H2S. These data suggest a possible model in which HIF-1 and SKN-1 act together to coordinate a transcriptional response to H2S that ultimately leads to alterations in the expression of genes involved in protein homeostasis.


H2S exposure leads to rapid and progressive changes in mRNA abundance

To investigate the transcriptional responses to H2S in C. elegans, we employed a microarray approach to identify mRNAs that were altered in abundance by exposure to H2S. In these experiments, we exposed synchronized cultures of C. elegans to 50 ppm H2S for 1, 12 or 48 h (schematized in Fig. 1A). We measured the response to 50 ppm H2S, as this is the same concentration of H2S that increases lifespan and thermotolerance [7]. Exposure to H2S was immediately prior to harvest. We previously showed that developmental rate is not affected by this concentration of H2S [7], ensuring that all animals were first-day adults when RNA was harvested. This experimental design enabled us to compare transcript abundance without confounding effects from comparing different developmental stages.

Figure 1
Exposure to H2S induces rapid and progressive changes in mRNA abundance.

The number of transcriptional changes we observed in animals exposed to H2S increased with the duration of exposure (Figure 1B). After 1 h exposure we detected significantly altered mRNA levels of 17 transcripts (adjusted p-value<0.05, see Experimental Procedures). All but one of these transcripts was more abundant in the animals exposed to H2S compared to untreated controls (Table 1), suggesting that the transcription of these genes was increased upon exposure to H2S. The effect of H2S on transcript abundance progressed with increased duration of exposure to H2S, both in magnitude of effect on mRNA level and number of gene products affected. After 12 h exposure to H2S, we observed 445 mRNA that were significantly changed, 259 (58%) of which were more abundant after exposure to H2S (Table S1). Nine of the 16 gene products that were increased after 1 h exposure to H2S were still significantly increased after 12 h in H2S. The scope and magnitude of transcript alterations was even more pronounced after 48 h in H2S, the time required to observe an increase in lifespan and thermotolerance [7]. We observed 5,089 genes that had significantly altered mRNA levels after 48 h exposure to H2S (Table S2), which represents 44% of the gene products included in the analysis. Of the significantly altered transcripts, 143 were increased by at least 5.5-fold and 126 were decreased by at least 5.5-fold. Together, these data indicate that there is a rapid and progressive induction of transcriptional activity upon exposure to H2S.

Table 1
Changes in mRNA abundance associated with exposure to H2S.

We used quantitative reverse transcript PCR (qRT-PCR) to validate changes in mRNA levels that we observed by microarray. We focused on transcripts that changed after brief exposure to H2S, reasoning that these early transcriptional changes represent the initial response to H2S and may be important to set up later, more progressive changes. When wild-type animals were grown on E. coli strain OP50, 6 of the 11 gene products that were predicted to be increased based on the microarray data were more abundant after 1 h exposure to H2S (Table 2). Several of the gene products trended toward higher expression, but did not reach significance in this assay. For many of these, we noticed that the level of transcript measured in the untreated sample was near background, which may have increased the variance in the measurements. These changes did reach significance when animals were grown on HT115(DE3) RNAi control food. In these conditions, 10 of the 11 gene products tested were upregulated after 1 h exposure to H2S (Table 2). In general, the magnitude of H2S-induced changes in transcript abundance was greater on RNAi food than on OP50. The source of this variation is unclear, but may hint at an effect of nutritional status on adaptation to H2S. Indeed, the HT115(DE3) bacterial food used for feeding RNAi has previously been shown to affect survival in hypoxia (DLM and MBR, unpublished observation and [15]), consistent with the idea that nutritional status can influence responses to environmental changes. The one gene product that was less abundant after exposure to H2S by microarray, C18H7.1, was not significantly altered after exposure to H2S in either nutrient condition in qRT-PCR measurements. Further validating these results, our microarray data corroborate previous studies that showed T05B4.2 (nhr-57) and K10H10.2 (cysl-2) are upregulated after short exposure to H2S [8]. We conclude that our microarray experiments identified transcripts that change in mRNA abundance associated with exposure to H2S.

Table 2
qRT-PCR Validation of changes in transcript abundance after 1 h H2S.

Function annotation of genes induced by exposure to H2S

Only 16 gene products were more abundant after exposure to H2S for 1 h, which precluded the use of bioinformatic analysis to measure enrichment of functional classes. However, we noticed that many of the genes on this list are annotated to be involved in cellular metabolic processes (Table 1). The most-highly expressed gene in response to H2S is a glutathione-S transferase, F37B1.8 (gst-19). Two other up-regulated genes are predicted to have a role in amino acid metabolism, including the rate-limiting enzyme in serine production, phosphoglycerate dehydrogenase (C31C9.2), and an enzyme with homology to cysteine synthase (cysl-2, K10H10.2). We also observed upregulation of nit-1 (ZK1058.6), a predicted carbon-nitrogen hydrolase. In addition to metabolic enzymes, exposure to H2S also resulted in the upregulation of 6 of the 8 nspe (nematode specific peptide, class E) genes. There is little known about these genes, other than they code for short peptides, 70–75 amino acids long, that are annotated to be integral to the membrane. Although the nspe transcripts were greatly increased in abundance after 1 h exposure to H2S, they were not significantly changed after 12 h exposure to H2S, suggesting that they were only transiently upregulated. Instead, after 12 h in H2S we found that 8 of the 12 nspa (nematode specific peptide, class A) transcripts were significantly upregulated (Table S1).

To evaluate functional categories of genes over-represented in the microarray dataset of transcripts changed after longer exposure to H2S, we employed the online Database for Annotation, Visualization and Integrated Discovery (DAVID, v6.7) functional annotation clustering tool [16], [17]. We focused on gene products that were most increased in response to H2S, as we hypothesize that these might be important for the phenotypic changes that we observe and likely to be most robust. This analysis showed that there were two highly enriched functional clusters in the genes increased by at least 3-fold after 12 h exposure to H2S (Table 3, Table S3). The most enriched cluster included gene ontology terms related to aging and stress resistance. This result is consistent with our previous observation that adaptation to H2S increases lifespan and thermotolerance [7].

Table 3
Functional annotation of gene products increased after exposure to H2S.

We observed that proteins with the F-box motif were enriched in the genes upregulated by H2S. F-box proteins were first identified as members of the SCF (Skp-cullin-F-box) ubiquitin ligase complex that targets proteins for ubiquitination and eventual degradation [18], [19]. F-box proteins were significantly enriched in genes upregulated after 12 h exposure to H2S (Table 3). These data raise the possibility that the response to H2S alters the stability of proteins regulated by the ubiquitin-proteosome system. When we extended this analysis to define the functional clusters in gene products upregulated by at least 5-fold after 48 h exposure to H2S, we observed an even greater enrichment of F-box containing proteins (Table 3). Moreover, the second most enriched cluster included genes associated with the BTB/POZ domain, another protein-protein interaction motif that has been associated with SCF function [20]. We did not observe an enrichment of genes involved in aging or stress resistance after 48 h exposure to H2S. Together, these data show that the transcriptional response to H2S includes many gene products involved in protein turnover mediated by the ubiquitin proteosome system.

The hif-1 transcription factor is required for H2S-induced transcriptional changes

HIF-1 is the C. elegans homologue of the highly conserved hypoxia inducible transcription factor, best known for its role in coordinating the transcriptional response to decreased O2 [9], [10], [21]. Recently, it has been demonstrated that C. elegans require hif-1 to survive exposure to H2S [8]. On exposure to H2S, the HIF-1 protein accumulates and is localized to the nucleus. At least some transcriptional responses to H2S result from the activation of HIF-1, including K10H10.2 (cysl-2), nhr-57, and sqrd-1 [8], [22]. These observations motivated us to consider the possibility that our microarray data might reveal other hif-1-dependent transcriptional responses to H2S.

We measured the abundance of mRNA in hif-1(ia04) mutant animals exposed to H2S for 1 h by qRT-PCR to evaluate if hif-1 is required for transcriptional changes that occur upon exposure to H2S. This short exposure to H2S was sufficient to induce transcriptional changes in wild-type, but was not lethal to the hif-1(ia04) mutant animals. In fact, after 1 h in H2S the hif-1(ia04) mutant animals were still moving normally [8]. We found that for 5 of 11 mRNAs tested, transcript levels in hif-1(ia04) mutant animals were significantly lower than wild-type (Figure 2). This includes the most highly-induced mRNAs. In the hif-1(ia04) mutant animals we observed very little change in any transcript abundance after exposure to H2S, both for messages that were highly induced as well as the lower-expressed transcripts where statistical significance of p<0.05 was not achieved. We conclude that hif-1 has a centrally important function in inducing transcriptional changes associated with exposure to H2S, particularly those changes that occur immediately upon exposure to H2S.

Figure 2
HIF-1 is required for early transcriptional responses to H2S.

Although the response to hypoxia and H2S both require the hif-1 transcription factor, there is little overlap between genes regulated in these two conditions. A previous report identified 68 hypoxia-induced, HIF-1-dependent transcriptional changes in C. elegans [15]. Of the 16 genes upregulated after 1 h exposure to H2S, only 3 (19%) are also regulated by hif-1 in hypoxia: nhr-57, rhy-1 and K10H10.2 (cysl-2) (Figure 2B, Table S4). The slight overlap between these data sets is greater than would be expected by chance (hypergeometric probability 0.001), consistent with our observation that that hif-1 is required for both responses. Similarly, of the 298 transcripts that are more abundant after 12 h exposure to H2S, 8 (3%) are also regulated by hypoxia (hypergeometric probability 0.006; Figure 2C, Table S4). Since hif-1(ia04) mutant animals die after 12 h exposure to H2S we could not determine which of these changes were hif-1-dependent. Thus, we included all hypoxia-induced genes, regardless of whether they require hif-1. These results show that there is statistically significant, though rather minimal similarity between transcriptional responses to H2S and hypoxia. Although our microarrays were not performed under exactly the same conditions as the previous hypoxia studies, these results suggest the interesting possibility that HIF-1 activates different spectrum of targets depending on whether it is activated by hypoxia or H2S. Consistent with this view, HIF-1 activity, as measured by an nhr-57::GFP transcriptional reporter, is in different tissues of animals exposed to hypoxia as compared to H2S [8].

Role of SKN-1 in response to H2S

We noted many Sdz genes were upregulated after 12 or 48 h exposure to H2S, including several with F-box domains (Tables S1 and S2). Sdz genes were named for their skn-1 dependent zygotic expression during embryogenesis between the 4- and 12- cell stage [23]. Several Sdz transcripts that were more abundant after exposure to H2S contained F-box domains (Table S3). skn-1 is a maternally-supplied factor required early in embryogenesis for specification of the EMS blastomere that gives rise to mesendodermal lineages [24], and acts postembryonically in the intestine to control the phase II response to oxidative stress [25] and in the two ASI neurons to control the effects of dietary restriction on lifespan [26]. The abundance of Sdz genes on the list led us to consider the possibility that skn-1 is also involved in the response to H2S.

To test the possibility that H2S-induced transcriptional changes depended on skn-1, we measured mRNA abundance by qRT-PCR from N2 wild-type animals raised on skn-1(RNAi) and exposed to H2S for 1 h as adults. These animals laid only dead eggs, demonstrating that SKN-1 levels were reduced below those required for early embryonic development. We found that 7 of the gene products showed skn-1-dependent changes in abundance after 1 h exposure to H2S (Figure 3). Three genes regulated in a skn-1-dependent manner in response to H2S, C31C9.2, K10H10.2 and ZK1058.6 (nit-1), were previously shown to be regulated by skn-1 postembryonically [27], as was embryonic expression of ZK1058.6 in the EMS lineage after the 4-cell stage and in E descendants after the 300-cell stage [23]. In addition to previously-identified skn-1 dependent transcripts, H2S-induced upregulation of Y38D10A.12 and T05B4.1 also required skn-1. The promoter region of all these genes contain core skn-1 consensus binding sequences, RTACT [27] (Figure 3A). We further observed that three gene products up-regulated by 1 h exposure to H2S were changed even more dramatically in the skn-1(RNAi) animals: K01H10.2, F37B1.8 (gst-19) and W07A12.7 (rhy-1). These data are consistent with reports that skn-1 can act to negatively regulate the expression of genes involved in the response to some oxidative stresses [27]. Our data show that skn-1 acts to both up and down regulate genes in response to H2S.

Figure 3
SKN-1 is essential for appropriate response to H2S.

There is minimal overlap between transcriptional responses to H2S and previously characterized, post-embryonic skn-1 dependent transcripts (Figure 3B, Table S5). We observed significant overlap between the set of genes changed in response to H2S and transcripts that require skn-1 for normal expression in unstressed conditions (16 genes in common, hypergeometric probability 0.006) [27], though 4 of these 16 genes were less abundant after exposure to H2S. Similarly, 10 H2S-regulated genes were identified as skn-1 depenendent targets in response to arsenic stress (hypergeometric probability 0.01) [27]. In contrast, we did not observe significant overlap with skn-1-dependent targets induced by exposure to t-butyl hydroperoxide (hypergeometric probability 0.06) [27], those downregulated in unstressed conditions (hypergeometric probability 0.13) [27], or messages regulated in response to hyperoxia (hypergeometric probability 0.15) [28]. Our observations suggest that skn-1-dependent transcriptional responses to H2S are somewhat distinct from skn-1 mediated responses to oxidative and xenobiotic stress. This is consistent with accumulating evidence that skn-1 transcriptional outputs are context dependent [27].

We monitored the viability of skn-1 mutant animals exposed to H2S to evaluate the functional significance of the skn-1-dependent transcriptional response (Figure 3C). Like wild-type animals, all skn-1(zu169)/nT1 control animals survived exposure to 50 ppm H2S (n = 24, 2 independent experiments). In contrast, none of the skn-1(zu169) homozygous mutant animals tested survived (n = 28 in same 2 experiments). To rule out the possibility that H2S-induced death was a result of the nT1 balancer chromosome, we crossed the balancer away from the zu169 allele. 24% (16/68) of the self-progeny from zu169/+ heterozygotes died when exposed to 50 ppm H2S. This is not statistically different the expected frequency of skn-1(zu169) homozygotes (25%; χ2 = 0.078, df = 1, p>0.05). Of the survivors, 14/15 were fertile, indicating that these animals were not homozygous for the skn-1 allele, which confers a maternal-effect lethal phenotype. We conclude that the skn-1(zu169) homozygous animals died when exposed to H2S. These data demonstrate that SKN-1 activity is essential to appropriately respond to H2S.


Our results indicate that hif-1 and skn-1 cooperate to orchestrate a progressive transcriptional response to H2S. Previous studies have demonstrated hif-1 dependent responses to H2S [8], [22]. We have extended this observation using an unbiased microarray approach that identified several new hif-1-dependent transcriptional changes upon exposure to H2S. In addition, we have identified skn-1 as another essential factor during exposure to H2S.

H2S protects mice from otherwise lethal whole-body hypoxia [2] and improves outcome in a variety of rodent models of ischemia/reperfusion injury [4], [29], [30]. The mammalian orthologue of HIF-1 has been implicated in protection against ischemia/reperfusion in mammals [31]. Thus, the observation that HIF-1 is activated by H2S suggests a mechanistic basis for the beneficial effects of H2S [8]. Curiously, our data suggest that there is little overlap between transcriptional targets of HIF-1 in hypoxia and H2S. These results may indicate that H2S does not mediate protection against ischemia simply by inducing a standard hypoxia response. Further understanding this conserved adaptive response to H2S will provide new insight into mechanisms that can improve homeostasis in changing conditions.

Our finding that SKN-1 plays a role in the response to H2S is consistent with a recent report that nuclear accumulation of NRF2, a mammalian homologue of SKN-1, is correlated with H2S-induced protection against from ischemia-induced heart failure [30]. SKN-1 controls the Phase II response to toxins and oxidative stress [25]. However, we do not favour the hypothesis that H2S activates the Phase II response. The canonical Phase II targets gst-4 nor gcs-1 were not dramatically induced by H2S, and the H2S-induced upregulation of another glutathione S-transferase, gst-19, was exaggerated in skn-1(RNAi) animals. Moreover, we observed little overlap in genes regulated by skn-1 in response to xenobiotic or oxidative stress and those that are changed in H2S. Instead, we found at least 7 of the gene products included in the F-box and BTB/POZ clusters (Table 3) were identified as skn-1-dependent zygotic transcripts [23]. Finally, we did not observe obvious accumulation of skn-1::GFP in the intestinal nuclei of animals exposed to H2S using epifluoresence microscopy (not shown), although we cannot rule out the possibility that the nuclear enrichment of GFP was below the detection limit in this experiment. We suggest that, during adaptation to H2S, skn-1 may play a role in remodeling the protein turnover machinery.

Protein homeostasis is increasingly appreciated for its importance to aging and age-associated decline [32]. We propose that one consequence of adaptation to H2S is to increase transcription of genes related to protein turnover by the ubiquitin ligase and proteasome system, including F-box and BTB/POZ domain proteins. In this model, the effect of H2S to increase lifespan and thermotolerance may be attributed, at least in part, from effects on protein homeostasis. Further understanding the mechanisms by which adaptation to H2S can improve homeostasis and influence lifespan may provide novel insights into the mechanisms that mediate the beneficial effects of H2S in mammals.

Materials and Methods

Nematode strains and culture

Strains used were N2 wild-type (Bristol), ZG31 hif-1(ia04), and EU35 skn-1(zu169)/nT1[unc-?(n754) let-?] (IV;V). EU35 and RNAi strains mentioned below were a gift from Dr. Jim Priess (Division of Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA). The ia04 mutation deletes the second through fourth exons of HIF-1 and is a predicted molecular null [9]. zu169 is an ochre mutation in an exon of skn-1 shared among all isoforms. The zu169 mutation is maternal effect lethal [33], abrogates paraquat-induced expression of gst-4 and gcs-1 in the intestine [25] and prevents increased lifespan in response to dietary restriction [26].

C. elegans were grown on nematode growth medium plates seeded with live Escherichia coli OP50 food (NGM/OP50 plates) as described previously [34]. All experiments and worm culturing were conducted at room temperature to avoid effects resulting from changing temperature. Exposure to H2S was in continuous flow H2S chambers that were created as previously described [7], by continuously diluting 5000 ppm H2S/balance N2 (Airgas, Seattle, WA USA) with house air to a final concentration of 50 ppm H2S. For viability assays, worms were exposed to H2S as fourth-stage larvae (L4) and scored for survival after 18–24 h. The effect of skn-1(zu169) mutations on viability in H2S was determined by picking sterile, non-Unc progeny of skn-1/nT1 animals, and compared to Unc heterozygous siblings. To cross the zu169 allele away from the nT1 balancer, N2 males were crossed with skn-1/nT1 mutant hermaphrodites. Non-Unc heterozygous F1 were allowed to produce F2 progeny, which were scored for sensitivity to H2S and the skn-1 maternal effect lethal phenotype.

skn-1(RNAi) animals were generated by feeding N2 from starved L1 on HT115(DE3) bacteria carrying either the skn-1 clone or empty vector control (L4440) from the Ahringer library [35]. RNAi by feeding was essentially as described [35]. Bacteria expressing the dsRNA was diluted from an overight culture grown in LB containing 25 mM carbenicillin and 10 mM tetracycline, regrown to OD600~0.6 in LB with 25 mM carbenicillin and then seeded onto NGM-lite plates that contained 3 mM isopropyl β-D-1-thiogalactopyranoside and 25 mM carbenicillin. RNAi plates were allowed to dry overnight, stored at 4 C, and used within 5 days of being seeded. skn-1(RNAi) adults laid only dead embryos.

RNA sample isolation

For microarray and qRT-PCR experiments, animals were synchronized as starved first-stage larvae (L1) after isolating embryos by hypochlorite treatment. For microarray analysis, 3000 starved L1 were distributed onto 15 cm NGM/OP50 plates, with each independent replicate performed on a different day. For quantitative RT-PCR (qRT-PCR), 1000 L1 larvae were distributed onto 10 cm NGM/OP50 plates. C. elegans were exposed to H2S on plates for the amount of time indicated immediately prior to harvest. All animals were harvested as first-day gravid adults (schematized in Figure 1A). For nematode harvest, plates were removed from H2S, the worms were immediately rinsed off the plates with distilled water, caught on a 43-micron nylon filter and collected by centrifugation. 100 µL of sedimented worms were added to 900 µL Trizol (Invitrogen, Carlsbad, CA), frozen in liquid nitrogen and stored at −70 C. Less than 2 min elapsed from when plates were removed from H2S until samples were frozen. Frozen samples were thawed, vortexed for 30 s, and the RNA was isolated following the protocol included with the Trizol manual, followed by isopropanol precipitation.

Microarray expression profiling and analysis

Each RNA sample was labelled, hybridized to a single-channel Nimblegen 4×72 K (build 160) expression array, and scanned following manufacturer's suggested protocols by the Fred Hutchinson Cancer Research Center's DNA Array Facility. Three biological replicates for each H2S-treated sample (1, 12 and 48 h exposure) and 5 biological replicates from untreated controls were used. Data were RMA normalized and probe-level data were summarized with the NimbleScan software. Genes with weak signal intensity across all groups and those with low variability across samples were excluded from further analysis. Each H2S-treated sample was statistically compared to a matched untreated control using the Bioconductor package limma [36]. The false discovery rate (FDR) method of Benjamini and Hochberg [37] was used to adjust p-values for multiple testing. An adjusted p-value≤0.05 was used to define differential expression. Results were annotated using WormBase WS190 ( Expression results and microarray raw intensity files, in compliance with MIAME guidelines, can be accessed through the Gene Expression Omnibus ( and are accessible through GEO series accession number GSE25199. Functional annotation clustering analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7 ( Gene list submitted for 12 h exposure to H2S included only the 91 gene products with logFC>1.6 (fold change>3). For the 48 h dataset, the 95 gene products with logFC>2.5 (fold change>5.5) were included. In both cases, the C. elegans background list in the database was used with default analysis parameters. Using only gene products included in the analysis after filtering did not alter the results. Annotation clusters that included at least one term with p<0.05 were considered to be functionally enriched clusters.

Hypergeometric probabilities were calculated including all 11,522 features included in the microarray analysis as the population, with successes and sample size as defined in the text ( The number of successes in each sample (overlap) was determined by manually comparing data from H2S-induced changes measured by our microarray experiments and hypoxia-induced genes [15] or skn-1-dependent transcripts [27], [28]. Probabilities less than 0.05 were considered significant. Core skn-1 consensus sites in the promoter region of candidate transcripts were defined manually, based on the published consensus RTCAT [27]. The promoter region was defined as 2 kb upstream of the start site. The ModENCODE database [38] was searched to determine if any of the 7 transcripts changed in a skn-1-dependent changes in response to H2S were shown to have SKN-1::GFP bound in the promoter.

Real-time RT-PCR (qRT-PCR)

Quantitative RT-PCR (qRT-PCR) was used to validate microarray measurements and determine if H2S-dependent changes occurred in hif-1(ia04) or skn-1(RNAi) animals. mRNA was isolated as described above, and cDNA was synthesized from 300 ng of RNA using the included random primers using the ProtoScript M-MuLV First Strand cDNA Synthesis Kit (New England Biolabs) according to manufacturer's suggested protocol. Primers to amplify cDNA targets were designed using Primer3 ( When possible, primer pairs spanned a small intron so that genomic and cDNA amplification products could be distinguished by agarose gel electrophoresis. Primer sequences are available upon request. Primers were tested to ensure amplification of the correct size genomic target, and then calibrated against serial dilutions of genomic DNA. qRT-PCR reactions were performed using an ep realplex2 S (Eppendorf). Each 10 µL reaction contained 5 µL 2X KAPA SYBR green Master Mix (Kapa Biosystems), 0.45 µL cDNA and 3 µL primers (10 µM each primer). Reactions were performed in duplicate and at least three independent biological replicates were included for each condition. Each experiment included primers that amplified only genomic DNA (negative controls to identify background signal levels) as well as 4 control targets (sir-2.1, tba-1, irs-2, and hil-1) that are not affected by H2S exposure for normalization. ΔCt for each gene product was calculated by subtracting Ct values from the geometric mean of the control targets [39]. ΔCt were averaged across the three experiments. Student's t-test was used to evaluate differences between ΔCt values of treated samples and untreated controls (EXCEL). For differences between genotypes (Figures 2 and and3),3), p-values were calculated by one-way Anova from summary statistics (mean, standard deviation, n) ( Reported fold-changes were calculated as 2∧−ΔΔCt [40], where ΔΔCt = ΔCt(H2S)−ΔCt(untreated). Error bars on graphs represent standard deviation, which was carried through the fold-change calculation using standard error propagation (reported as “variance”).

Supporting Information

Table S1

Transcripts that are significantly changed after 12 h exposure to H2S, listed in order of magnitude fold-change.


Table S2

Transcripts that are significantly changed after 48 h exposure to H2S, listed in order of magnitude fold-change.


Table S3

Genes included in functional annotation clusters.


Table S4

Transcripts altered by both 1 h exposure to H2S and hypoxia.


Table S5

Transcripts regulated by skn-1 in other studies that are altered by exposure to H2S.



Microarray experiments and analysis were performed with the expert assistance of Dr. Jeff Delrow (FHCRC DNA Array Facility), and Dr. Jerry Davison (FHCRC Computational Biology) provided valuable assistance with data analysis. Some nematode strains used in this work were provided by the Caenorhabditis Genetics Center, which is funded by the NIH National Center for Research Resources (NCRR).


Competing Interests: The authors have declared that no competing interests exist.

Funding: This work was supported by National Institutes of Health grant K99/R00 AG033050 to Dana L Miller. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Blackstone E, Morrison M, Roth MB. H2S induces a suspended animation-like state in mice. Science. 2005;308:518. [PubMed]
2. Blackstone E, Roth MB. Suspended animation-like state protects mice from lethal hypoxia. Shock. 2007;27:370–372. [PubMed]
3. Morrison ML, Blackwood JE, Lockett SL, Iwata A, Winn RK, et al. Surviving blood loss using hydrogen sulfide. J Trauma. 2008;65:183–188. [PubMed]
4. Elrod JW, Calvert JW, Morrison J, Doeller JE, Kraus DW, et al. Hydrogen sulfide attenuates myocardial ischemia-reperfusion injury by preservation of mitochondrial function. PNAS. 2007;104:15560–15565. [PubMed]
5. Simon F, Giudici R, Duy CN, Schelzig H, Oter S, et al. Hemodynamic and metabolic effects of hydrogen sulfide during porcine ischemia/reperfusion injury. Shock. 2008;30:359–364. [PubMed]
6. Jha S, Calvert JW, Duranski MR, Ramachandran A, Lefer DJ. Hydrogen sulfide attenuates hepatic ischemia-reperfusion injury: role of antioxidant and antiapoptotic signaling. Am J Physiol Heart Circ Physiol. 2008;295:H801–806. [PubMed]
7. Miller DL, Roth MB. Hydrogen sulfide increases thermotolerance and lifespan in Caenorhabditis elegans. PNAS. 2007;104:20618–20622. [PubMed]
8. Budde MW, Roth MB. Hydrogen sulfide increases hypoxia-inducible factor-1 activity independently of von Hippel-Lindau tumor suppressor-1 in C. elegans. Mol Biol Cell. 2010;21:212–217. [PMC free article] [PubMed]
9. Jiang H, Guo R, Powell-Coffman JA. The Caenorhabditis elegans hif-1 gene encodes a bHLH-PAS protein that is required for adaptation to hypoxia. PNAS. 2001;98:7916–7921. [PubMed]
10. Semenza GL. HIF-1: mediator of physiological and pathophysiological responses to hypoxia. J Appl Physiol. 2000;88:1474–1480. [PubMed]
11. Haigis MC, Yankner BA. The aging stress response. Mol Cell. 2010;40:333–344. [PMC free article] [PubMed]
12. Chen D, Thomas EL, Kapahi P. HIF-1 modulates dietary restriction-mediated lifespan extension via IRE-1 in Caenorhabditis elegans. PLoS Genet. 2009;5:e1000486. [PMC free article] [PubMed]
13. Mehta R, Steinkraus KA, Sutphin GL, Ramos FJ, Shamieh LS, et al. Proteasomal regulation of the hypoxic response modulates aging in C. elegans. Science. 2009;324:1196–1198. [PMC free article] [PubMed]
14. Zhang Y, Shao Z, Zhai Z, Shen C, Powell-Coffman JA. The HIF-1 hypoxia-inducible factor modulates lifespan in C. elegans. PLoS ONE. 2009;4:e6348. [PMC free article] [PubMed]
15. Shen C, Nettleton D, Jiang M, Kim SK, Powell-Coffman JA. Roles of the HIF-1 hypoxia-inducible factor during hypoxia response in Caenorhabditis elegans. J Biol Chem. 2005;280:20580–20588. [PubMed]
16. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. [PubMed]
17. Dennis G, Jr, Sherman BT, Hosack DA, Yang J, Gao W, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3. [PubMed]
18. Skowyra D, Craig KL, Tyers M, Elledge SJ, Harper JW. F-box proteins are receptors that recruit phosphorylated substrates to the SCF ubiquitin-ligase complex. Cell. 1997;91:209–219. [PubMed]
19. Bai C, Sen P, Hofmann K, Ma L, Goebl M, et al. SKP1 connects cell cycle regulators to the ubiquitin proteolysis machinery through a novel motif, the F-box. Cell. 1996;86:263–274. [PubMed]
20. Geyer R, Wee S, Anderson S, Yates J, Wolf DA. BTB/POZ domain proteins are putative substrate adaptors for cullin 3 ubiquitin ligases. Mol Cell. 2003;12:783–790. [PubMed]
21. Semenza GL. HIF-1 and human disease: one highly involved factor. Genes Dev. 2000;14:1983–1991. [PubMed]
22. Budde MW, Roth MB. The Response of Caenorhabditis elegans to Hydrogen Sulfide and Hydrogen Cyanide. Genetics. 2011 epub ahead of print PMID: 21840852. [PMC free article] [PubMed]
23. Robertson SM, Shetty P, Lin R. Identification of lineage-specific zygotic transcripts in early Caenorhabditis elegans embryos. Dev Biol. 2004;276:493–507. [PubMed]
24. Bowerman B, Eaton BA, Priess JR. skn-1, a maternally expressed gene required to specify the fate of ventral blastomeres in the early C. elegans embryo. Cell. 1992;68:1061–1075. [PubMed]
25. An JH, Blackwell TK. SKN-1 links C. elegans mesendodermal specification to a conserved oxidative stress response. Genes Dev. 2003;17:1882–1893. [PubMed]
26. Bishop NA, Guarente L. Two neurons mediate diet-restriction-induced longevity in C. elegans. Nature. 2007;447:545–549. [PubMed]
27. Oliveira RP, Porter Abate J, Dilks K, Landis J, Ashraf J, et al. Condition-adapted stress and longevity gene regulation by Caenorhabditis elegans SKN-1/Nrf. Aging Cell. 2009;8:524–541. [PMC free article] [PubMed]
28. Park SK, Tedesco PM, Johnson TE. Oxidative stress and longevity in Caenorhabditis elegans as mediated by SKN-1. Aging Cell. 2009;8:258–269. [PMC free article] [PubMed]
29. Nicholson CK, Calvert JW. Hydrogen sulfide and ischemia-reperfusion injury. Pharmacol Res. 2010;62:289–297. [PMC free article] [PubMed]
30. Calvert JW, Coetzee WA, Lefer DJ. Novel insights into hydrogen sulfide–mediated cytoprotection. Antioxid Redox Signal. 2010;12:1203–1217. [PMC free article] [PubMed]
31. Loor G, Schumacker PT. Role of hypoxia-inducible factor in cell survival during myocardial ischemia-reperfusion. Cell Death Differ. 2008;15:686–690. [PubMed]
32. Morimoto RI, Cuervo AM. Protein homeostasis and aging: taking care of proteins from the cradle to the grave. J Gerontol A Biol Sci Med Sci. 2009;64:167–170. [PMC free article] [PubMed]
33. Bowerman B, Tax FE, Thomas JH, Priess JR. Cell interactions involved in development of the bilaterally symmetrical intestinal valve cells during embryogenesis in Caenorhabditis elegans. Development. 1992;116:1113–1122. [PubMed]
34. Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77:71–94. [PubMed]
35. Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, et al. Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003;421:231–237. [PubMed]
36. Smyth GK. Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W, editors. Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. 2005. pp. 397–420.
37. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological) 195;57:289–300.
38. Celniker SE, Dillon LA, Gerstein MB, Gunsalus KC, Henikoff S, et al. Unlocking the secrets of the genome. Nature. 2009;459:927–930. [PMC free article] [PubMed]
39. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034. [PMC free article] [PubMed]
40. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3:1101–1108. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science