A spatial and temporal map of C. elegans gene expression

Functional Optimization in Distinct Tissues and Conditions Constrains the Rate of Protein Evolution

Understanding the main determinants of protein evolution is a fundamental challenge in biology. Despite many decades of active research, the molecular and cellular mechanisms underlying the substantial variability of evolutionary rates across cellular proteins are not currently well understood. It also remains unclear how protein molecular function is optimized in the context of multicellular species and why many proteins, such as enzymes, are only moderately efficient on average. Our analysis of genomics and functional datasets reveals in multiple organisms a strong inverse relationship between the optimality of protein molecular function and the rate of protein evolution. Furthermore, we find that highly expressed proteins tend to be substantially more functionally optimized. These results suggest that cellular expression costs lead to more pronounced functional optimization of abundant proteins and that the purifying selection to maintain high levels of functional optimality significantly slows protein evolution. We observe that in multicellular species both the rate of protein evolution and the degree of protein functional efficiency are primarily affected by expression in several distinct cell types and tissues, specifically, in developed neurons with upregulated synaptic processes in animals and in young and fast-growing tissues in plants. Overall, our analysis reveals how various constraints from the molecular, cellular, and species' levels of biological organization jointly affect the rate of protein evolution and the level of protein functional adaptation.

Distinct regulatory mechanisms by the nuclear Argonautes HRDE-1 and NRDE-3 in the soma of Caenorhabditis elegans

RNA interference is a conserved silencing mechanism that depends on the generation of small RNA molecules that disrupt synthesis of their corresponding transcripts. Nuclear RNA interference is a unique process that triggers regulation through epigenetic alterations to the genome. This pathway has been extensively characterized in Caenorhabditis elegans and involves the nuclear recruitment of H3K9 histone methyltransferases by the Argonautes HRDE-1 and NRDE-3. The coordinate regulation of genetic targets by H3K9 methylation and the nuclear Argonautes is highly complex and has been mainly described based on the small RNA populations that are involved. Recent studies have also linked the nuclear RNAi pathway to the compaction of the hermaphrodite X chromosomes during dosage compensation, a mechanism that balances genetic differences between the biological sexes by repressing X chromosomes in hermaphrodites. This chromosome-wide process provides an excellent opportunity to further investigate the relationship between H3K9 methylation and the nuclear Argonautes from the perspective of the transcriptome. Our work suggests that the nuclear RNAi and the H3K9 methylation pathways each contribute to the condensation of the X chromosomes during dosage compensation but the consequences on their transcriptional output are minimal. Instead, nuclear RNAi mutants exhibit global transcriptional differences, in which HRDE-1 and NRDE-3 affect expression of their native targets through different modes of regulation and different relationships to H3K9 methylation.

ARTICLE SUMMARY

This study examines the transcriptional consequences during the disruption of the nuclear RNAi silencing mechanism in C. elegans . Through microscopy and bioinformatic work, we demonstrate that although nuclear RNAi mutants exhibit significantly decondensed X chromosomes, chromosome-wide transcriptional de-repression is not detectable. Downstream analyses further explore the global influence of the nuclear RNAi pathway, indicating that the nuclear Argonautes HRDE-1 and NRDE-3 function through two distinct mechanisms.

The Molecular Basis of Differentiation Wave Activity in Embryogenesis

As development varies greatly across the tree of life, it may seem difficult to suggest a model that proposes a single mechanism for understanding collective cell behaviors and the coordination of tissue formation. Here we propose a mechanism called differentiation waves, which unify many disparate results involving developmental systems from across the tree of life. We demonstrate how a relatively simple model of differentiation proceeds not from function-related molecular mechanisms, but from so-called differentiation waves. A phenotypic model of differentiation waves is introduced, and its relation to molecular mechanisms is proposed. These waves contribute to a differentiation tree, which is an alternate way of viewing cell lineage and local action of the molecular factors. We construct a model of differentiation wave-related molecular mechanisms (genome, epigenome, and proteome) based on bioinformatic data from the nematode Caenorhabditis elegans. To validate this approach across different modes of development, we evaluate protein expression across different types of development by comparing Caenorhabditis elegans with several model organisms: fruit flies (Drosophila melanogaster), yeast (Saccharomyces cerevisiae), and mouse (Mus musculus). Inspired by gene regulatory networks, two Models of Interactive Contributions (fully-connected MICs and ordered MICs) are used to suggest potential genomic contributions to differentiation wave-related proteins. This, in turn, provides a framework for understanding differentiation and development.

Forward genetic screen of the C. elegans million mutation library reveals essential, cell-autonomous contributions of BBSome proteins to dopamine signaling

The nematode Caenorhabditis elegans is well known for its ability to support forward genetic screens to identify molecules involved in neuronal viability and signaling. The proteins involved in C. elegans dopamine (DA) regulation are highly conserved across evolution, with prior work demonstrating that the model can serve as an efficient platform to identify novel genes involved in disease-associated processes. To identify novel players in DA signaling, we took advantage of a recently developed library of pre-sequenced mutant nematodes arising from the million mutation project (MMP) to identify strains that display the DA-dependent swimming-induced-paralysis phenotype (Swip). Our screen identified novel mutations in the dopamine transporter encoding gene dat-1, whose loss was previously used to identify the Swip phenotype, as well as multiple genes with previously unknown connections to DA signaling. Here, we present our isolation and characterization of one of these genes, bbs-1, previously linked to the function of primary cilia in worms and higher organisms, including humans, and where loss-of-function mutations result in a human disorder known as Bardet-Biedl syndrome. Our studies of C. elegans BBS-1 protein, as well as other proteins that are known to be assembled into a higher order complex (the BBSome) reveal that functional or structural disruption of this complex leads to exaggerated C. elegans DA signaling to produce Swip via a cell-autonomous mechanism. We provide evidence that not only does the proper function of cilia in C. elegans DA neurons support normal swimming behavior, but also that bbs-1 maintains normal levels of DAT-1 trafficking or function via a RHO-1 and SWIP-13/MAPK-15 dependent pathway where mutants may contribute to Swip independent of altered ciliary function. Together, these studies demonstrate novel contributors to DA neuron function in the worm and demonstrate the utility and efficiency of forward genetic screens using the MMP library.

Transcriptomic analysis of the spatiotemporal axis of oogenesis and fertilization in C. elegans

Caenorhabditis elegans hermaphrodite presents a unique model to study the formation of oocytes. However, the size of the model animal and difficulties in retrieval of specific stages of the germline have obviated closer systematic studies of this process throughout the years. Here, we present a transcriptomic level analysis into the oogenesis of C. elegans hermaphrodites. We dissected a hermaphrodite gonad into seven sections corresponding to the mitotic distal region, the pachytene region, the diplotene region, the early diakinesis region and the 3 most proximal oocytes, and deeply sequenced the transcriptome of each of them along with that of the fertilized egg using a single-cell RNA-seq (scRNA-seq) protocol. We identified specific gene expression events as well as gene splicing events in finer detail along the gonad and provided novel insights into underlying mechanisms of the oogenesis process. Furthermore, through careful review of relevant research literature coupled with patterns observed in our analysis, we delineate transcripts that may serve functions in the interactions between the germline and cells of the somatic gonad. These results expand our knowledge of the transcriptomic space of the C. elegans germline and lay a foundation on which future studies of the germline can be based upon.

XOL-1 regulates developmental timing by modulating the H3K9 landscape in C. elegans early embryos

Sex determination in the nematode C. elegans is controlled by the master regulator XOL-1 during embryogenesis. Expression of xol-1 is dependent on the ratio of X chromosomes and autosomes, which differs between XX hermaphrodites and XO males. In males, xol-1 is highly expressed and in hermaphrodites, xol-1 is expressed at very low levels. XOL-1 activity is known to be critical for the proper development of C. elegans males, but its low expression was considered to be of minimal importance in the development of hermaphrodite embryos. Our study reveals that XOL-1 plays an important role as a regulator of developmental timing during hermaphrodite embryogenesis. Using a combination of imaging and bioinformatics techniques, we found that hermaphrodite embryos have an accelerated rate of cell division, as well as a more developmentally advanced transcriptional program when xol-1 is lost. Further analyses reveal that XOL-1 is responsible for regulating the timing of initiation of dosage compensation on the X chromosomes, and the appropriate expression of sex-biased transcriptional programs in hermaphrodites. We found that xol-1 mutant embryos overexpress the H3K9 methyltransferase MET-2 and have an altered H3K9me landscape. Some of these effects of the loss of xol-1 gene were reversed by the loss of met-2. These findings demonstrate that XOL-1 plays an important role as a developmental regulator in embryos of both sexes, and that MET-2 acts as a downstream effector of XOL-1 activity in hermaphrodites.

Mutations in nucleotide metabolism genes bypass proteasome defects in png-1/NGLY1-deficient Caenorhabditis elegans

The conserved SKN-1A/Nrf1 transcription factor regulates the expression of proteasome subunit genes and is essential for maintenance of adequate proteasome function in animal development, aging, and stress responses. Unusual among transcription factors, SKN-1A/Nrf1 is a glycoprotein synthesized in the endoplasmic reticulum (ER). N-glycosylated SKN-1A/Nrf1 exits the ER and is deglycosylated in the cytosol by the PNG-1/NGLY1 peptide:N-glycanase. Deglycosylation edits the protein sequence of SKN-1A/Nrf1 by converting N-glycosylated asparagine residues to aspartate, which is necessary for SKN-1A/Nrf1 transcriptional activation of proteasome subunit genes. Homozygous loss-of-function mutations in the peptide:N-glycanase (NGLY1) gene cause NGLY1 deficiency, a congenital disorder of deglycosylation. There are no effective treatments for NGLY1 deficiency. Since SKN-1A/Nrf1 is a major client of NGLY1, the resulting proteasome deficit contributes to NGLY1 disease. We sought to identify targets for mitigation of proteasome dysfunction in NGLY1 deficiency that might indicate new avenues for treatment. We isolated mutations that suppress the sensitivity to proteasome inhibitors caused by inactivation of the NGLY1 ortholog PNG-1 in Caenorhabditis elegans. We identified multiple suppressor mutations affecting 3 conserved genes: rsks-1, tald-1, and ent-4. We show that the suppressors act through a SKN-1/Nrf-independent mechanism and confer proteostasis benefits consistent with amelioration of proteasome dysfunction. ent-4 encodes an intestinal nucleoside/nucleotide transporter, and we show that restriction of nucleotide availability is beneficial, whereas a nucleotide-rich diet exacerbates proteasome dysfunction in PNG-1/NGLY1-deficient C. elegans. Our findings suggest that dietary or pharmacological interventions altering nucleotide availability have the potential to mitigate proteasome insufficiency in NGLY1 deficiency and other diseases associated with proteasome dysfunction.

Transcriptomic Analysis of the Spatiotemporal Axis of Oogenesis and Fertilization in C. elegans

The oocyte germline of the C. elegans hermaphrodite presents a unique model to study the formation of oocytes. However, the size of the model animal and difficulties in retrieval of specific stages of the germline have obviated closer systematic studies of this process throughout the years. Here, we present a transcriptomic level analysis into the oogenesis of C. elegans hermaphrodites. We dissected a hermaphrodite gonad into seven sections corresponding to the mitotic distal region, the pachytene, the diplotene, the early diakinesis region and the 3 most proximal oocytes, and deeply sequenced the transcriptome of each of them along with that of the fertilized egg using a single-cell RNA-seq protocol. We identified specific gene expression events as well as gene splicing events in finer detail along the oocyte germline and provided novel insights into underlying mechanisms of the oogenesis process. Furthermore, through careful review of relevant research literature coupled with patterns observed in our analysis, we attempt to delineate transcripts that may serve functions in the interaction between the germline and cells of the somatic gonad. These results expand our knowledge of the transcriptomic space of the C. elegans germline and lay a foundation on which future studies of the germline can be based upon.

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

The functional small RNA interactome reveals targets for the vancomycin-responsive sRNA RsaOI in vancomycin-tolerant Staphylococcus aureus

Small RNAs have been found to control a broad range of bacterial phenotypes including tolerance to antibiotics. Vancomycin tolerance in multidrug resistance Staphylococcus aureus is correlated with dysregulation of small RNAs although their contribution to antibiotic tolerance is poorly understood. RNA-RNA interactome profiling techniques are expanding our understanding of sRNA-mRNA interactions in bacteria; however, determining the function of these interactions for hundreds of sRNA-mRNA pairs is a major challenge. At steady-state, protein and mRNA abundances are often highly correlated and lower than expected protein abundance may indicate translational repression of an mRNA. To identify sRNA-mRNA interactions that regulate mRNA translation, we examined the correlation between gene transcript abundance, ribosome occupancy, and protein levels. We used the machine learning technique self-organizing maps (SOMs) to cluster genes with similar transcription and translation patterns and identified a cluster of mRNAs that appeared to be post-transcriptionally repressed. By integrating our clustering with sRNA-mRNA interactome data generated in vancomycin-tolerant S. aureus by RNase III-CLASH, we identified sRNAs that may be mediating translational repression. We have confirmed sRNA-dependant post-transcriptional repression of several mRNAs in this cluster. Two of these interactions are mediated by RsaOI, a sRNA that is highly upregulated by vancomycin. We demonstrate the regulation of HPr and the cell-wall autolysin Atl. These findings suggest that RsaOI coordinates carbon metabolism and cell wall turnover during vancomycin treatment.

IMPORTANCE

The emergence of multidrug-resistant Staphylococcus aureus (MRSA) is a major public health concern. Current treatment is dependent on the efficacy of last-line antibiotics like vancomycin. The most common cause of vancomycin treatment failure is strains with intermediate resistance or tolerance that arise through the acqusition of a diverse repertoire of point mutations. These strains have been shown to altered small RNA (sRNA) expression in response to antibiotic treatment. Here, we have used a technique termed RNase III-CLASH to capture sRNA interactions with their target mRNAs. To understand the function of these interactions, we have looked at RNA and protein abundance for mRNAs targeted by sRNAs. Messenger RNA and protein levels are generally well correlated and we use deviations from this correlation to infer post-transcriptional regulation and the function of individual sRNA-mRNA interactions. Using this approach we identify mRNA targets of the vancomycin-induced sRNA, RsaOI, that are repressed at the translational level. We find that RsaOI represses the cell wall autolysis Atl and carbon transporter HPr suggestion a link between vancomycin treatment and suppression of cell wall turnover and carbon metabolism.

Spatiotemporal analysis of mRNA-protein relationships enhances transcriptome-based developmental inference

Elucidating the complex relationships between mRNA and protein expression at high spatiotemporal resolution is critical for unraveling multilevel gene regulation and enhancing mRNA-based developmental analyses. In this study, we conduct a single-cell analysis of mRNA and protein expression of transcription factors throughout C. elegans embryogenesis. Initially, cellular co-presence of mRNA and protein is low, increasing to a medium-high level (73%) upon factoring in delayed protein synthesis and long-term protein persistence. These factors substantially affect mRNA-protein concordance, leading to potential inaccuracies in mRNA-reliant gene detection and specificity characterization. Building on the learned relationship, we infer protein presence from mRNA expression and demonstrate its utility in identifying tissue-specific genes and elucidating relationships between genes and cells. This approach facilitates identifying the role of sptf-1/SP7 in neuronal lineage development. Collectively, this study provides insights into gene expression dynamics during rapid embryogenesis and approaches for improving the efficacy of transcriptome-based developmental analyses.

Automated cell annotation in multi-cell images using an improved CRF_ID algorithm

Cell identification is an important yet difficult process in data analysis of biological images. Previously, we developed an automated cell identification method called CRF_ID and demonstrated its high performance in C. elegans whole-brain images (Chaudhary et al, 2021). However, because the method was optimized for whole-brain imaging, comparable performance could not be guaranteed for application in commonly used C. elegans multi-cell images that display a subpopulation of cells. Here, we present an advance CRF_ID 2.0 that expands the generalizability of the method to multi-cell imaging beyond whole-brain imaging. To illustrate the application of the advance, we show the characterization of CRF_ID 2.0 in multi-cell imaging and cell-specific gene expression analysis in C. elegans. This work demonstrates that high accuracy automated cell annotation in multi-cell imaging can expedite cell identification and reduce its subjectivity in C. elegans and potentially other biological images of various origins.

The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

Spaceflight Induces Strength Decline in Caenorhabditis elegans

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

VISTA: visualizing the spatial transcriptome of the C.elegans nervous system

Summary

Profiling the transcriptomes of single cells without sacrificing spatial information is a major goal of the field of spatial transcriptomics, but current technologies require tradeoffs between single-cell resolution and whole-transcriptome coverage. In one animal species, the nematode worm Caenorhabditis elegans, a comprehensive spatial transcriptome with single-cell resolution is attainable using existing datasets, thanks to the worm's invariant cell lineage and a series of recently generated single cell transcriptomes. Here we present VISTA, which leverages these datasets to provide a visualization of the worm spatial transcriptome, focusing specifically on the nervous system. VISTA allows users to input a query gene and visualize its expression across all neurons in the form of a "spatial heatmap" in which the color of a cell reports the expression level. Underlying gene expression values (in Transcripts Per Million) are displayed when an individual cell is selected. We provide examples of the utility of VISTA for identifying striking new gene expression patterns in specific neurons, and for resolving cellular identities of ambiguous expression patterns generated from in vivo reporter genes. The ability to easily obtain gene-level snapshots of the neuronal spatial transcriptome should facilitate studies on neuron-specific gene expression and regulation and provide a template for the high-resolution spatial transcriptomes the field hopes to obtain for various animal species in the future.

Availability and implementation

VISTA is freely available at the following URL: https://public.tableau.com/app/profile/smu.oit.data.insights/viz/VISTA_16814210566130/VISTA.

azyx-1 is a new gene that overlaps with zyxin and affects its translation in C. elegans, impacting muscular integrity and locomotion

Overlapping genes are widely prevalent; however, their expression and consequences are poorly understood. Here, we describe and functionally characterize a novel zyx-1 overlapping gene, azyx-1, with distinct regulatory functions in Caenorhabditis elegans. We observed conservation of alternative open reading frames (ORFs) overlapping the 5' region of zyxin family members in several animal species, and find shared sites of azyx-1 and zyxin proteoform expression in C. elegans. In line with a standard ribosome scanning model, our results support cis regulation of zyx-1 long isoform(s) by upstream initiating azyx-1a. Moreover, we report on a rare observation of trans regulation of zyx-1 by azyx-1, with evidence of increased ZYX-1 upon azyx-1 overexpression. Our results suggest a dual role for azyx-1 in influencing zyx-1 proteoform heterogeneity and highlight its impact on C. elegans muscular integrity and locomotion.

Transcriptome profiling of the Caenorhabditis elegans intestine reveals that ELT-2 negatively and positively regulates intestinal gene expression within the context of a gene regulatory network

ELT-2 is the major transcription factor (TF) required for Caenorhabditis elegans intestinal development. ELT-2 expression initiates in embryos to promote development and then persists after hatching through the larval and adult stages. Though the sites of ELT-2 binding are characterized and the transcriptional changes that result from ELT-2 depletion are known, an intestine-specific transcriptome profile spanning developmental time has been missing. We generated this dataset by performing Fluorescence Activated Cell Sorting on intestine cells at distinct developmental stages. We analyzed this dataset in conjunction with previously conducted ELT-2 studies to evaluate the role of ELT-2 in directing the intestinal gene regulatory network through development. We found that only 33% of intestine-enriched genes in the embryo were direct targets of ELT-2 but that number increased to 75% by the L3 stage. This suggests additional TFs promote intestinal transcription especially in the embryo. Furthermore, only half of ELT-2's direct target genes were dependent on ELT-2 for their proper expression levels, and an equal proportion of those responded to elt-2 depletion with over-expression as with under-expression. That is, ELT-2 can either activate or repress direct target genes. Additionally, we observed that ELT-2 repressed its own promoter, implicating new models for its autoregulation. Together, our results illustrate that ELT-2 impacts roughly 20-50% of intestine-specific genes, that ELT-2 both positively and negatively controls its direct targets, and that the current model of the intestinal regulatory network is incomplete as the factors responsible for directing the expression of many intestinal genes remain unknown.

Aminopeptidase MNP-1 triggers intestine protease production by activating daf-16 nuclear location to degrade pore-forming toxins in Caenorhabditis elegans

Pore-forming toxins (PFTs) are effective tools for pathogens infection. By disrupting epithelial barriers and killing immune cells, PFTs promotes the colonization and reproduction of pathogenic microorganisms in their host. In turn, the host triggers defense responses, such as endocytosis, exocytosis, or autophagy. Bacillus thuringiensis (Bt) bacteria produce PFT, known as crystal proteins (Cry) which damage the intestinal cells of insects or nematodes, eventually killing them. In insects, aminopeptidase N (APN) has been shown to act as an important receptor for Cry toxins. Here, using the nematode Caenorhabditis elegans as model, an extensive screening of APN gene family was performed to analyze the potential role of these proteins in the mode of action of Cry5Ba against the nematode. We found that one APN, MNP-1, participate in the toxin defense response, since the mnp-1(ok2434) mutant showed a Cry5Ba hypersensitive phenotype. Gene expression analysis in mnp-1(ok2434) mutant revealed the involvement of two protease genes, F19C6.4 and R03G8.6, that participate in Cry5Ba degradation. Finally, analysis of the transduction pathway involved in F19C6.4 and R03G8.6 expression revealed that upon Cry5Ba exposure, the worms up regulated both protease genes through the activation of the FOXO transcription factor DAF-16, which was translocated into the nucleus. The nuclear location of DAF-16 was found to be dependent on mnp-1 under Cry5Ba treatment. Our work provides evidence of new host responses against PFTs produced by an enteric pathogenic bacterium, resulting in activation of host intestinal proteases that degrade the PFT in the intestine.

VISTA: Visualizing the Spatial Transcriptome of the C. elegans Nervous System

Profiling the transcriptomes of single cells without sacrificing spatial information is a major goal of the field of spatial transcriptomics, but current technologies require tradeoffs between single-cell resolution and whole-transcriptome coverage. In one animal species, the nematode worm C. elegans, a comprehensive spatial transcriptome with single-cell resolution is attainable using existing datasets, thanks to the worm's invariant cell lineage and a series of recently-generated single cell transcriptomes. Here we present VISTA, which leverages these datasets to provide a visualization of the worm spatial transcriptome, focusing specifically on the nervous system. VISTA allows users to input a query gene and visualize its expression across all neurons in the form of a "spatial heatmap" in which the color of a cell reports the expression level. Underlying gene expression values (in Transcripts Per Million) are displayed when an individual cell is selected. We provide examples of the utility of VISTA for identifying striking new gene expression patterns in specific neurons, and for resolving cellular identities of ambiguous expression patterns generated from in vivo reporter genes. The ability to easily obtain gene-level snapshots of the neuronal spatial transcriptome should facilitate studies on neuron-specific gene expression and regulation, and provide a template for the high-resolution spatial transcriptomes the field hopes to obtain for various animal species in the future.

Temporal transitions in the postembryonic nervous system of the nematode Caenorhabditis elegans: Recent insights and open questions

After the generation, differentiation and integration into functional circuitry, post-mitotic neurons continue to change certain phenotypic properties throughout postnatal juvenile stages until an animal has reached a fully mature state in adulthood. We will discuss such changes in the context of the nervous system of the nematode C. elegans, focusing on recent descriptions of anatomical and molecular changes that accompany postembryonic maturation of neurons. We summarize the characterization of genetic timer mechanisms that control these temporal transitions or maturational changes, and discuss that many but not all of these transitions relate to sexual maturation of the animal. We describe how temporal, spatial and sex-determination pathways are intertwined to sculpt the emergence of cell-type specific maturation events. Finally, we lay out several unresolved questions that should be addressed to move the field forward, both in C. elegans and in vertebrates.

Utilization of Computer Classification Methods for Exposure Prediction and Gene Selection in Daphnia magna Toxicogenomics

Zinc (Zn) is an essential element that influences many cellular functions. Depending on bioavailability, Zn can cause both deficiency and toxicity. Zn bioavailability is influenced by water hardness. Therefore, water quality analysis for health-risk assessment should consider both Zn concentration and water hardness. However, exposure media selection for traditional toxicology tests are set to defined hardness levels and do not represent the diverse water chemistry compositions observed in nature. Moreover, these tests commonly use whole organism endpoints, such as survival and reproduction, which require high numbers of test animals and are labor intensive. Gene expression stands out as a promising alternative to provide insight into molecular events that can be used for risk assessment. In this work, we apply machine learning techniques to classify the Zn concentrations and water hardness from Daphnia magna gene expression by using quantitative PCR. A method for gene ranking was explored using techniques from game theory, namely, Shapley values. The results show that standard machine learning classifiers can classify both Zn concentration and water hardness simultaneously, and that Shapley values are a versatile and useful alternative for gene ranking that can provide insight about the importance of individual genes.

Systems-level transcriptional regulation of Caenorhabditis elegans metabolism

Metabolism is controlled to ensure organismal development and homeostasis. Several mechanisms regulate metabolism, including allosteric control and transcriptional regulation of metabolic enzymes and transporters. So far, metabolism regulation has mostly been described for individual genes and pathways, and the extent of transcriptional regulation of the entire metabolic network remains largely unknown. Here, we find that three-quarters of all metabolic genes are transcriptionally regulated in the nematode Caenorhabditis elegans. We find that many annotated metabolic pathways are coexpressed, and we use gene expression data and the iCEL1314 metabolic network model to define coregulated subpathways in an unbiased manner. Using a large gene expression compendium, we determine the conditions where subpathways exhibit strong coexpression. Finally, we develop "WormClust," a web application that enables a gene-by-gene query of genes to view their association with metabolic (sub)-pathways. Overall, this study sheds light on the ubiquity of transcriptional regulation of metabolism and provides a blueprint for similar studies in other organisms, including humans.

MJL-1 is a nuclear envelope protein required for homologous chromosome pairing and regulation of synapsis during meiosis in C. elegans

Interactions between chromosomes and LINC (linker of nucleoskeleton and cytoskeleton) complexes in the nuclear envelope (NE) promote homolog pairing and synapsis during meiosis. By tethering chromosomes to cytoskeletal motors, these connections lead to processive chromosome movements along the NE. This activity is usually mediated by telomeres, but in the nematode Caenorhabditis elegans, special chromosome regions called "pairing centers" (PCs) have acquired this meiotic function. Here, we identify a previously uncharacterized meiosis-specific NE protein, MJL-1 (MAJIN-Like-1), that is essential for interactions between PCs and LINC complexes in C. elegans. Mutations in MJL-1 eliminate active chromosome movements during meiosis, resulting in nonhomologous synapsis and impaired homolog pairing. Fission yeast and mice also require NE proteins to connect chromosomes to LINC complexes. Extensive similarities in the molecular architecture of meiotic chromosome-NE attachments across eukaryotes suggest a common origin and/or functions of this architecture during meiosis.

A new AMPK isoform mediates glucose-restriction induced longevity non-cell autonomously by promoting membrane fluidity

Dietary restriction (DR) delays aging and the onset of age-associated diseases. However, it is yet to be determined whether and how restriction of specific nutrients promote longevity. Previous genome-wide screens isolated several Escherichia coli mutants that extended lifespan of Caenorhabditis elegans. Here, using ¹H-NMR metabolite analyses and inter-species genetics, we demonstrate that E. coli mutants depleted of intracellular glucose extend C. elegans lifespans, serving as bona fide glucose-restricted (GR) diets. Unlike general DR, GR diets don't reduce the fecundity of animals, while still improving stress resistance and ameliorating neuro-degenerative pathologies of Aβ₄₂. Interestingly, AAK-2a, a new AMPK isoform, is necessary and sufficient for GR-induced longevity. AAK-2a functions exclusively in neurons to modulate GR-mediated longevity via neuropeptide signaling. Last, we find that GR/AAK-2a prolongs longevity through PAQR-2/NHR-49/Δ9 desaturases by promoting membrane fluidity in peripheral tissues. Together, our studies identify the molecular mechanisms underlying prolonged longevity by glucose specific restriction in the context of whole animals.

Transcriptomic Signature of the Simulated Microgravity Response in Caenorhabditis elegans and Comparison to Spaceflight Experiments

Given the growing interest in human exploration of space, it is crucial to identify the effects of space conditions on biological processes. Here, we analyze the transcriptomic response of Caenorhabditis elegans to simulated microgravity and observe the maintained transcriptomic response after returning to ground conditions for four, eight, and twelve days. We show that 75% of the simulated microgravity-induced changes on gene expression persist after returning to ground conditions for four days while most of these changes are reverted after twelve days. Our results from integrative RNA-seq and mass spectrometry analyses suggest that simulated microgravity affects longevity-regulating insulin/IGF-1 and sphingolipid signaling pathways. Finally, we identified 118 genes that are commonly differentially expressed in simulated microgravity- and space-exposed worms. Overall, this work provides insight into the effect of microgravity on biological systems during and after exposure.

Genetic polymorphism, constitutive expression and tissue localization of Dirofilaria immitis P-glycoprotein 11: a putative marker of macrocyclic lactone resistance

BACKGROUND

Dirofilaria immitis causes dirofilariosis, a potentially fatal condition in canids. Dirofilaria infections can be prevented with a macrocyclic lactone (ML) prophylactic regimen. However, some D. immitis isolates have become resistant to MLs. Genetic changes on the P-glycoprotein 11 gene, encoding an ABCB transporter, have been linked to the ML-resistant phenotypes and have been proposed as markers of drug resistance. However, nothing is known about the expression and the localization of this transporter in D. immitis, despite its strong link to ML-resistant phenotypes.

METHODS

We examined the clinically validated D. immitis P-glycoprotein 11 (DimPgp-11) single nucleotide polymorphism (SNP) via MiSeq analysis in three ML-susceptible isolates (Missouri, MP3 and Yazoo) and two ML-resistant isolates (JYD-34 and Metairie), and correlated the data with previously published MiSeq results of USA laboratory-maintained D. immitis isolates. The level of the expression of the DimPgp-11 messenger RNA transcript was analyzed by droplet digital PCR (ddPCR) and compared in the USA laboratory-maintained isolates, namely the ML-susceptible Missouri and Berkeley isolates, the putative ML-susceptible Georgia III and Big Head isolates and the ML-resistant isolate JYD-34. The immunolocalization of DimPgp-11 was visualized in the microfilaria (mf) life stage of the Missouri isolate using confocal microscopy.

RESULTS

The results confirmed that the SNP found on DimPgp-11 is differentially expressed in the USA laboratory-maintained isolates. The ML-susceptible isolates had an alternate allele frequency of between 0% and 15%, while it ranged between 17% and 56% in the ML-resistant isolates. The constitutive expression of DimPgp-11 was similar in the Berkeley, Georgia III and Big Head isolates, while it was significantly decreased in the ML-resistant JYD-34 isolate (P < 0.05), when compared to the ML-susceptible Missouri isolate. The DimPgp-11 protein was distinctly localized within the excretory-secretory (ES) duct, pore cells and the excretory cell and, more faintly, along the mf body wall.

CONCLUSION

Our data confirm that genetic polymorphism of DimPgp-11 is associated with ML resistance in USA laboratory-maintained D. imminits isolates. A link between DimPgp-11 and ML resistance in D. immitis is further supported by the lower protein expression in the ML-resistant JYD-34 isolate when compared with the ML-susceptible Missouri isolate. Interestingly, DimPgp-11 is strategically located surrounding the ES pore where it could play an active role in ML efflux.

Unconventional Myosins from Caenorhabditis elegans as a Probe to Study Human Orthologues

Unconventional myosins are a superfamily of actin-based motor proteins that perform a number of roles in fundamental cellular processes, including (but not limited to) intracellular trafficking, cell motility, endocytosis, exocytosis and cytokinesis. 40 myosins genes have been identified in humans, which belong to different 12 classes based on their domain structure and organisation. These genes are widely expressed in different tissues, and mutations leading to loss of function are associated with a wide variety of pathologies while over-expression often results in cancer. Caenorhabditis elegans (C. elegans) is a small, free-living, non-parasitic nematode. ~38% of the genome of C. elegans has predicted orthologues in the human genome, making it a valuable tool to study the function of human counterparts and human diseases. To date, 8 unconventional myosin genes have been identified in the nematode, from 6 different classes with high homology to human paralogues. The hum-1 and hum-5 (heavy chain of an unconventional myosin) genes encode myosin of class I, hum-2 of class V, hum-3 and hum-8 of class VI, hum-6 of class VII and hum-7 of class IX. The hum-4 gene encodes a high molecular mass myosin (307 kDa) that is one of the most highly divergent myosins and is a member of class XII. Mutations in many of the human orthologues are lethal, indicating their essential properties. However, a functional characterisation for many of these genes in C. elegans has not yet been performed. This article reviews the current knowledge of unconventional myosin genes in C. elegans and explores the potential use of the nematode to study the function and regulation of myosin motors to provide valuable insights into their role in diseases.

TGS1 impacts snRNA 3'-end processing, ameliorates survival motor neuron-dependent neurological phenotypes in vivo and prevents neurodegeneration

Trimethylguanosine synthase 1 (TGS1) is a highly conserved enzyme that converts the 5'-monomethylguanosine cap of small nuclear RNAs (snRNAs) to a trimethylguanosine cap. Here, we show that loss of TGS1 in Caenorhabditis elegans, Drosophila melanogaster and Danio rerio results in neurological phenotypes similar to those caused by survival motor neuron (SMN) deficiency. Importantly, expression of human TGS1 ameliorates the SMN-dependent neurological phenotypes in both flies and worms, revealing that TGS1 can partly counteract the effects of SMN deficiency. TGS1 loss in HeLa cells leads to the accumulation of immature U2 and U4atac snRNAs with long 3' tails that are often uridylated. snRNAs with defective 3' terminations also accumulate in Drosophila Tgs1 mutants. Consistent with defective snRNA maturation, TGS1 and SMN mutant cells also exhibit partially overlapping transcriptome alterations that include aberrantly spliced and readthrough transcripts. Together, these results identify a neuroprotective function for TGS1 and reinforce the view that defective snRNA maturation affects neuronal viability and function.

Channel-independent function of UNC-9/Innexin in spatial arrangement of GABAergic synapses in C. elegans

Precise synaptic connection of neurons with their targets is essential for the proper functioning of the nervous system. A plethora of signaling pathways act in concert to mediate the precise spatial arrangement of synaptic connections. Here we show a novel role for a gap junction protein in controlling tiled synaptic arrangement in the GABAergic motor neurons in Caenorhabditis elegans, in which their axons and synapses overlap minimally with their neighboring neurons within the same class. We found that while EGL-20/Wnt controls axonal tiling, their presynaptic tiling is mediated by a gap junction protein UNC-9/Innexin, that is localized at the presynaptic tiling border between neighboring dorsal D-type GABAergic motor neurons. Strikingly, the gap junction channel activity of UNC-9 is dispensable for its function in controlling tiled presynaptic patterning. While gap junctions are crucial for the proper functioning of the nervous system as channels, our finding uncovered the novel channel-independent role of UNC-9 in synapse patterning.

Germline immortality relies on TRIM32-mediated turnover of a maternal mRNA activator in C. elegans

How the germ line achieves a clean transition from maternal to zygotic gene expression control is a fundamental problem in sexually reproducing organisms. Whereas several mechanisms terminate the maternal program in the soma, this combined molecular reset and handover are poorly understood for primordial germ cells (PGCs). Here, we show that GRIF-1, a TRIM32-related and presumed E3 ubiquitin ligase in Caenorhabditis elegans, eliminates the maternal cytoplasmic poly(A) polymerase (cytoPAP) complex by targeting the germline-specific intrinsically disordered region of its enzymatic subunit, GLD-2, for proteasome-mediated degradation. Interference with cytoPAP turnover in PGCs causes frequent transgenerational sterility and, eventually, germline mortality. Hence, positively acting maternal RNA regulators are cleared via the proteasome system to avoid likely interference between maternal and zygotic gene expression programs to maintain transgenerational fertility and acquire germline immortality. This strategy is likely used in all animals that preform their immortal germ line via maternally inherited germplasm determinants.

Teaching the Applications of CRISPR/Cas9: Using the African Turquoise Killifish as a Novel Model of Aging and Age-Related Diseases

The development of genome editing technologies, including the novel CRISPR/Cas9 technique, has advanced scientific research concerning the contribution of genetics to disease through the creation of new model organisms. The subject of this review is a 2015 study done by Harel et al. from the journal Cell. This study is a prime example of using CRISPR/Cas9 to create a new model organism to accurately model the effects of aging and age-related diseases on a short-lived vertebrate. This study found that the African turquoise killifish is a reliable model to study the physiological process of aging due to its compressed lifespan. In addition, it provides a genotype-to-phenotype platform to study genes related to the hallmarks of aging and age-related diseases. This paper demonstrates this by showing that killifish deficient in the protein subunit of telomerase display telomerase-related pathologies faster than other established vertebrate models. From a teaching perspective, this paper could be used as a resource for educators to teach students about new technologies emerging in the field of neuroscience and the importance of model organisms. Specifically, for upper-level undergraduate students, this paper could serve as a real-world example of how scientific techniques such as CRISPR/Cas9 could be used to answer scientific questions. Further, it shows how these techniques could bring forward new model organisms better suited to answer the scientific questions being asked. Learning these techniques and being open minded to new approaches will be advantageous to students' future careers in science.

Using single-worm RNA sequencing to study C. elegans responses to pathogen infection

Background

Caenorhabditis elegans is an excellent research model whose populations have been used in many studies to address various biological questions. Although worm-to-worm phenotypic variations in isogenic populations have been persistently observed, they are not well understood and are often ignored or averaged out in studies, masking the impacts of such variations on data collection and interpretation. Single-worm RNA sequencing that profiles the transcriptomes of individual animals has the power to examine differences between individuals in a worm population, but this approach has been understudied. The integrity of the starting RNA, the quality of the library and sequence data, as well as the transcriptome-profiling effectiveness of single-worm RNA-seq remain unclear. Therefore, more studies are needed to improve this technique and its application in research.

Results

In this study, we aimed to develop a single-worm RNA-seq method that includes five steps: worm lysis and RNA extraction, cDNA synthesis, library preparation, sequencing, and sequence data analysis. We found that the mechanical lysis of worms using a Qiagen TissueLyser maintained RNA integrity and determined that the quality of our single-worm libraries was comparable to that of standard RNA-seq libraries based on assessments of a variety of parameters. Furthermore, analysis of pathogen infection-induced gene expression using single-worm RNA-seq identified a core set of genes and biological processes relating to the immune response and metabolism affected by infection. These results demonstrate the effectiveness of our single-worm RNA-seq method in transcriptome profiling and its usefulness in addressing biological questions.

Conclusions

We have developed a single-worm RNA-seq method to effectively profile gene expression in individual C. elegans and have applied this method to study C. elegans responses to pathogen infection. Key aspects of our single-worm RNA-seq libraries were comparable to those of standard RNA-seq libraries. The single-worm method captured the core set of, but not all, infection-affected genes and biological processes revealed by the standard method, indicating that there was gene regulation that is not shared by all individuals in a population. Our study suggests that combining single-worm and standard RNA-seq approaches will allow for detecting and distinguishing shared and individual-specific gene activities in isogenic populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08878-x.

Maternal H3K36 and H3K27 HMTs protect germline development via regulation of the transcription factor LIN-15B

Maternally synthesized products play critical roles in the development of offspring. A premier example is the Caenorhabditis elegans H3K36 methyltransferase MES-4, which is essential for germline survival and development in offspring. How maternal MES-4 protects the germline is not well understood, but its role in H3K36 methylation hinted that it may regulate gene expression in primordial germ cells (PGCs). We tested this hypothesis by profiling transcripts from nascent germlines (PGCs and their descendants) dissected from wild-type and mes-4 mutant (lacking maternal and zygotic MES-4) larvae. mes-4 nascent germlines displayed downregulation of some germline genes, upregulation of some somatic genes, and dramatic upregulation of hundreds of genes on the X chromosome. We demonstrated that upregulation of one or more genes on the X is the cause of germline death by generating and analyzing mes-4 mutants that inherited different endowments of X chromosome(s). Intriguingly, removal of the THAP transcription factor LIN-15B from mes-4 mutants reduced X misexpression and prevented germline death. lin-15B is X-linked and misexpressed in mes-4 PGCs, identifying it as a critical target for MES-4 repression. The above findings extend to the H3K27 methyltransferase MES-2/3/6, the C. elegans version of polycomb repressive complex 2. We propose that maternal MES-4 and PRC2 cooperate to protect germline survival by preventing synthesis of germline-toxic products encoded by genes on the X chromosome, including the key transcription factor LIN-15B.

Identification of healthspan-promoting genes in Caenorhabditis elegans based on a human GWAS study

To find drivers of healthy ageing, a genome-wide association study (GWAS) was performed in healthy and unhealthy older individuals. Healthy individuals were defined as free from cardiovascular disease, stroke, heart failure, major adverse cardiovascular event, diabetes, dementia, cancer, chronic obstructive pulmonary disease (COPD), asthma, rheumatism, Crohn’s disease, malabsorption or kidney disease. Six single nucleotide polymorphisms (SNPs) with unknown function associated with ten human genes were identified as candidate healthspan markers. Thirteen homologous or closely related genes were selected in the model organism C. elegans for evaluating healthspan after targeted RNAi-mediated knockdown using pathogen resistance, muscle integrity, chemotaxis index and the activity of known longevity and stress response pathways as healthspan reporters. In addition, lifespan was monitored in the RNAi-treated nematodes. RNAi knockdown of yap-1, wwp-1, paxt-1 and several acdh genes resulted in heterogeneous phenotypes regarding muscle integrity, pathogen resistance, chemotactic behaviour, and lifespan. Based on these observations, we hypothesize that their human homologues WWC2, CDKN2AIP and ACADS may play a role in health maintenance in the elderly.

Gut neuroendocrine signaling regulates synaptic assembly in C. elegans

Synaptic connections are essential to build a functional brain. How synapses are formed during development is a fundamental question in neuroscience. Recent studies provided evidence that the gut plays an important role in neuronal development through processing signals derived from gut microbes or nutrients. Defects in gut-brain communication can lead to various neurological disorders. Although the roles of the gut in communicating signals from its internal environment to the brain are well known, it remains unclear whether the gut plays a genetically encoded role in neuronal development. Using C. elegans as a model, we uncover that a Wnt-endocrine signaling pathway in the gut regulates synaptic development in the brain. A canonical Wnt signaling pathway promotes synapse formation through regulating the expression of the neuropeptides encoding gene nlp-40 in the gut, which functions through the neuronally expressed GPCR/AEX-2 receptor during development. Wnt-NLP-40-AEX-2 signaling likely acts to modulate neuronal activity. Our study reveals a genetic role of the gut in synaptic development and identifies a novel contribution of the gut-brain axis.

Tyraminergic corollary discharge filters reafferent perception in a chemosensory neuron

Interpreting sensory information requires its integration with the current behavior of the animal. However, how motor-related circuits influence sensory information processing is incompletely understood. Here, we report that current locomotor state directly modulates the activity of BAG CO₂ sensory neurons in Caenorhabditis elegans. By recording neuronal activity in animals freely navigating CO₂ landscapes, we found that during reverse crawling states, BAG activity is suppressed by tyraminergic corollary discharge signaling. We provide genetic evidence that tyramine released from the RIM reversal interneurons extrasynaptically activates the inhibitory chloride channel LGC-55 in BAG. Disrupting this pathway genetically leads to excessive behavioral responses to CO₂ stimuli. Moreover, we find that LGC-55 signaling cancels out perception of self-produced CO₂ and O₂ stimuli when animals reverse into their own gas plume in ethologically relevant aqueous environments. Our results show that sensorimotor integration involves corollary discharge signals directly modulating chemosensory neurons.

An extracellular matrix damage sensor signals through membrane-associated kinase DRL-1 to mediate cytoprotective responses in Caenorhabditis elegans

We and others previously identified circumferential bands of collagen named annular furrows as key components of a damage sensor in the cuticle of Caenorhabditis elegans that regulates cytoprotective genes. Mutation or loss of noncollagen secreted proteins OSM-7, OSM-8, and OSM-11 activate the same cytoprotective responses without obvious changes to the cuticle indicating that other extracellular proteins are involved. Here, we used RNAi screening to identify protein kinase DRL-1 as a key modulator of cytoprotective gene expression and stress resistance in furrow and extracellular OSM protein mutants. DRL-1 functions downstream from furrow disruption and is expressed in cells that induce cytoprotective genes. DRL-1 is not required for the expression of cytoprotective genes under basal or oxidative stress conditions consistent with specificity to extracellular signals. DRL-1 was previously shown to regulate longevity via a "Dietary Restriction-Like" state, but it functions downstream from furrow disruption by a distinct mechanism. The kinase domain of DRL-1 is related to mammalian MEKK3, and MEKK3 is recruited to a plasma membrane osmosensor complex by a scaffold protein. In C. elegans, DRL-1 contains an atypical hydrophobic C-terminus with predicted transmembrane domains and is constitutively expressed at or near the plasma membrane where it could function to receive extracellular damage signals for cells that mount cytoprotective responses.

AP-1 Recruits SMAP-1/SMAPs to the trans-Golgi Network to Promote Sorting in Polarized Epithelia

Coordinated AP-1 and clathrin coat assembly mediate secretory sorting on the trans-Golgi network (TGN) during conventional secretion. Here we found that SMAP-1/SMAPs deficiency caused the apical protein ERM-1 to accumulate on the basolateral side of the TGN. In contrast, the basolateral protein SLCF-1 appeared abnormally on the apical membrane. SMAP-1 colocalized with AP-1 on the TGN. The integrity of AP-1 is required for the subcellular presence of SMAP-1. Moreover, we found that the loss of SMAP-1 reduced clathrin-positive structures in the cytosol, suggesting that SMAP-1 has a regulatory role in clathrin assembly on the TGN. Functional experiments showed that overexpressing clathrin effectively alleviated exocytic defects due to the lack of SMAP-1, corroborating the role of SMAP-1 in promoting the assembly of clathrin on the TGN. Together, our results suggested that the AP-1 complex regulates the TGN localization of SMAP-1, promoting clathrin assembly to ensure polarized conventional secretion in C. elegans intestinal epithelia.

Sensory neuron transcriptomes reveal complex neuron-specific function and regulation of mec-2/Stomatin splicing

The function and identity of a cell is shaped by transcription factors controlling transcriptional networks, and further shaped by RNA binding proteins controlling post-transcriptional networks. To overcome limitations inherent to analysis of sparse single-cell post-transcriptional data, we leverage the invariant Caenorhabditis elegans cell lineage, isolating thousands of identical neuron types from thousands of isogenic individuals. The resulting deep transcriptomes facilitate splicing network analysis due to increased sequencing depth and uniformity. We focus on mechanosensory touch-neuron splicing regulated by MEC-8/RBPMS. We identify a small MEC-8-regulated network, where MEC-8 establishes touch-neuron isoforms differing from default isoforms found in other cells. MEC-8 establishes the canonical long mec-2/Stomatin isoform in touch neurons, but surprisingly the non-canonical short isoform predominates in other neurons, including olfactory neurons, and mec-2 is required for olfaction. Forced endogenous isoform-specific expression reveals that the short isoform functions in olfaction but not mechanosensation. The long isoform is functional in both processes. Remarkably, restoring the long isoform completely rescues mec-8 mutant mechanosensation, indicating a single MEC-8 touch-neuron target is phenotypically relevant. Within the long isoform we identify a cassette exon further diversifying mec-2 into long/extra-long isoforms. Neither is sufficient for mechanosensation. Both are simultaneously required, likely functioning as heteromers to mediate mechanosensation.

Phalloidin Staining of Actin Filaments for Visualization of Muscle Fibers in Caenorhabditis elegans

Advances in C. elegans research have allowed scientists to recapitulate different human disorders, from neurodegenerative diseases to muscle dysfunction, in these nematodes. Concomitantly, the interest in visualizing organs affected by these conditions has grown, leading to the establishment of different antibody- and dye-based staining protocols to verify tissue morphology. In particular, the quality of muscle tissue has been largely used in nematodes as a readout for fitness and healthspan. Phalloidin derivatives, which are commonly used to stain actin filaments in cells and tissues, have been implemented in the context of C. elegans research for visualization of muscle fibers. However, the majority of the phalloidin-based protocols depend on fixation steps using harmful compounds, preparation of specific buffers, and large amounts of worms. Herein, we implemented a safer and more flexible experimental procedure to stain actin filaments in C. elegans using phalloidin-based dyes. Lyophilization of the worms followed by their acetone permeabilization allows bypassing the fixation process while also providing the opportunity to suspend the experiment at different steps. Moreover, by using conventional buffers throughout our protocol, we avoid the additional preparation of solutions. Finally, our protocol requires a limited number of worms, making it suitable for slow-growing C. elegans strains. Overall, this protocol provides an efficient, fast, and safer method to stain actin filaments and visualize muscle fibers in C. elegans. Graphic abstract: Schematic overview of phalloidin staining in C. elegans for assessing muscle fiber morphology.

Gene expression profiling of epidermal cell types in C. elegans using Targeted DamID

The epidermis of Caenorhabditis elegans is an essential tissue for survival because it contributes to the formation of the cuticle barrier as well as facilitating developmental progression and animal growth. Most of the epidermis consists of the hyp7 hypodermal syncytium, the nuclei of which are largely generated by the seam cells, which exhibit stem cell-like behaviour during development. How seam cell progenitors differ transcriptionally from the differentiated hypodermis is poorly understood. Here, we introduce Targeted DamID (TaDa) in C. elegans as a method for identifying genes expressed within a tissue of interest without cell isolation. We show that TaDa signal enrichment profiles can be used to identify genes transcribed in the epidermis and use this method to resolve differences in gene expression between the seam cells and the hypodermis. Finally, we predict and functionally validate new transcription and chromatin factors acting in seam cell development. These findings provide insights into cell type-specific gene expression profiles likely associated with epidermal cell fate patterning.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Molecular Evolutionary Analysis of Nematode Zona Pellucida (ZP) Modules Reveals Disulfide-Bond Reshuffling and Standalone ZP-C Domains

Zona pellucida (ZP) modules mediate extracellular protein-protein interactions and contribute to important biological processes including syngamy and cellular morphogenesis. Although some biomedically relevant ZP modules are well studied, little is known about the protein family's broad-scale diversity and evolution. The increasing availability of sequenced genomes from "nonmodel" systems provides a valuable opportunity to address this issue and to use comparative approaches to gain new insights into ZP module biology. Here, through phylogenetic and structural exploration of ZP module diversity across the nematode phylum, I report evidence that speaks to two important aspects of ZP module biology. First, I show that ZP-C domains-which in some modules act as regulators of ZP-N domain-mediated polymerization activity, and which have never before been found in isolation-can indeed be found as standalone domains. These standalone ZP-C domain proteins originated in independent (paralogous) lineages prior to the diversification of extant nematodes, after which they evolved under strong stabilizing selection, suggesting the presence of ZP-N domain-independent functionality. Second, I provide a much-needed phylogenetic perspective on disulfide bond variability, uncovering evidence for both convergent evolution and disulfide-bond reshuffling. This result has implications for our evolutionary understanding and classification of ZP module structural diversity and highlights the usefulness of phylogenetics and diverse sampling for protein structural biology. All told, these findings set the stage for broad-scale (cross-phyla) evolutionary analysis of ZP modules and position Caenorhabditis elegans and other nematodes as important experimental systems for exploring the evolution of ZP modules and their constituent domains.

Comparative Transcriptomics across Nematode Life Cycles Reveal Gene Expression Conservation and Correlated Evolution in Adjacent Developmental Stages

Nematodes are highly abundant animals with diverse habitats and lifestyles. Some are free living whereas others parasitize animals or plants, and among the latter, infection abilities change across developmental stages to infect hosts and complete life cycles. To determine the relationship between transcriptome evolution and morphological divergences among nematodes, we compared 48 transcriptomes of different developmental stages across eight nematode species. The transcriptomes were clustered broadly into embryo, larva, and adult stages, with the developmental plastic stages were separated from common larval stages within the larval branch. This suggests that development was the major determining factor after lifestyle changes, such as parasitism, during transcriptome evolution. Such patterns were partly accounted for by tissue-specific genes—such as those in oocytes and the hypodermis—being expressed at different proportions. Although nematodes typically have 3–5 larval stages, the transcriptomes for these stages were found to be highly correlated within each species, suggesting high similarity among larval stages across species. For the Caenorhabditis elegans–Caenorhabditis briggsae and Strongyloides stercoralis–Strongyloides venezuelensis comparisons, we found that ∼50% of genes were expressed at multiple stages, whereas half of their orthologs were also expressed in multiple but different stages. Such frequent changes in expression have resulted in concerted transcriptome evolution across adjacent stages, thus generating species-specific transcriptomes over the course of nematode evolution. Our study provides a first insight into the evolution of nematode transcriptomes beyond embryonic development.

DAF-18/PTEN inhibits germline zygotic gene activation during primordial germ cell quiescence

Quiescence, an actively-maintained reversible state of cell cycle arrest, is not well understood. PTEN is one of the most frequently lost tumor suppressors in human cancers and regulates quiescence of stem cells and cancer cells. The sole PTEN ortholog in Caenorhabditis elegans is daf-18. In a C. elegans loss-of-function mutant for daf-18, primordial germ cells (PGCs) divide inappropriately in L1 larvae hatched into starvation conditions, in a TOR-dependent manner. Here, we further investigated the role of daf-18 in maintaining PGC quiescence in L1 starvation. We found that maternal or zygotic daf-18 is sufficient to maintain cell cycle quiescence, that daf-18 acts in the germ line and soma, and that daf-18 affects timing of PGC divisions in fed animals. Importantly, our results also implicate daf-18 in repression of germline zygotic gene activation, though not in germline fate specification. However, TOR is less important to germline zygotic gene expression, suggesting that in the absence of food, daf-18/PTEN prevents inappropriate germline zygotic gene activation and cell division by distinct mechanisms.

RTKN-1/Rhotekin shields endosome-associated F-actin from disassembly to ensure endocytic recycling

Cargo sorting and the subsequent membrane carrier formation require a properly organized endosomal actin network. To better understand the actin dynamics during endocytic recycling, we performed a genetic screen in C. elegans and identified RTKN-1/Rhotekin as a requisite to sustain endosome-associated actin integrity. Loss of RTKN-1 led to a prominent decrease in actin structures and basolateral recycling defects. Furthermore, we showed that the presence of RTKN-1 thwarts the actin disassembly competence of UNC-60A/cofilin. Consistently, in RTKN-1–deficient cells, UNC-60A knockdown replenished actin structures and alleviated the recycling defects. Notably, an intramolecular interaction within RTKN-1 could mediate the formation of oligomers. Overexpression of an RTKN-1 mutant form that lacks self-binding capacity failed to restore actin structures and recycling flow in rtkn-1 mutants. Finally, we demonstrated that SDPN-1/Syndapin acts to direct the recycling endosomal dwelling of RTKN-1 and promotes actin integrity there. Taken together, these findings consolidated the role of SDPN-1 in organizing the endosomal actin network architecture and introduced RTKN-1 as a novel regulatory protein involved in this process.

The Doublesex/Mab-3 domain transcription factor DMD-10 regulates ASH-dependent behavioral responses

The Doublesex/Mab-3 Domain transcription factor DMD-10 is expressed in several cell types in C. elegans, including in the nervous system. We sought to investigate whether DMD-10 is required for normal neuronal function using behavioral assays. We found that mutation of dmd-10 did not broadly affect behavior. dmd-10 mutants were normal in several behavioral assays including a body bends assay for locomotion, egg laying, chemotaxis and response to gentle touch to the body. dmd-10 mutants did have defects in nose-touch responsiveness, which requires the glutamate receptor GLR-1. However, using quantitative fluorescence microscopy to measure levels of a GLR-1::GFP fusion protein in the ventral nerve cord, we found no evidence supporting a difference in the number of GLR-1 synapses or in the amount of GLR-1 present in dmd-10 mutants. dmd-10 mutants did have decreased responsiveness to high osmolarity, which, along with nose-touch, is sensed by the polymodal sensory neuron ASH. Furthermore, mutation of dmd-10 impaired behavioral response to optogenetic activation of ASH, suggesting that dmd-10 promotes neuronal signaling in ASH downstream of sensory receptor activation. Together our results suggest that DMD-10 is important in regulating the frequency of multiple ASH-dependent behavioral responses.

VER/VEGF receptors regulate AMPA receptor surface levels and glutamatergic behavior

Several intracellular trafficking pathways contribute to the regulation of AMPA receptor (AMPAR) levels at synapses and the control of synaptic strength. While much has been learned about these intracellular trafficking pathways, a major challenge is to understand how extracellular factors, such as growth factors, neuropeptides and hormones, impinge on specific AMPAR trafficking pathways to alter synaptic function and behavior. Here, we identify the secreted ligand PVF-1 and its cognate VEGF receptor homologs, VER-1 and VER-4, as regulators of glutamate signaling in C. elegans. Loss of function mutations in ver-1, ver-4, or pvf-1, result in decreased cell surface levels of the AMPAR GLR-1 and defects in glutamatergic behavior. Rescue experiments indicate that PVF-1 is expressed and released from muscle, whereas the VERs function in GLR-1-expressing neurons to regulate surface levels of GLR-1 and glutamatergic behavior. Additionally, ver-4 is unable to rescue glutamatergic behavior in the absence of pvf-1, suggesting that VER function requires endogenous PVF-1. Inducible expression of a pvf-1 rescuing transgene suggests that PVF-1 can function in the mature nervous system to regulate GLR-1 signaling. Genetic double mutant analysis suggests that the VERs act together with the VPS-35/retromer recycling complex to promote cell surface levels of GLR-1. Our data support a genetic model whereby PVF-1/VER signaling acts with retromer to promote recycling and cell surface levels of GLR-1 to control behavior.

Sensation, behavior, and cognition all depend on the proper function of neuronal connections called synapses. Synapses that use the neurotransmitter glutamate to signal between nerve cells are the most abundant type in our brain. Presynaptic neurons release glutamate, which activates glutamate receptors on postsynaptic neurons. Dysfunction of glutamate synapses leads to several neurological disorders, and changing their strength–in part by altering glutamate receptors numbers on the surface of the postsynaptic cell—provides the cellular basis of learning and memory. Much remains to be learned about how factors released from other cell types affects synaptic communication. We took advantage of light-activated molecular switches engineered into specific sensory neurons of C. elegans worms to trigger a behavioral reflex that depends on glutamate synapses. Using this behavior, we identified proteins called VER-1 and VER-4 as important for glutamate synapse function. We found that worms missing these VER proteins or their activator PVF-1 have reduced levels of glutamate receptors at the postsynaptic surface and defects in glutamate-dependent behaviors. Our results suggest that inter-tissue cross-talk between muscle PVF-1 and neuronal VERs is important for controlling the number of glutamate receptors at the cell surface, robust neuronal communication and behavioral responses.

Conserved nicotine-activated neuroprotective pathways involve mitochondrial stress

Tobacco smoking is a risk factor for several human diseases. Conversely, smoking also reduces the prevalence of Parkinson's disease, whose hallmark is degeneration of substantia nigra dopaminergic neurons (DNs). We use C. elegans as a model to investigate whether tobacco-derived nicotine activates nicotinic acetylcholine receptors (nAChRs) to selectively protect DNs. Using this model, we demonstrate conserved functions of DN-expressed nAChRs. We find that DOP-2, a D3-receptor homolog; MCU-1, a mitochondrial calcium uniporter; PINK-1 (PTEN-induced kinase 1); and PDR-1 (Parkin) are required for nicotine-mediated protection of DNs. Together, our results support involvement of a calcium-modulated, mitochondrial stress-activated PINK1/Parkin-dependent pathway in nicotine-induced neuroprotection. This suggests that nicotine-selective protection of substantia nigra DNs is due to the confluence of two factors: first, their unique vulnerability to mitochondrial stress, which is mitigated by increased mitochondrial quality control due to PINK1 activation, and second, their specific expression of D3-receptors.