Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analysed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined four 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′UTR deletion as a mechanism of FGFR2 upregulation. Integration of non-coding 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.

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.

Many traditional pharmacopeias include and related plants, which contain nephrotoxins and mutagens in the form of aristolochic acids and similar compounds (collectively, AA). AA is implicated in multiple cancer types, sometimes with very high mutational burdens, especially in upper tract urothelial cancers (UTUCs). AA-associated kidney failure and UTUCs are prevalent in Taiwan, but AA's role in hepatocellular carcinomas (HCCs) there remains unexplored. Therefore, we sequenced the whole exomes of 98 HCCs from two hospitals in Taiwan and found that 78% showed the distinctive mutational signature of AA exposure, accounting for most of the nonsilent mutations in known cancer driver genes. We then searched for the AA signature in 1400 HCCs from diverse geographic regions. Consistent with exposure through known herbal medicines, 47% of Chinese HCCs showed the signature, albeit with lower mutation loads than in Taiwan. In addition, 29% of HCCs from Southeast Asia showed the signature. The AA signature was also detected in 13 and 2.7% of HCCs from Korea and Japan as well as in 4.8 and 1.7% of HCCs from North America and Europe, respectively, excluding one U.S. hospital where 22% of 87 "Asian" HCCs had the signature. Thus, AA exposure is geographically widespread. Asia, especially Taiwan, appears to be much more extensively affected, which is consistent with other evidence of patterns of AA exposure. We propose that additional measures aimed at primary prevention through avoidance of AA exposure and investigation of possible approaches to secondary prevention are warranted.

Mature T-cell lymphomas, including peripheral T-cell lymphoma (PTCL) and extranodal NK/T-cell lymphoma (NKTL), represent a heterogeneous group of non-Hodgkin lymphomas with dismal outcomes and limited treatment options. To determine the extent of involvement of the JAK/STAT pathway in this malignancy, we performed targeted capture sequencing of 188 genes in this pathway in 171 PTCL and NKTL cases. A total of 272 nonsynonymous somatic mutations in 101 genes were identified in 73% of the samples, including 258 single-nucleotide variants and 14 insertions or deletions. Recurrent mutations were most frequently located in and (15%), followed by and (6%) and (4%). A high prevalence of mutation (21%) was observed specifically in NKTL. Novel mutations (p.D427H, E616G, p.E616K, and p.E696K) were shown to increase STAT3 phosphorylation and transcriptional activity of in the absence of cytokine, in which p.E616K induced programmed cell death-ligand 1 (PD-L1) expression by robust binding of activated STAT3 to the PD-L1 gene promoter. Consistent with these findings, PD-L1 was overexpressed in NKTL cell lines harboring hotspot mutations, and similar findings were observed by the overexpression of p.E616K and p.E616G in the wild-type NKTL cell line. Conversely, STAT3 silencing and inhibition decreased PD-L1 expression in mutant NKTL cell lines. In NKTL tumors, STAT3 activation correlated significantly with PD-L1 expression. We demonstrated that STAT3 activation confers high PD-L1 expression, which may promote tumor immune evasion. The combination of PD-1/PD-L1 antibodies and STAT3 inhibitors might be a promising therapeutic approach for NKTL, and possibly PTCL.

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.