A compendium of Caenorhabditis elegans regulatory transcription factors: a resource for mapping transcription regulatory networks

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

A full-body transcription factor expression atlas with completely resolved cell identities in C. elegans

Invariant cell lineage in C. elegans enables spatiotemporal resolution of transcriptional regulatory mechanisms controlling the fate of each cell. Here, we develop RAPCAT (Robust-point-matching- And Piecewise-affine-based Cell Annotation Tool) to automate cell identity assignment in three-dimensional image stacks of L1 larvae and profile reporter expression of 620 transcription factors in every cell. Transcription factor profile-based clustering analysis defines 80 cell types distinct from conventional phenotypic cell types and identifies three general phenotypic modalities related to these classifications. First, transcription factors are broadly downregulated in quiescent stage Hermaphrodite Specific Neurons, suggesting stage- and cell type-specific variation in transcriptome size. Second, transcription factor expression is more closely associated with morphology than other phenotypic modalities in different pre- and post-differentiation developmental stages. Finally, embryonic cell lineages can be associated with specific transcription factor expression patterns and functions that persist throughout postembryonic life. This study presents a comprehensive transcription factor atlas for investigation of intra-cell type heterogeneity.

Mechanisms of lineage specification in Caenorhabditis elegans

The studies of cell fate and lineage specification are fundamental to our understanding of the development of multicellular organisms. Caenorhabditis elegans has been one of the premiere systems for studying cell fate specification mechanisms at single cell resolution, due to its transparent nature, the invariant cell lineage, and fixed number of somatic cells. We discuss the general themes and regulatory mechanisms that have emerged from these studies, with a focus on somatic lineages and cell fates. We next review the key factors and pathways that regulate the specification of discrete cells and lineages during embryogenesis and postembryonic development; we focus on transcription factors and include numerous lineage diagrams that depict the expression of key factors that specify embryonic founder cells and postembryonic blast cells, and the diverse somatic cell fates they generate. We end by discussing some future perspectives in cell and lineage specification.

Dynamic compartmentalization of the pro-invasive transcription factor NHR-67 reveals a role for Groucho in regulating a proliferative-invasive cellular switch in C. elegans

A growing body of evidence suggests that cell division and basement membrane invasion are mutually exclusive cellular behaviors. How cells switch between proliferative and invasive states is not well understood. Here, we investigated this dichotomy in vivo by examining two cell types in the developing Caenorhabditis elegans somatic gonad that derive from equipotent progenitors, but exhibit distinct cell behaviors: the post-mitotic, invasive anchor cell and the neighboring proliferative, non-invasive ventral uterine (VU) cells. We show that the fates of these cells post-specification are more plastic than previously appreciated and that levels of NHR-67 are important for discriminating between invasive and proliferative behavior. Transcription of NHR-67 is downregulated following post-translational degradation of its direct upstream regulator, HLH-2 (E/Daughterless) in VU cells. In the nuclei of VU cells, residual NHR-67 protein is compartmentalized into discrete punctae that are dynamic over the cell cycle and exhibit liquid-like properties. By screening for proteins that colocalize with NHR-67 punctae, we identified new regulators of uterine cell fate maintenance: homologs of the transcriptional co-repressor Groucho (UNC-37 and LSY-22), as well as the TCF/LEF homolog POP-1. We propose a model in which the association of NHR-67 with the Groucho/TCF complex suppresses the default invasive state in non-invasive cells, which complements transcriptional regulation to add robustness to the proliferative-invasive cellular switch in vivo.

MYCT1 in cancer development: Gene structure, regulation, and biological implications for diagnosis and treatment

Myc target 1 (MYCT1), located at 6q25.2, is a crucial player in cancer development. While widely distributed in cells, its subcellular localization varies across different cancer types. As a novel c-Myc target gene, MYCT1 is subject to regulation by multiple transcription factors. Studies have revealed aberrant expression of MYCT1 in various cancers, impacting pivotal biological processes such as proliferation, apoptosis, migration, genomic instability, and differentiation in cancer cells. Additionally, MYCT1 plays a critical role in modulating tumor angiogenesis and remodeling tumor immune responses within the tumor microenvironment. Despite certain debated functions, MYCT1 undeniably holds significance in cancer development. In this review, we comprehensively examine the relationship between MYCT1 and cancer, encompassing gene structure, regulation of gene expression, gene mutation, and biological function, with the aim of providing valuable insights for cancer diagnosis and treatment.

Whole-body gene expression atlas of an adult metazoan

Gene activity defines cell identity, drives intercellular communication, and underlies the functioning of multicellular organisms. We present the single-cell resolution atlas of gene activity of a fertile adult metazoan: Caenorhabditis elegans. This compendium comprises 180 distinct cell types and 19,657 expressed genes. We predict 7541 transcription factor expression profile associations likely responsible for defining cellular identity. We predict thousands of intercellular interactions across the C. elegans body and the ligand-receptor pairs that mediate them, some of which we experimentally validate. We identify 172 genes that show consistent expression across cell types, are involved in basic and essential functions, and are conserved across phyla; therefore, we present them as experimentally validated housekeeping genes. We developed the WormSeq application to explore these data. In addition to the integrated gene-to-systems biology, we present genome-scale single-cell resolution testable hypotheses that we anticipate will advance our understanding of the molecular mechanisms, underlying the functioning of a multicellular organism and the perturbations that lead to its malfunction.

Transcription factors regulating the fate and developmental potential of a multipotent progenitor in Caenorhabditis elegans

Multipotent stem and progenitor cells have the capacity to generate a limited array of related cell types. The Caenorhabditis elegans somatic gonadal precursors are multipotent progenitors that generate all 143 cells of the somatic gonad, including complex tissues and specialized signaling cells. To screen for candidate regulators of cell fate and multipotency, we identified transcription factor genes with higher expression in somatic gonadal precursors than in their differentiated sister, the head mesodermal cell. We used RNA interference or genetic mutants to reduce the function of 183 of these genes and examined the worms for defects in the somatic gonadal precursor cell fate or the ability to generate gonadal tissue types. We identify 8 genes that regulate somatic gonadal precursor fate, including the SWI/SNF chromatin remodeling complex gene swsn-3 and the Ci/GLI homolog tra-1, which is the terminal regulator of sex determination. Four genes are necessary for somatic gonadal precursors to generate the correct number and type of descendant cells. We show that the E2F homolog, efl-3, regulates the cell fate decision between distal tip cells and the sheath/spermathecal precursor. We find that the FACT complex gene hmg-4 is required for the generation of the correct number of somatic gonadal precursor descendants, and we define an earlier role for the nhr-25 nuclear hormone receptor-encoding gene, in addition to its previously described role in regulating the asymmetric division of somatic gonadal precursors. Overall, our data show that genes regulating cell fate are largely different from genes regulating developmental potential, demonstrating that these processes are genetically separable.

A developmental pathway for epithelial-to-motoneuron transformation in C. elegans

Motoneurons and motoneuron-like pancreatic β cells arise from radial glia and ductal cells, respectively, both tube-lining progenitors that share molecular regulators. To uncover programs underlying motoneuron formation, we studied a similar, cell-division-independent transformation of the C. elegans tube-lining Y cell into the PDA motoneuron. We find that lin-12/Notch acts through ngn-1/Ngn and its regulator hlh-16/Olig to control transformation timing. lin-12 loss blocks transformation, while lin-12(gf) promotes precocious PDA formation. Early basal expression of ngn-1/Ngn and hlh-16/Olig depends on sem-4/Sall and egl-5/Hox. Later, coincident with Y cell morphological changes, ngn-1/Ngn expression is upregulated in a sem-4/Sall and egl-5/Hox-dependent but hlh-16/Olig-independent manner. Subsequently, Y cell retrograde extension forms an anchored process priming PDA axon extension. Extension requires ngn-1-dependent expression of the cytoskeleton organizers UNC-119, UNC-44/ANK, and UNC-33/CRMP, which also activate PDA terminal-gene expression. Our findings uncover cell-division-independent regulatory events leading to motoneuron generation, suggesting a conserved pathway for epithelial-to-motoneuron/motoneuron-like cell differentiation.

Rashid et al. report on a conserved epithelial-to-motoneuron transformation pathway in C. elegans requiring ngn-1/Ngn and hlh-16/Olig. lin-12/Notch regulates transformation timing through these genes, while ngn-1/Ngn and hlh-16/Olig expression levels are regulated by sem-4/Sall and egl-5/Hox. Unexpectedly, the cytoskeleton organizers UNC-119, UNC-44, and UNC-33, which are ngn-1/Ngn targets, promote motoneuron terminal identity.

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.

Here, Diot et al. use the nematode Caenorhabditis elegans as a model to identify off-target toxicity mechanisms for tamoxifen, and find that these include fatty acid metabolism and cell death, which can be modulated by different bacterial species.

Division and Transmission: Malaria Parasite Development in the Mosquito

The malaria parasite life cycle alternates between two hosts: a vertebrate and the female Anopheles mosquito vector. Cell division, proliferation, and invasion are essential for parasite development, transmission, and survival. Most research has focused on Plasmodium development in the vertebrate, which causes disease; however, knowledge of malaria parasite development in the mosquito (the sexual and transmission stages) is now rapidly accumulating, gathered largely through investigation of the rodent malaria model, with Plasmodium berghei. In this review, we discuss the seminal genome-wide screens that have uncovered key regulators of cell proliferation, invasion, and transmission during Plasmodium sexual development. Our focus is on the roles of transcription factors, reversible protein phosphorylation, and molecular motors. We also emphasize the still-unanswered important questions around key pathways in cell division during the vector transmission stages and how they may be targeted in future studies.

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.

Transcription Factors That Control Behavior-Lessons From C. elegans

Behavior encompasses the physical and chemical response to external and internal stimuli. Neurons, each with their own specific molecular identities, act in concert to perceive and relay these stimuli to drive behavior. Generating behavioral responses requires neurons that have the correct morphological, synaptic, and molecular identities. Transcription factors drive the specific gene expression patterns that define these identities, controlling almost every phenomenon in a cell from development to homeostasis. Therefore, transcription factors play an important role in generating and regulating behavior. Here, we describe the transcription factors, the pathways they regulate, and the neurons that drive chemosensation, mechanosensation, thermosensation, osmolarity sensing, complex, and sex-specific behaviors in the animal model Caenorhabditis elegans. We also discuss the current limitations in our knowledge, particularly our minimal understanding of how transcription factors contribute to the adaptive behavioral responses that are necessary for organismal survival.

A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis

Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state-fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared with insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and extended haplotype homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency-derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbor a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.

Internetwork connectivity of molecular networks across species of life

Molecular interactions are studied as independent networks in systems biology. However, molecular networks do not exist independently of each other. In a network of networks approach (called multiplex), we study the joint organization of transcriptional regulatory network (TRN) and protein-protein interaction (PPI) network. We find that TRN and PPI are non-randomly coupled across five different eukaryotic species. Gene degrees in TRN (number of downstream genes) are positively correlated with protein degrees in PPI (number of interacting protein partners). Gene-gene and protein-protein interactions in TRN and PPI, respectively, also non-randomly overlap. These design principles are conserved across the five eukaryotic species. Robustness of the TRN-PPI multiplex is dependent on this coupling. Functionally important genes and proteins, such as essential, disease-related and those interacting with pathogen proteins, are preferentially situated in important parts of the human multiplex with highly overlapping interactions. We unveil the multiplex architecture of TRN and PPI. Multiplex architecture may thus define a general framework for studying molecular networks. This approach may uncover the building blocks of the hierarchical organization of molecular interactions.

Combining Auxin-Induced Degradation and RNAi Screening Identifies Novel Genes Involved in Lipid Bilayer Stress Sensing in Caenorhabditis elegans

Alteration of the lipid composition of biological membranes interferes with their function and can cause tissue damage by triggering apoptosis. Upon lipid bilayer stress, the endoplasmic reticulum mounts a stress response similar to the unfolded protein response. However, only a few genes are known to regulate lipid bilayer stress. We performed a suppressor screen that combined the auxin-inducible degradation (AID) system with conventional RNAi in C. elegans to identify members of the lipid bilayer stress response. AID-mediated degradation of the mediator MDT-15, a protein required for the upregulation of fatty acid desaturases, induced the activation of lipid bilayer stress-sensitive reporters. We screened through most C. elegans kinases and transcription factors by feeding RNAi. We discovered nine genes that suppressed the lipid bilayer stress response in C. elegans. These suppressor genes included drl-1/MAP3K3, gsk-3/GSK3, let-607/CREB3, ire-1/IRE1, and skn-1/NRF1,2,3. Our candidate suppressor genes suggest a network of transcription factors and the integration of multiple tissues for a centralized lipotoxicity response in the intestine. Thus, we demonstrated proof-of-concept for combining AID and RNAi as a new screening strategy and identified eight conserved genes that had not previously been implicated in the lipid bilayer stress response.

A Genetic Model of the Connectome

The connectomes of organisms of the same species show remarkable architectural and often local wiring similarity, raising the question: where and how is neuronal connectivity encoded? Here, we start from the hypothesis that the genetic identity of neurons guides synapse and gap-junction formation and show that such genetically driven wiring predicts the existence of specific biclique motifs in the connectome. We identify a family of large, statistically significant biclique subgraphs in the connectomes of three species and show that within many of the observed bicliques the neurons share statistically significant expression patterns and morphological characteristics, supporting our expectation of common genetic factors that drive the synapse formation within these subgraphs. The proposed connectome model offers a self-consistent framework to link the genetics of an organism to the reproducible architecture of its connectome, offering experimentally falsifiable predictions on the genetic factors that drive the formation of individual neuronal circuits.

Activation of the nucleotide excision repair pathway by crude oil exposure: A translational study from model organisms to the Hebei Spirit Oil Spill Cohort

In order to gain insight into the human health implications of the Hebei Spirit Oil Spill (HSOS), the mechanism of toxicity of the Iranian heavy crude (IHC), the main oil component in the HSOS was investigated in Caenorhabditis elegans and zebrafish. The identified mechanism was translated to humans using blood samples from Taean residents, who experienced HSOS with different levels of exposure to the spill. C. elegans TF RNAi screening with IHC oil revealed the nucleotide excision repair (NER) pathway as being significantly involved by oil exposure. To identify the main toxicity contributors within the chemical mixture of the crude oil, further studies were conducted on C. elegans by exposure to C3-naphthalene, an alkylated polycyclic aromatic hydrocarbon (PAH), which constitutes one of the major components of IHC oil. Increased expression of NER pathway genes was observed following exposure to the IHC oil, C3-naphthalene enriched fraction and C3-naphthalene. As the NER pathway is conserved in fish and humans, the same experiment was conducted in zebrafish, and the data were similar to what was seen in C. elegans. Increased expression of NER pathway genes was observed in human samples from the high exposure group, which suggests the involvement of the NER pathway in IHC oil exposure. Overall, the study suggests that IHC oil may cause bulk damage to DNA and activation of the NER system and Alkylated PAHs are the major contributor to DNA damage. Our study provides an innovative approach for studying translational toxicity testing from model organisms to human health.

Early life exposure of a biocide, CMIT/MIT causes metabolic toxicity via the O-GlcNAc transferase pathway in the nematode C. elegans

Unusual cases of fatal lung injury, later determined to be a result of exposure to chemicals used as humidifier disinfectants, were reported among Korean children from 2006 to 2011. This resulted in considerable study of the pulmonary toxicity of humidifier disinfectant chemicals to establish the causal relationship between exposure and lung disease. However, the systemic toxicity of the former and health effects other than lung disease are not fully understood. Here, we investigated the effect of 5-chloro-2-methyl-4-isothiazoline-3-one and 2-methyl-4-isothiazolin-3-one (CMIT/MIT), among the humidifier disinfectants used in the accidents, on the development of metabolic toxicity in the model organism, Caenorhabditis elegans using an exposure scenario comparison. We screened the potential of CMIT/MIT to induce metabolic toxicity using C. elegans oga-1(ok1207) and ogt-1(ok1474) mutants. We also performed a pathway analysis based on C. elegans transcription factor RNAi library screening to identify the underlying toxicity mechanisms. Finally, to understand the critical window of exposure for metabolic toxicity, responses to exposure during different periods in the life cycles of the worms were compared. We determined that CMIT/MIT could induce metabolic toxicity through O-linked N-acetylglucosamine transferase and early life seems to be the critical window for exposure for metabolic toxicity for this substance. The O-linked N-acetylglucosamine transferase pathway is conserved from worms to humans; our results thus insinuate that early-life exposure to CMIT/MIT could cause metabolic health problems during adult life in humans. We therefore suggest that a systemic toxicity approach should be considered to comprehensively understand the adverse health effects of humidifier disinfectant misuse.

Caenorhabditis elegans hub genes that respond to amyloid beta are homologs of genes involved in human Alzheimer's disease

The prominent characteristic of Alzheimer’s disease (AD) is the accumulation of amyloid beta (Abeta) proteins in the form of plaques that cause molecular and cellular alterations in the brain. Due to the paucity of brain samples of early-stage Abeta aggregation, animal models have been developed to study early events in AD. Caenorhabditis elegans is a genetically tractable animal model for AD. Here, we used transcriptomic data, network-based protein-protein interactions and weighted gene co-expression network analysis (WGCNA), to detect modules and their gene ontology in response to Abeta aggregation in C. elegans. Additionally, hub genes and their orthologues in human and mouse were identified to study their relation to AD. We also found several transcription factors (TFs) responding to Abeta accumulation. Our results show that Abeta expression in C. elegans relates to general processes such as molting cycle, locomotion, and larval development plus AD-associated processes, including protein phosphorylation, and G-protein coupled receptor-regulated pathways. We reveal that many hub genes and TFs including ttbk-2, daf-16, and unc-49 have human and mouse orthologues that are directly or potentially associated with AD and neural development. In conclusion, using systems biology we identified important genes and biological processes in C. elegans that respond to Abeta aggregation, which could be used as potential diagnostic or therapeutic targets. In addition, because of evolutionary relationship to AD in human, we suggest that C. elegans is a useful model for studying early molecular events in AD.

mRNA profiling reveals significant transcriptional differences between a multipotent progenitor and its differentiated sister

BACKGROUND

The two Caenorhabditis elegans somatic gonadal precursors (SGPs) are multipotent progenitors that generate all somatic tissues of the adult reproductive system. The sister cells of the SGPs are two head mesodermal cells (hmcs); one hmc dies by programmed cell death and the other terminally differentiates. Thus, a single cell division gives rise to one multipotent progenitor and one differentiated cell with identical lineage histories. We compared the transcriptomes of SGPs and hmcs in order to learn the determinants of multipotency and differentiation in this lineage.

RESULTS

We generated a strain that expressed fluorescent markers specifically in SGPs (ehn-3A::tdTomato) and hmcs (bgal-1::GFP). We dissociated cells from animals after the SGP/hmc cell division, but before the SGPs had further divided, and subjected the dissociated cells to fluorescence-activated cell sorting to collect isolated SGPs and hmcs. We analyzed the transcriptomes of these cells and found that 5912 transcripts were significantly differentially expressed, with at least two-fold change in expression, between the two cell types. The hmc-biased genes were enriched with those that are characteristic of neurons. The SGP-biased genes were enriched with those indicative of cell proliferation and development. We assessed the validity of our differentially expressed genes by examining existing reporters for five of the 10 genes with the most significantly biased expression in SGPs and found that two showed expression in SGPs. For one reporter that did not show expression in SGPs, we generated a GFP knock-in using CRISPR/Cas9. This reporter, in the native genomic context, was expressed in SGPs.

CONCLUSIONS

We found that the transcriptional profiles of SGPs and hmcs are strikingly different. The hmc-biased genes are enriched with those that encode synaptic transmission machinery, which strongly suggests that it has neuron-like signaling properties. In contrast, the SGP-biased genes are enriched with genes that encode factors involved in transcription and translation, as would be expected from a cell preparing to undergo proliferative divisions. Mediators of multipotency are likely to be among the genes differentially expressed in SGPs.

Options and Considerations When Using a Yeast One-Hybrid System

Comprehensive mapping of protein-DNA interactions is essential to uncover the mechanisms involved in gene regulation. However, the data generated has been sparse given the number of regulatory elements and transcription factors (TFs) encoded in the genomes of metazoan organisms. Yeast one-hybrid (Y1H) assays provide a powerful "DNA-centered" method, complementary to "TF-centered" methods such as chromatin immunoprecipitation, to identify the TFs that can bind a DNA sequence of interest. Here, we present different technical variations that should be considered when using a Y1H system, including the type of DNA sequence to test, source of TF clones, as well as types of vectors and screening format. Finally, we discuss limitations of the assay and future challenges.

Transcriptome analysis of adult Caenorhabditis elegans cells reveals tissue-specific gene and isoform expression

The biology and behavior of adults differ substantially from those of developing animals, and cell-specific information is critical for deciphering the biology of multicellular animals. Thus, adult tissue-specific transcriptomic data are critical for understanding molecular mechanisms that control their phenotypes. We used adult cell-specific isolation to identify the transcriptomes of C. elegans' four major tissues (or "tissue-ome"), identifying ubiquitously expressed and tissue-specific "enriched" genes. These data newly reveal the hypodermis' metabolic character, suggest potential worm-human tissue orthologies, and identify tissue-specific changes in the Insulin/IGF-1 signaling pathway. Tissue-specific alternative splicing analysis identified a large set of collagen isoforms. Finally, we developed a machine learning-based prediction tool for 76 sub-tissue cell types, which we used to predict cellular expression differences in IIS/FOXO signaling, stage-specific TGF-β activity, and basal vs. memory-induced CREB transcription. Together, these data provide a rich resource for understanding the biology governing multicellular adult animals.

Gateway-Compatible Yeast One-Hybrid and Two-Hybrid Assays

In the first section of this introduction, we provide background information for yeast two-hybrid (Y2H) assays that provide a genetic method for the identification and analysis of binary protein-protein interactions and that are complementary to biochemical methods such as immunoprecipitation. In the second section, we discuss yeast one-hybrid (Y1H) assays that provide a "gene-centered" (DNA-to-protein) genetic method to identify and study protein-DNA interactions between cis-regulatory elements and transcription factors (TFs). This method is complementary to "TF-centered" (protein-to-DNA) biochemical methods such as chromatin immunoprecipitation.

Evolutionary plasticity in the innate immune function of Akirin

Eukaryotic gene expression requires the coordinated action of transcription factors, chromatin remodelling complexes and RNA polymerase. The conserved nuclear protein Akirin plays a central role in immune gene expression in insects and mammals, linking the SWI/SNF chromatin-remodelling complex with the transcription factor NFκB. Although nematodes lack NFκB, Akirin is also indispensable for the expression of defence genes in the epidermis of Caenorhabditis elegans following natural fungal infection. Through a combination of reverse genetics and biochemistry, we discovered that in C. elegans Akirin has conserved its role of bridging chromatin-remodellers and transcription factors, but that the identity of its functional partners is different since it forms a physical complex with NuRD proteins and the POU-class transcription factor CEH-18. In addition to providing a substantial step forward in our understanding of innate immune gene regulation in C. elegans, our results give insight into the molecular evolution of lineage-specific signalling pathways.

The ModERN Resource: Genome-Wide Binding Profiles for Hundreds of Drosophila and Caenorhabditis elegans Transcription Factors

To develop a catalog of regulatory sites in two major model organisms, Drosophila melanogaster and Caenorhabditis elegans, the modERN (model organism Encyclopedia of Regulatory Networks) consortium has systematically assayed the binding sites of transcription factors (TFs). Combined with data produced by our predecessor, modENCODE (Model Organism ENCyclopedia Of DNA Elements), we now have data for 262 TFs identifying 1.23 M sites in the fly genome and 217 TFs identifying 0.67 M sites in the worm genome. Because sites from different TFs are often overlapping and tightly clustered, they fall into 91,011 and 59,150 regions in the fly and worm, respectively, and these binding sites span as little as 8.7 and 5.8 Mb in the two organisms. Clusters with large numbers of sites (so-called high occupancy target, or HOT regions) predominantly associate with broadly expressed genes, whereas clusters containing sites from just a few factors are associated with genes expressed in tissue-specific patterns. All of the strains expressing GFP-tagged TFs are available at the stock centers, and the chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center and also through a simple interface (http://epic.gs.washington.edu/modERN/) that facilitates rapid accessibility of processed data sets. These data will facilitate a vast number of scientific inquiries into the function of individual TFs in key developmental, metabolic, and defense and homeostatic regulatory pathways, as well as provide a broader perspective on how individual TFs work together in local networks and globally across the life spans of these two key model organisms.

A microRNA feedback loop regulates global microRNA abundance during aging

Expression levels of many microRNAs (miRNAs) change during aging, notably declining globally in a number of organisms and tissues across taxa. However, little is known about the mechanisms or the biological relevance for this change. We investigated the network of genes that controls miRNA transcription and processing during C. elegans aging. We found that miRNA biogenesis genes are highly networked with transcription factors and aging-associated miRNAs. In particular, miR-71, known to influence life span and itself up-regulated during aging, represses alg-1/Argonaute expression post-transcriptionally during aging. Increased ALG-1 abundance in mir-71 loss-of-function mutants led to globally increased miRNA expression. Interestingly, these mutants demonstrated widespread mRNA expression dysregulation and diminished levels of variability both in gene expression and in overall life span. Thus, the progressive molecular decline often thought to be the result of accumulated damage over an organism's life may be partially explained by a miRNA-directed mechanism of age-associated decline.

Systems biology of embryonic development: Prospects for a complete understanding of the Caenorhabditis elegans embryo

The convergence of developmental biology and modern genomics tools brings the potential for a comprehensive understanding of developmental systems. This is especially true for the Caenorhabditis elegans embryo because its small size, invariant developmental lineage, and powerful genetic and genomic tools provide the prospect of a cellular resolution understanding of messenger RNA (mRNA) expression and regulation across the organism. We describe here how a systems biology framework might allow large-scale determination of the embryonic regulatory relationships encoded in the C. elegans genome. This framework consists of two broad steps: (a) defining the "parts list"-all genes expressed in all cells at each time during development and (b) iterative steps of computational modeling and refinement of these models by experimental perturbation. Substantial progress has been made towards defining the parts list through imaging methods such as large-scale green fluorescent protein (GFP) reporter analysis. Imaging results are now being augmented by high-resolution transcriptome methods such as single-cell RNA sequencing, and it is likely the complete expression patterns of all genes across the embryo will be known within the next few years. In contrast, the modeling and perturbation experiments performed so far have focused largely on individual cell types or genes, and improved methods will be needed to expand them to the full genome and organism. This emerging comprehensive map of embryonic expression and regulatory function will provide a powerful resource for developmental biologists, and would also allow scientists to ask questions not accessible without a comprehensive picture. This article is categorized under: Invertebrate Organogenesis > Worms Technologies > Analysis of the Transcriptome Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

Construction of Arabidopsis Transcription Factor ORFeome Collections and Identification of Protein-DNA Interactions by High-Throughput Yeast One-Hybrid Screens

Identification of transcription factor (TF)-promoter interactions is key to understanding the basic molecular underpinnings of gene regulation. The complexity of gene regulation, however, is driven by the combined function of several TFs recruited to the promoter region, which often confounds the discovery of transcriptional regulatory mechanisms. Genome sequencing enabled the construction of TF-specific ORFeome clone collections that can be used to study TF function with unprecedented coverage. Among the recently developed methods, gene-centered yeast one-hybrid (Y1H) screens performed with these ORFeome collections provide a simple and reliable strategy to identify TF-promoter interactions. Here, we describe high-throughput cloning protocols used to generate a gold standard TF ORFeome collection for the model organism Arabidopsis thaliana. Furthermore, we outline the protocol to build a daughter clone collection suitable for the Y1H assay and a high-throughput Y1H screening procedure that enables rapid assessment of thousands TF-promoter interactions using a robotic platform. These protocols can be universally adopted to build ORFeome libraries and thus expand the usage of gene-centered Y1H screens or other alternative strategies for discovery and characterization of TF functions.

The Caenorhabditis elegans Female-Like State: Decoupling the Transcriptomic Effects of Aging and Sperm Status

Understanding genome and gene function in a whole organism requires us to fully comprehend the life cycle and the physiology of the organism in question. Caenorhabditis elegans XX animals are hermaphrodites that exhaust their sperm after 3 d of egg-laying. Even though C. elegans can live for many days after cessation of egg-laying, the molecular physiology of this state has not been as intensely studied as other parts of the life cycle, despite documented changes in behavior and metabolism. To study the effects of sperm depletion and aging of C. elegans during the first 6 d of adulthood, we measured the transcriptomes of first-day adult hermaphrodites and sixth-day sperm-depleted adults, and, at the same time points, mutant fog-2(lf) worms that have a feminized germline phenotype. We found that we could separate the effects of biological aging from sperm depletion. For a large subset of genes, young adult fog-2(lf) animals had the same gene expression changes as sperm-depleted sixth-day wild-type hermaphrodites, and these genes did not change expression when fog-2(lf) females reached the sixth day of adulthood. Taken together, this indicates that changing sperm status causes a change in the internal state of the worm, which we call the female-like state. Our data provide a high-quality picture of the changes that happen in global gene expression throughout the period of early aging in the worm.

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.

The Oxidative Stress Response in Caenorhabditis elegans Requires the GATA Transcription Factor ELT-3 and SKN-1/Nrf2

Cellular damage caused by reactive oxygen species is believed to be a major contributor to age-associated diseases. Previously, we characterized the Caenorhabditis elegans Brap2 ortholog (BRAP-2) and found that it is required to prevent larval arrest in response to elevated levels of oxidative stress. Here, we report that C. elegans brap-2 mutants display increased expression of SKN-1-dependent, phase II detoxification enzymes that is dependent on PMK-1 (a p38 MAPK C. elegans ortholog). An RNA-interference screen was conducted using a transcription factor library to identify genes required for increased expression of the SKN-1 target gst-4 in brap-2 mutants. We identified ELT-3, a member of the GATA transcription factor family, as a positive regulator of gst-4p::gfp expression. We found that ELT-3 interacts with SKN-1 to activate gst-4 transcription in vitro and that elt-3 is required for enhanced gst-4 expression in the brap-2(ok1492) mutant in vivo Furthermore, nematodes overexpressing SKN-1 required ELT-3 for life-span extension. Taken together, these results suggest a model where BRAP-2 acts as negative regulator of SKN-1 through inhibition of p38 MAPK activity, and that the GATA transcription factor ELT-3 is required along with SKN-1 for the phase II detoxification response in C. elegans.

Identification of an RNA Polymerase III Regulator Linked to Disease-Associated Protein Aggregation

Protein aggregation is associated with age-related neurodegenerative disorders, such as Alzheimer's and polyglutamine diseases. As a causal relationship between protein aggregation and neurodegeneration remains elusive, understanding the cellular mechanisms regulating protein aggregation will help develop future treatments. To identify such mechanisms, we conducted a forward genetic screen in a C. elegans model of polyglutamine aggregation and identified the protein MOAG-2/LIR-3 as a driver of protein aggregation. In the absence of polyglutamine, MOAG-2/LIR-3 regulates the RNA polymerase III-associated transcription of small non-coding RNAs. This regulation is lost in the presence of polyglutamine, which mislocalizes MOAG-2/LIR-3 from the nucleus to the cytosol. We then show biochemically that MOAG-2/LIR-3 can also catalyze the aggregation of polyglutamine-expanded huntingtin. These results suggest that polyglutamine can induce an aggregation-promoting activity of MOAG-2/LIR-3 in the cytosol. The concept that certain aggregation-prone proteins can convert other endogenous proteins into drivers of aggregation and toxicity adds to the understanding of how cellular homeostasis can be deteriorated in protein misfolding diseases.

Starvation-Induced Stress Response Is Critically Impacted by Ceramide Levels in Caenorhabditis elegans

Our understanding of the cellular mechanisms by which animals regulate their response to starvation is limited, despite the strong relevance of the problem to major human health issues. The L1 diapause of Caenorhabditis elegans, where first-stage larvae arrest in response to a food-less environment, is an excellent system to study this mechanism. We found, through genetic manipulation and lipid analysis, that biosynthesis of ceramide, particularly those with longer fatty acid side chains, critically impacts animal survival during L1 diapause. Genetic interaction analysis suggests that ceramide may act in both insulin-IGF-1 signaling (IIS)-dependent and IIS-independent pathways to affect starvation survival. Genetic and expression analyses indicate that ceramide is required for maintaining the proper expression of previously characterized starvation-responsive genes, genes that are regulated by the IIS pathway and tumor suppressor Rb, and genes responsive to pathogen. These findings provide an important insight into the roles of sphingolipid metabolism, not only in starvation response, but also in aging and food-response-related human health problems.

Gene-Centered Yeast One-Hybrid Assays

An important question when studying gene regulation is which transcription factors (TFs) interact with which cis-regulatory elements, such as promoters and enhancers. Addressing this issue in complex multicellular organisms is challenging as several hundreds of TFs and thousands of regulatory elements must be considered in the context of different tissues and physiological conditions. Yeast one-hybrid (Y1H) assays provide a powerful "gene-centered" method to identify the TFs that can bind a DNA sequence of interest. In this introduction, we describe the basic principles of the Y1H assay and its advantages and disadvantages and briefly discuss how it is complementary to "TF-centered" methods that identify protein-DNA interactions for a known protein of interest.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein–DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene‐centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two‐thirds of which are novel, as well as integrate gene expression, protein–protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.

DREISS: Using State-Space Models to Infer the Dynamics of Gene Expression Driven by External and Internal Regulatory Networks

Gene expression is controlled by the combinatorial effects of regulatory factors from different biological subsystems such as general transcription factors (TFs), cellular growth factors and microRNAs. A subsystem’s gene expression may be controlled by its internal regulatory factors, exclusively, or by external subsystems, or by both. It is thus useful to distinguish the degree to which a subsystem is regulated internally or externally–e.g., how non-conserved, species-specific TFs affect the expression of conserved, cross-species genes during evolution. We developed a computational method (DREISS, dreiss.gerteinlab.org) for analyzing the Dynamics of gene expression driven by Regulatory networks, both External and Internal based on State Space models. Given a subsystem, the “state” and “control” in the model refer to its own (internal) and another subsystem’s (external) gene expression levels. The state at a given time is determined by the state and control at a previous time. Because typical time-series data do not have enough samples to fully estimate the model’s parameters, DREISS uses dimensionality reduction, and identifies canonical temporal expression trajectories (e.g., degradation, growth and oscillation) representing the regulatory effects emanating from various subsystems. To demonstrate capabilities of DREISS, we study the regulatory effects of evolutionarily conserved vs. divergent TFs across distant species. In particular, we applied DREISS to the time-series gene expression datasets of C. elegans and D. melanogaster during their embryonic development. We analyzed the expression dynamics of the conserved, orthologous genes (orthologs), seeing the degree to which these can be accounted for by orthologous (internal) versus species-specific (external) TFs. We found that between two species, the orthologs have matched, internally driven expression patterns but very different externally driven ones. This is particularly true for genes with evolutionarily ancient functions (e.g. the ribosomal proteins), in contrast to those with more recently evolved functions (e.g., cell-cell communication). This suggests that despite striking morphological differences, some fundamental embryonic-developmental processes are still controlled by ancient regulatory systems.

The dynamics of a biological system can be controlled by its own internal mechanisms and external perturbations. To gain intuition on this, we may draw a comparison with a mass hanging from a spring. The mass will move naturally by itself but its dynamics is also affected by one’s pulling it. That is, the dynamics of the mass is governed by the effect of the external perturbations superimposed on the internal mechanism of the spring (i.e. Hooke’s law). Similarly, given a group of genes, their temporal gene expression dynamics can be controlled by both transcription factors inside the group and external regulatory factors. Therefore, it is useful to identify the expression dynamics that are exclusively controlled by internal or external factors and compare them across various systems. While state-space models have been widely used to decouple the internal and external effects in physical systems, such as the mass and spring, typical biological systems do not have enough time samples to infer all the model’s parameters, and applications of state-space models were not very effective in these instances. Hence, we developed a general-purpose computational method by integrating state-space models and dimensionality reduction to identify temporal gene expression patterns driven by internal and external regulatory networks. We applied our method to the embryonic developmental datasets in the worm and fly (and also in a human cancer context). We successfully identified the temporal expression dynamics of cross-species conserved genes that were driven by conserved and species-specific regulatory networks.

A cellular and regulatory map of the GABAergic nervous system of C. elegans

Neurotransmitter maps are important complements to anatomical maps and represent an invaluable resource to understand nervous system function and development. We report here a comprehensive map of neurons in the C. elegans nervous system that contain the neurotransmitter GABA, revealing twice as many GABA-positive neuron classes as previously reported. We define previously unknown glia-like cells that take up GABA, as well as 'GABA uptake neurons' which do not synthesize GABA but take it up from the extracellular environment, and we map the expression of previously uncharacterized ionotropic GABA receptors. We use the map of GABA-positive neurons for a comprehensive analysis of transcriptional regulators that define the GABA phenotype. We synthesize our findings of specification of GABAergic neurons with previous reports on the specification of glutamatergic and cholinergic neurons into a nervous system-wide regulatory map which defines neurotransmitter specification mechanisms for more than half of all neuron classes in C. elegans.

DOI:http://dx.doi.org/10.7554/eLife.17686.001

Big Data in Caenorhabditis elegans: quo vadis?

A clear definition of what constitutes “Big Data” is difficult to identify, but we find it most useful to define Big Data as a data collection that is complete. By this criterion, researchers on Caenorhabditis elegans have a long history of collecting Big Data, since the organism was selected with the idea of obtaining a complete biological description and understanding of development. The complete wiring diagram of the nervous system, the complete cell lineage, and the complete genome sequence provide a framework to phrase and test hypotheses. Given this history, it might be surprising that the number of “complete” data sets for this organism is actually rather small—not because of lack of effort, but because most types of biological experiments are not currently amenable to complete large-scale data collection. Many are also not inherently limited, so that it becomes difficult to even define completeness. At present, we only have partial data on mutated genes and their phenotypes, gene expression, and protein–protein interaction—important data for many biological questions. Big Data can point toward unexpected correlations, and these unexpected correlations can lead to novel investigations; however, Big Data cannot establish causation. As a result, there is much excitement about Big Data, but there is also a discussion on just what Big Data contributes to solving a biological problem. Because of its relative simplicity, C. elegans is an ideal test bed to explore this issue and at the same time determine what is necessary to build a multicellular organism from a single cell.

Cell-Autonomous and Non-Cell-Autonomous Regulation of a Feeding State-Dependent Chemoreceptor Gene via MEF-2 and bHLH Transcription Factors

Food and feeding-state dependent changes in chemoreceptor gene expression may allow Caenorhabditis elegans to modify their chemosensory behavior, but the mechanisms essential for these expression changes remain poorly characterized. We had previously shown that expression of a feeding state-dependent chemoreceptor gene, srh-234, in the ADL sensory neuron of C. elegans is regulated via the MEF-2 transcription factor. Here, we show that MEF-2 acts together with basic helix-loop-helix (bHLH) transcription factors to regulate srh-234 expression as a function of feeding state. We identify a cis-regulatory MEF2 binding site that is necessary and sufficient for the starvation-induced down regulation of srh-234 expression, while an E-box site known to bind bHLH factors is required to drive srh-234 expression in ADL. We show that HLH-2 (E/Daughterless), HLH-3 and HLH-4 (Achaete-scute homologs) act in ADL neurons to regulate srh-234 expression. We further demonstrate that the expression levels of srh-234 in ADL neurons are regulated remotely by MXL-3 (Max-like 3 homolog) and HLH-30 (TFEB ortholog) acting in the intestine, which is dependent on insulin signaling functioning specifically in ADL neurons. We also show that this intestine-to-neuron feeding-state regulation of srh-234 involves a subset of insulin-like peptides. These results combined suggest that chemoreceptor gene expression is regulated by both cell-autonomous and non-cell-autonomous transcriptional mechanisms mediated by MEF2 and bHLH factors, which may allow animals to fine-tune their chemosensory responses in response to changes in their feeding state.

Plasticity in chemoreceptor gene expression may be a simple strategy by which an animal can modulate its chemosensory responses in changing external and internal state conditions. However, the transcriptional mechanisms required for these chemoreceptor gene expression changes are poorly understood. Here, we describe the identification of a transcriptional module(s) consisting of MEF-2 and basic helix-loop-helix (bHLH) transcription factors and their cognate binding sites in Caenorhabditis elegans that act together in ADL sensory neurons to properly regulate expression of a feeding-state dependent chemoreceptor gene. We also showed that chemoreceptor gene expression in ADL neurons are regulated remotely by bHLH factors acting in the intestine through an insulin-mediated signaling pathway, implying a sensory neuron-gut interaction for modulating chemoreceptor gene expression as a function of feeding state. This work describes transcriptional mechanisms mediated by MEF-2 and bHLH factors by which the expression of individual chemoreceptor genes in C. elegans are changed in response to changes in feeding state conditions.

Mechanisms of transcription factor evolution in Metazoa

Transcriptions factors (TFs) are pivotal for the regulation of virtually all cellular processes, including growth and development. Expansions of TF families are causally linked to increases in organismal complexity. Here we study the evolutionary dynamics, genetic causes and functional implications of the five largest metazoan TF families. We find that family expansions dominate across the whole metazoan tree; however, some branches experience exceptional family-specific accelerated expansions. Additionally, we find that such expansions are often predated by modular domain rearrangements, which spur the expansion of a new sub-family by separating it from the rest of the TF family in terms of protein-protein interactions. This separation allows for radical shifts in the functional spectrum of a duplicated TF. We also find functional differentiation inside TF sub-families as changes in expression specificity. Furthermore, accelerated family expansions are facilitated by repeats of sequence motifs such as C2H2 zinc fingers. We quantify whole genome duplications and single gene duplications as sources of TF family expansions, implying that some, but not all, TF duplicates are preferentially retained. We conclude that trans-regulatory changes (domain rearrangements) are instrumental for fundamental functional innovations, that cis-regulatory changes (affecting expression) accomplish wide-spread fine tuning and both jointly contribute to the functional diversification of TFs.

Metabolic network rewiring of propionate flux compensates vitamin B12 deficiency in C. elegans

Metabolic network rewiring is the rerouting of metabolism through the use of alternate enzymes to adjust pathway flux and accomplish specific anabolic or catabolic objectives. Here, we report the first characterization of two parallel pathways for the breakdown of the short chain fatty acid propionate in Caenorhabditis elegans. Using genetic interaction mapping, gene co-expression analysis, pathway intermediate quantification and carbon tracing, we uncover a vitamin B12-independent propionate breakdown shunt that is transcriptionally activated on vitamin B12 deficient diets, or under genetic conditions mimicking the human diseases propionic- and methylmalonic acidemia, in which the canonical B12-dependent propionate breakdown pathway is blocked. Our study presents the first example of transcriptional vitamin-directed metabolic network rewiring to promote survival under vitamin deficiency. The ability to reroute propionate breakdown according to B12 availability may provide C. elegans with metabolic plasticity and thus a selective advantage on different diets in the wild.

Quantification of Cooperativity in Heterodimer-DNA Binding Improves the Accuracy of Binding Specificity Models

Many transcription factors (TFs) have the ability to cooperate on DNA elements as heterodimers. Despite the significance of TF heterodimerization for gene regulation, a quantitative understanding of cooperativity between various TF dimer partners and its impact on heterodimer DNA binding specificity models is still lacking. Here, we used a novel integrative approach, combining microfluidics-steered measurements of dimer-DNA assembly with mechanistic modeling of the implicated protein-protein-DNA interactions to quantitatively interrogate the cooperative DNA binding behavior of the adipogenic peroxisome proliferator-activated receptor γ (PPARγ):retinoid X receptor α (RXRα) heterodimer. Using the high throughput MITOMI (mechanically induced trapping of molecular interactions) platform, we derived equilibrium DNA binding data for PPARγ, RXRα, as well as the PPARγ:RXRα heterodimer to more than 300 target DNA sites and variants thereof. We then quantified cooperativity underlying heterodimer-DNA binding and derived an integrative heterodimer DNA binding constant. Using this cooperativity-inclusive constant, we were able to build a heterodimer-DNA binding specificity model that has superior predictive power than the one based on a regular one-site equilibrium. Our data further revealed that individual nucleotide substitutions within the target site affect the extent of cooperativity in PPARγ:RXRα-DNA binding. Our study therefore emphasizes the importance of assessing cooperativity when generating DNA binding specificity models for heterodimers.

Gene-diet interactions and aging in C. elegans

Diet is the most variable aspect of life history, as most individuals have a large diversity of food choices, varying in the type and amount that they ingest. In the short-term, diet can affect metabolism and energy levels. However, in the long run, the net deficiency or excess of calories from diet can influence the progression and severity of age-related diseases. An old and yet still debated question is: how do specific dietary choices impact health- and lifespan? It is clear that genetics can play a critical role - perhaps just as important as diet choices. For example, poor diet in combination with genetic susceptibility can lead to metabolic disorders, such as obesity and type 2 diabetes. Recent work in Caenorhabditis elegans has identified the existence of diet-gene pairs, where the consequence of mutating a specific gene is only realized on specific diets. Many core metabolic pathways are conserved from worm to human. Although only a handful of these diet-gene pairs has been characterized, there are potentially hundreds, if not thousands, of such interactions, which may explain the variability in the rates of aging in humans and the incidence and severity of age-related diseases.

Transcription Factor Activity Mapping of a Tissue-Specific in vivo Gene Regulatory Network

A wealth of physical interaction data between transcription factors (TFs) and DNA has been generated, but these interactions often do not have apparent regulatory consequences. Thus, equating physical interaction data with gene regulatory networks (GRNs) is problematic. Here, we comprehensively assay TF activity, rather than binding, to construct a network of gene regulatory interactions in the C. elegans intestine. By manually observing the in vivo tissue-specific knockdown of 921 TFs on a panel of 19 fluorescent transcriptional reporters, we identified a GRN of 411 interactions between 19 promoters and 177 TFs. This GRN shows only modest overlap with physical interactions, indicating that many regulatory interactions are indirect. We applied nested effects modeling to uncover information flow between TFs in the intestine that converges on a small set of physical TF-promoter interactions. We found numerous cell nonautonomous regulatory interactions, illustrating tissue-to-tissue communication. Altogether, our study illuminates the complexity of gene regulation in the context of a living animal.

The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis

Homeobox genes play crucial roles for the development of multicellular eukaryotes. We have generated a revised list of all homeobox genes for Caenorhabditis elegans and provide a nomenclature for the previously unnamed ones. We show that, out of 103 homeobox genes, 70 are co-orthologous to human homeobox genes. 14 are highly divergent, lacking an obvious ortholog even in other Caenorhabditis species. One of these homeobox genes encodes 12 homeodomains, while three other highly divergent homeobox genes encode a novel type of double homeodomain, termed HOCHOB. To understand how transcription factors regulate cell fate during development, precise spatio-temporal expression data need to be obtained. Using a new imaging framework that we developed, Endrov, we have generated spatio-temporal expression profiles during embryogenesis of over 60 homeobox genes, as well as a number of other developmental control genes using GFP reporters. We used dynamic feedback during recording to automatically adjust the camera exposure time in order to increase the dynamic range beyond the limitations of the camera. We have applied the new framework to examine homeobox gene expression patterns and provide an analysis of these patterns. The methods we developed to analyze and quantify expression data are not only suitable for C. elegans, but can be applied to other model systems or even to tissue culture systems.

Genome-Wide Identification of the Target Genes of AP2-O, a Plasmodium AP2-Family Transcription Factor

Stage-specific transcription is a fundamental biological process in the life cycle of the Plasmodium parasite. Proteins containing the AP2 DNA-binding domain are responsible for stage-specific transcriptional regulation and belong to the only known family of transcription factors in Plasmodium parasites. Comprehensive identification of their target genes will advance our understanding of the molecular basis of stage-specific transcriptional regulation and stage-specific parasite development. AP2-O is an AP2 family transcription factor that is expressed in the mosquito midgut-invading stage, called the ookinete, and is essential for normal morphogenesis of this stage. In this study, we identified the genome-wide target genes of AP2-O by chromatin immunoprecipitation-sequencing and elucidate how this AP2 family transcription factor contributes to the formation of this motile stage. The analysis revealed that AP2-O binds specifically to the upstream genomic regions of more than 500 genes, suggesting that approximately 10% of the parasite genome is directly regulated by AP2-O. These genes are involved in distinct biological processes such as morphogenesis, locomotion, midgut penetration, protection against mosquito immunity and preparation for subsequent oocyst development. This direct and global regulation by AP2-O provides a model for gene regulation in Plasmodium parasites and may explain how these parasites manage to control their complex life cycle using a small number of sequence-specific AP2 transcription factors.

Although malarial parasites have a complex life cycle, they harbor only 30 transcription factors in their genome. The majority of these transcription factors belong to a single family referred to as the AP2 family. Our previous study suggested that stage-specific AP2 family transcription factors have critical roles in maintaining the Plasmodium parasite life cycle. However, it remains fairly elusive as to how these transcription factors regulate each stage. AP2-O is an AP2 family transcription factor that is expressed during the mosquito midgut-invading stage, the ookinete, and is essential for normal development of this stage. In the present study, we identified the entire set of AP2-O target genes to elucidate how this AP2 family transcription factor contributes to the formation of this stage. Our results showed that AP2-O directly regulates 10% of the parasite genome and is involved in the whole process of mosquito midgut-invasion by ookinetes. The global and comprehensive regulation by the AP2 family transcription factor that we revealed provides a model for transcriptional regulation of this parasite and may explain how malarial parasites regulate their complex life cycle using a small number of sequence-specific transcription factors.

Cis- and trans-regulatory mechanisms of gene expression in the ASJ sensory neuron of Caenorhabditis elegans

The identity of a given cell type is determined by the expression of a set of genes sharing common cis-regulatory motifs and being regulated by shared transcription factors. Here, we identify cis and trans regulatory elements that drive gene expression in the bilateral sensory neuron ASJ, located in the head of the nematode Caenorhabditis elegans. For this purpose, we have dissected the promoters of the only two genes so far reported to be exclusively expressed in ASJ, trx-1 and ssu-1. We hereby identify the ASJ motif, a functional cis-regulatory bipartite promoter region composed of two individual 6 bp elements separated by a 3 bp linker. The first element is a 6 bp CG-rich sequence that presumably binds the Sp family member zinc-finger transcription factor SPTF-1. Interestingly, within the C. elegans nervous system SPTF-1 is also found to be expressed only in ASJ neurons where it regulates expression of other genes in these neurons and ASJ cell fate. The second element of the bipartite motif is a 6 bp AT-rich sequence that is predicted to potentially bind a transcription factor of the homeobox family. Together, our findings identify a specific promoter signature and SPTF-1 as a transcription factor that functions as a terminal selector gene to regulate gene expression in C. elegans ASJ sensory neurons.