Anna Gan

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.

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.

Exome-wide Sequencing Shows Low Mutation Rates and Identifies Novel Mutated Genes in Seminomas

BACKGROUND Testicular germ cell tumors are the most common cancer diagnosed in young men, and seminomas are the most common type of these cancers. There have been no exome-wide examinations of genes mutated in seminomas or of overall rates of nonsilent somatic mutations in these tumors. OBJECTIVE The objective was to analyze somatic mutations in seminomas to determine which genes are affected and to determine rates of nonsilent mutations. DESIGN, SETTING, AND PARTICIPANTS Eight seminomas and matched normal samples were surgically obtained from eight patients. INTERVENTION DNA was extracted from tissue samples and exome sequenced on massively parallel Illumina DNA sequencers. Single-nucleotide polymorphism chip-based copy number analysis was also performed to assess copy number alterations. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS The DNA sequencing read data were analyzed to detect somatic mutations including single-nucleotide substitutions and short insertions and deletions. The detected mutations were validated by independent sequencing and further checked for subclonality. RESULTS AND LIMITATIONS The rate of nonsynonymous somatic mutations averaged 0.31 mutations/Mb. We detected nonsilent somatic mutations in 96 genes that were not previously known to be mutated in seminomas, of which some may be driver mutations. Many of the mutations appear to have been present in subclonal populations. In addition, two genes, KIT and KRAS, were affected in two tumors each with mutations that were previously observed in other cancers and are presumably oncogenic. CONCLUSIONS Our study, the first report on exome sequencing of seminomas, detected somatic mutations in 96 new genes, several of which may be targetable drivers. Furthermore, our results show that seminoma mutation rates are five times higher than previously thought, but are nevertheless low compared to other common cancers. Similar low rates are seen in other cancers that also have excellent rates of remission achieved with chemotherapy. PATIENT SUMMARY We examined the DNA sequences of seminomas, the most common type of testicular germ cell cancer. Our study identified 96 new genes in which mutations occurred during seminoma development, some of which might contribute to cancer development or progression. The study also showed that the rates of DNA mutations during seminoma development are higher than previously thought, but still lower than for other common solid-organ cancers. Such low rates are also observed among other cancers that, like seminomas, show excellent rates of disease remission after chemotherapy.

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 implicate aristolochic acid in bladder cancer development

BackgroundAristolochic 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.MethodsHere, 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.ResultsThe 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.ConclusionsThis 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.

Multiregion ultra‐deep sequencing reveals early intermixing and variable levels of intratumoral heterogeneity in colorectal cancer

Intratumor heterogeneity (ITH) contributes to cancer progression and chemoresistance. We sought to comprehensively describe ITH of somatic mutations, copy number, and transcriptomic alterations involving clinically and biologically relevant gene pathways in colorectal cancer (CRC). We performed multiregion, high‐depth (384× on average) sequencing of 799 cancer‐associated genes in 24 spatially separated primary tumor and nonmalignant tissues from four treatment‐naïve CRC patients. We then used ultra‐deep sequencing (17 075× on average) to accurately verify the presence or absence of identified somatic mutations in each sector. We also digitally measured gene expression and copy number alterations using NanoString assays. We identified the subclonal point mutations and determined the mutational timing and phylogenetic relationships among spatially separated sectors of each tumor. Truncal mutations, those shared by all sectors in the tumor, affected the well‐described driver genes such as APC, TP53, and KRAS. With sequencing at 17 075×, we found that mutations first detected at a sequencing depth of 384× were in fact more widely shared among sectors than originally assessed. Interestingly, ultra‐deep sequencing also revealed some mutations that were present in all spatially dispersed sectors, but at subclonal levels. Ultra‐high‐depth validation sequencing, copy number analysis, and gene expression profiling provided a comprehensive and accurate genomic landscape of spatial heterogeneity in CRC. Ultra‐deep sequencing allowed more sensitive detection of somatic mutations and a more accurate assessment of ITH. By detecting the subclonal mutations with ultra‐deep sequencing, we traced the genomic histories of each tumor and the relative timing of mutational events. We found evidence of early mixing, in which the subclonal ancestral mutations intermixed across the sectors before the acquisition of subsequent nontruncal mutations. Our findings also indicate that different CRC patients display markedly variable ITH, suggesting that each patients tumor possesses a unique genomic history and spatial organization.