Mechanisms of resistance to oncogenic KRAS inhibition in pancreatic cancer.

KRAS inhibitors demonstrate clinical efficacy in pancreatic ductal adenocarcinoma (PDAC); however, resistance is common. Among patients with KRASG12C-mutant PDAC treated with adagrasib or sotorasib, mutations in PIK3CA and KRAS, and amplifications of KRASG12C, MYC, MET, EGFR, and CDK6 emerged at acquired resistance. In PDAC cell lines and organoid models treated with the KRASG12D inhibitor MRTX1133, epithelial-to-mesenchymal transition and PI3K-AKT-mTOR signaling associate with resistance to therapy. MRTX1133 treatment of the KrasLSL-G12D/+;Trp53LSL-R172H/+;p48-Cre (KPC) mouse model yielded deep tumor regressions, but drug resistance ultimately emerged, accompanied by amplifications of Kras, Yap1, Myc, and Cdk6/Abcb1a/b, and co-evolution of drug-resistant transcriptional programs. Moreover, in KPC and PDX models, mesenchymal and basal-like cell states displayed increased response to KRAS inhibition compared to the classical state. Combination treatment with KRASG12D inhibition and chemotherapy significantly improved tumor control in PDAC mouse models. Collectively, these data elucidate co-evolving resistance mechanisms to KRAS inhibition and support multiple combination therapy strategies.

Metastatic site influences driver gene function in pancreatic cancer.

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFß signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFß's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.

Plasticity-induced repression of Irf6 underlies acquired resistance to cancer immunotherapy in pancreatic ductal adenocarcinoma.

Acquired resistance to immunotherapy remains a critical yet incompletely understood biological mechanism. Here, using a mouse model of pancreatic ductal adenocarcinoma (PDAC) to study tumor relapse following immunotherapy-induced responses, we find that resistance is reproducibly associated with an epithelial-to-mesenchymal transition (EMT), with EMT-transcription factors ZEB1 and SNAIL functioning as master genetic and epigenetic regulators of this effect. Acquired resistance in this model is not due to immunosuppression in the tumor immune microenvironment, disruptions in the antigen presentation machinery, or altered expression of immune checkpoints. Rather, resistance is due to a tumor cell-intrinsic defect in T-cell killing. Molecularly, EMT leads to the epigenetic and transcriptional silencing of interferon regulatory factor 6 (Irf6), rendering tumor cells less sensitive to the pro-apoptotic effects of TNF-α. These findings indicate that acquired resistance to immunotherapy may be mediated by programs distinct from those governing primary resistance, including plasticity programs that render tumor cells impervious to T-cell killing.

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.

MACHETE identifies interferon-encompassing chromosome 9p21.3 deletions as mediators of immune evasion and metastasis.

The most prominent homozygous deletions in cancer affect chromosome 9p21.3 and eliminate CDKN2A/B tumor suppressors, disabling a cell-intrinsic barrier to tumorigenesis. Half of 9p21.3 deletions, however, also encompass a type I interferon (IFN) gene cluster; the consequences of this co-deletion remain unexplored. To functionally dissect 9p21.3 and other large genomic deletions, we developed a flexible deletion engineering strategy, MACHETE (molecular alteration of chromosomes with engineered tandem elements). Applying MACHETE to a syngeneic mouse model of pancreatic cancer, we found that co-deletion of the IFN cluster promoted immune evasion, metastasis and immunotherapy resistance. Mechanistically, IFN co-deletion disrupted type I IFN signaling in the tumor microenvironment, leading to marked changes in infiltrating immune cells and escape from CD8+ T-cell surveillance, effects largely driven by the poorly understood interferon epsilon. These results reveal a chromosomal deletion that disables both cell-intrinsic and cell-extrinsic tumor suppression and provide a framework for interrogating large deletions in cancer and beyond.

Measurable residual disease by flow cytometry in acute myeloid leukemia is prognostic, independent of genomic profiling.

Measurable residual disease (MRD) assessment provides a potent indicator of the efficacy of anti-leukemic therapy. It is unknown, however, whether integrating MRD with molecular profiling better identifies patients at risk of relapse. To investigate the clinical relevance of MRD in relation to a molecular-based prognostic schema, we measured MRD by flow cytometry in 189 AML patients enrolled in ECOG-ACRIN E1900 trial (NCT00049517) in morphologic complete remission (CR) (28.8 % of the original cohort) representing 44.4 % of CR patients. MRD positivity was defined as ≥ 0.1 % of leukemic bone marrow cells. Risk classification was based on standard cytogenetics, fluorescence-in-situ-hybridization, somatic gene analysis, and sparse whole genome sequencing for copy number ascertainment. At 84.6 months median follow-up of patients still alive at the time of analysis (range 47.0-120 months), multivariate analysis demonstrated that MRD status at CR (p = 0.001) and integrated molecular risk (p = 0.0004) independently predicted overall survival (OS). Among risk classes, MRD status significantly affected OS only in the favorable risk group (p = 0.002). Expression of CD25 (α-chain of the interleukin-2 receptor) by leukemic myeloblasts at diagnosis negatively affected OS independent of post-treatment MRD levels. These data suggest that integrating MRD with genetic profiling and pre-treatment CD25 expression may improve prognostication in AML.

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.

Mcm2 hypomorph leads to acute leukemia or hematopoietic stem cell failure, dependent on genetic context.

Minichromosome maintenance proteins (Mcm2-7) form a hexameric complex that unwinds DNA ahead of a replicative fork. The deficiency of Mcm proteins leads to replicative stress and consequent genomic instability. Mice with a germline insertion of a Cre cassette into the 3'UTR of the Mcm2 gene (designated Mcm2Cre ) have decreased Mcm2 expression and invariably develop precursor T-cell lymphoblastic leukemia/lymphoma (pre-T LBL), due to 100-1000 kb deletions involving important tumor suppressor genes. To determine whether mice that were protected from pre-T LBL would develop non-T-cell malignancies, we used two approaches. Mice engrafted with Mcm2Cre/Cre Lin- Sca-1+ Kit+ hematopoietic stem/progenitor cells did not develop hematologic malignancy; however, these mice died of hematopoietic stem cell failure by 6 months of age. Placing the Mcm2Cre allele onto an athymic nu/nu background completely prevented pre-T LBL and extended survival of these mice three-fold (median 296.5 vs. 80.5 days). Ultimately, most Mcm2Cre/Cre ;nu/nu mice developed B-cell precursor acute lymphoblastic leukemia (BCP-ALL). We identified recurrent deletions of 100-1000 kb that involved genes known or suspected to be involved in BCP-ALL, including Pax5, Nf1, Ikzf3, and Bcor. Moreover, whole-exome sequencing identified recurrent mutations of genes known to be involved in BCP-ALL progression, such as Jak1/Jak3, Ptpn11, and Kras. These findings demonstrate that an Mcm2Cre/Cre hypomorph can induce hematopoietic dysfunction via hematopoietic stem cell failure as well as a "deletor" phenotype affecting known or suspected tumor suppressor genes.

Human OTULIN haploinsufficiency impairs cell-intrinsic immunity to staphylococcal α-toxin.

The molecular basis of interindividual clinical variability upon infection with Staphylococcus aureus is unclear. We describe patients with haploinsufficiency for the linear deubiquitinase OTULIN, encoded by a gene on chromosome 5p. Patients suffer from episodes of life-threatening necrosis, typically triggered by S. aureus infection. The disorder is phenocopied in patients with the 5p- (Cri-du-Chat) chromosomal deletion syndrome. OTULIN haploinsufficiency causes an accumulation of linear ubiquitin in dermal fibroblasts, but tumor necrosis factor receptor-mediated nuclear factor κB signaling remains intact. Blood leukocyte subsets are unaffected. The OTULIN-dependent accumulation of caveolin-1 in dermal fibroblasts, but not leukocytes, facilitates the cytotoxic damage inflicted by the staphylococcal virulence factor α-toxin. Naturally elicited antibodies against α-toxin contribute to incomplete clinical penetrance. Human OTULIN haploinsufficiency underlies life-threatening staphylococcal disease by disrupting cell-intrinsic immunity to α-toxin in nonleukocytic cells.

BMP2/SMAD pathway activation in JAK2/p53-mutant megakaryocyte/erythroid progenitors promotes leukemic transformation.

Leukemic transformation (LT) of myeloproliferative neoplasm (MPN) has a dismal prognosis and is largely fatal. Mutational inactivation of TP53 is the most common somatic event in LT; however, the mechanisms by which TP53 mutations promote LT remain unresolved. Using an allelic series of mouse models of Jak2/Trp53 mutant MPN, we identify that only biallelic inactivation of Trp53 results in LT (to a pure erythroleukemia [PEL]). This PEL arises from the megakaryocyte-erythroid progenitor population. Importantly, the bone morphogenetic protein 2/SMAD pathway is aberrantly activated during LT and results in abnormal self-renewal of megakaryocyte-erythroid progenitors. Finally, we identify that Jak2/Trp53 mutant PEL is characterized by recurrent copy number alterations and DNA damage. Using a synthetic lethality strategy, by targeting active DNA repair pathways, we show that this PEL is highly sensitive to combination WEE1 and poly(ADP-ribose) polymerase inhibition. These observations yield new mechanistic insights into the process of p53 mutant LT and offer new, clinically translatable therapeutic approaches.

Disruption of the crypt niche promotes outgrowth of mutated colorectal tumor stem cells.

Recent data establish a logarithmic expansion of leucine rich repeat containing G protein coupled receptor 5-positive (Lgr5+) colonic epithelial stem cells (CESCs) in human colorectal cancer (CRC). Complementary studies using the murine 2-stage azoxymethane-dextran sulfate sodium (AOM-DSS) colitis-associated tumor model indicate early acquisition of Wnt pathway mutations drives CESC expansion during adenoma progression. Here, subdivision of the AOM-DSS model into in vivo and in vitro stages revealed DSS induced physical separation of CESCs from stem cell niche cells and basal lamina, a source of Wnt signals, within hours, disabling the stem cell program. While AOM delivery in vivo under non-adenoma-forming conditions yielded phenotypically normal mucosa and organoids derived thereof, niche injury ex vivo by progressive DSS dose escalation facilitated outgrowth of Wnt-independent dysplastic organoids. These organoids contained 10-fold increased Lgr5+ CESCs with gain-of-function Wnt mutations orthologous to human CRC driver mutations. We posit CRC originates by niche injury-induced outgrowth of normally suppressed mutated stem cells, consistent with models of adaptive oncogenesis.

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.

MYC Levels Regulate Metastatic Heterogeneity in Pancreatic Adenocarcinoma.

The degree of metastatic disease varies widely among patients with cancer and affects clinical outcomes. However, the biological and functional differences that drive the extent of metastasis are poorly understood. We analyzed primary tumors and paired metastases using a multifluorescent lineage-labeled mouse model of pancreatic ductal adenocarcinoma (PDAC)-a tumor type in which most patients present with metastases. Genomic and transcriptomic analysis revealed an association between metastatic burden and gene amplification or transcriptional upregulation of MYC and its downstream targets. Functional experiments showed that MYC promotes metastasis by recruiting tumor-associated macrophages, leading to greater bloodstream intravasation. Consistent with these findings, metastatic progression in human PDAC was associated with activation of MYC signaling pathways and enrichment for MYC amplifications specifically in metastatic patients. Collectively, these results implicate MYC activity as a major determinant of metastatic burden in advanced PDAC. SIGNIFICANCE: Here, we investigate metastatic variation seen clinically in patients with PDAC and murine PDAC tumors and identify MYC as a major driver of this heterogeneity.This article is highlighted in the In This Issue feature, p. 275.

Clonal populations of a human TNBC model display significant functional heterogeneity and divergent growth dynamics in distinct contexts.

Intratumoral heterogeneity has been described for various tumor types and models of human cancer, and can have profound effects on tumor progression and drug resistance. This study describes an in-depth analysis of molecular and functional heterogeneity among subclonal populations (SCPs) derived from a single triple-negative breast cancer cell line, including copy number analysis, whole-exome and RNA sequencing, proteome analysis, and barcode analysis of clonal dynamics, as well as functional assays. The SCPs were found to have multiple unique genetic alterations and displayed significant variation in anchorage independent growth and tumor forming ability. Analyses of clonal dynamics in SCP mixtures using DNA barcode technology revealed selection for distinct clonal populations in different in vitro and in vivo environmental contexts, demonstrating that in vitro propagation of cancer cell lines using different culture conditions can contribute to the establishment of unique strains. These analyses also revealed strong enrichment of a single SCP during the development of xenograft tumors in immune-compromised mice. This SCP displayed attenuated interferon signaling in vivo and reduced sensitivity to the antiproliferative effects of type I interferons. Reduction in interferon signaling was found to provide a selective advantage within the xenograft microenvironment specifically. In concordance with the previously described role of interferon signaling as tumor suppressor, these findings suggest that similar selective pressures may be operative in human cancer and patient-derived xenograft models.

High resolution copy number inference in cancer using short-molecule nanopore sequencing.

Genome copy number is an important source of genetic variation in health and disease. In cancer, Copy Number Alterations (CNAs) can be inferred from short-read sequencing data, enabling genomics-based precision oncology. Emerging Nanopore sequencing technologies offer the potential for broader clinical utility, for example in smaller hospitals, due to lower instrument cost, higher portability, and ease of use. Nonetheless, Nanopore sequencing devices are limited in the number of retrievable sequencing reads/molecules compared to short-read sequencing platforms, limiting CNA inference accuracy. To address this limitation, we targeted the sequencing of short-length DNA molecules loaded at optimized concentration in an effort to increase sequence read/molecule yield from a single nanopore run. We show that short-molecule nanopore sequencing reproducibly returns high read counts and allows high quality CNA inference. We demonstrate the clinical relevance of this approach by accurately inferring CNAs in acute myeloid leukemia samples. The data shows that, compared to traditional approaches such as chromosome analysis/cytogenetics, short molecule nanopore sequencing returns more sensitive, accurate copy number information in a cost effective and expeditious manner, including for multiplex samples. Our results provide a framework for short-molecule nanopore sequencing with applications in research and medicine, which includes but is not limited to, CNAs.

Leukemia Cell of Origin Influences Apoptotic Priming and Sensitivity to LSD1 Inhibition.

The cell of origin of oncogenic transformation is a determinant of therapeutic sensitivity, but the mechanisms governing cell-of-origin-driven differences in therapeutic response have not been delineated. Leukemias initiating in hematopoietic stem cells (HSC) are less sensitive to chemotherapy and highly express the transcription factor MECOM (EVI1) compared with leukemias derived from myeloid progenitors. Here, we compared leukemias initiated in either HSCs or myeloid progenitors to reveal a novel function for EVI1 in modulating p53 protein abundance and activity. HSC-derived leukemias exhibit decreased apoptotic priming, attenuated p53 transcriptional output, and resistance to lysine-specific demethylase 1 (LSD1) inhibitors in addition to classical genotoxic stresses. p53 loss of function in Evi1 lo progenitor-derived leukemias induces resistance to LSD1 inhibition, and EVI1hi leukemias are sensitized to LSD1 inhibition by venetoclax. Our findings demonstrate a role for EVI1 in p53 wild-type cancers in reducing p53 function and provide a strategy to circumvent drug resistance in chemoresistant EVI1 hi acute myeloid leukemia. SIGNIFICANCE: We demonstrate that the cell of origin of leukemia initiation influences p53 activity and dictates therapeutic sensitivity to pharmacologic LSD1 inhibitors via the transcription factor EVI1. We show that drug resistance could be overcome in HSC-derived leukemias by combining LSD1 inhibition with venetoclax.See related commentary by Gu et al., p. 1445.This article is highlighted in the In This Issue feature, p. 1426.

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.