Transcriptomic analysis of cutaneous squamous cell carcinoma reveals a multigene prognostic signature associated with metastasis.

BACKGROUND: Metastasis of cutaneous squamous cell carcinoma (cSCC) is uncommon. Current staging methods are reported to have sub-optimal performances in metastasis prediction. Accurate identification of patients with tumors at high risk of metastasis would have a significant impact on management. OBJECTIVE: To develop a robust and validated gene expression profile signature for predicting primary cSCC metastatic risk using an unbiased whole transcriptome discovery-driven approach. METHODS: Archival formalin-fixed paraffin-embedded primary cSCC with perilesional normal tissue from 237 immunocompetent patients (151 nonmetastasizing and 86 metastasizing) were collected retrospectively from four centers. TempO-seq was used to probe the whole transcriptome and machine learning algorithms were applied to derive predictive signatures, with a 3:1 split for training and testing datasets. RESULTS: A 20-gene prognostic model was developed and validated, with an accuracy of 86.0%, sensitivity of 85.7%, specificity of 86.1%, and positive predictive value of 78.3% in the testing set, providing more stable, accurate prediction than pathological staging systems. A linear predictor was also developed, significantly correlating with metastatic risk. LIMITATIONS: This was a retrospective 4-center study and larger prospective multicenter studies are now required. CONCLUSION: The 20-gene signature prediction is accurate, with the potential to be incorporated into clinical workflows for cSCC.

Whole-genome landscape of pancreatic neuroendocrine tumours.

The diagnosis of pancreatic neuroendocrine tumours (PanNETs) is increasing owing to more sensitive detection methods, and this increase is creating challenges for clinical management. We performed whole-genome sequencing of 102 primary PanNETs and defined the genomic events that characterize their pathogenesis. Here we describe the mutational signatures they harbour, including a deficiency in G:C > T:A base excision repair due to inactivation of MUTYH, which encodes a DNA glycosylase. Clinically sporadic PanNETs contain a larger-than-expected proportion of germline mutations, including previously unreported mutations in the DNA repair genes MUTYH, CHEK2 and BRCA2. Together with mutations in MEN1 and VHL, these mutations occur in 17% of patients. Somatic mutations, including point mutations and gene fusions, were commonly found in genes involved in four main pathways: chromatin remodelling, DNA damage repair, activation of mTOR signalling (including previously undescribed EWSR1 gene fusions), and telomere maintenance. In addition, our gene expression analyses identified a subgroup of tumours associated with hypoxia and HIF signalling.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours.

This corrects the article DOI: 10.1038/nature21063.

Targeting DNA Damage Response and Replication Stress in Pancreatic Cancer.

BACKGROUND & AIMS: Continuing recalcitrance to therapy cements pancreatic cancer (PC) as the most lethal malignancy, which is set to become the second leading cause of cancer death in our society. The study aim was to investigate the association between DNA damage response (DDR), replication stress, and novel therapeutic response in PC to develop a biomarker-driven therapeutic strategy targeting DDR and replication stress in PC. METHODS: We interrogated the transcriptome, genome, proteome, and functional characteristics of 61 novel PC patient-derived cell lines to define novel therapeutic strategies targeting DDR and replication stress. Validation was done in patient-derived xenografts and human PC organoids. RESULTS: Patient-derived cell lines faithfully recapitulate the epithelial component of pancreatic tumors, including previously described molecular subtypes. Biomarkers of DDR deficiency, including a novel signature of homologous recombination deficiency, cosegregates with response to platinum (P < .001) and PARP inhibitor therapy (P < .001) in vitro and in vivo. We generated a novel signature of replication stress that predicts response to ATR (P < .018) and WEE1 inhibitor (P < .029) treatment in both cell lines and human PC organoids. Replication stress was enriched in the squamous subtype of PC (P < .001) but was not associated with DDR deficiency. CONCLUSIONS: Replication stress and DDR deficiency are independent of each other, creating opportunities for therapy in DDR-proficient PC and after platinum therapy.

Dual targeting of p53 and c-MYC selectively eliminates leukaemic stem cells.

Chronic myeloid leukaemia (CML) arises after transformation of a haemopoietic stem cell (HSC) by the protein-tyrosine kinase BCR-ABL. Direct inhibition of BCR-ABL kinase has revolutionized disease management, but fails to eradicate leukaemic stem cells (LSCs), which maintain CML. LSCs are independent of BCR-ABL for survival, providing a rationale for identifying and targeting kinase-independent pathways. Here we show--using proteomics, transcriptomics and network analyses--that in human LSCs, aberrantly expressed proteins, in both imatinib-responder and non-responder patients, are modulated in concert with p53 (also known as TP53) and c-MYC regulation. Perturbation of both p53 and c-MYC, and not BCR-ABL itself, leads to synergistic cell kill, differentiation, and near elimination of transplantable human LSCs in mice, while sparing normal HSCs. This unbiased systems approach targeting connected nodes exemplifies a novel precision medicine strategy providing evidence that LSCs can be eradicated.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Whole-genome characterization of chemoresistant ovarian cancer.

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.

Whole genomes redefine the mutational landscape of pancreatic cancer.

Pancreatic cancer remains one of the most lethal of malignancies and a major health burden. We performed whole-genome sequencing and copy number variation (CNV) analysis of 100 pancreatic ductal adenocarcinomas (PDACs). Chromosomal rearrangements leading to gene disruption were prevalent, affecting genes known to be important in pancreatic cancer (TP53, SMAD4, CDKN2A, ARID1A and ROBO2) and new candidate drivers of pancreatic carcinogenesis (KDM6A and PREX2). Patterns of structural variation (variation in chromosomal structure) classified PDACs into 4 subtypes with potential clinical utility: the subtypes were termed stable, locally rearranged, scattered and unstable. A significant proportion harboured focal amplifications, many of which contained druggable oncogenes (ERBB2, MET, FGFR1, CDK6, PIK3R3 and PIK3CA), but at low individual patient prevalence. Genomic instability co-segregated with inactivation of DNA maintenance genes (BRCA1, BRCA2 or PALB2) and a mutational signature of DNA damage repair deficiency. Of 8 patients who received platinum therapy, 4 of 5 individuals with these measures of defective DNA maintenance responded.

Hypermutation In Pancreatic Cancer.

Pancreatic cancer is molecularly diverse, with few effective therapies. Increased mutation burden and defective DNA repair are associated with response to immune checkpoint inhibitors in several other cancer types. We interrogated 385 pancreatic cancer genomes to define hypermutation and its causes. Mutational signatures inferring defects in DNA repair were enriched in those with the highest mutation burdens. Mismatch repair deficiency was identified in 1% of tumors harboring different mechanisms of somatic inactivation of MLH1 and MSH2. Defining mutation load in individual pancreatic cancers and the optimal assay for patient selection may inform clinical trial design for immunotherapy in pancreatic cancer.

Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Integrated genomic and transcriptomic analysis of human brain metastases identifies alterations of potential clinical significance.

Treatment options for patients with brain metastases (BMs) have limited efficacy and the mortality rate is virtually 100%. Targeted therapy is critically under-utilized, and our understanding of mechanisms underpinning metastatic outgrowth in the brain is limited. To address these deficiencies, we investigated the genomic and transcriptomic landscapes of 36 BMs from breast, lung, melanoma and oesophageal cancers, using DNA copy-number analysis and exome- and RNA-sequencing. The key findings were as follows. (a) Identification of novel candidates with possible roles in BM development, including the significantly mutated genes DSC2, ST7, PIK3R1 and SMC5, and the DNA repair, ERBB-HER signalling, axon guidance and protein kinase-A signalling pathways. (b) Mutational signature analysis was applied to successfully identify the primary cancer type for two BMs with unknown origins. (c) Actionable genomic alterations were identified in 31/36 BMs (86%); in one case we retrospectively identified ERBB2 amplification representing apparent HER2 status conversion, then confirmed progressive enrichment for HER2-positivity across four consecutive metastatic deposits by IHC and SISH, resulting in the deployment of HER2-targeted therapy for the patient. (d) In the ERBB/HER pathway, ERBB2 expression correlated with ERBB3 (r(2) = 0.496; p < 0.0001) and HER3 and HER4 were frequently activated in an independent cohort of 167 archival BM from seven primary cancer types: 57.6% and 52.6% of cases were phospho-HER3(Y1222) or phospho-HER4(Y1162) membrane-positive, respectively. The HER3 ligands NRG1/2 were barely detectable by RNAseq, with NRG1 (8p12) genomic loss in 63.6% breast cancer-BMs, suggesting a microenvironmental source of ligand. In summary, this is the first study to characterize the genomic landscapes of BM. The data revealed novel candidates, potential clinical applications for genomic profiling of resectable BMs, and highlighted the possibility of therapeutically targeting HER3, which is broadly over-expressed and activated in BMs, independent of primary site and systemic therapy.

Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling.

The importance of epigenetic modifications such as DNA methylation in tumorigenesis is increasingly being appreciated. To define the genome-wide pattern of DNA methylation in pancreatic ductal adenocarcinomas (PDAC), we captured the methylation profiles of 167 untreated resected PDACs and compared them to a panel of 29 adjacent nontransformed pancreata using high-density arrays. A total of 11,634 CpG sites associated with 3,522 genes were significantly differentially methylated (DM) in PDAC and were capable of segregating PDAC from non-malignant pancreas, regardless of tumor cellularity. As expected, PDAC hypermethylation was most prevalent in the 5' region of genes (including the proximal promoter, 5'UTR and CpG islands). Approximately 33% DM genes showed significant inverse correlation with mRNA expression levels. Pathway analysis revealed an enrichment of aberrantly methylated genes involved in key molecular mechanisms important to PDAC: TGF-ß, WNT, integrin signaling, cell adhesion, stellate cell activation and axon guidance. Given the recent discovery that SLIT-ROBO mutations play a clinically important role in PDAC, the role of epigenetic perturbation of axon guidance was pursued in more detail. Bisulfite amplicon deep sequencing and qRT-PCR expression analyses confirmed recurrent perturbation of axon guidance pathway genes SLIT2, SLIT3, ROBO1, ROBO3, ITGA2 and MET and suggests epigenetic suppression of SLIT-ROBO signaling and up-regulation of MET and ITGA2 expression. Hypomethylation of MET and ITGA2 correlated with high gene expression, which was associated with poor survival. These data suggest that aberrant methylation plays an important role in pancreatic carcinogenesis affecting core signaling pathways with potential implications for the disease pathophysiology and therapy.

Efferocytosis reprograms the tumor microenvironment to promote pancreatic cancer liver metastasis.

Pancreatic ductal adenocarcinoma is a highly metastatic disease and macrophages support liver metastases. Efferocytosis, or engulfment of apoptotic cells by macrophages, is an essential process in tissue homeostasis and wound healing, but its role in metastasis is less well understood. Here, we found that the colonization of the hepatic metastatic site is accompanied by low-grade tissue injury and that efferocytosis-mediated clearance of parenchymal dead cells promotes macrophage reprogramming and liver metastasis. Mechanistically, progranulin expression in macrophages is necessary for efficient efferocytosis by controlling lysosomal acidification via cystic fibrosis transmembrane conductance regulator and the degradation of lysosomal cargo, resulting in LXRα/RXRα-mediated macrophage conversion and upregulation of arginase 1. Pharmacological blockade of efferocytosis or macrophage-specific genetic depletion of progranulin impairs macrophage conversion, improves CD8+ T cell functions, and reduces liver metastasis. Our findings reveal how hard-wired functions of macrophages in tissue repair contribute to liver metastasis and identify potential targets for prevention of pancreatic ductal adenocarcinoma liver metastasis.

Unraveling Heterogeneity in the Aging Hematopoietic Stem Cell Compartment: An Insight From Single-cell Approaches.

Specific cell types and, therefore, organs respond differently during aging. This is also true for the hematopoietic system, where it has been demonstrated that hematopoietic stem cells alter a variety of features, such as their metabolism, and accumulate DNA damage, which can lead to clonal outgrowth over time. In addition, profound changes in the bone marrow microenvironment upon aging lead to senescence in certain cell types such as mesenchymal stem cells and result in increased inflammation. This heterogeneity makes it difficult to pinpoint the molecular drivers of organismal aging gained from bulk approaches, such as RNA sequencing. A better understanding of the heterogeneity underlying the aging process in the hematopoietic compartment is, therefore, needed. With the advances of single-cell technologies in recent years, it is now possible to address fundamental questions of aging. In this review, we discuss how single-cell approaches can and indeed are already being used to understand changes observed during aging in the hematopoietic compartment. We will touch on established and novel methods for flow cytometric detection, single-cell culture approaches, and single-cell omics.

Metabolic profiling stratifies colorectal cancer and reveals adenosylhomocysteinase as a therapeutic target.

The genomic landscape of colorectal cancer (CRC) is shaped by inactivating mutations in tumour suppressors such as APC, and oncogenic mutations such as mutant KRAS. Here we used genetically engineered mouse models, and multimodal mass spectrometry-based metabolomics to study the impact of common genetic drivers of CRC on the metabolic landscape of the intestine. We show that untargeted metabolic profiling can be applied to stratify intestinal tissues according to underlying genetic alterations, and use mass spectrometry imaging to identify tumour, stromal and normal adjacent tissues. By identifying ions that drive variation between normal and transformed tissues, we found dysregulation of the methionine cycle to be a hallmark of APC-deficient CRC. Loss of Apc in the mouse intestine was found to be sufficient to drive expression of one of its enzymes, adenosylhomocysteinase (AHCY), which was also found to be transcriptionally upregulated in human CRC. Targeting of AHCY function impaired growth of APC-deficient organoids in vitro, and prevented the characteristic hyperproliferative/crypt progenitor phenotype driven by acute deletion of Apc in vivo, even in the context of mutant Kras. Finally, pharmacological inhibition of AHCY reduced intestinal tumour burden in ApcMin/+ mice indicating its potential as a metabolic drug target in CRC.

HNF4A and GATA6 Loss Reveals Therapeutically Actionable Subtypes in Pancreatic Cancer.

Pancreatic ductal adenocarcinoma (PDAC) can be divided into transcriptomic subtypes with two broad lineages referred to as classical (pancreatic) and squamous. We find that these two subtypes are driven by distinct metabolic phenotypes. Loss of genes that drive endodermal lineage specification, HNF4A and GATA6, switch metabolic profiles from classical (pancreatic) to predominantly squamous, with glycogen synthase kinase 3 beta (GSK3ß) a key regulator of glycolysis. Pharmacological inhibition of GSK3ß results in selective sensitivity in the squamous subtype; however, a subset of these squamous patient-derived cell lines (PDCLs) acquires rapid drug tolerance. Using chromatin accessibility maps, we demonstrate that the squamous subtype can be further classified using chromatin accessibility to predict responsiveness and tolerance to GSK3ß inhibitors. Our findings demonstrate that distinct patterns of chromatin accessibility can be used to identify patient subgroups that are indistinguishable by gene expression profiles, highlighting the utility of chromatin-based biomarkers for patient selection in the treatment of PDAC.

Epithelial NOTCH Signaling Rewires the Tumor Microenvironment of Colorectal Cancer to Drive Poor-Prognosis Subtypes and Metastasis.

The metastatic process of colorectal cancer (CRC) is not fully understood and effective therapies are lacking. We show that activation of NOTCH1 signaling in the murine intestinal epithelium leads to highly penetrant metastasis (100% metastasis; with >80% liver metastases) in KrasG12D-driven serrated cancer. Transcriptional profiling reveals that epithelial NOTCH1 signaling creates a tumor microenvironment (TME) reminiscent of poorly prognostic human CRC subtypes (CMS4 and CRIS-B), and drives metastasis through transforming growth factor (TGF) ß-dependent neutrophil recruitment. Importantly, inhibition of this recruitment with clinically relevant therapeutic agents blocks metastasis. We propose that NOTCH1 signaling is key to CRC progression and should be exploited clinically.

Copy number signatures and mutational processes in ovarian carcinoma.

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.