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.

Functional genomics identifies specific vulnerabilities in PTEN-deficient breast cancer.

BACKGROUND: Phosphatase and tensin homolog (PTEN) is one of the most frequently inactivated tumor suppressors in breast cancer. While PTEN itself is not considered a druggable target, PTEN synthetic-sick or synthetic-lethal (PTEN-SSL) genes are potential drug targets in PTEN-deficient breast cancers. Therefore, with the aim of identifying potential targets for precision breast cancer therapy, we sought to discover PTEN-SSL genes present in a broad spectrum of breast cancers. METHODS: To discover broad-spectrum PTEN-SSL genes in breast cancer, we used a multi-step approach that started with (1) a genome-wide short interfering RNA (siRNA) screen of ~ 21,000 genes in a pair of isogenic human mammary epithelial cell lines, followed by (2) a short hairpin RNA (shRNA) screen of ~ 1200 genes focused on hits from the first screen in a panel of 11 breast cancer cell lines; we then determined reproducibility of hits by (3) identification of overlaps between our results and reanalyzed data from 3 independent gene-essentiality screens, and finally, for selected candidate PTEN-SSL genes we (4) confirmed PTEN-SSL activity using either drug sensitivity experiments in a panel of 19 cell lines or mutual exclusivity analysis of publicly available pan-cancer somatic mutation data. RESULTS: The screens (steps 1 and 2) and the reproducibility analysis (step 3) identified six candidate broad-spectrum PTEN-SSL genes (PIK3CB, ADAMTS20, AP1M2, HMMR, STK11, and NUAK1). PIK3CB was previously identified as PTEN-SSL, while the other five genes represent novel PTEN-SSL candidates. Confirmation studies (step 4) provided additional evidence that NUAK1 and STK11 have PTEN-SSL patterns of activity. Consistent with PTEN-SSL status, inhibition of the NUAK1 protein kinase by the small molecule drug HTH-01-015 selectively impaired viability in multiple PTEN-deficient breast cancer cell lines, while mutations affecting STK11 and PTEN were largely mutually exclusive across large pan-cancer data sets. CONCLUSIONS: Six genes showed PTEN-SSL patterns of activity in a large proportion of PTEN-deficient breast cancer cell lines and are potential specific vulnerabilities in PTEN-deficient breast cancer. Furthermore, the NUAK1 PTEN-SSL vulnerability identified by RNA interference techniques can be recapitulated and exploited using the small molecule kinase inhibitor HTH-01-015. Thus, NUAK1 inhibition may be an effective strategy for precision treatment of PTEN-deficient breast tumors.

Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma.

Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analyzed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined 4 CCA clusters-fluke-positive CCAs (clusters 1/2) are enriched in ERBB2 amplifications and TP53 mutations; conversely, fluke-negative CCAs (clusters 3/4) exhibit high copy-number alterations and PD-1/PD-L2 expression, or epigenetic mutations (IDH1/2, BAP1) and FGFR/PRKA-related gene rearrangements. Whole-genome analysis highlighted FGFR2 3' untranslated region deletion as a mechanism of FGFR2 upregulation. Integration of noncoding promoter mutations with protein-DNA binding profiles demonstrates pervasive modulation of H3K27me3-associated sites in CCA. Clusters 1 and 4 exhibit distinct DNA hypermethylation patterns targeting either CpG islands or shores-mutation signature and subclonality analysis suggests that these reflect different mutational pathways. Our results exemplify how genetics, epigenetics, and environmental carcinogens can interplay across different geographies to generate distinct molecular subtypes of cancer.Significance: Integrated whole-genome and epigenomic analysis of CCA on an international scale identifies new CCA driver genes, noncoding promoter mutations, and structural variants. CCA molecular landscapes differ radically by etiology, underscoring how distinct cancer subtypes in the same organ may arise through different extrinsic and intrinsic carcinogenic processes. Cancer Discov; 7(10); 1116-35. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1047.

A formalin-fixed paraffin-embedded (FFPE)-based prognostic signature to predict metastasis in clinically low risk stage I/II microsatellite stable colorectal cancer.

Approximately 20% early-stage (I/II) colorectal cancer (CRC) patients develop metastases despite curative surgery. We aim to develop a formalin-fixed and paraffin-embedded (FFPE)-based predictor of metastases in early-stage, clinically-defined low risk, microsatellite-stable (MSS) CRC patients. We considered genome-wide mRNA and miRNA expression and mutation status of 20 genes assayed in 150 fresh-frozen tumours with known metastasis status. We selected 193 genes for further analysis using NanoString nCounter arrays on corresponding FFPE tumours. Neither mutation status nor miRNA expression improved the estimated prediction. The final predictor, ColoMet19, based on the top 19 genes' mRNA levels trained by Random Forest machine-learning strategy, had an estimated positive-predictive-value (PPV) of 0.66. We tested ColoMet19 on an independent test-set of 131 tumours and obtained a population-adjusted PPV of 0.67 indicating that early-stage CRC patients who tested positive have a 67% risk of developing metastases, substantially higher than the metastasis risk of 40% for node-positive (Stage III) patients who are generally treated with chemotherapy. Predicted-positive patients also had poorer metastasis-free survival (hazard ratios [HR] = 1.92, design-set; HR = 2.05, test-set). Thus, early-stage CRC patients who test positive may be considered for adjuvant therapy after surgery.

Transcription-associated mutation of lasR in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunistic pathogen which infects cystic fibrosis and cancer patients with compromised immune systems. LasR is a master regulator which controls the virulence of P. aeruginosa in response to bacterial cell-density and host signals. During infection, lasR is frequently mutated, conferring P. aeruginosa a growth advantage in hosts and enhances resistance to widely used antibiotics. However, the mechanistic basis of lasR mutation is not well understood. We have tested here the hypothesis that transcription strength is a contributory determinant of lasR mutagenesis. P. aeruginosa strains with different lasR transcription strengths were therefore engineered and the lasR mutations were monitored unbiasedly using next-generation sequencing technology. Our results suggest that the strength of transcription could be one of the deterministic factors that drive the mutagenesis of lasR in P. aeruginosa, shedding new insights into bacterial infection and antibiotic resistance.

Exome sequencing reveals recurrent REV3L mutations in cisplatin-resistant squamous cell carcinoma of head and neck.

Dacomitinib, an irreversible pan-HER inhibitor, had shown modest clinical activity in squamous cell carcinoma of head and neck (SCCHN) patients. Therefore, validated predictive biomarkers are required to identify patients most likely to benefit from this therapeutic option. To characterize the genetic landscape of cisplatin-treated SCCHN genomes and identify potential predictive biomarkers for dacomitinib sensitivity, we performed whole exome sequencing on 18 cisplatin-resistant metastatic SCCHN tumors and their matched germline DNA. Platinum-based chemotherapy elevated the mutation rates of SCCHN compared to chemotherapy-naïve SCCHNs. Cisplatin-treated SCCHN genomes uniquely exhibited a novel mutational signature characterized by C:G to A:T transversions at CCR sequence contexts that may have arisen due to error-prone translesional synthesis. Somatic mutations in REV3L, the gene encoding the catalytic subunit of DNA polymerase ζ involved in translesional synthesis, are significantly enriched in a subset of patients who derived extended clinical benefit to dacomitinib (P = 0.04). Functional assays showed that loss-of-function of REV3L dramatically enhanced the sensitivity of SCCHN cells to dacomitinib by the loss of both translesion synthesis and homologous recombination pathways. Our data suggest that the 'platinum' mutational signature and inactivation of REV3L may inform treatment options in patients of recurrent SCCHN.

SETD2 histone modifier loss in aggressive GI stromal tumours.

BACKGROUND: GI stromal tumours (GISTs) are clinically heterogenous exhibiting varying degrees of disease aggressiveness in individual patients. OBJECTIVES: We sought to identify genetic alterations associated with high-risk GIST, explore their molecular consequences, and test their utility as prognostic markers. DESIGNS: Exome sequencing of 18 GISTs was performed (9 patients with high-risk/metastatic and 5 patients with low/intermediate-risk), corresponding to 11 primary and 7 metastatic tumours. Candidate alterations were validated by prevalence screening in an independent patient cohort (n=120). Functional consequences of SETD2 mutations were investigated in primary tissues and cell lines. Transcriptomic profiles for 8 GISTs (4 SETD2 mutated, 4 SETD2 wild type) and DNA methylation profiles for 22 GISTs (10 SETD2 mutated, 12 SETD2 wild type) were analysed. Statistical associations between molecular, clinicopathological factors, and relapse-free survival were determined. RESULTS: High-risk GISTs harboured increased numbers of somatic mutations compared with low-risk GISTs (25.2 mutations/high-risk cases vs 6.8 mutations/low-risk cases; two sample t test p=3.1×10-5). Somatic alterations in the SETD2 histone modifier gene occurred in 3 out of 9 high-risk/metastatic cases but no low/intermediate-risk cases. Prevalence screening identified additional SETD2 mutations in 7 out of 80 high-risk/metastatic cases but no low/intermediate-risk cases (n=29). Combined, the frequency of SETD2 mutations was 11.2% (10/89) and 0% (0/34) in high-risk and low-risk GISTs respectively. SETD2 mutant GISTs exhibited decreased H3K36me3 expression while SETD2 silencing promoted DNA damage in GIST-T1 cells. In gastric GISTs, SETD2 mutations were associated with overexpression of HOXC cluster genes and a DNA methylation signature of hypomethylated heterochromatin. Gastric GISTs with SETD2 mutations, or GISTs with hypomethylated heterochromatin, showed significantly shorter relapse-free survival on univariate analysis (log rank p=4.1×10-5). CONCLUSIONS: Our data suggest that SETD2 is a novel GIST tumour suppressor gene associated with disease progression. Assessing SETD2 genetic status and SETD2-associated epigenomic phenotypes may guide risk stratification and provide insights into mechanisms of GIST clinical aggressiveness.

Genomic landscapes of breast fibroepithelial tumors.

Breast fibroepithelial tumors comprise a heterogeneous spectrum of pathological entities, from benign fibroadenomas to malignant phyllodes tumors. Although MED12 mutations have been frequently found in fibroadenomas and phyllodes tumors, the landscapes of genetic alterations across the fibroepithelial tumor spectrum remain unclear. Here, by performing exome sequencing of 22 phyllodes tumors followed by targeted sequencing of 100 breast fibroepithelial tumors, we observed three distinct somatic mutation patterns. First, we frequently observed MED12 and RARA mutations in both fibroadenomas and phyllodes tumors, emphasizing the importance of these mutations in fibroepithelial tumorigenesis. Second, phyllodes tumors exhibited mutations in FLNA, SETD2 and KMT2D, suggesting a role in driving phyllodes tumor development. Third, borderline and malignant phyllodes tumors harbored additional mutations in cancer-associated genes. RARA mutations exhibited clustering in the portion of the gene encoding the ligand-binding domain, functionally suppressed RARA-mediated transcriptional activation and enhanced RARA interactions with transcriptional co-repressors. This study provides insights into the molecular pathogenesis of breast fibroepithelial tumors, with potential clinical implications.

Mutational landscapes of tongue carcinoma reveal recurrent mutations in genes of therapeutic and prognostic relevance.

BACKGROUND: Carcinoma of the oral tongue (OTSCC) is the most common malignancy of the oral cavity, characterized by frequent recurrence and poor survival. The last three decades has witnessed a change in the OTSCC epidemiological profile, with increasing incidence in younger patients, females and never-smokers. Here, we sought to characterize the OTSCC genomic landscape and to determine factors that may delineate the genetic basis of this disease, inform prognosis and identify targets for therapeutic intervention. METHODS: Seventy-eight cases were subjected to whole-exome (n = 18) and targeted deep sequencing (n = 60). RESULTS: While the most common mutation was in TP53, the OTSCC genetic landscape differed from previously described cohorts of patients with head and neck tumors: OTSCCs demonstrated frequent mutations in DST and RNF213, while alterations in CDKN2A and NOTCH1 were significantly less frequent. Despite a lack of previously reported NOTCH1 mutations, integrated analysis showed enrichments of alterations affecting Notch signaling in OTSCC. Importantly, these Notch pathway alterations were prognostic on multivariate analyses. A high proportion of OTSCCs also presented with alterations in drug targetable and chromatin remodeling genes. Patients harboring mutations in actionable pathways were more likely to succumb from recurrent disease compared with those who did not, suggesting that the former should be considered for treatment with targeted compounds in future trials. CONCLUSIONS: Our study defines the Asian OTSCC mutational landscape, highlighting the key role of Notch signaling in oral tongue tumorigenesis. We also observed somatic mutations in multiple therapeutically relevant genes, which may represent candidate drug targets in this highly lethal tumor type.

MSIseq: Software for Assessing Microsatellite Instability from Catalogs of Somatic Mutations.

Microsatellite instability (MSI) is a form of hypermutation that occurs in some tumors due to defects in cellular DNA mismatch repair. MSI is characterized by frequent somatic mutations (i.e., cancer-specific mutations) that change the length of simple repeats (e.g., AAAAA…., GATAGATAGATA...). Clinical MSI tests evaluate the lengths of a handful of simple repeat sites, while next-generation sequencing can assay many more sites and offers a much more complete view of their somatic mutation frequencies. Using somatic mutation data from the exomes of a 361-tumor training set, we developed classifiers to determine MSI status based on four machine-learning frameworks. All frameworks had high accuracy, and after choosing one we determined that it had >98% concordance with clinical tests in a separate 163-tumor test set. Furthermore, this classifier retained high concordance even when classifying tumors based on subsets of whole-exome data. We have released a CRAN R package, MSIseq, based on this classifier. MSIseq is faster and simpler to use than software that requires large files of aligned sequenced reads. MSIseq will be useful for genomic studies in which clinical MSI test results are unavailable and for detecting possible misclassifications by clinical tests.

Mutation signatures implicate aristolochic acid in bladder cancer development.

BACKGROUND: Aristolochic acid (AA) is a natural compound found in many plants of the Aristolochia genus, and these plants are widely used in traditional medicines for numerous conditions and for weight loss. Previous work has connected AA-mutagenesis to upper-tract urothelial cell carcinomas and hepatocellular carcinomas. We hypothesize that AA may also contribute to bladder cancer. METHODS: Here, we investigated the involvement of AA-mutagenesis in bladder cancer by sequencing bladder tumor genomes from two patients with known exposure to AA. After detecting strong mutational signatures of AA exposure in these tumors, we exome-sequenced and analyzed an additional 11 bladder tumors and analyzed publicly available somatic mutation data from a further 336 bladder tumors. RESULTS: The somatic mutations in the bladder tumors from the two patients with known AA exposure showed overwhelming AA signatures. We also detected evidence of AA exposure in 1 out of 11 bladder tumors from Singapore and in 3 out of 99 bladder tumors from China. In addition, 1 out of 194 bladder tumors from North America showed a pattern of mutations that might have resulted from exposure to an unknown mutagen with a heretofore undescribed pattern of A > T mutations. Besides the signature of AA exposure, the bladder tumors also showed the CpG > TpG and activated-APOBEC signatures, which have been previously reported in bladder cancer. CONCLUSIONS: This study demonstrates the utility of inferring mutagenic exposures from somatic mutation spectra. Moreover, AA exposure in bladder cancer appears to be more pervasive in the East, where traditional herbal medicine is more widely used. More broadly, our results suggest that AA exposure is more extensive than previously thought both in terms of populations at risk and in terms of types of cancers involved. This appears to be an important public health issue that should be addressed by further investigation and by primary prevention through regulation and education. In addition to opportunities for primary prevention, knowledge of AA exposure would provide opportunities for secondary prevention in the form of intensified screening of patients with known or suspected AA exposure.

High-depth sequencing of over 750 genes supports linear progression of primary tumors and metastases in most patients with liver-limited metastatic colorectal cancer.

BACKGROUND: Colorectal cancer with metastases limited to the liver (liver-limited mCRC) is a distinct clinical subset characterized by possible cure with surgery. We performed high-depth sequencing of over 750 cancer-associated genes and copy number profiling in matched primary, metastasis and normal tissues to characterize genomic progression in 18 patients with liver-limited mCRC. RESULTS: High depth Illumina sequencing and use of three different variant callers enable comprehensive and accurate identification of somatic variants down to 2.5% variant allele frequency. We identify a median of 11 somatic single nucleotide variants (SNVs) per tumor. Across patients, a median of 79.3% of somatic SNVs present in the primary are present in the metastasis and 81.7% of all alterations present in the metastasis are present in the primary. Private alterations are found at lower allele frequencies; a different mutational signature characterized shared and private variants, suggesting distinct mutational processes. Using B-allele frequencies of heterozygous germline SNPs and copy number profiling, we find that broad regions of allelic imbalance and focal copy number changes, respectively, are generally shared between the primary tumor and metastasis. CONCLUSIONS: Our analyses point to high genomic concordance of primary tumor and metastasis, with a thick common trunk and smaller genomic branches in general support of the linear progression model in most patients with liver-limited mCRC. More extensive studies are warranted to further characterize genomic progression in this important clinical population.

Whole-exome sequencing studies of parathyroid carcinomas reveal novel PRUNE2 mutations, distinctive mutational spectra related to APOBEC-catalyzed DNA mutagenesis and mutational enrichment in kinases associated with cell migration and invasion.

CONTEXT: Cell division cycle 73 (CDC73), encoding the protein parafibromin, is the most prevalent mutated gene in familial and sporadic parathyroid carcinoma (PC). OBJECTIVE: To identify additional genetic abnormalities in PCs. DESIGN: Whole-exome sequencing was performed using DNA from seven pairs of matched PCs and one triplet containing double primary tumor and normal leukocyte. Somatic variants were confirmed using Sanger sequencing and recurrently mutated genes were assessed in 13 additional PCs as well as 40 parathyroid adenomas (PA). RESULTS: PC had an average of 51 somatic variants/tumor (range 3-176) with approximately 58% of variants occurring as nonsynonymous single nucleotide variants. The importance of CDC73 in PC is reinforced with a remarkable preferential amplification of the mutant CDC73 allele. Furthermore, recurrent germ line and somatic mutations in prune homolog 2 [Drosophila] (PRUNE2) were found in PC and computationally predicted to be deleterious; in addition, recurrent mutations in kinase genes related to cell migration and invasion were found. PRUNE2 showed recurrent mutations in 18% (4/22) of PCs with additional screening in 40 PAs revealing only one rare missense polymorphism (Asp1677Asn). For the first time, the mutational signature associated with apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC)-catalyzed cytosine-to-uracil deamination is found in a subset of PC. CONCLUSION: This study outlines the genetic landscape of PC and attempts to characterize the mutational processes shaping the PC genome.

Mutation signatures of carcinogen exposure: genome-wide detection and new opportunities for cancer prevention.

Exposure to environmental mutagens is an important cause of human cancer, and measures to reduce mutagenic and carcinogenic exposures have been highly successful at controlling cancer. Until recently, it has been possible to connect the chemical characteristics of mutagens to actual mutations observed in human tumors only indirectly. Now, next-generation sequencing technology enables us to observe in detail the DNA-sequence-level effects of well-known mutagens, such as ultraviolet radiation and tobacco smoke, as well as endogenous mutagenic processes, such as those involving activated DNA cytidine deaminases (APOBECs). We can also observe the effects of less well-known but potent mutagens, including those recently found to be present in some herbal remedies. Crucially, we can now tease apart the superimposed effects of several mutational exposures and processes and determine which ones occurred during the development of individual tumors. Here, we review advances in detecting these mutation signatures and discuss the implications for surveillance and prevention of cancer. The number of sequenced tumors from diverse cancer types and multiple geographic regions is growing explosively, and the genomes of these tumors will bear the signatures of even more diverse mutagenic exposures. Thus, we envision development of wide-ranging compendia of mutation signatures from tumors and a concerted effort to experimentally elucidate the signatures of a large number of mutagens. This information will be used to link signatures observed in tumors to the exposures responsible for them, which will offer unprecedented opportunities for prevention.

The NF1 gene revisited - from bench to bedside.

Neurofibromatosis type 1 (NF1) is a relatively common tumour predisposition syndrome related to germline aberrations of NF1, a tumour suppressor gene. The gene product neurofibromin is a negative regulator of the Ras cellular proliferation pathway, and also exerts tumour suppression via other mechanisms. Recent next-generation sequencing projects have revealed somatic NF1 aberrations in various sporadic tumours. NF1 plays a critical role in a wide range of tumours. NF1 alterations appear to be associated with resistance to therapy and adverse outcomes in several tumour types. Identification of a patient's germline or somatic NF1 aberrations can be challenging, as NF1 is one of the largest human genes, with a myriad of possible mutations. Epigenetic factors may also also contribute to inadequate levels of neurofibromin in cancer cells. Clinical trials of NF1-based therapeutic approaches are currently limited. Preclinical studies on neurofibromin-deficient malignancies have mainly been on malignant peripheral nerve sheath tumour cell lines or xenografts derived from NF1 patients. However, the emerging recognition of the role of NF1 in sporadic cancers may lead to the development of NF1-based treatments for other tumour types. Improved understanding of the implications of NF1 aberrations is critical for the development of novel therapeutic strategies.

Exome sequencing identifies highly recurrent MED12 somatic mutations in breast fibroadenoma.

Fibroadenomas are the most common breast tumors in women under 30 (refs. 1,2). Exome sequencing of eight fibroadenomas with matching whole-blood samples revealed recurrent somatic mutations solely in MED12, which encodes a Mediator complex subunit. Targeted sequencing of an additional 90 fibroadenomas confirmed highly frequent MED12 exon 2 mutations (58/98, 59%) that are probably somatic, with 71% of mutations occurring in codon 44. Using laser capture microdissection, we show that MED12 fibroadenoma mutations are present in stromal but not epithelial mammary cells. Expression profiling of MED12-mutated and wild-type fibroadenomas revealed that MED12 mutations are associated with dysregulated estrogen signaling and extracellular matrix organization. The fibroadenoma MED12 mutation spectrum is nearly identical to that of previously reported MED12 lesions in uterine leiomyoma but not those of other tumors. Benign tumors of the breast and uterus, both of which are key target tissues of estrogen, may thus share a common genetic basis underpinned by highly frequent and specific MED12 mutations.

Exome sequencing identifies distinct mutational patterns in liver fluke-related and non-infection-related bile duct cancers.

The impact of different carcinogenic exposures on the specific patterns of somatic mutation in human tumors remains unclear. To address this issue, we profiled 209 cholangiocarcinomas (CCAs) from Asia and Europe, including 108 cases caused by infection with the liver fluke Opisthorchis viverrini and 101 cases caused by non-O. viverrini-related etiologies. Whole-exome sequencing (n = 15) and prevalence screening (n = 194) identified recurrent somatic mutations in BAP1 and ARID1A, neither of which, to our knowledge, has previously been reported to be mutated in CCA. Comparisons between intrahepatic O. viverrini-related and non-O. viverrini-related CCAs demonstrated statistically significant differences in mutation patterns: BAP1, IDH1 and IDH2 were more frequently mutated in non-O. viverrini CCAs, whereas TP53 mutations showed the reciprocal pattern. Functional studies demonstrated tumor suppressive functions for BAP1 and ARID1A, establishing the role of chromatin modulators in CCA pathogenesis. These findings indicate that different causative etiologies may induce distinct somatic alterations, even within the same tumor type.

Genome-wide mutational signatures of aristolochic acid and its application as a screening tool.

Aristolochic acid (AA), a natural product of Aristolochia plants found in herbal remedies and health supplements, is a group 1 carcinogen that can cause nephrotoxicity and upper urinary tract urothelial cell carcinoma (UTUC). Whole-genome and exome analysis of nine AA-associated UTUCs revealed a strikingly high somatic mutation rate (150 mutations/Mb), exceeding smoking-associated lung cancer (8 mutations/Mb) and ultraviolet radiation-associated melanoma (111 mutations/Mb). The AA-UTUC mutational signature was characterized by A:T to T:A transversions at the sequence motif A[C|T]AGG, located primarily on nontranscribed strands. AA-induced mutations were also significantly enriched at splice sites, suggesting a role for splice-site mutations in UTUC pathogenesis. RNA sequencing of AA-UTUC confirmed a general up-regulation of nonsense-mediated decay machinery components and aberrant splicing events associated with splice-site mutations. We observed a high frequency of somatic mutations in chromatin modifiers, particularly KDM6A, in AA-UTUC, demonstrated the sufficiency of AA to induce renal dysplasia in mice, and reproduced the AA mutational signature in experimentally treated human renal tubular cells. Finally, exploring other malignancies that were not known to be associated with AA, we screened 93 hepatocellular carcinoma genomes/exomes and identified AA-like mutational signatures in 11. Our study highlights an unusual genome-wide AA mutational signature and the potential use of mutation signatures as "molecular fingerprints" for interrogating high-throughput cancer genome data to infer previous carcinogen exposures.

Exome sequencing of liver fluke-associated cholangiocarcinoma.

Opisthorchis viverrini-related cholangiocarcinoma (CCA), a fatal bile duct cancer, is a major public health concern in areas endemic for this parasite. We report here whole-exome sequencing of eight O. viverrini-related tumors and matched normal tissue. We identified and validated 206 somatic mutations in 187 genes using Sanger sequencing and selected 15 genes for mutation prevalence screening in an additional 46 individuals with CCA (cases). In addition to the known cancer-related genes TP53 (mutated in 44.4% of cases), KRAS (16.7%) and SMAD4 (16.7%), we identified somatic mutations in 10 newly implicated genes in 14.8-3.7% of cases. These included inactivating mutations in MLL3 (in 14.8% of cases), ROBO2 (9.3%), RNF43 (9.3%) and PEG3 (5.6%), and activating mutations in the GNAS oncogene (9.3%). These genes have functions that can be broadly grouped into three biological classes: (i) deactivation of histone modifiers, (ii) activation of G protein signaling and (iii) loss of genome stability. This study provides insight into the mutational landscape contributing to O. viverrini-related CCA.