Evolutionary histories of breast cancer and related clones.

Recent studies have documented frequent evolution of clones carrying common cancer mutations in apparently normal tissues, which are implicated in cancer development1-3. However, our knowledge is still missing with regard to what additional driver events take place in what order, before one or more of these clones in normal tissues ultimately evolve to cancer. Here, using phylogenetic analyses of multiple microdissected samples from both cancer and non-cancer lesions, we show unique evolutionary histories of breast cancers harbouring der(1;16), a common driver alteration found in roughly 20% of breast cancers. The approximate timing of early evolutionary events was estimated from the mutation rate measured in normal epithelial cells. In der(1;16)(+) cancers, the derivative chromosome was acquired from early puberty to late adolescence, followed by the emergence of a common ancestor by the patient's early 30s, from which both cancer and non-cancer clones evolved. Replacing the pre-existing mammary epithelium in the following years, these clones occupied a large area within the premenopausal breast tissues by the time of cancer diagnosis. Evolution of multiple independent cancer founders from the non-cancer ancestors was common, contributing to intratumour heterogeneity. The number of driver events did not correlate with histology, suggesting the role of local microenvironments and/or epigenetic driver events. A similar evolutionary pattern was also observed in another case evolving from an AKT1-mutated founder. Taken together, our findings provide new insight into how breast cancer evolves.

Two Aldehyde Clearance Systems Are Essential to Prevent Lethal Formaldehyde Accumulation in Mice and Humans.

Reactive aldehydes arise as by-products of metabolism and are normally cleared by multiple families of enzymes. We find that mice lacking two aldehyde detoxifying enzymes, mitochondrial ALDH2 and cytoplasmic ADH5, have greatly shortened lifespans and develop leukemia. Hematopoiesis is disrupted profoundly, with a reduction of hematopoietic stem cells and common lymphoid progenitors causing a severely depleted acquired immune system. We show that formaldehyde is a common substrate of ALDH2 and ADH5 and establish methods to quantify elevated blood formaldehyde and formaldehyde-DNA adducts in tissues. Bone-marrow-derived progenitors actively engage DNA repair but also imprint a formaldehyde-driven mutation signature similar to aging-associated human cancer mutation signatures. Furthermore, we identify analogous genetic defects in children causing a previously uncharacterized inherited bone marrow failure and pre-leukemic syndrome. Endogenous formaldehyde clearance alone is therefore critical for hematopoiesis and in limiting mutagenesis in somatic tissues.

Landscape of driver mutations and their clinical effects on Down syndrome-related myeloid neoplasms.

ABSTRACT: Transient abnormal myelopoiesis (TAM) is a common complication in newborns with Down syndrome (DS). It commonly progresses to myeloid leukemia (ML-DS) after spontaneous regression. In contrast to the favorable prognosis of primary ML-DS, patients with refractory/relapsed ML-DS have poor outcomes. However, the molecular basis for refractoriness and relapse and the full spectrum of driver mutations in ML-DS remain largely unknown. We conducted a genomic profiling study of 143 TAM, 204 ML-DS, and 34 non-DS acute megakaryoblastic leukemia cases, including 39 ML-DS cases analyzed by exome sequencing. Sixteen novel mutational targets were identified in ML-DS samples. Of these, inactivations of IRX1 (16.2%) and ZBTB7A (13.2%) were commonly implicated in the upregulation of the MYC pathway and were potential targets for ML-DS treatment with bromodomain-containing protein 4 inhibitors. Partial tandem duplications of RUNX1 on chromosome 21 were also found, specifically in ML-DS samples (13.7%), presenting its essential role in DS leukemia progression. Finally, in 177 patients with ML-DS treated following the same ML-DS protocol (the Japanese Pediatric Leukemia and Lymphoma Study Group acute myeloid leukemia -D05/D11), CDKN2A, TP53, ZBTB7A, and JAK2 alterations were associated with a poor prognosis. Patients with CDKN2A deletions (n = 7) or TP53 mutations (n = 4) had substantially lower 3-year event-free survival (28.6% vs 90.5%; P < .001; 25.0% vs 89.5%; P < .001) than those without these mutations. These findings considerably change the mutational landscape of ML-DS, provide new insights into the mechanisms of progression from TAM to ML-DS, and help identify new therapeutic targets and strategies for ML-DS.

Landscape and function of multiple mutations within individual oncogenes.

Sporadic reports have described cancer cases in which multiple driver mutations (MMs) occur in the same oncogene1,2. However, the overall landscape and relevance of MMs remain elusive. Here we carried out a pan-cancer analysis of 60,954 cancer samples, and identified 14 pan-cancer and 6 cancer-type-specific oncogenes in which MMs occur more frequently than expected: 9% of samples with at least one mutation in these genes harboured MMs. In various oncogenes, MMs are preferentially present in cis and show markedly different mutational patterns compared with single mutations in terms of type (missense mutations versus in-frame indels), position and amino-acid substitution, suggesting a cis-acting effect on mutational selection. MMs show an overrepresentation of functionally weak, infrequent mutations, which confer enhanced oncogenicity in combination. Cells with MMs in the PIK3CA and NOTCH1 genes exhibit stronger dependencies on the mutated genes themselves, enhanced downstream signalling activation and/or greater sensitivity to inhibitory drugs than those with single mutations. Together oncogenic MMs are a relatively common driver event, providing the underlying mechanism for clonal selection of suboptimal mutations that are individually rare but collectively account for a substantial proportion of oncogenic mutations.

Frequent mutations that converge on the NFKBIZ pathway in ulcerative colitis.

Chronic inflammation is accompanied by recurring cycles of tissue destruction and repair and is associated with an increased risk of cancer1-3. However, how such cycles affect the clonal composition of tissues, particularly in terms of cancer development, remains unknown. Here we show that in patients with ulcerative colitis, the inflamed intestine undergoes widespread remodelling by pervasive clones, many of which are positively selected by acquiring mutations that commonly involve the NFKBIZ, TRAF3IP2, ZC3H12A, PIGR and HNRNPF genes and are implicated in the downregulation of IL-17 and other pro-inflammatory signals. Mutational profiles vary substantially between colitis-associated cancer and non-dysplastic tissues in ulcerative colitis, which indicates that there are distinct mechanisms of positive selection in both tissues. In particular, mutations in NFKBIZ are highly prevalent in the epithelium of patients with ulcerative colitis but rarely found in both sporadic and colitis-associated cancer, indicating that NFKBIZ-mutant cells are selected against during colorectal carcinogenesis. In further support of this negative selection, we found that tumour formation was significantly attenuated in Nfkbiz-mutant mice and cell competition was compromised by disruption of NFKBIZ in human colorectal cancer cells. Our results highlight common and discrete mechanisms of clonal selection in inflammatory tissues, which reveal unexpected cancer vulnerabilities that could potentially be exploited for therapeutics in colorectal cancer.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

Germ line DDX41 mutations define a unique subtype of myeloid neoplasms.

Germ line DDX41 variants have been implicated in late-onset myeloid neoplasms (MNs). Despite an increasing number of publications, many important features of DDX41-mutated MNs remain to be elucidated. Here we performed a comprehensive characterization of DDX41-mutated MNs, enrolling a total of 346 patients with DDX41 pathogenic/likely-pathogenic (P/LP) germ line variants and/or somatic mutations from 9082 MN patients, together with 525 first-degree relatives of DDX41-mutated and wild-type (WT) patients. P/LP DDX41 germ line variants explained ∼80% of known germ line predisposition to MNs in adults. These risk variants were 10-fold more enriched in Japanese MN cases (n = 4461) compared with the general population of Japan (n = 20 238). This enrichment of DDX41 risk alleles was much more prominent in male than female (20.7 vs 5.0). P/LP DDX41 variants conferred a large risk of developing MNs, which was negligible until 40 years of age but rapidly increased to 49% by 90 years of age. Patients with myelodysplastic syndromes (MDS) along with a DDX41-mutation rapidly progressed to acute myeloid leukemia (AML), which was however, confined to those having truncating variants. Comutation patterns at diagnosis and at progression to AML were substantially different between DDX41-mutated and WT cases, in which none of the comutations affected clinical outcomes. Even TP53 mutations made no exceptions and their dismal effect, including multihit allelic status, on survival was almost completely mitigated by the presence of DDX41 mutations. Finally, outcomes were not affected by the conventional risk stratifications including the revised/molecular International Prognostic Scoring System. Our findings establish that MDS with DDX41-mutation defines a unique subtype of MNs that is distinct from other MNs.

Age-related remodelling of oesophageal epithelia by mutated cancer drivers.

Clonal expansion in aged normal tissues has been implicated in the development of cancer. However, the chronology and risk dependence of the expansion are poorly understood. Here we intensively sequence 682 micro-scale oesophageal samples and show, in physiologically normal oesophageal epithelia, the progressive age-related expansion of clones that carry mutations in driver genes (predominantly NOTCH1), which is substantially accelerated by alcohol consumption and by smoking. Driver-mutated clones emerge multifocally from early childhood and increase their number and size with ageing, and ultimately replace almost the entire oesophageal epithelium in the extremely elderly. Compared with mutations in oesophageal cancer, there is a marked overrepresentation of NOTCH1 and PPM1D mutations in physiologically normal oesophageal epithelia; these mutations can be acquired before late adolescence (as early as early infancy) and significantly increase in number with heavy smoking and drinking. The remodelling of the oesophageal epithelium by driver-mutated clones is an inevitable consequence of normal ageing, which-depending on lifestyle risks-may affect cancer development.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv.

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

A quantification method of somatic mutations in normal tissues and their accumulation in pediatric patients with chemotherapy.

Somatic mutations are accumulated in normal human tissues with aging and exposure to carcinogens. If we can accurately count any passenger mutations in any single DNA molecule, since their quantity is much larger than driver mutations, we can sensitively detect mutation accumulation in polyclonal normal tissues. Duplex sequencing, which tags both DNA strands in one DNA molecule, enables accurate count of such mutations, but requires a very large number of sequencing reads for each single sample of human-genome size. Here, we reduced the genome size to 1/90 using the BamHI restriction enzyme and established a cost-effective pipeline. The enzymatically cleaved and optimal sequencing (EcoSeq) method was able to count somatic mutations in a single DNA molecule with a sensitivity of as low as 3 × 10-8 per base pair (bp), as assessed by measuring artificially prepared mutations. Taking advantages of EcoSeq, we analyzed normal peripheral blood cells of pediatric sarcoma patients who received chemotherapy (n = 10) and those who did not (n = 10). The former had a mutation frequency of 31.2 ± 13.4 × 10-8 per base pair while the latter had 9.0 ± 4.5 × 10-8 per base pair (P < 0.001). The increase in mutation frequency was confirmed by analysis of the same patients before and after chemotherapy, and increased mutation frequencies persisted 46 to 64 mo after chemotherapy, indicating that the mutation accumulation constitutes a risk of secondary leukemia. EcoSeq has the potential to reveal accumulation of somatic mutations and exposure to environmental factors in any DNA samples and will contribute to cancer risk estimation.

Two novel high-risk adult B-cell acute lymphoblastic leukemia subtypes with high expression of CDX2 and IDH1/2 mutations.

The genetic basis of leukemogenesis in adults with B-cell acute lymphoblastic leukemia (B-ALL) is largely unclear, and its clinical outcome remains unsatisfactory. This study aimed to advance the understanding of biological characteristics, improve disease stratification, and identify molecular targets of adult B-ALL. Adolescents and young adults (AYA) (15 to 39 years old, n = 193) and adults (40 to 64 years old, n = 161) with Philadelphia chromosome-negative (Ph-) B-ALL were included in this study. Integrated transcriptomic and genetic analyses were used to classify the cohort into defined subtypes. Of the 323 cases included in the RNA sequencing analysis, 278 (86.1%) were classified into 18 subtypes. The ZNF384 subtype (22.6%) was the most prevalent, with 2 novel subtypes (CDX2-high and IDH1/2-mut) identified among cases not assigned to the established subtypes. The CDX2-high subtype (3.4%) was characterized by high expression of CDX2 and recurrent gain of chromosome 1q. The IDH1/2-mut subtype (1.9%) was defined by IDH1 R132C or IDH2 R140Q mutations with specific transcriptional and high-methylation profiles. Both subtypes showed poor prognosis and were considered inferior prognostic factors independent of clinical parameters. Comparison with a previously reported pediatric B-ALL cohort (n = 1003) showed that the frequencies of these subtypes were significantly higher in AYA/adults than in children. We delineated the genetic and transcriptomic landscape of adult B-ALL and identified 2 novel subtypes that predict poor disease outcomes. Our findings highlight the age-dependent distribution of subtypes, which partially accounts for the prognostic differences between adult and pediatric B-ALL.

Whole-genome landscape of adult T-cell leukemia/lymphoma.

Adult T-cell leukemia/lymphoma (ATL) is an aggressive neoplasm immunophenotypically resembling regulatory T cells, associated with human T-cell leukemia virus type-1. Here, we performed whole-genome sequencing (WGS) of 150 ATL cases to reveal the overarching landscape of genetic alterations in ATL. We discovered frequent (33%) loss-of-function alterations preferentially targeting the CIC long isoform, which were overlooked by previous exome-centric studies of various cancer types. Long but not short isoform-specific inactivation of Cic selectively increased CD4+CD25+Foxp3+ T cells in vivo. We also found recurrent (13%) 3'-truncations of REL, which induce transcriptional upregulation and generate gain-of-function proteins. More importantly, REL truncations are also common in diffuse large B-cell lymphoma, especially in germinal center B-cell-like subtype (12%). In the non-coding genome, we identified recurrent mutations in regulatory elements, particularly splice sites, of several driver genes. In addition, we characterized the different mutational processes operative in clustered hypermutation sites within and outside immunoglobulin/T-cell receptor genes and identified the mutational enrichment at the binding sites of host and viral transcription factors, suggesting their activities in ATL. By combining the analyses for coding and noncoding mutations, structural variations, and copy number alterations, we discovered 56 recurrently altered driver genes, including 11 novel ones. Finally, ATL cases were classified into 2 molecular groups with distinct clinical and genetic characteristics based on the driver alteration profile. Our findings not only help to improve diagnostic and therapeutic strategies in ATL, but also provide insights into T-cell biology and have implications for genome-wide cancer driver discovery.

UBTF-internal tandem duplication as a novel poor prognostic factor in pediatric acute myeloid leukemia.

The prognosis of pediatric acute myeloid leukemia (AML) has improved via stratification therapy. However, relapse or death occurs in 30%-40% of cases. Novel genetic factors for pediatric AML need to be elucidated to improve prognosis. We detected recurrent internal tandem duplication in upstream binding transcription factor (UBTF-ITD) in 1.2% (6/503) of Japanese pediatric patients with de novo AML. No UBTF-ITD was detected in 175 adult patients with AML or in 65 cell lines that included 15 AML, 39 acute lymphoblastic leukemia, five chronic myeloid leukemia, and six neuroblastoma cell lines. All UBTF-ITDs were found in exon 13 and shared a duplicated region. UBTF-ITD was more frequently detected in patients with trisomy 8, FLT3-ITD, WT1 mutation, and/or high PRDM16 expression (trisomy 8, 3/6; FLT3-ITD, 5/6; WT1 mutation, 2/6; and high PRDM16 expression, 6/6). Gene expression patterns of patients with UBTF-ITD were similar to those of patients with NUP98::NSD1 or FUS::ERG. Survival analysis of the AML-05 cohort revealed that patients with UBTF-ITD had worse outcomes than those without UBTF-ITD (3-year event-free survival, 20% vs. 55%; 3-year overall survival, 40% vs. 74%). Moreover, among the 27 patients with trisomy 8, all three patients with UBTF -ITD had a poor prognosis resulting in early events (relapse or non-complete remission) within 1 year. Our findings suggest that UBTF-ITD may be a novel and significant prognostic factor for pediatric patients with AML.

TP53 and RB1 alterations characterize poor prognostic subgroups in pediatric acute myeloid leukemia.

Pediatric acute myeloid leukemia (AML) is a poor prognostic subtype of pediatric leukemia. However, the detailed characteristics of many genetic abnormalities are yet to be established in this disease. Although TP53 and RB1 are established as representative tumor suppressor genes in various cancers, alterations of these two genes, especially RB1, have not been characterized in pediatric AML. We performed next-generation sequencing in 328 pediatric AML patients from the Japanese AML-05 trial to ascertain TP53 and RB1 alterations, and their prognostic implications. We identified seven patients with TP53 alterations (2.1%) and six patients with RB1 alterations (1.8%). These alterations were found in only patients without RUNX1::RUNX1T1, CBFB::MYH11, or KMT2A rearrangements. TP53 and RB1 were frequently co-deleted with their neighboring genes PRPF8 and ELF1, respectively. Patients with TP53 alterations had significantly lower 5-year overall survival (OS; 14.3% vs. 71.4%, p < 0.001) and lower 5-year event-free survival (EFS; 0% vs. 56.3%, p < 0.001); similarly, patients with RB1 had significantly lower 5-year OS (0% vs. 71.8%, p < 0.001) and lower 5-year EFS (0% vs. 56.0%, p < 0.001) when compared to patients without these alterations. In gene expression analyses, oxidative phosphorylation, glycolysis, and protein secretion were upregulated in patients with TP53 and/or RB1 alterations. Additionally, Kaplan-Meier analysis revealed that high expressions of SLC2A5, KCNAB2, and CD300LF were related to poor OS of non-core-binding factor AML patients (p < 0.001, p = 0.001, and p = 0.021, respectively). This study will contribute to the development of risk-stratified therapy and precision medicine in pediatric AML.

The Impact of FGFR3 Alterations on the Tumor Microenvironment and the Efficacy of Immune Checkpoint Inhibitors in Bladder Cancer.

BACKGROUND: Currently, only limited knowledge is available regarding the phenotypic association between fibroblast growth factor receptor 3 (FGFR3) alterations and the tumor microenvironment (TME) in bladder cancer (BLCA). METHODS: A multi-omics analysis on 389 BLCA and 35 adjacent normal tissues from a cohort of OMPU-NCC Consortium Japan was retrospectively performed by integrating the whole-exome and RNA-sequence dataset and clinicopathological record. A median follow-up duration of all BLCA cohort was 31 months. RESULTS: FGFR3 alterations (aFGFR3), including recurrent mutations and fusions, accounted for 44% of non-muscle invasive bladder cancer (NMIBC) and 15% of muscle-invasive bladder cancer (MIBC). Within MIBC, the consensus subtypes LumP was significantly more prevalent in aFGFR3, whereas the Ba/Sq subtype exhibited similarity between intact FGFR3 (iFGFR3) and aFGFR3 cases. We revealed that basal markers were significantly increased in MIBC/aFGFR3 compared to MIBC/iFGFR3. Transcriptome analysis highlighted TIM3 as the most upregulated immune-related gene in iFGFR3, with differential immune cell compositions observed between iFGFR3 and aFGFR3. Using EcoTyper, TME heterogeneity was discerned even within aFGFR cases, suggesting potential variations in the response to checkpoint inhibitors (CPIs). Among 72 patients treated with CPIs, the objective response rate (ORR) was comparable between iFGFR3 and aFGFR3 (20% vs 31%; p = 0.467). Strikingly, a significantly higher ORR was noted in LumP/aFGFR3 compared to LumP/iFGFR3 (50% vs 5%; p = 0.022). This trend was validated using data from the IMvigor210 trial. Additionally, several immune-related genes, including IDO1, CCL24, IL1RL1, LGALS4, and NCAM (CD56) were upregulated in LumP/iFGFR3 compared to LumP/aFGFR3 cases. CONCLUSIONS: Differential pathways influenced by aFGFR3 were observed between NMIBC and MIBC, highlighting the upregulation of both luminal and basal markers in MIBC/aFGFR3. Heterogeneous TME was identified within MIBC/aFGFR3, leading to differential outcomes for CPIs. Specifically, a favorable ORR in LumP/aFGFR3 and a poor ORR in LumP/iFGFR3 were observed. We propose TIM3 as a potential target for iFGFR3 (ORR: 20%) and several immune checkpoint genes, including IDO1 and CCL24, for LumP/iFGFR3 (ORR: 5%), indicating promising avenues for precision immunotherapy for BLCA.

Tertiary lymphoid structure and neutrophil-lymphocyte ratio coordinately predict outcome of pembrolizumab.

Emerging evidence suggests that the presence of tertiary lymphoid structures (TLS) and neutrophil-lymphocyte ratio (NLR) in peripheral blood is associated with the treatment response to checkpoint inhibitors (CPIs), whereas there is limited knowledge regarding whether these factors reciprocally impact the treatment outcomes of CPIs in metastatic urothelial carcinoma (mUC). Herein, we investigated treatment outcomes of platinum-refractory mUC patients (50 cases with whole-exome and transcriptome sequencing) treated with pembrolizumab. The pathological review identified 24% of cases of TLS in the specimens. There was no significant difference in the NLR between the TLS- and TLS+ groups (p = 0.153). In the lower NLR group, both overall survival and progression-free survival were significantly longer in patients with TLS than in those without TLS, whereas the favorable outcomes associated with TLS were not observed in patients in the higher NLR group. We explored transcriptomic differences in UC with TLS. The TLS was comparably observed between luminal (20%) and basal (25%) tumor subtypes (p = 0.736). Exploring putative immune-checkpoint genes revealed that ICOSLG (B7-H2) was significantly increased in tumors with lower NLR. KRT expression levels exhibited higher basal cell markers (KRT5 and KRT17) in the higher NLR group and lower differentiated cell markers (KRT8 and KRT18) in patients with TLS. In conclusion, the improved outcomes of pembrolizumab treatment in mUC are restricted to patients with lower NLR. Our findings begin to elucidate a distinct molecular pattern for the presence of TLS according to the NLR in peripheral blood.

Comprehensive analysis of indels in whole-genome microsatellite regions and microsatellite instability across 21 cancer types.

Microsatellites are repeats of 1- to 6-bp units, and approximately 10 million microsatellites have been identified across the human genome. Microsatellites are vulnerable to DNA mismatch errors and have thus been used to detect cancers with mismatch repair deficiency. To reveal the mutational landscape of microsatellite repeat regions at the genome level, we analyzed approximately 20.1 billion microsatellites in 2717 whole genomes of pan-cancer samples across 21 tissue types. First, we developed a new insertion and deletion caller (MIMcall) that takes into consideration the error patterns of different types of microsatellites. Among the 2717 pan-cancer samples, our analysis identified 31 samples, including colorectal, uterus, and stomach cancers, with a higher proportion of mutated microsatellite (≥0.03), which we defined as microsatellite instability (MSI) cancers of genome-wide level. Next, we found 20 highly mutated microsatellites that can be used to detect MSI cancers with high sensitivity. Third, we found that replication timing and DNA shape were significantly associated with mutation rates of microsatellites. Last, analysis of mutations in mismatch repair genes showed that somatic SNVs and short indels had larger functional impacts than germline mutations and structural variations. Our analysis provides a comprehensive picture of mutations in the microsatellite regions and reveals possible causes of mutations, as well as provides a useful marker set for MSI detection.

Long-read sequencing for non-small-cell lung cancer genomes.

Here, we report the application of a long-read sequencer, PromethION, for analyzing human cancer genomes. We first conducted whole-genome sequencing on lung cancer cell lines. We found that it is possible to genotype known cancerous mutations, such as point mutations. We also found that long-read sequencing is particularly useful for precisely identifying and characterizing structural aberrations, such as large deletions, gene fusions, and other chromosomal rearrangements. In addition, we identified several medium-sized structural aberrations consisting of complex combinations of local duplications, inversions, and microdeletions. These complex mutations occurred even in key cancer-related genes, such as STK11, NF1, SMARCA4, and PTEN The biological relevance of those mutations was further revealed by epigenome, transcriptome, and protein analyses of the affected signaling pathways. Such structural aberrations were also found in clinical lung adenocarcinoma specimens. Those structural aberrations were unlikely to be reliably detected by conventional short-read sequencing. Therefore, long-read sequencing may contribute to understanding the molecular etiology of patients for whom causative cancerous mutations remain unknown and therapeutic strategies are elusive.

Frequent genetic alterations in immune checkpoint-related genes in intravascular large B-cell lymphoma.

Intravascular large B-cell lymphoma (IVLBCL) is a unique type of extranodal lymphoma characterized by selective growth of tumor cells in small vessels without lymphadenopathy. Greater understanding of the molecular pathogenesis of IVLBCL is hampered by the paucity of lymphoma cells in biopsy specimens, creating a limitation in obtaining sufficient tumor materials. To uncover the genetic landscape of IVLBCL, we performed whole-exome sequencing (WES) of 21 patients with IVLBCL using plasma-derived cell-free DNA (cfDNA) (n = 18), patient-derived xenograft tumors (n = 4), and tumor DNA from bone marrow (BM) mononuclear cells (n = 2). The concentration of cfDNA in IVLBCL was significantly higher than that in diffuse large B-cell lymphoma (DLBCL) (P < .0001) and healthy donors (P = .0053), allowing us to perform WES; most mutations detected in BM tumor DNA were successfully captured in cfDNA and xenograft. IVLBCL showed a high frequency of genetic lesions characteristic of activated B-cell-type DLBCL, with the former showing conspicuously higher frequencies (compared with nodal DLBCL) of mutations in MYD88 (57%), CD79B (67%), SETD1B (57%), and HLA-B (57%). We also found that 8 IVLBCL (38%) harbored rearrangements of programmed cell death 1 ligand 1 and 2 (PD-L1/PD-L2) involving the 3' untranslated region; such rearrangements are implicated in immune evasion via PD-L1/PD-L2 overexpression. Our data demonstrate the utility of cfDNA and imply important roles for immune evasion in IVLBCL pathogenesis and PD-1/PD-L1/PD-L2 blockade in therapeutics for IVLBCL.