Automated identification of cell populations in flow cytometry data with transformers.

Acute Lymphoblastic Leukemia (ALL) is the most frequent hematologic malignancy in children and adolescents. A strong prognostic factor in ALL is given by the Minimal Residual Disease (MRD), which is a measure for the number of leukemic cells persistent in a patient. Manual MRD assessment from Multiparameter Flow Cytometry (FCM) data after treatment is time-consuming and subjective. In this work, we present an automated method to compute the MRD value directly from FCM data. We present a novel neural network approach based on the transformer architecture that learns to directly identify blast cells in a sample. We train our method in a supervised manner and evaluate it on publicly available ALL FCM data from three different clinical centers. Our method reaches a median F1 score of ≈0.94 when evaluated on 519 B-ALL samples and shows better results than existing methods on 4 different datasets.

An Extensive Quality Control and Quality Assurance (QC/QA) Program Significantly Improves Inter-Laboratory Concordance Rates of Flow-Cytometric Minimal Residual Disease Assessment in Acute Lymphoblastic Leukemia: An I-BFM-FLOW-Network Report.

Monitoring of minimal residual disease (MRD) by flow cytometry (FCM) is a powerful prognostic tool for predicting outcomes in acute lymphoblastic leukemia (ALL). To apply FCM-MRD in large, collaborative trials, dedicated laboratory staff must be educated to concordantly high levels of expertise and their performance quality should be continuously monitored. We sought to install a unique and comprehensive training and quality control (QC) program involving a large number of reference laboratories within the international Berlin-Frankfurt-Münster (I-BFM) consortium, in order to complement the standardization of the methodology with an educational component and persistent quality control measures. Our QC and quality assurance (QA) program is based on four major cornerstones: (i) a twinning maturation program, (ii) obligatory participation in external QA programs (spiked sample send around, United Kingdom National External Quality Assessment Service (UK NEQAS)), (iii) regular participation in list-mode-data (LMD) file ring trials (FCM data file send arounds), and (iv) surveys of independent data derived from trial results. We demonstrate that the training of laboratories using experienced twinning partners, along with continuous educational feedback significantly improves the performance of laboratories in detecting and quantifying MRD in pediatric ALL patients. Overall, our extensive education and quality control program improved inter-laboratory concordance rates of FCM-MRD assessments and ultimately led to a very high conformity of risk estimates in independent patient cohorts.

Risk factors and outcomes in children with high-risk B-cell precursor and T-cell relapsed acute lymphoblastic leukaemia: combined analysis of ALLR3 and ALL-REZ BFM 2002 clinical trials.

AIM: Outcomes of children with high-risk (HR) relapsed acute lymphoblastic leukaemia (ALL) (N = 393), recruited to ALLR3 and ALL-REZ BFM 2002 trials, were analysed. Minimal residual disease (MRD) was assessed after induction and at predetermined time points until haematopoietic stem cell transplantation (SCT). METHODS: Genetic analyses included karyotype, copy-number alterations and mutation analyses. Ten-year survivals were analysed using Kaplan-Meier and Cox models for multivariable analyses. RESULTS: Outcomes of patients were comparable in ALLR3 and ALL-REZ BFM 2002. The event-free survival of B-cell precursor (BCP) and T-cell ALL (T-ALL) was 22.6% and 26.2% (P = 0.94), respectively, and the overall survival (OS) was 32.6% and 28.2% (P = 0.11), respectively. Induction failures (38%) were associated with deletions of NR3C1 (P = 0.002) and BTG1 (P = 0.03) in BCP-ALL. The disease-free survival (DFS) and OS in patients with good vs poor MRD responses were 57.4% vs 22.6% (P < 0.0001) and 57.8% vs 32.0% (P = 0.0004), respectively. For BCP- and T-ALL, the post-SCT DFS and OS were 42.1% and 56.8% (P = 0.26) and 51.6% and 55.4% (P = 0.67), respectively. The cumulative incidences of post-SCT relapse for BCP- and T-ALL were 36.9% and 17.8% (P = 0.012) and of death were 10.7% and 25.5% (P = 0.013), respectively. Determinants of outcomes after SCT were acute graft versus host disease, pre-SCT MRD (≥10-3), HR cytogenetics and TP53 alterations in BCP-ALL. CONCLUSION: Improvements in outcomes for HR ALL relapses require novel compounds in induction therapy to improve remission rates and immune targeted therapy after induction to maintain remission after SCT. TRIAL REGISTRATION: ALLR3: NCT00967057; ALL REZ-BFM 2002: NCT00114348.

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

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

Reflection of neuroblastoma intratumor heterogeneity in the new OHC-NB1 disease model.

Accurate modeling of intratumor heterogeneity presents a bottleneck against drug testing. Flexibility in a preclinical platform is also desirable to support assessment of different endpoints. We established the model system, OHC-NB1, from a bone marrow metastasis from a patient diagnosed with MYCN-amplified neuroblastoma and performed whole-exome sequencing on the source metastasis and the different models and passages during model development (monolayer cell line, 3D spheroid culture and subcutaneous xenograft tumors propagated in mice). OHC-NB1 harbors a MYCN amplification in double minutes, 1p deletion, 17q gain and diploid karyotype, which persisted in all models. A total of 80-540 single-nucleotide variants (SNVs) was detected in each sample, and comparisons between the source metastasis and models identified 34 of 80 somatic SNVs to be propagated in the models. Clonal reconstruction using the combined copy number and SNV data revealed marked clonal heterogeneity in the originating metastasis, with four clones being reflected in the model systems. The set of OHC-NB1 models represents 43% of somatic SNVs and 23% of the cellularity in the originating metastasis with varying clonal compositions, indicating that heterogeneity is partially preserved in our model system.

Improving Stratification for Children With Late Bone Marrow B-Cell Acute Lymphoblastic Leukemia Relapses With Refined Response Classification and Integration of Genetics.

PURPOSE: Minimal residual disease (MRD) helps to accurately assess when children with late bone marrow relapses of B-cell precursor (BCP) acute lymphoblastic leukemia (ALL) will benefit from allogeneic hematopoietic stem-cell transplantation (allo-HSCT). More detailed dissection of MRD response heterogeneity and the specific genetic aberrations could improve current practice. PATIENTS AND METHODS: MRD was assessed after induction treatment and at different times during relapse treatment until allo-HSCT (indicated in poor responders to induction; MRD ≥ 10-3) for patients being treated for late BCP-ALL bone marrow relapses (n = 413; median follow-up, 9.4 years) in the ALL-REZ BFM 2002 trial/registry (ClinicalTrials.gov identifier: NCT00114348). RESULTS: Patients with both good (MRD < 10-3) and poor responses to induction treatment reached excellent event-free survival (EFS; 72% v 65%) and overall survival (OS; 82% v 74%). Patients with MRD of 10-2 or greater after induction had reduced EFS (56%), and their MRD persisted until allo-HSCT more frequently than it did in patients with MRD of 10-3 or greater to less than 10-2 (P = .037). Patients with 25% or more leukemic blasts after induction (early nonresponders) had the poorest prognosis (EFS, 22%). Interestingly, patients with MRD of 10-3 or greater before allo-HSCT (late nonresponders) still had an EFS of 50% and OS of 63%, which in principle justifies allo-HSCT in these patients. From a panel of selected candidate genes, TP53 alterations (frequency, 8%) were the only genetic alteration with independent prognostic value in any MRD-based response subgroup. CONCLUSION: After induction treatment, MRD-based treatment stratification resulted in excellent survival in patients with late relapsed BCP-ALL. Prognosis could be further improved in very poor responders by intensifying treatment directly after induction. TP53 alterations can be defined as a novel genetic high-risk marker in all MRD response groups in late relapsed BCP-ALL. Here we identified early and late nonresponders to be considered as events in future trials.

Automated Flow Cytometric MRD Assessment in Childhood Acute B- Lymphoblastic Leukemia Using Supervised Machine Learning.

Minimal residual disease (MRD) as measured by multiparameter flow cytometry (FCM) is an independent and strong prognostic factor in B-cell acute lymphoblastic leukemia (B-ALL). However, reliable flow cytometric detection of MRD strongly depends on operator skills and expert knowledge. Hence, an objective, automated tool for reliable FCM-MRD quantification, able to overcome the technical diversity and analytical subjectivity, would be most helpful. We developed a supervised machine learning approach using a combination of multiple Gaussian Mixture Models (GMM) as a parametric density model. The approach was used for finding the weights of a linear combination of multiple GMMs to represent new, "unseen" samples by an interpolation of stored samples. The experimental data set contained FCM-MRD data of 337 bone marrow samples collected at day 15 of induction therapy in three different laboratories from pediatric patients with B-ALL for which accurate, expert-set gates existed. We compared MRD quantification by our proposed GMM approach to operator assessments, its performance on data from different laboratories, as well as to other state-of-the-art automated read-out methods. Our proposed GMM-combination approach proved superior over support vector machines, deep neural networks, and a single GMM approach in terms of precision and average F 1 -scores. A high correlation of expert operator-based and automated MRD assessment was achieved with reliable automated MRD quantification (F 1 -scores >0.5 in more than 95% of samples) in the clinically relevant range. Although best performance was found, if test and training samples were from the same system (i.e., flow cytometer and staining panel; lowest median F 1 -score 0.92), cross-system performance remained high with a median F 1 -score above 0.85 in all settings. In conclusion, our proposed automated approach could potentially be used to assess FCM-MRD in B-ALL in an objective and standardized manner across different laboratories. © 2019 International Society for Advancement of Cytometry.

Aneuploidy in children with relapsed B-cell precursor acute lymphoblastic leukaemia: clinical importance of detecting a hypodiploid origin of relapse.

Aneuploidy is common in paediatric B-cell precursor acute lymphoblastic leukaemia (ALL). Specific subgroups, such as high hyperdiploidy (>50 chromosomes or DNA Index ≥1·16) and hypodiploidy (<45 chromosomes), predict outcome of patients after primary treatment. Whether aneuploidy has a prognostic value for relapsed disease is yet to be determined. Using DNA index and centromere screening by multiplex ligation-dependent probe amplification, we investigated aneuploidy in 413 children treated for first relapse of B-cell precursor ALL according to the ALL-REZ BFM 2002 protocol. Ten-year event-free survival of patients with high hyperdiploid relapses approached 70%, whereas it was only 40% in low hyperdiploid relapses. Three patients with apparent hyperdiploid relapse had TP53 mutations. In these cases, array-based allelotyping revealed a hypodiploid origin with absence of the hypodiploid founder clone (masked hypodiploidy). Collectively, patients with evident or masked hypodiploid relapses showed an extremely low event-free survival rate of 9%. Importantly, the current relapse risk stratification did not identify cases with masked hypodiploidy as high-risk patients, due to their favourable clinical presentation. In multivariate analysis, hypodiploidy proved to be an independent prognostic factor. This finding supports stratification of relapses with hypodiploid origin into high-risk arms in future trials or allocation of patients to alternative treatment approaches.

IKZF1plus Defines a New Minimal Residual Disease-Dependent Very-Poor Prognostic Profile in Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia.

Purpose Somatic deletions that affect the lymphoid transcription factor-coding gene IKZF1 have previously been reported as independently associated with a poor prognosis in pediatric B-cell precursor (BCP) acute lymphoblastic leukemia (ALL). We have now refined the prognostic strength of IKZF1 deletions by analyzing the effect of co-occurring deletions. Patients and Methods The analysis involved 991 patients with BCP ALL treated in the Associazione Italiana Ematologia ed Oncologia Pediatrica-Berlin-Frankfurt-Muenster (AIEOP-BFM) ALL 2000 trial with complete information for copy number alterations of IKZF1, PAX5, ETV6, RB1, BTG1, EBF1, CDKN2A, CDKN2B, Xp22.33/Yp11.31 (PAR1 region; CRLF2, CSF2RA, and IL3RA), and ERG; replication of findings involved 417 patients from the same trial. Results IKZF1 deletions that co-occurred with deletions in CDKN2A, CDKN2B, PAX5, or PAR1 in the absence of ERG deletion conferred the worst outcome and, consequently, were grouped as IKZF1plus. The IKZF1plus group comprised 6% of patients with BCP ALL, with a 5-year event-free survival of 53 ± 6% compared with 79 ± 5% in patients with IKZF1 deletion who did not fulfill the IKZF1plus definition and 87 ± 1% in patients who lacked an IKZF1 deletion ( P ≤ .001). Respective 5-year cumulative relapse incidence rates were 44 ± 6%, 11 ± 4%, and 10 ± 1% ( P ≤ .001). Results were confirmed in the replication cohort, and multivariable analyses demonstrated independence of IKZF1plus. The IKZF1plus prognostic effect differed dramatically in analyses stratified by minimal residual disease (MRD) levels after induction treatment: 5-year event-free survival for MRD standard-risk IKZF1plus patients was 94 ± 5% versus 40 ± 10% in MRD intermediate- and 30 ± 14% in high-risk IKZF1plus patients ( P ≤ .001). Corresponding 5-year cumulative incidence of relapse rates were 6 ± 6%, 60 ± 10%, and 60 ± 17% ( P ≤ .001). Conclusion IKZF1plus describes a new MRD-dependent very-poor prognostic profile in BCP ALL. Because current AIEOP-BFM treatment is largely ineffective for MRD-positive IKZF1plus patients, new experimental treatment approaches will be evaluated in our upcoming trial AIEOP-BFM ALL 2017.

Intragenic amplification of PAX5: a novel subgroup in B-cell precursor acute lymphoblastic leukemia?

Intragenic PAX5 amplification defines a novel, relapse-prone subtype of B-cell precursor acute lymphoblastic leukemia with a poor outcome.

Ras pathway mutations are prevalent in relapsed childhood acute lymphoblastic leukemia and confer sensitivity to MEK inhibition.

For most children who relapse with acute lymphoblastic leukemia (ALL), the prognosis is poor, and there is a need for novel therapies to improve outcome. We screened samples from children with B-lineage ALL entered into the ALL-REZ BFM 2002 clinical trial (www.clinicaltrials.gov, #NCT00114348) for somatic mutations activating the Ras pathway (KRAS, NRAS, FLT3, and PTPN11) and showed mutation to be highly prevalent (76 from 206). Clinically, they were associated with high-risk features including early relapse, central nervous system (CNS) involvement, and specifically for NRAS/KRAS mutations, chemoresistance. KRAS mutations were associated with a reduced overall survival. Mutation screening of the matched diagnostic samples found many to be wild type (WT); however, by using more sensitive allelic-specific assays, low-level mutated subpopulations were found in many cases, suggesting that they survived up-front therapy and subsequently emerged at relapse. Preclinical evaluation of the mitogen-activated protein kinase kinase 1/2 inhibitor selumetinib (AZD6244, ARRY-142886) showed significant differential sensitivity in Ras pathway-mutated ALL compared with WT cells both in vitro and in an orthotopic xenograft model engrafted with primary ALL; in the latter, reduced RAS-mutated CNS leukemia. Given these data, clinical evaluation of selumetinib may be warranted for Ras pathway-mutated relapsed ALL.

Small sizes and indolent evolutionary dynamics challenge the potential role of P2RY8-CRLF2-harboring clones as main relapse-driving force in childhood ALL.

The P2RY8-CRLF2 fusion defines a particular relapse-prone subset of childhood acute lymphoblastic leukemia (ALL) in Italian Association of Pediatric Hematology and Oncology Berlin-Frankfurt-Münster (AIEOP-BFM) 2000 protocols. To investigate whether and to what extent different clone sizes influence disease and relapse development, we quantified the genomic P2RY8-CRLF2 fusion product and correlated it with the corresponding CRLF2 expression levels in patients enrolled in the BFM-ALL 2000 protocol in Austria. Of 268 cases without recurrent chromosomal translocations and high hyperdiploidy, representing approximately 50% of all cases, 67 (25%) were P2RY8-CRLF2 positive. The respective clone sizes were ≥ 20% in 27% and < 20% in 73% of them. The cumulative incidence of relapse of the entire fusion-positive group was clone size independent and significantly higher than that of the fusion-negative group (35% ± 8% vs 13% ± 3%, P = .008) and primarily confined to the non-high-risk group. Of 22 P2RY8-CRLF2-positive diagnosis/relapse pairs, only 4/8 had the fusion-positive dominant clone conserved at relapse, whereas none of the original 14 fusion-positive small clones reappeared as the dominant relapse clone. We conclude that the majority of P2RY8-CRLF2-positive clones are small at diagnosis and virtually never generate a dominant relapse clone. Our findings therefore suggest that P2RY8-CRLF2-positive clones do not have the necessary proliferative or selective advantage to evolve into a disease-relevant relapse clone.

Analysis of binding sites on complement factor I using artificial N-linked glycosylation.

Factor I (FI) is a serine protease that inhibits all complement pathways by degrading activated complement components C3b and C4b. FI functions only in the presence of several cofactors, such as factor H, C4b-binding protein, complement receptor 1, and membrane cofactor protein. FI is composed of two chains linked by a disulfide bridge; the light chain comprises only the serine protease (SP) domain, whereas the heavy chain contains the FI membrane attack complex domain (FIMAC), CD5 domain, and low density lipoprotein receptor 1 (LDLr1) and LDLr2 domains. To better understand how FI inhibits complement, we used homology-based three-dimensional models of FI domains in an attempt to identify potential protein-protein interaction sites. Specific amino acids were then mutated to yield 20 recombinant mutants of FI carrying additional surface-exposed N-glycosylation sites that were expected to sterically hinder interactions. The Michaelis constant (K(m)) of all FI mutants toward a small substrate was not increased. We found that many mutations in the FIMAC and SP domains nearly abolished the ability of FI to degrade C4b and C3b in the fluid phase and on the surface, irrespective of the cofactor used. On the other hand, only a few alterations in the CD5 and LDLr1/2 domains impaired this activity. In conclusion, all analyzed cofactors form similar trimolecular complexes with FI and C3b/C4b, and the accessibility of FIMAC and SP domains is crucial for the function of FI.

Mutations and deletions of the TP53 gene predict nonresponse to treatment and poor outcome in first relapse of childhood acute lymphoblastic leukemia.

PURPOSE: In the clinical management of children with relapsed acute lymphoblastic leukemia (ALL), treatment resistance remains a major challenge. Alterations of the TP53 gene are frequently associated with resistance to chemotherapy, but their significance in relapsed childhood ALL has remained controversial because of small studies. PATIENTS AND METHODS: Therefore, we systematically studied 265 first-relapse patients enrolled in the German Acute Lymphoblastic Leukemia Relapse Berlin-Frankfurt-Mü nster 2002 (ALL-REZ BFM 2002) trial for sequence and copy number alterations of the TP53 gene by using direct sequencing and multiplex ligation-dependent probe amplification. RESULTS: We observed copy number and sequence alterations of TP53 in 12.4% (27 of 218) of patients with B-cell precursor ALL and 6.4% (three of 47) of patients with T-cell ALL relapse. Backtracking to initial ALL in 23 matched samples revealed that 54% of all TP53 alterations were gained at relapse. Within B-cell precursor ALL, TP53 alterations were consistently associated with nonresponse to chemotherapy (P < .001) and poor event-free survival (P < .001) and overall survival rates (P = .002). TP53 alterations also had a significant impact on survival within intermediate-risk (S2) and high-risk (S3/S4) relapse patients (P = .007 and P = .019, respectively). This prognostic significance of TP53 alterations was confirmed in multivariate analysis. Besides their clinical impact, TP53 alterations were associated with a higher fraction of leukemic cells in S/G(2)-M phase of the cell cycle at relapse diagnosis. CONCLUSION: Alterations of the TP53 gene are of particular importance in the relapse stage of childhood ALL, in which they independently predict high risk of treatment failure in a significant number of patients. Therefore, they will aid in future risk assessment of children with ALL relapse.

Clinical isolates of Streptococcus pneumoniae bind the complement inhibitor C4b-binding protein in a PspC allele-dependent fashion.

The complement system constitutes an important component of the innate immune system. To colonize their host and/or to cause disease, many pathogens have evolved strategies to avoid complement-mediated bacterial lysis and opsonophagocytosis. In this study, using a collection of 55 clinical isolates of Streptococcus pneumoniae, we demonstrate for the first time that pneumococci bind the complement inhibitor C4b-binding protein (C4BP). C4BP binding seems to be restricted to certain serotypes such as serotype 4, 6B, 7F, and 14, of which the strains of serotype 14 are the strongest binders. We show that bacteria-bound C4BP retains its functional activity and down-regulates the activation of the classical pathway. Thus, this major respiratory pathogen may escape immune recognition and eradication by the complement system. Furthermore, we show that C4BP binding varies between strains but is dependent on the expression of pneumococcal surface protein C, PspC of group 4. The study of the distribution of group 4 pspC locus shows that most of high-binder serotype 14 isolates harbor an allelic variant of group 4 pspC. Using PspC-negative mutant strains, we identified a new allelic variant of PspC (PspC4.4) as a major ligand for C4BP, revealing a new function for this important pneumococcal virulence factor. Thus pneumococci exploit host C4BP for complement evasion in a PspC allele-dependent manner.