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.

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.

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.

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.

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.

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.

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.

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.

Human Pluripotency Is Initiated and Preserved by a Unique Subset of Founder Cells.

The assembly of organized colonies is the earliest manifestation in the derivation or induction of pluripotency in vitro. However, the necessity and origin of this assemblance is unknown. Here, we identify human pluripotent founder cells (hPFCs) that initiate, as well as preserve and establish, pluripotent stem cell (PSC) cultures. PFCs are marked by N-cadherin expression (NCAD+) and reside exclusively at the colony boundary of primate PSCs. As demonstrated by functional analysis, hPFCs harbor the clonogenic capacity of PSC cultures and emerge prior to commitment events or phenotypes associated with pluripotent reprogramming. Comparative single-cell analysis with pre- and post-implantation primate embryos revealed hPFCs share hallmark properties with primitive endoderm (PrE) and can be regulated by non-canonical Wnt signaling. Uniquely informed by primate embryo organization in vivo, our study defines a subset of founder cells critical to the establishment pluripotent state.

Inadequate DNA Damage Repair Promotes Mammary Transdifferentiation, Leading to BRCA1 Breast Cancer.

Loss of BRCA1 p220 function often results in basal-like breast cancer (BLBC), but the underlying disease mechanism is largely opaque. In mammary epithelial cells (MECs), BRCA1 interacts with multiple proteins, including NUMB and HES1, to form complexes that participate in interstrand crosslink (ICL) DNA repair and MEC differentiation control. Unrepaired ICL damage results in aberrant transdifferentiation to a mesenchymal state of cultured, human basal-like MECs and to a basal/mesenchymal state in primary mouse luminal MECs. Loss of BRCA1, NUMB, or HES1 or chemically induced ICL damage in primary murine luminal MECs results in persistent DNA damage that triggers luminal to basal/mesenchymal transdifferentiation. In vivo single-cell analysis revealed a time-dependent evolution from normal luminal MECs to luminal progenitor-like tumor cells with basal/mesenchymal transdifferentiation during murine BRCA1 BLBC development. Growing DNA damage accompanied this malignant transformation.

Reversible, tunable epigenetic silencing of TCF1 generates flexibility in the T cell memory decision.

The immune system encodes information about the severity of a pathogenic threat in the quantity and type of memory cells it forms. This encoding emerges from lymphocyte decisions to maintain or lose self-renewal and memory potential during a challenge. By tracking CD8+ T cells at the single-cell and clonal lineage level using time-resolved transcriptomics, quantitative live imaging, and an acute infection model, we find that T cells will maintain or lose memory potential early after antigen recognition. However, following pathogen clearance, T cells may regain memory potential if initially lost. Mechanistically, this flexibility is implemented by a stochastic cis-epigenetic switch that tunably and reversibly silences the memory regulator, TCF1, in response to stimulation. Mathematical modeling shows how this flexibility allows memory T cell numbers to scale robustly with pathogen virulence and immune response magnitudes. We propose that flexibility and stochasticity in cellular decisions ensure optimal immune responses against diverse threats.

Identity Noise and Adipogenic Traits Characterize Dermal Fibroblast Aging.

During aging, stromal functions are thought to be impaired, but little is known whether this stems from changes of fibroblasts. Using population- and single-cell transcriptomics, as well as long-term lineage tracing, we studied whether murine dermal fibroblasts are altered during physiological aging under different dietary regimes that affect longevity. We show that the identity of old fibroblasts becomes undefined, with the fibroblast states present in young skin no longer clearly demarcated. In addition, old fibroblasts not only reduce the expression of genes involved in the formation of the extracellular matrix, but also gain adipogenic traits, paradoxically becoming more similar to neonatal pro-adipogenic fibroblasts. These alterations are sensitive to systemic metabolic changes: long-term caloric restriction reversibly prevents them, whereas a high-fat diet potentiates them. Our results therefore highlight loss of cell identity and the acquisition of adipogenic traits as a mechanism underlying cellular aging, which is influenced by systemic metabolism.

Asymmetric Expression of LincGET Biases Cell Fate in Two-Cell Mouse Embryos.

In early mammalian embryos, it remains unclear how the first cell fate bias is initially triggered and amplified toward cell fate segregation. Here, we report that a long noncoding RNA, LincGET, is transiently and asymmetrically expressed in the nucleus of two- to four-cell mouse embryos. Overexpression of LincGET in one of the two-cell blastomeres biases its progeny predominantly toward the inner cell mass (ICM) fate. Mechanistically, LincGET physically binds to CARM1 and promotes the nuclear localization of CARM1, which can further increase the level of H3 methylation at Arginine 26 (H3R26me), activate ICM-specific gene expression, upregulate transposons, and increase global chromatin accessibility. Simultaneous overexpression of LincGET and depletion of Carm1 no longer biased embryonic fate, indicating that the effect of LincGET in directing ICM lineage depends on CARM1. Thus, our data identify LincGET as one of the earliest known lineage regulators to bias cell fate in mammalian 2-cell embryos.

CARM1 and Paraspeckles Regulate Pre-implantation Mouse Embryo Development.

Nuclear architecture has never been carefully examined during early mammalian development at the stages leading to establishment of the embryonic and extra-embryonic lineages. Heterogeneous activity of the methyltransferase CARM1 during these stages results in differential methylation of histone H3R26 to modulate establishment of these two lineages. Here we show that CARM1 accumulates in nuclear granules at the 2- to 4-cell stage transition in the mouse embryo, with the majority corresponding to paraspeckles. The paraspeckle component Neat1 and its partner p54nrb are required for CARM1's association with paraspeckles and for H3R26 methylation. Conversely, CARM1 also influences paraspeckle organization. Depletion of Neat1 or p54nrb results in arrest at the 16- to 32-cell stage, with elevated expression of transcription factor Cdx2, promoting differentiation into the extra-embryonic lineage. This developmental arrest occurs at an earlier stage than following CARM1 depletion, indicating that paraspeckles act upstream of CARM1 but also have additional earlier roles in fate choice.

Systemic Human ILC Precursors Provide a Substrate for Tissue ILC Differentiation.

Innate lymphoid cells (ILCs) represent innate versions of T helper and cytotoxic T cells that differentiate from committed ILC precursors (ILCPs). How ILCPs give rise to mature tissue-resident ILCs remains unclear. Here, we identify circulating and tissue ILCPs in humans that fail to express the transcription factors and cytokine outputs of mature ILCs but have these signature loci in an epigenetically poised configuration. Human ILCPs robustly generate all ILC subsets in vitro and in vivo. While human ILCPs express low levels of retinoic acid receptor (RAR)-related orphan receptor C (RORC) transcripts, these cells are found in RORC-deficient patients and retain potential for EOMES+ natural killer (NK) cells, interferon gamma-positive (IFN-γ+) ILC1s, interleukin (IL)-13+ ILC2s, and for IL-22+, but not for IL-17A+ ILC3s. Our results support a model of tissue ILC differentiation ("ILC-poiesis"), whereby diverse ILC subsets are generated in situ from systemically distributed ILCPs in response to local environmental signals.

Fatty Acids Regulate Germline Sex Determination through ACS-4-Dependent Myristoylation.

Fat metabolism has been linked to fertility and reproductive adaptation in animals and humans, and environmental sex determination potentially plays a role in the process. To investigate the impact of fatty acids (FA) on sex determination and reproductive development, we examined and observed an impact of FA synthesis and mobilization by lipolysis in somatic tissues on oocyte fate in Caenorhabditis elegans. The subsequent genetic analysis identified ACS-4, an acyl-CoA synthetase and its FA-CoA product, as key germline factors that mediate the role of FA in promoting oocyte fate through protein myristoylation. Further tests indicated that ACS-4-dependent protein myristoylation perceives and translates the FA level into regulatory cues that modulate the activities of MPK-1/MAPK and key factors in the germline sex-determination pathway. These findings, including a similar role of ACS-4 in a male/female species, uncover a likely conserved mechanism by which FA, an environmental factor, regulates sex determination and reproductive development.

EGFR Ligands Differentially Stabilize Receptor Dimers to Specify Signaling Kinetics.

Epidermal growth factor receptor (EGFR) regulates many crucial cellular programs, with seven different activating ligands shaping cell signaling in distinct ways. Using crystallography and other approaches, we show how the EGFR ligands epiregulin (EREG) and epigen (EPGN) stabilize different dimeric conformations of the EGFR extracellular region. As a consequence, EREG or EPGN induce less stable EGFR dimers than EGF-making them partial agonists of EGFR dimerization. Unexpectedly, this weakened dimerization elicits more sustained EGFR signaling than seen with EGF, provoking responses in breast cancer cells associated with differentiation rather than proliferation. Our results reveal how responses to different EGFR ligands are defined by receptor dimerization strength and signaling dynamics. These findings have broad implications for understanding receptor tyrosine kinase (RTK) signaling specificity. Our results also suggest parallels between partial and/or biased agonism in RTKs and G-protein-coupled receptors, as well as new therapeutic opportunities for correcting RTK signaling output.