Mapping and analysis of Caenorhabditis elegans transcription factor sequence specificities

LET-381/FoxF and its target UNC-30/Pitx2 specify and maintain the molecular identity of C. elegans mesodermal glia that regulate motor behavior

While most glial cell types in the central nervous system (CNS) arise from neuroectodermal progenitors, some, like microglia, are mesodermally derived. To understand mesodermal glia development and function, we investigated C. elegans GLR glia, which envelop the brain neuropil and separate it from the circulatory system cavity. Transcriptome analysis shows that GLR glia combine astrocytic and endothelial characteristics, which are relegated to separate cell types in vertebrates. Combined fate acquisition is orchestrated by LET-381/FoxF, a fate-specification/maintenance transcription factor also expressed in glia and endothelia of other animals. Among LET-381/FoxF targets, the UNC-30/Pitx2 transcription factor controls GLR glia morphology and represses alternative mesodermal fates. LET-381 and UNC-30 co-expression in naive cells is sufficient for GLR glia gene expression. GLR glia inactivation by ablation or let-381 mutation disrupts locomotory behavior and promotes salt-induced paralysis, suggesting brain-neuropil activity dysregulation. Our studies uncover mechanisms of mesodermal glia development and show that like neuronal differentiation, glia differentiation requires autoregulatory terminal selector genes that define and maintain the glial fate.

Gene regulatory patterning codes in early cell fate specification of the C. elegans embryo

Pattern formation originates during embryogenesis by a series of symmetry-breaking steps throughout an expanding cell lineage. In Drosophila, classic work has shown that segmentation in the embryo is established by morphogens within a syncytium, and the subsequent action of the gap, pair-rule, and segment polarity genes. This classic model however does not translate directly to species that lack a syncytium - such as Caenorhabditis elegans - where cell fate is specified by cell-autonomous cell lineage programs and their inter-signaling. Previous single-cell RNA-Seq studies in C. elegans have analyzed cells from a mixed suspension of cells from many embryos to study late differentiation stages, or individual early stage embryos to study early gene expression in the embryo. To study the intermediate stages of early and late gastrulation (28- to 102-cells stages) missed by these approaches, here we determine the transcriptomes of the 1- to 102-cell stage to identify 119 embryonic cell states during cell fate specification, including 'equivalence-group' cell identities. We find that gene expression programs are modular according to the sub-cell lineages, each establishing a set of stripes by combinations of transcription factor gene expression across the anterior-posterior axis. In particular, expression of the homeodomain genes establishes a comprehensive lineage-specific positioning system throughout the embryo beginning at the 28-cell stage. Moreover, we find that genes that segment the entire embryo in Drosophila have orthologs in C. elegans that exhibit sub-lineage-specific expression. These results suggest that the C. elegans embryo is patterned by a juxtaposition of distinct lineage-specific gene regulatory programs each with a unique encoding of cell location and fate. This use of homologous gene regulatory patterning codes suggests a deep homology of cell fate specification programs across diverse modes of development.

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.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

LET-381/FoxF and UNC-30/Pitx2 control the development of C. elegans mesodermal glia that regulate motor behavior

While most CNS glia arise from neuroectodermal progenitors, some, like microglia, are mesodermally derived. To understand mesodermal glia development and function, we investigated C. elegans GLR glia, which ensheath the brain neuropil and separate it from the circulatory-system cavity. Transcriptome analysis suggests GLR glia merge astrocytic and endothelial characteristics relegated to separate cell types in vertebrates. Combined fate acquisition is orchestrated by LET-381/FoxF, a fate-specification/maintenance transcription factor expressed in glia and endothelia of other animals. Among LET-381/FoxF targets, UNC-30/Pitx2 transcription factor controls GLR glia morphology and represses alternative mesodermal fates. LET-381 and UNC-30 co-expression in naïve cells is sufficient for GLR glia gene expression. GLR glia inactivation by ablation or let-381 mutation disrupts locomotory behavior and induces salt hypersensitivity, suggesting brain-neuropil activity dysregulation. Our studies uncover mechanisms of mesodermal glia development and show that like neurons, glia differentiation requires autoregulatory terminal selector genes that define and maintain the glial fate.

Spaceflight Induces Strength Decline in Caenorhabditis elegans

Background: Understanding and countering the well-established negative health consequences of spaceflight remains a primary challenge preventing safe deep space exploration. Targeted/personalized therapeutics are at the forefront of space medicine strategies, and cross-species molecular signatures now define the 'typical' spaceflight response. However, a lack of direct genotype-phenotype associations currently limits the robustness and, therefore, the therapeutic utility of putative mechanisms underpinning pathological changes in flight. Methods: We employed the worm Caenorhabditis elegans as a validated model of space biology, combined with 'NemaFlex-S' microfluidic devices for assessing animal strength production as one of the most reproducible physiological responses to spaceflight. Wild-type and dys-1 (BZ33) strains (a Duchenne muscular dystrophy (DMD) model for comparing predisposed muscle weak animals) were cultured on the International Space Station in chemically defined media before loading second-generation gravid adults into NemaFlex-S devices to assess individual animal strength. These same cultures were then frozen on orbit before returning to Earth for next-generation sequencing transcriptomic analysis. Results: Neuromuscular strength was lower in flight versus ground controls (16.6% decline, p < 0.05), with dys-1 significantly more (23% less strength, p < 0.01) affected than wild types. The transcriptional gene ontology signatures characterizing both strains of weaker animals in flight strongly corroborate previous results across species, enriched for upregulated stress response pathways and downregulated mitochondrial and cytoskeletal processes. Functional gene cluster analysis extended this to implicate decreased neuronal function, including abnormal calcium handling and acetylcholine signaling, in space-induced strength declines under the predicted control of UNC-89 and DAF-19 transcription factors. Finally, gene modules specifically altered in dys-1 animals in flight again cluster to neuronal/neuromuscular pathways, suggesting strength loss in DMD comprises a strong neuronal component that predisposes these animals to exacerbated strength loss in space. Conclusions: Highly reproducible gene signatures are strongly associated with space-induced neuromuscular strength loss across species and neuronal changes in calcium/acetylcholine signaling require further study. These results promote targeted medical efforts towards and provide an in vivo model for safely sending animals and people into deep space in the near future.

Histone methyltransferase activity affects metabolism in human cells independently of transcriptional regulation

The N-terminal tails of eukaryotic histones are frequently posttranslationally modified. The role of these modifications in transcriptional regulation is well-documented. However, the extent to which the enzymatic processes of histone posttranslational modification might affect metabolic regulation is less clear. Here, we investigated how histone methylation might affect metabolism using metabolomics, proteomics, and RNA-seq data from cancer cell lines, primary tumour samples and healthy tissue samples. In cancer, the expression of histone methyltransferases (HMTs) was inversely correlated to the activity of NNMT, an enzyme previously characterised as a methyl sink that disposes of excess methyl groups carried by the universal methyl donor S-adenosyl methionine (SAM or AdoMet). In healthy tissues, histone methylation was inversely correlated to the levels of an alternative methyl sink, PEMT. These associations affected the levels of multiple histone marks on chromatin genome-wide but had no detectable impact on transcriptional regulation. We show that HMTs with a variety of different associations to transcription are co-regulated by the Retinoblastoma (Rb) tumour suppressor in human cells. Rb-mutant cancers show increased total HMT activity and down-regulation of NNMT. Together, our results suggest that the total activity of HMTs affects SAM metabolism, independent of transcriptional regulation.

Forkhead transcription factor FKH-8 cooperates with RFX in the direct regulation of sensory cilia in Caenorhabditis elegans

Cilia, either motile or non-motile (a.k.a primary or sensory), are complex evolutionarily conserved eukaryotic structures composed of hundreds of proteins required for their assembly, structure and function that are collectively known as the ciliome. Ciliome gene mutations underlie a group of pleiotropic genetic diseases known as ciliopathies. Proper cilium function requires the tight coregulation of ciliome gene transcription, which is only fragmentarily understood. RFX transcription factors (TF) have an evolutionarily conserved role in the direct activation of ciliome genes both in motile and non-motile cilia cell-types. In vertebrates, FoxJ1 and FoxN4 Forkhead (FKH) TFs work with RFX in the direct activation of ciliome genes, exclusively in motile cilia cell-types. No additional TFs have been described to act together with RFX in primary cilia cell-types in any organism. Here we describe FKH-8, a FKH TF, as a direct regulator of the sensory ciliome genes in Caenorhabditis elegans. FKH-8 is expressed in all ciliated neurons in C. elegans, binds the regulatory regions of ciliome genes, regulates ciliome gene expression, cilium morphology and a wide range of behaviors mediated by sensory ciliated neurons. FKH-8 and DAF-19 (C. elegans RFX) physically interact and synergistically regulate ciliome gene expression. C. elegans FKH-8 function can be replaced by mouse FOXJ1 and FOXN4 but not by other members of other mouse FKH subfamilies. In conclusion, RFX and FKH TF families act jointly as direct regulators of ciliome genes also in sensory ciliated cell types suggesting that this regulatory logic could be an ancient trait predating functional cilia sub-specialization.

Thousands of Pristionchus pacificus orphan genes were integrated into developmental networks that respond to diverse environmental microbiota

Adaptation of organisms to environmental change may be facilitated by the creation of new genes. New genes without homologs in other lineages are known as taxonomically-restricted orphan genes and may result from divergence or de novo formation. Previously, we have extensively characterized the evolution and origin of such orphan genes in the nematode model organism Pristionchus pacificus. Here, we employ large-scale transcriptomics to establish potential functional associations and to measure the degree of transcriptional plasticity among orphan genes. Specifically, we analyzed 24 RNA-seq samples from adult P. pacificus worms raised on 24 different monoxenic bacterial cultures. Based on coexpression analysis, we identified 28 large modules that harbor 3,727 diplogastrid-specific orphan genes and that respond dynamically to different bacteria. These coexpression modules have distinct regulatory architecture and also exhibit differential expression patterns across development suggesting a link between bacterial response networks and development. Phylostratigraphy revealed a considerably high number of family- and even species-specific orphan genes in certain coexpression modules. This suggests that new genes are not attached randomly to existing cellular networks and that integration can happen very fast. Integrative analysis of protein domains, gene expression and ortholog data facilitated the assignments of biological labels for 22 coexpression modules with one of the largest, fast-evolving module being associated with spermatogenesis. In summary, this work presents the first functional annotation for thousands of P. pacificus orphan genes and reveals insights into their integration into environmentally responsive gene networks.

Characterizing a standardized BioPart for PVQ-specific expression in C. elegans

Synthetic biology relies on standardized biological parts (BioParts), and we aim to identify cell-specific promoters for every class of neuron in C. elegans . Here, we characterize a short BioPart (P nlp-17 , 300 bp) for PVQ-specific expression. P nlp-17 ::mScarlet showed bright, persistent, and specific expression in hermaphrodite and male PVQ neurons from multicopy arrays and single-copy insertions starting from the comma stage. We generated standardized P nlp-17 cloning vectors with gfp and mScarlet compatible with single-copy or array expression for PVQ-specific transgene expression or identification. To facilitate gene synthesis, we have incorporated P nlp-17 as a standard BioPart in our online transgene design tool (www.wormbuilder.org/transgenebuilder).

Transcellular chaperone signaling is an intercellular stress-response distinct from the HSF-1-mediated heat shock response

Organismal proteostasis is maintained by intercellular signaling processes including cell nonautonomous stress responses such as transcellular chaperone signaling (TCS). When TCS is activated upon tissue-specific knockdown of hsp-90 in the Caenorhabditis elegans intestine, heat-inducible hsp-70 is induced in muscle cells at the permissive temperature resulting in increased heat stress resistance and lifespan extension. However, our understanding of the molecular mechanism and signaling factors mediating transcellular activation of hsp-70 expression from one tissue to another is still in its infancy. Here, we conducted a combinatorial approach using transcriptome RNA-Seq profiling and a forward genetic mutagenesis screen to elucidate how stress signaling from the intestine to the muscle is regulated. We find that the TCS-mediated "gut-to-muscle" induction of hsp-70 expression is suppressed by HSF-1 and instead relies on transcellular-X-cross-tissue (txt) genes. We identify a key role for the PDZ-domain guanylate cyclase txt-1 and the homeobox transcription factor ceh-58 as signaling hubs in the stress receiving muscle cells to initiate hsp-70 expression and facilitate TCS-mediated heat stress resistance and lifespan extension. Our results provide a new view on cell-nonautonomous regulation of "inter-tissue" stress responses in an organism that highlight a key role for the gut. Our data suggest that the HSF-1-mediated heat shock response is switched off upon TCS activation, in favor of an intercellular stress-signaling route to safeguard survival.

Pathogenic bacteria modulate pheromone response to promote mating

Pathogens generate ubiquitous selective pressures and host-pathogen interactions alter social behaviours in many animals1-4. However, very little is known about the neuronal mechanisms underlying pathogen-induced changes in social behaviour. Here we show that in adult Caenorhabditis elegans hermaphrodites, exposure to a bacterial pathogen (Pseudomonas aeruginosa) modulates sensory responses to pheromones by inducing the expression of the chemoreceptor STR-44 to promote mating. Under standard conditions, C. elegans hermaphrodites avoid a mixture of ascaroside pheromones to facilitate dispersal5-13. We find that exposure to the pathogenic Pseudomonas bacteria enables pheromone responses in AWA sensory neurons, which mediate attractive chemotaxis, to suppress the avoidance. Pathogen exposure induces str-44 expression in AWA neurons, a process regulated by a transcription factor zip-5 that also displays a pathogen-induced increase in expression in AWA. STR-44 acts as a pheromone receptor and its function in AWA neurons is required for pathogen-induced AWA pheromone response and suppression of pheromone avoidance. Furthermore, we show that C. elegans hermaphrodites, which reproduce mainly through self-fertilization, increase the rate of mating with males after pathogen exposure and that this increase requires str-44 in AWA neurons. Thus, our results uncover a causal mechanism for pathogen-induced social behaviour plasticity, which can promote genetic diversity and facilitate adaptation of the host animals.

Widespread employment of conserved C. elegans homeobox genes in neuronal identity specification

Homeobox genes are prominent regulators of neuronal identity, but the extent to which their function has been probed in animal nervous systems remains limited. In the nematode Caenorhabditis elegans, each individual neuron class is defined by the expression of unique combinations of homeobox genes, prompting the question of whether each neuron class indeed requires a homeobox gene for its proper identity specification. We present here progress in addressing this question by extending previous mutant analysis of homeobox gene family members and describing multiple examples of homeobox gene function in different parts of the C. elegans nervous system. To probe homeobox function, we make use of a number of reporter gene tools, including a novel multicolor reporter transgene, NeuroPAL, which permits simultaneous monitoring of the execution of multiple differentiation programs throughout the entire nervous system. Using these tools, we add to the previous characterization of homeobox gene function by identifying neuronal differentiation defects for 14 homeobox genes in 24 distinct neuron classes that are mostly unrelated by location, function and lineage history. 12 of these 24 neuron classes had no homeobox gene function ascribed to them before, while in the other 12 neuron classes, we extend the combinatorial code of transcription factors required for specifying terminal differentiation programs. Furthermore, we demonstrate that in a particular lineage, homeotic identity transformations occur upon loss of a homeobox gene and we show that these transformations are the result of changes in homeobox codes. Combining the present with past analyses, 113 of the 118 neuron classes of C. elegans are now known to require a homeobox gene for proper execution of terminal differentiation programs. Such broad deployment indicates that homeobox function in neuronal identity specification may be an ancestral feature of animal nervous systems.

A natural transdifferentiation event involving mitosis is empowered by integrating signaling inputs with conserved plasticity factors

Transdifferentiation, or direct cell reprogramming, is the conversion of one fully differentiated cell type into another. Whether core mechanisms are shared between natural transdifferentiation events when occurring with or without cell division is unclear. We have previously characterized the Y-to-PDA natural transdifferentiation in Caenorhabditis elegans, which occurs without cell division and requires orthologs of vertebrate reprogramming factors. Here, we identify a rectal-to-GABAergic transdifferentiation and show that cell division is required but not sufficient for conversion. We find shared mechanisms, including erasure of the initial identity, which requires the conserved reprogramming factors SEM-4/SALL, SOX-2, CEH-6/OCT, and EGL-5/HOX. We also find three additional and parallel roles of the Wnt signaling pathway: selection of a specific daughter, removal of the initial identity, and imposition of the precise final subtype identity. Our results support a model in which levels and antagonistic activities of SOX-2 and Wnt signaling provide a timer for the acquisition of final identity.

pop-1/TCF, ref-2/ZIC and T-box factors regulate the development of anterior cells in the C. elegans embryo

Patterning of the anterior-posterior axis is fundamental to animal development. The Wnt pathway plays a major role in this process by activating the expression of posterior genes in animals from worms to humans. This observation raises the question of whether the Wnt pathway or other regulators control the expression of the many anterior-expressed genes. We found that the expression of five anterior-specific genes in Caenorhabditis elegans embryos depends on the Wnt pathway effectors pop-1/TCF and sys-1/β-catenin. We focused further on one of these anterior genes, ref-2/ZIC, a conserved transcription factor expressed in multiple anterior lineages. Live imaging of ref-2 mutant embryos identified defects in cell division timing and position in anterior lineages. Cis-regulatory dissection identified three ref-2 transcriptional enhancers, one of which is necessary and sufficient for anterior-specific expression. This enhancer is activated by the T-box transcription factors TBX-37 and TBX-38, and surprisingly, concatemerized TBX-37/38 binding sites are sufficient to drive anterior-biased expression alone, despite the broad expression of TBX-37 and TBX-38. Taken together, our results highlight the diverse mechanisms used to regulate anterior expression patterns in the embryo.

The regulatory landscape of neurite development in Caenorhabditis elegans

Neuronal communication requires precise connectivity of neurite projections (axons and dendrites). Developing neurites express cell-surface receptors that interpret extracellular cues to enable correct guidance toward, and connection with, target cells. Spatiotemporal regulation of neurite guidance molecule expression by transcription factors (TFs) is critical for nervous system development and function. Here, we review how neurite development is regulated by TFs in the Caenorhabditis elegans nervous system. By collecting publicly available transcriptome and ChIP-sequencing data, we reveal gene expression dynamics during neurite development, providing insight into transcriptional mechanisms governing construction of the nervous system architecture.

The zinc-finger transcription factor LSL-1 is a major regulator of the germline transcriptional program in Caenorhabditis elegans

Specific gene transcriptional programs are required to ensure the proper proliferation and differentiation processes underlying the production of specialized cells during development. Gene activity is mainly regulated by the concerted action of transcription factors and chromatin proteins. In the nematode Caenorhabditis elegans, mechanisms that silence improper transcriptional programs in germline and somatic cells have been well studied, however, how are tissue-specific sets of genes turned on is less known. LSL-1 is herein defined as a novel crucial transcriptional regulator of germline genes in C. elegans. LSL-1 is first detected in the P4 blastomere and remains present at all stages of germline development, from primordial germ cell proliferation to the end of meiotic prophase. lsl-1 loss-of-function mutants exhibit many defects including meiotic prophase progression delay, a high level of germline apoptosis, and production of almost no functional gametes. Transcriptomic analysis and ChIP-seq data show that LSL-1 binds to promoters and acts as a transcriptional activator of germline genes involved in various processes, including homologous chromosome pairing, recombination, and genome stability. Furthermore, we show that LSL-1 functions by antagonizing the action of the heterochromatin proteins HPL-2/HP1 and LET-418/Mi2 known to be involved in the repression of germline genes in somatic cells. Based on our results, we propose LSL-1 to be a major regulator of the germline transcriptional program during development.

The anterior Hox gene ceh-13 and elt-1/GATA activate the posterior Hox genes nob-1 and php-3 to specify posterior lineages in the C. elegans embryo

Hox transcription factors play a conserved role in specifying positional identity during animal development, with posterior Hox genes typically repressing the expression of more anterior Hox genes. Here, we dissect the regulation of the posterior Hox genes nob-1 and php-3 in the nematode C. elegans. We show that nob-1 and php-3 are co-expressed in gastrulation-stage embryos in cells that previously expressed the anterior Hox gene ceh-13. This expression is controlled by several partially redundant transcriptional enhancers. These enhancers act in a ceh-13-dependant manner, providing a striking example of an anterior Hox gene positively regulating a posterior Hox gene. Several other regulators also act positively through nob-1/php-3 enhancers, including elt-1/GATA, ceh-20/ceh-40/Pbx, unc-62/Meis, pop-1/TCF, ceh-36/Otx, and unc-30/Pitx. We identified defects in both cell position and cell division patterns in ceh-13 and nob-1;php-3 mutants, suggesting that these factors regulate lineage identity in addition to positional identity. Together, our results highlight the complexity and flexibility of Hox gene regulation and function and the ability of developmental transcription factors to regulate different targets in different stages of development.

The C. elegans regulatory factor X (RFX) DAF-19M module: A shift from general ciliogenesis to cell-specific ciliary and behavioral specialization

Cilia are important for the interaction with environments and the proper function of tissues. While the basic structure of cilia is well conserved, ciliated cells have various functions. To understand the distinctive identities of ciliated cells, the identification of cell-specific proteins and its regulation is essential. Here, we report the mechanism that confers a specific identity on IL2 neurons in Caenorhabditis elegans, neurons important for the dauer larva-specific nictation behavior. We show that DAF-19M, an isoform of the sole C. elegans RFX transcription factor DAF-19, heads a regulatory subroutine, regulating target genes through an X-box motif variant under the control of terminal selector proteins UNC-86 and CFI-1 in IL2 neurons. Considering the conservation of DAF-19M module in IL2 neurons for nictation and in male-specific neurons for mating behavior, we propose the existence of an evolutionarily adaptable, hard-wired genetic module for distinct behaviors that share the feature "recognizing the environment."

Evolutionarily conserved transcription factors as regulators of longevity and targets for geroprotection

Aging is the single largest risk factor for many debilitating conditions, including heart diseases, stroke, cancer, diabetes, and neurodegenerative disorders. While far from understood in its full complexity, it is scientifically well-established that aging is influenced by genetic and environmental factors, and can be modulated by various interventions. One of aging's early hallmarks are aberrations in transcriptional networks, controlling for example metabolic homeostasis or the response to stress. Evidence in different model organisms abounds that a number of evolutionarily conserved transcription factors, which control such networks, can affect lifespan and healthspan across species. These transcription factors thus potentially represent conserved regulators of longevity and are emerging as important targets in the challenging quest to develop treatments to mitigate age-related diseases, and possibly even to slow aging itself. This review provides an overview of evolutionarily conserved transcription factors that impact longevity or age-related diseases in at least one multicellular model organism (nematodes, flies, or mice), and/or are tentatively linked to human aging. Discussed is the general evidence for transcriptional regulation of aging and disease, followed by a more detailed look at selected transcription factor families, the common metabolic pathways involved, and the targeting of transcription factors as a strategy for geroprotective interventions.

The enteric nervous system of the C. elegans pharynx is specified by the Sine oculis-like homeobox gene ceh-34

Overarching themes in the terminal differentiation of the enteric nervous system, an autonomously acting unit of animal nervous systems, have so far eluded discovery. We describe here the overall regulatory logic of enteric nervous system differentiation of the nematode Caenorhabditis elegans that resides within the foregut (pharynx) of the worm. A C. elegans homolog of the Drosophila Sine oculis homeobox gene, ceh-34, is expressed in all 14 classes of interconnected pharyngeal neurons from their birth throughout their life time, but in no other neuron type of the entire animal. Constitutive and temporally controlled ceh-34 removal shows that ceh-34 is required to initiate and maintain the neuron type-specific terminal differentiation program of all pharyngeal neuron classes, including their circuit assembly. Through additional genetic loss of function analysis, we show that within each pharyngeal neuron class, ceh-34 cooperates with different homeodomain transcription factors to individuate distinct pharyngeal neuron classes. Our analysis underscores the critical role of homeobox genes in neuronal identity specification and links them to the control of neuronal circuit assembly of the enteric nervous system. Together with the pharyngeal nervous system simplicity as well as its specification by a Sine oculis homolog, our findings invite speculations about the early evolution of nervous systems.

Ageing induces tissue‐specific transcriptomic changes in Caenorhabditis elegans

Ageing is a complex process with common and distinct features across tissues. Unveiling the underlying processes driving ageing in individual tissues is indispensable to decipher the mechanisms of organismal longevity. Caenorhabditis elegans is a well‐established model organism that has spearheaded ageing research with the discovery of numerous genetic pathways controlling its lifespan. However, it remains challenging to dissect the ageing of worm tissues due to the limited description of tissue pathology and access to tissue‐specific molecular changes during ageing. In this study, we isolated cells from five major tissues in young and old worms and profiled the age‐induced transcriptomic changes within these tissues. We observed a striking diversity of ageing across tissues and identified different sets of longevity regulators therein. In addition, we found novel tissue‐specific factors, including irx‐1 and myrf‐2, which control the integrity of the intestinal barrier and sarcomere structure during ageing respectively. This study demonstrates the complexity of ageing across worm tissues and highlights the power of tissue‐specific transcriptomic profiling during ageing, which can serve as a resource to the field.

Joint actions of diverse transcription factor families establish neuron-type identities and promote enhancer selectivity

To systematically investigate the complexity of neuron specification regulatory networks, we performed an RNA interference (RNAi) screen against all 875 transcription factors (TFs) encoded in Caenorhabditis elegans genome and searched for defects in nine different neuron types of the monoaminergic (MA) superclass and two cholinergic motoneurons. We identified 91 TF candidates to be required for correct generation of these neuron types, of which 28 were confirmed by mutant analysis. We found that correct reporter expression in each individual neuron type requires at least nine different TFs. Individual neuron types do not usually share TFs involved in their specification but share a common pattern of TFs belonging to the five most common TF families: homeodomain (HD), basic helix loop helix (bHLH), zinc finger (ZF), basic leucine zipper domain (bZIP), and nuclear hormone receptors (NHR). HD TF members are overrepresented, supporting a key role for this family in the establishment of neuronal identities. These five TF families are also prevalent when considering mutant alleles with previously reported neuronal phenotypes in C. elegans, Drosophila, and mouse. In addition, we studied terminal differentiation complexity focusing on the dopaminergic terminal regulatory program. We found two HD TFs (UNC-62 and VAB-3) that work together with known dopaminergic terminal selectors (AST-1, CEH-43, CEH-20). Combined TF binding sites for these five TFs constitute a cis-regulatory signature enriched in the regulatory regions of dopaminergic effector genes. Our results provide new insights on neuron-type regulatory programs in C. elegans that could help better understand neuron specification and evolution of neuron types.

Distinct structural bases for sequence-specific DNA binding by mammalian BEN domain proteins

The BEN domain is a recently recognized DNA binding module that is present in diverse metazoans and certain viruses. Several BEN domain factors are known as transcriptional repressors, but, overall, relatively little is known of how BEN factors identify their targets in humans. In particular, X-ray structures of BEN domain:DNA complexes are only known for Drosophila factors bearing a single BEN domain, which lack direct vertebrate orthologs. Here, we characterize several mammalian BEN domain (BD) factors, including from two NACC family BTB-BEN proteins and from BEND3, which has four BDs. In vitro selection data revealed sequence-specific binding activities of isolated BEN domains from all of these factors. We conducted detailed functional, genomic, and structural studies of BEND3. We show that BD4 is a major determinant for in vivo association and repression of endogenous BEND3 targets. We obtained a high-resolution structure of BEND3-BD4 bound to its preferred binding site, which reveals how BEND3 identifies cognate DNA targets and shows differences with one of its non-DNA-binding BEN domains (BD1). Finally, comparison with our previous invertebrate BEN structures, along with additional structural predictions using AlphaFold2 and RoseTTAFold, reveal distinct strategies for target DNA recognition by different types of BEN domain proteins. Together, these studies expand the DNA recognition activities of BEN factors and provide structural insights into sequence-specific DNA binding by mammalian BEN proteins.

DREAM represses distinct targets by cooperating with different THAP domain proteins

The DREAM (dimerization partner [DP], retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell-cycle and other genes, but its mechanism of action is unclear. Here, we demonstrate that two C. elegans THAP domain proteins, LIN-15B and LIN-36, co-localize with DREAM and function by different mechanisms for repression of distinct sets of targets. LIN-36 represses classical cell-cycle targets by promoting DREAM binding and gene body enrichment of H2A.Z, and we find that DREAM subunit EFL-1/E2F is specific for LIN-36 targets. In contrast, LIN-15B represses germline-specific targets in the soma by facilitating H3K9me2 promoter marking. We further find that LIN-36 and LIN-15B differently regulate DREAM binding. In humans, THAP proteins have been implicated in cell-cycle regulation by poorly understood mechanisms. We propose that THAP domain proteins are key mediators of Rb/DREAM function.

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.

In silico analysis of the transcriptional regulatory logic of neuronal identity specification throughout the C. elegans nervous system

The generation of the enormous diversity of neuronal cell types in a differentiating nervous system entails the activation of neuron type-specific gene batteries. To examine the regulatory logic that controls the expression of neuron type-specific gene batteries, we interrogate single cell expression profiles of all 118 neuron classes of the Caenorhabditis elegans nervous system for the presence of DNA binding motifs of 136 neuronally expressed C. elegans transcription factors. Using a phylogenetic footprinting pipeline, we identify cis-regulatory motif enrichments among neuron class-specific gene batteries and we identify cognate transcription factors for 117 of the 118 neuron classes. In addition to predicting novel regulators of neuronal identities, our nervous system-wide analysis at single cell resolution supports the hypothesis that many transcription factors directly co-regulate the cohort of effector genes that define a neuron type, thereby corroborating the concept of so-called terminal selectors of neuronal identity. Our analysis provides a blueprint for how individual components of an entire nervous system are genetically specified.

H3K9me selectively blocks transcription factor activity and ensures differentiated tissue integrity

The developmental role of histone H3K9 methylation (H3K9me), which typifies heterochromatin, remains unclear. In Caenorhabditis elegans, loss of H3K9me leads to a highly divergent upregulation of genes with tissue and developmental-stage specificity. During development H3K9me is lost from differentiated cell type-specific genes and gained at genes expressed in earlier developmental stages or other tissues. The continuous deposition of H3K9me2 by the SETDB1 homolog MET-2 after terminal differentiation is necessary to maintain repression. In differentiated tissues, H3K9me ensures silencing by restricting the activity of a defined set of transcription factors at promoters and enhancers. Increased chromatin accessibility following the loss of H3K9me is neither sufficient nor necessary to drive transcription. Increased ATAC-seq signal and gene expression correlate at a subset of loci positioned away from the nuclear envelope, while derepressed genes at the nuclear periphery remain poorly accessible despite being transcribed. In conclusion, H3K9me deposition can confer tissue-specific gene expression and maintain the integrity of terminally differentiated muscle by restricting transcription factor activity.

Evolution of Developmental GATA Factors in Nematodes

GATA transcription factors are found in animals, plants, and fungi. In animals, they have important developmental roles in controlling specification of cell identities and executing tissue-specific differentiation. The Phylum Nematoda is a diverse group of vermiform animals that inhabit ecological niches all over the world. Both free-living and parasitic species are known, including those that cause human infectious disease. To date, GATA factors in nematodes have been studied almost exclusively in the model system C. elegans and its close relatives. In this study, we use newly available sequences to identify GATA factors across the nematode phylum. We find that most species have fewer than six GATA factors, but some species have 10 or more. Comparisons of gene and protein structure suggest that there were at most two GATA factors at the base of the phylum, which expanded by duplication and modification to result in a core set of four factors. The high degree of structural similarity with the corresponding orthologues in C. elegans suggests that the nematode GATA factors share similar functions in development.

Gli3 utilizes Hand2 to synergistically regulate tissue-specific transcriptional networks

Despite a common understanding that Gli TFs are utilized to convey a Hh morphogen gradient, genetic analyses suggest craniofacial development does not completely fit this paradigm. Using the mouse model (Mus musculus), we demonstrated that rather than being driven by a Hh threshold, robust Gli3 transcriptional activity during skeletal and glossal development required interaction with the basic helix-loop-helix TF Hand2. Not only did genetic and expression data support a co-factorial relationship, but genomic analysis revealed that Gli3 and Hand2 were enriched at regulatory elements for genes essential for mandibular patterning and development. Interestingly, motif analysis at sites co-occupied by Gli3 and Hand2 uncovered mandibular-specific, low-affinity, 'divergent' Gli-binding motifs (dGBMs). Functional validation revealed these dGBMs conveyed synergistic activation of Gli targets essential for mandibular patterning and development. In summary, this work elucidates a novel, sequence-dependent mechanism for Gli transcriptional activity within the craniofacial complex that is independent of a graded Hh signal.

A compendium of DNA-binding specificities of transcription factors in Pseudomonas syringae

Pseudomonas syringae is a Gram-negative and model pathogenic bacterium that causes plant diseases worldwide. Here, we set out to identify binding motifs for all 301 annotated transcription factors (TFs) of P. syringae using HT-SELEX. We successfully identify binding motifs for 100 TFs. We map functional interactions between the TFs and their targets in virulence-associated pathways, and validate many of these interactions and functions using additional methods such as ChIP-seq, electrophoretic mobility shift assay (EMSA), RT-qPCR, and reporter assays. Our work identifies 25 virulence-associated master regulators, 14 of which had not been characterized as TFs before.

The authors set out to identify binding motifs for all 301 transcription factors of a plant pathogenic bacterium, Pseudomonas syringae, using HT-SELEX. They successfully identify binding motifs for 100 transcription factors, infer their binding sites on the genome, and validate the predicted interactions and functions.

Combinatorial Action of Temporally Segregated Transcription Factors

Combinatorial action of transcription factors (TFs) with partially overlapping expression is a widespread strategy to generate novel gene-expression patterns and, thus, cellular diversity. Known mechanisms underlying combinatorial activity require co-expression of TFs within the same cell. Here, we describe the mechanism by which two TFs that are never co-expressed generate a new, intersectional expression pattern in C. elegans embryos: lineage-specific priming of a gene by a transiently expressed TF generates a unique intersection with a second TF acting on the same gene four cell divisions later; the second TF is expressed in multiple cells but only activates transcription in those where priming occurred. Early induction of active transcription is necessary and sufficient to establish a competent state, maintained by broadly expressed regulators in the absence of the initial trigger. We uncover additional cells diversified through this mechanism. Our findings define a mechanism for combinatorial TF activity with important implications for generation of cell-type diversity.

•

Lineage-specific priming enables asymmetric gene expression in L/R neuron pairs

•

Transient, lineage-specific TFs prime a locus for later activation by a bilateral TF

•

An early active transcriptional state is necessary and sufficient for priming

•

Maintenance of asymmetric primed state occurs in a symmetric regulatory environment

How affinity of the ELT-2 GATA factor binding to cis-acting regulatory sites controls Caenorhabditis elegans intestinal gene transcription

We define a quantitative relationship between the affinity with which the intestine-specific GATA factor ELT-2 binds to cis-acting regulatory motifs and the resulting transcription of asp-1, a target gene representative of genes involved in Caenorhabditis elegans intestine differentiation. By establishing an experimental system that allows unknown parameters (e.g. the influence of chromatin) to effectively cancel out, we show that levels of asp-1 transcripts increase monotonically with increasing binding affinity of ELT-2 to variant promoter TGATAA sites. The shape of the response curve reveals that the product of the unbound ELT-2 concentration in vivo [i.e. (ELT-2free) or ELT-2 'activity'] and the largest ELT-XXTGATAAXX association constant (Kmax) lies between five and ten. We suggest that this (unitless) product [Kmax×(ELT-2free) or the equivalent product for any other transcription factor] provides an important quantitative descriptor of transcription-factor/regulatory-motif interaction in development, evolution and genetic disease. A more complicated model than simple binding affinity is necessary to explain the fact that ELT-2 appears to discriminate in vivo against equal-affinity binding sites that contain AGATAA instead of TGATAA.

The geometric influence on the Cys2His2 zinc finger domain and functional plasticity

The Cys₂His₂ zinc finger is the most common DNA-binding domain expanding in metazoans since the fungi human split. A proposed catalyst for this expansion is an arms race to silence transposable elements yet it remains poorly understood how this domain is able to evolve the required specificities. Likewise, models of its DNA binding specificity remain error prone due to a lack of understanding of how adjacent fingers influence each other's binding specificity. Here, we use a synthetic approach to exhaustively investigate binding geometry, one of the dominant influences on adjacent finger function. By screening over 28 billion protein–DNA interactions in various geometric contexts we find the plasticity of the most common natural geometry enables more functional amino acid combinations across all targets. Further, residues that define this geometry are enriched in genomes where zinc fingers are prevalent and specificity transitions would be limited in alternative geometries. Finally, these results demonstrate an exhaustive synthetic screen can produce an accurate model of domain function while providing mechanistic insight that may have assisted in the domains expansion.

Loss of an H3K9me anchor rescues laminopathy-linked changes in nuclear organization and muscle function in an Emery-Dreifuss muscular dystrophy model

In this study, Harr et al. use C. elegans to investigate the consequences of a missense mutation (Y45C) in lamin A (encoded by LMNA) found in the human Emery-Dreifuss muscular dystrophy (EDMD) syndrome. Using muscle-specific emerin Dam-ID and other in vivo approaches, the authors report that they were able to counteract the dominant muscle-specific defects provoked by LMNA mutation by the ablation of a lamin-associated H3K9me anchor, suggesting a novel therapeutic pathway for treating EDMD.

Mutations in the nuclear structural protein lamin A produce rare, tissue-specific diseases called laminopathies. The introduction of a human Emery-Dreifuss muscular dystrophy (EDMD)-inducing mutation into the C. elegans lamin (LMN-Y59C), recapitulates many muscular dystrophy phenotypes, and correlates with hyper-sequestration of a heterochromatic array at the nuclear periphery in muscle cells. Using muscle-specific emerin Dam-ID in worms, we monitored the effects of the mutation on endogenous chromatin. An increased contact with the nuclear periphery along chromosome arms, and an enhanced release of chromosomal centers, coincided with the disease phenotypes of reduced locomotion and compromised sarcomere integrity. The coupling of the LMN-Y59C mutation with the ablation of CEC-4, a chromodomain protein that anchors H3K9-methylated chromatin at the nuclear envelope (NE), suppressed the muscle-associated disease phenotypes. Deletion of cec-4 also rescued LMN-Y59C-linked alterations in chromatin organization and some changes in transcription. Sequences that changed position in the LMN-Y59C mutant, are enriched for E2F (EFL-2)-binding sites, consistent with previous studies suggesting that altered Rb-E2F interaction with lamin A may contribute to muscle dysfunction. In summary, we were able to counteract the dominant muscle-specific defects provoked by LMNA mutation by the ablation of a lamin-associated H3K9me anchor, suggesting a novel therapeutic pathway for EDMD.

Redox-mediated regulation of aging and healthspan by an evolutionarily conserved transcription factor HLH-2/Tcf3/E2A

Physiological aging is a complex process, influenced by a plethora of genetic and environmental factors. While being far from fully understood, a number of common aging hallmarks have been elucidated in recent years. Among these, transcriptomic alterations are hypothesized to represent a crucial early manifestation of aging. Accordingly, several transcription factors (TFs) have previously been identified as important modulators of lifespan in evolutionarily distant model organisms. Based on a set of TFs conserved between nematodes, zebrafish, mice, and humans, we here perform a RNA interference (RNAi) screen in C. elegans to discover evolutionarily conserved TFs impacting aging. We identify a basic helix-loop-helix TF, named HLH-2 in nematodes (Tcf3/E2A in mammals), to exert a pronounced lifespan-extending effect in C. elegans upon impairment. We further show that its impairment impacts cellular energy metabolism, increases parameters of healthy aging, and extends nematodal lifespan in a ROS-dependent manner. We then identify arginine kinases, orthologues of mammalian creatine kinases, as a target of HLH-2 transcriptional regulation, serving to mediate the healthspan-promoting effects observed upon impairment of hlh-2 expression. Consistently, HLH-2 is shown to epistatically interact with core components of known lifespan-regulating pathways, i.e. AAK-2/AMPK and LET-363/mTOR, as well as the aging-related TFs SKN-1/Nrf2 and HSF-1. Lastly, single-nucelotide polymorphisms (SNPs) in Tcf3/E2A are associated with exceptional longevity in humans. Together, these findings demonstrate that HLH-2 regulates energy metabolism via arginine kinases and thereby affects the aging phenotype dependent on ROS-signaling and established canonical effectors.

•

A C. elegans RNAi screen identifies conserved aging-related transcription factors.

•

Impairment of transcription factor hlh-2 has the most pronounced effect on lifespan.

•

C. elegans HLH-2 affects cellular energy homeostasis and redox signaling.

•

HLH-2 modulates arginine kinase to interact with downstream longevity pathways.

•

Polymorphisms (SNPs) in the hlh-2 orthologue Tcf3/E2A are linked to human longevity.

Brn3/POU-IV-type POU homeobox genes-Paradigmatic regulators of neuronal identity across phylogeny

One approach to understand the construction of complex systems is to investigate whether there are simple design principles that are commonly used in building such a system. In the context of nervous system development, one may ask whether the generation of its highly diverse sets of constituents, that is, distinct neuronal cell types, relies on genetic mechanisms that share specific common features. Specifically, are there common patterns in the function of regulatory genes across different neuron types and are those regulatory mechanisms not only used in different parts of one nervous system, but are they conserved across animal phylogeny? We address these questions here by focusing on one specific, highly conserved and well-studied regulatory factor, the POU homeodomain transcription factor UNC-86. Work over the last 30 years has revealed a common and paradigmatic theme of unc-86 function throughout most of the neuron types in which Caenorhabditis elegans unc-86 is expressed. Apart from its role in preventing lineage reiterations during development, UNC-86 operates in combination with distinct partner proteins to initiate and maintain terminal differentiation programs, by coregulating a vast array of functionally distinct identity determinants of specific neuron types. Mouse orthologs of unc-86, the Brn3 genes, have been shown to fulfill a similar function in initiating and maintaining neuronal identity in specific parts of the mouse brain and similar functions appear to be carried out by the sole Drosophila ortholog, Acj6. The terminal selector function of UNC-86 in many different neuron types provides a paradigm for neuronal identity regulation across phylogeny. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Invertebrate Organogenesis > Worms Nervous System Development > Vertebrates: Regional Development.

Diversification of the Caenorhabditis heat shock response by Helitron transposable elements

Heat Shock Factor 1 (HSF-1) is a key regulator of the heat shock response (HSR). Upon heat shock, HSF-1 binds well-conserved motifs, called Heat Shock Elements (HSEs), and drives expression of genes important for cellular protection during this stress. Remarkably, we found that substantial numbers of HSEs in multiple Caenorhabditis species reside within Helitrons, a type of DNA transposon. Consistent with Helitron-embedded HSEs being functional, upon heat shock they display increased HSF-1 and RNA polymerase II occupancy and up-regulation of nearby genes in C. elegans. Interestingly, we found that different genes appear to be incorporated into the HSR by species-specific Helitron insertions in C. elegans and C. briggsae and by strain-specific insertions among different wild isolates of C. elegans. Our studies uncover previously unidentified targets of HSF-1 and show that Helitron insertions are responsible for rewiring and diversifying the Caenorhabditis HSR.

Gene expression modulation by the linker of nucleoskeleton and cytoskeleton complex contributes to proteostasis

Cellular mechanisms that act in concert to maintain protein homeostasis (proteostasis) are vital for organismal functionality and survival. Nevertheless, subsets of aggregation-prone proteins form toxic aggregates (proteotoxicity) that in some cases, underlie the development of neurodegenerative diseases. Proteotoxic aggregates are often deposited in the vicinity of the nucleus, a process that is cytoskeleton-dependent. Accordingly, cytoskeletal dysfunction contributes to pathological hallmarks of various neurodegenerative diseases. Here, we asked whether the linker of nucleoskeleton and cytoskeleton (LINC) complex, which bridges these filaments across the nuclear envelope, is needed for the maintenance of proteostasis. Employing model nematodes, we discovered that knocking down LINC components impairs the ability of the worm to cope with proteotoxicity. Knocking down anc-1, which encodes a key component of the LINC complex, modulates the expression of transcription factors and E3 ubiquitin ligases, thereby affecting the rates of protein ubiquitination and impairing proteasome-mediated protein degradation. Our results establish a link between the LINC complex, protein degradation, and neurodegeneration-associated proteotoxicity.

A Myt1 family transcription factor defines neuronal fate by repressing non-neuronal genes

Cellular differentiation requires both activation of target cell transcriptional programs and repression of non-target cell programs. The Myt1 family of zinc finger transcription factors contributes to fibroblast to neuron reprogramming in vitro. Here, we show that ztf-11 (Zinc-finger Transcription Factor-11), the sole Caenorhabditis elegans Myt1 homolog, is required for neurogenesis in multiple neuronal lineages from previously differentiated epithelial cells, including a neuron generated by a developmental epithelial-to-neuronal transdifferentiation event. ztf-11 is exclusively expressed in all neuronal precursors with remarkable specificity at single-cell resolution. Loss of ztf-11 leads to upregulation of non-neuronal genes and reduced neurogenesis. Ectopic expression of ztf-11 in epidermal lineages is sufficient to produce additional neurons. ZTF-11 functions together with the MuvB corepressor complex to suppress the activation of non-neuronal genes in neurons. These results dovetail with the ability of Myt1l (Myt1-like) to drive neuronal transdifferentiation in vitro in vertebrate systems. Together, we identified an evolutionarily conserved mechanism to specify neuronal cell fate by repressing non-neuronal genes.

Self-sperm induce resistance to the detrimental effects of sexual encounters with males in hermaphroditic nematodes

Sexual interactions have a potent influence on health in several species, including mammals. Previous work in C. elegans identified strategies used by males to accelerate the demise of the opposite sex (hermaphrodites). But whether hermaphrodites evolved counter-strategies against males remains unknown. Here we discover that young C. elegans hermaphrodites are remarkably resistant to brief sexual encounters with males, whereas older hermaphrodites succumb prematurely. Surprisingly, it is not their youthfulness that protects young hermaphrodites, but the fact that they have self-sperm. The beneficial effect of self-sperm is mediated by a sperm-sensing pathway acting on the soma rather than by fertilization. Activation of this pathway in females triggers protection from the negative impact of males. Interestingly, the role of self-sperm in protecting against the detrimental effects of males evolved independently in hermaphroditic nematodes. Endogenous strategies to delay the negative effect of mating may represent a key evolutionary innovation to maximize reproductive success.

A nematode worm known as Caenorhabditis elegans is often used in the laboratory to study how animals grow and develop. There are two types of C. elegans worm: hermaphrodite individuals produce both female sex cells (eggs) and male sex cells (sperm), while male individuals only produce sperm.

The hermaphrodite worms are able to reproduce without mating with another worm, allowing populations of C. elegans to grow rapidly when they are living in favorable conditions. However, when the hermaphrodites do mate with males they tend to produce more offspring. These offspring are also usually healthier because they receive a mixture of genetic material from two different parents.

Although mating is beneficial for the survival of a species it can also harm an individual animal. Previous studies have shown that mating with male worms can accelerate aging of hermaphrodite worms and cause premature death. However, it remained unclear whether hermaphrodite worms have evolved any mechanisms to protect themselves after mating with a male.

To address this question, Booth et al. used genetic techniques to study the lifespans of hermaphrodite worms. The experiments found that the hermaphrodites’ own sperm (known as self-sperm) regulated a sperm-sensing signaling pathway that protected them from the negative impact of mating with males. Hermaphrodites with self-sperm that mated with males lived for a similar length of time as hermaphrodites that did not mate. On the other hand, hermaphrodites that did not have self-sperm (because they were older or had a genetic mutation) had shorter lifespans after mating than worms that did not mate. Modulating the sperm-sensing signaling pathway in worms that lacked self-sperm was sufficient to protect them from the negative effects of mating with males.

Further experiments found that the hermaphrodites of another nematode worm called C. briggsae – which evolved self-sperm independently of C. elegans – also protected themselves from the negative effects of mating with males in a similar way. This suggests that other animals may also have evolved similar mechanisms to protect themselves from harm when mating.

A separate study by Shi et al. has found that the beneficial effects of self-sperm are mediated by a pathway linked to longevity that also exists in mammals. The results of both investigations combined suggest possible avenues for future research into the complex relationship between health, longevity, and reproduction.

Similarity regression predicts evolution of transcription factor sequence specificity

Transcription factor (TF) binding specificities (motifs) are essential for the analysis of gene regulation. Accurate prediction of TF motifs is critical, because it is infeasible to assay all TFs in all sequenced eukaryotic genomes. There is ongoing controversy regarding the degree of motif diversification among related species that is, in part, because of uncertainty in motif prediction methods. Here we describe similarity regression, a significantly improved method for predicting motifs, which we use to update and expand the Cis-BP database. Similarity regression inherently quantifies TF motif evolution, and shows that previous claims of near-complete conservation of motifs between human and Drosophila are inflated, with nearly half of the motifs in each species absent from the other, largely due to extensive divergence in C2H2 zinc finger proteins. We conclude that diversification in DNA-binding motifs is pervasive, and present a new tool and updated resource to study TF diversity and gene regulation across eukaryotes.

Active chromatin marks drive spatial sequestration of heterochromatin in C. elegans nuclei

The execution of developmental programs of gene expression requires an accurate partitioning of the genome into subnuclear compartments, with active euchromatin enriched centrally and silent heterochromatin at the nuclear periphery¹. The existence of degenerative diseases linked to lamin A mutations suggests that perinuclear binding of chromatin contributes to cell-type integrity^2,3. The methylation of lysine 9 of histone H3 (H3K9me) characterizes heterochromatin and mediates both transcriptional repression and chromatin anchoring at the inner nuclear membrane⁴. In Caenorhabditis elegans embryos, chromodomain protein CEC-4 bound to the inner nuclear membrane tethers heterochromatin through H3K9me^3,5, whereas in differentiated tissues, a second heterochromatin-sequestering pathway is induced. Here we use an RNA interference screen in the cec-4 background and identify MRG-1 as a broadly expressed factor that is necessary for this second chromatin anchor in intestinal cells. However, MRG-1 is exclusively bound to euchromatin, suggesting that it acts indirectly. Heterochromatin detachment in double mrg-1; cec-4 mutants is rescued by depleting the histone acetyltransferase CBP-1/p300 or the transcription factor ATF-8, a member of the bZIP family (which is known to recruit CBP/p300). Overexpression of CBP-1 in cec-4 mutants is sufficient to delocalize heterochromatin in an ATF-8-dependent manner. CBP-1 and H3K27ac levels increase in heterochromatin upon mrg-1 knockdown, coincident with delocalization. This suggests that the spatial organization of chromatin in C. elegans is regulated both by the direct perinuclear attachment of silent chromatin, and by an active retention of CBP-1/p300 in euchromatin. The two pathways contribute differentially in embryos and larval tissues, with CBP-1 sequestration by MRG-1 having a major role in differentiated cells.

Spatiotemporal Gene Expression Analysis of the Caenorhabditis elegans Germline Uncovers a Syncytial Expression Switch

Developmental programs are executed by tightly controlled gene regulatory pathways. Here, we combined the unique sample retrieval capacity afforded by laser capture microscopy with analysis of mRNA abundance by CEL-Seq (cell expression by linear amplification and sequencing) to generate a spatiotemporal gene expression map of the Caenorhabditis elegans syncytial germline from adult hermaphrodites and males. We found that over 6000 genes exhibit spatiotemporally dynamic expression patterns throughout the hermaphrodite germline, with two dominant groups of genes exhibiting reciprocal shifts in expression at late pachytene during meiotic prophase I. We found a strong correlation between restricted spatiotemporal expression and known developmental and cellular processes, indicating that these gene expression changes may be an important driver of germ cell progression. Analysis of the male gonad revealed a shift in gene expression at early pachytene and upregulation of subsets of genes following the meiotic divisions, specifically in early and late spermatids, mostly transcribed from the X chromosome. We observed that while the X chromosome is silenced throughout the first half of the gonad, some genes escape this control and are highly expressed throughout the germline. Although we found a strong correlation between the expression of genes corresponding to CSR-1-interacting 22G-RNAs during germ cell progression, we also found that a large fraction of genes may bypass the need for CSR-1-mediated germline licensing. Taken together, these findings suggest the existence of mechanisms that enable a shift in gene expression during prophase I to promote germ cell progression.

Non-proteolytic activity of 19S proteasome subunit RPT-6 regulates GATA transcription during response to infection

GATA transcription factors play a crucial role in the regulation of immune functions across metazoans. In Caenorhabditis elegans, the GATA transcription factor ELT-2 is involved in the control of not only infections but also recovery after an infection. We identified RPT-6, part of the 19S proteasome subunit, as an ELT-2 binding partner that is required for the proper expression of genes required for both immunity against bacterial infections and recovery after infection. We found that the intact ATPase domain of RPT-6 is required for the interaction and that inhibition of rpt-6 affected the expression of ELT-2-controlled genes, preventing the appropriate immune response against Pseudomonas aeruginosa and recovery from infection by the pathogen. Further studies indicated that SKN-1, which is an Nrf transcription factor involved in the response to oxidative stress and infection, is activated by inhibition of rpt-6. Our results indicate that RPT-6 interacts with ELT-2 in vivo to control the expression of immune genes in a manner that is likely independent of the proteolytic activity of the proteasome.

The conserved GATA transcription factor ELT-2 plays an important role in the control of genes required for both defense and recovery from infection. We show that RPT-6, a component of the 19S subunit, physically interacts with ELT-2 in vivo, controlling the expression of ELT-2-dependent genes and the response of the nematode Caenorhabditis elegans to bacterial infection. The proteolytic activity of the proteasome has surfaced as a key regulator of gene expression, but our results provide evidence indicating that a non-canonical activity of the 26S proteasome subunit plays an important role in the control of gene expression during the response to bacterial infection.

FACT Sets a Barrier for Cell Fate Reprogramming in Caenorhabditis elegans and Human Cells

The chromatin regulator FACT (facilitates chromatin transcription) is essential for ensuring stable gene expression by promoting transcription. In a genetic screen using Caenorhabditis elegans, we identified that FACT maintains cell identities and acts as a barrier for transcription factor-mediated cell fate reprogramming. Strikingly, FACT's role as a barrier to cell fate conversion is conserved in humans as we show that FACT depletion enhances reprogramming of fibroblasts. Such activity is unexpected because FACT is known as a positive regulator of gene expression, and previously described reprogramming barriers typically repress gene expression. While FACT depletion in human fibroblasts results in decreased expression of many genes, a number of FACT-occupied genes, including reprogramming-promoting factors, show increased expression upon FACT depletion, suggesting a repressive function of FACT. Our findings identify FACT as a cellular reprogramming barrier in C. elegans and humans, revealing an evolutionarily conserved mechanism for cell fate protection.

BEESEM: estimation of binding energy models using HT-SELEX data

Motivation

Characterizing the binding specificities of transcription factors (TFs) is crucial to the study of gene expression regulation. Recently developed high-throughput experimental methods, including protein binding microarrays (PBM) and high-throughput SELEX (HT-SELEX), have enabled rapid measurements of the specificities for hundreds of TFs. However, few studies have developed efficient algorithms for estimating binding motifs based on HT-SELEX data. Also the simple method of constructing a position weight matrix (PWM) by comparing the frequency of the preferred sequence with single-nucleotide variants has the risk of generating motifs with higher information content than the true binding specificity.

Results

We developed an algorithm called BEESEM that builds on a comprehensive biophysical model of protein-DNA interactions, which is trained using the expectation maximization method. BEESEM is capable of selecting the optimal motif length and calculating the confidence intervals of estimated parameters. By comparing BEESEM with the published motifs estimated using the same HT-SELEX data, we demonstrate that BEESEM provides significant improvements. We also evaluate several motif discovery algorithms on independent PBM and ChIP-seq data. BEESEM provides significantly better fits to in vitro data, but its performance is similar to some other methods on in vivo data under the criterion of the area under the receiver operating characteristic curve (AUROC). This highlights the limitations of the purely rank-based AUROC criterion. Using quantitative binding data to assess models, however, demonstrates that BEESEM improves on prior models.

Availability and Implementation

Freely available on the web at http://stormo.wustl.edu/resources.html .

Contact

stormo@wustl.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.