The repertoire of mutational signatures in human cancer.

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Recurrent mutations in topoisomerase IIα cause a previously undescribed mutator phenotype in human cancers.

Topoisomerases nick and reseal DNA to relieve torsional stress associated with transcription and replication and to resolve structures such as knots and catenanes. Stabilization of the yeast Top2 cleavage intermediates is mutagenic in yeast, but whether this extends to higher eukaryotes is less clear. Chemotherapeutic topoisomerase poisons also elevate cleavage, resulting in mutagenesis. Here, we describe p.K743N mutations in human topoisomerase hTOP2α and link them to a previously undescribed mutator phenotype in cancer. Overexpression of the orthologous mutant protein in yeast generated a characteristic pattern of 2- to 4-base pair (bp) duplications resembling those in tumors with p.K743N. Using mutant strains and biochemical analysis, we determined the genetic requirements of this mutagenic process and showed that it results from trapping of the mutant yeast yTop2 cleavage complex. In addition to 2- to 4-bp duplications, hTOP2α p.K743N is also associated with deletions that are absent in yeast. We call the combined pattern of duplications and deletions ID_TOP2α. All seven tumors carrying the hTOP2α p.K743N mutation showed ID_TOP2α, while it was absent from all other tumors examined (n = 12,269). Each tumor with the ID_TOP2α signature had indels in several known cancer genes, which included frameshift mutations in tumor suppressors PTEN and TP53 and an activating insertion in BRAF. Sequence motifs found at ID_TOP2α mutations were present at 80% of indels in cancer-driver genes, suggesting that ID_TOP2α mutagenesis may contribute to tumorigenesis. The results reported here shed further light on the role of topoisomerase II in genome instability.

Comparison of the performance of two targeted metagenomic virus capture probe-based methods using reference control materials and clinical samples.

Viral enrichment by probe hybridization has been reported to significantly increase the sensitivity of viral metagenomics. This study compares the analytical performance of two targeted metagenomic virus capture probe-based methods: (i) SeqCap EZ HyperCap by Roche (ViroCap) and (ii) Twist Comprehensive Viral Research Panel workflow, for diagnostic use. Sensitivity, specificity, and limit of detection were analyzed using 25 synthetic viral sequences spiked in increasing proportions of human background DNA, eight clinical samples, and American Type Culture Collection (ATCC) Virome Virus Mix. Sensitivity and specificity were 95% and higher for both methods using the synthetic and reference controls as gold standard. Combining thresholds for viral sequence read counts and genome coverage [respectively 500 reads per million (RPM) and 10% coverage] resulted in optimal prediction of true positive results. Limits of detection were approximately 50-500 copies/mL for both methods as determined by ddPCR. Increasing proportions of spike-in cell-free human background sequences up to 99.999% (50 ng/mL) did not negatively affect viral detection, suggesting effective capture of viral sequences. These data show analytical performances in ranges applicable to clinical samples, for both probe hybridization metagenomic approaches. This study supports further steps toward more widespread use of viral metagenomics for pathogen detection, in clinical and surveillance settings using low biomass samples. IMPORTANCE: Viral metagenomics has been gradually applied for broad-spectrum pathogen detection of infectious diseases, surveillance of emerging diseases, and pathogen discovery. Viral enrichment by probe hybridization methods has been reported to significantly increase the sensitivity of viral metagenomics. During the past years, a specific hybridization panel distributed by Roche has been adopted in a broad range of different clinical and zoonotic settings. Recently, Twist Bioscience has released a new hybridization panel targeting human and animal viruses. This is the first report comparing the performance of viral metagenomic hybridization panels.

Enrichment of colibactin-associated mutational signatures in unexplained colorectal polyposis patients.

BACKGROUND: Colibactin, a genotoxin produced by polyketide synthase harboring (pks+) bacteria, induces double-strand breaks and chromosome aberrations. Consequently, enrichment of pks+Escherichia coli in colorectal cancer and polyposis suggests a possible carcinogenic effect in the large intestine. Additionally, specific colibactin-associated mutational signatures; SBS88 and ID18 in the Catalogue of Somatic Mutations in Cancer database, are detected in colorectal carcinomas. Previous research showed that a recurrent APC splice variant perfectly fits SBS88. METHODS: In this study, we explore the presence of colibactin-associated signatures and fecal pks in an unexplained polyposis cohort. Somatic targeted Next-Generation Sequencing (NGS) was performed for 379 patients. Additionally, for a subset of 29 patients, metagenomics was performed on feces and mutational signature analyses using Whole-Genome Sequencing (WGS) on Formalin-Fixed Paraffin Embedded (FFPE) colorectal tissue blocks. RESULTS: NGS showed somatic APC variants fitting SBS88 or ID18 in at least one colorectal adenoma or carcinoma in 29% of patients. Fecal metagenomic analyses revealed enriched presence of pks genes in patients with somatic variants fitting colibactin-associated signatures compared to patients without variants fitting colibactin-associated signatures. Also, mutational signature analyses showed enrichment of SBS88 and ID18 in patients with variants fitting these signatures in NGS compared to patients without. CONCLUSIONS: These findings further support colibactins ability to mutagenize colorectal mucosa and contribute to the development of colorectal adenomas and carcinomas explaining a relevant part of patients with unexplained polyposis.

Assessing pathogenicity of mismatch repair variants of uncertain significance by molecular tumor analysis.

Functional analyses are the main method to classify mismatch repair (MMR) gene variants of uncertain significance (VUSs). However, the pathogenicity remains unclear for many variants because of conflicting results between clinical, molecular, and functional data. In this study, we evaluated whether whole exome sequencing (WES) could add another layer of evidence to elucidate the pathogenicity of MMR variants with conflicting interpretations. WES was performed on formalin-fixed paraffin-embedded tumor tissue of eight patients with a constitutional MMR VUS (seven families), including eight colorectal and two endometrial carcinomas and one ovarian carcinoma. Cell-free CIMRA assays were performed to assign Odds of Pathogenicity to these VUSs. In four families, seven tumors showed MMR deficiency-associated mutational signatures, supporting the pathogenicity of the VUS. Moreover, somatic (second) MMR hits identified in the WES data were found to explain MMR staining patterns when the MMR staining was discordant with the reported germline MMR gene variant. In conclusion, WES did not significantly reclassify VUS in these cases but clarified some phenotypic aspects such as age of onset and explanations in case of discordant MMR stainings.

Integrative multiomics enhancer activity profiling identifies therapeutic vulnerabilities in cholangiocarcinoma of different etiologies.

OBJECTIVES: Cholangiocarcinoma (CCA) is a heterogeneous malignancy with high mortality and dismal prognosis, and an urgent clinical need for new therapies. Knowledge of the CCA epigenome is largely limited to aberrant DNA methylation. Dysregulation of enhancer activities has been identified to affect carcinogenesis and leveraged for new therapies but is uninvestigated in CCA. Our aim is to identify potential therapeutic targets in different subtypes of CCA through enhancer profiling. DESIGN: Integrative multiomics enhancer activity profiling of diverse CCA was performed. A panel of diverse CCA cell lines, patient-derived and cell line-derived xenografts were used to study identified enriched pathways and vulnerabilities. NanoString, multiplex immunohistochemistry staining and single-cell spatial transcriptomics were used to explore the immunogenicity of diverse CCA. RESULTS: We identified three distinct groups, associated with different etiologies and unique pathways. Drug inhibitors of identified pathways reduced tumour growth in in vitro and in vivo models. The first group (ESTRO), with mostly fluke-positive CCAs, displayed activation in estrogen signalling and were sensitive to MTOR inhibitors. Another group (OXPHO), with mostly BAP1 and IDH-mutant CCAs, displayed activated oxidative phosphorylation pathways, and were sensitive to oxidative phosphorylation inhibitors. Immune-related pathways were activated in the final group (IMMUN), made up of an immunogenic CCA subtype and CCA with aristolochic acid (AA) mutational signatures. Intratumour differences in AA mutation load were correlated to intratumour variation of different immune cell populations. CONCLUSION: Our study elucidates the mechanisms underlying enhancer dysregulation and deepens understanding of different tumourigenesis processes in distinct CCA subtypes, with potential significant therapeutics and clinical benefits.

Characterization of colibactin-associated mutational signature in an Asian oral squamous cell carcinoma and in other mucosal tumor types.

Mutational signatures can reveal the history of mutagenic processes that cells were exposed to before and during tumorigenesis. We expect that as-yet-undiscovered mutational processes will shed further light on mutagenesis leading to carcinogenesis. With this in mind, we analyzed the mutational spectra of 36 Asian oral squamous cell carcinomas. The mutational spectra of two samples from patients who presented with oral bacterial infections showed novel mutational signatures. One of these novel signatures, SBS_AnT, is characterized by a preponderance of thymine mutations, strong transcriptional strand bias, and enrichment for adenines in the 4 bp 5' of mutation sites. The mutational signature described in this manuscript was shown to be caused by colibactin, a bacterial mutagen produced by E. coli carrying the pks-island. Examination of publicly available sequencing data revealed SBS_AnT in 25 tumors from several mucosal tissue types, expanding the list of tissues in which this mutational signature is observed.

A tumor-associated splice-isoform of MAP2K7 drives dedifferentiation in MBNL1-low cancers via JNK activation.

Master splicing regulator MBNL1 shapes large transcriptomic changes that drive cellular differentiation during development. Here we demonstrate that MBNL1 is a suppressor of tumor dedifferentiation. We surveyed MBNL1 expression in matched tumor/normal pairs across The Cancer Genome Atlas and found that MBNL1 was down-regulated in several common cancers. Down-regulation of MBNL1 predicted poor overall survival in breast, lung, and stomach adenocarcinomas and increased relapse and distant metastasis in triple-negative breast cancer. Down-regulation of MBNL1 led to increased tumorigenic and stem/progenitor-like properties in vitro and in vivo. A discrete set of alternative splicing events (ASEs) are shared between MBNL1-low cancers and embryonic stem cells including a MAP2K7∆exon2 splice variant that leads to increased stem/progenitor-like properties via JNK activation. Accordingly, JNK inhibition is capable of reversing MAP2K7∆exon2-driven tumor dedifferentiation in MBNL1-low cancer cells. Our work elucidates an alternative-splicing mechanism that drives tumor dedifferentiation and identifies biomarkers that predict enhanced susceptibility to JNK inhibition.

The mutational landscape of early- and typical-onset oral tongue squamous cell carcinoma.

BACKGROUND: The incidence of oral tongue squamous cell carcinoma (OTSCC) is increasing among younger birth cohorts. The etiology of early-onset OTSCC (diagnosed before the age of 50 years) and cancer driver genes remain largely unknown. METHODS: The Sequencing Consortium of Oral Tongue Cancer was established through the pooling of somatic mutation data of oral tongue cancer specimens (n = 227 [107 early-onset cases]) from 7 studies and The Cancer Genome Atlas. Somatic mutations at microsatellite loci and Catalog of Somatic Mutations in Cancer mutation signatures were identified. Cancer driver genes were identified with the MutSigCV and WITER algorithms. Mutation comparisons between early- and typical-onset OTSCC were evaluated via linear regression with adjustments for patient-related factors. RESULTS: Two novel driver genes (ATXN1 and CDC42EP1) and 5 previously reported driver genes (TP53, CDKN2A, CASP8, NOTCH1, and FAT1) were identified. Six recurrent mutations were identified, with 4 occurring in TP53. Early-onset OTSCC had significantly fewer nonsilent mutations even after adjustments for tobacco use. No associations of microsatellite locus mutations and mutation signatures with the age of OTSCC onset were observed. CONCLUSIONS: This international, multicenter consortium is the largest study to characterize the somatic mutational landscape of OTSCC and the first to suggest differences by age of onset. This study validates multiple previously identified OTSCC driver genes and proposes 2 novel cancer driver genes. In analyses by age, early-onset OTSCC had a significantly smaller somatic mutational burden that was not explained by differences in tobacco use. LAY SUMMARY: This study identifies 7 specific areas in the human genetic code that could be responsible for promoting the development of tongue cancer. Tongue cancer in young patients (under the age of 50 years) has fewer overall changes to the genetic code in comparison with tongue cancer in older patients, but the authors do not think that this is due to differences in smoking rates between the 2 groups. The cause of increasing cases of tongue cancer in young patients remains unclear.

In-depth characterization of the cisplatin mutational signature in human cell lines and in esophageal and liver tumors.

Cisplatin reacts with DNA and thereby likely generates a characteristic pattern of somatic mutations, called a mutational signature. Despite widespread use of cisplatin in cancer treatment and its role in contributing to secondary malignancies, its mutational signature has not been delineated. We hypothesize that cisplatin's mutational signature can serve as a biomarker to identify cisplatin mutagenesis in suspected secondary malignancies. Knowledge of which tissues are at risk of developing cisplatin-induced secondary malignancies could lead to guidelines for noninvasive monitoring for secondary malignancies after cisplatin chemotherapy. We performed whole genome sequencing of 10 independent clones of cisplatin-exposed MCF-10A and HepG2 cells and delineated the patterns of single and dinucleotide mutations in terms of flanking sequence, transcription strand bias, and other characteristics. We used the mSigAct signature presence test and nonnegative matrix factorization to search for cisplatin mutagenesis in hepatocellular carcinomas and esophageal adenocarcinomas. All clones showed highly consistent patterns of single and dinucleotide substitutions. The proportion of dinucleotide substitutions was high: 8.1% of single nucleotide substitutions were part of dinucleotide substitutions, presumably due to cisplatin's propensity to form intra- and interstrand crosslinks between purine bases in DNA. We identified likely cisplatin exposure in nine hepatocellular carcinomas and three esophageal adenocarcinomas. All hepatocellular carcinomas for which clinical data were available and all esophageal cancers indeed had histories of cisplatin treatment. We experimentally delineated the single and dinucleotide mutational signature of cisplatin. This signature enabled us to detect previous cisplatin exposure in human hepatocellular carcinomas and esophageal adenocarcinomas with high confidence.

Genome-scale mutational signatures of aflatoxin in cells, mice, and human tumors.

Aflatoxin B1 (AFB1) is a mutagen and IARC (International Agency for Research on Cancer) Group 1 carcinogen that causes hepatocellular carcinoma (HCC). Here, we present the first whole-genome data on the mutational signatures of AFB1 exposure from a total of >40,000 mutations in four experimental systems: two different human cell lines, in liver tumors in wild-type mice, and in mice that carried a hepatitis B surface antigen transgene-this to model the multiplicative effects of aflatoxin exposure and hepatitis B in causing HCC. AFB1 mutational signatures from all four experimental systems were remarkably similar. We integrated the experimental mutational signatures with data from newly sequenced HCCs from Qidong County, China, a region of well-studied aflatoxin exposure. This indicated that COSMIC mutational signature 24, previously hypothesized to stem from aflatoxin exposure, indeed likely represents AFB1 exposure, possibly combined with other exposures. Among published somatic mutation data, we found evidence of AFB1 exposure in 0.7% of HCCs treated in North America, 1% of HCCs from Japan, but 16% of HCCs from Hong Kong. Thus, aflatoxin exposure apparently remains a substantial public health issue in some areas. This aspect of our study exemplifies the promise of future widespread resequencing of tumor genomes in providing new insights into the contribution of mutagenic exposures to cancer incidence.

HLA-G protein expression in colorectal cancer evaluated by immunohistochemistry and western blot analysis: Its expression characteristics remain enigmatic.

HLA-G protein expression could play a role in evasion of tumor immune surveillance. Accumulating evidence demonstrates that HLA-G is expressed in different types of malignancies, including colorectal cancer (CRC). The purpose of the current study was to further unravel whether HLA-G protein expression could play a role in immune evasion of CRC. Therefore, to firmly establish HLA-G protein expression, eight early passage human CRC cell lines and five human rectal cancer tissues were analyzed by western blot analysis. The results obtained by western blot analysis were compared with immunohistochemistry on tumor tissue sections of the same patient. Furthermore, multiple monoclonal antibodies (mAbs), 4H84, MEM-G/1 and 5A6G7, targeting HLA-G were used to unravel staining patterns. We showed that results obtained with immunohistochemistry did not correlate with protein expression detected by western blot analysis, using three different HLA-G targeting mAbs. Furthermore, with respect to the specificity of the mAbs employed, additional immune reactivity was detected using the mAbs MEM-G/1 and 5A6G7 in western blot analysis with K562 control cell lines overexpressing HLA-A2 or HLA-G, all tumor tissues and in two out of eight CRC cell lines. Based on the current study and our previously reported results, we conclude that claiming HLA-G plays a role in immune modulation of CRC seems premature, as results from anti-body based detection of HLA-G protein remain inconclusive. Until the time that detection of HLA-G is sensitive enough to detect all aspects of HLA-G expression in biological samples, rather than transfected cells or long time cultured cell lines, conclusions should be drawn with great care.

Evidence for genetic association between chromosome 1q loci and predisposition to colorectal neoplasia.

BACKGROUND: A substantial fraction of familial colorectal cancer (CRC) and polyposis heritability remains unexplained. This study aimed to identify predisposing loci in patients with these disorders. METHODS: Homozygosity mapping was performed using 222 563 SNPs in 302 index patients with various colorectal neoplasms and 3367 controls. Linkage analysis, exome and whole-genome sequencing were performed in a family affected by microsatellite stable CRCs. Candidate variants were genotyped in 10 554 cases and 21 480 controls. Gene expression was assessed at the mRNA and protein level. RESULTS: Homozygosity mapping revealed a disease-associated region at 1q32.3 which was part of the linkage region 1q32.2-42.2 identified in the CRC family. This includes a region previously associated with risk of CRC. Sequencing identified the p.Asp1432Glu variant in the MIA3 gene (known as TANGO1 or TANGO) and 472 additional rare, shared variants within the linkage region. In both cases and controls the population frequency was 0.02% for this MIA3 variant. The MIA3 mutant allele showed predominant mRNA expression in normal, cancer and precancerous tissues. Furthermore, immunohistochemistry revealed increased expression of MIA3 in adenomatous tissues. CONCLUSIONS: Taken together, our two independent strategies associate genetic variations in chromosome 1q loci and predisposition to familial CRC and polyps, which warrants further investigation.

Imprinted survival genes preclude loss of heterozygosity of chromosome 7 in cancer cells.

The genomes of a wide range of cancers, including colon, breast, and thyroid cancers, frequently show copy number gains of chromosome 7 and rarely show loss of heterozygosity. The molecular basis for this phenomenon is unknown. Strikingly, oncocytic follicular thyroid carcinomas can display an extreme genomic profile, with homozygosity of all chromosomes except for chromosome 7. The observation that homozygosity of chromosome 7 is never observed suggests that retention of heterozygosity is essential for cells. We hypothesized that cell survival genes are genetically imprinted on either of two copies of chromosome 7, which thwarts loss of heterozygosity at this chromosome in cancer cells. By employing a DNA methylation screen and gene expression analysis, we identified six imprinted genes that force retention of heterozygosity on chromosome 7. Subsequent knockdown of gene expression showed that CALCR, COPG2, GRB10, KLF14, MEST, and PEG10 were essential for cancer cell survival, resulting in reduced cell proliferation, G1 -phase arrest, and increased apoptosis. We propose that imprinted cell survival genes provide a genetic basis for retention of chromosome 7 heterozygosity in cancer cells. The monoallelically expressed cell survival genes identified in this study, and the cellular pathways that they are involved in, offer new therapeutic targets for the treatment of tumours showing retention of heterozygosity on chromosome 7. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Chemosensitivity of BRCA1-Mutated Ovarian Cancer Cells and Established Cytotoxic Agents.

OBJECTIVE: Serous adenocarcinomas that arise in patients with inherited mutations in the tumor suppressor genes BRCA1 and BRCA2 are initially well treatable with platinum/paclitaxel. For recurrent disease in patients with BRCA1 or BRCA2 mutations, olaparib treatment is available. To study additional therapeutic regimens, a better understanding of the cellular and molecular mechanisms of the tumors in in vitro models is important. METHODS/MATERIALS: From a high-grade serous ovarian tumor of a BRCA1 mutation carrier, we established 3 distinct cell line subclones, OVCA-TR3.1, -2, and -3. Immunohistochemical characterization, flow cytometric analyses, chemosensitivity, and somatic mutation profiling were performed. RESULTS: The cell lines expressed AE1/AE3, Pax8, WT-1, OC125, estrogen receptor (ER), and p53, comparable to the primary tumor. Synergism could be shown in the combination treatment eremophila-1-(10)-11(13)-dien-12,8ß-olide (EPD), with cisplatin, whereas combination with olaparib did not show synergism. Eremophila-1-(10)-11(13)-dien-12,8ß-olide, a sesquiterpene lactone, is a novel chemotherapeutic agent. The inherited BRCA1 c.2989_2990dupAA mutation was confirmed in the cell lines. Loss of heterozygosity of BRCA1 was detected in each cell line, as well as a homozygous TP53 c.722C>A mutation. Flow cytometry showed that all cell lines had a distinct DNA index. CONCLUSIONS: Three new isogenic ovarian cancer cell lines were developed from a patient with a germ line BRCA1 mutation. Chemosensitivity profiling of the cell lines showed high tolerance for olaparib. Treatment with EPD proved synergistic with cisplatin. The effects of EPD will be further investigated for future clinical efficacy.

Tumor LINE-1 Methylation Level in Association with Survival of Patients with Stage II Colon Cancer.

Genome-wide DNA hypomethylation is associated with a worse prognosis in early-stage colorectal cancer. To measure genome-wide DNA methylation levels, long interspersed nucleotide element (LINE-1) repeats are used as a surrogate marker. Cohort studies on the clinical impact of genome-wide DNA methylation level in patients with only early-stage colon cancer, are currently lacking. This study aimed to investigate the prognostic value of LINE-1 methylation in a stage II colon cancer cohort (n = 164). Manual needle microdissection of tumor areas was performed on formalin-fixed paraffin-embedded tumor tissue sections followed by DNA extraction. Bisulfite converted DNA was used to assess tumor LINE-1 methylation level by qPCR. Patients with LINE-1 hypomethylated tumors had a significantly worse overall survival compared to patients with a higher level of LINE-1 tumor DNA methylation (HR 1.68, 95% CI 1.03-2.75; p = 0.04). This effect was more prominent in patients aged over 65 years (HR 2.00, 95% CI 1.13-3.52; p = 0.02), although the test for age interaction was not significant. No significant effect on recurrence-free survival was observed. Based on these results, tumor LINE-1 hypomethylation is associated with a worse overall survival in stage II colon cancer. Whether the origin of this causation is cancer-specific or age-related can be debated.