Diseases caused by single-gene mutations can display substantial phenotypic variability, which may be due to genetic, environmental, or epigenetic modifiers. Here, we induce Gaucher disease (GD), a rare inherited metabolic disorder, by injecting 15 inbred mouse strains with a low dose of a chemical inhibitor of acid β-glucosidase, the enzyme defective in GD. Different mouse strains exhibit widely different lifespans, which is unrelated to levels of acid β-glucosidase’s substrate accumulation. Genome-wide association reveals a number of candidate risk loci, including a marker within Grin2b, which in combination with another marker allows us to predict the lifespan of additional mouse strains. An antagonist of the NMDA receptor (encoded by Grin2b) significantly increases the lifespan of GD mice that would otherwise have lived for a short time. Our data identify putative modifier genes that may be involved in determining GD severity, which might help elucidate phenotypic variability between patients with similar GD mutations.