Aging Biology

Research Paper

Genetic Dissection Reveals New Components of the NHR-49 Multi-Stress Resistance and Longevity Pathway

Authors

Kelsie R. S. Doering,1,2,3,4,5,* Glafira Ermakova,1,2,3,4 Brendil Sabatino,1,2,3 Rachel Cheng,1,2,3,6 Luxcia Kugathasan,1,2,3,7 Zoe DeBoer1,2,3 and Stefan Taubert1,2,3,4,*

1Centre for Molecular Medicine and Therapeutics, The University of British Columbia, Vancouver, BC, Canada

2Edwin S.H. Leong Centre for Healthy Aging, The University of British Columbia, Vancouver, BC, Canada

3British Columbia Children’s Hospital Research Institute, Vancouver, BC, Canada

4Department of Medical Genetics, The University of British Columbia, Vancouver, BC, Canada

5Department of Biological and Health Sciences, Kwantlen Polytechnic University, Surrey, BC, Canada

6Now at AbCellera Biologics Inc., Vancouver, BC, Canada

7Now at Department of Medicine, Division of Nephrology, Toronto General Hospital, University of Toronto, Toronto, ON, Canada

*Corresponding authors: Kelsie.doering@kpu.ca, taubert@cmmt.ubc.ca

DOI:https://doi.org/10.59368/agingbio.20250040

Received: 10/11/2024, Revised: 5/1/2025, Accepted: 5/8/2025, Published: 6/6/2025

Full Article |  PDF |  Supplementary 

Abstract

Cells, tissues, and organisms encounter many environmental and endogenous stresses, and the ability to mount specific responses to these stresses is critical to achieve homeostasis and survival. Stress also occurs in many pathological states and is therefore of great biomedical relevance, both as a cause or contributor to diseases and as a potential therapeutic target. Pathways that regulate the responses to oxidative stress, hypoxia, and starvation are evolutionarily conserved from the nematode Caenorhabditis elegans to humans and involve master regulators of gene expression. One such stress response pathway engages C. elegans NHR-49, an orthologue of mammalian PPARα, which is required to survive tBOOH-induced oxidative stress, hypoxia, and starvation, and is also a key effector in many genetic longevity pathways. However, in contrast to other stress response pathways such as SKN-1/Nrf2 and HIF-1/HIF signaling, many aspects of NHR-49 signaling remain poorly understood. To map new components of this pathway, we performed reverse genetic screens, uncovering 80 candidate kinases, transcription factors, and transcriptional coregulators required for NHR-49-controlled gene activation in tBOOH, hypoxia, and/or starvation. Further analysis revealed that the kinase hpk-1/HIPK is required in the nhr-49 pathway for hypoxia and oxidative stress survival, that pseudokinase nipi-3 is required with nhr-49 for hypoxia survival, and that transcription factor nhr-80 is required with nhr-49 for hypoxia and oxidative stress survival and for the normal lifespan of wild-type C. elegans. Collectively, our data identify numerous new potential regulators in an important stress response and thus substantially enhance our understanding of stress adaptation and longevity regulation.

Introduction

The ability to survive in the face of stress is fundamental for cells, tissues, and organisms in their natural surroundings, and altered stress responses are also a feature of many pathological states in humans, such as cancer. Common stresses include oxidative stress, hypoxia, and starvation. Oxidative stress occurs when reactive oxygen species (ROS), obligate and ubiquitous byproducts of aerobic respiration, accumulate within the cell to toxic levels. ROS are highly electrophilic and cause cellular damage by oxidizing DNA, proteins, and lipids, which promotes the activation of cytoprotective and damage repair pathways13. Hypoxia occurs when cellular oxygen levels are too low for normal physiological functions, resulting in major metabolic and developmental rewiring4,5. Finally, starvation causes major changes in cellular metabolism, including depletion of energy stores via lipolysis, autophagy, and fatty acid β-oxidation. Thus, organisms rewire their physiology and metabolism in the face of limited nutrient availability68.

The pathways that regulate the response to these stresses are evolutionarily conserved and typically require master regulators of gene expression to orchestrate changes in developmental and physiological gene programs. Oxidative stress activates evolutionarily conserved transcription factors of the nuclear factor erythroid 2-related factor 2 (Nrf2) family. In the animal model Caenorhabditis elegans, the Nrf2 family member SKiNhead-1 (SKN-1) is modified by several kinases and translocates to the nucleus, where it activates the expression of cellular detoxification and antioxidant genes, reestablishing homeostasis and promoting survival914. Similarly, transcription factors of the hypoxia-inducible factor (HIF) family play a critical role in hypoxia response in animals. For example, C. elegans typically activates the Egg-laying defective 9 (egl-9)–von Hippel–Lindau (vhl-1)–hif-1 pathway to promote metabolic changes in hypoxia, thus enabling survival5,15. Finally, the transcription factor EB (TFEB) family is critical for the response to starvation6,11,1619. In starving C. elegans, the TFEB ortholog HLH-30 translocates to the nucleus to regulate the expression of many lysosomal and autophagy genes6,8,17,20. In addition to their role in acute stress responses, these genes also play essential roles in maintaining the normal lifespan of wild-type C. elegans and the extended lifespan of many long-lived mutants21, highlighting the importance of stress response pathways in longevity.

Although the C. elegans SKN-1/Nrf2, HIF-1, and HLH-30/TFEB pathways are critical for the oxidative stress, hypoxia, and starvation response, respectively, other pathways of equal importance to stress survival have been discovered. Several studies show that the nuclear hormone receptor (NHR) transcription factor NHR-49, originally identified as a regulator of fatty acid metabolism, is necessary for the resistance to hypoxia, oxidative stress, infection with pathogens, and for the normal lifespan of wild-type and long-lived C. elegans strains2230. Notably, nhr-49 acts in parallel to hif-1 in hypoxia and regulates skn-1-independent gene responses in oxidative stress31,32, highlighting that several parallel pathways have evolved to provide C. elegans with complementary and optimal protection in the face of stress.

The activation and implementation of SKN-1-, HIF-1-, and HLH-30-driven signaling in C. elegans is complex and involves numerous factors, including kinases, transcriptional coregulators, and other transcription factors. For example, SKN-1 activity is controlled by post-translational modification by kinases such as P38 Map Kinase 1 (PMK-1), Glycogen Synthase Kinase 3 (GSK-3), AKT-1/-2, Serum- and Glucocorticoid-inducible Kinase 1 (SGK-1), and No Induction of Peptide after Drechmeria coniospora Infection 3 (NIPI-3)10,14,33,34, by the proteasome via WD Repeat protein 23 (WDR-23) and SKp1 Related 1 and 2 (SKR-1/-2)35,36, and by interactions with transcriptional coregulators such as cAMP response element-binding protein (CREB)-binding protein 1 (CBP-1) and Mediator subunit MDT-1537,38. Collectively these factors fine-tune SKN-1 activity in various settings related to stress, pathogen infection, and longevity.

Less is known about how the nhr-49 pathway is regulated, especially in conditions of stress. NHR-49 can dimerize with other NHRs3941, and some such partnerships (e.g., with NHR-66, -79, and -80) may regulate specific sets of lipid metabolism genes4244. However, how NHR-49 promotes certain stress response programs is unclear. Furthermore, although we and others showed that mdt-15 is a coregulator required for the induction of nhr-49-dependent stress response genes32, NHR-49 is likely to engage additional coregulators to achieve specialized stress responses. AMP-activated kinase 2 (aak-2) may regulate NHR-4945,46, and NHR-49 may also be controlled by fatty acids that could act as NHR-49 ligands47, but overall, little is known about how NHR-49 is activated by stress. Thus, a substantial gap in knowledge remains on whether and how transcription factors, coregulators, and kinases function in NHR-49-controlled stress responses.

To discover new genes acting with NHR-49 to regulate stress responses, we conducted a reverse genetic screen using a flavin-containing monooxygenase 2 transcriptional reporter (fmo-2p::gfp) as a readout. fmo-2p::gfp reporter induction depends on nhr-49 in oxidative stress caused by tert-butyl hydroperoxide (tBOOH), hypoxia, and starvation, but is independent of skn-1 in tBOOH and only weakly depends on hif-1 in hypoxia and hlh-30 in starvation31,32. Our screen identified genes unique to one stress, genes involved in two stresses, and a set of 29 genes required for the upregulation of fmo-2p::gfp in all three stresses. The latter set included the known NHR-49 coregulator mdt-15 and homeodomain interacting protein kinase 1 (hpk-1), which acts with nhr-49 in the hypoxia response31, validating our screen. Follow-up studies on hpk-1, the pseudokinase nipi-3, and the transcription factor nhr-80 revealed new potential partners of nhr-49-dependent stress responses with specific requirements in some, but not other, stresses. Our work has thus uncovered factors that may delineate stress-specific or general signaling branches of the NHR-49 stress response pathway.

Materials and Methods

Nematode strains and growth conditions

We cultured C. elegans strains on nematode growth media (NGM) plates using standard techniques. To avoid background effects, each mutant was crossed into our lab wild-type N2 strain at least six times. E. coli strain OP50 was the food source in all experiments except for RNAi experiments, for which we used E. coli strain HT115. All experiments were carried out at 20 °C. C. elegans strains used in this study are listed in Table 1. For synchronized C. elegans growths, embryos were isolated by standard sodium hypochlorite treatment48. Isolated embryos were allowed to hatch overnight on unseeded NGM plates until the population reached a synchronized halted development at L1 stage via short-term fasting (12–24 hours). Synchronized L1 stage larvae were then transferred to OP50-seeded plates and grown to the desired stage.

C. elegans strains used in this study.
StrainGenotypeReference
N2Wild type83
VE40eavEx20[fmo-2p::gfp + rol-6(su1006)] 32
STE117nhr-49(et13) I; eavEx20[fmo-2p::gfp + rol-6(su1006)] 32
STE68nhr-49(nr2041) I 28
EK273hpk-1(pk1393) X 84
STE132nhr-49(nr2041) I; hpk-1(pk1393) X 31
AGP25fglmEx5 (nhr-49p::nhr-49::gfp + myo-2p::mCherry) 27
STE142hpk-1(pk1393) X; glmEx5 (nhr-49p::nhr-49::gfp + myo-2p::mCherry) 31
AVS310artEx27[Phpk-1::hpk-1::GFP + pRF4(rol-6 (su1006))] 54
STE197nhr-49(nr2041) I; artEx27[Phpk-1::hpk-1::GFP + pRF4(rol-6 (su1006))] This study
IG544nipi-3(fr4) X 58
STE198nhr-49(nr2041) I; nipi-3(fr4) X This study
STE70nhr-80(tm1011) III 43
STE199nhr-49(nr2041) I;nhr-80(tm1011) This study

Feeding RNA interference

RNAi was performed on NGM plates supplemented with 25 μg/mL carbenicillin (BioBasic CDJ469), 1 mM IPTG (Santa Cruz CAS 367-93-1), and 12.5 μg/mL tetracycline (BioBasic TB0504; NGM-RNAi plates) and seeded with appropriate HT115 RNAi bacteria. RNAi sub-libraries covered 1,543 genes total, including 332 kinases, 768 transcription factors, and 443 coregulators40,49. RNAi clones were from the Ahringer and Vidal libraries50,51; for a complete list, see Supplementary Table 1H.

fmo-2p::gfp RNAi screens

Egg lays were performed for 3–4 hours on 12-well NGM-RNAi plates with VE40 eavEx20[fmo-2p::gfp + rol-6(su1006)] animals (hereafter fmo-2p::gfp animals) until there were approximately 60 eggs per well. Adults were then removed, and synchronized progeny (which include transgenic fmo-2p::gfp animals displaying the Roller phenotype and non-transgenic, non-Roller siblings) were grown for three days at 20 °C to the L4/young adult stage. At that time, animals were transferred to 10 mM tBOOH RNAi plates seeded with heat-inactivated HT115 bacteria for 5 hours for the tBOOH screen; transferred to 0.5% O2 for 4 hours followed by 1 hour in 21% O2 for the hypoxia screen; or washed at least five times with M9 to remove all excess bacteria and transferred to unseeded plates for 8 hours for the starvation screen. For the nhr-49(et13) screen, fmo-2p::gfp;nhr-49(et13) animals synchronized as described above were grown to the L4/young adult stage. Two control wells were always included during the screen: L4440 empty vector (EV) RNAi as a negative control and nhr-49 RNAi as a positive control (on fmo-2p::gfp or fmo-2p::gfp;nhr-49(et13) animals, as appropriate).

As a screen readout, we performed semi-quantitative scoring of fmo-2p::gfp fluorescence. Specifically, approximately 20 L4/young adult stage C. elegans were assessed for whole-body GFP fluorescence levels under a Leica M205FA fluorescence microscope, and fluorescence levels were scored as high, medium, low, or off for each animal in each well. In the tBOOH screen, wells with at least 40% of animals showing medium, low, or off GFP fluorescence were counted as “hits” (Supplementary Table 1A–C). For all subsequent screens, wells with at least 30% of animals showing low or off GFP fluorescence were counted as “hits” (Supplementary Table 1A, 1D–G). Following one repeat of all appropriate RNAi clones in each screen, all initial “hits” were re-screened. Only RNAi clones which yielded “hits” in all three repeats (original screen and two re-screens) were considered validated hits (see Supplementary Fig. 1A). Any culture that impaired animal development was subsequently removed from the total number of cultures tested for that screen, regardless of hit status. This included removal of 37 cultures from the tBOOH screen and removal of 2 cultures from the hypoxia screen (Supplementary Fig. 1A). All 80 RNAi cultures representing final hits were sequenced to verify that the correct target gene was targeted.

For the tBOOH screen, we used the high, medium, low, or off designations to calculate an “activity score” (AS) for each RNAi clone, as follows: AS = (0*[number of “off” animals]/[number of total transgenic animals]) + 1*[number of “low” animals]/[number of total transgenic animals]) + 2*[number of “medium” animals]/[number of total transgenic animals]) + 0*[number of “high” animals]/[number of total transgenic animals]). Thus, the minimum AS possible is 0, and the maximum possible AS is 3. We then used unpaired, two-tailed Student’s T-tests to calculate the significance of AS for each hit versus EV control (Supplementary Table 1C).

Analysis of fluorescent reporter lines via differential interference contrast (DIC) and fluorescence microscopy

To analyze fmo-2p::gfp fluorescence in the RNAi screens, C. elegans were grown and treated in the appropriate condition as described above. For NHR-49::GFP experiments, adult animals were transferred to 10 mM tBOOH plates for 3 hours. For microscopy, C. elegans were collected into M9 buffer containing 0.06% levamisole (Sigma L9756) for immobilization on 2% (w/v) agarose pads. We captured images at 10x magnification on a CoolSnap HQ camera (Photometrics) attached to a Zeiss Axioplan 2 compound microscope, followed by MetaMorph Imaging Software with Autoquant 3D digital deconvolution. All images for the same experiment were captured at the same exposure time. Images were analyzed using ImageJ software (https://imagej.nih.gov/ij/download.html), with fluorescence calculated by taking the difference of the background fluorescence (area beside the animal) from the mean whole animal fluorescence. For each experiment, at least three independent trials were performed with a minimum of 30 C. elegans per condition.

tBOOH survival assay

Synchronized C. elegans were grown to the late L4 stage and then transferred to plates containing 6 mM tBOOH seeded with heat-inactivated OP50. Dead and surviving animals were quantified twice a day for 96 hours or until all animals were dead, whichever came first. Death was assessed by prodding the animals with a platinum wire. Animals that crawled off the plate or died due to rupturing or internal hatching were censored. GraphPad Prism was used to generate survival curves and calculate statistical significance using the log-rank (Mantel–Cox) test, and statistics are listed in the accompanying tables.

Hypoxia sensitivity assays

Hypoxic conditions were maintained using continuous flow chambers, as previously described52. Compressed gas tanks (5000 ppm O2 balanced with N2) were certified as standard to be within 2% of indicated concentration from Praxair Canada (Delta, BC). Oxygen flow was regulated using Aalborg rotameters (Aalborg Instruments and Controls, Inc., Orangeburg, NY, USA). Hypoxic chambers (and room air controls) were maintained in a 20 °C incubator for the duration of the experiments.

For embryo survival assays, gravid first-day adult animals (picked as L4 the previous day) were allowed to lay eggs for 1–4 hours on plates seeded with 15 μL OP50 the previous day. Adults were removed, and eggs were exposed to 0.5% O2 for 24 hours. After being placed back into normoxia (21% O2) for 65 hours, animals were scored for developmental success, as defined as reaching at least L4 stage. For all normoxia comparison experiments, methods were as described above except plates were kept in room air (21% O2) for the duration (instead of being exposed to 0.5% O2).

L1 starvation recovery assay

Gravid day 1 adult C. elegans were picked onto agar plates seeded with E. coli (OP50), with 15–20 animals per plate and 2–3 plates per strain, and maintained at 20 °C. Four days later, a first sodium hypochlorite treatment was conducted to obtain synchronized eggs (see above). The bleached eggs were incubated on a rotator in 15 mL tube in S basal solution with antibiotics/antimycotics at 20 °C overnight. On the following day, hatched L1 stage animals were plated onto seeded plates and then kept in a 20 °C incubator. When the C. elegans reached day 1 adulthood (96 hours post first bleach), a second sodium hypochlorite was conducted, after which the bleached eggs were incubated on a rotator in 15 mL tube in S basal solution with antibiotics/antimycotics at 20 °C overnight. On the following day, hatched L1 animals were diluted to achieve a concentration of 1 C. elegans/μL in the 15 mL tubes. On the same day, 150 μL of C. elegans suspended in S Basal were plated onto plates seeded with OP50 for day 0. Recovery counts were performed 48 hours after plating to assess the percentage of animals that reached the L4 stage. Subsequent plating and survival counting were carried out on days 2 and 4, and continued every other day until day 26 or until all animals failed to recover. Statistical analysis and graphical representation of % recovery from starvation were done with GraphPad Prism using one-way ANOVA of the AUC of the individual recovery curves. To eliminate effects due to intrinsic developmental delays of some mutants, we normalized baseline viability of each strain to 100%.

Lifespan assay

Gravid day 1 adult C. elegans were picked onto agar plates seeded with E. coli OP50, with 15–20 animals per plate and 2–3 plates per strain, and maintained at 20 °C. Four days later, hypochlorite treatment was conducted to obtain synchronized eggs (see above). The bleached eggs were incubated on a rotator in 15 mL tube in S basal solution with antibiotics/antimycotics at 20 °C overnight. On the following day, hatched L1 stage animals were plated onto seeded plates and then kept in a 20 °C incubator. Ninety-six hours after sodium hypochlorite treatment, adult longevity experiments commenced by picking synchronized adult animals onto freshly seeded plates. Animals were scored as alive/dead/censored and moved to newly seeded plates every day for the first 10 days. Animals were scored every other day until day 30 or until all animals were dead. Statistical analysis and graphical representation of the C. elegans lifespan were done with GraphPad Prism using Kaplan-Meier survival analysis.

Results

Reverse genetic screen reveals 80 genes required to activate fmo-2 in oxidative stress, hypoxia, and/or starvation

To identify new genes in the NHR-49-regulated tBOOH oxidative stress response pathway, we prepared RNAi sub-libraries covering 1,543 genes, including 332 kinases, 768 transcription factors, and 443 coregulators (see the Materials and Methods section, Supplementary Table 1H). We then tested each RNAi clone for its requirement to induce the NHR-49-activated fmo-2p::gfp promoter reporter after exposing C. elegans to 10 mM tBOOH for 3 hours, that is, testing for nhr-49 phenocopy (Fig. 1A). Thirty-seven cultures caused animal development to become impaired and were subsequently removed from the total number of cultures tested for the tBOOH screen, bringing the total number of genes tested to 1,506. Screening all clones leading to healthy animals one time revealed 194 potential hits, including 35 kinases, 77 transcription factors, and 82 coregulators. These 194 initial hits were then screened again in the same conditions, resulting in 122 hits. These 122 hits were then screened one more time in the same conditions, which resulted in a final total of 80 hits, including 9 kinases, 16 transcription factors, and 55 coregulators (Table 2, Supplementary Table 1A–C, Supplementary Fig. 1A).


Reverse genetic RNAi screen reveals new genes that may function in the nhr-49 stress response pathway.

(A) Representative micrographs of fmo-2p::gfp adult animals fed empty vector (EV) or nhr-49 RNAi (negative and positive controls, respectively) or RNAi of hpk-1 (screen hit) or nhr-80 in non-stressed conditions, after 5 hours exposure to 10 mM tBOOH, after 4 hours exposure to 0.5% O2 followed by 1 hour recovery in 21% O2, and after 8 hours without food. (B) The Venn diagram shows the overlap of hits from RNAi screens for genes required to induce fmo-2p::gfp in oxidative stress, hypoxia, starvation, or in the nhr-49(et13) background (see Supplementary Table 1G for gene lists). Three genes that were hits in all screens are highlighted: mdt-15, a known nhr-49 coregulator32,41, and hpk-1 and nipi-3, which were further investigated in this study.

Eighty genes are required to induce the fmo-2p::gfp reporter in tBOOH oxidative stress.
WormBase Gene IDGene NameHuman HomologTypeHypoxia HitStarvation Hitet13 Hit
Kinases
C36B1.10gskl-2 GSK3BMAP kinaseXXX
F35C8.2sek-4 MAP2K6MAP kinaseXXX
W04B5.5 PDPK2serine/threonine kinase X 
C01C4.3 SBK1serine/threonine kinaseX  
F20B6.8hpk-1 HIPKserine/threonine kinaseXXX
M79.1abl-1 ABL1tyrosine kinase X 
K09A9.1nipi-3 TRIB1pseudokinaseXXX
T10D4.6  GTPase bindingXXX
B0024.6gcy-6 NPR1Guanylyl CyclaseX  
Transcription factors
C25A1.11aha-1 ARNTaryl hydrocarbon receptor associated proteinXX 
T05A7.4hmg-11 HMGA2high mobility group AT-hookX X
C28A5.4ceh-43 DLX1homeoboxX X
ZC247.3lin-11 LHX1homeobox X 
F31E8.3tab-1 BSXhomeoboxXX 
W01D2.1 RPL37ribosomal protein  X
T27A1.6mab-9 TBX20transcription factorX  
C01G12.1madf-5  transcription factorXXX
T22E7.2dhhc-11 ZDHHC11zinc-finger X 
F13D11.2hbl-1  zinc-finger  X
C27A12.3let-391 ZFP2zinc-fingerX  
K01H12.3nhr-38  zinc-fingerXX 
K12H6.1nhr-119  zinc-fingerX  
Y40B1A.4sptf-3 SP3zinc-fingerXX 
ZC123.3zfh-2 ZFHX3zinc-fingerX X
F52F12.6ztf-11 MYT1Lzinc-fingerX X
Coregulators
Y48B6A.11jmjd-2 KDM4demethylation   
F42A6.5  DNA repairXXX
C06E1.4glr-1 GRIA2/3glutamate receptor   
F33H2.1dog-1 BRIP1helicase XX
C26D10.2hel-1 DDX39Ahelicase  X
W08D2.7mtr-4 MTREXhelicaseXXX
Y48G8AL.6smg-2 UPF1helicase X 
F55F8.2 DDX24helicase X 
F35H10.1his-30 H2AXhistoneX X
F17E9.10his-32 H3Chistone XX
R06F6.1cdl-1 SLBPhistone binding proteinXX 
F54A5.1hmbx-1 HMBOX1homeobox   
C01H6.9hasp-1 HASPINkinase X 
F47A4.2dpy-22 MED12Lmediator complexXXX
Y54E5B.3let-49 MED7mediator complex  X
T09A5.6mdt-10 MED10mediator complex XX
R12B2.5mdt-15 MED15mediator complexXXX
C55B7.9mdt-18 MED18mediator complexXX 
Y17G7B.2ash-2 ASH2Lmethylation   
R06A4.7mes-2 EZH1/2methylationX  
F46A9.6mec-8 RBPMSmRNA processingX  
F37E3.1ncbp-1 NCBP1nuclear cap binding proteinXXX
F26A3.2ncbp-2 NCBP2nuclear cap binding protein XX
F20D12.2mcm-3AP MCM3APnuclear exportXXX
R12C12.2ran-5 RANBP3nuclear export X 
K07F5.13npp-1 NUP54nucleoporinX X
T01G9.4npp-2 NUP85nucleoporinXX 
Y54E5A.4npp-4  nucleoporinXXX
F56A3.3npp-6 NUP160nucleoporin XX
T19B4.2npp-7 POM121nucleoporinXXX
C03D6.4npp-14 NUP214nucleoporinXX 
C29E4.4npp-15 NUP133nucleoporinXX 
Y56A3A.17npp-16 NUP50nucleoporin XX
Y43F4B.4npp-18 SEH1Lnucleoporin  X
R06F6.5npp-19 NUP35nucleoporinXX 
Y77E11A.13npp-20 SEC13nucleoporin X 
F36A4.7ama-1 POLR2ARNA polymerase  X
C26E6.4rpb-2 POLR2BRNA polymeraseXXX
C41G7.4set-32  SET domain containing protein X 
C48B6.2 IMP3snoRNA X 
Y54E2A.3tac-1  spindle-associatedXX 
W02D3.11hrpf-1  splicing X 
C50C3.6prp-8 PRPF8splicing XX
W02B12.2rsp-2 SRSF5/6splicing   
C08B11.5sap-49 SF3B4splicing X 
C50D2.5sftb-6 SF3B6splicingX  
Y49E10.15snr-6 SNRPEsplicing  X
Y71F9B.4snr-7 SNRPGsplicing X 
Y116A8C.35uaf-2 U2AF1splicing XX
C23H3.3ascc-1 ASCC1transcription coregulatorXX 
K03H1.10cbp-2 CBP/p300transcription coregulator XX
C25H3.6mdt-26 TCEANC2transcription elongationXXX
F46F11.4ubl-5 UBL5ubiquitin-likeX X
ZK637.4   XXX
F01G12.5let-2 COL4Acollagen   
The table lists each hit, their predicted mammalian homolog, the type of protein it represents, and results of follow-up screens in hypoxia, starvation, and in the nhr-49(et13) background (“X” indicates hit). Hits are organized in groups of similar predicted molecular function/type. All 80 RNAi cultures were sequenced to verify that the correct target gene was knocked down.

In addition to its role in the oxidative stress response, nhr-49 is required to induce genes in hypoxia and for survival in and recovery from transient hypoxia stress and from starvation31,32. We therefore tested if the genes required for fmo-2p::gfp induction in oxidative stress are also required for its induction in hypoxia and starvation (Fig. 1A). Screening 78 candidates from the oxidative stress screen in hypoxia conditions revealed 49 potential hits, which were then screened twice more with the same assay, resulting in 44 validated hits, including 7 kinases, 12 transcription factors, and 25 coregulators (Table 2, Supplementary Fig. 1A, Supplementary Table 1A,D; note: two of 80 RNAi clone cultures identified in the tBOOH screen impaired animal development and were removed from the hypoxia screen analysis, resulting in 78 tested genes). Screening 80 candidates from the oxidative stress screen in starvation revealed 52 potential hits, which were then screened twice more with the same assay, resulting in 52 validated hits, including 7 kinases, 7 transcription factors, and 38 coregulators (Table 2, Supplementary Fig. 1A, Supplementary Table 1A,E). In total, of the genes studied in these screens, 13 genes only regulate fmo-2p::gfp activity in tBOOH; 15 genes regulate fmo-2p::gfp activity in tBOOH and hypoxia; 23 genes regulate fmo-2p::gfp activity in tBOOH and starvation; and 29 genes regulate fmo-2p::gfp activity in all three stresses, including 5 kinases, 5 transcription factors, and 19 coregulators (Supplementary Fig. 1B, Table 2, Supplementary Table 1G). Several of these 29 genes have known roles in stress responses and/or longevity, including the kinases hpk-1, nipi-3, and gskl-2, and the transcription factor aha-131,34,5358. As expected, the known NHR-49 binding transcriptional coregulator mdt-1532,41 was also in this list, validating our approach.

Forty genes are required for fmo-2 induction in the nhr-49(et-13) background

The stress conditions studied above likely engage numerous regulators that cooperatively activate fmo-2. To identify genes required for NHR-49 activity in the absence of stress, we crossed the fmo-2p::gfp reporter into a strain bearing the nhr-49(et13) gain-of-function allele. The et13 mutation affects the ligand binding domain of NHR-49, is sufficient to induce fmo-2 in the absence of stress, and promotes hypoxia resistance31,59. These properties make it a useful tool to study factors acting directly with NHR-49. Screening of the 80 genes required to induce fmo-2p::gfp in oxidative stress in the nhr-49(et13) background revealed 53 potential hits, which were then screened twice more with the same assay, resulting in 40 hits, including 5 kinases, 7 transcription factors, and 28 coregulators (Supplementary Fig. 1A, Table 2, Supplementary Table 1A,F). Comparing the list of genes required for fmo-2p::gfp reporter induction in the nhr-49(et13) background with the lists of genes required for induction in different stresses revealed that 17 genes are essential in all four conditions, including 5 kinases, 1 transcription factor, and 11 coregulators; these genes most closely phenocopy nhr-49’s requirement (Fig. 1B).

HPK-1 functions within the NHR-49-controlled hypoxia and oxidative stress response pathways

Our screens revealed that the kinase hpk-1 is required for fmo-2p::gfp induction in oxidative stress, hypoxia, starvation, and in the nhr-49(et13) background. We previously showed that hpk-1 functions upstream of NHR-49 in the hypoxia response pathway, with hpk-1(pk1393) mutant animals being sensitive to hypoxia and epistasis experiments using the nhr-49(nr2041);hpk-1(pk1393) double mutant showing that hpk-1 acts in the same pathway as nhr-49. Furthermore, an overexpressed NHR-49::GFP fusion protein is regulated post-translationally in hypoxia, and hpk-1 is required for increased NHR-49::GFP protein levels in low oxygen31. To determine if hpk-1 acts in the same pathway as nhr-49 to promote oxidative stress survival, wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) L4 animals were exposed to 6 mM tBOOH. All strains were sensitive to tBOOH compared to wild type (Fig. 2A). Epistasis analysis revealed that survival of the nhr-49(nr2041);hpk-1(pk1393) double mutant is not different from that of the single mutants, demonstrating that these two genes act in the same pathway to promote tBOOH resistance (Fig. 2A).


hpk-1 is required for resistance to oxidative stress and starvation and for normal lifespan, acting partly in the nhr-49 pathway.

(A) Survival plot of wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) L4 animals exposed to 6 mM tBOOH. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001 versus wild-type animals; ns versus nhr-49(nr2041);hpk-1(pk1393). Data and statistical analysis for additional repeats in Supplementary Table 2. (B) Fractional population recovery to L4 after L1 starvation (in percent, y-axis) over time (days of starvation, x-axis) of wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) animals. Error bars represent SD from at least three independent experiments. p-values were derived by calculating area under the curve (AUC) for each strain using one-way ANOVA with Dunnett’s multiple comparisons test. *p < 0.05, **p < 0.01, ***p < 0.001 (n = 3, >100 animals per time point, normalized to 100% viability for each strain). (C) Lifespan of wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) animals. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 3. ns = not significant.

To determine if hpk-1 acts in the same pathway as nhr-49 to promote recovery from starvation, wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) mutant L1 stage C. elegans were starved and tested for their ability to reach at least the L4 stage after recovering on food for 48 hours. Similar to nhr-49, hpk-1 animals were sensitive to starvation compared to wild type (Fig. 2B). Interestingly, epistasis analysis revealed that nhr-49(nr2041);hpk-1(pk1393) double mutant animals are significantly more starvation-sensitive than each single mutant, suggesting that these two genes act in separate pathways to promote recovery from starvation (Fig. 2B).

nhr-49 and hpk-1 are both required for normal lifespan in C. elegans27,28,53,54, but whether they interact genetically in this context is unknown. To test if hpk-1 acts in the same pathway as nhr-49 to promote lifespan, we analyzed population survival of wild-type, nhr-49(nr2041), hpk-1(pk1393), and nhr-49(nr2041);hpk-1(pk1393) animals (Fig. 2C). Consistent with previous findings, the nhr-49 and hpk-1 single mutants are short-lived compared to wild type. Interestingly, the nhr-49;hpk-1 double mutant displayed an additive defect, suggesting that these genes act in parallel pathways to promote normal C. elegans lifespan (Fig. 2C).

NHR-49 is regulated post-translationally during tBOOH-induced oxidative stress32 and in hypoxia31, and hpk-1 is required for the elevated levels of an overexpressed NHR-49::GFP protein in hypoxia31. To determine if hpk-1 is required for induction of this NHR-49::GFP protein during tBOOH-induced oxidative stress, we exposed C. elegans bearing the nhr-49p::nhr-49::gfp translational reporter (encoding overexpressed NHR-49::GFP)27 to 10 mM tBOOH for 3 hours in the wild-type and hpk-1(pk1393) mutant background. hpk-1 mutation reduced NHR-49::GFP levels in the wild-type background; in addition, tBOOH exposure significantly induced NHR-49::GFP levels in the wild-type background, but this induction was abrogated in the hpk-1(pk1393) background (Fig. 3A,B). Thus, hpk-1 is required for NHR-49::GFP protein induction during tBOOH oxidative stress, as seen in hypoxia31.


HPK-1 induces NHR-49::GFP in oxidative stress.

(A) Representative micrographs show nhr-49p::nhr-49::gfp and nhr-49p::nhr-49::gfp;hpk-1(pk1393) adult animals following 3 hours exposure to 10 mM tBOOH. (B) The graph shows the quantification of whole animal GFP levels in nhr-49p::nhr-49::gfp and nhr-49p::nhr-49::gfp;hpk-1(pk1393) animals following 3 hours exposure 10 mM tBOOH (three repeats totaling > 30 individual animals per strain). **p < 0.01, ***p < 0.001 (two-way ANOVA corrected for multiple comparisons using the Tukey method). (C) Survival plot of nhr-49p::nhr-49::gfp and nhr-49p::nhr-49::gfp;hpk-1(pk1393) and non-GFP sibling L4 animals exposed to 6 mM tBOOH. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 4. (D) Survival plot of hpk-1p::hpk-1::gfp and hpk-1p::hpk1::gfp;nhr-49(nr2041) and non-GFP sibling L4 animals exposed to 6 mM tBOOH. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 5. (E) Working model of how NHR-49 and HPK-1 interact in response to oxidative stress, hypoxia, starvation, and for promotion of normal lifespan. ns = not significant.

To determine if this requirement for hpk-1 in NHR-49 protein induction affected animal survival on tBOOH, we exposed NHR-49::GFP overexpressing animals in the wild-type and hpk-1 mutant backgrounds to 6 mM tBOOH. As negative controls, non-GFP siblings from each strain were used. Similar to previous findings, overexpressing NHR-49::GFP slightly enhanced the survival of C. elegans on tBOOH compared to the non-GFP sibling controls32 (Fig. 3C). Importantly, loss of hpk-1 abolished this extension (Fig. 3C), suggesting that hpk-1 is required for the increased stress resistance due to NHR-49::GFP overexpression.

In addition, we tested if stress resilience caused by HPK-1 overexpression was dependent on nhr-49. We studied C. elegans bearing an hpk-1p::hpk-1::gfp transgene (hereafter HPK-1::GFP)54 in wild-type and nhr-49(nr2041) mutant backgrounds. HPK-1::GFP overexpression increased the population survival of animals on tBOOH compared to non-GFP siblings (Fig. 3D). Critically, loss of nhr-49 abolished this enhanced stress resistance (Fig. 3D). Taken together, our results suggest that, like in hypoxia, HPK-1 acts together with NHR-49 in tBOOH-induced oxidative stress, possibly upstream of NHR-49, to induce NHR-49 protein levels, thus activating NHR-49 to promote animal survival (Fig. 3E).

NIPI-3 acts in the NHR-49 pathway to promote hypoxia survival

Our RNAi screens revealed the pseudokinase nipi-3 as a new gene required for fmo-2p::gfp induction in tBOOH-induced oxidative stress, hypoxia, starvation, and in the nhr-49(et13) background (Table 2, Supplementary Table 1A–F). To determine if nipi-3 is required to promote oxidative stress survival, wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) L4 animals were exposed to 6 mM tBOOH. Unexpectedly, loss of nipi-3 enhanced survival on tBOOH compared to wild-type (Fig. 4A,E). Epistasis analysis revealed that nipi-3 loss restored the oxidative stress sensitivity of the nhr-49 mutant to wild-type levels (Fig. 4A,E). This suggests that these two factors act in separate tBOOH survival pathways, with nhr-49 promoting and nipi-3 inhibiting survival.


C. elegans nipi-3 loss causes oxidative stress resistance and hypoxia sensitivity.

(A) The graph shows population survival of wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) L4 stage animals exposed to 6 mM tBOOH. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 6. (B) The graph shows the average population survival of wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) embryos following 24 hours exposure to 0.5% O2, then allowed to recover at 21% O2 for 65 hours, and counted as the ability to reach at least L4 stage (n = 4, totaling > 100 individual animals per strain). p-values were derived using ordinary one-way ANOVA corrected for multiple comparisons using the Tukey method. *p < 0.05, ****p < 0.0001. (C) Fractional population recovery to L4 after L1 starvation (in percent, y-axis) over time (days of starvation, x-axis) of wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) animals. Error bars represent SD from at least three independent experiments. p-values were derived by calculating AUC for each strain using one-way ANOVA with Dunnett’s multiple comparisons test. *p < 0.05, **p < 0.01 (n = 4, >100 animals per time point, normalized to 100% viability for each strain). (D) Lifespans of wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) animals. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 7. (E) Working model of how NHR-49 and NIPI-3 interact in response to oxidative stress and hypoxia. ns = not significant.

To determine if nipi-3 is required to promote hypoxia survival, wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) embryos were exposed to 0.5% O2 for 24 hours and then allowed to recover to the L4 or later stage after being placed back in normoxia (21% O2) for 65 hours. We found that nipi-3 is required for embryo survival and recovery in hypoxia (Fig. 4B,E). Epistasis analysis revealed that the nhr-49(nr2041);nipi-3(fr4) double mutant is not significantly more sensitive to hypoxia than the nhr-49 single mutants, suggesting that these two genes act in the same hypoxia response pathway (Fig. 4B,E).

To determine if nipi-3 acts in the same pathway as nhr-49 to promote recovery from starvation, wild-type, nhr-49(nr2041), nipi-3(fr4), and nhr-49(nr2041);nipi-3(fr4) L1 animals were starved and tested for their ability to reach the L4 stage after recovering on food for 48 hours. Unlike nhr-49, nipi-3 animals were not sensitive to starvation (Fig. 4C). Similarly, nipi-3 loss does not affect lifespan in the wild-type background on its own (Fig. 4D). Curiously, however, nhr-49(nr2041);nipi-3(fr4) double mutant animals live significantly shorter than nhr-49(nr2041) single mutants (Fig. 4D). This suggests a potential epistatic effect of nipi-3 in the absence of nhr-49.

NHR-80 acts in the NHR-49 pathway to promote oxidative stress and hypoxia survival and normal lifespan

NHR-49 likely dimerizes to exert its function, but NHR-49’s dimerization partner(s) in the control of stress responses are unknown. nhr-80 is required for the normal and extended lifespans of wild-type and glp-1 mutant C. elegans, respectively42,60, dimerizes with NHR-49, and both genes are required to activate fatty acid desaturase genes43, which are implicated in longevity13. Thus, although nhr-80 was not a hit in our tBOOH screen (Fig. 1A), its link to nhr-49 function made it an interesting candidate for further investigation. Indeed, we found that nhr-80 is required for fmo-2 induction in both hypoxia and starvation (Fig. 1A). To determine if nhr-80 is required to promote oxidative stress survival, wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) L4 animals were exposed to 6 mM tBOOH. All mutant animals were sensitive to tBOOH compared to wild-type animals, and nhr-49 mutant animals were more sensitive than nhr-80 mutant animals (Fig. 5A). Epistasis analysis revealed that there is no additive effect when both nhr-49 and nhr-80 are lost compared to the nhr-49 single mutant, suggesting that these two factors act in the same pathway in the tBOOH response (Fig. 5A,E). Thus, although nhr-80 is not required to induce fmo-2 in response to oxidative stress, it may regulate other important stress response genes to promote survival to tBOOH with nhr-49.


nhr-80 is required for C. elegans oxidative stress survival, hypoxia survival, and normal lifespan.

(A) Survival plot of wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) L4 animals exposed to 6 mM tBOOH. p-values were derived using the log-rank (Mantel-Cox) test. ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 8. (B) The graph shows the average population survival of wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) embryos following 24 hours exposure to 0.5% O2, then allowed to recover at 21% O2 for 65 hours, and counted as the ability to reach at least L4 stage (n = 4, totaling > 100 individual animals per strain). p-values were derived using ordinary one-way ANOVA corrected for multiple comparisons using the Tukey method. **p < 0.01, ***p < 0.001, ****p < 0.0001. (C) Fractional population recovery to L4 after L1 starvation (in percent, y-axis) over time (days of starvation, x-axis) of wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) animals. Error bars represent SD from at least three independent experiments. p-values were derived by calculating AUC for each strain using one-way ANOVA with Dunnett’s multiple comparisons test. *p < 0.05, **p < 0.01 (n = 4, >100 animals per time point, normalized to 100% viability for each strain). (D) Lifespan of wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041); nhr-80(tm1011) animals. p-values were derived using the log-rank (Mantel-Cox) test. *p < 0.05, **p < 0.01, ****p < 0.0001. Data and statistical analysis for additional repeats in Supplementary Table 9. (E) Working model of how NHR-49 and NHR-80 interact in response to oxidative stress, hypoxia, and for promotion of normal lifespan. ns = not significant.

To determine if nhr-80 is required to promote hypoxia survival, wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) embryos were exposed to 0.5% O2 for 24 hours and then allowed to recover to the L4 or later stage when placed back in normoxia for 65 hours. Similar to nhr-49, loss of nhr-80 rendered animals sensitive to hypoxia (Fig. 5B). Loss of both nhr-49 and nhr-80 showed similar survival to each of the single mutants, suggesting that these two genes act in the same pathway to promote hypoxia survival (Fig. 5B,E).

To determine if nhr-80 acts in the same pathway as nhr-49 to promote recovery from starvation, wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) L1 stage animals were starved and tested for ability to reach the L4 stage after recovering on food for 48 hours. Unlike nhr-49, nhr-80 animals were not sensitive to starvation, recovering like wild type (Fig. 5C).

Both nhr-49 and nhr-80 are required for normal lifespan in C. elegans, but whether they interact genetically in this context is unknown27,28,42,60. Thus, we carried out lifespan assays of wild-type, nhr-49(nr2041), nhr-80(tm1011), and nhr-49(nr2041);nhr-80(tm1011) animals (Fig. 5D). Consistent with previous findings, the nhr-49 and nhr-80 single mutants showed a shortened lifespan compared to the wild-type control, with nhr-49 mutant animals having a shorter lifespan than nhr-80. Interestingly, the nhr-49;nhr-80 double loss-of-function mutant had a similar lifespan to the nhr-49 single mutant animals, suggesting that these two genes act in the same genetic pathway to promote normal C. elegans lifespan (Fig. 5D,E).

Discussion

NHR-49 was first identified as an important regulator of lipid metabolism, through which it affects development and lifespan. More recently, additional roles of NHR-49 emerged in the control of stress responses, including in the response to tBOOH-induced oxidative stress, hypoxia, starvation, and innate immune responses after pathogen infection24. However, several aspects of NHR-49-mediated stress resilience remain poorly understood, including how it is activated in stress and how it achieves its effects on gene regulation, especially considering that it orchestrates different transcriptomic responses in different conditions. Here, we use a series of reverse genetic RNAi screens to identify 80 new genes that are required to induce the nhr-49-dependent reporter fmo-2p::gfp in at least one stress, including genes encoding for kinases, transcription factors, and coregulators. Functional follow-up studies of three such genes, the kinase hpk-1, the pseudokinase nipi-3, and the nuclear hormone receptor nhr-80, revealed that they feature nonequivalent contributions to NHR-49-driven stress resistance and longevity. Collectively, our work identifies numerous new genes, some of which may act with NHR-49 to implement condition-specific gene programs, representing a substantial advance in uncovering the composition of important cytoprotective and pro-longevity pathways.

fmo-2 is variably required for stress resistance

We used the fmo-2 promoter as a basis for our screen readout because it is highly induced by NHR-49 in various stresses23,3032,61. The FMO family is conserved and present in all eukaryotes examined so far62. Although S. cerevisiae has only one FMO gene, C. elegans, mice, and humans all contain five due to evolutionary duplication events, possibly due to selection pressure caused by increased xenobiotic exposure63. The five C. elegans fmo genes are paralogs and are most closely related to mammalian FMO564. Here, we exploit the strong induction of fmo-2 by using the transgenic fmo-2p::gfp strain as a transcriptional reporter in oxidative stress, hypoxia, starvation, and in the nhr-49(et13) background to identify potential genes acting within NHR-49-controlled stress response pathways. However, although fmo-2 expression is strongly induced in many stresses, its functional role is not entirely clear. Interestingly, FMOs regulate metabolism, and overexpression of any of five mouse Fmo genes render cells resistant to multiple stressors and affect cellular metabolism by increasing mitochondrial respiration and decreasing glycolytic activity62. In Fmo5 knockout mice, glycolysis and fatty acid metabolism were downregulated, and mice were leaner and exhibited higher energy expenditure, suggesting roles for Fmo5 in regulating metabolism and energy balance65,66. In addition, Fmo3 regulates cholesterol, glucose, and lipid metabolism in mice67. Interestingly, there is a significant overlap between metabolism genes regulated by Fmo3 and those regulated by Peroxisome proliferator-activated receptor alpha (PPARα), a mammalian NHR-49 homolog, suggesting that Fmo3 may regulate cellular metabolism via this pathway67. Thus, it is possible that the C. elegans fmo genes act with NHR-49 to regulate both cellular detoxification and metabolism. However, although it is highly induced by stress, loss of fmo-2 causes only weak sensitivity to starvation, doesn’t cause sensitivity to hypoxia, and causes resistance to tBOOH32. Additional C. elegans fmo genes are induced in stress, including fmo-1 and fmo-3 during starvation, and fmo-1, -3, and -4 during oxidative stress6870. Perhaps these genes act redundantly, with multiple fmo genes cooperating to achieve stress resistance. Functional analysis of combinatorial mutants in the fmo family genes may elucidate if this family of genes is important for stress resistance in C. elegans. In addition, we highlight that candidates from our regulatory screen will require detailed further analysis for functional roles in various stress and longevity contexts, as done herein for hpk-1, nipi-3, and nhr-80.

Several candidate nhr-49 pathway genes play known roles in stress responses and/or longevity

Notably, the 80 screen hits include many genes with known roles in stress responses or the regulation of animal life span, often associated with stress response programs. abl-1 is a tyrosine kinase required for apoptosis during oxidative stress71. Another kinase, gskl-2, becomes transcriptionally induced during dietary restriction by the transcription factor BLMP-1, which also regulates a HIF-1-independent hypoxia response pathway55,72. Whether these genes are in the same pathway in the hypoxia (or any other) response will be interesting to test in the future. In addition, many of the 55 coregulator hits are part of the transcription-export (TREX) complex, which couples transcription to nuclear mRNA export and ensures stress response gene transcripts are delivered to the cytoplasm for translation73. Interaction of NHR-49, and possibly of other cytoprotective transcription factors, with this machinery may represent an important mechanism of ensuring preferential export of stress resilience transcripts.

Homeodomain interacting protein kinase HPK-1 emerges as key contributor in NHR-49-regulated stress response pathways

Our screen identified the nuclear serine/threonine protein kinase hpk-1 as a regulator of fmo-2p::gfp in tBOOH, hypoxia, starvation, and in the nhr-49(et13) background. Expanding on our earlier work on nhr-49 and hpk-1 in hypoxia31, we show that hpk-1 is required in the nhr-49 pathway for survival on tBOOH as well as for the stress-induced increase in NHR-49::GFP levels (Figs. 2 and 3). These findings agree with previous reports showing that hpk-1 is required for C. elegans survival in the presence of the cytochrome oxidase inhibitor sodium azide, which causes oxidative stress53. Furthermore, hpk-1 regulates a suite of Phase I and Phase II detoxification genes in response to oxidative stress, including cytochrome P450s, short-chain dehydrogenases, and glutathione S transferases, which nhr-49 also regulates to promote redox homeostasis32,53. Interestingly, HPK-1 does not seem to be induced in tBOOH, but expression was only examined in the neurons and hypodermis54. It is therefore possible that HPK-1 acts in a different tissue, potentially the intestine, during this stress to exert its functions.

Notably, the function of HIPKs in redox responses is conserved in mammals, where HIPK2 is de-SUMOylated during ROS accumulation, which relieves it from proteasomal degradation. HIPK2 is then acetylated by p300/CBP, which increases cellular survival74. Interestingly, HIPK2 is regulated at the transcriptional level by Nrf2, a master regulator of redox homeostasis75, which induces expression of HIPK2 when cells are exposed to hydrogen peroxide. Downstream, HIPK2 promotes Nrf2 nuclear localization (potentially directly or indirectly) to promote expression of antioxidant genes such as NAD(P)H dehydrogenase [quinone] 1 (NQO1) and heme oxygenase 1 (HO-1)75. Future work may assess whether HIPK2 similarly regulates the NHR-49 homologs PPARα or hepatocyte nuclear factor 4 alpha (HNF4α) during stress in mammals.

Our data show that nhr-49 and hpk-1 function in the same response pathway to tBOOH and hypoxia, but whether or how they physically interact is not known. As HPK-1 phosphorylates transcription factors, such as p53 and CREB76,77, it will be interesting to determine if HPK-1 binds to and/or phosphorylates NHR-49.

The pseudokinase NIPI-3 interacts with the NHR-49 pathway differently in various stress responses

nipi-3 is required for fmo-2p::gfp induction in tBOOH, hypoxia, starvation, and in the nhr-49(et13) background. nipi-3 is a homolog of the highly conserved tribbles family of pseudokinases, which contain structural elements of canonical kinases but cannot phosphorylate substrates78. Instead, tribbles pseudokinases likely act as “decoys” to disrupt normal kinase function or adaptors that facilitate protein degradation79. Interestingly, nipi-3/TRIB3 has conserved roles in many stress response pathways, including starvation, hypoxia, and following pathogen infection34,58,80,81. In C. elegans, nipi-3 is required for survival upon pathogen infection, including by the fungus D. coniospora and the bacteria P. aeruginosa and E. faecalis; the latter two responses are also controlled by NHR-4923,26,34,57,58. nipi-3 functions upstream of SKN-1 during these infections, controlling activation of the MAPK pathway34,58. This function is conserved, as human TRIB3 functions upstream of MAPK in humans, too82. As the MAPK pathway appears to play a minor role in NHR-49-controlled gene expression32, future experiments may uncover if NIPI-3 acts upstream to regulate both NHR-49 and SKN-1 in these stress responses, or if NIPI-3 controls pathways in parallel to NHR-49.

NHR-80 may regulate a subset of NHR-49-dependent metabolism genes to selectively promote stress resistance

In vitro assays show that NHR-49 can interact with numerous other transcription factors, including many NHRs; in addition, NHR-49 also homodimerizes in these assays40,41,43. In vivo, nhr-49 and nhr-80 both regulate fatty acid desaturase genes to promote longevity28,60, suggesting that they interact genetically in this context. Overall, this raises the interesting hypothesis that NHR-49 engages in specific NHR dimerization partnerships to implement the regulation of specific gene sets and thus biological processes. In support of this model, NHR-79 appears to regulate peroxisomal genes together with NHR-4943,44, where as NHR-66 may cooperate with NHR-49 to repress sphingolipid metabolism genes. Our data indicate that nhr-80 and nhr-49 lie in one pathway, not only in the regulation of longevity but also in the hypoxia and oxidative stress response. In contrast, nhr-80 was dispensable for starvation resistance. As such, it is possible that, although nhr-80, like nhr-49, regulates lipid metabolism genes essential for longevity, such as the fatty acid desaturases, the nhr-49-dependent genes promoting starvation survival (e.g., acs-2) are not dependent on nhr-80. Further fine mapping of how various nhr genes contribute to nhr-49-dependent regulation may reveal additional specificity determinants in the future.

Acknowledgments

This work was supported by grants from the Canadian Institutes of Health Research (CIHR; PJT-153199, PJT-186144 to S.T.) and the Natural Sciences and Engineering Research Council of Canada (NSERC; RGPIN-2018-05133 and RGPIN-2024-06537 to S.T.). K.R.S.D. was supported by NSERC CGS-D, UBC 4-year fellowship (4YF), and BC Children’s Hospital Research Institute (BCCHR) scholarships, G.E. by CIHR CGS-D, UBC Medical Genetics Graduate Program, and UBC 4-year fellowship (4YF) scholarships, B.S. by a BCCHR Healthy Starts Summer Studentship, L.K. by an NSERC Undergraduate Summer Research Award (USRA), and S.T. by a Canada Research Chair and a BCCHR Investigator Grant Award Program (IGAP) award. Some strains were provided by the CGC, which is funded by the NIH Office of Research Infrastructure Programs (P40 OD010440).

Author Contributions

Design of research: K.R.S.D., G.E., S.T.; Performed experiments: K.R.S.D., G.E., B.S., R.C., L.K., Z.B.; Analyzed data: K.R.S.D., G.E.; Wrote the manuscript: K.R.S.D., G.E., S.T.; Edited the manuscript: All authors.

Conflict of Interest statement

The authors declare that they have no actual or apparent conflict of interest between authorship of this study and any other activities.

Supplementary Materials

Supplemental information can be found here: Supplementary.

References

Halliwell B. (1991). Reactive oxygen species in living systems: Source, biochemistry, and role in human disease. Am. J. Med. 91(3), S14S22. PMID: 1928205; doi: 10.1016/0002-9343(91)90279-7.

Shields H.J., Traa A., & Van Raamsdonk J.M. (2021). Beneficial and detrimental effects of reactive oxygen species on lifespan: A comprehensive review of comparative and experimental studies. Front. Cell. Dev. Biol. 9, 628157. PMID: 33644065; doi: 10.3389/fcell.2021.628157.

Sies H., Belousov V.V., Chandel N.S., Davies M.J., Jones D.P., Mann G.E., … Winterbourn C. (2022). Defining roles of specific reactive oxygen species (ROS) in cell biology and physiology. Nat. Rev. Mol. Cell Biol. 23, 499515. PMID: 35190722; doi: 10.1038/s41580-022-00456-z.

Lee P., Chandel N.S., & Simon M.C. (2020). Cellular adaptation to hypoxia through hypoxia inducible factors and beyond. Nat. Rev. Mol. Cell. Biol. 21, 268283. PMID: 32144406; doi: 10.1038/s41580-020-0227-y.

Powell-Coffman J.A. (2010). Hypoxia signaling and resistance in C. elegans. Trends Endocrinol. Metabol. TEM 21(7), 435440. PMID: 20335046; doi: 10.1016/j.tem.2010.02.006.

Harvald E.B., Sprenger R.R., Dall K.B., Ejsing C.S., Nielsen R., Mandrup S., … Færgeman N.J. (2017). Multi-omics analyses of starvation responses reveal a central role for lipoprotein metabolism in acute starvation survival in C. elegans. Cell Syst. 5(1), 3852. PMID: 28734827; doi: 10.1016/j.cels.2017.06.004.

Liu K. & Czaja M.J. (2013). Regulation of lipid stores and metabolism by lipophagy. Cell Death Differ 20, 311. PMID: 22595754; doi: 10.1038/cdd.2012.63.

Settembre C., De Cegli R., Mansueto G., Saha P.K., Vetrini F., Visvikis O., … Ballabio A. (2013). TFEB controls cellular lipid metabolism through a starvation-induced autoregulatory loop. Nat. Cell. Biol. 15, 647658. PMID: 23604321; doi: 10.1038/ncb2718.

An J.H. & Blackwell T.K. (2003). SKN-1 links C. elegans mesendodermal specification to a conserved oxidative stress response. Genes. Dev. 17, 18821893. PMID: 12869585; doi: 10.1101/gad.1107803.

An J.H., Vranas K., Lucke M., Inoue H., Hisamoto N., Matsumoto K., & Blackwell T.K. (2005). Regulation of the Caenorhabditis elegans oxidative stress defense protein SKN-1 by glycogen synthase kinase-3. Proc. Natl. Acad. Sci. U.S.A. 102(45), 1627516280. PMID: 16251270; doi: 10.1073/pnas.0508105102.

Blackwell T.K., Steinbaugh M.J., Hourihan J.M., Ewald C.Y., & Isik M. (2015). SKN-1/Nrf, stress responses, and aging in Caenorhabditis elegans . Free Radic. Biol. Med. 88(Part B), 290301. PMID: 26232625; doi: 10.1016/j.freeradbiomed.2015.06.008.

Oliveira R.P., Porter Abate J., Dilks K., Landis J., Ashraf J., Murphy C.T., & Blackwell T.K. (2009). Condition-adapted stress and longevity gene regulation by Caenorhabditis elegans SKN-1/Nrf. Aging Cell 8(5), 524541. PMID: 19575768; doi: 10.1111/j.1474-9726.2009.00501.x.

Steinbaugh M.J., Narasimhan S.D., Robida-Stubbs S., Moronetti Mazzeo L.E., Dreyfuss J.M., Hourihan J.M., … Blackwell T.K. (2015). Lipid-mediated regulation of SKN-1/Nrf in response to germ cell absence. Elife 4, e07836. PMID: 26196144; doi: 10.7554/eLife.07836.

Tullet J.M.A., Hertweck M., An J.H., Baker J., Hwang J.Y., Liu S., … Blackwell T.K. (2008). Direct inhibition of the longevity-promoting factor SKN-1 by insulin-like signaling in C. elegans. Cell 132(6), 10251038. PMID: 18358814; doi: 10.1016/j.cell.2008.01.030.

Epstein A.C., Gleadle J.M., McNeill L.A., Hewitson K.S., O’Rourke J., Mole D.R., … Ratcliffe P.J. (2001). C. elegans EGL-9 and mammalian homologs define a family of dioxygenases that regulate HIF by prolyl hydroxylation. Cell 107(1), 4354. PMID: 11595184; doi: 10.1016/s0092-8674(01)00507-4.

Choudhry H. & Harris A.L. (2018). Advances in Hypoxia-inducible factor biology. Cell Metab. 27(2), 281298. PMID: 29129785; doi: 10.1016/j.cmet.2017.10.005.

Lapierre L.R., De Magalhaes Filho C.D., McQuary P.R., Chu C.-C., Visvikis O., Chang J.T., … Hansen M. (2013). The TFEB orthologue HLH-30 regulates autophagy and modulates longevity in Caenorhabditis elegans . Nat. Commun. 4, 2267. PMID: 23925298; doi: 10.1038/ncomms3267.

Ma Q. (2013). Role of nrf2 in oxidative stress and toxicity. Annu. Rev. Pharmacol. Toxicol 53, 401426. PMID: 23294312; doi: 10.1146/annurev-pharmtox-011112-140320.

Wang L., Cui S., Ma L., Kong L., & Geng X. (2015). Current advances in the novel functions of hypoxia-inducible factor and prolyl hydroxylase in invertebrates. Insect. Mol. Biol. 24(6), 634648. PMID: 26387499; doi: 10.1111/imb.12189.

O’Rourke E.J. & Ruvkun G. (2013). MXL-3 and HLH-30 transcriptionally link lipolysis and autophagy to nutrient availability. Nat. Cell. Biol. 15, 668676. PMID: 23604316; doi: 10.1038/ncb2741.

Uno M. & Nishida E. (2016). Lifespan-regulating genes in C. elegans. NPJ Aging Mech. Dis. 2, 16010. PMID: 28721266; doi: 10.1038/npjamd.2016.10.

Burkewitz K., Morantte I., Weir H.J.M., Yeo R., Zhang Y., Huynh F.K., … Mair W.B. (2015). Neuronal CRTC-1 governs systemic mitochondrial metabolism and lifespan via a catecholamine signal. Cell 160(5), 842855. PMID: 25723162; doi: 10.1016/j.cell.2015.02.004.

Dasgupta M., Shashikanth M., Gupta A., Sandhu A., De A., Javed S., & Singh V. (2020). NHR-49 transcription factor regulates immunometabolic response and survival of Caenorhabditis elegans during enterococcus faecalis infection. Infect. Immun. 88(8), 101128. PMID: 32482643; doi: 10.1128/IAI.00130-20.

Doering K.R.S., Ermakova G., & Taubert S. (2023). Nuclear hormone receptor NHR-49 is an essential regulator of stress resilience and healthy aging in Caenorhabditis elegans . Front Physiol. 14, 1241591. PMID: 37645565; doi: 10.3389/fphys.2023.1241591.

Hu Q., D’Amora D.R., MacNeil L.T., Walhout A.J.M., & Kubiseski T.J. (2018). The Caenorhabditis elegans oxidative stress response requires the NHR-49 transcription factor. G3 (Bethesda) 8(12), 38573863. PMID: 30297383; doi: 10.1534/g3.118.200727.

Naim N., Amrit F.R.G., Ratnappan R., DelBuono N., Loose J.A., & Ghazi A. (2021). Cell nonautonomous roles of NHR-49 in promoting longevity and innate immunity. Aging Cell 20(7), e13413. PMID: 34156142; doi: 10.1111/acel.13413.

Ratnappan R., Amrit F.R.G., Chen S.-W., Gill H., Holden K., Ward J., … Ghazi A. (2014). Germline signals deploy NHR-49 to modulate fatty-acid β-oxidation and desaturation in somatic tissues of C. elegans. PLoS Genet. 10, e1004829. PMID: 25474470; doi: 10.1371/journal.pgen.1004829.

Van Gilst M.R., Hadjivassiliou H., Jolly A., & Yamamoto K.R. (2005). Nuclear hormone receptor NHR-49 controls fat consumption and fatty acid composition in C. elegans. PLoS Biol. 3, e53. PMID: 15719061; doi: 10.1371/journal.pbio.0030053.

Vozdek R., Long Y., & Ma D.K. (2018). The receptor tyrosine kinase HIR-1 coordinates HIF-independent responses to hypoxia and extracellular matrix injury. Sci. Signal 11(550), 117138. PMID: 30279166; doi: 10.1126/scisignal.aat0138.

Wani K.A., Goswamy D., Taubert S., Ratnappan R., Ghazi A., & Irazoqui J.E. (2021). NHR-49/PPAR-α and HLH-30/TFEB cooperate for C. elegans host defense via a flavin-containing monooxygenase. Elife 10, e62775. PMID: 33978570; doi: 10.7554/eLife.62775.

Doering K.R.S., Cheng X., Milburn L., Ratnappan R., Ghazi A., Miller D.L., & Taubert S. (2022). Nuclear hormone receptor NHR-49 acts in parallel with HIF-1 to promote hypoxia adaptation in Caenorhabditis elegans . Elife 11, e67911. PMID: 35285794; doi: 10.7554/eLife.67911.

Goh G.Y.S., Winter J.J., Bhanshali F., Doering K.R.S., Lai R., Lee K., … Taubert S. (2018). NHR-49/HNF4 integrates regulation of fatty acid metabolism with a protective transcriptional response to oxidative stress and fasting. Aging Cell 17(3), e12743. PMID: 29508513; doi: 10.1111/acel.12743.

Inoue H., Hisamoto N., An J.H., Oliveira R.P., Nishida E., Blackwell T.K., & Matsumoto K. (2005). The C. elegans p38 MAPK pathway regulates nuclear localization of the transcription factor SKN-1 in oxidative stress response. Genes. Dev. 19, 22782283. PMID: 16166371; doi: 10.1101/gad.1324805.

Wu C., Karakuzu O., & Garsin D.A. (2021). Tribbles pseudokinase NIPI-3 regulates intestinal immunity in Caenorhabditis elegans by controlling SKN-1/Nrf activity. Cell Rep. 36(7), 109529. PMID: 34407394; doi: 10.1016/j.celrep.2021.109529.

Choe K.P., Przybysz A.J., & Strange K. (2009). The WD40 repeat protein WDR-23 functions with the CUL4/DDB1 ubiquitin ligase to regulate nuclear abundance and activity of SKN-1 in Caenorhabditis elegans . Mol. Cell. Biol. 29(10), 27042715. PMID: 19273594; doi: 10.1128/MCB.01811-08.

Wu C.-W., Deonarine A., Przybysz A., Strange K., & Choe K.P. (2016). The Skp1 homologs SKR-1/2 are required for the Caenorhabditis elegans SKN-1 antioxidant/detoxification response independently of p38 MAPK. PLoS Genet. 12, e1006361. PMID: 27776126; doi: 10.1371/journal.pgen.1006361.

Ganner A., Gerber J., Ziegler A.-K., Li Y., Kandzia J., Matulenski T., … Neumann-Haefelin E. (2019). CBP-1/p300 acetyltransferase regulates SKN-1/Nrf cellular levels, nuclear localization, and activity in C. elegans. Exp. Gerontol 126, 110690. PMID: 31419472; doi: 10.1016/j.exger.2019.110690.

Goh G.Y.S., Martelli K.L., Parhar K.S., Kwong A.W.L., Wong M.A., Mah A., …, Taubert S. (2014). The conserved Mediator subunit MDT-15 is required for oxidative stress responses in Caenorhabditis elegans . Aging Cell 13(1), 7079. PMID: 23957350; doi: 10.1111/acel.12154.

Arda H.E., Taubert S., MacNeil L.T., Conine C.C., Tsuda B., Van Gilst M., … Walhout A.J.M. (2010). Functional modularity of nuclear hormone receptors in a Caenorhabditis elegans metabolic gene regulatory network. Mol. Syst. Biol. 6(1), 367. PMID: 20461074; doi: 10.1038/msb.2010.23.

Reece-Hoyes J.S., Pons C., Diallo A., Mori A., Shrestha S., Kadreppa S., … Walhout A.J.M. (2013). Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network. Molecular Cell. 51(1), 116127. PMID: 23791784; doi: 10.1016/j.molcel.2013.05.018.

Taubert S., Van Gilst M.R., Hansen M., & Yamamoto K.R. (2006). A Mediator subunit, MDT-15, integrates regulation of fatty acid metabolism by NHR-49-dependent and -independent pathways in C. elegans. Genes. Dev. 20, 11371149. PMID: 16651656; doi: 10.1101/gad.1395406.

Folick A., Oakley H.D., Yu Y., Armstrong E.H., Kumari M., Sanor L., … Wang M.C. (2015). Aging. Lysosomal signaling molecules regulate longevity in Caenorhabditis elegans . Science 347(6217), 8386. PMID: 25554789; doi: 10.1126/science.1258857.

Pathare P.P., Lin A., Bornfeldt K.E., Taubert S., & Van Gilst M.R. (2012). Coordinate regulation of lipid metabolism by novel nuclear receptor partnerships. PLoS Genet. 8, e1002645. PMID: 22511885; doi: 10.1371/journal.pgen.1002645.

Zeng L., Li X., Preusch C.B., He G.J., Xu N., Cheung T.H., … Mak H.Y. (2021). Nuclear receptors NHR-49 and NHR-79 promote peroxisome proliferation to compensate for aldehyde dehydrogenase deficiency in C. elegans . PLoS Genet. 17, e1009635. PMID: 34237064; doi: 10.1371/journal.pgen.1009635.

Jeong J.-H., Han J.-S., Jung Y., Lee S.-M., Park S.-H., Park M., … Kwon E.-S. (2023). A new AMPK isoform mediates glucose-restriction induced longevity non-cell autonomously by promoting membrane fluidity. Nat. Commun. 14, 288. PMID: 36653384; doi: 10.1038/s41467-023-35952-z.

Moreno-Arriola E., El Hafidi M., Ortega-Cuéllar D., & Carvajal K. (2016). AMP-activated protein Kinase regulates oxidative metabolism in Caenorhabditis elegans through the NHR-49 and MDT-15 transcriptional regulators. PLoS One 11, e0148089. PMID: 26824904; doi: 10.1371/journal.pone.0148089.

Fox B.W., Helf M.J., Burkhardt R.N., Artyukhin A.B., Curtis B.J., Palomino D.F., … Schroeder F.C. (2024). Evolutionarily related host and microbial pathways regulate fat desaturation in C. elegans. Nat. Commun. 15, 1520. PMID: 38374083; doi: 10.1038/s41467-024-45782-2.

Porta-de-la-Riva M., Fontrodona L., Villanueva A., & Cerón J. (2012). Basic Caenorhabditis elegans methods: synchronization and observation. J. Vis. Exp. e4019. PMID: 22710399; doi: 10.3791/4019.

Manning G. (2005). Genomic Overview of Protein Kinases. WormBook. PMID: 18050405; doi: 10.1895/wormbook.1.60.1.

Kamath R.S. & Ahringer J. (2003). Genome-wide RNAi screening in Caenorhabditis elegans . Methods 30(4), 313321. PMID: 12828945; doi: 10.1016/s1046-2023(03)00050-1.

Rual J.-F., Ceron J., Koreth J., Hao T., Nicot A.-S., Hirozane-Kishikawa T., … Vidal M. (2004). Toward improving Caenorhabditis elegans phenome mapping with an ORFeome-based RNAi library. Genome. Res. 14, 21622168. PMID: 15489339; doi: 10.1101/gr.2505604.

Fawcett E.M., Horsman J.W., & Miller D.L. (2012). Creating defined gaseous environments to study the effects of hypoxia on C. elegans. J. Vis. Exp. e4088. PMID: 22850348; doi: 10.3791/4088.

Berber S., Wood M., Llamosas E., Thaivalappil P., Lee K., Liao B.M., … Nicholas H.R. (2016). Homeodomain-Interacting Protein Kinase (HPK-1) regulates stress responses and ageing in C. elegans. Sci. Rep. 6, 19582. PMID: 26791749; doi: 10.1038/srep19582.

Das R., Melo J.A., Thondamal M., Morton E.A., Cornwell A.B., Crick B., … Samuelson A.V. (2017). The homeodomain-interacting protein kinase HPK-1 preserves protein homeostasis and longevity through master regulatory control of the HSF-1 chaperone network and TORC1-restricted autophagy in Caenorhabditis elegans . PLoS Genet 13, e1007038. PMID: 29036198; doi: 10.1371/journal.pgen.1007038.

Hu Q., Xu Y., Song M., Dai Y., Antebi A., & Shen Y. (2024). BLMP-1 is a critical temporal regulator of dietary-restriction-induced response in Caenorhabditis elegans . Cell Rep. 43(3), 113959. PMID: 38483903; doi: 10.1016/j.celrep.2024.113959.

Jiang H., Guo R., & Powell-Coffman J.A. (2001). The Caenorhabditis elegans hif-1 gene encodes a bHLH-PAS protein that is required for adaptation to hypoxia. Proc. Nat. Acad. Sci. U. S. A. 98(14), 79167921. PMID: 11427734; doi: 10.1073/pnas.141234698.

McEwan D.L., Feinbaum R.L., Stroustrup N., Haas W., Conery A.L., Anselmo A., … Ausubel F.M. (2016). Tribbles ortholog NIPI-3 and bZIP transcription factor CEBP-1 regulate a Caenorhabditis elegans intestinal immune surveillance pathway. BMC Biol. 14, 105. PMID: 27927200; doi: 10.1186/s12915-016-0334-6.

Pujol N., Cypowyj S., Ziegler K., Millet A., Astrain A., Goncharov A., … Ewbank J.J. (2008). Distinct innate immune responses to infection and wounding in the C. elegans epidermis. Curr. Biol. 18(7), 481489. PMID: 18394898; doi: 10.1016/j.cub.2008.02.079.

Lee K., Goh G.Y.S., Wong M.A., Klassen T.L., & Taubert S. (2016). Gain-of-function alleles in Caenorhabditis elegans nuclear hormone receptor nhr-49 are functionally distinct. PloS One 11, e0162708. PMID: 27618178; doi: 10.1371/journal.pone.0162708.

Goudeau J., Bellemin S., Toselli-Mollereau E., Shamalnasab M., Chen Y., & Aguilaniu H. (2011). Fatty acid desaturation links germ cell loss to longevity through NHR-80/HNF4 in C. elegans. PLoS Biol 9, e1000599. PMID: 21423649; doi: 10.1371/journal.pbio.1000599.

Leiser S.F., Miller H., Rossner R., Fletcher M., Leonard A., Primitivo M., … Kaeberlein M. (2015). Cell nonautonomous activation of flavin-containing monooxygenase promotes longevity and health span. Science (New York, N.Y.) 350, 13751378. PMID: 26586189; doi: 10.1126/science.aac9257.

Huang S., Howington M.B., Dobry C.J., Evans C.R., & Leiser S.F. (2021). Flavin-containing monooxygenases are conserved regulators of stress resistance and metabolism. Front. Cell Dev. Biol. 9, 630188. PMID: 33644069; doi: 10.3389/fcell.2021.630188.

Hao D.C., Chen S.L., Mu J., & Xiao P.G. (2009). Molecular phylogeny, long-term evolution, and functional divergence of flavin-containing monooxygenases. Genetica 137(2), 173187. PMID: 19579011; doi: 10.1007/s10709-009-9382-y.

Nicoll C.R., Bailleul G., Fiorentini F., Mascotti M.L., Fraaije M.W., & Mattevi A. (2020). Ancestral-sequence reconstruction unveils the structural basis of function in mammalian FMOs. Nat. Struct. Mol. Biol. 27, 1424. PMID: 31873300; doi: 10.1038/s41594-019-0347-2.

Gonzalez Malagon S.G., Melidoni A.N., Hernandez D., Omar B.A., Houseman L., Veeravalli S., … Shephard E.A. (2015). The phenotype of a knockout mouse identifies flavin-containing monooxygenase 5 (FMO5) as a regulator of metabolic ageing. Biochem. Pharmacol. 96(3), 267277. PMID: 26049045; doi: 10.1016/j.bcp.2015.05.013.

Scott F., Gonzalez Malagon S.G., O’Brien B.A., Fennema D., Veeravalli S., Coveney C.R., … Shephard E.A. (2017). Identification of flavin-containing monooxygenase 5 (FMO5) as a regulator of glucose homeostasis and a potential sensor of Gut Bacteria. Drug. Metab. Dispos. 45(9), 982989. PMID: 28646079; doi: 10.1124/dmd.117.076612.

Shih D.M., Wang Z., Lee R., Meng Y., Che N., Charugundla S., … Lusis A.J. (2015). Flavin containing monooxygenase 3 exerts broad effects on glucose and lipid metabolism and atherosclerosis. J. Lipid Res. 56(1), 2237. PMID: 25378658; doi: 10.1194/jlr.M051680.

Nhan J.D., Turner C.D., Anderson S.M., Yen C.-A., Dalton H.M., Cheesman H.K., … Curran S.P. (2019). Redirection of SKN-1 abates the negative metabolic outcomes of a perceived pathogen infection. Proc. Nat. Acad. Sci. U S A 116(44), 2232222330. PMID: 31611372; doi: 10.1073/pnas.1909666116.

Park S.-K., Tedesco P.M., & Johnson T.E. (2009). Oxidative stress and longevity in Caenorhabditis elegans as mediated by SKN-1. Aging Cell 8(3), 258269. PMID: 19627265; doi: 10.1111/j.1474-9726.2009.00473.x.

Uno M., Honjoh S., Matsuda M., Hoshikawa H., Kishimoto S., Yamamoto T., … Nishida E. (2013). A Fasting-Responsive Signaling Pathway that Extends Life Span in C. elegans. Cell Rep. 3(1), 7991. PMID: 23352664; doi: 10.1016/j.celrep.2012.12.018.

Salinas L.S., Maldonado E., & Navarro R.E. (2006). Stress-induced germ cell apoptosis by a p53 independent pathway in Caenorhabditis elegans . Cell. Death Differ. 13, 21292139. PMID: 16729024; doi: 10.1038/sj.cdd.4401976.

Padmanabha D., Padilla P.A., You Y.-J., & Baker K.D. (2015). A HIF-independent mediator of transcriptional responses to oxygen deprivation in Caenorhabditis elegans . Genetics 199(3), 739748. PMID: 25552276; doi: 10.1534/genetics.114.173989.

Heath C.G., Viphakone N., & Wilson S.A. (2016). The role of TREX in gene expression and disease. Biochem. J. 473(19), 29112935. PMID: 27679854; doi: 10.1042/BCJ20160010.

de la Vega L., Grishina I., Moreno R., Krüger M., Braun T., & Schmitz M.L. (2012). A redox-regulated SUMO/acetylation switch of HIPK2 controls the survival threshold to oxidative stress. Molecular Cell 46(4), 472483. PMID: 22503103; doi: 10.1016/j.molcel.2012.03.003.

Torrente L., Sanchez C., Moreno R., Chowdhry S., Cabello P., Isono K., … de la Vega L. (2017). Crosstalk between NRF2 and HIPK2 shapes cytoprotective responses. Oncogene 36, 62046212. PMID: 28692050; doi: 10.1038/onc.2017.221.

Hofmann T.G., Glas C., & Bitomsky N. (2013). HIPK2: A tumour suppressor that controls DNA damage-induced cell fate and cytokinesis. Bioessays 35(1), 5564. PMID: 23169233; doi: 10.1002/bies.201200060.

Sakamoto K., Huang B.-W., Iwasaki K., Hailemariam K., Ninomiya-Tsuji J., & Tsuji Y. (2010). Regulation of genotoxic stress response by homeodomain-interacting protein kinase 2 through phosphorylation of cyclic AMP response element-binding protein at serine 271. Mol. Biol. Cell. 21(16), 29662974. PMID: 20573984; doi: 10.1091/mbc.E10-01-0015.

Hegedus Z., Czibula A., & Kiss-Toth E. (2007). Tribbles: A family of kinase-like proteins with potent signalling regulatory function. Cellular Signal. 19(2), 238250. PMID: 16963228; doi: 10.1016/j.cellsig.2006.06.010.

Lohan F. & Keeshan K. (2013). The functionally diverse roles of tribbles. Biochem. Soc. Trans. 41(4), 10961100. PMID: 23863185; doi: 10.1042/BST20130105.

Bowers A.J., Scully S., & Boylan J.F. (2003). SKIP3, a novel Drosophila tribbles ortholog, is overexpressed in human tumors and is regulated by hypoxia. Oncogene 22, 28232835. PMID: 12743605; doi: 10.1038/sj.onc.1206367.

Wennemers M., Bussink J., Scheijen B., Nagtegaal I.D., van Laarhoven H.W.M., Raleigh J.A., … Span P.N. (2011). Tribbles homolog 3 denotes a poor prognosis in breast cancer and is involved in hypoxia response. Breast Cancer Res. 13, R82. PMID: 21864376; doi: 10.1186/bcr2934.

Kiss-Toth E., Bagstaff S.M., Sung H.Y., Jozsa V., Dempsey C., Caunt J.C., … Dower S.K. (2004). Human tribbles, a protein family controlling mitogen-activated protein kinase cascades. J. Biol. Chem. 279(41), 4270342708. PMID: 15299019; doi: 10.1074/jbc.M407732200.

Brenner S. (1974). The genetics of Caenorhabditis elegans. Genetics 77(1), 7194. PMID: 4366476; doi: 10.1093/genetics/77.1.71.

Raich W.B., Moorman C., Lacefield C.O., Lehrer J., Bartsch D., Plasterk R.H.A., … Hobert O. (2003). Characterization of Caenorhabditis elegans homologs of the down syndrome candidate gene DYRK1A. Genetics 163(2), 571580. PMID: 12618396; doi: 10.1093/genetics/163.2.571.