{'fraction_non_alphanumeric': 0.03792121721273219, 'fraction_numerical': 0.022472696723606835, 'mean_word_length': 5.025449943757031, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Autophagy and modular restructuring of metabolism control germline tumor differentiation and proliferation in C. elegans

ABSTRACT Autophagy can act either as a tumor suppressor or as a survival mechanism for established tumors. To understand how autophagy plays this dual role in cancer, in vivo models are required. By using a highly heterogeneous C. elegans germline tumor, we show that autophagy-related proteins are expressed in a specific subset of tumor cells, neurons. Inhibition of autophagy impairs neuronal differentiation and increases tumor cell number, resulting in a shorter life span of animals with tumors, while induction of autophagy extends their life span by impairing tumor proliferation. Fasting of animals with fully developed tumors leads to a doubling of their life span, which depends on modular changes in transcription including switches in transcription factor networks and mitochondrial metabolism. Hence, our results suggest that metabolic restructuring, cell-type specific regulation of autophagy and neuronal differentiation constitute central pathways preventing growth of heterogeneous tumors.


Introduction
Germ cell-derived tumors usually result from precocious differentiation of cells in the gonads or the embryo. 1 In C. elegans, induction of germline tumors can be achieved by depletion of gld-1, a member of the STAR KH-domain family of RNA binding proteins and a suppressor of translation. 2 Depletion results in tumors containing fully differentiated neurons, muscle and gut cells (Fig. 1A). 2 Their overproliferation leads to gonad swelling and invasion of other tissues by germline tumor cells, resulting in premature death of animals. 3,4 Mutations that generally extend C. elegans life span inhibit gld-1 tumor growth and prevent tumor-dependent life-span shortening, although the mechanism is not fully understood. 4 Interestingly, metabolic adjustments including upregulation of autophagy have been proposed for these long life-span mutants. [5][6][7] Autophagy is thought to play a dual role in cancer. It can either act as a tumor suppressor mechanism and prevent tumor initiation, or act as an adaptive response to cancer survival-associated stress and sustain cancer metabolism of established tumors. 8 Therefore, it is essential to understand the context-specific functions of autophagy in tumorigenesis. These will depend on tumor type, stage, and genetic alterations, among others. 9 A recent study showed that intratumor heterogeneity in the form of clonal diversity occurs over clinically relevant timeframes and can lead to clinically important phenotypic properties. 10 Hence, the heterogeneity of C. elegans germline tumors suggests that they could constitute an appealing model to study the role of autophagy in tumorigenesis.
Autophagy pathways have been shown to regulate key stages of C. elegans development. 11 During embryogenesis, autophagy selectively eliminates paternal mitochondria [12][13][14] and P granules 15 and contributes to clearance of cell corpses. 16 Autophagy also regulates cell size, 17 is essential for dauer development 7,18 and, together with apoptosis, is required for proper embryonic development. 19 In the former examples, regulation of autophagy has been mostly studied at the post-transcriptional level. Recently, the nuclear hormone receptor NHR-62, 20 the transcription factor HLH-30/TFEB 6,21 and the conserved transcription factor PHA-4/FOXA 5,22 have been identified as transcriptional regulators of autophagy. In addition, pathways such as XBP-1-mediated ER stress and ATFS-1-mediated mitochondrial stress also regulate autophagy at the transcriptional level in C. elegans. 23 However, transcriptional and post-transcriptional regulation of autophagy during differentiation remains to be explored.
We reasoned that the heterogeneity of C. elegans germline tumors, with cells expressing a variety of markers otherwise only found in differentiated somatic cells, calls for studying the role of autophagy in tumorigenesis. Thus, we analyzed expression of autophagy markers and performed a comprehensive transcriptome analysis of animals with germline tumors. Genes related to stress response and proteolysis, such as autophagyrelated genes, are upregulated in animals with tumors. In agreement, autophagy-related proteins are expressed in germline tumor neurons and inhibition of autophagy impairs neuronal differentiation. Conversely, autophagy inhibition leads to an increase in germ cell number and a rapid invasion of head and vulva by the tumor. In addition, depletion of autophagy-related genes shortens the life span of animals with germline tumors, while induction of autophagy extends their life span. In animals LGG-1 fluorescence intensity in gld-1 tumors in DMSO (control) and Baf-treated animals (n D 10 each), ÃÃÃ , P 0.001. (F) Representative maximum projections of GFP::mCherry::LGG-1 expressing animals. A control animal's gonad turn region (left) and a gld-1 RNAi-treated animal's germline tumor turn region (right) are shown. Boxed areas are shown enlarged on the right. Scale bar: 10 mm with fully developed tumors, fasting doubles their life span by triggering both an upregulation of autophagy-related genes and a remodeling of the metabolic transcriptome, including an upregulation of oxidative metabolism. Mitochondrial oxidative stress signaling and potentially also modular activation of gene networks governed by nuclear hormone receptors seem to be part of a positive feedback that leads to upregulation of autophagy. This suggests that even in mature tumors a straight-forward dietary intervention can inhibit tumor growth and that autophagy, metabolic restructuring and mitochondrial stress signaling constitute central antitumorigenic pathways in the C. elegans germline tumor model.

Results
Autophagy-related proteins are expressed in gld-1 germline tumors Since many tumors require drastic changes in metabolism and particularly in catabolism to grow, we asked whether autophagy is induced in germline tumors in C. elegans. To test this, we examined if autophagy-related proteins are expressed in the wild-type gonad and/or in the gld-1 germline tumor. We used transgenic animals that express fluorescence-tagged versions of the following proteins: GFP::LGG-1 (GABARAP ortholog), GFP::LGG-2 (LC3 ortholog), and BEC-1::GFP (BECN1 ortholog). Consistent with previous reports, 7,18,24 we observed that autophagy-related proteins are highly expressed in the nervous system, pharynx, intestine, hypodermis, vulva, spermatheca and somatic gonad of wild-type animals, but not in wild-type germline cells (Fig. 1B). However, cells within the germline tumor in gld-1 RNAi animals show high expression of the autophagy pathway proteins from the first day of adulthood throughout their life span ( Fig. 1C and D). Notably, these transgenes (with the exception of BEC-1::GFP) are not expressed in wild-type germlines and only become activated during zygotic gene activation in embryos (Fig. 1B) or somatic trans-differentiation in the tumor.
To test whether tumor formation leads to a general induction of autophagy rather than to a block of autophagosome maturation, we treated animals with germline tumors with bafilomycin A 1 (Baf) which blocks autophagosome-lysosome fusion. This leads to an increase in GFP::LGG-1 levels in germline tumors, consistent with the idea that autophagy is generally upregulated during tumor formation (Fig. 1E). Additionally, we also observed a massive increase of autophagosomes and lysosomes when using a germline-expressed transgene, Ppie-1:: GFP::mCherry::lgg-1, that allows to monitor autophagosome maturation due to quenching of GFP fluorescence and maintenance of mCherry fluorescence after acidification following autophagosome-lysosome fusion (Fig 1F).
Consistent with germline tumors containing cells expressing different fate markers, 2,4 we also observed heterogenetity in nuclear shape, size and nucleolar structure as well as cell shape ( Fig. 2A). Importantly, germline tumor cells show a heterogenous pattern of GFP::LGG-1 expression, with many cells showing elevated cytoplasmic expression and about 90% of them contain clearly discernible autophagosomes with rather uniform size of around 0.5 mm (Fig. 2B). Thus, unlike what has been reported for larval seam cells where only few autophagosomes can be detected, germline tumor cells show a distribution of autophagosome numbers with 4 or 5 autophagosomes per cell on average (Fig. 2B). Moreover, many cells expressing GFP::LGG-1 have elongated protrusions which are often connected to neighboring cells, thereby resembling neurons in their shape (Fig. 2C). By performing time-lapse microscopy of germline tumors, we also found that while most tumor cells show little intracellular autophagosome movement, in a few tumor cells with high levels of GFP::LGG-1, fast, random-walking autophagosomes can be observed ( Fig. 2D; Movie S1).
In gld-1 animals, changes in translation seem to be accompanied by transcriptional changes. Several gld-1 targets have been identified; [25][26][27] some of them are transcriptional regulators, like pal-1/Caudal, whose translation is repressed by GLD-1. 2 Additionally, GLD-1 can bind and protect its targets from nonsense-mediated mRNA decay (NMD), for instance the glucosamine 6-phosphate N-acetyltransferase gna-2 and the putative transposase Y75B12B.1. 28 To confirm our data using integrated reporters (Figs. 1 and 2) and to gain deeper understanding of the transcriptional regulation in gld-1 tumors, we performed a comprehensive analysis of the transcriptome in wild-type and gld-1 RNAi animals. Genome-wide transcriptome analysis of tumor-containing gld-1 RNAi animals (third d of adulthood) was carried out and compared with wild-type animals of the same age as well as young adults. First, we analyzed the transcripts of the 2 targets protected from NMD. These are strongly downregulated when comparing gld-1 tumor animals with staged gravid adults and young adults (5.2-and 2.6-fold for gna-2, 3.5-and 2.0-fold for Y75B12B.1, respectively). Second, we performed a gene ontology (GO) analysis, which showed that genes related to stress response pathways (e.g., DNA repair, ubiquitin-mediated proteolysis) and autophagy (9 genes) are upregulated in gld-1 RNAi animals compared to young adults ( Fig. 3A and B; Table S1), but not compared to wild-type adults of the same age (Table S2) since expression of these genes in the embryos seems to mask the differences between gld-1-depleted and older wild-type animals. All together, these findings are consistent with a transcriptional upregulation of autophagy-related genes followed by protein synthesis in tumor cells.
Autophagy-related proteins are preferentially expressed in germline tumor neurons gld-1 germline tumors do not only comprise germ cells but also ectopically differentiated somatic cells, such as muscle, neurons and intestinal cells ( Fig. 1A and S1A). 2 We observed that the morphology of cells expressing autophagy markers within the germline tumor resembles the morphology of neuronal cells ( Figs. 1 and 2C). Therefore, we used a C. elegans strain expressing a pan-neuronal marker (UNC-119::mCherry) and GFP:: LGG-1. Confocal microscopy of germline tumors induced by gld-1 RNAi, showed that LGG-1 is expressed in cells differentiating as neurons (Fig. 4A). Concomitantly, we examined differentiation of other somatic cell types in gld-1 germline tumors, by knocking down gld-1 in reporter strains expressing GFP or mCherry under the control of tissue-specific-neuronal (unc-119), muscle (myo-3) or pharyngeal (pha-4)-promoters.
Although gld-1-deficient germ cells reprogram mainly as neurons, we also observed-to a lesser extent-differentiation of muscle and pharynx cells (Fig. S1A). Consistently, transcriptome analysis showed that neuronal-specific genes are consistently upregulated in gld-1-depleted animals (Fig. S1B) and, to a lesser extent, muscle-and pharynx-specific genes are also upregulated (Fig. S1B). Hence, the concurrence of neuronal differentiation and high-level expression of autophagy pathway genes suggests that both processes might depend on each other.

Inhibition of autophagy impairs neuronal differentiation within germline tumors
A growing body of evidence suggests a crucial role for autophagy during neuronal differentiation. [29][30][31][32][33] To investigate whether modulation of autophagy affects neuronal differentiation within germline tumors, we depleted 2 mRNAs simultaneously. 34,35 We used dual depletion for gld-1 and an autophagy-related protein to induce both the formation of a tumor and to inhibit autophagy in the same animal. lgg-1 RNAi and lgg-1Cgld-1 RNAi lead to similar depletion of lgg-1; likewise, lgg-3 RNAi (Atg12 ortholog) and lgg-3Cgld-1 RNAi result in similar accumulation of the unmodified form of LGG-1 and the percentage of animals displaying tumors was similar to single depletion of gld-1 ( Fig. S2C and D). To assess neuronal differentiation, the tumor was induced in the pan-neuronal marker (UNC-119::mCherry) strain by gld-1 RNAi and maximum z-projections covering whole gonad arms were used to score the area occupied by UNC-119::mCherry-expressing cells. Concomitant downregulation of either lgg-3 or atg-7both required for autophagy (Fig. 3C)leads to a reduction of the area occupied by UNC-119::mCherry-expressing cells in the tumor by around 50% and 30%, respectively ( Fig. 4B and C). Therefore, downregulation of autophagy results in an impairment of neuronal differentiation within the germline tumor.

Autophagy affects tumor growth and apoptosis
A balance between self-renewal and differentiation of stem cells is critical for proliferation control and maintenance of a stem cell pool in mammalian tissues. 36 Hence, we tested whether autophagy downregulation impinges on this balance in gld-1 germline tumors. Contrary to the effect of autophagy depletion on neuronal differentiation, simultaneous knockdown of gld-1 and autophagy-related genes leads to a higher number of germ cells, indicative of increased tumor growth (Fig. 5A). In accordance, invasion of the head and vulva region by the tumor occurred earlier in autophagy-deficient animals than in controls (Tables S1 and S2). Moreover, quantitative analysis of gld-1 RNAi-induced tumors shows that nuclei density increases by »40% in the central and proximal part of the tumor in atg-7(bp411) animals ( Fig. 5B and C), and large and dense accumulations of nuclei can be observed under these conditions (Fig. 5B, red arrows; Fig. 5E). Also, a large number of chromosome threads that emerge from apoptotic corpses can be observed. These show a »3-fold increase in atg-7(bp411) animals ( Fig. 5C), particularly in the gonad turn ( Fig. 5D) but also in the proximal part of the tumor ( Fig. 5B and data not shown). However, upregulation of apoptosis cannot compensate for the strong upregulation of proliferation since the tumor apoptosis rate is very low (0.7 § 0.1% apoptotic cells in atg-7(bp411) vs. 0.3 § 0.1% in wild-type animals with gld-1 induced tumors).
Next, we asked whether induction of autophagy would have an opposite effect on tumor growth. Indeed, induction of autophagy by silencing the gene encoding the C. elegans TOR ortholog (let-363) (Fig. 3C) results in a significant reduction of the percentage of animals displaying invasion of the head and vulva region, even at the seventh day of adulthood, both in gld-1 RNAi and mutant animals ( Fig. 5A; Tables S3 and S4). Invasion of the vulva region often results in the appearance of a tumor prolapse that protrudes from the vulva ( Fig. S2C; Movies S2 and S3). Prolapses form by a flow of individual tumor cells and animals with tumor prolapses are motile and often shed their prolapse. In line with modulation of autophagy affecting tumor growth, inhibition of autophagy leads to an increase in the percentage of animals displaying a tumor prolapse both in gld-1 RNAi and mutant animals ( Fig. S2D and E). Conversely, let-363 RNAi results in a reduced percentage of gld-1 mutant animals with a tumor prolapse, consistent with impaired tumor growth (Fig. S2). Moreover, the inhibitory effect of let-363 RNAi is partially mediated through autophagy (Fig. 5E). Taken together, our data demonstrate that autophagy is crucial for curbing tumor proliferation.

Autophagy impinges on gld-1 animals' life span
As autophagy inhibition resulted in faster tumor growth, we asked whether it affects life span of gld-1 mutants. Overproliferation of gld-1 tumors is lethal and manifests in strongly shortened life spans when compared to wild-type animals without tumors ( Fig. 6A; as reported previously refs. 3 and 4). Interfering with autophagy by lgg-3 RNAI further shortens the life span of gld-1 RNAi and mutant animals by »12% but did not affect wild-type life span ( Fig. 6A and B). 5,22 The functionality of the RNAi approach is also obvious from the changes in the ratio between lipidated/unlipidated LGG-1 (Fig. 3C).
Since we observed that induction of autophagy, by let-363 knockdown, extends the life span of gld-1 mutants by »40% (Fig. 6B), we wondered whether inducing autophagy after tumors are already established would affect longevity of gld-1 animals. We triggered autophagy by uninterrupted fasting (Fig. 3C) and found that gld-1 worms fasting from d 2 of adulthood display a striking increase in life span by »100% (Fig. 7A). Remarkably, fasting inhibited tumor growth most likely by enhancing neuronal differentiation within the tumor (Fig. 7B), in line with a function for autophagy on tumor growth or differentiation balance.

Autophagy is required for fasting-induced life-span extension
Next, we analyzed whether autophagy is required for the fasting-induced life-span extension observed in gld-1 deficient animals. To this end, atg-18(gk378) and atg-7(bp411) autophagy mutants were fed with gld-1 RNAi till d 2 of adulthood and then either fasted or kept under feeding conditions throughout their life span. Fasting extends the life span of gld-1 RNAitreated wild-type animals by »100% (Fig. 7C). This is a substantially longer life-span extension than previously reported for fasting of wild-type animals without tumors. 37,38 On the other hand, fasting prolongs the life span of atg-18 (gk378) and atg-7(bp411) mutants by only »15% and »26%, respectively (Fig. 7C), showing that autophagy is a central pathway for fasting-induced life-span extension.

Fasting induces modular metabolic switches
Fasting does not exclusively induce autophagy, but it triggers changes at both transcriptional and post-transcriptional levels. 39 To acquire a broader picture, we performed a transcriptome analysis of gld-1 RNAi animals fasted for 48 h, from d 2 . Autophagosome paths are shown as black lines. Note that the movement of the bottom autophagosome is due to the sample shifting (see Movie S1). All panels in Fig. 2 show animals at the d 3 of adulthood.
of adulthood, and compared their transcriptome with fed gld-1 RNAi animals. Consistent with autophagy being required for fasting-induced life-span extension, we observed that autophagy-related genes, are upregulated in fasted gld-1 RNAi animals (Fig. 7D). In accordance with previous reports 40,41 fasting also induced upregulation of components of SCF-E3 ubiquitin ligase complexs, members of the P450 family of cytochromes, C-type lectins, secreted cysteine-rich, and saposin-like proteins (Table S5).
Functional clustering showed a modular reshuffle of metabolism in fasted, gld-1-depleted animals (Fig. 8A). A consistent upregulation of genes involved in fatty acid b¡oxidation, citric acid cycle, amino acid and propionate catabolism ( Fig. 8B and C; Table S5) was detected in fasted gld-1 RNAi animals. Conversely, fatty acid synthesis genes are downregulated in fasted animals (Fig. 8C). Regarding carbohydrate metabolism, upregulation of genes involved in glycolysis, gluconeogenesis and glycogenesis was detected (Fig. 8B). In line with a metabolic reorganization in fasted gld-1 animals, module 142 'Generation of precursor metabolites and energy, positive regulation of growth, ion transport' described by Vermeirssen and colleagues 42 is coregulated in fasted animals ( Fig. 8D; all components except cst-1 and dox-1, which are also found inversely regulated in other experiments 42 ). This module was identified using transcription profiles from various experiments by applying a reverse-engineering algorithm to extract ensembles of coexpressed genes. Module 142 is related to mitochondrial metabolism and contains transcripts of the mitochondrial phosphate carrier F01G4.6, the complex IV subunit F26E4.6, the ATP synthase subunits F58F12.1 and H28O16.1, the outer membrane voltage-gated anion channel vdac-1, the fructose bisphosphate aldolase aldo-2, and the peroxidase prdx-2. Taken together, the coregulation of several metabolic pathways and of module 142 during fasting leads to an upregulation of the major catabolic pathways and a downregulation of fatty acid synthesis, suggesting that mitochondrial respiration will increase under these conditions (Fig. 8E).
To corroborate whether changes on the transcript level correspond to metabolic alterations, we analyzed fat content of gld-1 RNAi animals by vital staining of lipid droplets with BODIPY 493/503. 43 In agreement with the observed transcriptional upregulation of genes involved in b-oxidation and downregulation of fatty acid synthesis, the number of lipid droplets in fasted animals is lower than in fed animals ( Fig. S3; lipid droplets are also observed in tumor prolapses). This observation sustains the hypothesis that metabolic adjustments reflect the observed transcriptional changes. We additionally examined N-acetyl-Dglucosamine (GlcNAc) levels, as we observed an upregulation of the hexosamine pathway in the tumor and a downregulation when the animals containing tumors were fasted (Fig. S4A). This pathway has been associated with life-span extension. 44 Again in agreement with the transcriptome data, we observe that O-linked GlcNAc-modified proteins are increased when tumors are formed and are not further increased under starvation since the transferase OGA-1 is downregulated (Fig. S4B).
Additionally, we analyzed the expression of genes comprised in a central C. elegans metabolic gene regulatory network. 45 This network was generated from experimental interactions between genes involved in lipid metabolism 46 and fasting responsive genes. 47 A substantial number of genes in this network are nuclear hormone receptors (NHRs), that can function as metabolic sensors and represent an expanded gene family (284 members in C. elegans vs. 48 in H. sapiens). Using topological overlap analysis, 5 modules could be identified in this network. 45 Two of these 5 modules, II and III, are coregulated in fasted gld-1-depleted animals (Fig. 8F), both containing several NHRs. Importantly, all components of module III that are directly linked to lipid catabolism are consistently upregulated during fasting, including nhr-28, -66, -70, -79, -109, -178, and -273, while transcription factors involved in developmental To establish whether these NHRs indeed play roles in tumor metabolism, we focused on NHR-178. For this NHR, it has been shown that depletion by RNAi results in increased lipid levels and that its expression is upregulated in the epidermis on food withdrawal. 45 Consistent with these results and a function of NHR-178 in tumor metabolism, we observed weak expression in germline tumors that increases »3-fold during starvation (Fig. 8G, red). This demonstrates that NHR-178 and most likely also additional NHRs of module III involved in metabolism function in adapting tumor metabolism to environmental conditions.
Although only moderately upregulated during tumorigenesis (1.1-fold) and further upregulated during fasting (1.1-fold), we also analyzed expression of the transcription factor NHR-49 in germline tumors by using an integrated reporter construct that drives histone H1 expression from an nhr-49 promoter. NHR-49 is a major regulator of metabolism and autophagy. 22,47 Notably, we detect expression of the nhr-49 reporter in germline tumors which increases »10-fold during fasting (Fig. 8G,  blue). Therefore, the global changes that we detect by RNAseq seem to underestimate the expression changes inside the tumor. Thus, far-ranging, coordinated changes in transcriptional networks seem to govern the response to fasting, including known (NHR-49, NHR-178, mitochondrial energy metabolism) and previously unknown regulators (Fig. 8F). These findings also support the notion that NHRs seem to act redundantly during metabolic adaptations. 45 Two other regulators, PHA-4 and HLH-30, have been previously proposed as regulators of metabolic changes and autophagy. To gain an even broader picture, we also analyzed their contribution to autophagy-related gene expression in the tumor. Expression of PHA-4 within the tumor is found only in a very low number of tumor cells (Fig. S1) and we could not detect coregulation of PHA-4 targets in the transcriptome data (Table S6). Hence, PHA-4 does not seem to be involved in autophagy regulation in the tumor. Regarding HLH-30, we examined whether HLH-30 depletion in the tumor affects LGG-1 expression or autophagosome formation. We did not observe major changes in GFP::LGG-1 levels in the tumor upon hlh-30 RNAi, although we observed a slight decrease in the number of autophagosomes and an increase in number of large autophagosomes, the latter statistically not being significant ( Fig. S5A to C). To examine the expression of HLH-30 in the tumor, we generated a strain expressing both HLH-30::GFP and UNC-119::mCherry transgenes. Very few cells express HLH-30 compared to the number of cells expressing UNC-119 and those expressing HLH-30 are not tumor neurons (Fig. S5D). This suggests that HLH-30 is not a key transcription factor in autophagy regulation either. In agreement, transcriptome analysis showed that HLH-30 targets are not coregulated in the tumor (Table S7) and lysosomal genes, previously shown to be influenced by HLH-30 levels, 6 are also not coregulated (Table S8).

ROS stress signaling affects tumor autophagy and differentiation
Our analysis reveals that all major pathways leading to mitochondrial catabolism are upregulated in fasted animals, including pathways that lead to anaplerotic reactions for the tricarboxylic acid cycle (Fig. 8E). Activation of these pathways under low glucose supply could lead to activation of AMP-activated protein kinase (aak-2), the main activator of general catabolism, and result in increased reactive oxygen species (ROS) production. Under these conditions, increased ROS levels have been suggested to act as a hormetic stressor, which induces defense mechanisms like phase II detoxification, transduced by a peroxiredoxin (prdx-2), thereby prolonging life span. 48,49 Remarkably, activators, transducers and downstream factors of mitochondrial stress signaling are upregulated in germline tumors animals under starvation, including prdx-2 (being part of module 142, Fig. 8D), aak-2, and many components of phase II detoxification (Fig. 9A). 49,50 To test a function of mitochondrial stress signaling in tumor survival during starvation, we used mutants in this signaling pathway, prdx-2(gk169) and aak-2(ok425). After induction of tumor formation by gld-1 RNAi, we performed life span analysis either under feeding (aak-2 and prdx-2; Figs. 9B and S6A) or starvation conditions (aak-2 only; Figs. 9C and S6B) with animals being additionally exposed either to DOG (2-deoxy-D-glucose, to increase ROS levels),  or an antioxidant (N-acetyL-cysteine, NAC). Notably, life span under feeding in animals with tumors depends on both aak-2 and prdx-2 function (Fig. 9B), whereas starvation-induced life-span extension is not affected in the aak-2 mutant (Fig. 9C), consistent with recent findings. 51 Moreover, glucose restriction by DOG only extends life span in wild-type animals but not in aak-2 mutants (Fig. 9C) and life span under fasting conditions was severely reduced in animals exposed to NAC, demonstrating that ROS production is crucial for starvation-induced life-span extension in animals with tumors (Fig. 9C). To understand whether ROS scavenging also influences tumor cell autophagy and differentiation, we first quantified the extent of neuronal differentiation in the tumor of fasted animals treated with NAC. NAC led to a severe reduction in the tumor area occupied by neurons (Fig. 9D) and exposure to NAC also leads to a severe reduction in the number of autophagosomes (Fig. 9E  and F). Collectively, these findings suggest that a reciprocal dependence exists between induction of autophagy and mitochondrial oxidative metabolism during germline tumorigenesis.

A metabolic model for tumor growth control in C. elegans
We show here that catabolic pathways including autophagy influence cell differentiation and proliferation in germline tumors. Regulation of tumor growth by autophagy, a pathway of central importance when animals rely on internal energy stores, has important implications for the life span of animals with tumors. Autophagy is required for life-span extension by various longevity pathways in C. elegans. 5,7,[52][53][54][55] Here, however, we show that inhibition of autophagy further decreases the life span of severely short-lived, tumor-bearing animals (Fig. 6). Strikingly, induction of autophagy by fasting animals with completely developed tumors from d 2 of adulthood extends their life span by »100% (Fig. 7A and C), fully rescuing the short life-span phenotype of these animals. Fasting from d 2 of adulthood was also reported to extend the life span of wildtype animals, but only by 40 to 50%. 37,38 Furthermore, our genome-wide transcriptome analysis uncovered a reshuffle of the metabolic transcriptome (Fig. 8). Our observations, together with previous reports, suggest a model (Fig. 9G) in which induction of autophagy seems to collaborate with other pathways to affect tumor proliferation through induction of cell-type specific effects: (1) stress-response genes are upregulated in animals with germline tumors (Fig. 3A); (2) long life-span mutations (eat-2, isp-1, clk-1, daf-2), known to induce autophagy, trigger tumor proliferation arrest; [4][5][6][7]22,38 and (3) autophagy promotes neuronal differentiation vs. proliferation within germline tumors while inhibition of autophagy leads to the opposite (Figs. 4C, 7B and 9D to F). Additionally, ROS produced by metabolic activity were also reported to be early activators of fasting-induced autophagy, 56 by oxidizing both AMPK, 57 leading to its activation, and Atg4, resulting in inhibition of its delipidating activity toward the LC3 and GABARAP subfamilies at the autophagosome formation site. 58 Therefore, it is possible that increasing oxidative metabolism and inducing autophagy collaboratively establish cytoprotection, being crucial for positive regulation of life span during fasting.

Strengthening the link between neurogenesis and autophagy
To date, a role for autophagy in cell differentiation during development in C. elegans is restricted to a microRNA-regulated case of bilaterally asymmetric neuron specification. 59 In mammals, accumulating evidence supports a role for autophagy in cell differentiation, particularly neuronal differentiation. 60 We suggest that autophagy might be similarly required for differentiation of neurons or the maintenance of the metabolic needs of this cell type in C. elegans gld-1 tumors since systemic inhibition of autophagy reduces neuronal differentiation ( Fig. 4B and C). In addition, autophagy proteins are strongly reexpressed in transdifferentiating neurons and not much in any other cell type within gld-1 tumors (Figs. 1, 2 and 4). Their regulation appears to be mainly transcriptional (Fig. 1C) with the exception of LGG-1/GABARAP and LGG-2/MAP1LC3 (LC3) whose mRNAs are direct GLD-1 targets 26,27 and their expression is probably regulated post-transcriptionally. Interestingly, in line with our observations, a role for autophagy in neuroendocrine differentiation of human prostate cancer cells has been recently reported, 61 suggesting that a conserved interdependency exists in C. elegans and mammalian tumors.
Moreover, we observed that fasting not only inhibits tumor growth but also induces neuronal differentiation within the tumor (Fig. 7B), in line with observations showing that dietary restriction in mice enhances neurogenesis. 62 Fasting and dietary restriction have shown positive effects in both cancer prevention and treatment in diverse mice models. 39 However, Pinkston and colleagues did not observe changes in the frequency of neuronal differentiation within the tumor in the daf-2 mutant, 4 pointing to the possibility that specific pathways inducing tumor growth arrest, likely also through autophagy, 7 might affect differentiation differently. It is therefore essential to better understand the molecular mechanisms that underlie tumor growth arrest upon activation of different catabolic pathways.

Triggering modular changes in metabolic networks
Our data show that fasted animals with germline tumors upregulate major catabolic pathways and downregulate fatty acid synthesis (Fig. 8), reminiscent of DAF-16-dependent metabolic restructuring. 63 However, DAF-16 targets 64 are not consistently regulated in fasted animals, suggesting that other regulatory mechanisms are responsible for the metabolic remodeling of fasted animals with germline tumors. Nonetheless, gene expression profiles of fasted animals with germline tumors partially overlap with those of daf-2 mutants 64 and fasted wild-type animals, 41 indicating that related pathways underlie longevity extension regardless of the presence of germline tumors. For instance, it has been suggested that the coregulated module 142 'Generation of precursor metabolites and energy, positive regulation of growth, ion transport' may be part of insulin-like signaling coordinating oxidative stress tolerance, metabolic programs and life span since all promoters of module 142 genes contain the DAF-16 consensus binding element. 42 At the organismal level, we observed coregulation of bioinformatically predicted gene modules ( Fig. 8D and F). 42,45 Remarkably, 2 modules that are part of a metabolic gene regulatory network 45 are coregulated in fasted animals (Fig. 8F). These modules are enriched for specific NHRs, the majority of which are consistently upregulated in fasted animals (Fig. 8F). NHRs are transcription factors that integrate intrinsic and environmental cues to coordinate transcriptional cascades that control reproduction, development, metabolism and homeostasis. 65 In particular, NHR-49 and NHR-62 have been shown to mediate the response to nutrient availability 20,66 and NHR-62 is required for dietary restriction-induced autophagy and longevity. 20 Recent studies 52 show that autophagy is also regulated by the NHRs FXR and PPARA (the ortholog of NHR-49) in mammals and that these factors might have antagonistic effects on transcription of autophagy genes. Additionally, NHR-49 has been shown to be involved in the control of fat consumption and the upregulation of mitochondrial b-oxidation 47,67 , which we find is also upregulated in our RNAseq analysis (Fig. 8C). Our findings also implicate NHR-49 in this regulation (Fig. 8G). Furthermore, additional NHRs (e.g., NHR-67, NHR-91 -the latter a likely ortholog of FXR, and the fasting-responsive NHR-178; Fig. 8F and G) seem to be involved in the response to fasting most likely as part of a transcription factor network. Further studies will be required to identify the NHRs directly affecting autophagy gene expression and which additional regulatory roles these NHRs have.
Altogether, our data suggest that spontaneous neuronal differentiation in germline tumors seems to rely on the induction of autophagy, which in turn depends on metabolic restructuring. Differentiation in this context apparently has a positive effect on animal survival. Thus, manipulation of metabolic restructuring seems to constitute a powerful way to drive cell differentiation and to interfere with tumor growth.

RNAi treatment
RNAi was performed by feeding according to standard protocols. 70 Clones used were prepared from cDNA or obtained from commercially available libraries (http://www.lifesciences. sourcebioscience.com/clone-products/mammalian/rnai/mirnarnai-resources.aspx). To knock down 2 target genes simultaneously, a gld-1 cDNA fragment was amplified, digested with XhoI and inserted into a vector (pL4440) containing already a cDNA fragment of the other gene, as described elsewhere. 34,35 Vectors were transformed into the HT115 E. coli strain. Synchronized L1 or L3 (for let-363 clones) larvae were used for RNAi feeding. The following geneÀspecific primers were used to generate RNAi clones: let-363 (forward: CATCGTAGAT-CTTCGGAATTCTTGGAGCAATC; reverse: GATACTAGA-TCTTCATCAACCACTGCAACCAT).

RNA sequencing
Synchronized L1 larvae were placed on L4440 control or gld-1 RNAi plates and kept at 20 C. Animals were either collected as young adults or picked as synchronized adults at the first d of adulthood into fresh RNAi plates. When performing gld-1 RNAi, only animals displaying a tumor at the first d of adulthood were picked. gld-1 RNAi animals were kept in gld-1 RNAi plates until d 3 of adulthood, when harvested. Control animals, kept in L4440 plates, were washed off the plates twice a day with M9 to remove the progeny and harvested at d 3 of adulthood. Fasted gld-1 RNAi animals were moved to unseeded plates at d 2 of adulthood for 48 h, before harvesting. Fed gld-1 RNAi control animals were moved to L4440 seeded plates at d 2 of adulthood for 48 h, before harvesting. A total of approximately 3000 animals pulled from 2 independent experiments were collected per condition. The animals were washed off the plates with M9, washed 4 additional times with M9 to remove remaining bacteria and the pellet of worms was frozen in liquid nitrogen. RNA isolation, cDNA libraries construction and sequencing were conducted by GATC Biotech AG.

RNAseq data analysis
3 0 -adaptor sequences were trimmed from the RNA sequencing reads using custom-built Perl script. Resulting reads were then aligned to Ensembl WS220 assembly using the TopHat. 71 Uniquely mapped reads were quantified by counting the overlap of reads with genes using the HtSeq python package (http:// www-huber.embl.de/users/anders/HTSeq). Differential expression analysis was performed using the DESeq software. 72 Count data were variance-stabilized to overcome the influence of bias due to difference in library size. DAVID 73 was used for gene ontology analysis of genes with 1.2 or higher fold-change between the experimental conditions. Heat maps of log 2 foldchanges in gene expression were produced using the heatmap2 package of R. 74 BODIPY 493/503 staining For vital lipid droplet staining, live animals were incubated in M9/BODIPY 493/503 (6.7 mg/ml; ThermoFisher Scientific, D-3922) solution for 20 min. 43 Animals were then washed 3 times with M9 and used immediately for microscopy.

DAPI staining
As described previously, 3 intact worms were fixed in cold methanol for 5 min at ¡20 C and incubated for 30 min in 100 ng/ mL DAPI (Sigma-Aldrich, D9542) in modified M9 (no Mg 2C ).

Bafilomycin A 1 (Baf) treatment
Animals were either injected with 50 mM Baf (LC Laboratories, B-1080) or DMSO (Carl Roth, A994.1) and imaged 3 h later or incubated in an M9 solution with 25 mM of Baf or DMSO for 6 h.
Animals stained with DAPI were imaged using a Leica DM IL epifluorescence microscope (Leica Microsystems, Wetzlar, Germany) equipped with a Leica HC PL FL 10£/0.3 objective and a Motic Moticam 1SP camera (Motic Deutschland GmbH, Wetzlar, Germany) with Motic acquisition software.

Quantifications
To quantify the area occupied by UNC-119 expressing cells per gonad arm, a qualitative scoring was performed blindly from maximum z-projections covering the whole gonad arms (0% to 25%, 25% to 50%, 50% to 75% or 75% to 100% of area occupied by the tumor scored 0, 1, 2 or 3 respectively) and the average score per gonad arm was calculated. The average score per condition was normalized by the control. Mean and standard error are shown including P values from Student t test statistics or from N-1 2 proportion statistics. The number of animals displaying head or vulva invasion (with clear disruption of the gonadal basement membrane) was counted in DAPI stained animals. Nuclei density and apoptotic nuclei number were counted from maximum projections of the superficial cell layer of the tumor. Autophagosomes were counted using fluorescence intensity cut-offs and defined circular windows in ImageJ.

Life-span analysis
Life-span experiments were carried at 20 C. Synchronized L1 larvae were placed on RNAi plates. Animals were picked as synchronized adults at the first day of adulthood into fresh RNAi plates. When gld-1 RNAi was used, only animals displaying a tumor at the first d of adulthood were picked. Worms were transferred to new plates every other day until the end of the reproductive period and every 4 or 5 d in the case of gld-1 animals or after the end of the reproductive period. Animals were fasted by transferring them to empty feeding plates at d 2 of adulthood; the respective control animals, kept under feeding conditions, were transferred to RNAi feeding plates with bacteria expressing the empty pL4440 vector. Animals were examined every day for touch-provoked movement, until death. Animals that crawled off the plate, displayed extruded internal organs during the first 3 d of adulthood, or died from internally hatched progeny were censored. Kaplan-Meier method and Log-Rank (Mantel-Cox) test were performed using OriginPro 8 or OASIS (Online Application for Survival Analysis). 76 In parallel with life-span analysis, the number of animals displaying tumor-extruded tissue through the vulva from d 3 of adulthood throughout their life span was quantified and the percentage of animals showing tumor prolapses calculated.

Protein extraction and western blot analysis
To prepare total worm extracts, worm pellets were resuspended in lysis buffer (62.5 mM Tris-HCl, pH 6.8, 10% glycerol, 2% SDS (Carl Roth, 0183.3) with glass beads (Carl Roth, N030.1), vortexed for 2 min and heated with shaking at 95 C for 5 min. The samples were centrifuged at 20,000 g for 5 min and supernatant fractions were collected. Protein concentration was determined by BCA assay (ThermoFisher Scientific, 23225).

Disclosure of potential conflicts of interest
No potential conflicts of interest were disclosed.