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.

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.

Human gene-centered transcription factor networks for enhancers and disease variants.

Gene regulatory networks (GRNs) comprising interactions between transcription factors (TFs) and regulatory loci control development and physiology. Numerous disease-associated mutations have been identified, the vast majority residing in non-coding regions of the genome. As current GRN mapping methods test one TF at a time and require the use of cells harboring the mutation(s) of interest, they are not suitable to identify TFs that bind to wild-type and mutant loci. Here, we use gene-centered yeast one-hybrid (eY1H) assays to interrogate binding of 1,086 human TFs to 246 enhancers, as well as to 109 non-coding disease mutations. We detect both loss and gain of TF interactions with mutant loci that are concordant with target gene expression changes. This work establishes eY1H assays as a powerful addition to the toolkit of mapping human GRNs and for the high-throughput characterization of genomic variants that are rapidly being identified by genome-wide association studies.

Widespread macromolecular interaction perturbations in human genetic disorders.

How disease-associated mutations impair protein activities in the context of biological networks remains mostly undetermined. Although a few renowned alleles are well characterized, functional information is missing for over 100,000 disease-associated variants. Here we functionally profile several thousand missense mutations across a spectrum of Mendelian disorders using various interaction assays. The majority of disease-associated alleles exhibit wild-type chaperone binding profiles, suggesting they preserve protein folding or stability. While common variants from healthy individuals rarely affect interactions, two-thirds of disease-associated alleles perturb protein-protein interactions, with half corresponding to "edgetic" alleles affecting only a subset of interactions while leaving most other interactions unperturbed. With transcription factors, many alleles that leave protein-protein interactions intact affect DNA binding. Different mutations in the same gene leading to different interaction profiles often result in distinct disease phenotypes. Thus disease-associated alleles that perturb distinct protein activities rather than grossly affecting folding and stability are relatively widespread.

Interspecies systems biology uncovers metabolites affecting C. elegans gene expression and life history traits.

Diet greatly influences gene expression and physiology. In mammals, elucidating the effects and mechanisms of individual nutrients is challenging due to the complexity of both the animal and its diet. Here, we used an interspecies systems biology approach with Caenorhabditis elegans and two of its bacterial diets, Escherichia coli and Comamonas aquatica, to identify metabolites that affect the animal's gene expression and physiology. We identify vitamin B12 as the major dilutable metabolite provided by Comamonas aq. that regulates gene expression, accelerates development, and reduces fertility but does not affect lifespan. We find that vitamin B12 has a dual role in the animal: it affects development and fertility via the methionine/S-Adenosylmethionine (SAM) cycle and breaks down the short-chain fatty acid propionic acid, preventing its toxic buildup. Our interspecies systems biology approach provides a paradigm for understanding complex interactions between diet and physiology.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

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.

Diet-induced developmental acceleration independent of TOR and insulin in C. elegans.

Dietary composition has major effects on physiology. Here, we show that developmental rate, reproduction, and lifespan are altered in C. elegans fed Comamonas DA1877 relative to those fed a standard E. coli OP50 diet. We identify a set of genes that change in expression in response to this diet and use the promoter of one of these (acdh-1) as a dietary sensor. Remarkably, the effects on transcription and development occur even when Comamonas DA1877 is diluted with another diet, suggesting that Comamonas DA1877 generates a signal that is sensed by the nematode. Surprisingly, the developmental effect is independent from TOR and insulin signaling. Rather, Comamonas DA1877 affects cyclic gene expression during molting, likely through the nuclear hormone receptor NHR-23. Altogether, our findings indicate that different bacteria elicit various responses via distinct mechanisms, which has implications for diseases such as obesity and the interactions between the human microbiome and intestinal cells.

Integration of metabolic and gene regulatory networks modulates the C. elegans dietary response.

Expression profiles are tailored according to dietary input. However, the networks that control dietary responses remain largely uncharacterized. Here, we combine forward and reverse genetic screens to delineate a network of 184 genes that affect the C. elegans dietary response to Comamonas DA1877 bacteria. We find that perturbation of a mitochondrial network composed of enzymes involved in amino acid metabolism and the TCA cycle affects the dietary response. In humans, mutations in the corresponding genes cause inborn diseases of amino acid metabolism, most of which are treated by dietary intervention. We identify several transcription factors (TFs) that mediate the changes in gene expression upon metabolic network perturbations. Altogether, our findings unveil a transcriptional response system that is poised to sense dietary cues and metabolic imbalances, illustrating extensive communication between metabolic networks in the mitochondria and gene regulatory networks in the nucleus.

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 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 multiparameter network reveals extensive divergence between C. elegans bHLH transcription factors.

Differences in expression, protein interactions, and DNA binding of paralogous transcription factors ("TF parameters") are thought to be important determinants of regulatory and biological specificity. However, both the extent of TF divergence and the relative contribution of individual TF parameters remain undetermined. We comprehensively identify dimerization partners, spatiotemporal expression patterns, and DNA-binding specificities for the C. elegans bHLH family of TFs, and model these data into an integrated network. This network displays both specificity and promiscuity, as some bHLH proteins, DNA sequences, and tissues are highly connected, whereas others are not. By comparing all bHLH TFs, we find extensive divergence and that all three parameters contribute equally to bHLH divergence. Our approach provides a framework for examining divergence for other protein families in C. elegans and in other complex multicellular organisms, including humans. Cross-species comparisons of integrated networks may provide further insights into molecular features underlying protein family evolution. For a video summary of this article, see the PaperFlick file available with the online Supplemental Data.

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.

Understanding Metabolic Regulation at a Systems Level: Metabolite Sensing, Mathematical Predictions, and Model Organisms.

Metabolic networks are extensively regulated to facilitate tissue-specific metabolic programs and robustly maintain homeostasis in response to dietary changes. Homeostatic metabolic regulation is achieved through metabolite sensing coupled to feedback regulation of metabolic enzyme activity or expression. With a wealth of transcriptomic, proteomic, and metabolomic data available for different cell types across various conditions, we are challenged with understanding global metabolic network regulation and the resulting metabolic outputs. Stoichiometric metabolic network modeling integrated with "omics" data has addressed this challenge by generating nonintuitive, testable hypotheses about metabolic flux rewiring. Model organism studies have also yielded novel insight into metabolic networks. This review covers three topics: the feedback loops inherent in metabolic regulatory networks, metabolic network modeling, and interspecies studies utilizing Caenorhabditis elegans and various bacterial diets that have revealed novel metabolic paradigms.

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.

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.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network.

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Functional Conservation of a Developmental Switch in Mammals since the Jurassic Age.

ThPOK is a "master regulator" of T lymphocyte lineage choice, whose presence or absence is sufficient to dictate development to the CD4 or CD8 lineages, respectively. Induction of ThPOK is transcriptionally regulated, via a lineage-specific silencer element, SilThPOK. Here, we take advantage of the available genome sequence data as well as site-specific gene targeting technology, to evaluate the functional conservation of ThPOK regulation across mammalian evolution, and assess the importance of motif grammar (order and orientation of TF binding sites) on SilThPOK function in vivo. We make three important points: First, the SilThPOK is present in marsupial and placental mammals, but is not found in available genome assemblies of nonmammalian vertebrates, indicating that it arose after divergence of mammals from other vertebrates. Secondly, by replacing the murine SilThPOK in situ with its marsupial equivalent using a knockin approach, we demonstrate that the marsupial SilThPOK supports correct CD4 T lymphocyte lineage-specification in mice. To our knowledge, this is the first in vivo demonstration of functional equivalency for a silencer element between marsupial and placental mammals using a definitive knockin approach. Finally, we show that alteration of the position/orientation of a highly conserved region within the murine SilThPOK is sufficient to destroy silencer activity in vivo, demonstrating that motif grammar of this "solid" synteny block is critical for silencer function. Dependence of SilThPOK function on motif grammar conserved since the mid-Jurassic age, 165 Ma, suggests that the SilThPOK operates as a silenceosome, by analogy with the previously proposed enhanceosome model.

Many transcription factors contribute to C. elegans growth and fat storage.

Reverse genetic screens by RNA interference (RNAi) in model organisms such as the nematode Caenorhabditis elegans have provided numerous insights into gene function, thereby connecting genotype to phenotype. However, genes that contribute only subtly are often missed because relatively large numbers of measurements and reliable quantification are required to overcome experimental and biological noise that may mask subtle phenotypic effects. Here, we address this challenge by focusing on two phenotypes in C. elegans: growth and fat storage. We carried out comprehensive RNAi knockdown of transcription factors (TFs), as these are known important regulators of biological processes during development and the maintenance of homeostasis. Microscopy images of TF knockdown animals stained with Oil Red O (ORO) were captured, and body size (proxy for growth) and ORO staining intensity (proxy for fat storage) were precisely quantified using a newly developed imaging tool we named IPPOME (Image Processing for Precise and Objective MEasurement). We found that a surprisingly large proportion of TFs contribute to growth and fat storage, but that most TFs have only subtle, yet significant effects. This study provides a blueprint for studies of other genes and phenotypes in C. elegans.