Several molecules can extend healthspan and lifespan across organisms. However, most are upstream signaling hubs or transcription factors orchestrating complex anti-aging programs. Therefore, these molecules point to but do not reveal the fundamental mechanisms driving longevity. Instead, downstream effectors that are necessary and sufficient to promote longevity across conditions or organisms may reveal the fundamental anti-aging drivers. Toward this goal, we searched for effectors acting downstream of the transcription factor EB (TFEB), known as HLH-30 in C. elegans, because TFEB/HLH-30 is necessary across anti-aging interventions and its overexpression is sufficient to extend C. elegans lifespan and reduce biomarkers of aging in mammals including humans. As a result, we present an alcohol-dehydrogenase-mediated anti-aging response (AMAR) that is essential for C. elegans longevity driven by HLH-30 overexpression, caloric restriction, mTOR inhibition, and insulin-signaling deficiency. The sole overexpression of ADH-1 is sufficient to activate AMAR, which extends healthspan and lifespan by reducing the levels of glycerol-an age-associated and aging-promoting alcohol. Adh1 overexpression is also sufficient to promote longevity in yeast, and adh-1 orthologs are induced in calorically restricted mice and humans, hinting at ADH-1 acting as an anti-aging effector across phyla.

Cell-cell interactions shape cellular function and ultimately organismal phenotype. Interacting cells can sense their mutual distance using combinations of ligand-receptor pairs, suggesting the existence of a spatial code, i.e., signals encoding spatial properties of cellular organization. However, this code driving and sustaining the spatial organization of cells remains to be elucidated. Here we present a computational framework to infer the spatial code underlying cell-cell interactions from the transcriptomes of the cell types across the whole body of a multicellular organism. As core of this framework, we introduce our tool cell2cell, which uses the coexpression of ligand-receptor pairs to compute the potential for intercellular interactions, and we test it across the Caenorhabditis elegans' body. Leveraging a 3D atlas of C. elegans' cells, we also implement a genetic algorithm to identify the ligand-receptor pairs most informative of the spatial organization of cells across the whole body. Validating the spatial code extracted with this strategy, the resulting intercellular distances are negatively correlated with the inferred cell-cell interactions. Furthermore, for selected cell-cell and ligand-receptor pairs, we experimentally confirm the communicatory behavior inferred with cell2cell and the genetic algorithm. Thus, our framework helps identify a code that predicts the spatial organization of cells across a whole-animal body.

We present a toolbox for high-throughput screening of image-based Caenorhabditis elegans phenotypes. The image analysis algorithms measure morphological phenotypes in individual worms and are effective for a variety of assays and imaging systems. This WormToolbox is available through the open-source CellProfiler project and enables objective scoring of whole-worm high-throughput image-based assays of C. elegans for the study of diverse biological pathways that are relevant to human disease.

DNA glycosylases catalyze the excision of chemically modified bases from DNA. Although most glycosylases are specific to a particular base, the 3-methyladenine (m3A) DNA glycosylases include both highly specific enzymes acting on a single modified base, and enzymes with broader specificity for alkylation-damaged DNA. Our structural understanding of these different enzymatic specificities is currently limited to crystal and NMR structures of the unliganded enzymes and complexes with abasic DNA inhibitors. Presented here are high-resolution crystal structures of the m3A DNA glycosylase from Helicobacter pylori (MagIII) in the unliganded form and bound to alkylated bases 3,9-dimethyladenine and 1,N6-ethenoadenine. These are the first structures of a nucleobase bound in the active site of a m3A glycosylase belonging to the helix-hairpin-helix superfamily. MagIII achieves its specificity for positively-charged m3A not by direct interactions with purine or methyl substituent atoms, but rather by stacking the base between two aromatic side chains in a pocket that excludes 7-methylguanine. We report base excision and DNA binding activities of MagIII active site mutants, together with a structural comparison of the HhH glycosylases.

The gut microbiota metabolizes drugs and alters their efficacy and toxicity. Diet alters drugs, the metabolism of the microbiota, and the host. However, whether diet-triggered metabolic changes in the microbiota can alter drug responses in the host has been largely unexplored. Here we show that dietary thymidine and serine enhance 5-fluoro 2'deoxyuridine (FUdR) toxicity in C. elegans through different microbial mechanisms. Thymidine promotes microbial conversion of the prodrug FUdR into toxic 5-fluorouridine-5'-monophosphate (FUMP), leading to enhanced host death associated with mitochondrial RNA and DNA depletion, and lethal activation of autophagy. By contrast, serine does not alter FUdR metabolism. Instead, serine alters E. coli's 1C-metabolism, reduces the provision of nucleotides to the host, and exacerbates DNA toxicity and host death without mitochondrial RNA or DNA depletion; moreover, autophagy promotes survival in this condition. This work implies that diet-microbe interactions can alter the host response to drugs without altering the drug or the host.

Obesity and its associated metabolic syndrome are a leading cause of morbidity and mortality. Given the disease's heavy burden on patients and the healthcare system, there has been increased interest in identifying pharmacological targets for the treatment and prevention of obesity. Towards this end, genome-wide association studies (GWAS) have identified hundreds of human genetic variants associated with obesity. The next challenge is to experimentally define which of these variants are causally linked to obesity, and could therefore become targets for the treatment or prevention of obesity. Here we employ high-throughput in vivo RNAi screening to test for causality 293 C. elegans orthologs of human obesity-candidate genes reported in GWAS. We RNAi screened these 293 genes in C. elegans subject to two different feeding regimens: (1) regular diet, and (2) high-fructose diet, which we developed and present here as an invertebrate model of diet-induced obesity (DIO). We report 14 genes that promote obesity and 3 genes that prevent DIO when silenced in C. elegans. Further, we show that knock-down of the 3 DIO genes not only prevents excessive fat accumulation in primary and ectopic fat depots but also improves the health and extends the lifespan of C. elegans overconsuming fructose. Importantly, the direction of the association between expression variants in these loci and obesity in mice and humans matches the phenotypic outcome of the loss-of-function of the C. elegans ortholog genes, supporting the notion that some of these genes would be causally linked to obesity across phylogeny. Therefore, in addition to defining causality for several genes so far merely correlated with obesity, this study demonstrates the value of model systems compatible with in vivo high-throughput genetic screening to causally link GWAS gene candidates to human diseases.

Fat accumulation is a complex phenotype affected by factors such as neuroendocrine signaling, feeding, activity, and reproductive output. Accordingly, the most informative screens for genes and compounds affecting fat accumulation would be those carried out in whole living animals. Caenorhabditis elegans is a well-established and effective model organism, especially for biological processes that involve organ systems and multicellular interactions, such as metabolism. Every cell in the transparent body of C. elegans is visible under a light microscope. Consequently, an accessible and reliable method to visualize worm lipid-droplet fat depots would make C. elegans the only metazoan in which genes affecting not only fat mass but also body fat distribution could be assessed at a genome-wide scale. Here we present a radical improvement in oil red O worm staining together with high-throughput image-based phenotyping. The three-step sample preparation method is robust, formaldehyde-free, and inexpensive, and requires only 15min of hands-on time to process a 96-well plate. Together with our free and user-friendly automated image analysis package, this method enables C. elegans sample preparation and phenotype scoring at a scale that is compatible with genome-wide screens. Thus we present a feasible approach to small-scale phenotyping and large-scale screening for genetic and/or chemical perturbations that lead to alterations in fat quantity and distribution in whole animals.

Diverse algorithms can integrate transcriptomics with genome-scale metabolic models (GEMs) to build context-specific metabolic models. These algorithms require identification of a list of high confidence (core) reactions from transcriptomics, but parameters related to identification of core reactions, such as thresholding of expression profiles, can significantly change model content. Importantly, current thresholding approaches are burdened with setting singular arbitrary thresholds for all genes; thus, resulting in removal of enzymes needed in small amounts and even many housekeeping genes. Here, we describe StanDep, a novel heuristic method for using transcriptomics to identify core reactions prior to building context-specific metabolic models. StanDep clusters gene expression data based on their expression pattern across different contexts and determines thresholds for each cluster using data-dependent statistics, specifically standard deviation and mean. To demonstrate the use of StanDep, we built hundreds of models for the NCI-60 cancer cell lines. These models successfully increased the inclusion of housekeeping reactions, which are often lost in models built using standard thresholding approaches. Further, StanDep also provided a transcriptomic explanation for inclusion of lowly expressed reactions that were otherwise only supported by model extraction methods. Our study also provides novel insights into how cells may deal with context-specific and ubiquitous functions. StanDep, as a MATLAB toolbox, is available at https://github.com/LewisLabUCSD/StanDep.

The concept of "housekeeping gene" has been used for four decades but remains loosely defined. Housekeeping genes are commonly described as "essential for cellular existence regardless of their specific function in the tissue or organism", and "stably expressed irrespective of tissue type, developmental stage, cell cycle state, or external signal". However, experimental support for the tenet that gene essentiality is linked to stable expression across cell types, conditions, and organisms has been limited. Here we use genome-scale functional genomic screens together with bulk and single-cell sequencing technologies to test this link and optimize a quantitative and experimentally validated definition of housekeeping gene. Using the optimized definition, we identify, characterize, and provide as resources, housekeeping gene lists extracted from several human datasets, and 10 other animal species that include primates, chicken, and C. elegans. We find that stably expressed genes are not necessarily essential, and that the individual genes that are essential and stably expressed can considerably differ across organisms; yet the pathways enriched among these genes are conserved. Further, the level of conservation of housekeeping genes across the analyzed organisms captures their taxonomic groups, showing evolutionary relevance for our definition. Therefore, we present a quantitative and experimentally supported definition of housekeeping genes that can contribute to better understanding of their unique biological and evolutionary characteristics.

Plasticity in multicellular organisms involves signaling pathways converting contexts-either natural environmental challenges or laboratory perturbations-into context-specific changes in gene expression. Congruently, the interactions between the signaling molecules and transcription factors (TF) regulating these responses are also context specific. However, when a target gene responds across contexts, the upstream TF identified in one context is often inferred to regulate it across contexts. Reconciling these stable TF-target gene pair inferences with the context-specific nature of homeostatic responses is therefore needed. The induction of the Caenorhabditis elegans genes lipl-3 and lipl-4 is observed in many genetic contexts and is essential to survival during fasting. We find DAF-16/FOXO mediating lipl-4 induction in all contexts tested; hence, lipl-4 regulation seems context independent and compatible with across-context inferences. In contrast, DAF-16-mediated regulation of lipl-3 is context specific. DAF-16 reduces the induction of lipl-3 during fasting, yet it promotes it during oxidative stress. Through discrete dynamic modeling and genetic epistasis, we define that DAF-16 represses HLH-30/TFEB-the main TF activating lipl-3 during fasting. Contrastingly, DAF-16 activates the stress-responsive TF HSF-1 during oxidative stress, which promotes C. elegans survival through induction of lipl-3 Furthermore, the TF MXL-3 contributes to the dominance of HSF-1 at the expense of HLH-30 during oxidative stress but not during fasting. This study shows how context-specific diverting of functional interactions within a molecular network allows cells to specifically respond to a large number of contexts with a limited number of molecular players, a mode of transcriptional regulation we name "contextualized transcription."