An excess in gamma-rays emanating from the galactic centre has recently been observed in the Fermi-LAT data. This signal can be interpreted as resulting from WIMP annihilation, with the spectrum well-fit by dark matter annihilating dominantly into either bottom-quark or Higgs pairs. Supersymmetric models provide a well-motivated framework to study the implications of this signal in these channels. With a neutralino dark matter candidate, the gamma-ray excess cannot be easily accommodated in the minimal supersymmetric model, which in any case requires tuning below the percent level to explain the observed Higgs mass. Instead we are naturally led to consider the next-to-minimal model with a singlet superfield. This not only allows for the annihilation channel into bottom-quark pairs to be implemented, but also provides new possibilities for annihilation into Higgs-pseudoscalar pairs. We show that the fit to the gamma-ray excess for the Higgs-pseudoscalar channel can be just as good as for annihilation into bottom-quark pairs. Moreover, in the parameter range of interest, the next-to-minimal supersymmetric model solves the mu-problem and can explain the 125 GeV Higgs mass with improved naturalness. We also consider an extension by adding a right-handed neutrino superfield with the right-handed sneutrino acting as a dark matter candidate. Interestingly, this allows for the annihilation into pseudoscalar pairs which also provide a good fit to the gamma-ray excess. Furthermore, in the case of a neutralino LSP, the late decay of a sneutrino NLSP can non-thermally produce the observed relic abundance. Finally, the WIMP annihilation into scalar pairs allows for the possibility of detecting the Higgs or pseudoscalar decay into two photons, providing a smoking-gun signal of the model.