  Extending the work by Bukhvostov, Frolov, Lipatov and Kuraev (BFLK) on the
renormalization of quasipartonic operators we derive a complete set of
two-particle renormalization group kernels that enter QCD evolution equations
to twist-four accuracy. It is shown that the 2->2 evolution kernels which
involve ``non-partonic'' components of field operators, and, most remarkably,
also 2->3 kernels do not require independent calculation and can be restored
from the known results for quasipartonic operators using conformal symmetry and
Lorentz transformations. The kernels are presented for the renormalization of
light-ray operators built of chiral fields in a particular basis such that the
conformal symmetry is manifest. The results can easily be recast in momentum
space, in the form of evolution equations for generalized parton distributions.
