  We rederive the small-$x$ evolution equations governing quark helicity
distribution in a proton using solely an operator-based approach. In our
previous works on the subject, the evolution equations were derived using a mix
of diagrammatic and operator-based methods. In this work, we re-derive the
double-logarithmic small-$x$ evolution equations for quark helicity in terms of
the "polarized Wilson lines", the operators consisting of light-cone Wilson
lines with one or two non-eikonal local operator insertions which bring in
helicity dependence. For the first time we give explicit and complete
expressions for the quark and gluon polarized Wilson line operators, including
insertions of both the gluon and quark sub-eikonal operators. We show that the
double-logarithmic small-$x$ evolution of the "polarized dipole amplitude"
operators, made out of regular light-cone Wilson lines along with the polarized
ones constructed here, reproduces the equations derived in our earlier works.
The method we present here can be used as a template for determining the
small-$x$ asymptotics of any transverse momentum-dependent (TMD) quark (or
gluon) parton distribution functions (PDFs), and is not limited to helicity.
