Does the Slavnov-Taylor identity still hold for scalar Yang-Mills?

I want to renormalize the minimally-coupled scalar Yang-Mills theory: $$\mathcal{L}_{YM\phi}=(D_\mu\phi)^\dagger(D^\mu\phi)-\frac{1}{4}F_{\mu\nu}^a{F^{\mu\nu}}^a-\frac{1}{2\xi}(\partial_\mu {A^\mu}^a)^2-\bar{c}^a\partial^\mu D_\mu^{ac}c^c$$ to find the $$\beta(g)$$ function, and I would like to do so using the vector-scalar-scalar vertex. In a Yang-Mills theory coupled to fermions, I can safely use the vector-fermion-fermion vertex knowing that choosing other vertices would yield the same solution because of the Slavnov-Taylor identity, obtained with BRST symmetry. Now, are there analogous "BRST" symmetry (and corresponding "Slavnov-Taylor" identity) for the scalar case?