In the BRST formalism of gauge theories, the Lautrup-Nakanishi field $B^a(x)$ appears as an auxiliary variable $$\mathcal{L}_\text{BRST}=-\frac{1}{4}F_{\mu\nu}^a F^{a\,\mu\nu}+\frac{1}{2}\xi B^a B^a + B^a\partial_\mu A^{a\,\mu}+\partial_\mu\bar\eta^a(D^\mu\eta)^a,$$ and in the superfield formalism of SUSY, the field $F(x)$ also appears as an auxiliary variable: $$\mathcal{L}_\text{SUSY}=\partial_\mu \phi\partial^\mu\phi+i\bar\psi^\dagger\bar\sigma^\mu\partial_\mu\psi+F^*F+\ldots\,.$$
It is very tempting to view $B^a(x)$ and $F$ as Lagrange multipliers since their equations of motion leads to constraints. But, these variables do not enter into the Lagrangian linearly, like a conventional Lagrange multiplier. Rather, they enter into the Lagrangian quadratically.
However, in Kugo and Ojima's paper Manifestly Covariant Canonical Formulation of the Yang-Mills Field Theores (1978), they refer to the $B^a(x)$ fields as the 'Lagrange Multiplier' fields (p.1882).
So my question is: Can these auxiliary fields be viewed as Lagrange multipliers? and in what ways do they behave differently/similar to the conventional Lagrange multipliers that enter into the function linearly?