Search for: All records

Abstract 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. 

Abstract 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. 

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. 

Abstract Oncogene-induced replication stress generates endogenous DNA damage that activates cGAS–STING-mediated signalling and tumour suppression^1–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 ofZBP1in 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. 

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. 

Cell cycle gene expression programs fuel proliferation and are universally dysregulated in cancer. The retinoblastoma (RB)-family of proteins, RB1, RBL1/p107, and RBL2/p130, coordinately represses cell cycle gene expression, inhibiting proliferation, and suppressing tumorigenesis. Phosphorylation of RB-family proteins by cyclin-dependent kinases is firmly established. Like phosphorylation, ubiquitination is essential to cell cycle control, and numerous proliferative regulators, tumor suppressors, and oncoproteins are ubiquitinated. However, little is known about the role of ubiquitin signaling in controlling RB-family proteins. A systems genetics analysis of CRISPR/Cas9 screens suggested the potential regulation of the RB-network by cyclin F, a substrate recognition receptor for the SCF family of E3 ligases. We demonstrate that RBL2/p130 is a direct substrate of SCF cyclin F . We map a cyclin F regulatory site to a flexible linker in the p130 pocket domain, and show that this site mediates binding, stability, and ubiquitination. Expression of a mutant version of p130, which cannot be ubiquitinated, severely impaired proliferative capacity and cell cycle progression. Consistently, we observed reduced expression of cell cycle gene transcripts, as well a reduced abundance of cell cycle proteins, analyzed by quantitative, iterative immunofluorescent imaging. These data suggest a key role for SCF cyclin F in the CDK-RB network and raise the possibility that aberrant p130 degradation could dysregulate the cell cycle in human cancers. 

This study demonstrates application of convolutional neural networks (CNNs) for the analysis of a unique image analysis problem in fluorescence microscopy. We employed the U-Net CNN architecture and trained a model to segment nuclear regions in images of a translocating biosensor—which alternates between the nucleus and cytoplasm—without the need for a constant nuclear marker. The modelprovided high-quality segmentation results that allowed us to accurately quantify the extent of cyclin-dependentkinase activity in a population of cells. We envision that the development of this kind of analysis tools will enable biologists to design live-cell fluorescence imaging experiments without the need for providing a constant marker for a subcellular region of interest. As a consequence, they willbe free to increase the number of biosensors measured in single cells or reduce the phototoxicity of cellular imaging.