I will present the full one-loop supersymmetric QCD corrections for gaugino annihilation and co-annihilation into light and heavy quarks in the Minimal Supersymmetric Standard Model (MSSM). I will show that these channels are phenomenologically relevant within the so-called phenomenological MSSM. Numerical results for the (co-)annihilation cross sections and the predicted neutralino relic density are presented. I will in particular show that the impact of including the radiative corrections on the cosmologically preferred region of the parameter space is larger than the current experimental uncertainty from Planck data. The presentation will include a few technical details as well as a discussion of the theoretical uncertainty due to scale variation in the calculation. Reference: arXiv:1404.2931 [hep-ph]