FAS Center for Systems Biology, Harvard University, Cambridge, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
Sandeep Choubey
FAS Center for Systems Biology, Harvard University, Cambridge, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
FAS Center for Systems Biology, Harvard University, Cambridge, United States; Biophysics Program, Harvard University, Cambridge, United States
Ling-Nan Zou
FAS Center for Systems Biology, Harvard University, Cambridge, United States
Adele Doyle
FAS Center for Systems Biology, Harvard University, Cambridge, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
Vilas Menon
Allen Institute for Brain Science, Seattle, United States
Ethan B Loew
FAS Center for Systems Biology, Harvard University, Cambridge, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
Anne-Rachel Krostag
Allen Institute for Brain Science, Seattle, United States
Refugio A Martinez
Allen Institute for Brain Science, Seattle, United States
Linda Madisen
Allen Institute for Brain Science, Seattle, United States
Boaz P Levi
Allen Institute for Brain Science, Seattle, United States
Sharad Ramanathan
FAS Center for Systems Biology, Harvard University, Cambridge, United States; Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States; Allen Institute for Brain Science, Seattle, United States; School of Engineering and Applied Sciences, Harvard University, Cambridge, United States; Harvard Stem Cell Institute, Harvard University, Cambridge, United States
The complexity of gene regulatory networks that lead multipotent cells to acquire different cell fates makes a quantitative understanding of differentiation challenging. Using a statistical framework to analyze single-cell transcriptomics data, we infer the gene expression dynamics of early mouse embryonic stem (mES) cell differentiation, uncovering discrete transitions across nine cell states. We validate the predicted transitions across discrete states using flow cytometry. Moreover, using live-cell microscopy, we show that individual cells undergo abrupt transitions from a naïve to primed pluripotent state. Using the inferred discrete cell states to build a probabilistic model for the underlying gene regulatory network, we further predict and experimentally verify that these states have unique response to perturbations, thus defining them functionally. Our study provides a framework to infer the dynamics of differentiation from single cell transcriptomics data and to build predictive models of the gene regulatory networks that drive the sequence of cell fate decisions during development.