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: DUOX2variant was associated with skeletal muscle function deficit, as defined by peak muscle power Zscore ≤ −2 SD (P= 4.5 × 10−5for genotyping model). Upon risk stratification analysis, common variants in the APOL3, COL12A1, and LY75genes were associated with Zscore ≤ −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.