In supersymmetric gauge theories, the vector potential is a part of a vector supermultiplet which is represented by a real superfield $V$. Expanded out in components, the Lagrangian for such a field is given by
$$\mathcal{L}=-\frac{1}{4g^2}F_{\mu\nu}^a F^{\mu\nu a}+\cdots,$$
where the extra terms are in terms of the extra superpartners in this multiplet.
Upon gauge fixing, one introduces Faddeev-Popov ghosts to counter the overcounting in the path integral. The gauge-fixed Lagrangian in, say, the Lorenz gauge takes the form
$$\mathcal{L}_{\text{g.f.}}=\text{Tr}\left[-\frac{1}{4g^2}F_{\mu\nu} F^{\mu\nu}+\overline{c}D_{\mu}\partial^{\mu}c+\cdots\right].$$
Thus, since the vector potential $A$ transforms non-trivially under SUSY transformations, but the ghosts do not, it is clear that the introduction of ghost fields has broken supersymmetry at ghost level $n_g\neq 0$.
Is there a gauge-fixing method which circumvents this and allows for intrinsically supersymmetric perturbative calculations (for instance introducing a super partner for the ghosts and absorbing the broken supersymmetry into a $\mathbb{Z}_2$-graded version of extended BRST symmetry)?
