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.

Vitamin B12 produced by gut bacteria modulates cholinergic signalling.

A growing body of evidence indicates that gut microbiota influence brain function and behaviour. However, the molecular basis of how gut bacteria modulate host nervous system function is largely unknown. Here we show that vitamin B12-producing bacteria that colonize the intestine can modulate excitatory cholinergic signalling and behaviour in the host Caenorhabditis elegans. Here we demonstrate that vitamin B12 reduces cholinergic signalling in the nervous system through rewiring of the methionine (Met)/S-adenosylmethionine cycle in the intestine. We identify a conserved metabolic crosstalk between the methionine/S-adenosylmethionine cycle and the choline-oxidation pathway. In addition, we show that metabolic rewiring of these pathways by vitamin B12 reduces cholinergic signalling by limiting the availability of free choline required by neurons to synthesize acetylcholine. Our study reveals a gut-brain communication pathway by which enteric bacteria modulate host behaviour and may affect neurological health.

A bacterial pathogen induces developmental slowing by high reactive oxygen species and mitochondrial dysfunction in Caenorhabditis elegans.

Host-pathogen interactions are complex by nature, and the host developmental stage increases this complexity. By utilizing Caenorhabditis elegans larvae as the host and the bacterium Pseudomonas aeruginosa as the pathogen, we investigated how a developing organism copes with pathogenic stress. By screening 36 P. aeruginosa isolates, we found that the CF18 strain causes a severe but reversible developmental delay via induction of reactive oxygen species (ROS) and mitochondrial dysfunction. While the larvae upregulate mitophagy, antimicrobial, and detoxification genes, mitochondrial unfolded protein response (UPRmt) genes are repressed. Either antioxidant or iron supplementation rescues the phenotypes. We examined the virulence factors of CF18 via transposon mutagenesis and RNA sequencing (RNA-seq). We found that non-phenazine toxins that are regulated by quorum sensing (QS) and the GacA/S system are responsible for developmental slowing. This study highlights the importance of ROS levels and mitochondrial health as determinants of developmental rate and how pathogens can attack these important features.

Disruption of sugar nucleotide clearance is a therapeutic vulnerability of cancer cells.

Identifying metabolic steps that are specifically required for the survival of cancer cells but are dispensable in normal cells remains a challenge1. Here we report a therapeutic vulnerability in a sugar nucleotide biosynthetic pathway that can be exploited in cancer cells with only a limited impact on normal cells. A systematic examination of conditionally essential metabolic enzymes revealed that UXS1, a Golgi enzyme that converts one sugar nucleotide (UDP-glucuronic acid, UDPGA) to another (UDP-xylose), is essential only in cells that express high levels of the enzyme immediately upstream of it, UGDH. This conditional relationship exists because UXS1 is required to prevent excess accumulation of UDPGA, which is produced by UGDH. UXS1 not only clears away UDPGA but also limits its production through negative feedback on UGDH. Excess UDPGA disrupts Golgi morphology and function, which impedes the trafficking of surface receptors such as EGFR to the plasma membrane and diminishes the signalling capacity of cells. UGDH expression is elevated in several cancers, including lung adenocarcinoma, and is further enhanced during chemoresistant selection. As a result, these cancer cells are selectively dependent on UXS1 for UDPGA detoxification, revealing a potential weakness in tumours with high levels of UGDH.

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.

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 regulate different promoters, and that both essential and intestinally expressed CFs have the greatest effects on promoter activity. 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.

Evolved bacterial resistance to the chemotherapy gemcitabine modulates its efficacy in co-cultured cancer cells.

Drug metabolism by the microbiome can influence anticancer treatment success. We previously suggested that chemotherapies with antimicrobial activity can select for adaptations in bacterial drug metabolism that can inadvertently influence the host's chemoresistance. We demonstrated that evolved resistance against fluoropyrimidine chemotherapy lowered its efficacy in worms feeding on drug-evolved bacteria (Rosener et al., 2020). Here, we examine a model system that captures local interactions that can occur in the tumor microenvironment. Gammaproteobacteria-colonizing pancreatic tumors can degrade the nucleoside-analog chemotherapy gemcitabine and, in doing so, can increase the tumor's chemoresistance. Using a genetic screen in Escherichia coli, we mapped all loss-of-function mutations conferring gemcitabine resistance. Surprisingly, we infer that one third of top resistance mutations increase or decrease bacterial drug breakdown and therefore can either lower or raise the gemcitabine load in the local environment. Experiments in three E. coli strains revealed that evolved adaptation converged to inactivation of the nucleoside permease NupC, an adaptation that increased the drug burden on co-cultured cancer cells. The two studies provide complementary insights on the potential impact of microbiome adaptation to chemotherapy by showing that bacteria-drug interactions can have local and systemic influence on drug activity.

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.

Bacterial diet modulates tamoxifen-induced death via host fatty acid metabolism.

Tamoxifen is a selective estrogen receptor (ER) modulator that is used to treat ER-positive breast cancer, but that at high doses kills both ER-positive and ER-negative breast cancer cells. We recapitulate this off-target effect in Caenorhabditis elegans, which does not have an ER ortholog. We find that different bacteria dramatically modulate tamoxifen toxicity in C. elegans, with a three-order of magnitude difference between animals fed Escherichia coli, Comamonas aquatica, and Bacillus subtilis. Remarkably, host fatty acid (FA) biosynthesis mitigates tamoxifen toxicity, and different bacteria provide the animal with different FAs, resulting in distinct FA profiles. Surprisingly these bacteria modulate tamoxifen toxicity by different death mechanisms, some of which are modulated by FA supplementation and others by antioxidants. Together, this work reveals a complex interplay between microbiota, FA metabolism and tamoxifen toxicity that may provide a blueprint for similar studies in more complex mammals.

C. elegans as a model for inter-individual variation in metabolism.

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors1-4. It is difficult to connect differences in metabolism to genomic variation and derive underlying molecular mechanisms in humans, owing to differences in diet and lifestyle, among others. Here we use the nematode Caenorhabditis elegans as a model to study inter-individual variation in metabolism. By comparing three wild strains and the commonly used N2 laboratory strain, we find differences in the abundances of both known metabolites and those that have not to our knowledge been previously described. The latter metabolites include conjugates between 3-hydroxypropionate (3HP) and several amino acids (3HP-AAs), which are much higher in abundance in one of the wild strains. 3HP is an intermediate in the propionate shunt pathway, which is activated when flux through the canonical, vitamin-B12-dependent propionate breakdown pathway is perturbed5. We show that increased accumulation of 3HP-AAs is caused by genetic variation in HPHD-1, for which 3HP is a substrate. Our results suggest that the production of 3HP-AAs represents a 'shunt-within-a-shunt' pathway to accommodate a reduction-of-function allele in hphd-1. This study provides a step towards the development of metabolic network models that capture individual-specific differences of metabolism and more closely represent the diversity that is found in entire species.

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.

The Plasmid pEX18Gm Indirectly Increases Caenorhabditis elegans Fecundity by Accelerating Bacterial Methionine Synthesis.

Plasmids are mostly found in bacteria as extrachromosomal genetic elements and are widely used in genetic engineering. Exploring the mechanisms of plasmid-host interaction can provide crucial information for the application of plasmids in genetic engineering. However, many studies have generally focused on the influence of plasmids on their bacterial hosts, and the effects of plasmids on bacteria-feeding animals have not been explored in detail. Here, we use a "plasmid-bacteria-Caenorhabditis elegans" model to explore the impact of plasmids on their host bacteria and bacterivorous nematodes. First, the phenotypic responses of C. elegans were observed by feeding Escherichia coli OP50 harboring different types of plasmids. We found that E. coli OP50 harboring plasmid pEX18Gm unexpectedly increases the fecundity of C. elegans. Subsequently, we found that the plasmid pEX18Gm indirectly affects C. elegans fecundity via bacterial metabolism. To explore the underlying regulatory mechanism, we performed bacterial RNA sequencing and performed in-depth analysis. We demonstrated that the plasmid pEX18Gm upregulates the transcription of methionine synthase gene metH in the bacteria, which results in an increase in methionine that supports C. elegans fecundity. Additionally, we found that a pEX18Gm-induced increase in C. elegans can occur in different bacterial species. Our findings highlight the plasmid-bacteria-C. elegans model to reveal the mechanism of plasmids' effects on their host and provide a new pattern for systematically studying the interaction between plasmids and multi-species.

Chromosome-Level Reference Genomes for Two Strains of Caenorhabditis briggsae: An Improved Platform for Comparative Genomics.

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome-conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Compartmentalization of metabolism between cell types in multicellular organisms: a computational perspective.

In multicellular organisms, metabolism is compartmentalized at many levels, including tissues and organs, different cell types, and subcellular compartments. Compartmentalization creates a coordinated homeostatic system where each compartment contributes to the production of energy and biomolecules the organism needs to carrying out specific metabolic tasks. Experimentally studying metabolic compartmentalization and metabolic interactions between cells and tissues in multicellular organisms is challenging at a systems level. However, recent progress in computational modeling provides an alternative approach to this problem. Here we discuss how integrating metabolic network modeling with omics data offers an opportunity to reveal metabolic states at the level of organs, tissues and, ultimately, individual cells. We review the current status of genome-scale metabolic network models in multicellular organisms, methods to study metabolic compartmentalization in silico, and insights gained from computational analyses. We also discuss outstanding challenges and provide perspectives for the future directions of the field.

Monomethyl branched-chain fatty acids are critical for Caenorhabitis elegans survival in elevated glucose conditions.

The maintenance of optimal membrane composition under basal and stress conditions is critical for the survival of an organism. High-glucose stress has been shown to perturb membrane properties by decreasing membrane fluidity, and the membrane sensor PAQR-2 is required to restore membrane integrity. However, the mechanisms required to respond to elevated dietary glucose are not fully established. In this study, we used a 13C stable isotope-enriched diet and mass spectrometry to better understand the impact of glucose on fatty acid dynamics in the membrane of Caenorhabditis elegans. We found a novel role for monomethyl branched-chain fatty acids (mmBCFAs) in mediating the ability of the nematodes to survive conditions of elevated dietary glucose. This requirement of mmBCFAs is unique to glucose stress and was not observed when the nematode was fed elevated dietary saturated fatty acid. In addition, when worms deficient in elo-5, the major biosynthesis enzyme of mmBCFAs, were fed Bacillus subtilis (a bacteria strain rich in mmBCFAs) in combination with high glucose, their survival rates were rescued to wild-type levels. Finally, the results suggest that mmBCFAs are part of the PAQR-2 signaling response during glucose stress. Taken together, we have identified a novel role for mmBCFAs in stress response in nematodes and have established these fatty acids as critical for adapting to elevated glucose.

Essential role of a ThPOK autoregulatory loop in the maintenance of mature CD4+ T cell identity and function.

The transcription factor ThPOK (encoded by the Zbtb7b gene) controls homeostasis and differentiation of mature helper T cells, while opposing their differentiation to CD4+ intraepithelial lymphocytes (IELs) in the intestinal mucosa. Thus CD4 IEL differentiation requires ThPOK transcriptional repression via reactivation of the ThPOK transcriptional silencer element (SilThPOK). In the present study, we describe a new autoregulatory loop whereby ThPOK binds to the SilThPOK to maintain its own long-term expression in CD4 T cells. Disruption of this loop in vivo prevents persistent ThPOK expression, leads to genome-wide changes in chromatin accessibility and derepresses the colonic regulatory T (Treg) cell gene expression signature. This promotes selective differentiation of naive CD4 T cells into GITRloPD-1loCD25lo (Triplelo) Treg cells and conversion to CD4+ IELs in the gut, thereby providing dominant protection from colitis. Hence, the ThPOK autoregulatory loop represents a key mechanism to physiologically control ThPOK expression and T cell differentiation in the gut, with potential therapeutic relevance.

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.

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.

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.

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.