Regulatory Logic of Pan-Neuronal Gene Expression in C. elegans

The Iroquois (Iro/Irx) homeobox genes are conserved Hox targets involved in motor neuron development

The Iroquois (Iro/Irx) homeobox genes encode transcription factors with fundamental roles in animal development. Despite their link to various congenital conditions in humans, our understanding of Iro/Irx gene expression, function, and regulation remains incomplete. Here, we conducted a systematic expression analysis of all six mouse Irx genes in the embryonic spinal cord. We found five Irx genes (Irx1, Irx2, Irx3, Irx5, and Irx6) to be confined mostly to ventral spinal domains, offering new molecular markers for specific groups of post-mitotic motor neurons (MNs). Further, we engineered Irx2, Irx5, and Irx6 mouse mutants and uncovered essential but distinct roles for Irx2 and Irx6 in MN development. Last, we found that the highly conserved regulators of MN development across species, the HOX proteins, directly control Irx gene expression both in mouse and C. elegans MNs, critically expanding the repertoire of HOX target genes in the developing nervous system. Altogether, our study provides important insights into Iro/Irx expression and function in the developing spinal cord, and uncovers an ancient gene regulatory relationship between HOX and Iro/Irx genes.

Deep Transcriptomics Reveals Cell-Specific Isoforms of Pan-Neuronal Genes

Profiling gene expression in single neurons using single-cell RNA-Seq is a powerful method for understanding the molecular diversity of the nervous system. Profiling alternative splicing in single neurons using these methods is more challenging, however, due to low capture efficiency and sensitivity. As a result, we know much less about splicing patterns and regulation across neurons than we do about gene expression. Here we leverage unique attributes of the C. elegans nervous system to investigate deep cell-specific transcriptomes complete with biological replicates generated by the CeNGEN consortium, enabling high-confidence assessment of splicing across neuron types even for lowly-expressed genes. Global splicing maps reveal several striking observations, including pan-neuronal genes that harbor cell-specific splice variants, abundant differential intron retention across neuron types, and a single neuron highly enriched for upstream alternative 3' splice sites. We develop an algorithm to identify unique cell-specific expression patterns and use it to discover both cell-specific isoforms and potential regulatory RNA binding proteins that establish these isoforms. Genetic interrogation of these RNA binding proteins in vivo identifies three distinct regulatory factors employed to establish unique splicing patterns in a single neuron. Finally, we develop a user-friendly platform for spatial transcriptomic visualization of these splicing patterns with single-neuron resolution.

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generate a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database. Single-cell RNA sequencing of 13,200 cells reveals that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. Extending C. elegans Neuronal Gene Expression Map and Network (CeNGEN) findings, all MN subclasses are delineated by distinct expression codes of either neuropeptide or transcription factor gene families. Strikingly, combinatorial codes of homeodomain transcription factor genes succinctly delineate adult MN diversity in both C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs and uncovers organizing principles and conserved molecular codes of adult MN diversity.

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.

Homeodomain proteins hierarchically specify neuronal diversity and synaptic connectivity

How our brain generates diverse neuron types that assemble into precise neural circuits remains unclear. Using Drosophila lamina neuron types (L1-L5), we show that the primary homeodomain transcription factor (HDTF) brain-specific homeobox (Bsh) is initiated in progenitors and maintained in L4/L5 neurons to adulthood. Bsh activates secondary HDTFs Ap (L4) and Pdm3 (L5) and specifies L4/L5 neuronal fates while repressing the HDTF Zfh1 to prevent ectopic L1/L3 fates (control: L1-L5; Bsh-knockdown: L1-L3), thereby generating lamina neuronal diversity for normal visual sensitivity. Subsequently, in L4 neurons, Bsh and Ap function in a feed-forward loop to activate the synapse recognition molecule DIP-β, thereby bridging neuronal fate decision to synaptic connectivity. Expression of a Bsh:Dam, specifically in L4, reveals Bsh binding to the DIP-β locus and additional candidate L4 functional identity genes. We propose that HDTFs function hierarchically to coordinate neuronal molecular identity, circuit formation, and function. Hierarchical HDTFs may represent a conserved mechanism for linking neuronal diversity to circuit assembly and function.

Nuclear factor Y is a pervasive regulator of neuronal gene expression

Nervous system function relies on the establishment of complex gene expression programs that provide neuron-type-specific and core pan-neuronal features. These complementary regulatory paradigms are controlled by terminal selector and parallel-acting transcription factors (TFs), respectively. Here, we identify the nuclear factor Y (NF-Y) TF as a pervasive direct and indirect regulator of both neuron-type-specific and pan-neuronal gene expression. Mapping global NF-Y targets reveals direct binding to the cis-regulatory regions of pan-neuronal genes and terminal selector TFs. We show that NFYA-1 controls pan-neuronal gene expression directly through binding to CCAAT boxes in target gene promoters and indirectly by regulating the expression of terminal selector TFs. Further, we find that NFYA-1 regulation of neuronal gene expression is important for neuronal activity and motor function. Thus, our research sheds light on how global neuronal gene expression programs are buffered through direct and indirect regulatory mechanisms.

Expression and function of Caenorhabditis elegans UNCP-18, a paralog of the SM protein UNC-18

Sec1/Munc18 (SM) proteins are important regulators of SNARE complex assembly during exocytosis throughout all major animal tissue types. However, expression of a founding member of the SM family, UNC-18, is mostly restricted to the nervous system of the nematode Caenorhabditis elegans, where it is important for synaptic transmission. Moreover, unc-18 null mutants do not display the lethality phenotype associated with (a) loss of all Drosophila and mouse orthologs of unc-18 and (b) with complete elimination of synaptic transmission in C. elegans. We investigated whether a previously uncharacterized unc-18 paralog, which we named uncp-18, may be able to explain the restricted expression and limited phenotypes of unc-18 null mutants. A reporter allele shows ubiquitous expression of uncp-18. Analysis of uncp-18 null mutants, unc-18 and uncp-18 double null mutants, as well as overexpression of uncp-18 in an unc-18 null mutant background, shows that these 2 genes can functionally compensate for one another and are redundantly required for embryonic viability. Our results indicate that the synaptic transmission defects of unc-18 null mutants cannot necessarily be interpreted as constituting a null phenotype for SM protein function at the synapse.

Distinct neurexin isoforms cooperate to initiate and maintain foraging activity

Neurexins are synaptic adhesion molecules that play diverse roles in synaptic development, function, maintenance, and plasticity. Neurexin genes have been associated with changes in human behavior, where variants in NRXN1 are associated with autism, schizophrenia, and Tourette syndrome. While NRXN1, NRXN2, and NRXN3 all encode major α and β isoforms, NRXN1 uniquely encodes a γ isoform, for which mechanistic roles in behavior have yet to be defined. Here, we show that both α and γ isoforms of neurexin/nrx-1 are required for the C. elegans behavioral response to food deprivation, a sustained period of hyperactivity upon food loss. We find that the γ isoform regulates initiation and the α isoform regulates maintenance of the behavioral response to food deprivation, demonstrating cooperative function of multiple nrx-1 isoforms in regulating a sustained behavior. The γ isoform alters monoamine signaling via octopamine, relies on specific expression of NRX-1 isoforms throughout the relevant circuit, and is independent of neuroligin/nlg-1, the canonical trans-synaptic partner of nrx-1. The α isoform regulates the pre-synaptic structure of the octopamine producing RIC neuron and its maintenance role is conditional on neuroligin/nlg-1. Collectively, these results demonstrate that neurexin isoforms can have separate behavioral roles and act cooperatively across neuronal circuits to modify behavior, highlighting the need to directly analyze and consider all isoforms when defining the contribution of neurexins to behavior.

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.

Electric shock causes a fleeing-like persistent behavioral response in the nematode Caenorhabditis elegans

Behavioral persistency reflects internal brain states, which are the foundations of multiple brain functions. However, experimental paradigms enabling genetic analyses of behavioral persistency and its associated brain functions have been limited. Here, we report novel persistent behavioral responses caused by electric stimuli in the nematode Caenorhabditis elegans. When the animals on bacterial food are stimulated by alternating current, their movement speed suddenly increases 2- to 3-fold, persisting for more than 1 minute even after a 5-second stimulation. Genetic analyses reveal that voltage-gated channels in the neurons are required for the response, possibly as the sensors, and neuropeptide signaling regulates the duration of the persistent response. Additional behavioral analyses implicate that the animal's response to electric shock is scalable and has a negative valence. These properties, along with persistence, have been recently regarded as essential features of emotion, suggesting that C. elegans response to electric shock may reflect a form of emotion, akin to fear.

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generated a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database (http://celegans.spinalcordatlas.org). Single-cell RNA-sequencing of 13,200 cells revealed that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. All subclasses are delineated by unique expression codes of either neuropeptide or transcription factor gene families. Strikingly, we found that combinatorial codes of homeodomain transcription factor genes define adult MN diversity both in C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs, and uncovers organizing principles and conserved molecular codes of adult MN diversity.

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.

Identification and characterization of transcribed enhancers during cerebellar development through enhancer RNA analysis

BACKGROUND

The development of the brain requires precise coordination of molecular processes across many cell-types. Underpinning these events are gene expression programs which require intricate regulation by non-coding regulatory sequences known as enhancers. In the context of the developing brain, transcribed enhancers (TEs) regulate temporally-specific expression of genes critical for cell identity and differentiation. Transcription of non-coding RNAs at active enhancer sequences, known as enhancer RNAs (eRNAs), is tightly associated with enhancer activity and has been correlated with target gene expression. TEs have been characterized in a multitude of developing tissues, however their regulatory role has yet to be described in the context of embryonic and early postnatal brain development. In this study, eRNA transcription was analyzed to identify TEs active during cerebellar development, as a proxy for the developing brain. Cap Analysis of Gene Expression followed by sequencing (CAGE-seq) was conducted at 12 stages throughout embryonic and early postnatal cerebellar development.

RESULTS

Temporal analysis of eRNA transcription identified clusters of TEs that peak in activity during either embryonic or postnatal times, highlighting their importance for temporally specific developmental events. Functional analysis of putative target genes identified molecular mechanisms under TE regulation revealing that TEs regulate genes involved in biological processes specific to neurons. We validate enhancer activity using in situ hybridization of eRNA expression from TEs predicted to regulate Nfib, a gene critical for cerebellar granule cell differentiation.

CONCLUSIONS

The results of this analysis provide a valuable dataset for the identification of cerebellar enhancers and provide insight into the molecular mechanisms critical for brain development under TE regulation. This dataset is shared with the community through an online resource ( https://goldowitzlab.shinyapps.io/trans-enh-app/ ).

Ancestral glycoprotein hormone-receptor pathway controls growth in C. elegans

In vertebrates, thyrostimulin is a highly conserved glycoprotein hormone that, besides thyroid stimulating hormone (TSH), is a potent ligand of the TSH receptor. Thyrostimulin is considered the most ancestral glycoprotein hormone and orthologs of its subunits, GPA2 and GPB5, are widely conserved across vertebrate and invertebrate animals. Unlike TSH, however, the functions of the thyrostimulin neuroendocrine system remain largely unexplored. Here, we identify a functional thyrostimulin-like signaling system in Caenorhabditis elegans. We show that orthologs of GPA2 and GPB5, together with thyrotropin-releasing hormone (TRH) related neuropeptides, constitute a neuroendocrine pathway that promotes growth in C. elegans. GPA2/GPB5 signaling is required for normal body size and acts through activation of the glycoprotein hormone receptor ortholog FSHR-1. C. elegans GPA2 and GPB5 increase cAMP signaling by FSHR-1 in vitro. Both subunits are expressed in enteric neurons and promote growth by signaling to their receptor in glial cells and the intestine. Impaired GPA2/GPB5 signaling causes bloating of the intestinal lumen. In addition, mutants lacking thyrostimulin-like signaling show an increased defecation cycle period. Our study suggests that the thyrostimulin GPA2/GPB5 pathway is an ancient enteric neuroendocrine system that regulates intestinal function in ecdysozoans, and may ancestrally have been involved in the control of organismal growth.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

C. elegans Gene Regulatory Alleles and Reporter Bashing Studies

Gene regulation has been studied in C. elegans for over 30 years. In this analysis of 102 publications, we find that most transcriptional cis-regulatory elements are located within 5,000 bp of the transcription start site. Over 75% of studies conclude that transcriptional elements and 5'UTRs activate-, while 3'UTRs repress gene expression. While gene regulatory mutations make up less than 0.8% of alleles in forward genetics screens, recent CRISPR-Cas approaches are increasing the number of tested mutations. This work provides a resource of known gene regulatory sequences in C.elegans .

Cell-specific effects of the sole C. elegans Daughterless/E protein homolog, HLH-2, on nervous system development

Are there common mechanisms of neurogenesis used throughout an entire nervous system? We explored to what extent canonical proneural class I/II bHLH complexes are responsible for neurogenesis throughout the entire Caenorhabditis elegans nervous system. Distinct, lineage-specific proneural class II bHLH factors are generally thought to operate via interaction with a common, class I bHLH subunit, encoded by Daughterless in flies, the E proteins in vertebrates and HLH-2 in C. elegans. To eliminate function of all proneuronal class I/II bHLH complexes, we therefore genetically removed maternal and zygotic hlh-2 gene activity. We observed broad effects on neurogenesis, but still detected normal neurogenesis in many distinct neuron-producing lineages of the central and peripheral nervous system. Moreover, we found that hlh-2 selectively affects some aspects of neuron differentiation while leaving others unaffected. Although our studies confirm the function of proneuronal class I/II bHLH complexes in many different lineages throughout a nervous system, we conclude that their function is not universal, but rather restricted by lineage, cell type and components of differentiation programs affected.

Coordinated control of neuronal differentiation and wiring by sustained transcription factors

The large diversity of cell types in nervous systems presents a challenge in identifying the genetic mechanisms that encode it. Here, we report that nearly 200 distinct neurons in the Drosophila visual system can each be defined by unique combinations of on average 10 continuously expressed transcription factors. We show that targeted modifications of this terminal selector code induce predictable conversions of neuronal fates that appear morphologically and transcriptionally complete. Cis-regulatory analysis of open chromatin links one of these genes to an upstream patterning factor that specifies neuronal fates in stem cells. Experimentally validated network models describe the synergistic regulation of downstream effectors by terminal selectors and ecdysone signaling during brain wiring. Our results provide a generalizable framework of how specific fates are implemented in postmitotic neurons.

O-GlcNAc transferase plays a non-catalytic role in C. elegans male fertility

Animal behavior is influenced by the competing drives to maintain energy and to reproduce. The balance between these evolutionary pressures and how nutrient signaling pathways intersect with mating remains unclear. The nutrient sensor O-GlcNAc transferase, which post-translationally modifies intracellular proteins with a single monosaccharide, is responsive to cellular nutrient status and regulates diverse biological processes. Though essential in most metazoans, O-GlcNAc transferase (ogt-1) is dispensable in Caenorhabditis elegans, allowing genetic analysis of its physiological roles. Compared to control, ogt-1 males had a four-fold reduction in mean offspring, with nearly two thirds producing zero progeny. Interestingly, we found that ogt-1 males transferred sperm less often, and virgin males had reduced sperm count. ogt-1 males were also less likely to engage in mate-searching and mate-response behaviors. Surprisingly, we found normal fertility for males with hypodermal expression of ogt-1 and for ogt-1 strains with catalytic-dead mutations. This suggests OGT-1 serves a non-catalytic function in the hypodermis impacting male fertility and mating behavior. This study builds upon research on the nutrient sensor O-GlcNAc transferase and demonstrates a role it plays in the interplay between the evolutionary drives for reproduction and survival.

Maintenance of neurotransmitter identity by Hox proteins through a homeostatic mechanism

Hox transcription factors play fundamental roles during early patterning, but they are also expressed continuously, from embryonic stages through adulthood, in the nervous system. However, the functional significance of their sustained expression remains unclear. In C. elegans motor neurons (MNs), we find that LIN-39 (Scr/Dfd/Hox4-5) is continuously required during post-embryonic life to maintain neurotransmitter identity, a core element of neuronal function. LIN-39 acts directly to co-regulate genes that define cholinergic identity (e.g., unc-17/VAChT, cho-1/ChT). We further show that LIN-39, MAB-5 (Antp/Hox6-8) and the transcription factor UNC-3 (Collier/Ebf) operate in a positive feedforward loop to ensure continuous and robust expression of cholinergic identity genes. Finally, we identify a two-component design principle for homeostatic control of Hox gene expression in adult MNs: Hox transcriptional autoregulation is counterbalanced by negative UNC-3 feedback. These findings uncover a noncanonical role for Hox proteins during post-embryonic life, critically broadening their functional repertoire from early patterning to the control of neurotransmitter identity.

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.

Sociosexual behavior requires both activating and repressive roles of Tfap2e/AP-2ε in vomeronasal sensory neurons

Neuronal identity dictates the position in an epithelium, and the ability to detect, process, and transmit specific signals to specified targets. Transcription factors (TFs) determine cellular identity via direct modulation of genetic transcription and recruiting chromatin modifiers. However, our understanding of the mechanisms that define neuronal identity and their magnitude remain a critical barrier to elucidate the etiology of congenital and neurodegenerative disorders. The rodent vomeronasal organ provides a unique system to examine in detail the molecular mechanisms underlying the differentiation and maturation of chemosensory neurons. Here, we demonstrated that the identity of postmitotic/maturing vomeronasal sensory neurons (VSNs), and vomeronasal-dependent behaviors can be reprogrammed through the rescue of Tfap2e/AP-2ε expression in the Tfap2eNull mice, and partially reprogrammed by inducing ectopic Tfap2e expression in mature apical VSNs. We suggest that the TF Tfap2e can reprogram VSNs bypassing cellular plasticity restrictions, and that it directly controls the expression of batteries of vomeronasal genes.

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.

Toward a more accurate 3D atlas of C. elegans neurons

Background

Determining cell identity in volumetric images of tagged neuronal nuclei is an ongoing challenge in contemporary neuroscience. Frequently, cell identity is determined by aligning and matching tags to an “atlas” of labeled neuronal positions and other identifying characteristics. Previous analyses of such C. elegans datasets have been hampered by the limited accuracy of such atlases, especially for neurons present in the ventral nerve cord, and also by time-consuming manual elements of the alignment process.

Results

We present a novel automated alignment method for sparse and incomplete point clouds of the sort resulting from typical C. elegans fluorescence microscopy datasets. This method involves a tunable learning parameter and a kernel that enforces biologically realistic deformation. We also present a pipeline for creating alignment atlases from datasets of the recently developed NeuroPAL transgene. In combination, these advances allow us to label neurons in volumetric images with confidence much higher than previous methods.

Conclusions

We release, to the best of our knowledge, the most complete full-body C. elegans 3D positional neuron atlas, incorporating positional variability derived from at least 7 animals per neuron, for the purposes of cell-type identity prediction for myriad applications (e.g., imaging neuronal activity, gene expression, and cell-fate).

The Q system for conditional gene expression is leaky and lacks dynamic range in C. elegans neurons

The Q system allows for conditional gene expression in several organisms, including C. elegans . We aimed to apply this system in C. elegans neurons to obtain temporally-resolved, tissue-specific expression of a fluorescent reporter. We report that, in our hands, there is undesired expression of the reporter in conditions where expression is supposed to be repressed. In addition, in this setup, the signal‑to‑noise ratio of the Q system is unfavorable. We conclude that the Q system is far from optimal in the C. elegans nervous system, and advise cautious use.

Open Frontiers in Neural Cell Type Investigations; Lessons From Caenorhabditis elegans and Beyond, Toward a Multimodal Integration

Nervous system cells, the building blocks of circuits, have been studied with ever-progressing resolution, yet neural circuits appear still resistant to schemes of reductionist classification. Due to their sheer numbers, complexity and diversity, their systematic study requires concrete classifications that can serve reduced dimensionality, reproducibility, and information integration. Conventional hierarchical schemes transformed through the history of neuroscience by prioritizing criteria of morphology, (electro)physiological activity, molecular content, and circuit function, influenced by prevailing methodologies of the time. Since the molecular biology revolution and the recent advents in transcriptomics, molecular profiling gains ground toward the classification of neurons and glial cell types. Yet, transcriptomics entails technical challenges and more importantly uncovers unforeseen spatiotemporal heterogeneity, in complex and simpler nervous systems. Cells change states dynamically in space and time, in response to stimuli or throughout their developmental trajectory. Mapping cell type and state heterogeneity uncovers uncharted terrains in neurons and especially in glial cell biology, that remains understudied in many aspects. Examining neurons and glial cells from the perspectives of molecular neuroscience, physiology, development and evolution highlights the advantage of multifaceted classification schemes. Among the amalgam of models contributing to neuroscience research, Caenorhabditis elegans combines nervous system anatomy, lineage, connectivity and molecular content, all mapped at single-cell resolution, and can provide valuable insights for the workflow and challenges of the multimodal integration of cell type features. This review reflects on concepts and practices of neuron and glial cells classification and how research, in C. elegans and beyond, guides nervous system experimentation through integrated multidimensional schemes. It highlights underlying principles, emerging themes, and open frontiers in the study of nervous system development, regulatory logic and evolution. It proposes unified platforms to allow integrated annotation of large-scale datasets, gene-function studies, published or unpublished findings and community feedback. Neuroscience is moving fast toward interdisciplinary, high-throughput approaches for combined mapping of the morphology, physiology, connectivity, molecular function, and the integration of information in multifaceted schemes. A closer look in mapped neural circuits and understudied terrains offers insights for the best implementation of these approaches.

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.

Robust regulatory architecture of pan-neuronal gene expression

Pan-neuronally expressed genes, such as genes involved in the synaptic vesicle cycle or in neuropeptide maturation, are critical for proper function of all neurons, but the transcriptional control mechanisms that direct such genes to all neurons of a nervous system remain poorly understood. We show here that six members of the CUT family of homeobox genes control pan-neuronal identity specification in Caenorhabditis elegans. Single CUT mutants show barely any effects on pan-neuronal gene expression or global nervous system function, but such effects become apparent and progressively worsen upon removal of additional CUT family members, indicating a critical role of gene dosage. Overexpression of each individual CUT gene rescued the phenotype of compound mutants, corroborating that gene dosage, rather than the activity of specific members of the gene family, is critical for CUT gene family function. Genome-wide binding profiles, as well as mutation of CUT homeodomain binding sites by CRISPR/Cas9 genome engineering show that CUT genes directly control the expression of pan-neuronal features. Moreover, CUT genes act in conjunction with neuron-type-specific transcription factors to control pan-neuronal gene expression. Our study, therefore, provides a previously missing key insight into how neuronal gene expression programs are specified and reveals a highly buffered and robust mechanism that controls the most critical functional features of all neuronal cell types.

An expansion of the non-coding genome and its regulatory potential underlies vertebrate neuronal diversity

Proper assembly and function of the nervous system requires the generation of a uniquely diverse population of neurons expressing a cell-type-specific combination of effector genes that collectively define neuronal morphology, connectivity, and function. How countless partially overlapping but cell-type-specific patterns of gene expression are controlled at the genomic level remains poorly understood. Here we show that neuronal genes are associated with highly complex gene regulatory systems composed of independent cell-type- and cell-stage-specific regulatory elements that reside in expanded non-coding genomic domains. Mapping enhancer-promoter interactions revealed that motor neuron enhancers are broadly distributed across the large chromatin domains. This distributed regulatory architecture is not a unique property of motor neurons but is employed throughout the nervous system. The number of regulatory elements increased dramatically during the transition from invertebrates to vertebrates, suggesting that acquisition of new enhancers might be a fundamental process underlying the evolutionary increase in cellular complexity.

Transcriptional regulation of neuronal identity

Neuronal diversity is an intrinsic feature of the nervous system. Transcription factors (TFs) are key regulators in the establishment of different neuronal identities; how are the actions of different TFs coordinated to orchestrate this diversity? Are there common features shared among the different neuron types of an organism or even among different animal groups? In this review, we provide a brief overview on common traits emerging on the transcriptional regulation of neuron type diversification with a special focus on the comparison between mouse and Caenorhabditis elegans model systems. In the first part, we describe general concepts on neuronal identity and transcriptional regulation of gene expression. In the second part of the review, TFs are classified in different categories according to their key roles at specific steps along the protracted process of neuronal specification and differentiation. The same TF categories can be identified both in mammals and nematodes. Importantly, TFs are very pleiotropic: Depending on the neuron type or the time in development, the same TF can fulfil functions belonging to different categories. Finally, we describe the key role of transcriptional repression at all steps controlling neuronal diversity and propose that acquisition of neuronal identities could be considered a metastable process.

A C. elegans model of C9orf72-associated ALS/FTD uncovers a conserved role for eIF2D in RAN translation

A hexanucleotide repeat expansion GGGGCC in the non-coding region of C9orf72 is the most common cause of inherited amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Toxic dipeptide repeats (DPRs) are synthesized from GGGGCC via repeat-associated non-AUG (RAN) translation. Here, we develop C. elegans models that express, either ubiquitously or exclusively in neurons, 75 GGGGCC repeats flanked by intronic C9orf72 sequence. The worms generate DPRs (poly-glycine-alanine [poly-GA], poly-glycine-proline [poly-GP]) and poly-glycine-arginine [poly-GR]), display neurodegeneration, and exhibit locomotor and lifespan defects. Mutation of a non-canonical translation-initiating codon (CUG) upstream of the repeats selectively reduces poly-GA steady-state levels and ameliorates disease, suggesting poly-GA is pathogenic. Importantly, loss-of-function mutations in the eukaryotic translation initiation factor 2D (eif-2D/eIF2D) reduce poly-GA and poly-GP levels, and increase lifespan in both C. elegans models. Our in vitro studies in mammalian cells yield similar results. Here, we show a conserved role for eif-2D/eIF2D in DPR 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 nervous system-specific subnuclear organelle in Caenorhabditis elegans

We describe here phase-separated subnuclear organelles in the nematode Caenorhabditis elegans, which we term NUN (NUclear Nervous system-specific) bodies. Unlike other previously described subnuclear organelles, NUN bodies are highly cell type specific. In fully mature animals, 4-10 NUN bodies are observed exclusively in the nucleus of neuronal, glial and neuron-like cells, but not in other somatic cell types. Based on co-localization and genetic loss of function studies, NUN bodies are not related to other previously described subnuclear organelles, such as nucleoli, splicing speckles, paraspeckles, Polycomb bodies, promyelocytic leukemia bodies, gems, stress-induced nuclear bodies, or clastosomes. NUN bodies form immediately after cell cycle exit, before other signs of overt neuronal differentiation and are unaffected by the genetic elimination of transcription factors that control many other aspects of neuronal identity. In one unusual neuron class, the canal-associated neurons, NUN bodies remodel during larval development, and this remodeling depends on the Prd-type homeobox gene ceh-10. In conclusion, we have characterized here a novel subnuclear organelle whose cell type specificity poses the intriguing question of what biochemical process in the nucleus makes all nervous system-associated cells different from cells outside the nervous system.

Tissue-specific targeting of DNA nanodevices in a multicellular living organism

Nucleic acid nanodevices present great potential as agents for logic-based therapeutic intervention as well as in basic biology. Often, however, the disease targets that need corrective action are localized in specific organs, and thus realizing the full potential of DNA nanodevices also requires ways to target them to specific cell types in vivo. Here, we show that by exploiting either endogenous or synthetic receptor-ligand interactions and leveraging the biological barriers presented by the organism, we can target extraneously introduced DNA nanodevices to specific cell types in Caenorhabditis elegans, with subcellular precision. The amenability of DNA nanostructures to tissue-specific targeting in vivo significantly expands their utility in biomedical applications and discovery biology.

A perspective on C. elegans neurodevelopment: from early visionaries to a booming neuroscience research

The formation of the nervous system and its striking complexity is a remarkable feat of development. C. elegans served as a unique model to dissect the molecular events in neurodevelopment, from its early visionaries to the current booming neuroscience community. Soon after being introduced as a model, C. elegans was mapped at the level of genes, cells, and synapses, providing the first metazoan with a complete cell lineage, sequenced genome, and connectome. Here, I summarize mechanisms underlying C. elegans neurodevelopment, from the generation and diversification of neural components to their navigation and connectivity. I point out recent noteworthy findings in the fields of glia biology, sex dimorphism and plasticity in neurodevelopment, highlighting how current research connects back to the pioneering studies by Brenner, Sulston and colleagues. Multifaceted investigations in model organisms, connecting genes to cell function and behavior, expand our mechanistic understanding of neurodevelopment while allowing us to formulate emerging questions for future discoveries.

Molecular topography of an entire nervous system

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

Piecemeal regulation of convergent neuronal lineages by bHLH transcription factors in Caenorhabditis elegans

Cells of the same type can be generated by distinct cellular lineages that originate in different parts of the developing embryo ('lineage convergence'). Several Caenorhabditis elegans neuron classes composed of left/right or radially symmetric class members display such lineage convergence. We show here that the C. elegans Atonal homolog lin-32 is differentially expressed in neuronal lineages that give rise to left/right or radially symmetric class members. Loss of lin-32 results in the selective loss of the expression of pan-neuronal markers and terminal selector-type transcription factors that confer neuron class-specific features. Another basic helix-loop-helix (bHLH) gene, the Achaete-Scute homolog hlh-14, is expressed in a mirror image pattern relative to lin-32 and is required to induce neuronal identity and terminal selector expression on the contralateral side of the animal. These findings demonstrate that distinct lineage histories converge via different bHLH factors at the level of induction of terminal selector identity determinants, which thus serve as integrators of distinct lineage histories. We also describe neuron-to-neuron identity transformations in lin-32 mutants, which we propose to also be the result of misregulation of terminal selector gene expression.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

Non-autonomous regulation of germline stem cell proliferation by somatic MPK-1/MAPK activity in C. elegans

Extracellular-signal-regulated kinase (ERK)/mitogen-activated protein kinase (MAPK) is a major positive regulator of cell proliferation, which is often upregulated in cancer. However, few studies have addressed ERK/MAPK regulation of proliferation within a complete organism. The Caenorhabditis elegans ERK/MAPK ortholog MPK-1 is best known for its control of somatic organogenesis and germline differentiation, but it also stimulates germline stem cell proliferation. Here, we show that the germline-specific MPK-1B isoform promotes germline differentiation but has no apparent role in germline stem cell proliferation. By contrast, the soma-specific MPK-1A isoform promotes germline stem cell proliferation non-autonomously. Indeed, MPK-1A functions in the intestine or somatic gonad to promote germline proliferation independent of its other known roles. We propose that a non-autonomous role of ERK/MAPK in stem cell proliferation may be conserved across species and various tissue types, with major clinical implications for cancer and other diseases.

The prevailing paradigm is that ERK/MAPK functions autonomously to promote cell proliferation upon mitogen stimulation. Robinson-Thiewes et al. now demonstrate that C. elegans ERK/MAPK acts within somatic tissues to non-autonomously promote the proliferation of germline stem cells. Germline ERK/MAPK is thus dispensable for germline stem cell proliferation.

DAF-16/FoxO and DAF-12/VDR control cellular plasticity both cell-autonomously and via interorgan signaling

Many cell types display the remarkable ability to alter their cellular phenotype in response to specific external or internal signals. Such phenotypic plasticity is apparent in the nematode Caenorhabditis elegans when adverse environmental conditions trigger entry into the dauer diapause stage. This entry is accompanied by structural, molecular, and functional remodeling of a number of distinct tissue types of the animal, including its nervous system. The transcription factor (TF) effectors of 3 different hormonal signaling systems, the insulin-responsive DAF-16/FoxO TF, the TGFβ-responsive DAF-3/SMAD TF, and the steroid nuclear hormone receptor, DAF-12/VDR, a homolog of the vitamin D receptor (VDR), were previously shown to be required for entering the dauer arrest stage, but their cellular and temporal focus of action for the underlying cellular remodeling processes remained incompletely understood. Through the generation of conditional alleles that allowed us to spatially and temporally control gene activity, we show here that all 3 TFs are not only required to initiate tissue remodeling upon entry into the dauer stage, as shown before, but are also continuously required to maintain the remodeled state. We show that DAF-3/SMAD is required in sensory neurons to promote and then maintain animal-wide tissue remodeling events. In contrast, DAF-16/FoxO or DAF-12/VDR act cell-autonomously to control anatomical, molecular, and behavioral remodeling events in specific cell types. Intriguingly, we also uncover non-cell autonomous function of DAF-16/FoxO and DAF-12/VDR in nervous system remodeling, indicating the presence of several insulin-dependent interorgan signaling axes. Our findings provide novel perspectives into how hormonal systems control tissue remodeling.

3DeeCellTracker, a deep learning-based pipeline for segmenting and tracking cells in 3D time lapse images

Despite recent improvements in microscope technologies, segmenting and tracking cells in three-dimensional time-lapse images (3D + T images) to extract their dynamic positions and activities remains a considerable bottleneck in the field. We developed a deep learning-based software pipeline, 3DeeCellTracker, by integrating multiple existing and new techniques including deep learning for tracking. With only one volume of training data, one initial correction, and a few parameter changes, 3DeeCellTracker successfully segmented and tracked ~100 cells in both semi-immobilized and 'straightened' freely moving worm's brain, in a naturally beating zebrafish heart, and ~1000 cells in a 3D cultured tumor spheroid. While these datasets were imaged with highly divergent optical systems, our method tracked 90-100% of the cells in most cases, which is comparable or superior to previous results. These results suggest that 3DeeCellTracker could pave the way for revealing dynamic cell activities in image datasets that have been difficult to analyze.

NvPOU4/Brain3 Functions as a Terminal Selector Gene in the Nervous System of the Cnidarian Nematostella vectensis

Terminal selectors are transcription factors that control the morphological, physiological, and molecular features that characterize distinct cell types. Here, we show that, in the sea anemone Nematostella vectensis, NvPOU4 is expressed in post-mitotic cells that give rise to a diverse set of neural cell types, including cnidocytes and NvElav1-expressing neurons. Morphological analyses of NvPOU4 mutants crossed to transgenic reporter lines show that the loss of NvPOU4 does not affect the initial specification of neural cells. Transcriptomes derived from the mutants and from different neural cell populations reveal that NvPOU4 is required for the execution of the terminal differentiation program of these neural cells. These findings suggest that POU4 genes have ancient functions as terminal selectors for morphologically and functionally disparate types of neurons and they provide experimental support for the relevance of terminal selectors for understanding the evolution of cell types.

NeuroPAL: A Multicolor Atlas for Whole-Brain Neuronal Identification in C. elegans

Comprehensively resolving neuronal identities in whole-brain images is a major challenge. We achieve this in C. elegans by engineering a multicolor transgene called NeuroPAL (a neuronal polychromatic atlas of landmarks). NeuroPAL worms share a stereotypical multicolor fluorescence map for the entire hermaphrodite nervous system that resolves all neuronal identities. Neurons labeled with NeuroPAL do not exhibit fluorescence in the green, cyan, or yellow emission channels, allowing the transgene to be used with numerous reporters of gene expression or neuronal dynamics. We showcase three applications that leverage NeuroPAL for nervous-system-wide neuronal identification. First, we determine the brainwide expression patterns of all metabotropic receptors for acetylcholine, GABA, and glutamate, completing a map of this communication network. Second, we uncover changes in cell fate caused by transcription factor mutations. Third, we record brainwide activity in response to attractive and repulsive chemosensory cues, characterizing multimodal coding for these stimuli.

Expression pattern determines regulatory logic

Large amounts of effort have been invested in trying to understand how a single genome is able to specify the identity of hundreds of cell types. Inspired by some aspects of Caenorhabditis elegans biology, we implemented an in silico evolutionary strategy to produce gene regulatory networks (GRNs) that drive cell-specific gene expression patterns, mimicking the process of terminal cell differentiation. Dynamics of the gene regulatory networks are governed by a thermodynamic model of gene expression, which uses DNA sequences and transcription factor degenerate position weight matrixes as input. In a version of the model, we included chromatin accessibility. Experimentally, it has been determined that cell-specific and broadly expressed genes are regulated differently. In our in silico evolved GRNs, broadly expressed genes are regulated very redundantly and the architecture of their cis-regulatory modules is different, in accordance to what has been found in C. elegans and also in other systems. Finally, we found differences in topological positions in GRNs between these two classes of genes, which help to explain why broadly expressed genes are so resilient to mutations. Overall, our results offer an explanatory hypothesis on why broadly expressed genes are regulated so redundantly compared to cell-specific genes, which can be extrapolated to phenomena such as ChIP-seq HOT regions.

The Conserved ASCL1/MASH-1 Ortholog HLH-3 Specifies Sex-Specific Ventral Cord Motor Neuron Fate in Caenorhabditis elegans

Neural specification is regulated by one or many transcription factors that control expression of effector genes that mediate function and determine neuronal type. Here we identify a novel role for one conserved proneural factor, the bHLH protein HLH-3, implicated in the specification of sex-specific ventral cord motor neurons in C. elegans. Proneural genes act in early stages of neurogenesis in early progenitors, but here, we demonstrate a later role for hlh-3. First, we document that differentiation of the ventral cord type C motor neuron class (VC) within their neuron class, is dynamic in time and space. Expression of VC class-specific and subclass-specific identity genes is distinct through development and is dependent on the VC position along the A-P axis and their proximity to the vulva. Our characterization of the expression of VC class and VC subclass-specific differentiation markers in the absence of hlh-3 function reveals that VC fate specification, differentiation, and morphology requires hlh-3 function. Finally, we conclude that hlh-3 cell-autonomously specifies VC cell fate.

Direct glia-to-neuron transdifferentiation gives rise to a pair of male-specific neurons that ensure nimble male mating

Sexually dimorphic behaviours require underlying differences in the nervous system between males and females. The extent to which nervous systems are sexually dimorphic and the cellular and molecular mechanisms that regulate these differences are only beginning to be understood. We reveal here a novel mechanism by which male-specific neurons are generated in Caenorhabditis elegans through the direct transdifferentiation of sex-shared glial cells. This glia-to-neuron cell fate switch occurs during male sexual maturation under the cell-autonomous control of the sex-determination pathway. We show that the neurons generated are cholinergic, peptidergic, and ciliated putative proprioceptors which integrate into male-specific circuits for copulation. These neurons ensure coordinated backward movement along the mate's body during mating. One step of the mating sequence regulated by these neurons is an alternative readjustment movement performed when intromission becomes difficult to achieve. Our findings reveal programmed transdifferentiation as a developmental mechanism underlying flexibility in innate behaviour.

cnd-1/NeuroD1 Functions with the Homeobox Gene ceh-5/Vax2 and Hox Gene ceh-13/labial To Specify Aspects of RME and DD Neuron Fate in Caenorhabditis elegans

Identifying the mechanisms behind neuronal fate specification are key to understanding normal neural development in addition to neurodevelopmental disorders such as autism and schizophrenia. In vivo cell fate specification is difficult to study in vertebrates. However, the nematode Caenorhabditis elegans, with its invariant cell lineage and simple nervous system of 302 neurons, is an ideal organism to explore the earliest stages of neural development. We used a comparative transcriptome approach to examine the role of cnd-1/NeuroD1 in C. elegans nervous system development and function. This basic helix-loop-helix transcription factor is deeply conserved across phyla and plays a crucial role in cell fate specification in both the vertebrate nervous system and pancreas. We find that cnd-1 controls expression of ceh-5, a Vax2-like homeobox class transcription factor, in the RME head motorneurons and PVQ tail interneurons. We also show that cnd-1 functions redundantly with the Hox gene ceh-13/labial in defining the fate of DD1 and DD2 embryonic ventral nerve cord motorneurons. These data highlight the utility of comparative transcriptomes for identifying transcription factor targets and understanding gene regulatory networks.

APP-Induced Patterned Neurodegeneration Is Exacerbated by APOE4 in Caenorhabditis elegans

Genetic and epidemiological studies have found that variations in the amyloid precursor protein (APP) and the apoliopoprotein E (APOE) genes represent major modifiers of the progressive neurodegeneration in Alzheimer's disease (AD). An extra copy of or gain-of-function mutations in APP correlate with early onset AD. Compared to the other variants (APOE2 and APOE3), the ε4 allele of APOE (APOE4) hastens and exacerbates early and late onset forms of AD. Convenient in vivo models to study how APP and APOE4 interact at the cellular and molecular level to influence neurodegeneration are lacking. Here, we show that the nematode C. elegans can model important aspects of AD including age-related, patterned neurodegeneration that is exacerbated by APOE4 Specifically, we found that APOE4, but not APOE3, acts with APP to hasten and expand the pattern of cholinergic neurodegeneration caused by APP Molecular mechanisms underlying how APP and APOE4 synergize to kill some neurons while leaving others unaffected may be uncovered using this convenient worm model of neurodegeneration.

Combining single-cell RNA-sequencing with a molecular atlas unveils new markers for Caenorhabditis elegans neuron classes

Single-cell RNA-sequencing (scRNA-seq) of the Caenorhabditis elegans nervous system offers the unique opportunity to obtain a partial expression profile for each neuron within a known connectome. Building on recent scRNA-seq data and on a molecular atlas describing the expression pattern of ∼800 genes at the single cell resolution, we designed an iterative clustering analysis aiming to match each cell-cluster to the ∼100 anatomically defined neuron classes of C. elegans. This heuristic approach successfully assigned 97 of the 118 neuron classes to a cluster. Sixty two clusters were assigned to a single neuron class and 15 clusters grouped neuron classes sharing close molecular signatures. Pseudotime analysis revealed a maturation process occurring in some neurons (e.g. PDA) during the L2 stage. Based on the molecular profiles of all identified neurons, we predicted cell fate regulators and experimentally validated unc-86 for the normal differentiation of RMG neurons. Furthermore, we observed that different classes of genes functionally diversify sensory neurons, interneurons and motorneurons. Finally, we designed 15 new neuron class-specific promoters validated in vivo. Amongst them, 10 represent the only specific promoter reported to this day, expanding the list of neurons amenable to genetic manipulations.