Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution

Small cells - big issues: biological implications and preclinical advancements in small cell lung cancer

Current treatment guidelines refer to small cell lung cancer (SCLC), one of the deadliest human malignancies, as a homogeneous disease. Accordingly, SCLC therapy comprises chemoradiation with or without immunotherapy. Meanwhile, recent studies have made significant advances in subclassifying SCLC based on the elevated expression of the transcription factors ASCL1, NEUROD1, and POU2F3, as well as on certain inflammatory characteristics. The role of the transcription regulator YAP1 in defining a unique SCLC subset remains to be established. Although preclinical analyses have described numerous subtype-specific characteristics and vulnerabilities, the so far non-existing clinical subtype distinction may be a contributor to negative clinical trial outcomes. This comprehensive review aims to provide a framework for the development of novel personalized therapeutic approaches by compiling the most recent discoveries achieved by preclinical SCLC research. We highlight the challenges faced due to limited access to patient material as well as the advances accomplished by implementing state-of-the-art models and methodologies.

Clonally heritable gene expression imparts a layer of diversity within cell types

Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information.

Functional analysis of cell plasticity using single-cell technologies

Metazoan organisms are heterocellular systems composed of hundreds of different cell types, which arise from an isogenic genome through differentiation. Cellular 'plasticity' further enables cells to alter their fate in response to exogenous cues and is involved in a variety of processes, such as wound healing, infection, and cancer. Recent advances in cellular model systems, high-dimensional single-cell technologies, and lineage tracing have sparked a renaissance in plasticity research. Here, we discuss the definition of cell plasticity, evaluate state-of-the-art model systems and techniques to study cell-fate dynamics, and explore the application of single-cell technologies to obtain functional insights into cell plasticity in healthy and diseased tissues. The integration of advanced biomimetic model systems, single-cell technologies, and high-throughput perturbation studies is enabling a new era of research into non-genetic plasticity in metazoan systems.

Cell type evolution reconstruction across species through cell phylogenies of single-cell RNA sequencing data

The origin and evolution of cell types has emerged as a key topic in evolutionary biology. Driven by rapidly accumulating single-cell datasets, recent attempts to infer cell type evolution have largely been limited to pairwise comparisons because we lack approaches to build cell phylogenies using model-based approaches. Here we approach the challenges of applying explicit phylogenetic methods to single-cell data by using principal components as phylogenetic characters. We infer a cell phylogeny from a large, comparative single-cell dataset of eye cells from five distantly related mammals. Robust cell type clades enable us to provide a phylogenetic, rather than phenetic, definition of cell type, allowing us to forgo marker genes and phylogenetically classify cells by topology. We further observe evolutionary relationships between diverse vessel endothelia and identify the myelinating and non-myelinating Schwann cells as sister cell types. Finally, we examine principal component loadings and describe the gene expression dynamics underlying the function and identity of cell type clades that have been conserved across the five species. A cell phylogeny provides a rigorous framework towards investigating the evolutionary history of cells and will be critical to interpret comparative single-cell datasets that aim to ask fundamental evolutionary questions.

Epigenetic regulation of mRNA mediates the phenotypic plasticity of cancer cells during metastasis and therapeutic resistance (Review)

Plasticity, the ability of cancer cells to transition between differentiation states without genomic alterations, has been recognized as a major source of intratumoral heterogeneity. It has a crucial role in cancer metastasis and treatment resistance. Thus, targeting plasticity holds tremendous promise. However, the molecular mechanisms of plasticity in cancer cells remain poorly understood. Several studies found that mRNA, which acts as a bridge linking the genetic information of DNA and protein, has an important role in translating genotypes into phenotypes. The present review provided an overview of the regulation of cancer cell plasticity occurring via changes in the transcription and editing of mRNAs. The role of the transcriptional regulation of mRNA in cancer cell plasticity was discussed, including DNA‑binding transcriptional factors, DNA methylation, histone modifications and enhancers. Furthermore, the role of mRNA editing in cancer cell plasticity was debated, including mRNA splicing and mRNA modification. In addition, the role of non‑coding (nc)RNAs in cancer plasticity was expounded, including microRNAs, long intergenic ncRNAs and circular RNAs. Finally, different strategies for targeting cancer cell plasticity to overcome metastasis and therapeutic resistance in cancer were discussed.

How important is EMT for cancer metastasis?

Epithelial-to-mesenchymal transition (EMT), a biological phenomenon of cellular plasticity initially reported in embryonic development, has been increasingly recognized for its importance in cancer progression and metastasis. Despite tremendous progress being made in the past 2 decades in our understanding of the molecular mechanism and functional importance of EMT in cancer, there are several mysteries around EMT that remain unresolved. In this Unsolved Mystery, we focus on the variety of EMT types in metastasis, cooperative and collective EMT behaviors, spatiotemporal characterization of EMT, and strategies of therapeutically targeting EMT. We also highlight new technical advances that will facilitate the efforts to elucidate the unsolved mysteries of EMT in metastasis.

Single-cell lineage tracing with endogenous markers

Resolving lineage relationships between cells in an organism provides key insights into the fate of individual cells and drives a fundamental understanding of the process of development and disease. A recent rapid increase in experimental and computational advances for detecting naturally occurring somatic nuclear and mitochondrial mutation at single-cell resolution has expanded lineage tracing from model organisms to humans. This review discusses the advantages and challenges of experimental and computational techniques for cell lineage tracing using somatic mutation as endogenous DNA barcodes to decipher the relationships between cells during development and tumour evolution. We outlook the advantages of spatial clonal evolution analysis and single-cell lineage tracing using endogenous genetic markers.

Data Storage Using DNA

The exponential growth of global data has outpaced the storage capacities of current technologies, necessitating innovative storage strategies. DNA, as a natural medium for preserving genetic information, has emerged as a highly promising candidate for next-generation storage medium. Storing data in DNA offers several advantages, including ultrahigh physical density and exceptional durability. Facilitated by significant advancements in various technologies, such as DNA synthesis, DNA sequencing, and DNA nanotechnology, remarkable progress has been made in the field of DNA data storage over the past decade. However, several challenges still need to be addressed to realize practical applications of DNA data storage. In this review, the processes and strategies of in vitro DNA data storage are first introduced, highlighting recent advancements. Next, a brief overview of in vivo DNA data storage is provided, with a focus on the various writing strategies developed to date. At last, the challenges encountered in each step of DNA data storage are summarized and promising techniques are discussed that hold great promise in overcoming these obstacles.

Cancer as a Disease of Development Gone Awry

In the 160 years since Rudolf Virchow first postulated that neoplasia arises by the same law that regulates embryonic development, scientists have come to recognize the striking overlap between the molecular and cellular programs used by cancers and embryos. Advances in cancer biology and molecular techniques have further highlighted the similarities between carcinogenesis and embryogenesis, where cellular growth, differentiation, motility, and intercellular cross talk are mediated by common drivers and regulatory networks. This review highlights the many connections linking cancer biology and developmental biology to provide a deeper understanding of how a tissue's developmental history may both enable and constrain cancer cell evolution.

Tumor-educated Gr1+CD11b+ cells drive breast cancer metastasis via OSM/IL-6/JAK-induced cancer cell plasticity

Cancer cell plasticity contributes to therapy resistance and metastasis, which represent the main causes of cancer-related death, including in breast cancer. The tumor microenvironment drives cancer cell plasticity and metastasis, and unraveling the underlying cues may provide novel strategies for managing metastatic disease. Using breast cancer experimental models and transcriptomic analyses, we show that stem cell antigen-1 positive (SCA1+) murine breast cancer cells enriched during tumor progression and metastasis had higher in vitro cancer stem cell-like properties, enhanced in vivo metastatic ability, and generated tumors rich in Gr1hiLy6G+CD11b+ cells. In turn, tumor-educated Gr1+CD11b+ (Tu-Gr1+CD11b+) cells rapidly and transiently converted low metastatic SCA1- cells into highly metastatic SCA1+ cells via secreted oncostatin M (OSM) and IL-6. JAK inhibition prevented OSM/IL-6-induced SCA1+ population enrichment, while OSM/IL-6 depletion suppressed Tu-Gr1+CD11b+-induced SCA1+ population enrichment in vitro and metastasis in vivo. Moreover, chemotherapy-selected highly metastatic 4T1 cells maintained high SCA1+ positivity through autocrine IL-6 production, and in vitro JAK inhibition blunted SCA1 positivity and metastatic capacity. Importantly, Tu-Gr1+CD11b+ cells invoked a gene signature in tumor cells predicting shorter overall survival (OS), relapse-free survival (RFS), and lung metastasis in breast cancer patients. Collectively, our data identified OSM/IL-6/JAK as a clinically relevant paracrine/autocrine axis instigating breast cancer cell plasticity and triggering metastasis.

Lineage Plasticity: The New Cancer Hallmark on the Block

Plasticity refers to the ability of cells to adopt a spectrum of states or phenotypes. In cancer, it is a critical contributor to tumor initiation, progression, invasiveness, and therapy resistance, and it has recently been recognized as an emerging cancer hallmark. Plasticity can occur as a result of cell-intrinsic factors (e.g., genetic, transcriptional, or epigenetic fluctuations), or through cell-extrinsic cues (e.g., signaling from components of the tumor microenvironment or selective pressure from therapy). Over the past decade, technological advances, analysis of patient samples, and studies in mouse model systems have led to a deeper understanding of how such plastic states come about. In this review, we discuss: (i) the definition of plasticity; (ii) methods to measure and quantify plasticity; (iii) the clinical relevance of plasticity; and (iv) therapeutic hypotheses to modulate plasticity in the clinic.

Charting Single Cell Lineage Dynamics and Mutation Networks via Homing CRISPR

Single cell lineage tracing, essential for unraveling cellular dynamics in disease evolution is critical for developing targeted therapies. CRISPR-Cas9, known for inducing permanent and cumulative mutations, is a cornerstone in lineage tracing. The novel homing guide RNA (hgRNA) technology enhances this by enabling dynamic retargeting and facilitating ongoing genetic modifications. Charting these mutations, especially through successive hgRNA edits, poses a significant challenge. Our solution, LINEMAP, is a computational framework designed to trace and map these mutations with precision. LINEMAP meticulously discerns mutation alleles at single-cell resolution and maps their complex interrelationships through a mutation evolution network. By utilizing a Markov Process model, we can predict mutation transition probabilities, revealing potential mutational routes and pathways. Our reconstruction algorithm, anchored in the Markov model's attributes, reconstructs cellular lineage pathways, shedding light on the cell's evolutionary journey to the minutiae of single-cell division. Our findings reveal an intricate network of mutation evolution paired with a predictive Markov model, advancing our capability to reconstruct single-cell lineage via hgRNA. This has substantial implications for advancing our understanding of biological mechanisms and propelling medical research forward.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

Characterizing chromosomal instability-driven cancer evolution and cell fitness at a glance

Chromosomal instability (CIN), an increased rate of chromosome segregation errors during mitosis, is a hallmark of cancer cells. CIN leads to karyotype differences between cells and thus large-scale heterogeneity among individual cancer cells; therefore, it plays an important role in cancer evolution. Studying CIN and its consequences is technically challenging, but various technologies have been developed to track karyotype dynamics during tumorigenesis, trace clonal lineages and link genomic changes to cancer phenotypes at single-cell resolution. These methods provide valuable insight not only into the role of CIN in cancer progression, but also into cancer cell fitness. In this Cell Science at a Glance article and the accompanying poster, we discuss the relationship between CIN, cancer cell fitness and evolution, and highlight techniques that can be used to study the relationship between these factors. To that end, we explore methods of assessing cancer cell fitness, particularly for chromosomally unstable cancer.

Characterizing Rare Dormant and Cycling Lineages Using the Watermelon System

Barcode-based lineage tracing approaches enable molecular characterization of clonal cell families. Barcodes that are expressed as mRNA can be used to deconvolve lineage identity from single-cell RNA sequencing transcriptional data. Here we describe the Watermelon system, which facilitates the simultaneous tracing of lineage, transcriptional, and proliferative state at a single cell level.

DNA Barcoding Technology for Lineage Recording and Tracing to Resolve Cell Fate Determination

In various biological contexts, cells receive signals and stimuli that prompt them to change their current state, leading to transitions into a future state. This change underlies the processes of development, tissue maintenance, immune response, and the pathogenesis of various diseases. Following the path of cells from their initial identity to their current state reveals how cells adapt to their surroundings and undergo transformations to attain adjusted cellular states. DNA-based molecular barcoding technology enables the documentation of a phylogenetic tree and the deterministic events of cell lineages, providing the mechanisms and timing of cell lineage commitment that can either promote homeostasis or lead to cellular dysregulation. This review comprehensively presents recently emerging molecular recording technologies that utilize CRISPR/Cas systems, base editing, recombination, and innate variable sequences in the genome. Detailing their underlying principles, applications, and constraints paves the way for the lineage tracing of every cell within complex biological systems, encompassing the hidden steps and intermediate states of organism development and disease progression.

Startle: A star homoplasy approach for CRISPR-Cas9 lineage tracing

CRISPR-Cas9-based genome editing combined with single-cell sequencing enables the tracing of the history of cell divisions, or cellular lineage, in tissues and whole organisms. Although standard phylogenetic approaches may be applied to reconstruct cellular lineage trees from this data, the unique features of the CRISPR-Cas9 editing process motivate the development of specialized models that describe the evolution of CRISPR-Cas9-induced mutations. Here, we introduce the "star homoplasy" evolutionary model that constrains a phylogenetic character to mutate at most once along a lineage, capturing the "non-modifiability" property of CRISPR-Cas9 mutations. We derive a combinatorial characterization of star homoplasy phylogenies and use this characterization to develop an algorithm, "Startle", that computes a maximum parsimony star homoplasy phylogeny. We demonstrate that Startle infers more accurate phylogenies on simulated lineage tracing data compared with existing methods and finds parsimonious phylogenies with fewer metastatic migrations on lineage tracing data from mouse metastatic lung adenocarcinoma.

ConvexML: Scalable and accurate inference of single-cell chronograms from CRISPR/Cas9 lineage tracing data

CRISPR/Cas9 gene editing technology has enabled lineage tracing for thousands of cells in vivo. However, most of the analysis of CRISPR/Cas9 lineage tracing data has so far been limited to the reconstruction of single-cell tree topologies, which depict lineage relationships between cells, but not the amount of time that has passed between ancestral cell states and the present. Time-resolved trees, known as chronograms, would allow one to study the evolutionary dynamics of cell populations at an unprecedented level of resolution. Indeed, time-resolved trees would reveal the timing of events on the tree, the relative fitness of subclones, and the dynamics underlying phenotypic changes in the cell population - among other important applications. In this work, we introduce the first scalable and accurate method to refine any given single-cell tree topology into a single-cell chronogram by estimating its branch lengths. To do this, we leverage a statistical model of CRISPR/Cas9 cutting with missing data, paired with a conservative version of maximum parsimony that reconstructs only the ancestral states that we are confident about. As part of our method, we propose a novel approach to represent and handle missing data - specifically, double-resection events - which greatly simplifies and speeds up branch length estimation without compromising quality. All this leads to a convex maximum likelihood estimation (MLE) problem that can be readily solved in seconds with off-the-shelf convex optimization solvers. To stabilize estimates in low-information regimes, we propose a simple penalized version of MLE using a minimum branch length and pseudocounts. We benchmark our method using simulations and show that it performs well on several tasks, outperforming more naive baselines. Our method, which we name 'ConvexML', is available through the cassiopeia open source Python package.

Global lineage evolution pattern of sars-cov-2 in Africa, America, Europe, and Asia: A comparative analysis of variant clusters and their relevance across continents

Objective

The objective of this study is to provide a comparative analysis of variant clusters and their relevance across Africa, America, Europe, and Asia, in order to understand the evolutionary patterns of the virus across different regions and to inform the development of targeted interventions and genomic surveillance eforts.

Methods

The study analyzed the global lineage evolution pattern of 74, 075 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from 32 countries across four continents, focusing on variant clusters and their relevance across regions. Variants were weighted according to their hierarchical level. The correlation between variants was visualized through Dimensionality reduction analysis and Pairwise Pearson's correlation. We presented a reconstructed phylogenetic tree based on correlation analysis and variant weights.

Results

The analysis revealed that each continent had distinct variant clusters and different evolutionary patterns. The Americas had two clustered variants before lineage divergence and a downstream confluence lineage, Europe had bifurcation into two global lineages with an early occurrence of certain cluster while Asia had a downstream confluence of two large lineages diverging by two distinct clusters. Based on the cluster patterns of shared variants of the SARS-CoV-2 virus, Africa demonstrated a relatively clear distinction among three distinct regions.

Conclusions

The study provides insights into the evolutionary patterns of SARS-CoV-2 and highlights the importance of international collaboration in tracking and responding to emerging variants. The study found that the global pandemic was driven by Omicron variants that evolved with significant differences between countries and regions, and with different patterns across continents.

A mouse model with high clonal barcode diversity for joint lineage, transcriptomic, and epigenomic profiling in single cells

Cellular lineage histories and their molecular states encode fundamental principles of tissue development and homeostasis. Current lineage-recording mouse models have insufficient barcode diversity and single-cell lineage coverage for profiling tissues composed of millions of cells. Here, we developed DARLIN, an inducible Cas9 barcoding mouse line that utilizes terminal deoxynucleotidyl transferase (TdT) and 30 CRISPR target sites. DARLIN is inducible, generates massive lineage barcodes across tissues, and enables the detection of edited barcodes in ∼70% of profiled single cells. Using DARLIN, we examined fate bias within developing hematopoietic stem cells (HSCs) and revealed unique features of HSC migration. Additionally, we established a protocol for joint transcriptomic and epigenomic single-cell measurements with DARLIN and found that cellular clonal memory is associated with genome-wide DNA methylation rather than gene expression or chromatin accessibility. DARLIN will enable the high-resolution study of lineage relationships and their molecular signatures in diverse tissues and physiological contexts.

Biophysical phenotype mixtures reveal advantages for tumor muscle invasion in vivo

Bladder, colon, gastric, prostate, and uterine cancers originate in organs surrounded by laminin-coated smooth muscle. In human prostate cancer, tumors that are organ confined, without extracapsular extension through muscle, have an overall cancer survival rate of up to 97% compared with 32% for metastatic disease. Our previous work modeling extracapsular extension reported the blocking of tumor invasion by mutation of a laminin-binding integrin called α6β1. Expression of the α6AA mutant resulted in a biophysical switch from cell-ECM (extracellular matrix) to cell-cell adhesion with drug sensitivity properties and an inability to invade muscle. Here we used different admixtures of α6AA and α6WT cells to test the cell heterogeneity requirements for muscle invasion. Time-lapse video microscopy revealed that tumor mixtures self-assembled into invasive networks in vitro, whereas α6AA cells assembled only as cohesive clusters. Invasion of α6AA cells into and through live muscle occurred using a 1:1 mixture of α6AA and α6WT cells. Electric cell-substrate impedance sensing measurements revealed that compared with α6AA cells, invasion-competent α6WT cells were 2.5-fold faster at closing a cell-ECM or cell-cell wound, respectively. Cell-ECM rebuilding kinetics show that an increased response occurred in mixtures since the response was eightfold greater compared with populations containing only one cell type. A synthetic cell adhesion cyclic peptide called MTI-101 completely blocked electric cell-substrate impedance sensing cell-ECM wound recovery that persisted in vitro up to 20 h after the wound. Treatment of tumor-bearing animals with 10 mg/kg MTI-101 weekly resulted in a fourfold decrease of muscle invasion by tumor and a decrease of the depth of invasion into muscle comparable to the α6AA cells. Taken together, these data suggest that mixed biophysical phenotypes of tumor cells within a population can provide functional advantages for tumor invasion into and through muscle that can be potentially inhibited by a synthetic cell adhesion molecule.

Marker-based CRISPR screening reveals a MED12-p63 interaction that activates basal identity in pancreatic ductal adenocarcinoma

The presence of basal lineage characteristics signifies hyper-aggressive human adenocarcinomas of the breast, bladder, and pancreas. However, the biochemical mechanisms that maintain this aberrant cell state are poorly understood. Here we performed marker-based genetic screens in search of factors needed to maintain basal identity in pancreatic ductal adenocarcinoma (PDAC). This approach revealed MED12 as a powerful regulator of the basal cell state in this disease. Using biochemical reconstitution and epigenomics, we show that MED12 carries out this function by bridging the transcription factor p63, a known master regulator of the basal lineage, with the Mediator complex to activate lineage-specific enhancer elements. Consistent with this finding, the growth of basal-like PDAC is hypersensitive to MED12 loss when compared to classical PDAC. Taken together, our comprehensive genetic screens have revealed a biochemical interaction that sustains basal identity in human cancer, which could serve as a target for tumor lineage-directed therapeutics.

Bayesian inference of relative fitness on high-throughput pooled competition assays

The tracking of lineage frequencies via DNA barcode sequencing enables the quantification of microbial fitness. However, experimental noise coming from biotic and abiotic sources complicates the computation of a reliable inference. We present a Bayesian pipeline to infer relative microbial fitness from high-throughput lineage tracking assays. Our model accounts for multiple sources of noise and propagates uncertainties throughout all parameters in a systematic way. Furthermore, using modern variational inference methods based on automatic differentiation, we are able to scale the inference to a large number of unique barcodes. We extend this core model to analyze multi-environment assays, replicate experiments, and barcodes linked to genotypes. On simulations, our method recovers known parameters within posterior credible intervals. This work provides a generalizable Bayesian framework to analyze lineage tracking experiments. The accompanying open-source software library enables the adoption of principled statistical methods in experimental evolution.

Cracking the pattern of tumor evolution based on single-cell copy number alterations

Copy number alterations (CNAs) are a key characteristic of tumor development and progression. The accumulation of various CNAs during tumor development plays a critical role in driving tumor evolution. Heterogeneous clones driven by distinct CNAs have different selective advantages, leading to differential patterns of tumor evolution that are essential for developing effective cancer therapies. Recent advances in single-cell sequencing technology have enabled genome-wide copy number profiling of tumor cell populations at single-cell resolution. This has made it possible to explore the evolutionary patterns of CNAs and accurately discover the mechanisms of intra-tumor heterogeneity. Here, we propose a two-step statistical approach that distinguishes neutral, linear, branching and punctuated evolutionary patterns for a tumor cell population based on single-cell copy number profiles. We assessed our approach using a variety of simulated and real single-cell genomic and transcriptomic datasets, demonstrating its high accuracy and robustness in predicting tumor evolutionary patterns. We applied our approach to single-cell DNA sequencing data from 20 breast cancer patients and observed that punctuated evolution is the dominant evolutionary pattern in breast cancer. Similar conclusions were drawn when applying the approach to single-cell RNA sequencing data obtained from 132 various cancer patients. Moreover, we found that differential immune cell infiltration is associated with specific evolutionary patterns. The source code of our study is available at https://github.com/FangWang-SYSU/PTEM.

Single-cell phylogenies reveal changes in the evolutionary rate within cancer and healthy tissues

Cell lineages accumulate somatic mutations during organismal development, potentially leading to pathological states. The rate of somatic evolution within a cell population can vary due to multiple factors, including selection, a change in the mutation rate, or differences in the microenvironment. Here, we developed a statistical test called the Poisson Tree (PT) test to detect varying evolutionary rates among cell lineages, leveraging the phylogenetic signal of single-cell DNA sequencing (scDNA-seq) data. We applied the PT test to 24 healthy and cancer samples, rejecting a constant evolutionary rate in 11 out of 15 cancer and five out of nine healthy scDNA-seq datasets. In six cancer datasets, we identified subclonal mutations in known driver genes that could explain the rate accelerations of particular cancer lineages. Our findings demonstrate the efficacy of scDNA-seq for studying somatic evolution and suggest that cell lineages often evolve at different rates within cancer and healthy tissues.

Machine learning based lineage tree reconstruction improved with knowledge of higher level relationships between cells and genomic barcodes

Tracking cells as they divide and progress through differentiation is a fundamental step in understanding many biological processes, such as the development of organisms and progression of diseases. In this study, we investigate a machine learning approach to reconstruct lineage trees in experimental systems based on mutating synthetic genomic barcodes. We refine previously proposed methodology by embedding information of higher level relationships between cells and single-cell barcode values into a feature space. We test performance of the algorithm on shallow trees (up to 100 cells) and deep trees (up to 10 000 cells). Our proposed algorithm can improve tree reconstruction accuracy in comparison to reconstructions based on a maximum parsimony method, but this comes at a higher computational time requirement.

Single-Cell Evolutionary Analysis Reveals Drivers of Plasticity and Mediators of Chemoresistance in Small Cell Lung Cancer

Small cell lung cancer (SCLC) is often a heterogeneous tumor, where dynamic regulation of key transcription factors can drive multiple populations of phenotypically different cells which contribute differentially to tumor dynamics. This tumor is characterized by a very low 2-year survival rate, high rates of metastasis, and rapid acquisition of chemoresistance. The heterogeneous nature of this tumor makes it difficult to study and to treat, as it is not clear how or when this heterogeneity arises. Here we describe temporal, single-cell analysis of SCLC to investigate tumor initiation and chemoresistance in both SCLC xenografts and an autochthonous SCLC model. We identify an early population of tumor cells with high expression of AP-1 network genes that are critical for tumor growth. Furthermore, we have identified and validated the cancer testis antigens (CTA) PAGE5 and GAGE2A as mediators of chemoresistance in human SCLC. CTAs have been successfully targeted in other tumor types and may be a promising avenue for targeted therapy in SCLC.

IMPLICATIONS

Understanding the evolutionary dynamics of SCLC can shed light on key mechanisms such as cellular plasticity, heterogeneity, and chemoresistance.

Phylogenetic inference from single-cell RNA-seq data

Tumors are comprised of subpopulations of cancer cells that harbor distinct genetic profiles and phenotypes that evolve over time and during treatment. By reconstructing the course of cancer evolution, we can understand the acquisition of the malignant properties that drive tumor progression. Unfortunately, recovering the evolutionary relationships of individual cancer cells linked to their phenotypes remains a difficult challenge. To address this need, we have developed PhylinSic, a method that reconstructs the phylogenetic relationships among cells linked to their gene expression profiles from single cell RNA-sequencing (scRNA-Seq) data. This method calls nucleotide bases using a probabilistic smoothing approach and then estimates a phylogenetic tree using a Bayesian modeling algorithm. We showed that PhylinSic identified evolutionary relationships underpinning drug selection and metastasis and was sensitive enough to identify subclones from genetic drift. We found that breast cancer tumors resistant to chemotherapies harbored multiple genetic lineages that independently acquired high K-Ras and β-catenin, suggesting that therapeutic strategies may need to control multiple lineages to be durable. These results demonstrated that PhylinSic can reconstruct evolution and link the genotypes and phenotypes of cells across monophyletic tumors using scRNA-Seq.

The Transcriptional and Epigenetic Landscape of Cancer Cell Lineage Plasticity

Lineage plasticity, a process whereby cells change their phenotype to take on a different molecular and/or histologic identity, is a key driver of cancer progression and therapy resistance. Although underlying genetic changes within the tumor can enhance lineage plasticity, it is predominantly a dynamic process controlled by transcriptional and epigenetic dysregulation. This review explores the transcriptional and epigenetic regulators of lineage plasticity and their interplay with other features of malignancy, such as dysregulated metabolism, the tumor microenvironment, and immune evasion. We also discuss strategies for the detection and treatment of highly plastic tumors.

SIGNIFICANCE

Lineage plasticity is a hallmark of cancer and a critical facilitator of other oncogenic features such as metastasis, therapy resistance, dysregulated metabolism, and immune evasion. It is essential that the molecular mechanisms of lineage plasticity are elucidated to enable the development of strategies to effectively target this phenomenon. In this review, we describe key transcriptional and epigenetic regulators of cancer cell plasticity, in the process highlighting therapeutic approaches that may be harnessed for patient benefit.

Modeling the effect of acquired resistance on cancer therapy outcomes

Adaptive therapy (AT) is an evolution-based treatment strategy that exploits cell-cell competition. Acquired resistance can change the competitive nature of cancer cells in a tumor, impacting AT outcomes. We aimed to determine if adaptive therapy can still be effective with cell's acquiring resistance. We developed an agent-based model for spatial tumor growth considering three different types of acquired resistance: random genetic mutations during cell division, drug-induced reversible (plastic) phenotypic changes, and drug-induced irreversible phenotypic changes. These three resistance mechanisms lead to different spatial distributions of resistant cells. To quantify the spatial distribution, we propose an extension of Ripley's K-function, Sampled Ripley's K-function (SRKF), which calculates the non-randomness of the resistance distribution over the tumor domain. Our model predicts that the emergent spatial distribution of resistance can determine the time to progression under both adaptive and continuous therapy (CT). Notably, a high rate of random genetic mutations leads to quicker progression under AT than CT due to the emergence of many small clumps of resistant cells. Drug-induced phenotypic changes accelerate tumor progression irrespective of the treatment strategy. Low-rate switching to a sensitive state reduces the benefits of AT compared to CT. Furthermore, we also demonstrated that drug-induced resistance necessitates aggressive treatment under CT, regardless of the presence of cancer-associated fibroblasts. However, there is an optimal dose that can most effectively delay tumor relapse under AT by suppressing resistance. In conclusion, this study demonstrates that diverse resistance mechanisms can shape the distribution of resistance and thus determine the efficacy of adaptive therapy.

Comprehensive spatiotemporal mapping of single-cell lineages in developing mouse brain by CRISPR-based barcoding

A fundamental interest in developmental neuroscience lies in the ability to map the complete single-cell lineages within the brain. To this end, we developed a CRISPR editing-based lineage-specific tracing (CREST) method for clonal tracing in Cre mice. We then used two complementary strategies based on CREST to map single-cell lineages in developing mouse ventral midbrain (vMB). By applying snapshotting CREST (snapCREST), we constructed a spatiotemporal lineage landscape of developing vMB and identified six progenitor archetypes that could represent the principal clonal fates of individual vMB progenitors and three distinct clonal lineages in the floor plate that specified glutamatergic, dopaminergic or both neurons. We further created pandaCREST (progenitor and derivative associating CREST) to associate the transcriptomes of progenitor cells in vivo with their differentiation potentials. We identified multiple origins of dopaminergic neurons and demonstrated that a transcriptome-defined progenitor type comprises heterogeneous progenitors, each with distinct clonal fates and molecular signatures. Therefore, the CREST method and strategies allow comprehensive single-cell lineage analysis that could offer new insights into the molecular programs underlying neural specification.

Phosphorylated nuclear DICER1 promotes open chromatin state and lineage plasticity of AT2 tumor cells in lung adenocarcinomas

KRAS/ERK pathway phosphorylates DICER1, causing its nuclear translocation, and phosphomimetic Dicer1 contributes to tumorigenesis in mice. Mechanisms through which phospho-DICER1 regulates tumor progression remain undefined. While DICER1 canonically regulates microRNAs (miRNA) and epithelial-to-mesenchymal transition (EMT), we found that phosphorylated nuclear DICER1 (phospho-nuclear DICER1) promotes late-stage tumor progression in mice with oncogenic Kras, independent of miRNAs and EMT. Instead, we observe that the murine AT2 tumor cells exhibit altered chromatin compaction, and cells from disorganized advanced tumors, but not localized tumors, express gastric genes. Collectively, this results in subpopulations of tumor cells transitioning from a restricted alveolar to a broader endodermal lineage state. In human LUADs, we observed expression of phospho-nuclear DICER1 in advanced tumors together with the expression of gastric genes. We define a multimeric chromatin-DICER1 complex composed of the Mediator complex subunit 12, CBX1, MACROH2A.1, and transcriptional regulators supporting the model that phospho-nuclear DICER1 leads to lineage reprogramming of AT2 tumor cells to mediate lung cancer progression.

Tumor-intrinsic LKB1-LIF signaling axis establishes a myeloid niche to promote immune evasion and tumor growth

Tumor mutations can influence the surrounding microenvironment leading to suppression of anti-tumor immune responses and thereby contributing to tumor progression and failure of cancer therapies. Here we use genetically engineered lung cancer mouse models and patient samples to dissect how LKB1 mutations accelerate tumor growth by reshaping the immune microenvironment. Comprehensive immune profiling of LKB1 -mutant vs wildtype tumors revealed dramatic changes in myeloid cells, specifically enrichment of Arg1 ⁺ interstitial macrophages and SiglecF ^Hi neutrophils. We discovered a novel mechanism whereby autocrine LIF signaling in Lkb1 -mutant tumors drives tumorigenesis by reprogramming myeloid cells in the immune microenvironment. Inhibiting LIF signaling in Lkb1 -mutant tumors, via gene targeting or with a neutralizing antibody, resulted in a striking reduction in Arg1 ⁺ interstitial macrophages and SiglecF ^Hi neutrophils, expansion of antigen specific T cells, and inhibition of tumor progression. Thus, targeting LIF signaling provides a new therapeutic approach to reverse the immunosuppressive microenvironment of LKB1 -mutant tumors.

Quantifying Cell-State Densities in Single-Cell Phenotypic Landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive cellular differentiation, regeneration, and disease. Here, we present Mellon, a novel computational algorithm for high-resolution estimation of cell-state densities from single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of various differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. Utilizing hematopoietic stem cell fate specification to B-cells as a case study, we present evidence implicating enhancer priming and the activation of master regulators in the emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during the inherently continuous developmental processes. Scalable and adaptable, Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide new insights into the regulatory mechanisms guiding cellular fate decisions.

High-Throughput Identification, Modeling, and Analysis of Cancer Driver Genes In Vivo

The vast number of genomic and molecular alterations in cancer pose a substantial challenge to uncovering the mechanisms of tumorigenesis and identifying therapeutic targets. High-throughput functional genomic methods in genetically engineered mouse models allow for rapid and systematic investigation of cancer driver genes. In this review, we discuss the basic concepts and tools for multiplexed investigation of functionally important cancer genes in vivo using autochthonous cancer models. Furthermore, we highlight emerging technical advances in the field, potential opportunities for future investigation, and outline a vision for integrating multiplexed genetic perturbations with detailed molecular analyses to advance our understanding of the genetic and molecular basis of cancer.

KRAS(G12D) drives lepidic adenocarcinoma through stem-cell reprogramming

Many cancers originate from stem or progenitor cells hijacked by somatic mutations that drive replication, exemplified by adenomatous transformation of pulmonary alveolar epithelial type II (AT2) cells1. Here we demonstrate a different scenario: expression of KRAS(G12D) in differentiated AT1 cells reprograms them slowly and asynchronously back into AT2 stem cells that go on to generate indolent tumours. Like human lepidic adenocarcinoma, the tumour cells slowly spread along alveolar walls in a non-destructive manner and have low ERK activity. We find that AT1 and AT2 cells act as distinct cells of origin and manifest divergent responses to concomitant WNT activation and KRAS(G12D) induction, which accelerates AT2-derived but inhibits AT1-derived adenoma proliferation. Augmentation of ERK activity in KRAS(G12D)-induced AT1 cells increases transformation efficiency, proliferation and progression from lepidic to mixed tumour histology. Overall, we have identified a new cell of origin for lung adenocarcinoma, the AT1 cell, which recapitulates features of human lepidic cancer. In so doing, we also uncover a capacity for oncogenic KRAS to reprogram a differentiated and quiescent cell back into its parent stem cell en route to adenomatous transformation. Our work further reveals that irrespective of a given cancer's current molecular profile and driver oncogene, the cell of origin exerts a pervasive and perduring influence on its subsequent behaviour.

Cellular barcoding: From developmental tracing to anti-tumor drug discovery

Clonal evolution has gained immense attention in explaining cancer cell status, history, and fate during cancer progression. Current single-cell or spatial transcriptome technologies have broadened our understanding of various mechanisms underlying cancer initiation, relapse, and drug resistance. However, technical challenges still hinder a better understanding of the dynamics of distinctive phenotypic states and abnormal trajectories from normal physiological transition to malignant stages. Cellular barcoding enabled lineage tracing on parallelly massive cells at single-cell resolution through different mechanisms lately, enabling new insights into exploring developmental trajectories, cancer progression, and targeted therapies. This review summarizes the latest noteworthy and robust strategies for different types of cellular barcodes. To introduce the major characteristics, advantages and limitations of these different strategies, this review will further guide in choosing or improving cellular barcoding technologies and their applications in cancer research.

Lineage motifs: developmental modules for control of cell type proportions

In multicellular organisms, cell types must be produced and maintained in appropriate proportions. One way this is achieved is through committed progenitor cells that produce specific sets of descendant cell types. However, cell fate commitment is probabilistic in most contexts, making it difficult to infer progenitor states and understand how they establish overall cell type proportions. Here, we introduce Lineage Motif Analysis (LMA), a method that recursively identifies statistically overrepresented patterns of cell fates on lineage trees as potential signatures of committed progenitor states. Applying LMA to published datasets reveals spatial and temporal organization of cell fate commitment in zebrafish and rat retina and early mouse embryo development. Comparative analysis of vertebrate species suggests that lineage motifs facilitate adaptive evolutionary variation of retinal cell type proportions. LMA thus provides insight into complex developmental processes by decomposing them into simpler underlying modules.

Conserved transcriptional connectivity of regulatory T cells in the tumor microenvironment informs new combination cancer therapy strategies

While regulatory T (Treg) cells are traditionally viewed as professional suppressors of antigen presenting cells and effector T cells in both autoimmunity and cancer, recent findings of distinct Treg cell functions in tissue maintenance suggest that their regulatory purview extends to a wider range of cells and is broader than previously assumed. To elucidate tumoral Treg cell 'connectivity' to diverse tumor-supporting accessory cell types, we explored immediate early changes in their single-cell transcriptomes upon punctual Treg cell depletion in experimental lung cancer and injury-induced inflammation. Before any notable T cell activation and inflammation, fibroblasts, endothelial and myeloid cells exhibited pronounced changes in their gene expression in both cancer and injury settings. Factor analysis revealed shared Treg cell-dependent gene programs, foremost, prominent upregulation of VEGF and CCR2 signaling-related genes upon Treg cell deprivation in either setting, as well as in Treg cell-poor versus Treg cell-rich human lung adenocarcinomas. Accordingly, punctual Treg cell depletion combined with short-term VEGF blockade showed markedly improved control of PD-1 blockade-resistant lung adenocarcinoma progression in mice compared to the corresponding monotherapies, highlighting a promising factor-based querying approach to elucidating new rational combination treatments of solid organ cancers.

Unraveling the Drivers of Tumorigenesis in the Context of Evolution: Theoretical Models and Bioinformatics Tools

Cancer originates from somatic cells that have accumulated mutations. These mutations alter the phenotype of the cells, allowing them to escape homeostatic regulation that maintains normal cell numbers. The emergence of malignancies is an evolutionary process in which the random accumulation of somatic mutations and sequential selection of dominant clones cause cancer cells to proliferate. The development of technologies such as high-throughput sequencing has provided a powerful means to measure subclonal evolutionary dynamics across space and time. Here, we review the patterns that may be observed in cancer evolution and the methods available for quantifying the evolutionary dynamics of cancer. An improved understanding of the evolutionary trajectories of cancer will enable us to explore the molecular mechanism of tumorigenesis and to design tailored treatment strategies.

Toward a systems-level probing of tumor clonality

Cancer has been described as a genetic disease that clonally evolves in the face of selective pressures imposed by cell-intrinsic and extrinsic factors. Although classical models based on genetic data predominantly propose Darwinian mechanisms of cancer evolution, recent single-cell profiling of cancers has described unprecedented heterogeneity in tumors providing support for alternative models of branched and neutral evolution through both genetic and non-genetic mechanisms. Emerging evidence points to a complex interplay between genetic, non-genetic, and extrinsic environmental factors in shaping the evolution of tumors. In this perspective, we briefly discuss the role of cell-intrinsic and extrinsic factors that shape clonal behaviors during tumor progression, metastasis, and drug resistance. Taking examples of pre-malignant states associated with hematological malignancies and esophageal cancer, we discuss recent paradigms of tumor evolution and prospective approaches to further enhance our understanding of this spatiotemporally regulated process.

m6 A-Dependent Modulation via IGF2BP3/MCM5/Notch Axis Promotes Partial EMT and LUAD Metastasis

The importance of mRNA N6-methyladenosine (m⁶ A) modification during tumor metastasis is controversial as it plays distinct roles in different biological contexts. Moreover, how cancer cell plasticity is shaped by m⁶ A modification is interesting but remains uncharacterized. Here, this work shows that m⁶ A reader insulin like growth factor 2 mRNA binding protein 3 (IGF2BP3) is remarkably upregulated in metastatic lung adenocarcinoma (LUAD) and indicates worse prognosis of patients. Interestingly, IGF2BP3 induces partial epithelial-mesenchymal-transition (EMT) and confers LUAD cells plasticity to metastasize through m⁶ A-dependent overactivation of Notch signaling. Mechanistically, IGF2BP3 recognized m⁶ A-modified minichromosome maintenance complex component (MCM5) mRNAs to prolong stability of them, subsequently upregulating MCM5 protein, which competitively inhibits SIRT1-mediated deacetylation of Notch1 intracellular domain (NICD1), stabilizes NICD1 protein and contributes to m⁶ A-dependent IGF2BP3-mediated cellular plasticity. Notably, a tight correlation of the IGF2BP3/MCM5/Notch axis is evidenced in clinical LUAD specimens. Therefore, this study elucidates a critical role of m⁶ A modification on LUAD cell plasticity in fostering tumor metastasis via the above axis, providing potential targets for metastatic LUAD.

Gene networks reveal stem-cell state convergence during preneoplasia and progression to malignancy in multistage skin carcinogenesis

Adult mammalian stem cells play critical roles in normal tissue homeostasis, as well as in tumor development, by contributing to cell heterogeneity, plasticity, and development of drug resistance. The relationship between different types of normal and cancer stem cells is highly controversial and poorly understood. Here, we carried out gene expression network analysis of normal and tumor samples from genetically heterogeneous mice to create network metagenes for visualization of stem-cell networks, rather than individual stem-cell markers, at the single-cell level during multistage carcinogenesis. We combined this approach with lineage tracing and single-cell RNASeq of stem cells and their progeny, identifying a previously unrecognized hierarchy in which Lgr6+ stem cells from tumors generate progeny that express a range of other stem-cell markers including Sox2, Pitx1, Foxa1, Klf5, and Cd44. Our data identify a convergence of multiple stem-cell and tumor-suppressor pathways in benign tumor cells expressing markers of lineage plasticity and oxidative stress. This same single-cell population expresses network metagenes corresponding to markers of cancer drug resistance in human tumors of the skin, lung and prostate. Treatment of mouse squamous carcinomas in vivo with the chemotherapeutic cis-platin resulted in elevated expression of the genes that mark this cell population. Our data have allowed us to create a simplified model of multistage carcinogenesis that identifies distinct stem-cell states at different stages of tumor progression, thereby identifying networks involved in lineage plasticity, drug resistance, and immune surveillance, providing a rich source of potential targets for cancer therapy.

Dissecting metastasis using preclinical models and methods

Metastasis has long been understood to lead to the overwhelming majority of cancer-related deaths. However, our understanding of the metastatic process, and thus our ability to prevent or eliminate metastases, remains frustratingly limited. This is largely due to the complexity of metastasis, which is a multistep process that likely differs across cancer types and is greatly influenced by many aspects of the in vivo microenvironment. In this Review, we discuss the key variables to consider when designing assays to study metastasis: which source of metastatic cancer cells to use and where to introduce them into mice to address different questions of metastasis biology. We also examine methods that are being used to interrogate specific steps of the metastatic cascade in mouse models, as well as emerging techniques that may shed new light on previously inscrutable aspects of metastasis. Finally, we explore approaches for developing and using anti-metastatic therapies, and how mouse models can be used to test them.

Recording of cellular physiological histories along optically readable self-assembling protein chains

Observing cellular physiological histories is key to understanding normal and disease-related processes. Here we describe expression recording islands-a fully genetically encoded approach that enables both continual digital recording of biological information within cells and subsequent high-throughput readout in fixed cells. The information is stored in growing intracellular protein chains made of self-assembling subunits, human-designed filament-forming proteins bearing different epitope tags that each correspond to a different cellular state or function (for example, gene expression downstream of neural activity or pharmacological exposure), allowing the physiological history to be read out along the ordered subunits of protein chains with conventional optical microscopy. We use expression recording islands to record gene expression timecourse downstream of specific pharmacological and physiological stimuli in cultured neurons and in living mouse brain, with a time resolution of a fraction of a day, over periods of days to weeks.

Modeling collective cell behavior in cancer: Perspectives from an interdisciplinary conversation

Collective cell behavior contributes to all stages of cancer progression. Understanding how collective behavior emerges through cell-cell interactions and decision-making will advance our understanding of cancer biology and provide new therapeutic approaches. Here, we summarize an interdisciplinary discussion on multicellular behavior in cancer, draw lessons from other scientific disciplines, and identify future directions.

Untangling Cellular Host-Pathogen Encounters at Infection Bottlenecks

Bacterial pathogens can invade the tissue and establish a protected intracellular niche at the site of invasion that can spread locally (e.g., microcolonies) or to systemic sites (e.g., granulomas). Invasion of the tissue and establishment of intracellular infection are rare events that are difficult to study in the in vivo setting but have critical clinical consequences, such as long-term carriage, reinfections, and emergence of antibiotic resistance. Here, I discuss Salmonella interactions with its host macrophage during early stages of infection and their critical role in determining infection outcome. The dynamics of host-pathogen interactions entail highly heterogenous host immunity, bacterial virulence, and metabolic cross talk, requiring in vivo analysis at single-cell resolution. I discuss models and single-cell approaches that provide a global understanding of the establishment of a protected intracellular niche within the tissue and the host-pathogen landscape at infection bottlenecks during early stages of infection. Studying cellular host-pathogen interactions in vivo can improve our knowledge of the trajectory of infection between the initial inoculation with a dose of pathogens and the appearance of symptoms of disease.

Metastasis

Most cancer-associated deaths occur due to metastasis, yet our understanding of metastasis as an evolving, heterogeneous, systemic disease and of how to effectively treat it is still emerging. Metastasis requires the acquisition of a succession of traits to disseminate, variably enter and exit dormancy, and colonize distant organs. The success of these events is driven by clonal selection, the potential of metastatic cells to dynamically transition into distinct states, and their ability to co-opt the immune environment. Here, we review the main principles of metastasis and highlight emerging opportunities to develop more effective therapies for metastatic cancer.

Accelerating the understanding of cancer biology through the lens of genomics

A core mission of cancer genomics is to comprehensively chart molecular underpinnings of cancer-driving events and to provide personalized therapeutic strategies. Primarily focused on cancer cells, cancer genomics studies have successfully uncovered many drivers for major cancer types. Since the emergence of cancer immune evasion as a critical cancer hallmark, the paradigm has been elevated to the holistic tumor ecosystem, with distinct cellular components and their functional states elucidated. We highlight the milestones of cancer genomics, depict the evolving path of the field, and discuss future directions in completing the understanding of the tumor ecosystem and in advancing therapeutic strategies.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.