Myelodysplastic syndromes (MDS) arise in older adults through stepwise acquisitions of multiple somatic mutations. Here, analyzing 1809 MDS patients, we infer clonal architecture by using a stringent, the single-cell sequencing validated PyClone bioanalytic pipeline, and assess the position of the mutations within the clonal architecture. All 3,971 mutations are grouped based on their rank in the deduced clonal hierarchy (dominant and secondary). We evaluated how they affect the resultant morphology, progression, survival and response to therapies. Mutations of SF3B1, U2AF1, and TP53 are more likely to be dominant, those of ASXL1, CBL, and KRAS are secondary. Among distinct combinations of dominant/secondary mutations we identified 37 significant relationships, of which 12 affect clinical phenotypes, 5 cooperatively associate with poor prognosis. They also predict response to hypomethylating therapies. The clonal hierarchy has distinct ranking and the resultant invariant combinations of dominant/secondary mutations yield novel insights into the specific clinical phenotype of MDS.