Willie Yu

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.

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

Genome-wide mutational signatures of the group 1 carcinogen aristolochic acid are observed in urothelial cancers and liver cancers from Asia. Carcinogen AAlert Aristolochic acid (AA) is a natural compound derived from plants in the Aristolochia genus. For centuries, Aristolochia has been used throughout Asia to treat a variety of ailments as a component of traditional Chinese medicine. In recent years, however, a more sinister side of this herb has come to light when it was linked to kidney damage and cancers of the urinary tract. Now, two studies by Poon et al. and Hoang et al. present a “molecular signature” of AA-induced DNA damage, which helps to explain the mutagenic effects of AA and may also be useful as a way to detect unsuspected AA exposure as a cause of cancer. The molecular signature seen in AA-associated tumors is characterized by a predominance of A:T-to-T:A transversions, a relatively unusual type of mutation that is infrequently seen in other types of cancer, including those caused by other carcinogens. These mutations concentrate at splice sites, causing the inappropriate inclusion or exclusion of entire exons in the resulting mRNA. The overall mutation rate is another notable feature of AA-associated cancers because it is several times higher than the rate of mutations caused by other carcinogens such as tobacco and ultraviolet light. In both studies, the authors also used the molecular signature to discover that AA was a likely cause of tumors previously attributed to other carcinogens. In one case, a urinary tract cancer that had been attributed to smoking and, in the other case, a liver cancer previously attributed to a chronic hepatitis infection were both identified as having the telltale signature of AA mutagenesis. The identification of a specific molecular signature for AA has both clinical and public health implications. For individual patients, the molecular signature could help physicians identify which tumors were caused by AA. Although this information cannot yet be used to optimize the treatment of individual patients, those who are diagnosed with AA-associated cancers could be monitored more closely for the appearance of additional tumors. Meanwhile, a better understanding of the mutagenic effects of AA should also help to strengthen public health efforts to decrease exposure to this carcinogenic herb. 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.

Janus Kinase 3–Activating Mutations Identified in Natural Killer/T-cell Lymphoma

UNLABELLED The molecular pathogenesis of natural killer/T-cell lymphoma (NKTCL) is not well understood. We conducted whole-exome sequencing and identified Janus kinase 3 (JAK3) somatic-activating mutations (A572V and A573V) in 2 of 4 patients with NKTCLs. Further validation of the prevalence of JAK3 mutations was determined by Sanger sequencing and high-resolution melt (HRM) analysis in an additional 61 cases. In total, 23 of 65 (35.4%) cases harbored JAK3 mutations. Functional characterization of the JAK3 mutations support its involvement in cytokine-independent JAK/STAT constitutive activation leading to increased cell growth. Moreover, treatment of both JAK3-mutant and wild-type NKTCL cell lines with a novel pan-JAK inhibitor, CP-690550, resulted in dose-dependent reduction of phosphorylated STAT5, reduced cell viability, and increased apoptosis. Hence, targeting the deregulated JAK/STAT pathway could be a promising therapy for patients with NKTCLs. SIGNIFICANCE Gene mutations causing NKTCL have not been fully identified. Through exome sequencing, we identified activating mutations of JAK3 that may play a significant role in the pathogenesis of NKTCLs. Our findings have important implications for the management of patients with NKTCLs.

Genetic and Structural Variation in the Gastric Cancer Kinome Revealed through Targeted Deep Sequencing

Genetic alterations in kinases have been linked to multiple human pathologies. To explore the landscape of kinase genetic variation in gastric cancer (GC), we used targeted, paired-end deep sequencing to analyze 532 protein and phosphoinositide kinases in 14 GC cell lines. We identified 10,604 single-nucleotide variants (SNV) in kinase exons including greater than 300 novel nonsynonymous SNVs. Family-wise analysis of the nonsynonymous SNVs revealed a significant enrichment in mitogen-activated protein kinase (MAPK)-related genes (P < 0.01), suggesting a preferential involvement of this kinase family in GC. A potential antioncogenic role for MAP2K4, a gene exhibiting recurrent alterations in 2 lines, was functionally supported by siRNA knockdown and overexpression studies in wild-type and MAP2K4 variant lines. The deep sequencing data also revealed novel, large-scale structural rearrangement events involving kinases including gene fusions involving CDK12 and the ERBB2 receptor tyrosine kinase in MKN7 cells. Integrating SNVs and copy number alterations, we identified Hs746T as a cell line exhibiting both splice-site mutations and genomic amplification of MET, resulting in MET protein overexpression. When applied to primary GCs, we identified somatic mutations in 8 kinases, 4 of which were recurrently altered in both primary tumors and cell lines (MAP3K6, STK31, FER, and CDKL5). These results demonstrate that how targeted deep sequencing approaches can deliver unprecedented multilevel characterization of a medically and pharmacologically relevant gene family. The catalog of kinome genetic variants assembled here may broaden our knowledge on kinases and provide useful information on genetic alterations in GC.

Next-generation sequencing of translocation renal cell carcinoma reveals novel RNA splicing partners and frequent mutations of chromatin-remodeling genes.

Purpose: MITF/TFE translocation renal cell carcinoma (TRCC) is a rare subtype of kidney cancer. Its incidence and the genome-wide characterization of its genetic origin have not been fully elucidated. Experimental Design: We performed RNA and exome sequencing on an exploratory set of TRCC (n = 7), and validated our findings using The Cancer Genome Atlas (TCGA) clear-cell RCC (ccRCC) dataset (n = 460). Results: Using the TCGA dataset, we identified seven TRCC (1.5%) cases and determined their genomic profile. We discovered three novel partners of MITF/TFE (LUC7L3, KHSRP, and KHDRBS2) that are involved in RNA splicing. TRCC displayed a unique gene expression signature as compared with other RCC types, and showed activation of MITF, the transforming growth factor β1 and the PI3K complex targets. Genes differentially spliced between TRCC and other RCC types were enriched for MITF and ID2 targets. Exome sequencing of TRCC revealed a distinct mutational spectrum as compared with ccRCC, with frequent mutations in chromatin-remodeling genes (six of eight cases, three of which were from the TCGA). In two cases, we identified mutations in INO80D, an ATP-dependent chromatin-remodeling gene, previously shown to control the amplitude of the S phase. Knockdown of INO80D decreased cell proliferation in a novel cell line bearing LUC7L3–TFE3 translocation. Conclusions: This genome-wide study defines the incidence of TRCC within a ccRCC-directed project and expands the genomic spectrum of TRCC by identifying novel MITF/TFE partners involved in RNA splicing and frequent mutations in chromatin-remodeling genes. Clin Cancer Res; 20(15); 4129–40. ©2014 AACR.

Regulatory crosstalk between lineage-survival oncogenes KLF5, GATA4 and GATA6 cooperatively promotes gastric cancer development.

Objective Gastric cancer (GC) is a deadly malignancy for which new therapeutic strategies are needed. Three transcription factors, KLF5, GATA4 and GATA6, have been previously reported to exhibit genomic amplification in GC. We sought to validate these findings, investigate how these factors function to promote GC, and identify potential treatment strategies for GCs harbouring these amplifications. Design KLF5, GATA4 and GATA6 copy number and gene expression was examined in multiple GC cohorts. Chromatin immunoprecipitation with DNA sequencing was used to identify KLF5/GATA4/GATA6 genomic binding sites in GC cell lines, and integrated with transcriptomics to highlight direct target genes. Phenotypical assays were conducted to assess the function of these factors in GC cell lines and xenografts in nude mice. Results KLF5, GATA4 and GATA6 amplifications were confirmed in independent GC cohorts. Although factor amplifications occurred in distinct sets of GCs, they exhibited significant mRNA coexpression in primary GCs, consistent with KLF5/GATA4/GATA6 cross-regulation. Chromatin immunoprecipitation with DNA sequencing revealed a large number of genomic sites co-occupied by KLF5 and GATA4/GATA6, primarily located at gene promoters and exhibiting higher binding strengths. KLF5 physically interacted with GATA factors, supporting KLF5/GATA4/GATA6 cooperative regulation on co-occupied genes. Depletion and overexpression of these factors, singly or in combination, reduced and promoted cancer proliferation, respectively, in vitro and in vivo. Among the KLF5/GATA4/GATA6 direct target genes relevant for cancer development, one target gene, HNF4α, was also required for GC proliferation and could be targeted by the antidiabetic drug metformin, revealing a therapeutic opportunity for KLF5/GATA4/GATA6 amplified GCs. Conclusions KLF5/GATA4/GATA6 may promote GC development by engaging in mutual crosstalk, collaborating to maintain a pro-oncogenic transcriptional regulatory network in GC cells.

Nucleus-Translocated ACSS2 Promotes Gene Transcription for Lysosomal Biogenesis and Autophagy

Overcoming metabolic stress is a critical step in tumor growth. Acetyl coenzyme A (acetyl-CoA) generated from glucose and acetate uptake is important for histone acetylation and gene expression. However, how acetyl-CoA is produced under nutritional stress is unclear. We demonstrate here that glucose deprivation results in AMP-activated protein kinase (AMPK)-mediated acetyl-CoA synthetase 2 (ACSS2) phosphorylation at S659, which exposed the nuclear localization signal of ACSS2 for importin α5 binding and nuclear translocation. In the nucleus, ACSS2 binds to transcription factor EB and translocates to lysosomal and autophagy gene promoter regions, where ACSS2 incorporates acetate generated from histone acetylation turnover to locally produce acetyl-CoA for histone H3 acetylation in these regions and promote lysosomal biogenesis, autophagy, cell survival, and brain tumorigenesis. In addition, ACSS2 S659 phosphorylation positively correlates with AMPK activity in glioma specimens and grades of glioma malignancy. These results underscore the significance of nuclear ACSS2-mediated histone acetylation in maintaining cell homeostasis and tumor development.

Aristolochic acids and their derivatives are widely implicated in liver cancers in Taiwan and throughout Asia

Mutational signatures reveal high burdens of aristolochic acid–related mutations in Asian liver cancers, with Taiwan most intensely affected. The dark side of an herbal medicine Aristolochic acid, an herbal compound found in many traditional medicines, had been previously linked to kidney failure, as well as cancers of the urinary tract. Because of these known toxicities, herbs containing this compound have been restricted or banned in some countries, but it is still available on the internet and in alternate formulations. By analyzing numerous samples from Taiwan and other countries in Asia and elsewhere, Ng et al. demonstrated the effects of aristolochic acid in hepatocellular carcinoma, a much more common tumor type. The authors showed that the use of this drug remains widespread in Asia and particularly in Taiwan, and that it appears to increase the risk of multiple different cancer types. Many traditional pharmacopeias include Aristolochia 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.

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.