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.

Caenorhabditis elegans ETR-1/CELF has broad effects on the muscle cell transcriptome, including genes that regulate translation and neuroblast migration.

Migration of neuroblasts and neurons from their birthplace is central to the formation of neural circuits and networks. ETR-1 is the Caenorhabditis elegans homolog of the CELF1 (CUGBP, ELAV-like family 1) RNA-processing factor involved in neuromuscular disorders. etr-1 regulates body wall muscle differentiation. Our previous work showed that etr-1 in muscle has a non-autonomous role in neuronal migration, suggesting that ETR-1 is involved in the production of a signal emanating from body wall muscle that controls neuroblast migration and that interacts with Wnt signaling. etr-1 is extensively alternatively-spliced, and we identified the viable etr-1(lq61) mutant, caused by a stop codon in alternatively-spliced exon 8 and only affecting etr-1 isoforms containing exon 8. We took advantage of viable etr-1(lq61) to identify potential RNA targets of ETR-1 in body wall muscle using a combination of fluorescence activated cell sorting (FACS) of body wall muscles from wild-type and etr-1(lq61) and subsequent RNA-seq. This analysis revealed genes whose splicing and transcript levels were controlled by ETR-1 exon 8 isoforms, and represented a broad spectrum of genes involved in muscle differentiation, myofilament lattice structure, and physiology. Genes with transcripts underrepresented in etr-1(lq61) included those involved in ribosome function and translation, similar to potential CELF1 targets identified in chick cardiomyocytes. This suggests that at least some targets of ETR-1 might be conserved in vertebrates, and that ETR-1 might generally stimulate translation in muscles. As proof-of-principle, a functional analysis of a subset of ETR-1 targets revealed genes involved in AQR and PQR neuronal migration. One such gene, lev-11/tropomyosin, requires ETR-1 for alternative splicing, and another, unc-52/perlecan, requires ETR-1 for the production of long isoforms containing 3' exons. In sum, these studies identified gene targets of ETR-1/CELF1 in muscles, which included genes involved in muscle development and physiology, and genes with novel roles in neuronal migration.

A spatial and temporal map of C. elegans gene expression.

The C. elegans genome has been completely sequenced, and the developmental anatomy of this model organism is described at single-cell resolution. Here we utilize strategies that exploit this precisely defined architecture to link gene expression to cell type. We obtained RNAs from specific cells and from each developmental stage using tissue-specific promoters to mark cells for isolation by FACS or for mRNA extraction by the mRNA-tagging method. We then generated gene expression profiles of more than 30 different cells and developmental stages using tiling arrays. Machine-learning-based analysis detected transcripts corresponding to established gene models and revealed novel transcriptionally active regions (TARs) in noncoding domains that comprise at least 10% of the total C. elegans genome. Our results show that about 75% of transcripts with detectable expression are differentially expressed among developmental stages and across cell types. Examination of known tissue- and cell-specific transcripts validates these data sets and suggests that newly identified TARs may exercise cell-specific functions. Additionally, we used self-organizing maps to define groups of coregulated transcripts and applied regulatory element analysis to identify known transcription factor- and miRNA-binding sites, as well as novel motifs that likely function to control subsets of these genes. By using cell-specific, whole-genome profiling strategies, we have detected a large number of novel transcripts and produced high-resolution gene expression maps that provide a basis for establishing the roles of individual genes in cellular differentiation.

Receptor-type guanylate cyclase is required for carbon dioxide sensation by Caenorhabditis elegans.

CO(2) is both a critical regulator of animal physiology and an important sensory cue for many animals for host detection, food location, and mate finding. The free-living soil nematode Caenorhabditis elegans shows CO(2) avoidance behavior, which requires a pair of ciliated sensory neurons, the BAG neurons. Using in vivo calcium imaging, we show that CO(2) specifically activates the BAG neurons and that the CO(2)-sensing function of BAG neurons requires TAX-2/TAX-4 cyclic nucleotide-gated ion channels and the receptor-type guanylate cyclase GCY-9. Our results delineate a molecular pathway for CO(2) sensing and suggest that activation of a receptor-type guanylate cyclase is an evolutionarily conserved mechanism by which animals detect environmental CO(2).

Patient-specific variants of NFU1/NFU-1 disrupt cholinergic signaling in a model of multiple mitochondrial dysfunctions syndrome 1.

Neuromuscular dysfunction is a common feature of mitochondrial diseases and frequently presents as ataxia, spasticity and/or dystonia, all of which can severely impact individuals with mitochondrial diseases. Dystonia is one of the most common symptoms of multiple mitochondrial dysfunctions syndrome 1 (MMDS1), a disease associated with mutations in the causative gene (NFU1) that impair iron-sulfur cluster biogenesis. We have generated Caenorhabditis elegans strains that recreated patient-specific point variants in the C. elegans ortholog (nfu-1) that result in allele-specific dysfunction. Each of these mutants, Gly147Arg and Gly166Cys, have altered acetylcholine signaling at neuromuscular junctions, but opposite effects on activity and motility. We found that the Gly147Arg variant was hypersensitive to acetylcholine and that knockdown of acetylcholine release rescued nearly all neuromuscular phenotypes of this variant. In contrast, we found that the Gly166Cys variant caused predominantly postsynaptic acetylcholine hypersensitivity due to an unclear mechanism. These results are important for understanding the neuromuscular conditions of MMDS1 patients and potential avenues for therapeutic intervention.

A head-to-head comparison of ribodepletion and polyA selection approaches for Caenorhabditis elegans low input RNA-sequencing libraries.

A recent and powerful technique is to obtain transcriptomes from rare cell populations, such as single neurons in Caenorhabditis elegans, by enriching dissociated cells using fluorescent sorting. However, these cell samples often have low yields of RNA that present challenges in library preparation. This can lead to PCR duplicates, noisy gene expression for lowly expressed genes, and other issues that limit endpoint analysis. Furthermore, some common resources, such as sequence-specific kits for removing ribosomal RNA, are not optimized for nonmammalian samples. To advance library construction for such challenging samples, we compared two approaches for building RNAseq libraries from less than 10 nanograms of C. elegans RNA: SMARTSeq V4 (Takara), a widely used kit for selecting poly-adenylated transcripts; and SoLo Ovation (Tecan Genomics), a newly developed ribodepletion-based approach. For ribodepletion, we used a custom kit of 200 probes designed to match C. elegans rRNA gene sequences. We found that SoLo Ovation, in combination with our custom C. elegans probe set for rRNA depletion, detects an expanded set of noncoding RNAs, shows reduced noise in lowly expressed genes, and more accurately counts expression of long genes. The approach described here should be broadly useful for similar efforts to analyze transcriptomics when RNA is limiting.

cAMP controls a trafficking mechanism that maintains the neuron specificity and subcellular placement of electrical synapses.

Electrical synapses are established between specific neurons and within distinct subcellular compartments, but the mechanisms that direct gap junction assembly in the nervous system are largely unknown. Here, we show that a developmental program tunes cAMP signaling to direct the neuron-specific assembly and placement of electrical synapses in the C. elegans motor circuit. We use live-cell imaging to visualize electrical synapses in vivo and an optogenetic assay to confirm that they are functional. In ventral A class (VA) motor neurons, the UNC-4 transcription factor blocks expression of cAMP antagonists that promote gap junction miswiring. In unc-4 mutants, VA electrical synapses are established with an alternative synaptic partner and are repositioned from the VA axon to soma. cAMP counters these effects by driving gap junction trafficking into the VA axon for electrical synapse assembly. Thus, our experiments establish that cAMP regulates gap junction trafficking for the biogenesis of functional electrical synapses.

Epidermal Growth Factor Signaling Promotes Sleep through a Combined Series and Parallel Neural Circuit.

Sleep requires sleep-active neurons that depolarize to inhibit wake circuits. Sleep-active neurons are under the control of homeostatic mechanisms that determine sleep need. However, little is known about the molecular and circuit mechanisms that translate sleep need into the depolarization of sleep-active neurons. During many stages and conditions in C. elegans, sleep requires a sleep-active neuron called RIS. Here, we defined the transcriptome of RIS and discovered that genes of the epidermal growth factor receptor (EGFR) signaling pathway are expressed in RIS. Because of cellular stress, EGFR directly activates RIS. Activation of EGFR signaling in the ALA neuron has previously been suggested to promote sleep independently of RIS. Unexpectedly, we found that ALA activation promotes RIS depolarization. Our results suggest that ALA is a drowsiness neuron with two separable functions: (1) it inhibits specific behaviors, such as feeding, independently of RIS, (2) and it activates RIS. Whereas ALA plays a strong role in surviving cellular stress, surprisingly, RIS does not. In summary, EGFR signaling can depolarize RIS by an indirect mechanism through activation of the ALA neuron that acts upstream of the sleep-active RIS neuron and through a direct mechanism using EGFR signaling in RIS. ALA-dependent drowsiness, rather than RIS-dependent sleep bouts, appears to be important for increasing survival after cellular stress, suggesting that different types of behavioral inhibition play different roles in restoring health. VIDEO ABSTRACT.

Food Sensation Modulates Locomotion by Dopamine and Neuropeptide Signaling in a Distributed Neuronal Network.

Finding food and remaining at a food source are crucial survival strategies. We show how neural circuits and signaling molecules regulate these food-related behaviors in Caenorhabditis elegans. In the absence of food, AVK interneurons release FLP-1 neuropeptides that inhibit motorneurons to regulate body posture and velocity, thereby promoting dispersal. Conversely, AVK photoinhibition promoted dwelling behavior. We identified FLP-1 receptors required for these effects in distinct motoneurons. The DVA interneuron antagonizes signaling from AVK by releasing cholecystokinin-like neuropeptides that potentiate cholinergic neurons, in response to dopaminergic neurons that sense food. Dopamine also acts directly on AVK via an inhibitory dopamine receptor. Both AVK and DVA couple to head motoneurons by electrical and chemical synapses to orchestrate either dispersal or dwelling behavior, thus integrating environmental and proprioceptive signals. Dopaminergic regulation of food-related behavior, via similar neuropeptides, may be conserved in mammals.

Inhibiting poly(ADP-ribosylation) improves axon regeneration.

The ability of a neuron to regenerate its axon after injury depends in part on its intrinsic regenerative potential. Here, we identify novel intrinsic regulators of axon regeneration: poly(ADP-ribose) glycohodrolases (PARGs) and poly(ADP-ribose) polymerases (PARPs). PARGs, which remove poly(ADP-ribose) from proteins, act in injured C. elegans GABA motor neurons to enhance axon regeneration. PARG expression is regulated by DLK signaling, and PARGs mediate DLK function in enhancing axon regeneration. Conversely, PARPs, which add poly(ADP-ribose) to proteins, inhibit axon regeneration of both C. elegans GABA neurons and mammalian cortical neurons. Furthermore, chemical PARP inhibitors improve axon regeneration when administered after injury. Our results indicate that regulation of poly(ADP-ribose) levels is a critical function of the DLK regeneration pathway, that poly-(ADP ribosylation) inhibits axon regeneration across species, and that chemical inhibition of PARPs can elicit axon regeneration.

Neuroendocrine modulation sustains the C. elegans forward motor state.

Neuromodulators shape neural circuit dynamics. Combining electron microscopy, genetics, transcriptome profiling, calcium imaging, and optogenetics, we discovered a peptidergic neuron that modulates C. elegans motor circuit dynamics. The Six/SO-family homeobox transcription factor UNC-39 governs lineage-specific neurogenesis to give rise to a neuron RID. RID bears the anatomic hallmarks of a specialized endocrine neuron: it harbors near-exclusive dense core vesicles that cluster periodically along the axon, and expresses multiple neuropeptides, including the FMRF-amide-related FLP-14. RID activity increases during forward movement. Ablating RID reduces the sustainability of forward movement, a phenotype partially recapitulated by removing FLP-14. Optogenetic depolarization of RID prolongs forward movement, an effect reduced in the absence of FLP-14. Together, these results establish the role of a neuroendocrine cell RID in sustaining a specific behavioral state in C. elegans.

Transcriptional Control of Synaptic Remodeling through Regulated Expression of an Immunoglobulin Superfamily Protein.

Neural circuits are actively remodeled during brain development, but the molecular mechanisms that trigger circuit refinement are poorly understood. Here, we describe a transcriptional program in C. elegans that regulates expression of an Ig domain protein, OIG-1, to control the timing of synaptic remodeling. DD GABAergic neurons reverse polarity during larval development by exchanging the locations of pre- and postsynaptic components. In newly born larvae, DDs receive cholinergic inputs in the dorsal nerve cord. These inputs are switched to the ventral side by the end of the first larval (L1) stage. VD class GABAergic neurons are generated in the late L1 and are postsynaptic to cholinergic neurons in the dorsal nerve cord but do not remodel. We investigated remodeling of the postsynaptic apparatus in DD and VD neurons using targeted expression of the acetylcholine receptor (AChR) subunit, ACR-12::GFP. We determined that OIG-1 antagonizes the relocation of ACR-12 from the dorsal side in L1 DD neurons. During the L1/L2 transition, OIG-1 is downregulated in DD neurons by the transcription factor IRX-1/Iroquois, allowing the repositioning of synaptic inputs to the ventral side. In VD class neurons, which normally do not remodel, the transcription factor UNC-55/COUP-TF turns off IRX-1, thus maintaining high levels of OIG-1 to block the removal of dorsally located ACR-12 receptors. OIG-1 is secreted from GABA neurons, but its anti-plasticity function is cell autonomous and may not require secretion. Our study provides a novel mechanism by which synaptic remodeling is set in motion through regulated expression of an Ig domain protein.

Isolation of specific neurons from C. elegans larvae for gene expression profiling.

BACKGROUND: The simple and well-described structure of the C. elegans nervous system offers an unprecedented opportunity to identify the genetic programs that define the connectivity and function of individual neurons and their circuits. A correspondingly precise gene expression map of C. elegans neurons would facilitate the application of genetic methods toward this goal. Here we describe a powerful new approach, SeqCeL (RNA-Seq of C. elegans cells) for producing gene expression profiles of specific larval C. elegans neurons. METHODS AND RESULTS: We have exploited available GFP reporter lines for FACS isolation of specific larval C. elegans neurons for RNA-Seq analysis. Our analysis showed that diverse classes of neurons are accessible to this approach. To demonstrate the applicability of this strategy to rare neuron types, we generated RNA-Seq profiles of the NSM serotonergic neurons that occur as a single bilateral pair of cells in the C. elegans pharynx. These data detected >1,000 NSM enriched transcripts, including the majority of previously known NSM-expressed genes. SIGNIFICANCE: This work offers a simple and robust protocol for expression profiling studies of post-embryonic C. elegans neurons and thus provides an important new method for identifying candidate genes for key roles in neuron-specific development and function.