{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,2,26]],"date-time":"2026-02-26T20:35:00Z","timestamp":1772138100329,"version":"3.50.1"},"reference-count":54,"publisher":"Oxford University Press (OUP)","issue":"13","license":[{"start":{"date-parts":[[2019,12,20]],"date-time":"2019-12-20T00:00:00Z","timestamp":1576800000000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by-nc\/4.0\/"}],"funder":[{"DOI":"10.13039\/100000002","name":"NIH","doi-asserted-by":"publisher","award":["R01GM081871"],"award-info":[{"award-number":["R01GM081871"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2020,7,1]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:sec>\n                    <jats:title>Motivation<\/jats:title>\n                    <jats:p>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.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Results<\/jats:title>\n                    <jats:p>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\u2019s 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\u2019s 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.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Availability and implementation<\/jats:title>\n                    <jats:p>Analysis software available upon request.<\/jats:p>\n                  <\/jats:sec>\n                  <jats:sec>\n                    <jats:title>Supplementary information<\/jats:title>\n                    <jats:p>Supplementary data are available at Bioinformatics online.<\/jats:p>\n                  <\/jats:sec>","DOI":"10.1093\/bioinformatics\/btz864","type":"journal-article","created":{"date-parts":[[2019,12,18]],"date-time":"2019-12-18T15:12:41Z","timestamp":1576681961000},"page":"4047-4057","source":"Crossref","is-referenced-by-count":7,"title":["Meta-analysis of\n                    <i>Caenorhabditis elegans<\/i>\n                    single-cell developmental data reveals multi-frequency oscillation in gene activation"],"prefix":"10.1093","volume":"36","author":[{"given":"Luke A D","family":"Hutchison","sequence":"first","affiliation":[{"name":"MIT Computer Science and AI Lab , Cambridge, MA 02139, USA"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-2724-7228","authenticated-orcid":false,"given":"Bonnie","family":"Berger","sequence":"additional","affiliation":[{"name":"MIT Computer Science and AI Lab , Cambridge, MA 02139, USA"}]},{"given":"Isaac S","family":"Kohane","sequence":"additional","affiliation":[{"name":"Harvard Medical School , Boston, MA 02115, USA"}]}],"member":"286","published-online":{"date-parts":[[2019,12,20]]},"reference":[{"key":"2023062312035057900_btz864-B1","doi-asserted-by":"crossref","first-page":"1161","DOI":"10.1242\/dev.105288","article-title":"BLIND ordering of large-scale transcriptomic developmental timecourses","volume":"141","author":"Anavy","year":"2014","journal-title":"Development"},{"key":"2023062312035057900_btz864-B2","doi-asserted-by":"crossref","first-page":"331","DOI":"10.1111\/j.1356-9597.2004.00725.x","article-title":"Caenorhabditis elegans T-box genes tbx-9 and tbx-8 are required for formation of hypodermis and body-wall muscle in embryogenesis","volume":"9","author":"Andachi","year":"2004","journal-title":"Genes Cells"},{"key":"2023062312035057900_btz864-B3","doi-asserted-by":"crossref","first-page":"400","DOI":"10.1038\/nature13497","article-title":"Regulatory analysis of the C. elegans genome with spatiotemporal resolution","volume":"512","author":"Araya","year":"2014","journal-title":"Nature"},{"key":"2023062312035057900_btz864-B4","doi-asserted-by":"crossref","first-page":"2707","DOI":"10.1073\/pnas.0511111103","article-title":"Automated cell lineage tracing in Caenorhabditis elegans","volume":"103","author":"Bao","year":"2006","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062312035057900_btz864-B5","doi-asserted-by":"crossref","first-page":"714","DOI":"10.1016\/j.cell.2014.04.005","article-title":"Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development","volume":"157","author":"Bendall","year":"2014","journal-title":"Cell"},{"key":"2023062312035057900_btz864-B6","doi-asserted-by":"crossref","first-page":"275","DOI":"10.1186\/1471-2105-7-275","article-title":"AceTree: a tool for visual analysis of Caenorhabditis elegans embryogenesis","volume":"7","author":"Boyle","year":"2006","journal-title":"BMC Bioinformatics"},{"key":"2023062312035057900_btz864-B7","doi-asserted-by":"crossref","first-page":"eaar5780","DOI":"10.1126\/science.aar5780","article-title":"The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution","volume":"360","author":"Briggs","year":"2018","journal-title":"Science"},{"key":"2023062312035057900_btz864-B8","doi-asserted-by":"crossref","first-page":"661","DOI":"10.1126\/science.aam8940","article-title":"Comprehensive single cell transcriptional profiling of a multicellular organism by combinatorial indexing","volume":"357","author":"Cao","year":"2017","journal-title":"Science"},{"key":"2023062312035057900_btz864-B9","doi-asserted-by":"crossref","first-page":"29","DOI":"10.1146\/annurev-biodatasci-080917-013452","article-title":"From tissues to cell types and back: single-cell gene expression analysis of tissue architecture","volume":"1","author":"Chen","year":"2018","journal-title":"Ann. Rev. Biomed. Data Sci"},{"key":"2023062312035057900_btz864-B10","doi-asserted-by":"crossref","first-page":"7345","DOI":"10.1128\/MCB.24.17.7345-7358.2004","article-title":"The Caenorhabditis elegans nuclear receptor gene nhr-25 regulates epidermal cell development","volume":"24","author":"Chen","year":"2004","journal-title":"Mol. Cell. Biol"},{"key":"2023062312035057900_btz864-B11","doi-asserted-by":"crossref","first-page":"185","DOI":"10.1016\/j.cels.2018.05.017","article-title":"Neural data visualization for scalable and generalizable single cell analysis","volume":"7","author":"Cho","year":"2018","journal-title":"Cell Syst"},{"key":"2023062312035057900_btz864-B12","doi-asserted-by":"crossref","first-page":"196","DOI":"10.1038\/254196a0","article-title":"Control of somite number during morphogenesis of a vertebrate, Xenopus laevis","volume":"254","author":"Cooke","year":"1975","journal-title":"Nature"},{"key":"2023062312035057900_btz864-B13","doi-asserted-by":"crossref","first-page":"72","DOI":"10.1002\/(SICI)1521-1878(200001)22:1<72::AID-BIES12>3.0.CO;2-S","article-title":"A clock-work somite","volume":"22","author":"Dale","year":"2000","journal-title":"Bioessays"},{"key":"2023062312035057900_btz864-B14","doi-asserted-by":"crossref","first-page":"1044","DOI":"10.1182\/blood.V120.21.1044.1044","article-title":"Single cell trajectory detection orders hallmarks of early human B cell development","volume":"120","author":"Davis","year":"2012","journal-title":"Blood"},{"key":"2023062312035057900_btz864-B15","doi-asserted-by":"crossref","first-page":"2863","DOI":"10.1242\/dev.127.13.2863","article-title":"Progressive restriction in fate potential by neural progenitors during cerebral cortical development","volume":"127","author":"Desai","year":"2000","journal-title":"Development"},{"key":"2023062312035057900_btz864-B16","doi-asserted-by":"crossref","first-page":"4612","DOI":"10.1158\/0008-5472.CAN-13-3338","article-title":"Using high-throughput transcriptomic data for prognosis: a critical overview and perspectives","volume":"74","author":"Domany","year":"2014","journal-title":"Cancer Res"},{"key":"2023062312035057900_btz864-B17","doi-asserted-by":"crossref","first-page":"eaar3131","DOI":"10.1126\/science.aar3131","article-title":"Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis","volume":"360","author":"Farrell","year":"2018","journal-title":"Science"},{"key":"2023062312035057900_btz864-B18","doi-asserted-by":"crossref","first-page":"4271","DOI":"10.1242\/dev.086256","article-title":"Quantitative semi-automated analysis of morphogenesis with single-cell resolution in complex embryos","volume":"139","author":"Giurumescu","year":"2012","journal-title":"Development"},{"key":"2023062312035057900_btz864-B19","first-page":"28","volume-title":"Handbook of Research on Machine Learning Applications and Trends: Algorithms, Methods, and Techniques","author":"Gorban","year":"2009"},{"key":"2023062312035057900_btz864-B20","doi-asserted-by":"crossref","first-page":"545","DOI":"10.1534\/genetics.105.044909","article-title":"Reevaluation of the role of the med-1 and med-2 genes in specifying the Caenorhabditis elegans endoderm","volume":"171","author":"Goszczynski","year":"2005","journal-title":"Genetics"},{"key":"2023062312035057900_btz864-B21","doi-asserted-by":"crossref","first-page":"e0126947","DOI":"10.1371\/journal.pone.0126947","article-title":"The homeobox genes of Caenorhabditis elegans and insights into their spatio-temporal expression dynamics during embryogenesis","volume":"10","author":"Hench","year":"2015","journal-title":"PLoS One"},{"key":"2023062312035057900_btz864-B22","doi-asserted-by":"crossref","first-page":"380","DOI":"10.1016\/j.molcel.2013.12.013","article-title":"Extensive oscillatory gene expression during C. elegans larval development","volume":"53","author":"Hendriks","year":"2014","journal-title":"Mol. Cell"},{"key":"2023062312035057900_btz864-B23","doi-asserted-by":"crossref","first-page":"685","DOI":"10.1038\/s41587-019-0113-3","article-title":"Panoramic stitching of heterogeneous single-cell transcriptomic data","volume":"37","author":"Hie","year":"2019","journal-title":"Nat. Biotechnol"},{"key":"2023062312035057900_btz864-B24","doi-asserted-by":"crossref","first-page":"1337","DOI":"10.1038\/ng.2763","article-title":"Dampening of expression oscillations by synchronous regulation of a microRNA and its target","volume":"45","author":"Kim","year":"2013","journal-title":"Nat. Genet"},{"key":"2023062312035057900_btz864-B25","doi-asserted-by":"crossref","first-page":"581","DOI":"10.1146\/annurev-cellbio-100913-013027","article-title":"Zygotic genome activation during the maternal-to-zygotic transition","volume":"30","author":"Lee","year":"2014","journal-title":"Annu. Rev. Cell Dev. Biol"},{"key":"2023062312035057900_btz864-B26","doi-asserted-by":"crossref","first-page":"623","DOI":"10.1016\/j.cell.2009.08.044","article-title":"Analysis of cell fate from single-cell gene expression profiles in C. elegans","volume":"139","author":"Liu","year":"2009","journal-title":"Cell"},{"key":"2023062312035057900_btz864-B27","doi-asserted-by":"crossref","first-page":"3722","DOI":"10.1242\/dev.086454","article-title":"Developmental fate and lineage commitment of singled mouse blastomeres","volume":"139","author":"Lorthongpanich","year":"2012","journal-title":"Development"},{"key":"2023062312035057900_btz864-B28","doi-asserted-by":"crossref","first-page":"15450","DOI":"10.1073\/pnas.1414856111","article-title":"Caenorhabditis elegans period homolog lin-42 regulates the timing of heterochronic miRNA expression","volume":"111","author":"McCulloch","year":"2014","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062312035057900_btz864-B29","doi-asserted-by":"crossref","first-page":"3266","DOI":"10.1242\/dev.096040","article-title":"Systematic quantification of developmental phenotypes at single-cell resolution during embryogenesis","volume":"140","author":"Moore","year":"2013","journal-title":"Development"},{"key":"2023062312035057900_btz864-B30","doi-asserted-by":"crossref","first-page":"1468","DOI":"10.1038\/nprot.2006.222","article-title":"The lineaging of fluorescently-labeled Caenorhabditis elegans embryos with StarryNite and AceTree","volume":"1","author":"Murray","year":"2006","journal-title":"Nat. Protoc"},{"key":"2023062312035057900_btz864-B31","doi-asserted-by":"crossref","first-page":"703","DOI":"10.1038\/nmeth.1228","article-title":"Automated analysis of embryonic gene expression with cellular resolution in C. elegans","volume":"5","author":"Murray","year":"2008","journal-title":"Nat. Methods"},{"key":"2023062312035057900_btz864-B32","doi-asserted-by":"crossref","first-page":"1282","DOI":"10.1101\/gr.131920.111","article-title":"Multidimensional regulation of gene expression in the C. elegans embryo","volume":"22","author":"Murray","year":"2012","journal-title":"Genome Res"},{"key":"2023062312035057900_btz864-B33","doi-asserted-by":"crossref","first-page":"3385","DOI":"10.1242\/dev.098012","article-title":"Gene transcription is coordinated with, but not dependent on, cell divisions during C. elegans embryonic fate specification","volume":"140","author":"Nair","year":"2013","journal-title":"Development"},{"key":"2023062312035057900_btz864-B34","doi-asserted-by":"crossref","first-page":"16","DOI":"10.1016\/j.ydbio.2009.01.044","article-title":"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","volume":"329","author":"Nicholas","year":"2009","journal-title":"Dev. Biol"},{"key":"2023062312035057900_btz864-B35","doi-asserted-by":"crossref","first-page":"639","DOI":"10.1016\/S0092-8674(00)80451-1","article-title":"Avian hairy gene expression identifies a molecular clock linked to vertebrate segmentation and somitogenesis","volume":"91","author":"Palmeirim","year":"1997","journal-title":"Cell"},{"key":"2023062312035057900_btz864-B36","doi-asserted-by":"crossref","first-page":"979","DOI":"10.1038\/nmeth.4402","article-title":"Reversed graph embedding resolves complex single-cell trajectories","volume":"14","author":"Qiu","year":"2017","journal-title":"Nat. Methods"},{"key":"2023062312035057900_btz864-B37","doi-asserted-by":"crossref","first-page":"777","DOI":"10.1038\/nbt.2282","article-title":"Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells","volume":"30","author":"Ramsk\u00f6ld","year":"2012","journal-title":"Nat. Biotechnol"},{"key":"2023062312035057900_btz864-B38","doi-asserted-by":"crossref","first-page":"901","DOI":"10.1038\/35002607","article-title":"The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans","volume":"403","author":"Reinhart","year":"2000","journal-title":"Nature"},{"key":"2023062312035057900_btz864-B39","doi-asserted-by":"crossref","first-page":"12","DOI":"10.1016\/j.ydbio.2012.11.034","article-title":"A quantitative model of normal Caenorhabditis elegans embryogenesis and its disruption after stress","volume":"374","author":"Richards","year":"2013","journal-title":"Dev. Biol"},{"key":"2023062312035057900_btz864-B40","doi-asserted-by":"crossref","first-page":"e106309","DOI":"10.1371\/journal.pone.0106309","article-title":"Uncoupling different characteristics of the C. elegans E lineage from differentiation of intestinal markers","volume":"9","author":"Robertson","year":"2014","journal-title":"PLoS One"},{"key":"2023062312035057900_btz864-B41","doi-asserted-by":"crossref","first-page":"156","DOI":"10.1111\/j.1432-0436.1982.tb01243.x","article-title":"Timing mechanisms in early embryonic development","volume":"22","author":"Satoh","year":"1982","journal-title":"Differentiation"},{"key":"2023062312035057900_btz864-B42","doi-asserted-by":"crossref","first-page":"33411","DOI":"10.1074\/jbc.M113.487975","article-title":"Collaborative regulation of development but independent control of metabolism by two epidermis-specific transcription factors in Caenorhabditis elegans","volume":"288","author":"Shao","year":"2013","journal-title":"J. Biol. Chem"},{"key":"2023062312035057900_btz864-B43","doi-asserted-by":"crossref","first-page":"82","DOI":"10.1093\/nar\/29.1.82","article-title":"WormBase: network access to the genome and biology of Caenorhabditis elegans","volume":"29","author":"Stein","year":"2001","journal-title":"Nucleic Acids Res"},{"key":"2023062312035057900_btz864-B44","doi-asserted-by":"crossref","first-page":"64","DOI":"10.1016\/0012-1606(83)90201-4","article-title":"The embryonic cell lineage of the nematode Caenorhabditis elegans","volume":"100","author":"Sulston","year":"1983","journal-title":"Dev. Biol"},{"key":"2023062312035057900_btz864-B45","doi-asserted-by":"crossref","first-page":"191","DOI":"10.1016\/j.ydbio.2016.03.002","article-title":"The paired-box protein PAX-3 regulates the choice between lateral and ventral epidermal cell fates in C. elegans","volume":"412","author":"Thompson","year":"2016","journal-title":"Dev. Biol"},{"key":"2023062312035057900_btz864-B46","doi-asserted-by":"crossref","first-page":"381","DOI":"10.1038\/nbt.2859","article-title":"The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells","volume":"32","author":"Trapnell","year":"2014","journal-title":"Nat. Biotechnol"},{"key":"2023062312035057900_btz864-B47","doi-asserted-by":"crossref","first-page":"371","DOI":"10.1038\/nature13173","article-title":"Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seq","volume":"509","author":"Treutlein","year":"2014","journal-title":"Nature"},{"key":"2023062312035057900_btz864-B48","doi-asserted-by":"crossref","first-page":"4499","DOI":"10.1073\/pnas.97.9.4499","article-title":"Caenorhabditis elegans embryonic axial patterning requires two recently discovered posterior-group Hox genes","volume":"97","author":"Van Auken","year":"2000","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"2023062312035057900_btz864-B49","doi-asserted-by":"crossref","first-page":"981","DOI":"10.1126\/science.aar4362","article-title":"Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo","volume":"360","author":"Wagner","year":"2018","journal-title":"Science"},{"key":"2023062312035057900_btz864-B50","year":"2006"},{"key":"2023062312035057900_btz864-B51","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1186\/s13059-017-1382-0","article-title":"SCANPY: large-scale single-cell gene expression data analysis","volume":"19","author":"Wolf","year":"2018","journal-title":"Genome Biol"},{"key":"2023062312035057900_btz864-B52","doi-asserted-by":"crossref","first-page":"1131","DOI":"10.1038\/nsmb.2660","article-title":"Single-cell RNA-Seq profiling of human preimplantation embryos and embryonic stem cells","volume":"20","author":"Yan","year":"2013","journal-title":"Nat. Struct. Mol. Biol"},{"key":"2023062312035057900_btz864-B53","doi-asserted-by":"crossref","first-page":"e1001089","DOI":"10.1371\/journal.pgen.1001089","article-title":"A negative regulatory loop between microRNA and Hox gene controls posterior identities in Caenorhabditis elegans","volume":"6","author":"Zhao","year":"2010","journal-title":"PLoS Genet"},{"key":"2023062312035057900_btz864-B54","doi-asserted-by":"crossref","first-page":"3809","DOI":"10.1101\/gad.12.24.3809","article-title":"Reprogramming of early embryonic blastomeres into endodermal progenitors by a Caenorhabditis elegans GATA factor","volume":"12","author":"Zhu","year":"1998","journal-title":"Genes Dev"}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/academic.oup.com\/bioinformatics\/advance-article-pdf\/doi\/10.1093\/bioinformatics\/btz864\/33049096\/btz864.pdf","content-type":"application\/pdf","content-version":"am","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/13\/4047\/50671438\/bioinformatics_36_13_4047.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/36\/13\/4047\/50671438\/bioinformatics_36_13_4047.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,6,24]],"date-time":"2023-06-24T17:43:04Z","timestamp":1687628584000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/36\/13\/4047\/5682412"}},"subtitle":[],"editor":[{"given":"Alfonso","family":"Valencia","sequence":"additional","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2019,12,20]]},"references-count":54,"journal-issue":{"issue":"13","published-print":{"date-parts":[[2020,7,1]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btz864","relation":{"has-review":[{"id-type":"doi","id":"10.3410\/f.737115629.793571384","asserted-by":"object"}]},"ISSN":["1367-4803","1367-4811"],"issn-type":[{"value":"1367-4803","type":"print"},{"value":"1367-4811","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2020,7]]},"published":{"date-parts":[[2019,12,20]]}}}