The $\Lambda$CDM simulations of Keller and Wadsley do not account for   the MOND mass-discrepancy-acceleration relation
  Keller and Wadsley (2016) have smugly suggested, recently, that the end of MOND may be in view. This is based on their claim that their highly-restricted sample of $\Lambda$CDM-simulated galaxies are "consistent" with the observed MOND mass-discrepancy-acceleration relation (MDAR), in particular, with its recent update by McGaugh et al. (2016), based on the SPARC sample. From this they extrapolate to "$\Lambda$CDM is fully consistent" with the MDAR. I explain why these simulated galaxies do not show that $\Lambda$CDM accounts for the MDAR. a. Their sample of simulated galaxies contains only 18 high-mass galaxies, within a narrow range of one order of magnitude in baryonic mass, at the very high end of the observed, SPARC sample, which spans 4.5 orders of magnitude in mass. More importantly, the simulated sample has none of the low-mass, low-acceleration galaxies -- abundant in SPARC -- which encapsulate the crux and the nontrivial aspects of the predicted and observed MDAR. The low-acceleration part of the simulated MDAR is achieved, rather trivially, from the flattish-velocity-curve regions of the simulated high-mass galaxies. b. Half of the simulated galaxies have "wrong" rotation curves that differ greatly from any observed ones. This, does not prevent these wrong galaxies from lying on the observed MDAR (for trivial reasons, again). They, in fact, define the high-acceleration branch of the simulated MDAR. c. To boot, even if $\Lambda$CDM were made "consistent" with the MDAR through the elaborate adjustments that the simulations allow, this would not obviate MOND, which predicts much more than the MDAR. 
Mordehai Milgrom
astro-ph.GA gr-qc