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

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

[1]  Bonnie Berger,et al.  Efficient integration of heterogeneous single-cell transcriptomes using Scanorama , 2019, Nature Biotechnology.

[2]  Garry P. Nolan,et al.  Single Cell Trajectory Detection Orders Hallmarks of Early Human B Cell Development , 2012 .

[3]  Paul W. Sternberg,et al.  WormBase: network access to the genome and biology of Caenorhabditis elegans , 2001, Nucleic Acids Res..

[4]  J. Cooke,et al.  Control of somite number during morphogenesis of a vertebrate, Xenopus laevis , 1975, Nature.

[5]  Jonathan Hodgkin,et al.  The C. elegans Hox gene egl-5 is required for correct development of the hermaphrodite hindgut and for the response to rectal infection by Microbacterium nematophilum. , 2009, Developmental biology.

[6]  B. Reinhart,et al.  The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans , 2000, Nature.

[7]  Ruiqiang Li,et al.  Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells , 2013, Nature Structural &Molecular Biology.

[8]  Ashley R Bonneau,et al.  Zygotic genome activation during the maternal-to-zygotic transition. , 2014, Annual review of cell and developmental biology.

[9]  Bonnie Berger,et al.  Neural Data Visualization for Scalable and Generalizable Single Cell Analysis , 2018, bioRxiv.

[10]  Leon Anavy,et al.  BLIND ordering of large-scale transcriptomic developmental timecourses , 2014, Development.

[11]  Thomas J. Nicholas,et al.  Automated analysis of embryonic gene expression with cellular resolution in C. elegans , 2008, Nature Methods.

[12]  O. Pourquié,et al.  A clock-work somite. , 2000, BioEssays : news and reviews in molecular, cellular and developmental biology.

[13]  Hao Wang,et al.  Collaborative Regulation of Development but Independent Control of Metabolism by Two Epidermis-specific Transcription Factors in Caenorhabditis elegans* , 2013, The Journal of Biological Chemistry.

[14]  Alexander N. Gorban,et al.  Principal Graphs and Manifolds , 2008, ArXiv.

[15]  Andrew C. Adey,et al.  Single-Cell Transcriptional Profiling of a Multicellular Organism , 2017 .

[16]  John Isaac Murray,et al.  Gene transcription is coordinated with, but not dependent on, cell divisions during C. elegans embryonic fate specification , 2013, Development.

[17]  A. Rougvie,et al.  Caenorhabditis elegans period homolog lin-42 regulates the timing of heterochronic miRNA expression , 2014, Proceedings of the National Academy of Sciences.

[18]  J. McGhee,et al.  Reevaluation of the Role of the med-1 and med-2 Genes in Specifying the Caenorhabditis elegans Endoderm , 2005, Genetics.

[19]  Joshua T. Burdick,et al.  A quantitative model of normal Caenorhabditis elegans embryogenesis and its disruption after stress. , 2013, Developmental biology.

[20]  R. Waterston,et al.  Multidimensional regulation of gene expression in the C. elegans embryo , 2012, Genome research.

[21]  Dimos Gaidatzis,et al.  Extensive oscillatory gene expression during C. elegans larval development. , 2014, Molecular cell.

[22]  Eugene W. Myers,et al.  Analysis of Cell Fate from Single-Cell Gene Expression Profiles in C. elegans , 2009, Cell.

[23]  John Isaac Murray,et al.  The lineaging of fluorescently-labeled Caenorhabditis elegans embryos with StarryNite and AceTree , 2006, Nature Protocols.

[24]  Allon M. Klein,et al.  Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo , 2018, Science.

[25]  Emilio Soria Olivas,et al.  Handbook of Research on Machine Learning Applications and Trends : Algorithms , Methods , and Techniques , 2009 .

[26]  R. Waterston,et al.  Automated cell lineage tracing in Caenorhabditis elegans. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Daniel L. Mace,et al.  Regulatory analysis of the C. elegans genome with spatiotemporal resolution , 2014, Nature.

[28]  Alexander van Oudenaarden,et al.  Dampening of expression oscillations by synchronous regulation of a microRNA and its target , 2013, Nature Genetics.

[29]  Nektarios Tavernarakis,et al.  The Homeobox Genes of Caenorhabditis elegans and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis , 2015, PloS one.

[30]  Allon M. Klein,et al.  The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution , 2018, Science.

[31]  N. Satoh,et al.  Timing mechanisms in early embryonic development. , 1982, Differentiation; research in biological diversity.

[32]  Eytan Domany,et al.  Using high-throughput transcriptomic data for prognosis: a critical overview and perspectives. , 2014, Cancer research.

[33]  Bonnie Berger,et al.  Generalizable and Scalable Visualization of Single-Cell Data Using Neural Networks. , 2018, Cell systems.

[34]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[35]  Andrew C. Adey,et al.  Comprehensive single cell transcriptional profiling of a multicellular organism by combinatorial indexing , 2017 .

[36]  Sean C. Bendall,et al.  Single-Cell Trajectory Detection Uncovers Progression and Regulatory Coordination in Human B Cell Development , 2014, Cell.

[37]  Vachiranee Limviphuvadh,et al.  Developmental fate and lineage commitment of singled mouse blastomeres , 2012, Development.

[38]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[39]  Fabian J Theis,et al.  SCANPY: large-scale single-cell gene expression data analysis , 2018, Genome Biology.

[40]  N. Neff,et al.  Reconstructing lineage hierarchies of the distal lung epithelium using single cell RNA-seq , 2014, Nature.

[41]  William B. Wood,et al.  A Negative Regulatory Loop between MicroRNA and Hox Gene Controls Posterior Identities in Caenorhabditis elegans , 2010, PLoS genetics.

[42]  S. Teichmann,et al.  From Tissues to Cell Types and Back: Single-Cell Gene Expression Analysis of Tissue Architecture , 2018, Annual Review of Biomedical Data Science.

[43]  Min Han,et al.  The Caenorhabditis elegans Nuclear Receptor Gene nhr-25 Regulates Epidermal Cell Development , 2004, Molecular and Cellular Biology.

[44]  A. Regev,et al.  Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis , 2018, Science.

[45]  W. B. Wood,et al.  Caenorhabditis elegans embryonic axial patterning requires two recently discovered posterior-group Hox genes. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Wouter Houthoofd,et al.  The embryonic cell lineage of the nematode Halicephalobus gingivalis (Nematoda: Cephalobina: Panagrolaimoidea) , 2007 .

[47]  Yoshiki Andachi,et al.  Caenorhabditis elegans T‐box genes tbx‐9 and tbx‐8 are required for formation of hypodermis and body‐wall muscle in embryogenesis , 2004, Genes to cells : devoted to molecular & cellular mechanisms.

[48]  Zhirong Bao,et al.  AceTree: a tool for visual analysis of Caenorhabditis elegans embryogenesis , 2006, BMC Bioinformatics.

[49]  Jessica S. Dymond,et al.  The Paired-box protein PAX-3 regulates the choice between lateral and ventral epidermal cell fates in C. elegans. , 2016, Developmental biology.

[50]  J. McGhee,et al.  Reprogramming of early embryonic blastomeres into endodermal progenitors by a Caenorhabditis elegans GATA factor. , 1998, Genes & development.

[51]  O. Pourquié,et al.  Avian hairy Gene Expression Identifies a Molecular Clock Linked to Vertebrate Segmentation and Somitogenesis , 1997, Cell.

[52]  Zhirong Bao,et al.  Systematic quantification of developmental phenotypes at single-cell resolution during embryogenesis , 2013, Development.

[53]  Bryan D. Bryson,et al.  Panoramic stitching of heterogeneous single-cell transcriptomic data , 2018, bioRxiv.

[54]  R. Lin,et al.  Uncoupling Different Characteristics of the C. elegans E Lineage from Differentiation of Intestinal Markers , 2014, PloS one.

[55]  S. Mcconnell,et al.  Progressive restriction in fate potential by neural progenitors during cerebral cortical development. , 2000, Development.

[56]  R. Sandberg,et al.  Full-Length mRNA-Seq from single cell levels of RNA and individual circulating tumor cells , 2012, Nature Biotechnology.

[57]  J. Sulston,et al.  The embryonic cell lineage of the nematode Caenorhabditis elegans. , 1983, Developmental biology.

[58]  Claudiu A. Giurumescu,et al.  Quantitative semi-automated analysis of morphogenesis with single-cell resolution in complex embryos , 2012, Development.