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

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

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.

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.

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.

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.

H3K9me2 orchestrates inheritance of spatial positioning of peripheral heterochromatin through mitosis.

Cell-type-specific 3D organization of the genome is unrecognizable during mitosis. It remains unclear how essential positional information is transmitted through cell division such that a daughter cell recapitulates the spatial genome organization of the parent. Lamina-associated domains (LADs) are regions of repressive heterochromatin positioned at the nuclear periphery that vary by cell type and contribute to cell-specific gene expression and identity. Here we show that histone 3 lysine 9 dimethylation (H3K9me2) is an evolutionarily conserved, specific mark of nuclear peripheral heterochromatin and that it is retained through mitosis. During mitosis, phosphorylation of histone 3 serine 10 temporarily shields the H3K9me2 mark allowing for dissociation of chromatin from the nuclear lamina. Using high-resolution 3D immuno-oligoFISH, we demonstrate that H3K9me2-enriched genomic regions, which are positioned at the nuclear lamina in interphase cells prior to mitosis, re-associate with the forming nuclear lamina before mitotic exit. The H3K9me2 modification of peripheral heterochromatin ensures that positional information is safeguarded through cell division such that individual LADs are re-established at the nuclear periphery in daughter nuclei. Thus, H3K9me2 acts as a 3D architectural mitotic guidepost. Our data establish a mechanism for epigenetic memory and inheritance of spatial organization of the genome.

Silencing the alternative.

The transcription factor ztf-11 promotes neuronal differentiation by repressing other cell fates in the nematode worm C. elegans.

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.

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.

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.

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

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.

Deconvolution of gene expression from cell populations across the C. elegans lineage.

BACKGROUND: Knowledge of when and in which cells each gene is expressed across multicellular organisms is critical in understanding both gene function and regulation of cell type diversity. However, methods for measuring expression typically involve a trade-off between imaging-based methods, which give the precise location of a limited number of genes, and higher throughput methods such as RNA-seq, which include all genes, but are more limited in their resolution to apply to many tissues. We propose an intermediate method, which estimates expression in individual cells, based on high-throughput measurements of expression from multiple overlapping groups of cells. This approach has particular benefits in organisms such as C. elegans where invariant developmental patterns make it possible to define these overlapping populations of cells at single-cell resolution. RESULT: We implement several methods to deconvolve the gene expression in individual cells from population-level data and determine the accuracy of these estimates on simulated data from the C. elegans embryo. CONCLUSION: These simulations suggest that a high-resolution map of expression in the C. elegans embryo may be possible with expression data from as few as 30 cell populations.

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

The invariant lineage of Caenorhabditis elegans has powerful potential for quantifying developmental variability in normal and stressed embryos. Previous studies of division timing by automated lineage tracing suggested that variability in cell cycle timing is low in younger embryos, but manual lineage tracing of specific lineages suggested that variability may increase for later divisions. We developed improved automated lineage tracing methods that allowroutine lineage tracing through the last round of embryonic cell divisions and we applied these methods to trace the lineage of 18 wild-type embryos. Cell cycle lengths, division axes and cell positions are remarkably consistent among these embryos at all stages, with only slight increase in variability later in development. The resulting quantitative 4-dimensional model of embryogenesis provides a powerful reference dataset to identify defects in mutants or in embryos that have experienced environmental perturbations. We also traced the lineages of embryos imaged at higher temperatures to quantify the decay in developmental robustness under temperature stress. Developmental variability increases modestly at 25°C compared with 22°C and dramatically at 26°C, and we identify homeotic transformations in a subset of embryos grown at 26°C. The deep lineage tracing methods provide a powerful tool for analysis of normal development, gene expression and mutants and we provide a graphical user interface to allow other researchers to explore the average behavior of arbitrary cells in a reference embryo.

Automated lineage and expression profiling in live Caenorhabditis elegans embryos.

Describing gene expression during animal development requires a way to quantitatively measure expression levels with cellular resolution and to describe how expression changes with time. Fluorescent protein reporters make it possible to measure expression dynamics in live cells by time-lapse microscopy, but it can be challenging to identify expressing cells in complex tissues and to compare expression across organisms. This protocol describes how to use automated lineage analysis to identify cells in Caenorhabditis elegans embryos expressing fluorescent reporters and how to quantify that expression with cellular resolution. Because C. elegans develops through an invariant pattern of cell divisions, every cell's identity and future fate can be predicted from its pattern of previous cell divisions. Automated analysis of images collected from embryos expressing a fluorescent histone transgene in all cells allows lineage tracing and cell identification. This provides a scaffold with which to describe expression of a second color reporter such as a fusion of a second fluorescent protein to a gene of interest or its regulatory sequences. These methods can also be used for analysis of reporter expression, cell division timing, and cell position in genetically perturbed embryos. The protocol describes how to prepare C. elegans strains containing nuclear-expressed fluorescent reporters, collect images of appropriate quality from embryos, perform automated lineage analysis, manually edit and curate the lineage, and, finally, extract and display reporter signals.

Multidimensional regulation of gene expression in the C. elegans embryo.

How cells adopt different expression patterns is a fundamental question of developmental biology. We quantitatively measured reporter expression of 127 genes, primarily transcription factors, in every cell and with high temporal resolution in C. elegans embryos. Embryonic cells are highly distinct in their gene expression; expression of the 127 genes studied here can distinguish nearly all pairs of cells, even between cells of the same tissue type. We observed recurrent lineage-regulated expression patterns for many genes in diverse contexts. These patterns are regulated in part by the TCF-LEF transcription factor POP-1. Other genes' reporters exhibited patterns correlated with tissue, position, and left-right asymmetry. Sequential patterns both within tissues and series of sublineages suggest regulatory pathways. Expression patterns often differ between embryonic and larval stages for the same genes, emphasizing the importance of profiling expression in different stages. This work greatly expands the number of genes in each of these categories and provides the first large-scale, digitally based, cellular resolution compendium of gene expression dynamics in live animals. The resulting data sets will be a useful resource for future research.

Automated analysis of embryonic gene expression with cellular resolution in C. elegans.

We describe a system that permits the automated analysis of reporter gene expression in Caenorhabditis elegans with cellular resolution continuously during embryogenesis. We demonstrate its utility by defining the expression patterns of reporters for several embryonically expressed transcription factors. The invariant cell lineage permits the automated alignment of multiple expression profiles, allowing direct comparison of the expression of different genes' reporters. We also used this system to monitor perturbations to normal development involving changes both in cell-division timing and in cell fate. Systematic application of this system could reveal the gene activity of each cell throughout development.

The lineaging of fluorescently-labeled Caenorhabditis elegans embryos with StarryNite and AceTree.

Lineage analysis of Caenorhabditis elegans is a powerful tool for characterizing developmental phenotypes and embryonic gene-expression patterns. We present a detailed protocol for the lineaging of embryos by computational analysis of 4D images of embryos that ubiquitously express histone-GFP (green fluorescent protein) fusion proteins through the 350 cell stage followed by manual editing. We describe how to optimize imaging settings for this purpose, the use of the lineage-extraction software, StarryNite, and the lineage-editing software, AceTree. In addition, we describe a useful polymer bead mounting technique for C. elegans embryos that has several advantages compared with the standard agar pad mounting technique. The protocol requires about 1 h of user time spread over 2 days to generate the raw lineage, and an additional 2 or 4 h to edit the lineage to the 194- or 350-cell stage, respectively.

Diverse and specific gene expression responses to stresses in cultured human cells.

We used cDNA microarrays in a systematic study of the gene expression responses of HeLa cells and primary human lung fibroblasts to heat shock, endoplasmic reticulum stress, oxidative stress, and crowding. Hierarchical clustering of the data revealed groups of genes with coherent biological themes, including genes that responded to specific stresses and others that responded to multiple types of stress. Fewer genes increased in expression after multiple stresses than in free-living yeasts, which have a large general stress response program. Most of the genes induced by multiple diverse stresses are involved in cell-cell communication and other processes specific to higher organisms. We found substantial differences between the stress responses of HeLa cells and primary fibroblasts. For example, many genes were induced by oxidative stress and dithiothreitol in fibroblasts but not HeLa cells; conversely, a group of transcription factors, including c-fos and c-jun, were induced by heat shock in HeLa cells but not in fibroblasts. The dataset is freely available for search and download at http://microarray-pubs.stanford.edu/human_stress/Home.shtml.