Busulfan (BU) is a key compound in conditioning myeloablative regimens for children undergoing hematopoietic stem cell transplantation (HSCT). There are wide interindividual differences in BU pharmacokinetics, which increase the risk of veno-occlusive disease, graft rejection and disease relapse. As BU is mainly metabolized by glutathione S-transferase (GST), it is hypothesized that functional polymorphisms in GST genes may explain in part the variability in BU pharmacokinetics. We analyzed polymorphisms in GSTA1 (C-69T, A-513G, G-631T, C-1142G), GSTM1 (deletion) and GSTP1 (A1578G, C2293T) genes in 28 children undergoing HSCT. All patients had individualized dosing based on pharmacokinetics after the first dose of intravenous BU. GSTM1-null individuals had higher drug exposure (P(Cmax)=0.008; P(AUC)=0.003; P(Css)=0.02) and lower clearance (P(CL)=0.001). Multivariate regression models showed that, other than the drug dose and age, the GSTM1 genotype was the best predictor of first-dose pharmacokinetic variability. GSTM1-null patients also received lower cumulative BU doses (P=0.02). No association was found between BU exposure and major GSTA1 or GSTP1 gene variants. In children, GSTM1 polymorphism seems to modify BU pharmacokinetics after intravenous drug administration.