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.

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.

Measuring and Modeling Single-Cell Heterogeneity and Fate Decision in Mouse Embryos.

Cellular heterogeneity is a property of any living system; however, its relationship with cellular fate decision remains an open question. Recent technological advances have enabled valuable insights, especially in complex systems such as the mouse embryo. In this review, we discuss recent studies that characterize cellular heterogeneity at different levels during mouse development, from the two-cell stage up to gastrulation. In addition to key experimental findings, we review mathematical modeling approaches that help researchers interpret these findings. Disentangling the role of heterogeneity in cell fate decision will likely rely on the refined integration of experiments, large-scale omics data, and mathematical modeling, complemented by the use of synthetic embryos and gastruloids as promising in vitro models.

A mathematical framework for measuring and tuning tempo in developmental gene regulatory networks.

Embryo development is a dynamic process governed by the regulation of timing and sequences of gene expression, which control the proper growth of the organism. Although many genetic programmes coordinating these sequences are common across species, the timescales of gene expression can vary significantly among different organisms. Currently, substantial experimental efforts are focused on identifying molecular mechanisms that control these temporal aspects. In contrast, the capacity of established mathematical models to incorporate tempo control while maintaining the same dynamical landscape remains less understood. Here, we address this gap by developing a mathematical framework that links the functionality of developmental programmes to the corresponding gene expression orbits (or landscapes). This unlocks the ability to find tempo differences as perturbations in the dynamical system that preserve its orbits. We demonstrate that this framework allows for the prediction of molecular mechanisms governing tempo, through both numerical and analytical methods. Our exploration includes two case studies: a generic network featuring coupled production and degradation, with a particular application to neural progenitor differentiation; and the repressilator. In the latter, we illustrate how altering the dimerisation rates of transcription factors can decouple the tempo from the shape of the resulting orbits. We conclude by highlighting how the identification of orthogonal molecular mechanisms for tempo control can inform the design of circuits with specific orbits and tempos.

CoBATCH for High-Throughput Single-Cell Epigenomic Profiling.

An efficient, generalizable method for genome-wide mapping of single-cell histone modifications or chromatin-binding proteins is lacking. Here, we develop CoBATCH, combinatorial barcoding and targeted chromatin release, for single-cell profiling of genomic distribution of chromatin-binding proteins in cell culture and tissue. Protein A in fusion to Tn5 transposase is enriched through specific antibodies to genomic regions, and Tn5 generates indexed chromatin fragments ready for library preparation and sequencing. Importantly, this strategy enables not only low-input epigenomic profiling in intact tissues but also measures scalable up to tens of thousands of single cells per experiment under both native and cross-linked conditions. CoBATCH produces ∼12,000 reads/cell with extremely low background. Mapping of endothelial cell lineages from ten embryonic mouse organs through CoBATCH allows for efficient deciphering of epigenetic heterogeneity of cell populations and cis-regulatory mechanisms. Thus, obviating specialized devices, CoBATCH is broadly applicable and easily deployable for single-cell profiling of protein-DNA interactions.

Resolving Cell Fate Decisions during Somatic Cell Reprogramming by Single-Cell RNA-Seq.

Somatic cells can be reprogrammed into induced pluripotent stem cells (iPSCs), which is a highly heterogeneous process. Here we report the cell fate continuum during somatic cell reprogramming at single-cell resolution. We first develop SOT to analyze cell fate continuum from Oct4/Sox2/Klf4- or OSK-mediated reprogramming and show that cells bifurcate into two categories, reprogramming potential (RP) or non-reprogramming (NR). We further show that Klf4 contributes to Cd34+/Fxyd5+/Psca+ keratinocyte-like NR fate and that IFN-γ impedes the final transition to chimera-competent pluripotency along the RP cells. We analyze more than 150,000 single cells from both OSK and chemical reprograming and identify additional NR/RP bifurcation points. Our work reveals a generic bifurcation model for cell fate decisions during somatic cell reprogramming that may be applicable to other systems and inspire further improvements for reprogramming.

Nutrient regulation of development and cell fate decisions.

Diet contributes to health at all stages of life, from embryonic development to old age. Nutrients, including vitamins, amino acids, lipids and sugars, have instructive roles in directing cell fate and function, maintaining stem cell populations, tissue homeostasis and alleviating the consequences of aging. This Review highlights recent findings that illuminate how common diets and specific nutrients impact cell fate decisions in healthy and disease contexts. We also draw attention to new models, technologies and resources that help to address outstanding questions in this emerging field and may lead to dietary approaches that promote healthy development and improve disease treatments.

The circadian clock mediates daily bursts of cell differentiation by periodically restricting cell-differentiation commitment.

Most mammalian cells have an intrinsic circadian clock that coordinates metabolic activity with the daily rest and wake cycle. The circadian clock is known to regulate cell differentiation, but how continuous daily oscillations of the internal clock can control a much longer, multiday differentiation process is not known. Here, we simultaneously monitor circadian clock and adipocyte-differentiation progression live in single cells. Strikingly, we find a bursting behavior in the cell population whereby individual preadipocytes commit to differentiate primarily during a 12-h window each day, corresponding to the time of rest. Daily gating occurs because cells irreversibly commit to differentiate within only a few hours, which is much faster than the rest phase and the overall multiday differentiation process. The daily bursts in differentiation commitment result from a differentiation-stimulus driven variable and slow increase in expression of PPARG, the master regulator of adipogenesis, overlaid with circadian boosts in PPARG expression driven by fast, clock-driven PPARG regulators such as CEBPA. Our finding of daily bursts in cell differentiation only during the circadian cycle phase corresponding to evening in humans is broadly relevant, given that most differentiating somatic cells are regulated by the circadian clock. Having a restricted time each day when differentiation occurs may open therapeutic strategies to use timed treatment relative to the clock to promote tissue regeneration.

Volumetric compression develops noise-driven single-cell heterogeneity.

Recent studies have revealed that extensive heterogeneity of biological systems arises through various routes ranging from intracellular chromosome segregation to spatiotemporally varying biochemical stimulations. However, the contribution of physical microenvironments to single-cell heterogeneity remains largely unexplored. Here, we show that a homogeneous population of non-small-cell lung carcinoma develops into heterogeneous subpopulations upon application of a homogeneous physical compression, as shown by single-cell transcriptome profiling. The generated subpopulations stochastically gain the signature genes associated with epithelial-mesenchymal transition (EMT; VIM, CDH1, EPCAM, ZEB1, and ZEB2) and cancer stem cells (MKI67, BIRC5, and KLF4), respectively. Trajectory analysis revealed two bifurcated paths as cells evolving upon the physical compression, along each path the corresponding signature genes (epithelial or mesenchymal) gradually increase. Furthermore, we show that compression increases gene expression noise, which interplays with regulatory network architecture and thus generates differential cell-fate outcomes. The experimental observations of both single-cell sequencing and single-molecule fluorescent in situ hybridization agrees well with our computational modeling of regulatory network in the EMT process. These results demonstrate a paradigm of how mechanical stimulations impact cell-fate determination by altering transcription dynamics; moreover, we show a distinct path that the ecology and evolution of cancer interplay with their physical microenvironments from the view of mechanobiology and systems biology, with insight into the origin of single-cell heterogeneity.

Bacteriophage self-counting in the presence of viral replication.

When host cells are in low abundance, temperate bacteriophages opt for dormant (lysogenic) infection. Phage lambda implements this strategy by increasing the frequency of lysogeny at higher multiplicity of infection (MOI). However, it remains unclear how the phage reliably counts infecting viral genomes even as their intracellular number increases because of replication. By combining theoretical modeling with single-cell measurements of viral copy number and gene expression, we find that instead of hindering lambda's decision, replication facilitates it. In a nonreplicating mutant, viral gene expression simply scales with MOI rather than diverging into lytic (virulent) and lysogenic trajectories. A similar pattern is followed during early infection by wild-type phage. However, later in the infection, the modulation of viral replication by the decision genes amplifies the initially modest gene expression differences into divergent trajectories. Replication thus ensures the optimal decision-lysis upon single-phage infection and lysogeny at higher MOI.

The role of flavonoids in the regulation of epithelial-mesenchymal transition in cancer: A review on targeting signaling pathways and metastasis.

One of the hallmarks of cancer is metastasis, a process that entails the spread of cancer cells to distant regions in the body, culminating in tumor formation in secondary organs. Importantly, the proinflammatory environment surrounding cancer cells further contributes to cancer cell transformation and extracellular matrix destruction. During metastasis, front-rear polarity and emergence of migratory and invasive features are manifestations of epithelial-mesenchymal transition (EMT). A variety of transcription factors (TFs) are implicated in the execution of EMT, the most prominent belonging to the Snail Family Transcriptional Repressor (SNAI) and Zinc Finger E-Box Binding Homeobox (ZEB) families of TFs. These TFs are regulated by interaction with specific microRNAs (miRNAs), as miR34 and miR200. Among the several secondary metabolites produced in plants, flavonoids constitute a major group of bioactive molecules, with several described effects including antioxidant, antiinflammatory, antidiabetic, antiobesogenic, and anticancer effects. This review scrutinizes the modulatory role of flavonoids on the activity of SNAI/ZEB TFs and on their regulatory miRNAs, miR-34, and miR-200. The modulatory role of flavonoids can attenuate mesenchymal features and stimulate epithelial features, thereby inhibiting and reversing EMT. Moreover, this modulation is concomitant with the attenuation of signaling pathways involved in diverse processes as cell proliferation, cell growth, cell cycle progression, apoptosis inhibition, morphogenesis, cell fate, cell migration, cell polarity, and wound healing. The antimetastatic potential of these versatile compounds is emerging and represents an opportunity for the synthesis of more specific and potent agents.

Ascl1 is required to specify a subset of ventromedial hypothalamic neurons.

Despite clear physiological roles, the ventromedial hypothalamus (VMH) developmental programs are poorly understood. Here, we asked whether the proneural gene achaete-scute homolog 1 (Ascl1) contributes to VMH development. Ascl1 transcripts were detected in embryonic day (E) 10.5 to postnatal day 0 VMH neural progenitors. The elimination of Ascl1 reduced the number of VMH neurons at E12.5 and E15.5, particularly within the VMH-central (VMHC) and -dorsomedial (VMHDM) subdomains, and resulted in a VMH cell fate change from glutamatergic to GABAergic. We observed a loss of Neurog3 expression in Ascl1-/- hypothalamic progenitors and an upregulation of Neurog3 when Ascl1 was overexpressed. We also demonstrated a glutamatergic to GABAergic fate switch in Neurog3-null mutant mice, suggesting that Ascl1 might act via Neurog3 to drive VMH cell fate decisions. We also showed a concomitant increase in expression of the central GABAergic fate determinant Dlx1/2 in the Ascl1-null hypothalamus. However, Ascl1 was not sufficient to induce an ectopic VMH fate when overexpressed outside the normal window of competency. Combined, Ascl1 is required but not sufficient to specify the neurotransmitter identity of VMH neurons, acting in a transcriptional cascade with Neurog3.

Oxidative stress as a key modulator of cell fate decision in osteoarthritis and osteoporosis: a narrative review.

During aging and after traumatic injuries, cartilage and bone cells are exposed to various pathophysiologic mediators, including reactive oxygen species (ROS), damage-associated molecular patterns, and proinflammatory cytokines. This detrimental environment triggers cellular stress and subsequent dysfunction, which not only contributes to the development of associated diseases, that is, osteoporosis and osteoarthritis, but also impairs regenerative processes. To counter ROS-mediated stress and reduce the overall tissue damage, cells possess diverse defense mechanisms. However, cellular antioxidative capacities are limited and thus ROS accumulation can lead to aberrant cell fate decisions, which have adverse effects on cartilage and bone homeostasis. In this narrative review, we address oxidative stress as a major driver of pathophysiologic processes in cartilage and bone, including senescence, misdirected differentiation, cell death, mitochondrial dysfunction, and impaired mitophagy by illustrating the consequences on tissue homeostasis and regeneration. Moreover, we elaborate cellular defense mechanisms, with a particular focus on oxidative stress response and mitophagy, and briefly discuss respective therapeutic strategies to improve cell and tissue protection.

Transcription factors regulate early T cell development via redeployment of other factors: Functional dynamics of constitutively required factors in cell fate decisions.

Establishment of cell lineage identity from multipotent progenitors is controlled by cooperative actions of lineage-specific and stably expressed transcription factors, combined with input from environmental signals. Lineage-specific master transcription factors activate and repress gene expression by recruiting consistently expressed transcription factors and chromatin modifiers to their target loci. Recent technical advances in genome-wide and multi-omics analysis have shed light on unexpected mechanisms that underlie more complicated actions of transcription factors in cell fate decisions. In this review, we discuss functional dynamics of stably expressed and continuously required factors, Notch and Runx family members, throughout developmental stages of early T cell development in the thymus. Pre- and post-commitment stage-specific transcription factors induce dynamic redeployment of Notch and Runx binding genomic regions. Thus, together with stage-specific transcription factors, shared transcription factors across distinct developmental stages regulate acquisition of T lineage identity.

Synchronization between Attractors: Genomic Mechanism of Cell-Fate Change.

Herein, we provide a brief overview of complex systems theory approaches to investigate the genomic mechanism of cell-fate changes. Cell trajectories across the epigenetic landscape, whether in development, environmental responses, or disease progression, are controlled by extensively coordinated genome-wide gene expression changes. The elucidation of the mechanisms underlying these coherent expression changes is of fundamental importance in cell biology and for paving the road to new therapeutic approaches. In previous studies, we pointed at dynamic criticality as a plausible characteristic of genome-wide transition dynamics guiding cell fate. Whole-genome expression develops an engine-like organization (genome engine) in order to establish an autonomous dynamical system, capable of both homeostasis and transition behaviors. A critical set of genes behaves as a critical point (CP) that serves as the organizing center of cell-fate change. When the system is pushed away from homeostasis, the state change that occurs at the CP makes local perturbation spread over the genome, demonstrating self-organized critical (SOC) control of genome expression. Oscillating-Mode genes (which normally keep genome expression on pace with microenvironment fluctuations), when in the presence of an effective perturbative stimulus, drive the dynamics of synchronization, and thus guide the cell-fate transition.

Genetic Control of Gonadal Sex Determination and Development.

Sex determination is the process by which the bipotential gonads develop as either testes or ovaries. With two distinct potential outcomes, the gonadal primordium offers a unique model for the study of cell fate specification and how distinct cell populations diverge from multipotent progenitors. This review focuses on recent advances in our understanding of the genetic programs and epigenetic mechanisms that regulate gonadal sex determination and the regulation of cell fate commitment in the bipotential gonads. We rely primarily on mouse data to illuminate the complex and dynamic genetic programs controlling cell fate decision and sex-specific cell differentiation during gonadal formation and gonadal sex determination.

Asymmetric expression of maternal mRNA governs first cell-fate decision.

The goal of preimplantation development is to establish the fates of the embryonic and extra-embryonic cells. However, when and how cell fates are determined during early mammalian embryonic development remains unclear. We report that the high mobility group (HMG) protein family member HMGA1 was distributed differentially in mouse two-cell blastomeres. Knockdown of Hmga1 expression in one of the two cells reduced the number of cells contributing to the inner cell mass (ICM), suggesting that differential distribution of HMGA1 in the blastomeres in two-cell mouse embryos affected the selection of embryonic cell lineages. Mechanistically, HMGA1 promotes the expression of the ICM-specific gene Sox2. The results of this study show that mouse embryos demonstrate heterogeneity as early as the two-cell stage, and that these differences are related to cell-fate differentiation in early mouse embryos.

Multistability and transitions between spatiotemporal patterns through versatile Notch-Hes signaling.

The Delta-Notch-Hes signaling pathway is involved in various developmental processes ranging from the formation of somites to the dynamic fine-grained patterns of cell types in developing or regenerating tissues. Such broad patterning capabilities rely in part on the versatile and tunable dynamics of the Notch-Hes feedback circuit eliciting both pulsatile and switching behaviors. This raises the theoretical issue of which specific spatiotemporal features emerge from lateral inhibition between cells that can display and transit between monostable, oscillatory and bistable regimes. To address this issue, I consider a discrete cell lattice model where intracellular dynamics is described by a phase-like variable and displays a typical cross-shaped phase diagram. Model analysis determines how the existence and stability of many spatially inhomogeneous and temporally synchronized states depends on key intracellular and intercellular parameters. It reveals a parameter-dependent multistability between those diverse spatiotemporal patterns, giving rise to tunable and robust developmental transition scenarios ensuring defect-free spatial patterns. Such broad repertoire and multistability of spatiotemporal patterns is corroborated with regulatory network modeling of the Delta-Notch-Hes pathway.

Arid3a is essential to execution of the first cell fate decision via direct embryonic and extraembryonic transcriptional regulation.

Despite their origin from the inner cell mass, embryonic stem (ES) cells undergo differentiation to the trophectoderm (TE) lineage by repression of the ES cell master regulator Oct4 or activation of the TE master regulator Caudal-type homeobox 2 (Cdx2). In contrast to the in-depth studies of ES cell self-renewal and pluripotency, few TE-specific regulators have been identified, thereby limiting our understanding of mechanisms underlying the first cell fate decision. Here we show that up-regulation and nuclear entry of AT-rich interactive domain 3a (Arid3a) drives TE-like transcriptional programs in ES cells, maintains trophoblast stem (TS) cell self-renewal, and promotes further trophoblastic differentiation both upstream and independent of Cdx2. Accordingly, Arid3a(-/-) mouse post-implantation placental development is severely impaired, resulting in early embryonic death. We provide evidence that Arid3a directly activates TE-specific and trophoblast lineage-specific genes while directly repressing pluripotency genes via differential regulation of epigenetic acetylation or deacetylation. Our results identify Arid3a as a critical regulator of TE and placental development through execution of the commitment and differentiation phases of the first cell fate decision.

Coordination of the AMPK, Akt, mTOR, and p53 Pathways under Glucose Starvation.

Glucose is a direct energy source for eukaryotic cells, and its deficiency elicits complex stress responses and diverse cellular outcomes. Although several signaling pathways involved have been identified, how they coordinately dictate the cell fate remains obscure. We propose a minimal network model for the cellular response to glucose restriction, characterizing the glucose uptake and signaling of the AMPK, Akt, mTOR, and p53 pathways. We demonstrate that in the presence of sufficient growth factors and amino acids, cells may undergo proliferation, senescence, or apoptosis, depending on the extracellular glucose level. AMPK is first activated upon glucose limitation, activating p53 to induce cell-cycle arrest; possibly, cells resume proliferation after timely glucose restoration. For long-term energy stress, cell senescence is maintained by low/intermediate levels of p53 and persistent activation of mTOR and Akt, or cells commit apoptosis when the proteins undergo biphasic dynamics, e.g., p53 switches from intermediate levels to high levels while mTOR and Akt become inactivated in the later phase. The biphasic dynamics of p53 are associated with flipping of two bistable switches. Appropriate mTOR levels are required for optimal cell-fate decision. This work suggests that senescence and apoptosis occur sequentially in glucose-depleted cells, and a theoretical framework is provided for exploring the cellular response to energy stress.