Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos</title>
<author><name sortKey="Goolam, Mubeen" sort="Goolam, Mubeen" uniqKey="Goolam M" first="Mubeen" last="Goolam">Mubeen Goolam</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Scialdone, Antonio" sort="Scialdone, Antonio" uniqKey="Scialdone A" first="Antonio" last="Scialdone">Antonio Scialdone</name>
<affiliation><nlm:aff id="aff2">European Bioinformatics Institute, European Molecular Biology Laboratory (EMBL-EBI), Wellcome Genome Campus, Cambridge CB10 1SD, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Graham, Sarah J L" sort="Graham, Sarah J L" uniqKey="Graham S" first="Sarah J. L." last="Graham">Sarah J. L. Graham</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Macaulay, Iain C" sort="Macaulay, Iain C" uniqKey="Macaulay I" first="Iain C." last="Macaulay">Iain C. Macaulay</name>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jedrusik, Agnieszka" sort="Jedrusik, Agnieszka" uniqKey="Jedrusik A" first="Agnieszka" last="Jedrusik">Agnieszka Jedrusik</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hupalowska, Anna" sort="Hupalowska, Anna" uniqKey="Hupalowska A" first="Anna" last="Hupalowska">Anna Hupalowska</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Voet, Thierry" sort="Voet, Thierry" uniqKey="Voet T" first="Thierry" last="Voet">Thierry Voet</name>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4">Laboratory of Reproductive Genomics, Department of Human Genetics, KU Leuven, Herestraat 49, 3000 Leuven, Belgium</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Marioni, John C" sort="Marioni, John C" uniqKey="Marioni J" first="John C." last="Marioni">John C. Marioni</name>
<affiliation><nlm:aff id="aff2">European Bioinformatics Institute, European Molecular Biology Laboratory (EMBL-EBI), Wellcome Genome Campus, Cambridge CB10 1SD, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5">Cancer Research UK-Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Cambridge CB2 0RE, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zernicka Goetz, Magdalena" sort="Zernicka Goetz, Magdalena" uniqKey="Zernicka Goetz M" first="Magdalena" last="Zernicka-Goetz">Magdalena Zernicka-Goetz</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">27015307</idno>
<idno type="pmc">4819611</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4819611</idno>
<idno type="RBID">PMC:4819611</idno>
<idno type="doi">10.1016/j.cell.2016.01.047</idno>
<date when="2016">2016</date>
<idno type="wicri:Area/Pmc/Corpus">000352</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos</title>
<author><name sortKey="Goolam, Mubeen" sort="Goolam, Mubeen" uniqKey="Goolam M" first="Mubeen" last="Goolam">Mubeen Goolam</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Scialdone, Antonio" sort="Scialdone, Antonio" uniqKey="Scialdone A" first="Antonio" last="Scialdone">Antonio Scialdone</name>
<affiliation><nlm:aff id="aff2">European Bioinformatics Institute, European Molecular Biology Laboratory (EMBL-EBI), Wellcome Genome Campus, Cambridge CB10 1SD, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Graham, Sarah J L" sort="Graham, Sarah J L" uniqKey="Graham S" first="Sarah J. L." last="Graham">Sarah J. L. Graham</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Macaulay, Iain C" sort="Macaulay, Iain C" uniqKey="Macaulay I" first="Iain C." last="Macaulay">Iain C. Macaulay</name>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jedrusik, Agnieszka" sort="Jedrusik, Agnieszka" uniqKey="Jedrusik A" first="Agnieszka" last="Jedrusik">Agnieszka Jedrusik</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hupalowska, Anna" sort="Hupalowska, Anna" uniqKey="Hupalowska A" first="Anna" last="Hupalowska">Anna Hupalowska</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Voet, Thierry" sort="Voet, Thierry" uniqKey="Voet T" first="Thierry" last="Voet">Thierry Voet</name>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4">Laboratory of Reproductive Genomics, Department of Human Genetics, KU Leuven, Herestraat 49, 3000 Leuven, Belgium</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Marioni, John C" sort="Marioni, John C" uniqKey="Marioni J" first="John C." last="Marioni">John C. Marioni</name>
<affiliation><nlm:aff id="aff2">European Bioinformatics Institute, European Molecular Biology Laboratory (EMBL-EBI), Wellcome Genome Campus, Cambridge CB10 1SD, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3">Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5">Cancer Research UK-Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Cambridge CB2 0RE, UK</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zernicka Goetz, Magdalena" sort="Zernicka Goetz, Magdalena" uniqKey="Zernicka Goetz M" first="Magdalena" last="Zernicka-Goetz">Magdalena Zernicka-Goetz</name>
<affiliation><nlm:aff id="aff1">Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Cell</title>
<idno type="ISSN">0092-8674</idno>
<idno type="eISSN">1097-4172</idno>
<imprint><date when="2016">2016</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><title>Summary</title>
<p>The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Anders, S" uniqKey="Anders S">S. Anders</name>
</author>
<author><name sortKey="Pyl, P T" uniqKey="Pyl P">P.T. Pyl</name>
</author>
<author><name sortKey="Huber, W" uniqKey="Huber W">W. Huber</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Biase, F H" uniqKey="Biase F">F.H. Biase</name>
</author>
<author><name sortKey="Cao, X" uniqKey="Cao X">X. Cao</name>
</author>
<author><name sortKey="Zhong, S" uniqKey="Zhong S">S. Zhong</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bischoff, M" uniqKey="Bischoff M">M. Bischoff</name>
</author>
<author><name sortKey="Parfitt, D E" uniqKey="Parfitt D">D.E. Parfitt</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brennecke, P" uniqKey="Brennecke P">P. Brennecke</name>
</author>
<author><name sortKey="Anders, S" uniqKey="Anders S">S. Anders</name>
</author>
<author><name sortKey="Kim, J K" uniqKey="Kim J">J.K. Kim</name>
</author>
<author><name sortKey="Kolodziejczyk, A A" uniqKey="Kolodziejczyk A">A.A. Kołodziejczyk</name>
</author>
<author><name sortKey="Zhang, X" uniqKey="Zhang X">X. Zhang</name>
</author>
<author><name sortKey="Proserpio, V" uniqKey="Proserpio V">V. Proserpio</name>
</author>
<author><name sortKey="Baying, B" uniqKey="Baying B">B. Baying</name>
</author>
<author><name sortKey="Benes, V" uniqKey="Benes V">V. Benes</name>
</author>
<author><name sortKey="Teichmann, S A" uniqKey="Teichmann S">S.A. Teichmann</name>
</author>
<author><name sortKey="Marioni, J C" uniqKey="Marioni J">J.C. Marioni</name>
</author>
<author><name sortKey="Heisler, M G" uniqKey="Heisler M">M.G. Heisler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buettner, F" uniqKey="Buettner F">F. Buettner</name>
</author>
<author><name sortKey="Natarajan, K N" uniqKey="Natarajan K">K.N. Natarajan</name>
</author>
<author><name sortKey="Casale, F P" uniqKey="Casale F">F.P. Casale</name>
</author>
<author><name sortKey="Proserpio, V" uniqKey="Proserpio V">V. Proserpio</name>
</author>
<author><name sortKey="Scialdone, A" uniqKey="Scialdone A">A. Scialdone</name>
</author>
<author><name sortKey="Theis, F J" uniqKey="Theis F">F.J. Theis</name>
</author>
<author><name sortKey="Teichmann, S A" uniqKey="Teichmann S">S.A. Teichmann</name>
</author>
<author><name sortKey="Marioni, J C" uniqKey="Marioni J">J.C. Marioni</name>
</author>
<author><name sortKey="Stegle, O" uniqKey="Stegle O">O. Stegle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burton, A" uniqKey="Burton A">A. Burton</name>
</author>
<author><name sortKey="Muller, J" uniqKey="Muller J">J. Muller</name>
</author>
<author><name sortKey="Tu, S" uniqKey="Tu S">S. Tu</name>
</author>
<author><name sortKey="Padilla Longoria, P" uniqKey="Padilla Longoria P">P. Padilla-Longoria</name>
</author>
<author><name sortKey="Guccione, E" uniqKey="Guccione E">E. Guccione</name>
</author>
<author><name sortKey="Torres Padilla, M E" uniqKey="Torres Padilla M">M.E. Torres-Padilla</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chakravarthy, H" uniqKey="Chakravarthy H">H. Chakravarthy</name>
</author>
<author><name sortKey="Ormsbee, B D" uniqKey="Ormsbee B">B.D. Ormsbee</name>
</author>
<author><name sortKey="Mallanna, S K" uniqKey="Mallanna S">S.K. Mallanna</name>
</author>
<author><name sortKey="Rizzino, A" uniqKey="Rizzino A">A. Rizzino</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deng, Q" uniqKey="Deng Q">Q. Deng</name>
</author>
<author><name sortKey="Ramskold, D" uniqKey="Ramskold D">D. Ramsköld</name>
</author>
<author><name sortKey="Reinius, B" uniqKey="Reinius B">B. Reinius</name>
</author>
<author><name sortKey="Sandberg, R" uniqKey="Sandberg R">R. Sandberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fleming, T P" uniqKey="Fleming T">T.P. Fleming</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gardner, R L" uniqKey="Gardner R">R.L. Gardner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goke, J" uniqKey="Goke J">J. Göke</name>
</author>
<author><name sortKey="Jung, M" uniqKey="Jung M">M. Jung</name>
</author>
<author><name sortKey="Behrens, S" uniqKey="Behrens S">S. Behrens</name>
</author>
<author><name sortKey="Chavez, L" uniqKey="Chavez L">L. Chavez</name>
</author>
<author><name sortKey="O Eeffe, S" uniqKey="O Eeffe S">S. O’Keeffe</name>
</author>
<author><name sortKey="Timmermann, B" uniqKey="Timmermann B">B. Timmermann</name>
</author>
<author><name sortKey="Lehrach, H" uniqKey="Lehrach H">H. Lehrach</name>
</author>
<author><name sortKey="Adjaye, J" uniqKey="Adjaye J">J. Adjaye</name>
</author>
<author><name sortKey="Vingron, M" uniqKey="Vingron M">M. Vingron</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Guo, G" uniqKey="Guo G">G. Guo</name>
</author>
<author><name sortKey="Huss, M" uniqKey="Huss M">M. Huss</name>
</author>
<author><name sortKey="Tong, G Q" uniqKey="Tong G">G.Q. Tong</name>
</author>
<author><name sortKey="Wang, C" uniqKey="Wang C">C. Wang</name>
</author>
<author><name sortKey="Li Sun, L" uniqKey="Li Sun L">L. Li Sun</name>
</author>
<author><name sortKey="Clarke, N D" uniqKey="Clarke N">N.D. Clarke</name>
</author>
<author><name sortKey="Robson, P" uniqKey="Robson P">P. Robson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hillman, N" uniqKey="Hillman N">N. Hillman</name>
</author>
<author><name sortKey="Sherman, M I" uniqKey="Sherman M">M.I. Sherman</name>
</author>
<author><name sortKey="Graham, C" uniqKey="Graham C">C. Graham</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jedrusik, A" uniqKey="Jedrusik A">A. Jedrusik</name>
</author>
<author><name sortKey="Parfitt, D E" uniqKey="Parfitt D">D.E. Parfitt</name>
</author>
<author><name sortKey="Guo, G" uniqKey="Guo G">G. Guo</name>
</author>
<author><name sortKey="Skamagki, M" uniqKey="Skamagki M">M. Skamagki</name>
</author>
<author><name sortKey="Grabarek, J B" uniqKey="Grabarek J">J.B. Grabarek</name>
</author>
<author><name sortKey="Johnson, M H" uniqKey="Johnson M">M.H. Johnson</name>
</author>
<author><name sortKey="Robson, P" uniqKey="Robson P">P. Robson</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Johnson, M H" uniqKey="Johnson M">M.H. Johnson</name>
</author>
<author><name sortKey="Ziomek, C A" uniqKey="Ziomek C">C.A. Ziomek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Johnson, W E" uniqKey="Johnson W">W.E. Johnson</name>
</author>
<author><name sortKey="Li, C" uniqKey="Li C">C. Li</name>
</author>
<author><name sortKey="Rabinovic, A" uniqKey="Rabinovic A">A. Rabinovic</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kiso, M" uniqKey="Kiso M">M. Kiso</name>
</author>
<author><name sortKey="Tanaka, S" uniqKey="Tanaka S">S. Tanaka</name>
</author>
<author><name sortKey="Saba, R" uniqKey="Saba R">R. Saba</name>
</author>
<author><name sortKey="Matsuda, S" uniqKey="Matsuda S">S. Matsuda</name>
</author>
<author><name sortKey="Shimizu, A" uniqKey="Shimizu A">A. Shimizu</name>
</author>
<author><name sortKey="Ohyama, M" uniqKey="Ohyama M">M. Ohyama</name>
</author>
<author><name sortKey="Okano, H J" uniqKey="Okano H">H.J. Okano</name>
</author>
<author><name sortKey="Shiroishi, T" uniqKey="Shiroishi T">T. Shiroishi</name>
</author>
<author><name sortKey="Okano, H" uniqKey="Okano H">H. Okano</name>
</author>
<author><name sortKey="Saga, Y" uniqKey="Saga Y">Y. Saga</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kuzmichev, A N" uniqKey="Kuzmichev A">A.N. Kuzmichev</name>
</author>
<author><name sortKey="Kim, S K" uniqKey="Kim S">S.K. Kim</name>
</author>
<author><name sortKey="D Lessio, A C" uniqKey="D Lessio A">A.C. D’Alessio</name>
</author>
<author><name sortKey="Chenoweth, J G" uniqKey="Chenoweth J">J.G. Chenoweth</name>
</author>
<author><name sortKey="Wittko, I M" uniqKey="Wittko I">I.M. Wittko</name>
</author>
<author><name sortKey="Campanati, L" uniqKey="Campanati L">L. Campanati</name>
</author>
<author><name sortKey="Mckay, R D" uniqKey="Mckay R">R.D. McKay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Langfelder, P" uniqKey="Langfelder P">P. Langfelder</name>
</author>
<author><name sortKey="Zhang, B" uniqKey="Zhang B">B. Zhang</name>
</author>
<author><name sortKey="Horvath, S" uniqKey="Horvath S">S. Horvath</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Loh, Y H" uniqKey="Loh Y">Y.H. Loh</name>
</author>
<author><name sortKey="Wu, Q" uniqKey="Wu Q">Q. Wu</name>
</author>
<author><name sortKey="Chew, J L" uniqKey="Chew J">J.L. Chew</name>
</author>
<author><name sortKey="Vega, V B" uniqKey="Vega V">V.B. Vega</name>
</author>
<author><name sortKey="Zhang, W" uniqKey="Zhang W">W. Zhang</name>
</author>
<author><name sortKey="Chen, X" uniqKey="Chen X">X. Chen</name>
</author>
<author><name sortKey="Bourque, G" uniqKey="Bourque G">G. Bourque</name>
</author>
<author><name sortKey="George, J" uniqKey="George J">J. George</name>
</author>
<author><name sortKey="Leong, B" uniqKey="Leong B">B. Leong</name>
</author>
<author><name sortKey="Liu, J" uniqKey="Liu J">J. Liu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Love, M I" uniqKey="Love M">M.I. Love</name>
</author>
<author><name sortKey="Huber, W" uniqKey="Huber W">W. Huber</name>
</author>
<author><name sortKey="Anders, S" uniqKey="Anders S">S. Anders</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mallanna, S K" uniqKey="Mallanna S">S.K. Mallanna</name>
</author>
<author><name sortKey="Ormsbee, B D" uniqKey="Ormsbee B">B.D. Ormsbee</name>
</author>
<author><name sortKey="Iacovino, M" uniqKey="Iacovino M">M. Iacovino</name>
</author>
<author><name sortKey="Gilmore, J M" uniqKey="Gilmore J">J.M. Gilmore</name>
</author>
<author><name sortKey="Cox, J L" uniqKey="Cox J">J.L. Cox</name>
</author>
<author><name sortKey="Kyba, M" uniqKey="Kyba M">M. Kyba</name>
</author>
<author><name sortKey="Washburn, M P" uniqKey="Washburn M">M.P. Washburn</name>
</author>
<author><name sortKey="Rizzino, A" uniqKey="Rizzino A">A. Rizzino</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Matoba, R" uniqKey="Matoba R">R. Matoba</name>
</author>
<author><name sortKey="Niwa, H" uniqKey="Niwa H">H. Niwa</name>
</author>
<author><name sortKey="Masui, S" uniqKey="Masui S">S. Masui</name>
</author>
<author><name sortKey="Ohtsuka, S" uniqKey="Ohtsuka S">S. Ohtsuka</name>
</author>
<author><name sortKey="Carter, M G" uniqKey="Carter M">M.G. Carter</name>
</author>
<author><name sortKey="Sharov, A A" uniqKey="Sharov A">A.A. Sharov</name>
</author>
<author><name sortKey="Ko, M S" uniqKey="Ko M">M.S. Ko</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mcdole, K" uniqKey="Mcdole K">K. McDole</name>
</author>
<author><name sortKey="Xiong, Y" uniqKey="Xiong Y">Y. Xiong</name>
</author>
<author><name sortKey="Iglesias, P A" uniqKey="Iglesias P">P.A. Iglesias</name>
</author>
<author><name sortKey="Zheng, Y" uniqKey="Zheng Y">Y. Zheng</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Morris, S A" uniqKey="Morris S">S.A. Morris</name>
</author>
<author><name sortKey="Teo, R T" uniqKey="Teo R">R.T. Teo</name>
</author>
<author><name sortKey="Li, H" uniqKey="Li H">H. Li</name>
</author>
<author><name sortKey="Robson, P" uniqKey="Robson P">P. Robson</name>
</author>
<author><name sortKey="Glover, D M" uniqKey="Glover D">D.M. Glover</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Morris, S A" uniqKey="Morris S">S.A. Morris</name>
</author>
<author><name sortKey="Guo, Y" uniqKey="Guo Y">Y. Guo</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parfitt, D E" uniqKey="Parfitt D">D.E. Parfitt</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Picelli, S" uniqKey="Picelli S">S. Picelli</name>
</author>
<author><name sortKey="Faridani, O R" uniqKey="Faridani O">O.R. Faridani</name>
</author>
<author><name sortKey="Bjorklund, A K" uniqKey="Bjorklund A">A.K. Björklund</name>
</author>
<author><name sortKey="Winberg, G" uniqKey="Winberg G">G. Winberg</name>
</author>
<author><name sortKey="Sagasser, S" uniqKey="Sagasser S">S. Sagasser</name>
</author>
<author><name sortKey="Sandberg, R" uniqKey="Sandberg R">R. Sandberg</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Piotrowska, K" uniqKey="Piotrowska K">K. Piotrowska</name>
</author>
<author><name sortKey="Wianny, F" uniqKey="Wianny F">F. Wianny</name>
</author>
<author><name sortKey="Pedersen, R A" uniqKey="Pedersen R">R.A. Pedersen</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Piotrowska Nitsche, K" uniqKey="Piotrowska Nitsche K">K. Piotrowska-Nitsche</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Piotrowska Nitsche, K" uniqKey="Piotrowska Nitsche K">K. Piotrowska-Nitsche</name>
</author>
<author><name sortKey="Perea Gomez, A" uniqKey="Perea Gomez A">A. Perea-Gomez</name>
</author>
<author><name sortKey="Haraguchi, S" uniqKey="Haraguchi S">S. Haraguchi</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Plachta, N" uniqKey="Plachta N">N. Plachta</name>
</author>
<author><name sortKey="Bollenbach, T" uniqKey="Bollenbach T">T. Bollenbach</name>
</author>
<author><name sortKey="Pease, S" uniqKey="Pease S">S. Pease</name>
</author>
<author><name sortKey="Fraser, S E" uniqKey="Fraser S">S.E. Fraser</name>
</author>
<author><name sortKey="Pantazis, P" uniqKey="Pantazis P">P. Pantazis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Samarage, C R" uniqKey="Samarage C">C.R. Samarage</name>
</author>
<author><name sortKey="White, M D" uniqKey="White M">M.D. White</name>
</author>
<author><name sortKey="Alvarez, Y D" uniqKey="Alvarez Y">Y.D. Álvarez</name>
</author>
<author><name sortKey="Fierro Gonzalez, J C" uniqKey="Fierro Gonzalez J">J.C. Fierro-González</name>
</author>
<author><name sortKey="Henon, Y" uniqKey="Henon Y">Y. Henon</name>
</author>
<author><name sortKey="Jesudason, E C" uniqKey="Jesudason E">E.C. Jesudason</name>
</author>
<author><name sortKey="Bissiere, S" uniqKey="Bissiere S">S. Bissiere</name>
</author>
<author><name sortKey="Fouras, A" uniqKey="Fouras A">A. Fouras</name>
</author>
<author><name sortKey="Plachta, N" uniqKey="Plachta N">N. Plachta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sharov, A A" uniqKey="Sharov A">A.A. Sharov</name>
</author>
<author><name sortKey="Masui, S" uniqKey="Masui S">S. Masui</name>
</author>
<author><name sortKey="Sharova, L V" uniqKey="Sharova L">L.V. Sharova</name>
</author>
<author><name sortKey="Piao, Y" uniqKey="Piao Y">Y. Piao</name>
</author>
<author><name sortKey="Aiba, K" uniqKey="Aiba K">K. Aiba</name>
</author>
<author><name sortKey="Matoba, R" uniqKey="Matoba R">R. Matoba</name>
</author>
<author><name sortKey="Xin, L" uniqKey="Xin L">L. Xin</name>
</author>
<author><name sortKey="Niwa, H" uniqKey="Niwa H">H. Niwa</name>
</author>
<author><name sortKey="Ko, M S" uniqKey="Ko M">M.S. Ko</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shi, J" uniqKey="Shi J">J. Shi</name>
</author>
<author><name sortKey="Chen, Q" uniqKey="Chen Q">Q. Chen</name>
</author>
<author><name sortKey="Li, X" uniqKey="Li X">X. Li</name>
</author>
<author><name sortKey="Zheng, X" uniqKey="Zheng X">X. Zheng</name>
</author>
<author><name sortKey="Zhang, Y" uniqKey="Zhang Y">Y. Zhang</name>
</author>
<author><name sortKey="Qiao, J" uniqKey="Qiao J">J. Qiao</name>
</author>
<author><name sortKey="Tang, F" uniqKey="Tang F">F. Tang</name>
</author>
<author><name sortKey="Tao, Y" uniqKey="Tao Y">Y. Tao</name>
</author>
<author><name sortKey="Zhou, Q" uniqKey="Zhou Q">Q. Zhou</name>
</author>
<author><name sortKey="Duan, E" uniqKey="Duan E">E. Duan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stegle, O" uniqKey="Stegle O">O. Stegle</name>
</author>
<author><name sortKey="Teichmann, S A" uniqKey="Teichmann S">S.A. Teichmann</name>
</author>
<author><name sortKey="Marioni, J C" uniqKey="Marioni J">J.C. Marioni</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Strumpf, D" uniqKey="Strumpf D">D. Strumpf</name>
</author>
<author><name sortKey="Mao, C A" uniqKey="Mao C">C.A. Mao</name>
</author>
<author><name sortKey="Yamanaka, Y" uniqKey="Yamanaka Y">Y. Yamanaka</name>
</author>
<author><name sortKey="Ralston, A" uniqKey="Ralston A">A. Ralston</name>
</author>
<author><name sortKey="Chawengsaksophak, K" uniqKey="Chawengsaksophak K">K. Chawengsaksophak</name>
</author>
<author><name sortKey="Beck, F" uniqKey="Beck F">F. Beck</name>
</author>
<author><name sortKey="Rossant, J" uniqKey="Rossant J">J. Rossant</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tabansky, I" uniqKey="Tabansky I">I. Tabansky</name>
</author>
<author><name sortKey="Lenarcic, A" uniqKey="Lenarcic A">A. Lenarcic</name>
</author>
<author><name sortKey="Draft, R W" uniqKey="Draft R">R.W. Draft</name>
</author>
<author><name sortKey="Loulier, K" uniqKey="Loulier K">K. Loulier</name>
</author>
<author><name sortKey="Keskin, D B" uniqKey="Keskin D">D.B. Keskin</name>
</author>
<author><name sortKey="Rosains, J" uniqKey="Rosains J">J. Rosains</name>
</author>
<author><name sortKey="Rivera Feliciano, J" uniqKey="Rivera Feliciano J">J. Rivera-Feliciano</name>
</author>
<author><name sortKey="Lichtman, J W" uniqKey="Lichtman J">J.W. Lichtman</name>
</author>
<author><name sortKey="Livet, J" uniqKey="Livet J">J. Livet</name>
</author>
<author><name sortKey="Stern, J N" uniqKey="Stern J">J.N. Stern</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Takaoka, K" uniqKey="Takaoka K">K. Takaoka</name>
</author>
<author><name sortKey="Hamada, H" uniqKey="Hamada H">H. Hamada</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tang, F" uniqKey="Tang F">F. Tang</name>
</author>
<author><name sortKey="Barbacioru, C" uniqKey="Barbacioru C">C. Barbacioru</name>
</author>
<author><name sortKey="Nordman, E" uniqKey="Nordman E">E. Nordman</name>
</author>
<author><name sortKey="Bao, S" uniqKey="Bao S">S. Bao</name>
</author>
<author><name sortKey="Lee, C" uniqKey="Lee C">C. Lee</name>
</author>
<author><name sortKey="Wang, X" uniqKey="Wang X">X. Wang</name>
</author>
<author><name sortKey="Tuch, B B" uniqKey="Tuch B">B.B. Tuch</name>
</author>
<author><name sortKey="Heard, E" uniqKey="Heard E">E. Heard</name>
</author>
<author><name sortKey="Lao, K" uniqKey="Lao K">K. Lao</name>
</author>
<author><name sortKey="Surani, M A" uniqKey="Surani M">M.A. Surani</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tarkowski, A K" uniqKey="Tarkowski A">A.K. Tarkowski</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Torres Padilla, M E" uniqKey="Torres Padilla M">M.E. Torres-Padilla</name>
</author>
<author><name sortKey="Parfitt, D E" uniqKey="Parfitt D">D.E. Parfitt</name>
</author>
<author><name sortKey="Kouzarides, T" uniqKey="Kouzarides T">T. Kouzarides</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watanabe, T" uniqKey="Watanabe T">T. Watanabe</name>
</author>
<author><name sortKey="Biggins, J S" uniqKey="Biggins J">J.S. Biggins</name>
</author>
<author><name sortKey="Tannan, N B" uniqKey="Tannan N">N.B. Tannan</name>
</author>
<author><name sortKey="Srinivas, S" uniqKey="Srinivas S">S. Srinivas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wu, Q" uniqKey="Wu Q">Q. Wu</name>
</author>
<author><name sortKey="Bruce, A W" uniqKey="Bruce A">A.W. Bruce</name>
</author>
<author><name sortKey="Jedrusik, A" uniqKey="Jedrusik A">A. Jedrusik</name>
</author>
<author><name sortKey="Ellis, P D" uniqKey="Ellis P">P.D. Ellis</name>
</author>
<author><name sortKey="Andrews, R M" uniqKey="Andrews R">R.M. Andrews</name>
</author>
<author><name sortKey="Langford, C F" uniqKey="Langford C">C.F. Langford</name>
</author>
<author><name sortKey="Glover, D M" uniqKey="Glover D">D.M. Glover</name>
</author>
<author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Xue, Z" uniqKey="Xue Z">Z. Xue</name>
</author>
<author><name sortKey="Huang, K" uniqKey="Huang K">K. Huang</name>
</author>
<author><name sortKey="Cai, C" uniqKey="Cai C">C. Cai</name>
</author>
<author><name sortKey="Cai, L" uniqKey="Cai L">L. Cai</name>
</author>
<author><name sortKey="Jiang, C Y" uniqKey="Jiang C">C.Y. Jiang</name>
</author>
<author><name sortKey="Feng, Y" uniqKey="Feng Y">Y. Feng</name>
</author>
<author><name sortKey="Liu, Z" uniqKey="Liu Z">Z. Liu</name>
</author>
<author><name sortKey="Zeng, Q" uniqKey="Zeng Q">Q. Zeng</name>
</author>
<author><name sortKey="Cheng, L" uniqKey="Cheng L">L. Cheng</name>
</author>
<author><name sortKey="Sun, Y E" uniqKey="Sun Y">Y.E. Sun</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yamanaka, Y" uniqKey="Yamanaka Y">Y. Yamanaka</name>
</author>
<author><name sortKey="Lanner, F" uniqKey="Lanner F">F. Lanner</name>
</author>
<author><name sortKey="Rossant, J" uniqKey="Rossant J">J. Rossant</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zernicka Goetz, M" uniqKey="Zernicka Goetz M">M. Zernicka-Goetz</name>
</author>
<author><name sortKey="Morris, S A" uniqKey="Morris S">S.A. Morris</name>
</author>
<author><name sortKey="Bruce, A W" uniqKey="Bruce A">A.W. Bruce</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Cell</journal-id>
<journal-id journal-id-type="iso-abbrev">Cell</journal-id>
<journal-title-group><journal-title>Cell</journal-title>
</journal-title-group>
<issn pub-type="ppub">0092-8674</issn>
<issn pub-type="epub">1097-4172</issn>
<publisher><publisher-name>Cell Press</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">27015307</article-id>
<article-id pub-id-type="pmc">4819611</article-id>
<article-id pub-id-type="publisher-id">S0092-8674(16)30061-7</article-id>
<article-id pub-id-type="doi">10.1016/j.cell.2016.01.047</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Goolam</surname>
<given-names>Mubeen</given-names>
</name>
<xref rid="aff1" ref-type="aff">1</xref>
<xref rid="fn1" ref-type="fn">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Scialdone</surname>
<given-names>Antonio</given-names>
</name>
<xref rid="aff2" ref-type="aff">2</xref>
<xref rid="aff3" ref-type="aff">3</xref>
<xref rid="fn1" ref-type="fn">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Graham</surname>
<given-names>Sarah J.L.</given-names>
</name>
<xref rid="aff1" ref-type="aff">1</xref>
<xref rid="fn1" ref-type="fn">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Macaulay</surname>
<given-names>Iain C.</given-names>
</name>
<xref rid="aff3" ref-type="aff">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Jedrusik</surname>
<given-names>Agnieszka</given-names>
</name>
<xref rid="aff1" ref-type="aff">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hupalowska</surname>
<given-names>Anna</given-names>
</name>
<xref rid="aff1" ref-type="aff">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Voet</surname>
<given-names>Thierry</given-names>
</name>
<xref rid="aff3" ref-type="aff">3</xref>
<xref rid="aff4" ref-type="aff">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Marioni</surname>
<given-names>John C.</given-names>
</name>
<email>marioni@ebi.ac.uk</email>
<xref rid="aff2" ref-type="aff">2</xref>
<xref rid="aff3" ref-type="aff">3</xref>
<xref rid="aff5" ref-type="aff">5</xref>
<xref rid="cor1" ref-type="corresp">∗</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Zernicka-Goetz</surname>
<given-names>Magdalena</given-names>
</name>
<email>mz205@cam.ac.uk</email>
<xref rid="aff1" ref-type="aff">1</xref>
<xref rid="cor2" ref-type="corresp">∗∗</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
Department of Physiology, Development & Neuroscience, University of Cambridge, Downing Street, Cambridge CB2 3EG, UK</aff>
<aff id="aff2"><label>2</label>
European Bioinformatics Institute, European Molecular Biology Laboratory (EMBL-EBI), Wellcome Genome Campus, Cambridge CB10 1SD, UK</aff>
<aff id="aff3"><label>3</label>
Wellcome Trust Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK</aff>
<aff id="aff4"><label>4</label>
Laboratory of Reproductive Genomics, Department of Human Genetics, KU Leuven, Herestraat 49, 3000 Leuven, Belgium</aff>
<aff id="aff5"><label>5</label>
Cancer Research UK-Cambridge Institute, University of Cambridge, Li Ka Shing Centre, Cambridge CB2 0RE, UK</aff>
<author-notes><corresp id="cor1"><label>∗</label>
Corresponding author <email>marioni@ebi.ac.uk</email>
</corresp>
<corresp id="cor2"><label>∗∗</label>
Corresponding author <email>mz205@cam.ac.uk</email>
</corresp>
<fn id="fn1"><label>6</label>
<p id="ntpara0010">Co-first author</p>
</fn>
</author-notes>
<pub-date pub-type="pmc-release"><day>24</day>
<month>3</month>
<year>2016</year>
</pub-date>
<pmc-comment> PMC Release delay is 0 months and 0 days and was based on .</pmc-comment>
<pub-date pub-type="ppub"><day>24</day>
<month>3</month>
<year>2016</year>
</pub-date>
<volume>165</volume>
<issue>1</issue>
<fpage>61</fpage>
<lpage>74</lpage>
<history><date date-type="received"><day>9</day>
<month>6</month>
<year>2015</year>
</date>
<date date-type="rev-recd"><day>10</day>
<month>11</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>22</day>
<month>1</month>
<year>2016</year>
</date>
</history>
<permissions><copyright-statement>© 2016 The Authors</copyright-statement>
<copyright-year>2016</copyright-year>
<license license-type="CC BY" xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).</license-p>
</license>
</permissions>
<related-article related-article-type="article-reference" id="d32e12" ext-link-type="doi" xlink:href="10.1016/j.cell.2016.02.032"></related-article>
<abstract><title>Summary</title>
<p>The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation.</p>
</abstract>
<abstract abstract-type="graphical"><title>Graphical Abstract</title>
<fig id="undfig1" position="anchor"><graphic xlink:href="fx1"></graphic>
</fig>
</abstract>
<abstract abstract-type="author-highlights"><title>Highlights</title>
<p><list list-type="simple"><list-item id="u0010"><label>•</label>
<p>Oct4 and Sox2 targets such as Sox21 are highly heterogeneous in the 4-cell embryo</p>
</list-item>
<list-item id="u0015"><label>•</label>
<p>Depleting Sox21 levels leads to upregulation of Cdx2 and extra-embryonic fate</p>
</list-item>
<list-item id="u0020"><label>•</label>
<p>Sox21 expression responds to the histone H3R26 methyltransferase CARM1</p>
</list-item>
<list-item id="u0025"><label>•</label>
<p>Heterogeneous gene expression initiates cell-fate decisions at the 4-cell stage</p>
</list-item>
</list>
</p>
</abstract>
<abstract abstract-type="teaser"><p>Heterogeneous gene expression in 4-cell stage in mouse embryos contributes to cell-fate decisions, with lower levels of Sox21 promoting extra-embryonic fate over pluripotency.</p>
</abstract>
</article-meta>
<notes><p id="misc0010">Published: March 24, 2016</p>
</notes>
</front>
<body><sec id="sec1"><title>Introduction</title>
<p>When in mammalian development cells first start to differ from each other and whether these first differences play any role in cell-fate specification remain key open questions. In many model systems, initiation of cell-fate specification stems from heterogeneity between the blastomeres of the early embryo, but whether this might also be the case in mammals remains unknown. The first cell-fate specification in the mammalian embryo leads to the separation of embryonic and extra-embryonic lineages. The embryonic lineage is pluripotent and will give rise to the fetus, while the extra-embryonic lineages will differentiate into supportive structures critical for embryo implantation and fetal development, the placenta, and yolk sac (<xref rid="bib39" ref-type="bibr">Takaoka and Hamada, 2012</xref>
, <xref rid="bib47" ref-type="bibr">Zernicka-Goetz et al., 2009</xref>
). How and when these lineages start to separate in morphologically homogenous cells has been very difficult to dissect in mammals. Historically, cells of the early mouse embryo were considered identical in their ability to give rise to embryonic or extra-embryonic lineages, due to the regulative ability of the embryo to compensate for alterations in cell arrangement (<xref rid="bib13" ref-type="bibr">Hillman et al., 1972</xref>
, <xref rid="bib41" ref-type="bibr">Tarkowski, 1959</xref>
). However, more recent evidence has suggested that cells as early as the 4-cell stage become heterogeneous, exhibiting differences in developmental fate and potential (<xref rid="bib3" ref-type="bibr">Bischoff et al., 2008</xref>
, <xref rid="bib31" ref-type="bibr">Piotrowska-Nitsche et al., 2005</xref>
, <xref rid="bib38" ref-type="bibr">Tabansky et al., 2013</xref>
) and in the activity of specific cell-fate regulators (<xref rid="bib6" ref-type="bibr">Burton et al., 2013</xref>
, <xref rid="bib32" ref-type="bibr">Plachta et al., 2011</xref>
, <xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
). This heterogeneity indicates the possibility that the breaking of embryo symmetry starts earlier than expected, prior to differences in cell position and polarity evident from the 16-cell-stage onward (<xref rid="bib9" ref-type="bibr">Fleming, 1987</xref>
, <xref rid="bib15" ref-type="bibr">Johnson and Ziomek, 1981</xref>
). However, finding causal links between this early heterogeneity and later lineage divergence has proved extremely difficult because the key evidence—differences in gene expression patterns between individual cells that regulate cell fate—has, until now, been hard to identify due to technical limitations.</p>
<p>High-throughput single-cell transcriptomics offers an unbiased approach for understanding the extent, basis, and function of gene expression variation between seemingly identical cells. So far, the focus of single-cell studies in the mouse embryo has been on gene expression patterns that characterize particular developmental stages or lineages within the blastocyst or mono versus bi-allelic gene expression (<xref rid="bib2" ref-type="bibr">Biase et al., 2014</xref>
, <xref rid="bib8" ref-type="bibr">Deng et al., 2014</xref>
, <xref rid="bib12" ref-type="bibr">Guo et al., 2010</xref>
, <xref rid="bib35" ref-type="bibr">Shi et al., 2015</xref>
, <xref rid="bib40" ref-type="bibr">Tang et al., 2011</xref>
, <xref rid="bib45" ref-type="bibr">Xue et al., 2013</xref>
), rather than on investigating the functional consequences of heterogeneity within the same embryo for cell-fate specification. Here, using single-cell transcriptomics, we determined the extent of transcriptional heterogeneities between individual cells in pre-implantation embryos and identified that target genes of the pluripotency master regulators Oct4 and Sox2 are highly heterogeneous at the 4-cell stage. We find that <italic>Sox21</italic>
, a Sox2 target that also displays proximal binding of Oct4 and that regulates expression of the master regulator of differentiation Cdx2, is one of the most highly heterogeneous genes at this stage. We also find that decreasing expression of Sox21 in individual cells results in premature and increased expression of Cdx2 and decreased contribution to the pluripotent lineage. We identify CARM1, which methylates histone H3 at arginine 26 (H3R26), as an upstream regulator of Sox21 expression. These results indicate that heterogeneity in gene expression patterns biases cell fate in the mouse embryo as early as the 4-cell stage.</p>
</sec>
<sec id="sec2"><title>Results</title>
<sec id="sec2.1"><title>Temporal and Spatial Gene Expression Patterns in Single Cells during Pre-implantation Development</title>
<p>To assess gene expression patterns throughout the pre-implantation development of mouse embryos, we isolated individual cells from successive stages and assayed their transcriptomes using the Smart-Seq2 single-cell RNA-sequencing protocol (<xref rid="bib28" ref-type="bibr">Picelli et al., 2014</xref>
). Transcriptomes were determined for all blastomeres of 28 embryos isolated at the 2-cell (n = 8), 4-cell (n = 16), and 8-cell (n = 4) stages, and for individual cells at the 16-cell (n = 6) and 32-cell (n = 6) stages, corresponding to the morula and blastocyst, respectively (<xref rid="fig1" ref-type="fig">Figure 1</xref>
A; <xref rid="mmc1" ref-type="supplementary-material">Table S1</xref>
). To assess the quality of the data we compared, across cells in each batch, three metrics: the fraction of reads mapped to endogenous RNA molecules, the number of genes with more than 10 reads per million, and the fraction of reads mapped to mitochondrial genes (<xref rid="bib36" ref-type="bibr">Stegle et al., 2015</xref>
). One 32-cell stage sample failed our quality check and was excluded from downstream analysis (<xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
; <xref rid="figs1" ref-type="fig">Figure S1</xref>
).</p>
<p>Extrinsic “spike-in” RNA-molecules were added to each cell’s lysate prior to cDNA conversion, amplification, and library preparation (<xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
). We applied principal component analysis (PCA) to the log-transformed normalized read counts from the spike-ins and, as expected, did not observe any structure in the data (<xref rid="figs2" ref-type="fig">Figure S2</xref>
). When we applied PCA to the log-transformed normalized read counts of highly variable genes from all cells analyzed, the cells grouped by developmental stage (<xref rid="fig1" ref-type="fig">Figure 1</xref>
B; <xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
), in accord with previous results (<xref rid="bib8" ref-type="bibr">Deng et al., 2014</xref>
, <xref rid="bib40" ref-type="bibr">Tang et al., 2011</xref>
, <xref rid="bib45" ref-type="bibr">Xue et al., 2013</xref>
).</p>
<p>To identify genes that are differentially expressed as development progresses, we analyzed temporal gene expression patterns throughout pre-implantation development. At a false discovery rate (FDR) cutoff of 0.05, 2,716 genes were differentially expressed between 2- and 4-cell, 3,222 between 4- and 8-cell, 2,030 between 8- and 16-cell, and 351 between 16- and 32-cell embryos. Consequently, our data suggest that the transition from the 4- to 8-cell stage is where most transcriptional changes occur. To quantitatively assess intra- and inter-embryo heterogeneity in gene expression, we focused on the stages where all individual cells could be collected from the same embryo (the 2-, 4-, and 8-cell stages). We computed the correlation of expression levels between all pairs of blastomeres within each complete embryo (intra-embryonic correlation) and compared this to the correlation of blastomeres between different embryos (inter-embryonic correlation). We found that intra-embryonic correlation was significantly higher than the inter-embryonic correlation for 2- and 4-cell embryos, while at the 8-cell stage no statistically significant difference could be identified (<xref rid="fig1" ref-type="fig">Figure 1</xref>
C; Welch’s t test, p values are 8.5 × 10<sup>−5</sup>
, 5.9 × 10<sup>−7</sup>
and 0.74 for the 2-, 4-, and 8-cell stage, respectively). Thus at the 2- and 4-cell stages inter-embryonic variability is higher than the transcriptional differences between cells of the same embryo. Conversely, at the 8-cell stage the differences between cells in the same embryo are comparable to those between cells from different embryos, suggesting a relatively higher degree of intra-embryonic heterogeneity. After correcting for inter-embryo differences in expression levels, we identified genes displaying significantly more intra-embryonic heterogeneity in expression than expected by chance (<xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
). These analyses revealed 659, 1,339, and 813 highly variable genes at 2-, 4-, and 8-cell embryos, respectively (<xref rid="fig1" ref-type="fig">Figures 1</xref>
D and <xref rid="figs3" ref-type="fig">S3</xref>
; <xref rid="mmc2" ref-type="supplementary-material">Data S1</xref>
). Of these variable genes, 47 showed heterogeneous expression at all three stages (<xref rid="mmc1" ref-type="supplementary-material">Table S2</xref>
), 82 are shared between the 2- and the 4-cell stage, 203 between the 4- and the 8-cell stage, and 31 between the 2- and the 8-cell stage (<xref rid="fig1" ref-type="fig">Figure 1</xref>
D).</p>
<p>Previous studies have shown that embryo development can differ depending on the orientation and order of cell divisions to the 4-cell stage (<xref rid="bib10" ref-type="bibr">Gardner, 2002</xref>
, <xref rid="bib30" ref-type="bibr">Piotrowska-Nitsche and Zernicka-Goetz, 2005</xref>
), leading us to characterize gene expression in embryos with different division patterns (<xref rid="figs4" ref-type="fig">Figure S4</xref>
). PCA analysis of 16 4-cell embryos revealed that embryos with distinct cell division patterns align approximately along the first principal component (<xref rid="fig1" ref-type="fig">Figure 1</xref>
E; <xref rid="figs4" ref-type="fig">Figure S4</xref>
). Thus, cell division orientation and order might be one of the factors contributing to gene expression patterns in the 4-cell embryo. However, here we focus on genes expressed highly heterogeneously between cells in all individual embryos, independent of their cell division pattern.</p>
</sec>
<sec id="sec2.2"><title><italic>Sox21</italic>
mRNA Expression Is Highly Variable at the 4-Cell Stage and Correlates with the Expression of Pluripotency-Related Genes</title>
<p>We reasoned that highly heterogeneous genes in the 4-cell embryo were of particular interest as cells at this stage can display differential fate (<xref rid="bib31" ref-type="bibr">Piotrowska-Nitsche et al., 2005</xref>
, <xref rid="bib3" ref-type="bibr">Bischoff et al., 2008</xref>
, <xref rid="bib32" ref-type="bibr">Plachta et al., 2011</xref>
, <xref rid="bib38" ref-type="bibr">Tabansky et al., 2013</xref>
) and potential (<xref rid="bib31" ref-type="bibr">Piotrowska-Nitsche et al., 2005</xref>
, <xref rid="bib26" ref-type="bibr">Morris et al., 2012</xref>
). One of the most highly heterogeneous genes in all embryos analyzed at the 4-cell stage is the gene encoding the transcription factor Sox21 (<xref rid="fig2" ref-type="fig">Figure 2</xref>
A), which is involved in regulating ES cell-fate downstream of Sox2 (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
, <xref rid="bib22" ref-type="bibr">Mallanna et al., 2010</xref>
). Sox21 has not been previously studied in the early mouse embryo but is known to inhibit expression of the trophectoderm (TE) master gene Cdx2 in ES cells and is important for reprogramming (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
). Furthermore, Sox21 expression is directly regulated by Sox2 (<xref rid="bib7" ref-type="bibr">Chakravarthy et al., 2011</xref>
, <xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
, <xref rid="bib22" ref-type="bibr">Mallanna et al., 2010</xref>
), and its regulatory region is bound by Oct4 (<xref rid="bib7" ref-type="bibr">Chakravarthy et al., 2011</xref>
, <xref rid="bib11" ref-type="bibr">Göke et al., 2011</xref>
), which has heterogeneous nuclear-cytoplasmic kinetics at the 4-cell stage (<xref rid="bib32" ref-type="bibr">Plachta et al., 2011</xref>
). Together, this suggests that the heterogeneous Sox21 expression may itself be regulated by heterogeneous Oct4 and/or Sox2 activity in the embryo. To explore whether Oct4 and/or Sox2 could drive heterogeneity in gene expression, we next determined whether their target genes are overrepresented in the “highly heterogeneous” group of genes, and whether their expression correlates with that of <italic>Sox21</italic>
, which would be expected if they are regulated by differential Oct4 and/or Sox2 activity in a cell-specific manner.</p>
<p>We applied Fisher’s exact test to a list of Oct4/Sox2 downstream target genes (<xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
) and found that they are indeed overrepresented in the set of highly variable genes at the 4-cell stage (p value 3 × 10<sup>−10</sup>
; <xref rid="fig2" ref-type="fig">Figure 2</xref>
A′). Furthermore, we found that a subset of highly variable Oct4/Sox2 target genes are correlated with Sox21 (<xref rid="fig2" ref-type="fig">Figures 2</xref>
B and <xref rid="figs5" ref-type="fig">S5</xref>
), suggesting coordinated expression of this module of genes at this stage. Interestingly, this module includes known pluripotency-related genes, such as <italic>Nanog</italic>
and <italic>Esrrb</italic>
, suggesting that <italic>Sox21</italic>
expression might be associated with a “more pluripotent” transcriptional state.</p>
<p>When we examined the mRNA expression of <italic>Sox21</italic>
through time, we observed a peak at the 4-cell stage followed by downregulation from the 8-cell stage onward (<xref rid="fig3" ref-type="fig">Figure 3</xref>
A). By contrast, the expression of <italic>Esrrb</italic>
and <italic>Nanog</italic>
increases at the 8- and 16-cell stages, respectively (<xref rid="fig3" ref-type="fig">Figure 3</xref>
A). The expression pattern of <italic>Sox21</italic>
mRNA was consistently highly variable, with a large difference in expression between the highest and lowest expressing cells across all 4-cell embryos examined (<xref rid="fig3" ref-type="fig">Figure 3</xref>
B; n = 64 cells, 16 embryos). Interestingly, in all cases there was at least one cell in each embryo with very low or no <italic>Sox21</italic>
mRNA (<xref rid="fig3" ref-type="fig">Figure 3</xref>
C). This cellular heterogeneity in Sox21 expression was confirmed at the protein level at both the late 4- and early 8-cell stages (<xref rid="fig3" ref-type="fig">Figure 3</xref>
D; n = 25 4-cell embryos and n = 23 8-cell embryos). These results, together with the role of Sox21 in regulating ES cell fate (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
), led us to investigate whether heterogeneous Sox21 expression might have functional consequences for cell-fate specification in the embryo.</p>
</sec>
<sec id="sec2.3"><title>Live-Cell Tracking Demonstrates Cells with Decreased Sox21 Expression Contribute More Extra-Embryonic than Pluripotent Progeny</title>
<p>To determine whether heterogeneity in Sox21 expression influences cell fate, we decreased Sox21 expression using a combination of three <italic>Sox21</italic>
small interfering RNAs (siRNAs), which we confirmed reduce Sox21 protein to undetectable levels (<xref rid="fig4" ref-type="fig">Figure 4</xref>
A). To generate Sox21 heterogeneity, we injected these siRNAs into single blastomeres at the late 2-cell stage alongside mRNA for GFP as a lineage marker to follow a cell’s developmental fate (<xref rid="fig4" ref-type="fig">Figure 4</xref>
B). We allowed the embryos to develop for 72 hr, until the late blastocyst stage, and scored the lineage contribution of each cell using molecular markers for each lineage and cell position within the embryo (<xref rid="fig4" ref-type="fig">Figure 4</xref>
C). Cells with decreased Sox21 expression were significantly more likely to initiate differentiation and contribute to the extra-embryonic TE lineage (p < 0.001) than cells injected with control siRNA, and correspondingly significantly less likely to contribute to pluripotent inner cell mass (ICM) lineages (p < 0.01, <xref rid="fig4" ref-type="fig">Figure 4</xref>
D, n = 75 and n = 38 embryos for <italic>Sox21</italic>
siRNA and control, respectively). When we injected each <italic>Sox21</italic>
siRNA individually, the same phenotype was observed (<xref rid="figs6" ref-type="fig">Figure S6</xref>
), indicating that the phenotype is specific to decreased Sox21 expression and not due to off-target effects.</p>
<p>Segregation of the TE and ICM occurs primarily as a result of asymmetric cell divisions at the 8–16 and 16–32 cell transitions, which contribute one daughter cell to the TE and one to the ICM (<xref rid="bib9" ref-type="bibr">Fleming, 1987</xref>
, <xref rid="bib15" ref-type="bibr">Johnson and Ziomek, 1981</xref>
, <xref rid="bib25" ref-type="bibr">Morris et al., 2010</xref>
, <xref rid="bib43" ref-type="bibr">Watanabe et al., 2014</xref>
). Inside cells can also be generated by apical constriction or cell engulfment (<xref rid="bib24" ref-type="bibr">McDole et al., 2011</xref>
, <xref rid="bib26" ref-type="bibr">Morris et al., 2012</xref>
, <xref rid="bib27" ref-type="bibr">Parfitt and Zernicka-Goetz, 2010</xref>
, <xref rid="bib33" ref-type="bibr">Samarage et al., 2015</xref>
, <xref rid="bib46" ref-type="bibr">Yamanaka et al., 2010</xref>
). To determine the cellular mechanism by which cells with decreased Sox21 expression contribute to the TE, we injected one blastomere of 2-cell embryos with three <italic>Sox21</italic>
siRNAs, or control siRNA, together with mRNA for a GFP-tagged membrane-bound protein (Gap43-GFP) as a marker, and filmed development of the embryos by time-lapse microscopy (<xref rid="fig4" ref-type="fig">Figures 4</xref>
E and 4F; n = 18 embryos, 288 cells for control siRNA and n = 16 embryos, 256 cells for <italic>Sox21</italic>
siRNA). Live-cell tracking of all Gap43-GFP-expressing cells revealed that, by the early blastocyst (32-cell) stage, control siRNA blastomeres had generated an average of 10.9 outside cells and 5.1 inside cells (<xref rid="fig4" ref-type="fig">Figure 4</xref>
G), with an average of 4.41 asymmetric divisions and 0.69 cells internalized by cell engulfment (<xref rid="fig4" ref-type="fig">Figures 4</xref>
H and 4I). In contrast, the <italic>Sox21</italic>
siRNA blastomeres had generated an average of 13.7 outside cells and only 2.3 inside cells, with significantly fewer asymmetric divisions (an average of 2.3) by the same stage (p < 0.001; <xref rid="fig4" ref-type="fig">Figures 4</xref>
G, 4H, and 4J). Therefore, cells with lower Sox21 expression undertake fewer asymmetric cell divisions and contribute more extra-embryonic rather than pluripotent progeny.</p>
</sec>
<sec id="sec2.4"><title>Decreased Expression of Sox21 Results in Premature and Increased Expression of Cdx2</title>
<p>Correct differentiation into extra-embryonic TE requires the transcription factor Cdx2 (<xref rid="bib14" ref-type="bibr">Jedrusik et al., 2008</xref>
, <xref rid="bib37" ref-type="bibr">Strumpf et al., 2005</xref>
). As Sox21 has been shown to bind to an enhancer of <italic>Cdx2</italic>
and repress its expression in ES cells (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
), we next investigated whether Sox21 might also regulate <italic>Cdx2</italic>
in the embryo. To this end, we first decreased expression of Sox21 at the zygote stage, by injecting a mixture of three siRNAs as before, and examined the developmental timing and level of <italic>Cdx2</italic>
expression. Decreased expression of Sox21 resulted in a significant increase in <italic>Cdx2</italic>
mRNA expression at the 8-cell stage, compared to control siRNA embryos (<xref rid="fig5" ref-type="fig">Figures 5</xref>
A and 5B; p < 0.05; n = 85 and n = 75 embryos for Sox21 siRNA and control, respectively, three biological replicates). When we repeated the experiment by downregulating Sox21 in the zygote but now examined Cdx2 protein at the early 8-cell stage, when Cdx2 protein is normally not yet detectable, we found that all embryos in which Sox21 expression was decreased already expressed Cdx2 protein (<xref rid="fig5" ref-type="fig">Figures 5</xref>
C and 5D; n = 96 cells [84 Cdx2-positive cells, 87.5%], 12 embryos for the experimental group and n = 88 cells [11 Cdx2-positive cells, 12.5%], 11 embryos for the control group).</p>
<p>Finally, to test the functional consequences of heterogeneous Sox21 expression on Cdx2 expression, we decreased Sox21 in just half of the embryo, from the 2-cell stage onward, using Gap43-RFP as a lineage marker (<xref rid="fig5" ref-type="fig">Figures 5</xref>
E and 5F). This resulted in a significant upregulation of Cdx2 in the half of the embryo in which cells had decreased Sox21 expression (<xref rid="fig5" ref-type="fig">Figure 5</xref>
G; p < 0.001; n = 104 cells, 13 embryos and n = 88 cells, 11 embryos for the experimental and control groups, respectively). Together, these results suggest that cells with lower Sox21 levels are the first to upregulate Cdx2 (mRNA and protein) and therefore are first to initiate development into the extra-embryonic TE lineage.</p>
</sec>
<sec id="sec2.5"><title>Sox21 Expression Is Regulated by CARM1 Activity</title>
<p>Since our results indicate that heterogeneity in Sox21 expression could bias cell fate in the early embryo, we sought to determine the upstream regulators of this heterogeneity. One factor we previously found to be heterogeneous at the 4-cell stage is histone H3R26 methylation mediated by CARM1, which itself is heterogeneously active at the 4-cell stage (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
). Furthermore, elevated CARM1 expression leads to increased histone H3R26me and contribution to the ICM (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
, <xref rid="bib44" ref-type="bibr">Wu et al., 2009</xref>
). To determine whether altering CARM1 activity affects the expression of Sox21, we first treated embryos with a CARM1-specific inhibitor from the 2-cell stage. To test its effectiveness, we determined the levels of histone H3R26me after 10 hr of treatment and observed it to be reduced to an undetectable level (<xref rid="fig6" ref-type="fig">Figure 6</xref>
A). Inhibition of CARM1 resulted in a complete loss of Sox21 expression at the 4-cell stage (<xref rid="fig6" ref-type="fig">Figure 6</xref>
B, n = 21 control; n = 15 CARM1 inhibition). To downregulate CARM1 through an alternative mechanism, we injected <italic>Carm1</italic>
siRNA into zygotes and cultured embryos until the blastocyst stage. Decreasing CARM1 expression resulted in a significant reduction in the number of pluripotent epiblast cells with a concurrent increase in the number of primitive endoderm cells and no effect on the TE lineage (<xref rid="fig6" ref-type="fig">Figures 6</xref>
C and 6D p < 0.001; n = 12 control; n = 8 <italic>Carm1</italic>
siRNA). To overexpress CARM1, we injected synthetic mRNA for CARM1 into zygotes. CARM1 upregulation led to a significant upregulation of <italic>Sox21</italic>
mRNA expression at the 8-stage, as well as increased expression of other pluripotency genes, including <italic>Sox2</italic>
and <italic>Nanog</italic>
(<xref rid="fig6" ref-type="fig">Figure 6</xref>
E; p < 0.01 and p < 0.05; n = 75 embryos and n = 66 embryos for CARM1 mRNA and control, respectively, three biological replicates). These results indicate that CARM1 activity regulates expression of a number of pluripotency regulators including Sox21 in the 4-cell embryo.</p>
</sec>
</sec>
<sec id="sec3"><title>Discussion</title>
<p>When heterogeneity in gene expression patterns first arises during the development of the mammalian embryo and how this might influence cell fate have remained, thus far, unknown. Approaching these questions has been difficult because of the inherent developmental flexibility of the embryo and because of previous technical limitations that prevented accurate quantification of gene expression levels in individual cells as development progresses. Although it has previously been thought that cells are homogenous until generation of inside and outside cells, more recent studies have opened a new possibility that cells at the 4-cell stage can already exhibit differences in cell fate (<xref rid="bib31" ref-type="bibr">Piotrowska-Nitsche et al., 2005</xref>
, <xref rid="bib38" ref-type="bibr">Tabansky et al., 2013</xref>
) and developmental potential (<xref rid="bib26" ref-type="bibr">Morris et al., 2012</xref>
, <xref rid="bib31" ref-type="bibr">Piotrowska-Nitsche et al., 2005</xref>
). This led to the suggestion that the differential activity of epigenetic regulators such as CARM1 (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
) or PRDM14 (<xref rid="bib6" ref-type="bibr">Burton et al., 2013</xref>
), or the differential behavior of transcription factors such as Oct4 (<xref rid="bib32" ref-type="bibr">Plachta et al., 2011</xref>
), at the 4-cell stage could be linked to differential cell fate and potential. Here, we sought to address this by first investigating and directly comparing global differences in gene expression between individual cells within the embryo during the first days of development. Focusing on the transcriptomes of each and every cell in 16 4-cell stage embryos uncovered consistent patterns of highly heterogeneous gene expression, with the most heterogeneously expressed genes being targets of Sox2 and Oct4, among which is the gene for the transcription factor Sox21.</p>
<p>The role of Sox21 in the early embryo has not been previously investigated, but studies in ES cells indicated its involvement in stem cell differentiation, the inhibition of Cdx2 expression, and cellular reprogramming (<xref rid="bib22" ref-type="bibr">Mallanna et al., 2010</xref>
, <xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
). The regulatory regions of the <italic>Sox21</italic>
gene are targets of Oct4 binding, and its expression is regulated by Sox2 (<xref rid="bib7" ref-type="bibr">Chakravarthy et al., 2011</xref>
, <xref rid="bib11" ref-type="bibr">Göke et al., 2011</xref>
, <xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
, <xref rid="bib22" ref-type="bibr">Mallanna et al., 2010</xref>
). Among other target genes of Oct4 and/or Sox2 showing similar heterogeneity of expression, we identified those for pluripotency transcription factors such as Nanog and Esrrb; however, their expression peaked later than that of <italic>Sox21</italic>
, suggesting that Sox21 could play a role very early in development. Indeed, our finding that expression of Sox21 represses premature differentiation in the embryo accords with the finding that Sox21 is essential for reprogramming in vitro (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
). However, <italic>Sox21</italic>
mutant mice are viable (<xref rid="bib17" ref-type="bibr">Kiso et al., 2009</xref>
). One possible explanation of this paradox could be that Sox21 is part of a functionally redundant mechanism typical of crucial biological processes and of importance in the regulative development of mammalian embryos. Heterogeneous expression of <italic>Sox21</italic>
, along with other co-regulated genes, could also set up a competitive advantage between blastomeres during early development without being essential. Thus, while embryos may not have an absolute requirement for Sox21, heterogeneous Sox21 expression levels could lead to biased cell-fate decisions. We tested this possibility by setting up clonal differences in Sox21 levels experimentally. This revealed that clones of cells with decreased Sox21 expression contribute more differentiated TE rather than pluripotent progeny via a decreased frequency of asymmetric cell divisions. It has been previously shown that the frequency of symmetric to asymmetric cell divisions can be altered by the expression of the TE-specific transcription factor Cdx2 (<xref rid="bib14" ref-type="bibr">Jedrusik et al., 2008</xref>
). In agreement, we find that decreased Sox21 expression leads to elevated expression of Cdx2 at the 8-cell stage. This can be explained by noting that Sox21 binds to the <italic>Cdx2</italic>
enhancer in stem cells to negatively regulate its expression (<xref rid="bib18" ref-type="bibr">Kuzmichev et al., 2012</xref>
). In addition, depletion of Sox21 in the whole embryo leads to earlier expression of Cdx2. Thus, in the case of uniform Sox21 depletion, there are no heterogeneities present and thus no bias in cell fate. These results suggest that Sox21 expression in the 4-cell embryo might help to safeguard against the premature initiation of differentiation.</p>
<p>How might the heterogeneous expression of Sox21 first arise? It is possible that the Sox21 heterogeneity we identify here is generated by random transcriptional differences; however, we wished to investigate whether this differential expression could be regulated by factors previously found to be heterogeneous at the 4-cell stage. Our best candidate for the heterogeneous regulation of Sox21 was methylation of histone H3R26 because we have previously found that CARM1, which mediates this modification, regulates the expression of pluripotency genes and is differentially active in the 4-cell embryo (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
). We find here that inhibiting CARM1 activity results in a loss of Sox21 expression in 4-cell embryos and, in contrast, its upregulation causes an increase in expression of <italic>Sox21</italic>
and the pluripotency genes <italic>Sox2</italic>
and <italic>Nanog</italic>
at the 8-cell stage. Therefore, differential activity of CARM1 in the 4-cell embryo is likely to be responsible for the differential expression of Sox21. We also show that by decreasing CARM1 expression we significantly reduce the number of pluripotent epiblast cells, indicating the critical role of CARM1 for maintaining the pluripotent state in the ICM. How CARM1 becomes heterogeneous remains an open question. The work presented here shows that CARM1 overexpression does not affect <italic>Oct4</italic>
mRNA expression, suggesting that the regulation of <italic>Sox21</italic>
by CARM1 is not achieved through transcriptional regulation of <italic>Oct4</italic>
. This, however, does not exclude the possibility that CARM1 is indirectly affecting Oct4 or Sox2 activity. It can be hypothesized that increased levels of CARM1, resulting in increased histone H3R26me, could influence Oct4 and/or Sox2 DNA binding and lead to the upregulation of Oct4 and/or Sox2 target genes (<xref rid="fig7" ref-type="fig">Figure 7</xref>
). Indeed, Oct4 expression is not heterogeneous at the 4-cell stage and its different DNA binding dynamics in individual cells might be due to differential accessibility to target genes (<xref rid="bib32" ref-type="bibr">Plachta et al., 2011</xref>
).</p>
<p>In conclusion, we suggest that highly heterogeneous expression of several Oct4 and/or Sox2 targets within the 4-cell embryo, regulated by differential activity of CARM1, creates heterogeneities that bias cell fate toward either an embryonic (pluripotent) or extra-embryonic (differentiating) fate. We further suggest that the highly heterogeneous expression of Sox21 is one of a number of mechanisms available to the embryo to direct cell fate. Our findings that several Oct4 and Sox2 target genes are overrepresented in the subset of highly variable genes at the 4-cell stage and that Sox21 is co-regulated with other transcription factors, such as Nanog and Esrrb, supports this hypothesis. Thus, if a single source of heterogeneity were to be removed, the embryo could compensate for this loss. In summary, our results indicate that heterogeneous gene expression as early as the 4-cell stage initiates cell-fate decisions.</p>
</sec>
<sec id="sec4"><title>Experimental Procedures</title>
<sec id="sec4.1"><title>Collection of Mouse Embryos and Cells</title>
<p>All animal experiments were performed in compliance with Home Office regulations. Embryos were recovered from superovulated F1 (C57Bl6xCBA) females into M2 media supplemented with 4% BSA, as described previously (<xref rid="bib29" ref-type="bibr">Piotrowska et al., 2001</xref>
). Individual cells were collected at late 2-cell (48 hr after hCG), late 4-cell (10 hr after the first 2-cell blastomere divided), late 8-cell (10 hr after the first 4-cell blastomere divided), 16-cell (78 hr after hCG), and 32-cell (86 hr after hCG) stages. Embryos for cell isolation at the 2-, 4,- and 8-cell stages were collected at the 2-cell stage and cultured in KSOM media (Millipore) at 37°C with 5% CO<sub>2</sub>
. Embryos for the 16- and 32-cell blastomeres were collected directly at the respective stages. The zona pellucida was removed using Tyrode’s solution (Sigma). Zona-free embryos were incubated for 5 min (for the 2-, 4-, and 8-cell stages) or 20 min (for the 16- and 32-cell stages) in Ca<sup>2+</sup>
and Mg<sup>2+</sup>
free M2 before disaggregation by careful pipetting. Each embryo was processed individually. Single blastomeres were placed into individual tubes containing 2.3 μl of 0.2% Triton X-100 (Sigma) supplemented with 1 U/μl RNAsIN (Ambion).</p>
</sec>
<sec id="sec4.2"><title>Scoring 2- to 4-Cell Division Pattern</title>
<p>At the late 2-cell stage, one cell was injected with rhodamine-dextran (Invitrogen), as described previously (<xref rid="bib3" ref-type="bibr">Bischoff et al., 2008</xref>
). Embryos were observed at 20 min intervals to categorize the division order and orientation, meridional (M) or equatorial (E), relative to the position of the second polar body. 4-cell embryos were sorted into ME, EM, EE, and MM groups. Embryos with absent or unattached polar body were discarded. Following disaggregation, the cells were inspected under an epifluorescent microscope and classified as originating from the first- or second-dividing 2-cell blastomere based on the rhodamine fluorescence.</p>
</sec>
<sec id="sec4.3"><title>RNA Sequencing and Mapping of Reads</title>
<p>mRNA from the single cells was amplified using the SMARTSeq2 protocol (<xref rid="bib28" ref-type="bibr">Picelli et al., 2014</xref>
), with the additional inclusion of ERCC spike-in control (1 μl of 1:250,000 [batches 1 and 2] or a 1:1,000,000 [batch 3 and batch 4] dilution of mix 1 [Ambion] per cell). Multiplex sequencing libraries were generated from amplified cDNA using Nextera XT (Illumina) and sequenced on a HiSeq 2500 running in rapid mode. Paired-end reads were mapped simultaneously to the <italic>M. musculus</italic>
genome (Ensembl v.38.77) and the ERCC sequences using GSNAP (v.2014-10-07) with default parameters. Htseq-count (<xref rid="bib1" ref-type="bibr">Anders et al., 2015</xref>
) was used to count the number of reads mapped to each gene (default options).</p>
</sec>
<sec id="sec4.4"><title>Quality Assessment of Cells</title>
<p>To assess the quality of the data, three metrics were used: the fraction of mapped reads, the number of genes with more than 10 reads per million (RPM), and the fraction of reads mapped to mitochondrial genes. <xref rid="figs1" ref-type="fig">Figures S1</xref>
A–S1C show these metrics as functions of the total number of reads for each cell. When PCA was carried out using these three quantities, cells that perform worse than the average according to all three criteria appear as outliers (<xref rid="figs1" ref-type="fig">Figure S1</xref>
D). One clear outlier was found (black arrow in all four panels) and was excluded from downstream analyses (sample name “32cell_F”, see <xref rid="mmc1" ref-type="supplementary-material">Table S1</xref>
for a list of samples).</p>
</sec>
<sec id="sec4.5"><title>Normalization of Read Counts and Analysis of Spike-ins</title>
<p>The data were divided into two groups according to whether they were mapped to ERCC spike-ins or to endogenous genes. For each of these groups a set of size factors (<xref rid="bib21" ref-type="bibr">Love et al., 2014</xref>
) was calculated and used to normalize the raw read counts. While the size factors calculated with ERCC spike-ins account for sequencing depth, the size factors calculated on endogenous genes also normalize for the amount of RNA obtained from each cell (<xref rid="bib4" ref-type="bibr">Brennecke et al., 2013</xref>
), which is highly correlated with cell cycle stage (<xref rid="bib5" ref-type="bibr">Buettner et al., 2015</xref>
).</p>
</sec>
<sec id="sec4.6"><title>Differential Expression Analysis</title>
<p>First, genes with an average expression level of less than 50 normalized read counts across the two groups of cells compared were removed. The Bioconductor package DESeq2 (<xref rid="bib21" ref-type="bibr">Love et al., 2014</xref>
) was then used to find sets of differentially expressed genes. An FDR of 0.05 was used as a threshold for significance. To avoid confounding effects due to embryos being collected in different batches, when comparing blastomeres from different stages, we used only embryos at 2-, 4-, 8-, 16-, and 32-cell stages from the same batch (batch 1, see <xref rid="mmc1" ref-type="supplementary-material">Table S1</xref>
for a list of samples).</p>
</sec>
<sec id="sec4.7"><title>Highly Variable Genes</title>
<p>To identify highly variable genes across all cells collected, we applied the method described in <xref rid="bib4" ref-type="bibr">Brennecke et al. (2013)</xref>
by fitting the squared coefficient of variation CV<sup>2</sup>
as a function of the mean normalized counts μ with the parameterization CV<sup>2</sup>
= a<sub>1</sub>
/μ + α<sub>0</sub>
. To minimize the skewing effect due to lowly expressed genes, only genes with a mean normalized count greater than 10 were used. Genes with an adjusted p value (Benjamini-Hochberg method) less than 0.1 were considered as significantly highly variable.</p>
</sec>
<sec id="sec4.8"><title>Adjusting for Batch Effects</title>
<p>The technical variation between samples in different batches was removed using the ComBat function in the R package “sva” with default options (<xref rid="bib16" ref-type="bibr">Johnson et al., 2007</xref>
). A log<sub>10</sub>
transformation was first applied to normalized read counts (by summing 1 to avoid infinities) and lowly expressed genes (less than 10 normalized read counts on average) were removed. ComBat was used to adjust for batch effect with the known batch covariate, controlling for sample stage and cell division pattern. All expression values that were 0 or negative after batch effect removal were set to 0. These batch-adjusted expression values were used to perform PCA on the highly variable genes. The results do not significantly change if different criteria are adopted to select genes (e.g., union of 1,000 or 5,000 most highly expressed genes, 3,000 genes with highest average expression).</p>
</sec>
<sec id="sec4.9"><title>Analysis of Inter- and Intra-embryonic Heterogeneity</title>
<p>To identify genes that show a significantly high variability across blastomeres in the same embryo, we first regressed out the embryo effect by fitting the linear model:<disp-formula id="ufd1"><mml:math id="M1" altimg="si1.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mi>g</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>B</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mi>ε</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
where <inline-formula><mml:math id="M2" altimg="si2.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mi>g</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
is the log-expression value of the gene <italic>i</italic>
in cell <italic>j</italic>
and embryo <italic>k</italic>
, <italic>A</italic>
<sub><italic>i</italic>
</sub>
is the average log-expression of gene <italic>i</italic>
across all cells, <italic>B</italic>
<sub><italic>ik</italic>
</sub>
represents the embryo effect and <italic>ε</italic>
<sub><italic>ijk</italic>
</sub>
is the fit residual. The expression level of the gene with the embryo effect removed, <inline-formula><mml:math id="M3" altimg="si3.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mover accent="true"><mml:mi>g</mml:mi>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
, is equal to <inline-formula><mml:math id="M4" altimg="si4.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mi>g</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msub><mml:mi>B</mml:mi>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
. We then applied the method described above to the <inline-formula><mml:math id="M5" altimg="si3.gif" overflow="scroll"><mml:mrow><mml:msub><mml:mrow><mml:mover accent="true"><mml:mi>g</mml:mi>
<mml:mo>ˆ</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
to identify highly variable genes.</p>
<p>The intra- and inter-embryonic transcriptional differences were quantified and compared between blastomeres by computing Spearman’s correlation coefficient between all pairs of blastomeres within the same embryo and in different embryos at 2-, 4-, and 8-cell stage. Genes with a mean normalized count less than 10 were removed. Only embryos collected in the same batch were used for the calculation of inter-embryonic correlations (<xref rid="mmc1" ref-type="supplementary-material">Table S1</xref>
for a list of samples).</p>
</sec>
<sec id="sec4.10"><title>Testing for Enrichment of Oct4 and Sox2 Targets among Highly Variable Genes</title>
<p>We downloaded a list of Oct4 target genes identified in three different publications (<xref rid="bib20" ref-type="bibr">Loh et al., 2006</xref>
, <xref rid="bib23" ref-type="bibr">Matoba et al., 2006</xref>
, <xref rid="bib34" ref-type="bibr">Sharov et al., 2008</xref>
) and selected only those that were targeted by Oct4 in at least two out of three publications. After removing all genes that were never expressed in our data, we tested the significance of the association between the Oct4 target genes and highly variable genes in 2-, 4-, and 8-cell embryos by Fisher’s exact test (p values are, respectively, 0.005, 10<sup>−9</sup>
and 4 × 10<sup>−6</sup>
). Analogously, we tested the enrichment of Sox2 target genes (<xref rid="bib34" ref-type="bibr">Sharov et al., 2008</xref>
) among the highly variable genes and found statistical significance at the 2- and 4-cell stages (p values 0.004 and 10<sup>−5</sup>
), but not at the 8-cell stage (p value 0.26). Finally, the enrichment test run on genes targeted by Oct4 and/or Sox2 gave the following p values for embryos at 2-, 4-, and 8-cell stage, respectively: 2 × 10<sup>−4</sup>
, 3 × 10<sup>−10</sup>
, and 3 × 10<sup>−5</sup>
.</p>
<p>We focused on genes that are highly variable at the 4-cell stage and are targeted by Oct4 and/or Sox2 (including <italic>Sox21</italic>
[<xref rid="bib7" ref-type="bibr">Chakravarthy et al., 2011</xref>
, <xref rid="bib11" ref-type="bibr">Göke et al., 2011</xref>
]) and searched for subsets of genes whose expression patterns are coordinated. To this aim, the matrix of pairwise Spearman correlation coefficients was computed and hierarchical clustering was carried out (R function “hclust,” default options). The Dynamic Tree Cut algorithm was used to identify clusters of genes (<xref rid="bib19" ref-type="bibr">Langfelder et al., 2008</xref>
) with the R function “cutreeDynamic” in the “dynamicTreeCut” package (default parameters).</p>
</sec>
<sec id="sec4.11"><title>Functional Assays</title>
<p>To downregulate Sox21, three <italic>Sox21</italic>
-specific siRNAs (QIAGEN) at a total (combined or individual) concentration of 12 μM were injected together with <italic>GFP</italic>
, <italic>Gap43-GFP</italic>
, or <italic>Gap43-RFP</italic>
mRNA (400 ng/ml). AllStars Negative Control siRNA (QIAGEN) was used as a control. The sequences of the siRNAs used are as follows: Sox21 siRNA 1: 5′-CCCGGTTTGTATGTACATAGA-3′, Sox21 siRNA 2: 5′-TACTCTGATTGTACTGTTGAA-3′, Sox21 siRNA 3: 5′-TTGTATGTACATAGATGTATA-3′. To inhibit CARM1 activity embryos were treated from the 2- to 4-cell stage with a chemical inhibitor against CARM1 (Millipore, catalog no. 217531). CARM1 inhibitor was dissolved in DMSO (final concentration of DMSO was 0.005%) and used at a concentration of 9 μM in KSOM. Control embryos were incubated in the equivalent DMSO concentration but in the absence of the inhibitor. To downregulate CARM1 expression, CARM1 stealth siRNAs (Life Technologies - AM16708) were injected at a concentration of 200 nM together with <italic>Gap43-RFP</italic>
mRNA (200 ng/ml). CARM1 mRNA was prepared and microinjected as previously described (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
) at a concentration of 500 ng/μl. Embryos were cultured in KSOM under paraffin oil at 37.5°C in a 5% CO<sub>2</sub>
.</p>
</sec>
<sec id="sec4.12"><title>Immunofluorescence</title>
<p>Primary antibodies used were goat anti-Sox17 (R&D Systems), rabbit anti-Nanog (Abcam), mouse anti-Cdx2 (Biogenex), rabbit anti-Histone H3 (Abcam), rabbit anti-Sox2 (Cell Signaling Technology), rabbit anti-Sox21 (R&D Systems), and rabbit anti-Histone H3 (symmetric di methyl R26; Abcam). ME embryos were used for Sox21 protein quantification at the 4- and 8-cell stages. Immunofluorescence was carried out as described previously (<xref rid="bib14" ref-type="bibr">Jedrusik et al., 2008</xref>
). Multichannel images were acquired using a Leica SP5 confocal microscope, and confocal z stacks were exported to IMAGEJ for image processing, intensity measurements and cell counting. Quantification of fluorescence intensity was done by normalizing to Histone H3 staining and layer-normalizing using the built-in ImageJ function. Intensity was measured in the normalized sections using the ImageJ measure function.</p>
</sec>
<sec id="sec4.13"><title>Live Imaging</title>
<p>Following injection with siRNA, embryos were observed on an inverted epifluorescent Zeiss Axiovert 200M microscope with a 20×/0.75 NA objective. Two-channel (green fluorescence and transmitted light) multisection images were acquired every 30 min with a Hamamatsu ORCA ER CCD camera on 16 focal planes every 6 μm, with an exposure of 10 ms for transmitted light and 50 ms for fluorescence. Simi Biocell software was used for cell tracking, as previously described (<xref rid="bib3" ref-type="bibr">Bischoff et al., 2008</xref>
). Cell identification and lineage rendering was done with MovIT software (from the BioEmergences platform).</p>
</sec>
<sec id="sec4.14"><title>qRT-PCR</title>
<p>Total RNA was extracted from 8-cell stage embryos using the Arcturus PicoPure RNA isolation kit (Arcturus Bioscience). qRT-PCR was carried out using SYBR Green in a StepOne Plus Real-time PCR machine (Applied Biosystems). Gapdh was used as an endogenous control. Relative mRNA concentrations were calculated using the ddCT method. The following primers were used: Gapdh forward, 5′-AGAGACGGCCGCATCTTC-3′ and reverse, 5′-CCCAATACGGCCAAATCCGT-3′; Cdx2 forward, 5′-AAACCTGTGCGAGTGGATG-3′ and reverse, 5′-TCTGTGTACACCACCCGGTA-3′; Sox21 forward, 5′-GCCGGTGACTCGTGTCTTTA-3′ and reverse, 5′-GAACGGCGGTCATCTCTCAT-3′; Nanog forward, 5′-GGTTGAAGACTAGCAATGGTCTGA-3′ and reverse, 5′-TGCAATGGATGCTGGGATACT −3′; Oct3/4 forward, 5′-TTGGGCTAGAGAAGGATGTGGTT-3′ and reverse, 5′-GGAAAAGGGACTGAGTAGAGTGTGG-3′; Sox2 forward, 5′-CCATCCACCCTTATGTATCCAAG-3′ and reverse, 5′-CGAAGGAAGTGGGTAAACAGCAC-3′.</p>
</sec>
</sec>
<sec id="sec5"><title>Author Contributions</title>
<p>M.G. carried out the majority of the experiments with help from S.J.L.G. and A.H. S.J.L.G., M.G., and A.J. scored cell divisions and collected individual cells. A.S. analyzed the single-cell transcriptomics data and performed all downstream analyses using these data. I.C.M. prepared the cDNA libraries and performed and analyzed the deep sequencing data. T.V. supervised the cDNA library preparation. J.C.M. oversaw the computational analysis and supervised the study. M.Z.-G. conceived the project, designed and oversaw experiments, and supervised the study. All authors contributed to writing the manuscript.</p>
</sec>
</body>
<back><ref-list><title>References</title>
<ref id="bib1"><element-citation publication-type="journal" id="sref1"><person-group person-group-type="author"><name><surname>Anders</surname>
<given-names>S.</given-names>
</name>
<name><surname>Pyl</surname>
<given-names>P.T.</given-names>
</name>
<name><surname>Huber</surname>
<given-names>W.</given-names>
</name>
</person-group>
<article-title>HTSeq—a Python framework to work with high-throughput sequencing data</article-title>
<source>Bioinformatics</source>
<volume>31</volume>
<year>2015</year>
<fpage>166</fpage>
<lpage>169</lpage>
<pub-id pub-id-type="pmid">25260700</pub-id>
</element-citation>
</ref>
<ref id="bib2"><element-citation publication-type="journal" id="sref2"><person-group person-group-type="author"><name><surname>Biase</surname>
<given-names>F.H.</given-names>
</name>
<name><surname>Cao</surname>
<given-names>X.</given-names>
</name>
<name><surname>Zhong</surname>
<given-names>S.</given-names>
</name>
</person-group>
<article-title>Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencing</article-title>
<source>Genome Res.</source>
<volume>24</volume>
<year>2014</year>
<fpage>1787</fpage>
<lpage>1796</lpage>
<pub-id pub-id-type="pmid">25096407</pub-id>
</element-citation>
</ref>
<ref id="bib3"><element-citation publication-type="journal" id="sref3"><person-group person-group-type="author"><name><surname>Bischoff</surname>
<given-names>M.</given-names>
</name>
<name><surname>Parfitt</surname>
<given-names>D.E.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Formation of the embryonic-abembryonic axis of the mouse blastocyst: relationships between orientation of early cleavage divisions and pattern of symmetric/asymmetric divisions</article-title>
<source>Development</source>
<volume>135</volume>
<year>2008</year>
<fpage>953</fpage>
<lpage>962</lpage>
<pub-id pub-id-type="pmid">18234722</pub-id>
</element-citation>
</ref>
<ref id="bib4"><element-citation publication-type="journal" id="sref4"><person-group person-group-type="author"><name><surname>Brennecke</surname>
<given-names>P.</given-names>
</name>
<name><surname>Anders</surname>
<given-names>S.</given-names>
</name>
<name><surname>Kim</surname>
<given-names>J.K.</given-names>
</name>
<name><surname>Kołodziejczyk</surname>
<given-names>A.A.</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name><surname>Proserpio</surname>
<given-names>V.</given-names>
</name>
<name><surname>Baying</surname>
<given-names>B.</given-names>
</name>
<name><surname>Benes</surname>
<given-names>V.</given-names>
</name>
<name><surname>Teichmann</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Marioni</surname>
<given-names>J.C.</given-names>
</name>
<name><surname>Heisler</surname>
<given-names>M.G.</given-names>
</name>
</person-group>
<article-title>Accounting for technical noise in single-cell RNA-seq experiments</article-title>
<source>Nat. Methods</source>
<volume>10</volume>
<year>2013</year>
<fpage>1093</fpage>
<lpage>1095</lpage>
<pub-id pub-id-type="pmid">24056876</pub-id>
</element-citation>
</ref>
<ref id="bib5"><element-citation publication-type="journal" id="sref5"><person-group person-group-type="author"><name><surname>Buettner</surname>
<given-names>F.</given-names>
</name>
<name><surname>Natarajan</surname>
<given-names>K.N.</given-names>
</name>
<name><surname>Casale</surname>
<given-names>F.P.</given-names>
</name>
<name><surname>Proserpio</surname>
<given-names>V.</given-names>
</name>
<name><surname>Scialdone</surname>
<given-names>A.</given-names>
</name>
<name><surname>Theis</surname>
<given-names>F.J.</given-names>
</name>
<name><surname>Teichmann</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Marioni</surname>
<given-names>J.C.</given-names>
</name>
<name><surname>Stegle</surname>
<given-names>O.</given-names>
</name>
</person-group>
<article-title>Computational analysis of cell-to-cell heterogeneity in single-cell RNA-sequencing data reveals hidden subpopulations of cells</article-title>
<source>Nat. Biotechnol.</source>
<volume>33</volume>
<year>2015</year>
<fpage>155</fpage>
<lpage>160</lpage>
<pub-id pub-id-type="pmid">25599176</pub-id>
</element-citation>
</ref>
<ref id="bib6"><element-citation publication-type="journal" id="sref6"><person-group person-group-type="author"><name><surname>Burton</surname>
<given-names>A.</given-names>
</name>
<name><surname>Muller</surname>
<given-names>J.</given-names>
</name>
<name><surname>Tu</surname>
<given-names>S.</given-names>
</name>
<name><surname>Padilla-Longoria</surname>
<given-names>P.</given-names>
</name>
<name><surname>Guccione</surname>
<given-names>E.</given-names>
</name>
<name><surname>Torres-Padilla</surname>
<given-names>M.E.</given-names>
</name>
</person-group>
<article-title>Single-cell profiling of epigenetic modifiers identifies PRDM14 as an inducer of cell fate in the mammalian embryo</article-title>
<source>Cell Rep.</source>
<volume>5</volume>
<year>2013</year>
<fpage>687</fpage>
<lpage>701</lpage>
<pub-id pub-id-type="pmid">24183668</pub-id>
</element-citation>
</ref>
<ref id="bib7"><element-citation publication-type="journal" id="sref7"><person-group person-group-type="author"><name><surname>Chakravarthy</surname>
<given-names>H.</given-names>
</name>
<name><surname>Ormsbee</surname>
<given-names>B.D.</given-names>
</name>
<name><surname>Mallanna</surname>
<given-names>S.K.</given-names>
</name>
<name><surname>Rizzino</surname>
<given-names>A.</given-names>
</name>
</person-group>
<article-title>Rapid activation of the bivalent gene Sox21 requires displacement of multiple layers of gene-silencing machinery</article-title>
<source>FASEB J.</source>
<volume>25</volume>
<year>2011</year>
<fpage>206</fpage>
<lpage>218</lpage>
<pub-id pub-id-type="pmid">20876214</pub-id>
</element-citation>
</ref>
<ref id="bib8"><element-citation publication-type="journal" id="sref8"><person-group person-group-type="author"><name><surname>Deng</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Ramsköld</surname>
<given-names>D.</given-names>
</name>
<name><surname>Reinius</surname>
<given-names>B.</given-names>
</name>
<name><surname>Sandberg</surname>
<given-names>R.</given-names>
</name>
</person-group>
<article-title>Single-cell RNA-seq reveals dynamic, random monoallelic gene expression in mammalian cells</article-title>
<source>Science</source>
<volume>343</volume>
<year>2014</year>
<fpage>193</fpage>
<lpage>196</lpage>
<pub-id pub-id-type="pmid">24408435</pub-id>
</element-citation>
</ref>
<ref id="bib9"><element-citation publication-type="journal" id="sref9"><person-group person-group-type="author"><name><surname>Fleming</surname>
<given-names>T.P.</given-names>
</name>
</person-group>
<article-title>A quantitative analysis of cell allocation to trophectoderm and inner cell mass in the mouse blastocyst</article-title>
<source>Dev. Biol.</source>
<volume>119</volume>
<year>1987</year>
<fpage>520</fpage>
<lpage>531</lpage>
<pub-id pub-id-type="pmid">3803716</pub-id>
</element-citation>
</ref>
<ref id="bib10"><element-citation publication-type="journal" id="sref10"><person-group person-group-type="author"><name><surname>Gardner</surname>
<given-names>R.L.</given-names>
</name>
</person-group>
<article-title>Experimental analysis of second cleavage in the mouse</article-title>
<source>Hum. Reprod.</source>
<volume>17</volume>
<year>2002</year>
<fpage>3178</fpage>
<lpage>3189</lpage>
<pub-id pub-id-type="pmid">12456621</pub-id>
</element-citation>
</ref>
<ref id="bib11"><element-citation publication-type="journal" id="sref11"><person-group person-group-type="author"><name><surname>Göke</surname>
<given-names>J.</given-names>
</name>
<name><surname>Jung</surname>
<given-names>M.</given-names>
</name>
<name><surname>Behrens</surname>
<given-names>S.</given-names>
</name>
<name><surname>Chavez</surname>
<given-names>L.</given-names>
</name>
<name><surname>O’Keeffe</surname>
<given-names>S.</given-names>
</name>
<name><surname>Timmermann</surname>
<given-names>B.</given-names>
</name>
<name><surname>Lehrach</surname>
<given-names>H.</given-names>
</name>
<name><surname>Adjaye</surname>
<given-names>J.</given-names>
</name>
<name><surname>Vingron</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Combinatorial binding in human and mouse embryonic stem cells identifies conserved enhancers active in early embryonic development</article-title>
<source>PLoS Comput. Biol.</source>
<volume>7</volume>
<year>2011</year>
<fpage>e1002304</fpage>
<pub-id pub-id-type="pmid">22215994</pub-id>
</element-citation>
</ref>
<ref id="bib12"><element-citation publication-type="journal" id="sref12"><person-group person-group-type="author"><name><surname>Guo</surname>
<given-names>G.</given-names>
</name>
<name><surname>Huss</surname>
<given-names>M.</given-names>
</name>
<name><surname>Tong</surname>
<given-names>G.Q.</given-names>
</name>
<name><surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name><surname>Li Sun</surname>
<given-names>L.</given-names>
</name>
<name><surname>Clarke</surname>
<given-names>N.D.</given-names>
</name>
<name><surname>Robson</surname>
<given-names>P.</given-names>
</name>
</person-group>
<article-title>Resolution of cell fate decisions revealed by single-cell gene expression analysis from zygote to blastocyst</article-title>
<source>Dev. Cell</source>
<volume>18</volume>
<year>2010</year>
<fpage>675</fpage>
<lpage>685</lpage>
<pub-id pub-id-type="pmid">20412781</pub-id>
</element-citation>
</ref>
<ref id="bib13"><element-citation publication-type="journal" id="sref13"><person-group person-group-type="author"><name><surname>Hillman</surname>
<given-names>N.</given-names>
</name>
<name><surname>Sherman</surname>
<given-names>M.I.</given-names>
</name>
<name><surname>Graham</surname>
<given-names>C.</given-names>
</name>
</person-group>
<article-title>The effect of spatial arrangement on cell determination during mouse development</article-title>
<source>J. Embryol. Exp. Morphol.</source>
<volume>28</volume>
<year>1972</year>
<fpage>263</fpage>
<lpage>278</lpage>
<pub-id pub-id-type="pmid">4674567</pub-id>
</element-citation>
</ref>
<ref id="bib14"><element-citation publication-type="journal" id="sref14"><person-group person-group-type="author"><name><surname>Jedrusik</surname>
<given-names>A.</given-names>
</name>
<name><surname>Parfitt</surname>
<given-names>D.E.</given-names>
</name>
<name><surname>Guo</surname>
<given-names>G.</given-names>
</name>
<name><surname>Skamagki</surname>
<given-names>M.</given-names>
</name>
<name><surname>Grabarek</surname>
<given-names>J.B.</given-names>
</name>
<name><surname>Johnson</surname>
<given-names>M.H.</given-names>
</name>
<name><surname>Robson</surname>
<given-names>P.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Role of Cdx2 and cell polarity in cell allocation and specification of trophectoderm and inner cell mass in the mouse embryo</article-title>
<source>Genes Dev.</source>
<volume>22</volume>
<year>2008</year>
<fpage>2692</fpage>
<lpage>2706</lpage>
<pub-id pub-id-type="pmid">18832072</pub-id>
</element-citation>
</ref>
<ref id="bib15"><element-citation publication-type="journal" id="sref15"><person-group person-group-type="author"><name><surname>Johnson</surname>
<given-names>M.H.</given-names>
</name>
<name><surname>Ziomek</surname>
<given-names>C.A.</given-names>
</name>
</person-group>
<article-title>The foundation of two distinct cell lineages within the mouse morula</article-title>
<source>Cell</source>
<volume>24</volume>
<year>1981</year>
<fpage>71</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="pmid">7237545</pub-id>
</element-citation>
</ref>
<ref id="bib16"><element-citation publication-type="journal" id="sref16"><person-group person-group-type="author"><name><surname>Johnson</surname>
<given-names>W.E.</given-names>
</name>
<name><surname>Li</surname>
<given-names>C.</given-names>
</name>
<name><surname>Rabinovic</surname>
<given-names>A.</given-names>
</name>
</person-group>
<article-title>Adjusting batch effects in microarray expression data using empirical Bayes methods</article-title>
<source>Biostatistics</source>
<volume>8</volume>
<year>2007</year>
<fpage>118</fpage>
<lpage>127</lpage>
<pub-id pub-id-type="pmid">16632515</pub-id>
</element-citation>
</ref>
<ref id="bib17"><element-citation publication-type="journal" id="sref17"><person-group person-group-type="author"><name><surname>Kiso</surname>
<given-names>M.</given-names>
</name>
<name><surname>Tanaka</surname>
<given-names>S.</given-names>
</name>
<name><surname>Saba</surname>
<given-names>R.</given-names>
</name>
<name><surname>Matsuda</surname>
<given-names>S.</given-names>
</name>
<name><surname>Shimizu</surname>
<given-names>A.</given-names>
</name>
<name><surname>Ohyama</surname>
<given-names>M.</given-names>
</name>
<name><surname>Okano</surname>
<given-names>H.J.</given-names>
</name>
<name><surname>Shiroishi</surname>
<given-names>T.</given-names>
</name>
<name><surname>Okano</surname>
<given-names>H.</given-names>
</name>
<name><surname>Saga</surname>
<given-names>Y.</given-names>
</name>
</person-group>
<article-title>The disruption of Sox21-mediated hair shaft cuticle differentiation causes cyclic alopecia in mice</article-title>
<source>Proc. Natl. Acad. Sci. USA</source>
<volume>106</volume>
<year>2009</year>
<fpage>9292</fpage>
<lpage>9297</lpage>
<pub-id pub-id-type="pmid">19470461</pub-id>
</element-citation>
</ref>
<ref id="bib18"><element-citation publication-type="journal" id="sref18"><person-group person-group-type="author"><name><surname>Kuzmichev</surname>
<given-names>A.N.</given-names>
</name>
<name><surname>Kim</surname>
<given-names>S.K.</given-names>
</name>
<name><surname>D’Alessio</surname>
<given-names>A.C.</given-names>
</name>
<name><surname>Chenoweth</surname>
<given-names>J.G.</given-names>
</name>
<name><surname>Wittko</surname>
<given-names>I.M.</given-names>
</name>
<name><surname>Campanati</surname>
<given-names>L.</given-names>
</name>
<name><surname>McKay</surname>
<given-names>R.D.</given-names>
</name>
</person-group>
<article-title>Sox2 acts through Sox21 to regulate transcription in pluripotent and differentiated cells</article-title>
<source>Curr. Biol.</source>
<volume>22</volume>
<year>2012</year>
<fpage>1705</fpage>
<lpage>1710</lpage>
<pub-id pub-id-type="pmid">22902753</pub-id>
</element-citation>
</ref>
<ref id="bib19"><element-citation publication-type="journal" id="sref19"><person-group person-group-type="author"><name><surname>Langfelder</surname>
<given-names>P.</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>B.</given-names>
</name>
<name><surname>Horvath</surname>
<given-names>S.</given-names>
</name>
</person-group>
<article-title>Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R</article-title>
<source>Bioinformatics</source>
<volume>24</volume>
<year>2008</year>
<fpage>719</fpage>
<lpage>720</lpage>
<pub-id pub-id-type="pmid">18024473</pub-id>
</element-citation>
</ref>
<ref id="bib20"><element-citation publication-type="journal" id="sref20"><person-group person-group-type="author"><name><surname>Loh</surname>
<given-names>Y.H.</given-names>
</name>
<name><surname>Wu</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Chew</surname>
<given-names>J.L.</given-names>
</name>
<name><surname>Vega</surname>
<given-names>V.B.</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name><surname>Chen</surname>
<given-names>X.</given-names>
</name>
<name><surname>Bourque</surname>
<given-names>G.</given-names>
</name>
<name><surname>George</surname>
<given-names>J.</given-names>
</name>
<name><surname>Leong</surname>
<given-names>B.</given-names>
</name>
<name><surname>Liu</surname>
<given-names>J.</given-names>
</name>
</person-group>
<article-title>The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells</article-title>
<source>Nat. Genet.</source>
<volume>38</volume>
<year>2006</year>
<fpage>431</fpage>
<lpage>440</lpage>
<pub-id pub-id-type="pmid">16518401</pub-id>
</element-citation>
</ref>
<ref id="bib21"><element-citation publication-type="journal" id="sref21"><person-group person-group-type="author"><name><surname>Love</surname>
<given-names>M.I.</given-names>
</name>
<name><surname>Huber</surname>
<given-names>W.</given-names>
</name>
<name><surname>Anders</surname>
<given-names>S.</given-names>
</name>
</person-group>
<article-title>Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2</article-title>
<source>Genome Biol.</source>
<volume>15</volume>
<year>2014</year>
<fpage>550</fpage>
<pub-id pub-id-type="pmid">25516281</pub-id>
</element-citation>
</ref>
<ref id="bib22"><element-citation publication-type="journal" id="sref22"><person-group person-group-type="author"><name><surname>Mallanna</surname>
<given-names>S.K.</given-names>
</name>
<name><surname>Ormsbee</surname>
<given-names>B.D.</given-names>
</name>
<name><surname>Iacovino</surname>
<given-names>M.</given-names>
</name>
<name><surname>Gilmore</surname>
<given-names>J.M.</given-names>
</name>
<name><surname>Cox</surname>
<given-names>J.L.</given-names>
</name>
<name><surname>Kyba</surname>
<given-names>M.</given-names>
</name>
<name><surname>Washburn</surname>
<given-names>M.P.</given-names>
</name>
<name><surname>Rizzino</surname>
<given-names>A.</given-names>
</name>
</person-group>
<article-title>Proteomic analysis of Sox2-associated proteins during early stages of mouse embryonic stem cell differentiation identifies Sox21 as a novel regulator of stem cell fate</article-title>
<source>Stem Cells</source>
<volume>28</volume>
<year>2010</year>
<fpage>1715</fpage>
<lpage>1727</lpage>
<pub-id pub-id-type="pmid">20687156</pub-id>
</element-citation>
</ref>
<ref id="bib23"><element-citation publication-type="journal" id="sref23"><person-group person-group-type="author"><name><surname>Matoba</surname>
<given-names>R.</given-names>
</name>
<name><surname>Niwa</surname>
<given-names>H.</given-names>
</name>
<name><surname>Masui</surname>
<given-names>S.</given-names>
</name>
<name><surname>Ohtsuka</surname>
<given-names>S.</given-names>
</name>
<name><surname>Carter</surname>
<given-names>M.G.</given-names>
</name>
<name><surname>Sharov</surname>
<given-names>A.A.</given-names>
</name>
<name><surname>Ko</surname>
<given-names>M.S.</given-names>
</name>
</person-group>
<article-title>Dissecting Oct3/4-regulated gene networks in embryonic stem cells by expression profiling</article-title>
<source>PLoS ONE</source>
<volume>1</volume>
<year>2006</year>
<fpage>e26</fpage>
<pub-id pub-id-type="pmid">17183653</pub-id>
</element-citation>
</ref>
<ref id="bib24"><element-citation publication-type="journal" id="sref24"><person-group person-group-type="author"><name><surname>McDole</surname>
<given-names>K.</given-names>
</name>
<name><surname>Xiong</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Iglesias</surname>
<given-names>P.A.</given-names>
</name>
<name><surname>Zheng</surname>
<given-names>Y.</given-names>
</name>
</person-group>
<article-title>Lineage mapping the pre-implantation mouse embryo by two-photon microscopy, new insights into the segregation of cell fates</article-title>
<source>Dev. Biol.</source>
<volume>355</volume>
<year>2011</year>
<fpage>239</fpage>
<lpage>249</lpage>
<pub-id pub-id-type="pmid">21539832</pub-id>
</element-citation>
</ref>
<ref id="bib25"><element-citation publication-type="journal" id="sref25"><person-group person-group-type="author"><name><surname>Morris</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Teo</surname>
<given-names>R.T.</given-names>
</name>
<name><surname>Li</surname>
<given-names>H.</given-names>
</name>
<name><surname>Robson</surname>
<given-names>P.</given-names>
</name>
<name><surname>Glover</surname>
<given-names>D.M.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Origin and formation of the first two distinct cell types of the inner cell mass in the mouse embryo</article-title>
<source>Proc. Natl. Acad. Sci. USA</source>
<volume>107</volume>
<year>2010</year>
<fpage>6364</fpage>
<lpage>6369</lpage>
<pub-id pub-id-type="pmid">20308546</pub-id>
</element-citation>
</ref>
<ref id="bib26"><element-citation publication-type="journal" id="sref26"><person-group person-group-type="author"><name><surname>Morris</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Developmental plasticity is bound by pluripotency and the Fgf and Wnt signaling pathways</article-title>
<source>Cell Rep.</source>
<volume>2</volume>
<year>2012</year>
<fpage>756</fpage>
<lpage>765</lpage>
<pub-id pub-id-type="pmid">23041313</pub-id>
</element-citation>
</ref>
<ref id="bib27"><element-citation publication-type="journal" id="sref27"><person-group person-group-type="author"><name><surname>Parfitt</surname>
<given-names>D.E.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Epigenetic modification affecting expression of cell polarity and cell fate genes to regulate lineage specification in the early mouse embryo</article-title>
<source>Mol. Biol. Cell</source>
<volume>21</volume>
<year>2010</year>
<fpage>2649</fpage>
<lpage>2660</lpage>
<pub-id pub-id-type="pmid">20554762</pub-id>
</element-citation>
</ref>
<ref id="bib28"><element-citation publication-type="journal" id="sref28"><person-group person-group-type="author"><name><surname>Picelli</surname>
<given-names>S.</given-names>
</name>
<name><surname>Faridani</surname>
<given-names>O.R.</given-names>
</name>
<name><surname>Björklund</surname>
<given-names>A.K.</given-names>
</name>
<name><surname>Winberg</surname>
<given-names>G.</given-names>
</name>
<name><surname>Sagasser</surname>
<given-names>S.</given-names>
</name>
<name><surname>Sandberg</surname>
<given-names>R.</given-names>
</name>
</person-group>
<article-title>Full-length RNA-seq from single cells using Smart-seq2</article-title>
<source>Nat. Protoc.</source>
<volume>9</volume>
<year>2014</year>
<fpage>171</fpage>
<lpage>181</lpage>
<pub-id pub-id-type="pmid">24385147</pub-id>
</element-citation>
</ref>
<ref id="bib29"><element-citation publication-type="journal" id="sref29"><person-group person-group-type="author"><name><surname>Piotrowska</surname>
<given-names>K.</given-names>
</name>
<name><surname>Wianny</surname>
<given-names>F.</given-names>
</name>
<name><surname>Pedersen</surname>
<given-names>R.A.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Blastomeres arising from the first cleavage division have distinguishable fates in normal mouse development</article-title>
<source>Development</source>
<volume>128</volume>
<year>2001</year>
<fpage>3739</fpage>
<lpage>3748</lpage>
<pub-id pub-id-type="pmid">11585800</pub-id>
</element-citation>
</ref>
<ref id="bib30"><element-citation publication-type="journal" id="sref30"><person-group person-group-type="author"><name><surname>Piotrowska-Nitsche</surname>
<given-names>K.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Spatial arrangement of individual 4-cell stage blastomeres and the order in which they are generated correlate with blastocyst pattern in the mouse embryo</article-title>
<source>Mech. Dev.</source>
<volume>122</volume>
<year>2005</year>
<fpage>487</fpage>
<lpage>500</lpage>
<pub-id pub-id-type="pmid">15804563</pub-id>
</element-citation>
</ref>
<ref id="bib31"><element-citation publication-type="journal" id="sref31"><person-group person-group-type="author"><name><surname>Piotrowska-Nitsche</surname>
<given-names>K.</given-names>
</name>
<name><surname>Perea-Gomez</surname>
<given-names>A.</given-names>
</name>
<name><surname>Haraguchi</surname>
<given-names>S.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Four-cell stage mouse blastomeres have different developmental properties</article-title>
<source>Development</source>
<volume>132</volume>
<year>2005</year>
<fpage>479</fpage>
<lpage>490</lpage>
<pub-id pub-id-type="pmid">15634695</pub-id>
</element-citation>
</ref>
<ref id="bib32"><element-citation publication-type="journal" id="sref32"><person-group person-group-type="author"><name><surname>Plachta</surname>
<given-names>N.</given-names>
</name>
<name><surname>Bollenbach</surname>
<given-names>T.</given-names>
</name>
<name><surname>Pease</surname>
<given-names>S.</given-names>
</name>
<name><surname>Fraser</surname>
<given-names>S.E.</given-names>
</name>
<name><surname>Pantazis</surname>
<given-names>P.</given-names>
</name>
</person-group>
<article-title>Oct4 kinetics predict cell lineage patterning in the early mammalian embryo</article-title>
<source>Nat. Cell Biol.</source>
<volume>13</volume>
<year>2011</year>
<fpage>117</fpage>
<lpage>123</lpage>
<pub-id pub-id-type="pmid">21258368</pub-id>
</element-citation>
</ref>
<ref id="bib33"><element-citation publication-type="journal" id="sref33"><person-group person-group-type="author"><name><surname>Samarage</surname>
<given-names>C.R.</given-names>
</name>
<name><surname>White</surname>
<given-names>M.D.</given-names>
</name>
<name><surname>Álvarez</surname>
<given-names>Y.D.</given-names>
</name>
<name><surname>Fierro-González</surname>
<given-names>J.C.</given-names>
</name>
<name><surname>Henon</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Jesudason</surname>
<given-names>E.C.</given-names>
</name>
<name><surname>Bissiere</surname>
<given-names>S.</given-names>
</name>
<name><surname>Fouras</surname>
<given-names>A.</given-names>
</name>
<name><surname>Plachta</surname>
<given-names>N.</given-names>
</name>
</person-group>
<article-title>Cortical tension allocates the first inner cells of the mammalian embryo</article-title>
<source>Dev. Cell</source>
<volume>34</volume>
<year>2015</year>
<fpage>435</fpage>
<lpage>447</lpage>
<pub-id pub-id-type="pmid">26279486</pub-id>
</element-citation>
</ref>
<ref id="bib34"><element-citation publication-type="journal" id="sref34"><person-group person-group-type="author"><name><surname>Sharov</surname>
<given-names>A.A.</given-names>
</name>
<name><surname>Masui</surname>
<given-names>S.</given-names>
</name>
<name><surname>Sharova</surname>
<given-names>L.V.</given-names>
</name>
<name><surname>Piao</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Aiba</surname>
<given-names>K.</given-names>
</name>
<name><surname>Matoba</surname>
<given-names>R.</given-names>
</name>
<name><surname>Xin</surname>
<given-names>L.</given-names>
</name>
<name><surname>Niwa</surname>
<given-names>H.</given-names>
</name>
<name><surname>Ko</surname>
<given-names>M.S.</given-names>
</name>
</person-group>
<article-title>Identification of Pou5f1, Sox2, and Nanog downstream target genes with statistical confidence by applying a novel algorithm to time course microarray and genome-wide chromatin immunoprecipitation data</article-title>
<source>BMC Genomics</source>
<volume>9</volume>
<year>2008</year>
<fpage>269</fpage>
<pub-id pub-id-type="pmid">18522731</pub-id>
</element-citation>
</ref>
<ref id="bib35"><element-citation publication-type="journal" id="sref35"><person-group person-group-type="author"><name><surname>Shi</surname>
<given-names>J.</given-names>
</name>
<name><surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Li</surname>
<given-names>X.</given-names>
</name>
<name><surname>Zheng</surname>
<given-names>X.</given-names>
</name>
<name><surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Qiao</surname>
<given-names>J.</given-names>
</name>
<name><surname>Tang</surname>
<given-names>F.</given-names>
</name>
<name><surname>Tao</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Zhou</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Duan</surname>
<given-names>E.</given-names>
</name>
</person-group>
<article-title>Dynamic transcriptional symmetry-breaking in pre-implantation mammalian embryo development revealed by single-cell RNA-seq</article-title>
<source>Development</source>
<volume>142</volume>
<year>2015</year>
<fpage>3468</fpage>
<lpage>3477</lpage>
<pub-id pub-id-type="pmid">26395495</pub-id>
</element-citation>
</ref>
<ref id="bib36"><element-citation publication-type="journal" id="sref36"><person-group person-group-type="author"><name><surname>Stegle</surname>
<given-names>O.</given-names>
</name>
<name><surname>Teichmann</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Marioni</surname>
<given-names>J.C.</given-names>
</name>
</person-group>
<article-title>Computational and analytical challenges in single-cell transcriptomics</article-title>
<source>Nat. Rev. Genet.</source>
<volume>16</volume>
<year>2015</year>
<fpage>133</fpage>
<lpage>145</lpage>
<pub-id pub-id-type="pmid">25628217</pub-id>
</element-citation>
</ref>
<ref id="bib37"><element-citation publication-type="journal" id="sref37"><person-group person-group-type="author"><name><surname>Strumpf</surname>
<given-names>D.</given-names>
</name>
<name><surname>Mao</surname>
<given-names>C.A.</given-names>
</name>
<name><surname>Yamanaka</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Ralston</surname>
<given-names>A.</given-names>
</name>
<name><surname>Chawengsaksophak</surname>
<given-names>K.</given-names>
</name>
<name><surname>Beck</surname>
<given-names>F.</given-names>
</name>
<name><surname>Rossant</surname>
<given-names>J.</given-names>
</name>
</person-group>
<article-title>Cdx2 is required for correct cell fate specification and differentiation of trophectoderm in the mouse blastocyst</article-title>
<source>Development</source>
<volume>132</volume>
<year>2005</year>
<fpage>2093</fpage>
<lpage>2102</lpage>
<pub-id pub-id-type="pmid">15788452</pub-id>
</element-citation>
</ref>
<ref id="bib38"><element-citation publication-type="journal" id="sref38"><person-group person-group-type="author"><name><surname>Tabansky</surname>
<given-names>I.</given-names>
</name>
<name><surname>Lenarcic</surname>
<given-names>A.</given-names>
</name>
<name><surname>Draft</surname>
<given-names>R.W.</given-names>
</name>
<name><surname>Loulier</surname>
<given-names>K.</given-names>
</name>
<name><surname>Keskin</surname>
<given-names>D.B.</given-names>
</name>
<name><surname>Rosains</surname>
<given-names>J.</given-names>
</name>
<name><surname>Rivera-Feliciano</surname>
<given-names>J.</given-names>
</name>
<name><surname>Lichtman</surname>
<given-names>J.W.</given-names>
</name>
<name><surname>Livet</surname>
<given-names>J.</given-names>
</name>
<name><surname>Stern</surname>
<given-names>J.N.</given-names>
</name>
</person-group>
<article-title>Developmental bias in cleavage-stage mouse blastomeres</article-title>
<source>Curr. Biol.</source>
<volume>23</volume>
<year>2013</year>
<fpage>21</fpage>
<lpage>31</lpage>
<pub-id pub-id-type="pmid">23177476</pub-id>
</element-citation>
</ref>
<ref id="bib39"><element-citation publication-type="journal" id="sref39"><person-group person-group-type="author"><name><surname>Takaoka</surname>
<given-names>K.</given-names>
</name>
<name><surname>Hamada</surname>
<given-names>H.</given-names>
</name>
</person-group>
<article-title>Cell fate decisions and axis determination in the early mouse embryo</article-title>
<source>Development</source>
<volume>139</volume>
<year>2012</year>
<fpage>3</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="pmid">22147950</pub-id>
</element-citation>
</ref>
<ref id="bib40"><element-citation publication-type="journal" id="sref40"><person-group person-group-type="author"><name><surname>Tang</surname>
<given-names>F.</given-names>
</name>
<name><surname>Barbacioru</surname>
<given-names>C.</given-names>
</name>
<name><surname>Nordman</surname>
<given-names>E.</given-names>
</name>
<name><surname>Bao</surname>
<given-names>S.</given-names>
</name>
<name><surname>Lee</surname>
<given-names>C.</given-names>
</name>
<name><surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name><surname>Tuch</surname>
<given-names>B.B.</given-names>
</name>
<name><surname>Heard</surname>
<given-names>E.</given-names>
</name>
<name><surname>Lao</surname>
<given-names>K.</given-names>
</name>
<name><surname>Surani</surname>
<given-names>M.A.</given-names>
</name>
</person-group>
<article-title>Deterministic and stochastic allele specific gene expression in single mouse blastomeres</article-title>
<source>PLoS ONE</source>
<volume>6</volume>
<year>2011</year>
<fpage>e21208</fpage>
<pub-id pub-id-type="pmid">21731673</pub-id>
</element-citation>
</ref>
<ref id="bib41"><element-citation publication-type="journal" id="sref41"><person-group person-group-type="author"><name><surname>Tarkowski</surname>
<given-names>A.K.</given-names>
</name>
</person-group>
<article-title>Experiments on the development of isolated blastomers of mouse eggs</article-title>
<source>Nature</source>
<volume>184</volume>
<year>1959</year>
<fpage>1286</fpage>
<lpage>1287</lpage>
<pub-id pub-id-type="pmid">13836947</pub-id>
</element-citation>
</ref>
<ref id="bib42"><element-citation publication-type="journal" id="sref42"><person-group person-group-type="author"><name><surname>Torres-Padilla</surname>
<given-names>M.E.</given-names>
</name>
<name><surname>Parfitt</surname>
<given-names>D.E.</given-names>
</name>
<name><surname>Kouzarides</surname>
<given-names>T.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>Histone arginine methylation regulates pluripotency in the early mouse embryo</article-title>
<source>Nature</source>
<volume>445</volume>
<year>2007</year>
<fpage>214</fpage>
<lpage>218</lpage>
<pub-id pub-id-type="pmid">17215844</pub-id>
</element-citation>
</ref>
<ref id="bib43"><element-citation publication-type="journal" id="sref43"><person-group person-group-type="author"><name><surname>Watanabe</surname>
<given-names>T.</given-names>
</name>
<name><surname>Biggins</surname>
<given-names>J.S.</given-names>
</name>
<name><surname>Tannan</surname>
<given-names>N.B.</given-names>
</name>
<name><surname>Srinivas</surname>
<given-names>S.</given-names>
</name>
</person-group>
<article-title>Limited predictive value of blastomere angle of division in trophectoderm and inner cell mass specification</article-title>
<source>Development</source>
<volume>141</volume>
<year>2014</year>
<fpage>2279</fpage>
<lpage>2288</lpage>
<pub-id pub-id-type="pmid">24866117</pub-id>
</element-citation>
</ref>
<ref id="bib44"><element-citation publication-type="journal" id="sref44"><person-group person-group-type="author"><name><surname>Wu</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Bruce</surname>
<given-names>A.W.</given-names>
</name>
<name><surname>Jedrusik</surname>
<given-names>A.</given-names>
</name>
<name><surname>Ellis</surname>
<given-names>P.D.</given-names>
</name>
<name><surname>Andrews</surname>
<given-names>R.M.</given-names>
</name>
<name><surname>Langford</surname>
<given-names>C.F.</given-names>
</name>
<name><surname>Glover</surname>
<given-names>D.M.</given-names>
</name>
<name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
</person-group>
<article-title>CARM1 is required in embryonic stem cells to maintain pluripotency and resist differentiation</article-title>
<source>Stem Cells</source>
<volume>27</volume>
<year>2009</year>
<fpage>2637</fpage>
<lpage>2645</lpage>
<pub-id pub-id-type="pmid">19544422</pub-id>
</element-citation>
</ref>
<ref id="bib45"><element-citation publication-type="journal" id="sref45"><person-group person-group-type="author"><name><surname>Xue</surname>
<given-names>Z.</given-names>
</name>
<name><surname>Huang</surname>
<given-names>K.</given-names>
</name>
<name><surname>Cai</surname>
<given-names>C.</given-names>
</name>
<name><surname>Cai</surname>
<given-names>L.</given-names>
</name>
<name><surname>Jiang</surname>
<given-names>C.Y.</given-names>
</name>
<name><surname>Feng</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Liu</surname>
<given-names>Z.</given-names>
</name>
<name><surname>Zeng</surname>
<given-names>Q.</given-names>
</name>
<name><surname>Cheng</surname>
<given-names>L.</given-names>
</name>
<name><surname>Sun</surname>
<given-names>Y.E.</given-names>
</name>
</person-group>
<article-title>Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing</article-title>
<source>Nature</source>
<volume>500</volume>
<year>2013</year>
<fpage>593</fpage>
<lpage>597</lpage>
<pub-id pub-id-type="pmid">23892778</pub-id>
</element-citation>
</ref>
<ref id="bib46"><element-citation publication-type="journal" id="sref46"><person-group person-group-type="author"><name><surname>Yamanaka</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Lanner</surname>
<given-names>F.</given-names>
</name>
<name><surname>Rossant</surname>
<given-names>J.</given-names>
</name>
</person-group>
<article-title>FGF signal-dependent segregation of primitive endoderm and epiblast in the mouse blastocyst</article-title>
<source>Development</source>
<volume>137</volume>
<year>2010</year>
<fpage>715</fpage>
<lpage>724</lpage>
<pub-id pub-id-type="pmid">20147376</pub-id>
</element-citation>
</ref>
<ref id="bib47"><element-citation publication-type="journal" id="sref47"><person-group person-group-type="author"><name><surname>Zernicka-Goetz</surname>
<given-names>M.</given-names>
</name>
<name><surname>Morris</surname>
<given-names>S.A.</given-names>
</name>
<name><surname>Bruce</surname>
<given-names>A.W.</given-names>
</name>
</person-group>
<article-title>Making a firm decision: multifaceted regulation of cell fate in the early mouse embryo</article-title>
<source>Nat. Rev. Genet.</source>
<volume>10</volume>
<year>2009</year>
<fpage>467</fpage>
<lpage>477</lpage>
<pub-id pub-id-type="pmid">19536196</pub-id>
</element-citation>
</ref>
</ref-list>
<sec id="app1"><title>Accession Numbers</title>
<p>The accession number for the sequencing data reported in this paper is ArrayExpress: <ext-link ext-link-type="uri" xlink:href="array-express:E-MTAB-3321" id="intref0010">E-MTAB-3321</ext-link>
.</p>
</sec>
<sec id="app3" sec-type="supplementary-material"><title>Supplemental Information</title>
<p><supplementary-material content-type="local-data" id="mmc1"><caption><title>Document S1. Tables S1 and S2</title>
</caption>
<media xlink:href="mmc1.pdf"></media>
</supplementary-material>
<supplementary-material content-type="local-data" id="mmc2"><caption><title>Data S1. Highly Variable Genes at the 2-, 4-, and 8-Cell Stages, Related to Figure 1</title>
</caption>
<media xlink:href="mmc2.xlsx"></media>
</supplementary-material>
<supplementary-material content-type="local-data" id="mmc3"><caption><title>Document S2. Article plus Supplemental Information</title>
</caption>
<media xlink:href="mmc3.pdf"></media>
</supplementary-material>
</p>
</sec>
<ack id="ack0010"><title>Acknowledgments</title>
<p>We are grateful to the Wellcome Trust, EMBL, and CRUK for supporting our work and our colleagues: M. Shahbazi, D. Glover, and F. Antonica for feedback and G. Recher for imaging assistance. T.V. also acknowledges KU Leuven SymBioSys (PFV/10/016).</p>
</ack>
<fn-group><fn id="d32e383"><p id="ccnp0005">This is an open access article under the CC BY license (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" id="ccintref0005">http://creativecommons.org/licenses/by/4.0/</ext-link>
).</p>
</fn>
<fn id="app2" fn-type="supplementary-material"><p>Supplemental Information includes six figures, two tables, and one data file and can be found with this article online at <ext-link ext-link-type="doi" xlink:href="10.1016/j.cell.2016.01.047" id="intref0015">http://dx.doi.org/10.1016/j.cell.2016.01.047</ext-link>
.</p>
</fn>
</fn-group>
</back>
<floats-group><fig id="fig1"><label>Figure 1</label>
<caption><p>Cellular Heterogeneity in Overall Gene Expression Patterns</p>
<p>(A) Scheme of development from the fertilized egg (zygote) to the blastocyst stage with three lineages: epiblast (EPI; blue), primitive endoderm (PE, red), and trophectoderm (TE, green). The mechanism(s) initiating these cell-fate decisions remains unknown (reviewed in <xref rid="bib47" ref-type="bibr">Zernicka-Goetz et al., 2009</xref>
).</p>
<p>(B) Principal component analysis of the gene expression patterns of single cells at different stages. The percentage of variance explained by each principal component is indicated in parentheses.</p>
<p>(C) Spearman correlation coefficients of blastomeres within the same embryo (intra-embryonic) and in different embryos (inter-embryonic) at the 2-, 4-, and 8-cell stage. Asterisks indicate statistically significant differences (Welch’s t test).</p>
<p>(D) Venn diagram with the number of genes having a significantly high degree of heterogeneous expression at 2-, 4-, and 8-cell stages.</p>
<p>(E) Principal component analysis of the transcriptomes of cells from embryos at the 4-cell stage with different division patterns: ME, six embryos; EM, six embryos; MM, three embryos; EE, one embryo. M, meridional; E, equatorial cell division.</p>
<p>See also <xref rid="figs1" ref-type="fig">Figures S1</xref>
, <xref rid="figs2" ref-type="fig">S2</xref>
, <xref rid="figs3" ref-type="fig">S3</xref>
, and <xref rid="figs4" ref-type="fig">S4</xref>
, <xref rid="mmc1" ref-type="supplementary-material">Tables S1</xref>
and <xref rid="mmc1" ref-type="supplementary-material">S2</xref>
, and <xref rid="mmc2" ref-type="supplementary-material">Data S1</xref>
.</p>
</caption>
<graphic xlink:href="gr1"></graphic>
</fig>
<fig id="fig2"><label>Figure 2</label>
<caption><p>Heterogeneous Expression of Oct4 and Sox2 Target Genes at the 4-Cell Stage</p>
<p>(A) To find highly variable genes, we used the parameterization CV<sup>2</sup>
= a<sub>1</sub>
/μ + α<sub>0</sub>
to fit the relationship between the square of the coefficient of variation, CV<sup>2</sup>
, and the average expression level μ (green continuous curve; see <xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
for further details). All highly variable genes (with an adjusted p value < 0.1) marked by red circles. Sox21, cyan circle.</p>
<p>(A′) The same information from (A) but with Oct4 and/or Sox2 targets indicated (diamonds mark shared Sox2 and Oct4 targets, squares mark Oct4 targets, and triangles mark Sox2 targets). Gray symbols mark Oct4 and/or Sox2 targets that are not variable and blue symbols mark Oct4 and/or Sox2 targets that are highly variable at the 4-cell stage. Sox21 is marked in gold.</p>
<p>(B) Heatmap showing the Spearman correlation coefficient among highly variable Oct4 and/or Sox2 targets. A dynamic tree cut algorithm identified two clusters of genes showing similar patterns of correlations (colored side bars; <xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
).</p>
<p>See also <xref rid="figs5" ref-type="fig">Figure S5</xref>
.</p>
</caption>
<graphic xlink:href="gr2"></graphic>
</fig>
<fig id="fig3"><label>Figure 3</label>
<caption><p>Sox21 mRNA and Protein Expression Is Heterogeneous and Peaks at the 4-Cell Stage</p>
<p>(A) Relative mRNA expression of Sox21, Nanog, and Esrrb.</p>
<p>(B) Average relative Sox21 mRNA levels in 4-cell embryos.</p>
<p>(C) Relative Sox21 mRNA expression in all individual 4-cell embryos.</p>
<p>(D) Immunofluorescence of Sox21 in 4-cell (n = 25) and 8-cell (n = 23) embryos.</p>
<p>Fluorescence quantified and normalized to the nucleus with the strongest staining per individual 4-cell or 8-cell embryo. Arrowheads indicate highest expressing cell. Asterisks indicated lowest expressing cells. Error bars represent SEM. Scale bars, 10 μm.</p>
</caption>
<graphic xlink:href="gr3"></graphic>
</fig>
<fig id="fig4"><label>Figure 4</label>
<caption><p>Decreasing Sox21 Expression Leads Cells to an Extra-Embryonic Fate</p>
<p>(A) Verification of <italic>Sox21</italic>
siRNA efficiency at the 4-cell stage (Control siRNA n = 20, <italic>Sox21</italic>
siRNA n = 21). The nuclear immunofluorescence signal is lost in embryos injected with <italic>Sox21</italic>
siRNA.</p>
<p>(B) Scheme of clonal <italic>Sox21</italic>
siRNA. One blastomere of 2-cell stage embryos injected with <italic>Sox21</italic>
siRNA or control siRNA, and GFP mRNA. Embryos cultured to the late blastocyst stage and the contribution of the injected cells’ progeny to each lineage analyzed.</p>
<p>(C) Confocal images of control (n = 38) and <italic>Sox21</italic>
(n = 75) siRNA embryos. Sox17 (primitive endoderm, PE) and Cdx2 (TE) used as lineage markers. Dotted lines mark the ICM.</p>
<p>(D) Contribution of <italic>Sox21</italic>
siRNA cells to TE, PE, and EPI, relative to control siRNA cells.</p>
<p>(E) Time-lapse study following <italic>Sox21</italic>
siRNA (n = 16 embryos, 256 cells) or control siRNA (n = 18 embryos, 288 cells). Gap43-GFP expression indicates blastomeres injected with either control or <italic>Sox21</italic>
siRNA. Fluorescence images are overlaid with cell-tracking spheres.</p>
<p>(F) Example slice through embryos showing inside and outside cell fate of either control or <italic>Sox21</italic>
siRNA blastomeres.</p>
<p>(G) Number of inside and outside cells contributed from injected blastomeres at the 32-cell stage.</p>
<p>(H) Number of asymmetric divisions injected blastomeres underwent from the 8- to 32-cell stage.</p>
<p>(I and J) Lineage trees from two representative embryos. All injected cells were traced to the 32-cell stage. A cell was defined as occupying an inside position when enclosed from the outside environment by neighboring cells. Inside or outside cell fate indicated.</p>
<p>Error bars represent SEM. Wilcoxon rank-sum test was used to test significance <sup>∗∗</sup>
p < 0.01, <sup>∗∗∗</sup>
p < 0.001. Scale bars represent 10 μm. See also <xref rid="figs6" ref-type="fig">Figure S6</xref>
.</p>
</caption>
<graphic xlink:href="gr4"></graphic>
</fig>
<fig id="fig5"><label>Figure 5</label>
<caption><p>Sox21 Depletion Prematurely Upregulates Cdx2 Expression</p>
<p>(A) Scheme of <italic>Sox21</italic>
siRNA experiment. Zygotes injected with <italic>Sox21</italic>
siRNA, or control siRNA, and isolated at the 8-cell stage for immunostaining or qRT-PCR.</p>
<p>(B) qRT-PCR of embryos injected with either control siRNA (n = 75 embryos, three biological replicates) or <italic>Sox21</italic>
siRNA (n = 85 embryos, three biological replicates) comparing mRNA expression of Cdx2 and Sox21 at the late 8-cell stage. Student’s t test was used to test significance <sup>∗</sup>
p < 0.05, <sup>∗∗∗</sup>
p < 0.001.</p>
<p>(C) Confocal images of Cdx2 and Histone H3 expression in control (n = 12) and <italic>Sox21</italic>
(n = 12) siRNA embryos at the early 8-cell stage.</p>
<p>(D) Quantification of the number of Cdx2 positive cells from (C). Student’s t test was used to test significance <sup>∗∗∗</sup>
p < 0.001.</p>
<p>(E) Scheme of <italic>Sox21</italic>
siRNA experiment. One blastomere of 2-cell stage embryos was injected with <italic>Sox21</italic>
siRNA or control siRNA and Gap43-RFP mRNA and the 8-cell embryos isolated.</p>
<p>(F) Confocal images of Cdx2 and Histone H3 expression in late 8-cell stage embryos after injection of one blastomere at the 2-cell stage with either control siRNA (n = 11) or <italic>Sox21</italic>
siRNA (n = 13). Gap43-RFP expression on the membrane identifies injected cells.</p>
<p>(G) Quantification of relative fluorescence intensity of Cdx2 staining from (F).</p>
<p>Wilcoxon rank-sum test was used to test significance <sup>∗∗∗</sup>
p < 0.001. Scale bars, 10 μm. Error bars represent SEM.</p>
</caption>
<graphic xlink:href="gr5"></graphic>
</fig>
<fig id="fig6"><label>Figure 6</label>
<caption><p>Sox21 Expression Is Regulated by CARM1 Activity</p>
<p>(A) Verification of CARM1 inhibition by staining for the presence of H3R26me (control n = 21, CARM1 inhibition n = 15). Scale bars, 20 μm.</p>
<p>(B) Confocal images of Sox21 expression in 4-cell embryos after being treated from the 2-cell stage with either DMSO control (n = 27) or CARM1 inhibitor (n = 24). Scale bar, 50 μm.</p>
<p>(C) Confocal images of control (n = 12) and <italic>Carm1</italic>
siRNA (n = 8) embryos. Sox17 (PE), Cdx2 (TE), and Nanog (EPI) used as lineage markers. Scale bars, 10 μm.</p>
<p>(D) Overall number of TE, PE, and EPI cells present in <italic>Carm1</italic>
siRNA-injected embryos, relative to control siRNA cells. Student’s t test was used to test significance <sup>∗∗∗</sup>
p < 0.001.</p>
<p>(E) qRT-PCR of embryos injected with either control mRNA (n = 66 embryos, three biological replicates) or CARM1 mRNA (n = 75 embryos, three biological replicates).</p>
<p>Embryos injected at the zygote stage and isolated for qRT-PCR at the late 8-cell stage Student’s t test was used to test significance <sup>∗</sup>
p < 0.05, <sup>∗∗</sup>
p < 0.01. Error bars represent SEM.</p>
</caption>
<graphic xlink:href="gr6"></graphic>
</fig>
<fig id="fig7"><label>Figure 7</label>
<caption><p>Model of How Cellular Heterogeneities at the 4-Cell Stage Regulate Cell Fate</p>
<p>(A–C) In 4-cell embryos, CARM1, which methylates histone H3R26, is differentially expressed (<xref rid="bib42" ref-type="bibr">Torres-Padilla et al., 2007</xref>
). We hypothesize that higher levels of histone H3R26me facilitate the binding to DNA of pluripotency regulators such as Oct4 and Sox2, resulting in increased transcription of pluripotency-related target genes, such as <italic>Sox21</italic>
, <italic>Nanog</italic>
, and <italic>Esrrb</italic>
, biasing these cells to contribute to the pluripotent lineage. Conversely, in cells with lower levels of histone H3R26me, pluripotency regulators are only able to bind to DNA for shorter periods of time, and therefore their target genes are not as highly expressed. These cells have lower levels of pluripotency and are thus more likely to initiate expression of differentiation genes, such as <italic>Cdx2</italic>
, and initiate development into the extra-embryonic TE.</p>
</caption>
<graphic xlink:href="gr7"></graphic>
</fig>
<fig id="figs1"><label>Figure S1</label>
<caption><p>Quality Control Analyses of Single-Cell Transcriptomes, Related to <xref rid="fig1" ref-type="fig">Figure 1</xref>
</p>
<p>(A–C) The fraction of mapped reads (panel A), the number of genes with more than 10 Reads per Million (panel B) and the fraction of mapped reads allocated to mitochondrial genes were plotted as a function of the total number of reads for all samples. Different colors mark different batches.</p>
<p>(D) Principal Component Analysis of the three metrics plotted in panels A-C. The percentage of variance explained by each principal component is indicated in parentheses. The black arrow in all four panels marks the outlier (sample name “32cell_F”) that was removed from all downstream analysis.</p>
</caption>
<graphic xlink:href="figs1"></graphic>
</fig>
<fig id="figs2"><label>Figure S2</label>
<caption><p>Principal Component Analysis of ERCC Spike-ins from the Samples Shown in <xref rid="fig1" ref-type="fig">Figure 1</xref>
B, Related to <xref rid="fig1" ref-type="fig">Figure 1</xref>
</p>
<p>Principal Component Analysis of the log-transformed counts of ERCC spike-ins added to each cells. The numbers in parentheses indicate the percentage of total variance explained by each principal component.</p>
</caption>
<graphic xlink:href="figs2"></graphic>
</fig>
<fig id="figs3"><label>Figure S3</label>
<caption><p>Highly Variable Genes in Embryos at 2- and 8-Cell Stage, Related to <xref rid="fig1" ref-type="fig">Figure 1</xref>
</p>
<p>Genes displaying a high variability in embryos at 2- (panel A) and 8-cell (panel B) stage were identified, as described in <xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
. All highly variable genes (with an adjusted p-value < 0.1) are marked by red circles or cyan diamonds. Diamonds mark Oct4/Sox2 target genes.</p>
</caption>
<graphic xlink:href="figs3"></graphic>
</fig>
<fig id="figs4"><label>Figure S4</label>
<caption><p>Gene Expression Differences of 4-Cell Stage Embryos with Different Division Patterns, Related to <xref rid="fig1" ref-type="fig">Figure 1</xref>
</p>
<p>(A) The four different 4-cell stage embryos defined by division pattern.</p>
<p>(B) Live-imaging of embryos to score 2-to 4-cell division pattern. C) 20 genes with the most positive and negative loadings on the first 2 principal components shown in <xref rid="fig1" ref-type="fig">Figure 1</xref>
C are plotted, with the arrow indicating the correlation with each principal component.</p>
</caption>
<graphic xlink:href="figs4"></graphic>
</fig>
<fig id="figs5"><label>Figure S5</label>
<caption><p>Heatmap Showing the Spearman Correlation Coefficient among Highly Variable Oct4 and/or Sox2 Targets, Related to <xref rid="fig2" ref-type="fig">Figure 2</xref>
</p>
<p>A dynamic tree cut algorithm identified two clusters of genes showing similar patterns of correlations (colored side bars; <xref rid="sec4" ref-type="sec">Experimental Procedures</xref>
).</p>
</caption>
<graphic xlink:href="figs5"></graphic>
</fig>
<fig id="figs6"><label>Figure S6</label>
<caption><p>Sox21 Depletion with Three Individual siRNAs Biases Cells toward an Extra-Embryonic Cell Fate, Related to <xref rid="fig4" ref-type="fig">Figure 4</xref>
</p>
<p>(A) Confocal images of control (n = 24) and Sox21 siRNA 1 (n = 12), Sox21 siRNA 2 (n = 14) and Sox21 siRNA 3 (n = 18) siRNA embryos. Sox17 (PE) and Cdx2 (TE) were used as cell lineage markers. Dotted lines mark the ICM.</p>
<p>(B) Contribution of Sox21 siRNA 1,2 and 3 injected cells to TE, PE and EPI, relative to control siRNA cells.</p>
</caption>
<graphic xlink:href="figs6"></graphic>
</fig>
</floats-group>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Belgique/explor/OpenAccessBelV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000352 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000352 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Belgique |area= OpenAccessBelV2 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.25. | ![]() |