We examine the formation of clusters of galaxies in numerical simulations of a QUMOND cosmogony with massive sterile neutrinos. Clusters formed in these exploratory simulations develop higher velocities than those found in {Lambda}CDM simulations. The bulk motions of clusters attain about 1000 km/s by low redshift, comparable to observations whereas {Lambda}CDM simulated clusters tend to fall short. Similarly, high pairwise velocities are common in cluster-cluster collisions like the Bullet cluster. There is also a propensity for the most massive clusters to be larger in QUMOND and to appear earlier than in {Lambda}CDM, potentially providing an explanation for pink elephants like El Gordo. However, it is not obvious that the cluster mass function can be recovered.