Exploring the interplay between colorectal cancer subtypes genomic variants and cellular morphology: A deep-learning approach.

Molecular subtypes of colorectal cancer (CRC) significantly influence treatment decisions. While convolutional neural networks (CNNs) have recently been introduced for automated CRC subtype identification using H&E stained histopathological images, the correlation between CRC subtype genomic variants and their corresponding cellular morphology expressed by their imaging phenotypes is yet to be fully explored. The goal of this study was to determine such correlations by incorporating genomic variants in CNN models for CRC subtype classification from H&E images. We utilized the publicly available TCGA-CRC-DX dataset, which comprises whole slide images from 360 CRC-diagnosed patients (260 for training and 100 for testing). This dataset also provides information on CRC subtype classifications and genomic variations. We trained CNN models for CRC subtype classification that account for potential correlation between genomic variations within CRC subtypes and their corresponding cellular morphology patterns. We assessed the interplay between CRC subtypes' genomic variations and cellular morphology patterns by evaluating the CRC subtype classification accuracy of the different models in a stratified 5-fold cross-validation experimental setup using the area under the ROC curve (AUROC) and average precision (AP) as the performance metrics. The CNN models that account for potential correlation between genomic variations within CRC subtypes and their cellular morphology pattern achieved superior accuracy compared to the baseline CNN classification model that does not account for genomic variations when using either single-nucleotide-polymorphism (SNP) molecular features (AUROC: 0.824±0.02 vs. 0.761±0.04, p<0.05, AP: 0.652±0.06 vs. 0.58±0.08) or CpG-Island methylation phenotype (CIMP) molecular features (AUROC: 0.834±0.01 vs. 0.787±0.03, p<0.05, AP: 0.687±0.02 vs. 0.64±0.05). Combining the CNN models account for variations in CIMP and SNP further improved classification accuracy (AUROC: 0.847±0.01 vs. 0.787±0.03, p = 0.01, AP: 0.68±0.02 vs. 0.64±0.05). The improved accuracy of CNN models for CRC subtype classification that account for potential correlation between genomic variations within CRC subtypes and their corresponding cellular morphology as expressed by H&E imaging phenotypes may elucidate the biological cues impacting cancer histopathological imaging phenotypes. Moreover, considering CRC subtypes genomic variations has the potential to improve the accuracy of deep-learning models in discerning cancer subtype from histopathological imaging data.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Therapy-induced APOBEC3A drives evolution of persistent cancer cells.

Acquired drug resistance to anticancer targeted therapies remains an unsolved clinical problem. Although many drivers of acquired drug resistance have been identified1-4, the underlying molecular mechanisms shaping tumour evolution during treatment are incompletely understood. Genomic profiling of patient tumours has implicated apolipoprotein B messenger RNA editing catalytic polypeptide-like (APOBEC) cytidine deaminases in tumour evolution; however, their role during therapy and the development of acquired drug resistance is undefined. Here we report that lung cancer targeted therapies commonly used in the clinic can induce cytidine deaminase APOBEC3A (A3A), leading to sustained mutagenesis in drug-tolerant cancer cells persisting during therapy. Therapy-induced A3A promotes the formation of double-strand DNA breaks, increasing genomic instability in drug-tolerant persisters. Deletion of A3A reduces APOBEC mutations and structural variations in persister cells and delays the development of drug resistance. APOBEC mutational signatures are enriched in tumours from patients with lung cancer who progressed after extended responses to targeted therapies. This study shows that induction of A3A in response to targeted therapies drives evolution of drug-tolerant persister cells, suggesting that suppression of A3A expression or activity may represent a potential therapeutic strategy in the prevention or delay of acquired resistance to lung cancer targeted therapy.

Efficacy of Nivolumab in Pediatric Cancers with High Mutation Burden and Mismatch Repair Deficiency.

PURPOSE: Checkpoint inhibitors have limited efficacy for children with unselected solid and brain tumors. We report the first prospective pediatric trial (NCT02992964) using nivolumab exclusively for refractory nonhematologic cancers harboring tumor mutation burden (TMB) ≥5 mutations/megabase (mut/Mb) and/or mismatch repair deficiency (MMRD). PATIENTS AND METHODS: Twenty patients were screened, and 10 were ultimately included in the response cohort of whom nine had TMB >10 mut/Mb (three initially eligible based on MMRD) and one patient had TMB between 5 and 10 mut/Mb. RESULTS: Delayed immune responses contributed to best overall response of 50%, improving on initial objective responses (20%) and leading to 2-year overall survival (OS) of 50% [95% confidence interval (CI), 27-93]. Four children, including three with refractory malignant gliomas are in complete remission at a median follow-up of 37 months (range, 32.4-60), culminating in 2-year OS of 43% (95% CI, 18.2-100). Biomarker analyses confirmed benefit in children with germline MMRD, microsatellite instability, higher activated and lower regulatory circulating T cells. Stochastic mutation accumulation driven by underlying germline MMRD impacted the tumor microenvironment, contributing to delayed responses. No benefit was observed in the single patient with an MMR-proficient tumor and TMB 7.4 mut/Mb. CONCLUSIONS: Nivolumab resulted in durable responses and prolonged survival for the first time in a pediatric trial of refractory hypermutated cancers including malignant gliomas. Novel biomarkers identified here need to be translated rapidly to clinical care to identify children who can benefit from checkpoint inhibitors, including upfront management of cancer. See related commentary by Mardis, p. 4701.

Genomic Microsatellite Signatures Identify Germline Mismatch Repair Deficiency and Risk of Cancer Onset.

PURPOSE: Diagnosis of Mismatch Repair Deficiency (MMRD) is crucial for tumor management and early detection in patients with the cancer predisposition syndrome constitutional mismatch repair deficiency (CMMRD). Current diagnostic tools are cumbersome and inconsistent both in childhood cancers and in determining germline MMRD. PATIENTS AND METHODS: We developed and analyzed a functional Low-pass Genomic Instability Characterization (LOGIC) assay to detect MMRD. The diagnostic performance of LOGIC was compared with that of current established assays including tumor mutational burden, immunohistochemistry, and the microsatellite instability panel. LOGIC was then applied to various normal tissues of patients with CMMRD with comprehensive clinical data including age of cancer presentation. RESULTS: Overall, LOGIC was 100% sensitive and specific in detecting MMRD in childhood cancers (N = 376). It was more sensitive than the microsatellite instability panel (14%, P = 4.3 × 10-12), immunohistochemistry (86%, P = 4.6 × 10-3), or tumor mutational burden (80%, P = 9.1 × 10-4). LOGIC was able to distinguish CMMRD from other cancer predisposition syndromes using blood and saliva DNA (P < .0001, n = 277). In normal cells, MMRDness scores differed between tissues (GI > blood > brain), increased over time in the same individual, and revealed genotype-phenotype associations within the mismatch repair genes. Importantly, increased MMRDness score was associated with younger age of first cancer presentation in individuals with CMMRD (P = 2.2 × 10-5). CONCLUSION: LOGIC was a robust tool for the diagnosis of MMRD in multiple cancer types and in normal tissues. LOGIC may inform therapeutic cancer decisions, provide rapid diagnosis of germline MMRD, and support tailored surveillance for individuals with CMMRD.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century.

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Di-genic inheritance of germline POLE and PMS2 pathogenic variants causes a unique condition associated with pediatric cancer predisposition.

Polymerase proofreading-associated polyposis (PPAP) and Lynch syndrome, caused by mutated POLE and mismatch repair (MMR) genes, respectively, are associated with adult-onset cancer. PPAP and MMR-deficient tumors are both hypermutated, and each has a unique mutational signature. We describe a 4.5-year-old boy with multiple café au lait spots who presented with metastatic Sonic Hedgehog-activated medulloblastoma, with partial response to intensive chemotherapy and immunotherapy. The tumor showed microsatellite stability, loss of PMS2 nuclear expression, and an exceptionally high tumor mutational burden of 276 Mut/Mb. Germline molecular analysis revealed an inherited heterozygous pathogenic POLE variant and a de novo heterozygous PMS2 pathogenic variant. The tumor featured the MMR, POLE, and POLE+MMR mutational signatures. This is the first description of a di-genic condition, which we named "POL-LYNCH syndrome," manifested by an aggressive ultra-mutant pediatric medulloblastoma with a unique genomic signature.

Radiation-related genomic profile of papillary thyroid carcinoma after the Chernobyl accident.

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid carcinoma (PTC) incidence in surrounding regions, particularly for radioactive iodine (131I)-exposed children. We analyzed genomic, transcriptomic, and epigenomic characteristics of 440 PTCs from Ukraine (from 359 individuals with estimated childhood 131I exposure and 81 unexposed children born after 1986). PTCs displayed radiation dose-dependent enrichment of fusion drivers, nearly all in the mitogen-activated protein kinase pathway, and increases in small deletions and simple/balanced structural variants that were clonal and bore hallmarks of nonhomologous end-joining repair. Radiation-related genomic alterations were more pronounced for individuals who were younger at exposure. Transcriptomic and epigenomic features were strongly associated with driver events but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth after environmental radiation exposure.

DNA Polymerase and Mismatch Repair Exert Distinct Microsatellite Instability Signatures in Normal and Malignant Human Cells.

Although replication repair deficiency, either by mismatch repair deficiency (MMRD) and/or loss of DNA polymerase proofreading, can cause hypermutation in cancer, microsatellite instability (MSI) is considered a hallmark of MMRD alone. By genome-wide analysis of tumors with germline and somatic deficiencies in replication repair, we reveal a novel association between loss of polymerase proofreading and MSI, especially when both components are lost. Analysis of indels in microsatellites (MS-indels) identified five distinct signatures (MS-sigs). MMRD MS-sigs are dominated by multibase losses, whereas mutant-polymerase MS-sigs contain primarily single-base gains. MS deletions in MMRD tumors depend on the original size of the MS and converge to a preferred length, providing mechanistic insight. Finally, we demonstrate that MS-sigs can be a powerful clinical tool for managing individuals with germline MMRD and replication repair-deficient cancers, as they can detect the replication repair deficiency in normal cells and predict their response to immunotherapy. SIGNIFICANCE: Exome- and genome-wide MSI analysis reveals novel signatures that are uniquely attributed to mismatch repair and DNA polymerase. This provides new mechanistic insight into MS maintenance and can be applied clinically for diagnosis of replication repair deficiency and immunotherapy response prediction.This article is highlighted in the In This Issue feature, p. 995.

Proteogenomic Landscape of Breast Cancer Tumorigenesis and Targeted Therapy.

The integration of mass spectrometry-based proteomics with next-generation DNA and RNA sequencing profiles tumors more comprehensively. Here this "proteogenomics" approach was applied to 122 treatment-naive primary breast cancers accrued to preserve post-translational modifications, including protein phosphorylation and acetylation. Proteogenomics challenged standard breast cancer diagnoses, provided detailed analysis of the ERBB2 amplicon, defined tumor subsets that could benefit from immune checkpoint therapy, and allowed more accurate assessment of Rb status for prediction of CDK4/6 inhibitor responsiveness. Phosphoproteomics profiles uncovered novel associations between tumor suppressor loss and targetable kinases. Acetylproteome analysis highlighted acetylation on key nuclear proteins involved in the DNA damage response and revealed cross-talk between cytoplasmic and mitochondrial acetylation and metabolism. Our results underscore the potential of proteogenomics for clinical investigation of breast cancer through more accurate annotation of targetable pathways and biological features of this remarkably heterogeneous malignancy.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes.

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

Single-cell RNA sequencing reveals compromised immune microenvironment in precursor stages of multiple myeloma.

Precursor states of Multiple Myeloma (MM) and its native tumor microenvironment need in-depth molecular characterization to better stratify and treat patients at risk. Using single-cell RNA sequencing of bone marrow cells from precursor stages, MGUS and smoldering myeloma (SMM), to full-blown MM alongside healthy donors, we demonstrate early immune changes during patient progression. We find NK cell abundance is frequently increased in early stages, and associated with altered chemokine receptor expression. As early as SMM, we show loss of GrK+ memory cytotoxic T-cells, and show their critical role in MM immunosurveillance in mouse models. Finally, we report MHC class II dysregulation in CD14+ monocytes, which results in T cell suppression in vitro. These results provide a comprehensive map of immune changes at play over the evolution of pre-malignant MM, which will help develop strategies for immune-based patient stratification.

Patient-Specific Tumor Growth Trajectories Determine Persistent and Resistant Cancer Cell Populations during Treatment with Targeted Therapies.

The importance of preexisting versus acquired drug resistance in patients with cancer treated with small-molecule tyrosine kinase inhibitors (TKI) remains controversial. The goal of this study is to provide a general estimate of the size and dynamics of a preexisting, drug-resistant tumor cell population versus a slow-growing persister population that is the precursor of acquired TKI resistance. We describe a general model of resistance development, including persister evolution and preexisting resistance, solely based on the macroscopic trajectory of tumor burden during treatment. We applied the model to 20 tumor volume trajectories of EGFR-mutant lung cancer patients treated with the TKI erlotinib. Under the assumption of only preexisting resistant cells or only persister evolution, it is not possible to explain the observed tumor trajectories with realistic parameter values. Assuming only persister evolution would require very high mutation induction rates, while only preexisting resistance would lead to very large preexisting populations of resistant cells at the initiation of treatment. However, combining preexisting resistance with persister populations can explain the observed tumor volume trajectories and yields an estimated preexisting resistant fraction varying from 10-4 to 10-1 at the time of treatment initiation for this study cohort. Our results also demonstrate that the growth rate of the resistant population is highly correlated to the time to tumor progression. These estimates of the size of the resistant and persistent tumor cell population during TKI treatment can inform combination treatment strategies such as multi-agent schedules or a combination of targeted agents and radiotherapy. SIGNIFICANCE: These findings quantify pre-existing resistance and persister cell populations, which are essential for the integration of targeted agents into the management of locally advanced disease and the timing of radiotherapy in metastatic patients.

WRN helicase is a synthetic lethal target in microsatellite unstable cancers.

Synthetic lethality-an interaction between two genetic events through which the co-occurrence of these two genetic events leads to cell death, but each event alone does not-can be exploited for cancer therapeutics1. DNA repair processes represent attractive synthetic lethal targets, because many cancers exhibit an impairment of a DNA repair pathway, which can lead to dependence on specific repair proteins2. The success of poly(ADP-ribose) polymerase 1 (PARP-1) inhibitors in cancers with deficiencies in homologous recombination highlights the potential of this approach3. Hypothesizing that other DNA repair defects would give rise to synthetic lethal relationships, we queried dependencies in cancers with microsatellite instability (MSI), which results from deficient DNA mismatch repair. Here we analysed data from large-scale silencing screens using CRISPR-Cas9-mediated knockout and RNA interference, and found that the RecQ DNA helicase WRN was selectively essential in MSI models in vitro and in vivo, yet dispensable in models of cancers that are microsatellite stable. Depletion of WRN induced double-stranded DNA breaks and promoted apoptosis and cell cycle arrest selectively in MSI models. MSI cancer models required the helicase activity of WRN, but not its exonuclease activity. These findings show that WRN is a synthetic lethal vulnerability and promising drug target for MSI cancers.

Quantification of somatic mutation flow across individual cell division events by lineage sequencing.

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon (POLE) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations.

Genetic and transcriptional evolution alters cancer cell line drug response.

Human cancer cell lines are the workhorse of cancer research. Although cell lines are known to evolve in culture, the extent of the resultant genetic and transcriptional heterogeneity and its functional consequences remain understudied. Here we use genomic analyses of 106 human cell lines grown in two laboratories to show extensive clonal diversity. Further comprehensive genomic characterization of 27 strains of the common breast cancer cell line MCF7 uncovered rapid genetic diversification. Similar results were obtained with multiple strains of 13 additional cell lines. Notably, genetic changes were associated with differential activation of gene expression programs and marked differences in cell morphology and proliferation. Barcoding experiments showed that cell line evolution occurs as a result of positive clonal selection that is highly sensitive to culture conditions. Analyses of single-cell-derived clones demonstrated that continuous instability quickly translates into heterogeneity of the cell line. When the 27 MCF7 strains were tested against 321 anti-cancer compounds, we uncovered considerably different drug responses: at least 75% of compounds that strongly inhibited some strains were completely inactive in others. This study documents the extent, origins and consequences of genetic variation within cell lines, and provides a framework for researchers to measure such variation in efforts to support maximally reproducible cancer research.

Analysis of somatic microsatellite indels identifies driver events in human tumors.

Microsatellites (MSs) are tracts of variable-length repeats of short DNA motifs that exhibit high rates of mutation in the form of insertions or deletions (indels) of the repeated motif. Despite their prevalence, the contribution of somatic MS indels to cancer has been largely unexplored, owing to difficulties in detecting them in short-read sequencing data. Here we present two tools: MSMuTect, for accurate detection of somatic MS indels, and MSMutSig, for identification of genes containing MS indels at a higher frequency than expected by chance. Applying MSMuTect to whole-exome data from 6,747 human tumors representing 20 tumor types, we identified >1,000 previously undescribed MS indels in cancer genes. Additionally, we demonstrate that the number and pattern of MS indels can accurately distinguish microsatellite-stable tumors from tumors with microsatellite instability, thus potentially improving classification of clinically relevant subgroups. Finally, we identified seven MS indel driver hotspots: four in known cancer genes (ACVR2A, RNF43, JAK1, and MSH3) and three in genes not previously implicated as cancer drivers (ESRP1, PRDM2, and DOCK3).