Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Spatially confined sub-tumor microenvironments in pancreatic cancer.

Intratumoral heterogeneity is a critical frontier in understanding how the tumor microenvironment (TME) propels malignant progression. Here, we deconvolute the human pancreatic TME through large-scale integration of histology-guided regional multiOMICs with clinical data and patient-derived preclinical models. We discover "subTMEs," histologically definable tissue states anchored in fibroblast plasticity, with regional relationships to tumor immunity, subtypes, differentiation, and treatment response. "Reactive" subTMEs rich in complex but functionally coordinated fibroblast communities were immune hot and inhabited by aggressive tumor cell phenotypes. The matrix-rich "deserted" subTMEs harbored fewer activated fibroblasts and tumor-suppressive features yet were markedly chemoprotective and enriched upon chemotherapy. SubTMEs originated in fibroblast differentiation trajectories, and transitory states were notable both in single-cell transcriptomics and in situ. The intratumoral co-occurrence of subTMEs produced patient-specific phenotypic and computationally predictable heterogeneity tightly linked to malignant biology. Therefore, heterogeneity within the plentiful, notorious pancreatic TME is not random but marks fundamental tissue organizational units.

Targeting Pin1 renders pancreatic cancer eradicable by synergizing with immunochemotherapy.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by notorious resistance to current therapies attributed to inherent tumor heterogeneity and highly desmoplastic and immunosuppressive tumor microenvironment (TME). Unique proline isomerase Pin1 regulates multiple cancer pathways, but its role in the TME and cancer immunotherapy is unknown. Here, we find that Pin1 is overexpressed both in cancer cells and cancer-associated fibroblasts (CAFs) and correlates with poor survival in PDAC patients. Targeting Pin1 using clinically available drugs induces complete elimination or sustained remissions of aggressive PDAC by synergizing with anti-PD-1 and gemcitabine in diverse model systems. Mechanistically, Pin1 drives the desmoplastic and immunosuppressive TME by acting on CAFs and induces lysosomal degradation of the PD-1 ligand PD-L1 and the gemcitabine transporter ENT1 in cancer cells, besides activating multiple cancer pathways. Thus, Pin1 inhibition simultaneously blocks multiple cancer pathways, disrupts the desmoplastic and immunosuppressive TME, and upregulates PD-L1 and ENT1, rendering PDAC eradicable by immunochemotherapy.

Neurons Release Serine to Support mRNA Translation in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) tumors have a nutrient-poor, desmoplastic, and highly innervated tumor microenvironment. Although neurons can release stimulatory factors to accelerate PDAC tumorigenesis, the metabolic contribution of peripheral axons has not been explored. We found that peripheral axons release serine (Ser) to support the growth of exogenous Ser (exSer)-dependent PDAC cells during Ser/Gly (glycine) deprivation. Ser deprivation resulted in ribosomal stalling on two of the six Ser codons, TCC and TCT, and allowed the selective translation and secretion of nerve growth factor (NGF) by PDAC cells to promote tumor innervation. Consistent with this, exSer-dependent PDAC tumors grew slower and displayed enhanced innervation in mice on a Ser/Gly-free diet. Blockade of compensatory neuronal innervation using LOXO-101, a Trk-NGF inhibitor, further decreased PDAC tumor growth. Our data indicate that axonal-cancer metabolic crosstalk is a critical adaptation to support PDAC growth in nutrient poor environments.

A Multiscale Map of the Stem Cell State in Pancreatic Adenocarcinoma.

Drug resistance and relapse remain key challenges in pancreatic cancer. Here, we have used RNA sequencing (RNA-seq), chromatin immunoprecipitation (ChIP)-seq, and genome-wide CRISPR analysis to map the molecular dependencies of pancreatic cancer stem cells, highly therapy-resistant cells that preferentially drive tumorigenesis and progression. This integrated genomic approach revealed an unexpected utilization of immuno-regulatory signals by pancreatic cancer epithelial cells. In particular, the nuclear hormone receptor retinoic-acid-receptor-related orphan receptor gamma (RORγ), known to drive inflammation and T cell differentiation, was upregulated during pancreatic cancer progression, and its genetic or pharmacologic inhibition led to a striking defect in pancreatic cancer growth and a marked improvement in survival. Further, a large-scale retrospective analysis in patients revealed that RORγ expression may predict pancreatic cancer aggressiveness, as it positively correlated with advanced disease and metastasis. Collectively, these data identify an orthogonal co-option of immuno-regulatory signals by pancreatic cancer stem cells, suggesting that autoimmune drugs should be evaluated as novel treatment strategies for pancreatic cancer patients.

Divergent Routes toward Wnt and R-spondin Niche Independency during Human Gastric Carcinogenesis.

Recent sequencing analyses have shed light on heterogeneous patterns of genomic aberrations in human gastric cancers (GCs). To explore how individual genetic events translate into cancer phenotypes, we established a biological library consisting of genetically engineered gastric organoids carrying various GC mutations and 37 patient-derived organoid lines, including rare genomically stable GCs. Phenotype analyses of GC organoids revealed divergent genetic and epigenetic routes to gain Wnt and R-spondin niche independency. An unbiased phenotype-based genetic screening identified a significant association between CDH1/TP53 compound mutations and the R-spondin independency that was functionally validated by CRISPR-based knockout. Xenografting of GC organoids further established the feasibility of Wnt-targeting therapy for Wnt-dependent GCs. Our results collectively demonstrate that multifaceted genetic abnormalities render human GCs independent of the stem cell niche and highlight the validity of the genotype-phenotype screening strategy in gaining deeper understanding of human cancers.

Conceptual and Experimental Tools to Understand Spatial Effects and Transport Phenomena in Nonlinear Biochemical Networks Illustrated with Patchy Switching.

Many biochemical systems are spatially heterogeneous and exhibit nonlinear behaviors, such as state switching in response to small changes in the local concentration of diffusible molecules. Systems as varied as blood clotting, intracellular calcium signaling, and tissue inflammation are all heavily influenced by the balance of rates of reaction and mass transport phenomena including flow and diffusion. Transport of signaling molecules is also affected by geometry and chemoselective confinement via matrix binding. In this review, we use a phenomenon referred to as patchy switching to illustrate the interplay of nonlinearities, transport phenomena, and spatial effects. Patchy switching describes a change in the state of a network when the local concentration of a diffusible molecule surpasses a critical threshold. Using patchy switching as an example, we describe conceptual tools from nonlinear dynamics and chemical engineering that make testable predictions and provide a unifying description of the myriad possible experimental observations. We describe experimental microfluidic and biochemical tools emerging to test conceptual predictions by controlling transport phenomena and spatial distribution of diffusible signals, and we highlight the unmet need for in vivo tools.

Innate Immune Landscape in Early Lung Adenocarcinoma by Paired Single-Cell Analyses.

To guide the design of immunotherapy strategies for patients with early stage lung tumors, we developed a multiscale immune profiling strategy to map the immune landscape of early lung adenocarcinoma lesions to search for tumor-driven immune changes. Utilizing a barcoding method that allows a simultaneous single-cell analysis of the tumor, non-involved lung, and blood cells, we provide a detailed immune cell atlas of early lung tumors. We show that stage I lung adenocarcinoma lesions already harbor significantly altered T cell and NK cell compartments. Moreover, we identified changes in tumor-infiltrating myeloid cell (TIM) subsets that likely compromise anti-tumor T cell immunity. Paired single-cell analyses thus offer valuable knowledge of tumor-driven immune changes, providing a powerful tool for the rational design of immune therapies. VIDEO ABSTRACT.

Enhancer Reprogramming Promotes Pancreatic Cancer Metastasis.

Pancreatic ductal adenocarcinoma (PDA) is one of the most lethal human malignancies, owing in part to its propensity for metastasis. Here, we used an organoid culture system to investigate how transcription and the enhancer landscape become altered during discrete stages of disease progression in a PDA mouse model. This approach revealed that the metastatic transition is accompanied by massive and recurrent alterations in enhancer activity. We implicate the pioneer factor FOXA1 as a driver of enhancer activation in this system, a mechanism that renders PDA cells more invasive and less anchorage-dependent for growth in vitro, as well as more metastatic in vivo. In this context, FOXA1-dependent enhancer reprogramming activates a transcriptional program of embryonic foregut endoderm. Collectively, our study implicates enhancer reprogramming, FOXA1 upregulation, and a retrograde developmental transition in PDA metastasis.

Myc Cooperates with Ras by Programming Inflammation and Immune Suppression.

The two oncogenes KRas and Myc cooperate to drive tumorigenesis, but the mechanism underlying this remains unclear. In a mouse lung model of KRasG12D-driven adenomas, we find that co-activation of Myc drives the immediate transition to highly proliferative and invasive adenocarcinomas marked by highly inflammatory, angiogenic, and immune-suppressed stroma. We identify epithelial-derived signaling molecules CCL9 and IL-23 as the principal instructing signals for stromal reprogramming. CCL9 mediates recruitment of macrophages, angiogenesis, and PD-L1-dependent expulsion of T and B cells. IL-23 orchestrates exclusion of adaptive T and B cells and innate immune NK cells. Co-blockade of both CCL9 and IL-23 abrogates Myc-induced tumor progression. Subsequent deactivation of Myc in established adenocarcinomas triggers immediate reversal of all stromal changes and tumor regression, which are independent of CD4+CD8+ T cells but substantially dependent on returning NK cells. We show that Myc extensively programs an immune suppressive stroma that is obligatory for tumor progression.

TGF-ß Tumor Suppression through a Lethal EMT.

TGF-ß signaling can be pro-tumorigenic or tumor suppressive. We investigated this duality in pancreatic ductal adenocarcinoma (PDA), which, with other gastrointestinal cancers, exhibits frequent inactivation of the TGF-ß mediator Smad4. We show that TGF-ß induces an epithelial-mesenchymal transition (EMT), generally considered a pro-tumorigenic event. However, in TGF-ß-sensitive PDA cells, EMT becomes lethal by converting TGF-ß-induced Sox4 from an enforcer of tumorigenesis into a promoter of apoptosis. This is the result of an EMT-linked remodeling of the cellular transcription factor landscape, including the repression of the gastrointestinal lineage-master regulator Klf5. Klf5 cooperates with Sox4 in oncogenesis and prevents Sox4-induced apoptosis. Smad4 is required for EMT but dispensable for Sox4 induction by TGF-ß. TGF-ß-induced Sox4 is thus geared to bolster progenitor identity, whereas simultaneous Smad4-dependent EMT strips Sox4 of an essential partner in oncogenesis. Our work demonstrates that TGF-ß tumor suppression functions through an EMT-mediated disruption of a lineage-specific transcriptional network.

Mutant KRAS Enhances Tumor Cell Fitness by Upregulating Stress Granules.

There is growing evidence that stress-coping mechanisms represent tumor cell vulnerabilities that may function as therapeutically beneficial targets. Recent work has delineated an integrated stress adaptation mechanism that is characterized by the formation of cytoplasmic mRNA and protein foci, termed stress granules (SGs). Here, we demonstrate that SGs are markedly elevated in mutant KRAS cells following exposure to stress-inducing stimuli. The upregulation of SGs by mutant KRAS is dependent on the production of the signaling lipid molecule 15-deoxy-delta 12,14 prostaglandin J2 (15-d-PGJ2) and confers cytoprotection against stress stimuli and chemotherapeutic agents. The secretion of 15-d-PGJ2 by mutant KRAS cells is sufficient to enhance SG formation and stress resistance in cancer cells that are wild-type for KRAS. Our findings identify a mutant KRAS-dependent cell non-autonomous mechanism that may afford the establishment of a stress-resistant niche that encompasses different tumor subclones. These results should inform the design of strategies to eradicate tumor cell communities.

Vitamin D receptor-mediated stromal reprogramming suppresses pancreatitis and enhances pancreatic cancer therapy.

The poor clinical outcome in pancreatic ductal adenocarcinoma (PDA) is attributed to intrinsic chemoresistance and a growth-permissive tumor microenvironment. Conversion of quiescent to activated pancreatic stellate cells (PSCs) drives the severe stromal reaction that characterizes PDA. Here, we reveal that the vitamin D receptor (VDR) is expressed in stroma from human pancreatic tumors and that treatment with the VDR ligand calcipotriol markedly reduced markers of inflammation and fibrosis in pancreatitis and human tumor stroma. We show that VDR acts as a master transcriptional regulator of PSCs to reprise the quiescent state, resulting in induced stromal remodeling, increased intratumoral gemcitabine, reduced tumor volume, and a 57% increase in survival compared to chemotherapy alone. This work describes a molecular strategy through which transcriptional reprogramming of tumor stroma enables chemotherapeutic response and suggests vitamin D priming as an adjunct in PDA therapy. PAPERFLICK:

Yap1 activation enables bypass of oncogenic Kras addiction in pancreatic cancer.

Activating mutations in KRAS are among the most frequent events in diverse human carcinomas and are particularly prominent in human pancreatic ductal adenocarcinoma (PDAC). An inducible Kras(G12D)-driven mouse model of PDAC has established a critical role for sustained Kras(G12D) expression in tumor maintenance, providing a model to determine the potential for and the underlying mechanisms of Kras(G12D)-independent PDAC recurrence. Here, we show that some tumors undergo spontaneous relapse and are devoid of Kras(G12D) expression and downstream canonical MAPK signaling and instead acquire amplification and overexpression of the transcriptional coactivator Yap1. Functional studies established the role of Yap1 and the transcriptional factor Tead2 in driving Kras(G12D)-independent tumor maintenance. The Yap1/Tead2 complex acts cooperatively with E2F transcription factors to activate a cell cycle and DNA replication program. Our studies, along with corroborating evidence from human PDAC models, portend a novel mechanism of escape from oncogenic Kras addiction in PDAC.

Extrachromosomal DNA in the cancerous transformation of Barrett's oesophagus.

Oncogene amplification on extrachromosomal DNA (ecDNA) drives the evolution of tumours and their resistance to treatment, and is associated with poor outcomes for patients with cancer1-6. At present, it is unclear whether ecDNA is a later manifestation of genomic instability, or whether it can be an early event in the transition from dysplasia to cancer. Here, to better understand the development of ecDNA, we analysed whole-genome sequencing (WGS) data from patients with oesophageal adenocarcinoma (EAC) or Barrett's oesophagus. These data included 206 biopsies in Barrett's oesophagus surveillance and EAC cohorts from Cambridge University. We also analysed WGS and histology data from biopsies that were collected across multiple regions at 2 time points from 80 patients in a case-control study at the Fred Hutchinson Cancer Center. In the Cambridge cohorts, the frequency of ecDNA increased between Barrett's-oesophagus-associated early-stage (24%) and late-stage (43%) EAC, suggesting that ecDNA is formed during cancer progression. In the cohort from the Fred Hutchinson Cancer Center, 33% of patients who developed EAC had at least one oesophageal biopsy with ecDNA before or at the diagnosis of EAC. In biopsies that were collected before cancer diagnosis, higher levels of ecDNA were present in samples from patients who later developed EAC than in samples from those who did not. We found that ecDNAs contained diverse collections of oncogenes and immunomodulatory genes. Furthermore, ecDNAs showed increases in copy number and structural complexity at more advanced stages of disease. Our findings show that ecDNA can develop early in the transition from high-grade dysplasia to cancer, and that ecDNAs progressively form and evolve under positive selection.

Punctuated evolution of prostate cancer genomes.

The analysis of exonic DNA from prostate cancers has identified recurrently mutated genes, but the spectrum of genome-wide alterations has not been profiled extensively in this disease. We sequenced the genomes of 57 prostate tumors and matched normal tissues to characterize somatic alterations and to study how they accumulate during oncogenesis and progression. By modeling the genesis of genomic rearrangements, we identified abundant DNA translocations and deletions that arise in a highly interdependent manner. This phenomenon, which we term "chromoplexy," frequently accounts for the dysregulation of prostate cancer genes and appears to disrupt multiple cancer genes coordinately. Our modeling suggests that chromoplexy may induce considerable genomic derangement over relatively few events in prostate cancer and other neoplasms, supporting a model of punctuated cancer evolution. By characterizing the clonal hierarchy of genomic lesions in prostate tumors, we charted a path of oncogenic events along which chromoplexy may drive prostate carcinogenesis.

Hijacking the neuronal NMDAR signaling circuit to promote tumor growth and invasion.

Glutamate and its receptor N-methyl-D-aspartate receptor (NMDAR) have been associated with cancer, although their functions are not fully understood. Herein, we implicate glutamate-driven NMDAR signaling in a mouse model of pancreatic neuroendocrine tumorigenesis (PNET) and in selected human cancers. NMDAR was upregulated at the periphery of PNET tumors, particularly invasive fronts. Moreover, elevated coexpression of NMDAR and glutamate exporters correlated with poor prognosis in cancer patients. Treatment of a tumor-derived cell line with NMDAR antagonists impaired cancer cell proliferation and invasion. Flow conditions mimicking interstitial fluid pressure induced autologous glutamate secretion, activating NMDAR and its downstream MEK-MAPK and CaMK effectors, thereby promoting invasiveness. Congruently, pharmacological inhibition of NMDAR in mice with PNET reduced tumor growth and invasiveness. Therefore, beyond its traditional role in neurons, NMDAR may be activated in human tumors by fluid flow consequent to higher interstitial pressure, inducing an autocrine glutamate signaling circuit with resultant stimulation of malignancy.

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.

Tumor Cell-Intrinsic Factors Underlie Heterogeneity of Immune Cell Infiltration and Response to Immunotherapy.

The biological and functional heterogeneity between tumors-both across and within cancer types-poses a challenge for immunotherapy. To understand the factors underlying tumor immune heterogeneity and immunotherapy sensitivity, we established a library of congenic tumor cell clones from an autochthonous mouse model of pancreatic adenocarcinoma. These clones generated tumors that recapitulated T cell-inflamed and non-T-cell-inflamed tumor microenvironments upon implantation in immunocompetent mice, with distinct patterns of infiltration by immune cell subsets. Co-injecting tumor cell clones revealed the non-T-cell-inflamed phenotype is dominant and that both quantitative and qualitative features of intratumoral CD8+ T cells determine response to therapy. Transcriptomic and epigenetic analyses revealed tumor-cell-intrinsic production of the chemokine CXCL1 as a determinant of the non-T-cell-inflamed microenvironment, and ablation of CXCL1 promoted T cell infiltration and sensitivity to a combination immunotherapy regimen. Thus, tumor cell-intrinsic factors shape the tumor immune microenvironment and influence the outcome of immunotherapy.

Mapping the hallmarks of lung adenocarcinoma with massively parallel sequencing.

Lung adenocarcinoma, the most common subtype of non-small cell lung cancer, is responsible for more than 500,000 deaths per year worldwide. Here, we report exome and genome sequences of 183 lung adenocarcinoma tumor/normal DNA pairs. These analyses revealed a mean exonic somatic mutation rate of 12.0 events/megabase and identified the majority of genes previously reported as significantly mutated in lung adenocarcinoma. In addition, we identified statistically recurrent somatic mutations in the splicing factor gene U2AF1 and truncating mutations affecting RBM10 and ARID1A. Analysis of nucleotide context-specific mutation signatures grouped the sample set into distinct clusters that correlated with smoking history and alterations of reported lung adenocarcinoma genes. Whole-genome sequence analysis revealed frequent structural rearrangements, including in-frame exonic alterations within EGFR and SIK2 kinases. The candidate genes identified in this study are attractive targets for biological characterization and therapeutic targeting of lung adenocarcinoma.