Printer Friendly

PAX5 Alterations in a Consecutive Childhood B-Cell Acute Lymphoblastic Leukemia Cohort Treated Using the ALL IC-BFM 2009 Protocol.

Author(s): Klementina Crepinšek [1,3]; Nika Klobucar [1]; Tine Tesovnik [1,3]; Robert Šket [1,3]; Barbara Jenko Bizjan [1,3]; Jernej Kovac [1,3]; Marko Kavcic [2,3]; Tomaž Prelog [2]; Lidija Kitanovski [2,3]; Janez Jazbec [2,3]; Maruša Debeljak (corresponding author) [1,3,*]

1. Introduction

Childhood B-cell acute lymphoblastic leukemia (B-ALL) remains a challenging hematological malignancy, necessitating ongoing research to identify novel markers for improved treatment outcomes. Among the genes garnering significant attention in recent years, PAX5 stands out as a critical regulator of B-cell development. Understanding its genetic background and underlying mechanisms has become paramount to unravelling the pathogenesis of B-ALL and exploring its clinical implications.

The PAX5 gene belongs to the PAX gene family, encompassing transcription factors that control several developmental processes. PAX5 exhibits a predominant expression in B-cell precursors and plays a central role in their commitment, differentiation, and maintenance [1]. The protein PAX5 contains 391 amino acids and is divided into the N-terminus domain, responsible for binding to specific DNA sequences, the octapeptide and partial homeodomain in the center, involved in protein–protein interactions, and the transactivation and inhibitory domains at the C-terminus, which regulate gene expression. The functional integrity of the C-terminal domain stands out as a pivotal factor in B-cell development. The disruption of its normal function has been found to significantly promote the pathogenesis of leukemias, particularly B-cell ALL [2,3].

Numerous abnormalities associated with PAX5 have been identified, including gene rearrangements, sequence variants, focal intragenic amplifications, and deletions [4,5]. Comprehensive genomic studies have highlighted PAX5 as the most commonly altered gene in B-ALL, with the alterations being present in approximately one-third of pediatric and adult cases [5,6,7]. Recent evidence suggests that alterations in PAX5 are not merely secondary events in B-ALL, but often serve as the primary drivers of leukemogenesis [7]. Moreover, PAX5 variants are also sometimes present in the germline, contributing to B-ALL susceptibility [8,9]. This highlights the significance of examining both somatic and germline alterations in PAX5, expanding our understanding of disease predisposition and the interplay between genetic and inherited factors.

In light of the diverse range of genetic alterations involving PAX5, two distinct B-ALL genetic subtypes have been defined: PAX5 P80R and PAX5alt. The PAX5 P80R subtype has demonstrated variable prognostic effects across different cohort studies, indicating the need for further investigation to accurately elucidate its clinical implications. Conversely, patients with the PAX5alt subtype exhibit consistently poorer outcomes, especially in the presence of co-occurring IKZF1 deletions [7,10].

With the identification and definition of the two novel genetic subtypes involving PAX5 and identifying PAX5 as one of the key players also in Ph-like ALL [7,11,12], there is a growing interest in investigating the various changes and underlying mechanisms, as well as their implications for patient prognoses. Therefore, in this study, we aimed to identify patients within our B-ALL cohort with altered PAX5 using a comprehensive analysis approach, characterize the types of genetic changes, determine their somatic or germline origin, and analyze the clinical outcomes associated with these alterations, ultimately aiming to improve patient outcomes and provide a better understanding of the role of PAX5 in leukemogenesis.

2. Materials and Methods

2.1. Patients and Materials

A consecutive cohort of 99 patients with B-ALL that were diagnosed between January 2012 and December 2020 and treated at the Children’s Hospital of the University Medical Centre Ljubljana according to the ALL IC-BFM 2009 protocol was included in our study. Diagnoses were established following standard clinical, cytomorphological, cytogenetic, molecular genetic, and immunological criteria. Minimal residual disease (MRD) was measured by flow cytometry according to the established laboratory protocol.

We collected bone marrow for 92 patients and peripheral blood samples for 7 patients as part of the diagnostic procedure before starting treatment. For 7 of these 99 patients, we also collected samples to test for the presence of germline variants, buccal swab samples from 5 patients, and peripheral blood samples during remission for 2 patients.

For 4 patients, there was not sufficient material available to perform the multiplex ligation-dependent probe amplification (MLPA) assay, and for an additional 3, the assay failed due to poor sample quality. Therefore, MLPA data analysis included 92 patients from the cohort. For RNA sequencing, there was not sufficient material for two patients, therefore the analysis was carried out for 97 samples.

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Republic of Slovenia National Medical Ethics Committee on 10 March 2011 and 15 March 2022 (reference number KME 51/03/11 and KME 0120-17/2022/6). Informed consent was obtained from all subjects involved in the study. Written informed consent was obtained from the patients or parents to publish this paper or any other anonymous data and results.

2.2. DNA and RNA Isolation

The DNA from bone marrow and peripheral blood samples was isolated using the FlexGene DNA kit according to the manufacturer’s instructions (Qiagen, Hilden, Germany). The DNA from buccal swab samples was isolated using the QIAamp DNA Mini Kit according to the manufacturer instructions (Qiagen, Hilden, Germany). The extracted DNA was stored at 4 °C.

The bone marrow and peripheral blood samples for RNA isolation were stored in the TRIzol[sup.®] reagent (Thermno Fisher Scientific, Waltham, MA, USA). After the addition of chloroform and phase separation, the Direct-zol™ RNA MiniPrep Kit was used for RNA extraction (Zymo, Irvine, CA, USA). The extracted RNA was stored at –70 °C.

2.3. RNA Sequencing and Data Analysis

Bone marrow and peripheral blood samples at the time of diagnosis for transcriptome analysis were prepared using the NEBNext Globin & rRNA Depletion Kit Human/Mouse/Rat (New England Biolabs, Ipswich, MA, USA) and NEBNext Ultra II Directional RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA). We sequenced the samples with 100 million reads per sample using a NovaSeq 6000 system (Illumina, San Diego, CA, USA).

The RNA expression was determined by bioinformatic analysis of the data based on the SnakePipes tool [13]. Shortly, the Cutadapt tool [14] was used to remove adapter sequences and low-quality sequences. We used the STAR tool [15] to align the reads to the human reference genome (GRCh38) and featureCounts [16] to count gene expression levels.

To determine the B-ALL genetic subtype for each patient, we used two R-based tools, the gene expression classifier tool ALLCatchR [17] and MD-ALL [18]. ALLCatchR uses only a gene expression read count as the input and distinguishes between 21 BCP-ALL molecular subtypes (BCL2/MYC, CDX2/UBTF, CEBP, DUX4, ETV6::RUNX1, ETV6::RUNX1-like, HLF, hyperdiploid, iAMP21, IKZF1 N159Y, KMT2A, low hypodiploid, MEF2D, near haploid, NUTM1, PAX5 P80R, PAX5alt, Ph-like, Ph-pos, TCF3::PBX1, ZNF384). MD-ALL takes the read counts, VCF files, and raw outputs of fusion callers (FusionCatcher and/or Cicero) as input files and classifies B-ALL into 26 subtypes (BCL2/MYC, CDX2/UBTF, CRLF2 (non-Ph-like), DUX4, ETV6::RUNX1, ETV6::RUNX1-like, HLF, hyperdiploid, iAMP21, IKZF1 N159Y, KMT2A, KMT2A-like, low hypodiploid, low hyperdiploid, MEF2D, near haploid, NUTM1, PAX5 P80R, PAX5::ETV6, PAX5alt, Ph, Ph-like, TCF3::PBX1, ZEB2/CEBP, ZNF384, ZNF384-like).

Gene fusions were detected from RNA sequencing data using nf-core’s pipeline rnafusion on default settings (version 2.3.4) [19]. This pipeline uses several tools for fusion detection (Arriba [20], Pizzly [21], STAR-fusion [22], Squid [23], and FusionCatcher [24]) and combines the results from all tools in a joint report for each sample. Additionally, we used CICERO [25] for fusion detection, as the use of this tool is needed for a complete subtype analysis using the MD-ALL tool. When needed, the data were manually inspected to confirm the presence of the fusion transcripts.

All sequences were aligned to the human reference genome (GRCh38) and variants were detected using the nf-core rnavar pipeline (version 1.0.0) [26] for RNA sequencing data. This pipeline aligns reads to the reference genome with STAR. Variant calling was carried out according to the Genome Analysis Toolkit (GATK) best practices for variant calling [27]. Annotation was carried out with Annovar [28]. We analyzed 549 cancer-related genes from the UCSF500 Cancer gene panel (accessed on 31 August 2023 on https://genomics.ucsf.edu/UCSF500) (see Supplementary Table S1). Mutations were required to have a read depth of at least 30 and a VAF > 10%. Population variants were excluded based on gnomAD and 1000G annotations (variants with allele frequency = 0.01 were kept for analysis). We carried out a further selection of the variants with the prediction tools UMD Predictor (probable pathogenic and pathogenic), PolyPhen 2 HumDiv (possibly damaging and probably damaging), Provean (damaging), and CADD score (>15).

2.4. Multiplex Ligation-Dependent Probe Amplification for CNV Detection

DNA was analyzed 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 (IKZF1, CDKN2A/B, PAX5, EBF1, ETV6, BTG1, and RB1), as well as in genes from the X/Y PAR1 region (CRLF2, CSF2RA, SHOX, IL3RA, and P2RY8). It also contains 13 reference probes that function as internal controls. MLPA reactions were carried out on a 96-well PCR thermocycler SimpliAmp Thermal Cycler (Applied Biosystems, Thermo Fisher, Waltham, MA, USA), and the products were separated by capillary electrophoresis on an ABI-3500 genetic analyzer (Applied Biosystems, Thermo Fisher, Waltham, MA, USA). The resulting peak intensities were analyzed using Coffalyser software (https://www.mrcholland.com/technology/software/coffalyser-net) (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.

2.5. Sanger Sequencing for Detection of Germline PAX5 Variants

The polymerase chain reaction (PCR) primers (see Supplementary Table S2) were designed according to the established laboratory protocol. The reference sequence for each amplicon was acquired from the ENSEMBL database. Optimal primers were selected by the Primer3 online tool (v. 0.4.0) [29] and analyzed for possible polymorphisms with the SNPcheck tool [30].

The primers from the PCR amplification reaction were removed in the first step with the use of the ExoSap-IT enzyme mix (Affymetrix, Santa Clara, CA, USA), followed by sequence reaction with the BigDye Terminator sequencing master mix (Applied Biosystems, Waltham, MA, USA) according to the established laboratory protocol. The sequencing was performed on the ABI-3500 genetic analyzer (Applied Biosystems, Thermo Fisher, Waltham, MA, USA). Sequences were aligned to the reference sequence (NG_033894.1) from the NCBI Reference Sequences and analyzed with the BLAST software 2.15.0 (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 31 August 2023).

2.6. Statistical Analysis

Comparisons of categorical values were carried out using the logistic regression and Fischer’s exact test, with the use of Bonferroni correction to adjust for multiple testing. For the comparison of numerical values, the Mann–Whitney U-test was used. Survival analysis was carried out with the Kaplan–Meier method, and Cox regression models were used for testing the significance of variables on the outcomes. To adjust for small sample sizes in survival analysis within individual genetic subtypes, we performed the permutations test on the log-rank statistic with the R package FHtest. The significance level for all the tests was 5% (p < 0.05 was considered to be statistically significant). An event was defined as failure to achieve remission, a relapse after remission, the development of a second malignant neoplasm, or death. Overall survival was defined as time to death. Both event-free survival (EFS) and overall survival (OS) were censored at last contact.

3. Results

3.1. PAX5alt and PAX5 P80R Patients

Among our cohort of 99 patients, we found eight patients with PAX5-associated genetic subtypes (Table 1). We identified six patients with the PAX5alt signature and one with the PAX5 P80R subtype with the tool ALLCatchR. The classifier MD-ALL identified one additional sample as PAX5alt which was classified as Ph-like by ALLCatchR. Manual inspection of the fusion transcripts and variants showed no compelling evidence for Ph-like genetic alterations. We did, however, identify a missense variant in PAX5 with computational prediction tools unanimously supporting a deleterious effect on the gene. Therefore, this sample was regarded as PAX5alt. For seven of these patients, cytogenetic results were available, and for all of them, RT-PCR data were available. All of these patients were previously classified as B-other as they did not exhibit any recurrent abnormalities (BCR::ABL1, ETV6::RUN1, TCF3::PBX1, KMT2A, hyperdiploidy, hypodiploidy). The PAX5alt subtype therefore accounted for 7.1% of our whole B-ALL cohort and 29.1% of the B-other group, while the one patient with PAX5 P80R represented 1.0% of the whole cohort and 4.2% of the B-other group. One sample with the BCR::ABL1 fusion gene was classified as PAX5alt by both prediction tools but was manually assigned to the BCR::ABL1 subgroup.

3.2. Genetic Changes Identified in PAX5alt and PAX5 P80R Patients

We found diverse genetic alterations in our patients with the PAX5alt gene expression profile (Table 1). We identified four patients with rearrangements in PAX5. Two patients carried the fusion gene PAX5::FKBP15, and two had the fusion gene PAX5::FOXP1. Furthermore, two patients carried a missense single nucleotide variant in PAX5. The first had a substitution in exon 2, c.101C>T, leading to an amino acid change from proline to leucine at position 34, which was classified as likely pathogenic according to the ACMG guidelines. The second carried a substitution in exon 6, c.725A>G, changing the glutamic acid at position 242 to glycine, classified as a variant of uncertain significance (VUS). One patient had a frameshift variant, a deletion of four nucleotides in exon 9, c.1045_1048del, leading to a stop codon after 53 amino acids, and a partial loss of the C-terminal domain. This variant was classified as likely pathogenic and explained the predicted PAX5alt subtype. Germline samples were available for these three patients, and we did not confirm the presence of these variants in their samples. We also identified one patient with the recurrent PAX5 P80R somatic variant without any CNVs that we investigated. We did, however, identify this patient as a compound heterozygote with three missense variants on the second allele that were not expressed. These were variants in exons 5 (c.596G>A, p.R199K), 7 (c.826G>A, p.D276N) and 10 (c.1102T>A, p.S368T), all classified as VUS.

MLPA results revealed that five of the seven PAX5alt patients carried additional deletions in PAX5, leading to the loss of heterozygosity, of which two patients had additional deletions present and were therefore classified as IKZF1[sup.plus]. One patient had an amplification in PAX5. All patients classified as PAX5alt also carried deletions in CDKN2A and CDKN2B.

Considering the other single nucleotide variants in these patients, we identified on average 17 additional cancer-related variants per patient in 43 different genes in patient samples with the PAX5alt and PAX5 P80R subtypes, the most commonly altered genes being ATM, ATRX, CYP2D6, CYP2D7, KDM5A, KMT2A, and NCOR1 (see Supplementary Table S3).

3.3. Clinical Characteristics and Outcomes for Patients with the PAX5alt and PAX5 P80R Subtypes

Most of the patients with the PAX5alt and PAX5 P80R subtypes were male (75%) and only one was initially classified as standard risk; however, the difference between this group and other B-ALL patients from our cohort was not statistically significant (Table 2). They also did not differ from other B-ALL cases in other clinical characteristics (WBC at diagnosis, blast cell count on day 8, MRD on days 15 and 33, CNS involvement, age at diagnosis). One patient (14.3%) in the PAX5alt group and seventeen patients (18.7%) in the rest of the B-ALL cohort experienced an event. When comparing the PAX5alt and PAX5 P80R patients to the rest of the cohort, event-free survival analysis also showed no differences between the groups (5-year EFS 88% vs. 84%, respectively). A similar trend was observed for overall survival, with one (14.3%) patient in the PAX5alt group and six (6.6%) patients in B-ALL other group that died (5-year OS 88% vs. 93%, respectively). However, when putting the EFS in the context of all of the B-ALL genetic subtypes in our cohort, we can see the patients with PAX5alt had intermediate outcomes compared to other subtypes (Figure 1).

Taking a closer look at the two patients with the PAX5::FKBP15 fusion, both were male and CNS negative, but presented with a high WBC at diagnosis. One presented with a complex karyotype with a small hypodiploid clone and was therefore classified as high risk. However, his MRD was negative on days 15 and 33, and the rest of the treatment was without complications despite the presence of the IKZF1[sup.plus] profile. On the other hand, the second patient was a good prednisone responder, but the day 15 MRD was positive, so he was classified as intermediate risk. MRD on day 33 was negative. MLPA results showed a partial deletion in PAX5. Later, this patient experienced a very early relapse with a monocytic switch after 13 months. This resulted in a rapid progression of the disease and death. Subsequent analysis showed that this patient carried a somatic pathogenic variant in PTPN11 (PTPN11:c.211T>C). We could not confirm the presence of this variant at the time of relapse as there was not sufficient material available.

3.4. PAX5 Variants in Other B-ALL Subtypes

One patient, who was otherwise classified as BCR::ABL1-positive, carried a secondary variant in PAX5 (c.197G>A, p.S66N, classified as VUS) alongside a deletion, leading to the loss of heterozygosity. This was the sample that was assigned to the PAX5alt signature by the two gene expression prediction tools. No other secondary PAX5 variants or gene fusions were identified in any of the patients that were not classified in the PAX5alt or PAX5 P80R subtypes.

3.5. PAX5 CNVs and Their Clinical Significance

For seven patients out of the entire cohort (7%), we did not have data from the MLPA analysis. Of the remaining 92, 31 (33.7%) carried CNVs in PAX5. In our cohort, neither deletions nor amplifications in PAX5 were notably predominant. Specifically, 17 of the patients with the CNVs were found to have deletions, while duplications were observed in 14 patients. The deletions were all heterozygous, the most common being the deletion of the whole gene, followed by different partial deletions. The most common were amplifications in exon 7 and/or 8 which were present in eight of the fourteen patients with duplications, followed by amplifications of the whole gene as a result of hyperdiploidy (four patients).

CNVs in PAX5 were significantly enriched in males compared to females (p = 0.015) but not in any specific genetic subtype despite more than half of the patients with the PAX5alt and BCR::ABL1 genetic profiles carrying them (Figure 2). Two subtypes stood out with a significantly lower incidence of PAX5 deletions, namely ETV6::RUNX1 (p = 0.011) and hyperdiploid cases (p = 0.005). The presence of any PAX5 CNVs was not shown to be associated with survival outcomes or other clinical parameters when looking at the whole cohort. We further investigated the duplications in PAX5 in hyperdiploid and ETV6::RUNX1-positive cases as they were abundantly present in these genetic subtypes (29.6% of hyperdiploid and 17% of ETV6::RUNX1 patients). Hyperdiploid cases carrying PAX5 duplications had significantly higher MRD values at the end of induction (p = 0.017) compared to those without the duplications. We also observed poorer event-free survival in hyperdiploid cases carrying PAX5 duplications compared to those without, with 5-year EFS of 75% vs. 100%, respectively (p = 0.024) (Figure 3); however, the same effect was not observed in patients with the ETV6::RUNX1 subtype (p = 0.371).

4. Discussion

Our aim was to investigate the frequency and types of genetic alterations in PAX5 in our cohort of 99 children with B-ALL. As expected, PAX5 was impacted in a high percentage of our patients (33.3%), with the genetic alterations ranging from CNVs and rearrangements to SNVs. The most common were CNVs, which were found in more than a third of the patients, followed by point mutations in 5.2%, and gene rearrangements in 4.1%.

The reports about the frequencies of the recurrent B-ALL genetic subtypes associated with PAX5 (PAX5alt and PAX5 P80R) show that they occur in 7.0–10.0% and 0.2–4.0% of childhood patients, respectively [7,31,32,33]. These frequencies were very similar in our cohort, with PAX5alt and PAX5 P80R representing 7.1% and 1.0% of all patients, respectively. With the current knowledge of these subtypes, there are conflicting reports regarding outcomes for these patients. While some reports state that both subtypes are associated with intermediate to favorable outcomes, the subject matter seems to be more diverse. In adult B-ALL cases, PAX5 P80R has been associated with favorable outcomes [34], while it has been shown to contribute to poorer outcomes in childhood cases [31]. The PAX5alt subtype is associated with intermediate outcomes [35], and patients with this subtype are reported to be classified as high risk more often than standard risk, according to National Cancer Institute criteria [36]. The event-free and overall survival of the patients with alterations in PAX5 did not differ from the rest of the patients in our cohort (Figure 1), with the 5-year EFS being 88% and 84%, respectively, and 5-year OS being 88% and 93%, respectively. When looking at the EFS within the distinct genetic subtypes, however, we notice that PAX5alt patients seem to have intermediate outcomes in our cohort as well, and they were in most cases classified as non-SR.

CNVs in PAX5 are present in around 30% of all B-ALL patients [5], the most common being heterozygous deletions leading to haploinsufficiency, as was also the case in our cohort. Confirming the high incidence of these CNVs, we also showed that they are more often found in males compared to females but are not significantly associated with clinical outcomes in patients, corroborating the findings that they are not considered predictive in the context of MRD levels and are not independently indicative of poor prognosis in childhood B-ALL [37]. We took a closer look at the CNVs within some genetic subgroups, and showed a trend of poorer survival in hyperdiploid cases carrying duplications in PAX5 compared to other hyperdiploid cases. This poorer survival did not seem to be influenced by other factors such as sex, age, or the presence of other CNVs; therefore, it could be indicative of a more complex interplay of genetic factors within a genetic subtype that is otherwise considered standard risk. However, as our results arise from a small cohort, additional research should be conducted and validated on a larger cohort with the possibility to better account for confounding variables. Additionally, many of the duplication cases we had in this cohort were duplications of exons 7 and/or 8, and they were mostly present in ETV6::RUNX1-positive and hyperdiploid cases, which has not been described in other study cohorts. There have been reports of recurrent intragenic amplifications of PAX5 involving exon 2 or exon 5 or both of these exons [38,39], which are found in patients that lack other major cytogenetic abnormalities and are associated with poorer outcomes. Interestingly, we did not identify any cases like this in our cohort.

Pathogenic variants in PAX5 are less common and are found in around 7% of the patients [5,40]. The consequence of these variants is the creation of a non-functional protein that loses its ability to bind to DNA or to regulate transcription. These variants are rarely germline; however, some cases have been described [8,9,41]. We identified PAX5 variants in five (5.2%) patients, of which one carried the recurring variant p.P80R, three were classified as PAX5alt, and one was BCR::ABL1-positive. Among the PAX5alt patients, two had a missense variant and one a frameshift deletion. One of the missense variants (c.101C>T), which affects the paired box domain, has previously been described as somatic in B-ALL patients [42], and one case with this variant has been reported in intestinal adenocarcinoma [43]. The other missense variant (c.725A>G) in the partial homeodomain and the frameshift deletion (c.1045_1048del) in the C-terminal domain have not yet been described in literature or reported in the COSMIC database. One patient had a secondary variant in this gene (c.197G>A) with the deletion of the second allele, while also carrying the BCR::ABL1 fusion gene. This variant has been reported in one sporadic case of B-ALL [9]. This sample was also classified as PAX5alt according to the gene expression profiling classifiers, indicating a strong influence of the PAX5 alterations on the global gene expression and processes in the cell. None of the variants we found were present in the germline samples we tested, confirming their somatic origin.

In previous studies, most patients with PAX5 P80R were compound heterozygotes [31], and this was also true in our case. In addition to the variant c.239C>G, this patient harbored three missense variants on the second allele, which was not expressed. Two variants (c.596G>A and c.826G>A) have previously been reported in malignant melanoma, bladder cancer, and adenocarcinoma [44,45], while the third (c.1102T>A) has not been described in the literature yet. The PAX5 p.P80R variant has also been shown to be a somatic and not a germline event [31,46], as we also confirmed for our patient who experienced no notable complications during the whole treatment period.

Translocations involving the PAX5 gene are observed in approximately 2.5% of pediatric and 1% of adult patients with B-ALL [4,7]. These translocations often result in the formation of chimeric genes that encode proteins preserving the DNA-binding paired box domain and the nuclear localization signal of PAX5, but with C-terminal domains derived from the fusion partners. Various partner genes, including transcription factors, structural proteins, and signal transducers, have been identified as fusion partners with the PAX5 gene [7,12]. We found PAX5 rearrangements in four patients, all classified as PAX5alt. In two patients, the fusion partner was FKBP15, and in the other two, PAX5 fused to FOXP1, both of which have already been described as fusion genes in B-ALL [5,10,47,48,49]. None of the Ph-like patients in our cohort carried any translocations involving this gene despite PAX5 translocations being very common in this subgroup of patients [7,11,12].

Among our patients with altered PAX5, we also report an interesting case of a patient with the PAX5::FKBP15 fusion gene who underwent a monocytic switch. Several cases of a switch from ALL to AML have been described, with the most frequent subtypes that undergo the switch being DUX4-rearranged, PAX5 P80R, and ZNF384-rearranged [50,51]. To our knowledge, this is the first described case of a PAX5alt patient undergoing a monocytic switch, further showing the complexity and interplay of genetic factors in this disease. This patient carried a somatic pathogenic variant in PTPN11. This variant, associated with myeloproliferative diseases, was already present at the time of the first diagnosis of ALL. PTPN11 encodes the Shp2 non-receptor protein–tyrosine phosphatase implicated in several signaling pathways. Activating mutations in PTPN11 are most commonly associated with juvenile myelomonocytic leukemia but are not as well defined in other neoplasms. They have been reported at a rate of 6.6% in AML [52,53] and some cases of high hyperdiploid ALL [54], with this particular variant mainly described in JMML and some AML cases. It is important to note that one of the main functions of PAX5 is maintaining non-B-cell genes repressed during B-cell differentiation; therefore, B-ALL tumors lacking PAX5 function may express lineage promiscuity markers such as monocytes without implying a lineage switch [55]. Nevertheless, identifying pathogenic variants in genes associated with myeloproliferative diseases might, therefore, aid in the identification of patients who are at high risk of a monocytic switch.

5. Conclusions

In summary, this study contributes valuable insights into the presence, frequency, and prognostic implications of diverse PAX5 alterations in B-ALL patients, underscoring the complex interplay of genetic factors and their impact on patient prognosis. Our findings reveal a high frequency of PAX5 alterations, the most common being CNVs, which may hold additional prognostic information within specific genetic subtypes, followed by single nucleotide variants and fusion genes that are in most cases the oncogenic drivers, particularly in PAX5alt and PAX5 P80R subtypes. Recognizing these alterations within routine diagnostic frameworks is crucial, especially in the anticipation of forthcoming investigations aimed at tailoring therapeutic approaches for these patients to improve clinical outcomes.

Author Contributions

Conceptualization, M.D. and J.J.; methodology, K.C., N.K. and T.T.; software, B.J.B. and K.C.; formal analysis, K.C. and R.Š.; investigation, K.C., N.K., T.T., M.K. and T.P.; resources, M.D., J.J., M.K., T.P., L.K., B.J.B. and J.K.; data curation, K.C., N.K., M.K. and T.P.; writing—original draft preparation, K.C.; writing—review and editing, M.D., J.J., M.K., T.P., L.K., T.T., R.Š., B.J.B. and J.K.; visualization, K.C.; supervision, M.D., L.K. and J.J.; project administration, M.D. and K.C.; funding acquisition, M.D., J.J. and M.K. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Republic of Slovenia National Medical Ethics Committee on 10 March 2011 and 15 March 2022 (reference number KME 51/03/11 and KME 0120-17/2022/6).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients or parents to publish this paper or any other anonymous data and results.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. Hardcopies of all data and results are also available in patients’ files and the laboratory that conducted the tests. The data are not publicly available due to privacy and ethical restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

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.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers16061164/s1, Table S1: UCSF500 Cancer Gene Panel; Table S2: Primer sequences for exons of PAX5; Table S3: Variants identified in PAX5alt and PAX5 P80R patients from RNAseq data.

References

1. C. Cobaleda; A. Schebesta; A. Delogu; M. Busslinger Pax5: The Guardian of B Cell Identity and Function., 2007, 8,pp. 463-470. DOI: https://doi.org/10.1038/ni1454.

2. S. Gruenbacher; M. Jaritz; L. Hill; M. Schäfer; M. Busslinger Essential Role of the Pax5 C-Terminal Domain in Controlling B Cell Commitment and Development., 2023, 220,p. e20230260. DOI: https://doi.org/10.1084/jem.20230260.

3. K. Okuyama; T. Strid; J. Kuruvilla; R. Somasundaram; S. Cristobal; E. Smith; M. Prasad; T. Fioretos; H. Lilljebjörn; S. Soneji et al. PAX5 Is Part of a Functional Transcription Factor Network Targeted in Lymphoid Leukemia., 2019, 15, e1008280. DOI: https://doi.org/10.1371/journal.pgen.1008280.

4. K. Nebral; D. Denk; A. Attarbaschi; M. König; G. Mann; O.A. Haas; S. Strehl Incidence and Diversity of PAX5 Fusion Genes in Childhood Acute Lymphoblastic Leukemia., 2009, 23,pp. 134-143. DOI: https://doi.org/10.1038/leu.2008.306.

5. C.G. Mullighan; S. Goorha; I. Radtke; C.B. Miller; E. Coustan-Smith; J.D. Dalton; K. Girtman; S. Mathew; J. Ma; S.B. Pounds et al. Genome-Wide Analysis of Genetic Alterations in Acute Lymphoblastic Leukaemia., 2007, 446,pp. 758-764. DOI: https://doi.org/10.1038/nature05690.

6. J. Familiades; M. Bousquet; M. Lafage-Pochitaloff; M.-C. Béné; K. Beldjord; J. De Vos; N. Dastugue; E. Coyaud; S. Struski; C. Quelen et al. PAX5 Mutations Occur Frequently in Adult B-Cell Progenitor Acute Lymphoblastic Leukemia and PAX5 Haploinsufficiency Is Associated with BCR-ABL1 and TCF3-PBX1 Fusion Genes: A GRAALL Study., 2009, 23,pp. 1989-1998. DOI: https://doi.org/10.1038/leu.2009.135.

7. Z. Gu; M.L. Churchman; K.G. Roberts; I. Moore; X. Zhou; J. Nakitandwe; K. Hagiwara; S. Pelletier; S. Gingras; H. Berns et al. PAX5-Driven Subtypes of B-Progenitor Acute Lymphoblastic Leukemia., 2019, 51,pp. 296-307. DOI: https://doi.org/10.1038/s41588-018-0315-5. PMID: https://www.ncbi.nlm.nih.gov/pubmed/30643249.

8. N. Duployez; L.A. Jamrog; V. Fregona; C. Hamelle; L. Fenwarth; S. Lejeune; N. Helevaut; S. Geffroy; A. Caillault; A. Marceau-Renaut et al. Germline PAX5 Mutation Predisposes to Familial B-Cell Precursor Acute Lymphoblastic Leukemia., 2021, 137,pp. 1424-1428. DOI: https://doi.org/10.1182/blood.2020005756. PMID: https://www.ncbi.nlm.nih.gov/pubmed/33036026.

9. S. Shah; K.A. Schrader; E. Waanders; A.E. Timms; J. Vijai; C. Miething; J. Wechsler; J. Yang; J. Hayes; R.J. Klein et al. A Recurrent Germline PAX5 Mutation Confers Susceptibility to Pre-B Cell Acute Lymphoblastic Leukemia., 2013, 45,pp. 1226-1231. DOI: https://doi.org/10.1038/ng.2754. PMID: https://www.ncbi.nlm.nih.gov/pubmed/24013638.

10. Z. Li; S.H.R. Lee; W.H.N. Chin; Y. Lu; N. Jiang; E.H. Lim; E. Coustan-Smith; K.H. Chiew; B.L.Z. Oh; G.S. Koh et al. Distinct Clinical Characteristics of DUX4- and PAX5-Altered Childhood B-Lymphoblastic Leukemia., 2021, 5,pp. 5226-5238. DOI: https://doi.org/10.1182/bloodadvances.2021004895. PMID: https://www.ncbi.nlm.nih.gov/pubmed/34547766.

11. D. Schinnerl; K. Fortschegger; M. Kauer; J.R. M Marchante; R. Kofler; M.L. Den Boer; S. Strehl The Role of the Janus-Faced Transcription Factor PAX5-JAK2 in Acute Lymphoblastic Leukemia., 2015, 125,pp. 1282-1291. DOI: https://doi.org/10.1182/blood-2014-04-570960.

12. G. Fazio; S. Bresolin; D. Silvestri; M. Quadri; C. Saitta; E. Vendramini; B. Buldini; C. Palmi; M. Bardini; A. Grioni et al. PAX5 Fusion Genes Are Frequent in Poor Risk Childhood Acute Lymphoblastic Leukaemia and Can Be Targeted with BIBF1120., 2022, 83, 104224. DOI: https://doi.org/10.1016/j.ebiom.2022.104224.

13. V. Bhardwaj; S. Heyne; K. Sikora; L. Rabbani; M. Rauer; F. Kilpert; A.S. Richter; D.P. Ryan; T. Manke SnakePipes: Facilitating Flexible, Scalable and Integrative Epigenomic Analysis., 2019, 35,pp. 4757-4759. DOI: https://doi.org/10.1093/bioinformatics/btz436.

14. M. Martin Cutadapt Removes Adapter Sequences From High-Throughput Sequencing Reads., 2011, 17,pp. 10-12. DOI: https://doi.org/10.14806/ej.17.1.200.

15. A. Dobin; C.A. Davis; F. Schlesinger; J. Drenkow; C. Zaleski; S. Jha; P. Batut; M. Chaisson; T.R. Gingeras STAR: Ultrafast Universal RNA-Seq Aligner., 2013, 29,pp. 15-21. DOI: https://doi.org/10.1093/bioinformatics/bts635. PMID: https://www.ncbi.nlm.nih.gov/pubmed/23104886.

16. Y. Liao; G.K. Smyth; W. Shi FeatureCounts: An Efficient General Purpose Program for Assigning Sequence Reads to Genomic Features., 2014, 30,pp. 923-930. DOI: https://doi.org/10.1093/bioinformatics/btt656. PMID: https://www.ncbi.nlm.nih.gov/pubmed/24227677.

17. T. Beder; B.T. Hansen; A.M. Hartmann; J. Zimmermann; E. Amelunxen; N. Wolgast; W. Walter; M. Zaliova; Ž. Antic; P. Chouvarine et al. The Gene Expression Classifier ALLCatchR Identifies B-Cell Precursor ALL Subtypes and Underlying Developmental Trajectories Across Age., 2023, 7,p. E939. DOI: https://doi.org/10.1097/HS9.0000000000000939. PMID: https://www.ncbi.nlm.nih.gov/pubmed/37645423.

18. Z. Hu; Z. Jia; J. Liu; A. Mao; H. Han; Z. Gu MD-ALL: An Integrative Platform for Molecular Diagnosis of B-Acute Lymphoblastic Leukemia., 2023, DOI: https://doi.org/10.3324/haematol.2023.283706. PMID: https://www.ncbi.nlm.nih.gov/pubmed/37981856.. early view

19. A. Renevey; M. Proks; M.U. Garcia; A. Peltzer; R.S. Praveen; O. Ziff; P. Ewels; M. De Smet; B. Moran; C. Saad et al. Nf-Core/Rnafusion: 2.3.4 (2.3.4)., 2023, DOI: https://doi.org/10.5281/zenodo.7866954.

20. S. Uhrig; J. Ellermann; T. Walther; P. Burkhardt; M. Fröhlich; B. Hutter; U.H. Toprak; O. Neumann; A. Stenzinger; C. Scholl et al. Accurate and Efficient Detection of Gene Fusions from RNA Sequencing Data., 2021, 31,pp. 448-460. DOI: https://doi.org/10.1101/gr.257246.119.

21. P. Melsted; S. Hateley; I.C. Joseph; H. Pimentel; N. Bray; L. Pachter Fusion Detection and Quantification by Pseudoalignment., 2017, DOI: https://doi.org/10.1101/166322.

22. B.J. Haas; A. Dobin; N. Stransky; B. Li; X. Yang; T. Tickle; A. Bankapur; C. Ganote; T.G. Doak; N. Pochet et al. STAR-Fusion: Fast and Accurate Fusion Transcript Detection from RNA-Seq., 2017, DOI: https://doi.org/10.1101/120295.

23. C. Ma; M. Shao; C. Kingsford SQUID: Transcriptomic Structural Variation Detection from RNA-Seq., 2018, 19, 52. DOI: https://doi.org/10.1186/s13059-018-1421-5.

24. D. Nicorici; M. Satalan; H. Edgren; S. Kangaspeska; A. Murumägi; O. Kallioniemi; S. Virtanen; O. Kilkku FusionCatcher—A Tool for Finding Somatic Fusion Genes in Paired-End RNA-Sequencing Data., 2014, 011650. DOI: https://doi.org/10.1101/011650.

25. L. Tian; Y. Li; M.N. Edmonson; X. Zhou; S. Newman; C. McLeod; A. Thrasher; Y. Liu; B. Tang; M.C. Rusch et al. CICERO: A Versatile Method for Detecting Complex and Diverse Driver Fusions Using Cancer RNA Sequencing Data., 2020, 21,pp. 1-18. DOI: https://doi.org/10.1186/s13059-020-02043-x. PMID: https://www.ncbi.nlm.nih.gov/pubmed/32466770.

26. R.S. Praveen; M.U. Garcia; K. Menden Nf-Core/Rnavar: Nf-Core/Rnavar v1.0.0—Tiny Firefly (1.0.0)., 2022, DOI: https://doi.org/10.5281/zenodo.6669637.

27. G.A. Van der Auwera; M.O. Carneiro; C. Hartl; R. Poplin; G. del Angel; A. Levy-Moonshine; T. Jordan; K. Shakir; D. Roazen; J. Thibault et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline., 2013, 43,pp. 11.10.1-11.10.33. DOI: https://doi.org/10.1002/0471250953.bi1110s43.

28. K. Wang; M. Li; H. Hakonarson ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data., 2010, 38,p. e164. DOI: https://doi.org/10.1093/nar/gkq603. PMID: https://www.ncbi.nlm.nih.gov/pubmed/20601685.

29. S. Rozen; H. Skaletsky Primer3 on the WWW for General Users and for Biologist Programmers., 2000, 132,pp. 365-386. DOI: https://doi.org/10.1385/1-59259-192-2:365.

30. B. Yu; C. Zhang In Silico PCR Analysis., 2011, 760,pp. 91-107. DOI: https://doi.org/10.1007/978-1-61779-176-5_6.

31. M. Jung; M. Schieck; W. Hofmann; M. Tauscher; J. Lentes; A. Bergmann; M. Stelter; A. Möricke; J. Alten; B. Schlegelberger et al. Frequency and Prognostic Impact of PAX5 p.P80R in Pediatric Acute Lymphoblastic Leukemia Patients Treated on an AIEOP-BFM Acute Lymphoblastic Leukemia Protocol., 2020, 59,pp. 667-671. DOI: https://doi.org/10.1002/gcc.22882.

32. N.A. Migita; P.Y. Jotta; N.P. do Nascimento; V.S. Vasconcelos; G.L. Centoducatte; K.B. Massirer; A.C. de Azevedo; S.R. Brandalise; J.A. Yunes Classification and Genetics of Pediatric B-Other Acute Lymphoblastic Leukemia by Targeted RNA Sequencing., 2023, 7,pp. 2957-2971. DOI: https://doi.org/10.1182/bloodadvances.2022009179. PMID: https://www.ncbi.nlm.nih.gov/pubmed/36848637.

33. J.F. Li; Y.T. Dai; H. Lilljebjörn; S.H. Shen; B.W. Cui; L. Bai; Y.F. Liu; M.X. Qian; Y. Kubota; H. Kiyoi et al. Transcriptional Landscape of B Cell Precursor Acute Lymphoblastic Leukemia Based on an International Study of 1223 Cases., 2018, 115,pp. E11711-E11720. DOI: https://doi.org/10.1073/pnas.1814397115.

34. M. Passet; N. Boissel; F. Sigaux; C. Saillard; M. Bargetzi; I. Ba; X. Thomas; C. Graux; Y. Chalandon; T. Leguay et al. PAX5 P80R Mutation Identifies a Novel Subtype of B-Cell Precursor Acute Lymphoblastic Leukemia with Favorable Outcome., 2019, 133,pp. 280-284. DOI: https://doi.org/10.1182/blood-2018-10-882142. PMID: https://www.ncbi.nlm.nih.gov/pubmed/30510083.

35. S. Jeha; J. Choi; K.G. Roberts; D. Pei; E. Coustan-Smith; H. Inaba; J.E. Rubnitz; R.C. Ribeiro; T.A. Gruber; S.C. Raimondi et al. Clinical Significance of Novel Subtypes of Acute Lymphoblastic Leukemia in the Context of Minimal Residual Disease-Directed Therapy., 2021, 2,pp. 326-337. DOI: https://doi.org/10.1158/2643-3230.BCD-20-0229. PMID: https://www.ncbi.nlm.nih.gov/pubmed/34250504.

36. I. Iacobucci; S. Kimura; C.G. Mullighan Clinical Medicine Biologic and Therapeutic Implications of Genomic Alterations in Acute Lymphoblastic Leukemia., 2021, 10, 3792. DOI: https://doi.org/10.3390/jcm10173792.

37. E.M.P. Steeghs; J.M. Boer; A.Q. Hoogkamer; A. Boeree; V. de Haas; H.A. de Groot-Kruseman; M.A. Horstmann; G. Escherich; R. Pieters; M.L. den Boer Copy Number Alterations in B-Cell Development Genes, Drug Resistance, and Clinical Outcome in Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia., 2019, 9,p. 4634. DOI: https://doi.org/10.1038/s41598-019-41078-4.

38. A.H. Pour Feizi; S. Zeinali; J. Toporski; R. Sheervalilou; S. Mehranfar Frequency and Correlation of Common Genes Copy Number Alterations in Childhood Acute Lymphoblastic Leukemia with Prognosis., 2020, 21,pp. 3493-3500. DOI: https://doi.org/10.31557/APJCP.2020.21.12.3493.

39. C. Schwab; K. Nebral; L. Chilton; C. Leschi; E. Waanders; J.M. Boer; M. Žaliová; R. Sutton; I.I. Öfverholm; K. Ohki et al. Intragenic Amplification of PAX5: A Novel Subgroup in B-Cell Precursor Acute Lymphoblastic Leukemia?., 2017, 1,pp. 1473-1477. DOI: https://doi.org/10.1182/bloodadvances.2017006734.

40. J. Li; Y. Dai; L. Wu; M. Zhang; W. Ouyang; J. Huang; S. Chen Emerging Molecular Subtypes and Therapeutic Targets in B-Cell Precursor Acute Lymphoblastic Leukemia., 2021, 15,pp. 347-371. DOI: https://doi.org/10.1007/s11684-020-0821-6.

41. F. Auer; F. Rüschendorf; M. Gombert; P. Husemann; S. Ginzel; S. Izraeli; M. Harit; M. Weintraub; O.Y. Weinstein; I. Lerer et al. Inherited Susceptibility to Pre B-ALL Caused by Germline Transmission of PAX5 c.547G>A., 2014, 28,pp. 1136-1138. DOI: https://doi.org/10.1038/leu.2013.363.

42. S.W. Brady; K.G. Roberts; Z. Gu; L. Shi; S. Pounds; D. Pei; C. Cheng; Y. Dai; M. Devidas; C. Qu et al. The Genomic Landscape of Pediatric Acute Lymphoblastic Leukemia., 2022, 54,pp. 1376-1389. DOI: https://doi.org/10.1038/s41588-022-01159-z. PMID: https://www.ncbi.nlm.nih.gov/pubmed/36050548.

43. T.-M. Kim; S.-H. Jung; M.S. Kim; I.-P. Baek; S.-W. Park; S.H. Lee; H.H. Lee; S.S. Kim; Y.-J. Chung; S.H. Lee The Mutational Burdens and Evolutionary Ages of Early Gastric Cancers Are Comparable to Those of Advanced Gastric Cancers., 2014, 234,pp. 365-374. DOI: https://doi.org/10.1002/path.4401.

44. M. Krauthammer; Y. Kong; A. Bacchiocchi; P. Evans; N. Pornputtapong; C. Wu; J.P. McCusker; S. Ma; E. Cheng; R. Straub et al. Exome Sequencing Identifies Recurrent Mutations in NF1 and RASopathy Genes in Sun-Exposed Melanomas., 2015, 47,pp. 996-1002. DOI: https://doi.org/10.1038/ng.3361. PMID: https://www.ncbi.nlm.nih.gov/pubmed/26214590.

45. E.M. Van Allen; K.W. Mouw; P. Kim; G. Iyer; N. Wagle; H. Al-Ahmadie; C. Zhu; I. Ostrovnaya; G.V. Kryukov; K.W. O’connor et al. Somatic ERCC2 Mutations Correlate with Cisplatin Sensitivity in Muscle-Invasive Urothelial Carcinoma., 2014, 4,pp. 1140-1153. DOI: https://doi.org/10.1158/2159-8290.CD-14-0623.

46. L. Bastian; M.P. Schroeder; C. Eckert; C. Schlee; J.O. Tanchez; S. Kämpf; D.L. Wagner; V. Schulze; K. Isaakidis; J. Lázaro-Navarro et al. PAX5 Biallelic Genomic Alterations Define a Novel Subgroup of B-Cell Precursor Acute Lymphoblastic Leukemia., 2019, 33,pp. 1895-1909. DOI: https://doi.org/10.1038/s41375-019-0430-z. PMID: https://www.ncbi.nlm.nih.gov/pubmed/30842609.

47. N. Put; D. Deeren; L. Michaux; P. Vandenberghe FOXP1 and PAX5 Are Rare but Recurrent Translocations Partners in Acute Lymphoblastic Leukemia., 2011, 204,pp. 462-464. DOI: https://doi.org/10.1016/j.cancergen.2011.07.006.

48. A. Oikonomou; L. Valsecchi; M. Quadri; T. Watrin; K. Scharov; S. Procopio; J.-W. Tu; M. Vogt; A.M. Savino; D. Silvestri et al. High-Throughput Screening as a Drug Repurposing Strategy for Poor Outcome Subgroups of Pediatric B-Cell Precursor Acute Lymphoblastic Leukemia., 2023, 217, 115809. DOI: https://doi.org/10.1016/j.bcp.2023.115809.

49. Y. Kaminskiy; V. Kuznetsova; A. Kudriaeva; E. Zmievskaya; E. Bulatov Neglected, yet Significant Role of FOXP1 in T-Cell Quiescence, Differentiation and Exhaustion., 2022, 13,p. 971045. DOI: https://doi.org/10.3389/fimmu.2022.971045.

50. M. Novakova; M. Zaliova; K. Fiser; B. Vakrmanova; L. Slamova; A. Musilova; M. Brüggemann; M. Ritgen; E. Fronkova; T. Kalina et al. DUX4r, ZNF384r and PAX5-P80R Mutated B-Cell Precursor Acute Lymphoblastic Leukemia Frequently Undergo Monocytic Switch., 2021, 106,pp. 2066-2075. DOI: https://doi.org/10.3324/haematol.2020.250423.

51. M. Novakova; B. Vakrmanova; L. Slamova; A. Musilova; M. Brüggemann; M. Ritgen; E. Fronkova; T. Kalina; J. Trka; J. Stary et al. Switching Towards Monocytic Lineage and Discordancy between Flow Cytometric and PCR Minimal Residual Disease Results Is a Hallmark Feature of DUX4 Rearranged B-Cell Precursor Acute Lymphoblastic Leukemia., 2018, 132,p. 2825. DOI: https://doi.org/10.1182/blood-2018-99-116575.

52. L. Chen; W. Chen; M. Mysliwski; J. Serio; J. Ropa; F.A. Abulwerdi; R.J. Chan; J.P. Patel; M.S. Tallman; E. Paietta et al. Mutated Ptpn11 Alters Leukemic Stem Cell Frequency and Reduces the Sensitivity of Acute Myeloid Leukemia Cells to Mcl1 Inhibition., 2015, 29,pp. 1290-1300. DOI: https://doi.org/10.1038/leu.2015.18. PMID: https://www.ncbi.nlm.nih.gov/pubmed/25650089.

53. M. Tartaglia; S. Martinelli; I. Iavarone; G. Cazzaniga; M. Spinelli; E. Giarin; V. Petrangeli; C. Carta; R. Masetti; M. Aricò et al. Somatic PTPN11 Mutations in Childhood Acute Myeloid Leukaemia., 2005, 129,pp. 333-339. DOI: https://doi.org/10.1111/j.1365-2141.2005.05457.x. PMID: https://www.ncbi.nlm.nih.gov/pubmed/15842656.

54. A.J. De Smith; J. Ojha; S.S. Francis; E. Sanders; A.A. Endicott; H.M. Hansen; I. Smirnov; A.M. Termuhlen; K.M. Walsh; C. Metayer et al. Clonal and Microclonal Mutational Heterogeneity in High Hyperdiploid Acute Lymphoblastic Leukemia., 2016, 7,pp. 72733-72745. DOI: https://doi.org/10.18632/oncotarget.12238.

55. A. Delogu; A. Schebesta; Q. Sun; K. Aschenbrenner; T. Perlot; M. Busslinger Gene Repression by Pax5 in B Cells Is Essential for Blood Cell Homeostasis and Is Reversed in Plasma Cells., 2006, 24,pp. 269-281. DOI: https://doi.org/10.1016/j.immuni.2006.01.012.

Figures and Tables

Figure 1: Event-free survival for patients with PAX5alt vs. other B-ALL subtypes in our cohort (PAX5alt subtype is marked in orange). [Please download the PDF to view the image]

Figure 2: Frequency of PAX5 CNVs in B-ALL genetic subtypes. [Please download the PDF to view the image]

Figure 3: Event-free survival for hyperdiploid patients with PAX5 duplications vs. those without PAX5 duplications. [Please download the PDF to view the image]

Table 1: Genetic and clinical characteristics of B-ALL patients with a PAX5-associated genetic subtype.
Patient Nr.
12345678


Predicted subtype


PAX5alt


PAX5alt


PAX5alt


PAX5alt


PAX5alt


PAX5alt


PAX5alt


PAX5 P80R


Genetic variants in PAX5


/


/


/


c.1045_1048del


c.101C>T


/


c.725A>G


[c.239C>G]; [c.596G>A, c.826G>A, c.1102T>A]


Protein change


/


/


/


p.P349Sfs*53


p.P34L


/


p.E242G


p.P80R, p.R199K, p.D276N, p.S368T


ACMG classification


/


/


/


LP


LP


/


VUS


LP, VUS, VUS, VUS


Fusion genes involving PAX5


PAX5::FKBP15


PAX5::FOXP1


PAX5::FOXP1


/


/


PAX5::FKBP15


/


/


CNVs


IKZF1 del, CDKN2A/2B del., PAX5 ex. 10 del. (IKZF1[sup.plus])


CDKN2A/2B del.


PAX5 del., CDKN2A/2B del.


PAX5 del., CDKN2A/2B del.


PAX5 ex. 8 dup., CDKN2A/2B del.


PAX5 ex. 8 and 10 del., CDKN2A/2B del.


IKZF1 del., PAX5 del., CDKN2A/2B del. (IKZF1[sup.plus])


/


Gender


male


female


female


male


male


male


male


male


Age at diagnosis (years)


4


7


4


1


19


3


3


8


WBC at diagnosis (×10[sup.9]/L)


22.6


8.7


1.1


17.2


5


32.2


372.3


1.98


CNS involvement


neg.


neg.


neg.


pos.


neg.


neg.


neg.


neg.


Day 8 prednisone response (blasts/mm[sup.3])


280


1750


45


78


0


0


1078


0


MRD day 15 (%)


0


50


0.019


0.6


0.01


1.17


0.43


0


MRD day 33 (%)


0


0


0


0.016


0


0


0


0


End of induction risk classification


HR


HR


SR


IR


IR


IR


IR


IR


Event


no


no


no


no


no


relapse with AML switch, death


no


no


(LP—likely pathogenic; VUS—variant of uncertain significance; ex.—exon; del.—deletion; dup.—duplication; neg.—negative; pos.—positive; SR—standard risk; IR—intermediate risk; HR—high risk; CNVs—copy number variations; WBC—white blood count; CNS—central nervous system; MRD—minimal residual disease).

Table 2: Patients’ clinical characteristics.


Patient Characteristics


PAX5alt or PAX5 P80R (n = 8)


Other (n = 91)


Gender


male


6


48


female


2


43


Age at diagnosis


?1 year


0


3


1–5 years


5


52


=6 years


3


36


WBC at diagnosis


?20 × 10[sup.9]/L


5


71


=20 × 10[sup.9]/L


3


20


CNS involvement


yes


1


6


no


7


83


unknown


0


2


Day 8 prednisone response


?1000 blasts/µL


6


79


=1000 blasts/µL


2


10


unknown


0


2


MRD day 15


?0.1%


4


31


0.1–10%


3


45


?10%


1


12


unknown


0


3


MRD day 33


?0.01%


7


66


=0.01%


1


22


unknown


0


3


End of induction risk classification


standard


1


16


intermediate


5


54


high


2


21


Author Affiliation(s):

[1] Clinical Institute for Special Laboratory Diagnostics, University Children’s Hospital, University Medical Centre Ljubljana, Vrazov trg 1, 1000 Ljubljana, Slovenia; [email protected] (K.C.);

[2] Department of Oncology and Haematology, University Children’s Hospital, University Medical Centre Ljubljana, Bohoriceva ulica 20, 1000 Ljubljana, Slovenia

[3] Faculty of Medicine, University of Ljubljana, Vrazov trg 2, 1000 Ljubljana, Slovenia

Author Note(s):

[*] Correspondence: [email protected]

DOI: 10.3390/cancers16061164
COPYRIGHT 2024 MDPI AG
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2024 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Crepinšek, Klementina; Klobucar, Nika; Tesovnik, Tine; Šket, Robert; Jenko Bizjan, Barbara; Kovac, J
Publication:Cancers
Geographic Code:4EXSL
Date:Mar 1, 2024
Words:7954
Previous Article:Pretreatment Circulating HPV16 DNA Viral Load Predicts Risk of Distant Metastasis in Patients with HPV16-Positive Oropharyngeal Cancer.
Next Article:Exploratory Assessment of Galectin-1, -3, and -9 in Non-Small Cell Lung Cancer.
Topics:

Terms of use | Privacy policy | Copyright © 2024 Farlex, Inc. | Feedback | For webmasters |