A quantitative model of normal Caenorhabditis elegans embryogenesis and its disruption after stress

Nuclear Structure, Size Regulation, and Role in Cell Migration

The nucleus serves as a pivotal regulatory and control hub in the cell, governing numerous aspects of cellular functions, including DNA replication, transcription, and RNA processing. Therefore, any deviations in nuclear morphology, structure, or organization can strongly affect cellular activities. In this review, we provide an updated perspective on the structure and function of nuclear components, focusing on the linker of nucleoskeleton and cytoskeleton complex, the nuclear envelope, the nuclear lamina, and chromatin. Additionally, nuclear size should be considered a fundamental parameter for the cellular state. Its regulation is tightly linked to environmental changes, development, and various diseases, including cancer. Hence, we also provide a concise overview of different mechanisms by which nuclear size is determined, the emerging role of the nucleus as a mechanical sensor, and the implications of altered nuclear morphology on the physiology of diseased cells.

Thermosensory behaviors of the free-living life stages of Strongyloides species support parasitism in tropical environments

Soil-transmitted parasitic nematodes infect over 1 billion people worldwide and are a common source of neglected disease. Strongyloides stercoralis is a potentially fatal skin-penetrating human parasite that is endemic to tropical and subtropical regions around the world. The complex life cycle of Strongyloides species is unique among human-parasitic nematodes in that it includes a single free-living generation featuring soil-dwelling, bacterivorous adults whose progeny all develop into infective larvae. The sensory behaviors that enable free-living Strongyloides adults to navigate and survive soil environments are unknown. S. stercoralis infective larvae display parasite-specific sensory-driven behaviors, including robust attraction to mammalian body heat. In contrast, the free-living model nematode Caenorhabditis elegans displays thermosensory behaviors that guide adult worms to stay within a physiologically permissive range of environmental temperatures. Do S. stercoralis and C. elegans free-living adults, which experience similar environmental stressors, display common thermal preferences? Here, we characterize the thermosensory behaviors of the free-living adults of S. stercoralis as well as those of the closely related rat parasite, Strongyloides ratti. We find that Strongyloides free-living adults are exclusively attracted to near-tropical temperatures, despite their inability to infect mammalian hosts. We further show that lifespan is shorter at higher temperatures for free-living Strongyloides adults, similar to the effect of temperature on C. elegans lifespan. However, we also find that the reproductive potential of the free-living life stage is enhanced at warmer temperatures, particularly for S. stercoralis. Together, our results reveal a novel role for thermotaxis to maximize the infectious capacity of obligate parasites and provide insight into the biological adaptations that may contribute to their endemicity in tropical climates.

Thermosensory behaviors of the free-living life stages of Strongyloides species support parasitism in tropical environments

Soil-transmitted parasitic nematodes infect over 1 billion people worldwide and are a common source of neglected disease. Strongyloides stercoralis is a potentially fatal skin-penetrating human parasite that is endemic to tropical and subtropical regions around the world. The complex life cycle of Strongyloides species is unique among human-parasitic nematodes in that it includes a single free-living generation featuring soil-dwelling, bacterivorous adults whose progeny all develop into infective larvae. The sensory behaviors that enable free-living Strongyloides adults to navigate and survive soil environments are unknown. S. stercoralis infective larvae display parasite-specific sensory-driven behaviors, including robust attraction to mammalian body heat. In contrast, the free-living model nematode Caenorhabditis elegans displays thermosensory behaviors that guide adult worms to stay within a physiologically permissive range of environmental temperatures. Do S. stercoralis and C. elegans free-living adults, which experience similar environmental stressors, display common thermal preferences? Here, we characterize the thermosensory behaviors of the free-living adults of S. stercoralis as well as those of the closely related rat parasite, Strongyloides ratti. We find that Strongyloides free-living adults are exclusively attracted to near-tropical temperatures, despite their inability to infect mammalian hosts. We further show that lifespan is shorter at higher temperatures for free-living Strongyloides adults, similar to the effect of temperature on C. elegans lifespan. However, we also find that the reproductive potential of the free-living life stage is enhanced at warmer temperatures, particularly for S. stercoralis. Together, our results reveal a novel role for thermotaxis to maximize the infectious capacity of obligate parasites and provide insight into the biological adaptations that may contribute to their endemicity in tropical climates.

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.

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.

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.

Delineating the 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.

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.

Single cell biology-a Keystone Symposia report

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

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.

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.

Loss of histone H3.3 results in DNA replication defects and altered origin dynamics in C. elegans

Histone H3.3 is a replication-independent variant of histone H3 with important roles in development, differentiation, and fertility. Here, we show that loss of H3.3 results in replication defects in Caenorhabditis elegans embryos at elevated temperatures. To characterize these defects, we adapt methods to determine replication timing, map replication origins, and examine replication fork progression. Our analysis of the spatiotemporal regulation of DNA replication shows that despite the very rapid embryonic cell cycle, the genome is replicated from early and late firing origins and is partitioned into domains of early and late replication. We find that under temperature stress conditions, additional replication origins become activated. Moreover, loss of H3.3 results in altered replication fork progression around origins, which is particularly evident at stress-activated origins. These replication defects are accompanied by replication checkpoint activation, a delayed cell cycle, and increased lethality in checkpoint-compromised embryos. Our comprehensive analysis of DNA replication in C. elegans reveals the genomic location of replication origins and the dynamics of their firing, and uncovers a role of H3.3 in the regulation of replication origins under stress conditions.

PhenoMIP: High-Throughput Phenotyping of Diverse Caenorhabditis elegans Populations via Molecular Inversion Probes

Whether generated within a lab setting or isolated from the wild, variant alleles continue to be an important resource for decoding gene function in model organisms such as Caenorhabditis elegans. With advances in massively parallel sequencing, multiple whole-genome sequenced (WGS) strain collections are now available to the research community. The Million Mutation Project (MMP) for instance, analyzed 2007 N2-derived, mutagenized strains. Individually, each strain averages ∼400 single nucleotide variants amounting to ∼80 protein-coding variants. The effects of these variants, however, remain largely uncharacterized and querying the breadth of these strains for phenotypic changes requires a method amenable to rapid and sensitive high-throughput analysis. Here we present a pooled competitive fitness approach to quantitatively phenotype subpopulations of sequenced collections via molecular inversion probes (PhenoMIP). We phenotyped the relative fitness of 217 mutant strains on multiple food sources and classified these into five categories. We also demonstrate on a subset of these strains, that their fitness defects can be genetically mapped. Overall, our results suggest that approximately 80% of MMP mutant strains may have a decreased fitness relative to the lab reference, N2. The costs of generating this form of analysis through WGS methods would be prohibitive while PhenoMIP analysis in this manner is accomplished at less than one-tenth of projected WGS costs. We propose methods for applying PhenoMIP to a broad range of population selection experiments in a cost-efficient manner that would be useful to the community at large.

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.

A portable, low-cost device for precise control of specimen temperature under stereomicroscopes

To facilitate precise and convenient control of biological sample temperature, we developed a low-cost device that can be used independently or with any stereomicroscope. The purpose of the device is to control the thermal environment during experimental intervals in which a specimen must be manipulated outside of an incubator, e.g. for dissection or slide-mounting in preparation for imaging. Sample temperatures can be both cooled to below and heated to above room temperatures, and stably maintained at a precision of +/- 0.1˚C. To demonstrate the utility of this device, we report improved characterization of the penetrance of a short-acting temperature-sensitive allele in C. elegans embryos, and identification of the upper temperature threshold for embryonic viability for six Caenorhabditis species. By controlling the temperature environment even as a specimen is manipulated, this device offers consistency and flexibility, reduces environmental noise, and enables precision timing in experiments requiring temperature shifts.

Systems Properties and Spatiotemporal Regulation of Cell Position Variability during Embryogenesis

An intriguing question in developmental biology is how do developmental processes achieve high reproducibility among individuals? An in-depth analysis of information contained in phenotypic variability provides an important perspective to address this question. In this work, we present a quantitative and functional analysis of cell position variability during Caenorhabditis elegans embryogenesis. We find that cell position variability is highly deterministic and regulated by intrinsic and extrinsic mechanisms. Positional variability is determined by cell lineage identity and is coupled to diverse developmental properties of cells, including embryonic localization, cell contact, and left-right symmetry. Temporal dynamics of cell position variability are highly concordant, and fate specification contributes to a systems-wide reduction of variability that could provide a buffering strategy. Positional variability is stringently regulated throughout embryogenesis and cell-cell junctions function to restrict variability. Our results provide insight into systems properties and spatiotemporal control of cellular variability during development.

Embryonic canalization and its limits-A view from temperature

Many animals are able to produce similar offspring over a range of environmental conditions. This property of the developmental process has been termed canalization-the channeling of developmental pathways to generate a stable outcome despite varying conditions. Temperature is one environmental parameter that has fundamental effects on cell physiology and biochemistry, yet developmental programs generally result in a stable phenotype under a range of temperatures. On the other hand, there are typically upper and lower temperature limits beyond which the developmental program is unable to produce normal offspring. This review summarizes data on how development is affected by temperature, particularly high temperature, in various animal species. It also brings together information on potential cell biological and developmental genetic factors that may be responsible for developmental stability in varying temperatures, and likely critical mechanisms that break down at high temperature. Also reviewed are possible means for studying temperature effects on embryogenesis and how to determine which factors are most critical at the high-temperature limits for normal development. Increased knowledge of these critical factors will point to the targets of selection under climate change, and more generally, how developmental robustness in varying environments is maintained.

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.

Sample path properties of the average generation of a Bellman-Harris process

Motivated by a recently proposed design for a DNA coded randomised algorithm that enables inference of the average generation of a collection of cells descendent from a common progenitor, here we establish strong convergence properties for the average generation of a super-critical Bellman-Harris process. We further extend those results to a two-type Bellman-Harris process where one type can give rise to the other, but not vice versa. These results further affirm the estimation method's potential utility by establishing its long run accuracy on individual sample-paths, and significantly expanding its remit to encompass cellular development that gives rise to differentiated offspring with distinct population dynamics.

Establishment of Signaling Interactions with Cellular Resolution for Every Cell Cycle of Embryogenesis

Intercellular signaling interactions play a key role in breaking fate symmetry during animal development. Identification of signaling interactions at cellular resolution is technically challenging, especially in a developing embryo. Here, we develop a platform that allows automated inference and validation of signaling interactions for every cell cycle of Caenorhabditis elegans embryogenesis. This is achieved by the generation of a systems-level cell contact map, which consists of 1114 highly confident intercellular contacts, by modeling analysis and is validated through cell membrane labeling coupled with cell lineage analysis. We apply the map to identify cell pairs between which a Notch signaling interaction takes place. By generating expression patterns for two ligands and two receptors of the Notch signaling pathway with cellular resolution using the automated expression profiling technique, we are able to refine existing and identify novel Notch interactions during C. elegans embryogenesis. Targeted cell ablation followed by cell lineage analysis demonstrates the roles of signaling interactions during cell division in breaking fate symmetry. Finally, we describe the development of a website that allows online access to the cell-cell contact map for mapping of other signaling interactions by the community. The platform can be adapted to establish cellular interactions from any other signaling pathway.

Automated C. elegans embryo alignments reveal brain neuropil position invariance despite lax cell body placement

The Caenorhabditis elegans cell lineage is nearly invariant. Whether this stereotyped cell-division pattern promotes reproducibility in cell shapes/positions is not generally known, as manual spatiotemporal cell-shape/position alignments are labor-intensive, and fully-automated methods are not described. Here, we report automated algorithms for spatiotemporal alignments of C. elegans embryos from pre-morphogenesis to motor-activity initiation. We use sparsely-labeled green-fluorescent nuclei and a pan-nuclear red-fluorescent reporter to register consecutive imaging time points and compare embryos. Using our method, we monitor early assembly of the nerve-ring (NR) brain neuropil. While NR pioneer neurons exhibit reproducible growth kinetics and axon positions, cell-body placements are variable. Thus, pioneer-neuron axon locations, but not cell-body positions, are under selection. We also show that anterior NR displacement in cam-1/ROR Wnt-receptor mutants is not an early NR assembly defect. Our results demonstrate the utility of automated spatiotemporal alignments of C. elegans embryos, and uncover key principles guiding nervous-system development in this animal.

Systems biology of embryonic development: Prospects for a complete understanding of the Caenorhabditis elegans embryo

The convergence of developmental biology and modern genomics tools brings the potential for a comprehensive understanding of developmental systems. This is especially true for the Caenorhabditis elegans embryo because its small size, invariant developmental lineage, and powerful genetic and genomic tools provide the prospect of a cellular resolution understanding of messenger RNA (mRNA) expression and regulation across the organism. We describe here how a systems biology framework might allow large-scale determination of the embryonic regulatory relationships encoded in the C. elegans genome. This framework consists of two broad steps: (a) defining the "parts list"-all genes expressed in all cells at each time during development and (b) iterative steps of computational modeling and refinement of these models by experimental perturbation. Substantial progress has been made towards defining the parts list through imaging methods such as large-scale green fluorescent protein (GFP) reporter analysis. Imaging results are now being augmented by high-resolution transcriptome methods such as single-cell RNA sequencing, and it is likely the complete expression patterns of all genes across the embryo will be known within the next few years. In contrast, the modeling and perturbation experiments performed so far have focused largely on individual cell types or genes, and improved methods will be needed to expand them to the full genome and organism. This emerging comprehensive map of embryonic expression and regulatory function will provide a powerful resource for developmental biologists, and would also allow scientists to ask questions not accessible without a comprehensive picture. This article is categorized under: Invertebrate Organogenesis > Worms Technologies > Analysis of the Transcriptome Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

Comparative proteome analysis between C . briggsae embryos and larvae reveals a role of chromatin modification proteins in embryonic cell division

Caenorhabditis briggsae has emerged as a model for comparative biology against model organism C. elegans. Most of its cell fate specifications are completed during embryogenesis whereas its cell growth is achieved mainly in larval stages. The molecular mechanism underlying the drastic developmental changes is poorly understood. To gain insights into the molecular changes between the two stages, we compared the proteomes between the two stages using iTRAQ. We identified a total of 2,791 proteins in the C. briggsae embryos and larvae, 247 of which undergo up- or down-regulation between the two stages. The proteins that are upregulated in the larval stages are enriched in the Gene Ontology categories of energy production, protein translation, and cytoskeleton; whereas those upregulated in the embryonic stage are enriched in the categories of chromatin dynamics and posttranslational modification, suggesting a more active chromatin modification in the embryos than in the larva. Perturbation of a subset of chromatin modifiers followed by cell lineage analysis suggests their roles in controlling cell division pace. Taken together, we demonstrate a general molecular switch from chromatin modification to metabolism during the transition from C. briggsae embryonic to its larval stages using iTRAQ approach. The switch might be conserved across metazoans.

Biologically constrained optimization based cell membrane segmentation in C. elegans embryos

BACKGROUND

Recent advances in bioimaging and automated analysis methods have enabled the large-scale systematic analysis of cellular dynamics during the embryonic development of Caenorhabditis elegans. Most of these analyses have focused on cell lineage tracing rather than cell shape dynamics. Cell shape analysis requires cell membrane segmentation, which is challenging because of insufficient resolution and image quality. This problem is currently solved by complicated segmentation methods requiring laborious and time consuming parameter adjustments.

RESULTS

Our new framework BCOMS (Biologically Constrained Optimization based cell Membrane Segmentation) automates the extraction of the cell shape of C. elegans embryos. Both the segmentation and evaluation processes are automated. To automate the evaluation, we solve an optimization problem under biological constraints. The performance of BCOMS was validated against a manually created ground truth of the 24-cell stage embryo. The average deviation of 25 cell shape features was 5.6%. The deviation was mainly caused by membranes parallel to the focal planes, which either contact the surfaces of adjacent cells or make no contact with other cells. Because segmentation of these membranes was difficult even by manual inspection, the automated segmentation was sufficiently accurate for cell shape analysis. As the number of manually created ground truths is necessarily limited, we compared the segmentation results between two adjacent time points. Across all cells and all cell cycles, the average deviation of the 25 cell shape features was 4.3%, smaller than that between the automated segmentation result and ground truth.

CONCLUSIONS

BCOMS automated the accurate extraction of cell shapes in developing C. elegans embryos. By replacing image processing parameters with easily adjustable biological constraints, BCOMS provides a user-friendly framework. The framework is also applicable to other model organisms. Creating the biological constraints is a critical step requiring collaboration between an experimentalist and a software developer.

Infrared laser-induced gene expression for tracking development and function of single C. elegans embryonic neurons

Visualizing neural-circuit assembly in vivo requires tracking growth of optically resolvable neurites. The Caenorhabditis elegans embryonic nervous system, comprising 222 neurons and 56 glia, is attractive for comprehensive studies of development; however, embryonic reporters are broadly expressed, making single-neurite tracking/manipulation challenging. We present a method, using an infrared laser, for reproducible heat-dependent gene expression in small sublineages (one to four cells) without radiation damage. We go beyond proof-of-principle, and use our system to label and track single neurons during early nervous-system assembly. We uncover a retrograde extension mechanism for axon growth, and reveal the aetiology of axon-guidance defects in sax-3/Robo and vab-1/EphR mutants. We also perform cell-specific rescues, determining DAF-6/patched-related site of action during sensory-organ development. Simultaneous ablation and labelling of cells using our system reveals roles for glia in dendrite extension. Our method can be applied to other optically/IR-transparent organisms, and opens the door to high-resolution systematic analyses of C. elegans morphogenesis.

Gut development in C. elegans

The midgut (intestine) of the nematode, C. elegans, is a tube consisting of 20 cells that arises from a single embryonic precursor. Owing to its comparatively simple anatomy and the advantages inherent to the C. elegans system, the gut has been used as a model for organogenesis for more than 25 years. In this review, the salient features of C. elegans gut development are described from the E progenitor through to the 20-cell intestine. The core gene regulatory network that drives specification of the gut, and other genes with roles in organogenesis, lumen morphogenesis and the cell cycle, are also described. Questions for future work are posed.

Site-specific recombinatorics: in situ cellular barcoding with the Cre Lox system

BACKGROUND

Cellular barcoding is a recently developed biotechnology tool that enables the familial identification of progeny of individual cells in vivo. In immunology, it has been used to track the burst-sizes of multiple distinct responding T cells over several adaptive immune responses. In the study of hematopoiesis, it revealed fate heterogeneity amongst phenotypically identical multipotent cells. Most existing approaches rely on ex vivo viral transduction of cells with barcodes followed by adoptive transfer into an animal, which works well for some systems, but precludes barcoding cells in their native environment such as those inside solid tissues.

RESULTS

With a view to overcoming this limitation, we propose a new design for a genetic barcoding construct based on the Cre Lox system that induces randomly created stable barcodes in cells in situ by exploiting inherent sequence distance constraints during site-specific recombination. We identify the cassette whose provably maximal code diversity is several orders of magnitude higher than what is attainable with previously considered Cre Lox barcoding approaches, exceeding the number of lymphocytes or hematopoietic progenitor cells in mice.

CONCLUSIONS

Its high diversity and in situ applicability, make the proposed Cre Lox based tagging system suitable for whole tissue or even whole animal barcoding. Moreover, it can be built using established technology.

Combinatorial decoding of the invariant C. elegans embryonic lineage in space and time

Understanding how a single cell, the zygote, can divide and differentiate to produce the diverse animal cell types is a central goal of developmental biology research. The model organism Caenorhabditis elegans provides a system that enables a truly comprehensive understanding of this process across all cells. Its invariant cell lineage makes it possible to identify all of the cells in each individual and compare them across organisms. Recently developed methods automate the process of cell identification, allowing high-throughput gene expression characterization and phenotyping at single cell resolution. In this Review, we summarize the sequences of events that pattern the lineage including establishment of founder cell identity, the signaling pathways that diversify embryonic fate, and the regulators involved in patterning within these founder lineages before cells adopt their terminal fates. We focus on insights that have emerged from automated approaches to lineage tracking, including insights into mechanisms of robustness, context-specific regulation of gene expression, and temporal coordination of differentiation. We suggest a model by which lineage history produces a combinatorial code of transcription factors that act, often redundantly, to ensure terminal fate.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Inferring average generation via division-linked labeling

For proliferating cells subject to both division and death, how can one estimate the average generation number of the living population without continuous observation or a division-diluting dye? In this paper we provide a method for cell systems such that at each division there is an unlikely, heritable one-way label change that has no impact other than to serve as a distinguishing marker. If the probability of label change per cell generation can be determined and the proportion of labeled cells at a given time point can be measured, we establish that the average generation number of living cells can be estimated. Crucially, the estimator does not depend on knowledge of the statistics of cell cycle, death rates or total cell numbers. We explore the estimator's features through comparison with physiologically parameterized stochastic simulations and extrapolations from published data, using it to suggest new experimental designs.

Experience Modulates the Reproductive Response to Heat Stress in C. elegans via Multiple Physiological Processes

Natural environments are considerably more variable than laboratory settings and often involve transient exposure to stressful conditions. To fully understand how organisms have evolved to respond to any given stress, prior experience must therefore be considered. We investigated the effects of individual and ancestral experience on C. elegans reproduction. We documented ways in which cultivation at 15°C or 25°C affects developmental time, lifetime fecundity, and reproductive performance after severe heat stress that exceeds the fertile range of the organism but is compatible with survival and future fecundity. We found that experience modulates multiple aspects of reproductive physiology, including the male and female germ lines and the interaction between them. These responses vary in their environmental sensitivity, suggesting the existence of complex mechanisms for coping with unpredictable and stressful environments.

Quantitative Differences in Nuclear β-catenin and TCF Pattern Embryonic Cells in C. elegans

The Wnt signaling pathway plays a conserved role during animal development in transcriptional regulation of distinct targets in different developmental contexts but it remains unclear whether quantitative differences in the nuclear localization of effector proteins TCF and β-catenin contribute to context-specific regulation. We investigated this question in Caenorhabditis elegans embryos by quantifying nuclear localization of fluorescently tagged SYS-1/β-catenin and POP-1/TCF and expression of Wnt ligands at cellular resolution by time-lapse microscopy and automated lineage tracing. We identified reproducible, quantitative differences that generate a subset of Wnt-signaled cells with a significantly higher nuclear concentration of the TCF/β-catenin activating complex. Specifically, β-catenin and TCF are preferentially enriched in nuclei of daughter cells whose parents also had high nuclear levels of that protein, a pattern that could influence developmental gene expression. Consistent with this, we found that expression of synthetic reporters of POP-1-dependent activation is biased towards cells that had high nuclear SYS-1 in consecutive divisions. We identified new genes whose embryonic expression patterns depend on pop-1. Most of these require POP-1 for either transcriptional activation or repression, and targets requiring POP-1 for activation are more likely to be expressed in the cells with high nuclear SYS-1 in consecutive divisions than those requiring POP-1 for repression. Taken together, these results indicate that SYS-1 and POP-1 levels are influenced by the parent cell’s SYS-1/POP-1 levels and this may provide an additional mechanism by which POP-1 regulates distinct targets in different developmental contexts.

The Wnt signaling pathway is active during the development of all multi-cellular animals and also improperly re-activated in many cancers. Here, we use time-lapse microscopy to quantify the nuclear localization of several proteins in response to Wnt signaling throughout early embryonic development in the nematode worm, C. elegans. We find that cells that received a Wnt signal in the previous division respond more strongly to a Wnt signal in the next division, in part by localizing more of the regulator β-catenin to the nucleus. This causes the relative enrichment of Wnt pathway proteins in the nuclei of repeatedly signaled cells, which we show likely impacts the activation of Wnt target genes. This represents a novel mechanism for the regulation of Wnt pathway targets in development and disease.

Systems-level quantification of division timing reveals a common genetic architecture controlling asynchrony and fate asymmetry

Coordination of cell division timing is crucial for proper cell fate specification and tissue growth. However, the differential regulation of cell division timing across or within cell types during metazoan development remains poorly understood. To elucidate the systems-level genetic architecture coordinating division timing, we performed a high-content screening for genes whose depletion produced a significant reduction in the asynchrony of division between sister cells (ADS) compared to that of wild-type during Caenorhabditis elegans embryogenesis. We quantified division timing using 3D time-lapse imaging followed by computer-aided lineage analysis. A total of 822 genes were selected for perturbation based on their conservation and known roles in development. Surprisingly, we find that cell fate determinants are not only essential for establishing fate asymmetry, but also are imperative for setting the ADS regardless of cellular context, indicating a common genetic architecture used by both cellular processes. The fate determinants demonstrate either coupled or separate regulation between the two processes. The temporal coordination appears to facilitate cell migration during fate specification or tissue growth. Our quantitative dataset with cellular resolution provides a resource for future analyses of the genetic control of spatial and temporal coordination during metazoan development.

Cerium oxide nanoparticle aggregates affect stress response and function in Caenorhabditis elegans

OBJECTIVE

The continual increase in production and disposal of nanomaterials raises concerns regarding the safety of nanoparticles on the environmental and human health. Recent studies suggest that cerium oxide (CeO2) nanoparticles may possess both harmful and beneficial effects on biological processes. The primary objective of this study is to evaluate how exposure to different concentrations (0.17-17.21 µg/mL) of aggregated CeO2 nanoparticles affects indices of whole animal stress and survivability in Caenorhabditis elegans.

METHODS

Caenorhabditis elegans were exposed to different concentrations of CeO2 nanoparticles and evaluated.

RESULTS

Our findings demonstrate that chronic exposure of CeO2 nanoparticle aggregates is associated with increased levels of reactive oxygen species and heat shock stress response (HSP-4) in Caenorhabditis elegans, but not mortality. Conversely, CeO2 aggregates promoted strain-dependent decreases in animal fertility, a decline in stress resistance as measured by thermotolerance, and shortened worm length.

CONCLUSION

The data obtained from this study reveal the sublethal toxic effects of CeO2 nanoparticle aggregates in Caenorhabditis elegans and contribute to our understanding of how exposure to CeO2 may affect the environment.

The Bicoid class homeodomain factors ceh-36/OTX and unc-30/PITX cooperate in C. elegans embryonic progenitor cells to regulate robust development

While many transcriptional regulators of pluripotent and terminally differentiated states have been identified, regulation of intermediate progenitor states is less well understood. Previous high throughput cellular resolution expression studies identified dozens of transcription factors with lineage-specific expression patterns in C. elegans embryos that could regulate progenitor identity. In this study we identified a broad embryonic role for the C. elegans OTX transcription factor ceh-36, which was previously shown to be required for the terminal specification of four neurons. ceh-36 is expressed in progenitors of over 30% of embryonic cells, yet is not required for embryonic viability. Quantitative phenotyping by computational analysis of time-lapse movies of ceh-36 mutant embryos identified cell cycle or cell migration defects in over 100 of these cells, but most defects were low-penetrance, suggesting redundancy. Expression of ceh-36 partially overlaps with that of the PITX transcription factor unc-30. unc-30 single mutants are viable but loss of both ceh-36 and unc-30 causes 100% lethality, and double mutants have significantly higher frequencies of cellular developmental defects in the cells where their expression normally overlaps. These factors are also required for robust expression of the downstream developmental regulator mls-2/HMX. This work provides the first example of genetic redundancy between the related yet evolutionarily distant OTX and PITX families of bicoid class homeodomain factors and demonstrates the power of quantitative developmental phenotyping in C. elegans to identify developmental regulators acting in progenitor cells.

Comprehensive single cell-resolution analysis of the role of chromatin regulators in early C. elegans embryogenesis

Chromatin regulators are widely expressed proteins with diverse roles in gene expression, nuclear organization, cell cycle regulation, pluripotency, physiology and development, and are frequently mutated in human diseases such as cancer. Their inhibition often results in pleiotropic effects that are difficult to study using conventional approaches. We have developed a semi-automated nuclear tracking algorithm to quantify the divisions, movements and positions of all nuclei during the early development of Caenorhabditis elegans and have used it to systematically study the effects of inhibiting chromatin regulators. The resulting high dimensional datasets revealed that inhibition of multiple regulators, including F55A3.3 (encoding FACT subunit SUPT16H), lin-53 (RBBP4/7), rba-1 (RBBP4/7), set-16 (MLL2/3), hda-1 (HDAC1/2), swsn-7 (ARID2), and let-526 (ARID1A/1B) affected cell cycle progression and caused chromosome segregation defects. In contrast, inhibition of cir-1 (CIR1) accelerated cell division timing in specific cells of the AB lineage. The inhibition of RNA polymerase II also accelerated these division timings, suggesting that normal gene expression is required to delay cell cycle progression in multiple lineages in the early embryo. Quantitative analyses of the dataset suggested the existence of at least two functionally distinct SWI/SNF chromatin remodeling complex activities in the early embryo, and identified a redundant requirement for the egl-27 and lin-40 MTA orthologs in the development of endoderm and mesoderm lineages. Moreover, our dataset also revealed a characteristic rearrangement of chromatin to the nuclear periphery upon the inhibition of multiple general regulators of gene expression. Our systematic, comprehensive and quantitative datasets illustrate the power of single cell-resolution quantitative tracking and high dimensional phenotyping to investigate gene function. Furthermore, the results provide an overview of the functions of essential chromatin regulators during the early development of an animal.

A novel cell nuclei segmentation method for 3D C. elegans embryonic time-lapse images

BACKGROUND

Recently a series of algorithms have been developed, providing automatic tools for tracing C. elegans embryonic cell lineage. In these algorithms, 3D images collected from a confocal laser scanning microscope were processed, the output of which is cell lineage with cell division history and cell positions with time. However, current image segmentation algorithms suffer from high error rate especially after 350-cell stage because of low signal-noise ratio as well as low resolution along the Z axis (0.5-1 microns). As a result, correction of the errors becomes a huge burden. These errors are mainly produced in the segmentation of nuclei. Thus development of a more accurate image segmentation algorithm will alleviate the hurdle for automated analysis of cell lineage.

RESULTS

This paper presents a new type of nuclei segmentation method embracing an bi-directional prediction procedure, which can greatly reduce the number of false negative errors, the most common errors in the previous segmentation. In this method, we first use a 2D region growing technique together with the level-set method to generate accurate 2D slices. Then a modified gradient method instead of the existing 3D local maximum method is adopted to detect all the 2D slices located in the nuclei center, each of which corresponds to one nucleus. Finally, the bi-directional prediction method based on the images before and after the current time point is introduced into the system to predict the nuclei in low quality parts of the images. The result of our method shows a notable improvement in the accuracy rate. For each nucleus, its precise location, volume and gene expression value (gray value) is also obtained, all of which will be useful in further downstream analyses.

CONCLUSIONS

The result of this research demonstrates the advantages of the bi-directional prediction method in the nuclei segmentation over that of StarryNite/MatLab StarryNite. Several other modifications adopted in our nuclei segmentation system are also discussed.

Efficient single-cell transgene induction in Caenorhabditis elegans using a pulsed infrared laser

The coupling of transgenes to heat shock promoters is a widely applied method for regulating gene expression. In C. elegans, gene induction can be controlled temporally through timing of heat shock and spatially via targeted rescue in heat shock mutants. Here, we present a method for evoking gene expression in arbitrary cells, with single-cell resolution. We use a focused pulsed infrared laser to locally induce a heat shock response in specific cells. Our method builds on and extends a previously reported method using a continuous-wave laser. In our technique, the pulsed laser illumination enables a much higher degree of spatial selectivity because of diffusion of heat between pulses. We apply our method to induce transient and long-term transgene expression in embryonic, larval, and adult cells. Our method allows highly selective spatiotemporal control of transgene expression and is a powerful tool for model organism biology.

Gene transcription is coordinated with, but not dependent on, cell divisions during C. elegans embryonic fate specification

Cell differentiation and proliferation are coordinated during animal development, but the link between them remains uncharacterized. To examine this relationship, we combined single-molecule RNA imaging with time-lapse microscopy to generate high-resolution measurements of transcriptional dynamics in Caenorhabditis elegans embryogenesis. We found that globally slowing the overall development rate of the embryo by altering temperature or by mutation resulted in cell proliferation and transcription slowing, but maintaining, their relative timings, suggesting that cell division may directly control transcription. However, using mutants with specific defects in cell cycle pathways that lead to abnormal lineages, we found that the order between cell divisions and expression onset can switch, showing that expression of developmental regulators is not strictly dependent on cell division. Delaying cell divisions resulted in only slight changes in absolute expression time, suggesting that expression and proliferation are independently entrained to a separate clock-like process. These changes in relative timing can change the number of cells expressing a gene at a given time, suggesting that timing may help determine which cells adopt particular transcriptional patterns. Our results place limits on the types of mechanisms that are used during normal development to ensure that division timing and fate specification occur at appropriate times.