Somatic mutations constitute key elements of the pathogenesis of myelodysplastic syndromes (MDS), a group of clonal hematologic neoplasms characterized by cytopenias, dysplasia and leukemic evolution. Whole exome sequencing followed by targeted deep sequencing in patients with MDS and related diseases has been performed cross-sectionally and serially. Bioinformatic analysis and confirmatory sequencing led to detection of in 1458 genes affected by somatic alterations, and identification of known and new driver events. For each patient, mutation spectrum as well as clonal hierarchy was determined and for each significantly mutated gene, its role in the clonal succession established. This approach allowed for a dynamic definition of MDS mutatome, including the spectrum of founding mutations and subsequent secondary mutational patterns. We demonstrate that certain founder events determine the mode and speed of disease progression, while secondary mutations may further modulate phenotypic features. Combinations of founder and secondary mutations further contribute to the phenotypic diversity but categorical grouping of cases based on the type of founder mutations may better define molecular subtypes of MDS and correlates with clinical parameters.