PSG and Other Candidate Genes as Potential Biomarkers of Therapy Resistance in B-ALL: Insights from Chromosomal Microarray Analysis and Machine Learning

Chromosomal microarray analysis (CMA) was performed for 40 patients with B-ALL undergoing treatment according to the ALL-2016 protocol to investigate the copy number alterations (CNAs) and copy neutral loss of heterozygosity (cnLOH) associated with minimal residual disease (MRD)-positive remission. Aberrations involving over 20,000 genes were identified, and a random forest approach was applied to isolate a subset of genes whose CNAs and cnLOH are significantly associated with poor therapeutic response. We have assembled the triple matched healthy population data and used that data as a reference, but not as a matched control. We identified a recurrent cluster of cnLOH in the 19q13.2-19q13.31 region, significantly enriched in MRD-positive patients (70% vs. 47% in the reference group vs. 16% in MRD-negative patients). This region includes the pregnancy-specific glycoprotein (<i>PSG</i>) gene family and the oncogene <i>ERF</i>, suggesting a potential role in leukemic persistence and treatment resistance. Additionally, we observed significant deletions involving 7p22.3 and 16q13, often as part of large-scale losses affecting almost the entire chromosomes 7 and 16, indicative of global chromosomal instability. These findings highlight specific genomic regions potentially involved in therapy resistance and may contribute to improved risk stratification in B-ALL. Our findings emphasize the value of high-resolution CMA in diagnostics and risk stratification and suggest that <i>PSG</i> genes and other candidate genes could serve as biomarkers for predicting treatment outcomes.