Next-to-leading order evolution of color dipoles

The small-x deep inelastic scattering in the saturation region is governed by the nonlinear evolution of Wilson-line operators. In the leading logarithmic approximation it is given by the Balitsky-Kovchegov equation for the evolution of color dipoles. In the next-to-leading order the Balitsky-Kovchegov equation gets contributions from quark and gluon loops as well as from the tree gluon diagrams with quadratic and cubic nonlinearities. We calculate the gluon contribution to the small-x evolution of Wilson lines (the quark part was obtained earlier)