Next Article in Journal
Long-Term Follow-Up of Patients with Conjunctival Lymphoma after Individualized Lens-Sparing Electron Radiotherapy: Results from a Longitudinal Study
Previous Article in Journal
Advances in Early Breast Cancer Risk Profiling: From Histopathology to Molecular Technologies
 
 
cancers-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prognostic Impact of Copy Number Alterations’ Profile and AID/RAG Signatures in Acute Lymphoblastic Leukemia (ALL) with BCR::ABL and without Recurrent Genetic Aberrations (NEG ALL) Treated with Intensive Chemotherapy

by
Marta Libura
1,*,
Karolina Karabin
1,
Paweł Tyrna
1,
Anna Czyż
2,
Hanna Makuch-Łasica
3,
Bożena Jaźwiec
2,
Monika Paluszewska
1,
Beata Piątkowska-Jakubas
4,
Magdalena Zawada
4,
Michał Gniot
5,
Joanna Trubicka
6,
Magdalena Szymańska
1,
Katarzyna Borg
3,
Marta Więsik
1,
Sylwia Czekalska
4,
Izabela Florek
4,
Maria Król
1,
Małgorzata Paszkowska-Kowalewska
1,
Lidia Gil
5,
Katarzyna Kapelko-Słowik
2,
Elżbieta Patkowska
3,
Agnieszka Tomaszewska
1,
Krzysztof Mądry
1,
Rafał Machowicz
1,
Tomasz Czerw
7,
Agnieszka Piekarska
8,
Magdalena Dutka
8,
Anna Kopińska
9,
Grzegorz Helbig
9,
Tomasz Gromek
10,
Krzysztof Lewandowski
5,
Marta Zacharczuk
2,
Anna Pastwińska
11,
Tomasz Wróbel
2,
Olga Haus
12,
Grzegorz Basak
1,
Jerzy Hołowiecki
7,
Przemysław Juszczyński
3,
Ewa Lech-Marańda
3,
Sebastian Giebel
7,† and
Wiesław Wiktor Jędrzejczak
1,†
add Show full author list remove Hide full author list
1
Department of Hematology, Transplantation and Internal Medicine, Medical University of Warsaw, 02-091 Warsaw, Poland
2
Department of Hematology, Blood Neoplasms and Bone Marrow Transplantation, Medical University of Wrocław, 50-137 Wrocław, Poland
3
Institute of Hematology and Transfusion Medicine, 02-776 Warsaw, Poland
4
Department of Hematology, Jagiellonian University Medical College, 31-008 Cracow, Poland
5
Department of Hematology and Bone Marrow Transplantation, Poznań University of Medical Sciences, 61-701 Poznań, Poland
6
Children’s Memorial Health Institute, 04-736 Warsaw, Poland
7
Maria Skłodowska-Curie National Research Institute of Oncology, Gliwice Branch, 44-102 Gliwice, Poland
8
Department of Hematology and Transplantology, Medical University of Gdańsk, 80-214 Gdańsk, Poland
9
Department of Hematology and Bone Marrow Transplantation, Medical University of Silesia, 40-032 Katowice, Poland
10
Department of Hematooncology and Bone Marrow Transplantation, Medical University of Lublin, 20-081 Lublin, Poland
11
Department of Tumor Biology and Genetics, Medical University of Warsaw, 02-106 Warsaw, Poland
12
Department of Clinical Genetics, Faculty of Medicine, Collegium Medicum in Bydgoszcz, Nicolaus Copernicus University, 87-100 Toruń, Poland
*
Author to whom correspondence should be addressed.
These authors jointly supervised this work: Sebastian Giebel and Wiesław Wiktor Jędrzejczak.
Cancers 2023, 15(22), 5431; https://doi.org/10.3390/cancers15225431
Submission received: 9 October 2023 / Revised: 7 November 2023 / Accepted: 9 November 2023 / Published: 15 November 2023
(This article belongs to the Section Cancer Pathophysiology)

Abstract

:

Simple Summary

Adult ALL is a highly aggressive blood cancer. Two classes of genetic aberrations are responsible for ALL: primary aberrations followed by secondary aberrations. Currently, primary aberrations are used for estimating patients’ risk in adult ALL. In this study, we reassessed the importance of primary and secondary copy number alterations (CNA) aberrations in intensively treated adult ALL patients in correlation to RAG/AID mutator enzyme expression. Primary aberrations alone specified the risk of 30% of patients. To define the prognosis of the remaining 70%, we identified high-risk and low-risk CNA profiles. We found the CNA profiles correlated with differential RAG/AID expression profiles. Furthermore, the outcome of CNAneg adult ALL was stratified by AID expression. Thus, we suggested mechanisms linking secondary aberrations with patients’ outcomes and mutator enzymes. Finally, we propose a revised version of risk stratification in adult ALL patients which incorporates primary and secondary genetic lesions.

Abstract

Adult acute lymphoblastic leukemia (ALL) is associated with poor outcomes. ALL is initiated by primary aberrations, but secondary genetic lesions are necessary for overt ALL. In this study, we reassessed the value of primary and secondary aberrations in intensively treated ALL patients in relation to mutator enzyme expression. RT-PCR, genomic PCR, and sequencing were applied to evaluate primary aberrations, while qPCR was used to measure the expression of RAG and AID mutator enzymes in 166 adult ALL patients. Secondary copy number alterations (CNA) were studied in 94 cases by MLPA assay. Primary aberrations alone stratified 30% of the patients (27% high-risk, 3% low-risk cases). The remaining 70% intermediate-risk patients included BCR::ABL1pos subgroup and ALL lacking identified genetic markers (NEG ALL). We identified three CNA profiles: high-risk bad-CNA (CNAhigh/IKZF1pos), low-risk good-CNA (all other CNAs), and intermediate-risk CNAneg. Furthermore, based on RAG/AID expression, we report possible mechanisms underlying the CNA profiles associated with poor outcome: AID stratified outcome in CNAneg, which accompanied most likely a particular profile of single nucleotide variations, while RAG in CNApos increased the odds for CNAhigh/IKZF1pos development. Finally, we integrated primary genetic aberrations with CNA to propose a revised risk stratification code, which allowed us to stratify 75% of BCR::ABL1pos and NEG patients.

Graphical Abstract

1. Introduction

B-cell acute lymphoblastic leukemia (B-ALL) is a biologically heterogeneous disease caused by specific genetic lesions leading to aberrant differentiation and proliferation of lymphoid progenitor cells [1]. Adult B-ALL is characterized by poor outcomes and a high rate of mortality when compared to children [2,3,4,5]. This is due to a particular genetic background and a higher rate of age-related comorbidities [6,7,8,9,10]. Although new therapies, as well as risk-adapted therapy protocols, led to a significant improvement in the outcome of adult ALL patients [11], the main reason for therapy failure in most cases is the emergence of new leukemic subclones leading to a relapse of the disease [12,13,14]. Thus, besides the standard risk group definition, the determination of the potential for leukemia’s evolution at the initial stage of the disease could be a key prognostication marker for those patients who would need further stratification.
A known hallmark of B-ALL development is the sequential acquisition of new genetic aberrations as a result of genetic instability and clonal evolution [13,15]. Chromosomal rearrangements and fusion genes, identified as primary aberrations, are the initiating events and key drivers of leukemogenesis, however, in most cases, they are insufficient for leukemia development. Secondary genetic events targeting B cell development genes are therefore an essential requirement for overt ALL [16,17]. These mutations, like copy number alterations (CNA; particularly gene deletions) and single-nucleotide variant (SNV) aberrations, have been frequently found within genes such as IKZF1, CDKN2A/B, and PAX5 and have been reported to cooperate with each other and with primary lesions, leading to a more aggressive phenotype. Importantly, except for rare genetic subtypes, e.g., MLL rearrangements, the progression of high-risk B-ALLs depends mostly on mutator mechanisms introducing new CNA/SNV-type aberrations [14,18]. Not surprisingly, disease relapse is associated with both clonal diversification and a higher level of CNA mutation burden [19,20,21].
Recently, two lymphocyte-specific mutator enzymes, RAG1/2 (recombination-activating gene) and AID (activation-induced cytidine deaminase), have been extensively studied in relation to genetic instability and oncogenicity in ALL. RAG1 and RAG2 genes encode recombinases, which introduce double-strand breaks into DNA during immunoglobulin/T-cell receptor genes’ rearrangement. However, aberrant targeting of RAG1/2 to non-IGH sites contributes to the development of driver mutations and the clonal evolution of ALL [22]. The AID gene encodes cytidine deaminase, which generates point mutations in immunoglobulin genes in a process known as somatic hypermutations in germinal B-cells. AID was documented to be a driver of oncogenic mutations in lymphomas [23,24]. However, recent studies have shown that aberrant activation of AID by infectious signals may also accelerate mutagenic processes leading to childhood ALL [13,25,26,27]. The above-presented mutational mechanisms have been well described in pediatric ALL, but the data are insufficient for adult ALL. Although recent studies suggest the involvement of RAG2-mediated aberrant recombinations in the evolution of adult B-ALL: t(9;22)/BCR::ABL1pos and BCR::ABL-like [28,29,30,31]; however, no comprehensive analyses on the expression of both enzymes across different genetic subgroups, in correlation to CNA status and clinical outcome, have been carried out in adult ALL so far.
Previous studies documented that both the mutation burden corresponding to clonal heterogeneity, and the type of the mutated gene (IKZF1, CDKN2A/B, or PAX5) may provide prognostic information [3,4,5,32,33,34]. However, controversial results have also been reported, particularly in adult ALL [2,35,36]. Although the factors that contribute to these discrepancies are not fully understood, they may be related to the complex nature of interactions between coexisting primary and/or secondary aberrations, different inclusion criteria (e.g., age, genetic background), or treatment protocols between independent trials. Because of these ambiguities, the fifth edition of WHO, as well as the ICC classification did not include a broad range of CNAs, as a diagnostic criterion for disease entities [37,38]. Importantly, most of these controversies apply to adult ALL. Whereas in pediatric ALL, the clinical significance of CNA has been well recognized and incorporated, e.g., into the integrated risk scoring system widely used in clinical trials, with the particular role of IKZF1 deletions coexisting with other concomitant CNAs. This CNA profile was first described by Stanulla et al. as IKZFplus and indicates a subgroup of patients with particularly poor outcomes (in our study, referred to as CNAhigh/IKZF1pos) [5,39,40]. As childhood ALL presents a cure rate of 85–90%, while most adult ALLs are poor-risk leukemias, it is important to understand better the possible role of mutational processes behind high-risk leukemias in adults.
Therefore, we carried out a detailed characterization of CNA profiles: CNA mutation burden, the type of mutated genes, and association with mutator enzyme AID/RAG1/2 expression—to establish their correlation with prognosis in an adult B-ALL population treated with a standard intensive protocol according to PALG. Based on the obtained results, we proposed a combined revision of the genetic risk classification integrating CNA data with well-established primary aberrations, which enabled further prognostic stratification of BCR::ABL1pos and NEG ALL patients. Parallelly, a comprehensive evaluation of RAG1/2 and AID expression signatures in correlation with CNA profiles allowed for further characterization of mutational processes behind selected ALL subtypes associated with poor prognosis.

2. Materials and Methods

2.1. Patients’ Characteristics, Treatment, Material Collection, and Detection of Fusion Genes and Mutations by Molecular Analysis

A total of 166 patients with B-ALL were enrolled in the present study. The patients were diagnosed, and 161 patients met the criteria to be treated according to the Polish Adult Leukemia Group (PALG)-ALL5 or PALG-ALL6 protocols between 2007 and 2017 in six Polish hematology centers. The remaining 5 patients were disqualified from intensive therapy and received palliative treatment; their outcomes were not analyzed in this study. Allogenic hematopoietic stem cell transplantation (alloHSCT) from either an HLA-matched sibling or an unrelated donor was performed in 77/161 (47%) patients out of the whole group.
The entire population was analyzed for common fusion genes: BCR::ABL1, TEL::AML1, E2A::PBX1, and MLL::AF4, using the RT-PCR protocol and primers designed according to the consensus of the European BIOMED-1 Concerted Action [41]. Subsequently, patients negative for the fusion genes listed above (102 patients) were screened for markers or surrogates of the BCR::ABL1-like phenotype (see details in Supplementary Materials). In parallel, cytogenetic analysis was performed on diagnostic bone marrow samples by PALG laboratories as previously described [42].

2.2. Treatment Protocols

The PALG ALL6 protocol encompassed all patients with newly diagnosed ALL. Patients were categorized for one of the first-line treatment options based on the subtype of the disease (ALL with or without the Philadelphia chromosome/BCR::ABL1 rearrangement) and their age (either up to 55 years old or older). Younger patients adhered to a “pediatric” approach consisting of a pre-treatment phase, intensive induction, and consolidation therapy, followed by extended maintenance. Additionally, the protocols incorporated intensive and prolonged central nervous system prophylaxis. AlloHSCT was considered for patients with BCR::ABL1-positive disease and those with post-remission positive minimal residual disease (MRD) assessed by the multiparameter flow cytometry (MFC) method (MRD ≥ 0.1% after induction and/or >0.01% after consolidation). Compared to PALG ALL6, PALG ALL5 was a similar MRD-driven regimen, except for a less individualized approach to the treatment of older patients over 55 years old. Further information on the treatment protocols is provided in Supplementary Materials.

2.3. Gene Expression Analysis of Mutator Enzymes by Quantitative Real-Time PCR

Gene expression of RAG1, RAG2, and AID mutator enzymes was determined in the entire B-ALL cohort by real-time polymerase chain reaction (qPCR; see details in the Supplementary Methods). The whole group was segregated into 2 subgroups based on RAG1, RAG2, and AID expression levels at diagnosis, i.e., higher or lower than the median. As in our setting, there was a highly positive correlation between RAG1 and RAG2 enzymes, we used only the RAG2 mRNA expression level for further analyses (Spearman’s rank correlation coefficient rho = 0.6; p < 0.05) [43].

2.4. CNA Detection: Multiplex Ligation-Dependent Probe Amplification (MLPA) and RT-PCR

CNAs were evaluated using a multiplex ligation-dependent probe amplification (MLPA) assay. The SALSA MLPA P335-B2 ALL-IKZF1 (MRC-Holland, Amsterdam, The Netherlands) was used according to the manufacturer’s protocol in a group of 94 cases whose genetic material was available. The deletions in the loci of the following 7 genes were scored as deleted or non-deleted: IKZF1, CDKN2A/B, PAX5, RB1, EBF1, ETV6, and BTG1 (see details in Supplementary Materials). Additionally, a cohort of 163 patients was evaluated for expression of 2 isoforms of IKZF1 mRNA (Δ3–6, Δ1–7; primers for RT-PCR protocol described by Iacobucci I. et al. [44].

2.5. Statistical Analysis

Complete remission (CR) rate and probability of overall survival (OS), relapse-free survival (RFS), and disease-free survival (DFS) were the study endpoints. The CR rate was defined according to previously published criteria [4]. DFS was defined as the time from the achievement of CR to hematological relapse, death, or the last follow-up. OS was defined as the time from diagnosis to death or the last follow-up.
Statistical differences between groups were tested with non-parametric tests. The Mann-Whitney U test was used for comparing continuous variables, the chi-squared test or Fisher exact test for categorical variables, and Spearman’s correlation test for correlating continuous data. The patients’ survival rates (OS, RFS, DFS) were compared with the log-rank test. In multivariate analyses, the hazard ratio was computed from the Cox proportional hazards model for survival rates (OS, RFS, DFS) and from the general linear model for categorical variables (CR). All statistical analyses were performed using the Statistica 13.3 software (TIBCO Software Inc., Palo Alto, CA, USA).

3. Results

3.1. Patients’ Characteristics; Frequency and Clinical Correlates of Primary Chromosomal Abnormalities

In the present study, 166 adults with B-ALL were investigated for the presence of fusion genes using the diagnostic workup according to the BIOMED-1 protocol. Comparisons of molecular results with cytogenetic data allowed for the identification of genetic subgroups according to the established primary chromosomal aberrations. The median age of the studied population was 37.5 years (range 18–70 years). According to the international criteria, 55 patients were classified into the very high-risk group, 94 into the high-risk group, and 17 into the standard-risk group. The demographic and clinical characteristics of the patients’ cohort are summarized in Table 1.
Out of the total 166 B-ALL patients analyzed, 55 (33%) had BCR::ABL1, 5 (3%) had KMT2A::AFF1 (MLL::AF4), 4 (2%) had TCF3::PBX1 (E2A::PBX1), and none had ETV6::RUNX1 (TEL::AML1) fusion genes (Figure 1). One hundred-two patients (62%) were found negative for the common fusion genes analyzed according to the BIOMED-1 protocol. A subsequent analysis of these 102 patients identified 16 cases accompanied by the markers or surrogates of the BCR::ABL1-like phenotype (10% of the total population). All of the BCR::ABL1-like ALL represented JAK::STAT2–class aberrations: 8 patients with overexpression of CRLF2 gene, 7 with P2RY8::CRLF2 rearrangement, 1 with CRLF2 gene point mutation, 6 with JAK2 (see details in Table S1).
When we matched the remaining group of 86 patients, who had no identified molecular aberrations, with cytogenetic data, 23 patients could be defined as a subgroup with poor risk genetic characteristics: 15 with complex karyotype, 6 hyperdiploid, 2 hypodiploid, and 1 with another aberration of poor risk. Out of the remaining 63 cases, 42 presented normal karyotypes, 1 had other aberrations of unknown significance, and 20 lacked the cytogenetic analysis. As this latter subgroup represented a “real world” population of ALL patients for whom the initial risk stratification is unknown, for the sake of further analyses in this manuscript, we grouped these patients under the label of “NEG ALL” (38%).

3.2. Clinical and Biological Characteristics of RAG2 and AID Mutator Enzymes’ Expression—In Correlation with Primary Genetic Subgroups

In order to evaluate the potential for RAG2- and AID-mediated genetic instability, all patients were assigned into subgroups based on both enzymes’ mRNA expression levels that were higher or lower than the median (see Table 1B and Table S2). Thus, using the integrated AID/RAG2 profile, we could identify 4 signatures of differential AID and RAG2 expression (referred to here as “sig. 1–4): sig. 1 AIDlow/RAG2low was found in 18.6% of ALL patients, who showed low expression levels of both enzymes; sig. 2 AIDhigh/RAG2low accompanied 32.7% of ALL population with high AID but low RAG2; sig. 3 AIDlow/RAG2high was found in 31.4% of ALL, that showed high RAG2 while little or no AID. Sig. 4 AIDhigh/RAG2high represented 17.3% of ALL who had both enzymes in abundance, suggesting co-synergic involvement of AID with RAG2. In summary, the highest frequency of sig. 3 (AIDlow/RAG2high) was found in the BCR::ABL1pos subgroup, while the lowest one was in the NEG B-ALL (43% vs. 22%, p = 0.005). In contrast, a higher incidence of AIDhigh/RAG2low (sig. 2) was reported for the NEG B-ALL subgroup when compared to the BCR::ABL1pos population (47% vs. 21.6%, p = 0.005). Detailed data on AID/RAG2 profiles and their distribution in correlation to primary genetic aberrations are presented in the Supplementary Results.

3.3. Clinical and Biological Characteristics of Secondary CNA Aberrations—In Correlation with Primary Genetic Subgroups

In the next step, we analyzed the characteristics of CNAs: CNA mutation burden and gene types in correlation with demographic data and B-ALL primary genetic aberrations in our study cohort (see Table 1C and Table S3). In summary, 1 or more deletions were observed in 66 out of 94 patients (70%), while 28 patients (30%) had no CNAs. The highest prevalence of CNAs was observed in the BCR::ABL1 and BCR::ABL1-like subgroups (80% and 83.3%), followed by an intermediate one in the NEG B-ALL (68.6%), and the lowest in MLL gene rearrangement (33%). The detailed data on secondary CNA characteristics and correlation with genetic subgroups are further presented in the Supplementary Results.

3.4. Correlation of Secondary CNA Aberrations with RAG2/AID Signatures

Furthermore, in order to establish a functional link, we correlated the integrated AID/RAG2 expression signatures with CNAs’ profiles (see Table 2 and Table S2). In summary, although we could confirm the presence of CNAs in each subgroup according to the integrated AID/RAG2 signature profile, we report that the abundance of RAG2 alone correlated rather with higher CNA levels and IKZF1 deletions, particularly in the BCR::ABL1pos context (p = 0.001), while a lower CNA mutation number or CNAneg was more frequently accompanied by AIDhigh/RAG2low, particularly in NEG ALL context, suggesting involvement of secondary mutations other than CNA, e.g., SNV (see Table S2 for details). Significantly, sig. 4, i.e., with parallel AID and RAG2 abundance, was totally absent in CNAneg patients, again emphasizing the link of both enzymes with CNA-type mutagenesis. Altogether, these results demonstrate a correlation of different CNA profiles with differential mutagenic properties linked to mutator enzymes or other mutagenic stimuli (AIDlow/RAG2low). The detailed description of correlation data between the integrated AID/RAG2 expression signatures with CNAs’ profiles is further provided in the Supplementary Results.

3.5. The Outcome of Intensively Treated B-ALL Patients in Relation to Established Primary Aberrations

In the next part of the study, we verified the clinical outcomes of primary and secondary aberrations in the adult ALL population treated according to PALG protocol. Among the 166 patients analyzed for primary aberrations using molecular testing, 161 were eligible for intensive induction therapy, and 152 out of 161 were evaluated for response to induction, while 9 patients died before the remission evaluation. In the whole cohort, complete remission (CR) was achieved in 131/152 (86%) of the study population, which is consistent with the previous observations [45]. The median follow-up of the study group was 37.9 months and the median survival reached 20.5 months, which is comparable to results in other studies of similar populations. The estimated 4-year overall survival (OS) for the whole population was 34%, with a standard error of ±4%.
The impact of primary chromosomal aberrations on survival is shown in Table 3 and Figure 2A,B. Our data confirmed the poor 4-year outcome for patients harboring MLL rearrangements (OS 53 ± 25%), complex karyotype (OS 20 ± 10%), hyperdiploid karyotype (OS 0%), and in the BCR::ABL1-like subgroup (OS 10 ± 9%). As the last subgroup represented JAK::STAT pathway abnormalities, the poor prognosis of these patients has already been described in other cohorts [46]. All the cases with E2A::PBX exhibited a good outcome. Interestingly, in our study cohort, a relatively good outcome was also documented for the BCR::ABL1pos subgroup (54 ± 8%). The remaining patients, for whom neither molecular nor cytogenetic abnormalities with established prognostic importance were detected, thus indicating obscure prognosis, were referred to here as NEG ALL and showed an intermediate risk (OS 32 ± 8%; Table 3).
In multivariate analysis, the BCR::ABL1 fusion gene independently predicted a lower risk of death, and hyperdiploid karyotype was independently associated with an increased death rate (OS for BCR::ABL1: p = 0.034 with HR: 0.55, 95% CI, 0.32–0.96; for hyperdiploid: p = 0.006 with HR: 3.91, 95% CI, 1.49–10.31). BCR::ABL1 fusion independently predicted a lower incidence of relapse, while MLL::AF4 independently associated with a higher risk of relapse (RFS for BCR::ABL1: p = 0.005 with HR: 0.34, 95% CI, 0.16–0.72; for MLL rearrangements: p = 0.013 with HR: 4.40, 95% CI, 1.37–14.13) after adjustment for WBC, age, and other primary aberrations in the total ALL series (Table 3B).
Subsequently, based on the obtained data, we have assigned B-ALL patients to four categories: (1) “bad” primary aberrations, which grouped patients with all markers of poor prognosis (i.e., MLL rearrangements, BCR::ABL-like aberrations, complex, hyperdiploid, and hypodiploid karyotype), (2) “good” ones (E2A::PBXpos), as well as “intermediate” ones (IM), which included 2 subgroups: (3) BCR::ABL1pos and (4) NEG ALL (Figure 2C,D; Table 3A). As the prognoses of the first and second categories of aberrations are well established (either “bad” or “good”), the prognoses of the 2 remaining IM subgroups: BCR::ABL1pos and NEG ALL needed further clarification. Thus, we asked the question of whether CNA markers may help to identify cases with a high and low risk of disease progression. As BCR::ABL1pos and NEG ALL genetic subgroups may constitute up to 70% of adult ALL, it was important for us to find relatively easily accessible prognostic markers like CNA for more relevant therapeutic decisions. In parallel, we evaluated the prognostic impact of CNA on the total ALL population as a historical group to compare with the results obtained by the others.

3.6. Prognostic Relevance of Secondary CNA Mutations in BCR::ABL1pos and NEG ALL

When we compared survival between CNApos and CNAneg patients, we did not observe any significant differences (Table 4A). As the CNApos subgroup was heterogeneous with respect to the mutation burden and the gene type, we stratified patients according to the number of CNAs and compared their survival. After several intermediate analyses, which are presented in detail in the Supplementary Results section, we found rationale to stratify patients according to genetic instability level: the first subgroup was labeled as “CNAlow” and included both CNAneg and 1CNA patients, while the second subgroup was described as “CNAhigh” and included patients who harbored ≥2CNAs (see Table 1C for details). Comparing these two groups, we observed significant differences in survival in the NEG cohort as well as in the whole series: patients with CNAhigh had lower OS and RFS rates than CNAlow patients (see Figure 3, Table 4A and Supplementary Results for details). Only a trend towards significance was observed in the BCR::ABL1pos subgroup. The data suggested that a high mutation burden together with the cooperation of concomitant gene mutations may play a predominant role in prognosis rather than a single CNA.
Thus, in the next step, we wanted to verify if any of the particular genes had a higher impact on survival than the others. We started with IKZF1, as it was the most frequent deletion in our study cohort, and found that its presence significantly correlated with the inferior outcomes of patients from NEG B-ALL. A trend toward statistical significance was observed for the BCR::ABL1pos subgroup (see Figure 4 and Table 4A). For a detailed description of the obtained data, see Supplementary Results. In contrast, gene deletions like CDKN2A/B and PAX5 had no prognostic impact in subgroup analysis or the entire series.
Since the total IKZF1pos population remained heterogeneous regarding the CNA burden, we assessed the difference between IKZF1 deletion alone and IKZF1 mutation coexisting with other CNAs (1CNA/IKZF1pos vs. CNAhigh/IKZF1pos). Of note, Stanulla et al. were the first to describe the IKZF1pos profile that co-occurred with deletions in CDNK2A/B, PAX5, or PAR1 in the absence of ERG deletion, referring to it as the IKZFplus subset, which defined poor prognostic subgroups among pediatric B-ALL [39]. In turn, in our study, the CNAhigh/IKZF1pos subgroup included all IKZF1pos patients accompanied by other gene deletions. As expected, the CNAhigh/IKZF1pos patients showed worse survival when compared to the 1CNA/IKZF1pos series. This difference was statistically significant in the BCR::ABL1pos cohort and the entire series, while only a trend towards significance was found in the NEG subpopulation, probably due to a small number of 1CNA/IKZF1pos patients. This suggests that the negative effect of IKZF1 deletions within the total IKZF1pos population was found exclusively in patients with CNAhigh/IKZF1pos profile (see Figure S1, Table S4 for further data). A detailed description of these results is provided in the Supplementary Results.
Then, we attempted to verify if IKZF1 had an adverse impact on other genes in the CNAhigh population. Interestingly, when we compared the outcome of CNAhigh/IKZF1pos vs. CNAhigh/IKZF1neg series, the presence of IKZF1 deletions conferred lower survival. Despite a relatively small number of cases in the subgroup analysis, these differences were statistically significant for NEG ALL, BCR::ABL1pos cohort, and whole ALL series. These data may suggest again that the negative outcomes observed in the total CNAhigh population may be mostly attributed to patients with CNAhigh/IKZF1pos (see Table S4 and Figure S2 for further results). A detailed description of these data is provided in the Supplementary Results section.
Furthermore, as some studies reported [47,48] on the negative prognostic impact of CDKN2A/B deletions when accompanied by other deletions, particularly in BCR::ABL1pos context, we compared the impact of IKZF1 deletions in the CNAhigh/CDNK2A/Bpos cohort. Interestingly, CDNK2A/Bpos patients accompanied by IKZF1 deletion showed worse outcomes when compared to the subgroup without IKZF1. Despite the relatively small number of patients, this effect was particularly evident in the BCR::ABL1pos series (p = 0.01), suggesting that the coexistence of a CDNK2A/B gene deletion alone may not affect the outcome of the CNAhigh population, but rather through cooperation with IKZF1 deletions (Figure S3). A detailed description of these data is presented in the Supplementary Results section.
In summary, these data document that CNAhigh/IKZF1pos status is the predominant CNA profile with a clear negative impact on the treatment outcome in our series. As the dominant negative Ik6 isoform of IKZF1 deletion was found in other study cohorts to be associated with the worst outcome, we compared the Ik6 prognostic power within the CNAhigh/IKZF1pos population [4,33]. Despite a small group of patients, we found a negative impact of Ik6 on relapse rate among the poor-risk CNAhigh/IKZF1pos subgroup.
Finally, we noticed that the prognostic impact of CNAhigh/IKZF1pos varied in different age groups; a statistically significant difference in outcomes was obtained for younger patients (<40 y.), while no difference was observed for the older population (Figure S6). This suggests age-related comorbidities confounding the power of genetic markers in elderly patients, which can contribute to discrepancies in the prognostic impact of CNAs between clinical trials.
Thus, based on the above CNA data, we proposed a revised coding for risk stratification of adult B-ALL by CNA profile, which comprised 3 key subgroups, accompanied by different clinical impact: (1) “bad”-CNAs comprising CNAhigh/IKZF1pos cases (OS 9 ± 6%; RFS 14 ± 9%); (2) relatively „good”-CNAs composed of 1CNA/IKZF1pos or 1-4CNA/IKZF1neg cases (OS 43 ± 9%; RFS 45 ± 9%); (3) the remaining CNAneg ALL with unresolved prognosis–labeled as “intermediate”-CNAs (OS 43 ± 10%; RFS 59 ± 14%; p = 0.003 for bad-CNA vs. good-CNA; see Figure 5, Table 4A).
The revised CNA-risk index was valid for the whole series, but most importantly it stratified the prognosis of the two intermediate subgroups: 29% of NEG and 44% of BCR::ABL1pos ALL were accompanied by bad-CNA, while 41% of NEG and 34% of BCR::ABL1pos were associated with good-CNA profile (Figure 5C–F; Table 4A). Though the bad-CNA profile had equally poor prognosis both in the NEG and BCR::ABL1pos subgroup, the good-CNA profile was associated with much better survival in the BCR::ABL1pos patients (OS 89 ± 10%; RFS 89 ± 10%) when compared to the NEG ALL (n = 14, OS 28 ± 13%; RFS 31 ± 17%). It is possible that at least some good-CNApos cases from NEG ALL were in fact BCR::ABL-like patients, which are frequently accompanied by CNAs with poor prognosis by themselves.
Finally, the results of the multivariate Cox regression analysis were consistent with the results of the univariate analysis both in the NEG and BCR::ABL subgroups. CNAhigh mutation burden, IKZF1 deletion as well as the bad-CNAhigh/IKZF1pos status were independent prognostic factors associated with worse OS and RFS in the NEG subgroup (OS p = 0.003 with HR: 4.95, 95% CI, 1.70–14.39 for CNAhigh; p = 0.011 with HR: 4.12, 95% CI, 1.37–12.33 for IKZF1pos; p = 0.007 with HR: 4.85, 95%CI, 1.56–15.09 for CNAhigh/IKZF1pos; Table 4B). In the BCR::ABL1pos series, CNAhigh/IKZF1pos had an independent prognostic impact on an increased risk of relapse (p = 0.030 with HR:7.65, 95% CI, 1.21–48.31), while good-CNAs, which lacked CNAhigh/IKZF1pos, independently predicted a lower risk of death (p = 0.04 with HR: 0.11, 95% CI, 0.01–0.90; Table 4B).

3.7. Combined Genetic Risk Classification by Revised Coding of Primary and Secondary Aberrations in Adult B-ALL

Therefore, taking into account the data shown above, we could propose a novel combined risk-adapted classification for adult B-ALL: by incorporating the stratification of NEG and BCR::ABL1pos subtypes by secondary CNA into the existing genetic classification based on primary aberrations with well-established prognostic effect (Figure 2C,D and Figure 5C–F) we could further refine prognosis of the entire adult ALL population. The entire cohort was divided into 3 subgroups, which were created by co-segregating different categories with comparable prognostic impact: intermediate-, low-, and high-risk features. The intermediate one (1) included intermediate-CNAneg cases from both NEG and BCR::ABL1pos subgroups as well as good-CNA cases from NEG ALL (1CNA/IKZF1pos or 1-4CNA/IKZF1pos) identifying 33% adult ALL from our series; (2) the low-risk subgroup comprised E2A::PBXpos cases and good-CNAs from BCR::ABL1pos series (1CNA/IKZF1pos and 1-4CNA/IKZF1neg) and contained 14% patients from our cohort; (3) while the 3rd one included both bad primary aberrations (MLL::AF4, BCR::ABL-like, complex, hyper/hypodiploid karyotype) and bad-CNAs (CNAhigh/IKZF1pos) from NEG or BCR::ABL1pos ALL, and was associated with the worst outcome in 53% of total adult ALL (see Figure 6).
These data indicate that more than half of adult ALL patients demonstrated poor outcomes, which has been revised compared to the frequency of 27% bad-risk ALL identified using only the well-established primary aberrations. The univariate and multivariate analyses for the combined revised risk stratification of adult ALL treated with PALG protocol are presented in Table 5.
Finally, Figure 7 summarizes the flowchart of revised risk classification based on both primary and secondary aberrations, leading to a hierarchical definition of genetic subcategories. The identified CNA-based risk subgroups can be used to refine patient treatment according to a more detailed description of their leukemia.

3.8. Prognostic Relevance of RAG2 and AID Expression in CNAneg and CNApos B-ALL Population

Interestingly, the three identified CNA profiles were associated not only with different outcomes but also with differential AID/RAG2 signatures, suggesting a possible link between disease progression and the causal mutagenic mechanisms. Importantly, markers of CNA-based mutagenesis could not stratify CNAneg cases, and these patients still constituted 31% and 20% of cases from the NEG and BCR::ABL1pos subtypes, respectively. Thus, we tested whether the AID/RAG2 signature was associated with disease progression in CNAneg patients. Unfortunately, in this study, we did not have a chance to investigate SNVs or larger genomic lesions, which probably could drive the leukemogenesis of this subgroup. However, as SNV mutagenesis is related to AID activity, we speculated that a high AID expression may correlate with disease progression in CNAneg cases. Indeed, despite a small number of cases, we report here a positive correlation between the AIDhigh/RAG2low expression signature and the incidence of adverse events in the CNAneg subtype; patients who relapsed in 71% (5/7) were accompanied by high AID expression when compared to 20% (3/15) of those who remained in CR in the entire B-ALL series (p = 0.03; Table S5A). Such correlation was reversed for the CNApos series: only 20% (5/25) of AIDhigh/RAG2low signature was found among patients who eventually relapsed (p = 0.09). Interestingly, opposite data were reported for RAG2 abundance in the CNAneg series: a negative correlation with relapses (0/7) while a positive one with remission status (8/14; p = 0.01).
These data prompted us to compare survival in the CNAneg subgroup according to both AID and RAG2 expression levels. Interestingly, high AID and low RAG2 expressors were associated with low RFS and OS in this subgroup (Figure S4, Table S5B). In multivariate analysis, high AID and low RAG2 independently predicted a higher incidence of relapse in the CNAneg subset (Table S5C). These data, although based on a limited series of patients, allowed us to speculate that genomic markers other than CNA, which associate with AIDhigh/RAG2low signature may contribute to the outcome of CNAneg patients. Thus, although we did not directly analyze the genetic background of CNAneg ALL, we can assume that CNAneg cases with AIDhigh/RAG2low signature, along with CNAhigh/IKZF1pos cases may also represent a profile of “bad” secondary mutations. This again may help in the selection of new poor- and low-risk subgroups among CNAneg patients, that can be used to refine patients’ treatment (about 17% of total ALL; see Figure 7).
Finally, in the same manner, we decided to verify whether both enzymes correlate with clinical outcomes in the CNApos series, although an inverse pattern of correlation was assumed. This time, due to larger subgroups, the analysis could be performed separately for the NEG ALL and BCR::ABL1pos series. Indeed, in the CNApos subgroup, AIDhigh was found to correlate with better survival in NEG ALL, while AIDlow/RAG2high profile was associated with poor outcomes in BCR::ABL1pos (Figure S5). Importantly, both AID and RAG2 did not impact the outcome when evaluated independently from CNA status in the whole series, as well as in genetic subgroups.

4. Discussion

In this study, we reassessed the value of primary and CNA-type secondary aberrations in adult ALL treated with standard remission induction protocol according to PALG in relation to mutator enzyme AID/RAG1/2 signatures. Our principal aim was to identify risk-defined subgroups and mutational processes underlying adult leukemias with poor outcomes.

4.1. Revised Risk Index Based on Primary Aberrations Characterized the Prognosis of 30% of ALL Patients; Rationale for Further Stratification of NEG and BCR::ABL1pos Subgroups

We have confirmed that established primary aberrations allowed 30% of adult ALL to be assigned to well-defined risk subgroups, either to the bad prognosis subgroup (28%; KMT2A::AFF1 (MLL::AF4), BCR::ABL1-like phenotype, complex karyotype, and hyper-/hypodiploid karyotype); or the good prognosis subgroup with TCF3::PBX1; E2A::PBX1 (2%). The remaining patients with BCR::ABL1 or lacking identified genetic markers (NEG ALL) represented an intermediate outcome in need of further stratification.
We proved a relatively good prognosis of BCR::ABL1pos patients [49,50,51,52,53], however, we questioned whether all BCR::ABL1pos patients may respond well to therapy. NEG ALL, which was negative for all investigated fusion genes and BCR::ABL1-like surrogate markers, was composed mainly of cases with normal karyotype (70%), or lacking cytogenetic analysis (30%). This subgroup thus represented a “real world” population of ALL patients with unknown prognoses. Importantly, although these two intermediate prognosis subgroups: NEG and BCR::ABL1pos ALL covered 70% of adult ALL, their treatment outcome could not be further improved just by using conventional markers for initial risk stratification. Classic prognostic factors like age, WBC, fusion genes, and MRD are used for risk groups’ assignment [54]. However, routine practice shows that recurrent disease is not restricted only to the high-risk group. Additionally, the application of high-resolution techniques revealed that the genomic landscape particularly of NEG ALL is complex and heterogeneous [1,7,35,55]. Yet, the full spectrum of diagnostic workup according to the recent WHO/ICC classification is not always available, with many single aberrations lacking proven prognostic significance [37,38]. In relation to poor treatment outcomes of adult ALL, all these facts present a challenge for genetic classification based on risk.
Therefore, we tested our ALL cohort for hidden genomic lesions: secondary copy number alterations (CNA), which cooperate with primary aberrations, and thus may impact the prognosis, and the choice of risk-adapted therapy. Their incidence reaches 60–70% across both childhood and adult ALL populations, representing an easily accessible marker for risk stratification [16,33]. However, as some results from different clinical trials were contradictory [2,3,4,36], we aimed to verify their clinical impact for intensively treated ALL patients according to PALG. As RAG and AID—known mutator enzymes—were documented to increase clonal heterogeneity [13,22], we decided to verify whether RAG2/AID expression signatures, may relate to CNA characteristics, their prognostic impact, and thus provide additional information on relapse potential.
Using the MLPA technique we detected CNAs in seven genes in 94 patients. Furthermore, RAG2 and AID expression signatures were measured on a broader population of 166 patients with available RNA material. The distribution and frequencies of CNAs were generally consistent with those from other study cohorts in the whole series, as well as across genetic subgroups, with the highest CNA frequency in BCR::ABL1pos (80% CNA) and BCR::ABL1-like subgroups (83% CNA), intermediate in NEG B-ALL (69% CNA), while the lowest in patients with MLL rearrangements [2,33]. This differential distribution already suggests a possible correlation with subtype-specific mutational processes and/or cooperation of concomitant gene aberrations.

4.2. RAG Is Associated with High CNA Mutation Burden, While AID Is Frequently Observed in CNAlow Patients

In order to elucidate a functional link between CNAs and the mutagenic process, we identified four different signatures of AID and RAG2 expression and found their correlation with particular CNA profiles. Importantly, the signature with the highest expression of both enzymes was totally absent in the CNAneg subgroup, re-emphasizing AID and RAG2 involvement in CNA formation. The highest CNA mutation burden (≥3), IKZF1 deletions, and particularly CNAhigh/IKZF1pos profile were associated with high RAG2 and low AID (AIDlow/RAG2high) signature. This suggests an increased RAG-mediated clonal diversity involved in the progression of these ALL subtypes. Interestingly, the highest incidence of the above CNAs as well as high RAG2 expression was observed in BCR::ABL1pos and BCR::ABL-like ALLs from our series. Increased RAG2 expression has been already described as a result of constitutive VDJ recombination activation by ABL and JAK-STAT kinase-activating lesions [30,31]. Additionally, our data indicate that IKZF1 deletion by itself activates RAG1/2 expression, thus increasing the risk for further deletions, which can lead to disastrous consequences [56].
In contrast to the above data, low CNA mutation burden (0–1), CNAs lacking IKZF1, and the absence of CNA (CNAneg) correlated with a higher representation of the AIDhigh/RAG2low signature, which signifies a possible role of AID-mediated mutagenesis underlying progression of these leukemias. AID is a member of the family of cytidine deaminases leading to substitution mutations, therefore its involvement in SNV-type rather than CNA-type mutagenesis may be expected. However, a co-synergic involvement of AID with RAG in CNA formation has been described as well [13]. Previous studies already reported higher AID expression in ALL lacking common primary aberrations (our NEG ALL subgroup) [57]. Here we expand these data by reporting prevalent AIDhigh/RAGlow signature not only in NEG ALL but also in the CNAneg series, suggesting its involvement in mutagenic processes other than CNA rearrangements in the absence of RAG.
Interestingly, in contrast to RAG2, AID was found to be also upregulated by environmental factors, like repeated exposure to inflammatory stimuli, paralleling chronic infections in childhood [13,22,25,27]. Recent studies also describe a depleted microbiome as an accelerator of ALL development, while metabolites of a healthy microbiome prevented leukemia through AID inhibition [58,59,60]. Importantly, this interplay of AID with inflammation may be a factor stimulating clonal heterogeneity, e.g., in the infectious environment of immunodeficient ALL patients. Thus, while RAG has a well-established role in pro-B cells, expression of AID represents a recently discovered threat to lymphoid precursor genome integrity, but also a link to other cancers associated with chronic inflammation [13,26,61,62,63].

4.3. Prognostic Impact of CNAs Aberrations in Correlation with Data from Other Studies

Most importantly, CNA profiles in correlation with various AID/RAG2 ratios—impacted treatment outcomes in our series. Thus, on the basis of the clinical data, we have proposed a revised CNA-based index stratifying prognosis of the intermediate subgroups from our series: NEG and BCR::ABL ALL. The CNA-based classification defined 3 subgroups with different CNA profiles and prognoses: bad-CNA (CNAhigh/IKZF1pos), relatively good-CNA (all other CNAs, lacking CNAhigh/IKZF1pos), and CNAneg with an intermediate prognosis.
We report here that bad-CNA: CNAhigh/IKZF1pos is the most important CNA profile identifying patients with the worst outcome both in NEG and BCR::ABL subgroups of intensively treated adult ALL, and with the high RAG2 expression as a causal driver of the mutagenic process. Additionally, despite a small group of patients, we report a particularly unfavorable outcome for the dominant negative Ik6 variant within the CNAhigh/IKZF1pos subgroup. Our data are consistent with the observations of many authors [4,33]. Stanulla et al. were the first to describe the CNAhigh/IKZF1pos profile as defining very poor prognosis in pediatric ALL (so-called IKZFplus subgroup) [39]. Since then, these data have been confirmed in several studies both in childhood and adult ALL [3,5,33,40,64]. However, opposite results have also been published. Discrepancies mainly concerned adult ALL cohorts, while in childhood ALL the prognostic impact of CNA profiles was consistently well documented [2,35]. This points to age-related factors frequently confounding the prognostic power of genetic lesions. The prevalence of IKZF1 deletions increases with age, thus studies omitting younger adults or adolescents will have a higher incidence of IKZF1pos patients with an overlapping impact of age-related comorbidities [6,7]. Importantly, we confirmed this hypothesis in our study, showing that the prognostic impact of CNAhigh/IKZF1pos was restricted to the younger population only (<40 y).
Another factor contributing to clinical discrepancies may relate to different group definitions based on genetic background. For example, many studies have chosen to analyze the impact of CNAhigh/IKZF1pos among BCR::ABL1neg patients [2,35]. In contrast, our NEG B-ALL population was not only negative for BCR::ABL1 fusion but also excluded all other known primary aberrations, e.g., the majority of BCR::ABL1-like cases.

4.4. Prognostic Impact of CNAs Aberrations in Particular Context of BCR::ABL1pos Subgroup

Importantly, the BCR::ABL1pos series from our cohort was stratified not only by bad-CNA but also by good-CNA profile with favorite prognosis. Thus, the presence of good-CNAs could contribute to a relatively good outcome of our BCR::ABL1pos series as a whole group. This observation is of importance, as some studies did not show any correlation between CNA and outcome for adult BCR::ABL1pos ALL in the final multivariate analyses [2,33,64]. This may be due to the complex nature of chemotherapy resistance in BCR::ABL1pos leukemias. According to the literature, most BCR::ABL1pos patients relapse due to kinase domain mutations, which represent SNV-type aberrations, while 30–40% of patients relapse without kinase domain mutations [31]. We assume that the latter population may correspond to our bad-CNA with CNAhigh/IKZF1pos profile, as IKZF1 itself was documented to confer poor responsiveness to ABL kinase inhibitor therapy in the absence of mutation [65]. Interestingly, 5/7 patients who relapsed in our BCR::ABL1pos series presented CNAhigh/IKZF1pos deletions. Although we did not look for ABL1 kinase domain point mutations in this subgroup, we assume that RAG-mediated CNA-type mutagenesis is more relevant in this subtype than SNV-type ABL1 lesions, in line with high RAG2 expression. However, overlap of CNAhigh/IKZF1pos with SNV may also be considered, in line with observations of Koptyra et al. suggesting involvement of ROS-mediated SNV-type mutagenesis within ABL point mutations hotspot [66]. In contrast, good-CNA BCR::ABL1pos patients from our series responded well to TKI treatment (see Figure 7).
Other studies report on CDKN2A/B deletions as a high-risk marker for adult BCR::ABL1pos ALL [47,48]. Although on a limited series, we document, that CDKN2A/B deletions confer an adverse prognostic impact only when accompanied by IKZF1 deletions. CDKN2A/B deletions were frequently accompanied by IKZF1 deletions in BCR::ABL1pos subgroups. Therefore, it is possible, that the poor prognosis of the CNAhigh/CDKN2A/Bpos subgroup could be attributed to the frequent coexistence of IKZF1 with CDKN2A/B deletions as a part of CNAhigh/IKZF1pos population.
In the NEG ALL subgroup, the clinical impact of good-CNA was labeled as “good”, although remained unresolved. This was mostly due to the unknown nature of other mutations, whose interplay with good-CNA could affect the prognostic power.

4.5. Prognostic Impact of AID and RAG Expression in CNAneg vs. CNApos Subgroups

Importantly, the evaluation of AID expression helped in the prognostic stratification of the CNAneg subgroup (30% of total ALL). High AID expression has already been reported by Swaminathan et al. as a strong indicator of poor overall survival [13]. However, as part of this evaluation was performed on a pediatric group, here we present data for adult patients, specifying that the prognostic significance of AID was reserved only for CNAneg patients. Although we do not know the mutation profile associated with the AIDhigh/RAG2low signature, our data signify distinct mutagenic processes accompanying disease progression in these high-risk CNAneg patients, probably related to SNV-type lesions. Their association with inflammatory signals is also possible.
In contrast, in the CNApos series high RAG2 marked patients with poor outcomes, which resulted from the presence of CNAhigh/IKZF1pos mutations. Thus, our data suggest that patients whose ALL involves increased RAG2 or AID-mediated clonal heterogeneity exhibit a more aggressive disease because RAG and AID increase the probability of acquiring high-risk mutations, distinct in the context of CNApos and CNAneg subgroups.

4.6. Combined Revised Risk Classification of Adult ALL: CNA Data with Primary Aberrations Reclassify Prognostic Index of Adult ALL

Finally, on the basis of an extensive analysis of genetic data, we proposed a combined version of risk classification for intensively treated adult ALL, incorporating the CNA data into the risk index for primary aberrations. This is in contrast to the previous study by Moorman et al., who proposed a revised UKALL genetic risk classification based on key established primary aberrations as a stronger prognostic marker than secondary deletions [2]. Although the treatment protocol may affect prognostic factors, our study provides evidence that, at least for intensive treatment protocols such as PALG, besides MLL and hyperdiploid aberrations, secondary CNA aberrations stratify prognosis stronger than the established conventional factors. For example, these aberrations may confer a totally different outcome in patients with the same primary aberration.
Unfortunately, we did not have a chance to correlate our results with MRD data for all patients. However, independently from the prognostic power of MRD by itself, one can imagine, that the blast cells possess the potential to evolve, which can change their genotype and phenotype. Thus, the negative MRD results assessed using techniques following initial profiles—may not always recapitulate the risk of relapse. Whereas initial CNA profiles inform on the risk of adverse events related to clonal evolution and heterogeneity.

5. Conclusions

In summary, we present a pragmatic approach to select risk-defined ALL subgroups based on a simple categorization of CNA profiles. Integrating CNA data with the risk index by primary aberrations, allowed us to stratify the outcome of about 70% of adult ALL: into the poor-risk (53%), and favorable-risk populations (14%). This resulted in the reclassification of an additional 37% of patients who changed their status from IM-risk patients to favorable- or poor-risk, which may result in more relevant treatment stratification (see Figure 7).
Importantly, specific CNA profiles correlated with different AID/RAG2 signatures. We propose that these signatures contribute to the high-risk behavior of selected ALL subsets: RAG as an active driver increasing the odds for CNAhigh/IKZF1pos development among CNApos, while the abundance of AID indicates poor risk subset among CNAneg ALL patients. This may be helpful to refine the treatment of CNAneg patients. These data may also encourage the notion of risk-reducing antimutagenic interventions in the course of treatment (e.g., reinforcing anti-inflammatory background).
Given the limitations of this study: its retrospective nature, the lack of a comprehensive SNV mutational profile at diagnosis, and relatively small subgroups, as well as scant data on mutational processes in adult ALL compared to the vast data in childhood leukemias—further investigations in this field are warranted.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15225431/s1, Table S1: Characteristics of detected markers or surrogates of BCR::ABL1-like B-ALL; Table S2: (A) AID/RAG expression signatures in correlation with the primary aberrations; (B) AID/RAG expression signatures in correlation with the distribution of CNA profiles; Table S3: (A) Frequency of CNA mutation burden and specific gene deletions in genetic subgroups of the studied cohort; (B) Correlation between CNA mutation burden and deletions in each gene; Table S4: Univariate analysis of patients’ outcome in subgroups according to CNA status; Table S5: (A) Frequency of patients with relapse and those remaining in remission—in subgroups according to mutator enzyme expression; (B) Univariate analysis of CNAneg patients’ outcomes according to mutator enzyme expression; (C) Multivariate analysis of CNAneg patients’ outcomes according to mutator enzyme expression, adjusted for age and WBC; Figure S1: Outcome of IKZF1-positive patients according to the number of CNAs. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in IKZF1-positive patients according to the number of copy number alterations (CNAs). CNAhigh—≥2 detected CNAs. (A,B)—overall survival and relapse-free survival in all IKZF1-mutated patients. (C,D)—overall survival and relapse-free survival in BCR::ABL1-positive IKZF1-positive patients only. n—number of patients, p—p-value; Figure S2: Outcome of CNAhigh patients according to IKZF1 mutation status. Kaplan-Meier estimates for the probability of (A) overall survival and (B) relapse-free survival in CNAhigh patients according to IKZF1 mutation status. CNAhigh—≥2 detected CNAs. n—number of patients, pp-value; Figure S3: Outcome of CNAhigh patients according to IKZF1 and CDKN2A/B mutation status. Kaplan-Meier estimates for the probability of (A) overall survival and (B) relapse-free survival in CNAhigh CDKN2A/Bpos patients according to IKZF1 mutation status. The graph only contains CDKN2A/Bpos patients. n—number of patients, pp-value; Figure S4: Outcome of CNAneg patients according to RAG2 and AID expression. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in CNAneg patients according to the AID and RAG2 expression levels. (A,B)—overall survival and relapse-free survival in CNAneg patients according to the AID expression level. (C,D)—overall survival and relapse-free survival in CNAneg patients according to the RAG2 expression level. CNAneg—0 detected copy number alterations, n—number of patients, pp-value; Figure S5: Outcome of CNApos patients according to RAG2 and AID expression. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in CNApos patients according to the AID and RAG2 expression levels. (A,B)—overall survival and relapse-free survival in NEG CNApos patients according to the AID expression level. (C,D)—overall survival and relapse-free survival in BCR::ABL1pos CNApos patients according to the RAG2 and AID expression levels. CNApos—at least 1 detected copy number alteration, n—number of patients, pp-value; Figure S6: Outcome of younger and older patients according to CNA profile. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in all studied patients according to the CNA profile, divided into younger (≤40 y.) and older (>40 y.) patients. (A,B)—overall survival and relapse-free survival in younger patients, grouped into CNAhigh / IKZF1pos and all other CNA profiles. (C,D)—overall survival and relapse-free survival in older patients, grouped into CNAhigh / IKZF1pos and all other CNA profiles. CNAhigh—at least 2 detected copy number alteration, n—number of patients, pp-value. References [44,46,67,68,69,70,71,72] are cited in the supplementary materials.

Author Contributions

Project design and supervision: M.L., K.K., J.H., E.L.-M., P.J., S.G. and W.W.J.; Genetic and laboratory testing: K.K., H.M.-Ł., B.J., M.Z. (Magdalena Zawada), M.G., J.T., K.B., M.W., S.C., I.F., K.L., M.Z. (Marta Zacharczuk), A.P. (Anna Pastwińska), M.K., M.P.-K. and O.H.; Patient management: A.C., M.P., B.P.-J., L.G., K.K.-S., E.P., A.T., K.M., R.M., T.C., A.P. (Agnieszka Piekarska), M.D., A.K., G.H., T.G., T.W., G.B., J.H., E.L.-M., S.G. and W.W.J.; Analysis and interpretation of clinical data: M.L., K.K., A.C., J.H. and S.G.; Performance and/or supervision of statistical studies: M.L., P.T. and S.G.; Manuscript writing M.L., K.K., P.T., A.C., M.S., R.M. and W.W.J.; Construction of manuscript and critical revision: A.C., O.H., E.L.-M., P.J., S.G. and W.W.J. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Science Centre in Poland (Narodowe Centrum Nauki), grant number 2015/19/N/NZ5/03495.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Bioethics Committee of the Medical University of Warsaw for all participating institutions (KB/69/2013, 19 March 2013).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roberts, K.G.; Mullighan, C.G. Genomics in Acute Lymphoblastic Leukaemia: Insights and Treatment Implications. Nat. Rev. Clin. Oncol. 2015, 12, 344–357. [Google Scholar] [CrossRef]
  2. Moorman, A.V.; Barretta, E.; Butler, E.R.; Ward, E.J.; Twentyman, K.; Kirkwood, A.A.; Enshaei, A.; Schwab, C.; Creasey, T.; Leongamornlert, D.; et al. Prognostic Impact of Chromosomal Abnormalities and Copy Number Alterations in Adult B-Cell Precursor Acute Lymphoblastic Leukaemia: A UKALL14 Study. Leukemia 2022, 36, 625–636. [Google Scholar] [CrossRef] [PubMed]
  3. Boer, J.M.; van der Veer, A.; Rizopoulos, D.; Fiocco, M.; Sonneveld, E.; de Groot-Kruseman, H.A.; Kuiper, R.P.; Hoogerbrugge, P.; Horstmann, M.; Zaliova, M.; et al. Prognostic Value of Rare IKZF1 Deletion in Childhood B-Cell Precursor Acute Lymphoblastic Leukemia: An International Collaborative Study. Leukemia 2016, 30, 32–38. [Google Scholar] [CrossRef] [PubMed]
  4. Beldjord, K.; Chevret, S.; Asnafi, V.; Huguet, F.; Boulland, M.L.; Leguay, T.; Thomas, X.; Cayuela, J.M.; Grardel, N.; Chalandon, Y.; et al. Oncogenetics and Minimal Residual Disease Are Independent Outcome Predictors in Adult Patients with Acute Lymphoblastic Leukemia. Blood 2014, 123, 3739–3749. [Google Scholar] [CrossRef]
  5. Moorman, A.V.; Enshaei, A.; Schwab, C.; Wade, R.; Chilton, L.; Elliott, A.; Richardson, S.; Hancock, J.; Kinsey, S.E.; Mitchell, C.D.; et al. A Novel Integrated Cytogenetic and Genomic Classification Refines Risk Stratification in Pediatric Acute Lymphoblastic Leukemia. Blood 2014, 124, 1434–1444. [Google Scholar] [CrossRef]
  6. Iacobucci, I.; Mullighan, C.G. Genetic Basis of Acute Lymphoblastic Leukemia. J. Clin. Oncol. 2017, 35, 975–983. [Google Scholar] [CrossRef] [PubMed]
  7. Creasey, T.; Barretta, E.; Ryan, S.L.; Butler, E.; Kirkwood, A.A.; Leongamornlert, D.; Papaemmanuil, E.; Patrick, P.; Clifton-Hadley, L.; Patel, B.; et al. Genetic and Genomic Analysis of Acute Lymphoblastic Leukemia in Older Adults Reveals a Distinct Profile of Abnormalities: Analysis of 210 Patients from the UKALL14 and UKALL60+ Clinical Trials. Haematologica 2022, 107, 2051–2063. [Google Scholar] [CrossRef] [PubMed]
  8. Jabbour, E.J.; Faderl, S.; Kantarjian, H.M. Adult Acute Lymphoblastic Leukemia. In Mayo Clinic Proceedings; Elsevier Ltd.: Amsterdam, The Netherlands, 2005; Volume 80, pp. 1517–1527. [Google Scholar]
  9. Larson, S.; Stock, W. Progress in the Treatment of Adults with Acute Lymphoblastic Leukemia. Curr. Opin. Hematol. 2008, 15, 400–407. [Google Scholar] [CrossRef]
  10. Roberts, K.G. Genetics and Prognosis of ALL in Children vs. Adults. Hematol. Am. Soc. Hematol. Educ. Program 2018, 201, 137–145. [Google Scholar] [CrossRef]
  11. Jabbour, E.; Short, N.J.; Jain, N.; Haddad, F.G.; Welch, M.A.; Ravandi, F.; Kantarjian, H. The Evolution of Acute Lymphoblastic Leukemia Research and Therapy at MD Anderson over Four Decades. J. Hematol. Oncol. 2023, 16, 22. [Google Scholar] [CrossRef]
  12. Sun, C.; Chang, L.; Zhu, X. Pathogenesis of ETV6/RUNX1-Positive Childhood Acute Lymphoblastic Leukemia and Mechanisms Underlying Its Relapse. Oncotarget 2017, 8, 35445–35459. [Google Scholar] [CrossRef] [PubMed]
  13. Swaminathan, S.; Klemm, L.; Park, E.; Papaemmanuil, E.; Ford, A.; Kweon, S.M.; Trageser, D.; Hasselfeld, B.; Henke, N.; Mooster, J.; et al. Mechanisms of Clonal Evolution in Childhood Acute Lymphoblastic Leukemia. Nat. Immunol. 2015, 16, 766–774. [Google Scholar] [CrossRef] [PubMed]
  14. Notta, F.; Mullighan, C.G.; Wang, J.C.Y.; Poeppl, A.; Doulatov, S.; Phillips, L.A.; Ma, J.; Minden, M.D.; Downing, J.R.; Dick, J.E. Evolution of Human BCR-ABL1 Lymphoblastic Leukaemia-Initiating Cells. Nature 2011, 469, 362–367. [Google Scholar] [CrossRef] [PubMed]
  15. Moorman, A.V. New and Emerging Prognostic and Predictive Genetic Biomarkers in B-Cell Precursor Acute Lymphoblastic Leukemia. Haematologica 2016, 101, 407–416. [Google Scholar] [CrossRef] [PubMed]
  16. Mullighan, C.G.; Goorha, S.; Radtke, I.; Miller, C.B.; Coustan-Smith, E.; Dalton, J.D.; Girtman, K.; Mathew, S.; Ma, J.; Pounds, S.B.; et al. Genome-Wide Analysis of Genetic Alterations in Acute Lymphoblastic Leukaemia. Nature 2007, 446, 758–764. [Google Scholar] [CrossRef] [PubMed]
  17. Iacobucci, I.; Storlazzi, C.T.; Cilloni, D.; Lonetti, A.; Ottaviani, E.; Soverini, S.; Astolfi, A.; Chiaretti, S.; Vitale, A.; Messa, F.; et al. Identification and Molecular Characterization of Recurrent Genomic Deletions on 7p12 in the IKZF1 Gene in a Large Cohort of BCR-ABL1-Positive Acute Lymphoblastic Leukemia Patients: On Behalf of Gruppo Italiano Malattie Ematologiche Dell’Adulto Acute Leukemia Working Party (GIMEMA AL WP). Blood 2009, 114, 2159–2167. [Google Scholar] [CrossRef]
  18. Anderson, K.; Lutz, C.; Van Delft, F.W.; Bateman, C.M.; Guo, Y.; Colman, S.M.; Kempski, H.; Moorman, A.V.; Titley, I.; Swansbury, J.; et al. Genetic Variegation of Clonal Architecture and Propagating Cells in Leukaemia. Nature 2011, 469, 356–361. [Google Scholar] [CrossRef]
  19. Patel, S.; Mason, C.C.; Glenn, M.J.; Paxton, C.N.; South, S.T.; Cessna, M.H.; Asch, J.; Cobain, E.F.; Bixby, D.L.; Smith, L.B.; et al. Genomic Analysis of Adult B-ALL Identifies Potential Markers of Shorter Survival. Leuk. Res. 2017, 56, 44–51. [Google Scholar] [CrossRef]
  20. Levy, G.; Kicinski, M.; Van der Straeten, J.; Uyttebroeck, A.; Ferster, A.; De Moerloose, B.; Dresse, M.F.; Chantrain, C.; Brichard, B.; Bakkus, M. Immunoglobulin Heavy Chain High-Throughput Sequencing in Pediatric B-Precursor Acute Lymphoblastic Leukemia: Is the Clonality of the Disease at Diagnosis Related to Its Prognosis? Front. Pediatr. 2022, 10, 874771. [Google Scholar] [CrossRef]
  21. Yang, J.J.; Bhojwani, D.; Yang, W.; Cai, X.; Stocco, G.; Crews, K.; Wang, J.; Morrison, D.; Devidas, M.; Hunger, S.P.; et al. Genome-Wide Copy Number Profiling Reveals Molecular Evolution from Diagnosis to Relapse in Childhood Acute Lymphoblastic Leukemia. Blood 2008, 112, 4178–4183. [Google Scholar] [CrossRef]
  22. Papaemmanuil, E.; Rapado, I.; Li, Y.; Potter, N.E.; Wedge, D.C.; Tubio, J.; Alexandrov, L.B.; Van Loo, P.; Cooke, S.L.; Marshall, J.; et al. RAG-Mediated Recombination Is the Predominant Driver of Oncogenic Rearrangement in ETV6-RUNX1 Acute Lymphoblastic Leukemia. Nat. Genet. 2014, 46, 116–125. [Google Scholar] [CrossRef] [PubMed]
  23. Muramatsu, M.; Sankaranand, V.S.; Anant, S.; Sugai, M.; Kinoshita, K.; Davidson, N.O.; Honjo, T. Specific Expression of Activation-Induced Cytidine Deaminase (AID), a Novel Member of the RNA-Editing Deaminase Family in Germinal Center B Cells. J. Biol. Chem. 1999, 274, 18470–18476. [Google Scholar] [CrossRef] [PubMed]
  24. Maul, R.W.; Gearhart, P.J. AID and Somatic Hypermutation. Adv. Immunol. 2010, 105, 159–191. [Google Scholar] [CrossRef] [PubMed]
  25. Gruber, T.A.; Chang, M.S.; Sposto, R.; Müschen, M. Activation-Induced Cytidine Deaminase Accelerates Clonal Evolution in BCR-ABL1-Driven B-Cell Lineage Acute Lymphoblastic Leukemia. Cancer Res. 2010, 70, 7411–7420. [Google Scholar] [CrossRef] [PubMed]
  26. Tsai, A.G.; Lu, H.; Raghavan, S.C.; Muschen, M.; Hsieh, C.L.; Lieber, M.R. Human Chromosomal Translocations at CpG Sites and a Theoretical Basis for Their Lineage and Stage Specificity. Cell 2008, 135, 1130–1142. [Google Scholar] [CrossRef]
  27. Zhang, P.; Qin, M.; Wang, Y.; Chen, X.; Miao, Y.; Yuan, M.; Zhou, W.; Li, D.; Wang, D.; Wang, M.; et al. Inflammation Accelerates BCR-ABL11 B-ALL Development through Upregulation of AID. Blood Adv. 2022, 6, 4060–4072. [Google Scholar] [CrossRef]
  28. Dong, Y.; Liu, F.; Wu, C.; Li, S.; Zhao, X.; Zhang, P.; Jiao, J.; Yu, X.; Ji, Y.; Zhang, M. Illegitimate RAG-Mediated Recombination Events Are Involved in IKZF1 Delta 3–6 Deletion in BCR-ABL1 Lymphoblastic Leukaemia. Clin. Exp. Immunol. 2016, 185, 320–331. [Google Scholar] [CrossRef]
  29. Mullighan, C.G.; Miller, C.B.; Radtke, I.; Phillips, L.A.; Dalton, J.; Ma, J.; White, D.; Hughes, T.P.; Le Beau, M.M.; Pui, C.H.; et al. BCR-ABL1 Lymphoblastic Leukaemia Is Characterized by the Deletion of Ikaros. Nature 2008, 453, 110–114. [Google Scholar] [CrossRef]
  30. Sanders, M.A.; Szabó, A.; Exalto, C.; Hoogenboezem, R.; Zeilemaker, A.; Koenders, J.E.; van Geel, P.; Schelen, A.; Beverloo, H.B.; Cornelissen, J.J.; et al. Extensive RAG-Mediated Rearrangements and Mutations in BCR-ABL1 and BCR-ABL1-like Adult Acute Lymphoblastic Leukemia. Blood 2016, 128, 4067. [Google Scholar] [CrossRef]
  31. Kim, J.C.; Chan-Seng-Yue, M.; Ge, S.; Zeng, A.G.X.; Ng, K.; Gan, O.I.; Garcia-Prat, L.; Flores-Figueroa, E.; Woo, T.; Zhang, A.X.W.; et al. Transcriptomic Classes of BCR-ABL1 Lymphoblastic Leukemia. Nat. Genet. 2023, 55, 1186–1197. [Google Scholar] [CrossRef]
  32. Mullighan, C.G.; Su, X.; Zhang, J.; Radtke, I.; Phillips, L.A.A.; Miller, C.B.; Ma, J.; Liu, W.; Cheng, C.; Schulman, B.A.; et al. Deletion of IKZF1 and Prognosis in Acute Lymphoblastic Leukemia. N. Engl. J. Med. 2009, 360, 470–480. [Google Scholar] [CrossRef] [PubMed]
  33. Ribera, J.; Morgades, M.; Zamora, L.; Montesinos, P.; Gõmez-Seguí, I.; Pratcorona, M.; Sarrà, J.; Guàrdia, R.; Nomdedeu, J.; Tormo, M.; et al. Prognostic Significance of Copy Number Alterations in Adolescent and Adult Patients with Precursor B Acute Lymphoblastic Leukemia Enrolled in PETHEMA Protocols. Cancer 2015, 121, 3809–3817. [Google Scholar] [CrossRef]
  34. Giebel, S.; Labopin, M.; Socié, G.; Beauvais, D.; Klein, S.; Wagner-Drouet, E.M.; Blaise, D.; Nguyen-Quoc, S.; Bourhis, J.H.; Thiebaut, A.; et al. Allogeneic HCT for Adults with B-Cell Precursor Acute Lymphoblastic Leukemia Harboring IKZF1 Gene Mutations. A Study by the Acute Leukemia Working Party of the EBMT. Bone Marrow Transpl. 2021, 56, 1047–1055. [Google Scholar] [CrossRef] [PubMed]
  35. Paietta, E.; Roberts, K.G.; Wang, V.; Gu, Z.; Buck, G.A.N.; Pei, D.; Cheng, C.; Levine, R.L.; Abdel-Wahab, O.; Cheng, Z.; et al. Molecular Classification Improves Risk Assessment in Adult BCR-ABL1– Negative B-ALL. Blood 2021, 138, 948–958. [Google Scholar] [CrossRef] [PubMed]
  36. Messina, M.; Chiaretti, S.; Fedullo, A.L.; Piciocchi, A.; Puzzolo, M.C.; Lauretti, A.; Gianfelici, V.; Apicella, V.; Fazi, P.; Te Kronnie, G.; et al. Clinical Significance of Recurrent Copy Number Aberrations in B-Lineage Acute Lymphoblastic Leukaemia without Recurrent Fusion Genes across Age Cohorts. Br. J. Haematol. 2017, 178, 583–587. [Google Scholar] [CrossRef]
  37. Duffield, A.S.; Mullighan, C.G.; Borowitz, M.J. International Consensus Classification of Acute Lymphoblastic Leukemia/Lymphoma. Virchows Arch. 2023, 482, 11–26. [Google Scholar] [CrossRef]
  38. Alaggio, R.; Amador, C.; Anagnostopoulos, I.; Attygalle, A.D.; de Araujo, I.B.O.; Berti, E.; Bhagat, G.; Borges, A.M.; Boyer, D.; Calaminici, M.; et al. The 5th Edition of the World Health Organization Classification of Haematolymphoid Tumours: Lymphoid Neoplasms. Leukemia 2022, 36, 1720–1748. [Google Scholar] [CrossRef]
  39. Stanulla, M.; Dagdan, E.; Zaliova, M.; Möricke, A.; Palmi, C.; Cazzaniga, G.; Eckert, C.; Te Kronnie, G.; Bourquin, J.-P.; Koehler, R.; et al. IKZF1 plus Defines a New Minimal Residual Disease-Dependent Very-Poor Prognostic Profile in Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia. J. Clin. Oncol. 2018, 36, 1240–1249. [Google Scholar] [CrossRef]
  40. Gupta, S.K.; Singh, M.; Chandrashekar, P.H.; Bakhshi, S.; Trehan, A.; Gupta, R.; Thakur, R.; Gajendra, S.; Sharma, P.; Sreedharanunni, S.; et al. Clinical and Prognostic Impact of Copy Number Alterations and Associated Risk Profiles in a Cohort of Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia Cases Treated under ICiCLe Protocol. Hemasphere 2022, 6, E782. [Google Scholar] [CrossRef]
  41. van Dongen, J.; Macintyre, E.; Gabert, J.; Delabesse, E.; Rossi, V.; Saglio, G.; Gottardi, E.; Rambaldi, A.; Dotti, G.; Griesinger, F.; et al. Standardized RT-PCR Analysis of Fusion Gene Transcripts from Chromosome Aberrations in Acute Leukemia for Detection of Minimal Residual Disease Report of the BIOMED-1 Concerted Action: Investigation of Minimal Residual Disease in Acute Leukemia. Leukemia 1999, 13, 1901–1928. [Google Scholar] [CrossRef]
  42. Soszynska, K.; Mucha, B.; Debski, R.; Skonieczka, K.; Duszenko, E.; Koltan, A.; Wysocki, M.; Haus, O. The Application of Conventional Cytogenetics, FISH, and RT-PCR to Detect Genetic Changes in 70 Children with ALL. Ann. Hematol. 2008, 87, 991–1002. [Google Scholar] [CrossRef] [PubMed]
  43. Oettinger, M.A.; Schatz, D.G.; Gorka, C.; Baltimore, D. RAG-1 and RAG-2, Adjacent Genes That Synergistically Activate V(D)J Recombination. Science 1990, 248, 1517–1523. [Google Scholar] [CrossRef] [PubMed]
  44. Iacobucci, I.; Lonetti, A.; Messa, F.; Cilloni, D.; Arruga, F.; Ottaviani, E.; Paolini, S.; Papayannidis, C.; Piccaluga, P.P.; Giannoulia, P.; et al. Expression of Spliced Oncogenic Ikaros Isoforms in Philadelphia-Positive Acute Lymphoblastic Leukemia Patients Treated with Tyrosine Kinase Inhibitors: Implications for a New Mechanism of Resistance. Blood 2008, 112, 3847–3855. [Google Scholar] [CrossRef] [PubMed]
  45. Huguet, F.; Chevret, S.; Leguay, T.; Thomas, X.; Boissel, N.; Escoffre-Barbe, M.; Chevallier, P.; Hunault, M.; Vey, N.; Bonmati, C.; et al. Intensified Therapy of Acute Lymphoblastic Leukemia in Adults: Report of the Randomized GRAALL-2005 Clinical Trial. J. Clin. Oncol. 2018, 36, 2514–2523. [Google Scholar] [CrossRef] [PubMed]
  46. Boer, J.M.; Koenders, J.E.; Van Der Holt, B.; Exalto, C.; Sanders, M.A.; Cornelissen, J.J.; Valk, P.J.M.; Den Boer, M.L.; Rijneveld, A.W. Expression Profiling of Adult Acute Lymphoblastic Leukemia Identifies a BCR-ABL1-like Subgroup Characterized by High Non-Response and Relapse Rates. Haematologica 2015, 100, 261–264. [Google Scholar] [CrossRef]
  47. Pfeifer, H.; Raum, K.; Markovic, S.; Nowak, V.; Fey, S.; Obländer, J.; Pressler, J.; Böhm, V.; Brüggemann, M.; Wunderle, L.; et al. Genomic CDKN2A/2B Deletions in Adult Ph1 All Are Adverse despite Allogeneic Stem Cell Transplantation. Blood 2018, 131, 1464–1475. [Google Scholar] [CrossRef]
  48. Zhang, W.; Kuang, P.; Liu, T. Prognostic Significance of CDKN2A/B Deletions in Acute Lymphoblastic Leukaemia: A Meta-Analysis. Ann. Med. 2019, 51, 28–40. [Google Scholar] [CrossRef]
  49. Jabbour, E.; Short, N.J.; Ravandi, F.; Huang, X.; Daver, N.; DiNardo, C.D.; Konopleva, M.; Pemmaraju, N.; Wierda, W.; Garcia-Manero, G.; et al. Combination of Hyper-CVAD with Ponatinib as First-Line Therapy for Patients with Philadelphia Chromosome-Positive Acute Lymphoblastic Leukaemia: Long-Term Follow-up of a Single-Centre, Phase 2 Study. Lancet Haematol. 2018, 5, 618–627. [Google Scholar] [CrossRef]
  50. Ravandi, F.; O’Brien, S.M.; Cortes, J.E.; Thomas, D.M.; Garris, R.; Faderl, S.; Burger, J.A.; Rytting, M.E.; Ferrajoli, A.; Wierda, W.G.; et al. Long-Term Follow-up of a Phase 2 Study of Chemotherapy plus Dasatinib for the Initial Treatment of Patients with Philadelphia Chromosome-Positive Acute Lymphoblastic Leukemia. Cancer 2015, 121, 4158–4164. [Google Scholar] [CrossRef]
  51. Bassan, R.; Rossi, G.; Pogliani, E.M.; Di Bona, E.; Angelucci, E.; Cavattoni, I.; Lambertenghi-Deliliers, G.; Mannelli, F.; Levis, A.; Ciceri, F.; et al. Chemotherapy-Phased Imatinib Pulses Improve Long-Term Outcome of Adult Patients with Philadelphia Chromosome-Positive Acute Lymphoblastic Leukemia: Northern Italy Leukemia Group Protocol 09/00. J. Clin. Oncol. 2010, 28, 3644–3652. [Google Scholar] [CrossRef]
  52. Daver, N.; Thomas, D.; Ravandi, F.; Cortes, J.; Garris, R.; Jabbour, E.; Garcia-Manero, G.; Borthakur, G.; Kadia, T.; Rytting, M.; et al. Final Report of a Phase II Study of Imatinib Mesylate with Hyper-CVAD for the Front-Line Treatment of Adult Patients with Philadelphia Chromosome-Positive Acute Lymphoblastic Leukemia. Haematologica 2015, 100, 653–661. [Google Scholar] [CrossRef]
  53. Short, N.J.; Kantarjian, H.M.; Ravandi, F.; Huang, X.; Daver, N.G.; DiNardo, C.D.; Konopleva, M.Y.; Pemmaraju, N.; Wierda, W.G.; Garcia-Manero, G.; et al. Long-Term Safety and Efficacy of Hyper-CVAD Plus Ponatinib As Frontline Therapy for Adults with Philadelphia Chromosome-Positive Acute Lymphoblastic Leukemia. Blood 2019, 134, 283. [Google Scholar] [CrossRef]
  54. Holowiecki, J.; Krawczyk-Kulis, M.; Giebel, S.; Jagoda, K.; Stella-Holowiecka, B.; Piatkowska-Jakubas, B.; Paluszewska, M.; Seferynska, I.; Lewandowski, K.; Kielbinski, M.; et al. Status of Minimal Residual Disease after Induction Predicts Outcome in Both Standard and High-Risk Ph-Negative Adult Acute Lymphoblastic Leukaemia. The Polish Adult Leukemia Group ALL 4-2002 MRD Study. Br. J. Haematol. 2008, 142, 227–237. [Google Scholar] [CrossRef] [PubMed]
  55. Zaliova, M.; Stuchly, J.; Winkowska, L.; Musilova, A.; Fiser, K.; Slamova, M.; Starkova, J.; Vaskova, M.; Hrusak, O.; Sramkova, L.; et al. Genomic Landscape of Pediatric B-Other Acute Lymphoblastic Leukemia in a Consecutive European Cohort. Haematologica 2019, 104, 1396–1406. [Google Scholar] [CrossRef]
  56. Han, Q.; Ma, J.; Gu, Y.; Song, H.; Kapadia, M.; Kawasawa, Y.I.; Dovat, S.; Song, C.; Ge, Z. RAG1 High Expression Associated with IKZF1 Dysfunction in Adult B-Cell Acute Lymphoblastic Leukemia. J. Cancer 2019, 10, 3842–3850. [Google Scholar] [CrossRef] [PubMed]
  57. Heinä niemi, M.; Vuorenmaa, T.; Teppo, S.; Kaikkonen, M.U.; Bouvy-Liivrand, M.; Mehtonen, J.; Niskanen, H.; Zachariadis, V.; Laukkanen, S.; Liuksiala, T.; et al. Transcription-Coupled Genetic Instability Marks Acute Lymphoblastic Leukemia Structural Variation Hotspots. eLife 2016, 6, e13087. [Google Scholar] [CrossRef]
  58. Peppas, I.; Ford, A.M.; Furness, C.L.; Greaves, M.F. Gut Microbiome Immaturity and Childhood Acute Lymphoblastic Leukaemia. Nat. Rev. Cancer 2023, 23, 565–576. [Google Scholar] [CrossRef]
  59. Vicente-Dueñas, C.; Janssen, S.; Oldenburg, M.; Auer, F.; González-Herrero, I.; Casado-García, A.; Isidro-Hernández, M.; Raboso-Gallego, J.; Westhoff, P.; Pandyra, A.A.; et al. An Intact Gut Microbiome Protects Genetically Predisposed Mice against Leukemia. Blood 2020, 136, 2003–2017. [Google Scholar] [CrossRef]
  60. Greaves, M. A Causal Mechanism for Childhood Acute Lymphoblastic Leukaemia. Nat. Rev. Cancer 2018, 18, 471–484. [Google Scholar] [CrossRef]
  61. Çakan, E.; Gunaydin, G. Activation Induced Cytidine Deaminase: An Old Friend with New Faces. Front. Immunol. 2022, 13, 965312. [Google Scholar] [CrossRef] [PubMed]
  62. Shimizu, T.; Marusawa, H.; Endo, Y.; Chiba, T. Inflammation-Mediated Genomic Instability: Roles of Activation-Induced Cytidine Deaminase in Carcinogenesis. Cancer Sci. 2012, 103, 1201–1206. [Google Scholar] [CrossRef]
  63. Kou, T.; Marusawa, H.; Kinoshita, K.; Endo, Y.; Okazaki, I.M.; Ueda, Y.; Kodama, Y.; Haga, H.; Ikai, I.; Chiba, T. Expression of Activation-Induced Cytidine Deaminase in Human Hepatocytes during Hepatocarcinogenesis. Int. J. Cancer 2007, 120, 469–476. [Google Scholar] [CrossRef]
  64. Fang, Q.; Song, Y.; Gong, X.; Wang, J.; Li, Q.; Liu, K.; Feng, Y.; Hao, Q.; Li, Y.; Wei, H.; et al. Gene Deletions and Prognostic Values in B-Linage Acute Lymphoblastic Leukemia. Front. Oncol. 2021, 11, 677034. [Google Scholar] [CrossRef]
  65. Churchman, M.L.; Low, J.; Qu, C.; Paietta, E.M.; Kasper, L.H.; Chang, Y.; Payne-Turner, D.; Althoff, M.J.; Song, G.; Chen, S.C.; et al. Efficacy of Retinoids in IKZF1-Mutated BCR-ABL1 Acute Lymphoblastic Leukemia. Cancer Cell 2015, 28, 343–356. [Google Scholar] [CrossRef]
  66. Koptyra, M.; Falinski, R.; Nowicki, M.O.; Stoklosa, T.; Majsterek, I.; Nieborowska-Skorska, M.; Blasiak, J.; Skorski, T. BCR/ABL Kinase Induces Self-Mutagenesis via Reactive Oxygen Species to Encode Imatinib Resistance. Blood 2006, 108, 319–327. [Google Scholar] [CrossRef]
  67. Shochat, C.; Tal, N.; Bandapalli, O.R.; Palmi, C.; Ganmore, I.; te Kronnie, G.; Cario, G.; Cazzaniga, G.; Kulozik, A.E.; Stanulla, M.; et al. Gain-of-Function Mutations in Interleukin-7 Receptor-α (IL7R) in Childhood Acute Lymphoblastic Leukemias. J. Exp. Med. 2011, 208, 901–908. [Google Scholar] [CrossRef]
  68. Holmfeldt, L.; Wei, L.; Diaz-Flores, E.; Walsh, M.; Zhang, J.; Ding, L.; Payne-Turner, D.; Churchman, M.; Andersson, A.; Chen, S.C.; et al. The Genomic Landscape of Hypodiploid Acute Lymphoblastic Leukemia. Nat. Genet. 2013, 45, 242–252. [Google Scholar] [CrossRef]
  69. Iacobucci, I.; Iraci, N.; Messina, M.; Lonetti, A.; Chiaretti, S.; Valli, E.; Ferrari, A.; Papayannidis, C.; Paoloni, F.; Vitale, A.; et al. IKAROS Deletions Dictate a Unique Gene Expression Signature in Patients with Adult B-Cell Acute Lymphoblastic Leukemia. PLoS ONE 2012, 7, e40934. [Google Scholar] [CrossRef]
  70. Hertzberg, L.; Vendramini, E.; Ganmore, I.; Cazzaniga, G.; Schmitz, M.; Chalker, J.; Shiloh, R.; Iacobucci, I.; Shochat, C.; Zeligson, S.; et al. Down Syndrome Acute Lymphoblastic Leukemia, a Highly Heterogeneous Disease in Which Aberrant Expression of CRLF2 Is Associated with Mutated JAK2: A Report from the International BFM Study Group. Blood 2010, 115, 1006–1017. [Google Scholar] [CrossRef]
  71. Rao, X.; Huang, X.; Zhou, Z.; Lin, X. An Improvement of the 2ˆ(−Delta Delta CT) Method for Quantitative Real-Time Polymerase Chain Reaction Data Analysis. Biostat. Bioinform. Biomath. 2013, 3, 71–85. [Google Scholar]
  72. Kwon, M.; Martínez-Laperche, C.; Infante, M.; Carretero, F.; Balsalobre, P.; Serrano, D.; Gayoso, J.; Pérez-Corral, A.; Anguita, J.; Díez-Martín, J.L.; et al. Evaluation of Minimal Residual Disease by Real-Time Quantitative PCR of Wilms’ Tumor 1 Expression in Patients with Acute Myelogenous Leukemia after Allogeneic Stem Cell Transplantation: Correlation with Flow Cytometry and Chimerism. Biol. Blood Marrow Transplant. 2012, 18, 1235–1242. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The frequency of primary and secondary aberrations in the patients’ cohort. (A) The frequency of primary aberrations in the patients’ cohort. The tables denote the frequencies of copy number alteration (CNA) mutation burden in the respective primary aberrations. (B) The frequency of genes with detected CNAs in the studied cohort. Some patients harbor more than one CNA, hence the sum exceeds 100%.
Figure 1. The frequency of primary and secondary aberrations in the patients’ cohort. (A) The frequency of primary aberrations in the patients’ cohort. The tables denote the frequencies of copy number alteration (CNA) mutation burden in the respective primary aberrations. (B) The frequency of genes with detected CNAs in the studied cohort. Some patients harbor more than one CNA, hence the sum exceeds 100%.
Cancers 15 05431 g001
Figure 2. Patients’ outcomes according to the primary aberrations. (A,B) Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the primary aberration. (C,D) Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the revised code of the primary aberration. n—number of patients, pp-value.
Figure 2. Patients’ outcomes according to the primary aberrations. (A,B) Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the primary aberration. (C,D) Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the revised code of the primary aberration. n—number of patients, pp-value.
Cancers 15 05431 g002
Figure 3. Patients’ outcomes according to the copy number alteration mutation burden. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the number of copy number alterations (CNAs). CNAlow—0–1 detected CNAs, CNAhigh—≥2 detected CNAs. (A,B)—Overall survival and relapse-free survival in all studied patients. (C,D)—Overall survival and relapse-free survival in NEG patients only. n—number of patients, pp-value.
Figure 3. Patients’ outcomes according to the copy number alteration mutation burden. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the number of copy number alterations (CNAs). CNAlow—0–1 detected CNAs, CNAhigh—≥2 detected CNAs. (A,B)—Overall survival and relapse-free survival in all studied patients. (C,D)—Overall survival and relapse-free survival in NEG patients only. n—number of patients, pp-value.
Cancers 15 05431 g003
Figure 4. Patients’ outcome according to IKZF1 mutation presence. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to IKZF1 mutation status. (A,B)—Overall survival and relapse-free survival in all studied patients. (C,D)—Overall survival and relapse-free survival in NEG patients only. n—number of patients, pp-value.
Figure 4. Patients’ outcome according to IKZF1 mutation presence. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to IKZF1 mutation status. (A,B)—Overall survival and relapse-free survival in all studied patients. (C,D)—Overall survival and relapse-free survival in NEG patients only. n—number of patients, pp-value.
Cancers 15 05431 g004
Figure 5. Patients’ outcomes according to the revised risk stratification by CNA profile. CNAhigh/IKZF1pos represent bad-CNA profile, 1CNA/IKZF1pos or CNAhigh/IKZF1neg represent good-CNA profile. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the revised code of secondary aberrations. (A,B)—all studied patients. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0566; good-CNA vs. bad-CNA p = 0.0026. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0176; good-CNA vs. bad-CNA p = 0.0061. (C,D)—NEG patients only. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0212; good-CNA vs. bad-CNA p = 0.0055. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0063; good-CNA vs. bad-CNA p = 0.0048. (E,F)—BCR::ABL1-positive patients only. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0152; good-CNA vs. bad-CNA p = 0.0032. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0799; good-CNA vs. bad-CNA p = 0.0167. n—number of patients.
Figure 5. Patients’ outcomes according to the revised risk stratification by CNA profile. CNAhigh/IKZF1pos represent bad-CNA profile, 1CNA/IKZF1pos or CNAhigh/IKZF1neg represent good-CNA profile. Kaplan-Meier estimates for the probability of overall survival and relapse-free survival in studied patients according to the revised code of secondary aberrations. (A,B)—all studied patients. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0566; good-CNA vs. bad-CNA p = 0.0026. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0176; good-CNA vs. bad-CNA p = 0.0061. (C,D)—NEG patients only. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0212; good-CNA vs. bad-CNA p = 0.0055. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0063; good-CNA vs. bad-CNA p = 0.0048. (E,F)—BCR::ABL1-positive patients only. OS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0152; good-CNA vs. bad-CNA p = 0.0032. RFS: CNAneg vs. good-CNA vs. bad-CNA p = 0.0799; good-CNA vs. bad-CNA p = 0.0167. n—number of patients.
Cancers 15 05431 g005
Figure 6. Patients’ outcomes according to the revised risk stratification by primary and secondary aberrations. Kaplan-Meier estimates for the probability of (A) overall survival and (B) relapse-free survival in studied patients according to the revised risk stratification code including the primary and the secondary aberration. n—number of patients, pp-value.
Figure 6. Patients’ outcomes according to the revised risk stratification by primary and secondary aberrations. Kaplan-Meier estimates for the probability of (A) overall survival and (B) relapse-free survival in studied patients according to the revised risk stratification code including the primary and the secondary aberration. n—number of patients, pp-value.
Cancers 15 05431 g006
Figure 7. Schematic hierarchical definition of ALL genetic subgroups.
Figure 7. Schematic hierarchical definition of ALL genetic subgroups.
Cancers 15 05431 g007
Table 1. Patients’ characteristics according to the primary aberrations, secondary aberrations and mutator enzyme expression. (A) Patients’ characteristics according to the primary aberration. (B) Patients’ characteristics according to RAG2 and AID expression. (C) Patients’ characteristics according to the number of CNAs and IKZF1 mutation status.
Table 1. Patients’ characteristics according to the primary aberrations, secondary aberrations and mutator enzyme expression. (A) Patients’ characteristics according to the primary aberration. (B) Patients’ characteristics according to RAG2 and AID expression. (C) Patients’ characteristics according to the number of CNAs and IKZF1 mutation status.
(A)
TotalNEGBCR::ABL1“Bad Primary”E2A::PBX
BCR::ABL1-likeMLL::AF4ComplexHyperdiploidAll “Bad Primary”
n = 166n = 63
(38%)
n = 55
(33%)
n = 16
(10%)
n = 5
(3%)
n = 15
(9%)
n = 6
(4%)
n = 44
(27%)
n = 4
(2%)
Age
≤37 years (%)83 (50%)33 (52%)24 (44%)6 (38%)5 (100%)8 (53%)3 (50%)23 (52%)3 (75%)
>37 years (%)83 (50%)30 (48%)31 (56%)10 (62%)0 (0%)7 (47%)3 (50%)21 (48%)1 (25%)
Gender
Female (%)57 (34%)15 (24%)24 (44%)7 (44%)0 (0%)5 (33%)3 (50%)16 (36%)2 (50%)
Male (%)109 (66%)48 (76%)31 (56%)9 (56%)5 (100%)10 (67%)3 (50%)28 (64%)2 (50%)
Immunological subtype
prepreB (%)22 (13%)11 (17%)3 (5%)1 (6%)3 (60%)3 (21%)1 (17%)8 (19%)0 (0%)
preB (%)42 (25%)17 (27%)13 (24%)4 (25%)2 (40%)4 (29%)1 (17%)11 (26%)1 (25%)
common (%)101 (62%)35 (56%)39 (71%)11 (69%)0 (0%)7 (50%)4 (67%)24 (56%)3 (75%)
Median of WBC, ×109/L17.013.526.028.7127.9 #9.24.610.924.3
Rate of CD20 positive, %18.014.028.031.00.535.025.025.09.5
Rate of CD52 positive, %18.018.022.043.08.712.05.514.762.0
Risk group according to PALG
very high risk (%)55 (33%)0 (0%)55 (100%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)
high risk (%)94 (57%)54 (86%)0 (0%)16 (100%)5 (100%)13 (87%)3 (50%)38 (86%)2 (50%)
standard risk (%)17 (10%)9 (14%)0 (0%)0 (0%)0 (0%)2 (13%)3 (50%)6 (14%)2 (50%)
MRD status after first induction
positive (%)45 (38%)13 (30%)19 (53%)5 (56%)0 (0%)6 (43%)1 (20%)13 (37%)0 (0%)
negative (%)73 (62%)31 (70%)17 (47%)4 (44%)5 (100%)8 (57%)4 (80%)22 (63%)3 (100%)
Induction regimen
PALG5 (%)15 (9%)4 (6%)6 (11%)0 (0%)0 (0%)2 (13%)1 (17%)3 (7%)2 (50%)
PALG6 (%)145 (88%)58 (92%)46 (84%)15 (100%)5 (100%)13 (87%)5 (83%)40 (93%)1 (25%)
PALG modified (%)5 (3%)1 (2%)3 (5%)0 (0%)0 (0%)0 (0%)0 (0%)0 (0%)1 (25%)
AlloHSCT in first complete response
yes (%)77 (47%)27 (44%)29 (53%)7 (44%)5 (100%)6 (40%)1 (17%)19 (43%)2 (50%)
no (%)88 (53%)35 (56%)26 (47%)9 (56%)0 (0%)9 (60%)5 (83%)25 (57%)2 (50%)
(B)
TotalRAG2 ExpressionAID Expression
LowHighpLowHighp
n = 166n = 83n = 78 n = 78n = 79
Median age, years37.540.037.00.4817 *35.541.00.2802 *
Gender 0.5191 0.1440
Female (%)5728 (34%)27 (35%)31 (40%)24 (30%)
Male (%)10955 (66%)51 (65%)47 (60%)55 (70%)
Primary aberration
NEG6337 (62%)23 (38%)0.0478 a22 (38%)36 (62%)0.0242 a
BCR::ABL15521 (39%)33 (61%)0.0224 b30 (57%)23 (43%)0.2156 b
BCR::ABL1-like168 (53%)7 (47%) 7 (50%)7 (50%)
MLL::AF452 (40%)3 (60%) 5 (100%)0 (0%)
Complex158 (53%)7 (47%) 7 (47%)8 (53%)
Hyperdiploid85 (83%)1 (17%) 3 (50%)3 (50%)
E2A::PBX41 (25%)3 (75%) 3 (75%)1 (25%)
Immunological subtype 0.4641 c 0.0420 c
prepreB (%)2213 (59%)9 (41%)14 (67%)7 (33%)
preB (%)4218 (46%)21(54%)14 (36%)25 (64%)
common (%)10151 (52%)48 (48%)50 (52%)46 (48%)
Median of WBC, ×109/L17.010.223.00.0559 *21.012.90.0589 *
Rate of CD20 positive, %18.022.020.00.5278 *11.030.00.3524 *
Rate of CD52 positive, %18.022.517.50.5950 *22.517.00.1320 *
Risk group according to PALG 0.0681 0.1691
very high risk (%)5521 (39%)33 (61%)30 (57%)23 (43%)
high risk (%)9452 (57%)39 (43%)38 (43%)50 (57%)
standard risk (%)1710 (63%)6 (37%)10 (63%)6 (37%)
MRD status after first induction 0.4792 0.1463
positive (%)4524 (55%)20 (45%)18 (43%)24 (57%)
negative (%)7336 (52%)33 (48%)38 (55%)31 (45%)
Induction regimen
PALG5 (%)158 (57%)6 (43%)4 (29%)10 (71%)
PALG6 (%)14573 (52%)68 (48%)71 (52%)66 (48%)
PALG modified (%)51 (20%)4 (80%)3 (60%)2 (40%)
AlloHSCT in first complete response
yes (%)7736 (49%)38(51%)34 (47%)38 (53%)
no (%)8847 (54%)40 (46%)44 (52%)41 (48%)
(C)
TotalCNA PresenceCNA Mutation BurdenIKZF1 Mutation StatusIKZF1pos Patients Only d
CNAnegCNApospCNAlow (0–1 CNA)CNAhigh (≥2 CNAs)pIKZFnegIKZFpos ep1 CNACNAhighp
n = 94n = 28
(30%)
n = 66
(70%)
n = 53
(56%)
n = 41
(44%)
n = 53
(52%)
n = 49
(48%)
n = 12
(29%)
n = 29
(71%)
Median age, years37.033.039.50.0227 *36.039.00.1128 *35.041.00.0152 *41.045.00.6466 *
Gender 0.5383 0.3658 0.3789 0.5365
Female (%)3611 (31%)25 (69%)19 (53%)17 (47%)21 (58%)17 (42%)4 (27%)11 (73%)
Male (%)58 17 (29%)41 (71%)34 (59%)24 (41%)32 (50%)32 (50%)8 (31%)18 (69%)
Primary aberration
NEG35 (37%)11 (31%)24 (69%)0.7887 a20 (57%)15 (43%)0.9089 a21 (60%)14 (40%)0.2402 a4 (29%)10 (71%)0.9437 a
BCR::ABL130 (32%)6 (20%)24 (80%)0.1554 b14 (47%)16 (53%)0.1934 b10 (27%)27 (73%)0.0001 b6 (30%)14 (70%)0.9200 b
BCR::ABL1-like12 (13%)2 (17%)10 (83%) 5 (42%)7 (58%) 6 (50%)6 (50%) 2 (33%)4 (67%)
MLL::AF43 (3%)2 (67%)1 (33%) 3 (100%)0 (0%) 3 (100%)0 (0%) 0 (0%)0 (0%)
Complex6 (6%)2 (33%)4 (67%) 3 (50%)3 (50%) 5 (71%)2 (29%) 0 (0%)1 (100%)
Hyperdiploid6 (6%)3 (50%)3 (50%) 6 (100%)0 (0%) 6 (100%)0 (0%) 0 (0%)0 (0%)
E2A::PBX2 (2%)2 (100%)0 (0%) 2 (100%)0 (0%) 2 (100%)0 (0%) 0 (0%)0 (0%)
Immunological subtype 0.2555 c 0.2765 c 0.8579 c 0.9534 c
prepreB (%)115 (45%)6 (55%)7 (64%)4 (46%)9 (82%)2 (18%)0 (0%)2 (100%)
preB (%)2610 (38%)16 (62%)17 (65%)9 (35%)16 (53%)14 (47%)3 (30%)7 (70%)
common (%)5713 (23%)44 (77%)29 (51%)28 (49%)28 (46%)33 (54%)9 (31%)20 (69%)
Median of WBC, ×109/L15.513.717.90.8178 *13.623.60.0939 *13.819.00.3394 *7.224.40.1188 *
Rate of CD20 positive, %18.011.018.00.7008 *10.020.00.2409 *18.020.00.9968 *6.219.00.1939 *
Rate of CD52 positive, %22.515.026.20.2507 *26.021.50.4730 *22.018.00.6487 *82.018.00.0502 *
Risk group acc. to PALG 0.0567 0.0561 0.0005  0.7786
very high risk (%)306 (20%)24 (80%)14 (47%)16 (53%)10 (27%)27 (73%)6 (30%)14 (70%)
high risk (%)5416 (30%)38 (70%)30 (56%)24 (44%)35 (64%)20 (36%)5 (26%)14 (74%)
standard risk (%)106 (60%)4 (40%)9 (90%)1 (10%)8 (80%)2 (20%)1 (50%)1 (50%)
MRD after 1st induction 0.0486  0.0003  0.0019  0.0149 
positive (%)255 (20%)20 (80%)9 (36%)16 (64%)9 (31%)20 (69%)2 (13%)14 (87%)
negative (%)4218 (43%)24 (57%)34 (81%)8 (19%)30 (68%)14 (32%)7 (58%)5 (42%)
Induction regimen
PALG5 (%)8 2 (25%)6 (75%)6 (75%)2 (25%)3 (38%)5 (63%)3 (60%)2 (40%)
PALG6 (%)8325 (30%)58 (70%)46 (55%)37 (45%)49 (54%)42 (46%)9 (26%)25 (74%)
PALG modified (%)20 (0%)2 (100%)0 (0%)2 (100%)0 (0%)2 (100%)0 (0%)2 (100%)
AlloHSCT in first CR
yes (%)40 14 (35%)26 (65%)29 (73%)11 (27%)25 (56%)20 (44%)8 (53%)7 (47%)
no (%)54 14 (26%)40 (74%)24 (44%)30 (56%)28 (49%)29 (51%)4 (15%)22 (85%)
Two hypodiploid cases are included in the “bad genetics” group but are not included in the preceding subgroups. computed by chi-squared or Fisher exact test. * computed by Mann-Whitney U test. # WBC in MLL::AF4 significantly higher than in NEG, Complex, and Hyperdiploid patients (ANOVA p = 0.0038, Tukey HSD test p = 0.0163, p = 0.0470 and p = 0.0271, respectively. Other comparisons not significant). a computed by chi-squared test: NEG vs. all other cases. b computed by chi-squared test: BCR::ABL1 vs. all other cases. c computed by chi-squared test: preB vs. all other cases. d “MLPA group”: this series includes IKZF1pos cases identified using MLPA (37) and RT-PCR (2) in a group of 94 patients, for whom other gene deletions were analyzed using MLPA; this series was used for evaluation of patients’ characteristic (see also Table S2). e “MLPA+RT-PCR group”: this series includes IKZF1pos cases identified using MPLA (37) and RT-PCR (10); This series was used for survival analyses. Abbreviations: RAG2, Recombination Activating Gene 2; AID, Activation Induced Cytidine Deaminase; WBC, White Blood Cell count; PALG, Polish Adult Leukemia Group; CNA, copy number alterations.
Table 2. Prevalence of markers of CNA-type genetic instability: mutational burden, gene type—in subgroups according to mutator enzyme expression.
Table 2. Prevalence of markers of CNA-type genetic instability: mutational burden, gene type—in subgroups according to mutator enzyme expression.
Total B-ALL
RAG2 ExpressionAID Expressionp 
Low
n = 43
High
n = 48
Low
n = 54
High
n = 36
CNA
CNAneg18 (67%)9 (33%)15 (54%)13 (46%)
CNApos25 (39%)39 (61%)36 (61%)23 (39%)RAG2: CNAneg vs. CNApos: 0.0160
1 CNA9 (39%)14 (61%)11 (50%)11 (50%)
CNAhigh16 (39%)25 (61%)25 (68%)12 (32%)RAG2: CNAneg vs. CNAhigh: 0.0257
IKZF1 deletion
IKZF1neg31 (60%)21 (40%)24 (46%)28 (54%)RAG2: 0.0109
IKZF1pos16 (34%)31 (66%)31 (72%)12 (28%)AID: 0.0108
CDKN2A/B deletion RAG2: 0.0013
CDKN2A/Bneg33 (56%)26 (44%)36 (63%)21 (37%)
CDKN2A/Bpos10 (31%)22 (69%)15 (50%)15 (50%)
PAX5 deletion NS
PAX5neg36 (47%)40 (53%)41 (55%)33 (45%)
PAX5pos7 (47%)8 (53%)10 (77%)3 (23%)
NEG B-ALL
RAG2 ExpressionAID Expressionp 
Low
n = 22
High
n = 12
Low
n = 12
High
n = 21
CNA RAG2: CNAneg vs. CNApos: 0.0464
CNAneg9 (90%)1 (10%)3 (27%)8 (73%)
CNApos13 (54%)11 (46%)9 (41%)13 (59%)
1 CNA4 (44%)5 (56%)2 (25%)6 (75%)
CNAhigh9 (60%)6 (40%)7 (50%)7 (50%)
IKZF1 deletion
IKZF1neg16 (80%)4 (20%)3 (15%)17 (85%)RAG2: 0.0257
IKZF1pos6 (43%)8 (57%)8 (67%)4 (33%)AID: 0.0029
BCR::ABL1pos B-ALL
RAG2 ExpressionAID Expressionp 
Low
n = 8
High
n = 21
Low
n = 22
High
n = 6
CNA NS
CNAneg3 (50%)3 (50%)4 (67%)2 (33%)
CNApos5 (22%)18 (78%)18 (82%)4 (18%)
1 CNA1 (14%)6 (86%)5 (71%)2 (29%)
CNAhigh4 (25%)12 (75%)13 (87%)2 (13%)
IKZF1 deletion NS
IKZF1neg4 (40%)6 (60%)6 (60%)4 (40%)
IKZF1pos7 (27%)19 (73%)20 (80%)5 (20%)
computed by chi-squared or Fisher exact test. Abbreviations: RAG2, Recombination Activating Gene 2; AID, Activation Induced Cytidine Deaminase; CNA, Copy Number Alterations.
Table 3. Patients’ outcomes according to the primary aberration. (A) Univariate analysis of patients’ outcomes according to the primary aberration. (B) Multivariate analysis of patients’ outcomes according to the primary aberrations and demographic data.
Table 3. Patients’ outcomes according to the primary aberration. (A) Univariate analysis of patients’ outcomes according to the primary aberration. (B) Multivariate analysis of patients’ outcomes according to the primary aberrations and demographic data.
(A)
End Point and VariablesTotal
n = 161
NEGBCR::ABL1
n = 54
“Bad Primary”E2A::PBX
n = 3
p
Total NEG *
n = 61
NK Only
n = 41
BCR::ABL1-like
n = 15
MLL::AF4
n = 5
Complex
n = 15
Hyperdiploid
n = 6
All “Bad Primary”
n = 43
CR NEG vs. BCR::ABL1: 0.5454
No. of patients131/15249/5734/4147/5410/125/511/144/532/383/3NEG vs. BCR::ABL1-like: 0.5555
(%)(86%)(86%)(83%)(87%)(83%)(100%)(79%)(80%)(78%)(100%)NEG vs. “Bad”: 0.5177
OS NEG vs. BCR::ABL1:0.0895 #
No. of patients161614154155156433NEG vs. BCR::ABL1-like: 0.1080 #
4-year rate ± SE34 ± 4%32 ± 8%25 ± 8%54 ± 8%10 ± 9% 53 ± 25%20 ± 10%0%16 ± 7%100%NEG vs. “Bad”: 0.0192 #
RFS NEG vs. BCR::ABL1: 0.0356 #
No. of patients129493446105106313NEG vs. BCR::ABL1-like: 0.3400 #
4-year rate ± SE42 ± 6%28 ± 9%19 ± 9%66 ± 8%18 ± 16%20 ± 18%24 ± 15%60 ± 22% 22 ± 10%100%NEG vs. “Bad”: 0.0365 #
DFS NEG vs. BCR::ABL1:0.0978 #
NEG vs. BCR::ABL1-like: 0.5055 #
NEG vs. “Bad”: 0.0319 #
No. of patients126463246105115313
4-year rate ± SE37 ± 5%25 ± 8%18 ± 9%58 ± 8%16 ± 15%20 ± 18%22 ± 13%40 ± 22%17 ± 8%100%
(B)
Primary Aberration
End Point and VariablesAgeWBCNEG
n = 61
BCR::ABL1
n = 54
All “Bad Primary”
n = 43
BCR::ABL1-like
n = 15
MLL::AF4
n = 5
Complex
n = 15
Hyperdiploid
n = 6
E2A::PBX
n = 3
CR
OR (95% CI)0.83 (0.70–0.98)1.08 (0.91–1.29)1.000.99 (0.71–1.39)0.90 (0.66–1.21)0.92 (0.64–1.31)1.07 (0.67–1.71)0.86 (0.61–1.21)0.86 (0.55–1.35)1.11 (0.65–1.92)
p *0.02940.3593(reference)0.97260.47820.62800.77380.37260.51560.6936
OS
HR (95% CI)1.04 (1.03–1.06)1.00 (1.00–1.00)1.000.55 (0.32–0.96)1.67 (1.03–2.73)1.39 (0.70–2.75)1.10 (0.24–4.99)1.62 (0.83–3.14)3.91 (1.49–10.31)2.62 (0.61–11.18)
p #<0.00010.8974(reference)0.03380.02610.34340.89840.15490.00580.1929
RFS -
HR (95% CI)1.05 (1.03–1.07)1.00 (1.00–1.01)1.000.34 (0.16–0.72)1.81 (0.98–3.37)1.21 (0.49–3.02)4.40 (1.37–14.13)1.72 (0.71–4.14)3.61 (0.81–16.18)
p #<0.00010.0929(reference)0.00490.05970.67820.01290.22590.0934
(A) * NEG includes patients negative for all molecular markers with NK, or other aberrations of unknown significance, or without cytogenetic data. Two hypodiploid cases are included in all the “bad genetics” groups but are not included in the preceding subgroups. computed by Fisher’s exact test. # computed by the log-rank test. Abbreviations: NK, normal karyotype; CR, complete remission; OS, overall survival; RFS, relapse-free survival; DFS, disease-free survival; SE, standard error. (B) * computed from General Regression Model. # computed by Cox proportional hazard regression. Abbreviations: WBC, white blood count; CR, complete remission; OS, overall survival; RFS, relapse-free survival; OR, odds ratio; HR, hazard ratio; 95% CI, 95% confidence interval.
Table 4. Patients’ outcomes according to the secondary aberrations. (A) Univariate analysis of patients’ outcomes according to the secondary aberrations. (B) Multivariate analysis of patients’ outcomes according to the secondary aberrations, adjusted for age and WBC.
Table 4. Patients’ outcomes according to the secondary aberrations. (A) Univariate analysis of patients’ outcomes according to the secondary aberrations. (B) Multivariate analysis of patients’ outcomes according to the secondary aberrations, adjusted for age and WBC.
(A)
End Point and VariablesTotalCNA Presence CNA Mutation Burden IKZF1 Mutation Status Bad CNA
(CNAhigh/IKZFpos)
Good CNA
(CNApos Other than Bad CNA)
CNAnegCNApospCNAlowCNAhighpIKZFnegIKZFpospCNAhigh
IKZFpos
All Other PatientspGood CNAAll Other Patientsp
TOTAL B-ALL
CR 0.5283 0.0422  0.3746 0.0514 0.0880
No. of patients73/8822/2651/6245/5028/3841/4838/4720/2853/6031/3442/54
(%)83%85%82%90%74%85%85%71%88%91%78%
OS 0.4921 # 0.0179 # 0.2734 # 0.0045 # 0.0427 #
No. of patients9227655240524728643755
4-year rate ± SE36 ± 6%43 ± 10%28 ± 6%42 ± 8%20 ± 7%41 ± 8%26 ± 7%9 ± 6%42 ± 7%43 ± 9%25 ± 6%
RFS 0.1460 # 0.0412 # 0.0320 # 0.0035 # 0.2096 #
No. of patients7322494427403719523041
4-year rate ± SE41 ± 7%43 ± 10%34 ± 8%53 ± 9%20 ± 10%48 ± 10%31 ± 9%14 ± 9%49 ± 8%45 ± 10%40 ± 9%
NEG B-ALL
CR 0.6317 0.1710 0.5419 0.3105 0.4581
No. of patients28/329/1019/2218/1910/1317/1911/137/921/2312/1316/19
(%)88%90%86%95%77%89%85%78%91%92%84%
OS 0.2117 # 0.0016 # 0.0055 # 0.0051 # 0.1513 #
No. of patients3410241915201410241420
4-year rate ± SE28 ± 9%53 ± 17%22 ± 9%45 ± 15%7 ± 6%46 ± 13%7 ± 7%0%40 ± 12%28 ± 13%25 ± 10%
RFS 0.3322 # 0.0171 # 0.0118 # 0.0037 # 0.1042 #
No. of patients28919181017117211216
4-year rate ± SE25 ± 10%42 ± 20%19 ± 11%43 ± 14%0%36 ± 15%10 ± 10%0%34 ± 14%31 ± 17%22 ± 12%
BCR::ABL1pos B-ALL
CR 0.3438 0.3954 0.5441 0.2613 0.0653
No. of patients24/304/620/2412/1412/168/1022/2610/1414/1610/1014/20
(%)80%67%83%86%75%80%85%71%88%100%70%
OS 0.3019 # 0.2129 # 0.6468 # 0.0406 # 0.0025 #
No. of patients296231415102513161019
4-year rate ± SE47 ± 10%33 ± 19%51 ± 11%62 ± 13%33 ± 13%60 ± 15%45 ± 11%20 ± 12%67 ± 12%89 ± 10%25 ± 11%
RFS 0.8073 # 0.1157 # 0.0875 # 0.0231 # 0.0404 #
No. of patients2341912118219141013
4-year rate ± SE65 ± 11%75 ± 22%63 ± 12%82 ± 12%44 ± 17%88 ± 12%48 ± 12%29 ± 17%85 ± 10%89 ± 10%44 ± 15%
(B)
End Point and VariablesCNAhighIKZFposBad CNA
(CNAhigh/IKZFpos)
Good CNA
(CNApos Other than Bad CNA)
TOTAL B-ALL
CR
OR (95% CI)0.83 (0.67–1.02)0.97 (0.78–1.21)0.83 (0.67–1.04)1.20 (0.97–1.48)
p *0.07910.81560.10540.0952
OS
HR (95% CI)1.81 (1.04–3.17)0.91 (0.51–1.61)1.92 (1.06–3.48)0.56 (0.31–1.01)
p 0.03740.73930.03180.0521
RFS
HR (95% CI)1.97 (0.99–3.92)1.48 (0.74–2.97)2.89 (1.37–6.10)0.57 (0.28–1.16)
p 0.05380.26950.00540.1189
NEG B-ALL
CR
OR (95% CI)0.68 (0.47–0.98)0.87 (0.58–1.31)0.74 (0.50–1.10)1.15 (0.77–1.71)
p *0.03970.49210.13050.4897
OS
HR (95% CI)4.95 (1.70–14.39)4.12 (1.37–12.33)4.85 (1.56–15.09)0.42 (0.14–1.20)
p 0.00330.01140.00650.1050
RFS
HR (95% CI)4.11 (1.30–12.99)4.59 (1.32–15.97)9.80 (2.32–41.35)0.32 (0.09–1.08)
p 0.01600.01650.00190.0660
BCR:ABL1pos B-ALL
CR
OR (95% CI)0.97 (0.68–1.40)1.16 (0.84–1.62)0.94 (0.65–1.38)1.28 (0.88–1.87)
p *0.88600.35950.74850.1871
OS
HR (95% CI)1.84 (0.60–5.63)0.88 (0.27–2.91)2.36 (0.76–7.27)0.11 (0.01–0.90)
p 0.28600.83260.13590.0399
RFS
HR (95% CI)6.29 (0.89–44.47)3.52 (0.40–31.35)7.65 (1.21–48.31)0.14 (0.02–1.22)
p 0.06540.25640.03040.0756
computed by chi-squared or Fisher’s exact test. # computed by log-rank test. * computed from General Linear Model. computed by Cox proportional hazard regression. Abbreviations: WBC, white blood count; CR, complete remission; OS, overall survival; RFS, relapse-free survival; SE, standard error; CNA, copy number alteration; HR, hazard ratio; OR, odds ratio; 95% CI, 95% confidence interval.
Table 5. Outcome of B-ALL patients according to the revised code for risk stratification.
Table 5. Outcome of B-ALL patients according to the revised code for risk stratification.
End Point and VariablesIntermediate RiskLow RiskHigh RiskTotal
CR
No. of patients24/2813/1344/5581/96
(%)86%100%80%85%
p *(reference)0.15140.5224
OS
No. of patients291358100
4-year rate ± SE41 ± 10%75 ± 13%13 ± 5%30 ± 5%
p (reference)0.01480.0011
RFS
No. of patients26134281
4-year rate ± SE39 ± 12%90 ± 9%23 ± 7%38 ± 7%
p (reference)0.02540.0037
DFS
No. of patients24134380
4-year rate ± SE36 ± 11%90 ± 9%17 ± 6%34 ± 6%
p (reference)0.01830.0010
Multivariate analysis adjusted for age and WBC
CR
OR (95% CI)1.000.92 (0.73–1.16)1.12 (0.89–1.40)
p (reference)0.32590.4713
OS
HR (95% CI)1.000.13 (0.02–0.99)2.12 (1.14–3.93)
p #(reference)0.04840.0176
RFS
HR (95% CI)1.000.16 (0.02–1.23)2.59 (1.25–5.33)
p #(reference)0.07870.0100
* computed from the chi-squared test. computed from log-rank test. computed from a General Linear Model. # computed by Cox proportional hazard regression. All p-values are calculated as a comparison to the intermediate risk group. Abbreviations: CR, complete remission; OS, overall survival; RFS, relapse-free survival; DFS, disease-free survival; SE, standard error; HR, hazard ratio; OR, odds ratio; 95% CI, 95% confidence interval. Intermediate risk: BCR::ABL1pos: CNAneg; NEG: CNAneg or 1 CNA/IKZF1pos or CNApos/IKZF1neg. Low risk: E2A::PBXpos; BCR::ABL1pos: 1 CNA/IKZF1pos or CNApos/IKZF1neg. High risk: MLL::AF4pos; Complex; Hyperdiploid; Hypodiploid; BCR::ABL1-like; CNAhigh/IKZF1pos.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Libura, M.; Karabin, K.; Tyrna, P.; Czyż, A.; Makuch-Łasica, H.; Jaźwiec, B.; Paluszewska, M.; Piątkowska-Jakubas, B.; Zawada, M.; Gniot, M.; et al. Prognostic Impact of Copy Number Alterations’ Profile and AID/RAG Signatures in Acute Lymphoblastic Leukemia (ALL) with BCR::ABL and without Recurrent Genetic Aberrations (NEG ALL) Treated with Intensive Chemotherapy. Cancers 2023, 15, 5431. https://doi.org/10.3390/cancers15225431

AMA Style

Libura M, Karabin K, Tyrna P, Czyż A, Makuch-Łasica H, Jaźwiec B, Paluszewska M, Piątkowska-Jakubas B, Zawada M, Gniot M, et al. Prognostic Impact of Copy Number Alterations’ Profile and AID/RAG Signatures in Acute Lymphoblastic Leukemia (ALL) with BCR::ABL and without Recurrent Genetic Aberrations (NEG ALL) Treated with Intensive Chemotherapy. Cancers. 2023; 15(22):5431. https://doi.org/10.3390/cancers15225431

Chicago/Turabian Style

Libura, Marta, Karolina Karabin, Paweł Tyrna, Anna Czyż, Hanna Makuch-Łasica, Bożena Jaźwiec, Monika Paluszewska, Beata Piątkowska-Jakubas, Magdalena Zawada, Michał Gniot, and et al. 2023. "Prognostic Impact of Copy Number Alterations’ Profile and AID/RAG Signatures in Acute Lymphoblastic Leukemia (ALL) with BCR::ABL and without Recurrent Genetic Aberrations (NEG ALL) Treated with Intensive Chemotherapy" Cancers 15, no. 22: 5431. https://doi.org/10.3390/cancers15225431

APA Style

Libura, M., Karabin, K., Tyrna, P., Czyż, A., Makuch-Łasica, H., Jaźwiec, B., Paluszewska, M., Piątkowska-Jakubas, B., Zawada, M., Gniot, M., Trubicka, J., Szymańska, M., Borg, K., Więsik, M., Czekalska, S., Florek, I., Król, M., Paszkowska-Kowalewska, M., Gil, L., ... Jędrzejczak, W. W. (2023). Prognostic Impact of Copy Number Alterations’ Profile and AID/RAG Signatures in Acute Lymphoblastic Leukemia (ALL) with BCR::ABL and without Recurrent Genetic Aberrations (NEG ALL) Treated with Intensive Chemotherapy. Cancers, 15(22), 5431. https://doi.org/10.3390/cancers15225431

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop