John I. Murray

Regulatory analysis of the C. elegans genome with spatiotemporal resolution

Discovering the structure and dynamics of transcriptional regulatory events in the genome with cellular and temporal resolution is crucial to understanding the regulatory underpinnings of development and disease. We determined the genomic distribution of binding sites for 92 transcription factors and regulatory proteins across multiple stages of Caenorhabditis elegans development by performing 241 ChIP-seq (chromatin immunoprecipitation followed by sequencing) experiments. Integration of regulatory binding and cellular-resolution expression data produced a spatiotemporally resolved metazoan transcription factor binding map. Using this map, we explore developmental regulatory circuits that encode combinatorial logic at the levels of co-binding and co-expression of transcription factors, characterizing the genomic coverage and clustering of regulatory binding, the binding preferences of, and biological processes regulated by, transcription factors, the global transcription factor co-associations and genomic subdomains that suggest shared patterns of regulation, and identifying key transcription factors and transcription factor co-associations for fate specification of individual lineages and cell types.

Notch and Ras promote sequential steps of excretory tube development in C. elegans

Receptor tyrosine kinases and Notch are crucial for tube formation and branching morphogenesis in many systems, but the specific cellular processes that require signaling are poorly understood. Here we describe sequential roles for Notch and Epidermal growth factor (EGF)-Ras-ERK signaling in the development of epithelial tube cells in the C. elegans excretory (renal-like) organ. This simple organ consists of three tandemly connected unicellular tubes: the excretory canal cell, duct and G1 pore. lin-12 and glp-1/Notch are required to generate the canal cell, which is a source of LIN-3/EGF ligand and physically attaches to the duct during de novo epithelialization and tubulogenesis. Canal cell asymmetry and let-60/Ras signaling influence which of two equivalent precursors will attach to the canal cell. Ras then specifies duct identity, inducing auto-fusion and a permanent epithelial character; the remaining precursor becomes the G1 pore, which eventually loses epithelial character and withdraws from the organ to become a neuroblast. Ras continues to promote subsequent aspects of duct morphogenesis and differentiation, and acts primarily through Raf-ERK and the transcriptional effectors LIN-1/Ets and EOR-1. These results reveal multiple genetically separable roles for Ras signaling in tube development, as well as similarities to Ras-mediated control of branching morphogenesis in more complex organs, including the mammalian kidney. The relative simplicity of the excretory system makes it an attractive model for addressing basic questions about how cells gain or lose epithelial character and organize into tubular networks.

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.

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.

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.

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.

SAVER: gene expression recovery for single-cell RNA sequencing

In single-cell RNA sequencing (scRNA-seq) studies, only a small fraction of the transcripts present in each cell are sequenced. This leads to unreliable quantification of genes with low or moderate expression, which hinders downstream analysis. To address this challenge, we developed SAVER (single-cell analysis via expression recovery), an expression recovery method for unique molecule index (UMI)-based scRNA-seq data that borrows information across genes and cells to provide accurate expression estimates for all genes.SAVER accurately recovers expression values in single-cell RNA-sequencing data to improve downstream analysis.

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.

The Nkx5/HMX homeodomain protein MLS-2 is required for proper tube cell shape in the C. elegans excretory system

Cells perform wide varieties of functions that are facilitated, in part, by adopting unique shapes. Many of the genes and pathways that promote cell fate specification have been elucidated. However, relatively few transcription factors have been identified that promote shape acquisition after fate specification. Here we show that the Nkx5/HMX homeodomain protein MLS-2 is required for cellular elongation and shape maintenance of two tubular epithelial cells in the C. elegans excretory system, the duct and pore cells. The Nkx5/HMX family is highly conserved from sea urchins to humans, with known roles in neuronal and glial development. MLS-2 is expressed in the duct and pore, and defects in mls-2 mutants first arise when the duct and pore normally adopt unique shapes. MLS-2 cooperates with the EGF-Ras-ERK pathway to turn on the LIN-48/Ovo transcription factor in the duct cell during morphogenesis. These results reveal a novel interaction between the Nkx5/HMX family and the EGF-Ras pathway and implicate a transcription factor, MLS-2, as a regulator of cell shape.

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. genesis, 2016.