Let $W^i_t$ and $W_t$ be pairwise independent Brownian motions for $i \in \{1, \dots , d\}$.
Let $X_t^i$ be $d$ independent Ornstein–Uhlenbeck processes for $i \in \{1, \dots , d\}$, i.e. each $X_t^i$ fulfills the SDE:
$$dX_t = \frac \gamma 2dW^i_t + \frac \alpha 2 X^i_t dt$$
To show: the following stochastic integrals are equally distributed:
$$\sum\int^T_0X^i_tdW^i_t \overset{d}= \int_0^T\sqrt{\sum (X^t_i)^2} dW_t$$
My attempt: if the integrands were deterministic, the Wiener integrals would be normally distributed with mean zero. However, the integrands are random.
One possibly useful observation is that the solution $X$ of the above Ornstein–Uhlenbeck equation is normally distributed with
$$X \sim \mathcal N \left(X_0e^{\,\alpha t / 2}, \frac {\gamma ^2 (e^{\alpha t} - 1)}{4\alpha} \right)$$
I thought I could apply Ito's isometry at some point, though in this situation it might only help me to calculate the variance of the random variable and not it's distribution.
Any help appreciated.