Senescence-Induced Vascular Remodeling Creates Therapeutic Vulnerabilities in Pancreas Cancer.

KRAS mutant pancreatic ductal adenocarcinoma (PDAC) is characterized by a desmoplastic response that promotes hypovascularity, immunosuppression, and resistance to chemo- and immunotherapies. We show that a combination of MEK and CDK4/6 inhibitors that target KRAS-directed oncogenic signaling can suppress PDAC proliferation through induction of retinoblastoma (RB) protein-mediated senescence. In preclinical mouse models of PDAC, this senescence-inducing therapy produces a senescence-associated secretory phenotype (SASP) that includes pro-angiogenic factors that promote tumor vascularization, which in turn enhances drug delivery and efficacy of cytotoxic gemcitabine chemotherapy. In addition, SASP-mediated endothelial cell activation stimulates the accumulation of CD8+ T cells into otherwise immunologically "cold" tumors, sensitizing tumors to PD-1 checkpoint blockade. Therefore, in PDAC models, therapy-induced senescence can establish emergent susceptibilities to otherwise ineffective chemo- and immunotherapies through SASP-dependent effects on the tumor vasculature and immune system.

p53 Represses the Mevalonate Pathway to Mediate Tumor Suppression.

There are still gaps in our understanding of the complex processes by which p53 suppresses tumorigenesis. Here we describe a novel role for p53 in suppressing the mevalonate pathway, which is responsible for biosynthesis of cholesterol and nonsterol isoprenoids. p53 blocks activation of SREBP-2, the master transcriptional regulator of this pathway, by transcriptionally inducing the ABCA1 cholesterol transporter gene. A mouse model of liver cancer reveals that downregulation of mevalonate pathway gene expression by p53 occurs in premalignant hepatocytes, when p53 is needed to actively suppress tumorigenesis. Furthermore, pharmacological or RNAi inhibition of the mevalonate pathway restricts the development of murine hepatocellular carcinomas driven by p53 loss. Like p53 loss, ablation of ABCA1 promotes murine liver tumorigenesis and is associated with increased SREBP-2 maturation. Our findings demonstrate that repression of the mevalonate pathway is a crucial component of p53-mediated liver tumor suppression and outline the mechanism by which this occurs.

STING inhibits the reactivation of dormant metastasis in lung adenocarcinoma.

Metastasis frequently develops from disseminated cancer cells that remain dormant after the apparently successful treatment of a primary tumour. These cells fluctuate between an immune-evasive quiescent state and a proliferative state liable to immune-mediated elimination1-6. Little is known about the clearing of reawakened metastatic cells and how this process could be therapeutically activated to eliminate residual disease in patients. Here we use models of indolent lung adenocarcinoma metastasis to identify cancer cell-intrinsic determinants of immune reactivity during exit from dormancy. Genetic screens of tumour-intrinsic immune regulators identified the stimulator of interferon genes (STING) pathway as a suppressor of metastatic outbreak. STING activity increases in metastatic progenitors that re-enter the cell cycle and is dampened by hypermethylation of the STING promoter and enhancer in breakthrough metastases or by chromatin repression in cells re-entering dormancy in response to TGFß. STING expression in cancer cells derived from spontaneous metastases suppresses their outgrowth. Systemic treatment of mice with STING agonists eliminates dormant metastasis and prevents spontaneous outbreaks in a T cell- and natural killer cell-dependent manner-these effects require cancer cell STING function. Thus, STING provides a checkpoint against the progression of dormant metastasis and a therapeutically actionable strategy for the prevention of disease relapse.

Ordered and deterministic cancer genome evolution after p53 loss.

Although p53 inactivation promotes genomic instability1 and presents a route to malignancy for more than half of all human cancers2,3, the patterns through which heterogenous TP53 (encoding human p53) mutant genomes emerge and influence tumorigenesis remain poorly understood. Here, in a mouse model of pancreatic ductal adenocarcinoma that reports sporadic p53 loss of heterozygosity before cancer onset, we find that malignant properties enabled by p53 inactivation are acquired through a predictable pattern of genome evolution. Single-cell sequencing and in situ genotyping of cells from the point of p53 inactivation through progression to frank cancer reveal that this deterministic behaviour involves four sequential phases-Trp53 (encoding mouse p53) loss of heterozygosity, accumulation of deletions, genome doubling, and the emergence of gains and amplifications-each associated with specific histological stages across the premalignant and malignant spectrum. Despite rampant heterogeneity, the deletion events that follow p53 inactivation target functionally relevant pathways that can shape genomic evolution and remain fixed as homogenous events in diverse malignant populations. Thus, loss of p53-the 'guardian of the genome'-is not merely a gateway to genetic chaos but, rather, can enable deterministic patterns of genome evolution that may point to new strategies for the treatment of TP53-mutant tumours.

Spatial organization of the mouse genome and its role in recurrent chromosomal translocations.

The extent to which the three-dimensional organization of the genome contributes to chromosomal translocations is an important question in cancer genomics. We generated a high-resolution Hi-C spatial organization map of the G1-arrested mouse pro-B cell genome and used high-throughput genome-wide translocation sequencing to map translocations from target DNA double-strand breaks (DSBs) within it. RAG endonuclease-cleaved antigen-receptor loci are dominant translocation partners for target DSBs regardless of genomic position, reflecting high-frequency DSBs at these loci and their colocalization in a fraction of cells. To directly assess spatial proximity contributions, we normalized genomic DSBs via ionizing radiation. Under these conditions, translocations were highly enriched in cis along single chromosomes containing target DSBs and within other chromosomes and subchromosomal domains in a manner directly related to pre-existing spatial proximity. By combining two high-throughput genomic methods in a genetically tractable system, we provide a new lens for viewing cancer genomes.

A gene-environment-induced epigenetic program initiates tumorigenesis.

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

Genome-wide translocation sequencing reveals mechanisms of chromosome breaks and rearrangements in B cells.

Whereas chromosomal translocations are common pathogenetic events in cancer, mechanisms that promote them are poorly understood. To elucidate translocation mechanisms in mammalian cells, we developed high-throughput, genome-wide translocation sequencing (HTGTS). We employed HTGTS to identify tens of thousands of independent translocation junctions involving fixed I-SceI meganuclease-generated DNA double-strand breaks (DSBs) within the c-myc oncogene or IgH locus of B lymphocytes induced for activation-induced cytidine deaminase (AID)-dependent IgH class switching. DSBs translocated widely across the genome but were preferentially targeted to transcribed chromosomal regions. Additionally, numerous AID-dependent and AID-independent hot spots were targeted, with the latter comprising mainly cryptic I-SceI targets. Comparison of translocation junctions with genome-wide nuclear run-ons revealed a marked association between transcription start sites and translocation targeting. The majority of translocation junctions were formed via end-joining with short microhomologies. Our findings have implications for diverse fields, including gene therapy and cancer genomics.

Phase and context shape the function of composite oncogenic mutations.

Cancers develop as a result of driver mutations1,2 that lead to clonal outgrowth and the evolution of disease3,4. The discovery and functional characterization of individual driver mutations are central aims of cancer research, and have elucidated myriad phenotypes5 and therapeutic vulnerabilities6. However, the serial genetic evolution of mutant cancer genes7,8 and the allelic context in which they arise is poorly understood in both common and rare cancer genes and tumour types. Here we find that nearly one in four human tumours contains a composite mutation of a cancer-associated gene, defined as two or more nonsynonymous somatic mutations in the same gene and tumour. Composite mutations are enriched in specific genes, have an elevated rate of use of less-common hotspot mutations acquired in a chronology driven in part by oncogenic fitness, and arise in an allelic configuration that reflects context-specific selective pressures. cis-acting composite mutations are hypermorphic in some genes in which dosage effects predominate (such as TERT), whereas they lead to selection of function in other genes (such as TP53). Collectively, composite mutations are driver alterations that arise from context- and allele-specific selective pressures that are dependent in part on gene and mutation function, and which lead to complex-often neomorphic-functions of biological and therapeutic importance.

Senolytic CAR T cells reverse senescence-associated pathologies.

Cellular senescence is characterized by stable cell-cycle arrest and a secretory program that modulates the tissue microenvironment1,2. Physiologically, senescence serves as a tumour-suppressive mechanism that prevents the expansion of premalignant cells3,4 and has a beneficial role in wound-healing responses5,6. Pathologically, the aberrant accumulation of senescent cells generates an inflammatory milieu that leads to chronic tissue damage and contributes to diseases such as liver and lung fibrosis, atherosclerosis, diabetes and osteoarthritis1,7. Accordingly, eliminating senescent cells from damaged tissues in mice ameliorates the symptoms of these pathologies and even promotes longevity1,2,8-10. Here we test the therapeutic concept that chimeric antigen receptor (CAR) T cells that target senescent cells can be effective senolytic agents. We identify the urokinase-type plasminogen activator receptor (uPAR)11 as a cell-surface protein that is broadly induced during senescence and show that uPAR-specific CAR T cells efficiently ablate senescent cells in vitro and in vivo. CAR T cells that target uPAR extend the survival of mice with lung adenocarcinoma that are treated with a senescence-inducing combination of drugs, and restore tissue homeostasis in mice in which liver fibrosis is induced chemically or by diet. These results establish the therapeutic potential of senolytic CAR T cells for senescence-associated diseases.

Senescence induction dictates response to chemo- and immunotherapy in preclinical models of ovarian cancer.

High-grade serous ovarian carcinoma (HGSOC) is a cancer with dismal prognosis due to the limited effectiveness of existing chemo- and immunotherapies. To elucidate mechanisms mediating sensitivity or resistance to these therapies, we developed a fast and flexible autochthonous mouse model based on somatic introduction of HGSOC-associated genetic alterations into the ovary of immunocompetent mice using tissue electroporation. Tumors arising in these mice recapitulate the metastatic patterns and histological, molecular, and treatment response features of the human disease. By leveraging these models, we show that the ability to undergo senescence underlies the clinically observed increase in sensitivity of homologous recombination (HR)-deficient HGSOC tumors to platinum-based chemotherapy. Further, cGas/STING-mediated activation of a restricted senescence-associated secretory phenotype (SASP) was sufficient to induce immune infiltration and sensitize HR-deficient tumors to immune checkpoint blockade. In sum, our study identifies senescence propensity as a predictor of therapy response and defines a limited SASP profile that appears sufficient to confer added vulnerability to concurrent immunotherapy and, more broadly, provides a blueprint for the implementation of electroporation-based mouse models to reveal mechanisms of oncogenesis and therapy response in HGSOC.

A preclinical platform for assessing antitumor effects and systemic toxicities of cancer drug targets.

Anticancer drug development campaigns often fail due to an incomplete understanding of the therapeutic index differentiating the efficacy of the agent against the cancer and its on-target toxicities to the host. To address this issue, we established a versatile preclinical platform in which genetically defined cancers are produced using somatic tissue engineering in transgenic mice harboring a doxycycline-inducible short hairpin RNA against the target of interest. In this system, target inhibition is achieved by the addition of doxycycline, enabling simultaneous assessment of efficacy and toxicity in the same animal. As proof of concept, we focused on CDK9­a cancer target whose clinical development has been hampered by compounds with poorly understood target specificity and unacceptable toxicities. We systematically compared phenotypes produced by genetic Cdk9 inhibition to those achieved using a recently developed highly specific small molecule CDK9 inhibitor and found that both perturbations led to robust antitumor responses. Remarkably, nontoxic levels of CDK9 inhibition could achieve significant treatment efficacy, and dose-dependent toxicities produced by prolonged CDK9 suppression were largely reversible upon Cdk9 restoration or drug withdrawal. Overall, these results establish a versatile in vivo target validation platform that can be employed for rapid triaging of therapeutic targets and lend support to efforts aimed at advancing CDK9 inhibitors for cancer therapy.

Notch activity characterizes a common hepatocellular carcinoma subtype with unique molecular and clinicopathologic features.

BACKGROUND & AIMS: The hepatocyte Notch pathway is a pathogenic factor in non-alcoholic steatohepatitis (NASH)-associated fibrosis, but its role in hepatocellular carcinoma (HCC) is less well defined. Herein, we aimed to characterize the molecular and clinical features of Notch-active human HCC, and to investigate the mechanisms by which Notch affects NASH-driven HCC. METHODS: Using a 14-gene Notch score, we stratified human HCCs from multiple comprehensively profiled datasets. We performed gene set enrichment analyses to compare Notch-active HCCs with published HCC subtype signatures. Next, we sorted Notch-active hepatocytes from Notch reporter mice for RNA sequencing and characterized Notch-active tumors in an HCC model combining a carcinogen and a NASH-inducing diet. We used genetic mouse models to manipulate hepatocyte Notch to investigate the sufficiency and necessity of Notch in NASH-driven tumorigenesis. RESULTS: Notch-active signatures were found in ~30% of human HCCs that transcriptionally resemble cholangiocarcinoma-like HCC, exhibiting a lack of activating CTNNB1 (ß-catenin) mutations and a generally poor prognosis. Endogenous Notch activation in hepatocytes is associated with repressed ß-catenin signaling and hepatic metabolic functions, in lieu of increased interactions with the extracellular matrix in NASH. Constitutive hepatocyte Notch activation is sufficient to induce ß-catenin-inactive HCC in mice with NASH. Notch and ß-catenin show a pattern of mutual exclusivity in carcinogen-induced HCC; in this mouse model, chronic blockade of Notch led to ß-catenin-dependent tumor development. CONCLUSIONS: Notch activity characterizes a distinct HCC molecular subtype with unique histology and prognosis. Sustained Notch signaling in chronic liver diseases can drive tumor formation without acquiring specific genomic driver mutations. LAY SUMMARY: The Notch signaling pathway is known to be involved in the pathogenesis of liver fibrosis. However, its role in liver cancer has not been well defined. Herein, we show that Notch activity is increased in a subset of liver cancers and is associated with poor outcomes. We also used a mouse model to show that aberrant Notch activity can drive cancer progression in obese mice.

Single-cell RNA-seq analysis identifies markers of resistance to targeted BRAF inhibitors in melanoma cell populations.

Single-cell RNA-seq's (scRNA-seq) unprecedented cellular resolution at a genome-wide scale enables us to address questions about cellular heterogeneity that are inaccessible using methods that average over bulk tissue extracts. However, scRNA-seq data sets also present additional challenges such as high transcript dropout rates, stochastic transcription events, and complex population substructures. Here, we present a single-cell RNA-seq analysis and klustering evaluation (SAKE), a robust method for scRNA-seq analysis that provides quantitative statistical metrics at each step of the analysis pipeline. Comparing SAKE to multiple single-cell analysis methods shows that most methods perform similarly across a wide range of cellular contexts, with SAKE outperforming these methods in the case of large complex populations. We next applied the SAKE algorithms to identify drug-resistant cellular populations as human melanoma cells respond to targeted BRAF inhibitors (BRAFi). Single-cell RNA-seq data from both the Fluidigm C1 and 10x Genomics platforms were analyzed with SAKE to dissect this problem at multiple scales. Data from both platforms indicate that BRAF inhibitor-resistant cells can emerge from rare populations already present before drug application, with SAKE identifying both novel and known markers of resistance. These experimentally validated markers of BRAFi resistance share overlap with previous analyses in different melanoma cell lines, demonstrating the generality of these findings and highlighting the utility of single-cell analysis to elucidate mechanisms of BRAFi resistance.

Orientation-specific joining of AID-initiated DNA breaks promotes antibody class switching.

During B-cell development, RAG endonuclease cleaves immunoglobulin heavy chain (IgH) V, D, and J gene segments and orchestrates their fusion as deletional events that assemble a V(D)J exon in the same transcriptional orientation as adjacent Cµ constant region exons. In mice, six additional sets of constant region exons (CHs) lie 100-200 kilobases downstream in the same transcriptional orientation as V(D)J and Cµ exons. Long repetitive switch (S) regions precede Cµ and downstream CHs. In mature B cells, class switch recombination (CSR) generates different antibody classes by replacing Cµ with a downstream CH (ref. 2). Activation-induced cytidine deaminase (AID) initiates CSR by promoting deamination lesions within Sµ and a downstream acceptor S region; these lesions are converted into DNA double-strand breaks (DSBs) by general DNA repair factors. Productive CSR must occur in a deletional orientation by joining the upstream end of an Sµ DSB to the downstream end of an acceptor S-region DSB. However, the relative frequency of deletional to inversional CSR junctions has not been measured. Thus, whether orientation-specific joining is a programmed mechanistic feature of CSR as it is for V(D)J recombination and, if so, how this is achieved is unknown. To address this question, we adapt high-throughput genome-wide translocation sequencing into a highly sensitive DSB end-joining assay and apply it to endogenous AID-initiated S-region DSBs in mouse B cells. We show that CSR is programmed to occur in a productive deletional orientation and does so via an unprecedented mechanism that involves in cis Igh organizational features in combination with frequent S-region DSBs initiated by AID. We further implicate ATM-dependent DSB-response factors in enforcing this mechanism and provide an explanation of why CSR is so reliant on the 53BP1 DSB-response factor.

Genetic and clinical correlates of entosis in pancreatic ductal adenocarcinoma.

Entosis is a type of regulated cell death that promotes cancer cell competition. Though several studies have revealed the molecular mechanisms that govern entosis, the clinical and genetic correlates of entosis in human tumors is less well understood. Here we reviewed entotic cell-in-cell (CIC) patterns in a large single institution sequencing cohort (MSK IMPACT clinical sequencing cohort) of more than 1600 human pancreatic ductal adenocarcinoma (PDAC) samples to identify the genetic and clinical correlates of this cellular feature. After case selection, 516 conventional PDACs and 21 ASCs entered this study and ~45,000 HPFs (median 80 HPFs per sample) were reviewed; 549 entotic-CICs were detected through our cohort. We observed that entotic-CIC occurred more frequently in liver metastasis compared with primary in PDAC. Moreover, poorly differentiated adenocarcinoma or adenosquamous carcinoma had more entotic-CIC than well or moderately differentiated adenocarcinoma. With respect to genetic features TP53 mutations, KRAS amplification, and MYC amplification were significantly associated with entosis in PDAC tissues. From a clinical standpoint entotic CICs were independently associated with a poor prognosis by multivariate Cox regression analysis when considering all cases or primary PDACs specifically. These results provide a contextual basis for understanding entosis in PDAC, a highly aggressive cancer for which molecular insights are needed to improve survival.

Clonal dynamics of native haematopoiesis.

It is currently thought that life-long blood cell production is driven by the action of a small number of multipotent haematopoietic stem cells. Evidence supporting this view has been largely acquired through the use of functional assays involving transplantation. However, whether these mechanisms also govern native non-transplant haematopoiesis is entirely unclear. Here we have established a novel experimental model in mice where cells can be uniquely and genetically labelled in situ to address this question. Using this approach, we have performed longitudinal analyses of clonal dynamics in adult mice that reveal unprecedented features of native haematopoiesis. In contrast to what occurs following transplantation, steady-state blood production is maintained by the successive recruitment of thousands of clones, each with a minimal contribution to mature progeny. Our results demonstrate that a large number of long-lived progenitors, rather than classically defined haematopoietic stem cells, are the main drivers of steady-state haematopoiesis during most of adulthood. Our results also have implications for understanding the cellular origin of haematopoietic disease.

Timing of treatment shapes the path to androgen receptor signaling inhibitor resistance in prostate cancer.

There is optimism that cancer drug resistance can be addressed through appropriate combination therapy, but success requires understanding the growing complexity of resistance mechanisms, including the evolution and population dynamics of drug-sensitive and drug-resistant clones over time. Using DNA barcoding to trace individual prostate tumor cells in vivo , we find that the evolutionary path to acquired resistance to androgen receptor signaling inhibition (ARSI) is dependent on the timing of treatment. In established tumors, resistance occurs through polyclonal adaptation of drug-sensitive clones, despite the presence of rare subclones with known, pre-existing ARSI resistance. Conversely, in an experimental setting designed to mimic minimal residual disease, resistance occurs through outgrowth of pre-existing resistant clones and not by adaptation. Despite these different evolutionary paths, the underlying mechanisms responsible for resistance are shared across the two evolutionary paths. Furthermore, mixing experiments reveal that the evolutionary path to adaptive resistance requires cooperativity between subclones. Thus, despite the presence of pre-existing ARSI-resistant subclones, acquired resistance in established tumors occurs primarily through cooperative, polyclonal adaptation of drug-sensitive cells. This tumor ecosystem model of resistance has new implications for developing effective combination therapy.

Somatic mouse models of gastric cancer reveal genotype-specific features of metastatic disease.

Metastatic gastric carcinoma is a highly lethal cancer that responds poorly to conventional and molecularly targeted therapies. Despite its clinical relevance, the mechanisms underlying the behavior and therapeutic response of this disease are poorly understood owing, in part, to a paucity of tractable models. Here we developed methods to somatically introduce different oncogenic lesions directly into the murine gastric epithelium. Genotypic configurations observed in patients produced metastatic gastric cancers that recapitulated the histological, molecular and clinical features of all nonviral molecular subtypes of the human disease. Applying this platform to both wild-type and immunodeficient mice revealed previously unappreciated links between the genotype, organotropism and immune surveillance of metastatic cells, which produced distinct patterns of metastasis that were mirrored in patients. Our results establish a highly portable platform for generating autochthonous cancer models with flexible genotypes and host backgrounds, which can unravel mechanisms of gastric tumorigenesis or test new therapeutic concepts.