MRE11 liberates cGAS from nucleosome sequestration during tumorigenesis.

Oncogene-induced replication stress generates endogenous DNA damage that activates cGAS-STING-mediated signalling and tumour suppression1-3. However, the precise mechanism of cGAS activation by endogenous DNA damage remains enigmatic, particularly given that high-affinity histone acidic patch (AP) binding constitutively inhibits cGAS by sterically hindering its activation by double-stranded DNA (dsDNA)4-10. Here we report that the DNA double-strand break sensor MRE11 suppresses mammary tumorigenesis through a pivotal role in regulating cGAS activation. We demonstrate that binding of the MRE11-RAD50-NBN complex to nucleosome fragments is necessary to displace cGAS from acidic-patch-mediated sequestration, which enables its mobilization and activation by dsDNA. MRE11 is therefore essential for cGAS activation in response to oncogenic stress, cytosolic dsDNA and ionizing radiation. Furthermore, MRE11-dependent cGAS activation promotes ZBP1-RIPK3-MLKL-mediated necroptosis, which is essential to suppress oncogenic proliferation and breast tumorigenesis. Notably, downregulation of ZBP1 in human triple-negative breast cancer is associated with increased genome instability, immune suppression and poor patient prognosis. These findings establish MRE11 as a crucial mediator that links DNA damage and cGAS activation, resulting in tumour suppression through ZBP1-dependent necroptosis.

Encoding and decoding cellular information through signaling dynamics.

A growing number of studies are revealing that cells can send and receive information by controlling the temporal behavior (dynamics) of their signaling molecules. In this Review, we discuss what is known about the dynamics of various signaling networks and their role in controlling cellular responses. We identify general principles that are emerging in the field, focusing specifically on how the identity and quantity of a stimulus is encoded in temporal patterns, how signaling dynamics influence cellular outcomes, and how specific dynamical patterns are both shaped and interpreted by the structure of molecular networks. We conclude by discussing potential functional roles for transmitting cellular information through the dynamics of signaling molecules and possible applications for the treatment of disease.

Cell cycle plasticity underlies fractional resistance to palbociclib in ER+/HER2- breast tumor cells.

The CDK4/6 inhibitor palbociclib blocks cell cycle progression in Estrogen receptor-positive, human epidermal growth factor 2 receptor-negative (ER+/HER2-) breast tumor cells. Despite the drug's success in improving patient outcomes, a small percentage of tumor cells continues to divide in the presence of palbociclib-a phenomenon we refer to as fractional resistance. It is critical to understand the cellular mechanisms underlying fractional resistance because the precise percentage of resistant cells in patient tissue is a strong predictor of clinical outcomes. Here, we hypothesize that fractional resistance arises from cell-to-cell differences in core cell cycle regulators that allow a subset of cells to escape CDK4/6 inhibitor therapy. We used multiplex, single-cell imaging to identify fractionally resistant cells in both cultured and primary breast tumor samples resected from patients. Resistant cells showed premature accumulation of multiple G1 regulators including E2F1, retinoblastoma protein, and CDK2, as well as enhanced sensitivity to pharmacological inhibition of CDK2 activity. Using trajectory inference approaches, we show how plasticity among cell cycle regulators gives rise to alternate cell cycle "paths" that allow individual tumor cells to escape palbociclib treatment. Understanding drivers of cell cycle plasticity, and how to eliminate resistant cell cycle paths, could lead to improved cancer therapies targeting fractionally resistant cells to improve patient outcomes.

Cellograph: a semi-supervised approach to analyzing multi-condition single-cell RNA-sequencing data using graph neural networks.

With the growing number of single-cell datasets collected under more complex experimental conditions, there is an opportunity to leverage single-cell variability to reveal deeper insights into how cells respond to perturbations. Many existing approaches rely on discretizing the data into clusters for differential gene expression (DGE), effectively ironing out any information unveiled by the single-cell variability across cell-types. In addition, DGE often assumes a statistical distribution that, if erroneous, can lead to false positive differentially expressed genes. Here, we present Cellograph: a semi-supervised framework that uses graph neural networks to quantify the effects of perturbations at single-cell granularity. Cellograph not only measures how prototypical cells are of each condition but also learns a latent space that is amenable to interpretable data visualization and clustering. The learned gene weight matrix from training reveals pertinent genes driving the differences between conditions. We demonstrate the utility of our approach on publicly-available datasets including cancer drug therapy, stem cell reprogramming, and organoid differentiation. Cellograph outperforms existing methods for quantifying the effects of experimental perturbations and offers a novel framework to analyze single-cell data using deep learning.

The consequences of differential origin licensing dynamics in distinct chromatin environments.

Eukaryotic chromosomes contain regions of varying accessibility, yet DNA replication factors must access all regions. The first replication step is loading MCM complexes to license replication origins during the G1 cell cycle phase. It is not yet known how mammalian MCM complexes are adequately distributed to both accessible euchromatin regions and less accessible heterochromatin regions. To address this question, we combined time-lapse live-cell imaging with immunofluorescence imaging of single human cells to quantify the relative rates of MCM loading in euchromatin and heterochromatin throughout G1. We report here that MCM loading in euchromatin is faster than that in heterochromatin in early G1, but surprisingly, heterochromatin loading accelerates relative to euchromatin loading in middle and late G1. This differential acceleration allows both chromatin types to begin S phase with similar concentrations of loaded MCM. The different loading dynamics require ORCA-dependent differences in origin recognition complex distribution. A consequence of heterochromatin licensing dynamics is that cells experiencing a truncated G1 phase from premature cyclin E expression enter S phase with underlicensed heterochromatin, and DNA damage accumulates preferentially in heterochromatin in the subsequent S/G2 phase. Thus, G1 length is critical for sufficient MCM loading, particularly in heterochromatin, to ensure complete genome duplication and to maintain genome stability.

In this study the authors have, for the first time, quantified DNA replication origin licensing dynamics and distribution in single cells at subnuclear resolution. The cell cycle and DNA replication fields have long appreciated that origin licensing is both absolutely essential for replication and that licensing is strictly confined to G1 phase. The biochemical process of origin licensing- which is the DNA loading of MCM complexes- is known in considerable detail. What has never been explored in any system, is the dynamics of origin licensing itself. Here the authors define the dynamics of human MCM loading at different times within G1 in both euchromatin and heterochromatin, and explore the consequences of those dynamics for genome stability.

The molecular architecture of cell cycle arrest.

The cellular decision governing the transition between proliferative and arrested states is crucial to the development and function of every tissue. While the molecular mechanisms that regulate the proliferative cell cycle are well established, we know comparatively little about what happens to cells as they diverge into cell cycle arrest. We performed hyperplexed imaging of 47 cell cycle effectors to obtain a map of the molecular architecture that governs cell cycle exit and progression into reversible ("quiescent") and irreversible ("senescent") arrest states. Using this map, we found multiple points of divergence from the proliferative cell cycle; identified stress-specific states of arrest; and resolved the molecular mechanisms governing these fate decisions, which we validated by single-cell, time-lapse imaging. Notably, we found that cells can exit into senescence from either G1 or G2; however, both subpopulations converge onto a single senescent state with a G1-like molecular signature. Cells can escape from this "irreversible" arrest state through the upregulation of G1 cyclins. This map provides a more comprehensive understanding of the overall organization of cell proliferation and arrest.

Sequential replication-coupled destruction at G1/S ensures genome stability.

Timely ubiquitin-mediated protein degradation is fundamental to cell cycle control, but the precise degradation order at each cell cycle phase transition is still unclear. We investigated the degradation order among substrates of a single human E3 ubiquitin ligase, CRL4(Cdt2), which mediates the S-phase degradation of key cell cycle proteins, including Cdt1, PR-Set7, and p21. Our analysis of synchronized cells and asynchronously proliferating live single cells revealed a consistent order of replication-coupled destruction during both S-phase entry and DNA repair; Cdt1 is destroyed first, whereas p21 destruction is always substantially later than that of Cdt1. These differences are attributable to the CRL4(Cdt2) targeting motif known as the PIP degron, which binds DNA-loaded proliferating cell nuclear antigen (PCNA(DNA)) and recruits CRL4(Cdt2). Fusing Cdt1's PIP degron to p21 causes p21 to be destroyed nearly concurrently with Cdt1 rather than consecutively. This accelerated degradation conferred by the Cdt1 PIP degron is accompanied by more effective Cdt2 recruitment by Cdt1 even though p21 has higher affinity for PCNA(DNA). Importantly, cells with artificially accelerated p21 degradation display evidence of stalled replication in mid-S phase and sensitivity to replication arrest. We therefore propose that sequential degradation ensures orderly S-phase progression to avoid replication stress and genome instability.

Inferring the structures of signaling motifs from paired dynamic traces of single cells.

Individual cells show variability in their signaling dynamics that often correlates with phenotypic responses, indicating that cell-to-cell variability is not merely noise but can have functional consequences. Based on this observation, we reasoned that cell-to-cell variability under the same treatment condition could be explained in part by a single signaling motif that maps different upstream signals into a corresponding set of downstream responses. If this assumption holds, then repeated measurements of upstream and downstream signaling dynamics in a population of cells could provide information about the underlying signaling motif for a given pathway, even when no prior knowledge of that motif exists. To test these two hypotheses, we developed a computer algorithm called MISC (Motif Inference from Single Cells) that infers the underlying signaling motif from paired time-series measurements from individual cells. When applied to measurements of transcription factor and reporter gene expression in the yeast stress response, MISC predicted signaling motifs that were consistent with previous mechanistic models of transcription. The ability to detect the underlying mechanism became less certain when a cell's upstream signal was randomly paired with another cell's downstream response, demonstrating how averaging time-series measurements across a population obscures information about the underlying signaling mechanism. In some cases, motif predictions improved as more cells were added to the analysis. These results provide evidence that mechanistic information about cellular signaling networks can be systematically extracted from the dynamical patterns of single cells.

Evidence that the human cell cycle is a series of uncoupled, memoryless phases.

The cell cycle is canonically described as a series of four consecutive phases: G1, S, G2, and M. In single cells, the duration of each phase varies, but the quantitative laws that govern phase durations are not well understood. Using time-lapse microscopy, we found that each phase duration follows an Erlang distribution and is statistically independent from other phases. We challenged this observation by perturbing phase durations through oncogene activation, inhibition of DNA synthesis, reduced temperature, and DNA damage. Despite large changes in durations in cell populations, phase durations remained uncoupled in individual cells. These results suggested that the independence of phase durations may arise from a large number of molecular factors that each exerts a minor influence on the rate of cell cycle progression. We tested this model by experimentally forcing phase coupling through inhibition of cyclin-dependent kinase 2 (CDK2) or overexpression of cyclin D. Our work provides an explanation for the historical observation that phase durations are both inherited and independent and suggests how cell cycle progression may be altered in disease states.

Decoding the insulin signal.

In this issue of Molecular Cell, Kubota et al. (2012) show how different temporal patterns of insulin are decoded by the AKT signaling network, providing both new mechanistic insights and physiological relevance.

Inheritance of OCT4 predetermines fate choice in human embryonic stem cells.

It is well known that clonal cells can make different fate decisions, but it is unclear whether these decisions are determined during, or before, a cell's own lifetime. Here, we engineered an endogenous fluorescent reporter for the pluripotency factor OCT4 to study the timing of differentiation decisions in human embryonic stem cells. By tracking single-cell OCT4 levels over multiple cell cycle generations, we found that the decision to differentiate is largely determined before the differentiation stimulus is presented and can be predicted by a cell's preexisting OCT4 signaling patterns. We further quantified how maternal OCT4 levels were transmitted to, and distributed between, daughter cells. As mother cells underwent division, newly established OCT4 levels in daughter cells rapidly became more predictive of final OCT4 expression status. These results imply that the choice between developmental cell fates can be largely predetermined at the time of cell birth through inheritance of a pluripotency factor.

An integrative transcriptomics approach identifies miR-503 as a candidate master regulator of the estrogen response in MCF-7 breast cancer cells.

Estrogen receptor α (ERα) is an important biomarker of breast cancer severity and a common therapeutic target. In response to estrogen, ERα stimulates a dynamic transcriptional program including both coding and noncoding RNAs. We generate a fine-scale map of expression dynamics by performing a temporal profiling of both messenger RNAs (mRNAs) and microRNAs (miRNAs) in MCF-7 cells (an ER+ model cell line for breast cancer) in response to estrogen stimulation. We identified three primary expression trends-transient, induced, and repressed-that were each enriched for genes with distinct cellular functions. Integrative analysis of mRNA and miRNA temporal expression profiles identified miR-503 as the strongest candidate master regulator of the estrogen response, in part through suppression of ZNF217-an oncogene that is frequently amplified in cancer. We confirmed experimentally that miR-503 directly targets ZNF217 and that overexpression of miR-503 suppresses MCF-7 cell proliferation. Moreover, the levels of ZNF217 and miR-503 are associated with opposite outcomes in breast cancer patient cohorts, with high expression of ZNF217 associated with poor survival and high expression of miR-503 associated with improved survival. Overall, these data indicate that miR-503 acts as a potent estrogen-induced candidate tumor suppressor miRNA that opposes cellular proliferation and has promise as a novel therapeutic for breast cancer. More generally, our work provides a systems-level framework for identifying functional interactions that shape the temporal dynamics of gene expression.

Computational analysis of signaling patterns in single cells.

Signaling proteins are flexible in both form and function. They can bind to multiple molecular partners and integrate diverse types of cellular information. When imaged by time-lapse microscopy, many signaling proteins show complex patterns of activity or localization that vary from cell to cell. This heterogeneity is so prevalent that it has spurred the development of new computational strategies to analyze single-cell signaling patterns. A collective observation from these analyses is that cells appear less heterogeneous when their responses are normalized to, or synchronized with, other single-cell measurements. In many cases, these transformed signaling patterns show distinct dynamical trends that correspond with predictable phenotypic outcomes. When signaling mechanisms are unclear, computational models can suggest putative molecular interactions that are experimentally testable. Thus, computational analysis of single-cell signaling has not only provided new ways to quantify the responses of individual cells, but has helped resolve longstanding questions surrounding many well-studied human signaling proteins including NF-κB, p53, ERK1/2, and CDK2. A number of specific challenges lie ahead for single-cell analysis such as quantifying the contribution of non-cell autonomous signaling as well as the characterization of protein signaling dynamics in vivo.

DELVE: feature selection for preserving biological trajectories in single-cell data.

Single-cell technologies can measure the expression of thousands of molecular features in individual cells undergoing dynamic biological processes. While examining cells along a computationally-ordered pseudotime trajectory can reveal how changes in gene or protein expression impact cell fate, identifying such dynamic features is challenging due to the inherent noise in single-cell data. Here, we present DELVE, an unsupervised feature selection method for identifying a representative subset of molecular features which robustly recapitulate cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effects of confounding sources of variation, and instead models cell states from dynamic gene or protein modules based on core regulatory complexes. Using simulations, single-cell RNA sequencing, and iterative immunofluorescence imaging data in the context of cell cycle and cellular differentiation, we demonstrate how DELVE selects features that better define cell-types and cell-type transitions. DELVE is available as an open-source python package: https://github.com/jranek/delve .

Feature selection for preserving biological trajectories in single-cell data.

Single-cell technologies can readily measure the expression of thousands of molecular features from individual cells undergoing dynamic biological processes, such as cellular differentiation, immune response, and disease progression. While examining cells along a computationally ordered pseudotime offers the potential to study how subtle changes in gene or protein expression impact cell fate decision-making, identifying characteristic features that drive continuous biological processes remains difficult to detect from unenriched and noisy single-cell data. Given that all profiled sources of feature variation contribute to the cell-to-cell distances that define an inferred cellular trajectory, including confounding sources of biological variation (e.g. cell cycle or metabolic state) or noisy and irrelevant features (e.g. measurements with low signal-to-noise ratio) can mask the underlying trajectory of study and hinder inference. Here, we present DELVE (dynamic selection of locally covarying features), an unsupervised feature selection method for identifying a representative subset of dynamically-expressed molecular features that recapitulates cellular trajectories. In contrast to previous work, DELVE uses a bottom-up approach to mitigate the effect of unwanted sources of variation confounding inference, and instead models cell states from dynamic feature modules that constitute core regulatory complexes. Using simulations, single-cell RNA sequencing data, and iterative immunofluorescence imaging data in the context of the cell cycle and cellular differentiation, we demonstrate that DELVE selects features that more accurately characterize cell populations and improve the recovery of cell type transitions. This feature selection framework provides an alternative approach for improving trajectory inference and uncovering co-variation amongst features along a biological trajectory. DELVE is implemented as an open-source python package and is publicly available at: https://github.com/jranek/delve.

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.

Integrating temporal single-cell gene expression modalities for trajectory inference and disease prediction.

BACKGROUND: Current methods for analyzing single-cell datasets have relied primarily on static gene expression measurements to characterize the molecular state of individual cells. However, capturing temporal changes in cell state is crucial for the interpretation of dynamic phenotypes such as the cell cycle, development, or disease progression. RNA velocity infers the direction and speed of transcriptional changes in individual cells, yet it is unclear how these temporal gene expression modalities may be leveraged for predictive modeling of cellular dynamics. RESULTS: Here, we present the first task-oriented benchmarking study that investigates integration of temporal sequencing modalities for dynamic cell state prediction. We benchmark ten integration approaches on ten datasets spanning different biological contexts, sequencing technologies, and species. We find that integrated data more accurately infers biological trajectories and achieves increased performance on classifying cells according to perturbation and disease states. Furthermore, we show that simple concatenation of spliced and unspliced molecules performs consistently well on classification tasks and can be used over more memory intensive and computationally expensive methods. CONCLUSIONS: This work illustrates how integrated temporal gene expression modalities may be leveraged for predicting cellular trajectories and sample-associated perturbation and disease phenotypes. Additionally, this study provides users with practical recommendations for task-specific integration of single-cell gene expression modalities.

The structure of the human cell cycle.

Understanding the organization of the cell cycle has been a longstanding goal in cell biology. We combined time-lapse microscopy, highly multiplexed single-cell imaging of 48 core cell cycle proteins, and manifold learning to render a visualization of the human cell cycle. This data-driven approach revealed the comprehensive "structure" of the cell cycle: a continuum of molecular states that cells occupy as they transition from one cell division to the next, or as they enter or exit cell cycle arrest. Paradoxically, progression deeper into cell cycle arrest was accompanied by increases in proliferative effectors such as CDKs and cyclins, which can drive cell cycle re-entry by overcoming p21 induction. The structure also revealed the molecular trajectories into senescence and the unique combination of molecular features that define this irreversibly arrested state. This approach will enable the comparison of alternative cell cycles during development, in response to environmental perturbation and in disease. A record of this paper's transparent peer review process is included in the supplemental information.