A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis

Cell lineage-resolved embryonic morphological map reveals signaling associated with cell fate and size asymmetry

How cells change shape is crucial for the development of tissues, organs and embryos. However, studying these shape changes in detail is challenging. Here we present a comprehensive real-time cellular map that covers over 95% of the cells formed during Caenorhabditis elegans embryogenesis, featuring nearly 400,000 3D cell regions. This map includes information on each cell's identity, lineage, fate, shape, volume, surface area, contact area, and gene expression profiles, all accessible through our user-friendly software and website. Our map allows for detailed analysis of key developmental processes, including dorsal intercalation, intestinal formation, and muscle assembly. We show how Notch and Wnt signaling pathways, along with mechanical forces from cell interactions, regulate cell fate decisions and size asymmetries. Our findings suggest that repeated Notch signaling drives size disparities in the large excretory cell, which functions like a kidney. This work sets the stage for in-depth studies of the mechanisms controlling cell fate differentiation and morphogenesis.

CelEst: a unified gene regulatory network for estimating transcription factor activities in C. elegans

Transcription factors (TFs) play a pivotal role in orchestrating critical intricate patterns of gene regulation. Although gene expression is complex, differential expression of hundreds of genes is often due to regulation by just a handful of TFs. Despite extensive efforts to elucidate TF-target regulatory relationships in Caenorhabditis elegans, existing experimental datasets cover distinct subsets of TFs and leave data integration challenging. Here, I introduce CelEst, a unified gene regulatory network designed to estimate the activity of 487 distinct C. elegans TFs-∼58% of the total-from gene expression data. To integrate data from ChIP-seq, DNA-binding motifs, and eY1H screens, optimal processing of each data type was benchmarked against a set of TF perturbation RNA-seq experiments. Moreover, I showcase how leveraging TF motif conservation in target promoters across genomes of related species can distinguish highly informative interactions, a strategy which can be applied to many model organisms. Integrated analyses of data from commonly studied conditions including heat shock, bacterial infection, and sex differences validates CelEst's performance and highlights overlooked TFs that likely play major roles in coordinating the transcriptional response to these conditions. CelEst can infer TF activity on a standard laptop computer within minutes. Furthermore, an R Shiny app with a step-by-step guide is provided for the community to perform rapid analysis with minimal coding required. I anticipate that widespread adoption of CelEsT will significantly enhance the interpretive power of transcriptomic experiments, both present and retrospective, thereby advancing our understanding of gene regulation in C. elegans and beyond.

EFL-3/E2F7 modulates Wnt signalling by repressing the Nemo-like kinase LIT-1 during asymmetric epidermal cell division in Caenorhabditis elegans

The E2F family of transcription factors is conserved in higher eukaryotes and plays pivotal roles in controlling gene expression during the cell cycle. Most canonical E2Fs associate with members of the Dimerisation Partner (DP) family to activate or repress target genes. However, atypical repressors, such as E2F7 and E2F8, lack DP interaction domains and their functions are less understood. We report here that EFL-3, the E2F7 homologue of Caenorhabditis elegans, regulates epidermal stem cell differentiation. We show that phenotypic defects in efl-3 mutants depend on the Nemo-like kinase LIT-1. EFL-3 represses lit-1 expression through direct binding to a lit-1 intronic element. Increased LIT-1 expression in efl-3 mutants reduces POP-1/TCF nuclear distribution, and consequently alters Wnt pathway activation. Our findings provide a mechanistic link between an atypical E2F family member and NLK during C. elegans asymmetric cell division, which may be conserved in other animals.

Lipid Metabolic Heterogeneity during Early Embryogenesis Revealed by Hyper-3D Stimulated Raman Imaging

Studying embryogenesis is fundamental to understanding developmental biology and reproductive medicine. Its process requires precise spatiotemporal regulations in which lipid metabolism plays a crucial role. However, the spatial dynamics of lipid species at the subcellular level remains obscure due to technical limitations. To address this challenge, we developed a hyperspectral 3D imaging and analysis method based on stimulated Raman scattering microscopy (hyper-3D SRS) to quantitatively assess lipid profiles in individual embryos through submicrometer resolution (x-y), 3D optical sectioning (z), and chemical bond-selective (Ω) imaging. Using hyper-3D SRS, individual lipid droplets (LDs) in single cells were identified and quantified. Our findings revealed that the LD profiles within a single embryo are not uniform, even as early as the 2-cell stage. Notably, we also discovered a dynamic relationship between the LD size and unsaturation degree as embryos develop, indicating diverse lipid metabolism during early development. Furthermore, abnormal LDs were observed in oocytes of a progeria mouse model, suggesting that LDs could serve as a potential biomarker for assessing oocyte/embryo quality. Overall, our results highlight the potential of hyper-3D SRS as a noninvasive method for studying lipid content, composition, and subcellular distribution in embryos. This technique provides valuable insights into lipid metabolism during embryonic development and has the potential for clinical applications in evaluating oocyte/embryo quality.

SEM-2/SoxC regulates multiple aspects of C. elegans postembryonic mesoderm development

Development of multicellular organisms requires well-orchestrated interplay between cell-intrinsic transcription factors and cell-cell signaling. One set of highly conserved transcription factors that plays diverse roles in development is the SoxC group. C. elegans contains a sole SoxC protein, SEM-2. SEM-2 is essential for embryonic development, and for specifying the sex myoblast (SM) fate in the postembryonic mesoderm, the M lineage. We have identified a novel partial loss-of-function sem-2 allele that has a proline to serine change in the C-terminal tail of the highly conserved DNA-binding domain. Detailed analyses of mutant animals harboring this point mutation uncovered new functions of SEM-2 in the M lineage. First, SEM-2 functions antagonistically with LET-381, the sole C. elegans FoxF/C forkhead transcription factor, to regulate dorsoventral patterning of the M lineage. Second, in addition to specifying the SM fate, SEM-2 is essential for the proliferation and diversification of the SM lineage. Finally, SEM-2 appears to directly regulate the expression of hlh-8, which encodes a basic helix-loop-helix Twist transcription factor and plays critical roles in proper patterning of the M lineage. Our data, along with previous studies, suggest an evolutionarily conserved relationship between SoxC and Twist proteins. Furthermore, our work identified new interactions in the gene regulatory network (GRN) underlying C. elegans postembryonic development and adds to the general understanding of the structure-function relationship of SoxC proteins.

Binding profiles for 961 Drosophila and C. elegans transcription factors reveal tissue-specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here, we present the culmination of the efforts of the modENCODE (model organism Encyclopedia of DNA Elements) and modERN (model organism Encyclopedia of Regulatory Networks) consortia to systematically assay TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). These data sets comprise 605 TFs identifying 3.6 M sites in the fly and 356 TFs identifying 0.9 M sites in the worm, and represent the majority of the regulatory space in each genome. We demonstrate that TFs associate with chromatin in clusters termed "metapeaks," that larger metapeaks have characteristics of high-occupancy target (HOT) regions, and that the importance of consensus sequence motifs bound by TFs depends on metapeak size and complexity. Combining ChIP-seq data with single-cell RNA-seq data in a machine-learning model identifies TFs with a prominent role in promoting target gene expression in specific cell types, even differentiating between parent-daughter cells during embryogenesis. These data are a rich resource for the community that should fuel and guide future investigations into TF function. To facilitate data accessibility and utility, all strains expressing green fluorescent protein (GFP)-tagged TFs are available at the stock centers for each organism. The chromatin immunoprecipitation sequencing data are available through the ENCODE Data Coordinating Center, GEO, and through a direct interface that provides rapid access to processed data sets and summary analyses, as well as widgets to probe the cell-type-specific TF-target relationships.

Analysis of Haemonchus embryos at single cell resolution identifies two eukaryotic elongation factors as intervention target candidates

Advances in single cell technologies are allowing investigations of a wide range of biological processes and pathways in animals, such as the multicellular model organism Caenorhabditis elegans - a free-living nematode. However, there has been limited application of such technology to related parasitic nematodes which cause major diseases of humans and animals worldwide. With no vaccines against the vast majority of parasitic nematodes and treatment failures due to drug resistance or inefficacy, new intervention targets are urgently needed, preferably informed by a deep understanding of these nematodes' cellular and molecular biology - which is presently lacking for most worms. Here, we created the first single cell atlas for an early developmental stage of Haemonchus contortus - a highly pathogenic, C. elegans-related parasitic nematode. We obtained and curated RNA sequence (snRNA-seq) data from single nuclei from embryonating eggs of H. contortus (150,000 droplets), and selected high-quality transcriptomic data for > 14,000 single nuclei for analysis, and identified 19 distinct clusters of cells. Guided by comparative analyses with C. elegans, we were able to reproducibly assign seven cell clusters to body wall muscle, hypodermis, neuronal, intestinal or seam cells, and identified eight genes that were transcribed in all cell clusters/types, three of which were inferred to be essential in H. contortus. Two of these genes (i.e. Hc-eef-1A and Hc-eef1G), coding for eukaryotic elongation factors (called Hc-eEF1A and Hc-eEF1G), were also demonstrated to be transcribed and expressed in all key developmental stages of H. contortus. Together with these findings, sequence- and structure-based comparative analyses indicated the potential of Hc-eEF1A and/or Hc-eEF1G as intervention targets within the protein biosynthesis machinery of H. contortus. Future work will focus on single cell studies of all key developmental stages and tissues of H. contortus, and on evaluating the suitability of the two elongation factor proteins as drug targets in H. contortus and related nematodes, with a view to finding new nematocidal drug candidates.

Automated cell lineage reconstruction using label-free 4D microscopy

Patterns of lineal descent play a critical role in the development of metazoan embryos. In eutelic organisms that generate a fixed number of somatic cells, invariance in the topology of their cell lineage provides a powerful opportunity to interrogate developmental events with empirical repeatability across individuals. Studies of embryonic development using the nematode Caenorhabditis elegans have been drivers of discovery. These studies have depended heavily on high-throughput lineage tracing enabled by 4D fluorescence microscopy and robust computer vision pipelines. For a range of applications, computer-aided yet manual lineage tracing using 4D label-free microscopy remains an essential tool. Deep learning approaches to cell detection and tracking in fluorescence microscopy have advanced significantly in recent years, yet solutions for automating cell detection and tracking in 3D label-free imaging of dense tissues and embryos remain inaccessible. Here, we describe embGAN, a deep learning pipeline that addresses the challenge of automated cell detection and tracking in label-free 3D time-lapse imaging. embGAN requires no manual data annotation for training, learns robust detections that exhibits a high degree of scale invariance, and generalizes well to images acquired in multiple labs on multiple instruments. We characterize embGAN's performance using lineage tracing in the C. elegans embryo as a benchmark. embGAN achieves near-state-of-the-art performance in cell detection and tracking, enabling high-throughput studies of cell lineage without the need for fluorescent reporters or transgenics.

Neurogenesis in Caenorhabditis elegans

Animals rely on their nervous systems to process sensory inputs, integrate these with internal signals, and produce behavioral outputs. This is enabled by the highly specialized morphologies and functions of neurons. Neuronal cells share multiple structural and physiological features, but they also come in a large diversity of types or classes that give the nervous system its broad range of functions and plasticity. This diversity, first recognized over a century ago, spurred classification efforts based on morphology, function, and molecular criteria. Caenorhabditis elegans, with its precisely mapped nervous system at the anatomical level, an extensive molecular description of most of its neurons, and its genetic amenability, has been a prime model for understanding how neurons develop and diversify at a mechanistic level. Here, we review the gene regulatory mechanisms driving neurogenesis and the diversification of neuron classes and subclasses in C. elegans. We discuss our current understanding of the specification of neuronal progenitors and their differentiation in terms of the transcription factors involved and ensuing changes in gene expression and chromatin landscape. The central theme that has emerged is that the identity of a neuron is defined by modules of gene batteries that are under control of parallel yet interconnected regulatory mechanisms. We focus on how, to achieve these terminal identities, cells integrate information along their developmental lineages. Moreover, we discuss how neurons are diversified postembryonically in a time-, genetic sex-, and activity-dependent manner. Finally, we discuss how the understanding of neuronal development can provide insights into the evolution of neuronal diversity.

Tissue-specific RNA-seq defines genes governing male tail tip morphogenesis in C. elegans

Caenorhabditis elegans males undergo sex-specific tail tip morphogenesis (TTM) under the control of the DM-domain transcription factor DMD-3. To find genes regulated by DMD-3, we performed RNA-seq of laser-dissected tail tips. We identified 564 genes differentially expressed (DE) in wild-type males versus dmd-3(-) males and hermaphrodites. The transcription profile of dmd-3(-) tail tips is similar to that in hermaphrodites. For validation, we analyzed transcriptional reporters for 49 genes and found male-specific or male-biased expression for 26 genes. Only 11 DE genes overlapped with genes found in a previous RNAi screen for defective TTM. GO enrichment analysis of DE genes finds upregulation of genes within the unfolded protein response pathway and downregulation of genes involved in cuticle maintenance. Of the DE genes, 40 are transcription factors, indicating that the gene network downstream of DMD-3 is complex and potentially modular. We propose modules of genes that act together in TTM and are co-regulated by DMD-3, among them the chondroitin synthesis pathway and the hypertonic stress response.

SEM-2/SoxC regulates multiple aspects of C. elegans postembryonic mesoderm development

Development of multicellular organisms requires well-orchestrated interplay between cell-intrinsic transcription factors and cell-cell signaling. One set of highly conserved transcription factors that plays diverse roles in development is the SoxC group. C. elegans contains a sole SoxC protein, SEM-2. SEM-2 is essential for embryonic development, and for specifying the sex myoblast (SM) fate in the postembryonic mesoderm, the M lineage. We have identified a novel partial loss-of-function sem-2 allele that has a proline to serine change in the C-terminal tail of the highly conserved DNA-binding domain. Detailed analyses of mutant animals harboring this point mutation uncovered new functions of SEM-2 in the M lineage. First, SEM-2 functions antagonistically with LET-381, the sole C. elegans FoxF/C forkhead transcription factor, to regulate dorsoventral patterning of the M lineage. Second, in addition to specifying the SM fate, SEM-2 is essential for the proliferation and diversification of the SM lineage. Finally, SEM-2 appears to directly regulate the expression of hlh-8, which encodes a basic helix-loop-helix Twist transcription factor and plays critical roles in proper patterning of the M lineage. Our data, along with previous studies, suggest an evolutionarily conserved relationship between SoxC and Twist proteins. Furthermore, our work identified new interactions in the gene regulatory network (GRN) underlying C. elegans postembryonic development and adds to the general understanding of the structure-function relationship of SoxC proteins.

Characterizing Cellular Physiological States with Three-Dimensional Shape Descriptors for Cell Membranes

The shape of a cell as defined by its membrane can be closely associated with its physiological state. For example, the irregular shapes of cancerous cells and elongated shapes of neuron cells often reflect specific functions, such as cell motility and cell communication. However, it remains unclear whether and which cell shape descriptors can characterize different cellular physiological states. In this study, 12 geometric shape descriptors for a three-dimensional (3D) object were collected from the previous literature and tested with a public dataset of ~400,000 independent 3D cell regions segmented based on fluorescent labeling of the cell membranes in Caenorhabditis elegans embryos. It is revealed that those shape descriptors can faithfully characterize cellular physiological states, including (1) cell division (cytokinesis), along with an abrupt increase in the elongation ratio; (2) a negative correlation of cell migration speed with cell sphericity; (3) cell lineage specification with symmetrically patterned cell shape changes; and (4) cell fate specification with differential gene expression and differential cell shapes. The descriptors established may be used to identify and predict the diverse physiological states in numerous cells, which could be used for not only studying developmental morphogenesis but also diagnosing human disease (e.g., the rapid detection of abnormal cells).

A lineage-resolved cartography of microRNA promoter activity in C. elegans empowers multidimensional developmental analysis

Elucidating the expression of microRNAs in developing single cells is critical for functional discovery. Here, we construct scCAMERA (single-cell cartography of microRNA expression based on reporter assay), utilizing promoter-driven fluorescent reporters in conjunction with imaging and lineage tracing. The cartography delineates the transcriptional activity of 54 conserved microRNAs in lineage-resolved single cells throughout C. elegans embryogenesis. The combinatorial expression of microRNAs partitions cells into fine clusters reflecting their function and anatomy. Notably, the expression of individual microRNAs exhibits high cell specificity and divergence among family members. Guided by cellular expression patterns, we identify developmental functions of specific microRNAs, including miR-1 in pharynx development and physiology, miR-232 in excretory canal morphogenesis by repressing NHR-25/NR5A, and a functional synergy between miR-232 and miR-234 in canal development, demonstrating the broad utility of scCAMERA. Furthermore, integrative analysis reveals that tissue-specific fate determinants activate microRNAs to repress protein production from leaky transcripts associated with alternative, especially neuronal, fates, thereby enhancing the fidelity of developmental fate differentiation. Collectively, our study offers rich opportunities for multidimensional expression-informed analysis of microRNA biology in metazoans.

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.

Gene regulatory patterning codes in early cell fate specification of the C. elegans embryo

Pattern formation originates during embryogenesis by a series of symmetry-breaking steps throughout an expanding cell lineage. In Drosophila, classic work has shown that segmentation in the embryo is established by morphogens within a syncytium, and the subsequent action of the gap, pair-rule, and segment polarity genes. This classic model however does not translate directly to species that lack a syncytium - such as Caenorhabditis elegans - where cell fate is specified by cell-autonomous cell lineage programs and their inter-signaling. Previous single-cell RNA-Seq studies in C. elegans have analyzed cells from a mixed suspension of cells from many embryos to study late differentiation stages, or individual early stage embryos to study early gene expression in the embryo. To study the intermediate stages of early and late gastrulation (28- to 102-cells stages) missed by these approaches, here we determine the transcriptomes of the 1- to 102-cell stage to identify 119 embryonic cell states during cell fate specification, including 'equivalence-group' cell identities. We find that gene expression programs are modular according to the sub-cell lineages, each establishing a set of stripes by combinations of transcription factor gene expression across the anterior-posterior axis. In particular, expression of the homeodomain genes establishes a comprehensive lineage-specific positioning system throughout the embryo beginning at the 28-cell stage. Moreover, we find that genes that segment the entire embryo in Drosophila have orthologs in C. elegans that exhibit sub-lineage-specific expression. These results suggest that the C. elegans embryo is patterned by a juxtaposition of distinct lineage-specific gene regulatory programs each with a unique encoding of cell location and fate. This use of homologous gene regulatory patterning codes suggests a deep homology of cell fate specification programs across diverse modes of development.

Automated Cell Lineage Reconstruction using Label-Free 4D Microscopy

Here we describe embGAN, a deep learning pipeline that addresses the challenge of automated cell detection and tracking in label-free 3D time lapse imaging. embGAN requires no manual data annotation for training, learns robust detections that exhibits a high degree of scale invariance and generalizes well to images acquired in multiple labs on multiple instruments.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

During embryonic development, cells undergo dynamic changes in gene expression that are required for appropriate cell fate specification. Although both transcription and mRNA degradation contribute to gene expression dynamics, patterns of mRNA decay are less well-understood. Here we directly measured spatiotemporally resolved mRNA decay rates transcriptome-wide throughout C. elegans embryogenesis by transcription inhibition followed by bulk and single-cell RNA-sequencing. This allowed us to calculate mRNA half-lives within specific cell types and developmental stages and identify differentially regulated mRNA decay throughout embryonic development. We identified transcript features that are correlated with mRNA stability and found that mRNA decay rates are associated with distinct peaks in gene expression over time. Moreover, we provide evidence that, on average, mRNA is more stable in the germline compared to in the soma and in later embryonic stages compared to in earlier stages. This work suggests that differential mRNA decay across cell states and time helps to shape developmental gene expression, and it provides a valuable resource for studies of mRNA turnover regulatory mechanisms.

Tissue-specific RNA-seq defines genes governing male tail tip morphogenesis in C. elegans

Caenorhabditis elegans males undergo sex-specific tail tip morphogenesis (TTM) under the control of the transcription factor DMD-3. To find genes regulated by DMD-3, We performed RNA-seq of laser-dissected tail tips. We identified 564 genes differentially expressed (DE) in wild-type males vs. dmd-3(-) males and hermaphrodites. The transcription profile of dmd-3(-) tail tips is similar to that in hermaphrodites. For validation, we analyzed transcriptional reporters for 49 genes and found male-specific or male-biased expression for 26 genes. Only 11 DE genes overlapped with genes found in a previous RNAi screen for defective TTM. GO enrichment analysis of DE genes finds upregulation of genes within the UPR (unfolded protein response) pathway and downregulation of genes involved in cuticle maintenance. Of the DE genes, 40 are transcription factors, indicating that the gene network downstream of DMD-3 is complex and potentially modular. We propose modules of genes that act together in TTM and are coregulated by DMD-3, among them the chondroitin synthesis pathway and the hypertonic stress response.

A full-body transcription factor expression atlas with completely resolved cell identities in C. elegans

Invariant cell lineage in C. elegans enables spatiotemporal resolution of transcriptional regulatory mechanisms controlling the fate of each cell. Here, we develop RAPCAT (Robust-point-matching- And Piecewise-affine-based Cell Annotation Tool) to automate cell identity assignment in three-dimensional image stacks of L1 larvae and profile reporter expression of 620 transcription factors in every cell. Transcription factor profile-based clustering analysis defines 80 cell types distinct from conventional phenotypic cell types and identifies three general phenotypic modalities related to these classifications. First, transcription factors are broadly downregulated in quiescent stage Hermaphrodite Specific Neurons, suggesting stage- and cell type-specific variation in transcriptome size. Second, transcription factor expression is more closely associated with morphology than other phenotypic modalities in different pre- and post-differentiation developmental stages. Finally, embryonic cell lineages can be associated with specific transcription factor expression patterns and functions that persist throughout postembryonic life. This study presents a comprehensive transcription factor atlas for investigation of intra-cell type heterogeneity.

On the origins of developmental robustness: modeling buffering mechanisms against cell-level noise

During development, cells are subject to stochastic fluctuations in their positions (i.e. cell-level noise) that can potentially lead to morphological noise (i.e. stochastic differences between morphologies that are expected to be equal, e.g. the right and left sides of bilateral organisms). In this study, we explore new and existing hypotheses on buffering mechanisms against cell-level noise. Many of these hypotheses focus on how the boundaries between territories of gene expression remain regular and well defined, despite cell-level noise and division. We study these hypotheses and how irregular territory boundaries lead to morphological noise. To determine the consistency of the different hypotheses, we use a general computational model of development: EmbryoMaker. EmbryoMaker can implement arbitrary gene networks regulating basic cell behaviors (contraction, adhesion, etc.), signaling and tissue biomechanics. We found that buffering mechanisms based on the orientation of cell divisions cannot lead to regular boundaries but that other buffering mechanisms can (homotypic adhesion, planar contraction, non-dividing boundaries, constant signaling and majority rule hypotheses). We also explore the effects of the shape and size of the territories on morphological noise.

Metabolic plasticity sustains the robustness of Caenorhabditis elegans embryogenesis

Embryogenesis is highly dependent on maternally loaded materials, particularly those used for energy production. Different environmental conditions and genetic backgrounds shape embryogenesis. The robustness of embryogenesis in response to extrinsic and intrinsic changes remains incompletely understood. By analyzing the levels of two major nutrients, glycogen and neutral lipids, we discovered stage-dependent usage of these two nutrients along with mitochondrial morphology changes during Caenorhabditis elegans embryogenesis. ATGL, the rate-limiting lipase in cellular lipolysis, is expressed and required in the hypodermis to regulate mitochondrial function and support embryogenesis. The embryonic lethality of atgl-1 mutants can be suppressed by reducing sinh-1/age-1-akt signaling, likely through modulating glucose metabolism to maintain sustainable glucose consumption. The embryonic lethality of atgl-1(xd314) is also affected by parental nutrition. Parental glucose and oleic acid supplements promote glycogen storage in atgl-1(xd314) embryos to compensate for the impaired lipolysis. The rescue by parental vitamin B12 supplement is likely through enhancing mitochondrial function in atgl-1 mutants. These findings reveal that metabolic plasticity contributes to the robustness of C. elegans embryogenesis.

Mechanisms of lineage specification in Caenorhabditis elegans

The studies of cell fate and lineage specification are fundamental to our understanding of the development of multicellular organisms. Caenorhabditis elegans has been one of the premiere systems for studying cell fate specification mechanisms at single cell resolution, due to its transparent nature, the invariant cell lineage, and fixed number of somatic cells. We discuss the general themes and regulatory mechanisms that have emerged from these studies, with a focus on somatic lineages and cell fates. We next review the key factors and pathways that regulate the specification of discrete cells and lineages during embryogenesis and postembryonic development; we focus on transcription factors and include numerous lineage diagrams that depict the expression of key factors that specify embryonic founder cells and postembryonic blast cells, and the diverse somatic cell fates they generate. We end by discussing some future perspectives in cell and lineage specification.

Novel metrics reveal new structure and unappreciated heterogeneity in Caenorhabditis elegans development

High throughput experimental approaches are increasingly allowing for the quantitative description of cellular and organismal phenotypes. Distilling these large volumes of complex data into meaningful measures that can drive biological insight remains a central challenge. In the quantitative study of development, for instance, one can resolve phenotypic measures for single cells onto their lineage history, enabling joint consideration of heritable signals and cell fate decisions. Most attempts to analyze this type of data, however, discard much of the information content contained within lineage trees. In this work we introduce a generalized metric, which we term the branch edit distance, that allows us to compare any two embryos based on phenotypic measurements in individual cells. This approach aligns those phenotypic measurements to the underlying lineage tree, providing a flexible and intuitive framework for quantitative comparisons between, for instance, Wild-Type (WT) and mutant developmental programs. We apply this novel metric to data on cell-cycle timing from over 1300 WT and RNAi-treated Caenorhabditis elegans embryos. Our new metric revealed surprising heterogeneity within this data set, including subtle batch effects in WT embryos and dramatic variability in RNAi-induced developmental phenotypes, all of which had been missed in previous analyses. Further investigation of these results suggests a novel, quantitative link between pathways that govern cell fate decisions and pathways that pattern cell cycle timing in the early embryo. Our work demonstrates that the branch edit distance we propose, and similar metrics like it, have the potential to revolutionize our quantitative understanding of organismal phenotype.

Dynamic compartmentalization of the pro-invasive transcription factor NHR-67 reveals a role for Groucho in regulating a proliferative-invasive cellular switch in C. elegans

A growing body of evidence suggests that cell division and basement membrane invasion are mutually exclusive cellular behaviors. How cells switch between proliferative and invasive states is not well understood. Here, we investigated this dichotomy in vivo by examining two cell types in the developing Caenorhabditis elegans somatic gonad that derive from equipotent progenitors, but exhibit distinct cell behaviors: the post-mitotic, invasive anchor cell and the neighboring proliferative, non-invasive ventral uterine (VU) cells. We show that the fates of these cells post-specification are more plastic than previously appreciated and that levels of NHR-67 are important for discriminating between invasive and proliferative behavior. Transcription of NHR-67 is downregulated following post-translational degradation of its direct upstream regulator, HLH-2 (E/Daughterless) in VU cells. In the nuclei of VU cells, residual NHR-67 protein is compartmentalized into discrete punctae that are dynamic over the cell cycle and exhibit liquid-like properties. By screening for proteins that colocalize with NHR-67 punctae, we identified new regulators of uterine cell fate maintenance: homologs of the transcriptional co-repressor Groucho (UNC-37 and LSY-22), as well as the TCF/LEF homolog POP-1. We propose a model in which the association of NHR-67 with the Groucho/TCF complex suppresses the default invasive state in non-invasive cells, which complements transcriptional regulation to add robustness to the proliferative-invasive cellular switch in vivo.

Recent methodological advances towards single-cell proteomics

Studying the central dogma at the single-cell level has gained increasing attention to reveal hidden cell lineages and functions that cannot be studied using traditional bulk analyses. Nonetheless, most single-cell studies exploiting genomic and transcriptomic levels fail to address information on proteins that are central to many important biological processes. Single-cell proteomics enables understanding of the functional status of individual cells and is particularly crucial when the specimen is composed of heterogeneous entities of cells. With the growing importance of this field, significant methodological advancements have emerged recently. These include miniaturized and automated sample preparation, multi-omics analyses, and combined analyses of multiple techniques such as mass spectrometry and microscopy. Moreover, artificial intelligence and single-molecule detection technologies have advanced throughput and improved sensitivity limitations, respectively, over conventional methods. In this review, we summarize cutting-edge methodologies for single-cell proteomics and relevant emerging technologies that have been reported in the last 5 years, and provide an outlook on this research field.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

Novel metrics reveal new structure and unappreciated heterogeneity in C. elegans development

High throughput experimental approaches are increasingly allowing for the quantitative description of cellular and organismal phenotypes. Distilling these large volumes of complex data into meaningful measures that can drive biological insight remains a central challenge. In the quantitative study of development, for instance, one can resolve phenotypic measures for single cells onto their lineage history, enabling joint consideration of heritable signals and cell fate decisions. Most attempts to analyze this type of data, however, discard much of the information content contained within lineage trees. In this work we introduce a generalized metric, which we term the branch distance, that allows us to compare any two embryos based on phenotypic measurements in individual cells. This approach aligns those phenotypic measurements to the underlying lineage tree, providing a flexible and intuitive framework for quantitative comparisons between, for instance, Wild-Type (WT) and mutant developmental programs. We apply this novel metric to data on cell-cycle timing from over 1300 WT and RNAi-treated Caenorhabditis elegans embryos. Our new metric revealed surprising heterogeneity within this data set, including subtle batch effects in WT embryos and dramatic variability in RNAi-induced developmental phenotypes, all of which had been missed in previous analyses. Further investigation of these results suggests a novel, quantitative link between pathways that govern cell fate decisions and pathways that pattern cell cycle timing in the early embryo. Our work demonstrates that the branch distance we propose, and similar metrics like it, have the potential to revolutionize our quantitative understanding of organismal phenotype.

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

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

Widespread employment of conserved C. elegans homeobox genes in neuronal identity specification

Homeobox genes are prominent regulators of neuronal identity, but the extent to which their function has been probed in animal nervous systems remains limited. In the nematode Caenorhabditis elegans, each individual neuron class is defined by the expression of unique combinations of homeobox genes, prompting the question of whether each neuron class indeed requires a homeobox gene for its proper identity specification. We present here progress in addressing this question by extending previous mutant analysis of homeobox gene family members and describing multiple examples of homeobox gene function in different parts of the C. elegans nervous system. To probe homeobox function, we make use of a number of reporter gene tools, including a novel multicolor reporter transgene, NeuroPAL, which permits simultaneous monitoring of the execution of multiple differentiation programs throughout the entire nervous system. Using these tools, we add to the previous characterization of homeobox gene function by identifying neuronal differentiation defects for 14 homeobox genes in 24 distinct neuron classes that are mostly unrelated by location, function and lineage history. 12 of these 24 neuron classes had no homeobox gene function ascribed to them before, while in the other 12 neuron classes, we extend the combinatorial code of transcription factors required for specifying terminal differentiation programs. Furthermore, we demonstrate that in a particular lineage, homeotic identity transformations occur upon loss of a homeobox gene and we show that these transformations are the result of changes in homeobox codes. Combining the present with past analyses, 113 of the 118 neuron classes of C. elegans are now known to require a homeobox gene for proper execution of terminal differentiation programs. Such broad deployment indicates that homeobox function in neuronal identity specification may be an ancestral feature of animal nervous systems.

Caenorhabditis elegans sine oculis/SIX-type homeobox genes act as homeotic switches to define neuronal subtype identities

The classification of neurons into distinct types reveals hierarchical taxonomic relationships that reflect the extent of similarity between neuronal cell types. At the base of such taxonomies are neuronal cells that are very similar to one another but differ in a small number of reproducible and select features. How are very similar members of a neuron class that share many features instructed to diversify into distinct subclasses? We show here that the six very similar members of the Caenorhabditis elegans IL2 sensory neuron class, which are all specified by a homeobox terminal selector, unc-86/BRN3, differentiate into two subtly distinct subclasses, a dorsoventral subclass and a lateral subclass, by the toggle switch-like action of the sine oculis/SIX homeobox gene unc-39. unc-39 is expressed only in the lateral IL2 neurons, and loss of unc-39 leads to a homeotic transformation of the lateral into the dorsoventral class; conversely, ectopic unc-39 expression converts the dorsoventral subclass into the lateral subclass. Hence, a terminal selector homeobox gene controls both class- as well as subclass-specific features, while a subordinate homeobox gene determines the ability of the class-specific homeobox gene to activate subtype-specific target genes. We find a similar regulatory mechanism operating in a distinct class of six motor neurons. Our findings underscore the importance of homeobox genes in neuronal identity control and invite speculations about homeotic identity transformations as potential drivers of evolutionary novelty during cell-type evolution in the brain.

pop-1/TCF, ref-2/ZIC and T-box factors regulate the development of anterior cells in the C. elegans embryo

Patterning of the anterior-posterior axis is fundamental to animal development. The Wnt pathway plays a major role in this process by activating the expression of posterior genes in animals from worms to humans. This observation raises the question of whether the Wnt pathway or other regulators control the expression of the many anterior-expressed genes. We found that the expression of five anterior-specific genes in Caenorhabditis elegans embryos depends on the Wnt pathway effectors pop-1/TCF and sys-1/β-catenin. We focused further on one of these anterior genes, ref-2/ZIC, a conserved transcription factor expressed in multiple anterior lineages. Live imaging of ref-2 mutant embryos identified defects in cell division timing and position in anterior lineages. Cis-regulatory dissection identified three ref-2 transcriptional enhancers, one of which is necessary and sufficient for anterior-specific expression. This enhancer is activated by the T-box transcription factors TBX-37 and TBX-38, and surprisingly, concatemerized TBX-37/38 binding sites are sufficient to drive anterior-biased expression alone, despite the broad expression of TBX-37 and TBX-38. Taken together, our results highlight the diverse mechanisms used to regulate anterior expression patterns in the embryo.

Defect-buffering cellular plasticity increases robustness of metazoan embryogenesis

Developmental processes are intrinsically robust so as to preserve a normal-like state in response to genetic and environmental fluctuations. However, the robustness and potential phenotypic plasticity of individual developing cells under genetic perturbations remain to be systematically evaluated. Using large-scale gene perturbation, live imaging, lineage tracing, and single-cell phenomics, we quantified the phenotypic landscape of C. elegans embryogenesis in >2,000 embryos following individual knockdown of over 750 conserved genes. We observed that cellular genetic systems are not sufficiently robust to single-gene perturbations across all cells; rather, gene knockdowns frequently induced cellular defects. Dynamic phenotypic analyses revealed many cellular defects to be transient, with cells exhibiting phenotypic plasticity that serves to alleviate, correct, and accommodate the defects. Moreover, potential developmentally related cell modules may buffer the phenotypic effects of individual cell position changes. Our findings reveal non-negligible contributions of cellular plasticity and multicellularity as compensatory strategies to increase developmental robustness.

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/).

The anterior Hox gene ceh-13 and elt-1/GATA activate the posterior Hox genes nob-1 and php-3 to specify posterior lineages in the C. elegans embryo

Hox transcription factors play a conserved role in specifying positional identity during animal development, with posterior Hox genes typically repressing the expression of more anterior Hox genes. Here, we dissect the regulation of the posterior Hox genes nob-1 and php-3 in the nematode C. elegans. We show that nob-1 and php-3 are co-expressed in gastrulation-stage embryos in cells that previously expressed the anterior Hox gene ceh-13. This expression is controlled by several partially redundant transcriptional enhancers. These enhancers act in a ceh-13-dependant manner, providing a striking example of an anterior Hox gene positively regulating a posterior Hox gene. Several other regulators also act positively through nob-1/php-3 enhancers, including elt-1/GATA, ceh-20/ceh-40/Pbx, unc-62/Meis, pop-1/TCF, ceh-36/Otx, and unc-30/Pitx. We identified defects in both cell position and cell division patterns in ceh-13 and nob-1;php-3 mutants, suggesting that these factors regulate lineage identity in addition to positional identity. Together, our results highlight the complexity and flexibility of Hox gene regulation and function and the ability of developmental transcription factors to regulate different targets in different stages of development.

The enteric nervous system of the C. elegans pharynx is specified by the Sine oculis-like homeobox gene ceh-34

Overarching themes in the terminal differentiation of the enteric nervous system, an autonomously acting unit of animal nervous systems, have so far eluded discovery. We describe here the overall regulatory logic of enteric nervous system differentiation of the nematode Caenorhabditis elegans that resides within the foregut (pharynx) of the worm. A C. elegans homolog of the Drosophila Sine oculis homeobox gene, ceh-34, is expressed in all 14 classes of interconnected pharyngeal neurons from their birth throughout their life time, but in no other neuron type of the entire animal. Constitutive and temporally controlled ceh-34 removal shows that ceh-34 is required to initiate and maintain the neuron type-specific terminal differentiation program of all pharyngeal neuron classes, including their circuit assembly. Through additional genetic loss of function analysis, we show that within each pharyngeal neuron class, ceh-34 cooperates with different homeodomain transcription factors to individuate distinct pharyngeal neuron classes. Our analysis underscores the critical role of homeobox genes in neuronal identity specification and links them to the control of neuronal circuit assembly of the enteric nervous system. Together with the pharyngeal nervous system simplicity as well as its specification by a Sine oculis homolog, our findings invite speculations about the early evolution of nervous systems.

Computable early Caenorhabditis elegans embryo with a phase field model

Morphogenesis is a precise and robust dynamic process during metazoan embryogenesis, consisting of both cell proliferation and cell migration. Despite the fact that much is known about specific regulations at molecular level, how cell proliferation and migration together drive the morphogenesis at cellular and organismic levels is not well understood. Using Caenorhabditis elegans as the model animal, we present a phase field model to compute early embryonic morphogenesis within a confined eggshell. With physical information about cell division obtained from three-dimensional time-lapse cellular imaging experiments, the model can precisely reproduce the early morphogenesis process as seen in vivo, including time evolution of location and morphology of each cell. Furthermore, the model can be used to reveal key cell-cell attractions critical to the development of C. elegans embryo. Our work demonstrates how genetic programming and physical forces collaborate to drive morphogenesis and provides a predictive model to decipher the underlying mechanism.

Embryonic development is a precise process involving cell division, cell-cell interaction, and cell migration. During the process, how each cell reaches its supposed location and be in contact with the right neighbors, and what roles genetic factors and physical forces play are important and fascinating questions. Using the worm Caenorhabditis elegans as a model system, we build a phase field model to simulate early morphogenesis. With a few physical inputs, the model can precisely reproduce the early morphological development of the worm. Such an accurate simulator can not only teach us how physical forces work together with genetic factors to shape up the complex process of development, but also make predictions, such as key cell-cell attractions critical in the process.

Transcriptional regulation of neuronal identity

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