Impact of Mutations in Subunit Genes of the Mammalian SWI/SNF Complex on Immunological Tumor Microenvironment.

BACKGROUND/AIM: Recently, inactivating somatic mutations of SWI/SNF chromatin-remodeling genes in cancers have been reported. However, few studies have been performed regarding the immunological analysis of the tumor microenvironment (TME) in chromatin remodeling complex gene-mutated tumors. In the present study, we identified cancer patients harboring various mammalian SWI/SNF complex mutations and investigated the immunological features in those mutated cancers. PATIENTS AND METHODS: Cancer patients harboring any type of chromatin remodeling complex gene mutation were selected and clinicopathological features were compared between chromatin remodeling complex gene expression-low and expression-high groups. Specifically, expression levels of immune response-associated genes and cancer-associated genes were compared between the SMARCA4 expression-low and expression-high groups using volcano plot analysis. RESULTS: Among cancers harboring PBRM1, SAMRACA4 and ARID2 gene mutations, T-cell marker and mature B-cell marker genes were up-regulated in the tumor. Specifically, T-cell effector genes (CD8B, CD40LG), central memory marker genes (CD27, CCR7) and mature B-cell marker genes (CD20, CD38, CD79 and IRF4) were up-regulated, and cancer-associated genes including MYB, MYC and AURKB genes were down-regulated in the SMARCA4 expression-low group. Remarkably, heatmap of gene expression and immunohistochemistry (IHC) data demonstrated that the tertiary lymphoid structure (TLS) gene signature of mature B cells was up-regulated in SMACA4 gene-mutated stomach cancers. CONCLUSION: These results suggest that immune tumor microenvironment status, such as mature B cell recruitment featuring the TLS gene signature and immune activation mediated by cancer signal down-regulation, might contribute to the classification of SMARCA4 gene-mutated tumors as immune checkpoint blockade therapy-sensitive target tumors.

An autopsy case of primary gliosarcoma with multiple extracranial metastases: pathology after administration of bevacizumab and genetic profile.

Gliosarcoma (GS), a morphological variant of glioblastoma, pathologically shows a biphasic pattern with gliomatous and sarcomatous components. It has been reported that GS has much higher metastatic capacity than glioblastoma. A few reports on the pathology of the extracranial metastasis of GS have shown that metastatic lesions had a sarcomatous component alone or a mixture of gliomatous and sarcomatous ones. Therefore, it is considered that GS tends to disseminate hematogenously due to its mesenchymal sarcomatous component. Herein, we report an autopsy case of GS with multiple extracranial metastases treated by craniotomy, radiotherapy, and bevacizumab. In this case, metastatic lesions at autopsy contained a gliomatous component alone, but no sarcomatous component. In addition, the sarcomatous component disappeared from the intracranial lesion at autopsy after the administration of bevacizumab. In this report, we discuss the clinical course and pathological findings at the initial state, recurrence, and autopsy, including the results of whole-genome analysis.

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.

A Primary Intraosseous Meningioma: A Rare Case of Malignancy with High Proliferative Ability.

Primary intraosseous meningioma (PIM) is a rare tumor that arises in the skull. Histopathologically, it is generally described as a slow-growing, benign lesion. However, on rare occasions, PIM presents as a malignancy with high proliferative ability, which requires maximal resection, adjuvant radiotherapy, and subsequent careful follow-up. Because of the rarity of such cases, they present a diagnostic challenge with unusual pathological findings. Herein, we report a case of a primary intraosseous anaplastic meningioma with extensive invasion inside and outside the skull, along with the results of whole-genome analysis. Histopathological diagnosis was a World Health Organization grade 3 anaplastic meningioma. In the literature, only two cases of anaplastic PIM have been reported, so its characteristics and treatment are poorly understood. Our patient was successfully treated with tumor resection, followed by intensity-modulated radiation therapy. Follow-up imaging studies revealed no recurrence or distant metastasis, including to lung, liver, and bone, at 8 months after the surgery.

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.

Generation of novel complete HLA class I monoallelic cell lines used in an MHC stabilization assay for neoantigen evaluation.

Immunogenic neoantigens derived from somatic mutations in cancer have been identified through clinical studies with the cloning of tumor-infiltrating T cells, and cancer driver gene mutation-derived epitopes have been reported; however, these are rare. At present, the validation of epitopes predicted in silico is difficult as human T-cell clonal diversity cannot be reproduced in vitro or in experimental animal models. To confirm the epitope peptides presented by human leukocyte antigen (HLA) class I molecules predicted in silico, biochemical methods such as major histocompatibility complex (MHC) stabilization assays and mass spectrometry-mediated identification have been developed based on HLA-A*02:01 monoallelic T2 cells and HLA-C*01:02 monoallelic LCL721.221 cells. Therefore, in the present study, to prevent confusion due to peptide cross-presentation among HLA molecules, HLA class I monoallelic B-cell clones were generated from the TISI cell line by knocking out HLA-ABC and TAP2, and knocking in HLA alleles. To explore cancer driver mutations as potential targets for immunotherapy, exome sequencing data from 5,143 patients with cancer enrolled in a comprehensive genome analysis project at the Shizuoka Cancer Center were used to identify somatic amino acid substituted mutations and the 50 most frequent mutations in five genes, TP53, EGFR, PIK3CA, KRAS and BRAF, were identified. Using NetMHC4.1, the present study predicted whether epitopes derived from these mutations are presented on major HLA-ABC alleles in Japanese individuals and synthesized 138 peptides for MHC stabilization assays. The authors also attempted to examine the candidate epitopes at physiological temperatures by using antibody clone G46-2.6, which can detect HLA-ABC, independent of ß2-microglobulin association. In the assays, although the peptide-induced HLA expression levels were associated with the predicted affinities, the respective HLA alleles exhibited varying degrees of responsiveness, and unexpectedly, p53-mutant epitopes with predicted weak affinities exhibited strong responses. These results suggested that MHC stabilization assays using completely monoallelic HLA-expressing B-cell lines are useful for evaluating the presentation of neoantigen epitopes.

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.

Characterization of the Immunological Status of Hypermutated Solid Tumors in the Cancer Genome Analysis Project HOPE.

BACKGROUND: Many reports demonstrate that a high tumor mutation burden (TMB-H) is closely associated with good prognosis of cancer. However, specific studies investigating the association of various TMB statuses with overall survival in patients with solid tumors are scarce. PATIENTS AND METHODS: In the present study, we investigated the association of TMB status with overall survival in 5,072 patients with cancer from the HOPE project and clarified the specific mechanism responsible for the good prognosis of the TMB-H group. All tumors were classified into one of four groups based on TMB: ultralow (UL), low (L), intermediate (I) and high (H). RESULTS: The TMB-H group had a better prognosis than the TMB-I and TMB-L groups, but not than the TMB-UL group. Analyzing the expression of 293 immune response-associated genes, 17 genes were up-regulated in the TMB-H group compared to the TMB-I and TNB-L groups, and two genes [CD274 and interferon-γ (IFNG)] were identified as good prognostic factors. Analysis of immune cell populations inside tumors demonstrated that the frequencies of exhausted CD8+ T-cells, activated effector CD8+ T-cells and natural killer cells were significantly higher in the TMB-H group. The T-cell receptor repertoire numbers and the diversity evenness score (DE50) were lower in the TMB-H group than in TMB-UL group; however, no association of the DE50 value with the binding or elution affinity of epitope peptides from neoantigens was found. CONCLUSION: One possible mechanism for the good prognosis of the TMB-UL group compared to the TMB-H group might be that the TMB-UL group features a balance between immunosuppression and immunostimulation.

Tumor cell enrichment by tissue suspension enables detection of mutations with low variant allele frequency and estimation of germline mutations.

Targeted sequencing offers an opportunity to select specific drugs for cancer patients based on alterations in their genome. However, accurate sequencing cannot be performed in cancers harboring diffuse tumor cells because of low tumor content. We performed tumor cell enrichment using tissue suspension of formalin-fixed, paraffin-embedded (FFPE) tissue sections with low tumor cell content. The enriched fractions were used to efficiently identify mutations by sequencing a target panel of cancer-related genes. Tumor-enriched and residual fractions were isolated from FFPE tissue sections of intestinal and diffuse gastric cancers harboring diffuse tumor cells and DNA of suitable quality was isolated for next-generation sequencing. Sequencing of a target panel of cancer-related genes using the tumor-enriched fraction increased the number of detectable mutations and variant allele frequency. Furthermore, mutation analysis of DNA isolated from tumor-enriched and residual fractions allowed us to estimate germline mutations without a blood reference. This approach of tumor cell enrichment will not only enhance the success rate of target panel sequencing, but can also improve the accuracy of detection of somatic mutations in archived specimens.

Identification of tumor microenvironment-associated immunological genes as potent prognostic markers in the cancer genome analysis project HOPE.

Project High-tech Omics-based Patient Evaluation (HOPE), which used whole-exome sequencing and gene expression profiling, was launched in 2014. A total of ~2,000 patients were enrolled until March 2016, and the survival time was observed up to July 2019. In our previous study, a tumor microenvironment immune type classification based on the expression levels of the programmed death-ligand 1 (PD-L1) and CD8B genes was performed based on four types: A, adaptive immune resistance; B, intrinsic induction; C, immunological ignorance; and D, tolerance. Type A (PD-L1+ and CD8B+) exhibited upregulated features of T helper 1 antitumor responses. In the present study, survival time analysis at 5 years revealed that patients in type A had a better prognosis than those in other categories [5 year survival rate (%); A (80.5) vs. B (73.9), C (73.4) and D (72.6), P=0.0005]. Based on the expression data of 293 immune response-associated genes, 62 specific genes were upregulated in the type A group. Among these genes, 18 specific genes, such as activated effector T-cell markers (CD8/CD40LG/GZMB), effector memory T-cell markers (PD-1/CD27/ICOS), chemokine markers (CXCL9/CXCL10) and activated dendritic cell markers (CD80/CD274/SLAMF1), were significantly associated with a good prognosis using overall survival time analysis. Finally, multivariate Cox proportional hazard regression analyses of overall survival demonstrated that four genes (GZMB, HAVCR2, CXCL9 and CD40LG) were independent prognostic markers, and GZMB, CXCL9 and CD40LG may contribute to the survival benefit of patients in the immune type A group.

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.

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.