We studied a cohort of 99 patients diagnosed with B-ALL from January 2012 to December 2020 and treated according to the ALL IC-BFM 2009 protocol. Eighty-eight bone marrow or peripheral blood samples were analysed for copy number variations (CNVs) using the SALSA MLPA P335 ALL-IKZF1 probemix.
At least one CNV was detected in more than 65% of analysed samples. The most frequently altered genes were
Our results are in concordance with the results obtained in larger cooperative clinical trials. Copy number variations analysis using the SALSA MLPA kit is a reliable tool for initial diagnostic approach in children with B-ALL, even in smaller institutions in low- and middle-income countries.
Key words
- B-acute lymphoblastic leukemia
- deletions
- MLPA
- pediatric
- copy number variations (CNVs)
Improvements in risk stratification and new therapeutic approaches have dramatically improved treatment outcomes in pediatric B-cell acute lymphoblastic leukemia (B-ALL). In developed countries, the overall survival for these patients is approaching 90%.1,2 Nevertheless, some genetic subtypes still imply poor outcomes, and 10–20% of patients experience a relapse that is often accompanied by treatment resistance and failure.3,4 Therefore, the need for new diagnostic and prognostic markers remains of paramount importance.
In the previous years, deletions in the
Recently, another, minimal residual disease (MRD) dependent prognostic profile
Deletions in
In total, 99 children with B-ALL that were treated at the University Children’s Hospital, University Medical Centre Ljubljana between January 2012 and December 2020 according to the ALL IC-BFM 2009 protocol were included in this study. Diagnoses were established following standard clinical, cytomorphological, and immunological criteria.
We obtained bone marrow samples for 92 patients as part of the diagnostic procedure before starting treatment. For 7 patients, bone marrow samples were not available, therefore, peripheral blood samples were obtained for the analysis. For four patients, there was no sufficient material available to perform the multiplex ligation-dependent probe amplification (MLPA) assay, 3 samples contained less than 40% of blasts and were excluded from analysis, and for an additional four, the assay failed due to poor sample quality. Therefore, data analysis was performed on 88 patient samples (82 bone marrow and 6 peripheral blood). The bone marrow samples contained 77.4 ± 16.7% of blast in average, and for the peripheral blood samples this value was 77.8 ± 16.0%. For survival analysis, patient samples from the year 2020 were excluded, due to the short follow-up period. Therefore, the survival analysis was carried out on 72 patients diagnosed between January 2012 and December 2019.
Informed consent was obtained from all subjects involved in the study, or their parents. The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the Republic of Slovenia (reference number KME 51/03/11).
Genomic DNA was extracted from bone marrow or peripheral blood samples using the FlexiGene DNA Kit 250 (QIAGEN®, Hilden, Germany) according to the manufacturer’s instructions. All samples were quantified using DS-11 FX+ Spectrophotometer (DeNovix, Wilmington, USA), and stored at 4°C. Before analysis, the DNA concentration was established at 25 ± 1 ng/μL for the MLPA assay.
DNA was analysed for copy number alterations using the SALSA MLPA P335 ALL-IKZF1 probemix, according to the manufacturer’s instructions (MRC Holland, Amsterdam, the Netherlands). The P335 probemix allows for the detection of deletions and duplications in B-cell differentiation and cell cycle control genes (
MLPA reactions were carried out on a 96-well PCR thermocycler SimpliAmp Thermal Cycler (Applied Biosystems, Thermo Fisher, Massachusetts, USA), and the products were separated by capillary electrophoresis on an ABI-3500 genetic analyser (Applied Biosystems, Thermo Fisher, Massachusetts, USA). The resulting peak intensities were analysed using Coffalyser software (MRC-Holland) which performed the intrasample and intersample normalization of the peaks with the manufacturer’s reference probes and normal control DNA, respectively. Values above 1.3 were considered as gain, between 1.3 and 0.75 normal, between 0.75 and 0.25 heterozygous loss, and below 0.25 homozygous loss.
All statistical analyses were performed using SPSS 22.0 software (IBM Corp., Armonk, NY, USA). Event-free survival (EFS; defined as the time between diagnosis and relapse or death) and overall survival (OS; defined as time between diagnosis and death or last follow-up) were analysed using the Kaplan-Meier method and the differences between multiple groups were analysed using the log-rank test. Multivariate analysis was performed using a Cox regression model, which was adjusted for other risk factors, namely sex, age at diagnosis, and risk group (HR
The cohort included 54 males and 45 females (N = 99), with median age 4 years (range from 1 day to 23 years). Of these, 75 harboured recurrent genetic abnormalities (28
Figure 1
Prevalence of ALL subtypes in the Slovenian pediatric B-ALL cohort.

The demographic and clinical characteristics of Slovenian B-ALL patients included in the study
Characteristic | |
---|---|
99 | |
54 (54.5%) | |
45 (45.5%) | |
28 (28.3%) | |
7 (7.1%) | |
4 (4.0%) | |
4 (4.0%) | |
Hyperdiploidy | 27 (27.3%) |
Hypodiploidy | 3 (3.0%) |
iAMP21 | 2 (2.0%) |
No abnormalities recurrent |
24 (24.2%) |
< 1 | 3 (3.0%) |
1–5 | 57 (57.6%) |
≥ 6 |
39 (39.4%) |
Standard risk | 17 (17.2%) |
Intermediate risk | 59 (59.6%) |
High risk | 23 (23.2%) |
< 0.1% | 35 (35.4%) |
0.1–10% | 48 (48.5%) |
> 10% | 13 (13.1%) |
Unknown | 3 (3.0%) |
< 0.01% | 73 (73.7%) |
0.01–1% | 20 (20.2%) |
> 1% | 3 (3.0%) |
Unknown | 3 (3.0%) |
Altogether, 92 patient samples and 5 controls were analysed by MLPA using the SALSA MLPA P335 ALL-IKZF1 probemix. Four patient samples failed the MLPA analysis. Out of the remaining 88 samples,
The comparison of patient characteristics depending on
Figure 2
Primary genetic alterations in patients with

Patients’ characteristics and response to treatment according to
Characteristic | |||
---|---|---|---|
72 | 5 | 11 | |
Male | 34 (47.2%) | 5 (100%) | 9 (81.8%) |
Female | 38 (52.8%) | 0 (0.0%) | 2 (18.2%) |
ETV6-RUNX1 | 24 (33.3%) | 0 (0.0%) | 1 (9.1%) |
BCR-ABL1 | 1 (1.4%) | 2 (40.0%) | 4 (36.4%) |
KMT2A rearrangements | 4 (5.6%) | 0 (0.0%) | 0 (0.0%) |
TCF3-PBX1 | 3 (4.2%) | 0 (0.0%) | 1 (9.1%) |
Hyperdiploidy | 20 (27.8%) | 2 (40.0%) | 1 (9.1%) |
Hypodiploidy | 2 (2.8%) | 0 (0.0%) | 1 (9.1%) |
iAMP21 | 1 (1.4%) | 0 (0.0%) | 0 (0.0%) |
No recurrent abnormalities | 18 (25.0%) | 1 (20.0%) | 3 (27.3%) |
< 1 | 3 (4.2%) | 0 (0.0%) | 0 (0.0%) |
1–5 | 41 (56.9%) | 2 (40.0%) | 5 (45.5%) |
≥ 6 | 28 (38.9%) | 3 (60.0%) | 6 (54.5%) |
Standard risk | 13 (18.1%) | 0 (0.0%) | 0 (0.0%) |
Intermediate risk | 46 (63.9%) | 2 (40.0%) | 4 (36.4%) |
High risk | 13 (18.1%) | 3 (60.0%) | 7 (63.6%) |
Day 15 | |||
< 0.1% | 30 (41.7%) | 0 (0.0%) | 1 (9.1%) |
0.1–10% | 33 (45.8%) | 2 (40.0%) | 6 (54.5%) |
> 10% | 6 (8.3%) | 3 (60.0%) | 4 (36.4%) |
Unknown | 3 (4.2%) | 0 (0.0%) | 0 (0.0%) |
< 0.01% | 56 (77.8%) | 0 (0.0%) | 8 (72.7%) |
0.01–1% | 13 (18.1%) | 3 (60.0%) | 2 (18.2%) |
> 1% | 2 (2.8%) | 1 (20.0%) | 0 (0.0%) |
Unknown | 1 (1.4%) | 1 (20.0%) | 1 (9.1%) |
The
In our cohort, 13 patients (14.8%) experienced an event (either relapse or death). Among these, four patients had no recurrent genetic alterations, three patients carried the
The SALSA MLPA P335 ALL-IKZF1 probemix can detect deletions or duplications in the following B-cell differentiation and cell cycle control genes:
Figure 3
The number of CNVs present in Slovenian B-ALL samples.

The most common CNVs were those in the
Figure 4
Frequency of copy number variations:

First, the patients were divided into two groups, a group with and a group without
The 5-year EFS and OS were also compared between groups with no
A multivariate Cox regression model was applied to this data to see, whether after adjusting for other relevant risk factors,
Figure 5

We further looked at the group of patients, classified as non-high risk (either SR or IR). In this group, six patients carried the
Figure 6

Various alterations in genes involved in cell differentiation and cell cycle regulation are a hallmark of B-ALL. The role of deletions in the
Our study was carried out on a smaller (88 patients), yet consecutive, unselected, and well-controlled population of pediatric patients with B-ALL. This is the first report about CNVs in cell differentiation and cell cycle regulation genes in Slovenian pediatric B-ALL patients. In our cohort,
The most common
In our cohort,
Previously published studies have shown a poorer event-free and overall survival of patients with
The CNV analysis of other genes showed that the most common alterations were in the
In our cohort, amplifications in the PAR1 region were observed quite frequently. Altogether, 17 patients (19.3%) had at least one gene amplification in this region. This is due to the fact that our cohort is unselected, and therefore also includes patients with hyperdiploidy. In hyperdiploidy, gains in the X chromosome are very common (present in 70% of hyperdiploid childhood B-ALL cases)38, and indeed, 14 out of our 17 patients with amplifications in PAR1 had a hyperdiploid karyotype. This karyotype is associated with a favourable outcome. However, the hyperdiploid patients in our cohort did not have a significantly better 5-year EFS, and the same was seen for the patients with amplifications in the PAR1 region that were identified with MLPA. Two patients carried deletions in this region (namely in
Despite the limitations of our study due to a lower number of analyzed samples and relatively short follow-up period, it produced results that are in concordance with the results obtained in larger cooperative clinical trials. We have shown that it is possible to provide comparable results regarding the presence of certain CNVs and their prognostic value in pediatric B-ALL patients even within a single-center experience. This study is only a starting point for the more comprehensive screening of patients diagnosed with B-ALL in Slovenia that we have planned for the future and will enable us to better evaluate and treat these patients.
Figure 1

Figure 2

Figure 3

Figure 4

Figure 5
![(A) Event-free survival in patients with or without IKZF1 deletions (5-year event-free survival [EFS] 54.8% vs. 85.9%, p = 0.016). (B) Overall survival in patients with or without IKZF1 deletions (5-year overall survival [OS] 81.5% vs. 93.0%, p = 0.295).](https://sciendo-parsed-data-feed.s3.eu-central-1.amazonaws.com/6062cc35282c524fbc6e2e22/j_raon-2021-0050_fig_005.jpg?X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Date=20220520T170258Z&X-Amz-SignedHeaders=host&X-Amz-Expires=18000&X-Amz-Credential=AKIA6AP2G7AKDOZOEZ7H%2F20220520%2Feu-central-1%2Fs3%2Faws4_request&X-Amz-Signature=61e6b4a5f03381ea18fc1afc2ada7d5d5749ff0f61fcd68977d9bd9679e8b3a3)
Figure 6

The demographic and clinical characteristics of Slovenian B-ALL patients included in the study
Characteristic | |
---|---|
99 | |
54 (54.5%) | |
45 (45.5%) | |
28 (28.3%) | |
7 (7.1%) | |
4 (4.0%) | |
4 (4.0%) | |
Hyperdiploidy | 27 (27.3%) |
Hypodiploidy | 3 (3.0%) |
iAMP21 | 2 (2.0%) |
No abnormalities recurrent |
24 (24.2%) |
< 1 | 3 (3.0%) |
1–5 | 57 (57.6%) |
≥ 6 |
39 (39.4%) |
Standard risk | 17 (17.2%) |
Intermediate risk | 59 (59.6%) |
High risk | 23 (23.2%) |
< 0.1% | 35 (35.4%) |
0.1–10% | 48 (48.5%) |
> 10% | 13 (13.1%) |
Unknown | 3 (3.0%) |
< 0.01% | 73 (73.7%) |
0.01–1% | 20 (20.2%) |
> 1% | 3 (3.0%) |
Unknown | 3 (3.0%) |
Patients’ characteristics and response to treatment according to IKZF1 deletion status in 91 Slovenian pediatric B-ALL patients
Characteristic | |||
---|---|---|---|
72 | 5 | 11 | |
Male | 34 (47.2%) | 5 (100%) | 9 (81.8%) |
Female | 38 (52.8%) | 0 (0.0%) | 2 (18.2%) |
ETV6-RUNX1 | 24 (33.3%) | 0 (0.0%) | 1 (9.1%) |
BCR-ABL1 | 1 (1.4%) | 2 (40.0%) | 4 (36.4%) |
KMT2A rearrangements | 4 (5.6%) | 0 (0.0%) | 0 (0.0%) |
TCF3-PBX1 | 3 (4.2%) | 0 (0.0%) | 1 (9.1%) |
Hyperdiploidy | 20 (27.8%) | 2 (40.0%) | 1 (9.1%) |
Hypodiploidy | 2 (2.8%) | 0 (0.0%) | 1 (9.1%) |
iAMP21 | 1 (1.4%) | 0 (0.0%) | 0 (0.0%) |
No recurrent abnormalities | 18 (25.0%) | 1 (20.0%) | 3 (27.3%) |
< 1 | 3 (4.2%) | 0 (0.0%) | 0 (0.0%) |
1–5 | 41 (56.9%) | 2 (40.0%) | 5 (45.5%) |
≥ 6 | 28 (38.9%) | 3 (60.0%) | 6 (54.5%) |
Standard risk | 13 (18.1%) | 0 (0.0%) | 0 (0.0%) |
Intermediate risk | 46 (63.9%) | 2 (40.0%) | 4 (36.4%) |
High risk | 13 (18.1%) | 3 (60.0%) | 7 (63.6%) |
Day 15 | |||
< 0.1% | 30 (41.7%) | 0 (0.0%) | 1 (9.1%) |
0.1–10% | 33 (45.8%) | 2 (40.0%) | 6 (54.5%) |
> 10% | 6 (8.3%) | 3 (60.0%) | 4 (36.4%) |
Unknown | 3 (4.2%) | 0 (0.0%) | 0 (0.0%) |
< 0.01% | 56 (77.8%) | 0 (0.0%) | 8 (72.7%) |
0.01–1% | 13 (18.1%) | 3 (60.0%) | 2 (18.2%) |
> 1% | 2 (2.8%) | 1 (20.0%) | 0 (0.0%) |
Unknown | 1 (1.4%) | 1 (20.0%) | 1 (9.1%) |
Advancements in the radiooncological treatment of high-risk prostate cancer: a quarter century of achievements Safety and efficacy of drug-eluting microspheres chemoembolization under cone beam computed tomography control in patients with early and intermediate stage hepatocellular carcinoma Treatment of skin tumors with intratumoral interleukin 12 gene electrotransfer in the head and neck region: a first-in-human clinical trial protocol