Whole genome sequencing provides comprehensive genetic testing in childhood B-cell acute lymphoblastic leukaemia

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Benefits for children with suspected cancer from routine whole-genome sequencing

Clinical whole-genome sequencing (WGS) has been shown to deliver potential benefits to children with cancer and to alter treatment in high-risk patient groups. It remains unknown whether offering WGS to every child with suspected cancer can change patient management. We collected WGS variant calls and clinical and diagnostic information from 281 children (282 tumors) across two English units (n = 152 from a hematology center, n = 130 from a solid tumor center) where WGS had become a routine test. Our key finding was that variants uniquely attributable to WGS changed the management in ~7% (20 out of 282) of cases while providing additional disease-relevant findings, beyond standard-of-care molecular tests, in 108 instances for 83 (29%) cases. Furthermore, WGS faithfully reproduced every standard-of-care molecular test (n = 738) and revealed several previously unknown genomic features of childhood tumors. We show that WGS can be delivered as part of routine clinical care to children with suspected cancer and can change clinical management by delivering unexpected genomic insights. Our experience portrays WGS as a clinically impactful assay for routine practice, providing opportunities for assay consolidation and for delivery of molecularly informed patient care.

Acute lymphoblastic leukaemia

Acute lymphoblastic leukaemia (ALL) is a haematological malignancy characterized by the uncontrolled proliferation of immature lymphoid cells. Over past decades, significant progress has been made in understanding the biology of ALL, resulting in remarkable improvements in its diagnosis, treatment and monitoring. Since the advent of chemotherapy, ALL has been the platform to test for innovative approaches applicable to cancer in general. For example, the advent of omics medicine has led to a deeper understanding of the molecular and genetic features that underpin ALL. Innovations in genomic profiling techniques have identified specific genetic alterations and mutations that drive ALL, inspiring new therapies. Targeted agents, such as tyrosine kinase inhibitors and immunotherapies, have shown promising results in subgroups of patients while minimizing adverse effects. Furthermore, the development of chimeric antigen receptor T cell therapy represents a breakthrough in ALL treatment, resulting in remarkable responses and potential long-term remissions. Advances are not limited to treatment modalities alone. Measurable residual disease monitoring and ex vivo drug response profiling screening have provided earlier detection of disease relapse and identification of exceptional responders, enabling clinicians to adjust treatment strategies for individual patients. Decades of supportive and prophylactic care have improved the management of treatment-related complications, enhancing the quality of life for patients with ALL.

MD-ALL: an integrative platform for molecular diagnosis of B-acute lymphoblastic leukemia

B-acute lymphoblastic leukemia (B-ALL) consists of dozens of subtypes defined by distinct gene expression profiles (GEP) and various genetic lesions. With the application of transcriptome sequencing (RNA sequencing [RNA-seq]), multiple novel subtypes have been identified, which lead to an advanced B-ALL classification and risk-stratification system. However, the complexity of analyzing RNA-seq data for B-ALL classification hinders the implementation of the new B-ALL taxonomy. Here, we introduce Molecular Diagnosis of Acute Lymphoblastic Leukemia (MD-ALL), an integrative platform featuring sensitive and accurate B-ALL classification based on GEP and sentinel genetic alterations from RNA-seq data. In this study, we systematically analyzed 2,955 B-ALL RNA-seq samples and generated a reference dataset representing all the reported B-ALL subtypes. Using multiple machine learning algorithms, we identified the feature genes and then established highly sensitive and accurate models for B-ALL classification using either bulk or single-cell RNA-seq data. Importantly, this platform integrates multiple aspects of key genetic lesions acquired from RNA-seq data, which include sequence mutations, large-scale copy number variations, and gene rearrangements, to perform comprehensive and definitive B-ALL classification. Through validation in a hold-out cohort of 974 samples, our models demonstrated superior performance for B-ALL classification compared with alternative tools. Moreover, to ensure accessibility and user-friendly navigation even for users with limited or no programming background, we developed an interactive graphical user interface for this MD-ALL platform, using the R Shiny package. In summary, MD-ALL is a user-friendly B-ALL classification platform designed to enable integrative, accurate, and comprehensive B-ALL subtype classification. MD-ALL is available from https://github.com/gu-lab20/MD-ALL.

Cas9-directed long-read sequencing to resolve optical genome mapping findings in leukemia diagnostics

Leukemias are genetically heterogeneous and diagnostics therefore includes various standard-of-care (SOC) techniques, including karyotyping, SNP-array and FISH. Optical genome mapping (OGM) may replace these as it detects different types of structural aberrations simultaneously and additionally detects much smaller aberrations (500 bp vs 5-10 Mb with karyotyping). However, its resolution may still be too low to define clinical relevance of aberrations when they are located between two OGM labels or when labels are not distinct enough. Here, we test the potential of Cas9-directed long-read sequencing (LRS) as an additional technique to resolve such potentially relevant new findings. From an internal Bionano implementation study we selected ten OGM calls that could not be validated with SOC methods. Per variant we designed crRNAs for Cas9 enrichment, prepared libraries and sequenced them on a MinION/GridION device. We could confirm all aberrations and, importantly, the actual breakpoints of the OGM calls were located between 0.2 and 5.5 kb of the OGM-estimated breakpoints, confirming the high reliability of OGM. Furthermore, we show examples of redefinition of aberrations between labels that enable judgment of clinical relevance. Our results suggest that Cas9-directed LRS can be a relevant and flexible secondary technique in diagnostic workflows including OGM.

Refining risk prediction in pediatric acute lymphoblastic leukemia through DNA methylation profiling

Acute lymphoblastic leukemia (ALL) is the most prevalent cancer in children, and despite considerable progress in treatment outcomes, relapses still pose significant risks of mortality and long-term complications. To address this challenge, we employed a supervised machine learning technique, specifically random survival forests, to predict the risk of relapse and mortality using array-based DNA methylation data from a cohort of 763 pediatric ALL patients treated in Nordic countries. The relapse risk predictor (RRP) was constructed based on 16 CpG sites, demonstrating c-indexes of 0.667 and 0.677 in the training and test sets, respectively. The mortality risk predictor (MRP), comprising 53 CpG sites, exhibited c-indexes of 0.751 and 0.754 in the training and test sets, respectively. To validate the prognostic value of the predictors, we further analyzed two independent cohorts of Canadian (n = 42) and Nordic (n = 384) ALL patients. The external validation confirmed our findings, with the RRP achieving a c-index of 0.667 in the Canadian cohort, and the RRP and MRP achieving c-indexes of 0.529 and 0.621, respectively, in an independent Nordic cohort. The precision of the RRP and MRP models improved when incorporating traditional risk group data, underscoring the potential for synergistic integration of clinical prognostic factors. The MRP model also enabled the definition of a risk group with high rates of relapse and mortality. Our results demonstrate the potential of DNA methylation as a prognostic factor and a tool to refine risk stratification in pediatric ALL. This may lead to personalized treatment strategies based on epigenetic profiling.

RNA-seq-based miRNA signature as an independent predictor of relapse in pediatric B-cell acute lymphoblastic leukemia

ABSTRACT

Aberrant micro-RNA (miRNA) expression profiles have been associated with disease progression and clinical outcome in pediatric cancers. However, few studies have analyzed genome-wide dysregulation of miRNAs and messenger RNAs (mRNAs) in pediatric B-cell precursor acute lymphoblastic leukemia (BCP-ALL). To identify novel prognostic factors, we comprehensively investigated miRNA and mRNA sequencing (miRNA-seq and mRNA-seq) data in pediatric BCP-ALL samples with poor outcome. We analyzed 180 patients, including 43 matched pairs at diagnosis and relapse. Consensus clustering of miRNA expression data revealed a distinct profile characterized by mainly downregulation of miRNAs (referred to as an miR-low cluster [MLC]). The MLC profile was not associated with any known genetic subgroups. Intriguingly, patients classified as MLC had significantly shorter event-free survival (median 21 vs 33 months; log-rank P = 3 ×10-5). Furthermore, this poor prognosis was retained even in hyperdiploid ALL. This poor prognostic MLC profiling was confirmed in the validation cohort. Notably, non-MLC profiling at diagnosis (n = 9 of 23; Fisher exact test, P = .039) often changed into MLC profiling at relapse for the same patient. Integrated analysis of miRNA-seq and mRNA-seq data revealed that the transcriptional profile of MLC was characterized by enrichment of MYC target and oxidative phosphorylation genes, reduced intron retention, and low expression of DICER1. Thus, our miRNA-mRNA integration approach yielded a truly unbiased molecular stratification of pediatric BCP-ALL cases based on a novel prognostic miRNA signature, which may lead to better clinical outcomes.

The role of next-generation sequencing in hematologic malignancies

Next-generation sequencing (NGS) allows high-throughput detection of molecular changes in tumors. Over the past 15 years, NGS has rapidly evolved from a promising research tool to a core component of the clinical laboratory. Sequencing of tumor cells provides an important step in detecting somatic driver mutations that not only characterize the disease but also influence treatment decisions. For patients with hematologic malignancies, NGS has been used for accurate classification and diagnosis based on genetic alterations. The recently revised World Health Organization classification and the European LeukemiaNet recommendations for acute myeloid leukemia consider genetic abnormalities as a top priority for diagnosis, prognostication, monitoring of measurable residual disease, and treatment choice. This review aims to present the role and utility of various NGS approaches for the diagnosis, treatment, and follow-up of hemato-oncology patients.

The landscape of cytogenetic and molecular genetic methods in diagnostics for hematologic neoplasia

Improvements made during the last decades in the management of patients with hematologic neoplasia have resulted in increase of overall survival. These advancements have become possible through progress in our understanding of genetic basis of different hematologic malignancies and their role in the current risk-adapted treatment protocols. In this review, we provide an overview of current cytogenetic and molecular genetic methods, commonly used in the genetic characterization of hematologic malignancies, describe the current developments in the cytogenetic and molecular diagnostics, and give an outlook into their future development. Furthermore, we give a brief overview of the most important public databases and guidelines for sequence variant interpretation.

Genomic technologies for detecting structural variations in hematologic malignancies

Genomic structural variations in myeloid, lymphoid, and plasma cell neoplasms can provide key diagnostic, prognostic, and therapeutic information while elucidating the underlying disease biology. Several molecular diagnostic approaches play a central role in evaluating hematological malignancies. Traditional cytogenetic diagnostic assays, such as chromosome banding and fluorescence in situ hybridization, are essential components of the current diagnostic workup that guide clinical care for most hematologic malignancies. However, each assay has inherent limitations, including limited resolution for detecting small structural variations and low coverage, and can only detect alterations in the target regions. Recently, the rapid expansion and increasing availability of novel and comprehensive genomic technologies have led to their use in clinical laboratories for clinical management and translational research. This review aims to describe the clinical relevance of structural variations in hematologic malignancies and introduce genomic technologies that may facilitate personalized tumor characterization and treatment.

IKZF1 Alterations and Therapeutic Targeting in B-Cell Acute Lymphoblastic Leukemia

IKZF1 encodes the transcription factor IKAROS, a zinc finger DNA-binding protein with a key role in lymphoid lineage development. IKAROS plays a critical role in the development of lineage-restricted mature lymphocytes. Deletions within IKZF1 in B-cell acute lymphoblastic leukemia (B-ALL) lead to a loss of normal IKAROS function, conferring leukemic stem cell properties, including self-renewal and subsequent uncontrolled growth. IKZF1 deletions are associated with treatment resistance and inferior outcomes. Early identification of IKZF1 deletions in B-ALL may inform the intensification of therapy and other potential treatment strategies to improve outcomes in this high-risk leukemia.

5-Aza-4'-thio-2'-deoxycytidine induces C>G transversions in a specific trinucleotide context and leads to acute lymphoid leukemia

DNA methyltransferase inhibitors (DNMTi), most commonly cytidine analogs, are compounds that are used clinically to decrease 5'-cytosine methylation, with the aim of re-expression of tumor suppressor genes. We used a murine pre-clinical model of myelodysplastic syndrome based on transplantation of cells expressing a NUP98::HOXD13 transgene to investigate 5-Aza-4'-thio-2'-deoxycytidine (Aza TdCyd or ATC), a thiol substituted DNMTi, as a potential therapy. We found that ATC treatment led to lymphoid leukemia in wild-type recipient cells; further study revealed that healthy mice treated with ATC also developed lymphoid leukemia. Whole exome sequencing revealed thousands of acquired mutations, almost all of which were C > G transversions in a previously unrecognized, specific 5'-NCG-3' context. These mutations involved dozens of genes well-known to be involved in human lymphoid leukemia, such as Notch1, Pten, Pax5, Trp53 , and Nf1 . Treatment of human cells in vitro showed thousands of acquired C > G transversions in a similar context. Deletion of Dck , the rate-limiting enzyme for the cytidine salvage pathway, eliminated C > G transversions. Taken together, these findings demonstrate that DNMTi can be potent mutagens in human and mouse cells, both in vitro and in vivo .

Cytogenetics in the management of B-cell acute lymphoblastic leukemia: Guidelines from the Groupe Francophone de Cytogénétique Hématologique (GFCH)

Cytogenetic analysis is mandatory at initial assessment of B-cell acute lymphoblastic leukemia (B-ALL) due to its diagnostic and prognostic value. Results from chromosome banding analysis and complementary FISH are taken into account in therapeutic protocols and further completed by other techniques (RT-PCR, SNP-array, MLPA, NGS, OGM). Indeed, new genomic entities have been identified by NGS, mostly RNA sequencing, such as Ph-like ALL that can benefit from targeted therapy. Here, we have attempted to establish cytogenetic guidelines by reviewing the most recent published data including the novel 5th World Health Organization and International Consensus Classifications. We also focused on newly described cytogenomic entities and indicate alternative diagnostic tools such as NGS technology, as its importance is vastly increasing in the diagnostic setting.

Reproducible Bioinformatics Analysis Workflows for Detecting IGH Gene Fusions in B-Cell Acute Lymphoblastic Leukaemia Patients

B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by diverse genomic alterations, the most frequent being gene fusions detected via transcriptomic analysis (mRNA-seq). Due to its hypervariable nature, gene fusions involving the Immunoglobulin Heavy Chain (IGH) locus can be difficult to detect with standard gene fusion calling algorithms and significant computational resources and analysis times are required. We aimed to optimize a gene fusion calling workflow to achieve best-case sensitivity for IGH gene fusion detection. Using Nextflow, we developed a simplified workflow containing the algorithms FusionCatcher, Arriba, and STAR-Fusion. We analysed samples from 35 patients harbouring IGH fusions (IGH::CRLF2 n = 17, IGH::DUX4 n = 15, IGH::EPOR n = 3) and assessed the detection rates for each caller, before optimizing the parameters to enhance sensitivity for IGH fusions. Initial results showed that FusionCatcher and Arriba outperformed STAR-Fusion (85-89% vs. 29% of IGH fusions reported). We found that extensive filtering in STAR-Fusion hindered IGH reporting. By adjusting specific filtering steps (e.g., read support, fusion fragments per million total reads), we achieved a 94% reporting rate for IGH fusions with STAR-Fusion. This analysis highlights the importance of filtering optimization for IGH gene fusion events, offering alternative workflows for difficult-to-detect high-risk B-ALL subtypes.

The genomic landscape of acute lymphoblastic leukemia with intrachromosomal amplification of chromosome 21

Intrachromosomal amplification of chromosome 21 defines a subtype of high-risk childhood acute lymphoblastic leukemia (iAMP21-ALL) characterized by copy number changes and complex rearrangements of chromosome 21. The genomic basis of iAMP21-ALL and the pathogenic role of the region of amplification of chromosome 21 to leukemogenesis remains incompletely understood. In this study, using integrated whole genome and transcriptome sequencing of 124 patients with iAMP21-ALL, including rare cases arising in the context of constitutional chromosomal aberrations, we identified subgroups of iAMP21-ALL based on the patterns of copy number alteration and structural variation. This large data set enabled formal delineation of a 7.8 Mb common region of amplification harboring 71 genes, 43 of which were differentially expressed compared with non-iAMP21-ALL ones, including multiple genes implicated in the pathogenesis of acute leukemia (CHAF1B, DYRK1A, ERG, HMGN1, and RUNX1). Using multimodal single-cell genomic profiling, including single-cell whole genome sequencing of 2 cases, we documented clonal heterogeneity and genomic evolution, demonstrating that the acquisition of the iAMP21 chromosome is an early event that may undergo progressive amplification during disease ontogeny. We show that UV-mutational signatures and high mutation load are characteristic secondary genetic features. Although the genomic alterations of chromosome 21 are variable, these integrated genomic analyses and demonstration of an extended common minimal region of amplification broaden the definition of iAMP21-ALL for more precise diagnosis using cytogenetic or genomic methods to inform clinical management.

Diagnostic utility of whole genome sequencing in adults with B-other acute lymphoblastic leukemia

Genomic profiling during the diagnosis of B-cell precursor acute lymphoblastic leukemia (BCP-ALL) in adults is used to guide disease classification, risk stratification, and treatment decisions. Patients for whom diagnostic screening fails to identify disease-defining or risk-stratifying lesions are classified as having B-other ALL. We screened a cohort of 652 BCP-ALL cases enrolled in UKALL14 to identify and perform whole genome sequencing (WGS) of paired tumor-normal samples. For 52 patients with B-other, we compared the WGS findings with data from clinical and research cytogenetics. WGS identified a cancer-associated event in 51 of 52 patients, including an established subtype defining genetic alterations that were previously missed with standard-of-care (SoC) genetics in 5 of them. Of the 47 true B-other ALL, we identified a recurrent driver in 87% (41). A complex karyotype via cytogenetics emerges as a heterogeneous group, including distinct genetic alterations associated with either favorable (DUX4-r) or poor outcomes (MEF2D-r and IGK::BCL2). For a subset of 31 cases, we integrated the findings from RNA sequencing (RNA-seq) analysis to include fusion gene detection and classification based on gene expression. Compared with RNA-seq, WGS was sufficient to detect and resolve recurrent genetic subtypes; however, RNA-seq can provide orthogonal validation of findings. In conclusion, we demonstrated that WGS can identify clinically relevant genetic abnormalities missed with SoC testing as well as identify leukemia driver events in virtually all cases of B-other ALL.

Multi-Omic Approaches to Classify, Predict, and Treat Acute Leukemias

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer, in which nearly 5% of the cases are diagnosed before the first year of age [...].