A Persistence Detector for Metabolic Network Rewiring in an Animal

Worm Perturb-Seq: massively parallel whole-animal RNAi and RNA-seq

Transcriptomes provide highly informative molecular phenotypes that, combined with gene perturbation, can connect genotype to phenotype. An ultimate goal is to perturb every gene and measure transcriptome changes, however, this is challenging, especially in whole animals. Here, we present 'Worm Perturb-Seq (WPS)', a method that provides high-resolution RNA-sequencing profiles for hundreds of replicate perturbations at a time in living animals. WPS introduces multiple experimental advances combining strengths of Caenhorhabditis elegans genetics and multiplexed RNA-sequencing with a novel analytical framework, EmpirDE. EmpirDE leverages the unique power of large transcriptomic datasets and improves statistical rigor by using gene-specific empirical null distributions to identify DEGs. We apply WPS to 103 nuclear hormone receptors (NHRs) and find a striking 'pairwise modularity' in which pairs of NHRs regulate shared target genes. We envision the advances of WPS to be useful not only for C. elegans, but broadly for other models, including human cells.

Systems-level design principles of metabolic rewiring in an animal

The regulation of metabolism is vital to any organism and can be achieved by transcriptionally activating or repressing metabolic genes1-3. Although many examples of transcriptional metabolic rewiring have been reported4, a systems-level study of how metabolism is rewired in response to metabolic perturbations is lacking in any animal. Here we apply Worm Perturb-Seq (WPS)-a high-throughput method combining whole-animal RNA-interference and RNA-sequencing5-to around 900 metabolic genes in the nematode Caenorhabditis elegans. We derive a metabolic gene regulatory network (mGRN) in which 385 perturbations are connected to 9,414 genes by more than 110,000 interactions. The mGRN has a highly modular structure in which 22 perturbation clusters connect to 44 gene expression programs. The mGRN reveals different modes of transcriptional rewiring from simple reaction and pathway compensation to rerouting and more complex network coordination. Using metabolic network modelling, we identify a design principle of transcriptional rewiring that we name the compensation-repression (CR) model. The CR model explains most transcriptional responses in metabolic genes and reveals a high level of compensation and repression in five core metabolic functions related to energy and biomass. We provide preliminary evidence that the CR model may also explain transcriptional metabolic rewiring in human cells.

A systems-level, semi-quantitative landscape of metabolic flux in C. elegans

Metabolic flux, or the rate of metabolic reactions, is one of the most fundamental metrics describing the status of metabolism in living organisms. However, measuring fluxes across the entire metabolic network remains nearly impossible, especially in multicellular organisms. Computational methods based on flux balance analysis have been used with genome-scale metabolic network models to predict network-level flux wiring1-6. However, such approaches have limited power because of the lack of experimental constraints. Here, we introduce a strategy that infers whole-animal metabolic flux wiring from transcriptional phenotypes in the nematode Caenorhabditis elegans. Using a large-scale Worm Perturb-Seq (WPS) dataset for roughly 900 metabolic genes7, we show that the transcriptional response to metabolic gene perturbations can be integrated with the metabolic network model to infer a highly constrained, semi-quantitative flux distribution. We discover several features of adult C. elegans metabolism, including cyclic flux through the pentose phosphate pathway, lack of de novo purine synthesis flux and the primary use of amino acids and bacterial RNA as a tricarboxylic acid cycle carbon source, all of which we validate by stable isotope tracing. Our strategy for inferring metabolic wiring based on transcriptional phenotypes should be applicable to a variety of systems, including human cells.

Enhanced flux potential analysis links changes in enzyme expression to metabolic flux

Algorithms that constrain metabolic network models with enzyme levels to predict metabolic activity assume that changes in enzyme levels are indicative of flux variations. However, metabolic flux can also be regulated by other mechanisms such as allostery and mass action. To systematically explore the relationship between fluctuations in enzyme expression and flux, we combine available yeast proteomic and fluxomic data to reveal that flux changes can be best predicted from changes in enzyme levels of pathways, rather than the whole network or only cognate reactions. We implement this principle in an 'enhanced flux potential analysis' (eFPA) algorithm that integrates enzyme expression data with metabolic network architecture to predict relative flux levels of reactions including those regulated by other mechanisms. Applied to human data, eFPA consistently predicts tissue metabolic function using either proteomic or transcriptomic data. Additionally, eFPA efficiently handles data sparsity and noisiness, generating robust flux predictions with single-cell gene expression data. Our approach outperforms alternatives by striking an optimal balance, evaluating enzyme expression at pathway level, rather than either single-reaction or whole-network levels.

Functional Specialization of S-Adenosylmethionine Synthases Links Phosphatidylcholine to Mitochondrial Function and Stress Survival

S-adenosylmethionine (SAM), produced by SAM synthases, is critical for various cellular regulatory pathways and the synthesis of diverse metabolites. Studies have often equated the effects of knocking down one synthase with broader SAM-dependent outcomes such as histone methylation or phosphatidylcholine (PC) production. Humans and many other organisms express multiple SAM synthases. Evidence in Caenorhabditis elegans , which possesses four SAM synthase genes, suggest that the enzymatic source of SAM impacts its function. For instance, loss of sams-1 leads to enhanced heat shock survival and increased lifespan, whereas reducing sams-4 adversely affects heat stress survival. Here, we show that SAMS-1 contributes to a variety of intermediary metabolic pathways, whereas SAMS-4 is more important to generate SAM for methylation reactions. We demonstrate that loss of sams-1 exerts age-dependent effects on nuclear-encoded mitochondrial gene expression, mitochondrial metabolites, and may induce mitophagy. We propose a mechanistic model where reduced SAM from SAMS-1 acts through PC to impact mitochondria, thereby enhancing survival during heat stress.

Worm Perturb-Seq: massively parallel whole-animal RNAi and RNA-seq

The transcriptome provides a highly informative molecular phenotype to connect genotype to phenotype and is most frequently measured by RNA-sequencing (RNA-seq). Therefore, an ultimate goal is to perturb every gene and measure changes in the transcriptome. However, this remains challenging, especially in intact organisms due to different experimental and computational challenges. Here, we present 'Worm Perturb-Seq (WPS)', which provides high-resolution RNA-seq profiles for hundreds of replicate perturbations at a time in a living animal. WPS introduces multiple experimental advances that combine strengths of bulk and single cell RNA-seq, and that further provides an analytical framework, EmpirDE, that leverages the unique power of the large WPS datasets. EmpirDE identifies differentially expressed genes (DEGs) by using gene-specific empirical null distributions, rather than control conditions alone, thereby systematically removing technical biases and improving statistical rigor. We applied WPS to 103 Caenhorhabditis elegans nuclear hormone receptors (NHRs) to delineate a Gene Regulatory Network (GRN) and found that this GRN presents a striking 'pairwise modularity' where pairs of NHRs regulate shared target genes. We envision that the experimental and analytical advances of WPS should be useful not only for C. elegans, but will be broadly applicable to other models, including human cells.

idh-1 neomorphic mutation confers sensitivity to vitamin B12 in Caenorhabditis elegans

In humans, a neomorphic isocitrate dehydrogenase mutation (idh-1neo) causes increased levels of cellular D-2-hydroxyglutarate (D-2HG), a proposed oncometabolite. However, the physiological effects of increased D-2HG and whether additional metabolic changes occur in the presence of an idh-1neo mutation are not well understood. We created a Caenorhabditis elegans model to study the effects of the idh-1neo mutation in a whole animal. Comparing the phenotypes exhibited by the idh-1neo to ∆dhgd-1 (D-2HG dehydrogenase) mutant animals, which also accumulate D-2HG, we identified a specific vitamin B12 diet-dependent vulnerability in idh-1neo mutant animals that leads to increased embryonic lethality. Through a genetic screen, we found that impairment of the glycine cleavage system, which generates one-carbon donor units, exacerbates this phenotype. In addition, supplementation with alternate sources of one-carbon donors suppresses the lethal phenotype. Our results indicate that the idh-1neo mutation imposes a heightened dependency on the one-carbon pool and provides a further understanding of how this oncogenic mutation rewires cellular metabolism.

Nutritional vitamin B12 regulates RAS/MAPK-mediated cell fate decisions through one-carbon metabolism

Vitamin B12 is an essential nutritional co-factor for the folate and methionine cycles, which together constitute one-carbon metabolism. Here, we show that dietary uptake of vitamin B12 modulates cell fate decisions controlled by the conserved RAS/MAPK signaling pathway in C. elegans. A bacterial diet rich in vitamin B12 increases vulval induction, germ cell apoptosis and oocyte differentiation. These effects are mediated by different one-carbon metabolites in a tissue-specific manner. Vitamin B12 enhances via the choline/phosphatidylcholine metabolism vulval induction by down-regulating fat biosynthesis genes and increasing H3K4 tri-methylation, which results in increased expression of RAS/MAPK target genes. Furthermore, the nucleoside metabolism and H3K4 tri-methylation positively regulate germ cell apoptosis and oocyte production. Using mammalian cells carrying different activated KRAS and BRAF alleles, we show that the effects of methionine on RAS/MAPK-regulated phenotype are conserved in mammals. Our findings suggest that the vitamin B12-dependent one-carbon metabolism is a limiting factor for diverse RAS/MAPK-induced cellular responses.

Differential utilization of vitamin B12-dependent and independent pathways for propionate metabolism across human cells

Propionic acid links the oxidation of branched-chain amino acids and odd-chain fatty acids to the TCA cycle. Gut microbes ferment complex fiber remnants, generating high concentrations of short chain fatty acids, acetate, propionate and butyrate, which are shared with the host as fuel sources. Analysis of vitamin B12-dependent propionate utilization in skin biopsy samples has been used to characterize and diagnose underlying inborn errors of cobalamin (or B12) metabolism. In these cells, the B12-dependent enzyme, methylmalonyl-CoA mutase (MMUT), plays a central role in funneling propionate to the TCA cycle intermediate, succinate. Our understanding of the fate of propionate in other cell types, specifically, the involvement of the β-oxidation-like and methylcitrate pathways, is limited. In this study, we have used [14C]-propionate tracing in combination with genetic ablation or inhibition of MMUT, to reveal the differential utilization of the B12-dependent and independent pathways for propionate metabolism in fibroblast versus colon cell lines. We demonstrate that itaconate can be used as a tool to investigate MMUT-dependent propionate metabolism in cultured cell lines. While MMUT gates the entry of propionate carbons into the TCA cycle in fibroblasts, colon-derived cell lines exhibit a quantitatively significant or exclusive reliance on the β-oxidation-like pathway. Lipidomics and metabolomics analyses reveal that propionate elicits pleiotropic changes, including an increase in odd-chain glycerophospholipids, and perturbations in the purine nucleotide cycle and arginine/nitric oxide metabolism. The metabolic rationale and the regulatory mechanisms underlying the differential reliance on propionate utilization pathways at a cellular, and possibly tissue level, warrant further elucidation.

The Caenorhabditis elegans proteome response to two protective Pseudomonas symbionts

The Caenorhabditis elegans natural microbiota isolates Pseudomonas lurida MYb11 and Pseudomonas fluorescens MYb115 protect the host against pathogens through distinct mechanisms. While P. lurida produces an antimicrobial compound and directly inhibits pathogen growth, P. fluorescens MYb115 protects the host without affecting pathogen growth. It is unknown how these two protective microbes affect host biological processes. We used a proteomics approach to elucidate the C. elegans response to MYb11 and MYb115. We found that both Pseudomonas isolates increase vitellogenin protein production in young adults, which confirms previous findings on the effect of microbiota on C. elegans reproductive timing. Moreover, the C. elegans responses to MYb11 and MYb115 exhibit common signatures with the response to other vitamin B12-producing bacteria, emphasizing the importance of vitamin B12 in C. elegans-microbe metabolic interactions. We further analyzed signatures in the C. elegans response specific to MYb11 or MYb115. We provide evidence for distinct modifications in lipid metabolism by both symbiotic microbes. We could identify the activation of host-pathogen defense responses as an MYb11-specific proteome signature and provide evidence that the intermediate filament protein IFB-2 is required for MYb115-mediated protection. These results indicate that MYb11 not only produces an antimicrobial compound but also activates host antimicrobial defenses, which together might increase resistance to infection. In contrast, MYb115 affects host processes such as lipid metabolism and cytoskeleton dynamics, which might increase host tolerance to infection. Overall, this study pinpoints proteins of interest that form the basis for additional exploration into the mechanisms underlying C. elegans microbiota-mediated protection from pathogen infection and other microbiota-mediated traits.IMPORTANCESymbiotic bacteria can defend their host against pathogen infection. While some protective symbionts directly interact with pathogenic bacteria, other protective symbionts elicit a response in the host that improves its own pathogen defenses. To better understand how a host responds to protective symbionts, we examined which host proteins are affected by two protective Pseudomonas bacteria in the model nematode Caenorhabditis elegans. We found that the C. elegans response to its protective symbionts is manifold, which was reflected in changes in proteins that are involved in metabolism, the immune system, and cell structure. This study provides a foundation for exploring the contribution of the host response to symbiont-mediated protection from pathogen infection.

Metabolic rescue of α-synuclein-induced neurodegeneration through propionate supplementation and intestine-neuron signaling in C. elegans

Microbial metabolites that can modulate neurodegeneration are promising therapeutic targets. Here, we found that the short-chain fatty acid propionate protects against α-synuclein-induced neuronal death and locomotion defects in a Caenorhabditis elegans model of Parkinson's disease (PD) through bidirectional regulation between the intestine and neurons. Both depletion of dietary vitamin B12, which induces propionate breakdown, and propionate supplementation suppress neurodegeneration and reverse PD-associated transcriptomic aberrations. Neuronal α-synuclein aggregation induces intestinal mitochondrial unfolded protein response (mitoUPR), which leads to reduced propionate levels that trigger transcriptional reprogramming in the intestine and cause defects in energy production. Weakened intestinal metabolism exacerbates neurodegeneration through interorgan signaling. Genetically enhancing propionate production or overexpressing metabolic regulators downstream of propionate in the intestine rescues neurodegeneration, which then relieves mitoUPR. Importantly, propionate supplementation suppresses neurodegeneration without reducing α-synuclein aggregation, demonstrating metabolic rescue of neuronal proteotoxicity downstream of protein aggregates. Our study highlights the involvement of small metabolites in the gut-brain interaction in neurodegenerative diseases.

idh-1 neomorphic mutation confers sensitivity to vitamin B12 via increased dependency on one-carbon metabolism in Caenorhabditis elegans

The isocitrate dehydrogenase neomorphic mutation ( idh-1neo ) generates increased levels of cellular D-2-hydroxyglutarate (D-2HG), a proposed oncometabolite. However, the physiological effects of increased D-2HG and whether additional metabolic changes occur in the presence of an idh-1neo mutation are not well understood. We created a C. elegans model to study the effects of the idh-1neo mutation in a whole animal. Comparing the phenotypes exhibited by the idh-1neo to Δdhgd-1 (D-2HG dehydrogenase) mutant animals, which also accumulate D-2HG, we identified a specific vitamin B12 diet-dependent vulnerability in idh-1neo mutant animals that leads to increased embryonic lethality. Through a genetic screen we found that impairment of the glycine cleavage system, which generates one-carbon donor units, exacerbates this phenotype. Additionally, supplementation with an alternate source of one-carbon donors suppresses the lethal phenotype. Our results indicate that the idh-1neo mutation imposes a heightened dependency on the one-carbon pool and provides a further understanding how this oncogenic mutation rewires cellular metabolism.

Starvation resistance in the nematode Pristionchus pacificus requires a conserved supplementary nuclear receptor

Nuclear hormone receptors (NHRs) are a deeply-conserved superfamily of metazoan transcription factors, which fine-tune the expression of their regulatory target genes in response to a plethora of sensory inputs. In nematodes, NHRs underwent an explosive expansion and many species have hundreds of nhr genes, most of which remain functionally uncharacterized. However, recent studies have reported that two sister receptors, Ppa-NHR-1 and Ppa-NHR-40, are crucial regulators of feeding-structure morphogenesis in the diplogastrid model nematode Pristionchus pacificus. In the present study, we functionally characterize Ppa-NHR-10, the sister paralog of Ppa-NHR-1 and Ppa-NHR-40, aiming to reveal whether it too regulates aspects of feeding-structure development. We used CRISPR/CAS9-mediated mutagenesis to create small frameshift mutations of this nuclear receptor gene and applied a combination of geometric morphometrics and unsupervised clustering to characterize potential mutant phenotypes. However, we found that Ppa-nhr-10 mutants do not show aberrant feeding-structure morphologies. Instead, multiple RNA-seq experiments revealed that many of the target genes of this receptor are involved in lipid catabolic processes. We hypothesized that their mis-regulation could affect the survival of mutant worms during starvation, where lipid catabolism is often essential. Indeed, using novel survival assays, we found that mutant worms show drastically decreased starvation resistance, both as young adults and as dauer larvae. We also characterized genome-wide changes to the transcriptional landscape in P. pacificus when exposed to 24 h of acute starvation, and found that Ppa-NHR-10 partially regulates some of these responses. Taken together, these results demonstrate that Ppa-NHR-10 is broadly required for starvation resistance and regulates different biological processes than its closest paralogs Ppa-NHR-1 and Ppa-NHR-40.

Abalone peptide increases stress resilience and cost-free longevity via SKN-1-governed transcriptional metabolic reprogramming in C. elegans

A major goal of healthy aging is to prevent declining resilience and increasing frailty, which are associated with many chronic diseases and deterioration of stress response. Here, we propose a loss-or-gain survival model, represented by the ratio of cumulative stress span to life span, to quantify stress resilience at organismal level. As a proof of concept, this is demonstrated by reduced survival resilience in Caenorhabditis elegans exposed to exogenous oxidative stress induced by paraquat or with endogenous proteotoxic stress caused by polyglutamine or amyloid-β aggregation. Based on this, we reveal that a hidden peptide ("cryptide")-AbaPep#07 (SETYELRK)-derived from abalone hemocyanin not only enhances survival resilience against paraquat-induced oxidative stress but also rescues proteotoxicity-mediated behavioral deficits in C. elegans, indicating its capacity against stress and neurodegeneration. Interestingly, AbaPep#07 is also found to increase cost-free longevity and age-related physical fitness in nematodes. We then demonstrate that AbaPep#07 can promote nuclear localization of SKN-1/Nrf, but not DAF-16/FOXO, transcription factor. In contrast to its effects in wild-type nematodes, AbaPep#07 cannot increase oxidative stress survival and physical motility in loss-of-function skn-1 mutant, suggesting an SKN-1/Nrf-dependent fashion of these effects. Further investigation reveals that AbaPep#07 can induce transcriptional activation of immune defense, lipid metabolism, and metabolic detoxification pathways, including many SKN-1/Nrf target genes. Together, our findings demonstrate that AbaPep#07 is able to boost stress resilience and reduce behavioral frailty via SKN-1/Nrf-governed transcriptional reprogramming, and provide an insight into the health-promoting potential of antioxidant cryptides as geroprotectors in aging and associated conditions.

Transcriptional suppression of sphingolipid catabolism controls pathogen resistance in C. elegans

Sphingolipids are required for diverse biological functions and are degraded by specific catabolic enzymes. However, the mechanisms that regulate sphingolipid catabolism are not known. Here we characterize a transcriptional axis that regulates sphingolipid breakdown to control resistance against bacterial infection. From an RNAi screen for transcriptional regulators of pathogen resistance in the nematode C. elegans, we identified the nuclear hormone receptor nhr-66, a ligand-gated transcription factor homologous to human hepatocyte nuclear factor 4. Tandem chromatin immunoprecipitation-sequencing and RNA sequencing experiments revealed that NHR-66 is a transcriptional repressor, which directly targets sphingolipid catabolism genes. Transcriptional de-repression of two sphingolipid catabolic enzymes in nhr-66 loss-of-function mutants drives the breakdown of sphingolipids, which enhances host susceptibility to infection with the bacterial pathogen Pseudomonas aeruginosa. These data define transcriptional control of sphingolipid catabolism in the regulation of cellular sphingolipids, a process that is necessary for pathogen resistance.

Nuclear hormone receptor NHR-49 is an essential regulator of stress resilience and healthy aging in Caenorhabditis elegans

The genome of Caenorhabditis elegans encodes 284 nuclear hormone receptor, which perform diverse functions in development and physiology. One of the best characterized of these is NHR-49, related in sequence and function to mammalian hepatocyte nuclear factor 4α and peroxisome proliferator-activated receptor α. Initially identified as regulator of lipid metabolism, including fatty acid catabolism and desaturation, additional important roles for NHR-49 have since emerged. It is an essential contributor to longevity in several genetic and environmental contexts, and also plays vital roles in the resistance to several stresses and innate immune response to infection with various bacterial pathogens. Here, we review how NHR-49 is integrated into pertinent signaling circuits and how it achieves its diverse functions. We also highlight areas for future investigation including identification of regulatory inputs that drive NHR-49 activity and identification of tissue-specific gene regulatory outputs. We anticipate that future work on this protein will provide information that could be useful for developing strategies to age-associated declines in health and age-related human diseases.

Optogenetic manipulation identifies the roles of ERK and AKT dynamics in controlling mouse embryonic stem cell exit from pluripotency

ERK and AKT signaling control pluripotent cell self-renewal versus differentiation. ERK pathway activity over time (i.e., dynamics) is heterogeneous between individual pluripotent cells, even in response to the same stimuli. To analyze potential functions of ERK and AKT dynamics in controlling mouse embryonic stem cell (ESC) fates, we developed ESC lines and experimental pipelines for the simultaneous long-term manipulation and quantification of ERK or AKT dynamics and cell fates. We show that ERK activity duration or amplitude or the type of ERK dynamics (e.g., transient, sustained, or oscillatory) alone does not influence exit from pluripotency, but the sum of activity over time does. Interestingly, cells retain memory of previous ERK pulses, with duration of memory retention dependent on duration of previous pulse length. FGF receptor/AKT dynamics counteract ERK-induced pluripotency exit. These findings improve our understanding of how cells integrate dynamics from multiple signaling pathways and translate them into cell fate cues.

A D-2-hydroxyglutarate dehydrogenase mutant reveals a critical role for ketone body metabolism in Caenorhabditis elegans development

In humans, mutations in D-2-hydroxyglutarate (D-2HG) dehydrogenase (D2HGDH) result in D-2HG accumulation, delayed development, seizures, and ataxia. While the mechanisms of 2HG-associated diseases have been studied extensively, the endogenous metabolism of D-2HG remains unclear in any organism. Here, we find that, in Caenorhabditis elegans, D-2HG is produced in the propionate shunt, which is transcriptionally activated when flux through the canonical, vitamin B12-dependent propionate breakdown pathway is perturbed. Loss of the D2HGDH ortholog, dhgd-1, results in embryonic lethality, mitochondrial defects, and the up-regulation of ketone body metabolism genes. Viability can be rescued by RNAi of hphd-1, which encodes the enzyme that produces D-2HG or by supplementing either vitamin B12 or the ketone bodies 3-hydroxybutyrate (3HB) and acetoacetate (AA). Altogether, our findings support a model in which C. elegans relies on ketone bodies for energy when vitamin B12 levels are low and in which a loss of dhgd-1 causes lethality by limiting ketone body production.

Systems-level transcriptional regulation of Caenorhabditis elegans metabolism

Metabolism is controlled to ensure organismal development and homeostasis. Several mechanisms regulate metabolism, including allosteric control and transcriptional regulation of metabolic enzymes and transporters. So far, metabolism regulation has mostly been described for individual genes and pathways, and the extent of transcriptional regulation of the entire metabolic network remains largely unknown. Here, we find that three-quarters of all metabolic genes are transcriptionally regulated in the nematode Caenorhabditis elegans. We find that many annotated metabolic pathways are coexpressed, and we use gene expression data and the iCEL1314 metabolic network model to define coregulated subpathways in an unbiased manner. Using a large gene expression compendium, we determine the conditions where subpathways exhibit strong coexpression. Finally, we develop "WormClust," a web application that enables a gene-by-gene query of genes to view their association with metabolic (sub)-pathways. Overall, this study sheds light on the ubiquity of transcriptional regulation of metabolism and provides a blueprint for similar studies in other organisms, including humans.

S-adenosylmethionine synthases specify distinct H3K4me3 populations and gene expression patterns during heat stress

Methylation is a widely occurring modification that requires the methyl donor S-adenosylmethionine (SAM) and acts in regulation of gene expression and other processes. SAM is synthesized from methionine, which is imported or generated through the 1-carbon cycle (1 CC). Alterations in 1 CC function have clear effects on lifespan and stress responses, but the wide distribution of this modification has made identification of specific mechanistic links difficult. Exploiting a dynamic stress-induced transcription model, we find that two SAM synthases in Caenorhabditis elegans, SAMS-1 and SAMS-4, contribute differently to modification of H3K4me3, gene expression and survival. We find that sams-4 enhances H3K4me3 in heat shocked animals lacking sams-1, however, sams-1 cannot compensate for sams-4, which is required to survive heat stress. This suggests that the regulatory functions of SAM depend on its enzymatic source and that provisioning of SAM may be an important regulatory step linking 1 CC function to phenotypes in aging and stress.

Mitochondrial form and function in hair cells

Hair cells (HCs) are specialised sensory receptors residing in the neurosensory epithelia of inner ear sense organs. The precise morphological and physiological properties of HCs allow us to perceive sound and interact with the world around us. Mitochondria play a significant role in normal HC function and are also intricately involved in HC death. They generate ATP essential for sustaining the activity of ion pumps, Ca²⁺ transporters and the integrity of the stereociliary bundle during transduction as well as regulating cytosolic calcium homoeostasis during synaptic transmission. Advances in imaging techniques have allowed us to study mitochondrial populations throughout the HC, and how they interact with other organelles. These analyses have identified distinct mitochondrial populations between the apical and basolateral portions of the HC, in which mitochondrial morphology appears determined by the physiological processes in the different cellular compartments. Studies in HCs across species show that ototoxic agents, ageing and noise damage directly impact mitochondrial structure and function resulting in HC death. Deciphering the molecular mechanisms underlying this mitochondrial sensitivity, and how their morphology relates to their function during HC death, requires that we first understand this relationship in the context of normal HC function.

A Transcriptional Cofactor Regulatory Network for the C. elegans Intestine

Chromatin modifiers and transcriptional cofactors (collectively referred to as CFs) work with DNA-binding transcription factors (TFs) to regulate gene expression. In multicellular eukaryotes, distinct tissues each execute their own gene expression program for accurate differentiation and subsequent functionality. While the function of TFs in differential gene expression has been studied in detail in many systems, the contribution of CFs has remained less explored. Here we uncovered the contributions of CFs to gene regulation in the Caenorhabditis elegans intestine. We first annotated 366 CFs encoded by the C. elegans genome and assembled a library of 335 RNAi clones. Using this library, we analyzed the effects of individually depleting these CFs on the expression of 19 fluorescent transcriptional reporters in the intestine and identified 216 regulatory interactions. We found that different CFs interact specifically with different promoters, and that both essential and intestinally expressed CFs exhibit the highest proportion of interactions. We did not find all members of CF complexes acting on the same set of reporters but instead found diversity in the promoter targets of each complex component. Finally, we found that previously identified activation mechanisms for the acdh-1 promoter use different CFs and TFs. Overall, we demonstrate that CFs function specifically rather than ubiquitously at intestinal promoters and provide an RNAi resource for reverse genetic screens.

Fatty acids derived from the probiotic Lacticaseibacillus rhamnosus HA-114 suppress age-dependent neurodegeneration

The human microbiota is believed to influence health. Microbiome dysbiosis may be linked to neurological conditions like Alzheimer's disease, amyotrophic lateral sclerosis, and Huntington's disease. We report the ability of a probiotic bacterial strain in halting neurodegeneration phenotypes. We show that Lacticaseibacillus rhamnosus HA-114 is neuroprotective in C. elegans models of amyotrophic lateral sclerosis and Huntington's disease. Our results show that neuroprotection from L. rhamnosus HA-114 is unique from other L. rhamnosus strains and resides in its fatty acid content. Neuroprotection by L. rhamnosus HA-114 requires acdh-1/ACADSB, kat-1/ACAT1 and elo-6/ELOVL3/6, which are associated with fatty acid metabolism and mitochondrial β-oxidation. Our data suggest that disrupted lipid metabolism contributes to neurodegeneration and that dietary intervention with L. rhamnosus HA-114 restores lipid homeostasis and energy balance through mitochondrial β-oxidation. Our findings encourage the exploration of L. rhamnosus HA-114 derived interventions to modify the progression of neurodegenerative diseases.

A metabolic regulatory network for the Caenorhabditis elegans intestine

Metabolic perturbations can affect gene expression, for instance to rewire metabolism. While numerous efforts have measured gene expression in response to individual metabolic perturbations, methods that determine all metabolic perturbations that affect the expression for a given gene or set of genes have not been available. Here, we use a gene-centered approach to derive a first-pass metabolic regulatory network for Caenorhabditis elegans by performing RNAi of more than 1,400 metabolic genes with a set of 19 promoter reporter strains that express a fluorescent protein in the animal's intestine. We find that metabolic perturbations generally increase promoter activity, which contrasts with transcription factor (TF) RNAi, which tends to repress promoter activity. We identify several TFs that modulate promoter activity in response to perturbations of the electron transport chain and explore complex genetic interactions among metabolic pathways. This work provides a blueprint for a systems-level understanding of how metabolism affects gene expression.

Dietary vitamin B12 regulates chemosensory receptor gene expression via the MEF2 transcription factor in Caenorhabditis elegans

Dynamic changes in chemoreceptor gene expression levels in sensory neurons are one strategy that an animal can use to modify their responses to dietary changes. However, the mechanisms underlying diet-dependent modulation of chemosensory gene expression are unclear. Here, we show that the expression of the srh-234 chemoreceptor gene localized in a single ADL sensory neuron type of Caenorhabditis elegans is downregulated when animals are fed a Comamonas aquatica bacterial diet, but not on an Escherichia coli diet. Remarkably, this diet-modulated effect on srh-234 expression is dependent on the micronutrient vitamin B12 endogenously produced by Comamonas aq. bacteria. Excess propionate and genetic perturbations in the canonical and shunt propionate breakdown pathways are able to override the repressive effects of vitamin B12 on srh-234 expression. The vitamin B12-mediated regulation of srh-234 expression levels in ADL requires the MEF-2 MADS domain transcription factor, providing a potential mechanism by which dietary vitamin B12 may transcriptionally tune individual chemoreceptor genes in a single sensory neuron type, which in turn may change animal responses to biologically relevant chemicals in their diet.

Integrated -omics approach reveals persistent DNA damage rewires lipid metabolism and histone hyperacetylation via MYS-1/Tip60

Although DNA damage is intricately linked to metabolism, the metabolic alterations that occur in response to DNA damage are not well understood. We use a DNA repair-deficient model of ERCC1-XPF in Caenorhabditis elegans to gain insights on how genotoxic stress drives aging. Using multi-omic approach, we discover that nuclear DNA damage promotes mitochondrial β-oxidation and drives a global loss of fat depots. This metabolic shift to β-oxidation generates acetyl-coenzyme A to promote histone hyperacetylation and an associated change in expression of immune-effector and cytochrome genes. We identify the histone acetyltransferase MYS-1, as a critical regulator of this metabolic-epigenetic axis. We show that in response to DNA damage, polyunsaturated fatty acids, especially arachidonic acid (AA) and AA-related lipid mediators, are elevated and this is dependent on mys-1. Together, these findings reveal that DNA damage alters the metabolic-epigenetic axis to drive an immune-like response that can promote age-associated decline.

Nematicidal effects of 2-methyl-aconitate isomerase from the phytopathogen Pseudomonas syringae MB03 on the model nematode Caenorhabditis elegans

The pathogenicity of a common phytopathogenic bacterium, Pseudomonas syringae, against animal model hosts, such as mice and Caenorhabditis elegans, has been recently revealed. However, most of the virulence determinants associated with pathogenesis remain elusive. In the current study, we performed predictive analysis of virulence factors against C. elegans in the genome of the wild-type P. syringae strain MB03. Nine predicted nematicidal proteins were expressed and purified in recombinant Escherichia coli strains and were evaluated to define their toxicity against C. elegans in liquid killing assays. Next, we focused on one essential 2-methyl citrate cycle protein, PrpF03, which showed the highest lethal activity against C. elegans compared to the other tested proteins with a half lethal concentration (LC₅₀) of 155.3 (123.4-176.6) µg mL^-1 and a half lethal time (LT₅₀) of 3.72 (1.64-4.85) days. Purified PrpF03 also caused adverse effects on the brood size, growth, and motility of C. elegans. Moreover, the PrpF03 protein exhibited pathological activity towards the intestinal tract of C. elegans. We surmise that the PrpF03 protein functions as a virulence factor when it blocks the average circulation of the 2-methyl citrate cycle of C. elegans by accumulating 2-methyl citrate in the gut of C. elegans, which damages and restrains the growth of intestinal tissues that ultimately kill C. elegans. The discovery of specific nematicidal activities of PrpF03 provides a better understanding of the mechanisms of phytopathogenic P. syringae against nematodes and could aid in developing nematode pest-controlling agents in agriculture.

WormPaths: Caenorhabditis elegans metabolic pathway annotation and visualization

In our group, we aim to understand metabolism in the nematode Caenorhabditis elegans and its relationships with gene expression, physiology, and the response to therapeutic drugs. Visualization of the metabolic pathways that comprise the metabolic network is extremely useful for interpreting a wide variety of experiments. Detailed annotated metabolic pathway maps for C. elegans are mostly limited to pan-organismal maps, many with incomplete or inaccurate pathway and enzyme annotations. Here, we present WormPaths, which is composed of two parts: (1) the careful manual annotation of metabolic genes into pathways, categories, and levels, and (2) 62 pathway maps that include metabolites, metabolite structures, genes, reactions, and pathway connections between maps. These maps are available on the WormFlux website. We show that WormPaths provides easy-to-navigate maps and that the different levels in WormPaths can be used for metabolic pathway enrichment analysis of transcriptomic data. In the future, we envision further developing these maps to be more interactive, analogous to road maps that are available on mobile devices.

WormPaths: Caenorhabditis elegans metabolic pathway annotation and visualization

In our group, we aim to understand metabolism in the nematode Caenorhabditis elegans and its relationships with gene expression, physiology and the response to therapeutic drugs. On March 15, 2020, a stay-at-home order was put into effect in the state of Massachusetts, USA, to flatten the curve of the spread of the novel SARS-CoV2 virus that causes COVID-19. For biomedical researchers in our state, this meant putting a hold on experiments for nine weeks until May 18, 2020. To keep the lab engaged and productive, and to enhance communication and collaboration, we embarked on an in-lab project that we all found important but that we never had the time for: the detailed annotation and drawing of C. elegans metabolic pathways. As a result, we present WormPaths, which is composed of two parts: 1) the careful manual annotation of metabolic genes into pathways, categories and levels, and 2) 66 pathway maps that include metabolites, metabolite structures, genes, reactions, and pathway connections between maps. These maps are available on our WormFlux website. We show that WormPaths provides easy-to-navigate maps and that the different levels in WormPaths can be used for metabolic pathway enrichment analysis of transcriptomic data. In the unfortunate event of additional lockdowns, we envision further developing these maps to be more interactive, with an analogy of road maps that are available on mobile devices.

Caenorhabditis elegans methionine/S-adenosylmethionine cycle activity is sensed and adjusted by a nuclear hormone receptor

Vitamin B12 is an essential micronutrient that functions in two metabolic pathways: the canonical propionate breakdown pathway and the methionine/S-adenosylmethionine (Met/SAM) cycle. In Caenorhabditis elegans, low vitamin B12, or genetic perturbation of the canonical propionate breakdown pathway results in propionate accumulation and the transcriptional activation of a propionate shunt pathway. This propionate-dependent mechanism requires nhr-10 and is referred to as ‘B12-mechanism-I’. Here, we report that vitamin B12 represses the expression of Met/SAM cycle genes by a propionate-independent mechanism we refer to as ‘B12-mechanism-II’. This mechanism is activated by perturbations in the Met/SAM cycle, genetically or due to low dietary vitamin B12. B12-mechanism-II requires nhr-114 to activate Met/SAM cycle gene expression, the vitamin B12 transporter, pmp-5, and adjust influx and efflux of the cycle by activating msra-1 and repressing cbs-1, respectively. Taken together, Met/SAM cycle activity is sensed and transcriptionally adjusted to be in a tight metabolic regime.

Modeling tissue-relevant Caenorhabditis elegans metabolism at network, pathway, reaction, and metabolite levels

Metabolism is a highly compartmentalized process that provides building blocks for biomass generation during development, homeostasis, and wound healing, and energy to support cellular and organismal processes. In metazoans, different cells and tissues specialize in different aspects of metabolism. However, studying the compartmentalization of metabolism in different cell types in a whole animal and for a particular stage of life is difficult. Here, we present MEtabolic models Reconciled with Gene Expression (MERGE), a computational pipeline that we used to predict tissue-relevant metabolic function at the network, pathway, reaction, and metabolite levels based on single-cell RNA-sequencing (scRNA-seq) data from the nematode Caenorhabditis elegans. Our analysis recapitulated known tissue functions in C. elegans, captured metabolic properties that are shared with similar tissues in human, and provided predictions for novel metabolic functions. MERGE is versatile and applicable to other systems. We envision this work as a starting point for the development of metabolic network models for individual cells as scRNA-seq continues to provide higher-resolution gene expression data.

Transcriptome Analysis of the Nematode Caenorhabditis elegans in Acidic Stress Environments

Ocean acidification and acid rain, caused by modern industries' fossil fuel burning, lead to a decrease in the living environmental pH, which results in a series of negative effects on many organisms. However, the underlying mechanisms of animals' response to acidic pH stress are largely unknown. In this study, we used the nematode Caenorhabditis elegans as an animal model to explore the regulatory mechanisms of organisms' response to pH decline. Two major stress-responsive pathways were found through transcriptome analysis in acidic stress environments. First, when the pH dropped from 6.33 to 4.33, the worms responded to the pH stress by upregulation of the col, nas, and dpy genes, which are required for cuticle synthesis and structure integrity. Second, when the pH continued to decrease from 4.33, the metabolism of xenobiotics by cytochrome P450 pathway genes (cyp, gst, ugt, and ABC transporters) played a major role in protecting the nematodes from the toxic substances probably produced by the more acidic environment. At the same time, the slowing down of cuticle synthesis might be due to its insufficient protective ability. Moreover, the systematic regulation pattern we found in nematodes might also be applied to other invertebrate and vertebrate animals to survive in the changing pH environments. Thus, our data might lay the foundation to identify the master gene(s) responding and adapting to acidic pH stress in further studies, and might also provide new solutions to improve assessment and monitoring of ecological restoration outcomes, or generate novel genotypes via genome editing for restoring in challenging environments especially in the context of acidic stress through global climate change.

A pathogen branched-chain amino acid catabolic pathway subverts host survival by impairing energy metabolism and the mitochondrial UPR

The mitochondrial unfolded protein response (UPRmt) is a stress-activated pathway promoting mitochondrial recovery and defense against infection. In C. elegans, the UPRmt is activated during infection with the pathogen Pseudomonas aeruginosa-but only transiently. As this may reflect a pathogenic strategy to target a pathway required for host survival, we conducted a P. aeruginosa genetic screen to uncover mechanisms associated with this temporary activation. Here, we find that loss of the P. aeruginosa acyl-CoA dehydrogenase FadE2 prolongs UPRmt activity and extends host survival. FadE2 shows substrate preferences for the coenzyme A intermediates produced during the breakdown of the branched-chain amino acids valine and leucine. Our data suggests that during infection, FadE2 restricts the supply of these catabolites to the host hindering host energy metabolism in addition to the UPRmt. Thus, a metabolic pathway in P. aeruginosa contributes to pathogenesis during infection through manipulation of host energy status and mitochondrial stress signaling potential.

Natural variation in a glucuronosyltransferase modulates propionate sensitivity in a C. elegans propionic acidemia model

Mutations in human metabolic genes can lead to rare diseases known as inborn errors of human metabolism. For instance, patients with loss-of-function mutations in either subunit of propionyl-CoA carboxylase suffer from propionic acidemia because they cannot catabolize propionate, leading to its harmful accumulation. Both the penetrance and expressivity of metabolic disorders can be modulated by genetic background. However, modifiers of these diseases are difficult to identify because of the lack of statistical power for rare diseases in human genetics. Here, we use a model of propionic acidemia in the nematode Caenorhabditis elegans to identify genetic modifiers of propionate sensitivity. Using genome-wide association (GWA) mapping across wild strains, we identify several genomic regions correlated with reduced propionate sensitivity. We find that natural variation in the putative glucuronosyltransferase GLCT-3, a homolog of human B3GAT, partly explains differences in propionate sensitivity in one of these genomic intervals. We demonstrate that loss-of-function alleles in glct-3 render the animals less sensitive to propionate. Additionally, we find that C. elegans has an expansion of the glct gene family, suggesting that the number of members of this family could influence sensitivity to excess propionate. Our findings demonstrate that natural variation in genes that are not directly associated with propionate breakdown can modulate propionate sensitivity. Our study provides a framework for using C. elegans to characterize the contributions of genetic background in models of human inborn errors in metabolism.

Niche-Selective Inhibition of Pathogenic Th17 Cells by Targeting Metabolic Redundancy

Targeting glycolysis has been considered therapeutically intractable owing to its essential housekeeping role. However, the context-dependent requirement for individual glycolytic steps has not been fully explored. We show that CRISPR-mediated targeting of glycolysis in T cells in mice results in global loss of Th17 cells, whereas deficiency of the glycolytic enzyme glucose phosphate isomerase (Gpi1) selectively eliminates inflammatory encephalitogenic and colitogenic Th17 cells, without substantially affecting homeostatic microbiota-specific Th17 cells. In homeostatic Th17 cells, partial blockade of glycolysis upon Gpi1 inactivation was compensated by pentose phosphate pathway flux and increased mitochondrial respiration. In contrast, inflammatory Th17 cells experience a hypoxic microenvironment known to limit mitochondrial respiration, which is incompatible with loss of Gpi1. Our study suggests that inhibiting glycolysis by targeting Gpi1 could be an effective therapeutic strategy with minimum toxicity for Th17-mediated autoimmune diseases, and, more generally, that metabolic redundancies can be exploited for selective targeting of disease processes.

Conserved nuclear hormone receptors controlling a novel plastic trait target fast-evolving genes expressed in a single cell

Environment shapes development through a phenomenon called developmental plasticity. Deciphering its genetic basis has potential to shed light on the origin of novel traits and adaptation to environmental change. However, molecular studies are scarce, and little is known about molecular mechanisms associated with plasticity. We investigated the gene regulatory network controlling predatory vs. non-predatory dimorphism in the nematode Pristionchus pacificus and found that it consists of genes of extremely different age classes. We isolated mutants in the conserved nuclear hormone receptor nhr-1 with previously unseen phenotypic effects. They disrupt mouth-form determination and result in animals combining features of both wild-type morphs. In contrast, mutants in another conserved nuclear hormone receptor nhr-40 display altered morph ratios, but no intermediate morphology. Despite divergent modes of control, NHR-1 and NHR-40 share transcriptional targets, which encode extracellular proteins that have no orthologs in Caenorhabditis elegans and result from lineage-specific expansions. An array of transcriptional reporters revealed co-expression of all tested targets in the same pharyngeal gland cell. Major morphological changes in this gland cell accompanied the evolution of teeth and predation, linking rapid gene turnover with morphological innovations. Thus, the origin of feeding plasticity involved novelty at the level of genes, cells and behavior.

WormCat: An Online Tool for Annotation and Visualization of Caenorhabditis elegans Genome-Scale Data

The emergence of large gene expression datasets has revealed the need for improved tools to identify enriched gene categories and visualize enrichment patterns. While gene ontogeny (GO) provides a valuable tool for gene set enrichment analysis, it has several limitations. First, it is difficult to graph multiple GO analyses for comparison. Second, genes from some model systems are not well represented. For example, ∼30% of Caenorhabditis elegans genes are missing from the analysis in commonly used databases. To allow categorization and visualization of enriched C. elegans gene sets in different types of genome-scale data, we developed WormCat, a web-based tool that uses a near-complete annotation of the C. elegans genome to identify coexpressed gene sets and scaled heat map for enrichment visualization. We tested the performance of WormCat using a variety of published transcriptomic datasets, and show that it reproduces major categories identified by GO. Importantly, we also found previously unidentified categories that are informative for interpreting phenotypes or predicting biological function. For example, we analyzed published RNA-seq data from C. elegans treated with combinations of lifespan-extending drugs, where one combination paradoxically shortened lifespan. Using WormCat, we identified sterol metabolism as a category that was not enriched in the single or double combinations, but emerged in a triple combination along with the lifespan shortening. Thus, WormCat identified a gene set with potential. phenotypic relevance not found with previous GO analysis. In conclusion, WormCat provides a powerful tool for the analysis and visualization of gene set enrichment in different types of C. elegans datasets.

The Computational Boundary of a "Self": Developmental Bioelectricity Drives Multicellularity and Scale-Free Cognition

All epistemic agents physically consist of parts that must somehow comprise an integrated cognitive self. Biological individuals consist of subunits (organs, cells, and molecular networks) that are themselves complex and competent in their own native contexts. How do coherent biological Individuals result from the activity of smaller sub-agents? To understand the evolution and function of metazoan creatures' bodies and minds, it is essential to conceptually explore the origin of multicellularity and the scaling of the basal cognition of individual cells into a coherent larger organism. In this article, I synthesize ideas in cognitive science, evolutionary biology, and developmental physiology toward a hypothesis about the origin of Individuality: "Scale-Free Cognition." I propose a fundamental definition of an Individual based on the ability to pursue goals at an appropriate level of scale and organization and suggest a formalism for defining and comparing the cognitive capacities of highly diverse types of agents. Any Self is demarcated by a computational surface - the spatio-temporal boundary of events that it can measure, model, and try to affect. This surface sets a functional boundary - a cognitive "light cone" which defines the scale and limits of its cognition. I hypothesize that higher level goal-directed activity and agency, resulting in larger cognitive boundaries, evolve from the primal homeostatic drive of living things to reduce stress - the difference between current conditions and life-optimal conditions. The mechanisms of developmental bioelectricity - the ability of all cells to form electrical networks that process information - suggest a plausible set of gradual evolutionary steps that naturally lead from physiological homeostasis in single cells to memory, prediction, and ultimately complex cognitive agents, via scale-up of the basic drive of infotaxis. Recent data on the molecular mechanisms of pre-neural bioelectricity suggest a model of how increasingly sophisticated cognitive functions emerge smoothly from cell-cell communication used to guide embryogenesis and regeneration. This set of hypotheses provides a novel perspective on numerous phenomena, such as cancer, and makes several unique, testable predictions for interdisciplinary research that have implications not only for evolutionary developmental biology but also for biomedicine and perhaps artificial intelligence and exobiology.