Modeling T Cell Fate.

Many of the pathways that underlie the diversification of naive T cells into effector and memory subsets, and the maintenance of these populations, remain controversial. In recent years a variety of experimental tools have been developed that allow us to follow the fates of cells and their descendants. In this review we describe how mathematical models provide a natural language for describing the growth, loss, and differentiation of cell populations. By encoding mechanistic descriptions of cell behavior, models can help us interpret these new datasets and reveal the rules underpinning T cell fate decisions, both at steady state and during immune responses.

Origin and Heterogeneity of Tissue Myeloid Cells: A Focus on GMP-Derived Monocytes and Neutrophils.

Myeloid cells are a significant proportion of leukocytes within tissues, comprising granulocytes, monocytes, dendritic cells, and macrophages. With the identification of various myeloid cells that perform separate but complementary functions during homeostasis and disease, our understanding of tissue myeloid cells has evolved significantly. Exciting findings from transcriptomics profiling and fate-mapping mouse models have facilitated the identification of their developmental origins, maturation, and tissue-specific specializations. This review highlights the current understanding of tissue myeloid cells and the contributing factors of functional heterogeneity to better comprehend the complex and dynamic immune interactions within the healthy or inflamed tissue. Specifically, we discuss the new understanding of the contributions of granulocyte-monocyte progenitor-derived phagocytes to tissue myeloid cell heterogeneity as well as the impact of niche-specific factors on monocyte and neutrophil phenotype and function. Lastly, we explore the developing paradigm of myeloid cell heterogeneity during inflammation and disease.

Microglia and Central Nervous System-Associated Macrophages-From Origin to Disease Modulation.

The immune system of the central nervous system (CNS) consists primarily of innate immune cells. These are highly specialized macrophages found either in the parenchyma, called microglia, or at the CNS interfaces, such as leptomeningeal, perivascular, and choroid plexus macrophages. While they were primarily thought of as phagocytes, their function extends well beyond simple removal of cell debris during development and diseases. Brain-resident innate immune cells were found to be plastic, long-lived, and host to an outstanding number of risk genes for multiple pathologies. As a result, they are now considered the most suitable targets for modulating CNS diseases. Additionally, recent single-cell technologies enhanced our molecular understanding of their origins, fates, interactomes, and functional cell statesduring health and perturbation. Here, we review the current state of our understanding and challenges of the myeloid cell biology in the CNS and treatment options for related diseases.

Time-resolved fate mapping identifies the intestinal upper crypt zone as an origin of Lgr5+ crypt base columnar cells.

In the prevailing model, Lgr5+ cells are the only intestinal stem cells (ISCs) that sustain homeostatic epithelial regeneration by upward migration of progeny through elusive upper crypt transit-amplifying (TA) intermediates. Here, we identify a proliferative upper crypt population marked by Fgfbp1, in the location of putative TA cells, that is transcriptionally distinct from Lgr5+ cells. Using a kinetic reporter for time-resolved fate mapping and Fgfbp1-CreERT2 lineage tracing, we establish that Fgfbp1+ cells are multi-potent and give rise to Lgr5+ cells, consistent with their ISC function. Fgfbp1+ cells also sustain epithelial regeneration following Lgr5+ cell depletion. We demonstrate that FGFBP1, produced by the upper crypt cells, is an essential factor for crypt proliferation and epithelial homeostasis. Our findings support a model in which tissue regeneration originates from upper crypt Fgfbp1+ cells that generate progeny propagating bi-directionally along the crypt-villus axis and serve as a source of Lgr5+ cells in the crypt base.

Automated profiling of gene function during embryonic development.

Systematic functional profiling of the gene set that directs embryonic development is an important challenge. To tackle this challenge, we used 4D imaging of C. elegans embryogenesis to capture the effects of 500 gene knockdowns and developed an automated approach to compare developmental phenotypes. The automated approach quantifies features-including germ layer cell numbers, tissue position, and tissue shape-to generate temporal curves whose parameterization yields numerical phenotypic signatures. In conjunction with a new similarity metric that operates across phenotypic space, these signatures enabled the generation of ranked lists of genes predicted to have similar functions, accessible in the PhenoBank web portal, for ∼25% of essential development genes. The approach identified new gene and pathway relationships in cell fate specification and morphogenesis and highlighted the utilization of specialized energy generation pathways during embryogenesis. Collectively, the effort establishes the foundation for comprehensive analysis of the gene set that builds a multicellular organism.

An oncogenic phenoscape of colonic stem cell polarization.

Cancer cells are regulated by oncogenic mutations and microenvironmental signals, yet these processes are often studied separately. To functionally map how cell-intrinsic and cell-extrinsic cues co-regulate cell fate, we performed a systematic single-cell analysis of 1,107 colonic organoid cultures regulated by (1) colorectal cancer (CRC) oncogenic mutations, (2) microenvironmental fibroblasts and macrophages, (3) stromal ligands, and (4) signaling inhibitors. Multiplexed single-cell analysis revealed a stepwise epithelial differentiation phenoscape dictated by combinations of oncogenes and stromal ligands, spanning from fibroblast-induced Clusterin (CLU)+ revival colonic stem cells (revCSCs) to oncogene-driven LRIG1+ hyper-proliferative CSCs (proCSCs). The transition from revCSCs to proCSCs is regulated by decreasing WNT3A and TGF-ß-driven YAP signaling and increasing KRASG12D or stromal EGF/Epiregulin-activated MAPK/PI3K flux. We find that APC loss and KRASG12D collaboratively limit access to revCSCs and disrupt stromal-epithelial communication-trapping epithelia in the proCSC fate. These results reveal that oncogenic mutations dominate homeostatic differentiation by obstructing cell-extrinsic regulation of cell-fate plasticity.

What Is a Plant Cell Type in the Age of Single-Cell Biology? It's Complicated.

One of the fundamental questions in developmental biology is how a cell is specified to differentiate as a specialized cell type. Traditionally, plant cell types were defined based on their function, location, morphology, and lineage. Currently, in the age of single-cell biology, researchers typically attempt to assign plant cells to cell types by clustering them based on their transcriptomes. However, because cells are dynamic entities that progress through the cell cycle and respond to signals, the transcriptome also reflects the state of the cell at a particular moment in time, raising questions about how to define a cell type. We suggest that these complexities and dynamics of cell states are of interest and further consider the roles signaling, stochasticity, cell cycle, and mechanical forces play in plant cell fate specification. Once established, cell identity must also be maintained. With the wealth of single-cell data coming out, the field is poised to elucidate both the complexity and dynamics of cell states.

The intrinsic and extrinsic effects of TET proteins during gastrulation.

Mice deficient for all ten-eleven translocation (TET) genes exhibit early gastrulation lethality. However, separating cause and effect in such embryonic failure is challenging. To isolate cell-autonomous effects of TET loss, we used temporal single-cell atlases from embryos with partial or complete mutant contributions. Strikingly, when developing within a wild-type embryo, Tet-mutant cells retain near-complete differentiation potential, whereas embryos solely comprising mutant cells are defective in epiblast to ectoderm transition with degenerated mesoderm potential. We map de-repressions of early epiblast factors (e.g., Dppa4 and Gdf3) and failure to activate multiple signaling from nascent mesoderm (Lefty, FGF, and Notch) as likely cell-intrinsic drivers of TET loss phenotypes. We further suggest loss of enhancer demethylation as the underlying mechanism. Collectively, our work demonstrates an unbiased approach for defining intrinsic and extrinsic embryonic gene function based on temporal differentiation atlases and disentangles the intracellular effects of the demethylation machinery from its broader tissue-level ramifications.

Embryonic origins of adult pluripotent stem cells.

Although adult pluripotent stem cells (aPSCs) are found in many animal lineages, mechanisms for their formation during embryogenesis are unknown. Here, we leveraged Hofstenia miamia, a regenerative worm that possesses collectively pluripotent aPSCs called neoblasts and produces manipulable embryos. Lineage tracing and functional experiments revealed that one pair of blastomeres gives rise to cells that resemble neoblasts in distribution, behavior, and gene expression. In Hofstenia, aPSCs include transcriptionally distinct subpopulations that express markers associated with differentiated tissues; our data suggest that despite their heterogeneity, aPSCs are derived from one lineage, not from multiple tissue-specific lineages during development. Next, we combined single-cell transcriptome profiling across development with neoblast cell-lineage tracing and identified a molecular trajectory for neoblast formation that includes transcription factors Hes, FoxO, and Tbx. This identification of a cellular mechanism and molecular trajectory for aPSC formation opens the door for in vivo studies of aPSC regulation and evolution.

Mapping transcriptomic vector fields of single cells.

Single-cell (sc)RNA-seq, together with RNA velocity and metabolic labeling, reveals cellular states and transitions at unprecedented resolution. Fully exploiting these data, however, requires kinetic models capable of unveiling governing regulatory functions. Here, we introduce an analytical framework dynamo (https://github.com/aristoteleo/dynamo-release), which infers absolute RNA velocity, reconstructs continuous vector fields that predict cell fates, employs differential geometry to extract underlying regulations, and ultimately predicts optimal reprogramming paths and perturbation outcomes. We highlight dynamo's power to overcome fundamental limitations of conventional splicing-based RNA velocity analyses to enable accurate velocity estimations on a metabolically labeled human hematopoiesis scRNA-seq dataset. Furthermore, differential geometry analyses reveal mechanisms driving early megakaryocyte appearance and elucidate asymmetrical regulation within the PU.1-GATA1 circuit. Leveraging the least-action-path method, dynamo accurately predicts drivers of numerous hematopoietic transitions. Finally, in silico perturbations predict cell-fate diversions induced by gene perturbations. Dynamo, thus, represents an important step in advancing quantitative and predictive theories of cell-state transitions.

A single-embryo, single-cell time-resolved model for mouse gastrulation.

Mouse embryonic development is a canonical model system for studying mammalian cell fate acquisition. Recently, single-cell atlases comprehensively charted embryonic transcriptional landscapes, yet inference of the coordinated dynamics of cells over such atlases remains challenging. Here, we introduce a temporal model for mouse gastrulation, consisting of data from 153 individually sampled embryos spanning 36 h of molecular diversification. Using algorithms and precise timing, we infer differentiation flows and lineage specification dynamics over the embryonic transcriptional manifold. Rapid transcriptional bifurcations characterize the commitment of early specialized node and blood cells. However, for most lineages, we observe combinatorial multi-furcation dynamics rather than hierarchical transcriptional transitions. In the mesoderm, dozens of transcription factors combinatorially regulate multifurcations, as we exemplify using time-matched chimeric embryos of Foxc1/Foxc2 mutants. Our study rejects the notion of differentiation being governed by a series of binary choices, providing an alternative quantitative model for cell fate acquisition.

Stem-like intestinal Th17 cells give rise to pathogenic effector T cells during autoimmunity.

While intestinal Th17 cells are critical for maintaining tissue homeostasis, recent studies have implicated their roles in the development of extra-intestinal autoimmune diseases including multiple sclerosis. However, the mechanisms by which tissue Th17 cells mediate these dichotomous functions remain unknown. Here, we characterized the heterogeneity, plasticity, and migratory phenotypes of tissue Th17 cells in vivo by combined fate mapping with profiling of the transcriptomes and TCR clonotypes of over 84,000 Th17 cells at homeostasis and during CNS autoimmune inflammation. Inter- and intra-organ single-cell analyses revealed a homeostatic, stem-like TCF1+ IL-17+ SLAMF6+ population that traffics to the intestine where it is maintained by the microbiota, providing a ready reservoir for the IL-23-driven generation of encephalitogenic GM-CSF+ IFN-γ+ CXCR6+ T cells. Our study defines a direct in vivo relationship between IL-17+ non-pathogenic and GM-CSF+ and IFN-γ+ pathogenic Th17 populations and provides a mechanism by which homeostatic intestinal Th17 cells direct extra-intestinal autoimmune disease.

Senescence: An Identity Crisis Originating from Deep Within the Nucleus.

Cellular senescence is implicated in a wide range of physiological and pathological conditions throughout an organism's entire lifetime. In particular, it has become evident that senescence plays a causative role in aging and age-associated disorders. This is not due simply to the loss of function of senescent cells. Instead, the substantial alterations of the cellular activities of senescent cells, especially the array of secretory factors, impact the surrounding tissues or even entire organisms. Such non-cell-autonomous functionality is largely coordinated by tissue-specific genes, constituting a cell fate-determining state. Senescence can be viewed as a gain-of-function phenotype or a process of cell identity shift. Cellular functionality or lineage-specific gene expression is tightly linked to the cell type-specific epigenetic landscape, reinforcing the heterogeneity of senescence across cell types. Here, we aim to define the senescence cellular functionality and epigenetic features that may contribute to the gain-of-function phenotype.

MorphoSeq: Full Single-Cell Transcriptome Dynamics Up to Gastrulation in a Chordate.

Single-cell RNA sequencing (scRNA-seq) provides a leap forward in resolving cellular diversity and developmental trajectories but fails to comprehensively delineate the spatial organization and precise cellular makeup of individual embryos. Here, we reconstruct from scRNA-seq and light sheet imaging data a canonical digital embryo that captures the genome-wide gene expression trajectory of every single cell at every cell division in the 18 lineages up to gastrulation in the ascidian Phallusia mammillata. By using high-coverage scRNA-seq, we devise a computational framework that stratifies single cells of individual embryos into cell types without prior knowledge. Unbiased transcriptome data analysis mapped each cell's physical position and lineage history, yielding the complete history of gene expression at the genome-wide level for every single cell in a developing embryo. A comparison of individual embryos reveals both extensive reproducibility between symmetric embryo sides and a large inter-embryonic variability due to small differences in embryogenesis timing.

TBL1XR1 Mutations Drive Extranodal Lymphoma by Inducing a Pro-tumorigenic Memory Fate.

The most aggressive B cell lymphomas frequently manifest extranodal distribution and carry somatic mutations in the poorly characterized gene TBL1XR1. Here, we show that TBL1XR1 mutations skew the humoral immune response toward generating abnormal immature memory B cells (MB), while impairing plasma cell differentiation. At the molecular level, TBL1XR1 mutants co-opt SMRT/HDAC3 repressor complexes toward binding the MB cell transcription factor (TF) BACH2 at the expense of the germinal center (GC) TF BCL6, leading to pre-memory transcriptional reprogramming and cell-fate bias. Upon antigen recall, TBL1XR1 mutant MB cells fail to differentiate into plasma cells and instead preferentially reenter new GC reactions, providing evidence for a cyclic reentry lymphomagenesis mechanism. Ultimately, TBL1XR1 alterations lead to a striking extranodal immunoblastic lymphoma phenotype that mimics the human disease. Both human and murine lymphomas feature expanded MB-like cell populations, consistent with a MB-cell origin and delineating an unforeseen pathway for malignant transformation of the immune system.

Patterns of Early p21 Dynamics Determine Proliferation-Senescence Cell Fate after Chemotherapy.

Chemotherapy is designed to induce cell death. However, at non-lethal doses, cancer cells can choose to remain proliferative or become senescent. The slow development of senescence makes studying this decision challenging. Here, by analyzing single-cell p21 dynamics before, during, and days after drug treatment, we link three distinct patterns of early p21 dynamics to final cell fate. Surprisingly, while high p21 expression is classically associated with senescence, we find the opposite at early times during drug treatment: most senescence-fated cells express much lower p21 levels than proliferation-fated cells. We demonstrate that these dynamics lead to a p21 "Goldilocks zone" for proliferation, in which modest increases of p21 expression can lead to an undesirable increase of cancer cell proliferation. Our study identifies a counter-intuitive role for early p21 dynamics in the cell-fate decision and pinpoints a source of proliferative cancer cells that can emerge after exposure to non-lethal doses of chemotherapy.

Lateral Inhibition in Cell Specification Mediated by Mechanical Signals Modulating TAZ Activity.

Cell fate specification by lateral inhibition typically involves contact signaling through the Delta-Notch signaling pathway. However, whether this is the only signaling mode mediating lateral inhibition remains unclear. Here we show that in zebrafish oogenesis, a group of cells within the granulosa cell layer at the oocyte animal pole acquire elevated levels of the transcriptional coactivator TAZ in their nuclei. One of these cells, the future micropyle precursor cell (MPC), accumulates increasingly high levels of nuclear TAZ and grows faster than its surrounding cells, mechanically compressing those cells, which ultimately lose TAZ from their nuclei. Strikingly, relieving neighbor-cell compression by MPC ablation or aspiration restores nuclear TAZ accumulation in neighboring cells, eventually leading to MPC re-specification from these cells. Conversely, MPC specification is defective in taz-/- follicles. These findings uncover a novel mode of lateral inhibition in cell fate specification based on mechanical signals controlling TAZ activity.

Effector TH17 Cells Give Rise to Long-Lived TRM Cells that Are Essential for an Immediate Response against Bacterial Infection.

Adaptive immunity provides life-long protection by generating central and effector memory T cells and the most recently described tissue resident memory T (TRM) cells. However, the cellular origin of CD4 TRM cells and their contribution to host defense remain elusive. Using IL-17A tracking-fate mouse models, we found that a significant fraction of lung CD4 TRM cells derive from IL-17A-producing effector (TH17) cells following immunization with heat-killed Klebsiella pneumonia (Kp). These exTH17 TRM cells are maintained in the lung by IL-7, produced by lymphatic endothelial cells. During a memory response, neither antibodies, γδ T cells, nor circulatory T cells are sufficient for the rapid host defense required to eliminate Kp. Conversely, using parabiosis and depletion studies, we demonstrated that exTH17 TRM cells play an important role in bacterial clearance. Thus, we delineate the origin and function of airway CD4 TRM cells during bacterial infection, offering novel strategies for targeted vaccine design.

Identification of Enteroendocrine Regulators by Real-Time Single-Cell Differentiation Mapping.

Homeostatic regulation of the intestinal enteroendocrine lineage hierarchy is a poorly understood process. We resolved transcriptional changes during enteroendocrine differentiation in real time at single-cell level using a novel knockin allele of Neurog3, the master regulator gene briefly expressed at the onset of enteroendocrine specification. A bi-fluorescent reporter, Neurog3Chrono, measures time from the onset of enteroendocrine differentiation and enables precise positioning of single-cell transcriptomes along an absolute time axis. This approach yielded a definitive description of the enteroendocrine hierarchy and its sub-lineages, uncovered differential kinetics between sub-lineages, and revealed time-dependent hormonal plasticity in enterochromaffin and L cells. The time-resolved map of transcriptional changes predicted multiple novel molecular regulators. Nine of these were validated by conditional knockout in mice or CRISPR modification in intestinal organoids. Six novel candidate regulators (Sox4, Rfx6, Tox3, Myt1, Runx1t1, and Zcchc12) yielded specific enteroendocrine phenotypes. Our time-resolved single-cell transcriptional map presents a rich resource to unravel enteroendocrine differentiation.

Optimal Decoding of Cellular Identities in a Genetic Network.

In developing organisms, spatially prescribed cell identities are thought to be determined by the expression levels of multiple genes. Quantitative tests of this idea, however, require a theoretical framework capable of exposing the rules and precision of cell specification over developmental time. We use the gap gene network in the early fly embryo as an example to show how expression levels of the four gap genes can be jointly decoded into an optimal specification of position with 1% accuracy. The decoder correctly predicts, with no free parameters, the dynamics of pair-rule expression patterns at different developmental time points and in various mutant backgrounds. Precise cellular identities are thus available at the earliest stages of development, contrasting the prevailing view of positional information being slowly refined across successive layers of the patterning network. Our results suggest that developmental enhancers closely approximate a mathematically optimal decoding strategy.