Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative neoplasm (MPN) of childhood with a poor prognosis. Mutations in NF1, NRAS, KRAS, PTPN11 and CBL occur in 85% of patients, yet there are currently no risk stratification algorithms capable of predicting which patients will be refractory to conventional treatment and therefore be candidates for experimental therapies. In addition, there have been few other molecular pathways identified aside from the Ras/MAPK pathway to serve as the basis for such novel therapeutic strategies. We therefore sought to genomically characterize serial samples from patients at diagnosis through relapse and transformation to acute myeloid leukemia in order to expand our knowledge of the mutational spectrum in JMML. We identified recurrent mutations in genes involved in signal transduction, gene splicing, the polycomb repressive complex 2 (PRC2) and transcription. Importantly, the number of somatic alterations present at diagnosis appears to be the major determinant of outcome.