Comprehensive genetic profiling reveals frequent alterations of driver genes on the X chromosome in extranodal NK/T-cell lymphoma.

Extranodal NK/T-cell lymphoma (ENKTCL) is an Epstein-Barr virus (EBV)-related neoplasm with male dominance and a poor prognosis. A better understanding of the genetic alterations and their functional roles in ENKTCL could help improve patient stratification and treatments. Here, we performed comprehensive genetic analysis of 177 ENKTCL cases to delineate the landscape of mutations, copy number alterations (CNAs), and structural variations, identifying 34 driver genes including six previously unappreciated ones, namely HLA-B, HLA-C, ROBO1, CD58, POT1, and MAP2K1. Among them, CD274 (24%) was the most frequently altered, followed by TP53 (20%), CDKN2A (19%), ARID1A (15%), HLA-A (15%), BCOR (14%), and MSN (14%). Chromosome X (chrX) losses were the most common arm-level CNAs in females (~40%), and alterations of four X-linked driver genes (MSN, BCOR, DDX3X, and KDM6A) were more frequent in males and females harboring chrX losses. Among X-linked drivers, MSN was the most recurrently altered, and its expression was lost in approximately one-third of cases using immunohistochemical analysis. Functional studies of human cell lines demonstrated that MSN disruption promoted cell proliferation and NF-κB activation. Moreover, MSN inactivation increased sensitivity to NF-κB inhibition in vitro and in vivo. In addition, recurrent deletions were observed at the origin of replication in the EBV genome (6%). Finally, by integrating the 34 drivers and 19 significant arm-level CNAs, non-negative matrix factorization and consensus clustering identified two molecular groups with different genetic features and prognosis irrespective of clinical prognostic factors. Together, these findings could help improve diagnostic and therapeutic strategies in ENKTCL.

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

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 AML-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 < 0.001; 25.0% vs. 89.5%, P < 0.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.

Assessing the efficacy of target adaptive sampling long-read sequencing through hereditary cancer patient genomes.

Innovations in sequencing technology have led to the discovery of novel mutations that cause inherited diseases. However, many patients with suspected genetic diseases remain undiagnosed. Long-read sequencing technologies are expected to significantly improve the diagnostic rate by overcoming the limitations of short-read sequencing. In addition, Oxford Nanopore Technologies (ONT) offers adaptive sampling and computationally driven target enrichment technology. This enables more affordable intensive analysis of target gene regions compared to standard non-selective long-read sequencing. In this study, we developed an efficient computational workflow for target adaptive sampling long-read sequencing (TAS-LRS) and evaluated it through application to 33 genomes collected from suspected hereditary cancer patients. Our workflow can identify single nucleotide variants with nearly the same accuracy as the short-read platform and elucidate complex forms of structural variations. We also newly identified several SINE-R/VNTR/Alu (SVA) elements affecting the APC gene in two patients with familial adenomatous polyposis, as well as their sites of origin. In addition, we demonstrated that off-target reads from adaptive sampling, which is typically discarded, can be effectively used to accurately genotype common single-nucleotide polymorphisms (SNPs) across the entire genome, enabling the calculation of a polygenic risk score. Furthermore, we identified allele-specific MLH1 promoter hypermethylation in a Lynch syndrome patient. In summary, our workflow with TAS-LRS can simultaneously capture monogenic risk variants including complex structural variations, polygenic background as well as epigenetic alterations, and will be an efficient platform for genetic disease research and diagnosis.

Inter- and intra-tumor heterogeneity of genetic and immune profiles in inherited renal cell carcinoma.

Patients with von Hippel-Lindau disease (vHL) are at risk of developing spatially and temporally multiple clear cell renal cell carcinomas (ccRCCs), which offers a valuable opportunity to analyze inter- and intra-tumor heterogeneity of genetic and immune profiles within the same patient. Here, we perform whole-exome and RNA sequencing, digital gene expression, and immunohistochemical analyses for 81 samples from 51 ccRCCs of 10 patients with vHL. Inherited ccRCCs are clonally independent and have less genomic alterations than sporadic ccRCCs. Hierarchical clustering of transcriptome profiles shows two clusters with distinct immune signatures: immune hot and cold clusters. Interestingly, not only samples from the same tumors but also different tumors from the same patients tend to show a similar immune signature, whereas samples from different patients frequently exhibit different signatures. Our findings reveal the genetic and immune landscape of inherited ccRCCs, demonstrating the relevance of host factors in shaping anti-tumor immunity.

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.

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.

Multiancestry genomic and transcriptomic analysis of gastric cancer.

Gastric cancer is among the most common malignancies worldwide, characterized by geographical, epidemiological and histological heterogeneity. Here, we report an extensive, multiancestral landscape of driver events in gastric cancer, involving 1,335 cases. Seventy-seven significantly mutated genes (SMGs) were identified, including ARHGAP5 and TRIM49C. We also identified subtype-specific drivers, including PIGR and SOX9, which were enriched in the diffuse subtype of the disease. SMGs also varied according to Epstein-Barr virus infection status and ancestry. Non-protein-truncating CDH1 mutations, which are characterized by in-frame splicing alterations, targeted localized extracellular domains and uniquely occurred in sporadic diffuse-type cases. In patients with gastric cancer with East Asian ancestry, our data suggested a link between alcohol consumption or metabolism and the development of RHOA mutations. Moreover, mutations with potential roles in immune evasion were identified. Overall, these data provide comprehensive insights into the molecular landscape of gastric cancer across various subtypes and ancestries.

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.

Systematic identification of intron retention associated variants from massive publicly available transcriptome sequencing data.

Many disease-associated genomic variants disrupt gene function through abnormal splicing. With the advancement of genomic medicine, identifying disease-associated splicing associated variants has become more important than ever. Most bioinformatics approaches to detect splicing associated variants require both genome and transcriptomic data. However, there are not many datasets where both of them are available. In this study, we develop a methodology to detect genomic variants that cause splicing changes (more specifically, intron retention), using transcriptome sequencing data alone. After evaluating its sensitivity and precision, we apply it to 230,988 transcriptome sequencing data from the publicly available repository and identified 27,049 intron retention associated variants (IRAVs). In addition, by exploring positional relationships with variants registered in existing disease databases, we extract 3,000 putative disease-associated IRAVs, which range from cancer drivers to variants linked with autosomal recessive disorders. The in-silico screening framework demonstrates the possibility of near-automatically acquiring medical knowledge, making the most of massively accumulated publicly available sequencing data. Collections of IRAVs identified in this study are available through IRAVDB ( https://iravdb.io/ ).

Multi-omics analysis defines highly refractory RAS burdened immature subgroup of infant acute lymphoblastic leukemia.

KMT2A-rearranged infant acute lymphoblastic leukemia (ALL) represents the most refractory type of childhood leukemia. To uncover the molecular heterogeneity of this disease, we perform RNA sequencing, methylation array analysis, whole exome and targeted deep sequencing on 84 infants with KMT2A-rearranged leukemia. Our multi-omics clustering followed by single-sample and single-cell inference of hematopoietic differentiation establishes five robust integrative clusters (ICs) with different master transcription factors, fusion partners and corresponding stages of B-lymphopoietic and early hemato-endothelial development: IRX-type differentiated (IC1), IRX-type undifferentiated (IC2), HOXA-type MLLT1 (IC3), HOXA-type MLLT3 (IC4), and HOXA-type AFF1 (IC5). Importantly, our deep mutational analysis reveals that the number of RAS pathway mutations predicts prognosis and that the most refractory subgroup of IC2 possesses 100% frequency and the heaviest burden of RAS pathway mutations. Our findings highlight the previously under-appreciated intra- and inter-patient heterogeneity of KMT2A-rearranged infant ALL and provide a rationale for the future development of genomics-guided risk stratification and individualized therapy.

Amplified EPOR/JAK2 Genes Define a Unique Subtype of Acute Erythroid Leukemia.

Acute erythroid leukemia (AEL) is a unique subtype of acute myeloid leukemia characterized by prominent erythroid proliferation whose molecular basis is poorly understood. To elucidate the underlying mechanism of erythroid proliferation, we analyzed 121 AEL using whole-genome, whole-exome, and/or targeted-capture sequencing, together with transcriptome analysis of 21 AEL samples. Combining publicly available sequencing data, we found a high frequency of gains and amplifications involving EPOR/JAK2 in TP53-mutated cases, particularly those having >80% erythroblasts designated as pure erythroid leukemia (10/13). These cases were frequently accompanied by gains and amplifications of ERG/ETS2 and associated with a very poor prognosis, even compared with other TP53-mutated AEL. In addition to activation of the STAT5 pathway, a common feature across all AEL cases, these AEL cases exhibited enhanced cell proliferation and heme metabolism and often showed high sensitivity to ruxolitinib in vitro and in xenograft models, highlighting a potential role of JAK2 inhibition in therapeutics of AEL. SIGNIFICANCE: This study reveals the major role of gains, amplifications, and mutations of EPOR and JAK2 in the pathogenesis of pure erythroleukemia. Their frequent response to ruxolitinib in patient-derived xenograft and cell culture models highlights a possible therapeutic role of JAK2 inhibition for erythroleukemia with EPOR/JAK2-involving lesions. This article is highlighted in the In This Issue feature, p. 369.

Genetic Analysis of Pheochromocytoma and Paraganglioma Complicating Cyanotic Congenital Heart Disease.

CONTEXT: Pheochromocytoma and paraganglioma (PPGL) may appear as a complication of cyanotic congenital heart disease (CCHD-PPGL) with frequent EPAS1 mutations, suggesting a close link between EPAS1 mutations and tissue hypoxia in CCHD-PPGL pathogenesis. OBJECTIVE: Our aim is to further investigate the role of EPAS1 mutations in the hypoxia-driven mechanism of CCHD-PPGL pathogenesis, particularly focusing on metachronous and/or multifocal CCHD-PPGL tumors. METHODS: We performed whole-exome sequencing (WES) for somatic and germline mutations in 15 PPGL samples from 7 CCHD patients, including 3 patients with metachronous and/or multifocal tumors, together with an adrenal medullary hyperplasia (AMH) sample. RESULTS: We detected EPAS1 mutations in 15 out of 16 PPGL/AMH samples from 7 cases. Conspicuously, all EPAS1 mutations in each of 3 cases with multifocal or metachronous tumors were mutually independent and typical examples of parallel evolution, which is suggestive of strong positive selection of EPAS1-mutated clones. Compared to 165 The Cancer Genome Atlas non-CCHD-PPGL samples, CCHD-PPGL/AMH samples were enriched for 11p deletions (13/16) and 2p amplifications (4/16). Of particular note, the multiple metachronous PPGL tumors with additional copy number abnormalities developed 18 to 23 years after the resolution of hypoxemia, suggesting that CCHD-induced hypoxic environments are critical for positive selection of EPAS1 mutants in early life, but may no longer be required for development of PPGL in later life. CONCLUSION: Our results highlight a key role of activated hypoxia-inducible factor 2α due to mutated EPAS1 in positive selection under hypoxic environments, although hypoxemia itself may not necessarily be required for the EPAS1-mutated clones to progress to PPGL.

The landscape of genetic aberrations in myxofibrosarcoma.

Myxofibrosarcoma (MFS) is a rare subtype of sarcoma, whose genetic basis is poorly understood. We analyzed 69 MFS cases using whole-genome (WGS), whole-exome (WES) and/or targeted-sequencing (TS). Newly sequenced genomic data were combined with additional deposited 116 MFS samples. WGS identified a high number of structural variations (SVs) per tumor most frequently affecting the TP53 and RB1 loci, 40% of tumors showed a BRCAness-associated mutation signature, and evidence of chromothripsis was found in all cases. Most frequently mutated/copy number altered genes affected known disease drivers such as TP53 (56.2%), CDKN2A/B (29.7%), RB1 (27.0%), ATRX (19.5%) and HDLBP (18.9%). Several previously unappreciated genetic aberrations including MUC17, FLG and ZNF780A were identified in more than 20% of patients. Longitudinal analysis of paired diagnosis and relapse time points revealed a 1.2-fold mutation number increase accompanied with substantial changes in clonal composition over time. Our study highlights the genetic complexity underlying sarcomagenesis of MFS.

Alteration of the immune environment in bone marrow from children with recurrent B cell precursor acute lymphoblastic leukemia.

Due to the considerable success of cancer immunotherapy for leukemia, the tumor immune environment has become a focus of intense research; however, there are few reports on the dynamics of the tumor immune environment in leukemia. Here, we analyzed the tumor immune environment in pediatric B cell precursor acute lymphoblastic leukemia by analyzing serial bone marrow samples from nine patients with primary and recurrent disease by mass cytometry using 39 immunophenotype markers, and transcriptome analysis. High-dimensional single-cell mass cytometry analysis elucidated a dynamic shift of T cells from naïve to effector subsets, and clarified that, during relapse, the tumor immune environment comprised a T helper 1-polarized immune profile, together with an increased number of effector regulatory T cells. These results were confirmed in a validation cohort using conventional flow cytometry. Furthermore, RNA transcriptome analysis identified the upregulation of immune-related pathways in B cell precursor acute lymphoblastic leukemia cells during relapse, suggesting interaction with the surrounding environment. In conclusion, a tumor immune environment characterized by a T helper 1-polarized immune profile, with an increased number of effector regulatory T cells, could contribute to the pathophysiology of recurrent B cell precursor acute lymphoblastic leukemia. This information could contribute to the development of effective immunotherapeutic approaches against B cell precursor acute lymphoblastic leukemia relapse.

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.

The Evolving Genomic Landscape of Esophageal Squamous Cell Carcinoma Under Chemoradiotherapy.

Esophageal squamous cell carcinoma (ESCC) often recurs after chemoradiotherapy, and the prognosis of ESCC after chemoradiotherapy has not improved over the past few decades. The mutation process in chemoradiotherapy-resistant clones and the functional relevance of genetic alterations remain unclear. To address these problems, we performed whole-exome sequencing of 52 tumor samples from 33 patients with ESCC who received radiotherapy combined with 5-fluorouracil/platinum. In multiregion analyses of pretreatment and locally recurrent lesions from five cases, most driver gene-altered clones remained under chemoradiotherapy selection pressure, while few driver gene alterations were acquired at recurrence. The mutation signatures of recurrent ESCC, including increased deletion frequency and platinum dose-dependent base substitution signatures, were substantially different from those of primary ESCC and reflected the iatrogenic impacts of chemoradiotherapy. Single-region analysis of 28 pretreatment tumors indicated that focal copy-number gain at the MYC locus was significantly associated with poor progression-free survival and overall survival after chemoradiotherapy. MYC gain remained throughout the chemoradiotherapy course and potentially contributes to intrinsic resistance to chemoradiotherapy. Consistent with these findings, MYC copy number and mRNA and protein levels in ESCC cell lines correlated positively with resistance to radiotherapy, and MYC knockdown improved sensitivity to radiotherapy. Overall, these data characterize the clonal evolution process induced by chemoradiotherapy and clinically relevant associations for genetic alterations in ESCC. These findings increase our understanding of therapeutic resistance and support the rationale for precision chemoradiotherapy. SIGNIFICANCE: Whole-exome sequencing reveals the genetic evolution of ESCC during chemoradiotherapy, highlighting MYC gain in pretreatment tumors as a potential marker of therapy resistance.

Combined landscape of single-nucleotide variants and copy number alterations in clonal hematopoiesis.

Clonal hematopoiesis (CH) in apparently healthy individuals is implicated in the development of hematological malignancies (HM) and cardiovascular diseases. Previous studies of CH analyzed either single-nucleotide variants and indels (SNVs/indels) or copy number alterations (CNAs), but not both. Here, using a combination of targeted sequencing of 23 CH-related genes and array-based CNA detection of blood-derived DNA, we have delineated the landscape of CH-related SNVs/indels and CNAs in 11,234 individuals without HM from the BioBank Japan cohort, including 672 individuals with subsequent HM development, and studied the effects of these somatic alterations on mortality from HM and cardiovascular disease, as well as on hematological and cardiovascular phenotypes. The total number of both types of CH-related lesions and their clone size positively correlated with blood count abnormalities and mortality from HM. CH-related SNVs/indels and CNAs exhibited statistically significant co-occurrence in the same individuals. In particular, co-occurrence of SNVs/indels and CNAs affecting DNMT3A, TET2, JAK2 and TP53 resulted in biallelic alterations of these genes and was associated with higher HM mortality. Co-occurrence of SNVs/indels and CNAs also modulated risks for cardiovascular mortality. These findings highlight the importance of detecting both SNVs/indels and CNAs in the evaluation of CH.

Molecular classification and diagnostics of upper urinary tract urothelial carcinoma.

Upper urinary tract urothelial carcinoma (UTUC) is one of the common urothelial cancers. Its molecular pathogenesis, however, is poorly understood, with no useful biomarkers available for accurate diagnosis and molecular classification. Through an integrated genetic study involving 199 UTUC samples, we delineate the landscape of genetic alterations in UTUC enabling genetic/molecular classification. According to the mutational status of TP53, MDM2, RAS, and FGFR3, UTUC is classified into five subtypes having discrete profiles of gene expression, tumor location/histology, and clinical outcome, which is largely recapitulated in an independent UTUC cohort. Sequencing of urine sediment-derived DNA has a high diagnostic value for UTUC with 82.2% sensitivity and 100% specificity. These results provide a solid basis for better diagnosis and management of UTUC.

Clonal evolution and clinical implications of genetic abnormalities in blastic transformation of chronic myeloid leukaemia.

Blast crisis (BC) predicts dismal outcomes in patients with chronic myeloid leukaemia (CML). Although additional genetic alterations play a central role in BC, the landscape and prognostic impact of these alterations remain elusive. Here, we comprehensively investigate genetic abnormalities in 136 BC and 148 chronic phase (CP) samples obtained from 216 CML patients using exome and targeted sequencing. One or more genetic abnormalities are found in 126 (92.6%) out of the 136 BC patients, including the RUNX1-ETS2 fusion and NBEAL2 mutations. The number of genetic alterations increase during the transition from CP to BC, which is markedly suppressed by tyrosine kinase inhibitors (TKIs). The lineage of the BC and prior use of TKIs correlate with distinct molecular profiles. Notably, genetic alterations, rather than clinical variables, contribute to a better prediction of BC prognosis. In conclusion, genetic abnormalities can help predict clinical outcomes and can guide clinical decisions in CML.