A central aspect of aging research concerns the question of when individuality in lifespan arises1. Here we show that a transient increase in reactive oxygen species (ROS), which occurs naturally during early development in a subpopulation of synchronized Caenorhabditis elegans, sets processes in motion that increase stress resistance, improve redox homeostasis and ultimately prolong lifespan in those animals. We find that these effects are linked to the global ROS-mediated decrease in developmental histone H3K4me3 levels. Studies in HeLa cells confirmed that global H3K4me3 levels are ROS-sensitive and that depletion of H3K4me3 levels increases stress resistance in mammalian cell cultures. In vitro studies identified SET1/MLL histone methyltransferases as redox sensitive units of the H3K4-trimethylating complex of proteins (COMPASS). Our findings implicate a link between early-life events, ROS-sensitive epigenetic marks, stress resistance and lifespan.
This is a preview of subscription content
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Get time limited or full article access on ReadCube.
All prices are NET prices.
Finch, C. E. & Tanzi, R. E. Genetics of aging. Science 278, 407–411 (1997).
Rea, S. L., Wu, D., Cypser, J. R., Vaupel, J. W. & Johnson, T. E. A stress-sensitive reporter predicts longevity in isogenic populations of Caenorhabditis elegans. Nat. Genet. 37, 894–898 (2005).
Holmström, K. M. & Finkel, T. Cellular mechanisms and physiological consequences of redox-dependent signalling. Nat. Rev. Mol. Cell. Biol. 15, 411–421 (2014).
Schulz, T. J. et al. Glucose restriction extends Caenorhabditis elegans life span by inducing mitochondrial respiration and increasing oxidative stress. Cell Metab. 6, 280–293 (2007).
Dillin, A. et al. Rates of behavior and aging specified by mitochondrial function during development. Science 298, 2398–2401 (2002).
Ristow, M. & Schmeisser, S. Extending life span by increasing oxidative stress. Free Radic. Biol. Med. 51, 327–336 (2011).
Knoefler, D. et al. Quantitative in vivo redox sensors uncover oxidative stress as an early event in life. Mol. Cell 47, 767–776 (2012).
Gutscher, M. et al. Real-time imaging of the intracellular glutathione redox potential. Nat. Methods 5, 553–559 (2008).
Labbadia, J. & Morimoto, R. I. Repression of the heat shock response is a programmed event at the onset of reproduction. Mol. Cell 59, 639–650 (2015).
Greer, E. L. et al. Members of the H3K4 trimethylation complex regulate lifespan in a germline-dependent manner in C. elegans. Nature 466, 383–387 (2010).
Shilatifard, A. The COMPASS family of histone H3K4 methylases: mechanisms of regulation in development and disease pathogenesis. Annu. Rev. Biochem. 81, 65–95 (2012).
Guenther, M. G., Levine, S. S., Boyer, L. A., Jaenisch, R. & Young, R. A. A chromatin landmark and transcription initiation at most promoters in human cells. Cell 130, 77–88 (2007).
Pu, M. et al. Unique patterns of trimethylation of histone H3 lysine 4 are prone to changes during aging in Caenorhabditis elegans somatic cells. PLoS Genet. 14, e1007466 (2018).
Wang, W. et al. SET-9 and SET-26 are H3K4me3 readers and play critical roles in germline development and longevity. eLife 7, e34970 (2018).
Dou, Y. et al. Regulation of MLL1 H3K4 methyltransferase activity by its core components. Nat. Struct. Mol. Biol. 13, 713–719 (2006).
Southall, S. M., Wong, P. S., Odho, Z., Roe, S. M. & Wilson, J. R. Structural basis for the requirement of additional factors for MLL1 SET domain activity and recognition of epigenetic marks. Mol. Cell 33, 181–191 (2009).
Shilatifard, A. Molecular implementation and physiological roles for histone H3 lysine 4 (H3K4) methylation. Curr. Opin. Cell Biol. 20, 341–348 (2008).
Li, T. & Kelly, W. G. A role for Set1/MLL-related components in epigenetic regulation of the Caenorhabditis elegans germ line. PLoS Genet. 7, e1001349 (2011).
Leichert, L. I. et al. Quantifying changes in the thiol redox proteome upon oxidative stress in vivo. Proc. Natl Acad. Sci. USA 105, 8197–8202 (2008).
Kenyon, C. J. The genetics of ageing. Nature 464, 504–512 (2010).
Yang, W. & Hekimi, S. A mitochondrial superoxide signal triggers increased longevity in Caenorhabditis elegans. PLoS Biol. 8, e1000556 (2010).
Weiner, A. et al. Systematic dissection of roles for chromatin regulators in a yeast stress response. PLoS Biol. 10, e1001369 (2012).
Hansen, M., Flatt, T. & Aguilaniu, H. Reproduction, fat metabolism, and life span: what is the connection? Cell Metab. 17, 10–19 (2013).
Han, S. et al. Mono-unsaturated fatty acids link H3K4me3 modifiers to C. elegans lifespan. Nature 544, 185–190 (2017).
Sun, L., Sadighi Akha, A. A., Miller, R. A. & Harper, J. M. Life-span extension in mice by preweaning food restriction and by methionine restriction in middle age. J. Gerontol. A 64A, 711–722 (2009).
Brenner, S. The genetics of Caenorhabditis elegans. Genetics 77, 71–94 (1974).
Stroustrup, N. et al. The Caenorhabditis elegans Lifespan Machine. Nat. Methods 10, 665–670 (2013).
Koopman, M. et al. A screening-based platform for the assessment of cellular respiration in Caenorhabditis elegans. Nat. Protocols 11, 1798–1816 (2016).
Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).
Trapnell, C., Pachter, L. & Salzberg, S. L. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25, 1105–1111 (2009).
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Ng, S. B. et al. Exome sequencing identifies MLL2 mutations as a cause of Kabuki syndrome. Nat. Genet. 42, 790–793 (2010).
Dorgan, K. M. et al. An enzyme-coupled continuous spectrophotometric assay for S-adenosylmethionine-dependent methyltransferases. Anal. Biochem. 350, 249–255 (2006).
Southall, S. M., Cronin, N. B. & Wilson, J. R. A novel route to product specificity in the Suv4-20 family of histone H4K20 methyltransferases. Nucleic Acids Res. 42, 661–671 (2014).
We thank M. Malinouski and T. Mullins for assistance with the reconfiguration of the Biosorter; G. Csankovszki for antibodies, C. elegans RNAi feeding clones and comments; B. Braeckman for the N2jrIs2[Prpl-17::Grx-1-roGFP2] strain; J. Nandakumar for HeLa (EM-2-11ht) cells; the Caenorhabditis Genetics Center (funded by National Institutes of Health Infrastructure Program P40 OD010440) for C. elegans strains; the DNA Sequencing Core (BRCF), R. Tagett, W. Wu and the Bioinformatics Core of University of Michigan for RNA sequencing and data analysis; K. Wan for protein purification; R. Sawarkar and J. Labbadia for important suggestions; Jakob laboratory members for comments on the manuscript and J. Bardwell for critically reading the manuscript. Mass spectrometry was performed by MS Bioworks. This work was supported by NIH grants GM122506 and AG046799 as well as the Priority Program SPP 1710 of the Deutsche Forschungsgemeinschaft (Schw823/3-2) to U.J., a NIH T32 Career Training in the Biology of Aging grant to D.B., a NIH T32 Career Training in the Biology of Aging grant and a Bright Focus ADR Fellowship (A2019250F) to B.J.O., and the National Natural Science Foundation of China (31470737) to Y.C.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Peer review information Nature thanks Anne Brunet, Michael Ristow and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Extended data figures and tables
Extended Data Fig. 1 In vivo readout of endogenous redox states at different stages during C. elegans lifespan and sorting parameters of oxidized and reduced subpopulations.
a, Microscopy analysis of the Grx1–roGFP2 ratio of individual N2jrIs2[Prpl-17::Grx1-roGFP2] worms (symbol) cultivated at 15 °C and imaged at the indicated time points. Data points that are not significantly different from each other (P > 0.05) are labelled with the same letter. Data are mean ± s.e.m; n, number of worms; one-way ANOVA with Tukey correction. b, The Grx–roGFP2 ratio (A405/A488) was calculated using the partial profiling feature (pp) configured to analyse extinction and emission data from 488 nm and 405 nm lasers that sequentially excited each worm. c, A population of N2jrIs2[Prpl-17::Grx1-roGFP2] at the L2 stage, separated based on opacity (extinction) and length (time of flight) was gated as R1.
a–d, Microscopy analysis of the Grx1–roGFP2 ratio of individual worms (dots) previously sorted into L2ox, L2mean and L2red subpopulations. n, number of worms; one-way ANOVA with Tukey correction. e–i, Survival curves of sorted L2ox, L2mean and L2red worms. For n numbers, P values (log-rank test) see Extended Data Table 2. Insets, Grx1–roGFP2 ratio of individual worms (dots), assessed by fluorescence microscopy after sorting. n, number of worms. For P values (two-sided unpaired t-test) see Extended Data Table 2.
a, Length measurements of L2ox and L2red worms (symbol) from nose to tail tip immediately after sorting. No significant difference; P = 0.4735 (unpaired two-sided t-test). b, Brood size of L2ox and L2red worms, measured at the indicated time points. n, number of worms. No significant difference within a single age; P = 0.6532 (two-way ANOVA). c–f, Basal respiration (c), maximal (d) and spare (e) respiratory capacity and basal rates of flux through glycolysis (f) of L2ox and L2red worms. n = 3 independent sorting experiments. ECAR, extracellular acidification rate; OCR, oxygen consumption rate. P = 0.9469 (c), P = 0.7784 (d), P = 0.7904 (e) and P = 0.7925 (f); two-sided unpaired t-test. g, Survival of L2ox, L2mean and L2red worms 20 h after heat shock. n = 5 independent sorting experiments; two-sided unpaired t-test. The connected data points represent data from the same sorting experiment. The survival of L2ox is set to 1. All data are mean ± s.e.m.
a, Steady-state transcript levels of selected oxidative stress-related genes in L2ox and L2red worms. n, number of independent sorting experiments; unpaired two-sided t-test. Data are mean ± s.e.m. b, Volcano plot showing fold changes versus P values for the transcriptomes of L2ox and L2red subpopulations. DEGs (P ≤ 0.05) are represented by red dots (see Methods for statistical definition of DEGs). Data were collected from four independent sorting experiments. c, GSEA of the 327 DEGs. Normalized enrichment scores (see Methods for calculation) are represented by the bar graph. Terms (for summary, see Supplementary Table 1) indicating origin, process or phenotype associated with genes known to have a role in the process are shown on the left. Some terms (*) have been merged and are represented as a single category bar for simplicity (for detailed values, see Supplementary Table 1). d, e, Percentage of DEGs identified in L2ox that intersect with H3K4me3 peak signals within their 5′ region (500 bp upstream and downstream from the transcription start site). The H3K4me3 chromatin immunoprecipitation (ChIP) datasets were generated from L3-staged N2 worms: ChIP chip, GEO entry GSE30789 in d and chromatin immunoprecipitation followed by sequencing (ChIP–seq), GEO entry GSE28770 in e, indicating that these marks are set during larval development. Hypergeometric probability: d, P = 0.064; e, P = 2.786 × 10−6. In f, g, Venn diagrams show the overlap among upregulated (f) or downregulated (g) gene sets in L2ox and downregulated or upregulated set-9(rw5) and set-26(tm2467) gene sets (GEO entry: GSE100623). See Supplementary Table 1 for datasets in d–g.
Extended Data Fig. 5 Redox sensitivity of in vivo H3K4m3e3 levels and in vitro histone methyltransferase complex activity.
a, Global H3K4me3 levels in the sorted L2ox and L2red worms. A representative western blot using antibodies against H3K4me3 is shown. b, c, Quantification of global H3K27ac (b) and H3K27me3 (c) levels by western blot. n = 3 independent sorting experiments. P = 0.3793 (b) and P = 0.0905 (c); unpaired two-sided t-test. Data represent mean ± s.e.m. d–f, Global H3K4me3 (d), ASH2L (e) and MLL1 (f) levels in HeLa cells before and after H2O2 treatment, as assessed by western blot. g, Time course of the in vitro methyltransferase reaction for core COMPASS subunits (SET domain of MLL1–WDR5–ASH2L–RBBP5). Reaction rates were derived from the first 20 min of the linear range. h–j, In vitro histone methyltransferase assays of core COMPASS subunits, consisting of purified GST–WDR5 (WDR5), GST–ASH2L (ASH2L), GST–RBBP5 (RBBP5) and either GST–MLL1 SET domain or untagged MLL1 SET domain (h), GST–SET1A SET domain (i) or GST–SET1B SET domain (j). Superscript ox indicates that the protein was pre-treated with either 1 mM (+) or 2 mM (++) H2O2 for 30 min before the activity assay. DTT was added after the H2O2 treatment. n = 3 independent experiments; one-way ANOVA with Sidak correction. Data are mean ± s.e.m. k, The MLL1 SET domain was treated with either 2 mM DTT, 2 mM H2O2 or 2 mM H2O2, followed by 4 mM DTT. Catalase was used to quench the H2O2. The proteins were denatured and thiols were modified with NEM before loading onto non-reducing SDS–PAGE to prevent non-specific thiol oxidation. The proteins were visualized by silver staining. M, marker. l, MLL1 SET domain treated with either 2 mM DTT, 2 mM H2O2 or 2 mM H2O2, followed by 4 mM DTT. All reduced protein thiols were then labelled with the 500-Da thiol-reactive compound AMS, causing a 500-Da mass decrease per oxidized thiol—detectable on reducing SDS–PAGE. M, protein marker. m, Cysteine oxidation state in MLL1 SET domain after treatment with either 2 mM DTT or 2 mM H2O2 followed by NEM labelling as assessed by LC–MS/MS. The peptide containing Cys3967 could not be detected. n, Schematic representation of the redox sensitivity of the MLL1 SET domain. For blot and gel source images, see Supplementary Figs. 1 and 3. o, Sequence alignment of the SET domain. All cysteines in MLL1 are shown in bold, and the five absolutely conserved cysteines are highlighted in yellow. Cysteines shown to be involved in zinc coordination are marked with an asterisk. NCBI protein BLAST and Clustal Omega Multiple Sequence Alignment, Clustal O (1.2.4) were used.
Extended Data Fig. 6 Effects of H3K4me3 downregulation on heat-shock response and endogenous redox state.
a, ASH-2 and SET-2 transcript levels of N2jrIs2[Prpl-17::Grx1-roGFP2] worms treated with ash-2 or set-2 RNAi for 2 generations. n = 6 (ash-2) and n = 2 (set-2) independent experiments; unpaired two-sided t-test. E.V., empty vector. b, ASH-2 protein levels in N2jrIs2[Prpl-17::Grx1-roGFP2] worms treated with control RNAi or ash-2 RNAi for two generations using western blot analysis. c, H3K4me3 levels in N2jrIs2[Prpl-17::Grx1-roGFP2] worms treated with control RNAi, ash-2 RNAi or set-2 RNAi for two generations. d, Transcript levels of selected heat-shock genes after heat-shock treatment of N2jrIs2[Prpl-17::Grx1-roGFP2] worms treated with the indicated RNAi. n = 3 independent experiments; one-way ANOVA with Bonferroni correction. e, Transcript levels of selected heat-shock genes in set-2 or wdr-5.1 mutants before and after heat-shock treatment. n = 3 independent experiments; one-way ANOVA with Bonferroni correction. f, ASH2L levels following ASH2L siRNA treatment of HeLa cells. n = 2 independent experiments. g; Grx1–roGFP2 ratios of L2 larval worms treated with ash-2 RNAi, set-2 RNAi, wdr-5.1 RNAi or the empty vector were measured using the BioSorter. n = 4 (ash-2, set-2) and n = 3 (wdr-5.1) independent sorting experiments; unpaired two-sided t-test. h, Representative survival curves of N2jrIs2[Prpl-17::Grx1-roGFP2] worms treated with ash-2 or set-2 RNAi for two generations and treated with 1 mM paraquat (PQ) for 10 h at the L2 larval stage. For n numbers, repetitions and statistics (log-rank), see Extended Data Table 4. Data in a, d–g, represent mean ± s.e.m. For blot source images, see Supplementary Figs. 1, 3.
Supplementary Figure 1-3: Full gel and blot images.
Supplementary Table 1: RNA-sequencing data including DEGs and GSEA analysis related to Extended Data Fig. 3b and c. Data used for Venn diagrams and pie charts in Fig. 3b, Extended Data Fig. 3d-g are also included.
About this article
Cite this article
Bazopoulou, D., Knoefler, D., Zheng, Y. et al. Developmental ROS individualizes organismal stress resistance and lifespan. Nature 576, 301–305 (2019). https://doi.org/10.1038/s41586-019-1814-y
Nature Reviews Molecular Cell Biology (2022)
Cell Death Discovery (2022)
Nature Chemical Biology (2022)
Cellular and Molecular Life Sciences (2022)
Global profiling of distinct cysteine redox forms reveals wide-ranging redox regulation in C. elegans
Nature Communications (2021)