Caenorhabditis elegans metabolic gene regulatory networks govern the cellular economy

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.

Non-canonical pattern recognition of a pathogen-derived metabolite by a nuclear hormone receptor identifies virulent bacteria in C. elegans

Distinguishing infectious pathogens from harmless microorganisms is essential for animal health. The mechanisms used to identify infectious microbes are not fully understood, particularly in metazoan hosts that eat bacteria as their food source. Here, we characterized a non-canonical pattern-recognition system in Caenorhabditis elegans (C. elegans) that assesses the relative threat of virulent Pseudomonas aeruginosa (P. aeruginosa) to activate innate immunity. We discovered that the innate immune response in C. elegans was triggered by phenazine-1-carboxamide (PCN), a toxic metabolite produced by pathogenic strains of P. aeruginosa. We identified the nuclear hormone receptor NHR-86/HNF4 as the PCN sensor in C. elegans and validated that PCN bound to the ligand-binding domain of NHR-86/HNF4. Activation of NHR-86/HNF4 by PCN directly engaged a transcriptional program in intestinal epithelial cells that protected against P. aeruginosa. Thus, a bacterial metabolite is a pattern of pathogenesis surveilled by nematodes to identify a pathogen in its bacterial diet.

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.

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.

Bacterial diets differentially alter lifespan and healthspan trajectories in C. elegans

Diet is one of the more variable aspects in life due to the variety of options that organisms are exposed to in their natural habitats. In the laboratory, C. elegans are raised on bacterial monocultures, traditionally the E. coli B strain OP50, and spontaneously occurring microbial contaminants are removed to limit experimental variability because diet-including the presence of contaminants-can exert a potent influence over animal physiology. In order to diversify the menu available to culture C. elegans in the lab, we have isolated and cultured three such microbes: Methylobacterium, Xanthomonas, and Sphingomonas. The nutritional composition of these bacterial foods is unique, and when fed to C. elegans, can differentially alter multiple life history traits including development, reproduction, and metabolism. In light of the influence each food source has on specific physiological attributes, we comprehensively assessed the impact of these bacteria on animal health and devised a blueprint for utilizing different food combinations over the lifespan, in order to promote longevity. The expansion of the bacterial food options to use in the laboratory will provide a critical tool to better understand the complexities of bacterial diets and subsequent changes in physiology and gene expression.

Hierarchical and Dynamic Regulation of Defense-Responsive Specialized Metabolism by WRKY and MYB Transcription Factors

The plant kingdom produces hundreds of thousands of specialized bioactive metabolites, some with pharmaceutical and biotechnological importance. Their biosynthesis and function have been studied for decades, but comparatively less is known about how transcription factors with overlapping functions and contrasting regulatory activities coordinately control the dynamics and output of plant specialized metabolism. Here, we performed temporal studies on pathogen-infected intact host plants with perturbed transcription factors. We identified WRKY33 as the condition-dependent master regulator and MYB51 as the dual functional regulator in a hierarchical gene network likely responsible for the gene expression dynamics and metabolic fluxes in the camalexin and 4-hydroxy-indole-3-carbonylnitrile (4OH-ICN) pathways. This network may have also facilitated the regulatory capture of the newly evolved 4OH-ICN pathway in Arabidopsis thaliana by the more-conserved transcription factor MYB51. It has long been held that the plasticity of plant specialized metabolism and the canalization of development should be differently regulated; our findings imply a common hierarchical regulatory architecture orchestrated by transcription factors for specialized metabolism and development, making it an attractive target for metabolic engineering.

Metabolic effects of manganese in the nematode Caenorhabditis elegans through DAergic pathway and transcription factors activation

Manganese (Mn) is an essential trace element for physiological functions since it acts as an enzymatic co-factor. Nevertheless, overexposure to Mn has been associated with a pathologic condition called manganism. Furthermore, Mn has been reported to affect lipid metabolism by mechanisms which have yet to be established. Herein, we used the nematode Caenorhabditis elegans to examine Mn's effects on the dopaminergic (DAergic) system and determine which transcription factors that regulate with lipid metabolism are affected by it. Worms were exposed to Mn for four hours in the presence of bacteria and in a liquid medium (85 mM NaCl). Mn increased fat storage as evidenced both by Oil Red O accumulation and triglyceride levels. In addition, metabolic activity was reduced as a reflection of decreased oxygen consumption caused by Mn. Mn also affected feeding behavior as evidenced by decreased pharyngeal pumping rate. DAergic neurons viability were not altered by Mn, however the dopamine levels were significantly reduced following Mn exposure. Furthermore, the expression of sbp-1 transcription factor and let-363 protein kinase responsible for lipid accumulation control was increased and decreased, respectively, by Mn. Altogether, our data suggest that Mn increases the fat storage in C. elegans, secondary to DAergic system alterations, under the control of SBP-1 and LET-363 proteins.

Forward and reverse genetics approaches to uncover metabolic aging pathways in Caenorhabditis elegans

The biological mechanisms of aging have been studied in depth and prominent findings in this field promote the development of new therapies for age-associated disorders. Various model organisms are used for research on aging; among these, the nematode Caenorhabditis elegans has been widely used and has provided valuable knowledge in determining the regulatory mechanisms driving the aging process. Many genes involved in lifespan regulation are associated with metabolic pathways and are influenced by genetic and environmental factors. In line with this, C. elegans provides a promising platform to study such gene by environment interactions, in either a reverse or forward genetics approach. In this review, we discuss longevity mechanisms related to metabolic networks that have been discovered in C. elegans. We also highlight the use of wild populations to study the complex genetic basis of natural variation for quantitative traits that mediate longevity.

The Natural Biotic Environment of Caenorhabditis elegans

Organisms evolve in response to their natural environment. Consideration of natural ecological parameters are thus of key importance for our understanding of an organism's biology. Curiously, the natural ecology of the model species Caenorhabditis elegans has long been neglected, even though this nematode has become one of the most intensively studied models in biological research. This lack of interest changed ∼10 yr ago. Since then, an increasing number of studies have focused on the nematode's natural ecology. Yet many unknowns still remain. Here, we provide an overview of the currently available information on the natural environment of C. elegans We focus on the biotic environment, which is usually less predictable and thus can create high selective constraints that are likely to have had a strong impact on C. elegans evolution. This nematode is particularly abundant in microbe-rich environments, especially rotting plant matter such as decomposing fruits and stems. In this environment, it is part of a complex interaction network, which is particularly shaped by a species-rich microbial community. These microbes can be food, part of a beneficial gut microbiome, parasites and pathogens, and possibly competitors. C. elegans is additionally confronted with predators; it interacts with vector organisms that facilitate dispersal to new habitats, and also with competitors for similar food environments, including competitors from congeneric and also the same species. Full appreciation of this nematode's biology warrants further exploration of its natural environment and subsequent integration of this information into the well-established laboratory-based research approaches.

Bacterial Metabolism Affects the C. elegans Response to Cancer Chemotherapeutics

The human microbiota greatly affects physiology and disease; however, the contribution of bacteria to the response to chemotherapeutic drugs remains poorly understood. Caenorhabditis elegans and its bacterial diet provide a powerful system to study host-bacteria interactions. Here, we use this system to study how bacteria affect the C. elegans response to chemotherapeutics. We find that different bacterial species can increase the response to one drug yet decrease the effect of another. We perform genetic screens in two bacterial species using three chemotherapeutic drugs: 5-fluorouracil (5-FU), 5-fluoro-2'-deoxyuridine (FUDR), and camptothecin (CPT). We find numerous bacterial nucleotide metabolism genes that affect drug efficacy in C. elegans. Surprisingly, we find that 5-FU and FUDR act through bacterial ribonucleotide metabolism to elicit their cytotoxic effects in C. elegans rather than by thymineless death or DNA damage. Our study provides a blueprint for characterizing the role of bacteria in the host response to chemotherapeutics.

Long-term time-lapse microscopy of C. elegans post-embryonic development

We present a microscopy technique that enables long-term time-lapse microscopy at single-cell resolution in moving and feeding Caenorhabditis elegans larvae. Time-lapse microscopy of C. elegans post-embryonic development is challenging, as larvae are highly motile. Moreover, immobilization generally leads to rapid developmental arrest. Instead, we confine larval movement to microchambers that contain bacteria as food, and use fast image acquisition and image analysis to follow the dynamics of cells inside individual larvae, as they move within each microchamber. This allows us to perform fluorescence microscopy of 10-20 animals in parallel with 20 min time resolution. We demonstrate the power of our approach by analysing the dynamics of cell division, cell migration and gene expression over the full ∼48 h of development from larva to adult. Our approach now makes it possible to study the behaviour of individual cells inside the body of a feeding and growing animal.

A Caenorhabditis elegans Genome-Scale Metabolic Network Model

Caenorhabditis elegans is a powerful model to study metabolism and how it relates to nutrition, gene expression, and life history traits. However, while numerous experimental techniques that enable perturbation of its diet and gene function are available, a high-quality metabolic network model has been lacking. Here, we reconstruct an initial version of the C. elegans metabolic network. This network model contains 1,273 genes, 623 enzymes, and 1,985 metabolic reactions and is referred to as iCEL1273. Using flux balance analysis, we show that iCEL1273 is capable of representing the conversion of bacterial biomass into C. elegans biomass during growth and enables the predictions of gene essentiality and other phenotypes. In addition, we demonstrate that gene expression data can be integrated with the model by comparing metabolic rewiring in dauer animals versus growing larvae. iCEL1273 is available at a dedicated website (wormflux.umassmed.edu) and will enable the unraveling of the mechanisms by which different macro- and micronutrients contribute to the animal's physiology.

Folic acid supplementation at lower doses increases oxidative stress resistance and longevity in Caenorhabditis elegans

Folic acid (FA) is an essential nutrient that the human body needs but cannot be synthesized on its own. Fortified foods and plant food sources such as green leafy vegetables, beans, fruits, and juices are good sources of FA to meet the daily requirements of the body. The aim was to evaluate the effect of dietary FA levels on the longevity of well-known experimental aging model Caenorhabditis elegans. Here, we show for first time that FA extends organism life span and causes a delay in aging. We observed that FA inhibits mechanistic target of rapamycin (mTOR) and insulin/insulin growth factor 1 (IGF-1) signaling pathways to control both oxidative stress levels and life span. The expression levels of stress- and life span-relevant gerontogenes, viz. daf-16, skn-1, and sir. 2.1, and oxidative enzymes, such as glutathione S-transferase 4 (GST-4) and superoxide dismutase 3 (SOD-3), were also found to be highly enhanced to attenuate the intracellular reactive oxygen species (ROS) damage and to delay the aging process. Our study promotes the use of FA to mitigate abiotic stresses and other aging-related ailments.

Cell-Nonautonomous Mechanisms Underlying Cellular and Organismal Aging

Cell-autonomous mechanisms underlying cellular and organismal aging in evolutionarily distant eukaryotes have been established; these mechanisms regulate longevity-defining processes within a single eukaryotic cell. Recent findings have provided valuable insight into cell-nonautonomous mechanisms modulating cellular and organismal aging in eukaryotes across phyla; these mechanisms involve a transmission of various longevity factors between different cells, tissues, and organisms. Herein, we review such cell-nonautonomous mechanisms of aging in eukaryotes. We discuss the following: (1) how low molecular weight transmissible longevity factors modulate aging and define longevity of cells in yeast populations cultured in liquid media or on solid surfaces, (2) how communications between proteostasis stress networks operating in neurons and nonneuronal somatic tissues define longevity of the nematode Caenorhabditis elegans by modulating the rates of aging in different tissues, and (3) how different bacterial species colonizing the gut lumen of C. elegans define nematode longevity by modulating the rate of organismal aging.

Hierarchical decomposition of dynamically evolving regulatory networks

Background

Gene regulatory networks describe the interplay between genes and their products. These networks control almost every biological activity in the cell through interactions. The hierarchy of genes in these networks as defined by their interactions gives important insights into how these functions are governed. Accurately determining the hierarchy of genes is however a computationally difficult problem. This problem is further complicated by the fact that an intrinsic characteristic of regulatory networks is that the wiring of interactions can change over time. Determining how the hierarchy in the gene regulatory networks changes with dynamically evolving network topology remains to be an unsolved challenge.

Results

In this study, we develop a new method, named D-HIDEN (Dynamic-HIerarchical DEcomposition of Networks) to find the hierarchy of the genes in dynamically evolving gene regulatory network topologies. Unlike earlier methods, which recompute the hierarchy from scratch when the network topology changes, our method adapts the hierarchy based on the wiring of the interactions only for the nodes which have the potential to move in the hierarchy.

Conclusions

We compare D-HIDEN to five currently available hierarchical decomposition methods on synthetic and real gene regulatory networks. Our experiments demonstrate that D-HIDEN significantly outperforms existing methods in running time, accuracy, or both. Furthermore, our method is robust against dynamic changes in hierarchy. Our experiments on human gene regulatory networks suggest that our method may be used to reconstruct hierarchy in gene regulatory networks.

RNAi Interrogation of Dietary Modulation of Development, Metabolism, Behavior, and Aging in C. elegans

Diet affects nearly every aspect of animal life such as development, metabolism, behavior, and aging, both directly by supplying nutrients and indirectly through gut microbiota. C. elegans feeds on bacteria, and like other animals, different bacterial diets induce distinct dietary responses in the worm. However, the lack of certain critical tools hampers the use of worms as a model for dietary signaling. Here, we genetically engineered the bacterial strain OP50, the standard laboratory diet for C. elegans, making it compatible for dsRNA production and delivery. Using this RNAi-compatible OP50 strain and the other bacterial strain HT115, we feed worms different diets while delivering RNAi to interrogate the genetic basis underlying diet-dependent differential modulation of development, metabolism, behavior, and aging. We show by RNAi that neuroendocrine and mTOR pathways are involved in mediating differential dietary responses. This genetic tool greatly facilitates the use of C. elegans as a model for dietary signaling.

Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation in Caenorhabditis elegans intestine and muscles

BACKGROUND

Tissue-specific RNA plasticity broadly impacts the development, tissue identity and adaptability of all organisms, but changes in composition, expression levels and its impact on gene regulation in different somatic tissues are largely unknown. Here we developed a new method, polyA-tagging and sequencing (PAT-Seq) to isolate high-quality tissue-specific mRNA from Caenorhabditis elegans intestine, pharynx and body muscle tissues and study changes in their tissue-specific transcriptomes and 3'UTRomes.

RESULTS

We have identified thousands of novel genes and isoforms differentially expressed between these three tissues. The intestine transcriptome is expansive, expressing over 30% of C. elegans mRNAs, while muscle transcriptomes are smaller but contain characteristic unique gene signatures. Active promoter regions in all three tissues reveal both known and novel enriched tissue-specific elements, along with putative transcription factors, suggesting novel tissue-specific modes of transcription initiation. We have precisely mapped approximately 20,000 tissue-specific polyadenylation sites and discovered that about 30% of transcripts in somatic cells use alternative polyadenylation in a tissue-specific manner, with their 3'UTR isoforms significantly enriched with microRNA targets.

CONCLUSIONS

For the first time, PAT-Seq allowed us to directly study tissue specific gene expression changes in an in vivo setting and compare these changes between three somatic tissues from the same organism at single-base resolution within the same experiment. We pinpoint precise tissue-specific transcriptome rearrangements and for the first time link tissue-specific alternative polyadenylation to miRNA regulation, suggesting novel and unexplored tissue-specific post-transcriptional regulatory networks in somatic cells.

Promoter-based integration in plant defense regulation

A key unanswered question in plant biology is how a plant regulates metabolism to maximize performance across an array of biotic and abiotic environmental stresses. In this study, we addressed the potential breadth of transcriptional regulation that can alter accumulation of the defensive glucosinolate metabolites in Arabidopsis (Arabidopsis thaliana). A systematic yeast one-hybrid study was used to identify hundreds of unique potential regulatory interactions with a nearly complete complement of 21 promoters for the aliphatic glucosinolate pathway. Conducting high-throughput phenotypic validation, we showed that >75% of tested transcription factor (TF) mutants significantly altered the accumulation of the defensive glucosinolates. These glucosinolate phenotypes were conditional upon the environment and tissue type, suggesting that these TFs may allow the plant to tune its defenses to the local environment. Furthermore, the pattern of TF/promoter interactions could partially explain mutant phenotypes. This work shows that defense chemistry within Arabidopsis has a highly intricate transcriptional regulatory system that may allow for the optimization of defense metabolite accumulation across a broad array of environments.