Impact of somatic mutations and transcriptomic alterations on cancer aneuploidy.

Aneuploidy has been recognized as one of hallmark of tumorigenesis since the early 20th century. Recent developments in structural variation analysis in the human genome have revealed the diversity of aneuploidy in cancer. However, the effects of gene mutation and expression in tumors on aneuploidy remain poorly understood. Here, we performed whole exome analysis of over 5,000 Japanese cancer cases and investigated the impact of somatic mutations and gene expression alterations on aneuploidy. First, we evaluated tumor content and genomic alterations that could influence aneuploidy. Next, we compared the aneuploidy frequency in 18 cancer types and observed that TP53 mutations were associated with the aneuploidy on specific chromosomes in colorectal and gastric cancers. Finally, we used expression analysis to isolate pathways involved in aneuploidy accumulation from tumors without TP53 mutations. Chromosomal instability and cell cycle aberration were associated with aneuploidy in TP53 wild-type tumors, and 26 commonly upregulated genes were identified in aneuploidy-high solid tumors without TP53 mutations. Among them, two cancer-related genes (CENPA and PBK) were involved in aneuploidy. Our integrated analysis revealed that both TP53 mutations and transcriptomic alterations independent of somatic mutations affect aneuploidy accumulation. Our findings will facilitate further understanding of diverse aneuploidies in the tumorigenesis.

Comparative analysis of tumor content estimation methods based on simu- lated tumor samples identified their impact on somatic variant detection in cancer whole genome sequencing.

Whole genome sequencing (WGS) in cancer genomics has become widespread with recent technological innovations, and the amount and types of information obtained from WGS are increasing rapidly. Appropriate interpretation of results is becoming increasingly important in clinical applications. This study aimed to evaluate the accuracy of tumor content estimation and its impact on somatic variant detection, using 100 simulated tumor samples covering 10-100% tumor content constructed from the sequencing data of cell line models. Extensive analysis revealed that the estimation results varied among computational analytical methods. Notably, there was a large discrepancy in low tumor content (≤ 30%). The reproducibility decreased in cases wherein chromosome-scale copy number changes were observed in normal cells. The minimum tumor content required to detect somatic alterations was estimated to be 10-30%. Identification of whole genome doubling was achieved with the lowest tumor content, followed by single nucleotide variation/insertion or deletion, structural variation, and copy number variation. Tumor content had a significantly higher impact on the false negatives than the false positives in variant calls. Results should be interpreted cautiously for samples wherein tumor content is a concern. These results can form the basis of developing important guidelines for evaluating cancer WGS.

Whole-genome and Epigenomic Landscapes of Malignant Gastrointestinal Stromal Tumors Harboring KIT Exon 11 557-558 Deletion Mutations.

Gastrointestinal stromal tumors (GIST) with KIT exon 11 deletions involving in codons 557-558 (KIT Δ557-558) exhibit higher proliferation rates and shorter disease-free survival times compared with GISTs with other KIT exon 11 mutations. We analyzed 30 GIST cases and observed genomic instability and global DNA hypomethylation only in high-risk malignant GISTs with KIT Δ557-558. Whole-genome sequencing revealed that the high-risk malignant GISTs with KIT Δ557-558 (12 cases) had more structural variations (SV), single-nucleotide variants, and insertions and deletions compared with the low-risk, less malignant GISTs with KIT Δ557-558 (six cases) and the high-risk (six cases) or low-risk (6 cases) GISTs with other KIT exon 11 mutations. The malignant GISTs with KIT Δ557-558 showed higher frequency and significance in copy number (CN) reduction on chromosome arms 9p and 22q, and 50% of them had LOH or CN-dependent expression reduction in CDKN2A. In addition, SVs with driver potential were detected in 75% of them, in which AKT3 and MGMT were recurrently identified. Genome-wide DNA methylation and gene expression analyses showed global intergenic DNA hypomethylation, SNAI2 upregulation, and higher expression signatures, including p53 inactivation and chromosomal instability, as characteristics of malignant GISTs with KIT Δ557-558 that distinguished them from other GISTs. These genomic and epigenomic profiling results revealed that KIT Δ557-558 mutations are associated with increased genomic instability in malignant GISTs. Significance: We present genomic and epigenomic insights into the malignant progression of GISTs with KIT exon 11 deletions involving in 557-558, demonstrating their unique chromosomal instability and global intergenic DNA hypomethylation.

Comparison of genetic susceptibility to lung adenocarcinoma and squamous cell carcinoma in Japanese patients using a novel panel for cancer-related drug-metabolizing enzyme genes.

The differences in genetic susceptibility to lung adenocarcinoma and squamous cell carcinoma remain unclear. We developed a customized, targeted gene sequencing panel for efficient and sensitive identification of germline variants, including whole-gene deletion types for cancer-related drug-metabolizing enzyme genes in lung adenocarcinoma and squamous cell carcinoma. The minor allele frequencies of the variants, confirmed as clinically significant in the Japanese population, did not differ significantly from those of normal participants listed in the public database. Genotype analysis comparing lung adenocarcinoma (n = 559) and squamous cell carcinoma (n = 151) indicated that the variants of DPYD (rs190771411, Fisher's exact test, P = 0.045; rs200562975, P = 0.045) and ALDH2 (rs568781254, P = 0.032) were associated with an increased risk of squamous cell carcinoma compared to adenocarcinoma. Conversely, whole-gene deletion of CYP2A6 was associated with adenocarcinoma but not squamous cell carcinoma. Notably, whole-gene deletion of CYP2A6 was confirmed in 22 patients with lung adenocarcinoma but not in any patients with squamous cell carcinoma. Most patients with whole-gene deletion of CYP2A6 were female non-smokers. The discovery of a whole-gene deletion of CYP2A6 in patients with lung adenocarcinoma may have an important role in clinical practice and advance our understanding of CYP2A6 germline variants and their association with carcinogenesis or their susceptibility to lung adenocarcinoma.

Development of two 410-cancer-gene panel tests for solid tumors and liquid biopsy based on genome data of 5,143 Japanese cancer patients.

Next-generation sequencing (NGS) is an integral part of precision medicine, and its power for detecting comprehensive genetic alterations may contribute to treatment decisions for patients with advanced, recurrent, or metastatic cancer. An NGS oncology panel developed in the U.S. and Europe, which targets cancer-related genes, has been approved in Japan, and testing is becoming more widespread in clinical oncology practice. However, these panels are based on cancer-related genes selected from cancer databases of Westerners. We aimed to develop an onco-panel for Japanese. We designed two High-tech Omics-based Patient Evaluation (HOPE) onco-panels: HOPE onco-panel Solid for solid tumors and HOPE onco-panel Liquid for liquid biopsy. These were based on genomic information of 5,143 cancer cases in the Japanese Cancer Genome Atlas (JCGA), a database of Japanese cancer cases. Their performance was confirmed using clinical data.

JCGA: the Japanese version of the Cancer Genome Atlas and its contribution to the interpretation of gene alterations detected in clinical cancer genome sequencing.

With the emergence of next-generation sequencing (NGS)-based cancer gene panel tests in routine oncological practice in Japan, an easily interpretable cancer genome database of Japanese patients in which mutational profiles are unaffected by racial differences is needed to improve the interpretation of the detected gene alterations. Considering this, we constructed the first Japanese cancer genome database, called the Japanese version of the Cancer Genome Atlas (JCGA), which includes multiple tumor types. The database includes whole-exome sequencing data from 4907 surgically resected primary tumor samples obtained from 4753 Japanese patients with cancer and graphically provides genome information on 460 cancer-associated genes, including the 336 genes that are included in two NGS-based cancer gene panel tests approved by the Pharmaceuticals and Medical Devices Agency. Moreover, most of the contents of this database are written in Japanese; this not only helps physicians explain the results of NGS-based cancer gene panel tests but also enables patients and their families to obtain further information regarding the detected gene alterations.

BMP4 and PHLDA1 are plausible drug-targetable candidate genes for KRAS G12A-, G12D-, and G12V-driven colorectal cancer.

Despite the frequent detection of KRAS driver mutations in patients with colorectal cancer (CRC), no effective treatments that target mutant KRAS proteins have been introduced into clinical practice. In this study, we identified potential effector molecules, based on differences in gene expression between CRC patients carrying wild-type KRAS (n = 390) and those carrying KRAS mutations in codon 12 (n = 240). CRC patients with wild-type KRAS harboring mutations in HRAS, NRAS, PIK3CA, PIK3CD, PIK3CG, RALGDS, BRAF, or ARAF were excluded from the analysis. At least 11 promising candidate molecules showed greater than two-fold change between the KRAS G12 mutant and wild-type and had a Benjamini-Hochberg-adjusted P value of less than 1E-08, evidence of significantly differential expression between these two groups. Among these 11 genes examined in cell lines transfected with KRAS G12 mutants, BMP4, PHLDA1, and GJB5 showed significantly higher expression level in KRAS G12A, G12D, and G12V transfected cells than in the wild-type transfected cells. We expect that this study will lead to the development of novel treatments that target signaling molecules functioning with KRAS G12-driven CRC.

Genomic profiling of multiple tissues in two patients with multiple endocrine neoplasia type 1.

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant tumor syndrome. This hereditary cancer is caused by germline variants in MEN1. Two patients with MEN1 were identified via whole exome sequencing and gene expression profile analysis, conducted for 5,063 patients with various types of cancers. We obtained multiple tumors from each patient; tumors derived from these two MEN1 patients had a loss of the normal MEN1 allele and frequently chromosomal copy number changes. Thus, we investigated whether structural variants were present in the MEN1 patient genomes. Whole-genome sequencing revealed no catastrophic rearrangements, and the tumor samples had very low somatic variants. The two patients had germline variants in MEN1 and some chromosomal copy number changes including on chromosome 11. The only pathogenic variant detected was the MEN1 germline variant, and chromosomal rearrangements led to tumorigenesis in somatic cells. Furthermore, the MEN1 tumor samples displayed a specific signature characterized by T:A>C:G transition. Studies of multiple tumors obtained from single patients are rare in hereditary cancer syndromes, and our results provide insights that the second hit of the tumor suppressor gene MEN1 may be caused by a gross genome rearrangement, not a small insertion and deletion, nor a change in epigenetic regulation.

Whole exome sequencing and deep sequencing of esophageal squamous cell carcinoma and adenocarcinoma in Japanese patients using the Japanese version of the Genome Atlas, JCGA.

BACKGROUND: Recent comprehensive mutation analyses have revealed a relatively small number of driver mutations in esophageal cancer, implicating a limited number of molecular targets, most of which are also implicated in squamous cell carcinoma. METHODS: In this study, we investigated genetic alterations in 44 esophageal squamous cell carcinomas (ESCC) and 8 adenocarcinomas (EAC) from Japanese patients as potential molecular targets, based on data from the Japanese version of The Genome Atlas (JCGA). RESULTS: Esophageal cancer was characterized by TP53 somatic mutations in ESCC (39/44, 88.6%) and EAC (5/8, 62.5%). In addition to TP53 mutations, somatic mutations in NFE2L2 (16/44, 36.4%), CDKN2A (7/44, 15.9%), and KMT2D (7/44, 15.9%) were more frequently detected in ESCC than in EAC. WRN-truncated type mutations that lead to genomic instability correlate with EAC, but not ESCC. ESCC samples were enriched in ALDH2-associated mutational signature 16 as well as the APOBEC signature. Patients with FAT2 mutations had significantly poorer overall survival compared with those with wild-type status at FAT2 (p < 0.05). Patients with EP300 or PTPRD mutations also had poor progression-free survival compared with respective wild-types (p < 0.05 or p < 0.001). CONCLUSIONS: These findings may facilitate future precision medicine approaches based on genomic profiling in ESCC and EAC.

Mutational concordance analysis provides supportive information for double cancer diagnosis.

BACKGROUND: Mutation analysis using next-generation sequencing highlights the features of tumors with somatic alterations. However, the mutation profile of double cancer remains unclear. Here, we analyzed tumors derived from the same patient using whole exome sequencing (WES) to investigate the coherence of somatic mutations in double cancer. METHODS: First, the tumor mutational burden (TMB) was investigated using WES of 5521 tumor specimens from a Japanese pan-cancer cohort. The frequencies of mutation concordance were then compared in these cancers. Finally, we calculated the expected value of mutational concordance fitting a Poisson distribution to determine the relationship between double and metastatic cancers. RESULTS: In all, 44, 58, and 121 paired samples were diagnosed as double cancer, multifocal lesions (derived from identical tissues), and metastasis, respectively. Our analysis revealed that common somatic mutations were almost entirely absent in double cancer, whereas primary tumors and metastatic foci harbored several identical alterations. Concordance of the mutation profile in the same patient reflects the tumor origin and development, suggesting the potential for identifying double cancer based on common somatic mutations. Furthermore, according to a Poisson distribution, double cancer could be discriminated based on paired samples from the same patient. The probability of double cancer with more than 10 mutations was ≤1 part-per-billion (ppb, 10- 9). In multifocal lesions, 74% of tumor pairs accumulated ≤10 common mutations, implying a difference in tumor origin within identical tissues. CONCLUSIONS: These findings indicate that counting common somatic mutations can indicate the differences in origin between tumors derived from the same patient. Our mutation coherence analysis can thus provide beneficial information for diagnosing double cancer.

Characterization of tumors with ultralow tumor mutational burden in Japanese cancer patients.

Tumor mutational burden analysis using whole-exome sequencing highlights features of tumors with various mutations or known driver alterations. Cancers with few changes in the exon regions have unclear characteristics, even though low-mutated tumors are often detected in pan-cancer analysis. In the present study, we analyzed tumors with low tumor mutational burden listed in the Japanese version of The Cancer Genome Atlas, a data set of 5020 primary solid tumors. Our analysis revealed that detection rates of known driver mutations and copy number variation were decreased in samples with tumor mutational burden below 1.0 (ultralow tumor), compared with those in samples with low tumor mutational burden (≤5 mutations/Mb). This trend was also observed in The Cancer Genome Atlas data set. In the ultralow tumor mutational burden tumors, expression analysis showed decreased TP53 inactivation and chromosomal instability. TP53 inactivation frequently correlated with PI3K/mTOR-related gene expression, implying suppression of the PI3K/mTOR pathway in ultralow tumor mutational burden tumors. In common with mutational burden, the T cell-inflamed gene expression profiling signature was a potential marker for prediction of an immune checkpoint inhibitor response, and some ultralow tumor mutational burden tumor populations highly expressed this signature. Our analysis focused on how these tumors could provide insight into tumors with low somatic alteration that are difficult to detect solely using whole-exome sequencing.

Japanese version of The Cancer Genome Atlas, JCGA, established using fresh frozen tumors obtained from 5143 cancer patients.

This study aimed to establish the Japanese Cancer Genome Atlas (JCGA) using data from fresh frozen tumor tissues obtained from 5143 Japanese cancer patients, including those with colorectal cancer (31.6%), lung cancer (16.5%), gastric cancer (10.8%) and other cancers (41.1%). The results are part of a single-center study called "High-tech Omics-based Patient Evaluation" or "Project HOPE" conducted at the Shizuoka Cancer Center, Japan. All DNA samples and most RNA samples were analyzed using whole-exome sequencing, cancer gene panel sequencing, fusion gene panel sequencing and microarray gene expression profiling, and the results were annotated using an analysis pipeline termed "Shizuoka Multi-omics Analysis Protocol" developed in-house. Somatic driver alterations were identified in 72.2% of samples in 362 genes (average, 2.3 driver events per sample). Actionable information on drugs that is applicable in the current clinical setting was associated with 11.3% of samples. When including those drugs that are used for investigative purposes, actionable information was assigned to 55.0% of samples. Germline analysis revealed pathogenic mutations in hereditary cancer genes in 9.2% of samples, among which 12.2% were confirmed as pathogenic mutations by confirmatory test. Pathogenic mutations associated with non-cancerous hereditary diseases were detected in 0.4% of samples. Tumor mutation burden (TMB) analysis revealed 5.4% of samples as having the hypermutator phenotype (TMB ≥ 20). Clonal hematopoiesis was observed in 8.4% of samples. Thus, the JCGA dataset and the analytical procedures constitute a fundamental resource for genomic medicine for Japanese cancer patients.

Driver gene alterations and activated signaling pathways toward malignant progression of gastrointestinal stromal tumors.

Mutually exclusive KIT and PDGFRA mutations are considered to be the earliest events in gastrointestinal stromal tumors (GIST), but insufficient for their malignant progression. Herein, we aimed to identify driver genes and signaling pathways relevant to GIST progression. We investigated genetic profiles of 707 driver genes, including mutations, gene fusions, copy number gain or loss, and gene expression for 65 clinical specimens of surgically dissected GIST, consisting of six metastatic tumors and 59 primary tumors from stomach, small intestine, rectum, and esophagus. Genetic alterations included oncogenic mutations and amplification-dependent expression enhancement for oncogenes (OG), and loss of heterozygosity (LOH) and expression reduction for tumor suppressor genes (TSG). We assigned activated OG and inactivated TSG to 27 signaling pathways, the activation of which was compared between malignant GIST (metastasis and high-risk GIST) and less malignant GIST (low- and very low-risk GIST). Integrative molecular profiling indicated that a greater incidence of genetic alterations of driver genes was detected in malignant GIST (96%, 22 of 23) than in less malignant GIST (73%, 24 of 33). Malignant GIST samples groups showed mutations, LOH, and aberrant expression dominantly in driver genes associated with signaling pathways of PI3K (PIK3CA, AKT1, and PTEN) and the cell cycle (RB1, CDK4, and CDKN1B). Additionally, we identified potential PI3K-related genes, the expression of which was upregulated (SNAI1 and TPX2) or downregulated (BANK1) in malignant GIST. Based on our observations, we propose that inhibition of PI3K pathway signals might potentially be an effective therapeutic strategy against malignant progression of GIST.

Mutational burden and signatures in 4000 Japanese cancers provide insights into tumorigenesis and response to therapy.

Tumor mutational burden (TMB) and mutational signatures reflect the process of mutation accumulation in cancer. However, the significance of these emerging characteristics remains unclear. In the present study, we used whole-exome sequencing to analyze the TMB and mutational signature in solid tumors of 4046 Japanese patients. Eight predominant signatures-microsatellite instability, smoking, POLE, APOBEC, UV, mismatch repair, double-strand break repair, and Signature 16-were observed in tumors with TMB higher than 1.0 mutation/Mb, whereas POLE and UV signatures only showed moderate correlation with TMB, suggesting the extensive accumulation of mutations due to defective POLE and UV exposure. The contribution ratio of Signature 16, which is associated with hepatocellular carcinoma in drinkers, was increased in hypopharynx cancer. Tumors with predominant microsatellite instability signature were potential candidates for treatment with immune checkpoint inhibitors such as pembrolizumab and were found in 2.8% of cases. Furthermore, based on microarray analysis, tumors with predominant signatures were classified into 2 subgroups depending on the expression of immune-related genes reflecting differences in the immune context of the tumor microenvironment. Tumor subpopulations differing in the content of infiltrating immune cells might respond differently to immunotherapeutics. An understanding of cancer characteristics based on TMB and mutational signatures could provide new insights into mutation-driven tumorigenesis.

Molecular profiling and sequential somatic mutation shift in hypermutator tumours harbouring POLE mutations.

Defective DNA polymerase ε (POLE) proofreading leads to extensive somatic mutations that exhibit biased mutational properties; however, the characteristics of POLE-mutated tumours remain unclear. In the present study, we describe a molecular profile using whole exome sequencing based on the transition of somatic mutations in 10 POLE-mutated solid tumours that were obtained from 2,042 Japanese patients. The bias of accumulated variations in these mutants was quantified to follow a pattern of somatic mutations, thereby classifying the sequential mutation shift into three periods. During the period prior to occurrence of the aberrant POLE, bare accumulation of mutations in cancer-related genes was observed, whereas PTEN was highly mutated in conjunction with or subsequent to the event, suggesting that POLE and PTEN mutations were responsible for the development of POLE-mutated tumours. Furthermore, homologous recombination was restored following the occurrence of PTEN mutations. Our strategy for estimation of the footprint of somatic mutations may provide new insight towards the understanding of mutation-driven tumourigenesis.

Tumor mutational burden analysis of 2,000 Japanese cancer genomes using whole exome and targeted gene panel sequencing.

Tumor mutational burden (TMB) is an emerging characteristic in cancer and has been associated with microsatellite instability, defective DNA replication/repair, and response to PD-1 and PD-L1 blockade immunotherapy. When estimating TMB, targeted panel sequencing is performed using a few hundred genes; however, a comparison of TMB results obtained with this platform and with whole exome sequencing (WES) has not been performed for various cancer types. In the present study, we compared TMB results using the above two platforms in 2,908 solid tumors that were obtained from Japanese patients. For next-generation sequencing, we used fresh-frozen tissue specimens. The Ion Proton System was employed to detect somatic mutations in the coding genome and to sequence an available cancer panel that targeted 409 genes. We then selected 2,040 samples with sufficient tumor cellularity for TMB analysis. In tumors with TMB-high (TMB ≥ 20 mutations/Mb), TMB derived from WES correlated well with the estimated TMB (eTMB) based on panel sequencing, whereas TMB in the remaining tumors showed a weak correlation. In particular, eTMB was overestimated in tumors with low-frequency mutations, resulting in the accumulation of EGFR mutations not being discriminated as a feature of lung cancer with low-frequency mutations. The eTMB in tumors harboring POLE mutations and microsatellite instability was not overestimated, suggesting that panel sequencing could accurately estimate TMB in tumors with high-frequency mutations such as hypermutator tumors. These results may provide helpful information for interpreting TMB results based on clinical sequencing using a targeted gene panel.

Germline and somatic genetic changes in multicentric tumors obtained from a patient with multiple endocrine neoplasia type 1.

Multiple endocrine neoplasia type 1 (MEN1) is a hereditary cancer syndrome caused by germline mutations of the MEN1 gene located in chromosome 11q13. In patients with MEN1, multicentric tumors develop in the involved organs; however, precise evaluation of genetic changes in these multicentric tumors has not been performed. In the present study, using whole-exome sequencing, we analyzed germline and somatic genetic changes in blood cells, two pancreatic endocrine tumors and one duodenal tumor obtained from a patient with MEN1 gastrinoma. We found that this patient possessed a novel germline mutation of the MEN1 gene [NM_137099.2:c.1505dupA (p.Lys502Lysfs); the localization was Chr11:64572134 on Assembly GRCh37], in which an adenine insertion in codon 502 of the MEN1 gene resulted in a frame shift and a premature stop codon. In terms of heterozygosity, the mutated allele was heterozygous in blood cells, hemizygous in the two pancreatic tumors and homozygous in the duodenal tumor. Immunohistochemical staining confirmed that only truncated menin protein accumulated in the nucleus of the tumor tissues. Further evaluation of tumor-specific somatic mutations in two pancreatic tumors did not detect single-nucleotide variations (SNVs) in 609 cancer-associated genes designated by the COSMIC cancer gene census, suggesting that the germline MEN1 mutation and resultant loss of heterozygosity played a major role in tumorigenesis. In the duodenal tumor, in addition to the germline MEN1 mutation, single-nucleotide variations in two cancer-associated genes were found. Further studies are required to clarify the role of these somatic single-nucleotide variations in the progression of MEN1 tumors.

Integrated analysis of gene expression and copy number identified potential cancer driver genes with amplification-dependent overexpression in 1,454 solid tumors.

Identification of driver genes contributes to the understanding of cancer etiology and is imperative for the development of individualized therapies. Gene amplification is a major event in oncogenesis. Driver genes with tumor-specific amplification-dependent overexpression can be therapeutic targets. In this study, we aimed to identify amplification-dependent driver genes in 1,454 solid tumors, across more than 15 cancer types, by integrative analysis of gene expression and copy number. Amplification-dependent overexpression of 64 known driver oncogenes were found in 587 tumors (40%); genes frequently observed were MYC (25%) and MET (18%) in colorectal cancer; SKP2 (21%) in lung squamous cell carcinoma; HIST1H3B (19%) and MYCN (13%) in liver cancer; KIT (57%) in gastrointestinal stromal tumors; and FOXL2 (12%) in squamous cell carcinoma across tissues. Genomic aberrations in 138 known cancer driver genes and 491 established fusion genes were found in 1,127 tumors (78%). Further analyses of 820 cancer-related genes revealed 16 as potential driver genes, with amplification-dependent overexpression restricted to the remaining 22% of samples (327 tumors) initially undetermined genetic drivers. Among them, AXL, which encodes a receptor tyrosine kinase, was recurrently overexpressed and amplified in sarcomas. Our studies of amplification-dependent overexpression identified potential drug targets in individual tumors.

Comprehensive characterization of genes associated with the TP53 signal transduction pathway in various tumors.

The TP53 signal transduction pathway is an attractive target for cancer treatments. In this study, we conducted a comprehensive molecular evaluation of 907 patients with cancer in Japan to identify genomic alterations in the TP53 pathway. TP53 mutations were frequently detected in many cancers, except melanoma, thymic tumors, gastrointestinal stromal tumors, and renal cancers. The frequencies of non-synonymous single nucleotide variants (SNVs) in the TP53 family members TP63 and TP73 were relatively low, although genes with increased frequencies of SNVs were as follows: PTEN (11.7%) in breast cancer, CDKN2A (11.1 and 9.6%) in pancreas and head and neck cancers, and ATM (18.0 and 11.1%) in liver and esophageal cancers. MDM2 expression was decreased or increased in patients with mutant or wild-type TP53, respectively. CDKN1A expression was increased with mutant TP53 in head and neck cancers. Moreover, TP63 overexpression was characteristically observed in squamous cell carcinomas of the lung, esophagus, and head and neck region. Additionally, overexpression of TP63 and TP73 was frequently observed in thymomas. Our results reveal a spectrum of genomic alterations in the TP53 pathway that is characteristic of many tumor types, and these data may be useful in the trials of targeted therapies.

Whole exome sequencing detects variants of genes that mediate response to anticancer drugs.

Certain interindividual differences affecting the efficacy of drug treatment and adverse drug reactions are caused by genetic variants, and their phenotypic effects differ among ethnic groups. In this study, we used whole exome sequencing (WES) systematically to identify germline mutations that influence the activities of drug-metabolizing enzymes, as well as that of a transporter. We analyzed DNA isolated from blood samples from 2,042 Japanese patients with diverse cancers. We identified sequence variants of CYP2B6 (rs3745274), CYP2C9 (rs1057910), CYP2C19 (rs4986893), CYP2C19 (rs4244285), TPMT (rs1142345), NAT2 (rs1799930), NAT2 (rs1799931), UGT1A1 (rs4148323), COMT (rs4680), ABCB1 (rs1045642), and CDA (rs60369023). Wider application of WES will help to determine the effects of mutations on the activities of proteins encoded by drug response genes, and the information gained will accelerate the development of personalized therapies for patients with cancer. Moreover, this knowledge may provide clues for preventing cancer before the onset of symptoms.