The origin of lineage correlations among single cells and the extent of heterogeneity in their intermitotic times (IMT) and apoptosis times (AT) remain incompletely understood. Here we developed single cell lineage-tracking experiments and computational algorithms to uncover correlations and heterogeneity in the IMT and AT of a colon cancer cell line before and during cisplatin treatment. These correlations could not be explained using simple protein production/degradation models. Sister cell fates were similar regardless of whether they divided before or after cisplatin administration and did not arise from proximity-related factors, suggesting fate determination early in a cell's lifetime. Based on these findings, we developed a theoretical model explaining how the observed correlation structure can arise from oscillatory mechanisms underlying cell fate control. Our model recapitulated the data only with very specific oscillation periods that fit measured circadian rhythms, thereby suggesting an important role of the circadian clock in controlling cellular fates.
Members of the KDM5 histone H3 lysine 4 demethylase family are associated with therapeutic resistance, including endocrine resistance in breast cancer, but the underlying mechanism is poorly defined. Here we show that genetic deletion of KDM5A/B or inhibition of KDM5 activity increases sensitivity to anti-estrogens by modulating estrogen receptor (ER) signaling and by decreasing cellular transcriptomic heterogeneity. Higher KDM5B expression levels are associated with higher transcriptomic heterogeneity and poor prognosis in ER+ breast tumors. Single-cell RNA sequencing, cellular barcoding, and mathematical modeling demonstrate that endocrine resistance is due to selection for pre-existing genetically distinct cells, while KDM5 inhibitor resistance is acquired. Our findings highlight the importance of cellular phenotypic heterogeneity in therapeutic resistance and identify KDM5A/B as key regulators of this process.
Necroptosis is a lytic programmed cell death mediated by the RIPK1-RIPK3-MLKL pathway. The loss of Receptor-interacting serine/threonine-protein kinase 3 (RIPK3) expression and necroptotic potential have been previously reported in several cancer cell lines; however, the extent of this loss across cancer types, as well as its mutational drivers, were unknown. Here, we show that RIPK3 expression loss occurs progressively during tumor growth both in patient tumor biopsies and tumor xenograft models. Using a cell-based necroptosis sensitivity screen of 941 cancer cell lines, we find that escape from necroptosis is prevalent across cancer types, with an incidence rate of 83%. Genome-wide bioinformatics analysis of this differential necroptosis sensitivity data in the context of differential gene expression and mutation data across the cell lines identified various factors that correlate with resistance to necroptosis and loss of RIPK3 expression, including oncogenes BRAF and AXL. Inhibition of these oncogenes can rescue the RIPK3 expression loss and regain of necroptosis sensitivity. This genome-wide analysis also identifies that the loss of RIPK3 expression is the primary factor correlating with escape from necroptosis. Thus, we conclude that necroptosis resistance of cancer cells is common and is oncogene driven, suggesting that escape from necroptosis could be a potential hallmark of cancer, similar to escape from apoptosis.
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.
Dissecting cellular differentiation hierarchies in the mammary gland is a prerequisite for understanding both normal development and malignant transformation during tumorigenesis and tumor cell-of-origin. To achieve these goals, several recent papers utilized single cell RNA-seq and lineage tracing to improve our understanding of the composition of the mammary epithelium at different developmental stages.
Immunotherapy has emerged as a major therapeutic modality in oncology. Currently, however, the majority of patients with cancer do not derive benefit from these treatments. Vascular abnormalities are a hallmark of most solid tumours and facilitate immune evasion. These abnormalities stem from elevated levels of proangiogenic factors, such as VEGF and angiopoietin 2 (ANG2); judicious use of drugs targeting these molecules can improve therapeutic responsiveness, partially owing to normalization of the abnormal tumour vasculature that can, in turn, increase the infiltration of immune effector cells into tumours and convert the intrinsically immunosuppressive tumour microenvironment (TME) to an immunosupportive one. Immunotherapy relies on the accumulation and activity of immune effector cells within the TME, and immune responses and vascular normalization seem to be reciprocally regulated. Thus, combining antiangiogenic therapies and immunotherapies might increase the effectiveness of immunotherapy and diminish the risk of immune-related adverse effects. In this Perspective, we outline the roles of VEGF and ANG2 in tumour immune evasion and progression, and discuss the evidence indicating that antiangiogenic agents can normalize the TME. We also suggest ways that antiangiogenic agents can be combined with immune-checkpoint inhibitors to potentially improve patient outcomes, and highlight avenues of future research.
The architecture of normal and diseased tissues strongly influences the development and progression of disease as well as responsiveness and resistance to therapy. We describe a tissue-based cyclic immunofluorescence (t-CyCIF) method for highly multiplexed immuno-fluorescence imaging of formalin-fixed, paraffin-embedded (FFPE) specimens mounted on glass slides, the most widely used specimens for histopathological diagnosis of cancer and other diseases. t-CyCIF generates up to 60-plex images using an iterative process (a cycle) in which conventional low-plex fluorescence images are repeatedly collected from the same sample and then assembled into a high-dimensional representation. t-CyCIF requires no specialized instruments or reagents and is compatible with super-resolution imaging; we demonstrate its application to quantifying signal transduction cascades, tumor antigens and immune markers in diverse tissues and tumors. The simplicity and adaptability of t-CyCIF makes it an effective method for pre-clinical and clinical research and a natural complement to single-cell genomics.
In injured tissues, regeneration is often associated with cell fate plasticity, in that cells deviate from their normal lineage paths. It is becoming increasingly clear that this plasticity often creates alternative strategies to restore damaged or lost cells. Alternatively, cell fate plasticity is also part and parcel of pathologic tissue transformations that accompany disease. In this Perspective, we summarize a few illustrative examples of physiologic and aberrant cellular plasticity. Then, we speculate on how one could enhance endogenous plasticity to promote regeneration and reverse pathologic plasticity, perhaps inspiring interest in a new class of therapies targeting cell fate modulation.
Background: Craniopharyngiomas are neoplasms of the sellar/parasellar region that are classified into adamantinomatous craniopharyngioma (ACP) and papillary craniopharyngioma (PCP) subtypes. Surgical resection of craniopharyngiomas is challenging, and recurrence is common, frequently leading to profound morbidity. BRAF V600E mutations render PCP susceptible to BRAF/MEK inhibitors, but effective targeted therapies are needed for ACP. We explored the feasibility of targeting the programmed cell death protein 1/programmed death-ligand 1 (PD-1/PD-L1) immune checkpoint pathway in ACP and PCP. Methods: We mapped and quantified PD-L1 and PD-1 expression in ACP and PCP resections using immunohistochemistry, immunofluorescence, and RNA in situ hybridization. We used tissue-based cyclic immunofluorescence to map the spatial distribution of immune cells and characterize cell cycle and signaling pathways in ACP tumor cells which intrinsically express PD-1. Results: All ACP (15 ± 14% of cells, n = 23, average ± SD) and PCP (35 ± 22% of cells, n = 18) resections expressed PD-L1. In ACP, PD-L1 was predominantly expressed by tumor cells comprising the cyst lining. In PCP, PD-L1 was highly expressed by tumor cells surrounding the stromal fibrovascular cores. ACP also exhibited tumor cell-intrinsic PD-1 expression in whorled epithelial cells with nuclear-localized beta-catenin. These cells exhibited evidence of elevated mammalian target of rapamycin (mTOR) and mitogen-activated protein kinase (MAPK) signaling. Profiling of immune populations in ACP and PCP showed a modest density of CD8+ T cells. Conclusions: ACP exhibit PD-L1 expression in the tumor cyst lining and intrinsic PD-1 expression in cells proposed to comprise an oncogenic stem-like population. In PCP, proliferative tumor cells express PD-L1 in a continuous band at the stromal-epithelial interface. Targeting PD-L1 and/or PD-1 in both subtypes of craniopharyngioma might therefore be an effective therapeutic strategy.
Triple-negative breast cancer (TNBC) is an aggressive subtype characterized by extensive intratumoral heterogeneity. To investigate the underlying biology, we conducted single-cell RNA-sequencing (scRNA-seq) of >1500 cells from six primary TNBC. Here, we show that intercellular heterogeneity of gene expression programs within each tumor is variable and largely correlates with clonality of inferred genomic copy number changes, suggesting that genotype drives the gene expression phenotype of individual subpopulations. Clustering of gene expression profiles identified distinct subgroups of malignant cells shared by multiple tumors, including a single subpopulation associated with multiple signatures of treatment resistance and metastasis, and characterized functionally by activation of glycosphingolipid metabolism and associated innate immunity pathways. A novel signature defining this subpopulation predicts long-term outcomes for TNBC patients in a large cohort. Collectively, this analysis reveals the functional heterogeneity and its association with genomic evolution in TNBC, and uncovers unanticipated biological principles dictating poor outcomes in this disease.
Cells respond to DNA damage by activating complex signaling networks that decide cell fate, promoting not only DNA damage repair and survival but also cell death. We have developed a multiscale computational model that quantitatively links chemotherapy-induced DNA damage response signaling to cell fate. The computational model was trained and calibrated on extensive data from U2OS osteosarcoma cells, including the cell cycle distribution of the initial cell population, signaling data measured by Western blotting, and cell fate data in response to chemotherapy treatment measured by time-lapse microscopy. The resulting mechanistic model predicted the cellular responses to chemotherapy alone and in combination with targeted inhibitors of the DNA damage response pathway, which we confirmed experimentally. Computational models such as the one presented here can be used to understand the molecular basis that defines the complex interplay between cell survival and cell death and to rationally identify chemotherapy-potentiating drug combinations.
Cancer cells subvert host immune surveillance by altering immune checkpoint (IC) proteins. Some Epstein-Barr virus (EBV)-associated tumors have higher Programmed Cell Death Ligand, PD-L1 expression. However, it is not known how EBV alters ICs in the context of its preferred host, the B lymphocyte and in derived lymphomas. Here, we found that latency III-expressing Burkitt lymphoma (BL), diffuse large B-cell lymphomas (DLBCL) or their EBNA2-transfected derivatives express high PD-L1. In a DLBCL model, EBNA2 but not LMP1 is sufficient to induce PD-L1. Latency III-expressing DLBCL biopsies showed high levels of PD-L1. The PD-L1 targeting oncosuppressor microRNA miR-34a was downregulated in EBNA2-transfected lymphoma cells. We identified early B-cell factor 1 (EBF1) as a repressor of miR-34a transcription. Short hairpin RNA (shRNA)-mediated knockdown of EBF1 was sufficient to induce miR-34a transcription, which in turn reduced PD-L1. MiR-34a reconstitution in EBNA2-transfected DLBCL reduced PD-L1 expression and increased its immunogenicity in mixed lymphocyte reactions (MLR) and in three-dimensional biomimetic microfluidic chips. Given the importance of PD-L1 inhibition in immunotherapy and miR-34a dysregulation in cancers, our findings may have important implications for combinatorial immunotherapy, which include IC inhibiting antibodies and miR-34a, for EBV-associated cancers.
Cancer cell survival is dependent on oxidative-stress defenses against reactive oxygen species (ROS) that accumulate during tumorigenesis. Here, we show a non-canonical oxidative-stress defense mechanism through TRPA1, a neuronal redox-sensing Ca2+-influx channel. In TRPA1-enriched breast and lung cancer spheroids, TRPA1 is critical for survival of inner cells that exhibit ROS accumulation. Moreover, TRPA1 promotes resistance to ROS-producing chemotherapies, and TRPA1 inhibition suppresses xenograft tumor growth and enhances chemosensitivity. TRPA1 does not affect redox status but upregulates Ca2+-dependent anti-apoptotic pathways. NRF2, an oxidant-defense transcription factor, directly controls TRPA1 expression, thus providing an orthogonal mechanism for protection against oxidative stress together with canonical ROS-neutralizing mechanisms. These findings reveal an oxidative-stress defense program involving TRPA1 that could be exploited for targeted cancer therapies.
Resistance to chemotherapy plays a significant role in cancer mortality. To identify genetic units affecting sensitivity to cytarabine, the mainstay of treatment for acute myeloid leukemia (AML), we developed a comprehensive and integrated genome-wide platform based on a dual protein-coding and non-coding integrated CRISPRa screening (DICaS). Putative resistance genes were initially identified using pharmacogenetic data from 760 human pan-cancer cell lines. Subsequently, genome scale functional characterization of both coding and long non-coding RNA (lncRNA) genes by CRISPR activation was performed. For lncRNA functional assessment, we developed a CRISPR activation of lncRNA (CaLR) strategy, targeting 14,701 lncRNA genes. Computational and functional analysis identified novel cell-cycle, survival/apoptosis, and cancer signaling genes. Furthermore, transcriptional activation of the GAS6-AS2 lncRNA, identified in our analysis, leads to hyperactivation of the GAS6/TAM pathway, a resistance mechanism in multiple cancers including AML. Thus, DICaS represents a novel and powerful approach to identify integrated coding and non-coding pathways of therapeutic relevance.
The PI3K-AKT kinase signaling pathway is frequently deregulated in human cancers, particularly breast cancer, where amplification and somatic mutations of PIK3CA occur with high frequency in patients. Numerous small-molecule inhibitors targeting both PI3K and AKT are under clinical evaluation, but dose-limiting toxicities and the emergence of resistance limit therapeutic efficacy. Various resistance mechanisms to PI3K inhibitors have been identified, including de novo mutations, feedback activation of AKT, or cross-talk pathways. We found a previously unknown resistance mechanism to PI3K pathway inhibition that results in AKT rebound activation. In a subset of triple-negative breast cancer cell lines, treatment with a PI3K inhibitor or depletion of PIK3CA expression ultimately promoted AKT reactivation in a manner dependent on the E3 ubiquitin ligase Skp2, the kinases IGF-1R (insulin-like growth factor 1 receptor) and PDK-1 (phosphoinositide-dependent kinase-1), and the cell growth and metabolism-regulating complex mTORC2 (mechanistic target of rapamycin complex 2), but was independent of PI3K activity or PIP3 production. Resistance to PI3K inhibitors correlated with the increased abundance of Skp2, ubiquitylation of AKT, cell proliferation in culture, and xenograft tumor growth in mice. These findings reveal a ubiquitin signaling feedback mechanism by which PI3K inhibitor resistance may emerge in aggressive breast cancer cells.
Triple negative breast cancers (TNBC) are genetically characterized by aberrations in TP53 and a low rate of activating point mutations in common oncogenes, rendering it challenging in applying targeted therapies. We performed whole exome sequencing (WES) and RNAseq to identify somatic genetic alterations in mouse models of TNBCs driven by loss of Trp53 alone or in combination with Brca1. Amplifications or translocations that resulted in elevated oncoprotein expressions or oncoprotein-containing fusions, respectively, as well as frame-shift mutations of tumor suppressors were identified in approximately 50% of the tumors evaluated. While the spectrum of sporadic genetic alterations was diverse, the majority had in common the ability to activate the MAPK/PI3K pathways. Importantly, we demonstrated that approved or experimental drugs efficiently induce tumor regression specifically in tumors harboring somatic aberrations of the drug target. Our study suggests that the combination of WES and RNAseq on human TNBC will lead to the identification of actionable therapeutic targets for precision medicine guided TNBC treatment.
The precursor homocysteine is metabolized either through the methionine cycle to produce methionine or through the transsulfuration pathway to synthesize cysteine. Alternatively, cysteine can be obtained through uptake of its oxidized form, cystine. Many cancer cells exhibit methionine dependency such that their proliferation is impaired in growth media in which methionine is replaced by homocysteine. We showed that oncogenic PIK3CA and decreased expression of SLC7A11, a gene that encodes a cystine transporter also known as xCT, correlated with increased methionine dependency in breast cancer cells. Oncogenic PIK3CA was sufficient to confer methionine dependency to mammary epithelial cells, partly by decreasing cystine uptake through the transcriptional and posttranslational inhibition of xCT. Manipulation of xCT activity altered the proliferation of breast cancer cells in methionine-deficient, homocysteine-containing media, suggesting that it functionally contributed to methionine dependency. We propose that concurrent with decreased cystine uptake through xCT, PIK3CA mutant cells use homocysteine through the transsulfuration pathway to synthesize cysteine. Consequently, less homocysteine is available to produce methionine, contributing to methionine dependency. These results indicate that oncogenic PIK3CA alters methionine and cysteine utilization, partly by inhibiting xCT to contribute to the methionine dependency phenotype in breast cancer cells.
Combination cancer therapies aim to improve the probability and magnitude of therapeutic responses and reduce the likelihood of acquired resistance in an individual patient. However, drugs are tested in clinical trials on genetically diverse patient populations. We show here that patient-to-patient variability and independent drug action are sufficient to explain the superiority of many FDA-approved drug combinations in the absence of drug synergy or additivity. This is also true for combinations tested in patient-derived tumor xenografts. In a combination exhibiting independent drug action, each patient benefits solely from the drug to which his or her tumor is most sensitive, with no added benefit from other drugs. Even when drug combinations exhibit additivity or synergy in pre-clinical models, patient-to-patient variability and low cross-resistance make independent action the dominant mechanism in clinical populations. This insight represents a different way to interpret trial data and a different way to design combination therapies.
The tumor-suppressing transcription factor p53 is highly conserved at the protein level and plays a key role in the DNA damage response. One important aspect of p53 regulation is its dynamics in response to DNA damage, which include oscillations. Here, we observe that, while the qualitative oscillatory nature of p53 dynamics is conserved across cell lines derived from human, monkey, dog, mouse, and rat, the oscillation period is variable. Specifically, rodent cells exhibit rapid p53 oscillations, whereas dog, monkey, and human cells show slower oscillations. Computational modeling and experiments identify stronger negative feedback between p53 and MDM2 as the driver of faster oscillations in rodents, suggesting that the period of oscillation is a network-level property. In total, our study shows that despite highly conserved signaling, the quantitative features of p53 oscillations can diverge across evolution. We caution that strong amino acid conservation of proteins and transcriptional network similarity do not necessarily imply conservation of time dynamics.