  It is presented the general method that allows to formulate 4D $SU(N)$ Yang -
Mills theory in terms of only local gauge invariant variables. For the case
N=2, that is discussed in details, this gauge invariant formulation appears to
be very similar to $ R^2 $-gravity.
