Cell-specific microarray profiling experiments reveal a comprehensive picture of gene expression in the C. elegans nervous system

MAB-5/Hox regulates the Q neuroblast transcriptome, including cwn-1/Wnt, to mediate posterior migration in Caenorhabditis elegans

Neurogenesis involves the precisely coordinated action of genetic programs controlling large-scale neuronal fate specification down to terminal events of neuronal differentiation. The Q neuroblasts in Caenorhabditis elegans, QL on the left and QR on the right, divide, differentiate, and migrate in a similar pattern to produce three neurons each. However, QL on the left migrates posteriorly, and QR on the right migrates anteriorly. The MAB-5/Hox transcription factor is necessary and sufficient for posterior Q lineage migration and is normally expressed only in the QL lineage. To define genes controlled by MAB-5 in the Q cells, fluorescence-activated cell sorting was utilized to isolate populations of Q cells at a time in early L1 larvae when MAB-5 first becomes active. Sorted Q cells from wild-type, mab-5 loss-of-function (lof), and mab-5 gain-of-function (gof) mutants were subject to RNA-seq and differential expression analysis. Genes enriched in Q cells included those involved in cell division, DNA replication, and DNA repair, consist with the neuroblast stem cell identity of the Q cells at this stage. Genes affected by mab-5 included those involved in neurogenesis, neural development, and interaction with the extracellular matrix. cwn-1, which encodes a Wnt signaling molecule, showed a paired response to mab-5 in the Q cells: cwn-1 expression was reduced in mab-5(lof) and increased in mab-5(gof), suggesting that MAB-5 is required for cwn-1 expression in Q cells. MAB-5 is required to prevent anterior migration of the Q lineage while it transcriptionally reprograms the Q lineage for posterior migration. Functional genetic analysis revealed that CWN-1 is required downstream of MAB-5 to inhibit anterior migration of the QL lineage, likely in parallel to EGL-20/Wnt in a noncanonical Wnt pathway. In sum, work here describes a Q cell transcriptome, and a set of genes regulated by MAB-5 in the QL lineage. One of these genes, cwn-1, acts downstream of mab-5 in QL migration, indicating that this gene set includes other genes utilized by MAB-5 to facilitate posterior neuroblast migration.

Sink or swim: Does a worm paralysis phenotype hold clues to neurodegenerative disease?

Receiving a neurodegenerative disease (NDD) diagnosis, such as Alzheimer's disease, Parkinson's disease, Huntington's disease, or amyotrophic lateral sclerosis, is devastating, particularly given the limited options for treatment. Advances in genetic technologies have allowed for efficient modeling of NDDs in animals and brought hope for new disease-modifying medications. The complexity of the mammalian brain and the costs and time needed to identify and develop therapeutic leads limits progress. Modeling NDDs in invertebrates, such as the fruit fly Drosophila melanogaster and the nematode Caenorhabditis elegans, offers orders of magnitude increases in speed of genetic analysis and manipulation, and can be pursued at substantially reduced cost, providing an important, platform complement and inform research with mammalian NDD models. In this review, we describe how our efforts to exploit C. elegans for the study of neural signaling and health led to the discovery of a paralytic phenotype (swimming-induced paralysis) associated with altered dopamine signaling and, surprisingly, to the discovery of a novel gene and pathway whose dysfunction in glial cells triggers neurodegeneration. Research to date on swip-10 and its putative mammalian ortholog MBLAC1, suggests that a tandem analysis will offer insights into NDD mechanisms and insights into novel, disease-modifying therapeutics.

MAB-5/Hox regulates the Q neuroblast transcriptome, including cwn-1/Wnt, to mediate posterior migration in Caenorhabditis elegans

Neurogenesis involves the precisely-coordinated action of genetic programs controlling large-scale neuronal fate specification down to terminal events of neuronal differentiation. The Q neuroblasts in C. elegans, QL on the left and QR on the right, divide, differentiate, and migrate in a similar pattern to produce three neurons each. However, QL on the left migrates posteriorly, and QR on the right migrates anteriorly. The MAB-5/Hox transcription factor is necessary and sufficient for posterior Q lineage migration, and is normally expressed only in the QL lineage. To define genes controlled by MAB-5 in the Q cells, fluorescence-activated cell sorting was utilized to isolate populations of Q cells at a time in early L1 larvae when MAB-5 first becomes active. Sorted Q cells from wild-type, mab-5 loss-of-function (lof), and mab-5 gain-of-function (gof) mutants were subject to RNA-seq and differential expression analysis. Genes enriched in Q cells included those involved in cell division, DNA replication, and DNA repair, consist with the neuroblast stem cell identity of the Q cells at this stage. Genes affected by mab-5 included those involved in neurogenesis, neural development, and interaction with the extracellular matrix. cwn-1, which encodes a Wnt signaling molecule, showed a paired response to mab-5 in the Q cells: cwn-1 expression was reduced in mab-5(lof) and increased in mab-5(gof), suggesting that MAB-5 is required for cwn-1 expression in Q cells. MAB-5 is required to prevent anterior migration of the Q lineage while it transcriptionally reprograms the Q lineage for posterior migration. Functional genetic analysis revealed that CWN-1 is required downstream of MAB-5 to inhibit anterior migration of the QL lineage, likely in parallel to EGL-20/Wnt in a non-canonical Wnt pathway. In sum, work here describes a Q cell transcriptome, and a set of genes regulated by MAB-5 in the QL lineage. One of these genes, cwn-1, acts downstream of mab-5 in QL migration, indicating that this gene set includes other genes utilized by MAB-5 to facilitate posterior neuroblast migration.

3'sialyllactose and 6'sialyllactose enhance performance in endurance-type exercise through metabolic adaptation

Human milk oligosaccharides (HMOs) belong to a group of multifunctional glycans that are abundantly present in human breast milk. While health effects of neutral oligosaccharides have been investigated extensively, a lot remains unknown regarding health effects of acidic oligosaccharides, such as the two sialyllactoses (SLs), 3'sialyllactose (3'SL), and 6'sialyllactose (6'SL). We utilized Caenorhabditis elegans (C. elegans) to investigate the effects of SLs on exercise performance. Using swimming as an endurance-type exercise, we found that SLs decrease exhaustion, signifying an increase in endurance that is strongest for 6'SL. Through an unbiased metabolomics approach, we identified changes in energy metabolism that correlated with endurance performance. Further investigation suggested that these metabolic changes were related to adaptations of muscle mitochondria that facilitated a shift from beta oxidation to glycogenolysis during exercise. We found that the effect of SLs on endurance performance required AMPK- (aak-1/aak-2) and adenosine receptor (ador-1) signaling. We propose a model where SLs alter the metabolic status in the gut, causing a signal from the intestine to the nervous system toward muscle cells, where metabolic adaptation increases exercise performance. Together, our results underline the potential of SLs in exercise-associated health and contribute to our understanding of the molecular processes involved in nutritionally-induced health benefits.

The Caenorhabditis elegans Shugoshin regulates TAC-1 in cilia

The conserved Shugoshin (SGO) protein family is essential for mediating proper chromosome segregation from yeast to humans but has also been implicated in diverse roles outside of the nucleus. SGO's roles include inhibiting incorrect spindle attachment in the kinetochore, regulating the spindle assembly checkpoint (SAC), and ensuring centriole cohesion in the centrosome, all functions that involve different microtubule scaffolding structures in the cell. In Caenorhabditis elegans, a species with holocentric chromosomes, SGO-1 is not required for cohesin protection or spindle attachment but appears important for licensing meiotic recombination. Here we provide the first functional evidence that in C. elegans, Shugoshin functions in another extranuclear, microtubule-based structure, the primary cilium. We identify the centrosomal and microtubule-regulating transforming acidic coiled-coil protein, TACC/TAC-1, which also localizes to the basal body, as an SGO-1 binding protein. Genetic analyses indicate that TAC-1 activity must be maintained below a threshold at the ciliary base for correct cilia function, and that SGO-1 likely participates in constraining TAC-1 to the basal body by influencing the function of the transition zone 'ciliary gate'. This research expands our understanding of cellular functions of Shugoshin proteins and contributes to the growing examples of overlap between kinetochore, centrosome and cilia proteomes.

Nervous system-related gene regulatory networks and functional evolution of ETS proteins across species

The ETS domain transcription factor family is one of the major transcription factor superfamilies that play regulatory roles in development, cell growth, and cancer progression. Although different functions of ETS member proteins in the nervous system have been demonstrated in various studies, their role in neuronal cell differentiation and the evolutionary conservation of its target genes have not yet been extensively studied. In this study, we focused on the regulatory role of ETS transcription factors in neuronal differentiation and their functional evolution by comparative transcriptomics. In order to investigate the regulatory role of ETS transcription factors in neuronal differentiation across species, transcriptional profiles of ETS members and their target genes were investigated by comparing differentially expressed genes and gene regulatory networks, which were analyzed using human, gorilla, mouse, fruit fly and worm transcriptomics datasets. Bioinformatics approaches to examine the evolutionary conservation of ETS transcription factors during neuronal differentiation have shown that ETS member proteins regulate genes associated with neuronal differentiation, nervous system development, axon, and synaptic regulation in different organisms. This study is a comparative transcriptomic study of ETS transcription factors in terms of neuronal differentiation using a gene regulatory network inference algorithm. Overall, a comparison of gene regulation networks revealed that ETS members are indeed evolutionarily conserved in the regulation of neuronal differentiation. Nonetheless, ETS, PEA3, and ELF subfamilies were found to be relatively more active transcription factors in the transcriptional regulation of neuronal differentiation.

Laser Microdissection for Species-Agnostic Single-Tissue Applications

Single-cell methodologies have revolutionized the analysis of the transcriptomes of specific cell types. However, they often require species-specific genetic "toolkits," such as promoters driving tissue-specific expression of fluorescent proteins. Further, protocols that disrupt tissues to isolate individual cells remove cells from their native environment (e.g., signaling from neighbors) and may result in stress responses or other differences from native gene expression states. In the present protocol, laser microdissection (LMD) is optimized to isolate individual nematode tail tips for the study of gene expression during male tail tip morphogenesis. LMD allows the isolation of a portion of the animal without the need for cellular disruption or species-specific toolkits and is thus applicable to any species. Subsequently, single-cell RNA-seq library preparation protocols such as CEL-Seq2 can be applied to LMD-isolated single tissues and analyzed using standard pipelines, given that a well-annotated genome or transcriptome is available for the species. Such data can be used to establish how conserved or different the transcriptomes are that underlie the development of that tissue in different species. Limitations include the ability to cut out the tissue of interest and the sample size. A power analysis shows that as few as 70 tail tips per condition are required for 80% power. Tight synchronization of development is needed to obtain this number of animals at the same developmental stage. Thus, a method to synchronize animals at 1 h intervals is also described.

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.

Form and Function of the Vertebrate and Invertebrate Blood-Brain Barriers

The need to protect neural tissue from toxins or other substances is as old as neural tissue itself. Early recognition of this need has led to more than a century of investigation of the blood-brain barrier (BBB). Many aspects of this important neuroprotective barrier have now been well established, including its cellular architecture and barrier and transport functions. Unsurprisingly, most research has had a human orientation, using mammalian and other animal models to develop translational research findings. However, cell layers forming a barrier between vascular spaces and neural tissues are found broadly throughout the invertebrates as well as in all vertebrates. Unfortunately, previous scenarios for the evolution of the BBB typically adopt a classic, now discredited 'scala naturae' approach, which inaccurately describes a putative evolutionary progression of the mammalian BBB from simple invertebrates to mammals. In fact, BBB-like structures have evolved independently numerous times, complicating simplistic views of the evolution of the BBB as a linear process. Here, we review BBBs in their various forms in both invertebrates and vertebrates, with an emphasis on the function, evolution, and conditional relevance of popular animal models such as the fruit fly and the zebrafish to mammalian BBB research.

Reverse complementary matches simultaneously promote both back-splicing and exon-skipping

Background

Circular RNAs (circRNAs) play diverse roles in different biological and physiological environments and are always expressed in a tissue-specific manner. Especially, circRNAs are enriched in the brain tissues of almost all investigated species, including humans, mice, Drosophila, etc. Although circRNAs were found in C. elegans, the neuron-specific circRNA data is not available yet. Exon-skipping is found to be correlated to circRNA formation, but the mechanisms that link them together are not clear.

Results

Here, through large-scale neuron isolation from the first larval (L1) stage of C. elegans followed by RNA sequencing with ribosomal RNA depletion, the neuronal circRNA data in C. elegans were obtained. Hundreds of novel circRNAs were annotated with high accuracy. circRNAs were highly expressed in the neurons of C. elegans and were positively correlated to the levels of their cognate linear mRNAs. Disruption of reverse complementary match (RCM) sequences in circRNA flanking introns effectively abolished circRNA formation. In the zip-2 gene, deletion of either upstream or downstream RCMs almost eliminated the production of both the circular and the skipped transcript. Interestingly, the 13-nt RCM in zip-2 is highly conserved across five nematode ortholog genes, which show conserved exon-skipping patterns. Finally, through in vivo one-by-one mutagenesis of all the splicing sites and branch points required for exon-skipping and back-splicing in the zip-2 gene, I showed that back-splicing still happened without exon-skipping, and vice versa.

Conclusions

Through protocol optimization, total RNA obtained from sorted neurons is increased to hundreds of nanograms. circRNAs highly expressed in the neurons of C. elegans are more likely to be derived from genes also highly expressed in the neurons. RCMs are abundant in circRNA flanking introns, and RCM-deletion is an efficient way to knockout circRNAs. More importantly, these RCMs are not only required for back-splicing but also promote the skipping of exon(s) to be circularized. Finally, RCMs in circRNA flanking introns can directly promote both exon-skipping and back-splicing, providing a new explanation for the correlation between them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07910-w.

The phosphatidylinositol transfer protein PITP-1 facilitates fast recovery of eating behavior after hypoxia in the nematode Caenorhabditis elegans

Among the fascinating adaptations to limiting oxygen conditions (hypoxia) is the suppression of food intake and weight loss. In humans, this phenomenon is called high-altitude anorexia and is observed in people suffering from acute mountain syndrome. The high-altitude anorexia appears to be conserved in evolution and has been seen in species across the animal kingdom. However, the mechanism underlying the recovery of eating behavior after hypoxia is still not known. Here, we show that the phosphatidylinositol transfer protein PITP-1 is essential for the fast recovery of eating behavior after hypoxia in the nematode Caenorhabditis elegans. Unlike the neuroglobin GLB-5 that accelerates the recovery of eating behavior through its function in the oxygen (O₂ )-sensing neurons, PITP-1 appears to act downstream, in neurons that express the mod-1 serotonin receptor. Indeed, pitp-1 mutants display wild-type-like O₂ -evoked-calcium responses in the URX O₂ -sensing neuron. Intriguingly, loss-of-function of protein kinase C 1 (PKC-1) rescues pitp-1 mutants' recovery after hypoxia. Increased diacylglycerol (DAG), which activates PKC-1, attenuates the recovery of wild-type worms. Together, these data suggest that PITP-1 enables rapid recovery of eating behavior after hypoxia by limiting DAG's availability, thereby limiting PKC activity in mod-1-expressing neurons.

Impairment of the neurotrophic signaling hub B-Raf contributes to motoneuron degeneration in spinal muscular atrophy

Spinal muscular atrophy (SMA) is a motoneuron disease caused by deletions of the Survival of Motoneuron 1 gene (SMN1) and low SMN protein levels. SMN restoration is the concept behind a number of recently approved drugs which result in impressive yet limited effects. Since SMN has already been enhanced in treated patients, complementary SMN-independent approaches are needed. Previously, a number of altered signaling pathways which regulate motoneuron degeneration have been identified as candidate targets. However, signaling pathways form networks, and their connectivity is still unknown in SMA. Here, we used presymptomatic SMA mice to elucidate the network of altered signaling in SMA. The SMA network is structured in two clusters with AKT and 14-3-3 ζ/δ in their centers. Both clusters are connected by B-Raf as a major signaling hub. The direct interaction of B-Raf with 14-3-3 ζ/δ is important for an efficient neurotrophic activation of the MEK/ERK pathway and crucial for motoneuron survival. Further analyses in SMA mice revealed that both proteins were down-regulated in motoneurons and the spinal cord with B-Raf being reduced at presymptomatic stages. Primary fibroblasts and iPSC-derived motoneurons from SMA patients both showed the same pattern of down-regulation. This mechanism is conserved across species since a Caenorhabditis elegans SMA model showed less expression of the B-Raf homolog lin-45 Accordingly, motoneuron survival was rescued by a cell autonomous lin-45 expression in a C. elegans SMA model resulting in improved motor functions. This rescue was effective even after the onset of motoneuron degeneration and mediated by the MEK/ERK pathway.

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

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

Tissue-Specific Transcription Footprinting Using RNA PoI DamID (RAPID) in Caenorhabditis elegans

Differential gene expression across cell types underlies development and cell physiology in multicellular organisms. Caenorhabditis elegans is a powerful, extensively used model to address these biological questions. A remaining bottleneck relates to the difficulty to obtain comprehensive tissue-specific gene transcription data, since available methods are still challenging to execute and/or require large worm populations. Here, we introduce the RNA Polymerase DamID (RAPID) approach, in which the Dam methyltransferase is fused to a ubiquitous RNA polymerase subunit to create transcriptional footprints via methyl marks on the DNA of transcribed genes. To validate the method, we determined the polymerase footprints in whole animals, in sorted embryonic blastomeres and in different tissues from intact young adults by driving tissue-specific Dam fusion expression. We obtained meaningful transcriptional footprints in line with RNA-sequencing (RNA-seq) studies in whole animals or specific tissues. To challenge the sensitivity of RAPID and demonstrate its utility to determine novel tissue-specific transcriptional profiles, we determined the transcriptional footprints of the pair of XXX neuroendocrine cells, representing 0.2% of the somatic cell content of the animals. We identified 3901 candidate genes with putatively active transcription in XXX cells, including the few previously known markers for these cells. Using transcriptional reporters for a subset of new hits, we confirmed that the majority of them were expressed in XXX cells and identified novel XXX-specific markers. Taken together, our work establishes RAPID as a valid method for the determination of RNA polymerase footprints in specific tissues of C. elegans without the need for cell sorting or RNA tagging.

Clozapine, nimodipine and endosulfan differentially suppress behavioral defects caused by gain-of-function mutations in a two-pore domain K+ channel (UNC-58)

Two-pore domain K⁺ channels (K2Ps) regulate the resting membrane potential in excitable cells and determine ease of depolarization. Gain-of-function (gf) mutations in one of these channels (unc-58) in C. elegans switch it to a Na⁺ conductance channel and cause tremors, paralysis and other defects. We hypothesized that it should be possible to identify drugs that corrected these defects in unc-58(gf) mutant animals by blocking or modulating the over-active channels. We examined dispersal of animals on food because the absence of effective forward locomotion is the most obvious defect. In addition, we quantified egg release over 24 h. Starting with a known inhibitor of mammalian K2Ps and directed structure-based screening, we evaluated numerous drugs in these assays. Loratadine, which inhibits human KCNK18, significantly improved movement as did methiothepin. We confirmed that endosulfan, a GABA-A receptor antagonist, corrected locomotion in the unc-58(gf) strains. Based on structural similarities to other hits, we found that clozapine, loxapine and amoxapine potently suppressed abnormal phenotypes. Curiously, nimodipine, a Ca⁺⁺-channel blocker, dramatically improved movement and egg laying in unc-58(e665), but not unc-58(n495) animals. Molecular modeling provided initial insights into a possible basis for this difference based on the location of the e665 and n495 mutations. This research may lead to identification of novel K2P modulators and potential leads for drug discovery.

Neuronal transcriptome analyses reveal novel neuropeptide modulators of excitation and inhibition imbalance in C. elegans

Neuropeptides are secreted molecules that have conserved roles modulating many processes, including mood, reproduction, and feeding. Dysregulation of neuropeptide signaling is also implicated in neurological disorders such as epilepsy. However, much is unknown about the mechanisms regulating specific neuropeptides to mediate behavior. Here, we report that the expression levels of dozens of neuropeptides are up-regulated in response to circuit activity imbalance in C. elegans. acr-2 encodes a homolog of human nicotinic receptors, and functions in the cholinergic motoneurons. A hyperactive mutation, acr-2(gf), causes an activity imbalance in the motor circuit. We performed cell-type specific transcriptomic analysis and identified genes differentially expressed in acr-2(gf), compared to wild type. The most over-represented class of genes are neuropeptides, with insulin-like-peptides (ILPs) the most affected. Moreover, up-regulation of neuropeptides occurs in motoneurons, as well as sensory neurons. In particular, the induced expression of the ILP ins-29 occurs in the BAG neurons, which were previously shown to function in gas-sensing. We also show that this up-regulation of ins-29 in acr-2(gf) animals is activity-dependent. Our genetic and molecular analyses support cooperative effects for ILPs and other neuropeptides in promoting motor circuit activity in the acr-2(gf) background. Together, this data reveals that a major transcriptional response to motor circuit dysregulation is in up-regulation of multiple neuropeptides, and suggests that BAG sensory neurons can respond to intrinsic activity states to feedback on the motor circuit.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Sentryn and SAD Kinase Link the Guided Transport and Capture of Dense Core Vesicles in Caenorhabditis elegans

Dense core vesicles (DCVs) can transmit signals by releasing neuropeptides from specialized synaptic regions called active zones. DCVs reach the active zone by motorized transport through a long axon. A reverse motor frequently interrupts progress by taking DCVs in the opposite direction. "Guided transport" refers to the mechanism by which outward movements ultimately dominate to bring DCVs to the synaptic region. After guided transport, DCVs alter their interactions with motors and enter a "captured" state. The mechanisms of guided transport and capture of DCVs are unknown. Here, we discovered two proteins that contribute to both processes in Caenorhabditis elegans SAD kinase and a novel conserved protein we named Sentryn are the first proteins found to promote DCV capture. By imaging DCVs moving in various regions of single identified neurons in living animals, we found that DCV guided transport and capture are linked through SAD kinase, Sentryn, and Liprin-α. These proteins act together to regulate DCV motorized transport in a region-specific manner. Between the cell body and the synaptic region, they promote forward transport. In the synaptic region, where all three proteins are highly enriched at active zones, they promote DCV pausing by inhibiting transport in both directions. These three proteins appear to be part of a special subset of active zone-enriched proteins because other active zone proteins do not share their unique functions.

Revisiting Neuronal Cell Type Classification in Caenorhabditis elegans

We revisit the classification of neuronal cell types in the nervous system of the nematode Caenorhabditis elegans. Based on anatomy and synaptic connectivity patterns, the 302 neurons of the nervous system of the hermaphrodite were categorized into 118 neuron classes more than 30 years ago. Analysis of all presently available neuronal gene expression patterns reveals a remarkable congruence of anatomical and molecular classification and further suggests subclassification schemes. Transcription factor expression profiles alone are sufficient to uniquely classify more than 90% of all neuron classes in the C. elegans nervous system. Neuron classification in C. elegans may be paradigmatic for neuron classification schemes in vertebrate nervous systems.

JAK/STAT and TGF-ß activation as potential adverse outcome pathway of TiO2NPs phototoxicity in Caenorhabditis elegans

Titanium dioxide nanoparticles (TiO₂NPs) are widely used nanoparticles, whose catalytic activity is mainly due to photoactivation. In this study, the toxicity of TiO₂NPs was investigated on the nematode Caenorhabditis elegans, with and without UV activation. Comparative analyses across the four treatments revealed that UV-activated TiO₂NPs led to significant reproductive toxicity through oxidative stress. To understand the underlying molecular mechanism, transcriptomics and metabolomics analyses were conducted, followed by whole-genome network-based pathway analyses. Differential expression analysis from microarray data revealed only 4 DEGs by exposure to TiO₂NPs alone, compared to 3,625 and 3,286 DEGs by UV alone and UV-activated TiO₂NPs, respectively. Pathway analyses suggested the possible involvement of the JAK/STAT and TGF-ß pathways in the phototoxicity of TiO₂NPs, which correlated with the observation of increased gene expression of those pathways. Comparative analysis of C. elegans response across UV activation and TiO₂NPs exposure was performed using loss-of-function mutants of genes in these pathways. Results indicated that the JAK/STAT pathway was specific to TiO₂NPs, whereas the TGF-ß pathway was specific to UV. Interestingly, crosstalk between these pathways was confirmed by further mutant analysis. We consider that these findings will contribute to understand the molecular mechanisms of toxicity of TiO₂NPs in the natural environment.

An Elongin-Cullin-SOCS Box Complex Regulates Stress-Induced Serotonergic Neuromodulation

Neuromodulatory cells transduce environmental information into long-lasting behavioral responses. However, the mechanisms governing how neuronal cells influence behavioral plasticity are difficult to characterize. Here, we adapted the translating ribosome affinity purification (TRAP) approach in C. elegans to profile ribosome-associated mRNAs from three major tissues and the neuromodulatory dopaminergic and serotonergic cells. We identified elc-2, an Elongin C ortholog, specifically expressed in stress-sensing amphid neuron dual ciliated sensory ending (ADF) serotonergic sensory neurons, and we found that it plays a role in mediating a long-lasting change in serotonin-dependent feeding behavior induced by heat stress. We demonstrate that ELC-2 and the von Hippel-Lindau protein VHL-1, components of an Elongin-Cullin-SOCS box (ECS) E3 ubiquitin ligase, modulate this behavior after experiencing stress. Also, heat stress induces a transient redistribution of ELC-2, becoming more nuclearly enriched. Together, our results demonstrate dynamic regulation of an E3 ligase and a role for an ECS complex in neuromodulation and control of lasting behavioral states.

SAP97 Binding Partner CRIPT Promotes Dendrite Growth In Vitro and In Vivo

The dendritic tree is a key determinant of neuronal information processing. In the motor system, the dendritic tree of spinal cord neurons undergoes dramatic remodeling in an activity-dependent manner during early postnatal life. This leads to the proper segmental spinal cord connectivity that subserves normal locomotor behavior. One molecular system driving the establishment of dendrite architecture of mammalian motor neurons relies on AMPA receptors (AMPA-Rs) assembled with the GluA1 subunit, and this occurs in an NMDA receptor (NMDA-R)-independent manner. The dendrite growth promoting activity of GluA1-containing AMPA-Rs depends on its intracellular binding partner, SAP97, and SAP97's PDZ3 domain. We show here that cysteine-rich interactor of PDZ3 (CRIPT) is a bona fide SAP97 PDZ3-domain binding partner, localizes to synapses with GluA1 and SAP97 along the dendritic tree, and is a determinant of the dendritic growth of mammalian spinal cord neurons. We further show that CRIPT has a well-conserved ortholog in the nematode, Caenorhabditis elegans, and animals lacking CRIPT display decreased dendrite branching of the well-studied PVD neuron in vivo. The lack of CRIPT leads to a selective defect in touch perception, and this is rescued by expression of wild-type (WT) human CRIPT (hCRIPT) in the nervous system. This work brings new light into the molecular machinery that drives dendritic growth during development and may prove relevant to the promotion of nervous system plasticity following insult.

Interplay of neuronal and non-neuronal genes regulates intestinal DAF-16-mediated immune response during Fusarium infection of Caenorhabditis elegans

Although precisely controlled innate immune response is governed by conserved cellular events in phylogenetically diverse hosts, the underlying molecular mechanisms by which this process is regulated against a multi-host pathogen remain unknown. Fusarium oxysporum is a model multi-host pathogen, known to be associated with neuronal stress in humans and vascular wilt in plants. The interaction between innate immune and neuronal pathways is the basis of many diverse biological responses. How these processes are coordinated in response to fungal disease is not well understood. Here, we show that F. oxysporum f. sp. ciceri causes neuronal stress and intestinal disintegration, ultimately leading to the death of Caenorhabditis elegans. To explore the regulatory framework of Fusarium-associated disease, we analysed the gene expression during infection, integrated temporal gene expression, and network analysis with genetic inactivation data in Caenorhabditis elegans. We identified 1024 genes showing significant changes in expression (corrected P-values <0.05) in response to Fusarium infection. Co-expression network analysis of our data identified prognostic genes related to disease progression. These genes were dynamically expressed in various neuronal and non-neuronal tissues exhibiting diverse biological functions, including cellular homeostasis, organ patterning, stress response, and lipid metabolism. The RNA-seq analysis further identified shared and unique signalling pathways regulated by DAF-16/FOXO and SIR-2.1 linking neuronal stress, which facilitates negative regulation of intestinal innate immunity. Genetic analysis revealed that GCY-5 in ASE functions upstream of DAF-16, whereas ASI-specific SRD-1 regulates behavioural immunity. Overall, our results indicate that a ubiquitous response occurs during Fusarium infection mediated by highly conserved regulatory components and pathways, which can be exploited further for the identification of disease-responsive genes conserved among animals and plants. Finally, this study provided a novel insight into cross-species immune signalling and may facilitate the discovery of cellular therapeutic targets for Fusarium-associated disease.

Dorsal and Ventral Hippocampus Differentiate in Functional Pathways and Differentially Associate with Neurological Disease-Related Genes during Postnatal Development

The dorsal and ventral regions of the hippocampus are important in cognitive and emotional processing, respectively. Various approaches have revealed the differential molecular and structural characteristics, and functional roles of the hippocampus. Recent RNA sequencing (RNA-seq) technology has enriched our understanding of the hippocampus by elucidating more detailed information on gene expression patterns. However, no RNA-seq–based study on gene profiles in the developing hippocampus has been reported. Using RNA-seq–based bioinformatic analysis in conjunction with quantitative real-time polymerase chain reaction analysis and a comparison of in situ hybridization data obtained from the Allen Brain Atlas, we provide a thorough analysis of differentially expressed genes in the dorsal and ventral hippocampus at specific developmental ages representing the postnatally maturing hippocampus. Genes associated with particular functional pathways and marker genes for particular neurological diseases were found to be distinctively segregated within either the dorsal or ventral hippocampus at specific or at all developmental ages examined. We also report novel molecular markers enriched in the dorsal or ventral hippocampus. Taken together, this study provides insights into the molecular mechanisms underlying physiological functions linked to the dorsal or ventral hippocampus. The information provided in the study also contributes to a better understanding of brain functions and serves as a resource for future studies on the pathophysiology of dorsal and ventral hippocampal functions.

Characterization of gene expression associated with the adaptation of the nematode C. elegans to hypoxia and reoxygenation stress reveals an unexpected function of the neuroglobin GLB-5 in innate immunity

Oxygen (O₂) is a double-edged sword to cells, for while it is vital for energy production in all aerobic animals and insufficient O₂ (hypoxia) can lead to cell death, the reoxygenation of hypoxic tissues may trigger the generation of reactive oxygen species (ROS) that can destroy any biological molecule. Indeed, both hypoxia and hypoxia-reoxygenation (H/R) stress are harmful, and may play a critical role in the pathophysiology of many human diseases, such as myocardial ischemia and stroke. Therefore, understanding how animals adapt to hypoxia and H/R stress is critical for developing better treatments for these diseases. Previous studies showed that the neuroglobin GLB-5(Haw) is essential for the fast recovery of the nematode Caenorhabditis elegans (C. elegans) from H/R stress. Here, we characterize the changes in neuronal gene expression during the adaptation of worms to hypoxia and recovery from H/R stress. Our analysis shows that innate immunity genes are differentially expressed during both adaptation to hypoxia and recovery from H/R stress. Moreover, we reveal that the prolyl hydroxylase EGL-9, a known regulator of both adaptation to hypoxia and the innate immune response, inhibits the fast recovery from H/R stress through its activity in the O₂-sensing neurons AQR, PQR, and URX. Finally, we show that GLB-5(Haw) acts in AQR, PQR, and URX to increase the tolerance of worms to Pseudomonas aeruginosa pathogenesis. Together, our studies suggest that innate immunity and recovery from H/R stress are regulated by overlapping signaling pathways.

Cell type-specific transcriptome profiling in C. elegans using the Translating Ribosome Affinity Purification technique

Organs and specific cell types execute specialized functions in multicellular organisms, in large part through customized gene expression signatures. Thus, profiling the transcriptomes of specific cell and tissue types remains an important tool for understanding how cells become specialized. Methodological approaches to detect gene expression differences have utilized samples from whole animals, dissected tissues, and more recently single cells. Despite these advances, there is still a challenge and a need in most laboratories to implement less invasive yet powerful cell-type specific transcriptome profiling methods. Here, we describe the use of the Translating Ribosome Affinity Purification (TRAP) method for C. elegans to detect cell type-specific gene expression patterns at the level of translating mRNAs. In TRAP, a ribosomal protein is fused to a tag (GFP) and is expressed under cell type-specific promoters to mark genetically defined cell types in vivo. Affinity purification of lysates of animals expressing the tag enriches for ribosome-associated mRNAs of the targeted tissue. The purified mRNAs are used for making cDNA libraries subjected to high-throughput sequencing to obtain genome-wide profiles of transcripts from the targeted cell type. The ease of exposing C. elegans to diverse stimuli, coupled with available cell type specific promoters, makes TRAP a useful approach to enable the discovery of molecular components in response to external or genetic perturbations.

Genome-Wide Temporal Expression Profiling in Caenorhabditis elegans Identifies a Core Gene Set Related to Long-Term Memory

The identification of genes related to encoding, storage, and retrieval of memories is a major interest in neuroscience. In the current study, we analyzed the temporal gene expression changes in a neuronal mRNA pool during an olfactory long-term associative memory (LTAM) in Caenorhabditis elegans hermaphrodites. Here, we identified a core set of 712 (538 upregulated and 174 downregulated) genes that follows three distinct temporal peaks demonstrating multiple gene regulation waves in LTAM. Compared with the previously published positive LTAM gene set (Lakhina et al., 2015), 50% of the identified upregulated genes here overlap with the previous dataset, possibly representing stimulus-independent memory-related genes. On the other hand, the remaining genes were not previously identified in positive associative memory and may specifically regulate aversive LTAM. Our results suggest a multistep gene activation process during the formation and retrieval of long-term memory and define general memory-implicated genes as well as conditioning-type-dependent gene sets.SIGNIFICANCE STATEMENT The identification of genes regulating different steps of memory is of major interest in neuroscience. Identification of common memory genes across different learning paradigms and the temporal activation of the genes are poorly studied. Here, we investigated the temporal aspects of Caenorhabditis elegans gene expression changes using aversive olfactory associative long-term memory (LTAM) and identified three major gene activation waves. Like in previous studies, aversive LTAM is also CREB dependent, and CREB activity is necessary immediately after training. Finally, we define a list of memory paradigm-independent core gene sets as well as conditioning-dependent genes.

Ca2+ permeability and Na+ conductance in cellular toxicity caused by hyperactive DEG/ENaC channels

Hyperactivated DEG/ENaC channels cause neuronal death mediated by intracellular Ca²⁺ overload. Mammalian ASIC1a channels and MEC-4(d) neurotoxic channels in Caenorhabditis elegans both conduct Na⁺ and Ca²⁺, raising the possibility that direct Ca²⁺ influx through these channels contributes to intracellular Ca²⁺ overload. However, we showed that the homologous C. elegans DEG/ENaC channel UNC-8(d) is not Ca²⁺ permeable, yet it is neurotoxic, suggesting that Na⁺ influx is sufficient to induce cell death. Interestingly, UNC-8(d) shows small currents due to extracellular Ca²⁺ block in the Xenopus oocyte expression system. Thus, MEC-4(d) and UNC-8(d) differ both in current amplitude and Ca²⁺ permeability. Given that these two channels show a striking difference in toxicity, we wondered how Na⁺ conductance vs. Ca²⁺ permeability contributes to cell death. To address this question, we built an UNC-8/MEC-4 chimeric channel that retains the calcium permeability of MEC-4 and characterized its properties in Xenopus oocytes. Our data support the hypothesis that for Ca²⁺-permeable DEG/ENaC channels, both Ca²⁺ permeability and Na⁺ conductance contribute to toxicity. However, for Ca²⁺-impermeable DEG/ENaCs (e.g., UNC-8), our evidence shows that constitutive Na⁺ conductance is sufficient to induce toxicity, and that this effect is enhanced as current amplitude increases. Our work further refines the contribution of different channel properties to cellular toxicity induced by hyperactive DEG/ENaC channels.

Analysis of C. elegans muscle transcriptome using trans-splicing-based RNA tagging (SRT)

Current approaches to profiling tissue-specific gene expression in C. elegans require delicate manipulation and are difficult under certain conditions, e.g. from dauer or aging worms. We have developed an easy and robust method for tissue-specific RNA-seq by taking advantage of the endogenous trans-splicing process. In this method, transgenic worms are generated in which a spliced leader (SL) RNA gene is fused with a sequence tag and driven by a tissue-specific promoter. Only in the tissue of interest, the tagged SL RNA gene is transcribed and then trans-spliced onto mRNAs. The tag allows enrichment and sequencing of mRNAs from that tissue only. As a proof of principle, we profiled the muscle transcriptome, which showed high coverage and efficient enrichment of muscle specific genes, with low background noise. To demonstrate the robustness of our method, we profiled muscle gene expression in dauer larvae and aging worms, revealing gene expression changes consistent with the physiology of these stages. The resulting muscle transcriptome also revealed 461 novel RNA transcripts, likely muscle-expressed long non-coding RNAs. In summary, the splicing-based RNA tagging (SRT) method provides a convenient and robust tool to profile trans-spliced genes and identify novel transcripts in a tissue-specific manner, with a low false positive rate.

Overlapping cell population expression profiling and regulatory inference in C. elegans

Background

Understanding gene expression across the diverse metazoan cell types during development is critical to understanding their function and regulation. However, most cell types have not been assayed for expression genome-wide.

Results

We applied a novel approach we term “Profiling of Overlapping Populations of cells (POP-Seq)” to assay differential expression across all embryonic cells in the nematode Caenorhabditis elegans. In this approach, we use RNA-seq to define the transcriptome of diverse partially overlapping FACS-sorted cell populations. This identified thousands of transcripts differentially expressed across embryonic cells. Hierarchical clustering analysis identified over 100 sets of coexpressed genes corresponding to distinct patterns of cell type specific expression. We identified thousands of candidate regulators of these clusters based on enrichment of transcription factor motifs and experimentally determined binding sites.

Conclusions

Our analysis provides new insight into embryonic gene regulation, and provides a resource for improving our knowledge of tissue-specific expression and its regulation throughout C. elegans development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2482-z) contains supplementary material, which is available to authorized users.

Developmental Defects of Caenorhabditis elegans Lacking Branched-chain α-Ketoacid Dehydrogenase Are Mainly Caused by Monomethyl Branched-chain Fatty Acid Deficiency

Branched-chain α-ketoacid dehydrogenase (BCKDH) catalyzes the critical step in the branched-chain amino acid (BCAA) catabolic pathway and has been the focus of extensive studies. Mutations in the complex disrupt many fundamental metabolic pathways and cause multiple human diseases including maple syrup urine disease (MSUD), autism, and other related neurological disorders. BCKDH may also be required for the synthesis of monomethyl branched-chain fatty acids (mmBCFAs) from BCAAs. The pathology of MSUD has been attributed mainly to BCAA accumulation, but the role of mmBCFA has not been evaluated. Here we show that disrupting BCKDH in Caenorhabditis elegans causes mmBCFA deficiency, in addition to BCAA accumulation. Worms with deficiency in BCKDH function manifest larval arrest and embryonic lethal phenotypes, and mmBCFA supplementation suppressed both without correcting BCAA levels. The majority of developmental defects caused by BCKDH deficiency may thus be attributed to lacking mmBCFAs in worms. Tissue-specific analysis shows that restoration of BCKDH function in multiple tissues can rescue the defects, but is especially effective in neurons. Taken together, we conclude that mmBCFA deficiency is largely responsible for the developmental defects in the worm and conceivably might also be a critical contributor to the pathology of human MSUD.

Getting Down to Specifics: Profiling Gene Expression and Protein-DNA Interactions in a Cell Type-Specific Manner

The majority of multicellular organisms are comprised of an extraordinary range of cell types, with different properties and gene expression profiles. Understanding what makes each cell type unique and how their individual characteristics are attributed are key questions for both developmental and neurobiologists alike. The brain is an excellent example of the cellular diversity expressed in the majority of eukaryotes. The mouse brain comprises of approximately 75 million neurons varying in morphology, electrophysiology, and preferences for synaptic partners. A powerful process in beginning to pick apart the mechanisms that specify individual characteristics of the cell, as well as their fate, is to profile gene expression patterns, chromatin states, and transcriptional networks in a cell type-specific manner, i.e., only profiling the cells of interest in a particular tissue. Depending on the organism, the questions being investigated, and the material available, certain cell type-specific profiling methods are more suitable than others. This chapter reviews the approaches presently available for selecting and isolating specific cell types and evaluates their key features.

Expression of five acetylcholine receptor subunit genes in Brugia malayi adult worms

Acetylcholine receptors (AChRs) are required for body movement in parasitic nematodes and are targets of “classical” anthelmintic drugs such as levamisole and pyrantel and of newer drugs such as tribendimidine and derquantel. While neurotransmission explains the effects of these drugs on nematode movement, their effects on parasite reproduction are unexplained. The levamisole AChR type (L-AChRs) in Caenorhabditis elegans is comprised of five subunits: Cel-UNC-29, Cel-UNC-38, Cel-UNC-63, Cel-LEV-1 and Cel-LEV-8. The genome of the filarial parasite Brugia malayi contains nine AChRs subunits including orthologues of Cel-unc-29, Cel-unc-38, and Cel-unc-63. We performed in situ hybridization with RNA probes to localize the expression of five AChR genes (Bm1_35890-Bma-unc-29, Bm1_20330-Bma-unc-38, Bm1_38195-Bma-unc-63, Bm1_48815-Bma-acr-26 and Bm1_40515-Bma-acr-12) in B. malayi adult worms. Four of these genes had similar expression patterns with signals in body muscle, developing embryos, spermatogonia, uterine wall adjacent to stretched microfilariae, wall of Vas deferens, and lateral cord. Three L-AChR subunit genes (Bma-unc-29, Bma-unc-38 and Bma-unc-63) were expressed in body muscle, which is a known target of levamisole. Bma-acr-12 was co-expressed with these levamisole subunit genes in muscle, and this suggests that its protein product may form receptors with other alpha subunits. Bma-acr-26 was expressed in male muscle but not in female muscle. Strong expression signals of these genes in early embryos and gametes in uterus and testis suggest that AChRs may have a role in nervous system development of embryogenesis and spermatogenesis. This would be consistent with embryotoxic effects of drugs that target these receptors in filarial worms. Our data show that the expression of these receptor genes is tightly regulated with regard to localization in adult worms and developmental stage in embryos and gametes. These results may help to explain the broad effects of drugs that target AChRs in filarial worms.

•

Expression patterns of Brugia malayi AChR subunit genes studied by in situ hybridization.

•

All genes highly expressed in developing embryos and sperm precursors.

•

Highly expressed in the walls of uterus and Vas deferens with mature offspring.

•

Four of five genes expressed in body muscle of adult worms.

•

Expression patterns shed new light on the action of anthelmintics in filarial parasites.

Identification of Wnt Pathway Target Genes Regulating the Division and Differentiation of Larval Seam Cells and Vulval Precursor Cells in Caenorhabditis elegans

The evolutionarily conserved Wnt/β-catenin signaling pathway plays a fundamental role during metazoan development, regulating numerous processes including cell fate specification, cell migration, and stem cell renewal. Wnt ligand binding leads to stabilization of the transcriptional effector β-catenin and upregulation of target gene expression to mediate a cellular response. During larval development of the nematode Caenorhabditis elegans, Wnt/β-catenin pathways act in fate specification of two hypodermal cell types, the ventral vulval precursor cells (VPCs) and the lateral seam cells. Because little is known about targets of the Wnt signaling pathways acting during larval VPC and seam cell differentiation, we sought to identify genes regulated by Wnt signaling in these two hypodermal cell types. We conditionally activated Wnt signaling in larval animals and performed cell type-specific "mRNA tagging" to enrich for VPC and seam cell-specific mRNAs, and then used microarray analysis to examine gene expression compared to control animals. Two hundred thirty-nine genes activated in response to Wnt signaling were identified, and we characterized 50 genes further. The majority of these genes are expressed in seam and/or vulval lineages during normal development, and reduction of function for nine genes caused defects in the proper division, fate specification, fate execution, or differentiation of seam cells and vulval cells. Therefore, the combination of these techniques was successful at identifying potential cell type-specific Wnt pathway target genes from a small number of cells and at increasing our knowledge of the specification and behavior of these C. elegans larval hypodermal cells.

Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation in Caenorhabditis elegans intestine and muscles

BACKGROUND

Tissue-specific RNA plasticity broadly impacts the development, tissue identity and adaptability of all organisms, but changes in composition, expression levels and its impact on gene regulation in different somatic tissues are largely unknown. Here we developed a new method, polyA-tagging and sequencing (PAT-Seq) to isolate high-quality tissue-specific mRNA from Caenorhabditis elegans intestine, pharynx and body muscle tissues and study changes in their tissue-specific transcriptomes and 3'UTRomes.

RESULTS

We have identified thousands of novel genes and isoforms differentially expressed between these three tissues. The intestine transcriptome is expansive, expressing over 30% of C. elegans mRNAs, while muscle transcriptomes are smaller but contain characteristic unique gene signatures. Active promoter regions in all three tissues reveal both known and novel enriched tissue-specific elements, along with putative transcription factors, suggesting novel tissue-specific modes of transcription initiation. We have precisely mapped approximately 20,000 tissue-specific polyadenylation sites and discovered that about 30% of transcripts in somatic cells use alternative polyadenylation in a tissue-specific manner, with their 3'UTR isoforms significantly enriched with microRNA targets.

CONCLUSIONS

For the first time, PAT-Seq allowed us to directly study tissue specific gene expression changes in an in vivo setting and compare these changes between three somatic tissues from the same organism at single-base resolution within the same experiment. We pinpoint precise tissue-specific transcriptome rearrangements and for the first time link tissue-specific alternative polyadenylation to miRNA regulation, suggesting novel and unexplored tissue-specific post-transcriptional regulatory networks in somatic cells.

The C. elegans embryonic fate specification factor EGL-18 (GATA) is reutilized downstream of Wnt signaling to maintain a population of larval progenitor cells

In metazoans, stem cells in developing and adult tissues can divide asymmetrically to give rise to a daughter that differentiates and a daughter that retains the progenitor fate. Although the short-lived nematode C. elegans does not possess adult somatic stem cells, the lateral hypodermal seam cells behave in a similar manner: they divide once per larval stage to generate an anterior daughter that adopts a non-dividing differentiated fate and a posterior daughter that retains the seam fate and the ability to divide further. Wnt signaling pathway is known to regulate the asymmetry of these divisions and maintain the progenitor cell fate in one daughter, but how activation of the Wnt pathway accomplished this was unknown. We describe here our recent work that identified the GATA transcription factor EGL-18 as a downstream target of Wnt signaling necessary for maintenance of a progenitor population of larval seam cells. EGL-18 was previously shown to act in the initial specification of the seam cells in the embryo. Thus the acquisition of a Wnt-responsive cis-regulatory module allows an embryonic fate specification factor to be reutilized later in life downstream of a different regulator (Wnt signaling) to maintain a progenitor cell population. These results support the use of seam cell development in C. elegans as a simple model system for studying stem and progenitor cell biology.

Cell type-specific affinity purification of nuclei for chromatin profiling in whole animals

Analyzing cell differentiation during development in a complex organism requires the analysis of expression and chromatin profiles in individual cell types. Our laboratory has developed a simple and generally applicable strategy to purify specific cell types from whole organisms for simultaneous analysis of chromatin and expression. The method, termed INTACT for Isolation of Nuclei TAgged in specific Cell Types, depends on the expression of an affinity-tagged nuclear envelope protein in the cell type of interest. These nuclei can be affinity-purified from the total pool of nuclei and used as a source for RNA and chromatin. The method serves as a simple and scalable alternative to FACS sorting or laser capture microscopy to circumvent the need for expensive equipment and specialized skills. This chapter provides detailed protocols for the cell-type specific purification of nuclei from Caenorhabditis elegans.

High-throughput optical quantification of mechanosensory habituation reveals neurons encoding memory in Caenorhabditis elegans

A major goal of neuroscience studies is to identify the neurons and molecules responsible for memory. Mechanosensory habituation in Caenorhabditis elegans is a simple form of learning and memory, in which a circuit of several sensory neurons and interneurons governs behavior. However, despite the usefulness of this paradigm, there are hardly any systems for rapid and accurate behavioral genetic analysis. Here, we developed a multiplexed optical system to genetically analyze C. elegans mechanosensory habituation, and identified two interneurons involved in memory formation. The system automatically trains large populations of animals and simultaneously quantifies the behaviors of various strains by optically discriminating between transgenic and nontransgenic animals. Biochemical and cell-specific behavioral analyses indicated that phosphorylation of cyclic AMP response element-binding protein (CREB), a factor known to regulate memory allocation, was facilitated during training and this phosphorylation in AVA and AVD interneurons was required for habituation. These interneurons are a potential target for cell-specific exploration of the molecular substrates of memory.

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.

Freedom of expression: cell-type-specific gene profiling

Cell fate and behavior are results of differential gene regulation, making techniques to profile gene expression in specific cell types highly desirable. Many methods now enable investigation at the DNA, RNA and protein level. This review introduces the most recent and popular techniques, and discusses key issues influencing the choice between these such as ease, cost and applicability of information gained. Interdisciplinary collaborations will no doubt contribute further advances, including not just in single cell type but single-cell expression profiling.

DAF-16 target identification in C. elegans: past, present and future

In C. elegans, mutations in the conserved insulin/IGF-1 signaling (IIS) pathway lead to a robust extension in lifespan, improved late life health, and protection from age-related disease. These effects are mediated by the FoxO transcription factor DAF-16 which lies downstream of the IIS kinase cascade. Identifying and functionally testing DAF-16 target genes has been a focal point of ageing research for the last 10 years. Here, I review the recent advances in identifying and understanding IIS/DAF-16 targets. These studies continue to reveal the intricate nature of the IIS/DAF-16 gene regulation network and are helping us to understand the mechanisms that control lifespan. Ageing and age related disease is an area of intense public interest, and the biochemical characterization of the genes involved will be critical for identifying drugs to improve the health of our ageing population.

Simplified method for cell-specific gene expression analysis in Caenorhabditis elegans

In the neural circuit functional identities of individual neurons are mainly specified by their differential gene expression patterns. Unveiling functional roles of each neuron requires cell-specific interrogation of neural circuitry in the context of gene expressions. The mRNA tagging strategy in Caenorhabditis elegans is a powerful technique, in which cell-specific transcripts can be isolated by co-immunoprecipitating the complexes of mRNAs and epitope-tagged poly(A) binding protein (3× FLAG-PAB-1), expressed in target neurons. However, the conventional protocol requires laborious and time-consuming procedures; chromosomal integration of gene encoding 3× FLAG-PAB-1 and bleaching of obtained integrant animals for the isolation of huge amounts of synchronized animals. In this paper, we have presented a simplified methodology for cell-specific mRNA tagging analysis in C. elegans. We show that mRNA tagging was achieved using transgenic animals expressing 3× FLAG-PAB-1 as an extrachromosomal array under the control of the flp-18 promoter, without the chromosomal integration procedure. Furthermore, we successfully isolated cell-specific mRNAs from adult transgenic animals synchronously grown from eggs laid by gravid adults during a time window of 3h. This simplification facilitates the implementation of cell-specific gene expression analysis of C. elegans, which contributes to the understanding of neural circuitry at a cell-specific resolution.

Proteomic profiling of cardiac tissue by isolation of nuclei tagged in specific cell types (INTACT)

The proper dissection of the molecular mechanisms governing the specification and differentiation of specific cell types requires isolation of pure cell populations from heterogeneous tissues and whole organisms. Here, we describe a method for purification of nuclei from defined cell or tissue types in vertebrate embryos using INTACT (isolation of nuclei tagged in specific cell types). This method, previously developed in plants, flies and worms, utilizes in vivo tagging of the nuclear envelope with biotin and the subsequent affinity purification of the labeled nuclei. In this study we successfully purified nuclei of cardiac and skeletal muscle from Xenopus using this strategy. We went on to demonstrate the utility of this approach by coupling the INTACT approach with liquid chromatography-tandem mass spectrometry (LC-MS/MS) proteomic methodologies to profile proteins expressed in the nuclei of developing hearts. From these studies we have identified the Xenopus orthologs of 12 human proteins encoded by genes, which when mutated in human lead to congenital heart disease. Thus, by combining these technologies we are able to identify tissue-specific proteins that are expressed and required for normal vertebrate organ development.

Genome-wide analysis links emerin to neuromuscular junction activity in Caenorhabditis elegans

Background

Laminopathies are diseases characterized by defects in nuclear envelope structure. A well-known example is Emery-Dreifuss muscular dystrophy, which is caused by mutations in the human lamin A/C and emerin genes. While most nuclear envelope proteins are ubiquitously expressed, laminopathies often affect only a subset of tissues. The molecular mechanisms underlying these tissue-specific manifestations remain elusive. We hypothesize that different functional subclasses of genes might be differentially affected by defects in specific nuclear envelope components.

Results

Here we determine genome-wide DNA association profiles of two nuclear envelope components, lamin/LMN-1 and emerin/EMR-1 in adult Caenorhabditis elegans. Although both proteins bind to transcriptionally inactive regions of the genome, EMR-1 is enriched at genes involved in muscle and neuronal function. Deletion of either EMR-1 or LEM-2, another integral envelope protein, causes local changes in nuclear architecture as evidenced by altered association between DNA and LMN-1. Transcriptome analyses reveal that EMR-1 and LEM-2 are associated with gene repression, particularly of genes implicated in muscle and nervous system function. We demonstrate that emr-1, but not lem-2, mutants are sensitive to the cholinesterase inhibitor aldicarb, indicating altered activity at neuromuscular junctions.

Conclusions

We identify a class of elements that bind EMR-1 but do not associate with LMN-1, and these are enriched for muscle and neuronal genes. Our data support a redundant function of EMR-1 and LEM-2 in chromatin anchoring to the nuclear envelope and gene repression. We demonstrate a specific role of EMR-1 in neuromuscular junction activity that may contribute to Emery-Dreifuss muscular dystrophy in humans.

The dauer hypothesis and the evolution of parasitism: 20 years on and still going strong

How any complex trait has evolved is a fascinating question, yet the evolution of parasitism among the nematodes is arguably one of the most arresting. How did free-living nematodes cross that seemingly insurmountable evolutionary chasm between soil dwelling and survival inside another organism? Which of the many finely honed responses to the varied and harsh environments of free-living nematodes provided the material upon which natural selection could act? Although several complementary theories explain this phenomenon, I will focus on the dauer hypothesis. The dauer hypothesis posits that the arrested third-stage dauer larvae of free-living nematodes such as Caenorhabditis elegans are, due to their many physiological similarities with infective third-stage larvae of parasitic nematodes, a pre-adaptation to parasitism. If so, then a logical extension of this hypothesis is that the molecular pathways which control entry into and recovery from dauer formation by free-living nematodes in response to environmental cues have been co-opted to control the processes of infective larval arrest and activation in parasitic nematodes. The molecular machinery that controls dauer entry and exit is present in a wide range of parasitic nematodes. However, the developmental outputs of the different pathways are both conserved and divergent, not only between populations of C. elegans or between C. elegans and parasitic nematodes but also between different species of parasitic nematodes. Thus the picture that emerges is more nuanced than originally predicted and may provide insights into the evolution of such an interesting and complex trait.