  We give a simple proof of L^p boundedness of iterated commutators of Riesz
transforms and a product BMO function. We use a representation of the Riesz
transforms by means of simple dyadic operators - dyadic shifts - which in turn
reduces the estimate quickly to paraproduct estimates.
