Lineage tracing reveals the phylodynamics, plasticity, and paths of tumor evolution.

Tumor evolution is driven by the progressive acquisition of genetic and epigenetic alterations that enable uncontrolled growth and expansion to neighboring and distal tissues. The study of phylogenetic relationships between cancer cells provides key insights into these processes. Here, we introduced an evolving lineage-tracing system with a single-cell RNA-seq readout into a mouse model of Kras;Trp53(KP)-driven lung adenocarcinoma and tracked tumor evolution from single-transformed cells to metastatic tumors at unprecedented resolution. We found that the loss of the initial, stable alveolar-type2-like state was accompanied by a transient increase in plasticity. This was followed by the adoption of distinct transcriptional programs that enable rapid expansion and, ultimately, clonal sweep of stable subclones capable of metastasizing. Finally, tumors develop through stereotypical evolutionary trajectories, and perturbing additional tumor suppressors accelerates progression by creating novel trajectories. Our study elucidates the hierarchical nature of tumor evolution and, more broadly, enables in-depth studies of tumor progression.

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.

Conventional type I dendritic cells maintain a reservoir of proliferative tumor-antigen specific TCF-1+ CD8+ T cells in tumor-draining lymph nodes.

In tumors, a subset of CD8+ T cells expressing the transcription factor TCF-1 drives the response to immune checkpoint blockade. We examined the mechanisms that maintain these cells in an autochthonous model of lung adenocarcinoma. Longitudinal sampling and single-cell sequencing of tumor-antigen specific TCF-1+ CD8+ T cells revealed that while intratumoral TCF-1+ CD8+ T cells acquired dysfunctional features and decreased in number as tumors progressed, TCF-1+ CD8+ T cell frequency in the tumor draining LN (dLN) remained stable. Two discrete intratumoral TCF-1+ CD8+ T cell subsets developed over time-a proliferative SlamF6+ subset and a non-cycling SlamF6- subset. Blocking dLN egress decreased the frequency of intratumoral SlamF6+ TCF-1+ CD8+ T cells. Conventional type I dendritic cell (cDC1) in dLN decreased in number with tumor progression, and Flt3L+anti-CD40 treatment recovered SlamF6+ T cell frequencies and decreased tumor burden. Thus, cDC1s in tumor dLN maintain a reservoir of TCF-1+ CD8+ T cells and their decrease contributes to failed anti-tumor immunity.

Deciphering the immunopeptidome in vivo reveals new tumour antigens.

Immunosurveillance of cancer requires the presentation of peptide antigens on major histocompatibility complex class I (MHC-I) molecules1-5. Current approaches to profiling of MHC-I-associated peptides, collectively known as the immunopeptidome, are limited to in vitro investigation or bulk tumour lysates, which limits our understanding of cancer-specific patterns of antigen presentation in vivo6. To overcome these limitations, we engineered an inducible affinity tag into the mouse MHC-I gene (H2-K1) and targeted this allele to the KrasLSL-G12D/+Trp53fl/fl mouse model (KP/KbStrep)7. This approach enabled us to precisely isolate MHC-I peptides from autochthonous pancreatic ductal adenocarcinoma and from lung adenocarcinoma (LUAD) in vivo. In addition, we profiled the LUAD immunopeptidome from the alveolar type 2 cell of origin up to late-stage disease. Differential peptide presentation in LUAD was not predictable by mRNA expression or translation efficiency and is probably driven by post-translational mechanisms. Vaccination with peptides presented by LUAD in vivo induced CD8+ T cell responses in naive mice and tumour-bearing mice. Many peptides specific to LUAD, including immunogenic peptides, exhibited minimal expression of the cognate mRNA, which prompts the reconsideration of antigen prediction pipelines that triage peptides according to transcript abundance8. Beyond cancer, the KbStrep allele is compatible with other Cre-driver lines to explore antigen presentation in vivo in the pursuit of understanding basic immunology, infectious disease and autoimmunity.

Mutant p53 disrupts mammary tissue architecture via the mevalonate pathway.

p53 is a frequent target for mutation in human tumors, and mutant p53 proteins can actively contribute to tumorigenesis. We employed a three-dimensional culture model in which nonmalignant breast epithelial cells form spheroids reminiscent of acinar structures found in vivo, whereas breast cancer cells display highly disorganized morphology. We found that mutant p53 depletion is sufficient to phenotypically revert breast cancer cells to a more acinar-like morphology. Genome-wide expression analysis identified the mevalonate pathway as significantly upregulated by mutant p53. Statins and sterol biosynthesis intermediates reveal that this pathway is both necessary and sufficient for the phenotypic effects of mutant p53 on breast tissue architecture. Mutant p53 associates with sterol gene promoters at least partly via SREBP transcription factors. Finally, p53 mutation correlates with highly expressed sterol biosynthesis genes in human breast tumors. These findings implicate the mevalonate pathway as a therapeutic target for tumors bearing mutations in p53.

Mutant p53 cooperates with the SWI/SNF chromatin remodeling complex to regulate VEGFR2 in breast cancer cells.

Mutant p53 impacts the expression of numerous genes at the level of transcription to mediate oncogenesis. We identified vascular endothelial growth factor receptor 2 (VEGFR2), the primary functional VEGF receptor that mediates endothelial cell vascularization, as a mutant p53 transcriptional target in multiple breast cancer cell lines. Up-regulation of VEGFR2 mediates the role of mutant p53 in increasing cellular growth in two-dimensional (2D) and three-dimensional (3D) culture conditions. Mutant p53 binds near the VEGFR2 promoter transcriptional start site and plays a role in maintaining an open conformation at that location. Relatedly, mutant p53 interacts with the SWI/SNF complex, which is required for remodeling the VEGFR2 promoter. By both querying individual genes regulated by mutant p53 and performing RNA sequencing, the results indicate that >40% of all mutant p53-regulated gene expression is mediated by SWI/SNF. We surmise that mutant p53 impacts transcription of VEGFR2 as well as myriad other genes by promoter remodeling through interaction with and likely regulation of the SWI/SNF chromatin remodeling complex. Therefore, not only might mutant p53-expressing tumors be susceptible to anti VEGF therapies, impacting SWI/SNF tumor suppressor function in mutant p53 tumors may also have therapeutic potential.

Accurate and sensitive quantification of protein-DNA binding affinity.

Transcription factors (TFs) control gene expression by binding to genomic DNA in a sequence-specific manner. Mutations in TF binding sites are increasingly found to be associated with human disease, yet we currently lack robust methods to predict these sites. Here, we developed a versatile maximum likelihood framework named No Read Left Behind (NRLB) that infers a biophysical model of protein-DNA recognition across the full affinity range from a library of in vitro selected DNA binding sites. NRLB predicts human Max homodimer binding in near-perfect agreement with existing low-throughput measurements. It can capture the specificity of the p53 tetramer and distinguish multiple binding modes within a single sample. Additionally, we confirm that newly identified low-affinity enhancer binding sites are functional in vivo, and that their contribution to gene expression matches their predicted affinity. Our results establish a powerful paradigm for identifying protein binding sites and interpreting gene regulatory sequences in eukaryotic genomes.

Mutant p53: one name, many proteins.

There is now strong evidence that mutation not only abrogates p53 tumor-suppressive functions, but in some instances can also endow mutant proteins with novel activities. Such neomorphic p53 proteins are capable of dramatically altering tumor cell behavior, primarily through their interactions with other cellular proteins and regulation of cancer cell transcriptional programs. Different missense mutations in p53 may confer unique activities and thereby offer insight into the mutagenic events that drive tumor progression. Here we review mechanisms by which mutant p53 exerts its cellular effects, with a particular focus on the burgeoning mutant p53 transcriptome, and discuss the biological and clinical consequences of mutant p53 gain of function.

A prime editor mouse to model a broad spectrum of somatic mutations in vivo.

Genetically engineered mouse models only capture a small fraction of the genetic lesions that drive human cancer. Current CRISPR-Cas9 models can expand this fraction but are limited by their reliance on error-prone DNA repair. Here we develop a system for in vivo prime editing by encoding a Cre-inducible prime editor in the mouse germline. This model allows rapid, precise engineering of a wide range of mutations in cell lines and organoids derived from primary tissues, including a clinically relevant Kras mutation associated with drug resistance and Trp53 hotspot mutations commonly observed in pancreatic cancer. With this system, we demonstrate somatic prime editing in vivo using lipid nanoparticles, and we model lung and pancreatic cancer through viral delivery of prime editing guide RNAs or orthotopic transplantation of prime-edited organoids. We believe that this approach will accelerate functional studies of cancer-associated mutations and complex genetic combinations that are challenging to construct with traditional models.

Getting a handle on KRAS inhibitor resistance with hapten-mediated anti-tumor immunity.

Covalent inhibitors of oncogenic KRASG12C have demonstrated impressive clinical responses; however, therapeutic resistance has been commonly observed. In this issue, Zhang and colleagues demonstrate that small molecule KRASG12C inhibitors can generate haptenated major histocompatibility complex (MHC) class I:peptide complexes, which represent attractive targets for immune-based therapies to combat pharmacologic resistance.

Emerging Role of Targeted Therapy in Metastatic Pancreatic Adenocarcinoma.

The aggressive biology of pancreatic ductal adenocarcinoma (PDAC), along with its limited sensitivity to many systemic therapies, presents a major challenge in the management of patients with metastatic PDAC. Over the past decade, the incorporation of combinatorial cytotoxic chemotherapy regimens has improved patient outcomes. Despite these advances, resistance to cytotoxic chemotherapy inevitably occurs, and there is a great need for effective therapies. A major focus of research has been to identify molecularly defined subpopulations of patients with PDAC who may benefit from targeted therapies that are matched to their molecular profile. Recent successes include the demonstration of the efficacy of maintenance PARP inhibition in PDAC tumors harboring deleterious BRCA1, BRCA2, and PALB2 alterations. In addition, while therapeutic targeting of KRAS was long thought to be infeasible, emerging data on the efficacy of KRAS G12C inhibitors have increased optimism about next-generation KRAS-directed therapies in PDAC. Meanwhile, KRAS wild-type PDAC encompasses a unique molecular subpopulation of PDAC that is enriched for targetable genetic alterations, such as oncogenic BRAF alterations, mismatch repair deficiency, and FGFR2, ALK, NTRK, ROS1, NRG1, and RET rearrangements. As more molecularly targeted therapies are developed, precision medicine has the potential to revolutionize the treatment of patients with metastatic PDAC.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Measuring kinetics and metastatic propensity of CTCs by blood exchange between mice.

Existing preclinical methods for acquiring dissemination kinetics of rare circulating tumor cells (CTCs) en route to forming metastases have not been capable of providing a direct measure of CTC intravasation rate and subsequent half-life in the circulation. Here, we demonstrate an approach for measuring endogenous CTC kinetics by continuously exchanging CTC-containing blood over several hours between un-anesthetized, tumor-bearing mice and healthy, tumor-free counterparts. By tracking CTC transfer rates, we extrapolated half-life times in the circulation of between 40 and 260 s and intravasation rates between 60 and 107,000 CTCs/hour in mouse models of small-cell lung cancer (SCLC), pancreatic ductal adenocarcinoma (PDAC), and non-small cell lung cancer (NSCLC). Additionally, direct transfer of only 1-2% of daily-shed CTCs using our blood-exchange technique from late-stage, SCLC-bearing mice generated macrometastases in healthy recipient mice. We envision that our technique will help further elucidate the role of CTCs and the rate-limiting steps in metastasis.

The CD155/TIGIT axis promotes and maintains immune evasion in neoantigen-expressing pancreatic cancer.

The CD155/TIGIT axis can be co-opted during immune evasion in chronic viral infections and cancer. Pancreatic adenocarcinoma (PDAC) is a highly lethal malignancy, and immune-based strategies to combat this disease have been largely unsuccessful to date. We corroborate prior reports that a substantial portion of PDAC harbors predicted high-affinity MHC class I-restricted neoepitopes and extend these findings to advanced/metastatic disease. Using multiple preclinical models of neoantigen-expressing PDAC, we demonstrate that intratumoral neoantigen-specific CD8+ T cells adopt multiple states of dysfunction, resembling those in tumor-infiltrating lymphocytes of PDAC patients. Mechanistically, genetic and/or pharmacologic modulation of the CD155/TIGIT axis was sufficient to promote immune evasion in autochthonous neoantigen-expressing PDAC. Finally, we demonstrate that the CD155/TIGIT axis is critical in maintaining immune evasion in PDAC and uncover a combination immunotherapy (TIGIT/PD-1 co-blockade plus CD40 agonism) that elicits profound anti-tumor responses in preclinical models, now poised for clinical evaluation.

Quantitative Analysis of the DNA Methylation Sensitivity of Transcription Factor Complexes.

Although DNA modifications play an important role in gene regulation, the underlying mechanisms remain elusive. We developed EpiSELEX-seq to probe the sensitivity of transcription factor binding to DNA modification in vitro using massively parallel sequencing. Feature-based modeling quantifies the effect of cytosine methylation (5mC) on binding free energy in a position-specific manner. Application to the human bZIP proteins ATF4 and C/EBPß and three different Pbx-Hox complexes shows that 5mCpG can both increase and decrease affinity, depending on where the modification occurs within the protein-DNA interface. The TF paralogs tested vary in their methylation sensitivity, for which we provide a structural rationale. We show that 5mCpG can also enhance in vitro p53 binding and provide evidence for increased in vivo p53 occupancy at methylated binding sites, correlating with primed enhancer histone marks. Our results establish a powerful strategy for dissecting the epigenomic modulation of protein-DNA interactions and their role in gene regulation.

Stability of checkpoint kinase 2 is regulated via phosphorylation at serine 456.

Checkpoint kinase 2 (Chk2), a DNA damage-activated protein kinase, is phosphorylated at Thr-68 by ataxia telangiectasia mutated leading to its activation by phosphorylation at several additional sites. Using mass spectrometry we identified a new Chk2 phosphorylation site at Ser-456. We show that phosphorylation of Ser-456 plays a role in the regulation of Chk2 stability particularly after DNA damage. Mutation of Ser-456 to alanine results in hyperubiquitination of Chk2 and dramatically reduced Chk2 stability. Furthermore, cells expressing S456A Chk2 show a reduction in the apoptotic response to DNA damage. These findings suggest a mechanism for stabilization of Chk2 in response to DNA damage via phosphorylation at Ser-456 and proteasome-dependent turnover of Chk2 protein via dephosphorylation of the same residue.