Contrasting roles of histone 3 lysine 27 demethylases in acute lymphoblastic leukaemia.

T-cell acute lymphoblastic leukaemia (T-ALL) is a haematological malignancy with a dismal overall prognosis, including a relapse rate of up to 25%, mainly because of the lack of non-cytotoxic targeted therapy options. Drugs that target the function of key epigenetic factors have been approved in the context of haematopoietic disorders, and mutations that affect chromatin modulators in a variety of leukaemias have recently been identified; however, 'epigenetic' drugs are not currently used for T-ALL treatment. Recently, we described that the polycomb repressive complex 2 (PRC2) has a tumour-suppressor role in T-ALL. Here we delineated the role of the histone 3 lysine 27 (H3K27) demethylases JMJD3 and UTX in T-ALL. We show that JMJD3 is essential for the initiation and maintenance of T-ALL, as it controls important oncogenic gene targets by modulating H3K27 methylation. By contrast, we found that UTX functions as a tumour suppressor and is frequently genetically inactivated in T-ALL. Moreover, we demonstrated that the small molecule inhibitor GSKJ4 (ref. 5) affects T-ALL growth, by targeting JMJD3 activity. These findings show that two proteins with a similar enzymatic function can have opposing roles in the context of the same disease, paving the way for treating haematopoietic malignancies with a new category of epigenetic inhibitors.

CONSERTING: integrating copy-number analysis with structural-variation detection.

We developed Copy Number Segmentation by Regression Tree in Next Generation Sequencing (CONSERTING), an algorithm for detecting somatic copy-number alteration (CNA) using whole-genome sequencing (WGS) data. CONSERTING performs iterative analysis of segmentation on the basis of changes in read depth and the detection of localized structural variations, with high accuracy and sensitivity. Analysis of 43 cancer genomes from both pediatric and adult patients revealed novel oncogenic CNAs, complex rearrangements and subclonal CNAs missed by alternative approaches.

The genetic basis of early T-cell precursor acute lymphoblastic leukaemia.

Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL.

CREBBP mutations in relapsed acute lymphoblastic leukaemia.

Relapsed acute lymphoblastic leukaemia (ALL) is a leading cause of death due to disease in young people, but the biological determinants of treatment failure remain poorly understood. Recent genome-wide profiling of structural DNA alterations in ALL have identified multiple submicroscopic somatic mutations targeting key cellular pathways, and have demonstrated substantial evolution in genetic alterations from diagnosis to relapse. However, DNA sequence mutations in ALL have not been analysed in detail. To identify novel mutations in relapsed ALL, we resequenced 300 genes in matched diagnosis and relapse samples from 23 patients with ALL. This identified 52 somatic non-synonymous mutations in 32 genes, many of which were novel, including the transcriptional coactivators CREBBP and NCOR1, the transcription factors ERG, SPI1, TCF4 and TCF7L2, components of the Ras signalling pathway, histone genes, genes involved in histone modification (CREBBP and CTCF), and genes previously shown to be targets of recurring DNA copy number alteration in ALL. Analysis of an extended cohort of 71 diagnosis-relapse cases and 270 acute leukaemia cases that did not relapse found that 18.3% of relapse cases had sequence or deletion mutations of CREBBP, which encodes the transcriptional coactivator and histone acetyltransferase CREB-binding protein (CREBBP, also known as CBP). The mutations were either present at diagnosis or acquired at relapse, and resulted in truncated alleles or deleterious substitutions in conserved residues of the histone acetyltransferase domain. Functionally, the mutations impaired histone acetylation and transcriptional regulation of CREBBP targets, including glucocorticoid responsive genes. Several mutations acquired at relapse were detected in subclones at diagnosis, suggesting that the mutations may confer resistance to therapy. These results extend the landscape of genetic alterations in leukaemia, and identify mutations targeting transcriptional and epigenetic regulation as a mechanism of resistance in ALL.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis.

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

CREST maps somatic structural variation in cancer genomes with base-pair resolution.

We developed 'clipping reveals structure' (CREST), an algorithm that uses next-generation sequencing reads with partial alignments to a reference genome to directly map structural variations at the nucleotide level of resolution. Application of CREST to whole-genome sequencing data from five pediatric T-lineage acute lymphoblastic leukemias (T-ALLs) and a human melanoma cell line, COLO-829, identified 160 somatic structural variations. Experimental validation exceeded 80%, demonstrating that CREST had a high predictive accuracy.

Cell of origin strongly influences genetic selection in a mouse model of T-ALL.

Identifying the normal cell from which a tumor originates is crucial to understanding the etiology of that cancer. However, retrospective identification of the cell of origin in cancer is challenging because of the accumulation of genetic and epigenetic changes in tumor cells. The biologic state of the cell of origin likely influences the genetic events that drive transformation. We directly tested this hypothesis by performing a Sleeping Beauty transposon mutagenesis screen in which common insertion sites were identified in tumors that were produced by mutagenesis of cells at varying time points throughout the T lineage. Mutation and gene expression data derived from these tumors were then compared with data obtained from a panel of 84 human T-cell acute lymphoblastic leukemia samples, including copy number alterations and gene expression profiles. This revealed that altering the cell of origin produces tumors that model distinct subtypes of human T-cell acute lymphoblastic leukemia, suggesting that even subtle changes in the cell of origin dramatically affect genetic selection in tumors. These findings have broad implications for the genetic analysis of human cancers as well as the production of mouse models of cancer.

The BCL11B tumor suppressor is mutated across the major molecular subtypes of T-cell acute lymphoblastic leukemia.

The BCL11B transcription factor is required for normal T-cell development, and has recently been implicated in the pathogenesis of T-cell acute lymphoblastic leukemia (T-ALL) induced by TLX overexpression or Atm deficiency. To comprehensively assess the contribution of BCL11B inactivation to human T-ALL, we performed DNA copy number and sequencing analyses of T-ALL diagnostic specimens, revealing monoallelic BCL11B deletions or missense mutations in 9% (n = 10 of 117) of cases. Structural homology modeling revealed that several of the BCL11B mutations disrupted the structure of zinc finger domains required for this transcription factor to bind DNA. BCL11B haploinsufficiency occurred across each of the major molecular subtypes of T-ALL, including early T-cell precursor, HOXA-positive, LEF1-inactivated, and TAL1-positive subtypes, which have differentiation arrest at diverse stages of thymocyte development. Our findings provide compelling evidence that BCL11B is a haploinsufficient tumor suppressor that collaborates with all major T-ALL oncogenic lesions in human thymocyte transformation.

A hyper-mutant of the unusual sigma70-Pr promoter bypasses synergistic ppGpp/DksA co-stimulation.

The activities of promoters can be temporally and conditionally regulated by mechanisms other than classical DNA-binding repressors and activators. One example is the inherently weak σ(70)-dependent Pr promoter that ultimately controls catabolism of phenolic compounds. The activity of Pr is up-regulated through the joint action of ppGpp and DksA that enhance the performance of RNA polymerase at this promoter. Here, we report a mutagenesis analysis that revealed substantial differences between Pr and other ppGpp/DksA co-stimulated promoters. In vitro transcription and RNA polymerase binding assays show that it is the T at the -11 position of the extremely suboptimal -10 element of Pr that underlies both poor binding of σ(70)-RNAP and a slow rate of open complex formation--the process that is accelerated by ppGpp and DksA. Our findings support the idea that collaborative action of ppGpp and DksA lowers the rate-limiting transition energy required for conversion between intermediates on the road to open complex formation.

Proteogenomic analysis of acute myeloid leukemia associates relapsed disease with reprogrammed energy metabolism both in adults and children.

Despite improvement of current treatment strategies and novel targeted drugs, relapse and treatment resistance largely determine the outcome for acute myeloid leukemia (AML) patients. To identify the underlying molecular characteristics, numerous studies have been aimed to decipher the genomic- and transcriptomic landscape of AML. Nevertheless, further molecular changes allowing malignant cells to escape treatment remain to be elucidated. Mass spectrometry is a powerful tool enabling detailed insights into proteomic changes that could explain AML relapse and resistance. Here, we investigated AML samples from 47 adult and 22 pediatric patients at serial time-points during disease progression using mass spectrometry-based in-depth proteomics. We show that the proteomic profile at relapse is enriched for mitochondrial ribosomal proteins and subunits of the respiratory chain complex, indicative of reprogrammed energy metabolism from diagnosis to relapse. Further, higher levels of granzymes and lower levels of the anti-inflammatory protein CR1/CD35 suggest an inflammatory signature promoting disease progression. Finally, through a proteogenomic approach, we detected novel peptides, which present a promising repertoire in the search for biomarkers and tumor-specific druggable targets. Altogether, this study highlights the importance of proteomic studies in holistic approaches to improve treatment and survival of AML patients.

Interpretable machine learning identifies paediatric Systemic Lupus Erythematosus subtypes based on gene expression data.

Transcriptomic analyses are commonly used to identify differentially expressed genes between patients and controls, or within individuals across disease courses. These methods, whilst effective, cannot encompass the combinatorial effects of genes driving disease. We applied rule-based machine learning (RBML) models and rule networks (RN) to an existing paediatric Systemic Lupus Erythematosus (SLE) blood expression dataset, with the goal of developing gene networks to separate low and high disease activity (DA1 and DA3). The resultant model had an 81% accuracy to distinguish between DA1 and DA3, with unsupervised hierarchical clustering revealing additional subgroups indicative of the immune axis involved or state of disease flare. These subgroups correlated with clinical variables, suggesting that the gene sets identified may further the understanding of gene networks that act in concert to drive disease progression. This included roles for genes (i) induced by interferons (IFI35 and OTOF), (ii) key to SLE cell types (KLRB1 encoding CD161), or (iii) with roles in autophagy and NF-κB pathway responses (CKAP4). As demonstrated here, RBML approaches have the potential to reveal novel gene patterns from within a heterogeneous disease, facilitating patient clinical and therapeutic stratification.

Transcriptomic analysis reveals proinflammatory signatures associated with acute myeloid leukemia progression.

Numerous studies have been performed over the last decade to exploit the complexity of genomic and transcriptomic lesions driving the initiation of acute myeloid leukemia (AML). These studies have helped improve risk classification and treatment options. Detailed molecular characterization of longitudinal AML samples is sparse, however; meanwhile, relapse and therapy resistance represent the main challenges in AML care. To this end, we performed transcriptome-wide RNA sequencing of longitudinal diagnosis, relapse, and/or primary resistant samples from 47 adult and 23 pediatric AML patients with known mutational background. Gene expression analysis revealed the association of short event-free survival with overexpression of GLI2 and IL1R1, as well as downregulation of ST18. Moreover, CR1 downregulation and DPEP1 upregulation were associated with AML relapse both in adults and children. Finally, machine learning-based and network-based analysis identified overexpressed CD6 and downregulated INSR as highly copredictive genes depicting important relapse-associated characteristics among adult patients with AML. Our findings highlight the importance of a tumor-promoting inflammatory environment in leukemia progression, as indicated by several of the herein identified differentially expressed genes. Together, this knowledge provides the foundation for novel personalized drug targets and has the potential to maximize the benefit of current treatments to improve cure rates in AML.

What Is Abnormal in Normal Karyotype Acute Myeloid Leukemia in Children? Analysis of the Mutational Landscape and Prognosis of the TARGET-AML Cohort.

Normal karyotype acute myeloid leukemia (NK-AML) constitutes 20-25% of pediatric AML and detailed molecular analysis is essential to unravel the genetic background of this group. Using publicly available sequencing data from the TARGET-AML initiative, we investigated the mutational landscape of NK-AML in comparison with abnormal karyotype AML (AK-AML). In 164 (97.6%) of 168 independent NK-AML samples, at least one somatic protein-coding mutation was identified using whole-genome or targeted capture sequencing. We identified a unique mutational landscape of NK-AML characterized by a higher prevalence of mutated CEBPA, FLT3, GATA2, NPM1, PTPN11, TET2, and WT1 and a lower prevalence of mutated KIT, KRAS, and NRAS compared with AK-AML. Mutated CEBPA often co-occurred with mutated GATA2, whereas mutated FLT3 co-occurred with mutated WT1 and NPM1. In multivariate regression analysis, we identified younger age, WBC count ≥50 × 109/L, FLT3-internal tandem duplications, and mutated WT1 as independent predictors of adverse prognosis and mutated NPM1 and GATA2 as independent predictors of favorable prognosis in NK-AML. In conclusion, NK-AML in children is characterized by a unique mutational landscape which impacts the disease outcome.

Measurable Residual Disease Monitoring of SPAG6, ST18, PRAME, and XAGE1A Expression in Peripheral Blood May Detect Imminent Relapse in Childhood Acute Myeloid Leukemia.

Overexpressed genes may be useful for monitoring of measurable residual disease (MRD) in patients with childhood acute myeloid leukemia (AML) without a leukemia-specific target. The normal expression of five leukemia-associated genes (SPAG6, ST18, MSLN, PRAME, XAGE1A) was defined in children without hematologic disease (n = 53) and children with suspected infection (n = 90). Gene expression at AML diagnosis (n=50) and during follow-up (n = 21) was compared with child-specific reference values. At diagnosis, 34/50 children (68%) had high expression of at least one of the five genes, and so did 16/31 children (52%) without a leukemia-specific target. Gene expression was quantified in 110 peripheral blood (PB) samples (median, five samples/patient; range, 1 to 10) during follow-up in 21 patients with high expression at diagnosis. All nine patients with PB sampling performed within 100 days of disease recurrence displayed overexpression of SPAG6, ST18, PRAME, or XAGE1A at a median of 2 months (range, 0.6 to 9.6 months) before hematologic relapse, whereas MSLN did not reach expression above normal prior to hematologic relapse. Only 1 of 130 (0.8%) follow-up analyses performed in 10 patients in continuous complete remission had transient expression above normal. SPAG6, ST18, PRAME, and XAGE1A expression in PB may predict relapse in childhood AML patients and facilitate MRD monitoring in most patients without a leukemia-specific target.

Genomic characterization of relapsed acute myeloid leukemia reveals novel putative therapeutic targets.

Relapse is the leading cause of death of adult and pediatric patients with acute myeloid leukemia (AML). Numerous studies have helped to elucidate the complex mutational landscape at diagnosis of AML, leading to improved risk stratification and new therapeutic options. However, multi-whole-genome studies of adult and pediatric AML at relapse are necessary for further advances. To this end, we performed whole-genome and whole-exome sequencing analyses of longitudinal diagnosis, relapse, and/or primary resistant specimens from 48 adult and 25 pediatric patients with AML. We identified mutations recurrently gained at relapse in ARID1A and CSF1R, both of which represent potentially actionable therapeutic alternatives. Further, we report specific differences in the mutational spectrum between adult vs pediatric relapsed AML, with MGA and H3F3A p.Lys28Met mutations recurrently found at relapse in adults, whereas internal tandem duplications in UBTF were identified solely in children. Finally, our study revealed recurrent mutations in IKZF1, KANSL1, and NIPBL at relapse. All of the mentioned genes have either never been reported at diagnosis in de novo AML or have been reported at low frequency, suggesting important roles for these alterations predominantly in disease progression and/or resistance to therapy. Our findings shed further light on the complexity of relapsed AML and identified previously unappreciated alterations that may lead to improved outcomes through personalized medicine.

Generation of human acute lymphoblastic leukemia xenografts for use in oncology drug discovery.

The establishment of reproducible mouse models of acute lymphoblastic leukemia (ALL) is necessary to provide in vivo therapeutic test systems that recapitulate human ALL, and for amplification of limited amounts of primary tumor material. A popular assay is the primary xenograft model that utilizes immunocompromised mice. The protocol includes injection of primary patient tumor specimens into mice with subsequent serial passaging of the tumors by retransplants of cells harvested from the mouse bone marrow and spleen. The tumors generated are then used for genomic profiling, ex vivo compound testing, mechanistic studies and retransplantation. Detailed in this unit are procedures for the establishment and maintenance of primary ALL xenograft panels for use in basic research and translational studies.

NALP3 inflammasome upregulation and CASP1 cleavage of the glucocorticoid receptor cause glucocorticoid resistance in leukemia cells.

Glucocorticoids are universally used in the treatment of acute lymphoblastic leukemia (ALL), and resistance to glucocorticoids in leukemia cells confers poor prognosis. To elucidate mechanisms of glucocorticoid resistance, we determined the prednisolone sensitivity of primary leukemia cells from 444 patients newly diagnosed with ALL and found significantly higher expression of CASP1 (encoding caspase 1) and its activator NLRP3 in glucocorticoid-resistant leukemia cells, resulting from significantly lower somatic methylation of the CASP1 and NLRP3 promoters. Overexpression of CASP1 resulted in cleavage of the glucocorticoid receptor, diminished the glucocorticoid-induced transcriptional response and increased glucocorticoid resistance. Knockdown or inhibition of CASP1 significantly increased glucocorticoid receptor levels and mitigated glucocorticoid resistance in CASP1-overexpressing ALL. Our findings establish a new mechanism by which the NLRP3-CASP1 inflammasome modulates cellular levels of the glucocorticoid receptor and diminishes cell sensitivity to glucocorticoids. The broad impact on the glucocorticoid transcriptional response suggests that this mechanism could also modify glucocorticoid effects in other diseases.

The landscape of somatic mutations in infant MLL-rearranged acute lymphoblastic leukemias.

Infant acute lymphoblastic leukemia (ALL) with MLL rearrangements (MLL-R) represents a distinct leukemia with a poor prognosis. To define its mutational landscape, we performed whole-genome, exome, RNA and targeted DNA sequencing on 65 infants (47 MLL-R and 18 non-MLL-R cases) and 20 older children (MLL-R cases) with leukemia. Our data show that infant MLL-R ALL has one of the lowest frequencies of somatic mutations of any sequenced cancer, with the predominant leukemic clone carrying a mean of 1.3 non-silent mutations. Despite this paucity of mutations, we detected activating mutations in kinase-PI3K-RAS signaling pathway components in 47% of cases. Surprisingly, these mutations were often subclonal and were frequently lost at relapse. In contrast to infant cases, MLL-R leukemia in older children had more somatic mutations (mean of 6.5 mutations/case versus 1.3 mutations/case, P = 7.15 × 10(-5)) and had frequent mutations (45%) in epigenetic regulators, a category of genes that, with the exception of MLL, was rarely mutated in infant MLL-R ALL.