AceTree: a tool for visual analysis of Caenorhabditis elegans embryogenesis

Deep learning-based enhancement of fluorescence labeling for accurate cell lineage tracing during embryogenesis

MOTIVATION

Automated cell lineage tracing throughout embryogenesis plays a key role in the study of regulatory control of cell fate differentiation, morphogenesis and organogenesis in the development of animals, including nematode Caenorhabditis elegans. However, automated cell lineage tracing suffers from an exponential increase in errors at late embryo because of the dense distribution of cells, relatively low signal-to-noise ratio (SNR) and imbalanced intensity profiles of fluorescence images, which demands a huge amount of human effort to manually correct the errors. The existing image enhancement methods are not sensitive enough to deal with the challenges posed by the crowdedness and low signal-to-noise ratio. An alternative method is urgently needed to assist the existing detection methods in improving their detection and tracing accuracy, thereby reducing the huge burden for manual curation.

RESULTS

We developed a new method, termed as DELICATE, that dramatically improves the accuracy of automated cell lineage tracing especially during the stage post 350 cells of C. elegans embryo. DELICATE works by increasing the local SNR and improving the evenness of nuclei fluorescence intensity across cells especially in the late embryos. The method both dramatically reduces the segmentation errors by StarryNite and the time required for manually correcting tracing errors up to 550-cell stage, allowing the generation of accurate cell lineage at large-scale with a user-friendly software/interface.

AVAILABILITY AND IMPLEMENTATION

All images and data are available at https://doi.org/10.6084/m9.figshare.26778475.v1. The code and user-friendly software are available at https://github.com/plcx/NucApp-develop.

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

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

Ectopic transcription due to inappropriately inherited histone methylation may interfere with the ongoing function of terminally differentiated cells

Many human neurodevelopmental disorders are caused by de novo mutations in histone modifying enzymes. These patients have craniofacial defects, developmental delay, intellectual disability and behavioral abnormalities, but it remains unclear how the mutations lead to such developmental defects. Here we take advantage of the invariant C. elegans lineage along with a unique double mutant in the H3K4me1/2 demethylase SPR-5/LSD1/KDM1A and the H3K9 methyltransferase MET-2/SETDB1 to address this question. We demonstrate that spr-5; met-2 double mutant worms have a severe chemotaxis defect that is dependent upon the ectopic expression of germline genes in somatic tissues. In addition, by performing single-cell RNAseq, we find that germline genes begin to be ectopically expression widely in spr-5; met-2 embryos. However, surprisingly we found that spr-5; met-2 mutants have no somatic lineage defects prior to the 200-cell stage of embryogenesis. This suggests that the altered chemotaxis behavior may be due to ongoing defect in terminally differentiated cells rather than a defect in development. To test this directly, we used RNAi to shut off the ectopic expression of germline genes in L2 spr-5; met-2 larvae, which have a fully formed nervous system. Remarkably, we find that shutting off the ectopic germline expression rescues normal chemotaxis behavior in the same adult worms that previously had a chemotaxis defect at the L2 stage. This suggests that ongoing ectopic transcription can block normal behavior in a fully intact nervous system. These data raise the possibility that intellectual disability and altered behavior in neurodevelopmental syndromes, caused by mutations in histone modifying enzymes, could be due to ongoing ectopic transcription and may be reversible.

A compound PCP scheme underlies sequential rosettes-based cell intercalation

The formation of sequential rosettes is a type of collective cell behavior recently discovered in the Caenorhabditis elegans embryo that mediates directional cell migration through sequential formation and resolution of multicellular rosettes involving the migrating cell and its neighboring cells along the way. Here, we show that a planar cell polarity (PCP)-based polarity scheme regulates sequential rosettes, which is distinct from the known mode of PCP regulation in multicellular rosettes during the process of convergent extension. Specifically, non-muscle myosin (NMY) localization and edge contraction are perpendicular to that of Van Gogh as opposed to colocalizing with Van Gogh. Further analyses suggest a two-component polarity scheme: one being the canonical PCP pathway with MIG-1/Frizzled and VANG-1/Van Gogh localized to the vertical edges, the other being MIG-1/Frizzled and NMY-2 localized to the midline/contracting edges. The NMY-2 localization and contraction of the midline edges also required LAT-1/Latrophilin, an adhesion G protein-coupled receptor that has not been shown to regulate multicellular rosettes. Our results establish a distinct mode of PCP-mediated cell intercalation and shed light on the versatile nature of the PCP pathway.

Lineage Tracing and Single-Cell RNA-seq in C. elegans to Analyze Transgenerational Epigenetic Phenotypes Inherited from Germ Cells

The last several years have seen an increasing number of examples of transgenerational epigenetic inheritance, in which phenotypes are inherited for three or more generations without changes to the underlying DNA sequence. One model system that has been particularly useful for studying transgenerational epigenetic inheritance is C. elegans. Their short generation time and hermaphroditic reproduction have allowed multiple transgenerational phenotypes to be identified, including aging, fertility, and behavior. However, it is still not clear how transgenerational epigenetic inheritance from the germline affects embryogenesis. Fortunately, the C. elegans embryo has a unique property that makes it ideal for addressing this question: they develop via an invariant lineage, with each cell undergoing stereotypical cell divisions to adopt the same cell fate in every individual embryo. Because of this invariant cell lineage, automated lineage tracing and single-cell RNA-seq can be employed to determine how transgenerational epigenetic inheritance from the germline affects developmental timing and cell fate. Unfortunately, difficulties with these techniques have severely limited their adoption in the community. Here, we provide a practical guide to automated lineage tracing coupled with single-cell RNA-seq to facilitate their use in studying transgenerational epigenetic inheritance in C. elegans embryos.

An Exact Hypergraph Matching algorithm for posture identification in embryonic C. elegans

The nematode Caenorhabditis elegans (C. elegans) is a model organism used frequently in developmental biology and neurobiology [White, (1986), Sulston, (1983), Chisholm, (2016) and Rapti, (2020)]. The C. elegans embryo can be used for cell tracking studies to understand how cell movement drives the development of specific embryonic tissues. Analyses in late-stage development are complicated by bouts of rapid twitching motions which invalidate traditional cell tracking approaches. However, the embryo possesses a small set of cells which may be identified, thereby defining the coiled embryo's posture [Christensen, 2015]. The posture serves as a frame of reference, facilitating cell tracking even in the presence of twitching. Posture identification is nevertheless challenging due to the complete repositioning of the embryo between sampled images. Current approaches to posture identification rely on time-consuming manual efforts by trained users which limits the efficiency of subsequent cell tracking. Here, we cast posture identification as a point-set matching task in which coordinates of seam cell nuclei are identified to jointly recover the posture. Most point-set matching methods comprise coherent point transformations that use low order objective functions [Zhou, (2016) and Zhang, (2019)]. Hypergraphs, an extension of traditional graphs, allow more intricate modeling of relationships between objects, yet existing hypergraphical point-set matching methods are limited to heuristic algorithms which do not easily scale to handle higher degree hypergraphs [Duchenne, (2010), Chertok, (2010) and Lee, (2011)]. Our algorithm, Exact Hypergraph Matching (EHGM), adapts the classical branch-and-bound paradigm to dynamically identify a globally optimal correspondence between point-sets under an arbitrarily intricate hypergraphical model. EHGM with hypergraphical models inspired by C. elegans embryo shape identified posture more accurately (56%) than established point-set matching methods (27%), correctly identifying twice as many sampled postures as a leading graphical approach. Posterior region seeding empowered EHGM to correctly identify 78% of postures while reducing runtime, demonstrating the efficacy of the method on a cutting-edge problem in developmental biology.

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.

Cross-modality synthesis of EM time series and live fluorescence imaging

Analyses across imaging modalities allow the integration of complementary spatiotemporal information about brain development, structure, and function. However, systematic atlasing across modalities is limited by challenges to effective image alignment. We combine highly spatially resolved electron microscopy (EM) and highly temporally resolved time-lapse fluorescence microscopy (FM) to examine the emergence of a complex nervous system in Caenorhabditis elegans embryogenesis. We generate an EM time series at four classic developmental stages and create a landmark-based co-optimization algorithm for cross-modality image alignment, which handles developmental heterochrony among datasets to achieve accurate single-cell level alignment. Synthesis based on the EM series and time-lapse FM series carrying different cell-specific markers reveals critical dynamic behaviors across scales of identifiable individual cells in the emergence of the primary neuropil, the nerve ring, as well as a major sensory organ, the amphid. Our study paves the way for systematic cross-modality data synthesis in C. elegans and demonstrates a powerful approach that may be applied broadly.

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.

Open Frontiers in Neural Cell Type Investigations; Lessons From Caenorhabditis elegans and Beyond, Toward a Multimodal Integration

Nervous system cells, the building blocks of circuits, have been studied with ever-progressing resolution, yet neural circuits appear still resistant to schemes of reductionist classification. Due to their sheer numbers, complexity and diversity, their systematic study requires concrete classifications that can serve reduced dimensionality, reproducibility, and information integration. Conventional hierarchical schemes transformed through the history of neuroscience by prioritizing criteria of morphology, (electro)physiological activity, molecular content, and circuit function, influenced by prevailing methodologies of the time. Since the molecular biology revolution and the recent advents in transcriptomics, molecular profiling gains ground toward the classification of neurons and glial cell types. Yet, transcriptomics entails technical challenges and more importantly uncovers unforeseen spatiotemporal heterogeneity, in complex and simpler nervous systems. Cells change states dynamically in space and time, in response to stimuli or throughout their developmental trajectory. Mapping cell type and state heterogeneity uncovers uncharted terrains in neurons and especially in glial cell biology, that remains understudied in many aspects. Examining neurons and glial cells from the perspectives of molecular neuroscience, physiology, development and evolution highlights the advantage of multifaceted classification schemes. Among the amalgam of models contributing to neuroscience research, Caenorhabditis elegans combines nervous system anatomy, lineage, connectivity and molecular content, all mapped at single-cell resolution, and can provide valuable insights for the workflow and challenges of the multimodal integration of cell type features. This review reflects on concepts and practices of neuron and glial cells classification and how research, in C. elegans and beyond, guides nervous system experimentation through integrated multidimensional schemes. It highlights underlying principles, emerging themes, and open frontiers in the study of nervous system development, regulatory logic and evolution. It proposes unified platforms to allow integrated annotation of large-scale datasets, gene-function studies, published or unpublished findings and community feedback. Neuroscience is moving fast toward interdisciplinary, high-throughput approaches for combined mapping of the morphology, physiology, connectivity, molecular function, and the integration of information in multifaceted schemes. A closer look in mapped neural circuits and understudied terrains offers insights for the best implementation of these approaches.

Differential adhesion regulates neurite placement via a retrograde zippering mechanism

During development, neurites and synapses segregate into specific neighborhoods or layers within nerve bundles. The developmental programs guiding placement of neurites in specific layers, and hence their incorporation into specific circuits, are not well understood. We implement novel imaging methods and quantitative models to document the embryonic development of the C. elegans brain neuropil, and discover that differential adhesion mechanisms control precise placement of single neurites onto specific layers. Differential adhesion is orchestrated via developmentally-regulated expression of the IgCAM SYG-1, and its partner ligand SYG-2. Changes in SYG-1 expression across neuropil layers result in changes in adhesive forces, which sort SYG-2-expressing neurons. Sorting to layers occurs, not via outgrowth from the neurite tip, but via an alternate mechanism of retrograde zippering, involving interactions between neurite shafts. Our study indicates that biophysical principles from differential adhesion govern neurite placement and synaptic specificity in vivo in developing neuropil bundles.

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.

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

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

Structural and developmental principles of neuropil assembly in C. elegans

Neuropil is a fundamental form of tissue organization within the brain¹, in which densely packed neurons synaptically interconnect into precise circuit architecture^2,3. However, the structural and developmental principles that govern this nanoscale precision remain largely unknown^4,5. Here we use an iterative data coarse-graining algorithm termed 'diffusion condensation'⁶ to identify nested circuit structures within the Caenorhabditis elegans neuropil, which is known as the nerve ring. We show that the nerve ring neuropil is largely organized into four strata that are composed of related behavioural circuits. The stratified architecture of the neuropil is a geometrical representation of the functional segregation of sensory information and motor outputs, with specific sensory organs and muscle quadrants mapping onto particular neuropil strata. We identify groups of neurons with unique morphologies that integrate information across strata and that create neural structures that cage the strata within the nerve ring. We use high resolution light-sheet microscopy^7,8 coupled with lineage-tracing and cell-tracking algorithms^9,10 to resolve the developmental sequence and reveal principles of cell position, migration and outgrowth that guide stratified neuropil organization. Our results uncover conserved structural design principles that underlie the architecture and function of the nerve ring neuropil, and reveal a temporal progression of outgrowth-based on pioneer neurons-that guides the hierarchical development of the layered neuropil. Our findings provide a systematic blueprint for using structural and developmental approaches to understand neuropil organization within the brain.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

Analysis software available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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

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

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

Cadherin preserves cohesion across involuting tissues during C. elegans neurulation

The internalization of the central nervous system, termed neurulation in vertebrates, is a critical step in embryogenesis. Open questions remain regarding how force propels coordinated tissue movement during the process, and little is known as to how internalization happens in invertebrates. We show that in C. elegans morphogenesis, apical constriction in the retracting pharynx drives involution of the adjacent neuroectoderm. HMR-1/cadherin mediates this process via inter-tissue attachment, as well as cohesion within the neuroectoderm. Our results demonstrate that HMR-1 is capable of mediating embryo-wide reorganization driven by a centrally located force generator, and indicate a non-canonical use of cadherin on the basal side of an epithelium that may apply to vertebrate neurulation. Additionally, we highlight shared morphology and gene expression in tissues driving involution, which suggests that neuroectoderm involution in C. elegans is potentially homologous with vertebrate neurulation and thus may help elucidate the evolutionary origin of the brain.

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 Tug-of-War between Cell Shape and Polarity Controls Division Orientation to Ensure Robust Patterning in the Mouse Blastocyst

Oriented cell division patterns tissues by modulating cell position and fate. While cell geometry, junctions, cortical tension, and polarity are known to control division orientation, relatively little is known about how these are coordinated to ensure robust patterning. Here, we systematically characterize cell division, volume, and shape changes during mouse pre-implantation development by in toto live imaging. The analysis leads us to a model in which the apical domain competes with cell shape to determine division orientation. Two key predictions of the model are verified experimentally: when outside cells of the 16-cell embryo are released from cell shape asymmetry, the axis of division is guided by the apical domain. Conversely, orientation cues from the apical domain can be overcome by applied shape asymmetry in the 8-cell embryo. We propose that such interplay between cell shape and polarity in controlling division orientation ensures robust patterning of the blastocyst and possibly other tissues.

•

Cell division, volume, and shape changes are characterized by in toto embryo imaging

•

Cell shape and the apical domain compete to determine division orientation

•

Two key predictions of the model are verified experimentally

•

The tug-of-war mechanism ensures robust cell allocation and patterning

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

Motivation

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

Results

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

Availability and implementation

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

Supplementary information

Supplementary data are available at Bioinformatics online.

Isotropic Light-Sheet Microscopy and Automated Cell Lineage Analyses to Catalogue Caenorhabditis elegans Embryogenesis with Subcellular Resolution

Caenorhabditis elegans (C. elegans) stands out as the only organism in which the challenge of understanding the cellular origins of an entire nervous system can be observed, with single cell resolution, in vivo. Here, we present an integrated protocol for the examination of neurodevelopment in C. elegans embryos. Our protocol combines imaging, lineaging and neuroanatomical tracing of single cells in developing embryos. We achieve long-term, four-dimensional (4D) imaging of living C. elegans embryos with nearly isotropic spatial resolution through the use of Dual-view Inverted Selective Plane Illumination Microscopy (diSPIM). Nuclei and neuronal structures in the nematode embryos are imaged and isotropically fused to yield images with resolution of ~330 nm in all three dimensions. These minute-by-minute high-resolution 4D data sets are then analyzed to correlate definitive cell-lineage identities with gene expression and morphological dynamics at single-cell and subcellular levels of detail. Our protocol is structured to enable modular implementation of each of the described steps and enhance studies on embryogenesis, gene expression, or neurodevelopment.

3DMMS: robust 3D Membrane Morphological Segmentation of C. elegans embryo

BACKGROUND

Understanding the cellular architecture is a fundamental problem in various biological studies. C. elegans is widely used as a model organism in these studies because of its unique fate determinations. In recent years, researchers have worked extensively on C. elegans to excavate the regulations of genes and proteins on cell mobility and communication. Although various algorithms have been proposed to analyze nucleus, cell shape features are not yet well recorded. This paper proposes a method to systematically analyze three-dimensional morphological cellular features.

RESULTS

Three-dimensional Membrane Morphological Segmentation (3DMMS) makes use of several novel techniques, such as statistical intensity normalization, and region filters, to pre-process the cell images. We then segment membrane stacks based on watershed algorithms. 3DMMS achieves high robustness and precision over different time points (development stages). It is compared with two state-of-the-art algorithms, RACE and BCOMS. Quantitative analysis shows 3DMMS performs best with the average Dice ratio of 97.7% at six time points. In addition, 3DMMS also provides time series of internal and external shape features of C. elegans.

CONCLUSION

We have developed the 3DMMS based technique for embryonic shape reconstruction at the single-cell level. With cells accurately segmented, 3DMMS makes it possible to study cellular shapes and bridge morphological features and biological expression in embryo research.

Embryo timelapses can be compiled and quantified to understand canonical histone dynamics across multiple cell cycles

In the last decade, computational processing of big datasets has facilitated the analyses of unprecedented quantities of biological data. Thus, automation and big data analysis have been revolutionary in detecting and quantifying subtle phenotypes in cell biological contexts. Analyzing similar quantities of data in larger and more complicated biological systems such as developing embryos has been more challenging due to experimental limitations on both ensemble data collection and analysis. These challenges include photosensitivity of living samples and of the fluorescently tagged proteins under study, collectively limiting the number of images that can be acquired in a single timelapse series. Here we present a streamlined workflow to quantify dynamics of fluorescently labeled proteins over the course of several cell cycles in early embryos, taking advantage of the stereotypical nature of early development that is inherent for many organisms. We benchmark this pipeline studying a fluorescently labeled histone during early embryonic development of the nematode Caenorhabditis elegans. Our strategy allowed us to overcome biological and experimental variation among our timelapse series and quantify nuclear accumulation rate, chromatin incorporation, and turnover/stability of canonical histones. We find that histone proteins are broadly stable in early C. elegans development. Thus, changes in genome regulation occurring in early development do not manifest in gross changes in histone metabolism. Our method enabling characterization of cumulative protein dynamics over several cell cycles of developmental time with high temporal resolution can be applied to expand our understanding of diverse cellular and developmental processes.

AceTree: a major update and case study in the long term maintenance of open-source scientific software

BACKGROUND

AceTree, a software application first released in 2006, facilitates exploration, curation and editing of tracked C. elegans nuclei in 4-dimensional (4D) fluorescence microscopy datasets. Since its initial release, AceTree has been continuously used to interact with, edit and interpret C. elegans lineage data. In its 11 year lifetime, AceTree has been periodically updated to meet the technical and research demands of its community of users. This paper presents the newest iteration of AceTree which contains extensive updates, demonstrates the new applicability of AceTree in other developmental contexts, and presents its evolutionary software development paradigm as a viable model for maintaining scientific software.

RESULTS

Large scale updates have been made to the user interface for an improved user experience. Tools have been grouped according to functionality and obsolete methods have been removed. Internal requirements have been changed that enable greater flexibility of use both in C. elegans contexts and in other model organisms. Additionally, the original 3-dimensional (3D) viewing window has been completely reimplemented. The new window provides a new suite of tools for data exploration.

CONCLUSION

By responding to technical advancements and research demands, AceTree has remained a useful tool for scientific research for over a decade. The updates made to the codebase have extended AceTree's applicability beyond its initial use in C. elegans and enabled its usage with other model organisms. The evolution of AceTree demonstrates a viable model for maintaining scientific software over long periods of time.

An In Toto Approach to Dissecting Cellular Interactions in Complex Tissues

Single-cell measurements have broadened our understanding of heterogeneity in biology, yet have been limited to mostly observational studies of normal or globally perturbed systems. Typically, perturbations are utilized in an open-ended approach wherein an endpoint is assayed during or after the biological event of interest. Here we describe ShootingStar, a platform for perturbation analysis in vivo, which combines live imaging, real-time image analysis, and automated optical perturbations. ShootingStar builds a quantitative record of the state of the sample being analyzed, which is used to automate the identification of target cells for perturbation, as well as to validate the impacts of the perturbation. We used ShootingStar to dissect the cellular basis of development, morphogenesis, and polarity in the lateral line of Danio rerio and the embryo of Caenorhabditis elegans. ShootingStar can be extended to diverse optical manipulations and enables more robust and informative single-cell perturbations in complex tissues.

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.

The Apical Domain Is Required and Sufficient for the First Lineage Segregation in the Mouse Embryo

Mammalian development begins with segregation of the extra-embryonic trophectoderm from the embryonic lineage in the blastocyst. While cell polarity and adhesion play key roles, the decisive cue driving this lineage segregation remains elusive. Here, to study symmetry breaking, we use a reduced system in which isolated blastomeres recapitulate the first lineage segregation. We find that in the 8-cell stage embryo, the apical domain recruits a spindle pole to ensure its differential distribution upon division. Daughter cells that inherit the apical domain adopt trophectoderm fate. However, the fate of apolar daughter cells depends on whether their position within the embryo facilitates apical domain formation by Cdh1-independent cell contact. Finally, we develop methods for transplanting apical domains and show that acquisition of this domain is not only required but also sufficient for the first lineage segregation. Thus, we provide mechanistic understanding that reconciles previous models for symmetry breaking in mouse development.

•

A reduced system was established to study symmetry breaking in mouse development

•

8-cell stage blastomeres acquire the capacity to self-organize the apical domain

•

The apical domain is required and sufficient for the first lineage segregation

•

Contact asymmetry specifies cell fate, leading to self-organized embryo patterning

ImaEdge: a platform for the quantitative analysis of cortical proteins spatiotemporal dynamics during cell polarization

Cell polarity involves the compartmentalization of the cell cortex. The establishment of cortical compartments arises from the spatial bias in the activity and concentration of cortical proteins. The mechanistic dissection of cell polarity requires the accurate detection of dynamic changes in cortical proteins, but the fluctuations of cell shape and the inhomogeneous distributions of cortical proteins greatly complicate the quantitative extraction of their global and local changes during cell polarization. To address these problems, we introduce an open-source software package, ImaEdge, which automates the segmentation of the cortex from time-lapse movies, and enables quantitative extraction of cortical protein intensities. We demonstrate that ImaEdge enables efficient and rigorous analysis of the dynamic evolution of cortical PAR proteins during C. elegans embryogenesis. It is also capable of accurate tracking of varying levels of transgene expression and discontinuous signals of the actomyosin cytoskeleton during multiple rounds of cell division. ImaEdge provides a unique resource for the quantitative studies of cortical polarization, with the potential for application to many types of polarized cells.

The AP-2 Transcription Factor APTF-2 Is Required for Neuroblast and Epidermal Morphogenesis in Caenorhabditis elegans Embryogenesis

The evolutionarily conserved family of AP-2 transcription factors (TF) regulates proliferation, differentiation, and apoptosis. Mutations in human AP-2 TF have been linked with bronchio-occular-facial syndrome and Char Syndrome, congenital birth defects characterized by craniofacial deformities and patent ductus arteriosus, respectively. How mutations in AP-2 TF cause the disease phenotypes is not well understood. Here, we characterize the aptf-2(qm27) allele in Caenorhabditis elegans, which carries a point mutation in the conserved DNA binding region of AP-2 TF. We show that compromised APTF-2 activity leads to defects in dorsal intercalation, aberrant ventral enclosure and elongation defects, ultimately culminating in the formation of morphologically deformed larvae or complete arrest during epidermal morphogenesis. Using cell lineaging, we demonstrate that APTF-2 regulates the timing of cell division, primarily in ABarp, D and C cell lineages to control the number of neuroblasts, muscle and epidermal cells. Live imaging revealed nuclear enrichment of APTF-2 in lineages affected by the qm27 mutation preceding the relevant morphogenetic events. Finally, we found that another AP-2 TF, APTF-4, is also essential for epidermal morphogenesis, in a similar yet independent manner. Thus, our study provides novel insight on the cellular-level functions of an AP-2 transcription factor in development.

Mutations in the evolutionarily conserved family of AP-2 transcription factors are associated with multiple birth defects in Char syndrome and Brancio-oculo-facial syndrome. These DNA-binding proteins are known to regulate the proliferation, differentiation and death of specific cells during embryonic development but how point mutations in the AP-2 DNA-binding domain lead to these diseases during development is currently unknown. We have identified a mutation in one of the AP-2 orthologs of the nematode Caenorhabditis elegans, APTF-2, which falls in the same mutation hotspot as in human Char syndrome and Brancio-oculo-facial syndrome patients. Compromised APTF-2 activity in C. elegans results in embryonic lethality and embryos that survive to hatching displays body morphological defects, reminiscent of the aforementioned human diseases. Using time-lapse microscopy, we found that misregulation of cell division in the skin, muscle and neuronal cell lineages is the primary cause of developmental arrest. Our study provides insight into the regulation of cell division timing by AP-2 transcription factors and provides a model to study human diseases associated with AP-2 mutations.

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.

Digital development: a database of cell lineage differentiation in C. elegans with lineage phenotypes, cell-specific gene functions and a multiscale model

Developmental systems biology is poised to exploit large-scale data from two approaches: genomics and live imaging. The combination of the two offers the opportunity to map gene functions and gene networks in vivo at single-cell resolution using cell tracking and quantification of cellular phenotypes. Here we present Digital Development (http://www.digital-development.org), a database of cell lineage differentiation with curated phenotypes, cell-specific gene functions and a multiscale model. The database stores data from recent systematic studies of cell lineage differentiation in the C. elegans embryo containing ∼ 200 conserved genes, 1400 perturbed cell lineages and 600,000 digitized single cells. Users can conveniently browse, search and download four categories of phenotypic and functional information from an intuitive web interface. This information includes lineage differentiation phenotypes, cell-specific gene functions, differentiation landscapes and fate choices, and a multiscale model of lineage differentiation. Digital Development provides a comprehensive, curated, multidimensional database for developmental biology. The scale, resolution and richness of biological information presented here facilitate exploration of gene-specific and systems-level mechanisms of lineage differentiation in Metazoans.

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.

WormGUIDES: an interactive single cell developmental atlas and tool for collaborative multidimensional data exploration

Background

Imaging and image analysis advances are yielding increasingly complete and complicated records of cellular events in tissues and whole embryos. The ability to follow hundreds to thousands of cells at the individual level demands a spatio-temporal data infrastructure: tools to assemble and collate knowledge about development spatially in a manner analogous to geographic information systems (GIS). Just as GIS indexes items or events based on their spatio-temporal or 4D location on the Earth these tools would organize knowledge based on location within the tissues or embryos. Developmental processes are highly context-specific, but the complexity of the 4D environment in which they unfold is a barrier to assembling an understanding of any particular process from diverse sources of information. In the same way that GIS aids the understanding and use of geo-located large data sets, software can, with a proper frame of reference, allow large biological data sets to be understood spatially. Intuitive tools are needed to navigate the spatial structure of complex tissue, collate large data sets and existing knowledge with this spatial structure and help users derive hypotheses about developmental mechanisms.

Results

Toward this goal we have developed WormGUIDES, a mobile application that presents a 4D developmental atlas for Caenorhabditis elegans. The WormGUIDES mobile app enables users to navigate a 3D model depicting the nuclear positions of all cells in the developing embryo. The identity of each cell can be queried with a tap, and community databases searched for available information about that cell. Information about ancestry, fate and gene expression can be used to label cells and craft customized visualizations that highlight cells as potential players in an event of interest. Scenes are easily saved, shared and published to other WormGUIDES users. The mobile app is available for Android and iOS platforms.

Conclusion

WormGUIDES provides an important tool for examining developmental processes and developing mechanistic hypotheses about their control. Critically, it provides the typical end user with an intuitive interface for developing and sharing custom visualizations of developmental processes. Equally important, because users can select cells based on their position and search for information about them, the app also serves as a spatially organized index into the large body of knowledge available to the C. elegans community online. Moreover, the app can be used to create and publish the result of exploration: interactive content that brings other researchers and students directly to the spatio-temporal point of insight. Ultimately the app will incorporate a detailed time lapse record of cell shape, beginning with neurons. This will add the key ability to navigate and understand the developmental events that result in the coordinated and precise emergence of anatomy, particularly the wiring of the nervous system.

The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis

Homeobox genes play crucial roles for the development of multicellular eukaryotes. We have generated a revised list of all homeobox genes for Caenorhabditis elegans and provide a nomenclature for the previously unnamed ones. We show that, out of 103 homeobox genes, 70 are co-orthologous to human homeobox genes. 14 are highly divergent, lacking an obvious ortholog even in other Caenorhabditis species. One of these homeobox genes encodes 12 homeodomains, while three other highly divergent homeobox genes encode a novel type of double homeodomain, termed HOCHOB. To understand how transcription factors regulate cell fate during development, precise spatio-temporal expression data need to be obtained. Using a new imaging framework that we developed, Endrov, we have generated spatio-temporal expression profiles during embryogenesis of over 60 homeobox genes, as well as a number of other developmental control genes using GFP reporters. We used dynamic feedback during recording to automatically adjust the camera exposure time in order to increase the dynamic range beyond the limitations of the camera. We have applied the new framework to examine homeobox gene expression patterns and provide an analysis of these patterns. The methods we developed to analyze and quantify expression data are not only suitable for C. elegans, but can be applied to other model systems or even to tissue culture systems.

Single Cell Quantification of Reporter Gene Expression in Live Adult Caenorhabditis elegans Reveals Reproducible Cell-Specific Expression Patterns and Underlying Biological Variation

In multicellular organisms such as Caenorhabditis elegans, differences in complex phenotypes such as lifespan correlate with the level of expression of particular engineered reporter genes. In single celled organisms, quantitative understanding of responses to extracellular signals and of cell-to-cell variation in responses has depended on precise measurement of reporter gene expression. Here, we developed microscope-based methods to quantify reporter gene expression in cells of Caenorhabditis elegans with low measurement error. We then quantified expression in strains that carried different configurations of P_hsp-16.2-fluorescent-protein reporters, in whole animals, and in all 20 cells of the intestine tissue, which is responsible for most of the fluorescent signal. Some animals bore more recently developed single copy P_hsp-16.2 reporters integrated at defined chromosomal sites, others, “classical” multicopy reporter gene arrays integrated at random sites. At the level of whole animals, variation in gene expression was similar: strains with single copy reporters showed the same amount of animal-to-animal variation as strains with multicopy reporters. At the level of cells, in animals with single copy reporters, the pattern of expression in cells within the tissue was highly stereotyped. In animals with multicopy reporters, the cell-specific expression pattern was also stereotyped, but distinct, and somewhat more variable. Our methods are rapid and gentle enough to allow quantification of expression in the same cells of an animal at different times during adult life. They should allow investigators to use changes in reporter expression in single cells in tissues as quantitative phenotypes, and link those to molecular differences. Moreover, by diminishing measurement error, they should make possible dissection of the causes of the remaining, real, variation in expression. Understanding such variation should help reveal its contribution to differences in complex phenotypic outcomes in multicellular organisms.

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.

E3 ubiquitin ligases promote progression of differentiation during C. elegans embryogenesis

Regulated choice between cell fate maintenance and differentiation provides decision points in development to progress toward more restricted cell fates or to maintain the current one. Caenorhabditis elegans embryogenesis follows an invariant cell lineage where cell fate is generally more restricted upon each cell division. EMS is a progenitor cell in the four-cell embryo that gives rise to the endomesoderm. We recently found that when ubiquitin-mediated protein degradation is compromised, the anterior daughter of EMS, namely MS, reiterates the EMS fate. This observation demonstrates an essential function of ubiquitin-mediated protein degradation in driving the progression of EMS-to-MS differentiation. Here we report a genome-wide screen of the ubiquitin pathway and extensive lineage analyses. The results suggest a broad role of E3 ligases in driving differentiation progression. First, we identified three substrate-binding proteins for two Cullin-RING ubiquitin ligase (CRL) E3 complexes that promote the progression from the EMS fate to MS, namely LIN-23/β-TrCP and FBXB-3 for the CRL1/SCF complex and ZYG-11/ZYG-11B for the CRL2 complex. Genetic analyses suggest these E3 ligases function through a multifunctional protein OMA-1 and the endomesoderm lineage specifier SKN-1 to drive differentiation. Second, we found that depletion of components of the CRL1/SCF complex induces fate reiteration in all major founder cell lineages. These data suggest that regulated choice between self-renewal and differentiation is widespread during C. elegans embryogenesis as in organisms with regulative development, and ubiquitin-mediated protein degradation drives the choice towards differentiation. Finally, bioinformatic analysis of time series gene expression data showed that expression of E3 genes is transiently enriched during time windows of developmental stage transitions. Transcription factors show similar enrichment, but not other classes of regulatory genes. Based on these findings we propose that ubiquitin-mediated protein degradation, like many transcription factors, function broadly as regulators driving developmental progression during embryogenesis in C. elegans.

Biological Dynamics Markup Language (BDML): an open format for representing quantitative biological dynamics data

Motivation: Recent progress in live-cell imaging and modeling techniques has resulted in generation of a large amount of quantitative data (from experimental measurements and computer simulations) on spatiotemporal dynamics of biological objects such as molecules, cells and organisms. Although many research groups have independently dedicated their efforts to developing software tools for visualizing and analyzing these data, these tools are often not compatible with each other because of different data formats.

Results: We developed an open unified format, Biological Dynamics Markup Language (BDML; current version: 0.2), which provides a basic framework for representing quantitative biological dynamics data for objects ranging from molecules to cells to organisms. BDML is based on Extensible Markup Language (XML). Its advantages are machine and human readability and extensibility. BDML will improve the efficiency of development and evaluation of software tools for data visualization and analysis.

Availability and implementation: A specification and a schema file for BDML are freely available online at http://ssbd.qbic.riken.jp/bdml/.

Contact:sonami@riken.jp

Supplementary Information:Supplementary data are available at Bioinformatics online.

Irises: A practical tool for image-based analysis of cellular DNA content

The DNA content of nuclei is a valuable measure of cell cycle status. Irises is a software tool to facilitate systematic in situ determination of DNA content for cell cycle analysis at single-nucleus resolution within complex tissues. We demonstrate the utility of the tool with analysis of DNA content in germline nuclei of C. elegans. Compared with results obtained by manual analysis, we find the tool greatly facilitates analysis by improving speed at least 5-fold while maintaining accuracy. The source code and instruction manual (including installation for both Mac and PC) are provided.

De novo inference of systems-level mechanistic models of development from live-imaging-based phenotype analysis

Elucidation of complex phenotypes for mechanistic insights presents a significant challenge in systems biology. We report a strategy to automatically infer mechanistic models of cell fate differentiation based on live-imaging data. We use cell lineage tracing and combinations of tissue-specific marker expression to assay progenitor cell fate and detect fate changes upon genetic perturbation. Based on the cellular phenotypes, we further construct a model for how fate differentiation progresses in progenitor cells and predict cell-specific gene modules and cell-to-cell signaling events that regulate the series of fate choices. We validate our approach in C. elegans embryogenesis by perturbing 20 genes in over 300 embryos. The result not only recapitulates current knowledge but also provides insights into gene function and regulated fate choice, including an unexpected self-renewal. Our study provides a powerful approach for automated and quantitative interpretation of complex in vivo information.

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.

Towards 3D in silico modeling of the sea urchin embryonic development

Embryogenesis is a dynamic process with an intrinsic variability whose understanding requires the integration of molecular, genetic, and cellular dynamics. Biological circuits function over time at the level of single cells and require a precise analysis of the topology, temporality, and probability of events. Integrative developmental biology is currently looking for the appropriate strategies to capture the intrinsic properties of biological systems. The "-omic" approaches require disruption of the function of the biological circuit; they provide static information, with low temporal resolution and usually with population averaging that masks fast or variable features at the cellular scale and in a single individual. This data should be correlated with cell behavior as cells are the integrators of biological activity. Cellular dynamics are captured by the in vivo microscopy observation of live organisms. This can be used to reconstruct the 3D + time cell lineage tree to serve as the basis for modeling the organism's multiscale dynamics. We discuss here the progress that has been made in this direction, starting with the reconstruction over time of three-dimensional digital embryos from in toto time-lapse imaging. Digital specimens provide the means for a quantitative description of the development of model organisms that can be stored, shared, and compared. They open the way to in silico experimentation and to a more theoretical approach to biological processes. We show, with some unpublished results, how the proposed methodology can be applied to sea urchin species that have been model organisms in the field of classical embryology and modern developmental biology for over a century.

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.

A high-fidelity cell lineage tracing method for obtaining systematic spatiotemporal gene expression patterns in Caenorhabditis elegans

Advances in microscopy and fluorescent reporters have allowed us to detect the onset of gene expression on a cell-by-cell basis in a systemic fashion. This information, however, is often encoded in large repositories of images, and developing ways to extract this spatiotemporal expression data is a difficult problem that often uses complex domain-specific methods for each individual data set. We present a more unified approach that incorporates general previous information into a hierarchical probabilistic model to extract spatiotemporal gene expression from 4D confocal microscopy images of developing Caenorhabditis elegans embryos. This approach reduces the overall error rate of our automated lineage tracing pipeline by 3.8-fold, allowing us to routinely follow the C. elegans lineage to later stages of development, where individual neuronal subspecification becomes apparent. Unlike previous methods that often use custom approaches that are organism specific, our method uses generalized linear models and extensions of standard reversible jump Markov chain Monte Carlo methods that can be readily extended to other organisms for a variety of biological inference problems relating to cell fate specification. This modeling approach is flexible and provides tractable avenues for incorporating additional previous information into the model for similar difficult high-fidelity/low error tolerance image analysis problems for systematically applied genomic experiments.

The secretory pathway calcium ATPase PMR-1/SPCA1 has essential roles in cell migration during Caenorhabditis elegans embryonic development

Maintaining levels of calcium in the cytosol is important for many cellular events, including cell migration, where localized regions of high calcium are required to regulate cytoskeletal dynamics, contractility, and adhesion. Studies show inositol-trisphosphate receptors (IP3R) and ryanodine receptors (RyR), which release calcium into the cytosol, are important regulators of cell migration. Similarly, proteins that return calcium to secretory stores are likely to be important for cell migration. The secretory protein calcium ATPase (SPCA) is a Golgi-localized protein that transports calcium from the cytosol into secretory stores. SPCA has established roles in protein processing, metal homeostasis, and inositol-trisphosphate signaling. Defects in the human SPCA1/ATP2C1 gene cause Hailey-Hailey disease (MIM# 169600), a genodermatosis characterized by cutaneous blisters and fissures as well as keratinocyte cell adhesion defects. We have determined that PMR-1, the Caenorhabditis elegans ortholog of SPCA1, plays an essential role in embryogenesis. Pmr-1 strains isolated from genetic screens show terminal phenotypes, such as ventral and anterior enclosure failures, body morphogenesis defects, and an unattached pharynx, which are caused by earlier defects during gastrulation. In Pmr-1 embryos, migration rates are significantly reduced for cells moving along the embryo surface, such as ventral neuroblasts, C-derived, and anterior-most blastomeres. Gene interaction experiments show changing the activity of itr-1/IP3R and unc-68/RyR modulates levels of embryonic lethality in Pmr-1 strains, indicating pmr-1 acts with these calcium channels to regulate cell migration. This analysis reveals novel genes involved in C. elegans cell migration, as well as a new role in cell migration for the highly conserved SPCA gene family.

During growth or regeneration after damage, skin cells migrate from basal to superficial layers, forming tight attachments that protect an individual from environmental assaults. Proteins that remove calcium from the cell cytosol into secretory stores, where it is available for future release, play a key role in skin cell integrity. Defects in these secretory pathway calcium ATPase (SPCA) channels in humans cause Hailey-Hailey disease, a chronic disorder marked by skin lesions in areas of high-stress. Our study of the SPCA gene pmr-1 in Caenorhabditis elegans indicates the gene is essential for viability. Embryos with defective PMR-1 die with cell attachment defects superficially similar to those of Hailey-Hailey disease patients. To better understand this phenotype, we tracked the position of individual cells during development of pmr-1 mutant embryos. This analysis revealed that the cell attachment defects are caused by primary failures in cell migration. We also identified other calcium channel proteins involved in this process, indicating proper regulation of calcium is crucial for cell migration in C. elegans. If SPCA proteins act similarly in humans, this research will lead to better understanding of the molecules important for skin cell regeneration, as well as help to explain the defects observed in Hailey-Hailey disease patients.

Chapter 17: bioimage informatics for systems pharmacology

Recent advances in automated high-resolution fluorescence microscopy and robotic handling have made the systematic and cost effective study of diverse morphological changes within a large population of cells possible under a variety of perturbations, e.g., drugs, compounds, metal catalysts, RNA interference (RNAi). Cell population-based studies deviate from conventional microscopy studies on a few cells, and could provide stronger statistical power for drawing experimental observations and conclusions. However, it is challenging to manually extract and quantify phenotypic changes from the large amounts of complex image data generated. Thus, bioimage informatics approaches are needed to rapidly and objectively quantify and analyze the image data. This paper provides an overview of the bioimage informatics challenges and approaches in image-based studies for drug and target discovery. The concepts and capabilities of image-based screening are first illustrated by a few practical examples investigating different kinds of phenotypic changes caEditorsused by drugs, compounds, or RNAi. The bioimage analysis approaches, including object detection, segmentation, and tracking, are then described. Subsequently, the quantitative features, phenotype identification, and multidimensional profile analysis for profiling the effects of drugs and targets are summarized. Moreover, a number of publicly available software packages for bioimage informatics are listed for further reference. It is expected that this review will help readers, including those without bioimage informatics expertise, understand the capabilities, approaches, and tools of bioimage informatics and apply them to advance their own studies.