The genomic basis of childhood T-lineage acute lymphoblastic leukaemia.

T-lineage acute lymphoblastic leukaemia (T-ALL) is a high-risk tumour1 that has eluded comprehensive genomic characterization, which is partly due to the high frequency of noncoding genomic alterations that result in oncogene deregulation2,3. Here we report an integrated analysis of genome and transcriptome sequencing of tumour and remission samples from more than 1,300 uniformly treated children with T-ALL, coupled with epigenomic and single-cell analyses of malignant and normal T cell precursors. This approach identified 15 subtypes with distinct genomic drivers, gene expression patterns, developmental states and outcomes. Analyses of chromatin topology revealed multiple mechanisms of enhancer deregulation that involve enhancers and genes in a subtype-specific manner, thereby demonstrating widespread involvement of the noncoding genome. We show that the immunophenotypically described, high-risk entity of early T cell precursor ALL is superseded by a broader category of 'early T cell precursor-like' leukaemia. This category has a variable immunophenotype and diverse genomic alterations of a core set of genes that encode regulators of hematopoietic stem cell development. Using multivariable outcome models, we show that genetic subtypes, driver and concomitant genetic alterations independently predict treatment failure and survival. These findings provide a roadmap for the classification, risk stratification and mechanistic understanding of this disease.

G-quadruplexes are a source of vulnerability in BRCA2 deficient granule cell progenitors and medulloblastoma.

Biallelic pathogenic variants in the essential DNA repair gene BRCA2 causes Fanconi anemia, complementation group FA-D1. Patients in this group are highly prone to develop embryonal tumors, most commonly medulloblastoma arising from the cerebellar granule cell progenitors (GCPs). GCPs undergo high proliferation in the postnatal cerebellum under SHH activation, but the type of DNA lesions that require the function of the BRCA2 to prevent tumorigenesis remains unknown. To identify such lesions, we assessed both GCP neurodevelopment and tumor formation using a mouse model with deletion of exons three and four of Brca2 in the central nervous system, coupled with global Trp53 loss. Brca2 Δex3-4 ;Trp53 -/- animals developed SHH subgroup medulloblastomas with complete penetrance. Whole-genome sequencing of the tumors identified structural variants with breakpoints enriched in areas overlapping G-quadruplexes (G4s). Brca2-deficient GCPs exhibited decreased replication speed in the presence of the G4-stabilizer pyridostatin. Pif1 helicase, which resolves G4s during replication, was highly upregulated in tumors, and Pif1 knockout in primary MB tumor cells resulted in increased genome instability upon pyridostatin treatment. These data suggest that G4s may represent sites prone to replication stalling in highly proliferative GCPs and without BRCA2, G4s become a source of genome instability. Tumor cells upregulate G4-resolving helicases to facilitate rapid proliferation through G4s highlighting PIF1 helicase as a potential therapeutic target for treatment of BRCA2-deficient medulloblastomas.

The Fanconi anemia pathway induces chromothripsis and ecDNA-driven cancer drug resistance.

Chromothripsis describes the catastrophic shattering of mis-segregated chromosomes trapped within micronuclei. Although micronuclei accumulate DNA double-strand breaks and replication defects throughout interphase, how chromosomes undergo shattering remains unresolved. Using CRISPR-Cas9 screens, we identify a non-canonical role of the Fanconi anemia (FA) pathway as a driver of chromothripsis. Inactivation of the FA pathway suppresses chromosome shattering during mitosis without impacting interphase-associated defects within micronuclei. Mono-ubiquitination of FANCI-FANCD2 by the FA core complex promotes its mitotic engagement with under-replicated micronuclear chromosomes. The structure-selective SLX4-XPF-ERCC1 endonuclease subsequently induces large-scale nucleolytic cleavage of persistent DNA replication intermediates, which stimulates POLD3-dependent mitotic DNA synthesis to prime shattered fragments for reassembly in the ensuing cell cycle. Notably, FA-pathway-induced chromothripsis generates complex genomic rearrangements and extrachromosomal DNA that confer acquired resistance to anti-cancer therapies. Our findings demonstrate how pathological activation of a central DNA repair mechanism paradoxically triggers cancer genome evolution through chromothripsis.

Widespread chromatin context-dependencies of DNA double-strand break repair proteins.

DNA double-strand breaks are repaired by multiple pathways, including non-homologous end-joining (NHEJ) and microhomology-mediated end-joining (MMEJ). The balance of these pathways is dependent on the local chromatin context, but the underlying mechanisms are poorly understood. By combining knockout screening with a dual MMEJ:NHEJ reporter inserted in 19 different chromatin environments, we identified dozens of DNA repair proteins that modulate pathway balance dependent on the local chromatin state. Proteins that favor NHEJ mostly synergize with euchromatin, while proteins that favor MMEJ generally synergize with distinct types of heterochromatin. Examples of the former are BRCA2 and POLL, and of the latter the FANC complex and ATM. Moreover, in a diversity of human cancer types, loss of several of these proteins alters the distribution of pathway-specific mutations between heterochromatin and euchromatin. Together, these results uncover a complex network of proteins that regulate MMEJ:NHEJ balance in a chromatin context-dependent manner.

An IL-1ß-driven neutrophil-stromal cell axis fosters a BAFF-rich protumor microenvironment in individuals with multiple myeloma.

Human bone marrow permanently harbors high numbers of neutrophils, and a tumor-supportive bias of these cells could significantly impact bone marrow-confined malignancies. In individuals with multiple myeloma, the bone marrow is characterized by inflammatory stromal cells with the potential to influence neutrophils. We investigated myeloma-associated alterations in human marrow neutrophils and the impact of stromal inflammation on neutrophil function. Mature neutrophils in myeloma marrow are activated and tumor supportive and transcribe increased levels of IL1B and myeloma cell survival factor TNFSF13B (BAFF). Interactions with inflammatory stromal cells induce neutrophil activation, including BAFF secretion, in a STAT3-dependent manner, and once activated, neutrophils gain the ability to reciprocally induce stromal activation. After first-line myeloid-depleting antimyeloma treatment, human bone marrow retains residual stromal inflammation, and newly formed neutrophils are reactivated. Combined, we identify a neutrophil-stromal cell feed-forward loop driving tumor-supportive inflammation that persists after treatment and warrants novel strategies to target both stromal and immune microenvironments in multiple myeloma.

A Single-Cell Taxonomy Predicts Inflammatory Niche Remodeling to Drive Tissue Failure and Outcome in Human AML.

Cancer initiation is orchestrated by an interplay between tumor-initiating cells and their stromal/immune environment. Here, by adapted single-cell RNA sequencing, we decipher the predicted signaling between tissue-resident hematopoietic stem/progenitor cells (HSPC) and their neoplastic counterparts with their native niches in the human bone marrow. LEPR+ stromal cells are identified as central regulators of hematopoiesis through predicted interactions with all cells in the marrow. Inflammatory niche remodeling and the resulting deprivation of critical HSPC regulatory factors are predicted to repress high-output hematopoietic stem cell subsets in NPM1-mutated acute myeloid leukemia (AML), with relative resistance of clonal cells. Stromal gene signatures reflective of niche remodeling are associated with reduced relapse rates and favorable outcomes after chemotherapy across all genetic risk categories. Elucidation of the intercellular signaling defining human AML, thus, predicts that inflammatory remodeling of stem cell niches drives tissue repression and clonal selection but may pose a vulnerability for relapse-initiating cells in the context of chemotherapeutic treatment. SIGNIFICANCE: Tumor-promoting inflammation is considered an enabling characteristic of tumorigenesis, but mechanisms remain incompletely understood. By deciphering the predicted signaling between tissue-resident stem cells and their neoplastic counterparts with their environment, we identify inflammatory remodeling of stromal niches as a determinant of normal tissue repression and clinical outcomes in human AML. See related commentary by Lisi-Vega and Méndez-Ferrer, p. 349. This article is featured in Selected Articles from This Issue, p. 337.

Hematopoietic Cell Autonomous Disruption of Hematopoiesis in a Germline Loss-of-function Mouse Model of RUNX1-FPD.

RUNX1 familial platelet disorder (RUNX1-FPD) is a hematopoietic disorder caused by germline loss-of-function mutations in the RUNX1 gene and characterized by thrombocytopathy, thrombocytopenia, and an increased risk of developing hematologic malignancies, mostly of myeloid origin. Disease pathophysiology has remained incompletely understood, in part because of a shortage of in vivo models recapitulating the germline RUNX1 loss of function found in humans, precluding the study of potential contributions of non-hematopoietic cells to disease pathogenesis. Here, we studied mice harboring a germline hypomorphic mutation of one Runx1 allele with a loss-of-function mutation in the other Runx1 allele (Runx1 L148A/- mice), which display many hematologic characteristics found in human RUNX1-FPD patients. Runx1 L148A/- mice displayed robust and pronounced thrombocytopenia and myeloid-biased hematopoiesis, associated with an HSC intrinsic reconstitution defect in lymphopoiesis and expansion of myeloid progenitor cell pools. We demonstrate that specific deletion of Runx1 from bone marrow stromal cells in Prrx1-cre;Runx1 fl/fl mice did not recapitulate these abnormalities, indicating that the hematopoietic abnormalities are intrinsic to the hematopoietic lineage, and arguing against a driving role of the bone marrow microenvironment. In conclusion, we report a RUNX1-FPD mouse model faithfully recapitulating key characteristics of human disease. Findings do not support a driving role of ancillary, non-hematopoietic cells in the disruption of hematopoiesis under homeostatic conditions.

APOBEC mutagenesis is a common process in normal human small intestine.

APOBEC mutational signatures SBS2 and SBS13 are common in many human cancer types. However, there is an incomplete understanding of its stimulus, when it occurs in the progression from normal to cancer cell and the APOBEC enzymes responsible. Here we whole-genome sequenced 342 microdissected normal epithelial crypts from the small intestines of 39 individuals and found that SBS2/SBS13 mutations were present in 17% of crypts, more frequent than most other normal tissues. Crypts with SBS2/SBS13 often had immediate crypt neighbors without SBS2/SBS13, suggesting that the underlying cause of SBS2/SBS13 is cell-intrinsic. APOBEC mutagenesis occurred in an episodic manner throughout the human lifespan, including in young children. APOBEC1 mRNA levels were very high in the small intestine epithelium, but low in the large intestine epithelium and other tissues. The results suggest that the high levels of SBS2/SBS13 in the small intestine are collateral damage from APOBEC1 fulfilling its physiological function of editing APOB mRNA.

Gata2-regulated Gfi1b expression controls endothelial programming during endothelial-to-hematopoietic transition.

The first hematopoietic stem cells (HSCs) are formed through endothelial-to-hematopoietic transition (EHT) during embryonic development. The transcription factor GATA2 is a crucial regulator of EHT and HSC function throughout life. Because patients with GATA2 haploinsufficiency have inborn mutations, prenatal defects are likely to influence disease development. In mice, Gata2 haploinsufficiency (Gata2+/-) reduces the number and functionality of embryonic hematopoietic stem and progenitor cells (HSPCs) generated through EHT. However, the embryonic HSPC pool is heterogeneous and the mechanisms underlying this defect in Gata2+/- embryos remain unclear. Here, we investigated whether Gata2 haploinsufficiency selectively affects a cellular subset undergoing EHT. We showed that Gata2+/- HSPCs initiate, but cannot fully activate, hematopoietic programming during EHT. In addition, due to the reduced activity of the endothelial repressor Gfi1b, Gata2+/- HSPCs cannot repress endothelial identity to complete maturation. Finally, we showed that hematopoietic-specific induction of gfi1b could restore HSC production in gata2b-null (gata2b-/-) zebrafish embryos. This study illustrates the pivotal role of Gata2 in the regulation of the transcriptional network governing HSPC identity throughout the EHT.

Prognostic Value of FLT3-Internal Tandem Duplication Residual Disease in Acute Myeloid Leukemia.

PURPOSE: The applicability of FLT3-internal tandem duplications (FLT3-ITD) for assessing measurable residual disease (MRD) in acute myeloid leukemia (AML) in complete remission (CR) has been hampered by patient-specific duplications and potential instability of FLT3-ITD during relapse. Here, we comprehensively investigated the impact of next-generation sequencing (NGS)-based FLT3-ITD MRD detection on treatment outcome in a cohort of patients with newly diagnosed AML in relation to established prognostic factors at diagnosis and other MRD measurements, ie, mutant NPM1 and multiparameter flow cytometry. METHODS: In 161 patients with de novo FLT3-ITD AML, NGS was performed at diagnosis and in CR after intensive remission induction treatment. FLT3-ITD MRD status was correlated with the cumulative incidence of relapse and overall survival (OS). RESULTS: NGS-based FLT3-ITD MRD was present in 47 of 161 (29%) patients with AML. Presence of FLT3-ITD MRD was associated with increased risk of relapse (4-year cumulative incidence of relapse, 75% FLT3-ITD MRD v 33% no FLT3-ITD MRD; P < .001) and inferior OS (4-year OS, 31% FLT3-ITD MRD v 57% no FLT3-ITD MRD; P < .001). In multivariate analysis, detection of FLT3-ITD MRD in CR confers independent prognostic significance for relapse (hazard ratio, 3.55; P < .001) and OS (hazard ratio 2.51; P = .002). Strikingly, FLT3-ITD MRD exceeds the prognostic value of most generally accepted clinical and molecular prognostic factors, including the FLT3-ITD allelic ratio at diagnosis and MRD assessment by NGS-based mutant NPM1 detection or multiparameter flow cytometry. CONCLUSION: NGS-based detection of FLT3-ITD MRD in CR identifies patients with AML with profound risk of relapse and death that outcompetes the significance of most established prognostic factors at diagnosis and during therapy, and furnishes support for FLT3-ITD as a clinically relevant biomarker for dynamic disease risk assessment in AML.

Myeloid cells promote interferon signaling-associated deterioration of the hematopoietic system.

Innate and adaptive immune cells participate in the homeostatic regulation of hematopoietic stem cells (HSCs). Here, we interrogate the contribution of myeloid cells, the most abundant cell type in the mammalian bone marrow, in a clinically relevant mouse model of neutropenia. Long-term genetic depletion of neutrophils and eosinophils results in activation of multipotent progenitors but preservation of HSCs. Depletion of myeloid cells abrogates HSC expansion, loss of serial repopulation and lymphoid reconstitution capacity and remodeling of HSC niches, features previously associated with hematopoietic aging. This is associated with mitigation of interferon signaling in both HSCs and their niches via reduction of NK cell number and activation. These data implicate myeloid cells in the functional decline of hematopoiesis, associated with activation of interferon signaling via a putative neutrophil-NK cell axis. Innate immunity may thus come at the cost of system deterioration through enhanced chronic inflammatory signaling to stem cells and their niches.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Clonal diversification and histogenesis of malignant germ cell tumours.

Germ cell tumours (GCTs) are a collection of benign and malignant neoplasms derived from primordial germ cells. They are uniquely able to recapitulate embryonic and extraembryonic tissues, which carries prognostic and therapeutic significance. The developmental pathways underpinning GCT initiation and histogenesis are incompletely understood. Here, we study the relationship of histogenesis and clonal diversification in GCTs by analysing the genomes and transcriptomes of 547 microdissected histological units. We find no correlation between genomic and histological heterogeneity. However, we identify unifying features including the retention of fetal developmental transcripts across tissues, expression changes on chromosome 12p, and a conserved somatic evolutionary sequence of whole genome duplication followed by clonal diversification. While this pattern is preserved across all GCTs, the developmental timing of the duplication varies between prepubertal and postpubertal cases. In addition, tumours of younger children exhibit distinct substitution signatures which may lend themselves as potential biomarkers for risk stratification. Our findings portray the extensive diversification of GCT tissues and genetic subclones as randomly distributed, while identifying overarching transcriptional and genomic features.

Inherited MUTYH mutations cause elevated somatic mutation rates and distinctive mutational signatures in normal human cells.

Cellular DNA damage caused by reactive oxygen species is repaired by the base excision repair (BER) pathway which includes the DNA glycosylase MUTYH. Inherited biallelic MUTYH mutations cause predisposition to colorectal adenomas and carcinoma. However, the mechanistic progression from germline MUTYH mutations to MUTYH-Associated Polyposis (MAP) is incompletely understood. Here, we sequence normal tissue DNAs from 10 individuals with MAP. Somatic base substitution mutation rates in intestinal epithelial cells were elevated 2 to 4-fold in all individuals, except for one showing a 31-fold increase, and were also increased in other tissues. The increased mutation burdens were of multiple mutational signatures characterised by C > A changes. Different mutation rates and signatures between individuals are likely due to different MUTYH mutations or additional inherited mutations in other BER pathway genes. The elevated base substitution rate in normal cells likely accounts for the predisposition to neoplasia in MAP. Despite ubiquitously elevated mutation rates, individuals with MAP do not display overt evidence of premature ageing. Thus, accumulation of somatic mutations may not be sufficient to cause the global organismal functional decline of ageing.

Clonal evolution after treatment pressure in multiple myeloma: heterogenous genomic aberrations and transcriptomic convergence.

We investigated genomic and transcriptomic changes in paired tumor samples of 29 in-house multiple myeloma (MM) patients and 28 patients from the MMRF CoMMpass study before and after treatment. A change in clonal composition was found in 46/57 (82%) of patients, and single-nucleotide variants (SNVs) increased from median 67 to 86. The highest increase in prevalence of genetic aberrations was found in RAS genes (60% to 72%), amp1q21 (18% to 35%), and TP53 (9% to 18%). The SBS-MM1 mutation signature was detected both in patients receiving high and low dose melphalan. A total of 2589 genes were differentially expressed between early and late samples (FDR < 0.05). Gene set enrichment analysis (GSEA) showed increased expression of E2F, MYC, and glycolysis pathways and a decreased expression in TNF-NFkB and TGFbeta pathways in late compared to early stage. Single sample GSEA (ssGSEA) scores of differentially expressed pathways revealed that these changes were most evident in end-stage disease. Increased expression of several potentially targetable genes was found at late disease stages, including cancer-testis antigens, XPO1 and ABC transporters. Our study demonstrates a transcriptomic convergence of pathways supporting increased proliferation and metabolism during disease progression in MM.

Mutational landscape of normal epithelial cells in Lynch Syndrome patients.

Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole genome sequence 152 crypts from normal and neoplastic epithelial tissues from 10 LS patients. In normal tissues the repertoire of mutational processes and mutation rates is similar to that found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and MMR deficiency-associated mutational signatures is identified, which may represent a very early stage of LS pathogenesis. Phylogenetic trees of tumour crypts indicate that the most recent ancestor cell of each tumour is already MMR deficient and has experienced multiple cycles of clonal evolution. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

Somatic mutation rates scale with lifespan across mammals.

The rates and patterns of somatic mutation in normal tissues are largely unknown outside of humans1-7. Comparative analyses can shed light on the diversity of mutagenesis across species, and on long-standing hypotheses about the evolution of somatic mutation rates and their role in cancer and ageing. Here we performed whole-genome sequencing of 208 intestinal crypts from 56 individuals to study the landscape of somatic mutation across 16 mammalian species. We found that somatic mutagenesis was dominated by seemingly endogenous mutational processes in all species, including 5-methylcytosine deamination and oxidative damage. With some differences, mutational signatures in other species resembled those described in humans8, although the relative contribution of each signature varied across species. Notably, the somatic mutation rate per year varied greatly across species and exhibited a strong inverse relationship with species lifespan, with no other life-history trait studied showing a comparable association. Despite widely different life histories among the species we examined-including variation of around 30-fold in lifespan and around 40,000-fold in body mass-the somatic mutation burden at the end of lifespan varied only by a factor of around 3. These data unveil common mutational processes across mammals, and suggest that somatic mutation rates are evolutionarily constrained and may be a contributing factor in ageing.

Molecular characterization of mutant TP53 acute myeloid leukemia and high-risk myelodysplastic syndrome.

Substantial heterogeneity within mutant TP53 acute myeloid leukemia (AML) and myelodysplastic syndrome with excess of blast (MDS-EB) precludes the exact assessment of prognostic impact for individual patients. We performed in-depth clinical and molecular analysis of mutant TP53 AML and MDS-EB to dissect the molecular characteristics in detail and determine its impact on survival. We performed next-generation sequencing on 2200 AML/MDS-EB specimens and assessed the TP53 mutant allelic status (mono- or bi-allelic), the number of TP53 mutations, mutant TP53 clone size, concurrent mutations, cytogenetics, and mutant TP53 molecular minimal residual disease and studied the associations of these characteristics with overall survival. TP53 mutations were detected in 230 (10.5%) patients with AML/MDS-EB with a median variant allele frequency of 47%. Bi-allelic mutant TP53 status was observed in 174 (76%) patients. Multiple TP53 mutations were found in 49 (21%) patients. Concurrent mutations were detected in 113 (49%) patients. No significant difference in any of the aforementioned molecular characteristics of mutant TP53 was detected between AML and MDS-EB. Patients with mutant TP53 have a poor outcome (2-year overall survival, 12.8%); however, no survival difference between AML and MDS-EB was observed. Importantly, none of the molecular characteristics were significantly associated with survival in mutant TP53 AML/MDS-EB. In most patients, TP53 mutations remained detectable in complete remission by deep sequencing (73%). Detection of residual mutant TP53 was not associated with survival. Mutant TP53 AML and MDS-EB do not differ with respect to molecular characteristics and survival. Therefore, mutant TP53 AML/MDS-EB should be considered a distinct molecular disease entity.