Spatial Transcriptomics of C. elegans Males and Hermaphrodites Identifies Sex-Specific Differences in Gene Expression Patterns

Combined matrix-assisted laser desorption/ionisation-mass spectrometry imaging with liquid chromatography-tandem mass spectrometry for observing spatial distribution of lipids in whole Caenorhabditis elegans

RATIONALE

Matrix-assisted laser desorption/ionisation-mass spectrometry imaging (MALDI-MSI) is a powerful label-free technique for biomolecule detection (e.g., lipids), within tissue sections across various biological species. However, despite its utility in many applications, the nematode Caenorhabditis elegans is not routinely used in combination with MALDI-MSI. The lack of studies exploring spatial distribution of biomolecules in nematodes is likely due to challenges with sample preparation.

METHODS

This study developed a sample preparation method for whole intact nematodes, evaluated using cryosectioning of nematodes embedded in a 10% gelatine solution to obtain longitudinal cross sections. The slices were then subjected to MALDI-MSI, using a RapifleX Tissuetyper in positive and negative polarities. Samples were also prepared for liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis using an Exploris 480 coupled to a HPLC Vanquish system to confirm the MALDI-MSI results.

RESULTS

An optimised embedding method was developed for longitudinal cross-sectioning of individual worms. To obtain longitudinal cross sections, nematodes were frozen at -80°C so that all worms were rod shaped. Then, the samples were defrosted and transferred to a 10% gelatine matrix in a cryomold; the worms aligned, and the whole cryomold submerged in liquid nitrogen. Using MALDI-MSI, we were able to observe the distribution of lipids within C. elegans, with clear differences in their spatial distribution at a resolution of 5 μm. To confirm the lipids from MALDI-MSI, age-matched nematodes were subjected to LC-MS/MS. Here, 520 lipids were identified using LC-MS/MS, indicating overlap with MALDI-MSI data.

CONCLUSIONS

This optimised sample preparation technique enabled (un)targeted analysis of spatially distributed lipids within individual nematodes. The possibility to detect other biomolecules using this method thus laid the basis for prospective preclinical and toxicological studies on C. elegans.

Male-specific behavioral and transcriptomic changes in aging C. elegans neurons

Aging is a complex biological process with sexually dimorphic aspects. Although cognitive aging of Caenorhabditis elegans hermaphrodites has been studied, less is known about cognitive decline in males. We found that cognitive aging has both sex-shared and sex-dimorphic characteristics, and we identified neuron-specific age-associated sex-differential targets. In addition to sex-shared neuronal aging genes, males differentially downregulate mitochondrial metabolic genes and upregulate GPCR genes with age, while the X chromosome exhibits increased gene expression in hermaphrodites and altered dosage compensation complex expression with age, indicating possible X chromosome dysregulation that contributes to sexual dimorphism in cognitive aging. Finally, the sex-differentially expressed gene hrg-7, an aspartic-type endopeptidase, regulates male cognitive aging but does not affect hermaphrodites' behaviors. These results suggest that males and hermaphrodites exhibit different age-related neuronal changes. This study will strengthen our understanding of sex-specific vulnerability and resilience and identify pathways to target with treatments that could benefit both sexes.

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.

Evolutionarily related host and microbial pathways regulate fat desaturation in C. elegans

Fatty acid desaturation is central to metazoan lipid metabolism and provides building blocks of membrane lipids and precursors of diverse signaling molecules. Nutritional conditions and associated microbiota regulate desaturase expression, but the underlying mechanisms have remained unclear. Here, we show that endogenous and microbiota-dependent small molecule signals promote lipid desaturation via the nuclear receptor NHR-49/PPARα in C. elegans. Untargeted metabolomics of a β-oxidation mutant, acdh-11, in which expression of the stearoyl-CoA desaturase FAT-7/SCD1 is constitutively increased, revealed accumulation of a β-cyclopropyl fatty acid, becyp#1, that potently activates fat-7 expression via NHR-49. Biosynthesis of becyp#1 is strictly dependent on expression of cyclopropane synthase by associated bacteria, e.g., E. coli. Screening for structurally related endogenous metabolites revealed a β-methyl fatty acid, bemeth#1, which mimics the activity of microbiota-dependent becyp#1 but is derived from a methyltransferase, fcmt-1, that is conserved across Nematoda and likely originates from bacterial cyclopropane synthase via ancient horizontal gene transfer. Activation of fat-7 expression by these structurally similar metabolites is controlled by distinct mechanisms, as microbiota-dependent becyp#1 is metabolized by a dedicated β-oxidation pathway, while the endogenous bemeth#1 is metabolized via α-oxidation. Collectively, we demonstrate that evolutionarily related biosynthetic pathways in metazoan host and associated microbiota converge on NHR-49/PPARα to regulate fat desaturation.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

Spatial Transcriptomics: Emerging Technologies in Tissue Gene Expression Profiling

In this Perspective, we discuss the current status and advances in spatial transcriptomics technologies, which allow high-resolution mapping of gene expression in intact cell and tissue samples. Spatial transcriptomics enables the creation of high-resolution maps of gene expression patterns within their native spatial context, adding an extra layer of information to the bulk sequencing data. Spatial transcriptomics has expanded significantly in recent years and is making a notable impact on a range of fields, including tissue architecture, developmental biology, cancer, and neurodegenerative and infectious diseases. The latest advancements in spatial transcriptomics have resulted in the development of highly multiplexed methods, transcriptomic-wide analysis, and single-cell resolution utilizing diverse technological approaches. In this Perspective, we provide a detailed analysis of the molecular foundations behind the main spatial transcriptomics technologies, including methods based on microdissection, in situ sequencing, single-molecule FISH, spatial capturing, selection of regions of interest, and single-cell or nuclei dissociation. We contextualize the detection and capturing efficiency, strengths, limitations, tissue compatibility, and applications of these techniques as well as provide information on data analysis. In addition, this Perspective discusses future directions and potential applications of spatial transcriptomics, highlighting the importance of the continued development to promote widespread adoption of these techniques within the research community.

Miniature spatial transcriptomics for studying parasite-endosymbiont relationships at the micro scale

Several important human infectious diseases are caused by microscale-sized parasitic nematodes like filarial worms. Filarial worms have their own spatial tissue organization; to uncover this tissue structure, we need methods that can spatially resolve these miniature specimens. Most filarial worms evolved a mutualistic association with endosymbiotic bacteria Wolbachia. However, the mechanisms underlying the dependency of filarial worms on the fitness of these bacteria remain unknown. As Wolbachia is essential for the development, reproduction, and survival of filarial worms, we spatially explored how Wolbachia interacts with the worm's reproductive system by performing a spatial characterization using Spatial Transcriptomics (ST) across a posterior region containing reproductive tissue and developing embryos of adult female Brugia malayi worms. We provide a proof-of-concept for miniature-ST to explore spatial gene expression patterns in small sample types, demonstrating the method's ability to uncover nuanced tissue region expression patterns, observe the spatial localization of key B. malayi - Wolbachia pathway genes, and co-localize the B. malayi spatial transcriptome in Wolbachia tissue regions, also under antibiotic treatment. We envision our approach will open up new avenues for the study of infectious diseases caused by micro-scale parasitic worms.

Male gonad-enriched microRNAs function to control sperm production in C. elegans

Germ cell development and gamete production in animals require small RNA pathways. While studies indicate that microRNAs (miRNAs) are necessary for normal sperm production and function, the specific roles for individual miRNAs are largely unknown. Here, we use small RNA sequencing of dissected gonads and functional analysis of new loss of function alleles to identify functions for miRNAs in the control of fecundity and sperm production in Caenorhabditis elegans males and hermaphrodites. We describe a set of 29 male gonad-enriched miRNAs and identify a set of 3 individual miRNAs (mir-58.1, mir-83, and mir-235) and a miRNA cluster (mir-4807-4810.1) that are required for optimal sperm production at 20°C and 5 additional miRNAs (mir-49, mir-57, mir-261, and mir-357/358) that are required for sperm production at 25°C. We observed defects in meiotic progression in mir-58.1, mir-83, mir-235, and mir-4807-4810.1 mutants that may contribute to the reduced number of sperm. Further, analysis of multiple mutants of these miRNAs suggested complex genetic interactions between these miRNAs for sperm production. This study provides insights on the regulatory roles of miRNAs that promote optimal sperm production and fecundity in males and hermaphrodites.

Evolutionarily related host and microbial pathways regulate fat desaturation

Fatty acid desaturation is central to metazoan lipid metabolism and provides building blocks of membrane lipids and precursors of diverse signaling molecules. Nutritional conditions and associated microbiota regulate desaturase expression1-4, but the underlying mechanisms have remained unclear. Here, we show that endogenous and microbiota-dependent small molecule signals promote lipid desaturation via the nuclear receptor NHR-49/PPARα in C. elegans. Untargeted metabolomics of a β-oxidation mutant, acdh-11, in which expression of the stearoyl-CoA desaturase FAT-7/SCD1 is constitutively increased, revealed accumulation of a β-cyclopropyl fatty acid, becyp#1, that potently activates fat-7 expression via NHR-49. Biosynthesis of becyp#1 is strictly dependent on expression of cyclopropane synthase by associated bacteria, e.g., E. coli. Screening for structurally related endogenous metabolites revealed a β-methyl fatty acid, bemeth#1, whose activity mimics that of microbiota-dependent becyp#1, but is derived from a methyltransferase, fcmt-1, that is conserved across Nematoda and likely originates from bacterial cyclopropane synthase via ancient horizontal gene transfer. Activation of fat-7 expression by these structurally similar metabolites is controlled by distinct mechanisms, as microbiota-dependent becyp#1 is metabolized by a dedicated β-oxidation pathway, while the endogenous bemeth#1 is metabolized via α-oxidation. Collectively, we demonstrate that evolutionarily related biosynthetic pathways in metazoan host and associated microbiota converge on NHR-49/PPARα to regulate fat desaturation.

Caenorhabditis elegans in microgravity: An omics perspective

The application of omics to study Caenorhabditis elegans (C. elegans) in the context of spaceflight is increasing, illuminating the wide-ranging biological impacts of spaceflight on physiology. In this review, we highlight the application of omics, including transcriptomics, genomics, proteomics, multi-omics, and integrated omics in the study of spaceflown C. elegans, and discuss the impact, use, and future direction of this branch of research. We highlight the variety of molecular alterations that occur in response to spaceflight, most notably changes in metabolic and neuromuscular gene regulation. These transcriptional features are reproducible and evident across many spaceflown species (e.g., mice and astronauts), supporting the use of C. elegans as a model organism to study spaceflight physiology with translational capital. Integrating tissue-specific, spatial, and multi-omics approaches, which quantitatively link molecular responses to phenotypic adaptations, will facilitate the identification of candidate regulatory molecules for therapeutic intervention and thus represents the next frontiers in C. elegans space omics research.

A Caenorhabditis elegans Male Pheromone Feminizes Germline Gene Expression in Hermaphrodites and Imposes Life-History Costs

Sex pheromones not only improve the reproductive success of the recipients, but also impose costs, such as a reduced life span. The underlying mechanisms largely remain to be elucidated. Here, we show that even a brief exposure to physiological amounts of the dominant Caenorhabditis elegans male pheromone, ascr#10, alters the expression of thousands of genes in hermaphrodites. The most dramatic effect on the transcriptome is the upregulation of genes expressed during oogenesis and the downregulation of genes associated with male gametogenesis. This result reveals a way in which social signals help to resolve the inherent conflict between spermatogenesis and oogenesis in a simultaneous hermaphrodite, presumably to optimally align reproductive function with the presence of potential mating partners. We also found that exposure to ascr#10 increased the risk of persistent intestinal infections in hermaphrodites due to pathological pharyngeal hypertrophy. Thus, our study reveals ways in which the male pheromone can not only have beneficial effects on the recipients' reproduction, but also cause harmful consequences that reduce life span.

A C. elegans male pheromone feminizes germline gene expression in hermaphrodites and imposes life-history costs

Sex pheromones improve reproductive success, but also impose costs. Here we show that even brief exposure to physiological amounts of the dominant C. elegans male pheromone, ascr#10, alters the expression of thousands of genes in hermaphrodites. The most dramatic effect on the transcriptome was the upregulation of genes expressed during oogenesis and downregulation of genes associated with male gametogenesis. Among the detrimental effects of ascr#10 on hermaphrodites is the increased risk of persistent infections caused by pathological pharyngeal hypertrophy. Our results reveal a way in which social signals help to resolve the inherent conflict between spermatogenesis and oogenesis in a simultaneous hermaphrodite, presumably to optimally align reproductive function to the presence of potential mating partners. They also show that the beneficial effects of the pheromone are accompanied by harmful consequences that reduce lifespan.

Analysis of the C. elegans Germline Stem Cell Pool

The Caenorhabditis elegans germline is an excellent model for studying the genetic and molecular regulation of stem cell self-renewal and progression of cells from a stem cell state to a differentiated state. The germline tissue is organized in an assembly line with the germline stem cell (GSC) pool at one end and differentiated gametes at the other. A simple mesenchymal niche caps the GSC pool and maintains GSCs in an undifferentiated state by signaling through the conserved Notch pathway. Notch signaling activates transcription of the key GSC regulators lst-1 and sygl-1 proteins in a gradient through the GSC pool. LST-1 and SYGL-1 proteins work with PUF RNA regulators in a self-renewal hub to maintain the GSC pool. In this chapter, we present methods for characterizing the C. elegans GSC pool and early stages of germ cell differentiation. The methods include examination of germlines in living and fixed worms, cell cycle analysis, and analysis of markers. We also discuss assays to separate mutant phenotypes that affect the stem cell vs. differentiation decision from those that affect germ cell processes more generally.

A widely distributed family of eukaryotic and bacterial deubiquitinases related to herpesviral large tegument proteins

Distinct families of eukaryotic deubiquitinases (DUBs) are regulators of ubiquitin signaling. Here, we report on the presence of an additional DUB class broadly distributed in eukaryotes and several bacteria. The only described members of this family are the large tegument proteins of herpesviruses, which are attached to the outside of the viral capsid. By using a bioinformatics screen, we have identified distant homologs of this VTD (Viral tegument-like DUB) family in vertebrate transposons, fungi, insects, nematodes, cnidaria, protists and bacteria. While some VTD activities resemble viral tegument DUBs in that they favor K48-linked ubiquitin chains, other members are highly specific for K6- or K63-linked ubiquitin chains. The crystal structures of K48- and K6-specific members reveal considerable differences in ubiquitin recognition. The VTD family likely evolved from non-DUB proteases and spread through transposons, many of which became 'domesticated', giving rise to the Drosophila male sterile (3)76Ca gene and several nematode genes with male-specific expression.

Computational modeling and analysis of the morphogenetic domain signaling networks regulating C. elegans embryogenesis

Caenorhabditis elegans, often referred to as the ‘roundworm’, provides a powerful model for studying cell autonomous and cell–cell interactions through the direct observation of embryonic development in vivo. By leveraging the precisely mapped cell lineage at single cell resolution, we are able to study at a systems level how early embryonic cells communicate across morphogenetic domains for the coordinated processes of gene expressions and collective cellular behaviors that regulate tissue morphogenesis. In this study, we developed a computational framework for the exploration of the morphogenetic domain cell signaling networks that may regulate C. elegans gastrulation and embryonic organogenesis. We demonstrated its utility by producing the following results, i) established a virtual reference model of developing C. elegans embryos through the spatiotemporal alignment of individual embryo cell nuclear imaging samples; ii) integrated the single cell spatiotemporal gene expression profile with the established virtual embryo model by data pooling; iii) trained a Machine Learning model (Random Forest Regression), which predicts accurately the spatial positions of the cells given their gene expression profiles for a given developmental time (e.g. total cell number of the embryo); iv) enabled virtual 4-dimensional tomographic graphical modeling of single cell data; v) inferred the biology signaling pathways that act in each of morphogenetic domains by meta-data analysis. It is intriguing that the morphogenetic domain cell signaling network seems to involve some crosstalk of multiple biology signaling pathways during the formation of tissue boundary pattern. Lastly, we developed the Software tool ‘Embryo aligner version 1.0’ and provided it as an Open Source program to the research community for virtual embryo modeling, and phenotype perturbation analyses (https://github.com/csniuben/embryo_aligner/wiki and https://bioinfo89.github.io/C.elegansEmbryonicOrganogenesisweb/).

Isolation, profiling, and tracking of extracellular vesicle cargo in Caenorhabditis elegans

Extracellular vesicles (EVs) may mediate intercellular communication by carrying protein and RNA cargo. The composition, biology, and roles of EVs in physiology and pathology have been primarily studied in the context of biofluids and in cultured mammalian cells. The experimental tractability of C. elegans makes for a powerful in vivo animal system to identify and study EV cargo from its cellular source. We developed an innovative method to label, track, and profile EVs using genetically encoded, fluorescent-tagged EV cargo and conducted a large-scale isolation and proteomic profiling. Nucleic acid binding proteins (∼200) are overrepresented in our dataset. By integrating our EV proteomic dataset with single-cell transcriptomic data, we identified and validated ciliary EV cargo: CD9-like tetraspanin (TSP-6), ectonucleotide pyrophosphatase/phosphodiesterase (ENPP-1), minichromosome maintenance protein (MCM-3), and double-stranded RNA transporter SID-2. C. elegans EVs also harbor RNA, suggesting that EVs may play a role in extracellular RNA-based communication.

Spatial transcriptomics reveals antiparasitic targets associated with essential behaviors in the human parasite Brugia malayi

Lymphatic filariasis (LF) is a chronic debilitating neglected tropical disease (NTD) caused by mosquito-transmitted nematodes that afflicts over 60 million people. Control of LF relies on routine mass drug administration with antiparasitics that clear circulating larval parasites but are ineffective against adults. The development of effective adulticides is hampered by a poor understanding of the processes and tissues driving parasite survival in the host. The adult filariae head region contains essential tissues that control parasite feeding, sensory, secretory, and reproductive behaviors, which express promising molecular substrates for the development of antifilarial drugs, vaccines, and diagnostics. We have adapted spatial transcriptomic approaches to map gene expression patterns across these prioritized but historically intractable head tissues. Spatial and tissue-resolved data reveal distinct biases in the origins of known drug targets and secreted antigens. These data were used to identify potential new drug and vaccine targets, including putative hidden antigens expressed in the alimentary canal, and to spatially associate receptor subunits belonging to druggable families. Spatial transcriptomic approaches provide a powerful resource to aid gene function inference and seed antiparasitic discovery pipelines across helminths of relevance to human and animal health.

High Polymorphism in the Dmrt2a Gene Is Incompletely Sex-Linked in Spotted Scat, Scatophagus argus

Unlike mammals and birds, many fishes have young sex chromosomes, providing excellent models to study sex chromosome differentiation at early stages. Previous studies showed that spotted scat possesses an XX-XY sex determination system. The X has a complete Dmrt3 copy (termed normal) and a truncated copy of Dmrt1 (called Dmrt1b), while the Y has the opposite (normal Dmrt1, which is male-specific, and a truncated Dmrt3 called Dmrt3△-Y). Dmrt1 is the candidate sex determination gene, while the differentiation of other sex-linked genes remains unknown. The spotted scat has proven to be a good model to study the evolution of sex chromosomes in vertebrates. Herein, we sequenced a neighbor gene of this family, Dmrt2, positioned farther from Dmrt1 and closer to Dmrt3 in the spotted scat, and analyzed its sequence variation and expression profiles. The physical locations of the three genes span across an estimated size of >40 kb. The open reading frames of Dmrt2a and its paralog Dmrt2b are 1578 bp and 1311 bp, encoding peptides of 525 and 436 amino acid residues, respectively. Dmrt2a is positioned close to Dmrt3 but farther from Dmrt1 on the same chromosome, while Dmrt2b is not. Sequence analysis revealed several mutations; insertions, and deletions (indels) on Dmrt2a non-coding regions and single-nucleotide polymorphisms (SNPs) on the Dmrt2a transcript. These indels and SNPs are sex-linked and showed high male heterogeneity but do not affect gene translation. The markers designed to span the mutation sites tested on four different populations showed varied concordance with the genetic sexes. Dmrt2a is transcribed solely in the gonads and gills, while Dmrt2b exists in the gonads, hypothalamus, gills, heart, and spleen. The Dmrt2a and Dmrt2b transcripts are profoundly expressed in the male gonads. Analyses of the transcriptome data from five other fish species (Hainan medaka (Oryzias curvinotus), silver sillago (Sillago sihama), Nile tilapia (Oreochromis niloticus), Hong Kong catfish (Clarias fuscus), and spot-fin porcupine fish (Diodon hystrix)) revealed testes-biased expression of Dmrt1 in all, similar to spotted scat. Additionally, the expression of Dmrt2a is higher in the testes than the ovaries in spotted scat and Hainan medaka. The Dmrt2a transcript was not altered in the coding regions as found in Dmrt1 and Dmrt3 in spotted scat. This could be due to the functional importance of Dmrt2a in development. Another possibility is that because Dmrt2a is positioned farther from Dmrt1 and the chromosome is still young, meaning it is only a matter of time before it differentiates. This study undeniably will aid in understanding the functional divergence of the sex-linked genes in fish.

Balancing Selection of the Intracellular Pathogen Response in Natural Caenorhabditis elegans Populations

Genetic variation in host populations may lead to differential viral susceptibilities. Here, we investigate the role of natural genetic variation in the Intracellular Pathogen Response (IPR), an important antiviral pathway in the model organism Caenorhabditis elegans against Orsay virus (OrV). The IPR involves transcriptional activity of 80 genes including the pals-genes. We examine the genetic variation in the pals-family for traces of selection and explore the molecular and phenotypic effects of having distinct pals-gene alleles. Genetic analysis of 330 global C. elegans strains reveals that genetic diversity within the IPR-related pals-genes can be categorized in a few haplotypes worldwide. Importantly, two key IPR regulators, pals-22 and pals-25, are in a genomic region carrying signatures of balancing selection, suggesting that different evolutionary strategies exist in IPR regulation. We infected eleven C. elegans strains that represent three distinct pals-22 pals-25 haplotypes with Orsay virus to determine their susceptibility. For two of these strains, N2 and CB4856, the transcriptional response to infection was also measured. The results indicate that pals-22 pals-25 haplotype shapes the defense against OrV and host genetic variation can result in constitutive activation of IPR genes. Our work presents evidence for balancing genetic selection of immunity genes in C. elegans and provides a novel perspective on the functional diversity that can develop within a main antiviral response in natural host populations.

piRNAs initiate transcriptional silencing of spermatogenic genes during C. elegans germline development

Eukaryotic genomes harbor invading transposable elements that are silenced by PIWI-interacting RNAs (piRNAs) to maintain genome integrity in animal germ cells. However, whether piRNAs also regulate endogenous gene expression programs remains unclear. Here, we show that C. elegans piRNAs trigger the transcriptional silencing of hundreds of spermatogenic genes during spermatogenesis, promoting sperm differentiation and function. This silencing signal requires piRNA-dependent small RNA biogenesis and loading into downstream nuclear effectors, which correlates with the dynamic reorganization of two distinct perinuclear biomolecular condensates present in germ cells. In addition, the silencing capacity of piRNAs is temporally counteracted by the Argonaute CSR-1, which targets and licenses spermatogenic gene transcription. The spatial and temporal overlap between these opposing small RNA pathways contributes to setting up the timing of the spermatogenic differentiation program. Thus, our work identifies a prominent role for piRNAs as direct regulators of endogenous transcriptional programs during germline development and gamete differentiation.

Neuronal HSF-1 coordinates the propagation of fat desaturation across tissues to enable adaptation to high temperatures in C. elegans

To survive elevated temperatures, ectotherms adjust the fluidity of membranes by fine-tuning lipid desaturation levels in a process previously described to be cell autonomous. We have discovered that, in Caenorhabditis elegans, neuronal heat shock factor 1 (HSF-1), the conserved master regulator of the heat shock response (HSR), causes extensive fat remodeling in peripheral tissues. These changes include a decrease in fat desaturase and acid lipase expression in the intestine and a global shift in the saturation levels of plasma membrane's phospholipids. The observed remodeling of plasma membrane is in line with ectothermic adaptive responses and gives worms a cumulative advantage to warm temperatures. We have determined that at least 6 TAX-2/TAX-4 cyclic guanosine monophosphate (cGMP) gated channel expressing sensory neurons, and transforming growth factor ß (TGF-β)/bone morphogenetic protein (BMP) are required for signaling across tissues to modulate fat desaturation. We also find neuronal hsf-1 is not only sufficient but also partially necessary to control the fat remodeling response and for survival at warm temperatures. This is the first study to show that a thermostat-based mechanism can cell nonautonomously coordinate membrane saturation and composition across tissues in a multicellular animal.

Bisection of the X chromosome disrupts the initiation of chromosome silencing during meiosis in Caenorhabditis elegans

During meiosis, gene expression is silenced in aberrantly unsynapsed chromatin and in heterogametic sex chromosomes. Initiation of sex chromosome silencing is disrupted in meiocytes with sex chromosome-autosome translocations. To determine whether this is due to aberrant synapsis or loss of continuity of sex chromosomes, we engineered Caenorhabditis elegans nematodes with non-translocated, bisected X chromosomes. In early meiocytes of mutant males and hermaphrodites, X segments are enriched with euchromatin assembly markers and active RNA polymerase II staining, indicating active transcription. Analysis of RNA-seq data showed that genes from the X chromosome are upregulated in gonads of mutant worms. Contrary to previous models, which predicted that any unsynapsed chromatin is silenced during meiosis, our data indicate that unsynapsed X segments are transcribed. Therefore, our results suggest that sex chromosome chromatin has a unique character that facilitates its meiotic expression when its continuity is lost, regardless of whether or not it is synapsed.

Spatial tissue profiling by imaging-free molecular tomography

Several techniques are currently being developed for spatially resolved omics profiling, but each new method requires the setup of specific detection strategies or specialized instrumentation. Here we describe an imaging-free framework to localize high-throughput readouts within a tissue by cutting the sample into thin strips in a way that allows subsequent image reconstruction. We implemented this framework to transform a low-input RNA sequencing protocol into an imaging-free spatial transcriptomics technique (called STRP-seq) and validated it by profiling the spatial transcriptome of the mouse brain. We applied the technique to the brain of the Australian bearded dragon, Pogona vitticeps. Our results reveal the molecular anatomy of the telencephalon of this lizard, providing evidence for a marked regionalization of the reptilian pallium and subpallium. We expect that STRP-seq can be used to derive spatially resolved data from a range of other omics techniques.

Spatial transcriptomics of the nematode Caenorhabditis elegans using RNA tomography

RNA tomography or tomo-seq combines mRNA sequencing and cryo-sectioning to spatially resolve gene expression. We have adapted this method for the nematode Caenorhabditis elegans to generate anteroposterior gene expression maps at near-cellular resolution. Here, we provide a detailed overview of the method and present two approaches: one that includes RNA isolation for maximum sensitivity and one that is suitable for partial automatization and is therefore less time-consuming. For complete details on the use and execution of this protocol, please refer to Ebbing et al. (2018).

Sexual Dimorphism through the Lens of Genome Manipulation, Forward Genetics, and Spatiotemporal Sequencing

Sexual reproduction often leads to selection that favors the evolution of sex-limited traits or sex-specific variation for shared traits. These sexual dimorphisms manifest due to sex-specific genetic architectures and sex-biased gene expression across development, yet the molecular mechanisms underlying these patterns are largely unknown. The first step is to understand how sexual dimorphisms arise across the genotype-phenotype-fitness map. The emergence of "4D genome technologies" allows for efficient, high-throughput, and cost-effective manipulation and observations of this process. Studies of sexual dimorphism will benefit from combining these technological advances (e.g., precision genome editing, inducible transgenic systems, and single-cell RNA sequencing) with clever experiments inspired by classic designs (e.g., bulked segregant analysis, experimental evolution, and pedigree tracing). This perspective poses a synthetic view of how manipulative approaches coupled with cutting-edge observational methods and evolutionary theory are poised to uncover the molecular genetic basis of sexual dimorphism with unprecedented resolution. We outline hypothesis-driven experimental paradigms for identifying genetic mechanisms of sexual dimorphism among tissues, across development, and over evolutionary time.

Spatial Transcriptomics of Nematodes Identifies Sperm Cells as a Source of Genomic Novelty and Rapid Evolution

Divergence of gene function and expression during development can give rise to phenotypic differences at the level of cells, tissues, organs, and ultimately whole organisms. To gain insights into the evolution of gene expression and novel genes at spatial resolution, we compared the spatially resolved transcriptomes of two distantly related nematodes, Caenorhabditis elegans and Pristionchus pacificus, that diverged 60-90 Ma. The spatial transcriptomes of adult worms show little evidence for strong conservation at the level of single genes. Instead, regional expression is largely driven by recent duplication and emergence of novel genes. Estimation of gene ages across anatomical structures revealed an enrichment of novel genes in sperm-related regions. This provides first evidence in nematodes for the "out of testis" hypothesis that has been previously postulated based on studies in Drosophila and mammals. "Out of testis" genes represent a mix of products of pervasive transcription as well as fast evolving members of ancient gene families. Strikingly, numerous novel genes have known functions during meiosis in Caenorhabditis elegans indicating that even universal processes such as meiosis may be targets of rapid evolution. Our study highlights the importance of novel genes in generating phenotypic diversity and explicitly characterizes gene origination in sperm-related regions. Furthermore, it proposes new functions for previously uncharacterized genes and establishes the spatial transcriptome of Pristionchus pacificus as a catalog for future studies on the evolution of gene expression and function.

The conserved molting/circadian rhythm regulator NHR-23/NR1F1 serves as an essential co-regulator of C. elegans spermatogenesis

In sexually reproducing metazoans, spermatogenesis is the process by which uncommitted germ cells give rise to haploid sperm. Work in model systems has revealed mechanisms controlling commitment to the sperm fate, but how this fate is subsequently executed remains less clear. While studying the well-established role of the conserved nuclear hormone receptor transcription factor, NHR-23/NR1F1, in regulating C. elegans molting, we discovered that NHR-23/NR1F1 is also constitutively expressed in developing primary spermatocytes and is a critical regulator of spermatogenesis. In this novel role, NHR-23/NR1F1 functions downstream of the canonical sex-determination pathway. Degron-mediated depletion of NHR-23/NR1F1 within hermaphrodite or male germlines causes sterility due to an absence of functional sperm, as depleted animals produce arrested primary spermatocytes rather than haploid sperm. These spermatocytes arrest in prometaphase I and fail to either progress to anaphase or attempt spermatid-residual body partitioning. They make sperm-specific membranous organelles but fail to assemble their major sperm protein into fibrous bodies. NHR-23/NR1F1 appears to function independently of the known SPE-44 gene regulatory network, revealing the existence of an NHR-23/NR1F1-mediated module that regulates the spermatogenesis program.

The molecular underpinnings of fertility: Genetic approaches in Caenorhabditis elegans

The study of mutations that impact fertility has a catch-22. Fertility mutants are often lost since they cannot simply be propagated and maintained. This has hindered progress in understanding the genetics of fertility. In mice, several molecules are found to be required for the interactions between the sperm and egg, with JUNO and IZUMO1 being the only known receptor pair on the egg and sperm surface, respectively. In Caenorhabditis elegans, a total of 12 proteins on the sperm or oocyte have been identified to mediate gamete interactions. Majority of these genes were identified through mutants isolated from genetic screens. In this review, we summarize the several key screening strategies that led to the identification of fertility mutants in C. elegans and provide a perspective about future research using genetic approaches. Recently, advancements in new technologies such as high-throughput sequencing and Crispr-based genome editing tools have accelerated the molecular, cell biological, and mechanistic analysis of fertility genes. We review how these valuable tools advance our understanding of the molecular underpinnings of fertilization. We draw parallels of the molecular mechanisms of fertilization between worms and mammals and argue that our work in C. elegans complements fertility research in humans and other species.

Roles for the RNA polymerase III regulator MAFR-1 in regulating sperm quality in Caenorhabditis elegans

The negative regulator of RNA polymerase (pol) III mafr-1 has been shown to affect RNA pol III transcript abundance, lipid biosynthesis and storage, progeny output, and lifespan. We deleted mafr-1 from the Caenorhabditis elegans genome and found that animals lacking mafr-1 replicated many phenotypes from previous RNAi-based studies and discovered a new sperm-specific role. Utilizing a yeast two-hybrid assay, we discovered several novel interactors of MAFR-1 that are expressed in a sperm- and germline-enriched manner. In support of a role for MAFR-1 in the male germline, we found mafr-1 null males have smaller spermatids that are less capable in competition for fertilization; a phenotype that was dependent on RNA pol III activity. Restoration of MAFR-1 expression specifically in the germline rescued the spermatid-related phenotypes, suggesting a cell autonomous role for MAFR-1 in nematode male fertility. Based on the high degree of conservation of Maf1 activity across species, our study may inform similar roles for Maf1 and RNA pol III in mammalian male fertility.

Speciation and the developmental alarm clock

New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.

Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides that are transcribed from non-coding loci yet undergo biosynthesis similar to coding mRNAs. The disproportional number of lincRNAs expressed in testes suggests that lincRNAs are important during gametogenesis, but experimental evidence has implicated very few lincRNAs in this process. We took advantage of the relatively limited number of lincRNAs in the genome of the nematode Caenorhabditis elegans to systematically analyse the functions of lincRNAs during meiosis. We deleted six lincRNA genes that are highly and dynamically expressed in the C. elegans gonad and tested the effects on central meiotic processes. Surprisingly, whereas the lincRNA deletions did not strongly impact fertility, germline apoptosis, crossovers, or synapsis, linc-4 was required for somatic growth. Slower growth was observed in linc-4-deletion mutants and in worms depleted of linc-4 using RNAi, indicating that linc-4 transcripts are required for this post-embryonic process. Unexpectedly, analysis of worms depleted of linc-4 in soma versus germline showed that the somatic role stems from linc-4 expression in germline cells. This unique feature suggests that some lincRNAs, like some small non-coding RNAs, are required for germ-soma interactions.

SpatialDB: a database for spatially resolved transcriptomes

Spatially resolved transcriptomic techniques allow the characterization of spatial organization of cells in tissues, which revolutionize the studies of tissue function and disease pathology. New strategies for detecting spatial gene expression patterns are emerging, and spatially resolved transcriptomic data are accumulating rapidly. However, it is not convenient for biologists to exploit these data due to the diversity of strategies and complexity in data analysis. Here, we present SpatialDB, the first manually curated database for spatially resolved transcriptomic techniques and datasets. The current version of SpatialDB contains 24 datasets (305 sub-datasets) from 5 species generated by 8 spatially resolved transcriptomic techniques. SpatialDB provides a user-friendly web interface for visualization and comparison of spatially resolved transcriptomic data. To further explore these data, SpatialDB also provides spatially variable genes and their functional enrichment annotation. SpatialDB offers a repository for research community to investigate the spatial cellular structure of tissues, and may bring new insights into understanding the cellular microenvironment in disease. SpatialDB is freely available at https://www.spatialomics.org/SpatialDB.

SLC17A6/7/8 Vesicular Glutamate Transporter Homologs in Nematodes

Members of the superfamily of solute carrier (SLC) transmembrane proteins transport diverse substrates across distinct cellular membranes. Three SLC protein families transport distinct neurotransmitters into synaptic vesicles to enable synaptic transmission in the nervous system. Among them is the SLC17A6/7/8 family of vesicular glutamate transporters, which endows specific neuronal cell types with the ability to use glutamate as a neurotransmitter. The genome of the nematode Caenorhabditis elegans encodes three SLC17A6/7/8 family members, one of which, eat-4/VGLUT, has been shown to be involved in glutamatergic neurotransmission. Here, we describe our analysis of the two remaining, previously uncharacterized SLC17A6/7/8 family members, vglu-2 and vglu-3 These two genes directly neighbor one another and are the result of a recent gene duplication event in C. elegans, but not in other Caenorhabditis species. Compared to EAT-4, the VGLU-2 and VGLU-3 protein sequences display a more distant similarity to canonical, vertebrate VGLUT proteins. We tagged both genomic loci with gfp and detected no expression of vglu-3 at any stage of development in any cell type of both C. elegans sexes. In contrast, vglu-2::gfp is dynamically expressed in a restricted set of distinct cell types. Within the nervous system, vglu-2::gfp is exclusively expressed in a single interneuron class, AIA, where it localizes to vesicular structures in the soma, but not along the axon, suggesting that VGLU-2 may not be involved in synaptic transport of glutamate. Nevertheless, vglu-2 mutants are partly defective in the function of the AIA neuron in olfactory behavior. Outside the nervous system, VGLU-2 is expressed in collagen secreting skin cells where VGLU-2 most prominently localizes to early endosomes, and to a lesser degree to apical clathrin-coated pits, the trans-Golgi network, and late endosomes. On early endosomes, VGLU-2 colocalizes most strongly with the recycling promoting factor SNX-1, a retromer component. Loss of vglu-2 affects the permeability of the collagen-containing cuticle of the worm, and based on the function of a vertebrate VGLUT1 protein in osteoclasts, we speculate that vglu-2 may have a role in collagen trafficking in the skin. We conclude that C. elegans SLC17A6/7/8 family members have diverse functions within and outside the nervous system.

Defects in mating behavior and tail morphology are the primary cause of sterility in Caenorhabditis elegans males at high temperature

Reproduction is a fundamental imperative of all forms of life. For all the advantages sexual reproduction confers, it has a deeply conserved flaw: it is temperature sensitive. As temperatures rise, fertility decreases. Across species, male fertility is particularly sensitive to elevated temperature. Previously, we have shown in the model nematode Caenorhabditis elegans that all males are fertile at 20°C, but almost all males have lost fertility at 27°C. Male fertility is dependent on the production of functional sperm, successful mating and transfer of sperm, and successful fertilization post-mating. To determine how male fertility is impacted by elevated temperature, we analyzed these aspects of male reproduction at 27°C in three wild-type strains of C. elegans: JU1171, LKC34 and N2. We found no effect of elevated temperature on the number of immature non-motile spermatids formed. There was only a weak effect of elevated temperature on sperm activation. In stark contrast, there was a strong effect of elevated temperature on male mating behavior, male tail morphology and sperm transfer such that males very rarely completed mating successfully when exposed to 27°C. Therefore, we propose a model where elevated temperature reduces male fertility as a result of the negative impacts of temperature on the somatic tissues necessary for mating. Loss of successful mating at elevated temperature overrides any effects that temperature may have on the germline or sperm cells.

A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution

Caenorhabditis elegans is an animal with few cells but a wide diversity of cell types. In this study, we characterize the molecular basis for their specification by profiling the transcriptomes of 86,024 single embryonic cells. We identify 502 terminal and preterminal cell types, mapping most single-cell transcriptomes to their exact position in C. elegans' invariant lineage. Using these annotations, we find that (i) the correlation between a cell's lineage and its transcriptome increases from middle to late gastrulation, then falls substantially as cells in the nervous system and pharynx adopt their terminal fates; (ii) multilineage priming contributes to the differentiation of sister cells at dozens of lineage branches; and (iii) most distinct lineages that produce the same anatomical cell type converge to a homogenous transcriptomic state.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Insulin-like peptides and the mTOR-TFEB pathway protect Caenorhabditis elegans hermaphrodites from mating-induced death

Lifespan is shortened by mating, but these deleterious effects must be delayed long enough for successful reproduction. Susceptibility to brief mating-induced death is caused by the loss of protection upon self-sperm depletion. Self-sperm maintains the expression of a DAF-2 insulin-like antagonist, INS-37, which promotes the nuclear localization of intestinal HLH-30/TFEB, a key pro-longevity regulator. Mating induces the agonist INS-8, promoting HLH-30 nuclear exit and subsequent death. In opposition to the protective role of HLH-30 and DAF-16/FOXO, TOR/LET-363 and the IIS-regulated Zn-finger transcription factor PQM-1 promote seminal-fluid-induced killing. Self-sperm maintenance of nuclear HLH-30/TFEB allows hermaphrodites to resist mating-induced death until self-sperm are exhausted, increasing the chances that mothers will survive through reproduction. Mothers combat males' hijacking of their IIS pathway by expressing an insulin antagonist that keeps her healthy through the activity of pro-longevity factors, as long as she has her own sperm to utilize.

Spatiotemporal m(i)RNA Architecture and 3' UTR Regulation in the C. elegans Germline

In animal germlines, regulation of cell proliferation and differentiation is particularly important but poorly understood. Here, using a cryo-cut approach, we mapped RNA expression along the Caenorhabditis elegans germline and, using mutants, dissected gene regulatory mechanisms that control spatiotemporal expression. We detected, at near single-cell resolution, >10,000 mRNAs, >300 miRNAs, and numerous unannotated miRNAs. Most RNAs were organized in distinct spatial patterns. Germline-specific miRNAs and their targets were co-localized. Moreover, we observed differential 3' UTR isoform usage for hundreds of mRNAs. In tumorous gld-2 gld-1 mutants, gene expression was strongly perturbed. In particular, differential 3' UTR usage was significantly impaired. We propose that PIE-1, a transcriptional repressor, functions to maintain spatial gene expression. Our data also suggest that cpsf-4 and fipp-1 control differential 3' UTR usage for hundreds of genes. Finally, we constructed a "virtual gonad" enabling "virtual in situ hybridizations" and access to all data (https://shiny.mdc-berlin.de/spacegerm/).