A UV complete model where the Dark Matter (DM) particle interacts with gluons via a colored scalar mediator provides a viable phenomenological model that can be tested at hadron colliders. While monojet signatures are relevant for Collider searches, zero-jet processes would mean complete annihilation of Standard Model (SM) particles to DM particles, which contribute to the relic-density of DM. We look at the DM annihilation to SM colored particles which, at leading order, is a loop induced process in our model. We compute two-loop amplitudes in QCD that contribute to the process. We consider both gluon and quark channels for annihilation. Decomposing the amplitude in terms of form factors and making use of the projector technique, scalar Feynman integrals are obtained. Further, with the help of the IBP identities, an analytical expression for amplitude is obtained in terms of Master Integrals. The amplitude is made UV finite by Counterterm Renormalization. Our results can be used to predict DM pair production at hadron colliders and DM annihilation to SM colored particles at next-to-leading order in QCD.