A modular master regulator landscape controls cancer transcriptional identity.

Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.

Proteogenomic characterization of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

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.

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.

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.

SIRT6 Suppresses Pancreatic Cancer through Control of Lin28b.

Chromatin remodeling proteins are frequently dysregulated in human cancer, yet little is known about how they control tumorigenesis. Here, we uncover an epigenetic program mediated by the NAD(+)-dependent histone deacetylase Sirtuin 6 (SIRT6) that is critical for suppression of pancreatic ductal adenocarcinoma (PDAC), one of the most lethal malignancies. SIRT6 inactivation accelerates PDAC progression and metastasis via upregulation of Lin28b, a negative regulator of the let-7 microRNA. SIRT6 loss results in histone hyperacetylation at the Lin28b promoter, Myc recruitment, and pronounced induction of Lin28b and downstream let-7 target genes, HMGA2, IGF2BP1, and IGF2BP3. This epigenetic program defines a distinct subset with a poor prognosis, representing 30%-40% of human PDAC, characterized by reduced SIRT6 expression and an exquisite dependence on Lin28b for tumor growth. Thus, we identify SIRT6 as an important PDAC tumor suppressor and uncover the Lin28b pathway as a potential therapeutic target in a molecularly defined PDAC subset. PAPERCLIP.

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.

SMYD5 methylation of rpL40 links ribosomal output to gastric cancer.

Dysregulated transcription due to disruption in histone lysine methylation dynamics is an established contributor to tumorigenesis1,2. However, whether analogous pathologic epigenetic mechanisms act directly on the ribosome to advance oncogenesis is unclear. Here we find that trimethylation of the core ribosomal protein L40 (rpL40) at lysine 22 (rpL40K22me3) by the lysine methyltransferase SMYD5 regulates mRNA translation output to promote malignant progression of gastric adenocarcinoma (GAC) with lethal peritoneal ascites. A biochemical-proteomics strategy identifies the monoubiquitin fusion protein partner rpL40 (ref. 3) as the principal physiological substrate of SMYD5 across diverse samples. Inhibiting the SMYD5-rpL40K22me3 axis in GAC cell lines reprogrammes protein synthesis to attenuate oncogenic gene expression signatures. SMYD5 and rpL40K22me3 are upregulated in samples from patients with GAC and negatively correlate with clinical outcomes. SMYD5 ablation in vivo in familial and sporadic mouse models of malignant GAC blocks metastatic disease, including peritoneal carcinomatosis. Suppressing SMYD5 methylation of rpL40 inhibits human cancer cell and patient-derived GAC xenograft growth and renders them hypersensitive to inhibitors of PI3K and mTOR. Finally, combining SMYD5 depletion with PI3K-mTOR inhibition and chimeric antigen receptor T cell administration cures an otherwise lethal in vivo mouse model of aggressive GAC-derived peritoneal carcinomatosis. Together, our work uncovers a ribosome-based epigenetic mechanism that facilitates the evolution of malignant GAC and proposes SMYD5 targeting as part of a potential combination therapy to treat this cancer.

CRISPR-Cas9 knockin mice for genome editing and cancer modeling.

CRISPR-Cas9 is a versatile genome editing technology for studying the functions of genetic elements. To broadly enable the application of Cas9 in vivo, we established a Cre-dependent Cas9 knockin mouse. We demonstrated in vivo as well as ex vivo genome editing using adeno-associated virus (AAV)-, lentivirus-, or particle-mediated delivery of guide RNA in neurons, immune cells, and endothelial cells. Using these mice, we simultaneously modeled the dynamics of KRAS, p53, and LKB1, the top three significantly mutated genes in lung adenocarcinoma. Delivery of a single AAV vector in the lung generated loss-of-function mutations in p53 and Lkb1, as well as homology-directed repair-mediated Kras(G12D) mutations, leading to macroscopic tumors of adenocarcinoma pathology. Together, these results suggest that Cas9 mice empower a wide range of biological and disease modeling applications.

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.

Current treatment and recent progress in gastric cancer.

Gastric cancer is not a top-10 malignancy in the United States but represents one of the most common causes of cancer death worldwide. Biological differences between tumors from Eastern and Western countries add to the complexity of identifying standard-of-care therapy based on international trials. Systemic chemotherapy, radiotherapy, surgery, immunotherapy, and targeted therapy all have proven efficacy in gastric adenocarcinoma; therefore, multidisciplinary treatment is paramount to treatment selection. Triplet chemotherapy for resectable gastric cancer is now accepted and could represent a plateau of standard cytotoxic chemotherapy for localized disease. Classification of gastric cancer based on molecular subtypes is providing an opportunity for personalized therapy. Biomarkers, in particular microsatellite instability (MSI), programmed cell death ligand 1 (PD-L1), human epidermal growth factor receptor 2 (HER2), tumor mutation burden, and Epstein-Barr virus, are increasingly driving systemic therapy approaches and allowing for the identification of populations most likely to benefit from immunotherapy and targeted therapy. Significant research opportunities remain for the less differentiated histologic subtypes of gastric adenocarcinoma and those without markers of immunotherapy activity.

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.

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.

Zmat3 Is a Key Splicing Regulator in the p53 Tumor Suppression Program.

Although TP53 is the most commonly mutated gene in human cancers, the p53-dependent transcriptional programs mediating tumor suppression remain incompletely understood. Here, to uncover critical components downstream of p53 in tumor suppression, we perform unbiased RNAi and CRISPR-Cas9-based genetic screens in vivo. These screens converge upon the p53-inducible gene Zmat3, encoding an RNA-binding protein, and we demonstrate that ZMAT3 is an important tumor suppressor downstream of p53 in mouse KrasG12D-driven lung and liver cancers and human carcinomas. Integrative analysis of the ZMAT3 RNA-binding landscape and transcriptomic profiling reveals that ZMAT3 directly modulates exon inclusion in transcripts encoding proteins of diverse functions, including the p53 inhibitors MDM4 and MDM2, splicing regulators, and components of varied cellular processes. Interestingly, these exons are enriched in NMD signals, and, accordingly, ZMAT3 broadly affects target transcript stability. Collectively, these studies reveal ZMAT3 as a novel RNA-splicing and homeostasis regulator and a key component of p53-mediated tumor suppression.

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.

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.

Single-cell and single-nucleus RNA-sequencing from paired normal-adenocarcinoma lung samples provide both common and discordant biological insights.

Whether single-cell RNA-sequencing (scRNA-seq) captures the same biological information as single-nucleus RNA-sequencing (snRNA-seq) remains uncertain and likely to be context-dependent. Herein, a head-to-head comparison was performed in matched normal-adenocarcinoma human lung samples to assess biological insights derived from scRNA-seq versus snRNA-seq and better understand the cellular transition that occurs from normal to tumoral tissue. Here, the transcriptome of 160,621 cells/nuclei was obtained. In non-tumor lung, cell type proportions varied widely between scRNA-seq and snRNA-seq with a predominance of immune cells in the former (81.5%) and epithelial cells (69.9%) in the later. Similar results were observed in adenocarcinomas, in addition to an overall increase in cell type heterogeneity and a greater prevalence of copy number variants in cells of epithelial origin, which suggests malignant assignment. The cell type transition that occurs from normal lung tissue to adenocarcinoma was not always concordant whether cells or nuclei were examined. As expected, large differential expression of the whole-cell and nuclear transcriptome was observed, but cell-type specific changes of paired normal and tumor lung samples revealed a set of common genes in the cells and nuclei involved in cancer-related pathways. In addition, we showed that the ligand-receptor interactome landscape of lung adenocarcinoma was largely different whether cells or nuclei were evaluated. Immune cell depletion in fresh specimens partly mitigated the difference in cell type composition observed between cells and nuclei. However, the extra manipulations affected cell viability and amplified the transcriptional signatures associated with stress responses. In conclusion, research applications focussing on mapping the immune landscape of lung adenocarcinoma benefit from scRNA-seq in fresh samples, whereas snRNA-seq of frozen samples provide a low-cost alternative to profile more epithelial and cancer cells, and yield cell type proportions that more closely match tissue content.