With the latest Planck results the dark matter relic density is determined to an unprecedented precision. In order to reduce current theoretical uncertainties in the dark matter relic density prediction, we have calculated next-to-leading order SUSY-QCD corrections to neutralino (co)annihilation processes including Coulomb enhancement effects. We demonstrate that these corrections can have significant impact on the cosmologically favoured MSSM parameter space and are thus of general interest for parameter studies and global fits.