Error-prone repair of stalled replication forks drives mutagenesis and loss of heterozygosity in haploinsufficient BRCA1 cells.

Germline mutations in the BRCA genes are associated with a higher risk of carcinogenesis, which is linked to an increased mutation rate and loss of the second unaffected BRCA allele (loss of heterozygosity, LOH). However, the mechanisms triggering mutagenesis are not clearly understood. The BRCA genes contain high numbers of repetitive DNA sequences. We detected replication forks stalling, DNA breaks, and deletions at these sites in haploinsufficient BRCA cells, thus identifying the BRCA genes as fragile sites. Next, we found that stalled forks are repaired by error-prone pathways, such as microhomology-mediated break-induced replication (MMBIR) in haploinsufficient BRCA1 breast epithelial cells. We detected MMBIR mutations in BRCA1 tumor cells and noticed deletions-insertions (>50 bp) at the BRCA1 genes in BRCA1 patients. Altogether, these results suggest that under stress, error-prone repair of stalled forks is upregulated and induces mutations, including complex genomic rearrangements at the BRCA genes (LOH), in haploinsufficient BRCA1 cells.

A nondegenerate code of deleterious variants in Mendelian loci contributes to complex disease risk.

Although countless highly penetrant variants have been associated with Mendelian disorders, the genetic etiologies underlying complex diseases remain largely unresolved. By mining the medical records of over 110 million patients, we examine the extent to which Mendelian variation contributes to complex disease risk. We detect thousands of associations between Mendelian and complex diseases, revealing a nondegenerate, phenotypic code that links each complex disorder to a unique collection of Mendelian loci. Using genome-wide association results, we demonstrate that common variants associated with complex diseases are enriched in the genes indicated by this "Mendelian code." Finally, we detect hundreds of comorbidity associations among Mendelian disorders, and we use probabilistic genetic modeling to demonstrate that Mendelian variants likely contribute nonadditively to the risk for a subset of complex diseases. Overall, this study illustrates a complementary approach for mapping complex disease loci and provides unique predictions concerning the etiologies of specific diseases.

Super-enhancer hypermutation alters oncogene expression in B cell lymphoma.

Diffuse large B cell lymphoma (DLBCL) is the most common B cell non-Hodgkin lymphoma and remains incurable in around 40% of patients. Efforts to sequence the coding genome identified several genes and pathways that are altered in this disease, including potential therapeutic targets1-5. However, the non-coding genome of DLBCL remains largely unexplored. Here we show that active super-enhancers are highly and specifically hypermutated in 92% of samples from individuals with DLBCL, display signatures of activation-induced cytidine deaminase activity, and are linked to genes that encode B cell developmental regulators and oncogenes. As evidence of oncogenic relevance, we show that the hypermutated super-enhancers linked to the BCL6, BCL2 and CXCR4 proto-oncogenes prevent the binding and transcriptional downregulation of the corresponding target gene by transcriptional repressors, including BLIMP1 (targeting BCL6) and the steroid receptor NR3C1 (targeting BCL2 and CXCR4). Genetic correction of selected mutations restored repressor DNA binding, downregulated target gene expression and led to the counter-selection of cells containing corrected alleles, indicating an oncogenic dependency on the super-enhancer mutations. This pervasive super-enhancer mutational mechanism reveals a major set of genetic lesions deregulating gene expression, which expands the involvement of known oncogenes in DLBCL pathogenesis and identifies new deregulated gene targets of therapeutic relevance.

Clonal evolution mechanisms in NT5C2 mutant-relapsed acute lymphoblastic leukaemia.

Relapsed acute lymphoblastic leukaemia (ALL) is associated with resistance to chemotherapy and poor prognosis. Gain-of-function mutations in the 5'-nucleotidase, cytosolic II (NT5C2) gene induce resistance to 6-mercaptopurine and are selectively present in relapsed ALL. Yet, the mechanisms involved in NT5C2 mutation-driven clonal evolution during the initiation of leukaemia, disease progression and relapse remain unknown. Here we use a conditional-and-inducible leukaemia model to demonstrate that expression of NT5C2(R367Q), a highly prevalent relapsed-ALL NT5C2 mutation, induces resistance to chemotherapy with 6-mercaptopurine at the cost of impaired leukaemia cell growth and leukaemia-initiating cell activity. The loss-of-fitness phenotype of NT5C2+/R367Q mutant cells is associated with excess export of purines to the extracellular space and depletion of the intracellular purine-nucleotide pool. Consequently, blocking guanosine synthesis by inhibition of inosine-5'-monophosphate dehydrogenase (IMPDH) induced increased cytotoxicity against NT5C2-mutant leukaemia lymphoblasts. These results identify the fitness cost of NT5C2 mutation and resistance to chemotherapy as key evolutionary drivers that shape clonal evolution in relapsed ALL and support a role for IMPDH inhibition in the treatment of ALL.

A novel and highly effective mitochondrial uncoupling drug in T-cell leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignancy. Despite recent advances in treatments with intensified chemotherapy regimens, relapse rates and associated morbidities remain high. In this context, metabolic dependencies have emerged as a druggable opportunity for the treatment of leukemia. Here, we tested the antileukemic effects of MB1-47, a newly developed mitochondrial uncoupling compound. MB1-47 treatment in T-ALL cells robustly inhibited cell proliferation via both cytostatic and cytotoxic effects as a result of compromised mitochondrial energy and metabolite depletion, which severely impaired nucleotide biosynthesis. Mechanistically, acute treatment with MB1-47 in primary leukemias promoted adenosine monophosphate-activated serine/threonine protein kinase (AMPK) activation and downregulation of mammalian target of rapamycin (mTOR) signaling, stalling anabolic pathways that support leukemic cell survival. Indeed, MB1-47 treatment in mice harboring either murine NOTCH1-induced primary leukemias or human T-ALL patient-derived xenografts (PDXs) led to potent antileukemic effects with a significant extension in survival without overlapping toxicities. Overall, our findings demonstrate a critical role for mitochondrial oxidative phosphorylation in T-ALL and uncover MB1-47-driven mitochondrial uncoupling as a novel therapeutic strategy for the treatment of this disease.

Mutational signatures reveal ternary relationships between homologous recombination repair, APOBEC, and mismatch repair in gynecological cancers.

BACKGROUND: Revealing the impacts of endogenous and exogenous mutagenesis processes is essential for understanding the etiology of somatic genomic alterations and designing precise prognostication and treatment strategies for cancer. DNA repair deficiency is one of the main sources of endogenous mutagenesis and is increasingly recognized as a target for cancer therapeutics. The role and prevalence of mechanisms that underly different forms of DNA repair deficiencies and their interactions remain to be elucidated in gynecological malignancies. METHODS: We analyzed 1231 exomes and 268 whole-genomes from three major gynecological malignancies including uterine corpus endometrial carcinoma (UCEC) as well as ovarian and cervical cancers. We also analyzed data from 134 related cell lines. We extracted and compared de novo and refitted mutational signature profiles using complementary and confirmatory approaches and performed interaction analysis to detect co-occurring and mutually exclusive signatures. RESULTS: We found an inverse relationship between homologous recombination deficiency (HRd) and mismatch repair deficiency (MMRd). Moreover, APOBEC co-occurred with HRd but was mutually exclusive with MMRd. UCEC tumors were dominated by MMRd, yet a subset of them manifested the HRd and APOBEC signatures. Conversely, ovarian tumors were dominated by HRd, while a subset represented MMRd and APOBEC. In contrast to both, cervical tumors were dominated by APOBEC with a small subsets showing the POLE, HRd, and MMRd signatures. Although the type, prevalence, and heterogeneity of mutational signatures varied across the tumor types, the patterns of co-occurrence and exclusivity were consistently observed in all. Notably, mutational signatures in gynecological tumor cell lines reflected those detected in primary tumors. CONCLUSIONS: Taken together, these analyses indicate that application of mutation signature analysis not only advances our understanding of mutational processes and their interactions, but also it has the potential to stratify patients that could benefit from treatments available for tumors harboring distinct mutational signatures and to improve clinical decision-making for gynecological malignancies.

Subclonal NT5C2 mutations are associated with poor outcomes after relapse of pediatric acute lymphoblastic leukemia.

Activating mutations in cytosolic 5'-nucleotidase II (NT5C2) are considered to drive relapse formation in acute lymphoblastic leukemia (ALL) by conferring purine analog resistance. To examine the clinical effects of NT5C2 mutations in relapsed ALL, we analyzed NT5C2 in 455 relapsed B-cell precursor ALL patients treated within the ALL-REZ BFM 2002 relapse trial using sequencing and sensitive allele-specific real-time polymerase chain reaction. We detected 110 NT5C2 mutations in 75 (16.5%) of 455 B-cell precursor ALL relapses. Two-thirds of relapses harbored subclonal mutations and only one-third harbored clonal mutations. Event-free survival after relapse was inferior in patients with relapses with clonal and subclonal NT5C2 mutations compared with those without (19% and 25% vs 53%, P < .001). However, subclonal, but not clonal, NT5C2 mutations were associated with reduced event-free survival in multivariable analysis (hazard ratio, 1.89; 95% confidence interval, 1.28-2.69; P = .001) and with an increased rate of nonresponse to relapse treatment (subclonal 32%, clonal 12%, wild type 9%, P < .001). Nevertheless, 27 (82%) of 33 subclonal NT5C2 mutations became undetectable at the time of nonresponse or second relapse, and in 10 (71%) of 14 patients subclonal NT5C2 mutations were undetectable already after relapse induction treatment. These results show that subclonal NT5C2 mutations define relapses associated with high risk of treatment failure in patients and at the same time emphasize that their role in outcome is complex and goes beyond mutant NT5C2 acting as a targetable driver during relapse progression. Sensitive, prospective identification of NT5C2 mutations is warranted to improve the understanding and treatment of this aggressive ALL relapse subtype.

Tissue- and development-stage-specific mRNA and heterogeneous CNV signatures of human ribosomal proteins in normal and cancer samples.

We give results from a detailed analysis of human Ribosomal Protein (RP) levels in normal and cancer samples and cell lines from large mRNA, copy number variation and ribosome profiling datasets. After normalizing total RP mRNA levels per sample, we find highly consistent tissue specific RP mRNA signatures in normal and tumor samples. Multiple RP mRNA-subtypes exist in several cancers, with significant survival and genomic differences. Some RP mRNA variations among subtypes correlate with copy number loss of RP genes. In kidney cancer, RP subtypes map to molecular subtypes related to cell-of-origin. Pan-cancer analysis of TCGA data showed widespread single/double copy loss of RP genes, without significantly affecting survival. In several cancer cell lines, CRISPR-Cas9 knockout of RP genes did not affect cell viability. Matched RP ribosome profiling and mRNA data in humans and rodents stratified by tissue and development stage and were strongly correlated, showing that RP translation rates were proportional to mRNA levels. In a small dataset of human adult and fetal tissues, RP protein levels showed development stage and tissue specific heterogeneity of RP levels. Our results suggest that heterogeneous RP levels play a significant functional role in cellular physiology, in both normal and disease states.

JAK/STAT of all trades: linking inflammation with cancer development, tumor progression and therapy resistance.

Inflammation is at the forefront of carcinogenesis, tumor progression and resistance to therapy. The Janus kinase (JAK)/signal transducer and activator of transcription (STAT) signaling axis is a central pathway that mediates the cellular response to inflammation and contributes to carcinogenesis. The JAK/STAT pathway coordinates intercellular communication between tumor cells and their immune microenvironment, and JAK/STAT activation leads to the expression of a variety of proteins involved in cell proliferation, cell survival, stemness, self-renewal, evasion of immunosurveillance mechanisms and overall tumor progression. Activation of JAK/STAT signaling also mediates resistance to radiation therapy or cytotoxic agents and modulates tumor cell responses to molecularly targeted and immune modulating drugs. Despite extensive research focused on understanding its signaling mechanisms and downstream phenotypic and functional consequences in hematological disorders, the importance of JAK/STAT signaling in solid tumor initiation and progression has been underappreciated. We highlight the role of chronic inflammation in cancer, the epidemiological evidence for contribution of JAK/STAT to carcinogenesis, the current cancer prevention measures involving JAK/STAT inhibition and the impact of JAK/STAT signaling activity on cancer development, progression and treatment resistance. We also discuss recent therapeutic advances in targeting key factors within the JAK/STAT pathway with single agents and the use of these agents in combination with other targeted therapies and immune checkpoint inhibitors.

Gene expression of adipokines and adipokine receptors in the tumor microenvironment: associations of lower expression with more aggressive breast tumor features.

PURPOSE: Limited epidemiologic data are available on the expression of adipokines leptin (LEP) and adiponectin (ADIPOQ) and adipokine receptors (LEPR, ADIPOR1, ADIPOR2) in the breast tumor microenvironment (TME). The associations of gene expression of these biomarkers with tumor clinicopathology are not well understood. METHODS: NanoString multiplexed assays were used to assess the gene expression levels of LEP, LEPR, ADIPOQ, ADIPOR1, and ADIPOR2 within tumor tissues among 162 Black and 55 White women with newly diagnosed breast cancer. Multivariate mixed effects models were used to estimate associations of gene expression with breast tumor clinicopathology (overall and separately among Blacks). RESULTS: Black race was associated with lower gene expression of LEPR (P = 0.002) and ADIPOR1 (P = 0.01). Lower LEP, LEPR, and ADIPOQ gene expression were associated with higher tumor grade (P = 0.0007, P < 0.0001, and P < 0.0001, respectively) and larger tumor size (P < 0.0001, P = 0.0005, and P < 0.0001, respectively). Lower ADIPOQ expression was associated with ER- status (P = 0.0005), and HER2-enriched (HER2-E; P = 0.0003) and triple-negative (TN; P = 0.002) subtypes. Lower ADIPOR2 expression was associated with Ki67+ status (P = 0.0002), ER- status (P < 0.0001), PR- status (P < 0.0001), and TN subtype (P = 0.0002). Associations of lower adipokine and adipokine receptor gene expression with ER-, HER2-E, and TN subtypes were confirmed using data from The Cancer Genome Atlas (P-values < 0.005). CONCLUSION: These findings suggest that lower expression of ADIPOQ, ADIPOR2, LEP, and LEPR in the breast TME might be indicators of more aggressive breast cancer phenotypes. Validation of these findings are warranted to elucidate the role of the adipokines and adipokine receptors in long-term breast cancer prognosis.

All-FIT: allele-frequency-based imputation of tumor purity from high-depth sequencing data.

SUMMARY: Clinical sequencing aims to identify somatic mutations in cancer cells for accurate diagnosis and treatment. However, most widely used clinical assays lack patient-matched control DNA and additional analysis is needed to distinguish somatic and unfiltered germline variants. Such computational analyses require accurate assessment of tumor cell content in individual specimens. Histological estimates often do not corroborate with results from computational methods that are primarily designed for normal-tumor matched data and can be confounded by genomic heterogeneity and presence of sub-clonal mutations. Allele-frequency-based imputation of tumor (All-FIT) is an iterative weighted least square method to estimate specimen tumor purity based on the allele frequencies of variants detected in high-depth, targeted, clinical sequencing data. Using simulated and clinical data, we demonstrate All-FIT's accuracy and improved performance against leading computational approaches, highlighting the importance of interpreting purity estimates based on expected biology of tumors. AVAILABILITY AND IMPLEMENTATION: Freely available at http://software.khiabanian-lab.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Leukaemogenesis induced by an activating ß-catenin mutation in osteoblasts.

Cells of the osteoblast lineage affect the homing and the number of long-term repopulating haematopoietic stem cells, haematopoietic stem cell mobilization and lineage determination and B cell lymphopoiesis. Osteoblasts were recently implicated in pre-leukaemic conditions in mice. However, a single genetic change in osteoblasts that can induce leukaemogenesis has not been shown. Here we show that an activating mutation of ß-catenin in mouse osteoblasts alters the differentiation potential of myeloid and lymphoid progenitors leading to development of acute myeloid leukaemia with common chromosomal aberrations and cell autonomous progression. Activated ß-catenin stimulates expression of the Notch ligand jagged 1 in osteoblasts. Subsequent activation of Notch signalling in haematopoietic stem cell progenitors induces the malignant changes. Genetic or pharmacological inhibition of Notch signalling ameliorates acute myeloid leukaemia and demonstrates the pathogenic role of the Notch pathway. In 38% of patients with myelodysplastic syndromes or acute myeloid leukaemia, increased ß-catenin signalling and nuclear accumulation was identified in osteoblasts and these patients showed increased Notch signalling in haematopoietic cells. These findings demonstrate that genetic alterations in osteoblasts can induce acute myeloid leukaemia, identify molecular signals leading to this transformation and suggest a potential novel pharmacotherapeutic approach to acute myeloid leukaemia.

Activating mutations and translocations in the guanine exchange factor VAV1 in peripheral T-cell lymphomas.

Peripheral T-cell lymphomas (PTCLs) are a heterogeneous group of non-Hodgkin lymphomas frequently associated with poor prognosis and for which genetic mechanisms of transformation remain incompletely understood. Using RNA sequencing and targeted sequencing, here we identify a recurrent in-frame deletion (VAV1 Δ778-786) generated by a focal deletion-driven alternative splicing mechanism as well as novel VAV1 gene fusions (VAV1-THAP4, VAV1-MYO1F, and VAV1-S100A7) in PTCL. Mechanistically these genetic lesions result in increased activation of VAV1 catalytic-dependent (MAPK, JNK) and non-catalytic-dependent (nuclear factor of activated T cells, NFAT) VAV1 effector pathways. These results support a driver oncogenic role for VAV1 signaling in the pathogenesis of PTCL.

Mutational landscape, clonal evolution patterns, and role of RAS mutations in relapsed acute lymphoblastic leukemia.

Although multiagent combination chemotherapy is curative in a significant fraction of childhood acute lymphoblastic leukemia (ALL) patients, 20% of cases relapse and most die because of chemorefractory disease. Here we used whole-exome and whole-genome sequencing to analyze the mutational landscape at relapse in pediatric ALL cases. These analyses identified numerous relapse-associated mutated genes intertwined in chemotherapy resistance-related protein complexes. In this context, RAS-MAPK pathway-activating mutations in the neuroblastoma RAS viral oncogene homolog (NRAS), kirsten rat sarcoma viral oncogene homolog (KRAS), and protein tyrosine phosphatase, nonreceptor type 11 (PTPN11) genes were present in 24 of 55 (44%) cases in our series. Interestingly, some leukemias showed retention or emergence of RAS mutant clones at relapse, whereas in others RAS mutant clones present at diagnosis were replaced by RAS wild-type populations, supporting a role for both positive and negative selection evolutionary pressures in clonal evolution of RAS-mutant leukemia. Consistently, functional dissection of mouse and human wild-type and mutant RAS isogenic leukemia cells demonstrated induction of methotrexate resistance but also improved the response to vincristine in mutant RAS-expressing lymphoblasts. These results highlight the central role of chemotherapy-driven selection as a central mechanism of leukemia clonal evolution in relapsed ALL, and demonstrate a previously unrecognized dual role of RAS mutations as drivers of both sensitivity and resistance to chemotherapy.

MERIT reveals the impact of genomic context on sequencing error rate in ultra-deep applications.

BACKGROUND: Rapid progress in high-throughput sequencing (HTS) and the development of novel library preparation methods have improved the sensitivity of detecting mutations in heterogeneous samples, specifically in high-depth (> 500×) clinical applications. However, HTS methods are bounded by their technical and theoretical limitations and sequencing errors cannot be completely eliminated. Comprehensive quantification of the background noise can highlight both the efficiency and the limitations of any HTS methodology, and help differentiate true mutations at low abundance from artifacts. RESULTS: We introduce MERIT (Mutation Error Rate Inference Toolkit), designed for in-depth quantification of erroneous substitutions and small insertions and deletions. MERIT incorporates an all-inclusive variant caller and considers genomic context, including the nucleotides immediately at 5 'and 3 ', thereby establishing error rates for 96 possible substitutions as well as four single-base and 16 double-base indels. We applied MERIT to ultra-deep sequencing data (1,300,000 ×) obtained from the amplification of multiple clinically relevant loci, and showed a significant relationship between error rates and genomic contexts. In addition to observing significant difference between transversion and transition rates, we identified variations of more than 100-fold within each error type at high sequencing depths. For instance, T >G transversions in trinucleotide GTCs occurred 133.5 ± 65.9 more often than those in ATAs. Similarly, C >T transitions in GCGs were observed at 73.8 ± 10.5 higher rate than those in TCTs. We also devised an in silico approach to determine the optimal sequencing depth, where errors occur at rates similar to those of expected true mutations. Our analyses showed that increasing sequencing depth might improve sensitivity for detecting some mutations based on their genomic context. For example, T >G rate of error in GTCs did not change when sequenced beyond 10,000 ×; in contrast, T >G rate in TTAs consistently improved even at above 500,000 ×. CONCLUSIONS: Our results demonstrate significant variation in nucleotide misincorporation rates, and suggest that genomic context should be considered for comprehensive profiling of specimen-specific and sequencing artifacts in high-depth assays. This data provide strong evidence against assigning a single allele frequency threshold to call mutations, for it can result in substantial false positive as well as false negative variants, with important clinical consequences.

The genetics of nodal marginal zone lymphoma.

Nodal marginal zone lymphoma (NMZL) is a rare, indolent B-cell tumor that is distinguished from splenic marginal zone lymphoma (SMZL) by the different pattern of dissemination. NMZL still lacks distinct markers and remains orphan of specific cancer gene lesions. By combining whole-exome sequencing, targeted sequencing of tumor-related genes, whole-transcriptome sequencing, and high-resolution single nucleotide polymorphism array analysis, we aimed at disclosing the pathways that are molecularly deregulated in NMZL and we compare the molecular profile of NMZL with that of SMZL. These analyses identified a distinctive pattern of nonsilent somatic lesions in NMZL. In 35 NMZL patients, 41 genes were found recurrently affected in ≥3 (9%) cases, including highly prevalent molecular lesions of MLL2 (also known as KMT2D; 34%), PTPRD (20%), NOTCH2 (20%), and KLF2 (17%). Mutations of PTPRD, a receptor-type protein tyrosine phosphatase regulating cell growth, were enriched in NMZL across mature B-cell tumors, functionally caused the loss of the phosphatase activity of PTPRD, and were associated with cell-cycle transcriptional program deregulation and increased proliferation index in NMZL. Although NMZL shared with SMZL a common mutation profile, NMZL harbored PTPRD lesions that were otherwise absent in SMZL. Collectively, these findings provide new insights into the genetics of NMZL, identify PTPRD lesions as a novel marker for this lymphoma across mature B-cell tumors, and support the distinction of NMZL as an independent clinicopathologic entity within the current lymphoma classification.

Lactobacillus rhamnosus GG probiotic enteric regimen does not appreciably alter the gut microbiome or provide protection against GVHD after allogeneic hematopoietic stem cell transplantation.

Graft-versus-host disease (GVHD) is a major adverse effect associated with allogeneic stem cell transplant. Previous studies in mice indicated that administration of the probiotic Lactobacillus rhamnosus GG can reduce the incidence of GVHD after hematopoietic stem cell transplant. Here we report results from the first randomized probiotic enteric regimen trial in which allogenic hematopoietic stem cell patients were supplemented with Lactobacillus rhamnosus GG. Gut microbiome analysis confirmed a previously reported gut microbiome association with GVHD. However, the clinical trial was terminated when interim analysis did not detect an appreciable probiotic-related change in the gut microbiome or incidence of GVHD. Additional studies are necessary to determine whether probiotics can alter the incidence of GVHD after allogeneic stem cell transplant.

Genetic lesions associated with chronic lymphocytic leukemia chemo-refractoriness.

Fludarabine refractoriness (FR) represents an unsolved clinical problem of chronic lymphocytic leukemia (CLL) management. Although next-generation sequencing studies have led to the identification of a number of genes frequently mutated in FR-CLL, a comprehensive evaluation of the FR-CLL genome has not been reported. Toward this end, we studied 10 FR-CLLs by combining whole-exome sequencing and copy number aberration (CNA) analysis, which showed an average of 16.3 somatic mutations and 4 CNAs per sample. Screening of recurrently mutated genes in 48 additional FR-CLLs revealed that ~70% of FR-CLLs carry ≥1 mutation in genes previously associated with CLL clinical course, including TP53 (27.5%), NOTCH1 (24.1%), SF3B1 (18.9%), and BIRC3 (15.5%). In addition, this analysis showed that 10.3% of FR-CLL cases display mutations of the FAT1 gene, which encodes for a cadherin-like protein that negatively regulates Wnt signaling, consistent with a tumor suppressor role. The frequency of FAT1-mutated cases was significantly higher in FR-CLL than in unselected CLLs at diagnosis (10.3% vs 1.1%, P = .004), suggesting a role in the development of a high-risk phenotype. These findings have general implications for the mechanisms leading to FR and point to Wnt signaling as a potential therapeutic target in FR-CLL.

Clinical impact of small TP53 mutated subclones in chronic lymphocytic leukemia.

TP53 mutations are strong predictors of poor survival and refractoriness in chronic lymphocytic leukemia (CLL) and have direct implications for disease management. Clinical information on TP53 mutations is limited to lesions represented in >20% leukemic cells. Here, we tested the clinical impact and prediction of chemorefractoriness of very small TP53 mutated subclones. The TP53 gene underwent ultra-deep-next generation sequencing (NGS) in 309 newly diagnosed CLL. A robust bioinformatic algorithm was established for the highly sensitive detection of few TP53 mutated cells (down to 3 out of ∼1000 wild-type cells). Minor subclones were validated by independent approaches. Ultra-deep-NGS identified small TP53 mutated subclones in 28/309 (9%) untreated CLL that, due to their very low abundance (median allele frequency: 2.1%), were missed by Sanger sequencing. Patients harboring small TP53 mutated subclones showed the same clinical phenotype and poor survival (hazard ratio = 2.01; P = .0250) as those of patients carrying clonal TP53 lesions. By longitudinal analysis, small TP53 mutated subclones identified before treatment became the predominant population at the time of CLL relapse and anticipated the development of chemorefractoriness. This study provides a proof-of-principle that very minor leukemia subclones detected at diagnosis are an important driver of the subsequent disease course.