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.

Paediatric B lymphoblastic leukaemia with hyperdiploidy and a false-positive KMT2A fluorescence in situ hybridization result.

The dramatic improvement in the event-free survival of paediatric B-lymphoblastic leukaemia (B-ALL) has led to risk-stratified treatment. Through a combination of clinical features, cytogenetic abnormalities and assessment of treatment response, patients are stratified to receive different intensities of therapy. The presence of high hyperdiploidy (>50 chromosomes) is considered a favourable genetic feature. Conversely, KMT2A fusion genes in B-ALL are associated with a poor prognosis, resulting in intensification of treatment. We present a seven-year-old female with B-ALL, a high hyperdiploid karyotype (56 chromosomes) and KMT2A rearrangement detected on FISH, but with no productive fusion identified. Single nucleotide polymorphism (SNP) array suggested the KMT2A rearrangement was due to chromosome 11 chromothripsis. Subsequent targeted RNA fusion panel and whole transcriptomic sequencing (mRNA-seq) did not detect an expressed KMT2A fusion. Differential expression analyses of the mRNA-seq data led to clustering of this case with other hyperdiploid cases, consistent with the hyperdiploid cytogenetic results. Given the additional intensity and potential toxicity of high-risk treatment, unusual findings by chromosome analysis, FISH and/or chromosomal microarray should prompt consideration of testing for a KMT2A fusion by another method to avoid misclassification.

Case report: Rare case of donor cell-derived T-cell acute lymphoblastic leukaemia in a female patient after receiving an allo-transplant from her male sibling.

Donor-derived haematological neoplasms, in which recipients present with haematological malignancies that have evolved from transplant donor stem cells, have previously been described for myelodysplastic syndrome, myeloproliferative neoplasms, acute myeloid leukaemia and less often, leukaemias of lymphoid origin. Here we describe a rare and complex case of donor-derived T-cell acute lymphoblastic leukaemia with a relatively short disease latency of less than 4 years. Through genomic and in vitro analyses, we identified novel mutations in NOTCH1 as well as a novel activating mutation in STAT5B; the latter targetable with the clinically available drugs, venetoclax and ruxolitinib.

Gain of chromosome 21 increases the propensity for P2RY8::CRLF2 acute lymphoblastic leukemia via increased HMGN1 expression.

Acute lymphoblastic leukemia (ALL) patients with a gain of chromosome 21, intrachromosomal amplification of chromosome 21 (iAMP21), or Down syndrome (DS), have increased expression of genes in the DS critical region (DSCR) of chromosome 21, including the high-mobility group nucleosome-binding protein 1, HMGN1. Children with DS are predisposed to develop hematologic malignancies, providing insight into the role of chromosome 21 in the development of leukemias. A 320-kb deletion in the pseudoautosomal region of the X/Y chromosome in leukemic cells, resulting in a gene fusion between the purinergic receptor and cytokine receptor-like factor-2 (P2Y Receptor Family Member 8 (P2RY8)::CRLF2), is a common feature in ~60% of DS-ALL and ~40% of iAMP21 patients, suggesting a link between chromosome 21 and P2RY8::CRLF2. In an Australian cohort of pediatric B-ALL patients with P2RY8::CRLF2 (n = 38), eight patients harbored gain of chromosome 21 (+21), and two patients had iAMP21, resulting in a significantly increased HMGN1 expression. An inducible CRISPR/Cas9 system was used to model P2RY8::CRLF2 and investigate its cooperation with HMGN1. This model was then used to validate HMGN1 as an influencing factor for P2RY8::CRLF2 development. Using Cas9 to cleave the DNA at the pseudoautosomal region without directed repair, cells expressing HMGN1 favored repair, resulting in P2RY8::CRLF2 generation, compared with cells without HMGN1. CRISPR/Cas9 P2RY8::CRLF2 cells expressing HMGN1 exhibit increased proliferation, thymic stromal lymphopoietin receptor (TSLPR) expression, and JAK/STAT signaling, consistent with cells from patients with P2RY8::CRLF2. Our patient expression data and unique CRISPR/Cas9 modeling, when taken together, suggest that HMGN1 increases the propensity for P2RY8::CRLF2 development. This has important implications for patients with DS, +21, or iAMP21.

Case Report: Rare IKZF1 Gene Fusions Identified in Neonate with Congenital KMT2A-Rearranged Acute Lymphoblastic Leukemia.

Chromosomal rearrangements involving the KMT2A gene occur frequently in acute lymphoblastic leukaemia (ALL). KMT2A-rearranged ALL (KMT2Ar ALL) has poor long-term survival rates and is the most common ALL subtype in infants less than 1 year of age. KMT2Ar ALL frequently occurs with additional chromosomal abnormalities including disruption of the IKZF1 gene, usually by exon deletion. Typically, KMT2Ar ALL in infants is accompanied by a limited number of cooperative le-sions. Here we report a case of aggressive infant KMT2Ar ALL harbouring additional rare IKZF1 gene fusions. Comprehensive genomic and transcriptomic analyses were performed on sequential samples. This report highlights the genomic complexity of this particular disease and describes the novel gene fusions IKZF1::TUT1 and KDM2A::IKZF1.

RaScALL: Rapid (Ra) screening (Sc) of RNA-seq data for prognostically significant genomic alterations in acute lymphoblastic leukaemia (ALL).

RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, integrating RNA-seq in a clinical setting requires rapid detection and accurate reporting of clinically relevant alterations. Here we present RaScALL, an implementation of the k-mer based variant detection tool km, capable of identifying more than 100 prognostically significant lesions observed in ALL, including gene fusions, single nucleotide variants and focal gene deletions. We compared genomic alterations detected by RaScALL and those reported by alignment-based de novo variant detection tools in a study cohort of 180 Australian patient samples. Results were validated using 100 patient samples from a published North American cohort. RaScALL demonstrated a high degree of accuracy for reporting subtype defining genomic alterations. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% of reported cases in the study cohort (n = 164) and 95% of samples (n = 63) in the validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n = 12) and IKZF1 p.N159Y (n = 4). Intragenic IKZF1 deletions resulting in aberrant transcript isoforms were also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches. The application of RaScALL enables rapid identification and reporting of previously identified genomic alterations of known clinical relevance.

Multi-Cohort Transcriptomic Subtyping of B-Cell Acute Lymphoblastic Leukemia.

RNA sequencing provides a snapshot of the functional consequences of genomic lesions that drive acute lymphoblastic leukemia (ALL). The aims of this study were to elucidate diagnostic associations (via machine learning) between mRNA-seq profiles, independently verify ALL lesions and develop easy-to-interpret transcriptome-wide biomarkers for ALL subtyping in the clinical setting. A training dataset of 1279 ALL patients from six North American cohorts was used for developing machine learning models. Results were validated in 767 patients from Australia with a quality control dataset across 31 tissues from 1160 non-ALL donors. A novel batch correction method was introduced and applied to adjust for cohort differences. Out of 18,503 genes with usable expression, 11,830 (64%) were confounded by cohort effects and excluded. Six ALL subtypes (ETV6::RUNX1, KMT2A, DUX4, PAX5 P80R, TCF3::PBX1, ZNF384) that covered 32% of patients were robustly detected by mRNA-seq (positive predictive value ≥ 87%). Five other frequent subtypes (CRLF2, hypodiploid, hyperdiploid, PAX5 alterations and Ph-positive) were distinguishable in 40% of patients at lower accuracy (52% ≤ positive predictive value ≤ 73%). Based on these findings, we introduce the Allspice R package to predict ALL subtypes and driver genes from unadjusted mRNA-seq read counts as encountered in real-world settings. Two examples of Allspice applied to previously unseen ALL patient samples with atypical lesions are included.

Case Report: Precision Medicine Target Revealed by In Vitro Modeling of Relapsed, Refractory Acute Lymphoblastic Leukemia From a Child With Neurofibromatosis.

Children with neurofibromatosis have a higher risk of developing juvenile myelomonocytic leukemia and acute myeloid leukemia, but rarely develop B-cell acute lymphoblastic leukemia (B-ALL). Through in-vitro modeling, a novel NF1 p.L2467 frameshift (fs) mutation identified in a relapsed/refractory Ph-like B-ALL patient with neurofibromatosis demonstrated cytokine independence and increased RAS signaling, indicative of leukemic transformation. Furthermore, these cells were sensitive to the MEK inhibitors trametinib and mirdametinib. Bi-allelic NF1 loss of function may be a contributing factor to relapse and with sensitivity to MEK inhibitors, suggests a novel precision medicine target in the setting of neurofibromatosis patients with B-ALL.

Enhancer retargeting of CDX2 and UBTF::ATXN7L3 define a subtype of high-risk B-progenitor acute lymphoblastic leukemia.

Transcriptome sequencing has identified multiple subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL) of prognostic significance, but a minority of cases lack a known genetic driver. Here, we used integrated whole-genome (WGS) and -transcriptome sequencing (RNA-seq), enhancer mapping, and chromatin topology analysis to identify previously unrecognized genomic drivers in B-ALL. Newly diagnosed (n = 3221) and relapsed (n = 177) B-ALL cases with tumor RNA-seq were studied. WGS was performed to detect mutations, structural variants, and copy number alterations. Integrated analysis of histone 3 lysine 27 acetylation and chromatin looping was performed using HiChIP. We identified a subset of 17 newly diagnosed and 5 relapsed B-ALL cases with a distinct gene expression profile and 2 universal and unique genomic alterations resulting from aberrant recombination-activating gene activation: a focal deletion downstream of PAN3 at 13q12.2 resulting in CDX2 deregulation by the PAN3 enhancer and a focal deletion of exons 18-21 of UBTF at 17q21.31 resulting in a chimeric fusion, UBTF::ATXN7L3. A subset of cases also had rearrangement and increased expression of the PAX5 gene, which is otherwise uncommon in B-ALL. Patients were more commonly female and young adult with median age 35 (range,12-70 years). The immunophenotype was characterized by CD10 negativity and immunoglobulin M positivity. Among 16 patients with known clinical response, 9 (56.3%) had high-risk features including relapse (n = 4) or minimal residual disease >1% at the end of remission induction (n = 5). CDX2-deregulated, UBTF::ATXN7L3 rearranged (CDX2/UBTF) B-ALL is a high-risk subtype of leukemia in young adults for which novel therapeutic approaches are required.

Identification of a novel GOLGA4-JAK2 fusion gene in B-cell acute lymphoblastic leukaemia.

Rearrangements of Janus kinase 2 (JAK2r) form a subtype of acute lymphoblastic leukaemia (ALL) associated with poor patient outcomes. We present a high-risk case of B-cell ALL (B-ALL) where retrospective mRNA sequencing identified a novel GOLGA4-JAK2 fusion gene. Expression of GOLGA4-JAK2 in murine pro-B cells promoted factor-independent growth, implicating GOLGA4-JAK2 as an oncogenic driver. Cells expressing GOLGA4-JAK2 demonstrated constitutive activation of JAK/STAT signalling and were sensitive to JAK inhibitors. This study contributes to the diverse collection of JAK2 fusion genes identified in B-ALL and supports the incorporation of JAK inhibitors into treatment strategies to improve outcomes for this subtype.

Acquired JAK2 mutations confer resistance to JAK inhibitors in cell models of acute lymphoblastic leukemia.

Ruxolitinib (rux) Phase II clinical trials are underway for the treatment of high-risk JAK2-rearranged (JAK2r) B-cell acute lymphoblastic leukemia (B-ALL). Treatment resistance to targeted inhibitors in other settings is common; elucidating potential mechanisms of rux resistance in JAK2r B-ALL will enable development of therapeutic strategies to overcome or avert resistance. We generated a murine pro-B cell model of ATF7IP-JAK2 with acquired resistance to multiple type-I JAK inhibitors. Resistance was associated with mutations within the JAK2 ATP/rux binding site, including a JAK2 p.G993A mutation. Using in vitro models of JAK2r B-ALL, JAK2 p.G993A conferred resistance to six type-I JAK inhibitors and the type-II JAK inhibitor, CHZ-868. Using computational modeling, we postulate that JAK2 p.G993A enabled JAK2 activation in the presence of drug binding through a unique resistance mechanism that modulates the mobility of the conserved JAK2 activation loop. This study highlights the importance of monitoring mutation emergence and may inform future drug design and the development of therapeutic strategies for this high-risk patient cohort.

DUX Hunting-Clinical Features and Diagnostic Challenges Associated with DUX4-Rearranged Leukaemia.

DUX4-rearrangement (DUX4r) is a recently discovered recurrent genomic lesion reported in 4-7% of childhood B cell acute lymphoblastic leukaemia (B-ALL) cases. This subtype has favourable outcomes, especially in children and adolescents treated with intensive chemotherapy. The fusion most commonly links the hypervariable IGH gene to DUX4 a gene located within the D4Z4 macrosatellite repeat on chromosome 4, with a homologous polymorphic repeat on chromosome 10. DUX4r is cryptic to most standard diagnostic techniques, and difficult to identify even with next generation sequencing assays. This review summarises the clinical features and molecular genetics of DUX4r B-ALL and proposes prospective new diagnostic methods.