Multidimensional regulation of gene expression in the C. elegans embryo

The Molecular Basis of Differentiation Wave Activity in Embryogenesis Bradly Alicea

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

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 mathematical framework for understanding the spontaneous emergence of complexity applicable to growing multicellular systems

In embryonic development and organogenesis, cells sharing identical genetic codes acquire diverse gene expression states in a highly reproducible spatial distribution, crucial for multicellular formation and quantifiable through positional information. To understand the spontaneous growth of complexity, we constructed a one-dimensional division-decision model, simulating the growth of cells with identical genetic networks from a single cell. Our findings highlight the pivotal role of cell division in providing positional cues, escorting the system toward states rich in information. Moreover, we pinpointed lateral inhibition as a critical mechanism translating spatial contacts into gene expression. Our model demonstrates that the spatial arrangement resulting from cell division, combined with cell lineages, imparts positional information, specifying multiple cell states with increased complexity-illustrated through examples in C.elegans. This study constitutes a foundational step in comprehending developmental intricacies, paving the way for future quantitative formulations to construct synthetic multicellular patterns.

Biclustering analysis on tree-shaped time-series single cell gene expression data of Caenorhabditis elegans

BACKGROUND

In recent years, gene clustering analysis has become a widely used tool for studying gene functions, efficiently categorizing genes with similar expression patterns to aid in identifying gene functions. Caenorhabditis elegans is commonly used in embryonic research due to its consistent cell lineage from fertilized egg to adulthood. Biologists use 4D confocal imaging to observe gene expression dynamics at the single-cell level. However, on one hand, the observed tree-shaped time-series datasets have characteristics such as non-pairwise data points between different individuals. On the other hand, the influence of cell type heterogeneity should also be considered during clustering, aiming to obtain more biologically significant clustering results.

RESULTS

A biclustering model is proposed for tree-shaped single-cell gene expression data of Caenorhabditis elegans. Detailedly, a tree-shaped piecewise polynomial function is first employed to fit non-pairwise gene expression time series data. Then, four factors are considered in the objective function, including Pearson correlation coefficients capturing gene correlations, p-values from the Kolmogorov-Smirnov test measuring the similarity between cells, as well as gene expression size and bicluster overlapping size. After that, Genetic Algorithm is utilized to optimize the function.

CONCLUSION

The results on the small-scale dataset analysis validate the feasibility and effectiveness of our model and are superior to existing classical biclustering models. Besides, gene enrichment analysis is employed to assess the results on the complete real dataset analysis, confirming that the discovered biclustering results hold significant biological relevance.

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.

LET-381/FoxF and its target UNC-30/Pitx2 specify and maintain the molecular identity of C. elegans mesodermal glia that regulate motor behavior

While most glial cell types in the central nervous system (CNS) arise from neuroectodermal progenitors, some, like microglia, are mesodermally derived. To understand mesodermal glia development and function, we investigated C. elegans GLR glia, which envelop the brain neuropil and separate it from the circulatory system cavity. Transcriptome analysis shows that GLR glia combine astrocytic and endothelial characteristics, which are relegated to separate cell types in vertebrates. Combined fate acquisition is orchestrated by LET-381/FoxF, a fate-specification/maintenance transcription factor also expressed in glia and endothelia of other animals. Among LET-381/FoxF targets, the UNC-30/Pitx2 transcription factor controls GLR glia morphology and represses alternative mesodermal fates. LET-381 and UNC-30 co-expression in naive cells is sufficient for GLR glia gene expression. GLR glia inactivation by ablation or let-381 mutation disrupts locomotory behavior and promotes salt-induced paralysis, suggesting brain-neuropil activity dysregulation. Our studies uncover mechanisms of mesodermal glia development and show that like neuronal differentiation, glia differentiation requires autoregulatory terminal selector genes that define and maintain the glial fate.

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.

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.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

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.

Mapping of centriolar proteins onto the post-embryonic lineage of C. elegans

Centrioles, together with the surrounding peri-centriolar material (PCM), constitute the centrosome, a major microtubule-organizing center of animal cells. Despite being critical in many cells for signaling, motility and division, centrioles can be eliminated in some systems, including in the vast majority of differentiating cells during embryogenesis in Caenorhabditis elegans. Whether the cells retaining centrioles in the resulting L1 larvae do so because they lack an activity that eliminates centrioles in the other cells is not known. Moreover, the extent to which centrioles and PCM remain present in later stages of worm development, when all cells but those of the germ line are terminally differentiated, is not known. Here, by fusing cells that lack centrioles with cells that retain them, we established that L1 larvae do not possess a diffusible elimination activity sufficient to remove centrioles. Moreover, analyzing PCM core proteins in L1 larval cells that retain centrioles, we found that some such proteins, but not all, are present as well. Furthermore, we uncovered that foci of centriolar proteins remain present in specific terminally differentiated cells of adult hermaphrodites and males, in particular in the somatic gonad. Correlating the time at which cells were born with the fate of their centrioles revealed that it is not cell age, but instead cell fate, that determines whether and when centrioles are eliminated. Overall, our work maps the localization of centriolar and PCM core proteins in the post-embryonic C. elegans lineage, thereby providing an essential blueprint for uncovering mechanisms modulating their presence and function.

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

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

Extensive programmed centriole elimination unveiled in C. elegans embryos

Centrioles are critical for fundamental cellular processes, including signaling, motility, and division. The extent to which centrioles are present after cell cycle exit in a developing organism is not known. The stereotypical lineage of Caenorhabditis elegans makes it uniquely well-suited to investigate this question. Using notably lattice light-sheet microscopy, correlative light electron microscopy, and lineage assignment, we found that ~88% of cells lose centrioles during embryogenesis. Our analysis reveals that centriole elimination is stereotyped, occurring invariably at a given time in a given cell type. Moreover, we established that experimentally altering cell fate results in corresponding changes in centriole fate. Overall, we uncovered the existence of an extensive centriole elimination program, which we anticipate to be paradigmatic for a broad understanding of centriole fate regulation.

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.

An atonal homolog, lin-32, regulates hypodermal morphogenesis in Caenorhabditis elegans

The transcription factor atonal contributes to patterning and cell fate determination in specialized epithelial cells in various animals, but its function in hypodermis is unknown. Here, we analyzed the atonal homolog lin-32 in C. elegans to clarify whether atonal acts in hypodermal development. The lin-32 null mutants exhibited bulges and cavities in their head, which were prevented by LIN-32 expression. Fluorescent protein was expressed in hypodermis cells at the embryonic stage by the lin-32 promoter. These results certify that atonal plays an essential role in the development of a broader range of tissues as hypodermis than initially thought.

Systematic analysis of cell morphodynamics in C. elegans early embryogenesis

The invariant cell lineage of Caenorhabditis elegans allows unambiguous assignment of the identity for each cell, which offers a unique opportunity to study developmental dynamics such as the timing of cell division, dynamics of gene expression, and cell fate decisions at single-cell resolution. However, little is known about cell morphodynamics, including the extent to which they are variable between individuals, mainly due to the lack of sufficient amount and quality of quantified data. In this study, we systematically quantified the cell morphodynamics in 52 C. elegans embryos from the two-cell stage to mid-gastrulation at the high spatiotemporal resolution, 0.5 μm thickness of optical sections, and 30-second intervals of recordings. Our data allowed systematic analyses of the morphological features. We analyzed sphericity dynamics and found a significant increase at the end of metaphase in every cell, indicating the universality of the mitotic cell rounding. Concomitant with the rounding, the volume also increased in most but not all cells, suggesting less universality of the mitotic swelling. Combining all features showed that cell morphodynamics was unique for each cell type. The cells before the onset of gastrulation could be distinguished from all the other cell types. Quantification of reproducibility in cell-cell contact revealed that variability in division timings and cell arrangements produced variability in contacts between the embryos. However, the area of such contacts occupied less than 5% of the total area, suggesting the high reproducibility of spatial occupancies and adjacency relationships of the cells. By comparing the morphodynamics of identical cells between the embryos, we observed diversity in the variability between cells and found it was determined by multiple factors, including cell lineage, cell generation, and cell-cell contact. We compared the variabilities of cell morphodynamics and cell-cell contacts with those in ascidian Phallusia mammillata embryos. The variabilities were larger in C. elegans, despite smaller differences in embryo size and number of cells at each developmental stage.

V-ATPase/TORC1-mediated ATFS-1 translation directs mitochondrial UPR activation in C. elegans

To adapt mitochondrial function to the ever-changing intra- and extracellular environment, multiple mitochondrial stress response (MSR) pathways, including the mitochondrial unfolded protein response (UPRmt), have evolved. However, how the mitochondrial stress signal is sensed and relayed to UPRmt transcription factors, such as ATFS-1 in Caenorhabditis elegans, remains largely unknown. Here, we show that a panel of vacuolar H+-ATPase (v-ATPase) subunits and the target of rapamycin complex 1 (TORC1) activity are essential for the cytosolic relay of mitochondrial stress to ATFS-1 and for the induction of the UPRmt. Mechanistically, mitochondrial stress stimulates v-ATPase/Rheb-dependent TORC1 activation, subsequently promoting ATFS-1 translation. Increased translation of ATFS-1 upon mitochondrial stress furthermore relies on a set of ribosomal components but is independent of GCN-2/PEK-1 signaling. Finally, the v-ATPase and ribosomal subunits are required for mitochondrial surveillance and mitochondrial stress-induced longevity. These results reveal a v-ATPase-TORC1-ATFS-1 signaling pathway that links mitochondrial stress to the UPRmt through intimate crosstalks between multiple organelles.

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.

A metabolic regulatory network for the Caenorhabditis elegans intestine

Metabolic perturbations can affect gene expression, for instance to rewire metabolism. While numerous efforts have measured gene expression in response to individual metabolic perturbations, methods that determine all metabolic perturbations that affect the expression for a given gene or set of genes have not been available. Here, we use a gene-centered approach to derive a first-pass metabolic regulatory network for Caenorhabditis elegans by performing RNAi of more than 1,400 metabolic genes with a set of 19 promoter reporter strains that express a fluorescent protein in the animal's intestine. We find that metabolic perturbations generally increase promoter activity, which contrasts with transcription factor (TF) RNAi, which tends to repress promoter activity. We identify several TFs that modulate promoter activity in response to perturbations of the electron transport chain and explore complex genetic interactions among metabolic pathways. This work provides a blueprint for a systems-level understanding of how metabolism affects gene expression.

Delineating mechanisms and design principles of Caenorhabditis elegans embryogenesis using in toto high-resolution imaging data and computational modeling

The nematode (roundworm) Caenorhabditis elegans is one of the most popular animal models for the study of developmental biology, as its invariant development and transparent body enable in toto cellular-resolution fluorescence microscopy imaging of developmental processes at 1-min intervals. This has led to the development of various computational tools for the systematic and automated analysis of imaging data to delineate the molecular and cellular processes throughout the embryogenesis of C. elegans, such as those associated with cell lineage, cell migration, cell morphology, and gene activity. In this review, we first introduce C. elegans embryogenesis and the development of techniques for tracking cell lineage and reconstructing cell morphology during this process. We then contrast the developmental modes of C. elegans and the customized technologies used for studying them with the ones of other animal models, highlighting its advantage for studying embryogenesis with exceptional spatial and temporal resolution. This is followed by an examination of the physical models that have been devised—based on accurate determinations of developmental processes afforded by analyses of imaging data—to interpret the early embryonic development of C. elegans from subcellular to intercellular levels of multiple cells, which focus on two key processes: cell polarization and morphogenesis. We subsequently discuss how quantitative data-based theoretical modeling has improved our understanding of the mechanisms of C. elegans embryogenesis. We conclude by summarizing the challenges associated with the acquisition of C. elegans embryogenesis data, the construction of algorithms to analyze them, and the theoretical interpretation.

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.

mNG-tagged mls-2 knock-in alleles in C. elegans

The Caenorhabditis elegans HMX/NKX MLS-2 transcription factor was previously shown to play sequential roles in AWC general identity and the stochastic AWC^ON/AWC^OFF subtype choice during embryogenesis. Here we analyze the expression pattern of endogenous mls-2 during AWC development using mNeonGreen (mNG) knock-in strains. Similar to transgenic GFP::MLS-2, functional mNG::MLS-2 knock-in displayed nuclear localization in AWC precursor cells but was not observed in AWC during the later embryonic stage. These results suggest that the expression of mls-2 is below the detectable level and/or the stability of MLS-2 protein is very low in AWC cells.

The transcription factor unc-130/FOXD3/4 contributes to the biphasic calcium response required to optimize avoidance behavior

The central neural network optimizes avoidance behavior depending on the nociceptive stimulation intensity and is essential for survival. How the property of hub neurons that enables the selection of behaviors is genetically defined is not well understood. We show that the transcription factor unc-130, a human FOXD3/4 ortholog, is required to optimize avoidance behavior depending on stimulus strength in Caenorhabditis elegans. unc-130 is necessary for both ON responses (calcium decreases) and OFF responses (calcium increases) in AIBs, central neurons of avoidance optimization. Ablation of predicted upstream inhibitory neurons reduces the frequency of turn behavior, suggesting that optimization needs both calcium responses. At the molecular level, unc-130 upregulates the expression of at least three genes: nca-2, a homolog of the vertebrate cation leak channel NALCN; glr-1, an AMPA-type glutamate receptor; and eat-4, a hypothetical L-glutamate transmembrane transporter in the central neurons of optimization. unc-130 shows more limited regulation in optimizing behavior than an atonal homolog lin-32, and unc-130 and lin-32 appear to act in parallel molecular pathways. Our findings suggest that unc-130 is required for the establishment of some AIB identities to optimize avoidance behavior.

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.

Lineage-specific control of convergent differentiation by a Forkhead repressor

During convergent differentiation, multiple developmental lineages produce a highly similar or identical cell type. However, few molecular players that drive convergent differentiation are known. Here, we show that the C. elegans Forkhead transcription factor UNC-130 is required in only one of three convergent lineages that produce the same glial cell type. UNC-130 acts transiently as a repressor in progenitors and newly-born terminal cells to allow the proper specification of cells related by lineage rather than by cell type or function. Specification defects correlate with UNC-130:DNA binding, and UNC-130 can be functionally replaced by its human homolog, the neural crest lineage determinant FoxD3. We propose that, in contrast to terminal selectors that activate cell type-specific transcriptional programs in terminally differentiating cells, UNC-130 acts early and specifically in one convergent lineage to produce a cell type that also arises from molecularly distinct progenitors in other lineages.

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

Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state-fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.

Piecemeal regulation of convergent neuronal lineages by bHLH transcription factors in Caenorhabditis elegans

Cells of the same type can be generated by distinct cellular lineages that originate in different parts of the developing embryo ('lineage convergence'). Several Caenorhabditis elegans neuron classes composed of left/right or radially symmetric class members display such lineage convergence. We show here that the C. elegans Atonal homolog lin-32 is differentially expressed in neuronal lineages that give rise to left/right or radially symmetric class members. Loss of lin-32 results in the selective loss of the expression of pan-neuronal markers and terminal selector-type transcription factors that confer neuron class-specific features. Another basic helix-loop-helix (bHLH) gene, the Achaete-Scute homolog hlh-14, is expressed in a mirror image pattern relative to lin-32 and is required to induce neuronal identity and terminal selector expression on the contralateral side of the animal. These findings demonstrate that distinct lineage histories converge via different bHLH factors at the level of induction of terminal selector identity determinants, which thus serve as integrators of distinct lineage histories. We also describe neuron-to-neuron identity transformations in lin-32 mutants, which we propose to also be the result of misregulation of terminal selector gene expression.

Single-cell dynamics of chromatin activity during cell lineage differentiation in Caenorhabditis elegans embryos

Elucidating the chromatin dynamics that orchestrate embryogenesis is a fundamental question in developmental biology. Here, we exploit position effects on expression as an indicator of chromatin activity and infer the chromatin activity landscape in every lineaged cell during Caenorhabditis elegans early embryogenesis. Systems-level analyses reveal that chromatin activity distinguishes cellular states and correlates with fate patterning in the early embryos. As cell lineage unfolds, chromatin activity diversifies in a lineage-dependent manner, with switch-like changes accompanying anterior-posterior fate asymmetry and characteristic landscapes being established in different cell lineages. Upon tissue differentiation, cellular chromatin from distinct lineages converges according to tissue types but retains stable memories of lineage history, contributing to intra-tissue cell heterogeneity. However, the chromatin landscapes of cells organized in a left-right symmetric pattern are predetermined to be analogous in early progenitors so as to pre-set equivalent states. Finally, genome-wide analysis identifies many regions exhibiting concordant chromatin activity changes that mediate the co-regulation of functionally related genes during differentiation. Collectively, our study reveals the developmental and genomic dynamics of chromatin activity at the single-cell level.

Physically asymmetric division of the C. elegans zygote ensures invariably successful embryogenesis

Asymmetric divisions that yield daughter cells of different sizes are frequent during early embryogenesis, but the importance of such a physical difference for successful development remains poorly understood. Here, we investigated this question using the first division of Caenorhabditis elegans embryos, which yields a large AB cell and a small P1 cell. We equalized AB and P1 sizes using acute genetic inactivation or optogenetic manipulation of the spindle positioning protein LIN-5. We uncovered that only some embryos tolerated equalization, and that there was a size asymmetry threshold for viability. Cell lineage analysis of equalized embryos revealed an array of defects, including faster cell cycle progression in P1 descendants, as well as defects in cell positioning, division orientation, and cell fate. Moreover, equalized embryos were more susceptible to external compression. Overall, we conclude that unequal first cleavage is essential for invariably successful embryonic development of C. elegans.

Meta-analysis of Caenorhabditis elegans single-cell developmental data reveals multi-frequency oscillation in gene activation

MOTIVATION

The advent of in vivo automated techniques for single-cell lineaging, sequencing and analysis of gene expression has begun to dramatically increase our understanding of organismal development. We applied novel meta-analysis and visualization techniques to the EPIC single-cell-resolution developmental gene expression dataset for Caenorhabditis elegans from Bao, Murray, Waterston et al. to gain insights into regulatory mechanisms governing the timing of development.

RESULTS

Our meta-analysis of the EPIC dataset revealed that a simple linear combination of the expression levels of the developmental genes is strongly correlated with the developmental age of the organism, irrespective of the cell division rate of different cell lineages. We uncovered a pattern of collective sinusoidal oscillation in gene activation, in multiple dominant frequencies and in multiple orthogonal axes of gene expression, pointing to the existence of a coordinated, multi-frequency global timing mechanism. We developed a novel method based on Fisher's Discriminant Analysis to identify gene expression weightings that maximally separate traits of interest, and found that remarkably, simple linear gene expression weightings are capable of producing sinusoidal oscillations of any frequency and phase, adding to the growing body of evidence that oscillatory mechanisms likely play an important role in the timing of development. We cross-linked EPIC with gene ontology and anatomy ontology terms, employing Fisher's Discriminant Analysis methods to identify previously unknown positive and negative genetic contributions to developmental processes and cell phenotypes. This meta-analysis demonstrates new evidence for direct linear and/or sinusoidal mechanisms regulating the timing of development. We uncovered a number of previously unknown positive and negative correlations between developmental genes and developmental processes or cell phenotypes. Our results highlight both the continued relevance of the EPIC technique, and the value of meta-analysis of previously published results. The presented analysis and visualization techniques are broadly applicable across developmental and systems biology.

AVAILABILITY AND IMPLEMENTATION

Analysis software available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Establishment of a morphological atlas of the Caenorhabditis elegans embryo using deep-learning-based 4D segmentation

The invariant development and transparent body of the nematode Caenorhabditis elegans enables complete delineation of cell lineages throughout development. Despite extensive studies of cell division, cell migration and cell fate differentiation, cell morphology during development has not yet been systematically characterized in any metazoan, including C. elegans. This knowledge gap substantially hampers many studies in both developmental and cell biology. Here we report an automatic pipeline, CShaper, which combines automated segmentation of fluorescently labeled membranes with automated cell lineage tracing. We apply this pipeline to quantify morphological parameters of densely packed cells in 17 developing C. elegans embryos. Consequently, we generate a time-lapse 3D atlas of cell morphology for the C. elegans embryo from the 4- to 350-cell stages, including cell shape, volume, surface area, migration, nucleus position and cell-cell contact with resolved cell identities. We anticipate that CShaper and the morphological atlas will stimulate and enhance further studies in the fields of developmental biology, cell biology and biomechanics.

The systematic characterization of C. elegans morphology during development has yet to be performed. Here, the authors produce a 3D atlas of C. elegans morphology from 17 embryos and 54 developmental stages, using an automated pipeline, CShaper (combining segmentation of fluorescently labeled membranes with automated cell lineage tracing).

Single-cell RNA profiling links ncRNAs to spatiotemporal gene expression during C. elegans embryogenesis

Recent studies show that non-coding RNAs (ncRNAs) can regulate the expression of protein-coding genes and play important roles in mammalian development. Previous studies have revealed that during C. elegans (Caenorhabditis elegans) embryo development, numerous genes in each cell are spatiotemporally regulated, causing the cell to differentiate into distinct cell types and tissues. We ask whether ncRNAs participate in the spatiotemporal regulation of genes in different types of cells and tissues during the embryogenesis of C. elegans. Here, by using marker-free full-length high-depth single-cell RNA sequencing (scRNA-seq) technique, we sequence the whole transcriptomes from 1031 embryonic cells of C. elegans and detect 20,431 protein-coding genes, including 22 cell-type-specific protein-coding markers, and 9843 ncRNAs including 11 cell-type-specific ncRNA markers. We induce a ncRNAs-based clustering strategy as a complementary strategy to the protein-coding gene-based clustering strategy for single-cell classification. We identify 94 ncRNAs that have never been reported to regulate gene expressions, are co-expressed with 1208 protein-coding genes in cell type specific and/or embryo time specific manners. Our findings suggest that these ncRNAs could potentially influence the spatiotemporal expression of the corresponding genes during the embryogenesis of C. elegans.

cnd-1/NeuroD1 Functions with the Homeobox Gene ceh-5/Vax2 and Hox Gene ceh-13/labial To Specify Aspects of RME and DD Neuron Fate in Caenorhabditis elegans

Identifying the mechanisms behind neuronal fate specification are key to understanding normal neural development in addition to neurodevelopmental disorders such as autism and schizophrenia. In vivo cell fate specification is difficult to study in vertebrates. However, the nematode Caenorhabditis elegans, with its invariant cell lineage and simple nervous system of 302 neurons, is an ideal organism to explore the earliest stages of neural development. We used a comparative transcriptome approach to examine the role of cnd-1/NeuroD1 in C. elegans nervous system development and function. This basic helix-loop-helix transcription factor is deeply conserved across phyla and plays a crucial role in cell fate specification in both the vertebrate nervous system and pancreas. We find that cnd-1 controls expression of ceh-5, a Vax2-like homeobox class transcription factor, in the RME head motorneurons and PVQ tail interneurons. We also show that cnd-1 functions redundantly with the Hox gene ceh-13/labial in defining the fate of DD1 and DD2 embryonic ventral nerve cord motorneurons. These data highlight the utility of comparative transcriptomes for identifying transcription factor targets and understanding gene regulatory networks.

Regulation of Gliogenesis by lin-32/Atoh1 in Caenorhabditis elegans

The regulation of gliogenesis is a fundamental process for nervous system development, as the appropriate glial number and identity is required for a functional nervous system. To investigate the molecular mechanisms involved in gliogenesis, we used C. elegans as a model and identified the function of the proneural gene lin-32/Atoh1 in gliogenesis. We found that lin-32 functions during embryonic development to negatively regulate the number of AMsh glia. The ectopic AMsh cells at least partially arise from cells originally fated to become CEPsh glia, suggesting that lin-32 is involved in the specification of specific glial subtypes. Moreover, we show that lin-32 acts in parallel with cnd-1/ NeuroD1 and ngn-1/ Neurog1 in negatively regulating an AMsh glia fate. Furthermore, expression of murine Atoh1 fully rescues lin-32 mutant phenotypes, suggesting lin-32/Atoh1 may have a conserved role in glial specification.

Multilevel regulation of muscle-specific transcription factor hlh-1 during Caenorhabditis elegans embryogenesis

hlh-1 is a myogenic transcription factor required for body-wall muscle specification during embryogenesis in Caenorhabditis elegans. Despite its well-known role in muscle specification, comprehensive regulatory control upstream of hlh-1 remains poorly defined. Here, we first established a statistical reference for the spatiotemporal expression of hlh-1 at single-cell resolution up to the second last round of divisions for most of the cell lineages (from 4- to 350-cell stage) using 13 wild-type embryos. We next generated lineal expression of hlh-1 after RNA interference (RNAi) perturbation of 65 genes, which were selected based on their degree of conservation, mutant phenotypes, and known roles in development. We then compared the expression profiles between wild-type and RNAi embryos by clustering according to their lineal expression patterns using mean-shift and density-based clustering algorithms, which not only confirmed the roles of existing genes but also uncovered the potential functions of novel genes in muscle specification at multiple levels, including cellular, lineal, and embryonic levels. By combining the public data on protein-protein interactions, protein-DNA interactions, and genetic interactions with our RNAi data, we inferred regulatory pathways upstream of hlh-1 that function globally or locally. This work not only revealed diverse and multilevel regulatory mechanisms coordinating muscle differentiation during C. elegans embryogenesis but also laid a foundation for further characterizing the regulatory pathways controlling muscle specification at the cellular, lineal (local), or embryonic (global) level.

Alignment of Cell Lineage Trees Elucidates Genetic Programs for the Development and Evolution of Cell Types

A full understanding of the developmental process requires fine-scale characterization of cell divisions and cell types, which are naturally organized as the developmental cell lineage tree (CLT). Technological breakthroughs facilitated determination of more CLTs, but complete comprehension of the data remains difficult without quantitative comparison among CLTs. We hereby quantified phenotypic similarity between CLTs using a novel computational method that exhaustively searches for optimal correspondence between individual cells meanwhile retaining their topological relationships. The revealed CLT similarities allowed us to infer functional similarity at the transcriptome level, identify cell fate transformations, predict functional relationships between mutants, and find evolutionary correspondence between cell types of different species. By allowing quantitative comparison between CLTs, our work is expected to greatly enhance the interpretability of relevant data and help answer the myriad of questions surrounding the developmental process.

•

Align cell lineage trees (CLTs) to search/quantify their phenotypic similarities

•

Aligning worm CLTs captured known genetic/developmental programs

•

Similarities between knockdown CLTs revealed functional relationships between genes

•

CLT alignments between species gave insight on the evolution of cell types

ngn-1/neurogenin Activates Transcription of Multiple Terminal Selector Transcription Factors in the Caenorhabditis elegans Nervous System

Proper nervous system development is required for an organism's survival and function. Defects in neurogenesis have been linked to neurodevelopmental disorders such as schizophrenia and autism. Understanding the gene regulatory networks that orchestrate neural development, specifically cascades of proneural transcription factors, can better elucidate which genes are most important during early neurogenesis. Neurogenins are a family of deeply conserved factors shown to be both necessary and sufficient for the development of neural subtypes. However, the immediate downstream targets of neurogenin are not well characterized. The objective of this study was to further elucidate the role of ngn-1/neurogenin in nervous system development and to identify its downstream transcriptional targets, using the nematode Caenorhabditis elegans as a model for this work. We found that ngn-1 is required for axon outgrowth, nerve ring architecture, and neuronal cell fate specification. We also showed that ngn-1 may have roles in neuroblast migration and epithelial integrity during embryonic development. Using RNA sequencing and comparative transcriptome analysis, we identified eight transcription factors (hlh-34/NPAS1, unc-42/PROP1, ceh-17/PHOX2A, lim-4/LHX6, fax-1/NR2E3, lin-11/LHX1, tlp-1/ZNF503, and nhr-23/RORB) whose transcription is activated, either directly or indirectly, by ngn-1 Our results show that ngn-1 has a role in transcribing known terminal regulators that establish and maintain cell fate of differentiated neural subtypes and confirms that ngn-1 functions as a proneural transcription factor in C. elegans neurogenesis.

Inferring Causal Gene Regulatory Networks from Coupled Single-Cell Expression Dynamics Using Scribe

Here, we present Scribe (https://github.com/aristoteleo/Scribe-py), a toolkit for detecting and visualizing causal regulatory interactions between genes and explore the potential for single-cell experiments to power network reconstruction. Scribe employs restricted directed information to determine causality by estimating the strength of information transferred from a potential regulator to its downstream target. We apply Scribe and other leading approaches for causal network reconstruction to several types of single-cell measurements and show that there is a dramatic drop in performance for "pseudotime"-ordered single-cell data compared with true time-series data. We demonstrate that performing causal inference requires temporal coupling between measurements. We show that methods such as "RNA velocity" restore some degree of coupling through an analysis of chromaffin cell fate commitment. These analyses highlight a shortcoming in experimental and computational methods for analyzing gene regulation at single-cell resolution and suggest ways of overcoming it.

Dendrites with specialized glial attachments develop by retrograde extension using SAX-7 and GRDN-1

Dendrites develop elaborate morphologies in concert with surrounding glia, but the molecules that coordinate dendrite and glial morphogenesis are mostly unknown. C. elegans offers a powerful model for identifying such factors. Previous work in this system examined dendrites and glia that develop within epithelia, similar to mammalian sense organs. Here, we focus on the neurons BAG and URX, which are not part of an epithelium but instead form membranous attachments to a single glial cell at the nose, reminiscent of dendrite-glia contacts in the mammalian brain. We show that these dendrites develop by retrograde extension, in which the nascent dendrite endings anchor to the presumptive nose and then extend by stretching during embryo elongation. Using forward genetic screens, we find that dendrite development requires the adhesion protein SAX-7/L1CAM and the cytoplasmic protein GRDN-1/CCDC88C to anchor dendrite endings at the nose. SAX-7 acts in neurons and glia, while GRDN-1 acts in glia to non-autonomously promote dendrite extension. Thus, this work shows how glial factors can help to shape dendrites, and identifies a novel molecular mechanism for dendrite growth by retrograde extension.

Deep reinforcement learning of cell movement in the early stage of C.elegans embryogenesis

Motivation

Cell movement in the early phase of Caenorhabditis elegans development is regulated by a highly complex process in which a set of rules and connections are formulated at distinct scales. Previous efforts have demonstrated that agent-based, multi-scale modeling systems can integrate physical and biological rules and provide new avenues to study developmental systems. However, the application of these systems to model cell movement is still challenging and requires a comprehensive understanding of regulatory networks at the right scales. Recent developments in deep learning and reinforcement learning provide an unprecedented opportunity to explore cell movement using 3D time-lapse microscopy images.

Results

We present a deep reinforcement learning approach within an agent-based modeling system to characterize cell movement in the embryonic development of C.elegans. Our modeling system captures the complexity of cell movement patterns in the embryo and overcomes the local optimization problem encountered by traditional rule-based, agent-based modeling that uses greedy algorithms. We tested our model with two real developmental processes: the anterior movement of the Cpaaa cell via intercalation and the rearrangement of the superficial left-right asymmetry. In the first case, the model results suggested that Cpaaa's intercalation is an active directional cell movement caused by the continuous effects from a longer distance (farther than the length of two adjacent cells), as opposed to a passive movement caused by neighbor cell movements. In the second case, a leader-follower mechanism well explained the collective cell movement pattern in the asymmetry rearrangement. These results showed that our approach to introduce deep reinforcement learning into agent-based modeling can test regulatory mechanisms by exploring cell migration paths in a reverse engineering perspective. This model opens new doors to explore the large datasets generated by live imaging.

Availability and implementation

Source code is available at https://github.com/zwang84/drl4cellmovement.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

A Myt1 family transcription factor defines neuronal fate by repressing non-neuronal genes

Cellular differentiation requires both activation of target cell transcriptional programs and repression of non-target cell programs. The Myt1 family of zinc finger transcription factors contributes to fibroblast to neuron reprogramming in vitro. Here, we show that ztf-11 (Zinc-finger Transcription Factor-11), the sole Caenorhabditis elegans Myt1 homolog, is required for neurogenesis in multiple neuronal lineages from previously differentiated epithelial cells, including a neuron generated by a developmental epithelial-to-neuronal transdifferentiation event. ztf-11 is exclusively expressed in all neuronal precursors with remarkable specificity at single-cell resolution. Loss of ztf-11 leads to upregulation of non-neuronal genes and reduced neurogenesis. Ectopic expression of ztf-11 in epidermal lineages is sufficient to produce additional neurons. ZTF-11 functions together with the MuvB corepressor complex to suppress the activation of non-neuronal genes in neurons. These results dovetail with the ability of Myt1l (Myt1-like) to drive neuronal transdifferentiation in vitro in vertebrate systems. Together, we identified an evolutionarily conserved mechanism to specify neuronal cell fate by repressing non-neuronal genes.

Digging deeper: methodologies for high-content phenotyping in Caenorhabditis elegans

Deep phenotyping is an emerging conceptual paradigm and experimental approach aimed at measuring and linking many aspects of a phenotype to understand its underlying biology. To date, deep phenotyping has been applied mostly in cultured cells and used less in multicellular organisms. However, in the past decade, it has increasingly been recognized that deep phenotyping could lead to a better understanding of how genetics, environment and stochasticity affect the development, physiology and behavior of an organism. The nematode Caenorhabditis elegans is an invaluable model system for studying how genes affect a phenotypic trait, and new technologies have taken advantage of the worm's physical attributes to increase the throughput and informational content of experiments. Coupling of these technical advancements with computational and analytical tools has enabled a boom in deep-phenotyping studies of C. elegans. In this Review, we highlight how these new technologies and tools are digging into the biological origins of complex, multidimensional phenotypes.

The C. elegans embryonic transcriptome with tissue, time, and alternative splicing resolution

We have used RNA-seq in Caenorhabditis elegans to produce transcription profiles for seven specific embryonic cell populations from gastrulation to the onset of terminal differentiation. The expression data for these seven cell populations, covering major cell lineages and tissues in the worm, reveal the complex and dynamic changes in gene expression, both spatially and temporally. Also, within genes, start sites and exon usage can be highly differential, producing transcripts that are specific to developmental periods or cell lineages. We have also found evidence of novel exons and introns, as well as differential usage of SL1 and SL2 splice leaders. By combining this data set with the modERN ChIP-seq resource, we are able to support and predict gene regulatory relationships. The detailed information on differences and similarities between gene expression in cell lineages and tissues should be of great value to the community and provides a framework for the investigation of expression in individual cells.

Progression of Meiosis Is Coordinated by the Level and Location of MAPK Activation Via OGR-2 in Caenorhabditis elegans

During meiosis, a series of evolutionarily conserved events allow for reductional chromosome division, which is required for sexual reproduction. Although individual meiotic processes have been extensively studied, we currently know far less about how meiosis is regulated and coordinated. In the Caenorhabditis elegans gonad, mitogen-activated protein kinase (MAPK) signaling drives oogenesis while undergoing spatial activation and deactivation waves. However, it is currently unclear how MAPK activation is governed and how it facilitates the progression of oogenesis. Here, we show that the oocyte and germline-related 2 (ogr-2) gene affects proper progression of oogenesis. Complete deletion of ogr-2 results in delayed meiotic entry and late spatial onset of double-strand break repair. Elevated levels of apoptosis are observed in this mutant, independent of the meiotic canonical checkpoints; however, they are dependent on the MAPK terminal member MPK-1/ERK. MPK-1 activation is elevated in diplotene in ogr-2 mutants and its aberrant spatial activation correlates with stages where meiotic progression defects are evident. Deletion of ogr-2 significantly reduces the expression of lip-1, a phosphatase reported to repress MPK-1, which is consistent with OGR-2 localization at chromatin in germ cells. We suggest that OGR-2 modulates the expression of lip-1 to promote the timely progression of meiosis through MPK-1 spatial deactivation.

MRG-1 is required for both chromatin-based transcriptional silencing and genomic integrity of primordial germ cells in Caenorhabditis elegans

In Caenorhabditis elegans, germline cells remain transcriptionally silenced during embryogenesis. The transcriptional silencing is achieved by two different mechanisms: One is the inhibition of RNA polymerase II in P2-P4 cells at the establishment stage, and another is chromatin-based silencing in two primordial germ cells (PGCs) at the maintenance stage; however, the molecular mechanism underlying chromatin-based silencing is less understood. We investigated the role of the chromodomain protein MRG-1, which is an essential maternal factor for germline development, in transcriptional silencing in PGCs. PGCs lacking maternal MRG-1 showed increased levels of two histone modifications (H3K4me2 and H4K16ac), which are epigenetic markers for active transcription, and precocious activation of germline promoters. Loss of MES-4, a H3K36 methyltransferase, also caused similar derepression of the germline genes in PGCs, suggesting that both MRG-1 and MES-4 function in chromatin-based silencing in PGCs. In addition, the mrg-1 null mutant showed abnormal chromosome structures and a decrease in homologous recombinase RAD-51 foci in PGCs, but the mes-4 null mutant did not show such phenotypes. Taken together, we propose that MRG-1 has two distinct functions: chromatin-based transcriptional silencing and preserving genomic integrity at the maintenance stage of PGCs.