To explore the distinct genotypic and phenotypic states of melanoma tumors we applied single-cell RNA-seq to 4,645 single cells isolated from 19 patients, profiling malignant, immune, stromal and endothelial cells. Malignant cells within the same tumor displayed transcriptional heterogeneity associated with the cell cycle, spatial context, and a drug resistance program. In particular, all tumors harbored malignant cells from two distinct transcriptional cell states, such that “MITF-high” tumors also contained “AXL-high” tumor cells. Single-cell analyses suggested distinct tumor micro-environmental patterns, including cell-to-cell interactions. Analysis of tumor-infiltrating T cells revealed exhaustion programs, their connection to T cell activation and to clonal expansion, and their variability across patients. Overall, we begin to unravel the cellular ecosystem of tumors and how single cell genomics offers insights with implications for both targeted and immune therapies.

Microsystems create new opportunities for the spatial and temporal control of cell growth and stimuli by combining surfaces that mimic complex biochemistries and geometries of the extracellular matrix with microfluidic channels that regulate transport of fluids and soluble factors. Further integration with bioanalytic microsystems results in multifunctional platforms for basic biological insights into cells and tissues, as well as for cell-based sensors with biochemical, biomedical and environmental functions. Highly integrated microdevices show great promise for basic biomedical and pharmaceutical research, and robust and portable point-of-care devices could be used in clinical settings, in both the developed and the developing world.

In microorganisms, noise in gene expression gives rise to cell-to-cell variability in protein concentrations1–7. In mammalian cells, protein levels also vary8–10 and individual cells differ widely in responsiveness to uniform physiological stimuli11–15. In the case of apoptosis mediated by TRAIL (TNF related apoptosis-inducing ligand) it is common for some cells in a clonal population to die while others survive – a striking divergence in cell fate. Among cells that die, the time between TRAIL exposure and caspase activation is highly variable. Here we image sister cells expressing reporters of caspase activation and mitochondrial outer membrane permeabilisation (MOMP) following exposure to TRAIL. We show that naturally occurring differences in the levels or states of proteins regulating receptor-mediated apoptosis are the primary causes of cell-to-cell variability in the timing and probability of death. Protein state is transmitted from mother to daughter, giving rise to transient heritability in fate, but protein synthesis promotes rapid divergence so that sister cells soon become no more similar to each other than pairs of cells chosen at random. Our results have implications for understanding “fractional killing” of tumor cells following exposure to chemotherapy, and for variability in mammalian signal transduction in general.

SUMMARY Immune checkpoint inhibitors (ICIs) produce durable responses in some melanoma patients, but many patients derive no clinical benefit, and the molecular underpinnings of such resistance remain elusive. Here, we leveraged single-cell RNA sequencing (scRNA-seq) from 33 melanoma tumors and computational analyses to interrogate malignant cell states that promote immune evasion. We identified a resistance program expressed by malignant cells that is associated with T cell exclusion and immune evasion. The program is expressed prior to immunotherapy, characterizes cold niches in situ, and predicts clinical responses to anti-PD-1 therapy in an independent cohort of 112 melanoma patients. CDK4/6-inhibition represses this program in individual malignant cells, induces senescence, and reduces melanoma tumor outgrowth in mouse models in vivo when given in combination with immunotherapy. Our study provides a high-resolution landscape of ICI-resistant cell states, identifies clinically predictive signatures, and suggests new therapeutic strategies to overcome immunotherapy resistance.

Drug sensitivity and resistance are conventionally quantified by IC50 or Emax values, but these metrics are highly sensitive to the number of divisions taking place over the course of a response assay. The dependency of IC50 and Emax on division rate creates artefactual correlations between genotype and drug sensitivity while obscuring valuable biological insights and interfering with biomarker discovery. We derive alternative drug response metrics that are insensitive to division number. These are based on estimating the magnitude of drug-induced growth rate inhibition (GR) using endpoint or time-course assays. We show that GR50 and GRmax are superior to conventional metrics for assessing the effects of drugs in dividing cells. Moreover, adopting GR metrics requires only modest changes in experimental protocols. We expect GR metrics to improve the study of cell signaling and growth using drugs, discovery of drug response biomarkers, and identification of drugs effective on specific patient-derived tumor cells.

The mitotic checkpoint protein hsMad2 is required to arrest cells in mitosis when chromosomes are unattached to the mitotic spindle. The presence of a single, lagging chromosome is sufficient to activate the checkpoint, producing a delay at the metaphase-anaphase transition until the last spindle attachment is made. Complete loss of the mitotic checkpoint results in embryonic lethality owing to chromosome mis-segregation in various organisms. Whether partial loss of checkpoint control leads to more subtle rates of chromosome instability compatible with cell viability remains unknown. Here we report that deletion of one MAD2 allele results in a defective mitotic checkpoint in both human cancer cells and murine primary embryonic fibroblasts. Checkpoint-defective cells show premature sister-chromatid separation in the presence of spindle inhibitors and an elevated rate of chromosome mis-segregation events in the absence of these agents. Furthermore, Mad2+/- mice develop lung tumours at high rates after long latencies, implicating defects in the mitotic checkpoint in tumorigenesis.

SUMMARY Crosstalk and complexity within signaling pathways, and their perturbation by oncogenes, limits component-by-component approaches to understanding human disease. Network analysis of how normal and oncogenic signaling can be re-wired by drugs may provide opportunities to target tumors with high specificity and efficacy. Using targeted inhibition of oncogenic signaling pathways combined with DNA damaging chemotherapy, we report that time-staggered EGFR inhibition, but not simultaneous co-administration, dramatically sensitizes a subset of triple-negative breast cancer cells to genotoxic drugs. Systems-level analysis—using high-density time-dependent measurements of signaling networks, gene expression profiles, and cell phenotypic responses in combination with mathematical modeling— revealed an approach for altering the intrinsic state of the cell through dynamic re-wiring of oncogenic signaling pathways. This process converts these cells to a less tumorigenic state that is more susceptible to DNA damage-induced cell death by re-activation of an extrinsic apoptotic pathway whose function is suppressed in the oncogene-addicted state.

No abstract