I am stuck at a computation and I would appreciate any help. $U$ is the pion matrix in chiral perturbation theory

$$U=e^{i\sigma_a\phi_a/f}$$

where $\sigma_a$ are Pauli matrices, $\phi_a$ are three real scalar fields and $f$ is just a constant with mass dimensions. It is well known that this matrix is unitary, that is

$$U^{\dagger}U=I$$

Now comes the question. I want to compute this

$$Tr(\partial_{\mu}U^{\dagger}\partial_{\nu}U\partial^{\mu}U^{\dagger}\partial^{\nu}U-\partial_{\mu}U^{\dagger}\partial_{\nu}U\partial^{\nu}U^{\dagger}\partial^{\mu}U)$$

where $Tr$ denotes trace. I know the that answer is proportional to

$$Tr(\partial_{\mu}U^{\dagger}\partial_{\nu}U)Tr(\partial^{\mu}U^{\dagger}\partial^{\nu}U)$$

(if you want to know how I know this, it is claimed in page 9 of this http://arxiv.org/pdf/hep-ph/9403202v2.pdf paper). In any case, I am very puzzled here because a trace of something gives a product of traces. So, how does this follow?

EDIT::

I still haven't solved it but I think that equation (D.7) of this review http://arxiv.org/pdf/hep-ph/0210398v1.pdf might might help to solve it