Genetic deletion of JAM-C in pre-leukemic cells rewires leukemic stem cell gene expression program in AML.

The leukemic stem cell (LSC) score LSC-17 based on a stemness-related gene expression signature is an indicator of poor disease outcome in acute myeloid leukemia (AML). However, our understanding of the relationships between LSC and pre-leukemic cells is still incomplete. In particular, it is not known whether "niche-anchoring" of pre-leukemic cell affects disease evolution. To address this issue, we conditionally inactivated the adhesion molecule JAM-C expressed by haematopoietic stem cells (HSC) and LSC in an inducible iMLL-AF9-driven AML mouse model. Deletion of Jam3 (encoding JAM-C) before induction of the leukemia-initiating iMLL-AF9 fusion resulted in a shift from long term to short term-HSC expansion, without affecting disease initiation and progression. In vitro experiments showed that JAM-C controlled leukemic cell nesting irrespective of the bone marrow stromal cells used. RNA sequencing performed on leukemic HSC isolated from diseased mice revealed that genes upregulated in Jam3-deficient animals belonged to Activation Protein-1 (AP-1) and TNF-/NFB pathways. Human orthologs of dysregulated genes allowed to identify a score based on AP-1/TNF-a gene expression that was distinct and complementary from LSC-17 score. Sub-stratification of AML patients with LSC-17 and AP-1/TNF-genes signature defined four groups with median survival ranging from below one year to a median not reached after 8 years. Finally, coculture experiments showed that AP-1 activation in leukemic cells was dependent on the nature of stromal cells. Altogether, our results identify the AP-1/TNF- gene signature as a proxy of LSC anchoring in specific bone marrow niches which improves the prognosis value of the LSC-17 score. NCT02320656.

Immunotherapeutic targeting of surfaceome heterogeneity in AML.

Immunotherapy remains underexploited in acute myeloid leukemia (AML) compared to other hematological malignancies. Currently, gemtuzumab ozogamicin is the only therapeutic antibody approved for this disease. Here, to identify potential targets for immunotherapeutic intervention, we analyze the surface proteome of 100 genetically diverse primary human AML specimens for the identification of cell surface proteins and conduct single-cell transcriptome analyses on a subset of these specimens to assess antigen expression at the sub-population level. Through this comprehensive effort, we successfully identify numerous antigens and markers preferentially expressed by primitive AML cells. Many identified antigens are targeted by therapeutic antibodies currently under clinical evaluation for various cancer types, highlighting the potential therapeutic value of the approach. Importantly, this initiative uncovers AML heterogeneity at the surfaceome level, identifies several antigens and potential primitive cell markers characterizing AML subgroups, and positions immunotherapy as a promising approach to target AML subgroup specificities.

SF3B1 mutations provide genetic vulnerability to copper ionophores in human acute myeloid leukemia.

In a phenotypical screen of 56 acute myeloid leukemia (AML) patient samples and using a library of 10,000 compounds, we identified a hit with increased sensitivity toward SF3B1-mutated and adverse risk AMLs. Through structure-activity relationship studies, this hit was optimized into a potent, specific, and nongenotoxic molecule called UM4118. We demonstrated that UM4118 acts as a copper ionophore that initiates a mitochondrial-based noncanonical form of cell death known as cuproptosis. CRISPR-Cas9 loss-of-function screen further revealed that iron-sulfur cluster (ISC) deficiency enhances copper-mediated cell death. Specifically, we found that loss of the mitochondrial ISC transporter ABCB7 is synthetic lethal to UM4118. ABCB7 is misspliced and down-regulated in SF3B1-mutated leukemia, creating a vulnerability to copper ionophores. Accordingly, ABCB7 overexpression partially rescued SF3B1-mutated cells to copper overload. Together, our work provides mechanistic insights that link ISC deficiency to cuproptosis, as exemplified by the high sensitivity of SF3B1-mutated AMLs. We thus propose SF3B1 mutations as a biomarker for future copper ionophore-based therapies.

DELE1 haploinsufficiency causes resistance to mitochondrial stress-induced apoptosis in monosomy 5/del(5q) AML.

Monosomy 5 and deletions of the chromosome 5q (-5/del(5q)) are recurrent events in de novo adult acute myeloid leukemia (AML), reaching up to 40% of cases in secondary AML. These chromosome anomalies are associated with TP53 mutations and with very poor prognosis. Using the large Leucegene genomic and transcriptomic dataset composed of 48 -5/del(5q) patient specimens and 367 control AML, we identified DELE1 - located in the common deleted region - as the most consistently downregulated gene in these leukemias. DELE1 encodes a mitochondrial protein recently characterized as the relay of mitochondrial stress to the cytosol through a newly defined OMA1-DELE1-HRI pathway which ultimately leads to the activation of ATF4, the master transcription factor of the integrated stress response. Here, we showed that the partial loss of DELE1 expression observed in -5/del(5q) patients was sufficient to significantly reduce the sensitivity to mitochondrial stress in AML cells. Overall, our results suggest that DELE1 haploinsufficiency could represent a new driver mechanism in -5/del(5q) AML.

HMGA2 expression defines a subset of human AML with immature transcriptional signature and vulnerability to G2/M inhibition.

High-mobility group AT-hook 2 (HMGA2) is a nonhistone chromatin-binding protein that is normally expressed in stem cells of various tissues and aberrantly detected in several tumor types. We recently observed that one-fourth of human acute myeloid leukemia (AML) specimens express HMGA2, which associates with a very poor prognosis. We present results indicating that HMGA2+ AMLs share a distinct transcriptional signature representing an immature phenotype. Using single-cell analyses, we showed that HMGA2 is expressed in CD34+ subsets of stem cells and early progenitors, whether normal or derived from AML specimens. Of interest, we found that one of the strongest gene expression signatures associated with HMGA2 in AML is the upregulation of G2/M checkpoint genes. Whole-genome CRISPR/Cas9 screening in HMGA2 overexpressing cells further revealed a synthetic lethal interaction with several G2/M checkpoint genes. Accordingly, small molecules that target G2/M proteins were preferentially active in vitro and in vivo on HMGA2+ AML specimens. Together, our findings suggest that HMGA2 is a key functional determinant in AML and is associated with stem cell features, G2/M status, and related drug sensitivity.

Cut-like homeobox 1 (CUX1) tumor suppressor gene haploinsufficiency induces apoptosis evasion to sustain myeloid leukemia.

While oncogenes promote tumorigenesis, they also induce deleterious cellular stresses, such as apoptosis, that cancer cells must combat by coopting adaptive responses. Whether tumor suppressor gene haploinsufficiency leads to such phenomena and their mechanistic basis is unclear. Here, we demonstrate that elevated levels of the anti-apoptotic factor, CASP8 and FADD-like apoptosis regulator (CFLAR), promotes apoptosis evasion in acute myeloid leukemia (AML) cells haploinsufficient for the cut-like homeobox 1 (CUX1) transcription factor, whose loss is associated with dismal clinical prognosis. Genome-wide CRISPR/Cas9 screening identifies CFLAR as a selective, acquired vulnerability in CUX1-deficient AML, which can be mimicked therapeutically using inhibitor of apoptosis (IAP) antagonists in murine and human AML cells. Mechanistically, CUX1 deficiency directly alleviates CUX1 repression of the CFLAR promoter to drive CFLAR expression and leukemia survival. These data establish how haploinsufficiency of a tumor suppressor is sufficient to induce advantageous anti-apoptosis cell survival pathways and concurrently nominate CFLAR as potential therapeutic target in these poor-prognosis leukemias.

UM171 Preserves Epigenetic Marks that Are Reduced in Ex Vivo Culture of Human HSCs via Potentiation of the CLR3-KBTBD4 Complex.

Human hematopoietic stem cells (HSCs) exhibit attrition of their self-renewal capacity when cultured ex vivo, a process that is partially reversed upon treatment with epigenetic modifiers, most notably inhibitors of histone deacetylases (HDACs) or lysine-specific demethylase LSD1. A recent study showed that the human HSC self-renewal agonist UM171 modulates the CoREST complex, leading to LSD1 degradation, whose inhibition mimics the activity of UM171. The mechanism underlying the UM171-mediated loss of CoREST function remains undetermined. We now report that UM171 potentiates the activity of a CULLIN3-E3 ubiquitin ligase (CRL3) complex whose target specificity is dictated by the poorly characterized Kelch/BTB domain protein KBTBD4. CRL3KBTBD4 targets components of the LSD1/RCOR1 corepressor complex for proteasomal degradation, hence re-establishing H3K4me2 and H3K27ac epigenetic marks, which are rapidly decreased upon ex vivo culture of human HSCs.

High frequency of germline RUNX1 mutations in patients with RUNX1-mutated AML.

RUNX1 is mutated in ∼10% of adult acute myeloid leukemia (AML). Although most RUNX1 mutations in this disease are believed to be acquired, they can also be germline. Indeed, germline RUNX1 mutations result in the well-described autosomal-dominant familial platelet disorder with predisposition to hematologic malignancies (RUNX1-FPD, FPD/AML, FPDMM); ∼44% of affected individuals progress to AML or myelodysplastic syndromes. Using the Leucegene RUNX1 AML patient group, we sought to investigate the proportion of germline vs acquired RUNX1 mutations in this cohort. Our results showed that 30% of RUNX1 mutations in our AML cohort are germline. Molecular profiling revealed higher frequencies of NRAS mutations and other mutations known to activate various signaling pathways in these patients with RUNX1 germline-mutated AML. Moreover, 2 patients (mother and son) had co-occurrence of RUNX1 and CEBPA germline mutations, with variable AML disease onset at 59 and 27 years, respectively. Together, these data suggest a higher than anticipated frequency of germline RUNX1 mutations in the Leucegene cohort and further highlight the importance of testing for RUNX1 mutations in instances in which allogeneic stem cell transplantation using a related donor is envisioned.

Genetic characterization of ABT-199 sensitivity in human AML.

Acute myeloid leukemias (AML) with mutations in the NPM1 gene (NPM1c+) represent a large AML subgroup with varying response to conventional treatment, highlighting the need to develop targeted therapeutic strategies for this disease. We screened a library of clinical drugs on a cohort of primary human AML specimens and identified the BCL2 inhibitor ABT-199 as a selective agent against NPM1c+ AML. Mutational analysis of ABT-199-sensitive and -resistant specimens identified mutations in NPM1, RAD21, and IDH1/IDH2 as predictors of ABT-199 sensitivity. Comparative transcriptome analysis further uncovered BCL2A1 as a potential mediator of ABT-199 resistance in AML. In line with our observation that RAD21 mutation confers sensitivity to ABT-199, we provide functional evidence that reducing RAD21 levels can sensitize AML cells to BCL2 inhibition. Moreover, we demonstrate that ABT-199 is able to produce selective anti-AML activity in vivo toward AML with mutations associated with compound sensitivity in PDX models. Overall, this study delineates the contribution of several genetic events to the response to ABT-199 and provides a rationale for the development of targeted therapies for NPM1c+ AML.

Human models of NUP98-KDM5A megakaryocytic leukemia in mice contribute to uncovering new biomarkers and therapeutic vulnerabilities.

Acute megakaryoblastic leukemia (AMKL) represents ∼10% of pediatric acute myeloid leukemia cases and typically affects young children (<3 years of age). It remains plagued with extremely poor treatment outcomes (<40% cure rates), mostly due to primary chemotherapy refractory disease and/or early relapse. Recurrent and mutually exclusive chimeric fusion oncogenes have been detected in 60% to 70% of cases and include nucleoporin 98 (NUP98) gene rearrangements, most commonly NUP98-KDM5A. Human models of NUP98-KDM5A-driven AMKL capable of faithfully recapitulating the disease have been lacking, and patient samples are rare, further limiting biomarkers and drug discovery. To overcome these impediments, we overexpressed NUP98-KDM5A in human cord blood hematopoietic stem and progenitor cells using a lentiviral-based approach to create physiopathologically relevant disease models. The NUP98-KDM5A fusion oncogene was a potent inducer of maturation arrest, sustaining long-term proliferative and progenitor capacities of engineered cells in optimized culture conditions. Adoptive transfer of NUP98-KDM5A-transformed cells into immunodeficient mice led to multiple subtypes of leukemia, including AMKL, that phenocopy human disease phenotypically and molecularly. The integrative molecular characterization of synthetic and patient NUP98-KDM5A AMKL samples revealed SELP, MPIG6B, and NEO1 as distinctive and novel disease biomarkers. Transcriptomic and proteomic analyses pointed to upregulation of the JAK-STAT signaling pathway in the model AMKL. Both synthetic models and patient-derived xenografts of NUP98-rearranged AMKL showed in vitro therapeutic vulnerability to ruxolitinib, a clinically approved JAK2 inhibitor. Overall, synthetic human AMKL models contribute to defining functional dependencies of rare genotypes of high-fatality pediatric leukemia, which lack effective and rationally designed treatments.

UM171 induces a homeostatic inflammatory-detoxification response supporting human HSC self-renewal.

Elucidation of the molecular cues required to balance adult stem cell self-renewal and differentiation is critical for advancing cellular therapies. Herein, we report that the hematopoietic stem cell (HSC) self-renewal agonist UM171 triggers a balanced pro- and anti-inflammatory/detoxification network that relies on NFKB activation and protein C receptor-dependent ROS detoxification, respectively. We demonstrate that within this network, EPCR serves as a critical protective component as its deletion hypersensitizes primitive hematopoietic cells to pro-inflammatory signals and ROS accumulation resulting in compromised stem cell function. Conversely, abrogation of the pro-inflammatory activity of UM171 through treatment with dexamethasone, cAMP elevating agents or NFkB inhibitors abolishes EPCR upregulation and HSC expansion. Together, these results show that UM171 stimulates ex vivo HSC expansion by establishing a critical balance between key pro- and anti-inflammatory mediators of self-renewal.

Mubritinib Targets the Electron Transport Chain Complex I and Reveals the Landscape of OXPHOS Dependency in Acute Myeloid Leukemia.

To identify therapeutic targets in acute myeloid leukemia (AML), we chemically interrogated 200 sequenced primary specimens. Mubritinib, a known ERBB2 inhibitor, elicited strong anti-leukemic effects in vitro and in vivo. In the context of AML, mubritinib functions through ubiquinone-dependent inhibition of electron transport chain (ETC) complex I activity. Resistance to mubritinib characterized normal CD34+ hematopoietic cells and chemotherapy-sensitive AMLs, which displayed transcriptomic hallmarks of hypoxia. Conversely, sensitivity correlated with mitochondrial function-related gene expression levels and characterized a large subset of chemotherapy-resistant AMLs with oxidative phosphorylation (OXPHOS) hyperactivity. Altogether, our work thus identifies an ETC complex I inhibitor and reveals the genetic landscape of OXPHOS dependency in AML.

Genetic risk factors for VIPN in childhood acute lymphoblastic leukemia patients identified using whole-exome sequencing.

AIM: To identify genetic markers associated with vincristine-induced peripheral neuropathy (VIPN) in childhood acute lymphoblastic leukemia. PATIENTS & METHODS: Whole-exome sequencing data were combined with exome-wide association study to identify predicted-functional germline variants associated with high-grade VIPN. Genotyping was then performed for top-ranked signals (n = 237), followed by validation in independent replication group (n = 405). RESULTS: Minor alleles of rs2781377/SYNE2 (p = 0.01) and rs10513762/MRPL47 (p = 0.01) showed increased risk, whereas that of rs3803357/BAHD1 had a protective effect (p = 0.007). Using a genetic model based on weighted genetic risk scores, an additive effect of combining these loci was observed (p = 0.003). The addition of rs1135989/ACTG1 further enhanced model performance (p = 0.0001). CONCLUSION: Variants in SYNE2, MRPL47 and BAHD1 genes are putative new risk factors for VIPN in childhood acute lymphoblastic leukemia.

Mutational dynamics of early and late relapsed childhood ALL: rapid clonal expansion and long-term dormancy.

Childhood acute lymphoblastic leukemia (cALL) is the most frequent pediatric cancer. Refractory/relapsed cALL presents a survival rate of ∼45% and is still one of the leading causes of death by disease among children. Mechanisms, such as clonal competition and evolutionary adaptation, govern treatment resistance. However, the underlying clonal dynamics leading to multiple relapses and differentiating early (<36 months postdiagnosis) from late relapse events remain elusive. Here, we use an integrative genome-based analysis combined with serial sampling of relapsed tumors (from primary tumor to ≤4 relapse events) from 19 pre-B-cell cALL patients (8 early and 11 late relapses) to assess the fitness of somatic mutations and infer their ancestral relationships. By quantifying both general clonal dynamics and newly acquired subclonal diversity, we show that 2 distinct evolutionary patterns govern early and late relapse: on one hand, a highly dynamic pattern, sustained by a putative defect of DNA repair processes, illustrating the quick emergence of fitter clones, and on the other hand, a quasi-inert evolution pattern, suggesting the escape from dormancy of leukemia stem cells likely spared from initial cytoreductive therapy. These results offer new insights into cALL relapse mechanisms and highlight the pressing need for adapted treatment strategies to circumvent resistance mechanisms.

Characterization of the microDNA through the response to chemotherapeutics in lymphoblastoid cell lines.

Recently, a new class of extrachromosomal circular DNA, called microDNA, was identified. They are on average 100 to 400 bp long and are derived from unique non-repetitive genomic regions with high gene density. MicroDNAs are thought to arise from DNA breaks associated with RNA metabolism or replication slippage. Given the paucity of information on this entirely novel phenomenon, we aimed to get an additional insight into microDNA features by performing the microDNA analysis in 20 independent human lymphoblastoid cell lines (LCLs) prior and after treatment with chemotherapeutic drugs. The results showed non-random genesis of microDNA clusters from the active regions of the genome. The size periodicity of 190 bp was observed, which matches DNA fragmentation typical for apoptotic cells. The chemotherapeutic drug-induced apoptosis of LCLs increased both number and size of clusters further suggesting that part of microDNAs could result from the programmed cell death. Interestingly, proportion of identified microDNA sequences has common loci of origin when compared between cell line experiments. While compatible with the original observation that microDNAs originate from a normal physiological process, obtained results imply complementary source of its production. Furthermore, non-random genesis of microDNAs depicted by redundancy between samples makes these entities possible candidates for new biomarker generation.

Whole-exome sequencing identified genetic risk factors for asparaginase-related complications in childhood ALL patients.

Allergy, pancreatitis and thrombosis are common side-effects of childhood acute lymphoblastic leukemia (ALL) treatment that are associated with the use of asparaginase (ASNase), a key component in most ALL treatment protocols. Starting with predicted functional germline variants obtained through whole-exome sequencing (WES) data of the Quebec childhood ALL cohort we performed exome-wide association studies with ASNase-related toxicities. A subset of top-ranking variants was further confirmed by genotyping (N=302) followed by validation in an independent replication group (N=282); except for thrombosis which was not available for that dataset. SNPs in 12 genes were associated with ASNase complications in discovery cohort including 3 that were associated with allergy, 3 with pancreatitis and 6 with thrombosis. The risk was further increased through combined SNPs effect (p≤0.002), suggesting synergistic interactions between the SNPs identified in each of the studied toxicities. Interestingly, rs3809849 in the MYBBP1A gene was associated with allergy (p= 0.0006), pancreatitis (p=0.002), thrombosis (p=0.02), event-free survival (p=0.02) and overall survival (p=0.003). Furthermore, rs11556218 in IL16 and rs34708521 in SPEF2 were both associated with thrombosis (p=0.01 and p=0.03, respectively) and pancreatitis (p=0.02). The association of SNPs in MYBBP1A, SPEF2 and IL16 geneswith pancreatitis was replicated in the validation cohort (p ≤0.05) as well as in combined cohort (p=0.0003, p=0.008 and p=0.02, respectively). The synergistic effect of combining risk loci had the highest power to predict the development of pancreatitis in both cohorts and was further potentiated in the combined cohort (p=1x10-8).The present work demonstrates that using WES data is a successful "hypothesis-free" strategy for identifying significant genetic markers modulating the effect of the treatment in childhood ALL.

SNooPer: a machine learning-based method for somatic variant identification from low-pass next-generation sequencing.

BACKGROUND: Next-generation sequencing (NGS) allows unbiased, in-depth interrogation of cancer genomes. Many somatic variant callers have been developed yet accurate ascertainment of somatic variants remains a considerable challenge as evidenced by the varying mutation call rates and low concordance among callers. Statistical model-based algorithms that are currently available perform well under ideal scenarios, such as high sequencing depth, homogeneous tumor samples, high somatic variant allele frequency (VAF), but show limited performance with sub-optimal data such as low-pass whole-exome/genome sequencing data. While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. RESULTS: For these reasons, we developed SNooPer, a versatile machine learning approach that uses Random Forest classification models to accurately call somatic variants in low-depth sequencing data. SNooPer uses a subset of variant positions from the sequencing output for which the class, true variation or sequencing error, is known to train the data-specific model. Here, using a real dataset of 40 childhood acute lymphoblastic leukemia patients, we show how the SNooPer algorithm is not affected by low coverage or low VAFs, and can be used to reduce overall sequencing costs while maintaining high specificity and sensitivity to somatic variant calling. When compared to three benchmarked somatic callers, SNooPer demonstrated the best overall performance. CONCLUSIONS: While the goal of any cancer sequencing project is to identify a relevant, and limited, set of somatic variants for further sequence/functional validation, the inherently complex nature of cancer genomes combined with technical issues directly related to sequencing and alignment can affect either the specificity and/or sensitivity of most callers. The flexibility of SNooPer's random forest protects against technical bias and systematic errors, and is appealing in that it does not rely on user-defined parameters. The code and user guide can be downloaded at https://sourceforge.net/projects/snooper/ .

Genomic characterization of pediatric T-cell acute lymphoblastic leukemia reveals novel recurrent driver mutations.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematologic malignancy with variable prognosis. It represents 15% of diagnosed pediatric ALL cases and has a threefold higher incidence among males. Many recurrent alterations have been identified and help define molecular subgroups of T-ALL, however the full range of events involved in driving transformation remain to be defined. Using an integrative approach combining genomic and transcriptomic data, we molecularly characterized 30 pediatric T-ALLs and identified common recurrent T-ALL targets such as FBXW7, JAK1, JAK3, PHF6, KDM6A and NOTCH1 as well as novel candidate T-ALL driver mutations including the p.R35L missense mutation in splicesome factor U2AF1 found in 3 patients and loss of function mutations in the X-linked tumor suppressor genes MED12 (frameshit mutation p.V167fs, splice site mutation g.chrX:70339329T>C, missense mutation p.R1989H) and USP9X (nonsense mutation p.Q117*). In vitro functional studies further supported the putative role of these novel T-ALL genes in driving transformation. U2AF1 p.R35L was shown to induce aberrant splicing of downstream target genes, and shRNA knockdown of MED12 and USP9X was shown to confer resistance to apoptosis following T-ALL relevant chemotherapy drug treatment in Jurkat leukemia cells. Interestingly, nearly 60% of novel candidate driver events were identified among immature T-ALL cases, highlighting the underlying genomic complexity of pediatric T-ALL, and the need for larger integrative studies to decipher the mechanisms that contribute to its various subtypes and provide opportunities to refine patient stratification and treatment.

CLIC5: a novel ETV6 target gene in childhood acute lymphoblastic leukemia.

The most common rearrangement in childhood precursor B-cell acute lymphoblastic leukemia is the t(12;21)(p13;q22) translocation resulting in the ETV6-AML1 fusion gene. A frequent concomitant event is the loss of the residual ETV6 allele suggesting a critical role for the ETV6 transcriptional repressor in the etiology of this cancer. However, the precise mechanism through which loss of functional ETV6 contributes to disease pathogenesis is still unclear. To investigate the impact of ETV6 loss on the transcriptional network and to identify new transcriptional targets of ETV6, we used whole transcriptome analysis of both pre-B leukemic cell lines and patients combined with chromatin immunoprecipitation. Using this integrative approach, we identified 4 novel direct ETV6 target genes: CLIC5, BIRC7, ANGPTL2 and WBP1L To further evaluate the role of chloride intracellular channel protein CLIC5 in leukemogenesis, we generated cell lines overexpressing CLIC5 and demonstrated an increased resistance to hydrogen peroxide-induced apoptosis. We further described the implications of CLIC5's ion channel activity in lysosomal-mediated cell death, possibly by modulating the function of the transferrin receptor with which it colocalizes intracellularly. For the first time, we showed that loss of ETV6 leads to significant overexpression of CLIC5, which in turn leads to decreased lysosome-mediated apoptosis. Our data suggest that heightened CLIC5 activity could promote a permissive environment for oxidative stress-induced DNA damage accumulation, and thereby contribute to leukemogenesis.

A novel somatic mutation in ACD induces telomere lengthening and apoptosis resistance in leukemia cells.

BACKGROUND: The identification of oncogenic driver mutations has largely relied on the assumption that genes that exhibit more mutations than expected by chance are more likely to play an active role in tumorigenesis. Major cancer sequencing initiatives have therefore focused on recurrent mutations that are more likely to be drivers. However, in specific genetic contexts, low frequency mutations may also be capable of participating in oncogenic processes. Reliable strategies for identifying these rare or even patient-specific (private) mutations are needed in order to elucidate more personalized approaches to cancer diagnosis and treatment. METHODS: Here we performed whole-exome sequencing on three cases of childhood pre-B acute lymphoblastic leukemia (cALL), representing three cytogenetically-defined subgroups (high hyperdiploid, t(12;21) translocation, and cytogenetically normal). We applied a data reduction strategy to identify both common and rare/private somatic events with high functional potential. Top-ranked candidate mutations were subsequently validated at high sequencing depth on an independent platform and in vitro expression assays were performed to evaluate the impact of identified mutations on cell growth and survival. RESULTS: We identified 6 putatively damaging non-synonymous somatic mutations among the three cALL patients. Three of these mutations were well-characterized common cALL mutations involved in constitutive activation of the mitogen-activated protein kinase pathway (FLT3 p.D835Y, NRAS p.G13D, BRAF p.G466A). The remaining three patient-specific mutations (ACD p.G223V, DOT1L p.V114F, HCFC1 p.Y103H) were novel mutations previously undescribed in public cancer databases. Cytotoxicity assays demonstrated a protective effect of the ACD p.G223V mutation against apoptosis in leukemia cells. ACD plays a key role in protecting telomeres and recruiting telomerase. Using a telomere restriction fragment assay, we also showed that this novel mutation in ACD leads to increased telomere length in leukemia cells. CONCLUSION: This study identified ACD as a novel gene involved in cALL and points to a functional role for ACD in enhancing leukemia cell survival. These results highlight the importance of rare/private somatic mutations in understanding cALL etiology, even within well-characterized molecular subgroups.