ABSTRACT
BACKGROUND: A substantial number of survivors of childhood acute lymphoblastic leukemia (ALL) suffer from treatment-related late adverse effects. While multiple studies have identified the effects of chemotherapeutics and radiation therapy on musculoskeletal outcomes, few have investigated their associations with genetic factors. METHODS: Here we analyzed musculoskeletal complications in relation to common and rare genetic variants derived through whole-exome sequencing of the PETALE cohort. Top-ranking associations were further assessed through stratified and multivariate analyses. RESULTS: DUOX2 variant was associated with skeletal muscle function deficit, as defined by peak muscle power Z score ≤ -2 SD (P = 4.5 × 10-5 for genotyping model). Upon risk stratification analysis, common variants in the APOL3, COL12A1, and LY75 genes were associated with Z score ≤ -2 SD at the cross-sectional area (CSA) at 4% radial length and lumbar bone mineral density (BMD) in high-risk patients (P ≤ 0.01). The modulation of the effect by risk group was driven by the interaction of the genotype with cumulative glucocorticoid dose. Identified variants remained significant throughout multivariate analyses incorporating non-genetic factors of the studied cohort. CONCLUSION: This exploratory study identified novel genetic variants associated with long-term musculoskeletal impairments in childhood ALL survivors. Replication in an independent cohort is needed to confirm the association found in this study.