In the search for genes that define critical steps of relapse in pediatric T-cell acute lymphoblastic leukemia (T-ALL) and can serve as prognostic markers, we performed targeted sequencing of 313 leukemia-related genes in 214 patients: 67 samples collected at the time of relapse and 147 at initial diagnosis. As relapse-specific genetic events, we identified activating mutations in NT5C2 (P=0.0001, Fisher’s exact test), inactivation of TP53 (P=0.0007, Fisher’s exact test) and duplication of chr17:q11.2-24.3 (P=0.0068, Fisher’s exact test) in 32/67 of T-ALL relapse samples. Alterations of TP53 were frequently homozygous events, which significantly correlated with higher rates of copy number alterations in other genes compared with wild-type TP53 (P=0.0004, Mann–Whitney’s test). We subsequently focused on mutations with prognostic impact and identified genes governing DNA integrity (TP53, n=8; USP7, n=4; MSH6, n=4), having key roles in the RAS signaling pathway (KRAS, NRAS, n=8), as well as IL7R (n=4) and CNOT3 (n=4) to be exclusively mutated in fatal relapses. These markers recognize 24/49 patients with a second event. In 17 of these patients with mostly refractory relapse and dire need for efficient treatment, we identified candidate targets for personalized therapy with p53 reactivating compounds, MEK inhibitors or JAK/STAT-inhibitors that may be incorporated in future treatment strategies.

The treatment of acute lymphoblastic leukemia (ALL) and osteosarcoma (OSC) is very effective: the vast majority of patients recover and survive for decades. However, they still need to face serious adverse effects of chemotherapy. One of these is cardiotoxicity which may lead to progressive heart failure in the long term. Cardiotoxicity is contributed mainly to the use of anthracyclines and might have genetic risk factors. Our goal was to test the association between left ventricular function and genetic variations of candidate genes. Echocardiography data from medical records of 622 pediatric ALL and 39 OSC patients were collected from the period 1989–2015. Fractional shortening (FS) and ejection fraction (EF) were determined, 70 single nucleotide polymorphisms (SNPs) in 26 genes were genotyped. Multivariate logistic regression and multi-adjusted general linear model were performed to investigate the influence of genetic polymorphisms on the left ventricular parameters. Bayesian network based Bayesian multilevel analysis of relevance (BN-BMLA) method was applied to test for the potential interaction of the studied cofactors and SNPs. Our results indicate that variations in ABCC2, CYP3A5, NQO1, SLC22A6 and SLC28A3 genes might influence the left ventricular parameters. CYP3A5 rs4646450 TT was 17% among ALL cases with FS lower than 28, and 3% in ALL patients without pathological FS (p = 5.60E-03; OR = 6.94 (1.76–27.39)). SLC28A3 rs7853758 AA was 12% in ALL cases population, while only 1% among controls (p = 6.50E-03; OR = 11.56 (1.98–67.45)). Patients with ABCC2 rs3740066 GG genotype had lower FS during the acute phase of therapy and 5–10 years after treatment (p = 7.38E-03, p = 7.11E-04, respectively). NQO1 rs1043470 rare T allele was associated with lower left ventricular function in the acute phase and 5–10 years after the diagnosis (p = 4.28E-03 and 5.82E-03, respectively), and SLC22A6 gene rs6591722 AA genotype was associated with lower mean FS (p = 1.71E-03), 5–10 years after the diagnosis. Genetic variants in transporters and metabolic enzymes might modulate the individual risk to cardiac toxicity after chemotherapy.

Hypersensitivity reactions are the most frequent dose-limiting adverse reactions to Escherichia coli-derived asparaginase in pediatric acute lymphoblastic leukemia patients. The aim of the present study was to identify associations between sequence-based Human Leukocyte Antigen Class II region alleles and asparaginase hypersensitivity in a Hungarian acute lymphoblastic leukemia population. Four-digit typing of HLA-DRB1 and HLA-DQB1 loci was performed in 359 pediatric acute lymphoblastic leukemia patients by using next-generation sequencing method. Based on genotypic data of the two loci, haplotype reconstruction was carried out. In order to investigate the possible role of the HLA-DQ complex, the HLA-DQA1 alleles were also inferred. Multivariate logistic regression analysis and a Bayesian network-based approach were applied to identify relevant genetic risk factors of asparaginase hypersensitivity. Patients with HLA-DRB1*07:01 and HLA-DQB1*02:02 alleles had significantly higher risk to develop asparaginase hypersensitivity compared to noncarriers (p=4.56x10^-5; OR=2.86 (1.73-4.75) and p=1.85x10^-4; OR=2.99 (1.68-5.31); n=359, respectively). After haplotype reconstruction the HLA-DRB1*07:01-HLA-DQB1*02:02 haplotype was associated with an increased risk. After inferring the HLA-DQA1 alleles the HLA-DRB1*07:01–HLA-DQA1*02:01–HLA-DQB1*02:02 haplotype was associated with the highest risk of asparaginase hypersensitivity (p=1.22x10^-5; OR=5.00 (2.43-10.29); n=257). Significantly fewer proportion of T-cell acute lymphoblastic leukemia patients carried the HLA-DQB1*02:02 allele and the associated haplotype than did pre-B-cell acute lymphoblastic leukemia patients (6.5%; vs. 19.2%, respectively; p=0.047). In conclusion, we identified a haplotype in the Human Leukocyte Antigen Class II region associated with a higher risk of asparaginase hypersensitivity. Our results confirm that variations in HLA-D region might influence the development of asparaginase hypersensitivity.

Background: Refractory central nervous system (CNS) involvement is among the major causes of therapy failure in childhood acute leukemia. Applying contemporary diagnostic methods, CNS disease is often underdiagnosed. To explore more sensitive and less invasive CNS status indicators, we examined microRNA (miR) expressions and extracellular vesicle (EV) characteristics. Methods: In an acute lymphoblastic leukemia (ALL) discovery cohort, 47 miRs were screened using Custom TaqMan Advanced Low-Density Array gene expression cards. As a validation step, a candidate miR family was further scrutinized with TaqMan Advanced miRNA Assays on serial cerebrospinal fluid (CSF), bone marrow (BM) and peripheral blood samples with different acute leukemia subtypes. Furthermore, small EV-rich fractions were isolated from CSF and the samples were processed for immunoelectron microscopy with anti-CD63 and anti-CD81 antibodies, simultaneously. Results: Regarding the discovery study, principal component analysis identified the role of miR-181-family (miR-181a-5p, miR-181b-5p, miR-181c-5p) in clustering CNS-positive (CNS⁺) and CNS-negative (CNS^‒) CSF samples. We were able to validate miR-181a expression differences: it was about 52 times higher in CSF samples of CNS⁺ ALL patients compared to CNS^‒ cases (n = 8 vs. n = 10, ΔFC = 52.30, p = 1.5E−4), and CNS⁺ precursor B cell subgroup also had ninefold higher miR-181a levels in their BM (p = 0.04). The sensitivity of CSF miR-181a measurement in ALL highly exceeded those of conventional cytospin in the initial diagnosis of CNS leukemia (90% vs. 54.5%). Pellet resulting from ultracentrifugation of CNS⁺ CSF samples of ALL patients showed atypical CD63⁻/CD81⁻ small EVs in high density by immunoelectron microscopy. Conclusions: After validating in extensive cohorts, quantification of miR-181a or a specific EV subtype might provide novel tools to monitor CNS disease course and further adjust CNS-directed therapy in pediatric ALL.

Background: Treatment stratification based on bone marrow minimal residual disease (MRD) at set time points has resulted in considerably improved survival in pediatric acute lymphoblastic leukemia (ALL). Treatment response is assessed using bone marrow samples. MicroRNAs (miRs) easily traffic among fluid spaces and are more stable than most other RNA classes. We examined the role of circulating miRs as putative less invasive MRD biomarkers. Methods: In an exploratory experiment, expression of 46 preselected miRs was studied in platelet-free blood plasma samples of 15 de novo, 5 relapsed ALL patients and 10 controls by Custom TaqMan Array Advanced MicroRNA Card. Based on their high expression in ALL compared to controls, and on the reduction observed along the induction therapy, four miRs were selected for further analyses: miR-128-3p, -181a-5p, -181b-5p and 222-3p. Their expression was measured by qPCR at 4 time points in 27 de novo ALL patients treated in the ALL IC-BFM 2009 study. Results: The expression of all 4 miRs significantly decreased over the first week of therapy (miR-128-3p: log₂ fold change − 2.86; adjusted p 3.6 × 10⁻⁷; miR-181b-5p: log₂ fold change − 1.75; adjusted p 1.48 × 10⁻²; miR-181a-5p: log₂ fold change -1.33; adjusted p 3.12 × 10⁻²; miR-222-3p: log₂ fold change − 1.25; adjusted p 1.66 × 10⁻²). However, no significant further reduction in miR expression was found after the 8th day of therapy. Measured drop in expression of 2 miRs at day 8 strongly correlated with day 15 bone marrow flow cytometry MRD results (miR-128-3p: Pearson’s r = 0.88, adjusted p = 2.71 × 10⁻⁴; miR-222-3p: r = 0.81, adjusted p = 2.99 × 10⁻³). Conclusion: In conclusion, these circulating miRs might act as biomarkers of residual leukemia. MiR-128-3p and miR-222-3p in blood predict day 15 flow cytometry MRD results 7 days earlier. Although, their sensitivity falls behind that of bone marrow flow cytometry MRD at day 15.

L-asparaginase (ASP) is a key element in the treatment of paediatric acute lymphoblastic leukaemia (ALL). However, hypersensitivity reactions (HSRs) to ASP are major challenges in paediatric patients. Our aim was to investigate genetic variants that may influence the risk to Escherichia coli-derived ASP hypersensitivity. Sample and clinical data collection was carried out from 576 paediatric ALL patients who were treated according to protocols from the Berlin—Frankfurt—Münster Study Group. A total of 20 single nucleotide polymorphisms (SNPs) in GRIA1 and GALNT10 genes were genotyped. Patients with GRIA1 rs4958351 AA/AG genotype showed significantly reduced risk to ASP hypersensitivity compared to patients with GG genotype in the T-cell ALL subgroup (OR = 0.05 (0.01–0.26); p = 4.70E-04), while no such association was found in pre-B-cell ALL. In the medium risk group two SNPs of GRIA1 (rs2055083 and rs707176) were associated significantly with the occurrence of ASP hypersensitivity (OR = 0.21 (0.09–0.53); p = 8.48E-04 and OR = 3.02 (1.36–6.73); p = 6.76E-03, respectively). Evaluating the genders separately, however, the association of rs707176 with ASP HSRs was confined only to females. Our results suggest that genetic variants of GRIA1 might influence the risk to ASP hypersensitivity, but subgroups of patients can differ significantly in this respect.