Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Resource
  • Published:

A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis

Nature Methods volume 18, pages 893–902 (2021)Cite this article

Subjects

An Author Correction to this article was published on 08 December 2021

This article has been updated

Abstract

Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state–fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.

This is a preview of subscription content, access via your institution

Access options

Buy this article

  • Purchase on Springer Link
  • Instant access to full article PDF
Buy now

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: A developmental protein expression atlas describing hundreds of TFs in nearly all lineage-mapped embryonic cells.
Fig. 2: Spatiotemporal patterning of embryogenesis by lineage-, tissue- and time-specific TFs.
Fig. 3: Spatiotemporal TF cascades of cell lineage differentiation.
Fig. 4: Expression-guided discovery of new regulators of spatiotemporal fate patterning.
Fig. 5: Lineage-dependent diversity of cell states within functionally related cells.
Fig. 6: Systematic undirectional state transition during cell differentiation.

Similar content being viewed by others

Data availability

All data generated or analyzed in this study are included in the article (Extended Data files and Supplementary Tables). The 4D live imaging data generated and analyzed that support the findings of this study are deposited in BioStudies at https://www.ebi.ac.uk/biostudies/studies/S-BIAD56/. The unprocessed and processed quantitative expression of all TFs in all traced cells at each time point are deposited in Zenodo at https://doi.org/10.5281/zenodo.4737593 (ref. 163). All processed data generated in this study and their visualization are also available at http://dulab.genetics.ac.cn/TF-atlas/. Source data are provided with this paper.

Change history

References

  1. Sulston, J. E. & Horvitz, H. R. Post-embryonic cell lineages of the nematode, Caenorhabditis elegans. Dev. Biol. 56, 110–156 (1977).

    Article  CAS  PubMed  Google Scholar 

  2. Sulston, J. E., Schierenberg, E., White, J. G. & Thomson, J. N. The embryonic cell lineage of the nematode Caenorhabditis elegans. Dev. Biol. 100, 64–119 (1983).

    Article  CAS  PubMed  Google Scholar 

  3. White, J. G., Southgate, E., Thomson, J. N. & Brenner, S. The structure of the nervous system of the nematode Caenorhabditis elegans. Philos. Trans. R. Soc. Lond. B Biol. Sci. 314, 1–340 (1986).

    Article  CAS  PubMed  Google Scholar 

  4. Albertson, D. G. & Thomson, J. N. The pharynx of Caenorhabditis elegans. Philos. Trans. R. Soc. Lond. B Biol. Sci. 275, 299–325 (1976).

    Article  CAS  PubMed  Google Scholar 

  5. Marx, V. A dream of single-cell proteomics. Nat. Methods 16, 809–812 (2019).

    Article  CAS  PubMed  Google Scholar 

  6. Thul, P. J. et al. A subcellular map of the human proteome. Science 356, eaal3321 (2017).

    Article  PubMed  Google Scholar 

  7. Chalfie, M., Tu, Y., Euskirchen, G., Ward, W. W. & Prasher, D. C. Green fluorescent protein as a marker for gene expression. Science 263, 802–805 (1994).

    Article  CAS  PubMed  Google Scholar 

  8. Du, Z., Santella, A., He, F., Tiongson, M. & Bao, Z. De novo inference of systems-level mechanistic models of development from live-imaging-based phenotype analysis. Cell 156, 359–372 (2014).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Murray, J. I. et al. Automated analysis of embryonic gene expression with cellular resolution in C. elegans. Nat. Methods 5, 703–709 (2008).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Long, F., Peng, H., Liu, X., Kim, S. K. & Myers, E. A 3D digital atlas of C. elegans and its application to single-cell analyses. Nat. Methods 6, 667–672 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Liu, X. et al. Analysis of cell fate from single-cell gene expression profiles in C. elegans. Cell 139, 623–633 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Murray, J. I. et al. Multidimensional regulation of gene expression in the C. elegans embryo. Genome Res. 22, 1282–1294 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Reilly, M. B., Cros, C., Varol, E., Yemini, E. & Hobert, O. Unique homeobox codes delineate all the neuron classes of C. elegans. Nature 584, 595–601 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Sarov, M. et al. A genome-scale resource for in vivo tag-based protein function exploration in C. elegans. Cell 150, 855–866 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Araya, C. L. et al. Regulatory analysis of the C. elegans genome with spatiotemporal resolution. Nature 512, 400–405 (2014).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Mace, D. L., Weisdepp, P., Gevirtzman, L., Boyle, T. & Waterston, R. H. A high-fidelity cell lineage tracing method for obtaining systematic spatiotemporal gene expression patterns in Caenorhabditis elegans. G3 3, 851–863 (2013).

    Article  PubMed  PubMed Central  Google Scholar 

  17. Packer, J. S. et al. A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution. Science 365, eaax1971 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Boeck, M. E. et al. The time-resolved transcriptome of C. elegans. Genome Res. 26, 1441–1450 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hashimshony, T., Feder, M., Levin, M., Hall, B. K. & Yanai, I. Spatiotemporal transcriptomics reveals the evolutionary history of the endoderm germ layer. Nature 519, 219–222 (2015).

    Article  CAS  PubMed  Google Scholar 

  20. Grun, D. et al. Conservation of mRNA and protein expression during development of C. elegans. Cell Rep. 6, 565–577 (2014).

    Article  PubMed  Google Scholar 

  21. Liu, Y., Beyer, A. & Aebersold, R. On the dependency of cellular protein levels on mRNA abundance. Cell 165, 535–550 (2016).

    Article  CAS  PubMed  Google Scholar 

  22. Taniguchi, Y. et al. Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science 329, 533–538 (2010).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Labouesse, M. & Mango, S. E. Patterning the C. elegans embryo: moving beyond the cell lineage. Trends Genet. 15, 307–313 (1999).

    Article  CAS  PubMed  Google Scholar 

  24. Good, K. et al. The T-box transcription factors TBX-37 and TBX-38 link GLP-1/Notch signaling to mesoderm induction in C. elegans embryos. Development 131, 1967–1978 (2004).

    Article  CAS  PubMed  Google Scholar 

  25. Broitman-Maduro, G. et al. The NK-2 class homeodomain factor CEH-51 and the T-box factor TBX-35 have overlapping function in C. elegans mesoderm development. Development 136, 2735–2746 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Hunter, C. P. & Kenyon, C. Spatial and temporal controls target pal-1 blastomere-specification activity to a single blastomere lineage in C. elegans embryos. Cell 87, 217–226 (1996).

    Article  CAS  PubMed  Google Scholar 

  27. Mizumoto, K. & Sawa, H. Two betas or not two betas: regulation of asymmetric division by beta-catenin. Trends Cell Biol. 17, 465–473 (2007).

    Article  CAS  PubMed  Google Scholar 

  28. Lin, R., Hill, R. J. & Priess, J. R. POP-1 and anterior–posterior fate decisions in C. elegans embryos. Cell 92, 229–239 (1998).

    Article  CAS  PubMed  Google Scholar 

  29. Warner, A. D., Gevirtzman, L., Hillier, L. W., Ewing, B. & Waterston, R. H. The C. elegans embryonic transcriptome with tissue, time and alternative splicing resolution. Genome Res. 29, 1036–1045 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Kudron, M. M. et al. The ModERN Resource: genome-wide binding profiles for hundreds of Drosophila and Caenorhabditis elegans transcription factors. Genetics 208, 937–949 (2018).

    Article  CAS  PubMed  Google Scholar 

  31. Page, B. D., Zhang, W., Steward, K., Blumenthal, T. & Priess, J. R. ELT-1, a GATA-like transcription factor, is required for epidermal cell fates in Caenorhabditis elegans embryos. Genes Dev. 11, 1651–1661 (1997).

    Article  CAS  PubMed  Google Scholar 

  32. Vidal, B. et al. C. elegans SoxB genes are dispensable for embryonic neurogenesis but required for terminal differentiation of specific neuron types. Development 142, 2464–2477 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Maduro, M. F. Cell fate specification in the C. elegans embryo. Dev. Dyn. 239, 1315–1329 (2010).

    CAS  PubMed  Google Scholar 

  34. Yanai, I. et al. Pairing of competitive and topologically distinct regulatory modules enhances patterned gene expression. Mol. Syst. Biol. 4, 163 (2008).

    Article  PubMed  PubMed Central  Google Scholar 

  35. Gaudet, J. & Mango, S. E. Regulation of organogenesis by the Caenorhabditis elegans FoxA protein PHA-4. Science 295, 821–825 (2002).

    Article  CAS  PubMed  Google Scholar 

  36. Mango, S.E. The C. elegans pharynx: a model for organogenesis. WormBook, 1–26 (2007).

  37. Gilleard, J. S. & McGhee, J. D. Activation of hypodermal differentiation in the Caenorhabditis elegans embryo by GATA transcription factors ELT-1 and ELT-3. Mol. Cell. Biol. 21, 2533–2544 (2001).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Zhao, C. & Emmons, S. W. A transcription factor controlling development of peripheral sense organs in C. elegans. Nature 373, 74–78 (1995).

    Article  CAS  PubMed  Google Scholar 

  39. Fukushige, T., Brodigan, T. M., Schriefer, L. A., Waterston, R. H. & Krause, M. Defining the transcriptional redundancy of early bodywall muscle development in C. elegans: evidence for a unified theory of animal muscle development. Genes Dev. 20, 3395–3406 (2006).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Reinke, A. W., Mak, R., Troemel, E. R. & Bennett, E. J. In vivo mapping of tissue- and subcellular-specific proteomes in Caenorhabditis elegans. Sci. Adv. 3, e1602426 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  41. Gieseler, K., Qadota, H. & Benian, G. M. Development, structure, and maintenance of C. elegans body wall muscle. WormBook 2017, 1–59 (2017).

    Article  PubMed  Google Scholar 

  42. McInnes, L., Healy, J. & Melville, J. UMAP: Uniform manifold approximation and projection for dimension reduction. Preprint at https://arxiv.org/abs/1802.03426v2 (2018).

  43. Ahringer, J. Posterior patterning by the Caenorhabditis elegans even-skipped homolog vab-7. Genes Dev. 10, 1120–1130 (1996).

    Article  CAS  PubMed  Google Scholar 

  44. Tintori, S. C., Osborne Nishimura, E., Golden, P., Lieb, J. D. & Goldstein, B. A transcriptional lineage of the early C. elegans embryo. Dev. Cell 38, 430–444 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Buszczak, M., Signer, R. A. & Morrison, S. J. Cellular differences in protein synthesis regulate tissue homeostasis. Cell 159, 242–251 (2014).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Buenrostro, J. D. et al. Integrated single-cell analysis maps the continuous regulatory landscape of human hematopoietic differentiation. Cell 173, 1535–1548 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Velten, L. et al. Human haematopoietic stem cell lineage commitment is a continuous process. Nat. Cell Biol. 19, 271–281 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Yu, Y. C. et al. Preferential electrical coupling regulates neocortical lineage-dependent microcircuit assembly. Nature 486, 113–117 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Li, Y. et al. Clonally related visual cortical neurons show similar stimulus feature selectivity. Nature 486, 118–121 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Weinreb, C., Rodriguez-Fraticelli, A., Camargo, F. D. & Klein, A. M. Lineage tracing on transcriptional landscapes links state to fate during differentiation. Science 367, eaaw3381 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Reece-Hoyes, J. S. et al. A compendium of Caenorhabditis elegans regulatory transcription factors: a resource for mapping transcription regulatory networks. Genome Biol. 6, R110 (2005).

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kim, W., Underwood, R. S., Greenwald, I. & Shaye, D. D. OrthoList 2: a new comparative genomic analysis of human and Caenorhabditis elegans genes. Genetics 210, 445–461 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Harris, T. W. et al. WormBase: a modern model organism information resource. Nucleic Acids Res. 48, D762–D767 (2020).

    CAS  PubMed  Google Scholar 

  54. Li, X. et al. Systems properties and spatiotemporal regulation of cell position variability during embryogenesis. Cell Rep. 26, 313–321 (2019).

    Article  Google Scholar 

  55. Yu, B. et al. Convergent transcriptional programs regulate cAMP levels in C. elegans GABAergic motor neurons. Dev. Cell 43, 212–226 (2017).

    Article  Google Scholar 

  56. Kerk, S. Y., Kratsios, P., Hart, M., Mourao, R. & Hobert, O. Diversification of C. elegans motor neuron identity via selective effector gene repression. Neuron 93, 80–98 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Paix, A. et al. Scalable and versatile genome editing using linear DNAs with microhomology to Cas9 sites in Caenorhabditis elegans. Genetics 198, 1347–1356 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  58. Lloret-Fernandez, C. et al. A transcription factor collective defines the HSN serotonergic neuron regulatory landscape. eLife 7, e32785 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  59. Walton, T. et al. The Bicoid class homeodomain factors ceh-36/OTX and unc-30/PITX cooperate in C. elegans embryonic progenitor cells to regulate robust development. PLoS Genet. 11, e1005003 (2015).

    Article  PubMed  PubMed Central  Google Scholar 

  60. Hostettler, L. et al. The bright fluorescent protein mNeonGreen facilitates protein expression analysis in vivo. G3 7, 607–615 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Dickinson, D. J., Ward, J. D., Reiner, D. J. & Goldstein, B. Engineering the Caenorhabditis elegans genome using Cas9-triggered homologous recombination. Nat. Methods 10, 1028–1034 (2013).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bao, Z. & Murray, J. I. Mounting Caenorhabditis elegans embryos for live imaging of embryogenesis. Cold Spring Harb. Protoc. 2011, pdb.prot065599 (2011).

    Article  PubMed  Google Scholar 

  63. Tenenhaus, C., Schubert, C. & Seydoux, G. Genetic requirements for PIE-1 localization and inhibition of gene expression in the embryonic germ lineage of Caenorhabditis elegans. Dev. Biol. 200, 212–224 (1998).

    Article  CAS  PubMed  Google Scholar 

  64. Mello, C. C. et al. The PIE-1 protein and germline specification in C. elegans embryos. Nature 382, 710–712 (1996).

    Article  CAS  PubMed  Google Scholar 

  65. Santella, A., Du, Z. & Bao, Z. A semi-local neighborhood-based framework for probabilistic cell lineage tracing. BMC Bioinformatics 15, 217 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  66. Santella, A., Du, Z., Nowotschin, S., Hadjantonakis, A. K. & Bao, Z. A hybrid blob-slice model for accurate and efficient detection of fluorescence-labeled nuclei in 3D. BMC Bioinformatics 11, 580 (2010).

    Article  PubMed  PubMed Central  Google Scholar 

  67. Katzman, B., Tang, D., Santella, A. & Bao, Z. AceTree: a major update and case study in the long-term maintenance of open-source scientific software. BMC Bioinformatics 19, 121 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  68. Du, Z. et al. The regulatory landscape of lineage differentiation in a metazoan embryo. Dev. Cell 34, 592–607 (2015).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Kervrann, C., Legland, D. & Pardini, L. Robust incremental compensation of the light attenuation with depth in 3D fluorescence microscopy. J. Microsc. 214, 297–314 (2004).

    Article  CAS  PubMed  Google Scholar 

  70. Bao, Z., Zhao, Z., Boyle, T. J., Murray, J. I. & Waterston, R. H. Control of cell cycle timing during C. elegans embryogenesis. Dev. Biol. 318, 65–72 (2008).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Moore, J. L., Du, Z. & Bao, Z. Systematic quantification of developmental phenotypes at single-cell resolution during embryogenesis. Development 140, 3266–3274 (2013).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Qin, H. & Powell-Coffman, J. A. The Caenorhabditis elegans aryl hydrocarbon receptor, AHR-1, regulates neuronal development. Dev. Biol. 270, 64–75 (2004).

    Article  CAS  PubMed  Google Scholar 

  73. Ceol, C. J. & Horvitz, H. R. dpl-1 DP and efl-1 E2F act with lin-35 Rb to antagonize Ras signaling in C. elegans vulval development. Mol. Cell 7, 461–473 (2001).

    Article  CAS  PubMed  Google Scholar 

  74. Sherwood, D. R., Butler, J. A., Kramer, J. M. & Sternberg, P. W. FOS-1 promotes basement-membrane removal during anchor-cell invasion in C. elegans. Cell 121, 951–962 (2005).

    Article  CAS  PubMed  Google Scholar 

  75. Hope, I. A. PES-1 is expressed during early embryogenesis in Caenorhabditis elegans and has homology to the fork head family of transcription factors. Development 120, 505–514 (1994).

    Article  CAS  PubMed  Google Scholar 

  76. Ray, P., Schnabel, R. & Okkema, P. G. Behavioral and synaptic defects in C. elegans lacking the NK-2 homeobox gene ceh-28. Dev. Neurobiol. 68, 421–433 (2008).

    Article  CAS  PubMed  Google Scholar 

  77. Huang, T. F. et al. BLMP-1/Blimp-1 regulates the spatiotemporal cell migration pattern in C. elegans. PLoS Genet. 10, e1004428 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  78. Chen, L., Krause, M., Draper, B., Weintraub, H. & Fire, A. Body-wall muscle formation in Caenorhabditis elegans embryos that lack the MyoD homolog hlh-1. Science 256, 240–243 (1992).

    Article  CAS  PubMed  Google Scholar 

  79. Krause, M. et al. A C. elegans E/daughterless bHLH protein marks neuronal but not striated muscle development. Development 124, 2179–2189 (1997).

    Article  CAS  PubMed  Google Scholar 

  80. Jia, Y., Xie, G., McDermott, J. B. & Aamodt, E. The C. elegans gene pag-3 is homologous to the zinc finger proto-oncogene gfi-1. Development 124, 2063–2073 (1997).

    Article  CAS  PubMed  Google Scholar 

  81. Mathies, L. D., Henderson, S. T. & Kimble, J. The C. elegans hand gene controls embryogenesis and early gonadogenesis. Development 130, 2881–2892 (2003).

    Article  CAS  PubMed  Google Scholar 

  82. Baran, R., Aronoff, R. & Garriga, G. The C. elegans homeodomain gene unc-42 regulates chemosensory and glutamate receptor expression. Development 126, 2241–2251 (1999).

    Article  CAS  PubMed  Google Scholar 

  83. Yuan, J., Tirabassi, R. S., Bush, A. B. & Cole, M. D. The C. elegans MDL-1 and MXL-1 proteins can functionally substitute for vertebrate MAD and MAX. Oncogene 17, 1109–1118 (1998).

    Article  CAS  PubMed  Google Scholar 

  84. Pickett, C. L., Breen, K. T. & Ayer, D. E. A C. elegans Myc-like network cooperates with semaphorin and Wnt signaling pathways to control cell migration. Dev. Biol. 310, 226–239 (2007).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Hallam, S., Singer, E., Waring, D. & Jin, Y. The C. elegans NeuroD homolog cnd-1 functions in multiple aspects of motor neuron fate specification. Development 127, 4239–4252 (2000).

    Article  CAS  PubMed  Google Scholar 

  86. Svendsen, P. C. & McGhee, J. D. The C. elegans neuronally expressed homeobox gene ceh-10 is closely related to genes expressed in the vertebrate eye. Development 121, 1253–1262 (1995).

    Article  CAS  PubMed  Google Scholar 

  87. Clark, S. G. & Chiu, C. C. elegans ZAG-1, a Zn-finger-homeodomain protein, regulates axonal development and neuronal differentiation. Development 130, 3781–3794 (2003).

    Article  CAS  PubMed  Google Scholar 

  88. Aspock, G., Ruvkun, G. & Burglin, T. R. The Caenorhabditis elegans ems class homeobox gene ceh-2 is required for M3 pharynx motoneuron function. Development 130, 3369–3378 (2003).

    Article  PubMed  Google Scholar 

  89. Woollard, A. & Hodgkin, J. The Caenorhabditis elegans fate-determining gene mab-9 encodes a T-box protein required to pattern the posterior hindgut. Genes Dev. 14, 596–603 (2000).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Hobert, O., Tessmar, K. & Ruvkun, G. The Caenorhabditis elegans lim-6 LIM homeobox gene regulates neurite outgrowth and function of particular GABAergic neurons. Development 126, 1547–1562 (1999).

    Article  CAS  PubMed  Google Scholar 

  91. Labouesse, M., Hartwieg, E. & Horvitz, H. R. The Caenorhabditis elegans LIN-26 protein is required to specify and/or maintain all non-neuronal ectodermal cell fates. Development 122, 2579–2588 (1996).

    Article  CAS  PubMed  Google Scholar 

  92. Okkema, P. G. & Fire, A. The Caenorhabditis elegans NK-2 class homeoprotein CEH-22 is involved in combinatorial activation of gene expression in pharyngeal muscle. Development 120, 2175–2186 (1994).

    Article  CAS  PubMed  Google Scholar 

  93. Dozier, C., Kagoshima, H., Niklaus, G., Cassata, G. & Burglin, T. R. The Caenorhabditis elegans Six/sine oculis class homeobox gene ceh-32 is required for head morphogenesis. Dev. Biol. 236, 289–303 (2001).

    Article  CAS  PubMed  Google Scholar 

  94. Andachi, Y. Caenorhabditis elegans T-box genes tbx-9 and tbx-8 are required for formation of hypodermis and body-wall muscle in embryogenesis. Genes Cells 9, 331–344 (2004).

    Article  CAS  PubMed  Google Scholar 

  95. Cassata, G. et al. ceh-16/engrailed patterns the embryonic epidermis of Caenorhabditis elegans. Development 132, 739–749 (2005).

    Article  CAS  PubMed  Google Scholar 

  96. Zheng, X., Chung, S., Tanabe, T. & Sze, J. Y. Cell-type-specific regulation of serotonergic identity by the C. elegans LIM-homeodomain factor LIM-4. Dev. Biol. 286, 618–628 (2005).

    Article  CAS  PubMed  Google Scholar 

  97. Krause, M., Fire, A., Harrison, S. W., Priess, J. & Weintraub, H. CeMyoD accumulation defines the body wall muscle cell fate during C. elegans embryogenesis. Cell 63, 907–919 (1990).

    Article  CAS  PubMed  Google Scholar 

  98. Kostrouchova, M., Krause, M., Kostrouch, Z. & Rall, J. E. CHR3: a Caenorhabditis elegans orphan nuclear hormone receptor required for proper epidermal development and molting. Development 125, 1617–1626 (1998).

    Article  CAS  PubMed  Google Scholar 

  99. Haynes, C. M., Petrova, K., Benedetti, C., Yang, Y. & Ron, D. ClpP mediates activation of a mitochondrial unfolded protein response in C. elegans. Dev. Cell 13, 467–480 (2007).

    Article  CAS  PubMed  Google Scholar 

  100. Asahina, M. et al. The conserved nuclear receptor Ftz-F1 is required for embryogenesis, moulting and reproduction in Caenorhabditis elegans. Genes Cells 5, 711–723 (2000).

    Article  CAS  PubMed  Google Scholar 

  101. Fakhouri, T. H., Stevenson, J., Chisholm, A. D. & Mango, S. E. Dynamic chromatin organization during foregut development mediated by the organ selector gene PHA-4/FoxA. PLoS Genet. 6, e1001060 (2010).

    Article  PubMed  PubMed Central  Google Scholar 

  102. Gilleard, J. S., Shafi, Y., Barry, J. D. & McGhee, J. D. ELT-3: a Caenorhabditis elegans GATA factor expressed in the embryonic epidermis during morphogenesis. Dev. Biol. 208, 265–280 (1999).

    Article  CAS  PubMed  Google Scholar 

  103. Koh, K. & Rothman, J. H. ELT-5 and ELT-6 are required continuously to regulate epidermal seam cell differentiation and cell fusion in C. elegans. Development 128, 2867–2880 (2001).

    Article  CAS  PubMed  Google Scholar 

  104. Zhu, J. W. et al. end-1 encodes an apparent GATA factor that specifies the endoderm precursor in Caenorhabditis elegans embryos. Genes Dev. 11, 2883–2896 (1997).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  105. Suggs, B. Z., Latham, A. L., Dawes, A. T. & Chamberlin, H. M. FACT complex gene duplicates exhibit redundant and non-redundant functions in C. elegans. Dev. Biol. 444, 71–82 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Azzaria, M., Goszczynski, B., Chung, M. A., Kalb, J. M. & McGhee, J. D. A fork head/HNF-3 homolog expressed in the pharynx and intestine of the Caenorhabditis elegans embryo. Dev. Biol. 178, 289–303 (1996).

    Article  CAS  PubMed  Google Scholar 

  107. Hope, I. A., Mounsey, A., Bauer, P. & Aslam, S. The forkhead gene family of Caenorhabditis elegans. Gene 304, 43–55 (2003).

    Article  CAS  PubMed  Google Scholar 

  108. Hobert, O. & Westphal, H. Functions of LIM-homeobox genes. Trends Genet. 16, 75–83 (2000).

    Article  CAS  PubMed  Google Scholar 

  109. Fukushige, T., Hawkins, M. G. & McGhee, J. D. The GATA-factor elt-2 is essential for formation of the Caenorhabditis elegans intestine. Dev. Biol. 198, 286–302 (1998).

    CAS  PubMed  Google Scholar 

  110. Maduro, M. F. et al. Genetic redundancy in endoderm specification within the genus Caenorhabditis. Dev. Biol. 284, 509–522 (2005).

    Article  CAS  PubMed  Google Scholar 

  111. Ratnappan, R. et al. Germline signals deploy NHR-49 to modulate fatty-acid beta-oxidation and desaturation in somatic tissues of C. elegans. PLoS Genet. 10, e1004829 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  112. Pflugrad, A., Meir, J. Y., Barnes, T. M. & Miller, D. M. 3rd The Groucho-like transcription factor UNC-37 functions with the neural specificity gene unc-4 to govern motor neuron identity in C. elegans. Development 124, 1699–1709 (1997).

    Article  CAS  PubMed  Google Scholar 

  113. Pujol, N., Torregrossa, P., Ewbank, J. J. & Brunet, J. F. The homeodomain protein CePHOX2/CEH-17 controls anteroposterior axonal growth in C. elegans. Development 127, 3361–3371 (2000).

    Article  CAS  PubMed  Google Scholar 

  114. Reece-Hoyes, J. S. et al. Insight into transcription factor gene duplication from Caenorhabditis elegans promoterome-driven expression patterns. BMC Genomics 8, 27 (2007).

    Article  PubMed  PubMed Central  Google Scholar 

  115. Cassata, G. et al. The LIM homeobox gene ceh-14 confers thermosensory function to the AFD neurons in Caenorhabditis elegans. Neuron 25, 587–597 (2000).

    Article  CAS  PubMed  Google Scholar 

  116. Nimmo, R., Antebi, A. & Woollard, A. mab-2 encodes RNT-1, a C. elegans Runx homologue essential for controlling cell proliferation in a stem cell-like developmental lineage. Development 132, 5043–5054 (2005).

    Article  CAS  PubMed  Google Scholar 

  117. Maduro, M. F. & Rothman, J. H. Making worm guts: the gene regulatory network of the Caenorhabditis elegans endoderm. Dev. Biol. 246, 68–85 (2002).

    Article  CAS  PubMed  Google Scholar 

  118. Bowerman, B., Draper, B. W., Mello, C. C. & Priess, J. R. The maternal gene skn-1 encodes a protein that is distributed unequally in early C. elegans embryos. Cell 74, 443–452 (1993).

    Article  CAS  PubMed  Google Scholar 

  119. Unhavaithaya, Y. et al. MEP-1 and a homolog of the NURD complex component Mi-2 act together to maintain germline-soma distinctions in C. elegans. Cell 111, 991–1002 (2002).

    Article  CAS  PubMed  Google Scholar 

  120. Lee, C. S., Lu, T. & Seydoux, G. Nanos promotes epigenetic reprograming of the germline by down-regulation of the THAP transcription factor LIN-15B. eLife 6, e30201 (2017).

    Article  PubMed  PubMed Central  Google Scholar 

  121. Gissendanner, C. R. & Sluder, A. E. nhr-25, the Caenorhabditis elegans ortholog of ftz-f1, is required for epidermal and somatic gonad development. Dev. Biol. 221, 259–272 (2000).

    Article  CAS  PubMed  Google Scholar 

  122. Broitman-Maduro, G., Maduro, M. F. & Rothman, J. H. The noncanonical binding site of the MED-1 GATA factor defines differentially regulated target genes in the C. elegans mesendoderm. Dev. Cell 8, 427–433 (2005).

    Article  CAS  PubMed  Google Scholar 

  123. Johnston, R. J. Jr. & Hobert, O. A novel C. elegans zinc finger transcription factor, lsy-2, required for the cell-type-specific expression of the lsy-6 microRNA. Development 132, 5451–5460 (2005).

    Article  CAS  PubMed  Google Scholar 

  124. Buckley, M. S., Chau, J., Hoppe, P. E. & Coulter, D. E. odd-skipped homologs function during gut development in C. elegans. Dev. Genes Evol. 214, 10–18 (2004).

    Article  PubMed  Google Scholar 

  125. Thompson, K. W. et al. The paired-box protein PAX-3 regulates the choice between lateral and ventral epidermal cell fates in C. elegans. Dev. Biol. 412, 191–207 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  126. Chisholm, A. D. & Horvitz, H. R. Patterning of the Caenorhabditis elegans head region by the Pax-6 family member vab-3. Nature 377, 52–55 (1995).

    Article  CAS  PubMed  Google Scholar 

  127. Morck, C., Rauthan, M., Wagberg, F. & Pilon, M. pha-2 encodes the C. elegans ortholog of the homeodomain protein HEX and is required for the formation of the pharyngeal isthmus. Dev. Biol. 272, 403–418 (2004).

    Article  CAS  PubMed  Google Scholar 

  128. Kalb, J. M. et al. pha-4 is Ce-fkh-1, a fork head/HNF-3 alpha,beta,gamma homolog that functions in organogenesis of the C. elegans pharynx. Development 125, 2171–2180 (1998).

    Article  CAS  PubMed  Google Scholar 

  129. Horner, M. A. et al. pha-4, an HNF-3 homolog, specifies pharyngeal organ identity in Caenorhabditis elegans. Genes Dev. 12, 1947–1952 (1998).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  130. Tepper, R. G. et al. PQM-1 complements DAF-16 as a key transcriptional regulator of DAF-2-mediated development and longevity. Cell 154, 676–690 (2013).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  131. Cassata, G. et al. Rapid expression screening of Caenorhabditis elegans homeobox open reading frames using a two-step polymerase chain reaction promoter-GFP reporter construction technique. Gene 212, 127–135 (1998).

    Article  CAS  PubMed  Google Scholar 

  132. Neves, A. & Priess, J. R. The REF-1 family of bHLH transcription factors pattern C. elegans embryos through Notch-dependent and Notch-independent pathways. Dev. Cell 8, 867–879 (2005).

    Article  CAS  PubMed  Google Scholar 

  133. Burglin, T. R. & Ruvkun, G. Regulation of ectodermal and excretory function by the C. elegans POU homeobox gene ceh-6. Development 128, 779–790 (2001).

    Article  CAS  PubMed  Google Scholar 

  134. Pocock, R., Ahringer, J., Mitsch, M., Maxwell, S. & Woollard, A. A regulatory network of T-box genes and the even-skipped homologue vab-7 controls patterning and morphogenesis in C. elegans. Development 131, 2373–2385 (2004).

    Article  CAS  PubMed  Google Scholar 

  135. Zhu, J. W., Fukushige, T., McGhee, J. D. & Rothman, J. H. Reprogramming of early embryonic blastomeres into endodermal progenitors by a Caenorhabditis elegans GATA factor. Genes Dev. 12, 3809–3814 (1998).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  136. Maduro, M. F., Meneghini, M. D., Bowerman, B., Broitman-Maduro, G. & Rothman, J. H. Restriction of mesendoderm to a single blastomere by the combined action of SKN-1 and a GSK-3beta homolog is mediated by MED-1 and -2 in C. elegans. Mol. Cell 7, 475–485 (2001).

    Article  CAS  PubMed  Google Scholar 

  137. Doyle, T. G., Wen, C. & Greenwald, I. SEL-8, a nuclear protein required for LIN-12 and GLP-1 signaling in Caenorhabditis elegans. Proc. Natl Acad. Sci. USA 97, 7877–7881 (2000).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  138. Roy Chowdhuri, S., Crum, T., Woollard, A., Aslam, S. & Okkema, P. G. The T-box factor TBX-2 and the SUMO conjugating enzyme UBC-9 are required for ABa-derived pharyngeal muscle in C. elegans. Dev. Biol. 295, 664–677 (2006).

    Article  PubMed  Google Scholar 

  139. Deng, H. et al. Transcription factor NFY globally represses the expression of the C. elegans Hox gene abdominal-B homolog egl-5. Dev. Biol. 308, 583–592 (2007).

    Article  CAS  PubMed  Google Scholar 

  140. Finney, M. & Ruvkun, G. The unc-86 gene product couples cell lineage and cell identity in C. elegans. Cell 63, 895–905 (1990).

    Article  CAS  PubMed  Google Scholar 

  141. Masoudi, N. et al. Unconventional function of an Achaete-Scute homolog as a terminal selector of nociceptive neuron identity. PLoS Biol. 16, e2004979 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

  142. Wacker, I., Schwarz, V., Hedgecock, E. M. & Hutter, H. zag-1, a Zn-finger homeodomain transcription factor controlling neuronal differentiation and axon outgrowth in C. elegans. Development 130, 3795–3805 (2003).

    Article  CAS  PubMed  Google Scholar 

  143. Brunschwig, K. et al. Anterior organization of the Caenorhabditis elegans embryo by the labial-like Hox gene ceh-13. Development 126, 1537–1546 (1999).

    Article  CAS  PubMed  Google Scholar 

  144. Gladden, J. M. & Meyer, B. J. A ONECUT homeodomain protein communicates X chromosome dose to specify Caenorhabditis elegans sexual fate by repressing a sex switch gene. Genetics 177, 1621–1637 (2007).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  145. Jiang, Y., Horner, V. & Liu, J. The HMX homeodomain protein MLS-2 regulates cleavage orientation, cell proliferation and cell fate specification in the C. elegans postembryonic mesoderm. Development 132, 4119–4130 (2005).

    Article  CAS  PubMed  Google Scholar 

  146. Powell, J. R., Jow, M. M. & Meyer, B. J. The T-box transcription factor SEA-1 is an autosomal element of the X:A signal that determines C. elegans sex. Dev. Cell 9, 339–349 (2005).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  147. Prasad, B. C. et al. unc-3, a gene required for axonal guidance in Caenorhabditis elegans, encodes a member of the O/E family of transcription factors. Development 125, 1561–1568 (1998).

    Article  CAS  PubMed  Google Scholar 

  148. Sarafi-Reinach, T. R. & Sengupta, P. The forkhead domain gene unc-130 generates chemosensory neuron diversity in C. elegans. Genes Dev. 14, 2472–2485 (2000).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  149. Lee, J. et al. A Myt1 family transcription factor defines neuronal fate by repressing non-neuronal genes. eLife 8, e46703 (2019).

    Article  PubMed  PubMed Central  Google Scholar 

  150. Sarin, S., Antonio, C., Tursun, B. & Hobert, O. The C. elegans tailless/TLX transcription factor nhr-67 controls neuronal identity and left/right asymmetric fate diversification. Development 136, 2933–2944 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  151. Maduro, M. F. et al. MED GATA factors promote robust development of the C. elegans endoderm. Dev. Biol. 404, 66–79 (2015).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  152. Miyahara, K., Suzuki, N., Ishihara, T., Tsuchiya, E. & Katsura, I. TBX2/TBX3 transcriptional factor homologue controls olfactory adaptation in Caenorhabditis elegans. J. Neurobiol. 58, 392–402 (2004).

    Article  CAS  PubMed  Google Scholar 

  153. Raharjo, W. H., Ghai, V., Dineen, A., Bastiani, M. & Gaudet, J. Cell architecture: surrounding muscle cells shape gland cell morphology in the Caenorhabditis elegans pharynx. Genetics 189, 885–897 (2011).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  154. Chen, L., Krause, M., Sepanski, M. & Fire, A. The Caenorhabditis elegans MYOD homologue HLH-1 is essential for proper muscle function and complete morphogenesis. Development 120, 1631–1641 (1994).

    Article  CAS  PubMed  Google Scholar 

  155. Dowen, R. H., Breen, P. C., Tullius, T., Conery, A. L. & Ruvkun, G. A microRNA program in the C. elegans hypodermis couples to intestinal mTORC2/PQM-1 signaling to modulate fat transport. Genes Dev. 30, 1515–1528 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  156. Edgar, L. G., Carr, S., Wang, H. & Wood, W. B. Zygotic expression of the caudal homolog pal-1 is required for posterior patterning in Caenorhabditis elegans embryogenesis. Dev. Biol. 229, 71–88 (2001).

    Article  CAS  PubMed  Google Scholar 

  157. Kamath, R. S., Martinez-Campos, M., Zipperlen, P., Fraser, A. G. & Ahringer, J. Effectiveness of specific RNA-mediated interference through ingested double-stranded RNA in Caenorhabditis elegans. Genome Biol. 2, RESEARCH0002 (2001).

    CAS  PubMed  Google Scholar 

  158. Moskowitz, I. P. & Rothman, J. H. lin-12 and glp-1 are required zygotically for early embryonic cellular interactions and are regulated by maternal GLP-1 signaling in Caenorhabditis elegans. Development 122, 4105–4117 (1996).

    Article  CAS  PubMed  Google Scholar 

  159. Bowerman, B., Eaton, B. A. & Priess, J. R. skn-1, a maternally expressed gene required to specify the fate of ventral blastomeres in the early C. elegans embryo. Cell 68, 1061–1075 (1992).

    Article  CAS  PubMed  Google Scholar 

  160. Draper, B. W., Mello, C. C., Bowerman, B., Hardin, J. & Priess, J. R. MEX-3 is a KH domain protein that regulates blastomere identity in early C. elegans embryos. Cell 87, 205–216 (1996).

    Article  CAS  PubMed  Google Scholar 

  161. Schubert, C. M., Lin, R., de Vries, C. J., Plasterk, R. H. & Priess, J. R. MEX-5 and MEX-6 function to establish soma/germline asymmetry in early C. elegans embryos. Mol. Cell 5, 671–682 (2000).

    Article  CAS  PubMed  Google Scholar 

  162. Angeles-Albores, D., Lee, R., Chan, J. & Sternberg, P. Two new functions in the WormBase Enrichment Suite. MicroPubl. Biol. 2018 (2018).

  163. Zhao, Z., Du, Z. A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis. Zenodo https://doi.org/10.5281/zenodo.4737593 (2021).

  164. Smith, P. A. & Mango, S. E. Role of T-box gene tbx-2 for anterior foregut muscle development in C. elegans. Dev. Biol. 302, 25–39 (2007).

    Article  CAS  PubMed  Google Scholar 

  165. Van Auken, K. et al. Roles of the Homothorax/Meis/Prep homolog UNC-62 and the Exd/Pbx homologs CEH-20 and CEH-40 in C. elegans embryogenesis. Development 129, 5255–5268 (2002).

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We thank M. -Q. Dong, National Institute of Biological Sciences, for providing CRISPR–Cas9 reagents, and Caenorhabditis Genetics Center for providing some strains. This work was supported by grants from the ‘Strategic Priority Research Program’ of the Chinese Academy of Sciences to Z.D. (XDB19000000), the National Natural Science Foundation of China to Z.D. (31771598, 32061143010 and 31722035) and to X.M. (31900578), and the State Key Laboratory of Molecular Developmental Biology, China to Z.D.

Author information

Author notes

  1. These authors contributed equally: Xuehua Ma, Zhiguang Zhao, Long Xiao.

Authors and Affiliations

  1. State Key Laboratory of Molecular Developmental Biology, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing, China

    Xuehua Ma, Zhiguang Zhao, Long Xiao, Weina Xu, Yahui Kou, Yangyang Wang & Zhuo Du

  2. University of Chinese Academy of Sciences, Beijing, China

    Zhiguang Zhao, Long Xiao, Weina Xu, Yahui Kou & Zhuo Du

  3. National Institute of Biological Sciences, Beijing, China

    Yanping Zhang & Gang Wu

Authors
  1. Xuehua Ma
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhiguang Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Long Xiao
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Weina Xu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yahui Kou
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yanping Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Gang Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Yangyang Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Zhuo Du
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Z.D., X.M. and Z.Z. conceived the project and designed the study; X.M., L.X., Z.D., W.X., Y.K., Y.Z., G.W. and Y.W. conducted the experiments and generated the data. Z.Z. and Z.D. performed the data analysis. Z.D. wrote the manuscript with input from all authors.

Corresponding author

Correspondence to Zhuo Du.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Methods thanks David Matus, John I. Murray, Marta Shahbazi and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Madhura Mukhopadhyay was the primary editor on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Generation of protein-fusion fluorescent reporter of TFs.

a, Selection of 290 high-confidence TFs as target genes with human homologs or worm-specific TFs with potential interest in cell differentiation. Right: TF classes based on typical domains. Classes with less than three TFs are not included. b, Maximum projections of 3D images showing cellular protein expression of representative TFs. 4D imaging of each strain was repeated in ≥2 embryos, showing similar expression patterns. Scale bar, 10 µm. c, Fosmid-based protein-fusion reporters. Two hundred two protein-fusion reporters of 185 TFs have been generated by the ModENCODE projects in which GFP (green) was fused to a TF in a fosmid by recombineering, and the fosmid was then integrated into the C. elegans genome14,15. d, Tagging endogenous TFs with mNeonGreen/GFP by CRISPR/Cas9-mediated gene knock-in. A homologous recombination (HR) template sequence containing the coding sequence of mNeonGreen/GFP (green) flanked by 800-1,000 bp sequence homologous to the endogenous TF locus was used to guide HR following cleavage by Cas9 at the N- or C-terminals of the TF gene. (e) Flowchart showing the sources of 291 protein-fusion reporters for 266 TFs.

Extended Data Fig. 2 High quality of cellular protein expression of TFs.

a, Distribution of the Pearson correlation coefficient (R) for experimental replicates of TF expression assays. b, Representative examples showing the consistency of cellular TF expression between replicates. c, Comparison of cellular TF expression between different fosmid-based protein-fusion reporter strains that tag the same TF. d, Comparison of cellular TF expression between fosmid-based (top) and knockin-based (bottom) reporter strains that tag the same TF. e, Comparison of the measured TF expression to that previously reported for 30 TFs selected from the literature in cell tracks (a series of temporally ordered mother-daughter cells) leading to terminal cells. f, Pie charts comparing the distribution of the status of lineage identity (ID) assignment of all embryonic cells (n = 1341) between scRNA-seq and protein expression data. g, Comparison of the precision and sensitivity of mRNA and protein expression results at different expression cut-offs. The true positive (TP), false-positive (FP), and false-negative (FN) rates are calculated to determine the sensitivity and precision as TP/(TP + FN) and TP/(TP + FP), respectively. h, Correlation between mRNA and protein expression in identical cells. ρ and P-value were calculated by a two-sided Spearman’s rank correlation, n = 253,159.

Source data

Extended Data Fig. 3 Spatiotemporal patterning of cell states by lineage-, tissue-, and time-specific TFs.

a, Tree visualization of binarized expression (green) of representative lineage-specific TFs in lineaged cells. b, Representative example showing single-cell expression of three anterior-posterior asymmetric TFs. Purple and blue bars shown on the bottom of the cell lineage trees indicate the lineages produced by anteriorly and posteriorly localized progenitor cells following a division, respectively. c, Enrichment of the overlap between tissue-specific TFs identified in this study and by29. Statistics: One-tailed hypergeometric test. d, Protein expression of representative tissue-specific TFs in lineaged cells. Barcodes on the bottom highlight cells of specific tissue types. e, Tree visualization of the cellular expression patterns of TFs before (top) and after (bottom) switching progenitor cell fate by RNAi against specific fate specifiers. Arrows indicate fate transformations. Similar results were observed in ≥2 embryos. f, Relative enrichment of tissue-specific genes in the targets of tissue-specific TFs. In each tissue, relative enrichment was measured as the ratio of the frequency of genes specifically expressed in this tissue to that of genes specifically expressed in other tissues. Tissue-specific genes identified by29 were used. g, Tree visualization of the number of transiently-expressed TFs in each lineaged cell. The instances in which a TF is expressed before degradation were used to calculate the number. h, Left: single-cell expression of two transiently-expressed TFs. Right: changes in TF expression levels before and after degradation in indicated cells. i, Heatmap showing the contribution of transient expression to the net difference in TF expression between mother and daughter cells (Mother-daughter cell pair number n = 1,200). j, Percent of distinguishable mother-daughter states under various criteria before and after considering transient TF expression. Statistics: Two-tailed Wilcoxon signed-rank test. k, Distribution of the fraction of specifically expressed TFs in all TFs that distinguish each cell from individual other cells. l, Each curve shows the fractions of cells (averaged from 1000 times of simulations) whose TF expression is distinct from all other cells under different definitions (number of TFs showing differential binary expression) after removing a certain number of specifically-expressed TFs.

Source data

Extended Data Fig. 4 Construction of spatiotemporal TF cascades.

a, Tree visualization of the cell lineage leading to corresponding progenitors (colored circles) of defined tissue types. b, Spatiotemporal modules of cell lineage differentiation. Each column represents the development of a tissue progenitor cell lineage (spatial module). Each row represents the classification of the regulation of each progenitor cell lineage into three temporal modules (color gradient). c, Schematic of calculating the similarity in single-cell expression of a given gene. d, Single-cell expression similarity between 30 pairs of TFs known to function in the same pathway24,34,127,132,151,164,165. e, Comparison of expression similarity between genes in the same pathway (n = 30) and between randomly selected gene pairs (n = 30). The boxplot shows the median with IQR and whiskers extending to 1.5× the IQR; outliers are plotted as points. Statistics: Two-tailed Wilcoxon signed-rank test. f, Tree visualization of cellular expression of PHA-4 and CEH-27 in each lineaged cell. The barcode at the bottom indicates the pharyngeal cells. g, Changes in cellular CEH-27 protein expression in pha-4(zu225) embryos. Above showing cellular expression of CEH-27 protein in lineaged cells of wild-type embryos with vertical lines on the bottom indicating cells in which loss or reduction of CEH-27 expression was detected in pha-4(zu225) embryos. Numbers indicate the penetrance of expression changes.

Extended Data Fig. 5 ELT-1 regulates LIN-32 expression in multiple cell lineages.

a, Single-cell protein expression of ELT-1 and LIN-32. b, Top: Neuronal-specific TFs targeted by ELT-1. Bottom: ELT-1 binding pattern near the lin-32 gene. c, Comparison of single-cell mNeonGreen::LIN-32 presence in all analyzed elt-1(ok1002)/+ and elt-1(ok1002) embryos in lineaged cells. d, Representative micrographs showing loss and gain of mNeonGreen::LIN-32 expression and in cells from various lineages. Similar results were observed in ≥3 embryos. Scale bar, 2 μm. The center Z planes of corresponding cells were identified using mCherry signals and cropped from embryos. Figure organization is as in Fig. 4b. e, Representative micrographs showing loss and gain of cell death in elt-1(ok1002) embryos. Similar results were observed in ≥4 embryos. Scale bar, 2 μm. Programmed cell death was determined by a characteristic sequence of morphological changes in the nucleus of the corresponding cell. The center Z planes of corresponding cells were identified using mCherry signals and cropped from embryos. Figure organization is as in Fig. 4d.

Extended Data Fig. 6 ELT-1 regulates SOX-2 expression in multiple cell lineages.

a, Single-cell protein expression of ELT-1 and SOX-2::mNeonGreen. The color-coded barcode on the bottom highlights the neuronal and skin cells. b, Left: loss and gain of SOX-2:: mNeonGreen expression in corresponding cell lineages of elt-1(ok1002) embryos. Similar results were observed in ≥2 embryos. Numbers at the bottom indicate the penetrance of each phenotype. c, Comparison of single-cell SOX-2::mNeonGreen presence in all analyzed elt-1(ok1002)/+ and elt-1(ok1002) embryos in corresponding cell lineages.

Extended Data Fig. 7 M03D4.4 converges muscle differentiation in the body wall and pharynx.

a, Pan-muscle specificity of M03D4.4. Left: cell types of M03D4.4-expressing cells (apoptotic cells were not considered). Right: coverage of M03D4.4-expressing cells in different muscle cell types. b, Single-cell protein expression of M03D4.4::GFP in cell lineages before and after switching fates of specific progenitor cells. The homeotic transformation of ABala, ABalp, ABara, and ABarp cell fates to C (ABa-to-4C) was induced by mex-3(RNAi). The homeotic transformation of ABa cell fate to EMS (ABa-to-EMS) was induced by mex-5(RNAi). Similar results were observed in ≥2 embryos. c, Comparison of single-cell protein expression of M03D4.4 and other known regulatory TFs (PAL-1, HND-1, UNC-120, HLH-1) in the C and D lineages. d, Comparison of single-cell protein expression of PHA-4, TBX-2, and M03D4.4 in ABalpa, ABara, and MS lineages that produce pharyngeal muscle cells. For c, d, similar results were observed in ≥2 embryos for each reporter strain. e, Inferred regulatory relationship between M03D4.4 and other regulators. f, Comparison of single-cell M03D4.4 protein expression in corresponding cell lineages between all analyzed wild-type and pal-1(RNAi) embryos. g, Comparison of single-cell M03D4.4 protein expression in corresponding cell lineages between all analyzed wild-type and pha-4(zu225) embryos. h, Enrichment of proteins preferentially expressed in pharyngeal and body muscle cells40 in the ChIP-seq targets of M03D4.4. Statistics: One-tailed hypergeometric test. i, Enrichment of genes involved in body muscle development and function41 in the ChIP-seq targets of M03D4.4. Statistics: One-tailed hypergeometric test. j, Enrichment of RNAi phenotypes in M03D4.4 target genes. Statistics: One-tailed hypergeometric test, Benjamini-Hochberg corrected P-value. k, Generation of a M03D4.4 loss of function allele. A 1390-bp deletion that removes exons 4-6 and most of exon 7 was induced by CRISPR/Cas9-mediated gene editing.

Source data

Extended Data Fig. 8 Lineage-dependent cell state diversity within each tissue.

a, UMAP plot of neuronal (left, n = 41 for sensory neuron; 54 for inter neuron; 52 for motor neuron) and pharyngeal (right, n = 15 for pharyngeal neuron; 37 for pharyngeal muscle) subtypes based on TF expression. Parameters: min_dist = 0.8, n_neighbors = 5. b, Distribution of tissue types (colors) across terminal cells. c, Comparison of intra-tissue state divergences between all cells and between cells of each tissue progenitor cell lineage (cell pair numbers from left to right: Neu, n = 17,391; 300; 55; 21; 21; 55; 210; 3; 21; 66; 55; 36; 21; 45; 66; 36; 15. Pha, n = 2,415; 36; 55; 6; 91; 6; 6; 3; 10; 10; 10. Ski, n = 4,005; 6; 6; 28; 6; 6; 6; 6; 6; 6; 6; 6; 6; 6. Mus, n = 2,346; 120; 120; 120; 6; 6). The boxplots show the median with IQR and whiskers extending to 1.5 × the IQR; outliers are plotted as points. Statistics: Two-tailed Mann–Whitney U test. d, Changes in intra-tissue state divergence in terminal cells as a function of cell lineage distance in corresponding tissues. State divergence at certain cell lineage distance was not included because of a small number of eligible cells (n < 3). Cell pair numbers from left to right: Pha n = 24; 48; 73; 133; 220; 220; 156; 416; 1,125. Ski n = 32; 61; 82; 124; 297; 13; 254; 773; 126. Mus n = 31; 52; 94; 204; 22; 355; 128; 720. Int n = 4; 8; 18; 36. The boxplots show the median with IQR and whiskers extending to 1.5× the IQR; outliers are plotted as points. e, Changes in intra-sub-tissue state divergence in terminal neuronal and pharyngeal cells as a function of cell lineage distance. Neuronal and pharyngeal cells were divided into subtypes according to known functional subdivisions. Cell pair numbers from left to right: Motor neuron n = 10; 19; 36; 57; 42; 166; 69; 420; 507. Sensory neuron n = 4; 7; 15; 14; 32; 14; 99; 147; 336. Inter neuron n = 7; 18; 31; 75; 116; 119; 268; 589. Pharyngeal muscle n = 5; 13; 22; 27; 62; 45; 80; 70; 342. Pharyngeal neuron n = 5; 11; 15; 16; 50. The boxplots show the median with IQR and whiskers extending to 1.5× the IQR; outliers are plotted as points. f, Tissue and lineage compositions of 20 de novo identified state clusters. Individual tissue progenitor cell lineages were used as the unit to represent the lineage compositions of cells in each state cluster with hyphens linking left-right symmetric lineages. g, Lineage distribution of cells from representative clusters. Connected lines indicate left-right symmetric lineages.

Extended Data Fig. 9 Lineage-restricted expression of tissue-specific TFs.

a, Comparison of tissue coverage of each lineage-specific TFs in cells of all five broad tissue types. b, Expression of tissue-specific TFs (rows) in each clonal cell lineage of corresponding tissues (columns). c, Heatmap showing pair-wise state divergence (range 0-1) between all multipotent progenitor cells (n = 140) calculated using lineage-specific TFs. Cells are ordered first by lineage and then by generation. Cells before the 26-cell stage were not included due to a small number of expression TFs.

Extended Data Fig. 10 State transitions during embryogenesis.

a,b, Comparison of state divergence of cells at each cell generation to the initial and terminal state for cells that differentiate into different tissue types (a) and for cells that are from different cell lineages (b). Comparison numbers from left to right: n = 1.303; 490; 550; 403; 48 for tissue; n = 553; 543; 471; 533; 295; 48; 267; 84 for lineage. c, Heatmap showing the difference between SDm-t and SDd-t calculated by using single-cell transcriptomes following cell divisions (cell track number n = 382, ordered by lineage). d, Comparison of SDm-d at different developmental stages. Mother-daughter cell pair numbers from left to right: n = 20; 50; 102; 176; 334; 343. The boxplot shows the median with IQR and whiskers extending to 1.5× the IQR; outliers are plotted as points. Each dot is the result of a mother-daughter cell comparison. e, Correlation between SDm-d and the difference between SDm-t and SDd-t (each dot represents a comparison, n = 1,938). R and P-value were calculated by a two-sided Pearson’s correlation. Results for the last round of cell division were not included because, by definition, SDm-d equals SDm-t − SDd-t in these cells.

Supplementary information

Supplementary Information

Supplementary Notes 1–6, Supplementary Figs. 1–5 and references.

Reporting Summary

Supplementary Table 1

Selection of TFs. The table lists the names, sources, conservation and related information of all 609 high-confidence TFs encoded in the C. elegans genome (compiled from the wTF 2.0 database and WormBase followed by manual curation) and the selection criteria of 290 TFs as the target of this study.

Supplementary Table 2

Reporter verification and all strains used in this study. a, The sources and verification of fluorescent reporters. b, The strain names and associated genotypes of all dual-fluorescent reporters of TFs. c, The other strains and associated reagents used in this study.

Supplementary Table 3

Cellular expression pattern of all TFs and quality controls. a, Expression levels of 291 protein-fusion TF reporters (for 266 TFs) in 1,204 lineaged cells. TF reporters marked in yellow indicate those that exhibit phenotypic abnormalities of cell proliferation or cell position (sheets 3 and 4). b, The quantitative expression of all TFs in multiple embryos (≥2) for the same reporter strain up to the 350-cell stage. c, Cells in each TF reporter strain that exhibits reproducible defects in cell cycle length. d, Cells in each TF reporter strain that exhibits reproducible cell position defects. e, Previous RNA-seq expression of the TFs that were identified as exhibiting sporadic expression in this study. f, The results of comparing the expression patterns of TFs reported in this study to those reported in the literature. We only considered TFs for which expression patterns in specific lineages, tissue types, cells and developmental stages were previously described. g, The correlation coefficient of single-cell protein expression of TFs between this and previous studies. h, The results of cell-by-cell expression comparisons for 30 TFs with well-documented cellular expression patterns.

Supplementary Table 4

TFs exhibiting various expression specificity. a, All identified lineage-specific TFs, enriched lineages and enrichment scores. b, All identified A–P asymmetric TFs and related bias scores. c, All identified tissue-specific TFs, enriched tissue types and enrichment scores. d, The enriched RNAi phenotypes for the ChIP–seq target genes of each tissue-specific TF. Yellow highlights those phenotypes that are related to the function of corresponding tissue. e, All TFs with transient expression. ‘1’ denotes a TF that is expressed in a mother cell but absent in at least one of its daughters.

Supplementary Table 5

The spatiotemporal cascade of TFs. a, The assignment of TFs into spatiotemporal modules. b, All relationships between TFs inferred based on expression similarity and previous ChIP–seq experiments. c, The inferred TF cascade for each tissue progenitor cell lineage.

Supplementary Table 6

Embryonic cell differentiation phenotypes in TF mutants. a, Changes in cellular expression of CEH-27 in pha-4(zu225) embryos. The table lists cellular expression levels of CEH-27 in all analyzed wild-type and pha-4(zu225) embryos. b, Changes in cellular expression of LIN-32 and SOX-2 in elt-1(ok1002) embryos. The table lists cellular expression levels of LIN-32 and SOX-2 in all analyzed elt-1(ok1002)/+ and elt-1(ok1002) embryos. c, Changes in cellular expression of M03D4.4 in pal-1(RNAi) and pha-4(zu225) embryos. The Table lists changes in cellular expression levels of M03D4.4 in all analyzed embryos upon pal-1 and pha-4 perturbations.

Supplementary Table 7

State clusters of terminal cells. The table lists the lineage and functional name, tissue type and clustering results of all terminally differentiated cells.

Supplementary Table 8

State transitions during cell lineage differentiation. a, State divergences of progenitor cells at all stages to the initial and terminal states. b, State divergences between mother and daughter cells following each cell division.

Source data

Source Data Fig. 2

Statistical source data.

Source Data Fig. 3

Statistical source data.

Source Data Fig. 4

Statistical source data.

Source Data Fig. 5

Statistical source data.

Source Data Fig. 6

Statistical source data.

Source Data Extended Data Fig. 2

Statistical source data.

Source Data Extended Data Fig. 3

Statistical source data.

Source Data Extended Data Fig. 7

Statistical source data.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, X., Zhao, Z., Xiao, L. et al. A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis. Nat Methods 18, 893–902 (2021). https://doi.org/10.1038/s41592-021-01216-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41592-021-01216-1

This article is cited by

Search

Advanced search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing