We give an overview of recent developments in the computation of the anomalous dimension matrix of composite operators in non-forward kinematics. The elements of this matrix determine the scale dependence of non-perturbative parton distributions, such as GPDs, and hence constitute important input for phenomenological studies of exclusive processes like deeply-virtual Compton scattering. Particular emphasis will be put on a recently developed method that exploits consistency relations for the anomalous dimension matrix which follow from the renormalization structure of the operators.