NHR-8 and P-glycoproteins uncouple xenobiotic resistance from longevity in chemosensory C. elegans mutants
Abstract
Longevity is often associated with stress resistance, but whether they are causally linked is incompletely understood. Here we investigate chemosensory-defective Caenorhabditis elegans mutants that are long-lived and stress resistant. We find that mutants in the intraflagellar transport protein gene osm-3 were significantly protected from tunicamycin-induced ER stress. While osm-3 lifespan extension is dependent on the key longevity factor DAF-16/FOXO, tunicamycin resistance was not. osm-3 mutants are protected from bacterial pathogens, which is pmk-1 p38 MAP kinase dependent, while TM resistance was pmk-1 independent. Expression of P-glycoprotein (PGP) xenobiotic detoxification genes was elevated in osm-3 mutants and their knockdown or inhibition with verapamil suppressed tunicamycin resistance. The nuclear hormone receptor nhr-8 was necessary to regulate a subset of PGPs. We thus identify a cell-nonautonomous regulation of xenobiotic detoxification and show that separate pathways are engaged to mediate longevity, pathogen resistance, and xenobiotic detoxification in osm-3 mutants.
Introduction
Chemosensation is a genetically tractable phenotype in various organisms. In Caenorhabditis elegans, many mutants with defective chemosensation have been identified. Sensory phenotypes are complex in nature, and many of the classical chemosensory mutants were originally characterized by their behavioral phenotypes. Odr mutants, for example, have an abnormal odorant response, while Osm (osmotic avoidance abnormal) mutants do not avoid high-salt environments (Bargmann, 2006). Mutations in components of neuronal G-protein coupled receptor (GPCR) signaling, such as tax-4 and odr-1, cause atypical chemosensory behavior (Bargmann, 2006). These mutants are not only characterized by a failure to adequately respond to their environment, but show additional phenotypes linked to various life traits including pathogen resistance, increased lifespan, and drug detoxification (Gaglia et al., 2012; Apfeld and Kenyon, 1999; Dent et al., 2000). Many chemosensory mutants are long-lived, and this phenotype depends on the DAF-16/FOXO transcription factor that is regulated by the insulin signaling pathway (Apfeld and Kenyon, 1999).
Beyond mutations in GPCRs, variants in genes involved in the development of amphid sensory neurons can lead to chemosensory defects. One example are mutations in intraflagellar transport (IFT) proteins that prevent full cilia development of amphid sensory neurons (Inglis et al., 2007). Alterations in the IFT genes osm-3 and daf-10 disrupt the development of neuronal dendrites that project toward the tip of the nose, where they are exposed to the outer environment (Inglis et al., 2007). Thus, IFT mutants are usually characterized by defects in the chemical perception of their environment. Furthermore, developmental defects in lumen formation of the amphid head channel in daf-6 mutants prevent the direct contact of amphid sensory neurons with the outside (Perens and Shaham, 2005), resulting in chemosensory defects.
In wild-type (WT) worms, amphid and phasmid sensory neurons can be visualized using the lipophilic dye DiI (1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine perchlorate), which fluorescently stains them bright red. The dye is passively taken up by sensory neurons when they are fully developed and in contact with the environment. Dye filling defective (Dyf) C. elegans are defined by an inability to take up DiI into their sensory neurons (Inglis et al., 2007). The Dyf phenotype arises in daf-6 animals and also in many IFT mutants. Interestingly, Dyf mutants have unique stress resistance phenotypes that act via different signaling pathways. One such Dyf mutant, daf-10(m79), has a unique resistance to pathogenic bacteria that was proposed to be downstream of sensory input (Gaglia et al., 2012). Furthermore, pairs of amphid head neurons, such as ASH and ASJ that often become disrupted in Dyf mutants, help coordinate the innate immune response to bacterial stress after infection (Meisel et al., 2014). Longevity is often associated with activation of physiological stress response pathways and has been shown to be regulated by the insulin signaling pathway in Dyf C. elegans (Apfeld and Kenyon, 1999).
Conserved mechanisms of longevity and stress resistance have additionally been linked to hyposmia in higher model organisms. In Drosophila melanogaster, loss of the putative chemoreceptor Or83b has been shown to significantly increase lifespan (Libert et al., 2007). Moreover, mice with ablated olfactory sensory neurons have a unique metabolic signature that protects them from deleterious effects of a high-fat diet (Riera et al., 2017).
While the olfactory machinery is seemingly more complex in higher organisms, C. elegans have unique phenotypes associated with the loss or alteration of their sensory neurons. First, unbiased screens for drug resistance have independently identified the Dyf phenotype as a common feature in drug-resistant C. elegans (Fujii et al., 2004; Ménez et al., 2016; Collins et al., 2008). Selection for resistance to the anthelmintic drug ivermectin was shown to enrich for Dyf mutants in several forward genetic screens (Ménez et al., 2016; Page, 2018; Dent et al., 2000). Follow-up work identified P-glycoproteins (PGPs), a class of ATP binding cassette (ABC) transporters, as the drivers behind ivermectin resistance in Dyf mutants (Ardelli and Prichard, 2013). Second, some chemosensory mutations, such as in daf-10, lead to resistance to pathogenic bacteria (Gaglia et al., 2012).
Here, we characterize a long-lived Dyf mutant, osm-3, that was previously identified in a forward genetic screen for resistance to the potent ER stressor tunicamycin (TM). From this screen, we reported the role of the hexosamine biosynthetic pathway in TM resistance and longevity (Denzel et al., 2014). Our aim was to uncover the mechanism by which Dyf mutants are resistant to TM. We found that the TM resistance of osm-3 mutants is independent from its longevity and pathogen resistance phenotypes. Furthermore, we showed that PGPs mediate TM resistance in Dyf animals, and we found that the nuclear hormone receptor nhr-8 mediates osm-3 drug resistance. The ‘Green Theory’ of aging proposes that the accumulation of toxic endo- and xenobiotic debris contributes to physiological decline with age (Afschar et al., 2016; Gems and McElwee, 2005). Our work highlights drug resistance as a unique characteristic of Dyf mutants and functionally uncouples it from known stress response pathways that have been implicated in stress resistance and longevity of chemosensory mutants.
Results
Dye filling defective C. elegans mutants are resistant to TM
Tunicamycin (TM) is commonly used to induce endoplasmic reticulum (ER) stress by inhibiting the addition of N-glycans to nascent polypeptides (Parodi, 2000; Heifetz et al., 1979). In C. elegans, where the ER machinery is conserved, treatment with TM is toxic and induces the ER unfolded protein response. In fact, on plates containing TM at concentrations above 4 µg/mL, newly hatched larvae die at early developmental stages (Shen et al., 2001). Previously, we had carried out an ethyl methane-sulfonate mutagenesis screen to identify long-lived mutants using TM resistance in C. elegans as a proxy phenotype for longevity (Denzel et al., 2014). The largest cohort of classifiable mutants presented the Dyf phenotype (Figure 1A).
Among the Dyf mutants from the TM resistance screen, we found new alleles of genes previously linked to ciliary development. The mutant allele we selected for our investigation, osm-3(dh441)IV (this allele will be referred to hereafter as osm-3), carries a premature stop in the eighth exon, which prevents normal cilia development causing the Dyf phenotype (Figure 1A,B). osm-3 is a member of the kinesin family involved in axonal development and transport, and mutants are known to show a Dyf phenotype (Inglis et al., 2007). While WT worms fail to develop on TM at concentrations above 4 µg/mL, osm-3 mutants fully develop at concentrations at least up to 10 µg/mL TM (Figure 1C; Figure 1—source data 1). Loss of osm-3, as well as other Dyf mutations, promote longevity in C. elegans (Apfeld and Kenyon, 1999). Demographic lifespan analysis of osm-3 indeed confirmed extended lifespan in our point mutant (Figure 1D and Supplementary file 1). osm-3 mutants, and Dyf mutants in general, have never been described as TM resistant especially at concentrations as high as 10 µg/mL.
To further explore whether the Dyf phenotype is linked to TM resistance, we tested other Dyf mutants that were previously described as long-lived (Apfeld and Kenyon, 1999). To our surprise, all long-lived Dyf mutants that we tested proved to be significantly TM resistant at the WT lethal dose of 10 µg/mL TM (Figure 1E, Figure 1—figure supplement 1A–G, Figure 1—figure supplement 1—source data 1–7). In contrast, the long-lived chemosensory odr-1(n1936)X, odr-3(n2150)V, and tax-4(p678)III mutants that do not display a Dyf phenotype were only slightly TM resistant (Figure 1E, Figure 1—figure supplement 1H–J, Figure 1—figure supplement 1—source data 8–10). Our data thus suggest that, among the larger class of chemosensory defective mutants, lifespan extension is uncoupled from TM resistance.
Tunicamycin-Induced ER stress response is blunted in Osm-3 mutants
TM treatment activates the ER unfolded protein response (UPRER); we thus characterized the overall transcriptional response of osm-3 mutants upon TM treatment. We performed transcriptomic analyses of osm-3 and WT worms after 6 hr of TM treatment. Notably, the gene ontology (GO) terms related to ER stress and ER protein folding showed significant upregulation in WT animals, but not in osm-3 mutants (Figure 2A). RT-qPCR analysis of UPR target genes further confirmed that there was no UPRER induction in TM-treated osm-3 mutants compared to the WT (Figure 2—figure supplement 1, Figure 2—figure supplement 1—source data 1). A hallmark of UPRER activation in C. elegans is the upregulation of the molecular chaperone HSP-4/BiP (Shen et al., 2001). Upon TM treatment of osm-3 mutants carrying an hsp-4::GFP reporter construct, there was no significant increase in the GFP signal compared to the untreated control, while in the WT, the GFP levels were significantly increased (Figure 2B,C, Figure 2—source data 1). This observation corroborates the results from the transcriptome analysis.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig2-data1-v2.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig2-data2-v2.xlsx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig2-data3-v2.xlsx
One potential explanation for a suppressed UPRER response in osm-3 mutants might be a general defect in the stress signaling pathway. To address this possibility, we sought to activate the UPR by alternate means. The ER membrane protein SEL-1/HRD3 is a cofactor of the HRD1 ubiquitin ligase complex involved in ER-associated degradation (ERAD) (Hampton et al., 1996). Using sel-1 RNAi to inhibit ERAD, we observed a robust activation of the hsp-4::GFP reporter, suggesting a functional UPRER response upon ER stress in osm-3 mutants (Figure 2D,E, Figure 2—source data 2). We further used DTT treatment to induce ER stress and found normal levels of hsp-4::GFP induction and xbp-1 splicing (Figure 2F, Figure 2—source data 3) in osm-3 mutants. Together, these observations rule out that the reduced response to TM observed in osm-3 mutants was due to a defect in UPR signaling.
TM resistance in Osm-3 mutants is independent from Daf-16 or the PMK-1/p38 MAPK pathway
The lifespan extension observed in chemosensory-defective C. elegans as well as in Drosophila has been shown to be at least partially insulin signaling dependent (Apfeld and Kenyon, 1999; Libert et al., 2007). Therefore, we performed a demographic lifespan analysis to determine the role of the insulin signaling pathway in the lifespan extension of osm-3 mutants. Indeed, the osm-3 lifespan extension was fully daf-16 dependent, as the lifespan of the osm-3; daf-16 double mutants was identical to the daf-16 lifespan (Figure 3A, Supplementary file 1). In stark difference to osm-3 survival, the osm-3; daf-16 double mutant remained resistant to TM while daf-16 single mutants do not develop at 10 µg/mL TM (Figure 3B, Figure 3—source data 1). A previous study by Henis-Korenblit et al., 2010 had shown elevated hsp-4::GFP levels during daf-16 knockdown in a daf-2 mutant. We therefore asked whether the knockdown of daf-16 would sensitize osm-3 mutants to TM. Interestingly, knockdown of daf-16 by RNAi did not elevate the hsp-4::GFP response of osm-3 mutants after TM treatment (Figure 3—figure supplement 1, Figure 3—figure supplement 1—source data 1). Given the links between the insulin signaling pathway and ER stress signaling (Henis-Korenblit et al., 2010; Kyriakakis et al., 2017; Matai et al., 2019; Labbadia and Morimoto, 2014), we were surprised to find that daf-2 mutants failed to develop on 10 µg/mL TM, while osm-3 mutants were fully developed (Figure 3—figure supplement 1, Figure 3—figure supplement 1—source data 2). Together, this evidence suggests that the osm-3 lifespan extension and TM drug resistance are uncoupled and act via independent pathways.
In C. elegans, Pseudomonas aeruginosa PA14 is a pathogen that is frequently used to study innate immunity. The PMK-1 p38 mitogen-activated protein kinase (MAPK) pathway as well as UPRER targets have been implicated in the innate immune response of C. elegans during PA14 infection (Richardson et al., 2010; Haskins et al., 2008). Given the role of ER homeostasis in innate immunity and the TM resistance of osm-3 mutants, we speculated that they would be more robust on PA14 than WT. Indeed, osm-3 mutants were resistant to PA14, and this increased pathogen resistance was lost upon pmk-1 knockdown (Figure 3C). With this in mind, we tested whether pmk-1 mutation could suppress the osm-3 TM resistance. The osm-3; pmk-1 double mutant was fully resistant to TM, while the pmk-1 single mutant was sensitive (Figure 3D, Figure 3—source data 2), demonstrating that TM resistance and pmk-1 mediated pathogen resistance are uncoupled in osm-3 mutants. Taken together, insulin signaling and the MAPK pathway, two major stress response pathways in C. elegans that have links to ER protein quality control, are not required for TM resistance in osm-3 mutants.
NHR-8 and PGPs regulate xenobiotic detoxification and TM resistance in osm-3 mutants
Given the absence of ER stress in TM-treated osm-3 mutants along with the genetic separation of TM resistance from insulin or MAPK signaling, we hypothesized that TM might be cleared from worm tissues through xenobiotic detoxification. PGPs are a conserved family of ATP binding cassette (ABC) transporters that locate to the cell membrane (Sangster, 1994) and they are involved in toxin clearance. C. elegans have 15 pgp genes. Using quantitative PCR, we found a significant upregulation of pgp-3, pgp-5, pgp-11, and pgp-13 in osm-3 mutants (Figure 4A). The nuclear hormone receptor NHR-8 has been linked to xenobiotic detoxification and pgp regulation in C. elegans (Ménez et al., 2019; Lindblom et al., 2001). Indeed, we found that the upregulation of pgp-3, pgp-5, pgp-11 and pgp-13, in osm-3 mutants, was suppressed in the osm-3; nhr-8 double mutant (Figure 4A, Figure 4—source data 1). PGP upregulation suggested a role for PGPs in the TM resistance of osm-3 mutants. Therefore, we tested if RNAi-mediated PGP knockdown in osm-3 mutants would suppress TM resistance. We focused on the upregulated pgp-5, pgp-11 and pgp-13, and on pgp-8 and pgp-12. Interestingly, knockdown of pgp-5, pgp-11, and pgp-12 partially reduced the TM resistance of the osm-3 mutants (Figure 4B, Figure 4—source data 2). Next, we used verapamil to inhibit PGPs. Verapamil has been used to specifically inhibit PGP activity, re-sensitizing worms to the anthelmintic compound ivermectin (Ménez et al., 2016). Indeed, 1 nM verapamil significantly suppressed development of osm-3 C. elegans in the presence of TM, showing no effect on controls without TM (Figure 4C, Figure 4—source data 3). The partial effects of pgp knockdown (Figure 4B) and verapamil treatment (Figure 4C) suggest that xenobiotic detoxification is achieved by multiple PGPs, which are likely to be transcriptionally orchestrated by nhr-8. This led us to ask whether osm-3 TM resistance was dependent on NHR-8. Strikingly, the developmental TM assay revealed that the resistance of osm-3 mutants was largely suppressed in the osm-3; nhr-8(hd117)IV double mutant (Figure 4D. Figure 4—source data 4), suggesting that NHR-8 is necessary for osm-3 TM resistance. As suggested by a role of PGPs in ivermectin resistance, we found that osm-3 mutants fully developed to adults in the presence of 6 µg/mL ivermectin that is lethal to WT animals (Figure 4E, Figure 4—source data 5). Ivermectin resistance of osm-3 mutants was suppressed by loss of NHR-8. Dyf C. elegans mutants have been independently identified in drug resistance screens (Fujii et al., 2004; Ménez et al., 2016; Dent et al., 2000). Consistent with these findings, we observed that osm-3 mutants were resistant to 200 nM paraquat, as they fully develop at a concentration that is toxic to WT C. elegans (Figure 4F, Figure 4—source data 6). Consistent with a role in xenobiotic detoxification, NHR-8 was required for the paraquat resistance of the osm-3 mutants (Figure 4F).
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data1-v2.xlsx
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data2-v2.xlsx
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data3-v2.xlsx
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data4-v2.xlsx
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data5-v2.xlsx
-
Figure 4—source data 6
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig4-data6-v2.xlsx
We next reasoned that resistance of osm-3 mutants should be limited to chemical stressors that can be eliminated by xenobiotic detoxification. We thus tested heat stress and hydrogen peroxide (H2O2) induced oxidative stress. Paraquat and hydrogen peroxide are both oxidative stressors, but paraquat acts indirectly via the mitochondria (Castello et al., 2007), while hydrogen peroxide is a primary cause for oxidative damage. When challenged with heat or hydrogen peroxide, osm-3 mutants showed no difference to WT animals, suggesting that osm-3 mutants are not resistant to these non-xenobiotic stressors (Figure 4—figure supplement 1A, B, Figure 4—figure supplement 1—source data 1 and Figure 4—figure supplement 1—source data 2).
Given the involvement of PGP-dependent xenobiotic detoxification in TM resistance of osm-3 mutants, we wondered about its role in the longevity of osm-3 mutants. Notably, verapamil did not shorten osm-3 lifespan (Figure 4G). Moreover, verapamil did not affect WT lifespan (Figure 4G and Supplementary file 1) and had no visible effect on WT development (Figure 4—figure supplement 1C, Figure 4—figure supplement 1—source data 3). Finally, we tested a possible role of NHR-8 in osm-3 longevity. While nhr-8 mutation shortens lifespan, the osm-3 mutation extended both WT and nhr-8 survival to a similar degree (Figure 4H). Together, these data support the conclusion that the xenobiotic stress resistance and longevity phenotypes of osm-3 mutants are mediated by independent pathways.
Sterol signals modulate TM resistance
We aimed to further understand the systemic signals that control xenobiotic detoxification of TM. The NHR-8 ligand remains unknown, but is likely a sterol molecule derived from cholesterol metabolism (Magner et al., 2013). We thus tested the role of cholesterol in TM resistance and found that omitting cholesterol from the culture plates for one generation leads to TM hypersensitivity in both WT and osm-3 animals (Figure 5A, Figure 5—source data 1). To specifically assess a role of cholesterol-derived signals in TM resistance, we used the hsp-4::GFP reporter and found that high cholesterol significantly reduced TM-induced ER stress (Figure 5B, Figure 5—source data 2). Developmental cholesterol dose–response assays in the presence of TM confirmed that WT TM resistance is cholesterol dependent. In nhr-8 and osm-3; nhr-8 double mutants, however, cholesterol treatment failed to protect from TM toxicity (Figure 5C, Figure 5—source data 3). This demonstrates a key role of NHR-8 in mediating TM resistance downstream of sterol signals.
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig5-data1-v2.xlsx
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig5-data2-v2.xlsx
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/53174/elife-53174-fig5-data3-v2.xlsx
Discussion
In this study, we found that osm-3 mutants as well as other long-lived Dyf mutants are resistant to the ER toxin TM. Instead of an activated stress signaling status that might explain the TM resistance, osm-3 Dyf mutants show no UPRER induction upon TM treatment. Consistent with previous findings (Apfeld and Kenyon, 1999), the lifespan extension of osm-3 was daf-16 dependent. Despite the established link between the insulin signaling pathway and stress resistance, the TM resistance in osm-3 mutants was not insulin signaling dependent. Furthermore, osm-3 mutants were pathogen resistant, and this resistance was fully dependent on the PMK-1/p38 MAPK pathway. We were surprised that neither of the two pathways weakened the TM resistance phenotype. Dyf C. elegans drug resistance is not specific to TM, as we also observed resistance to paraquat and ivermectin. While paraquat and ivermectin resistance have been previously reported in Dyf mutants (Dent et al., 2000; Fujii et al., 2004; Ménez et al., 2016), no studies have demonstrated similar resistance mechanisms in Dyf C. elegans mutants to TM. Our findings link drug resistance to increased PGP expression in Dyf worms, suggesting a causal role of xenobiotic detoxification in the drug resistance phenotype. Importantly, we found that the regulation of a subset of PGPs is mediated by the nuclear hormone receptor NHR-8. Genetic or pharmacological PGP inhibition, as well as loss of nhr-8 significantly resensitized osm-3 mutants to TM, ivermectin, and paraquat. Potential NHR-8 ligands are derived from cholesterol (Magner et al., 2013), and we observed that WT development on lower TM concentrations was strongly cholesterol and NHR-8 dependent.
Given that improved protein homeostasis is one of the cellular hallmarks of longevity (López-Otín et al., 2013), we presumed that selection for TM resistance would serve as a proxy phenotype for longevity given that TM specifically targets ER protein folding. More specifically, we expected that high fidelity ER protein quality control would be the driver of longevity in these mutants. Previous findings describing increased hexosamine biosynthetic pathway flux (Denzel et al., 2014) or constitutive XBP-1 activation Taylor and Dillin, 2013 have linked protein homeostasis to longevity and stress resistance (Denzel and Antebi, 2015); therefore, the connection between TM resistance and longevity seemed apparent. Given that Dyf mutants are resistant to pathogenic bacteria through enhanced ER homeostasis, we speculated that TM resistance might serve as a proxy phenotype for pathogen resistance. Contrary to that, we found that TM resistance is a consequence of xenobiotic detoxification via PGPs, and NHR-8 regulated control of PGPs was not involved in the longevity of osm-3 mutants. pmk-1, and daf-16 were required for pathogen resistance and longevity, respectively, but not for the drug resistance of osm-3 mutants. Together, these findings demonstrate that the three phenotypes of osm-3 mutants (i.e. xenobiotic resistance, longevity, and protection from pathogens) are independent and controlled by separable genetic pathways. Our observations thus support a model in which the chemosensation defect activates a number of parallel pathways that are responsible for the distinct phenotypes of osm-3 mutants.
Adding another dimension to the cell-nonautonomous regulation of innate immunity (Aballay, 2013), we propose that the nervous system regulates systemic xenobiotic detoxification in C. elegans, likely through a sterol ligand. Our data suggest that mutants with defective amphid neurons have increased PGP expression. Several PGPs have been shown to be expressed in non-neuronal tissue (Lincke et al., 1993; Sheps et al., 2004). Moreover, we found that PGP expression is regulated by nhr-8, which has also been shown to be expressed in the intestine of C. elegans (Magner et al., 2013). Combining this previous information with our finding that Dyf mutants are resistant to several uniquely toxic drugs, we conclude that neuronal signaling controls drug resistance through nuclear hormone signaling. Based on our results, we cannot make conclusions about the sterol metabolism involved in the production of a putative NHR-8 ligand and we cannot make direct conclusions about PGP regulation by sterols. Given our observation that supplementation with cholesterol is itself sufficient for TM resistance in WT worms, it is now important to investigate the link between neuronal states and the specific cholesterol-derived signal that drives these drug resistance phenotypes.
PGPs are likely regulated in C. elegans to help them combat toxins found in their natural habitat. Toxic metabolites are undoubtedly common in the soil where C. elegans are found. The drugs used in our study, ivermectin (a derivative of avermectin) (Burg et al., 1979) and TM, were first discovered as antibiotics that are produced by soil bacteria (Takatsuki et al., 1971). Because C. elegans have no adaptive immune response, having PGPs as part of their innate immune system allows them to clear toxic molecules they may encounter in the wild. While Dyf mutants likely do not occur in the wild, our data nonetheless demonstrate a link between environmental sensing and drug resistance.
In humans, PGPs have been implicated in drug resistant malignancies (Lehne, 2000). Using C. elegans as a tool to study the cross talk between tissues, one might be able to better understand how extracellular signaling drives ABC transporter expression in chemotherapy-resistant cancer. Further studies in cell culture might also characterize TM, as well as other toxic metabolites, as substrates for C. elegans PGPs.
Materials and methods
Worm maintenance
Request a detailed protocolC. elegans nematodes were cultured using standard husbandry methods at 20 °C on nematode growth media (NGM) agar plates seeded with E. coli strain OP50 unless otherwise stated (Brenner, 1974).
Dye filling assay in C. elegans
Request a detailed protocolThe dye filling assay was performed on a synchronized population of day one adults. Forty to 60 adult worms were placed in M9 containing 10 µg/mL DiI. They were left at room temperature for 2 hr in the staining medium, then transferred to NGM plates seeded with E. coli OP50. The worms were then allowed to crawl on the food for 1 hr to allow for the dye to be cleared from the gut. The worms were scored for dye filling using a Leica SPX-8 confocal microscope or a Leica fluorescent stereo-microscope. Images were taken using the confocal microscope.
Drug resistance assays
Request a detailed protocolNGM plates containing either TM, paraquat, or ivermectin were used in developmental assays to test for resistance. To this end, we performed a 4 hr egg lay on NGM plates at room temperature, then transferred the eggs to control and drug containing plates, and recorded the number of eggs per plate. After incubation for 4 days at 20°C, plates were scored by counting the total number of L4 larvae or adults on each plate. Paraquat plates were prepared by adding 1 M methyl viologen dichloride (paraquat) directly onto seeded NGM plates for a final concentration of 200 µM and allowed to dry immediately before use. Ivermectin plates were prepared by adding ivermectin directly to the NGM medium at a final concentration of 6 μg/mL before pouring the plates. Modified NGM plates containing no cholesterol were made by using standard NGM without the addition of cholesterol. NGM agar plates containing 5, 25, or 50 µg/mL Cholesterol were used for development when specifically stated. These modified NGM plates were then supplemented with drugs as described above. All of the drug development assays were performed using E. coli OP50 or E. coli HT115 bacteria. For the TM RNAi experiments (Figure 4B), animals were cultured for one generation on pgp-5, pgp-8, pgp-11, pgp-12, or pgp-13 RNAi supplemented with DMSO. The offspring was placed on RNAi plates supplemented with 10 µg/mL TM.
ER stress quantification
Request a detailed protocolSynchronized day one adults were transferred to control and TM containing NGM plates seeded with OP50. After 6 hr of stress induction, hsp-4::GFP levels were measured by large particle flow cytometry in both the WT and osm-3(dh441)IV background using the Union Biometrica Biosorter and imaged using a Leica fluorescent stereo-microscope. For RNAi experiments, animals were raised on control and daf-16 or sel-1 RNAi. Synchronized day one adults were transferred to vehicle or TM containing plates seeded with control or experimental RNAi. After 6 hr of stress, GFP levels were measured using the Biosorter. GFP values were normalized to time of flight. At least 100 animals were analyzed per condition. For RNA isolation, day one adults were treated with 10 mM DTT in S Basal buffer containing E. coli OP50 for 2 hr and then snap frozen.
Quantitative PCR
Request a detailed protocolFor ER stress induction day one adults were washed from their plates and transferred to either control or TM containing plates, where they were incubated for 6 hr. The animals were then washed off using M9 and snap frozen in trizol. Unstressed synchronized animals were collected at day 1 of adult hood for RNA extraction. RNA was prepared using Zymo Research Direct-zol RNA Microprep kit. SYBR green was used to perform quantitative PCR (RT-qPCR). See Supplementary file 2 for list of qPCR primers used in this study. Primer efficiency was tested and was always above 90% and below 110%.
Lifespan analysis
Request a detailed protocolAdult lifespan analysis was performed at 20°C on mutant and WT C. elegans. The animals were synchronized in a 4-hr egg lay. Animals were scored as dead or alive every second day until the end of the experiment. Animals were transferred every day for the first 7 days. Statistical analysis on the Kaplan–Meier survival curves was performed using Microsoft Excel. See Supplementary file 1 for lifespan statistics. Each survival experiment was performed at least twice.
Hydrogen peroxide survival assay
Request a detailed protocolOne molar hydrogen peroxide (H2O2) was added to unseeded NGM plates to a final concentration of 1 µM and allowed to dry for several minutes. Synchronized day one adults were then transferred onto the plates and incubated at 25°C. The animals were then scored every 2 hr for survival. Twenty to 50 animals were used in each experimental condition.
Pathogenic bacteria survival assay
Request a detailed protocolThe P. aeruginosa strain PA14 was grown in LB media and seeded onto high peptone NGM plates and incubated for 12 hr at 37 °C immediately before the start of the experiment. Day one animals were first transferred to unseeded NGM plates to crawl around the plate and remove any excess E. coli OP50 off of their bodies. They were then transferred to the PA14 containing plates and incubated at 25°C. The animals were scored every 12 hr until all animals were dead. At least 75 animals were used in each experimental condition.
Heat stress survival assay
Request a detailed protocolDay one synchronized C. elegans were transferred to fresh NGM plates seeded with E. coli OP50. These plates were transferred to a 35°C incubator, where they were evenly distributed throughout the incubator to ensure even heat exposure. The plates were scored for live worms every 2 hr until all of the worms were dead. Twenty to 50 animals were used in each experimental condition.
Sample collection and RNA purification for sequencing
Request a detailed protocolFor RNA sequencing, we used day one adults that were hatched within 1 hr. To achieve this synchronization, we washed all adults and larvae of plates with M9 and allowed the eggs on hatch for 1 hr. The freshly hatched L1 larvae were then washed off and transferred to fresh plates seeded with OP50 and incubated at 20°C until they developed to adults. At day 1 of adulthood, animals were washed from their plates and transferred to either control or TM containing plates, where they were incubated for 6 hr. The animals were then washed off using M9 and snap frozen in trizol. All biological replicates were collected and prepared using a snaking collection method to reduce batch effects. Total RNA was purified using Zymo Research Direct-zol RNA Microprep kit.
RNA-Seq library preparation
Request a detailed protocolRNA quality was assessed using Agilent’s Bioanalyzer platform, and samples with RIN >9 were used for library construction. 2.5 µg of total RNA was subjected to ribosomal RNA (rRNA) depletion using the Illumina’s Ribo-zero Gold kit (Illumina), according to the manufacturer’s protocol. Strand specific RNA-seq libraries were then constructed using the SMARTer Stranded RNA-Seq HT Kit (Clontech #634839), according to the manufacturer’s protocol. Based on rRNA-depleted input amount, 15 cycles of amplification were performed to generate RNA-seq libraries. Paired-end 150 bp reads were sent for sequencing on the Illumina HiSeq-Xten platform at the Novogene Corporation (USA). The resulting data was then analyzed with a standardized RNA-seq data analysis pipeline (described below).
RNA-seq analysis pipeline cDNA sequences of worm genes were obtained through ENSEMBL Biomart for the WBCel235 build of the C. elegans genome (Ensemble release v94; accessed 2018-12-03). Trimmed reads were mapped to this reference using kallisto 0.43.0–1 and the –fr-stranded option (Bray et al., 2016). All subsequent analyses were performed in the R statistical software (https://cran.r-project.org/). Read counts were imported into R and summarized at the gene level. Differential gene expression was estimated using the ‘DESeq2’ R package (DESeq2 1.16.1). DAVID analysis (version 6.8) was performed to identify significantly enriched GO terms.
Developmental verapamil resistance assay and lifespan analysis
Request a detailed protocolPlates containing 10 µg/mL TM were supplemented with 1 nM verapamil by adding verapamil solved in DMSO onto the plates. Plates were allowed to dry for 6 hr before eggs were transferred to them to determine developmental resistance to TM. Plates were then scored after 4 days to determine the relative number of L4 larvae or adults. Lifespan assays were performed as above; however, the animals were transferred to fresh verapamil containing plates every second day for the whole lifespan assay. Twenty to 50 animals were used in each experimental condition.
Data availability
Raw sequencing data were deposited to the NCBI Gene Expression Omnibus (GEO) under the accession number GSE144675.
-
NCBI Gene Expression OmnibusID GSE144675. NHR-8 regulated P-glycoproteins uncouple xenobiotic stress resistance from longevity in chemosensory C. elegans mutants.
References
-
Inhibition of p-glycoprotein enhances sensitivity of Caenorhabditis elegans to ivermectinVeterinary Parasitology 191:264–275.https://doi.org/10.1016/j.vetpar.2012.09.021
-
Near-optimal probabilistic RNA-seq quantificationNature Biotechnology 34:525–527.https://doi.org/10.1038/nbt.3519
-
Avermectins, new family of potent anthelmintic agents: producing organism and fermentationAntimicrobial Agents and Chemotherapy 15:361–367.https://doi.org/10.1128/AAC.15.3.361
-
Mitochondria are a major source of paraquat-induced reactive oxygen species production in the brainThe Journal of Biological Chemistry 282:14186–14193.https://doi.org/10.1074/jbc.M700827200
-
Hexosamine pathway and (ER) protein quality controlCurrent Opinion in Cell Biology 33:14–18.https://doi.org/10.1016/j.ceb.2014.10.001
-
Mutations in chemosensory cilia cause resistance to paraquat in nematode Caenorhabditis elegansThe Journal of Biological Chemistry 279:20277–20282.https://doi.org/10.1074/jbc.M313119200
-
Broad spectrum detoxification: the major longevity assurance process regulated by insulin/IGF-1 signaling?Mechanisms of Ageing and Development 126:381–387.https://doi.org/10.1016/j.mad.2004.09.001
-
Proteostasis and longevity: when does aging really begin?F1000prime Reports 6:7.https://doi.org/10.12703/P6-7
-
Acquired tolerance to ivermectin and moxidectin after drug selection pressure in the nematode Caenorhabditis elegansAntimicrobial Agents and Chemotherapy 60:4809–4819.https://doi.org/10.1128/AAC.00713-16
-
The sensory amphidial structures of Caenorhabditis elegans are involved in macrocyclic lactone uptake and anthelmintic resistanceInternational Journal for Parasitology 48:1035–1042.https://doi.org/10.1016/j.ijpara.2018.06.003
-
Role of n-oligosaccharide endoplasmic reticulum processing reactions in glycoprotein folding and degradationThe Biochemical Journal 348:1–13.
-
The sense of smell impacts metabolic health and obesityCell Metabolism 26:198–211.https://doi.org/10.1016/j.cmet.2017.06.015
-
P-glycoproteins in nematodesParasitology Today 10:319–322.https://doi.org/10.1016/0169-4758(94)90088-4
-
Tunicamycin, a new antibiotic IThe Journal of Antibiotics 24:215–223.https://doi.org/10.7164/antibiotics.24.215
Article and author information
Author details
Funding
H2020 European Research Council (ERC-StG 640254)
- Martin S Denzel
Max-Planck-Gesellschaft
- Gabriel A Guerrero
German Research Foundation (73111208)
- Martin S Denzel
German Research Foundation (SFB 829)
- Martin S Denzel
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank all Denzel lab members for lively and helpful discussions. We thank Adam Antebi for valuable comments and for C. elegans strains. We thank Seung-Jae Lee for valuable comments on the manuscript. We thank the Caenorhabditis Genetics Center (CGC) for strains. We thank the members of the bioinformatics core facility and the FACS and imaging core facility at MPI AGE. BAB was supported by a seed grant from the NAVIGAGE foundation, a generous gift from the Hanson-Thorell Family and the Kathleen Gilmore award in aging biology. MSD was supported by ERC-StG 640254, by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Projektnummer 73111208 – SFB 829, and by the Max Planck Society.
Copyright
© 2021, Guerrero et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,532
- views
-
- 221
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Genetics and Genomics
- Immunology and Inflammation
Systemic lupus erythematosus (SLE) is an autoimmune disease, the pathophysiology and genetic basis of which are incompletely understood. Using a forward genetic screen in multiplex families with SLE, we identified an association between SLE and compound heterozygous deleterious variants in the non-receptor tyrosine kinases (NRTKs) ACK1 and BRK. Experimental blockade of ACK1 or BRK increased circulating autoantibodies in vivo in mice and exacerbated glomerular IgG deposits in an SLE mouse model. Mechanistically, NRTKs regulate activation, migration, and proliferation of immune cells. We found that the patients’ ACK1 and BRK variants impair efferocytosis, the MERTK-mediated anti-inflammatory response to apoptotic cells, in human induced pluripotent stem cell (hiPSC)-derived macrophages, which may contribute to SLE pathogenesis. Overall, our data suggest that ACK1 and BRK deficiencies are associated with human SLE and impair efferocytosis in macrophages.
-
- Genetics and Genomics
- Microbiology and Infectious Disease
The sustained success of Mycobacterium tuberculosis as a pathogen arises from its ability to persist within macrophages for extended periods and its limited responsiveness to antibiotics. Furthermore, the high incidence of resistance to the few available antituberculosis drugs is a significant concern, especially since the driving forces of the emergence of drug resistance are not clear. Drug-resistant strains of Mycobacterium tuberculosis can emerge through de novo mutations, however, mycobacterial mutation rates are low. To unravel the effects of antibiotic pressure on genome stability, we determined the genetic variability, phenotypic tolerance, DNA repair system activation, and dNTP pool upon treatment with current antibiotics using Mycobacterium smegmatis. Whole-genome sequencing revealed no significant increase in mutation rates after prolonged exposure to first-line antibiotics. However, the phenotypic fluctuation assay indicated rapid adaptation to antibiotics mediated by non-genetic factors. The upregulation of DNA repair genes, measured using qPCR, suggests that genomic integrity may be maintained through the activation of specific DNA repair pathways. Our results, indicating that antibiotic exposure does not result in de novo adaptive mutagenesis under laboratory conditions, do not lend support to the model suggesting antibiotic resistance development through drug pressure-induced microevolution.